ABAQUS/Explicit中的準(zhǔn)靜態(tài)分析

第十講ABAQUS/analysis_ABAQUS/Explicit中的準(zhǔn)靜態(tài)分析概述簡(jiǎn)介載荷速率能量平衡質(zhì)量縮放總結(jié)ABAQUS/analysis_ABAQUS/Explicit中的準(zhǔn)靜態(tài)分析簡(jiǎn)介ABAQUS/analysis_ABAQUS/Explicit中的準(zhǔn)靜態(tài)分析簡(jiǎn)介顯式方法是真正的動(dòng)力學(xué)過(guò)程。它最初用于模擬高速碰撞問(wèn)題。它用于求解結(jié)構(gòu)的動(dòng)力平衡狀態(tài)。在求解過(guò)程中，慣性力起到?jīng)Q定性的作用。非平衡力以應(yīng)力波的方式在相鄰單元之間傳播。穩(wěn)定時(shí)間增量一般較小。ABAQUS/analysis_ABAQUS/Explicit中的準(zhǔn)靜態(tài)分析簡(jiǎn)介顯式動(dòng)力學(xué)方法還可以模擬準(zhǔn)靜態(tài)問(wèn)題，比如金屬成型過(guò)程，但是需要特殊的考慮：如果以自然時(shí)間周期計(jì)算，用顯式動(dòng)力學(xué)方法求解準(zhǔn)靜態(tài)問(wèn)題是不切實(shí)際的。一般需要上百萬(wàn)的時(shí)間增量。對(duì)稱L型截面的軋制圓柱鋼坯的鐓鍛ABAQUS/analysis_ABAQUS/Explicit中的準(zhǔn)靜態(tài)分析簡(jiǎn)介為節(jié)省計(jì)算時(shí)間，可以在模擬過(guò)程中人為的增加軋制過(guò)程的速度。在增加軋制速度之后，靜平衡問(wèn)題演化為動(dòng)平衡問(wèn)題。慣性力的影響將會(huì)增加。準(zhǔn)靜態(tài)分析的目標(biāo)就是：在慣性力的影響較小的前提之下，盡量縮短計(jì)算的時(shí)間周期。ABAQUS/analysis_ABAQUS/Explicit中的準(zhǔn)靜態(tài)分析載荷速率ABAQUS/analysis_ABAQUS/Explicit中的準(zhǔn)靜態(tài)分析載荷速率在模擬過(guò)程中，人為的增加準(zhǔn)靜態(tài)成型過(guò)程的速度是必要的，它可以讓求解過(guò)程更經(jīng)濟(jì)。但是，在不使結(jié)果退化的前提下，究竟可以把速度提高多少呢？比如，金屬成型過(guò)程中，典型的工具速度大約為1m/s的數(shù)量級(jí)。這個(gè)速度與金屬中的典型波速相比是非常小的（鋼中的波速為5000m/s）。一般推薦的載荷速率為材料中波速的1%。ABAQUS/analysis_ABAQUS/Explicit中的準(zhǔn)靜態(tài)分析載荷速率推薦的方法：以不同的速率多次模擬（比如，工具的速度為100,50,5m/s）。因?yàn)橐缘偷妮d荷速率進(jìn)行分析的時(shí)間比較長(zhǎng)，所以以從高的載荷速度到低的載荷速度進(jìn)行分析。檢查結(jié)果（變形形狀、應(yīng)力、應(yīng)變、能量），分析不同載荷速率對(duì)結(jié)果的影響。在顯式板金成型模擬過(guò)程中，過(guò)大的工具速度將抑制起皺現(xiàn)象，并激起非真實(shí)的局部拉伸。在屈曲成型過(guò)程中，過(guò)大的工具速度將引起“噴注”效應(yīng)—水動(dòng)力學(xué)響應(yīng)（下頁(yè)中有圖形）。ABAQUS/analysis_ABAQUS/Explicit中的準(zhǔn)靜態(tài)分析載荷速率工具速度對(duì)變形形狀的影響噴注工具速度=500m/s

工具速度=10m/s噴注考慮下面的屈曲成型過(guò)程（軸對(duì)稱模型180

