第6章結(jié)構(gòu)動(dòng)力問(wèn)題的有限元法(2學(xué)時(shí))

1、有限元基礎(chǔ)第六章 結(jié)構(gòu)動(dòng)力問(wèn)題的有限元法劉艷芳Email: Add: 工程訓(xùn)練中心西配樓406Tel： 8233 8121劉艷芳，交通運(yùn)輸工程系， 第六章 結(jié)構(gòu)動(dòng)力問(wèn)題的有限元法章節(jié)內(nèi)容 介紹如何應(yīng)用有限元法求解結(jié)構(gòu)的動(dòng)力問(wèn)題？動(dòng)力學(xué)基本方程、質(zhì)量矩陣、阻尼矩陣；結(jié)構(gòu)的自由振動(dòng)和結(jié)構(gòu)的動(dòng)力響應(yīng)問(wèn)題；動(dòng)力學(xué)問(wèn)題的解法。章節(jié)目錄6.1 概述6.2 結(jié)構(gòu)動(dòng)力方程6.3 結(jié)構(gòu)的自由振動(dòng)劉艷芳，交通運(yùn)輸工程系， 6.1 概述n靜力分析假設(shè)載荷和結(jié)構(gòu)的響應(yīng)隨時(shí)間的變化非常緩慢。結(jié)構(gòu)的響應(yīng)：位移、應(yīng)變、應(yīng)力和力，與時(shí)間無(wú)關(guān)。靜力分析用于計(jì)算固定不變的載荷對(duì)結(jié)構(gòu)或部件的影響。n動(dòng)力分析考慮載荷和結(jié)構(gòu)的響應(yīng)隨時(shí)

2、間的變化而變化。結(jié)構(gòu)的響應(yīng)：位移、應(yīng)變、應(yīng)力和力是時(shí)間的函數(shù)；除此之外，還有速度、加速度。動(dòng)力分析用于求解隨時(shí)間變化的載荷對(duì)結(jié)構(gòu)或部件的影響。劉艷芳，交通運(yùn)輸工程系， 6.1 概述n工程意義結(jié)構(gòu)振動(dòng)是研究機(jī)械設(shè)備運(yùn)動(dòng)和力學(xué)問(wèn)題的重要基礎(chǔ)。振動(dòng)問(wèn)題引起的機(jī)械故障率高達(dá)6070，特別是運(yùn)動(dòng)機(jī)械。機(jī)械系統(tǒng)向高參數(shù)化發(fā)展，振動(dòng)和噪聲問(wèn)題日益突出。事故典型：n1940年，美國(guó)的Tacoma Narrows吊橋在中速風(fēng)載下，橋身發(fā)生嚴(yán)重的扭轉(zhuǎn)振動(dòng)和垂直振動(dòng)而導(dǎo)致崩塌。n1972年，日本海南電廠的一臺(tái)66萬(wàn)KW汽輪發(fā)電機(jī)組，在試車中因?yàn)榘l(fā)生異常振動(dòng)而全機(jī)毀壞，長(zhǎng)達(dá)51m的主軸斷裂飛射，聯(lián)軸節(jié)及汽輪機(jī)葉片竟然

