淺談結(jié)構(gòu)動(dòng)力有限元

1、得分：南京林業(yè)大學(xué)研克生課程論丈20132014學(xué)年第一學(xué)期課程號(hào)：63332課程名稱(chēng)：結(jié)構(gòu)動(dòng)力學(xué)論文題目：動(dòng)力有限元法在結(jié)構(gòu)工程中的應(yīng)用學(xué)科專(zhuān)業(yè)：結(jié)構(gòu)工程學(xué)號(hào)：3130194姓名：馬予煒任課教師：孫敦本二o三年十二月動(dòng)力有限元法在結(jié)構(gòu)工程中的應(yīng)用摘要：有限元法是口前應(yīng)用最為廣泛的一種離散化數(shù)值方法，其基木思想就是人 為地將連續(xù)體結(jié)構(gòu)分為有限個(gè)單元，規(guī)定每個(gè)單元所共有的一組變形形式，稱(chēng)之 為單元位移模式或插值函數(shù)。該方法在工程中冇著廣泛的應(yīng)用，比如：橋梁，建 筑上部和建筑基礎(chǔ)等。關(guān)鍵詞：有限元；動(dòng)力；位移；ansys1動(dòng)力有限元法簡(jiǎn)介冇限元方法在20世紀(jì)50年代起源于航空工程中飛機(jī)結(jié)構(gòu)的矩陣分

2、析。冇限 單元法發(fā)展至今，不斷開(kāi)拓新的領(lǐng)域，其范圍已經(jīng)ft!桿件結(jié)構(gòu)問(wèn)題擴(kuò)展到了彈性 力學(xué)乃至塑性力學(xué)問(wèn)題，由平面問(wèn)題擴(kuò)展到空間問(wèn)題，出靜力學(xué)問(wèn)題擴(kuò)展到動(dòng)力 學(xué)問(wèn)題和穩(wěn)定性問(wèn)題，由固體力學(xué)問(wèn)題擴(kuò)展到流體力學(xué)、熱力學(xué)和電磁學(xué)等問(wèn)題。結(jié)構(gòu)動(dòng)力學(xué)問(wèn)題的有限元法，實(shí)質(zhì)就是將一個(gè)彈性連續(xù)體的振動(dòng)問(wèn)題，離散 為一個(gè)以冇限個(gè)節(jié)點(diǎn)位移為廣義坐標(biāo)的多自曲度系統(tǒng)的振動(dòng)問(wèn)題。其基本原理和 分析方法類(lèi)同靜力學(xué)的冇限元法,按桿梁、薄板等不同結(jié)構(gòu)進(jìn)行分析。不同的是， 應(yīng)用振動(dòng)理論建立動(dòng)力學(xué)方程時(shí)，在單元分析中除需形成剛度矩陣外，還需形成 質(zhì)量矩陣，阻尼矩陣；在整體分析中，不僅求動(dòng)力響應(yīng)，還有求解特征值問(wèn)題（結(jié) 構(gòu)振動(dòng)的固

3、冇頻率及相應(yīng)的振動(dòng)型（或模態(tài)）。在動(dòng)力分析中，由于慣性力和阻 尼力岀現(xiàn)在平衡方程中，因此引入了質(zhì)量矩陣和阻尼矩陣，最后得到求解方程不 是代數(shù)方程組，而是常微分方程組。用有限元方法分析的步驟可概括為以下幾點(diǎn)：1單元剖分：把連續(xù)彈性體分割成許多個(gè)有限大小的單元，并為單元和節(jié)點(diǎn) 編號(hào)。2. 單元特征分析：對(duì)每個(gè)冇限單元，以節(jié)點(diǎn)位移為基本未知量，設(shè)選一 個(gè)單元位移函數(shù)，之后：（1）用節(jié)點(diǎn)位移表示單位位移，/ = "（2）通過(guò)幾何方程用節(jié)點(diǎn)位移表示單元應(yīng)變，£ = ；（3）通過(guò)物理方程用節(jié)點(diǎn)位移表示單元應(yīng)力，k = ga（4）通過(guò)虛功方程用節(jié)點(diǎn)位移表示節(jié)點(diǎn)力，fe =k e f得出單