的截面）。當(dāng)工具速度非常大時(shí)，產(chǎn)生高度局部化的變形（噴注）。ABAQUS/analysis_ABAQUS/Explicit中的準(zhǔn)靜態(tài)分析載荷速率例子：板金右圖為汽車(chē)門(mén)標(biāo)準(zhǔn)門(mén)梁的撞擊測(cè)試的簡(jiǎn)化模型。圓梁在每個(gè)端點(diǎn)固定，接觸剛體圓柱后變形。測(cè)試為準(zhǔn)靜態(tài)的。剛體圓柱與可變形梁的碰撞ABAQUS/analysis_ABAQUS/Explicit中的準(zhǔn)靜態(tài)分析速度25m/s:好的全局結(jié)果0.1mV速度400m/s:局部效應(yīng)V載荷速率如果碰撞速度很高，400

m/sec，變形高度局部化，梁沒(méi)有結(jié)構(gòu)響應(yīng)。靜態(tài)測(cè)試中的主要響應(yīng)為梁的一階模態(tài)。該模態(tài)的頻率用于預(yù)計(jì)碰撞速度。一階頻率大約為250

Hz。碰撞在4微秒內(nèi)完成。利用25

m/sec的速度碰撞，圓柱在4微秒內(nèi)向梁推動(dòng)0.1

m。ABAQUS/analysis_ABAQUS/Explicit中的準(zhǔn)靜態(tài)分析載荷速率為什么25m/sec的速度是合適的？一階頻率(f)大約為

250Hz。相應(yīng)的周期為

t=0.004秒。在此周期內(nèi)，剛體圓柱被推向梁

d=0.1m。這樣，估計(jì)的速度

v

為

v=d/t=0.1/0.004=25m/sec。金屬波速為5000m/sec，所以碰撞速度25m/sec為波速的

0.5%

。碰撞速度應(yīng)該小于材料波速的1%。在分析步內(nèi)，以光滑緩坡的方式，把碰撞速度從零增加到所施加的碰撞速度可以得到更加精確的解。ABAQUS/analysis_ABAQUS/Explicit中的準(zhǔn)靜態(tài)分析載荷速率SMOOTHSTEP幅值曲線通過(guò)逐步施加載荷的方式，可以提高準(zhǔn)靜態(tài)解的精度：工具中的常速度條件導(dǎo)致金屬毛坯遭受突然的沖擊載荷。這將導(dǎo)致應(yīng)力波穿過(guò)毛坯，產(chǎn)生不希望發(fā)生的結(jié)果。以斜坡的方式，使工具速度從零逐步增加，將減小這些不利的效應(yīng)。出于同樣的原因，在分析結(jié)束移除工具的過(guò)程中，也以斜坡的方式將工具速度逐步降為零。ABAQUS/analysis_ABAQUS/Explicit中的準(zhǔn)靜態(tài)分析載荷速率SMOOTHSTEP幅值定義兩個(gè)幅值之間以5階多項(xiàng)式過(guò)渡。比如，在過(guò)渡開(kāi)始和結(jié)束時(shí)一階和二階時(shí)間導(dǎo)數(shù)為零。在使用SMOOTHSTEP定義位移時(shí)間歷程時(shí)，每個(gè)指定的幅值處的速度和加速度為零。*AMPLITUDE,NAME=SSTEP,DEFINITION=SMOOTHSTEP0.0,0.0,1.0E-5,1.0*BOUNDARY,TYPE=DISPLACEMENT,AMP=SSTEP12,2,2,2.5ABAQUS/analysis_ABAQUS/Explicit中的準(zhǔn)靜態(tài)分析能量平衡ABAQUS/analysis_ABAQUS/Explicit中的準(zhǔn)靜態(tài)分析能量平衡能量平衡方程可以用于幫助評(píng)估計(jì)算結(jié)果是否為合理的準(zhǔn)靜態(tài)響應(yīng)。在ABAQUS/Explicit中，能量平衡可以寫(xiě)為

其中EI

為內(nèi)能

（包括彈性應(yīng)變能、塑性應(yīng)變能和與沙漏控制相關(guān)的偽能）。

EV

為粘性機(jī)制耗散的能量。EFD

