有限元分析與應(yīng)用第一章

有限元分析與應(yīng)用第一章第一頁(yè)，共二十九頁(yè)，2022年，8月28日第二頁(yè)，共二十九頁(yè)，2022年，8月28日（1）工程問(wèn)題（2）數(shù)值方法（3）有限元方法與ANSYS簡(jiǎn)介（4）有限元方法的基本步驟（5）直接公式法（6）最小總勢(shì)能公式（7）加權(quán)余數(shù)法（8）結(jié)果的驗(yàn)證（9）理解問(wèn)題本章的內(nèi)容第三頁(yè)，共二十九頁(yè)，2022年，8月28日

工程問(wèn)題一般是物理問(wèn)題的數(shù)學(xué)模型。數(shù)學(xué)模型是帶有邊界條件和初始問(wèn)題的微分方程，微分方程是通過(guò)對(duì)系統(tǒng)應(yīng)用自然的基本定律和原理推導(dǎo)出來(lái)的，如波動(dòng)方程等。這些微分方程代表了某種物理規(guī)律或平衡。（數(shù)理方程）在可能的情況下，由給定的條件求解微分方程可以得到系統(tǒng)的精確行為。在任何給定的工程問(wèn)題中，存在兩種影響系統(tǒng)行為的參數(shù)。一種表示給定系統(tǒng)自然行為的參數(shù)，例如彈性模量、熱傳導(dǎo)因子和粘度等。另一種是系統(tǒng)存在產(chǎn)生擾動(dòng)的參數(shù)，如外力、力矩、介質(zhì)的溫度差和流體的壓力差等。在有限元建模中，理解參數(shù)在剛度或傳導(dǎo)矩陣以及負(fù)荷矩陣中的作用是非常重要的。系統(tǒng)特性總是在剛度矩陣、傳導(dǎo)矩陣中得到體現(xiàn)，而擾動(dòng)參數(shù)總是出現(xiàn)在負(fù)荷矩陣中。工程問(wèn)題第四頁(yè)，共二十九頁(yè)，2022年，8月28日在許多實(shí)際工程問(wèn)題中，由于微分方程的復(fù)雜性或邊界條件和初始條件的難以確定性，得不到系統(tǒng)的精確解。為此我們借助數(shù)值方法來(lái)近似。解析解表明了系統(tǒng)在任何點(diǎn)上的精確行為，而數(shù)值解只在稱為節(jié)點(diǎn)的離散點(diǎn)上近似于解析解。數(shù)值法的第一步都是離散化。這一過(guò)程將系統(tǒng)分成一些單元和節(jié)點(diǎn)，然后對(duì)每一單元或節(jié)點(diǎn)建立代數(shù)方程組。這種方法假設(shè)代表每個(gè)單元的近似函數(shù)是連續(xù)的。假設(shè)單元間的邊界是連續(xù)的，通過(guò)組合各單元的解產(chǎn)生系統(tǒng)的完全解。數(shù)值方法第五頁(yè)，共二十九頁(yè)，2022年，8月28日

有限元方法是用于求解工程中各類問(wèn)題的數(shù)值方法，應(yīng)力分析中穩(wěn)態(tài)的、瞬態(tài)的、線性的或非線性的問(wèn)題以及熱傳導(dǎo)、流體流動(dòng)和電磁學(xué)中的問(wèn)題都可以用有限元方法進(jìn)行分析解決。現(xiàn)代有限元方法的20世紀(jì)早期開(kāi)始，20世紀(jì)50年代，boeing公司采用三角元對(duì)機(jī)翼進(jìn)行建模，推動(dòng)了有限元方法的應(yīng)用。到20世紀(jì)60年代，人們接受了“有限元”這個(gè)詞。ANSYS是一個(gè)通用的有限元計(jì)算機(jī)程序，其代碼長(zhǎng)度超過(guò)10萬(wàn)行。應(yīng)用ANSYS可以進(jìn)行靜態(tài)、動(dòng)態(tài)、熱傳導(dǎo)、流體流動(dòng)和電磁學(xué)等分析。在過(guò)去的20多年里，ANSYS是主要的有限元分析程序?，F(xiàn)在ANSYS被廣泛應(yīng)用在如航天、汽車、電子、核科學(xué)等領(lǐng)域。有限元方法與ANSYS簡(jiǎn)介第六頁(yè)，共二十九頁(yè)，2022年，8月28日有限元方法的基本步驟

