abaqus剪切的案例
有限元理論基礎及Abaqus內部實現方式研究系列15： 殼的剪切應力
我們關注CAE中的結構有限元，所以主要選擇了商用結構有限元軟件中文檔相對較完備的Abaqus來研究內部實現方式，同時對某些問題也會涉及其它的Nastran/Ansys等商軟。為了理解方便有很多問題在數學上其實并不嚴謹，同時由于水平有限可能有許多的理論錯誤，歡迎交流討論，也期待有更多的合作機會。
iSolver介紹視頻：
http://www.yqgqt.org.cn/college/video/c12884
==第15篇：殼的剪切應力 ==
自編有限元應力的校核除了Mises等合力外，也應該校核各個應力分量。材料力學中六個應力分量如下：
其中Tau11，Tau22，Tau33為正應力，Tau12,13,23為三個剪切應力，對殼來說，Tau33=0，Tau12為面內剪應力，Tau13,23即為本文所說的橫向剪切應力。
最近在做iSolver殼的應力分量和Abaqus比對時，發現Abaqus的橫向剪切應力和預想的不一致。iSolver按照常用的殼的理論得到的剪切應力是個與厚度無關的常量，但Abaqus的橫向剪切應力分量TSHR13，TSHR23，在各個截面方向積分點section point不一樣。
花了點時間細致的研究了一下，猜測Abaqus中剪切應力TSHR13、23是真實應力，但有限元理論和iSolver中計算的是板殼近似理論中平均剪切應力。本章將介紹殼單元中實際的和板殼近似理論中的剪切應力，也猜測了Abaqus的內部實現流程，最后通過一個算例來驗算Abaqus中的真實的剪切應力，并通過iSolver來計算板殼理論的平均剪切應力。
1.1 殼的真實的剪切應力
剪應力是材料由于抗拒面之間的滑動而產生的沿表面方向的應力。殼的中間層存在剪切應力，這個可以通過下面簡單的例子驗證。
展開 有限元理論基礎及Abaqus內部實現方式研究系列3：S4殼單元剪切自鎖和沙漏控制
圖1：剪切自鎖
圖2：沙漏
===S4殼單元剪切自鎖和沙漏控制研究總結===
完全積分單元才有剪切自鎖，雖然Abaqus的S4單元是完全積分，但內部已經做了修正完全消除了剪切自鎖，所以不需要用戶做任何設置。
減縮積分單元才有沙漏現象，Abaqus的S4R默認增加一個人工的沙漏剛度來控制沙漏現象，如果發現結果還是不理想，那么需要采用其它建模方法才能控制沙漏了。
Abaqus針對剪切自鎖和沙漏控制做的修正如下：
單元類型
元素
修正情況
剪切自鎖
S4
薄膜剛度
假設應變方法修正
面外彎曲剛度
無
面外橫向剪切剛度
假設應變方法修正
S4R
所有項
無
沙漏控制
S4
薄膜剛度
無
面外彎曲剛度
無
面外橫向剪切剛度
沙漏控制，和Belytchko公式不一致，暫時沒研究
S4R
薄膜剛度
和Belytchko公式一致，因子取0.005G
面外彎曲剛度
和Belytchko公式一致，因子取0.00375G*4
面外橫向剪切剛度
沙漏控制，和Belytchko公式不一致，暫時沒研究
詳細研究方法，見附件：
有限元理論基礎及Abaqus內部實現方式研究系列3：S4殼單元剪切自鎖和沙漏控制(SnowWave02 20171018).pdf
展開 【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幫助文檔中，鋁合金三點彎曲的案例（延性損傷+剪切損傷）
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為厚殼單元的一種，在設置單元屬性時會有明確的說明：
現在，測試一下薄殼與厚殼計算結果的差異有多大。
展開 新手最近學UMAT,VUAMT的學習筆記（1），請大神多多指點。
微信公眾號關注這個橙子不太冷，一只ABAQUS菜鳥的學習之路
看了幾天幫助文檔，
發現學習UMAT、VUMAT
需要一些基礎知識，
沒有相關的知識一個簡單的
子程序都看不懂，
至少需要有力學基礎，數學基礎等，
下面介紹下目前遇到的一些基礎知識。
愛因斯坦求和約定
某項中某指標重復出現兩次，
該指標稱為啞標，例如，
拉丁指標i、j、k...表示三維，
希臘指標α、β、γ...表示二維。
在應力應變方程中經常遇到，
第一次看我沒太注意，到后面看子程序，
發現程序中表達的應力應變方程和
之前的對不上，才考慮到可能是個函數。
拉美常數
一階拉美常數λ
表示材料的壓縮性
二階拉美常數μ表示剪切模量，
大致與G相當
ABAQUS中的默認慣例
應力存儲順序，
Standard中
Explicit中
ABAQUS把剪切應變作為工程剪切應變γ
以上就是目前看幫助文檔中
遇到的一些基礎問題，
通過上網查閱所了解的。
先說到這，下次再說說
幫助文檔中的例子。
還請各位大神多多指點
展開 abaqus顯式動力學在煤層頂板垮落中的應用
abaqus中的剪切損傷模型+DP塑性準則可以很好的模擬煤層頂板的垮落破壞，但計算時容易產生不收斂問題。
abaqus中的顯式動力學不存在不存在收斂性問題。
所以，剪切損傷模型+DP塑性準則+顯式動力學=完美的煤層頂板仿真模擬。
模擬效果如下：
視頻連接：
仿真任務說明書：
計算任務描述：模擬在煤層巷道開挖后，煤層頂板的垮落情況。
仿真計算所采用的設備基本情況：華碩四核筆記本電腦，酷睿i5處理器，2G內存。
計算模型的處理技術：將實際的三維地質情況，假定為平面應變彈性力學問題；對煤層頂板的力學行為假定為彈性-DP塑性-剪切損傷。
計算機的耗時情況：5min。
仿真結果：在煤層巷道開挖后，煤層頂板可以完全垮落，放頂采煤效果100
展開 
有限元理論基礎及Abaqus內部實現方式研究系列7：C3D8R六面體單元的剛度矩陣
http://www.yqgqt.org.cn/content/post/343905
第三篇：S4殼單元的剪切自鎖和沙漏控制。介紹Abaqus的S4單元如何來消除剪切自鎖以及S4R如何來抑制沙漏的。
http://www.yqgqt.org.cn/content/post/350865
第四篇：非線性問題的求解。介紹Abaqus在非線性分析中采用的數值計算的求解方法。
http://www.yqgqt.org.cn/content/post/360565
第五篇：單元正確性驗證。介紹有限元單元正確性的驗證方法，通過多個實例比較自研結構求解器程序iSolver與Abaqus的分析結果，從而說明整個正確性驗證的過程和iSolver結果的正確性。
https://www.yqgqt.org.cn/content/post/373743
第六篇：General梁單元的剛度矩陣。介紹梁單元的基礎理論和Abaqus中General梁單元的剛度矩陣的修正方式，采用這些修正方式可以得到和Abaqus梁單元完全一致的剛度矩陣。
https://www.yqgqt.org.cn/content/post/403932
展開 【加密技術】Abaqus inp文件加密方法_教學
對materials_enc.inp文件解密
Abaqus通過隨機字符的移位運算對文件進行加密，對于同一個被加密的文件，每次加密得到的密文都是不同的。但是，只要在執行加密命令時，沒有使用include_only選項，密文都可以通過密碼解密，還原成原始文件。
Abaqus文件加密與解密示例
同樣，也可以跳過解密直接運行求解文件，這時需要使用*include將加密后的材料數據inp文件和密碼包含在內。
跳過解密直接求解
如文章開頭所述，假設我們想要對方只能將密文包含在求解文件中進行分析、查看計算結果，但無法解密、獲得材料或連接器參數，則可以在執行加密的命令行中添加include_only選項。
示例：
abaqus encrypt input=materials.inp include_only output=materials_enc.inp password=USim8888
對于采用這種方式加密的文件，非授權用戶即使獲取了密碼也無法將其還原。
參考資料：
Abaqus Documentation
【螺栓斷裂】Abaqus韌性損傷與剪切損傷準則-工程應用案例，獲取inp文件：鏈接。
展開 【技術鄰雙十一來啦】即日起至11月13日，千套CAE/CAD視頻六折起，錯過等一年
-方塊粘結滑移模擬（cohesive element和cohesive contact...
6
6
ABAQUS-鋼筋混泥土梁CDP（混泥土損傷塑性）模擬
6
12
ABAQUS-三維復合材料裂紋擴展失效模擬（XFEM）
6
18
ABAQUS-股骨沖擊骨折模擬
6
6
ABAQUS-滾石沖擊鐵絲網模擬
6
9
ABAQUS-超彈性橡膠襯套壓縮模擬
6
6
ABAQUS-高速剪切載荷下的動態裂紋擴展（XFEM）
6
9
ABAQUS-纖維纏繞模擬
6
6
ABAQUS-螺紋連接水壓滲透模擬
6
6
ABAQUS-拉壓循環載荷試驗模擬
6
6
ABAQUS-鋼球淬火熱傳模擬
6
6
ABAQUS-復合膠黏層材料重物沖擊破壞模擬
6
9
ABAQUS-復合膠黏層材料拉伸破壞模擬
6
12
ABAQUS-球殼內爆（TNT）模擬
6
6
ABAQUS-鐘擺撞擊鋁板模擬
6
6
ABAQUS-橫梁撞擊模擬wmv
6
6
ABAQUS-車輪軌道模擬
6
6
ABAQUS-壓彎裂紋擴展模擬（XFEM-無預置裂紋）
6
12
ABAQUS-銷軸與套筒過盈分析
6
6
ABAQUS-內環外環過盈分析（2D）
6
6
ABAQUS-水滴滴落到薄板模擬（SPH方法）
6
6
ABAQUS-薄板彎曲成型模擬
6
6
ABAQUS-可樂罐跌落模擬（CEL）
6
6
ABAQUS-沖擊損傷模擬（J-C模型）
6
9
ABAQUS-子 彈穿過裝水水管模擬（SPH）
6
12
ABAQUS-鋁筒非線性屈曲分析
展開 