有限元方法及國內(nèi)外研究現(xiàn)狀

現(xiàn)代機(jī)械設(shè)計(jì)理論與方法有限元方法學(xué)院：機(jī)械工程學(xué)院日期：2012年目錄摘要3關(guān)鍵詞3Abstract3KeyWords41有限元方法的國內(nèi)外研究現(xiàn)狀及應(yīng)用實(shí)例41.1有限元的開展趨勢(shì)41.2有限元的應(yīng)用實(shí)例42有限元方法的分析過程52.1有限元分析的三個(gè)階段52.2有限元分析的七個(gè)步驟52.3有限元軟件的分析過程63參考文獻(xiàn)8有限元方法摘要：有限元方法法的根本概念是用較簡(jiǎn)單的問題代替復(fù)雜問題后再求解。有限元法的根本思想是先化整為零﹑再積零為整，也就是把一個(gè)連續(xù)體分割成有限個(gè)單元；即把一個(gè)結(jié)構(gòu)看成由假設(shè)干通過節(jié)點(diǎn)相連的單元組成的整體，先進(jìn)行單元分析，然后再把這些單元組合起來代表原來的結(jié)構(gòu)進(jìn)行整體分析。關(guān)鍵詞：有限元方法；單元；節(jié)點(diǎn)FiniteElementMethodAbstract：Thebasicconceptsofthefiniteelementmethodissolvingcomplexproblemswithasimplequestioninstead.Thebasicideaofthefiniteelementmethodisdismembered,andthenplotthepartsintoawhole,thatisdividedacontinuumintoafinitenumberofunit;thatistoregardastructureasawholeconnectedbymanynodes，firsttoanalysisunit，thenanalysistheoverallcombinedbytheseunits，whichrepresentstheoriginalstructure.KeyWords：finiteelementmethod；unit；node1有限元方法的國內(nèi)外研究現(xiàn)狀及應(yīng)用實(shí)例“有限單元法”這一名稱是克拉夫〔Clough〕在1960年首先引用的。它是隨著電子計(jì)算機(jī)的開展而迅速開展起來的一種現(xiàn)代計(jì)算方法。它雖然是50年代首先在連續(xù)體力學(xué)領(lǐng)域—飛機(jī)結(jié)構(gòu)靜、動(dòng)態(tài)特性分析中應(yīng)用過的一種有效的數(shù)值分析方法，但是，由于它所依據(jù)理論的普遍性，已經(jīng)能夠成功地用來求解其它工程領(lǐng)域中的許多問題[1]。隨后很快廣泛的應(yīng)用于求應(yīng)力場(chǎng)、位移場(chǎng)、電磁場(chǎng)、溫度場(chǎng)、流體場(chǎng)等連續(xù)性問題。涉及了很多的工程學(xué)科，如機(jī)械設(shè)計(jì)、聲學(xué)、電磁學(xué)、巖土力學(xué)、流體力學(xué)等。在機(jī)械工程領(lǐng)域，有限元被廣泛的應(yīng)用于機(jī)構(gòu)、振動(dòng)和傳熱問題上。1.1有限元的開展趨勢(shì)縱觀當(dāng)今國際上有限元軟件的開展情況，可看出有限元軟件的一些開展趨勢(shì)：與CAD軟件的無縫集成；更為強(qiáng)大的網(wǎng)格處理能力；由求解線性問題開展到求解非線性問題；由單一結(jié)構(gòu)場(chǎng)求解開展到耦合場(chǎng)問題的求解；程序面向用戶開放性等。1.2有限元的應(yīng)用實(shí)例地鐵振動(dòng)預(yù)測(cè)的周期性有限元—邊界元耦合模型Jones利用有限元—無限元耦合二維模型計(jì)算了鐵路隧道內(nèi)及周圍土體的動(dòng)力響應(yīng)[2]。針對(duì)地鐵列車運(yùn)行引起的隧道結(jié)構(gòu)和自由場(chǎng)中的振動(dòng)響應(yīng)問題，提出了一個(gè)在頻率—波數(shù)域內(nèi)的三維周期性有限元—邊界元耦合的數(shù)值模型，此模型中隧道結(jié)構(gòu)采用有限元法計(jì)算，自由場(chǎng)采用邊界元法模擬。此模型采用Floquet變換，利用隧道和自由場(chǎng)在隧道軸線方向上的周期性，把無限長(zhǎng)的隧道及自由場(chǎng)的網(wǎng)格劃分限制在一個(gè)根本元內(nèi)，這樣使動(dòng)力學(xué)數(shù)值計(jì)算的效率大大提高。利用此模型計(jì)算了在隧道底板上施加固定單位諧振荷載情況下隧道—自由場(chǎng)相互作用系統(tǒng)的動(dòng)力響應(yīng)，結(jié)果說明此模型可應(yīng)用于地鐵列車運(yùn)行引起的隧道和自由場(chǎng)中的動(dòng)力響應(yīng)預(yù)測(cè)[3]。永磁直線同步電機(jī)解析分析及有限元驗(yàn)證對(duì)永磁直線同步電動(dòng)機(jī)具有邊端效應(yīng)、氣隙不均勻、氣隙磁場(chǎng)分布復(fù)雜、永久磁極幾何尺寸難以確定的問題，給出了單段永磁直線同步電動(dòng)機(jī)氣隙磁場(chǎng)的二維解析分析，討論了電機(jī)的各主要尺寸參數(shù)對(duì)氣隙磁場(chǎng)和性能的影響，提出了永久磁極的設(shè)計(jì)原那么和計(jì)算公式，分析了考慮飽和影響、計(jì)及邊端效應(yīng)時(shí)分段式永磁直線同步電動(dòng)機(jī)的非線性、不對(duì)稱、變化的自感和互感參數(shù)的解析計(jì)算方法，得到了定子繞組自感和互感系數(shù)隨動(dòng)子位置不同時(shí)的變化曲線。有限元數(shù)值計(jì)算結(jié)果說明，導(dǎo)出的分段式永磁直線同步電動(dòng)機(jī)永久磁極的設(shè)計(jì)計(jì)算公式以及電機(jī)變電感參數(shù)解析計(jì)算公式是適用的[4]?；陧憫?yīng)面的橋梁有限元模型修正采用試驗(yàn)設(shè)計(jì)和回歸分析方法，以顯式的響應(yīng)面模型逼近特征量與設(shè)計(jì)參數(shù)間復(fù)雜的隱式函數(shù)關(guān)系，得到簡(jiǎn)化的結(jié)構(gòu)模型(Meta-model)，給出有限元模型修正過程。針對(duì)復(fù)雜的土木工程結(jié)構(gòu)，討論樣本選擇、修正參數(shù)選取以及如何從眾多因素中較合理地建立結(jié)構(gòu)的響應(yīng)面模型。用數(shù)值模擬算例和六跨連續(xù)梁橋環(huán)境振動(dòng)試驗(yàn)結(jié)果，實(shí)現(xiàn)基于響應(yīng)面模型的土木工程結(jié)構(gòu)有限元模型修正，并與傳統(tǒng)的基于靈敏度方法直接對(duì)結(jié)構(gòu)有限元模型修正結(jié)果進(jìn)行比擬。結(jié)果說明，基于響應(yīng)面方法的有限元模型修正和驗(yàn)證，能顯著提高修正的效率，修正過程計(jì)算簡(jiǎn)潔、迭代收斂快，避開每次迭代都需要進(jìn)行有限元計(jì)算，易于工程實(shí)際應(yīng)用[5]。2有限元方法的分析過程2.1有限元分析的三個(gè)階段前處理階段：將整體結(jié)構(gòu)或其一局部簡(jiǎn)化為理想的數(shù)學(xué)力學(xué)模型，用離散化的單元代替連續(xù)實(shí)體結(jié)構(gòu)或求解區(qū)域；主要包括：定義分析的類型、添加材料屬性、添加載荷和約束、網(wǎng)格的劃分等……分析求解階段：運(yùn)用有限元法對(duì)結(jié)構(gòu)離散模型進(jìn)行分析計(jì)算，這個(gè)過程是由計(jì)算機(jī)來完成的。后處理階段：對(duì)計(jì)算結(jié)果進(jìn)行分析、整理和歸納。2.2有限元分析的七個(gè)步驟1、結(jié)構(gòu)力學(xué)模型的簡(jiǎn)化從實(shí)際的問題中抽象出力學(xué)模型，對(duì)實(shí)際問題的邊界條件、約束條件和外載荷條件進(jìn)行簡(jiǎn)化。抽象簡(jiǎn)化出來的力學(xué)模型應(yīng)該能盡可能的反映真實(shí)的實(shí)際問題，合理的模型既能保證計(jì)算結(jié)構(gòu)的精度，又不會(huì)帶來結(jié)構(gòu)上的過分復(fù)雜〔模型的建立在有限元分析的過程中是一個(gè)比擬重要的階段〕。2、結(jié)構(gòu)的離散化所謂結(jié)構(gòu)的離散化就是將連續(xù)的結(jié)構(gòu)體劃分為有限個(gè)單元體以代替原來的結(jié)構(gòu)。這個(gè)過程也就是網(wǎng)格的劃分，網(wǎng)格劃分是建立有限元模型的中心工作，模型的合理性很大程度上可以通過所劃分的網(wǎng)格形式反映出來。3、位移模式的選擇為了能用節(jié)點(diǎn)位移表示單元體的位移、應(yīng)變和應(yīng)力，在分析時(shí)就要對(duì)單元體位移的分布進(jìn)行一定的假設(shè)，假設(shè)位移是坐標(biāo)的某種簡(jiǎn)單函數(shù)，這個(gè)函數(shù)就是位移模式或者位移函數(shù)。一般情況下選擇多項(xiàng)式作為位移模式，因?yàn)樗械墓饣瘮?shù)的局部都可以用多項(xiàng)式逼近。4、分析單元的力學(xué)特性單元特性的分析包括以下三局部的內(nèi)容：〔1〕利用幾何方程，用位移模式導(dǎo)出用節(jié)點(diǎn)位移表示單元應(yīng)變的關(guān)系式?！?〕利用物理方程，由應(yīng)變表達(dá)式導(dǎo)出用節(jié)點(diǎn)位移表示單元應(yīng)力的關(guān)系式。〔3〕利用虛功原理建立作用于單元上的節(jié)點(diǎn)力和節(jié)點(diǎn)位移之間的關(guān)系式，也就是單元的剛度方程。5、計(jì)算等效節(jié)點(diǎn)力彈性體經(jīng)過離散后，假定力是通過節(jié)點(diǎn)在單元體之間進(jìn)行傳遞的。但是，實(shí)際的連續(xù)體，力是通過單元的公共邊界進(jìn)行傳遞的。因此，作用在單元上的各種力就需要等效移植到節(jié)點(diǎn)上去，也就是用等效的節(jié)點(diǎn)力來代替單元上的力。移植的方法是按照虛功等效原那么進(jìn)行的。6、集合所有單元的剛度方程，建立這個(gè)結(jié)構(gòu)的平衡方程這個(gè)過程包括兩個(gè)方面的內(nèi)容：一是由各個(gè)單元的剛度矩陣集合成整個(gè)物體的整體剛度矩陣；二是將作用于各單元的等效節(jié)點(diǎn)力列陣合成總的載荷列陣。集合所依據(jù)的原那么是要求相鄰的單元在公共節(jié)點(diǎn)處的位移相等。7、求解未知節(jié)點(diǎn)的位移和計(jì)算單元應(yīng)力由集合起來的平衡方程組，解出未知位移。在線彈性平衡問題中，可以根據(jù)方程組的具體特點(diǎn)選擇適宜的計(jì)算方法。對(duì)于非線性問題，那么通過一系列的步驟，并逐步修正剛度矩陣或載荷矩陣，才能獲得解答。最后，利用物理方程和求出的節(jié)點(diǎn)位移，計(jì)算各單元的應(yīng)力，并加以整理得出所要的結(jié)果。2.3有限元軟件的分析過程用Solidworks軟件中的simulation插件可以解決一些簡(jiǎn)單的有限元問題。數(shù)學(xué)模型的建立從實(shí)際問題中抽象出數(shù)學(xué)模型，是一個(gè)長(zhǎng)、寬、高分別為1000mm、100mm、50mm接下來，我們添加靜態(tài)算例，因?yàn)樗惺艿妮d荷是靜態(tài)的。添加約束和載荷將實(shí)體的左端面固定，并且在最右端添加一個(gè)豎直向下的5000N的力。如下列圖所示：將實(shí)體進(jìn)行網(wǎng)格劃分網(wǎng)格劃分后的結(jié)果如下：運(yùn)行求解應(yīng)力截圖位移截圖應(yīng)變截圖結(jié)果分析對(duì)結(jié)果進(jìn)行分析可知，與實(shí)際的工程應(yīng)用實(shí)例結(jié)果一致。3參考文獻(xiàn)[1]任中全，寇子明，趙燦等．現(xiàn)代機(jī)械設(shè)計(jì)理論與方法[M]．北京：煤炭工業(yè)出版社，2000：24-25.[2]JonesCJC，WangA，DawnTM．Modellingpropagationofvibrationfromrailwaytunnels[A]．ComputationalAcousticsanditsEnvironmentalApplications[c]．Southampton，April