預(yù)處理階段1建立求解域并將之離散化成有限元，即將問(wèn)題分解成節(jié)點(diǎn)和單元。2假設(shè)代表單元物理行為的形函數(shù)，即假設(shè)代表單元解的近似連續(xù)函數(shù)。3對(duì)單元建立方程。4將單元組合成總體的問(wèn)題，構(gòu)造總體剛度矩陣。5應(yīng)用邊界條件、初值條件和負(fù)荷。求解階段6求解線性或非線性的方程組，以得到節(jié)點(diǎn)的值，例如得到不同節(jié)點(diǎn)的位移量或熱傳遞問(wèn)題中不同節(jié)點(diǎn)的溫度值。后處理階段7得到其他重要的信息。如主應(yīng)力、熱量值等。第七頁(yè)，共二十九頁(yè)，2022年，8月28日1.5直接公式法例1.1桿的一端固定，另一端承受負(fù)荷P，桿的厚度t，長(zhǎng)度L。桿的彈性模量用E表示。假設(shè)應(yīng)用的負(fù)荷比桿的重量大的多。將問(wèn)題域離散成有限的單元，首先將問(wèn)題分解成節(jié)點(diǎn)和單元。我們用5節(jié)點(diǎn)4單元的模型代替桿，如圖所示，分的越細(xì)越精確。桿的模型中有四個(gè)獨(dú)立的分段，每個(gè)分段都有統(tǒng)一的橫截面積。預(yù)處理階段第八頁(yè)，共二十九頁(yè)，2022年，8月28日第九頁(yè)，共二十九頁(yè)，2022年，8月28日假設(shè)近似單元行為的近似解考慮一個(gè)帶有統(tǒng)一橫截面A的實(shí)體的偏轉(zhuǎn)量，橫截面的長(zhǎng)度為l，承受的外力為F，如圖所示。

實(shí)體的平均應(yīng)力為實(shí)體的平均應(yīng)變?yōu)閼?yīng)變與應(yīng)力的關(guān)系為這里E是彈性模量。第十頁(yè)，共二十九頁(yè)，2022年，8月28日這樣有

令

，則

這相當(dāng)于彈簧。我們認(rèn)為桿由四個(gè)彈性系數(shù)不同的彈簧組成，對(duì)每個(gè)彈簧來(lái)說(shuō)，

和分別是和處的節(jié)點(diǎn)的橫截面積，l是單元的長(zhǎng)度。第十一頁(yè)，共二十九頁(yè)，2022年，8月28日靜力學(xué)要求每個(gè)節(jié)點(diǎn)上的應(yīng)力和為零，所以節(jié)點(diǎn)1：節(jié)點(diǎn)2：節(jié)點(diǎn)3：節(jié)點(diǎn)4：節(jié)點(diǎn)5：第十二頁(yè)，共二十九頁(yè)，2022年，8月28日重組方程組，得到

第十三頁(yè)，共二十九頁(yè)，2022年，8月28日寫成矩陣的形式為=在負(fù)荷矩陣中，將反作用力和負(fù)荷區(qū)分開(kāi)來(lái)是很重要的，于是=第十四頁(yè)，共二十九頁(yè)，2022年，8月28日上式表示成

上式中，R為反作用力，K為剛度矩陣，u為位移矩陣，F(xiàn)為負(fù)荷矩陣。由于桿的頂端是固定的，其位移為零，上面的方程簡(jiǎn)化為求解上面的方程，得到節(jié)點(diǎn)的位移量。第十五頁(yè)，共二十九頁(yè)，2022年，8月28日3．對(duì)單元建立方程每個(gè)單元有兩個(gè)節(jié)點(diǎn)，每個(gè)節(jié)點(diǎn)相應(yīng)一個(gè)位移量，因此我們需要對(duì)每個(gè)單元建立兩個(gè)方程，這些方程必須和節(jié)點(diǎn)的位移量和單元的剛度有關(guān)，考慮單元內(nèi)部傳遞的力的關(guān)系。第十六頁(yè)，共二十九頁(yè)，2022年，8月28日節(jié)點(diǎn)滿足的方程為寫成矩陣形式

第十七頁(yè)，共二十九頁(yè)，2022年，8月28日4.將單元組合起來(lái)表示整體問(wèn)題單元1的剛度矩陣為它在總體剛度矩陣中的位置為類似的，對(duì)于單元2、3、4，我們有第十八頁(yè)，共二十九頁(yè)，2022年，8月28日第十九頁(yè)，共二十九頁(yè)，2022年，8月28日總體剛度矩陣為利用單元建立方程和利用節(jié)點(diǎn)建立方程，得到的結(jié)果一樣。第二十頁(yè)，共二十九頁(yè)，2022年，8月28日5.應(yīng)用邊界條件和負(fù)荷桿的頂端是固定的，即有邊界條件。在節(jié)點(diǎn)5處加外力P，于是矩陣中的第一行必須包含一個(gè)1和四個(gè)0，以讀取給定的邊界條件在固體力學(xué)的問(wèn)題中，有限元公式一般有如下的形式：

第二十一頁(yè)，共二十九頁(yè)，2022年，8月28日6.求解階段我們以例題來(lái)說(shuō)明它。假設(shè)鋁975845715585單元節(jié)點(diǎn)平均截面面積（in2）長(zhǎng)度（in）彈性模量單元?jiǎng)偠认禂?shù)1120.2343752.52230.2031252.53340.1718752.54450.1406252.5第二十二頁(yè)，共二十九頁(yè)，2022年，8月28日桿的橫截面面積的變化可以由下式來(lái)表示這樣，每個(gè)橫截面上的面積為每個(gè)單元的對(duì)等剛度系數(shù)為：第二十三頁(yè)，共二十九頁(yè)，2022年，8月28日單元矩陣為：第二十四頁(yè)，共二十九頁(yè)，2022年，8月28日總體剛度矩陣為應(yīng)用邊界條件和負(fù)荷，得到第二十五頁(yè)，共二十九頁(yè)，2022年，8月28日第二行中，系數(shù)－975乘以的結(jié)果為零，我們得到如下矩