3、穿透廠房飛落到百米之外。振動(dòng)分析和振動(dòng)設(shè)計(jì)成為產(chǎn)品設(shè)計(jì)中的一個(gè)關(guān)鍵環(huán)節(jié)。劉艷芳，交通運(yùn)輸工程系， 6.1 概述n動(dòng)力問(wèn)題的類型瞬態(tài)分析：用于分析結(jié)構(gòu)在隨時(shí)間變化的載荷作用下的動(dòng)力響應(yīng)模態(tài)分析：用于分析結(jié)構(gòu)的固有頻率和模態(tài)（振型）諧響應(yīng)分析：用于分析結(jié)構(gòu)在隨時(shí)間正弦變化的載荷作用下的響應(yīng)響應(yīng)譜分析：用于求解結(jié)構(gòu)在沖擊載荷下的響應(yīng)隨機(jī)振動(dòng)分析：用于分析結(jié)構(gòu)對(duì)隨機(jī)激勵(lì)的響應(yīng)顯式動(dòng)力分析：用于計(jì)算高度非線性動(dòng)力學(xué)和復(fù)雜的接觸問(wèn)題劉艷芳，交通運(yùn)輸工程系， 6.2 結(jié)構(gòu)動(dòng)力方程n單自由度振動(dòng)系統(tǒng) 設(shè)物體的質(zhì)量為m，阻尼系數(shù)為c，彈簧剛度為k，受到隨時(shí)間變化的激振力F(t)作用。 則其動(dòng)力方程為：)(tFk

4、zzczm 加速度速度位移劉艷芳，交通運(yùn)輸工程系， 6.2 結(jié)構(gòu)動(dòng)力方程舉例：?jiǎn)钨|(zhì)量彈簧阻尼振動(dòng)器zkzk劉艷芳，交通運(yùn)輸工程系， 6.2 結(jié)構(gòu)動(dòng)力方程舉例：?jiǎn)钨|(zhì)量彈簧阻尼振動(dòng)器解析解：劉艷芳，交通運(yùn)輸工程系， 6.2 結(jié)構(gòu)動(dòng)力方程n有限元?jiǎng)恿Ψ治龌静襟E單元離散（網(wǎng)格劃分）單元特性分析n單元位移模式n單元節(jié)點(diǎn)載荷n單元?jiǎng)恿ζ胶夥匠?.整體分析劉艷芳，交通運(yùn)輸工程系， 6.2 結(jié)構(gòu)動(dòng)力方程n單元位移模式假定位移模式為：?jiǎn)卧魏瘮?shù)矩陣N僅是點(diǎn)的坐標(biāo)的函數(shù)，與時(shí)間無(wú)關(guān)；單元節(jié)點(diǎn)位移e既是點(diǎn)的坐標(biāo)的函數(shù)，又是時(shí)間的函數(shù)；單元內(nèi)任一點(diǎn)的位移f既是點(diǎn)的坐標(biāo)的函數(shù)，又是時(shí)間的函數(shù)。單元的速度：?jiǎn)卧募铀俣?/p>

5、：顯然，單元的速度和加速度：既是點(diǎn)的坐標(biāo)的函數(shù)，又是時(shí)間的函數(shù)。劉艷芳，交通運(yùn)輸工程系， 6.2 結(jié)構(gòu)動(dòng)力方程n單元節(jié)點(diǎn)載荷載荷：來(lái)自于外界的激振力、慣性力和阻尼力n激振力：無(wú)論是集中載荷還是分布載荷，都按照靜力分析的虛功等效原則移置到單元的節(jié)點(diǎn)上。n單位體積的慣性力： 為材料的密度為材料的密度n單位體積的阻尼力： 為材料的阻尼系數(shù)為材料的阻尼系數(shù)n單元質(zhì)量矩陣和單元阻尼矩陣都是對(duì)稱矩陣。單元質(zhì)量矩陣單元阻尼矩陣劉艷芳，交通運(yùn)輸工程系， 6.2 結(jié)構(gòu)動(dòng)力方程n單元平衡方程 或者寫成：動(dòng)力分析與靜力分析的區(qū)別！單元?jiǎng)偠染仃嚭挽o力分析中的完全一樣。劉艷芳，交通運(yùn)輸工程系， 6.2 結(jié)構(gòu)動(dòng)力方程n整

6、體平衡方程整體質(zhì)量矩陣、整體阻尼矩陣和整體剛度矩陣都是對(duì)稱矩陣。節(jié)點(diǎn)載荷列陣節(jié)點(diǎn)位移列陣節(jié)點(diǎn)加速度列陣節(jié)點(diǎn)速度列陣整體剛度矩陣整體質(zhì)量矩陣整體阻尼矩陣劉艷芳，交通運(yùn)輸工程系， 6.2 結(jié)構(gòu)動(dòng)力方程n動(dòng)力方程的簡(jiǎn)化與靜力分析相比，結(jié)構(gòu)有限元?jiǎng)恿Ψ治龅挠?jì)算量非常大。簡(jiǎn)化途徑：n單元質(zhì)量矩陣（一致質(zhì)量矩陣）n舉例：平面3節(jié)點(diǎn)三角形單元桿單元梁?jiǎn)卧?.組裝后的整體質(zhì)量矩陣為帶狀對(duì)稱矩陣（稀疏矩陣）劉艷芳，交通運(yùn)輸工程系， 6.2 結(jié)構(gòu)動(dòng)力方程n動(dòng)力方程的簡(jiǎn)化與靜力分析相比，結(jié)構(gòu)有限元?jiǎng)恿Ψ治龅挠?jì)算量非常大。簡(jiǎn)化途徑：n單元質(zhì)量矩陣（一致質(zhì)量矩陣）n簡(jiǎn)化辦法：假設(shè)單元的質(zhì)量集中在各節(jié)點(diǎn)上。 因此，某一節(jié)

