新手!什么叫顯示動力學,什么叫隱式動力學分析

1、新手！什么叫顯示動力學，什么叫隱式動力學分析我是一個ANSYS新手，最近寫論文需要仿真封口機封罐的過程，于是就開始學起了LS-dyna,書上老說顯示動力學分析，我不理解其義，于是問一下什么是顯示分析，什么是隱式分析？ yuminhust2005 發(fā)表于 2009-10-16 11:39i=s 本帖最后由 yuminhust2005 于 2009-10-16 11:41 編輯 /i"-SA!Ui9k;x'Rw 說得簡單點，1l:x0t xc0q顯式分析，下一步的計算結(jié)果只和前面的計算結(jié)果有關。有條件收斂，要求時間步較小。通常做動力分析用這種方法。,r5F1LMI隱式分析，下一步

2、的計算結(jié)果不僅和前面的結(jié)果有關，而且和下一步的結(jié)果有關，通過迭代得到。無條件收斂。通常做靜力分析用這種方法。 sammy388 發(fā)表于 2010-4-10 19:01版主總結(jié)的極是，在此再補充一點color=red（以下內(nèi)容全部引自中國機械CAD論壇，貼出來幫助新成員以及對此問題迷惑的人）/color_I&N,ql Z:_*u2f 1、顯式算法基于動力學方程，因此無需迭代；而靜態(tài)隱式算法基于虛功原理，一般需要迭代計算*y+vk4eK2S2、顯式算法t6A0iS/Z$c5    最大優(yōu)點是有較好的穩(wěn)定性。gB9U:Rh2d     動態(tài)顯式算法采用

3、動力學方程的一些差分格式（如廣泛使用的中心差分法、線性加速度法、Newmark法和wilson法等），不用直接求解切線剛度，不需要進行平衡迭代，計算速度快，時間步長只要取的足夠小，一般不存在收斂性問題。因此需要的內(nèi)存也比隱式算法要少。并且數(shù)值計算過程可以很容易地進行并行計算，程序編制也相對簡單。但顯式算法要求質(zhì)量矩陣為對角矩陣，而且只有在單元級計算盡可能少時速度優(yōu)勢才能發(fā)揮, 因而往往采用減縮積分方法，容易激發(fā)沙漏模式，影響應力和應變的計算精度。(m&_qXv.jZ     靜態(tài)顯式法基于率形式的平衡方程組與Euler向前差分法，不需要迭代求解。由于平衡方程式僅在率

4、形式上得到滿足，所以得出的結(jié)果會慢慢偏離正確值。為了減少相關誤差，必須每步使用很小的增量。eP&MV;tLu 3、隱式算法 )LdD B"Qv;aJ P     隱式算法中，在每一增量步內(nèi)都需要對靜態(tài)平衡方程進行迭代求解，并且每次迭代都需要求解大型的線性方程組，這以過程需要占用相當數(shù)量的計算資源、磁盤空間和內(nèi)存。該算法中的增量步可以比較大，至少可以比顯式算法大得多，但是實際運算中上要受到迭代次數(shù)及非線性程度的限制，需要取一個合理值。6 _;Cc4、求解時間eJY C%i$y     使用顯式方法，計算成本消耗與單元數(shù)量成正比，并且大致與

5、最小單元的尺寸成反比ka'h*s0Ae:z    應用隱式方法，經(jīng)驗表明對于許多問題的計算成本大致與自由度數(shù)目的平方成正比1vN-?9T0:|     因此如果網(wǎng)格是相對均勻的，隨著模型尺寸的增長，顯式方法表明比隱式方法更加節(jié)省計算成本 *B.b:ac KV8H 隱式求解法*BvF7UP9S    將沖壓成型過程的計算作為動態(tài)問題來處理后，就涉及到時間域的數(shù)值積分方法問題。在80年代中期以前，人們基本上使用牛曼法進行時間域的積分。根據(jù)牛曼法，位移、速度和加速度有著如下的關系：上面式子中 ， 分別為當前時刻和前一時刻的位移， 和

