斷裂能
關注創建者：老司機1 創建時間：2019-10-05
斷裂能的實例教程
絕大部分論文僅引用了這四個公式就去做漸進失效分析了，是不完整的，這僅僅能預測損傷的起始，還缺乏損傷起始以后的剛度退化。
與損傷起始相關的是6個強度值，定義如下，6個參數分別是：縱向拉伸強度XT，縱向壓縮強度XC，橫向拉伸強度YT，橫向壓縮強度YC，縱向剪切強度SL，橫向剪切強度ST。
圖2
03 基于能量的損傷演化
定義損傷演化參數的界面如下，輸入四種失效模式對應的斷裂能:縱向拉伸斷裂能Gftc,縱向壓縮斷裂能Gfcc,橫向拉伸斷裂能Gmtc,橫向壓縮斷裂能Gmcc。(斷裂能，即臨界應變能釋放率，是裂紋擴展單位面積所需要的能量，一般單位取為N/m，J/m2，N/mm或mJ/mm2。
展開 Type選擇判斷損傷的演化的方式：斷裂能和位移； softing指定損傷后的軟化方式；
混合模式行為：對以位移為判斷損傷依據時，只有兩種：直接指定損傷時的位移和表格形式;對于與斷裂能為判斷損傷依據時，還有BK定律和power law兩種，均需要指定法向和切向的斷裂能，其中BK定律中power給定有效數值即可，power law模式需要指定指數power。各項同性時，三種斷裂能相等，對指數不敏感。
按照下圖定義損傷的演化：
即能量釋放率（最大位移）達到給定的斷裂能（位移）認為完全損傷，裂紋將繼續向前擴展。
3）定義截面屬性
選擇材料和響應，本文選擇牽引開裂響應，使用節點坐標確定膠粘單元初始厚度。對于牽引開裂響應，默認厚度為1，其他默認為基于節點坐標的幾何厚度。如果采用指定為單位厚度，則創建材料時彈性模量需要進行縮放。采用單位厚度，則名義應變與開裂值相同，因為對于牽引開裂準則，彈性模量當罰剛度來理解，彈性響應通過En控制。
4）采用通用靜力學分析步，固定分析時間，初始增量步0.05，總時間1，左端完全固定，右下邊施加向右的位移0.05.
5)劃分網格
膠粘單元部分需采用結構或掃略技術，且形狀為四邊形，選擇cohesive單元族。最后單元類型為COH2D4（四節點的2維膠粘單元）。
6)提交計算，查看結國，截取裂紋尖端名義應變圖如下：
可以看到圖中裂紋尖端兩點變形前后后距離差為5.73036e-3,初始距離1.99986e-3,計算得到名義應變與圖中名義應變為2.86532相等。
展開 “我們需要了解鞋子斷裂的原因，因此，作為設計師，我們可以將材料發揮到極致，制作出更高更細的更復雜的鞋跟，并確信它們可以隨時穿著，并且經久耐用，”Tacchificio Villa Cortese 的研發經理DavideCarminati 說。
通過Marc虛擬模擬鞋子的性能測試，使我們可以嘗試更多的設計來調整和平衡工藝和性能，讓我們對產品的質量充滿信心，并在生產單個零件之前靈活地進行更多實驗以制造出更好的產品。
從技術角度來看，高跟鞋可能面臨幾個關鍵問題：
? 靜態負載下的損壞或斷裂，例如用戶在站立不動時
? 鞋跟固定在鞋子上時發生的斷裂
? 易受金屬嵌件損壞；作用于鞋跟的塑料部件，因為插入物會導致應力集中現象
? 由于疲勞導致的斷裂，無論是由于缺陷導致破損，還是由于穿鞋時可能發生的粘彈性效應
仿真使 Tacchificio Villa Cortese 了解如何改進他們的產品，以及如何優化鞋跟的材料。在使用 Marc 之前，他們無法在原型制作之前確定產品形狀（幾何）中的問題，有時需要對多個原型進行測試以確定正確的幾何形狀。
通過仿真，他們可以設計優化鞋跟的幾何形狀，在項目初期平衡所有的設計性能，并且避免了不好的設計。
對于當今的高級時裝公司而言，可持續性是一個重要趨勢。即使根據定義，高級時尚比快速時尚更具可持續性，但仍有更多工作要做：引入合適的材料和工藝變得越來越緊迫。
此外，在瞬息萬變且競爭激烈的時尚行業，按時為評論家和時裝秀提供新款式是一項大生意。
展開 緒論
斷裂是混凝土材料破壞的主要模式。可靠、高效的混凝土斷裂模型在橋梁、隧道、大壩等土木工程結構的安全評估中發揮著重要作用。對混凝土斷裂的研究,尤其對其裂紋萌生和擴展的研究,引起了國內外學者越來越多的關注。混凝土斷裂的數值模擬與斷裂理論、物理試驗相互印證與補充,并隨著科技發展不斷地提高著混凝土斷裂問題模擬的準確性。近年發展起來的斷裂相場法,通過場變量的自動演化獲取裂紋路徑,可方便地模擬出裂紋的動態擴展過程。因此本案列將采用基于<a href="/major/<a href="/major/ABAQUS的斷裂相場模型實現對混凝土斷裂問題的模擬分析并探討該模型的工程實際適用性
理論基礎
相場法是一種以經典熱、動力學理論為基礎,由耦合的非線性的力平衡方程和相場梯度型演化方程組合而成的唯象方法。該方法引入一組場變量來描述結構的相變過程。與銳界面法中場變量的不連續性相反的是,相場法中場變量在界面區域具有連續性,可以用來描述材料初始時和完全破壞之間的平滑過渡。相場變量能分成保守的場變量與非保守的場變量兩種,總量在物體結構演化中保持不變的為保守的場變量,如原子和電荷的濃度場;總量在物體結構演化中為不守恒的并從0到1變化的是非保守的場變量,如馬氏相變。
Frankfort和Marigo基于能量最小化原理提出了Griffith理論的變分形式。描述斷裂的相場法中材料勢能分為兩部分，彈性應變能和表面能，分別對應于完好相和斷裂相。Griffith理論的泛函形式可以表達為:
其中是對稱的小應變張量，代表裂紋面，Ω為求解區域。斷裂問題系統自由能由彈性應變能(等號右邊第一項)和斷裂表面能(等號右邊第二項)構成，裂紋的擴展受自由能最小化原理控制。通過求能量泛函的極值可以獲得材料系統的控制方程。
采用有限寬度的彌散區域來近似表征離散裂紋面，如圖 1所示。
展開 1 引言
本部分介紹來自于《斷裂相場法》書籍。
“1998年Francfort和Marigo根據Griffith脆性斷裂理論，提出了一種斷裂力學變分原理，他們以結構內可能的位移場和裂紋面作為自變量，將變形能與斷裂面之和定義為結構總能量，并且認為真實的位移場與裂紋面使得該總能量最小。然而在數值模擬中將離散的裂紋面作為未知量來求解是非常困難的。因此2000年Bourdin等提出了一種相場模型，其中引入了一個連續的標量場，即相場，來近似地描述裂紋。相場值為1和0分別代表材料完全破壞和完好兩種極限狀態，而它們之間的值代表了一種損傷狀態，并且裂紋的彌散程度由相場特征寬度來控制，其值越大彌散寬度越大，反之則越小。然后通過一個與相場相關的裂紋面密度泛函來重構結構內的斷裂能，并將因損傷而退化的變形能與重構的斷裂能代入Francfort-Marigo變分原理就得到了相場模型的基本列式。相場模型中的自變量為兩個連續變化的場，即位移場和相場，因此它可以很方便地由不同數值方法實現。直觀來看，相場模型將一個結構內裂紋萌生與演化問題，轉化為了一個多場耦合情況下求最小能量的優化問題，因此它可以用于直接求解（例如分叉、交叉、融合、扭結等）復雜斷裂問題，而不需要額外的裂紋路徑追蹤方法。”
2 理論
將系統的總勢能表示為如下兩項：
式中第一項能量為：
考慮損傷帶來的退化，彈性能的表達式為：
式中
k為一個小值，用于防止數值不穩定現象。另一項斷裂能為：
因此代入具體表達式可將系統總勢能表達為：
對上述能量進行一階變分可得：
即可得弱形式方程為：
具體外力虛功為：
式中本構方程為：
該弱形式方程是后續推導有限元方程的基礎。同時，通過弱形式方程也可推導得到強形式的控制方程，即位移場和相場的控制方程。
展開 
斷裂能的相關專題、標簽、搜索
斷裂能的最新內容
由于模擬過程會牽扯到鍵的形成和斷裂，因此不能對鍵長進行約束，shake要設置為0。
Abaqus的強大之處在于其引入了斷裂能正則化機制*DAMAGE EVOLUTION, TYPE=ENERGY。這一機制徹底改變了傳統材料卡片隨網格尺寸變小而急劇變“脆”的網格敏感性缺陷，使得能量耗散成為一個相對客觀的物理不變量。
材料系統與參數</strong></p><p class="ql-align-justify"> 文獻中的材料為 T700/M21，其層內力學參數、強度值、斷裂能以及界面參數均見于文獻表 3（亦為插件預置值）。層間損傷演化為二次應力準則與 B?K 混合模式能量準則。這些參數構成插件中 T700 材料數據庫的核心。
</p><p>三、與實驗對比</p><p><img src="https://img.jishulink.com/msimage/202601/9395db4701eba46fb3afcc7e0ef50e51.png"></p><p>圖5兩種方法的實驗數據對比：（a）界面剪切剛度；（b）極限剪應力；（c）界面斷裂能</p><p>圖5中與實驗結果的對比主要用于驗證本文提出的反演方法在考慮熱效應條件下對實際試驗數據的適用性和準確性
同時也開始了苦苦探求：這個結構在斷裂之前到底能堅持多久？
直到19世紀中期，德國工程師奧古斯特·沃勒通過試驗發現：應力幅度越小，材料能承受的循環次數就越多。
進一步，他以應力幅度（S，Stress）為縱軸，以循環次數（N，Number of Cycles）為橫軸，繪制了著名的S-N曲線。
一百多年后的今天，S-N曲線依然是結構疲勞分析最基礎和最常用的方法。
</p><p><img src="https://img.jishulink.com/msimage/202507/a4d0a8badde65f417c7a863d5df20510.png"></p><p>損傷演化有兩種定義模式，一種是基于斷裂能，另外一種是基于相對位移displacement。
其核心功能包含三部分：首先基于正交各向異性彈性本構更新應力，通過材料屬性計算剛度矩陣并響應應變增量；其次實現彈塑性修正，采用J2流動理論判斷屈服狀態，通過牛頓迭代求解塑性變形并更新應力；最后建立漸進損傷模型，分別針對纖維方向（拉伸/壓縮失效）和基體方向（通過180°平面搜索臨界斷裂面）定義損傷初始判據，結合斷裂能與特征長度控制損傷演化過程。
平面應力脆性斷裂相場AT2模型10個月前
另一項斷裂能為：
因此代入具體表達式可將系統總勢能表達為：
對上述能量進行一階變分可得：
即可得弱形式方程為：
具體外力虛功為：
式中本構方程為：
該弱形式方程是后續推導有限元方程的基礎。同時，通過弱形式方程也可推導得到強形式的控制方程，即位移場和相場的控制方程。
對于混凝土，需明確其彈性模量、泊松比、抗拉強度、密度、斷裂能等力學參數，以及導熱系數、比熱容等熱物理參數（若涉及熱 - 力耦合分析）。
對于纖維材料，要確定其彈性模量、抗拉強度、密度、直徑和長度等參數。此外，還需定義纖維與混凝土之間的界面屬性，如界面粘結強度、界面摩擦系數等，以準確模擬纖維與混凝土之間的相互作用。?
3、 XFEM 原理應用
基于 XFEM 的基本原理，在模型中引入裂縫單元。
3.15
3.29
3.17
3.11
1.0%-ASFRC
4.21
3.69
3.50
3.36
斷裂能
