abaqus 罰剛度
關(guān)注創(chuàng)建者：王靖雯 創(chuàng)建時間：2023-02-27
abaqus 罰剛度的視頻教程
詳解abaqus接觸分析之二：abaqus中基于面的接觸詳
詳細(xì)講解了abaqus中基于面的接觸行為： 接觸模型 接觸離散化（N-S S-S的接觸原理，優(yōu)劣勢，選用建議） 接觸施加方法（罰剛度、拉格朗日法的使用方法和注意事項(xiàng)，已經(jīng)如何通過他們來控制穿透和收斂） 體間相互滑移（有限滑移和小滑移的原理和注意事項(xiàng)） 等。
￥25 1小時10分鐘 2143播放
查看
abaqus 罰剛度的實(shí)例教程
1.3.3.2 仿真結(jié)果
Abaqus中設(shè)置接觸，切向無摩擦，法向采用Penalty算法，且設(shè)置stiffness=10：
分析后得到的位移如下：
U1=1.047e-1，U2=8.968e-2，明顯和Penalty方法理論值有較大差異，要么Abaqus的Stiffness的含義和理論的罰剛度不同，要么就是接觸分析的迭代導(dǎo)致了和理論的差異，哪位大神要了解，歡迎討論。
如果Stiffness=100，分析后如下，顯然U1、U2更加接近0.1024這個實(shí)際問題的理論解。
1.4 視頻講解和操作驗(yàn)證演示
如果覺得上面的文字太復(fù)雜，也可以看一下視頻的簡要講解和軟件演示，地址如下：
https://www.yqgqt.org.cn/college/video/c12884 20理論系列文章35-接觸求解算法理論.mp4
1.5 總結(jié)
本文從接觸分析時如何將接觸的約束加入到有限元基本方程出發(fā)介紹了商軟中最常用也是最基本的兩個算法Lagrange和Penalty方法，可以看到Penalty方法中stiffness剛度的選取不同對最終的求解具有較大的影響，我們想要對標(biāo)Abaqus實(shí)現(xiàn)Abaqus完全一致的接觸分析結(jié)果，但苦于找不到Abaqus個罰剛度的選取規(guī)則，如果有哪位大神對這個問題比較了解，還請不吝賜教。
以往的系列文章：
1.6.1 ========第一階段========
第一篇：S4殼單元剛度矩陣研究。
http://www.yqgqt.org.cn/content/post/338859
第二篇：S4殼單元質(zhì)量矩陣研究。
http://www.yqgqt.org.cn/content/post/343905
第三篇：S4殼單元的剪切自鎖和沙漏控制。
展開 在ABAQUS/Standard的模擬中，在粘結(jié)（剪應(yīng)力小于μP）和滑移兩種狀態(tài)之間的不連續(xù)性可能導(dǎo)致收斂問題。因此，只有當(dāng)摩擦力對模型的響應(yīng)有顯著影響時，才應(yīng)該在模型中包含摩擦。如果在有摩擦的接觸模擬中出現(xiàn)了收斂問題，首先應(yīng)該嘗試的診斷和修改問題的方法之一就是在無摩擦的情況下重新運(yùn)算。一般情況下，ABAQUS/Explicit引入摩擦并不會引起附加的計算困難。
模擬理想的摩擦行為可能是非常困難的，因此在默認(rèn)的大多數(shù)情況下，ABAQUS使用一個允許“彈性滑動”的罰摩擦公式，“彈性滑動”是在粘結(jié)的接觸面之間所發(fā)生的小量的相對運(yùn)動。ABAQUS自動地選擇罰剛度，因此，這個允許的“彈性滑動”是單元特征長度的很小一部分。罰摩擦公式適用于大多數(shù)問題，包括在大部分金屬成形問題中的應(yīng)用。
必須包含 理想的粘結(jié)-滑動摩擦行為的問題中，在ABAQUS/Standard中使用 Lagrange摩擦公式，在ABAQUS/Explicit中使用 動力學(xué)摩擦公式。在計算機(jī)資源的消耗上，Lagrange摩擦公式更加昂貴，因?yàn)閷τ诿總€采用摩擦接觸的表面節(jié)點(diǎn)，ABAQUS/Standard應(yīng)用附加的變量，求解速度也是更慢的，一般需要附加的迭代。
ABAQUS/Explicit中摩擦約束的動力學(xué)施加方法是 基于預(yù)測/修正算法。在預(yù)測構(gòu)形中，應(yīng)用與節(jié)點(diǎn)相關(guān)的質(zhì)量、節(jié)點(diǎn)滑動的距離和時間增量來計算用于保持另一側(cè)表面上節(jié)點(diǎn)位置所需要的力。如果在節(jié)點(diǎn)上應(yīng)用這個力計算得到的切應(yīng)力大于 τ
crit，則表面是在滑動，并施加了一個相應(yīng)于 τ
crit的力。在任何情況下，對于在處于接觸中的從屬節(jié)點(diǎn)與主控表面的節(jié)點(diǎn)上，這個力將導(dǎo)致沿表面切向的加速度修正。
展開 
abaqus 罰剛度的相關(guān)專題、標(biāo)簽、搜索
abaqus 罰剛度的最新內(nèi)容
</p><p class="ql-align-justify"><strong>接觸與約束建模（若涉及）</strong></p><p>2D/3D 接觸、摩擦、粘著/分離判定、主從面、罰項(xiàng)與拉格朗日乘子等實(shí)現(xiàn)。</p><p>接觸探測、接觸對的激活/去激活規(guī)則，以及對接觸剛度的處理。
界面力的大小與發(fā)生的貫穿相對位移成正比，即
(4)
其中：—罰剛度系數(shù)，僅取決于材料特性。
圖2 罰函數(shù)耦合算法
1.4.2.3 研究現(xiàn)狀
CEL算法在解決流固耦合問題時具有強(qiáng)大的優(yōu)勢，目前已有廣泛的應(yīng)用。左輝等人對比了簡單氣體壓縮過程的仿真結(jié)果和解析結(jié)果，成功證明了Abaqus中理想氣體狀態(tài)方程的可靠性與正確性，并應(yīng)用于產(chǎn)品分析。
以2節(jié)點(diǎn)單元為例，由和的插值表達(dá)式，將得到：
當(dāng)梁的高度越來越小時，通過罰函數(shù)迫使約束條件得以實(shí)現(xiàn)，隱含著上式右端的常數(shù)項(xiàng)和一次項(xiàng)系數(shù)分別等于零，即：
第一個式子表示單元內(nèi)梁的平均轉(zhuǎn)動應(yīng)等于梁法向位移的變化率，也即表示梁中點(diǎn)的直線法假設(shè)，在計算中是可以實(shí)現(xiàn)的。
在顯示動力學(xué)接觸設(shè)置中，abaqus軟件默認(rèn)的是動態(tài)接觸算法（Kinematic contact algorithm）,當(dāng)開啟多線程時就會由于求解速度過高而產(chǎn)生計算的不穩(wěn)定性，而該算法的接觸約束嚴(yán)格性很高，因此當(dāng)遇到求解不穩(wěn)定時就會產(chǎn)生報錯從而導(dǎo)致計算終止。由于罰函數(shù)法的接觸約束嚴(yán)格性要低于動態(tài)接觸算法，因此改為罰函數(shù)法（penalty contact method）即可。
考慮到?jīng)_擊過程存在摩擦，將接 觸面法向定義為“硬”接觸，切向摩擦采用罰函數(shù) 算法，摩擦因數(shù)設(shè)置為 0.1[15-16]。
隱式分析法在接觸條件復(fù)雜的模型時，計算結(jié)果不易收斂，而顯示分析不但能夠保證計算結(jié)果的收斂性，還能很好地模擬準(zhǔn)靜態(tài)工況，所以采用Abaqus/Explicit顯示分析進(jìn)行仿真模擬。
(2)切削刀具的強(qiáng)度和剛度遠(yuǎn)高于巖石材料，將刀具假設(shè)為剛體，且在切削過程中不發(fā)生磨損。
(3)不考慮溫度對切削過程的影響。
案例—天然橡膠隔震支座豎(軸)向力學(xué)性能
【JY】Abaqus6.14-4如何關(guān)聯(lián)fortran？
結(jié)論
本文借助ABAQUS有限元分析軟件實(shí)現(xiàn)了對某型USM系列軌道橡膠減振墊的數(shù)值分析，得到了形變、應(yīng)力、應(yīng)變和剛度等一系列結(jié)果，并對其中結(jié)構(gòu)不合理的地方提出了改進(jìn)建議。
實(shí)際上，輥?zhàn)訉﹄姵匕臄D壓速度很緩慢，擠壓的過程中可以看成無數(shù)個微小時間間隔的靜態(tài)過程，因此電池包的擠壓仿真用Abaqus做準(zhǔn)靜態(tài)仿真更準(zhǔn)確。由于水平有限和硬件不足，僅僅以低速碰撞替代準(zhǔn)靜態(tài)擠壓，為電池包的擠壓仿真探索思路。
在本案例中用到的卡片和關(guān)鍵字總結(jié)如下。
