非線性有限元法(5)-7

4方程解法4.1顯示方法（ExplicitMethods）在計(jì)算力學(xué)和計(jì)算物理中應(yīng)用最廣泛的顯式方法是中心差分法定義時(shí)間步增量為

1

4方程解法速度的中心差分公式為

積分公式為加速度及積分公式

2

4方程解法在第n時(shí)間步的有限元半離散方程為廣義約束方程代半離散方程到加速度積分公式中這樣，在任意時(shí)間步n，已知位移，通過順序計(jì)算應(yīng)變-位移方程、由或形式表示的本構(gòu)方程和外力，可確定節(jié)點(diǎn)力，上式右端全部賦值，可獲得，由位移的積分公式，可確定。

3

4方程解法特點(diǎn)a)M矩陣是常數(shù)陣，在整個(gè)求解過程中不變;b)如果M為對角陣，實(shí)現(xiàn)節(jié)點(diǎn)速度和節(jié)點(diǎn)位移的更新可以不用求解任何方程；c)因此，在顯示方法中，離散的動量方程一般都用團(tuán)聚質(zhì)量矩陣?！陲@示方法中，對離散的動量方程的時(shí)間積分不需求解任何方程

4

4方程解法流程圖

5

4方程解法流程圖（續(xù)）

6

4方程解法優(yōu)點(diǎn)－－方法簡單，不求解方程

顯式時(shí)間積分是很強(qiáng)健的（Robust），很少因?yàn)閿?shù)值失敗而終止缺點(diǎn)－－

顯式方法是條件穩(wěn)定的

如果時(shí)間步長超過一個(gè)臨界值，其計(jì)算結(jié)果會增長至無窮臨界時(shí)間步長

對于采用率無關(guān)材料的常應(yīng)變單元的網(wǎng)格，穩(wěn)定時(shí)間步長為

線性系統(tǒng)的最大頻率，單元e的特征長度，單元e的當(dāng)前波速，是折減系數(shù)。

7

4方程解法有限差分的臨界步長問題稱為Courant條件1928年Courant，F(xiàn)riedrichs和Lewy首先發(fā)表這一結(jié)果。臨界步長隨網(wǎng)格細(xì)劃和材料剛度的增加而減小對于彈－塑性問題，塑性響應(yīng)會減慢波速，但并不能增加臨界時(shí)間步長，因?yàn)?，在?shù)值計(jì)算中常常發(fā)生卸載從單元時(shí)間步長確定網(wǎng)格時(shí)間步長

8

4方程解法4.2平衡解法和隱式時(shí)間積分（Implicittimeintegration）離散的動量方程可寫為Newmark方程

更新的位移和速度為

的意義為人工粘性或人工阻尼。

9

4方程解法

10

4方程解法更新的加速度表示為這樣離散問題成為：這是在節(jié)點(diǎn)位移的非線性方程組。問題成為

尋找使服從于

11

4方程解法

Newton－Raphson法首先研究一個(gè)未知量的問題Newton方法是一個(gè)求解迭代過程。首先將殘數(shù)進(jìn)行Taylor展開，只保留線性項(xiàng)稱為非線性方程的線性化模型。線性模型是非線性殘差函數(shù)得正切。12

4方程解法

對于位移增量，求解線性模型，得更新值為持續(xù)迭代直到獲得理想得精度。13

4方程解法

對應(yīng)的載荷位移曲線14

4方程解法

有n個(gè)未知量的Newton法

線性化后方程為A稱為Jacobian矩陣，在計(jì)算力學(xué)中A稱為切線剛度矩陣。求解步驟：（1）解線性方程組（2）迭加到前一步得迭代（3）檢驗(yàn)其收斂性，若收斂得到解；若沒有滿足準(zhǔn)則，重新構(gòu)造A，

返回（1）15

4方程解法

16

4方程解法

17

4方程解法

收斂準(zhǔn)則

Newton法的迭代是否終止是由收斂準(zhǔn)則決定的。這是關(guān)于的收斂。一般有三種收斂準(zhǔn)則：（1）根據(jù)參數(shù)r的量級的準(zhǔn)則（2）根據(jù)位移增量的量級準(zhǔn)則（3）能量誤差的準(zhǔn)則18

4方程解法

19初應(yīng)力法

4方程解法

20修正牛頓法

4方程解法

214.3Wilsonq

4方程解法

22

4方程解法

23K.J.Bathe,FiniteElementProcedures,Prentice-Hall,Inc.,1996

4方程解法

24

4方程解法

4.4弧長法（Arc-lengthmethod）

25

4方程解法

單自由度示意圖

26

4方程解法

在跟蹤結(jié)構(gòu)得后屈曲（post-buckling）曲線時(shí)，或跟蹤分支時(shí)，Newton-Raphson法出現(xiàn)困難，弧長法是一個(gè)好的替代解法。在跟蹤分支時(shí)，載荷參數(shù)從零開始，并逐步增加。對于參數(shù)的每一個(gè)增量，計(jì)算平衡解答，即

在弧長法中，除了逐漸增加的載荷參數(shù)外，在位移-載荷參數(shù)空間中的弧長度量也要增加，這是通過在平衡方程中增量參數(shù)化方程實(shí)現(xiàn)的：27

4方程解法

方程組中增加一個(gè)未知數(shù)，并增加了一個(gè)參數(shù)化方程：使用Newton線性化方程組28

4方程解法

29計(jì)算后屈曲的例子

4方程解法