為摩擦耗散的能量。EKE

為動(dòng)能。EW

為外力功ETOT

為系統(tǒng)的總能量。ABAQUS/analysis_ABAQUS/Explicit中的準(zhǔn)靜態(tài)分析能量平衡考慮單軸拉伸試樣的拉伸測(cè)試如果實(shí)際測(cè)試是準(zhǔn)靜態(tài)的，拉伸試件的外力功等于試件的內(nèi)能。單軸拉伸測(cè)試ABAQUS/analysis_ABAQUS/Explicit中的準(zhǔn)靜態(tài)分析能量平衡準(zhǔn)靜態(tài)測(cè)試的能量歷程顯示在右圖中：慣性力是可以忽略的。測(cè)試試件的材料速度是很小的。動(dòng)能是可以忽略的。當(dāng)測(cè)試的速度增加以后：試件的響應(yīng)偏離靜態(tài)、趨于動(dòng)態(tài)。因此，材料速度和動(dòng)能更加明顯。準(zhǔn)靜態(tài)拉伸測(cè)試的能量歷程ABAQUS/analysis_ABAQUS/Explicit中的準(zhǔn)靜態(tài)分析能量平衡因此，能量檢查為ABAQUS/Explicit金屬成型過(guò)程的結(jié)果是否反應(yīng)了準(zhǔn)靜態(tài)解提供了另外的評(píng)估方法。在主要的成型過(guò)程中，變形材料的動(dòng)能不可以超過(guò)內(nèi)能的一小部分。這個(gè)小部分一般為1–5%。因?yàn)槊髟诋a(chǎn)生顯著變形之前將被移動(dòng)，因此在成型過(guò)程的早期，一般很難達(dá)到這個(gè)值。使用光滑幅值曲線將改進(jìn)早期響應(yīng)。

不關(guān)心工具的動(dòng)能。從模型的全部動(dòng)能中減掉工具的動(dòng)能，或者限制變形組件的能量輸出。ABAQUS/analysis_ABAQUS/Explicit中的準(zhǔn)靜態(tài)分析能量平衡例子：圓柱杯的深沖壓右圖為有限元模型的1/4模型。在所有的接觸面中定義摩擦：沖頭和毛坯：m

=0.25。砧和毛坯：m

=0.125。毛坯夾具和毛坯：m

=0。通過(guò)為毛坯夾具施加22.87KN的向下力模擬深沖壓過(guò)程，并將沖頭向下移動(dòng)36mm。圓柱杯深沖壓的初始構(gòu)型ABAQUS/analysis_ABAQUS/Explicit中的準(zhǔn)靜態(tài)分析能量平衡嘗試三種不同的沖壓速度：3m/s30m/s150m/s下表中總結(jié)了每個(gè)圓柱杯深沖壓的計(jì)算費(fèi)用：沖壓速度

(m/s)時(shí)間增量正則化的CPU時(shí)間3(1X)279291.030(10X)27040.097150(50X)5290.019ABAQUS/analysis_ABAQUS/Explicit中的準(zhǔn)靜態(tài)分析能量平衡最終構(gòu)型中毛坯厚度云圖過(guò)大的沖壓速度導(dǎo)致結(jié)果與實(shí)際的物理現(xiàn)象不符。盡管計(jì)算費(fèi)用相差10倍，以30m/s和3m/s的速度沖壓的結(jié)果非常接近。Vpunch=3m/sVpunch=30m/sVpunch=150m/sABAQUS/analysis_ABAQUS/Explicit中的準(zhǔn)靜態(tài)分析能量平衡比較內(nèi)能和動(dòng)能當(dāng)沖壓速度為150m/s時(shí)，毛坯中的動(dòng)能與內(nèi)能相比占很大的比例。當(dāng)沖壓速度為3m/s和30m/s時(shí)，在成型過(guò)程中，動(dòng)能與內(nèi)能相比只占很小的一部分。ABAQUS/analysis_ABAQUS/Explicit中的準(zhǔn)靜態(tài)分析質(zhì)量縮放ABAQUS/analysis_ABAQUS/Explicit中的準(zhǔn)靜態(tài)分析質(zhì)量縮放人為的增加成型速度可以提高解的經(jīng)濟(jì)性。同時(shí)，材料應(yīng)變率以同樣的速度增加。如果材料對(duì)于應(yīng)變率是不敏感的，這是不相關(guān)的。如果模型中考慮應(yīng)變率敏感性，將導(dǎo)致錯(cuò)誤的結(jié)果。如果考慮率相關(guān)性，一般需要用自然時(shí)間周期模擬成型過(guò)程?？梢酝ㄟ^(guò)質(zhì)量縮放實(shí)現(xiàn)這樣的功能。ABAQUS/analysis_ABAQUS/Explicit中的準(zhǔn)靜態(tài)分析質(zhì)量縮放顯式動(dòng)力學(xué)過(guò)程穩(wěn)定極限的估計(jì)公式為

