abaqus 剪應(yīng)力
關(guān)注創(chuàng)建者：王靖雯 創(chuàng)建時間：2023-02-27
abaqus 剪應(yīng)力的實例教程
本文在GTN模型中引入剪應(yīng)力的影響，編寫了相關(guān)的VUMAT子程序。
GTN模型的屈服函數(shù)可以用下式表示
其中q1，q2是模型參數(shù)，取q1=1.5，q2=1,σ0為等效應(yīng)力，p為靜水應(yīng)力，q為Mises等效應(yīng)力；f為空洞的體積分?jǐn)?shù)。
p和q可以通過徑向返回算法得到
應(yīng)變控制的孔洞形核系數(shù)
GTN模型可以通過以下4個方程進行描述
Nahshon and Hutchinson考慮了剪應(yīng)力對模型的影響
于是孔隙體積分?jǐn)?shù)的演化可以通過下式描述
仿真計算得到的結(jié)果如下圖所示
最后，歡迎大家關(guān)注我們的微信公眾號聯(lián)系我們。
展開 切線剛度法程序設(shè)計
算法設(shè)計
1：定義程序需要用到的常數(shù)和變量
2：讀取ABAQUS定義的材料常數(shù)和狀態(tài)變量(這里只定義了一個狀態(tài)變量)，材料常數(shù)為，彈性模量E，泊松比Mu，屈服應(yīng)力Yield0，參數(shù)A,B,C,并且計算出剪切模量G,狀態(tài)變量為等效塑性應(yīng)變EQPLAS
3：讀取應(yīng)力分量，計算平均應(yīng)力，應(yīng)力偏量以及Mises等效應(yīng)力
平均應(yīng)力：
應(yīng)力偏量：
Mises等效應(yīng)力：
4：根據(jù)3計算的Mises等效應(yīng)力和2讀取的屈服應(yīng)力Yield0比較，如果Mises等效應(yīng)力小于屈服應(yīng)力，表明此時材料未屈服，那么轉(zhuǎn)到5，否則轉(zhuǎn)到6
5：雅可比矩陣,初始化為0，然后計算彈性矩陣，按照彈性理論更新應(yīng)力
6：雅可比矩陣,初始化為0
1計算切線模量H
注意到當(dāng)?shù)刃苄詰?yīng)時對應(yīng)于本構(gòu)關(guān)系的屈服點，此時的H不能通過上式計算，可以取此時的H為彈性模量
2計算w
根據(jù)前一章推導(dǎo)的公式
3計算等效塑性應(yīng)變增量DEQPLAS，并更新狀態(tài)變量
根據(jù)前一章推導(dǎo)的公式：
帶入后可得等效塑性應(yīng)變增量DEQPLAS
然后EQPLAS=DEQPLAS+EQPLAS更新狀態(tài)變量
4計算雅可比矩陣
首先初始化默認(rèn)為0，然后用下式計算雅可比矩陣
注意：abaqus的剪應(yīng)力方向跟彈性力學(xué)規(guī)定的方向不一致，所以上式的最后兩行的應(yīng)交換
5更新應(yīng)力
算法流程圖
圖5-5 切線剛度法算法流程圖
展開 
abaqus 剪應(yīng)力的相關(guān)專題、標(biāo)簽、搜索
abaqus 剪應(yīng)力的最新內(nèi)容
</p><p class="ql-align-justify">同時，在層合板內(nèi)每一相鄰實體層與 Cohesive 層界面之間，自動建立 Surface?to?Surface Penalty 面面接觸對（法向硬接觸、切向罰摩擦系數(shù) 0.3），實現(xiàn)層間正應(yīng)力與剪應(yīng)力的真實傳遞。該設(shè)置還原了文獻(xiàn)中有限厚度模型對最大中心位移和接觸時間更為準(zhǔn)確的預(yù)測能力。
TRIAX：應(yīng)力三軸度，可用于評估斷裂行為，尤其適用于延性損傷模型，用戶手冊定義如下：
編輯
跳轉(zhuǎn)
2.主應(yīng)力
主應(yīng)力是指該點應(yīng)力狀態(tài)中僅存在正應(yīng)力、無剪應(yīng)力的特殊方向上的應(yīng)力值。
https://www.yqgqt.org.cn/content/post/1266640
第二十七篇：Abaqus內(nèi)部計算和顯示的應(yīng)變。
而切向粘結(jié)則通過剪應(yīng)力分量S13/S23獲得，如圖所示。
組建剛度矩陣K，abaqus自己處理
2. 載荷列陣F，abaqus自己處理
3. x=K^-1*F，所有節(jié)點的位移 abaqus自己處理
4. 根據(jù)x，求每個單元的應(yīng)變增量值 abaqus自己處理
5. 根據(jù)單元應(yīng)變和單元剛度矩陣，求應(yīng)力。
【相關(guān)閱讀】
【JY】Abaqus殼單元概述與應(yīng)用（一）
【JY】Abaqus 三維應(yīng)力單元解析、選擇與應(yīng)用指南
【JY】Abaqus“殼”單元概述與應(yīng)用（二）——固體殼單元
傳統(tǒng)固體殼單元在處理幾何非線性、材料非線性及復(fù)雜邊界條件時，存在諸多難以克服的缺陷，這促使研究者探索新的單元構(gòu)造方法。
對于需要精確模擬三維應(yīng)力狀態(tài)的厚壁復(fù)合材料結(jié)構(gòu)，C3D8I 單元可能更合適，但需注意網(wǎng)格質(zhì)量。
在復(fù)合材料層合板三點彎曲試驗的基準(zhǔn)測試中，CSS8 單元能夠準(zhǔn)確模擬層間應(yīng)力分布，而 SC8R 單元則僅僅只是在計算效率上表現(xiàn)更佳。對于三層復(fù)合材料板，CSS8 單元可以準(zhǔn)確捕捉到中間層的剪應(yīng)力分布，而 SC8R 單元在相同情況下可能會出現(xiàn)中間層無剪應(yīng)力的現(xiàn)象。
，dmc分別為表征纖維拉伸，纖維壓縮，基體拉伸和基體壓縮這四種損傷模式的損傷變量，E 和μ 分別為彈性模量和泊松比，Smt和Smc 分別為復(fù)合材料基體拉剪耦合和壓剪耦合系數(shù)；定義等效應(yīng)力和等效位移：
其中l(wèi)c為特征單元長度。
<p>1、 引言</p><p>雙鋼板 - 混凝土組合結(jié)構(gòu)的抗剪性能與傳統(tǒng)鋼筋混凝土結(jié)構(gòu)存在顯著差異。該結(jié)構(gòu)通過拉結(jié)筋和栓釘實現(xiàn)鋼板與混凝土的連接，在剪力作用下易產(chǎn)生界面滑移，導(dǎo)致試件剛度與承載力下降。本案例聚焦于論文第 4 章雙鋼板 - 混凝土組合梁的建模復(fù)現(xiàn)，旨在通過 ABAQUS 有限元分析軟件，對組合梁抗剪性能進行數(shù)值模擬。需特別說明的是，本次復(fù)現(xiàn)僅涵蓋建模過程教學(xué)，不涉及曲線擬合內(nèi)容。
（滑移系剪切應(yīng)變率）
SDV91~102：resolvedstress_slip（滑移系臨界分切應(yīng)力）
SDV103~114：accumulatedshear_slip（滑移系累積剪應(yīng)變）
3.其他說明
（1）在終端輸入
abq_addUserOutput -m 114 JiaHe1.inp
可以直接修改inp文件的狀態(tài)變量個數(shù)。
