abaqus 簡單剪切的案例
平面拉伸（純剪切），簡單剪切
foam, rubber, hyperelastic, hyperfoam, low density foam
https://www.sciencedirect.com/topics/engineering/pure-shear
https://en.wikipedia.org/wiki/Simple_shear#Simple_shear_stress%E2%80%93strain_relation
Simple shear is a deformation in which parallel planes in a material remain parallel and maintain a constant distance, while translating relative to each other.
https://www.youtube.com/watch?app=desktop&v=srgMoMoHU1A
展開 基于黃永剛umat的FCC多晶在簡單剪切變形下的有限元分析------案例二 ￥99
基于黃永剛umat的FCC多晶在簡單剪切變形下的有限元分析
1，使用neper軟件生成包含500個晶粒的幾何模型
2，賦予相應的材料參數（基于腳本完成材料的批量賦值）
3，施加相應的邊界條件（X0方向完全固定，X1方向施加Z1方向0.15mm的位移）
4，結果與后處理
5，模型文件見附錄
6，提供材料屬性賦予腳本
【ABAQUS筆記】什么是剪切閉鎖？剪切閉鎖會導致什么？
完全積分的二階單元每個方向上有3個積分點，如下圖：
探究元素的階數(一階與二階)和網格密度對結果精度的影響
采用了幾種不同單元和網格密度,在Abaqus/Standard模擬懸臂梁問題：
各個模擬的末端撓度位移與理論值3.09 mm的比值如下：
對于線性的，完全集成的單元，單元再密都不準。以上撓度模擬不準是因為剪切鎖緊引起的，這是所有完全積分的一階實體單元都存在的問題==
什么是剪切閉鎖？（shear lock）
剪切鎖緊會導致單元在彎曲時過于僵硬。
根據材料力學，一個材料微元承受彎矩M時的變形如下：
單元變形的應力情況：
這種不正確的剪切應力的產生是因為線性單元的邊無法彎曲。剪切應力的存在意味著單元應變導致剪切變形，而不是預期的彎曲變形，因此整體撓度較小，也就是說單元剛度太大了
剪切鎖定只影響完全積分的一階單元在彎曲載荷下的單元性能。
對于二階單元，剪切閉鎖不是問題，因為二階單元的邊是可以彎曲的，但它依舊有細微的剪切閉鎖效應。
所以，ABAQUS文檔建議：
只有相當確定載荷會在模型中產生很小的彎曲時，才推薦使用完全積分的一階單元。完全積分的二階單元在復雜的應力狀態也可能shear lock。
展開 ABAQUS中的單元選擇-理解剪切自鎖和沙漏
為了限制沙漏現象的擴展，ABAQUS引進了“防沙漏剛度”Hourglass stiffness，一般情況下采用默認值即可，如果確有需要可在圖1中的Hourglass control選項中設置。
圖4
四、小結
如果模型中有比較明顯的彎曲現象，為避免出現剪切自鎖現象，優先選擇二階單元，或者采用縮減積分方案（網格需要更細，通常厚度方向4層以上）。
來源： ABAQUS在巖土工程中的應用
abaqus剪切修正GTN模型的VUMAT子程序開發
這里對相應的算法進行簡要說明：
NH-GTN模型
屈服函數：
其中等效孔洞體積分數定義為：
孔洞體積分數包含新孔隙形核，原有空隙生長以及剪切相關的等效體積分數增加：
形核，生長，剪切相關體積分數的演化遵循：
其中：
剪切效應的修正，考慮應力狀態的影響
參數的物理含義如下
通過將文獻中的數值算法編程實現在VUMAT子程序中，可以用來實現對延性金屬材料在不同應力狀態下的損傷演化進行合理的數值預測，應用于金屬成型領域（沖壓，軋制，擠壓等）
預測修正后的模型應該在簡單拉伸情況下于abaqus自帶的GTN模型保持相同的損傷和其他狀態變量的分布，并在剪切情況中損傷發展顯著高于abaqus自帶的模型（自帶的模型忽略了剪切效應）。（為了進行對比使用于自帶的本構相同的硬化方式，模擬中使用了相同的質量縮放，但質量縮放容易產生數值振蕩，模擬的拉伸曲線存在波動。）
初步模擬結果：
拉伸情況（abaqus-VUMAT）
應力分情況
孔洞體積分數
剪切模型（abaqus-VUMAT）
不同變形時刻的應力分布
T=0.1s
局部放大圖
T=0.5s
局部放大圖
T=0.6s
局部放大圖
可以看到模型在拉伸預測中與原始模型保持一致，而在剪切修正后損傷發展顯著快于原始模型，利用作者提出的方法可以應用于復雜應力狀態下金屬材料的損傷分析，相關參數部分參考文獻，其中Kw=3.T1=0.2，T2=0.7.模擬結果符合文獻所提出方法的基本趨勢。
最后，如果有相關需要歡迎通過公眾號“320科技工作室”聯系我們。
展開 abaqus離散元做直接剪切試驗
需要源文件的可以抖音私信我哦
ABAQUS UMAT-混凝土受拉狀態下塑性損傷模型的簡單實現 ￥600
本文利用ABAQUS UMAT子程序，簡單實現了混凝土受拉狀態下的破壞。本構模型的實現算法摘抄自DeBorst的書籍《Nonlinear Finite Element Analysis of Solids and Structures》，基本如下：
為了簡化模型，筆者將書中損傷部分做了簡化，不再采用損傷屈服面進行判定。損傷影子w的計算直接由塑性等效應變確定。
在ABAQUS中建立100*100*100的立方體塊，試件的底部固定，頂部反復加載-卸載，通過UMAT得到的模擬結果如下：
Abaqus幫助文檔中，鋁合金三點彎曲的案例（延性損傷+剪切損傷）
threepointbending_alextrusion.rar
文檔.pdf
【螺栓斷裂】Abaqus韌性損傷與剪切損傷準則---{ 問題答疑 +工程案例 + 模型文件 } ￥99.9
Abaqus中韌性金屬失效分析需要定義c點的損傷初始化準則，以及cd段的損傷演化（損傷后材料剛度退化路徑）。材料軟化后可持續承載，直到達到d點，材料失效，失去承載能力。
圖1-韌性金屬的全載荷區間應力-應變曲線
圖2-韌性金屬的損傷準則
ABAQUS為韌性金屬提供不同的損傷初始化準則，大致分為兩種類型：
金屬裂紋的損傷初始化準則，包括韌性準則（ductile damage、Johnson-Cook damage）和剪切準則（shear damage）。也就是圖2中紅框內的三個準則，它們都屬于金屬承載后產生裂紋的準則。
金屬板的徑縮不穩定損傷初始化準則，包括幾種成形極限圖，用于評估鈑金件的可成形性。也就是紅框外的幾個準則，不在本文討論范圍。
圖3-漸進損傷失效分類【摘自Abaqus材料本構模型導圖，完整版鏈接】
····································常見問題解答····································
······Q1: 韌性準則和剪切準則有何不同？
······A1: 韌性金屬開裂有兩種主要機理，基于唯象觀察，仿真模擬這兩種機理時用到不同的損傷起始準則（hooputra2004）：
機理1，由于內部（微裂紋）的成核、生長和孔隙的聚集產生的韌性斷裂，這種情況下ductile damage、Johnson-Cook damage兩種韌性準則是適用的，常見于拉伸工況。
圖4-機理1韌性斷裂
機理2，由于剪力帶局部化產生的剪切斷裂，這時shear damage比較適合，常見于剪切工況。
展開 Abaqus中考慮橫向剪切的復合材料厚殼單元分析
在前面的兩篇文章中，已經對Abaqus復合材料殼單元分析模型的傳統建模方法和快捷建模方法做了詳細的介紹。熟悉Abaqus復合材料分析的人都知道，在采用二維Lamina復合材料模型配合殼單元進行分析時，材料參數中除了輸入兩個方向模量E1,E2，面內泊松比及面內剪切模量G12之外，還要額外輸入兩個面外剪切模量G13和G23，如下圖所示。
這里的面外剪切模量G13和G23就是用來考慮橫向剪切變形的。
一般，針對薄板結構（跨厚比大于20），通常做以下假設（Kirchhoof假設）：
(1)平行于中面 的各層互不擠壓：即垂直于中面法向的正應力很小，可以忽略。
(2)直法線假設：變形前垂直于中面的直線段，在變形后仍保持是直線，且仍垂直于變形后的中面。
(3)撓度沿板厚度方向的變化可以忽略，即統一厚度各點的撓度都 等于中面的撓度
(4)板的中面無伸縮和剪切變形
根據上述假設，板的橫向變形為零，相當于垂直于中面的各個面內剪切模量無窮大。薄板理論的假設在求解薄板問題時，精度足以滿足工程計算要求。
但對于中厚板或者厚板、集中力作用點附近、薄板邊界周圍以及開孔周圍，上述理論將不再適用，誤差大甚至會導致錯誤的結果，因此為了解決此類問題，便有了考慮剪切變形的中厚板理論。
那么在Abaqus分析中怎么考慮橫向剪切變形的影響呢？Abaqus默認的復合材料模型定義及截面屬性定義中是已經考慮了橫向剪切的，軟件會自動計算橫向剪切剛度。
而薄板、中厚板的區分在于單元類型，如下圖所示，S8R5為薄殼單元的一種，S8R為厚殼單元的一種，在設置單元屬性時會有明確的說明：
現在，測試一下薄殼與厚殼計算結果的差異有多大。
展開 有限元理論基礎及Abaqus內部實現方式研究系列3：S4殼單元剪切自鎖和沙漏控制
商用軟件對外就是一個黑盒子，除了開發人員，使用人員只能在黑盒子外猜測內部實現方式，一方面我們查閱Abaqus軟件手冊得到修正方法的說明，另一方面我們自己編程實現簡單的結構有限元求解器，通過自研求解器和Abaqus的結果比較結合理論手冊如同管中窺豹一般來研究Abaqus的修正方法，從而猜測商用有限元軟件的內部計算方法。為了理解方便有很多問題在數學上其實并不嚴謹，同時由于水平有限也可能有許多的理論錯誤，歡迎交流討論。
==以往的系列文章==
第一篇：S4殼單元剛度矩陣研究 http://www.yqgqt.org.cn/content/post/338859
研究基于Mindlin厚殼理論的S4殼單元的剛度矩陣在Abaqus中的實現方式
第二篇：S4殼單元質量矩陣研究 http://www.yqgqt.org.cn/content/post/343905
研究一致質量矩陣和集中質量矩陣在Abaqus的S4殼單元和Nastran的Quad4殼單元中的實現方式
==第三篇：S4殼單元剪切自鎖和沙漏控制==
商用有限元軟件的健壯性體現在對各種特殊情況，求解過程和解的正確性依然能得到保證，而這些特殊情況在自編程序中如果沒有考慮到，那么結果就可能相差極大。其中剪切自鎖和沙漏現象是最常見的會影響正確性的兩個特殊情況。這兩者具有相似性，所以我們在本文中一起研究Abaqus中線性殼單元S4針對這兩種情況下的內部實現方式。剪切自鎖和沙漏現象影響的是剛度矩陣和應力，我們研究方式是在自編程序iSolver中根據成熟的消除剪切自鎖和沙漏控制的理論實現剛度矩陣的修正，通過比較同一模型的Abaqus的剛度矩陣結果，結合幫助文檔猜測Abaqus軟件單元消除剪切自鎖和控制沙漏的內部實現方法。
展開 
有限元理論基礎及Abaqus內部實現方式研究系列15： 殼的剪切應力
1.4.4 平均剪切應力結果
Abaqus中無法輸出S13和S23，也就是平均剪切應力，為了便于大家理解S13這個值大概是什么量級，我們采用iSolver，在內部計算得到積分點上的剪切應力得到S13=8.1e5，而只有采用平均面積得到的V/(b*h)=2e6的2/5左右，就算加上剪切修正因子5/6也不是簡單的按平均面積得到的結果，猜測是因為在這種情況下，殼只在兩個節點處加載荷，而不是在一個自由邊上均勻加載，而平均面積計算的剪切應力是假定均勻加載得到，所以不一致。當然，有限元中采用這個積分點應力來求節點力，得到的節點力依然和外力是平衡的。
==總結==
本章介紹了殼單元中實際的和板殼近似理論中的剪切應力，也簡單猜測了一下Abaqus的內部實現流程，最后通過一個算例來驗算Abaqus中的真實的剪切應力。同時還有兩個疑問也希望與大家討論，得到大家的指點。
（1）我們暫時沒有找到方法來驗證用于Abaqus流程的剪切應力是取的平均應力還是真實應力，因為Abaqus的S13，S23沒有找到方法輸出，不知道誰了解怎么輸出殼的這兩個應力？
（2）按照下面的公式，Mises力應該包括正應力和剪切應力，但查看殼的Mises應力結果可以看出，殼的Mises力沒有計入Tau13，Tau23，不知道為什么這樣取Mises力？這樣用來校核殼的Mises應力達到屈服是否會有問題？
如果有任何其它疑問或者項目合作意向，也歡迎聯系我們：
snowwave02 From www.yqgqt.org.cn
email: snowwave02@qq.com
以往的系列文章：
第一篇：S4殼單元剛度矩陣研究。介紹Abaqus的S4剛度矩陣在普通厚殼理論上的修正。
展開 ABAQUS UEL - Embedded crack model 的簡單實現和應用 ￥1500
關于Embedded crack model的理論基礎和UEL實現可以參考，
https://zhuanlan.zhihu.com/p/496005714
https://zhuanlan.zhihu.com/p/486789331
本文介紹一種簡單實現 Embedded crack model 的方法以及該方法在有限元模擬中的表現（以三點受彎模型為例）。如上面兩篇連接里介紹的，在單元階段，開裂面上的應力可以寫為，
在材料階段，開裂面上的應力可以寫為，
兩個階段的應力需要相等，
這就形成了一個非線性式子，其中未知的開裂相對位移 e 就是我們想要求解的。對于非線性函數，我們可以采用牛頓拉弗森迭代，建立一個殘差方程，
殘差方程對相對位移的求導為，
算法實現如下，
考慮到在主應力平面中剪切力為 0，所以將 fs 的值設為 0，只對受拉破壞行為進行定義。
將以上概念通過UEL在Abaqus中實現。建立一個帶有預設缺口的粱模型。該模型的邊界條件為簡支，在粱中心處設置向下的位移。
展開 abaqus簡單分析完整操作步驟 ￥1
有任何疑問，歡迎聯系交流,QQ：1317425016.
Abaqus+Fluent簡單CFD實例詳解
1、實例簡介
2、仿真結果
3、詳細步驟
Abaqus+Fluent簡單CFD實例詳解1.rar
Abaqus+Fluent簡單CFD實例詳解2.rar
