幾何非線性問題的有限元法

1、1第八章幾何非線性問題的有限元法8.1引言前面各章所討論的問題都是在小變形假設的前提下進行的，即假定物體所發(fā)生的位移遠小于物體自身的幾何尺寸，應變遠小于1。在此前提下，建立物體或微元體的平衡條件時可以不考慮物體的位置和形狀（簡稱位形）的變化，因此在分析中不必區(qū)別變形前后位形的差別，且應變可用一階無窮小的線性應變表達。實際上，上述假設有時是不成立的。即使實際應變可能是小的，且不超過材料的彈性極限，但如果需要精確地確定位移，就必須考慮幾何非

2、線性，即平衡方程應該相對于變形后的位置得出，而幾何關系應該計及二次項。例如平板大撓度理論中，由于考慮了中面內(nèi)的薄膜應力，求得的撓度比小撓度理論的結果有顯著的減低。再如在結構穩(wěn)定性問題中，當載荷達到一定數(shù)值后，撓度比線性解答予示的結果更劇烈地增加，并且確實存在承載能力隨繼續(xù)變形而減低的現(xiàn)象。在冷卻塔、薄壁結構及其它比較細長的結構中，幾何非線性分析都顯得十分重要。幾何非線性問題可以分為以下幾種類型：（1）大位移小應變問題。一般工程結構所遇到

3、的幾何非線性問題大多屬于這一類。例如高層建筑或高聳構筑物以及大跨度網(wǎng)殼等結構的分析常需要考慮到結構大位移的影響。（2）大位移大應變問題，如金屬壓力加工中所遇到的問題就屬于這一類型。（3）結構的變形引起外載荷大小、方向或邊界支承條件的變化等。結構的平衡實際上是在結構發(fā)生變形之后達到的，對于幾何非線性問題來說，平衡方程必須建立在結構變形之后的狀態(tài)上。為了描述結構的變形需要設置一定的參考系統(tǒng)。一種做法是讓單元的局部坐標系始終固定在結構發(fā)生變形

4、之前的位置，以結構變形前的原始位形作為基本的參考位形，這種分析方法稱作總體的拉格朗日（Ｌagrange）列式法；另一種做法是讓單元的局部坐標系跟隨結構一起發(fā)生變位，分析過程中參考位形是不斷被更新的，這種分析方法稱作更新的拉格朗日列式法。本章首先對幾何非線性問題作一般性討論，從中導出經(jīng)典的線性屈曲問題的公式；然后建立平板大撓度問題和殼體的大位移（及大轉(zhuǎn)動）分析的有限方法公式；接著還給出了大應變及大位移的一般公式，最后還詳細討論了桿系結構幾

5、何非線性問題的有關公式。在討論中我們采用總體的拉格朗日列式法，但對桿系結構，為應用方便我們給出了兩種列式法的公式。8.2一般性討論8.2.1理論基礎無論是對于何種幾何非線性問題，虛功原理總是成立的。由虛功原理，單元的虛功方程可以寫成如下的形式（8.1）????????0???????veeTeeTFdv???其中為單元節(jié)點力向量，為單元的虛應變，為節(jié)點虛位移向量。??F??e????e??增量形式的應變一位移關系可表示為（8.2）???

6、???eedBd???3（8.10）??????????eeeeD00)(???????上式中和分別為單元材料中可能存在的初應變和初應力。??e0???e0?將式（8.2）代入式（8.9）就可以得到應力增量與單元節(jié)點位移增量之間的關系（8.11）????????eedBDd???將式（8.4）代入式（8.11）后得（8.12）??????????eLedBBDd??)(0??于是，式（8.8）左端中的第二項便可表示為（8.13）????

7、?????????????????????????????????????????????veLvTLTLvvLTTevTddvBDBdvBDBdvBDBdvBDBdvdB??))((0000若記(8.14)????????dvBDBkTv000????是與單元節(jié)點位移無關，它就是一般線性分析時的單元剛度矩陣。式（8.13）右端第??0k二層括號內(nèi)的項可記為（8.15）??????????????????????????vLTLTLLT

8、LdvBDBBDBBDBk)(00稱為單元的初位移矩陣或大位移矩陣，表示單元位置的變動對單元剛度矩陣的影響。??Lk現(xiàn)在再來看式（8.8）左端的第一項?？紤]到式（8.4）的關系并注意到與節(jié)點位移??0B無關，因此對節(jié)點位移的微分等于零，對于一個確定的有限元分析模型，式（8.8）左端的第一項可一般地寫成（8.16）????????????eevTLevTdkdvBddvBd????????????上式中稱為單元的初應力矩陣或幾何剛度矩陣，