復(fù)雜表面問題的有限元計算與分析

華沛

（河南城建學(xué)院，河南 平頂山 467000）

復(fù)雜表面問題的有限元計算與分析

華沛

（河南城建學(xué)院，河南 平頂山 467000）

本文主要通過對復(fù)雜接觸表面問題、流固耦合方程組以及邊界條件切觸不確定性的分析詳細探討了復(fù)雜表面問題的有限元計算的問題，希望能夠?qū)ξ覈蒲袦y試中的碰撞安全測試的最優(yōu)接觸問題的計算方法的研究有所幫助.

復(fù)雜接觸表面；流固耦合方程組；邊界條件切觸不確定性；有限元

近年來，隨著我國經(jīng)濟的高速發(fā)展，我國的在科學(xué)技術(shù)領(lǐng)域的研究越來越精確.有限元算法（FEA）在科研測試中也開始變得越來月重要.就目前來看，有限元計算法在接觸表面問題的數(shù)值計算中開始獲得較大的應(yīng)用，其中較為典型就是在三維立方元模型高性能前處理計算中應(yīng)用.三維立方元模型的計算以算法作為分類的標準，主要可以分為隱式以及顯式，其中隱式算法主要就是要對立方程組的聯(lián)解，在數(shù)學(xué)研究上來看，也即是應(yīng)該拉格朗日乘子法，且會在計算中常有難以收斂的情況發(fā)生，改重算法的優(yōu)點在于能夠?qū)Ρ砻娼佑|節(jié)點邊界條件有著精確的計算.對于顯式算法來說.在計算中有點是不需要對聯(lián)立方程求解，缺點在于受穩(wěn)定性限制.總之，就目前對復(fù)雜表面問題的有限元計算分析來看，立方元相關(guān)的計算是主要的研究方向，同時也是力學(xué)與經(jīng)濟學(xué)相關(guān)領(lǐng)域的隨機模型重要的參考對象.

1 復(fù)雜接觸表面問題

對于復(fù)雜接觸表面問題主要可以從兩個方面進行研究分析，也即是從確定性模型之上的隨機邊界條件分析以及蜂窩多孔復(fù)合材料的應(yīng)用和算法結(jié)果，下面對此進行簡單的分析探討：

1.1 確定性模型之上的隨機邊界條件分析

確定性模型之上的隨機邊界條件分析，主要可以分為三個步驟，第一步就是對于涉及時間變量的初邊值問題的分析，對與時間變量初邊值的處理基本就是依靠有限元方法求解的基本思想，也就是分開處理涉及到的時間變量以及與空間變量，使用有限元法的離散化對空間變量進行處理計算，從而得到與初值相關(guān)的常微分方程組.第二步驟，對第一步驟中有關(guān)的時間變量的處理，對于時間變量的處理是不同于空間變量的，時間變量的處理主要是應(yīng)用差分方法的離散化，從而得到線性代數(shù)方程組.第三步，利用確定性模型的工程試驗相關(guān)的數(shù)據(jù)庫，綜合第一步和第二步所得的時間和空間相關(guān)變量的初邊值可以確定離散系統(tǒng)隨機初邊值條件，進而對該離散系統(tǒng)進行數(shù)學(xué)化的運算，從而能夠取得相應(yīng)研究對象，比如像應(yīng)力、塑性變形率以及常見的碰撞加速度控制等等相關(guān)的各種分布函數(shù)的高精度解.圖1即為碰撞研究對象中應(yīng)力（VMS）與碰撞角度（IA）的聯(lián)合分布頻率函數(shù)的圖解.

圖1 應(yīng)力（VMS）與碰撞角度（IA）的聯(lián)合分布頻率隨時間的變化

1.2 蜂窩多孔復(fù)合材料的應(yīng)用和算法結(jié)果

關(guān)于有限元計算方法在蜂窩多孔復(fù)合材料（如圖2所示的力學(xué)結(jié)構(gòu)圖）這種典型復(fù)雜表面問題中應(yīng)用計算來說，重要的問題是對動力方程中相關(guān)的多尺度算法的計算規(guī)模的考慮，主要方法那就是在蜂窩多孔材料區(qū)域耦合方程中，關(guān)于對時間變量的處理時，加強對積分變換的應(yīng)用，同時對帶有參數(shù)關(guān)于蜂窩多孔材料的新耦合方程進行合理的多尺度分析，以求得簡單易于計算多尺度算法.不僅如此，在對蜂窩夾層結(jié)構(gòu)的表面問題進行計算時，可以應(yīng)用三維網(wǎng)格進行數(shù)值進行計算，以達到對蜂窩多孔復(fù)合材料表面問題對有限元的應(yīng)用，從而求出結(jié)果的目地.

圖2 蜂窩多孔復(fù)合材料區(qū)域上耦合彈塑性方程的算法結(jié)果

2 流固耦合方程組

在復(fù)雜表面問題的有限元計算中，流固耦合方程組也是整個計算中重要的環(huán)節(jié)，流固耦合方程組是對復(fù)雜表面性質(zhì)的描述的基礎(chǔ)，也是計算的方法之一，下面就從流固耦合方程組的離散系統(tǒng)和復(fù)雜接觸表面的計算方法這兩方面進行簡單分析：

2.1 流固耦合方程組

