搜索Abaqus有限元分析——動力學顯式有限元方法
2013-08-14 by:Abaqus軟件應用培訓中心 來源:仿真在線
Abaqus有限元分析——動力學顯式有限元方法
這一節(jié)包括ABAQUS/Explicit求解器的算法描述,在隱式和顯式時間積分之間進行比較,并討論了顯式方法的優(yōu)越性。
9.2.1 顯式時間積分
ABAQUS/Explicit應用中心差分方法對運動方程進行顯示的時間積分,應用一個增量步的動力學條件計算下一個增量步的動力學條件。在增量步開始時,程序求解動力學平衡方程,表示為用節(jié)點質量矩陣M乘以節(jié)點加速度等于節(jié)點的合力(在所施加的外力P與單元內力I之間的差值):
在當前增量步開始時(t時刻),計算加速度為:
由于顯式算法總是采用一個對角的、或者集中的質量矩陣,所以求解加速度并不復雜;不必同時求解聯(lián)立方程。任何節(jié)點的加速度是完全取決于節(jié)點質量和作用在節(jié)點上的合力,使得節(jié)點計算的成本非常低。
對加速度在時間上進行積分采用中心差分方法,在計算速度的變化時假定加速度為常數。應用這個速度的變化值加上前一個增量步中點的速度來確定當前增量步中點的速度:
速度對時間的積分并加上在增量步開始時的位移以確定增量步結束時的位移:
這樣,在增量步開始時提供了滿足動力學平衡條件的加速度。得到了加速度,在時間上“顯式地”前推速度和位移。所謂“顯式”是指在增量步結束時的狀態(tài)僅依賴于該增量步開始時的位移、速度和加速度。這種方法精確地積分常值的加速度。為了使該方法產生精確的結果,時間增量必須相當小,這樣在增量步中加速度幾乎為常數。由于時間增量步必須很小,一個典型的分析需要成千上萬個增量步。幸運的是,因為不必同時求解聯(lián)立方程組,所以每一個增量步的計算成本很低。大部分的計算成本消耗在單元的計算上,以此確定作用在節(jié)點上的單元內力。單元的計算包括確定單元應變和應用材料本構關系(單元剛度)確定單元應力,從而進一步地計算內力。
這里給出了顯式動力學方法的總結:
1. 節(jié)點計算
a. 動力學平衡方程
b. 對時間顯式積分
2. 單元計算
a. 根據應變速率,計算單元應變增量
b. 根據本構關系計算應力
c. 集成節(jié)點內力
3. 設置時間 t為,返回到步驟1。
9.2.2 比較隱式和顯式時間積分程序
對于隱式和顯式時間積分程序,都是以所施加的外力P、單元內力I和節(jié)點加速度的形式定義平衡:
其中M是質量矩陣。兩個程序求解節(jié)點加速度,并應用同樣的單元計算以獲得單元內力。兩個程序之間最大的不同在于求解節(jié)點加速度的方式上。在隱式程序中,通過直接求解的方法求解一組線性方程組,與應用顯式方法節(jié)點計算的相對較低成本比較,求解這組方程組的計算成本要高得多。
在完全Newton迭代求解方法的基礎上,ABAQUS/Standard使用自動增量步。在時刻增量步結束時,Newton方法尋求滿足動力學平衡方程,并計算出同一時刻的位移。由于隱式算法是無條件穩(wěn)定的,所以時間增量比應用于顯式方法的時間增量相對地大一些。對于非線性問題,每一個典型的增量步需要經過幾次迭代才能獲得滿足給定容許誤差的解答。每次Newton迭代都會得到對于位移增量的修正值。每次迭代需要求解的一組瞬時方程為
對于較大的模型,這是一個昂貴的計算過程。有效剛度矩陣是關于本次迭代的切向剛度矩陣和質量矩陣的線性組合。直到一些量滿足了給定的容許誤差才結束迭代,如力殘差、位移修正值等。對于一個光滑的非線性響應,Newton方法以二次速率收斂,描述如下:
迭代 相對誤差
1 1
2 10-2
3 10-4
. .
. .
. .
然而,如果模型包含高度的非連續(xù)過程,如接觸和滑動摩擦,則有可能失去二次收斂,而是可能需要大量的迭代過程。為了滿足平衡條件,減小時間增量的值可能是必要的。在極端情況下,在隱式分析中的求解時間增量值可能與在顯式分析中的典型穩(wěn)定時間增量值在同一量級上,但是仍然承擔著隱式迭代的高昂求解成本。在某些情況下,應用隱式方法甚至可能不會收斂。
在隱式分析中,每一次迭代都需要求解大型的線性方程組,這一過程需要占用相當數量的計算資源、磁盤空間和內存。對于大型問題,對這些方程求解器的需求優(yōu)于對單元和材料的計算的需求,對于在ABAQUS/Explicit中的分析這是類似的。隨著問題尺度的增加,對方程求解器的需求迅速地增加,因此在實踐中,隱式分析的最大尺度常常取決于給定計算機中的磁盤空間的大小和可用內存的數量,而不是取決于需要的計算時間。
9.2.3 顯式時間積分方法的優(yōu)越性
顯式方法特別地適用于求解高速動力學事件,它需要許多小的時間增量來獲得高精度的解答。如果事件持續(xù)的時間是非常短,則可能得到高效率的解答。
在顯式方法中可以很容易地模擬接觸條件和其它一些極度不連續(xù)的情況,并且能夠一個節(jié)點一個節(jié)點地求解而不必迭代。為了平衡在接觸時的外力和內力,可以調整節(jié)點加速度。
顯式方法最顯著的特點是沒有在隱式方法中所需要的整體切線剛度矩陣。由于是顯式地前推模型的狀態(tài),所以不需要迭代和收斂準則。
相關標簽搜索：Abaqus有限元分析——動力學顯式有限元方法 蘇州Abaqus分析培訓 寧波Abaqus有限元 Ansys有限元培訓 Ansys workbench培訓 ansys視頻教程 ansys workbench教程 ansys APDL經典教程 ansys資料下載 ansys技術咨詢 ansys基礎知識 ansys代做 Fluent、CFX流體分析