7、點(diǎn)的加速度不引起其它節(jié)點(diǎn)的慣性力。形成的單元質(zhì)量矩陣為對(duì)角矩陣，組裝后的整體質(zhì)量矩陣也是對(duì)角矩陣。1.簡(jiǎn)化后的質(zhì)量矩陣稱為 集中質(zhì)量矩陣。劉艷芳，交通運(yùn)輸工程系， 6.2 結(jié)構(gòu)動(dòng)力方程n動(dòng)力方程的簡(jiǎn)化與靜力分析相比，結(jié)構(gòu)有限元?jiǎng)恿Ψ治龅挠?jì)算量非常大。簡(jiǎn)化途徑：n阻尼矩陣2.當(dāng)密度和阻尼系數(shù)都為常數(shù)時(shí)，單元阻尼矩陣和單元質(zhì)量矩陣是成比例關(guān)系的。劉艷芳，交通運(yùn)輸工程系， 6.2 結(jié)構(gòu)動(dòng)力方程n動(dòng)力方程的簡(jiǎn)化與靜力分析相比，結(jié)構(gòu)有限元?jiǎng)恿Ψ治龅挠?jì)算量非常大。簡(jiǎn)化途徑：n阻尼矩陣n簡(jiǎn)化辦法(1) 用單元質(zhì)量矩陣m和單元?jiǎng)偠染仃噆的線性組合來(lái)表示單元阻尼矩陣。和與結(jié)構(gòu)的固有頻率和阻尼比有關(guān)。劉艷芳，交通

8、運(yùn)輸工程系， 6.2 結(jié)構(gòu)動(dòng)力方程n動(dòng)力方程的簡(jiǎn)化與靜力分析相比，結(jié)構(gòu)有限元?jiǎng)恿Ψ治龅挠?jì)算量非常大。簡(jiǎn)化途徑：n阻尼矩陣n簡(jiǎn)化辦法(2) 假設(shè)單元阻尼矩陣正比于單元?jiǎng)偠染仃?。劉艷芳，交通運(yùn)輸工程系， 6.3 結(jié)構(gòu)的自由振動(dòng)n彈性體振動(dòng)的有限元法系統(tǒng)的振動(dòng)特性：固有頻率和固有振型經(jīng)驗(yàn)證明，阻尼對(duì)結(jié)構(gòu)的固有頻率和振型影響很小，可以忽略。 因此，在沒有激振力的情況下，物體的動(dòng)力平衡方程（無(wú)阻尼振動(dòng)方程）為： 解析解：各節(jié)點(diǎn)的振幅向量（固有振型）固有頻率相位角劉艷芳，交通運(yùn)輸工程系， 6.3 結(jié)構(gòu)的自由振動(dòng)n無(wú)阻尼自由振動(dòng)方程各節(jié)點(diǎn)的振幅不可能全為零。 因此，可得結(jié)構(gòu)固有頻率的方程：一般解劉艷芳，交通運(yùn)輸工程系， 6.4 結(jié)構(gòu)的動(dòng)力響應(yīng)n結(jié)構(gòu)系統(tǒng)的動(dòng)力響應(yīng)根據(jù)結(jié)構(gòu)無(wú)阻尼自由振動(dòng)的頻率和振型，可以求解結(jié)構(gòu)的振動(dòng)方程，求出結(jié)構(gòu)各節(jié)點(diǎn)在強(qiáng)迫振動(dòng)時(shí)的瞬態(tài)位移，并進(jìn)而求得各單元在振動(dòng)時(shí)的瞬態(tài)應(yīng)力。n求解方法振型疊加法：適用于線性系統(tǒng)基本思想：將n階自由度系統(tǒng)的動(dòng)力方程，經(jīng)振型模態(tài)矩陣變換，化為互不耦合的n個(gè)單自由度問(wèn)題，進(jìn)行逐個(gè)求解后，再疊加得到動(dòng)力響應(yīng)的結(jié)果。逐步積分法：適用于非線性問(wèn)題的求解基本思想：把時(shí)間離散化，由初始狀態(tài)開始逐步求出每一個(gè)時(shí)間間隔上的狀態(tài)向量（位移、速度、加速度），最后求出的狀態(tài)向量就是結(jié)構(gòu)系統(tǒng)動(dòng)力響應(yīng)的解。 劉