其中

Le

是最小的特征單元長(zhǎng)度，

cd

是材料的膨脹波速。泊松比為零線彈性材料的膨脹波速為

其中

E

為楊氏模量，

為材料密度。ABAQUS/analysis_ABAQUS/Explicit中的準(zhǔn)靜態(tài)分析質(zhì)量縮放如果以f2（輸入的Factor值）的方式人為的增加材料密度：膨脹波速以f的方式減小。穩(wěn)定時(shí)間增量以f的方式增加。通過(guò)質(zhì)量縮放的方式人為的增加穩(wěn)定時(shí)間，使得可以以自然時(shí)間周期分析成型過(guò)程。人為的增加工具速度之后，質(zhì)量縮放對(duì)慣性效應(yīng)具有同樣的影響。過(guò)多的質(zhì)量縮放將導(dǎo)致非真實(shí)的解。如果質(zhì)量縮放用于完全的動(dòng)態(tài)條件下，總質(zhì)量的變化應(yīng)該盡量?。ㄐ∮?%）。ABAQUS/analysis_ABAQUS/Explicit中的準(zhǔn)靜態(tài)分析質(zhì)量縮放利用*FIXEDMASSSCALING選項(xiàng)，可以進(jìn)行質(zhì)量縮放。*FIXEDMASSSCALING在分析步開(kāi)始時(shí)施加質(zhì)量縮放。句法：*FIXEDMASSSCALING,ELSET=name,FACTOR=f

2

在相應(yīng)的單元集中，每個(gè)單元的密度以

f

2的方式增加，因此以f的方式增加穩(wěn)定時(shí)間增量。ABAQUS/analysis_ABAQUS/Explicit中的準(zhǔn)靜態(tài)分析質(zhì)量縮放質(zhì)量縮放例子：右圖為低碳鋼平面應(yīng)變?cè)嚰睦鞙y(cè)試。由于對(duì)稱的原因，只選取模型的1/4。單軸拉伸測(cè)試ABAQUS/analysis_ABAQUS/Explicit中的準(zhǔn)靜態(tài)分析質(zhì)量縮放圖形顯示了三種分析的不同結(jié)果

（PEEQ云圖）。右圖中左邊的結(jié)果和中間的結(jié)果幾乎相同。中間的結(jié)果與左邊的結(jié)果相比，只需要1/5的計(jì)算時(shí)間。與原始的靜態(tài)解相比，右邊的結(jié)果基本是沒(méi)有意義的。25質(zhì)量縮放因子 100001變化實(shí)際應(yīng)變速率ABAQUS/analysis_ABAQUS/Explicit中的準(zhǔn)靜態(tài)分析總結(jié)ABAQUS/analysis_ABAQUS/Explicit中的準(zhǔn)靜態(tài)分析總結(jié)過(guò)大的載荷速率將產(chǎn)生帶有顯著慣性效應(yīng)的結(jié)果。一般的建議為限制載荷施加的速率，比如，工具速度小于材料波速的1%。以斜坡的方式，從零增加載荷速率，也可以提高準(zhǔn)靜態(tài)響應(yīng)的準(zhǔn)確性。使用SMOOTHSTEP幅值定義質(zhì)量縮放可以用于率相關(guān)材料行為，允許以自然時(shí)間周期進(jìn)行模擬。能量平衡可以用于評(píng)估計(jì)算結(jié)果