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