6、 為當前時刻和前一時刻的速度， 和 為當前時刻和前一時刻的加速度，和為兩個待定參數(shù)。由上式可知，在牛曼法中任一時刻的位移、速度和加速度都相互關聯(lián)，這就使得運動方程的求解變成一系列相互關聯(lián)的非線性方程的求解。這個求解過程必須通過迭代和求解聯(lián)立方程組才能實現(xiàn)。這就是通常所說的隱式求解法。隱式求解法可能遇到兩個問題。一是迭代過程不一定收斂；二是聯(lián)立方程組可能出現(xiàn)病態(tài)而無確定的解。隱式求解法的最大優(yōu)點是它具有無條件穩(wěn)定性，即時間步長可以任意大。(RS0yXI 顯式求解法b:FGMj6F    如果采用中心差分法來進行動態(tài)問題的時域積分，則有如下位移、速度和加速度關系b D5U4LB

7、  由上式可以看出，當前時刻的位移只與前一時刻的加速度和位移有關，這就意味著當前時刻的位移求解無需迭代過程。另外，只要將運動方程中的質(zhì)量矩陣和阻尼矩陣對角化，前一時刻的加速度求解無需解聯(lián)立方程組，從而使問題大大簡化，這就是所謂的顯式求解法。顯式求解法的優(yōu)點是它即沒有收斂性問題，也不需求解聯(lián)立方程組，其缺點是時間步長受到數(shù)值積分穩(wěn)定性的限制，不能超過系統(tǒng)的臨界時間步長。由于沖壓成型過程具有很強的非線性，從解的精度考慮，時間步長也不能太大，這就在很大程度上彌補了顯式求解法的缺陷。Pmqg.F%7uq5HE     在80年代中期以前顯式算法主要用于高速碰撞的

8、仿真計算，效果很好。自80年代后期被越來越廣泛地用于沖壓成型過程的仿真，目前在這方面的應用效果已超過隱式算法。顯式算法在沖壓成型過程的仿真中獲得成功應用的關鍵，在于它不像隱式算法那樣有解的收斂性問題。1w7|4d#Y$E nJ7e 顯式算法和隱式算法，有時也稱為顯式解法和隱式解法，是計算力學中常見的兩個概念，但是它們并沒有普遍認可的定義，下面只是我的一些理解。先看看一般對兩種方法的理解和比較，2e6N3awv9E=g+i-Zp                    

9、   顯式算法            隱式算法M%b zfs-.(!Dv3WV y   (01)適用問題        動力學（動態(tài)）      靜力學（靜態(tài)）k wx4J XDyP   (02)阻尼            人工阻尼      

10、60;     數(shù)值阻尼ii"hx.E(q-4G1%Nh3? :wW   (03)每步求解方法    矩陣乘法            線性方程組zpaf c   (04)大矩陣（總剛）  否                  是+ynI2p J*j   (05

11、)數(shù)據(jù)存貯量      小                  大C1L zfIg;Kp,W  (06)每步計算速度    快                  慢_3&wJP   (07)迭代收斂性      無

12、                  有#sMw#EOe3W  (08)確定解          有確定解            可能是病態(tài)無確定解)F.Rd9g?-AUz)S'w zr5E   (09)時步穩(wěn)定性      有條件  

13、0;           無條件HycL iM  (10)時間步          小                  大.cQ2X_2Y)G&b  (11)計算精度        低    

14、60;            高7p x4yo:y-P8g;tG =dcO(t/b (01)是明顯不對的，只是對兩種方法的初級理解，(02)也是同樣。下面要詳細討論這兩點。/o*x'G5b|(03)是每一步求解的方法，(04)(05)(06)(07)(08)是由(03)所決定的，它們不是兩種方法的基本特點。同樣，(09)是時間步選擇的方法，(10)(11)是由(09)所決定的。s|S(G;i2S通過(03)(09)可以得到兩種方法的計算特點，顯式算法是每一步求解為矩陣乘法，時間步選擇為條件穩(wěn)定；隱式算法是每

15、一步求解為線性方程組求解，時間步選擇為無條件穩(wěn)定。;p$n:G%m sk2Z7K 下面主要分析兩種方法的應用范圍。Qv1MR2L7c,x0swY a)在求解動力學問題時，將方程在空間上采用有限元法（或其他方法）進行離散后，變?yōu)槌Ｎ⒎址匠探MM.u+C.u+Ku=f。求解這種方程的其中兩種方法為，中心差分法和Newmark法。采用中心差分法解決動力學問題被稱為顯式算法，采用Newmark法解決動力學問題被稱為隱式算法。!c7(w_i9y:J b)在求解動力學問題時，離散元法（也有其他方法）主要有兩種思想：動態(tài)松弛法（向后時步迭代），靜態(tài)松弛法（每一步要平衡）。動態(tài)松弛法是顯式算法，靜態(tài)松弛法是隱式

16、算法。其中沖壓成型就是動態(tài)松弛法的主要例子。kx:F*fC$j57A JH#Th*Lc)在求解靜力學問題時，有時候?qū)⑵淇醋鲃恿W問題來處理而采用動態(tài)松弛法，這是顯式算法。Flac就是主要例子。5iI(GO-d 最后總結(jié)，e+qfz:C U*w =B"W2sN'J5r.X0P                      顯式算法            隱式算法&#

17、39;Hj)T7I+oB-O S5oN5cF  (01)每步求解方法    矩陣乘法            線性方程組wT xjb   (02)時步穩(wěn)定性      有條件              無條件4F*W&D)v9kJ -B+ Ssh/S2U  (03)適用問題 

18、60;     動力中心差分法    動力Newmark法b:Wo$ X e4j                       動力動態(tài)松弛法    動力靜態(tài)松弛法y Ne(K-i#a-K                      靜力動態(tài)松弛法%LbAf X3Fc;r6 =Me2X0m-d&