流固耦合方程組對于復(fù)雜接觸表面問題的來說，主要是起到媒介的作用，也就是可以通過流固耦合方程組進行描述的同時，還能夠通過流固耦合方程組進行計算，換從研究的角度來說，主要是表現(xiàn)在下面的兩個方面：在工程材料科學(xué)中的相關(guān)標準變量的基礎(chǔ)之上，以對非牛頓P－T／T方程所表現(xiàn)的處理特性為切入點，分析復(fù)雜接觸表面的外延以及簡單剪切速率的阻力；以Cauchy守恒方程為基礎(chǔ)，通過對方程的計算可以得出，由于宏觀與微觀應(yīng)力場中因素τ分布的變化，產(chǎn)生的彈塑性材料大變形的基本情況，也即是研究中通常所說的切變致稀），其中關(guān)于移動邊界層在時間點數(shù)值結(jié)果如圖3所示.

圖3 三維碰撞模型中接觸表面移動邊界層，部分時間節(jié)點上的結(jié)果

2.2 復(fù)雜接觸表面的計算方法

通過上文對復(fù)雜接觸表面的基本情況以及性質(zhì)的分析可以得知，在對復(fù)雜表面問題處理是可以以非牛頓耦合方程模型為基礎(chǔ)進行探究，具體是可以通過非牛頓耦合方程模型隨機邊界層的特點來分析研究，在此研究為基礎(chǔ)，對相關(guān)研究的任意元素的計算進行數(shù)值仿真.具體的來說也就是應(yīng)用人工邊界法數(shù)值來求解非線性外問題，從而達到對復(fù)雜接觸表面年的計算.但是，該種計算方法應(yīng)該注意的是，人工邊界法數(shù)值求解非線性外問題，其基礎(chǔ)是在得到人工邊界上相對于該非線性問題相關(guān)的準確近似的邊界條件.該邊界條件一般來說就是未知函數(shù)及其微商相對于與計算相關(guān)的人工邊界的非線性積分方程.

3 邊界條件切觸不確定性

對于復(fù)雜接觸表面處理的有限元的計算與分析中，邊界條件切觸的不確定性也是整個計算中的重點難點，下面就從邊界條件切觸不確定性的邊界層特征函數(shù)和邊界層的外解漸進分析以及內(nèi)解漸進分析進行闡述：

邊界層的特征函數(shù).在復(fù)雜表面問題中利用邊界層理論對有限元接觸邊界相關(guān)的非線性特征函數(shù)基進行探討，這是對復(fù)雜表面問題進行有限元計算的關(guān)鍵點之一.簡單的來說也即是以邊界層外解與內(nèi)解逼近為基礎(chǔ)，通過對非線性邊界條件的匹配方法.從而得到跳躍條件特征值譜P與參數(shù)R，G，N及χ的關(guān)系，也即是得到了復(fù)雜表面有限元后處理的優(yōu)化基函數(shù)

外解的漸近分析.邊界條件切觸性對于漸近分析，主要就是針對對外解問題微分方程（方程一）而言，通過該方程的分析，可得到相關(guān)的彈性方程以及耦合方程組在復(fù)雜接觸表面層厚中當δ→0時，將會喪失方程或方程組的橢圓性，這也即是要求構(gòu)造新的空間，去討論光滑解相關(guān)的收斂性.因此，就外解的漸近分析來說，邊界層外穩(wěn)態(tài)約束場是好的逼近.方程如下：

內(nèi)解的漸近分析.對于內(nèi)解的漸近分析，主要是對邊界層內(nèi)解的進行計算，也即是通過對非線性相關(guān)的耦合程進行線性變換，同時對傅立葉進行調(diào)和分析，在進新的參數(shù)的情況下，將邊界相關(guān)的三維流體動力學(xué)方程轉(zhuǎn)換為一維復(fù)空間的階微分方程，方程如下：

通過對方程的分析與計算，從而可以確定復(fù)雜接觸表面的非線性有限元特征函數(shù)基，進而能夠保證在復(fù)雜表面有限元相關(guān)計算中的的穩(wěn)定與收斂.

總而言之，在對復(fù)雜表面問題的有限元計算與分析時，計算的方法和基礎(chǔ)就在于去采用流固耦合非牛頓流體方程，使用該方程對復(fù)雜表面的初邊值求解三維層結(jié)構(gòu)特性，同時可以利用變分原理的攝動問題的有限元方法，在已經(jīng)研發(fā)使用的高性能軟件平臺對所得數(shù)據(jù)實現(xiàn)挖掘處理，進而得到有效的結(jié)果.此外，復(fù)雜接觸表面的有限元計算，還可以利用基于漸近攝動方法，已經(jīng)較為成熟邊界層理論進行研究，從而得到可解的微分方程特征的函數(shù)空間，這也可以用來優(yōu)化有限元基函數(shù)的解，達到對復(fù)雜表面問題有限元計算的目地.

〔1〕林群，周俊明，陳竤燾.橢圓形方程四面體線元的超逼近與外推[J].數(shù)學(xué)實踐與認識，2009，39（15）：200－208.

〔2〕林群，謝和虎，羅福生，等.Stokes方程非協(xié)調(diào)混合元的特征值下界[J].數(shù)學(xué)的實踐與認識，2010，40（19）：157－168.

〔3〕李開泰，黃愛香，黃慶懷.有限元方法及其應(yīng)用[M].北京：科學(xué)出版社，2005.

O241.82

A

1673-260X（2014）03-0006-02