4、元?jiǎng)?度矩陣。3. 總體結(jié)構(gòu)合成（1）分析整理各單元?jiǎng)偠染仃嚕ㄟ^(guò)節(jié)點(diǎn)的平衡方程形成節(jié)點(diǎn)載荷列陣、 合成總體剛度矩陣，建立以節(jié)點(diǎn)位移為未知量的、以總體剛度矩陣為系數(shù)的線性代數(shù)方程:ka = f°(2) 對(duì)線性代數(shù)方程組進(jìn)行邊界條件處理，求解節(jié)點(diǎn)位移。進(jìn)而由 b = gsf可求得單元應(yīng)力。4. 建立有限單元集合體的運(yùn)動(dòng)方程 mu+cu+ku = p(f)式屮，阻尼矩陣c可按rayleigh阻尼假設(shè)形成。2動(dòng)力有限元法在結(jié)構(gòu)工程中的應(yīng)用2.1問(wèn)題的描述某一框架一筒體結(jié)構(gòu)，總高度為54m, 18層，層高3m,結(jié)構(gòu)平面圖如圖1 所示。處于某一 7級(jí)近源地震荷載作用下，支坐所處的場(chǎng)地為iii類(lèi)

5、。一.八| 簡(jiǎn)體勰rsfc驟1外環(huán)梁1111111框架柱、171圖1結(jié)構(gòu)平面示意圖結(jié)構(gòu)主要承重構(gòu)件的截面尺寸及混凝土強(qiáng)度標(biāo)號(hào)如表1所示。在此模型 中，鋼筋混凝土的密度統(tǒng)一取2700kg/m3,彈性模量按混凝土的彈性模量取值, 泊松比取0.2。表1模型參數(shù)表構(gòu)件截面尺寸/m混凝土彈性模量e/mpa框架柱1.1x1.1c403.25el0外環(huán)梁0.4x0.6c403.25el0內(nèi)框架梁0.5 x 0.8c403.25el0次梁0.3x0.5c403.25el0筒體墻肢0.3c403.25el0樓屋面板0.2c303.0el0外圍墻體0.2c303.0el0豎向荷載僅考慮結(jié)構(gòu)自重，水平荷載僅考慮地震

6、作用。2.2模態(tài)分析模型建立完成后，進(jìn)入ansys后處理模塊p0st1進(jìn)行模態(tài)分析。表2為模型前10階振型自振頻率。前三階的振型圖如圖3圖5所示。 表2模型前十階口振頻率settime/freqload stepsubstepcumulative11.557311122.507712232.83941334585314455.835915566.041116676.345017786.466718896.56441991()6.881811()1()由模型的振型圖可以看出，結(jié)構(gòu)以水平振型為主，一階振型為y方向的平動(dòng)，二階振型為繞z軸的轉(zhuǎn)動(dòng)，三階振型也為繞z軸的轉(zhuǎn)動(dòng)。elementsurot、a

7、nsysr1a.5dec 27 201308:43:51圖3.階振羽圖4二階振型圖5三階振烈結(jié)論動(dòng)力學(xué)的有限元法同靜力學(xué)問(wèn)題，是把物體離散為有限個(gè)單元體，考慮單元的 慣性力和阻尼力等動(dòng)力因素的特性。動(dòng)力有限元法精度的一個(gè)較為全面的對(duì)比分 析應(yīng)包描動(dòng)力時(shí)程計(jì)算結(jié)果的分析。從自振頻率和精確頻率的對(duì)比分析可以發(fā)現(xiàn):(1) 有限元分析低階自振頻率的精度高于高階頻率。(2) 增加有限元數(shù)量可以很快提高分析精度。(3) 在采用數(shù)目相同的有限元法的情況下，一致質(zhì)量法的結(jié)果優(yōu)于集中質(zhì) 量法，但一致質(zhì)量法要花費(fèi)更多的時(shí)間來(lái)解決特征值問(wèn)題，因?yàn)閱卧臄?shù)量相同 時(shí)，一致質(zhì)量法的動(dòng)力口由度比集中質(zhì)量法的多一倍。動(dòng)力學(xué)的有限元法如果能求出各單元的物理特性，也就可以求出這個(gè)組合結(jié) 構(gòu)的物理特性。因此，有限元法就是對(duì)單元力學(xué)性質(zhì)的分析。它將在工程中得到 廣泛的應(yīng)用。參考文獻(xiàn)：1母恩爲(wèi)，陳國(guó)平，基于ansys的高層框架結(jié)構(gòu)地震響應(yīng)分析j,四川建筑第28卷6期王新敏,ansys工程結(jié)