abaqus材料斷裂的案例
ABAQUS材料斷裂與失效-XFEM|VCCT|COHESIVE|疲勞|侵蝕
ABAQUS材料斷裂與失效（XFEM|VCCT|COHESIVE|疲勞|侵蝕）篇
上個系列我們講了有關ABAQUS混凝土塑性損傷模型（CDP）的相關內容，得到了一些不錯的反饋；因此想趁著這股熱勁，繼續錄制一套課程；最終選擇了材料的斷裂與失效仿真。在很久以前，通用有限元軟件中，只有ABAQUS提供了擴展有限元方法（XFEM），一種不依賴于網格邊界的裂縫仿真方法。當時我還在學校，由于課題的原因，接觸了很多巖石斷裂和材料失效的實驗和仿真任務，從而也就使用上ABAQUS。
09年，想學好這款軟件并非易事，一方面是身邊沒有直接可以求助的前輩，因為大家都在用ANSYS；另一方面是當時市面上的相關書籍非常少（當時出的書都看過），加之自己英文不好，而且不怎么看幫助文檔，所以學習斷裂仿真就更加不易。后來有幸得到一份PPT課件：Modeling with Fracture and Fracture with Abaqus，是達索在2009年出的一套斷裂仿真的培訓教程，內容非常系統，不僅介紹了相關的斷裂理論、仿真方法和仿真技巧，還有實操案例。
達索09年斷裂課程內容目錄如下（公眾號中回復：“斷裂與失效”，可獲得此課件）：
雖然有好的資料在手，但看這里面五花八門的知識點，確實看著頭大，當時也是看了個一知半解，完成任務后，就放在一邊沒有去細細專研。
隨著時光推移，使用ABAQUS已經有八個年頭，仿真水平大幅提高，不知不覺也做了ABAQUS二次開發的工作，漸漸喜歡使用它來解決問題，自己沒課題，就幫別人做，也積累起一些經驗；但人的精力是有限的，所以以后想盡量做專一點，又不想將以前所學的知識輕易拋棄，因此打算將所專所長一一錄制成視頻教程，供大家學習，讓初學者少走一些彎路，視頻設置收費，是因為自己家庭條件并不好，收取少量費用還請大家理解。
展開 使用Abaqus求解金屬材料斷裂破壞實例
本文簡單介紹使用Abaqus計算帶有漸進損傷破壞參數的韌性金屬模型，圖 1為典型材料漸進損傷曲線，其中A點為漸進損傷起始點，AB段為材料損傷過程，點B為材料完全失效點。
圖 2為Abaqus漸進損傷破壞相關參數，Fracture strain為破壞應變、stress triaxiality為應力三軸度、strain rate為破壞應變率、displacement at failure為漸進損傷失效位移。
算例:
該模型分為兩部分，上端為限位座，限位座兩螺栓孔為固定約束，下端為限位塊，限位塊整個為剛性體，剛性參考點處施加強制位移，兩部分接觸位置定義接觸關系。
下表為整個模型的計算結果
使用abaqus求解金屬材料斷裂破壞實例.pdf
展開 ABAQUS多晶體材料斷裂模型
多晶體材料的斷裂研究有助于深入了解材料在微觀尺度下的力學行為，包括裂紋如何形成、擴展以及停止，這對于發展和完善固體力學和斷裂力學理論至關重要。本案例介紹在ABAQUS內基于Voronoi建立多晶體材料晶粒及晶界模型，并進行多晶材料的斷裂模擬。
多晶材料晶粒及晶界模型采用CAD Voronoi V3 多圖層版生成，插件可將不同組分的晶粒在CAD內進行分圖層繪制，可控制晶粒占比參數，以精確建立多晶體模型。
在AutoCAD內將不同成分的晶粒分別另存為dxf格式文件，并導入到ABAQUS建立草圖，利用草圖建立多組晶粒及晶界部件，本案例中，共建立了五種不同的晶粒。
新建荷載施加裝置，并與多晶體模型裝配為整體，同時對不同組分的晶粒及晶界設置材料。由于本案例研究多組分晶粒模型的斷裂情況，因此不同組分的晶粒設置了不同的損傷破壞材料參數。
設置加載塊及支座與試件間的接觸。
編輯
跳轉
將下部支座固定，上部施加豎向位移，完成載荷的設置。
進行網格劃分。
建立作業提交計算并查看多晶模型的開裂結果。
展開 材料的理論斷裂強度 附晶體材料強度與斷裂微觀理論下載
材料力學低碳鋼拉伸試驗中，材料的變形分為四個階段：彈性階段、屈服流動階段、強化階段和徑縮斷裂階段，如圖1，其中當材料經過d點后，材料很快發生斷裂，該點對應的應力σb即為強度極限。但這只是實驗觀察到的現象，它與材料的理論斷裂值還有很大的區別。
假設材料的斷裂是由于原子間距被拉的太遠，超過了極限從而發生的斷裂。我們知道，原子之間的力與原子間的距離存在一定的關系，當原子靠的特別近的時候，原子間存在排斥力，當原子離的比較遠的時候，原子間存在相互吸引力，在某一距離下，原子間的作用力為0，即平衡位置。
現在我們來考慮原子間的力與應力的關系，根據應力的定義
顯然，曲線上的最大值σm即代表原子間的最大結合力——理論斷裂強度，即在理論上認為材料應力超過σm時將被拉斷。作為一級近似，該曲線可用正弦曲線表示。
而實際上，對于純鐵的抗拉強度是只有170～270MPa左右，我們熟知的Q235鋼，其抗拉極限為375～460MPa，Q345鋼的抗拉強度約是490-620MPa，遠遠低于材料的理論斷裂強度。主要原因在于公式(11)表示的是理想材料的斷裂強度，也就是說材料中沒有任何的缺陷。但這是不可能的，材料在冶金、鑄造、加工等過程中難免會產生一些初始缺陷，造成應力集中從而大大降低了材料的強度缺陷。
下載地址：晶體材料強度與斷裂微觀理論
展開 
Abaqus復合材料3D Hashin失效準則，脆性斷裂-Vumat
對于纖維增強復合材料的模擬，在ABAQUS中，集成了二維Hashin失效準則與多種損傷演化準則，但缺少三維的復合材料本構模型。
參考已有的3Dhashin失效準則編寫復合材料脆性斷裂子程序。
首先介紹該子程序的使用方法
1.在ABAQUS中建立三維復合材料模型，這里建立一個簡單的方塊。1,2方向分別表示絲束的方向，3方向表示垂直于1，2的方向，也就是面外方向。
2.建立材料屬性（圖片中材料參數為假設值）
表1 16個參數對應含義
1
2
3
4
5
6
7
8
E11
E22
E33
G12
G13
G23
U12
U13
9
10
11
12
13
14
15
16
U23
1方向拉伸強度
1方向壓縮強度
2方向拉伸強度
2方向壓縮強度
12方向剪切強度
13方向剪切強度
23方向剪切強度
3.建立顯示Explicit計算時間步，在場變量中勾選輸出 SDV和 STATUS.
4.劃分網格，賦給Explicit 3D stress單元類型，邊界條件根據需要設定即可。此處劃分為一個單元，單向加載。建立Job，提交模型前在Job中選擇該子程序，進行計算。
5.查看結果，滿足失效準則后無承載，單元被刪除。
子程序輸出的state1-6為儲存的應變（順序為11 22 33 12 23 13），state7為單元刪除變量，state8-11為Hashin失效判斷系數（0~1）。
接下來簡要介紹該子程序的相關理論
彈性階段總應力與總彈性應變之間的關系為
式中，σ是柯西應力，S0是柔度矩陣，ε是彈性應變。
展開 基于cohesive的三維脆性材料斷裂（abaqus cohesive單元本構） ￥999
1、根據論文《Three-dimensional modeling of fracture in quasi-brittle materials using plasticity and cohesive finite elements》DOI:https://doi.org/10.1007/s10704-021-00514-1 編寫的cohesive單元本構
2、適用于三維模型
3、包含umat以及vumat
4、umat適用范圍小，多個cohesive單元一般采用vumat進行計算
5、軟化曲線為Hordijk和bilinear
ABAQUS損傷斷裂 (例1) 金屬切割或沙柳切割斷裂 ￥26.67
1)該模型模擬了材料在旋轉切割下的損傷斷裂全過程，模型考慮了材料的彈性變形，塑性應變，損傷破壞的標準，損傷演化及斷裂的全過程，并考慮了溫度的影響；
2)模型可用于模擬沙柳切割過程，金屬切割過程及材料的損傷斷裂過程。
材料的疲勞損傷與斷裂 ￥5
材料的疲勞損傷與斷裂
復合材料的斷裂和韌性
復合材料的斷裂和韌性.ppt
材料脆性斷裂有限元模擬的UEL子程序實現方法
一、引言
相場斷裂模型是描述當裂紋尖端狀態達到臨界能量釋放率時發生裂紋擴展現象的工具，被廣泛應用于材料斷裂過程的模擬研究。UEL (User Element Subroutine) 子程序允許用戶自定義單元的切線剛度矩陣及節點力向量，在實現相場斷裂模型時具有靈活性與便利性。UMAT (User Material Subroutine)子程序可以供用戶自定義材料的本構模型，同時可彌補UEL子程序無法可視化的缺陷。本文中的計算結合UEL子程序與UMAT子程序，采用雙層模型進行交互計算材料的斷裂過程。
展開 MARC模擬金屬材料的斷裂破壞
我最近在用MARC模擬金屬材料在楔橫軋過程中的軋制過程的斷裂和破壞，做了一個簡單的二維的模擬，我把文件傳到附件中，有興趣的哥們可以下著看看，有MARC的高手希望能指導一下，我現在只是在初級階段，我的目的是研究三維的模擬。
crack_job1.rar
復合材料疊層結構的拉伸斷裂仿真 ￥800
本案例基于COMSOL軟件中的固體力學模塊的損傷模型模擬了一復合疊層結構在受到兩端拉伸作用下的拉伸變形過程以及斷裂帶生成過程，模擬結果如圖所示：
感興趣的朋友，歡迎合作交流！
在MARC中分析脆性材料斷裂問題
我用的材料是顆粒增強復合樹脂，延伸率只有2％，屬于脆性材料
A：我也正在用MARC做斷裂分析,不知道你算的是什么樣的問題,如果是MODE I 斷裂模型的話,我推薦你使用能量法,MARC會為你算出J積分值的,而且對于脆性材料,其J積分值正好就是能量釋放率.我做的問題是納米涂層結構在MIXED MODE下的斷裂分析,由于摩擦的存在,J積分可能不再保持其積分路徑無關性,所以采用直接法來計算應力強度因子的,也就是用應力或位移來算SIF.聽說MARC中可以在尖端處采用特異單元來擬出應力特異性,但我不會用,所以只好通過REFINE來細分尖端區域(最小單元長度在10E-3量級就差不多了(在MODE I下.采用能量法計算的話,不用細分尖端點).然后通過應力(或位移)結果求K值,并用外插法來求出尖端點處的K值.隨裂紋成長,K和斷裂長度的關系,我用了個笨辦法.就是預先知道斷裂擴展方向的基礎上,在每個裂紋長度下計算出K值,就可以得出K和a的關系曲線了.不過很麻煩,我一共建了14個模型算的,光建模就要花費了很多時間,但結果還算可以.
我再給你介紹一種方法,供你參考.我研究室有個老頭,10年前用MARC算過纖維復合材料的單纖維push-out問題,他在纖維和基體的界面處,用用戶自定義TYING進行約束,然后在UFORM子程序中給出TYING拖開的條件(他是通過實驗測出最大剪應力以后,以此應力作為拖開基準的,不過好象不準),這樣,就可以實現在加載過程中裂紋自動擴展.但大前提仍然是事先預測好擴展方向.我本來想試試,不過我這的MARC好象安裝有問題,子程序的compiler總是不行.
展開 復合材料纖維方向斷裂韌性的測量方法
近期的材料力學性能數據庫中收錄的部分材料中有一些是有纖維方向斷裂韌性的，歡迎大家點擊“材料庫”菜單查閱。
文章轉載自“復合材料力學”微信公眾號，原文鏈接如下：
復合材料纖維方向斷裂韌性的測量方法
激光非接觸發測量材料形變、斷裂
最高精度，線性度0.001%到0.1%，分辨率0.5nm到0.1mm
最大量程，130um-2000mm，最遠可測距離1mm到4000mm
最小尺寸，直徑6mm
最高采樣速度，2kHz到400kHz
最高可耐溫度，2200℃超高溫表面可測
應用
在線檢測
? 平整度監控
? 涂膠高度測量
? 翹曲度監控
位移測量
? 主軸跳動
? 仿生肌肉
形貌測量
? 沖壓\磨損形貌
? 板材厚度
? 材料熱變形
? 鋼軌形狀
振動測試
? 振動臺試驗\風洞試驗
? 沖擊試驗
? 模態分析
上海思信科學儀器有限公司面向全國各大高校、科研單位提供檢測及實驗用高精密儀器。
主營產品包括：激光位移傳感器、色散共焦位移計、高速攝像機、紅外熱像儀、激光測振儀、光學形變測量儀、激光剪切散斑干涉儀；日本YAMATO實驗室通用設備、YAMAOT等離子刻蝕/清洗機、YAMAOT等離子灰化裝置、YAMAOT噴霧干燥機等；各種顯微鏡、內窺鏡。
電話：021-31177311
E-mail：sparkshi@think-foucus.com
展開 