應力強度因子計算的案例
應用三維有限單元法計算應力強度因子
應用三維有限單元法計算應力強度因子
來源：中國機械工程 作者：林曉斌
摘要 描述了兩種基于有限單元計算面形裂紋應力強度因子的方法，建議了一種創造三維有限單元網格的途徑。計算方法的精度通過和其它解析解或數值解的比較得到了說明。
關鍵詞 應力強度因子 有限元 損傷容限設計 斷裂評定
無論在損傷容限設計還是在缺陷評定階段，工程師們需要知道正在分析的構件中裂紋的應力強度因子，因為判斷含裂紋構件的斷裂，或者計算剩余疲勞壽命大多依賴于這一參量。因此，在斷裂力學發展中，如何求取應力強度因子一直是一個重要的課題。當前已有許多方法可用來計算應力強度因子，較為典型的有解析法、邊界配位法、有限單元法、邊界元素法、體力法、權函數法和線彈簧模型。利用這些方法，大量的應力強度因子解已經獲得，已出版的應力強度手冊［1］中收編了許多典型的解。盡管如此，工程師們仍然會感到自己所需要的應力強度因子解很難找到，這是因為要解決的工程問題往往是一些受復雜載荷的構件，包含的裂紋也往往是一些不規則裂紋。
本文簡單介紹了兩種基于三維有限單元法計算面形裂紋應力強度因子的方法。有限單元法已經成為工程設計分析領域中一個強有力的計算工具，它能模擬非常復雜的構件。基于有限元的應力強度因子計算方法，自然也將具有卓越的工程能力。除了計算方法的介紹以外，還將簡單描述一種簡化網格的生成方法。最后提供了一些所得到的典型應力強度因子解，并和大家熟知的解進行了比較，以說明本文所描述的方法的可靠性。
展開 斷裂力學—有限寬板含雙邊裂紋的應力強度因子計算 ￥19.89
σ
σ
圖 C-3-4（取 1/4 研究）
應力強度因子的理論解： KI ? ???a ??
1、 計算不同載荷下的應力強度因子的數值解和理論解，列表給出各種情況的數值解和理論解，并畫比較圖。
厚度t ? 1mm ， a ? 5mm,b ? 80mm,c ? 200mm ，?? 25MPa , 26MPa , 27MPa , 28MPa ,
29MPa
2、 計算不同板寬時應力強度因子的數值解和理論解，列表給出各種情況的數值解和理論解，并畫比較圖。
載荷?? 10MPa ， a ? 11mm, c ? 200mm 。 每個人都計算b / a ? 4、4.2、4.4、4.6、4.8 時的應力強度因子。
3、 計算不同板長時應力強度因子的數值解和理論解，列表給出各種情況的數值解和理論解，并畫比較圖。
載荷?? 10MPa ， a ? 6mm,b=76mm 。每個人都計算c / b ? 2.2、2.4、2.6、2.8、3 時的應力強度因子。
展開 ABAQUS xfem裂紋擴展計算應力強度因子
Delete or suppressthe interaction, or change the crack to allow growth.輸入文件未生成，作業未提交分析
想請教一下，這個應該怎么改，未勾選允許裂紋生長，想要輸出應力強度因子，但是一值提交不上
在ANSYS中計算裂縫應力強度因子的技巧
在ANSYS中計算裂縫應力強度因子的技巧
裂縫應力強度因子用ANSYS中怎么求呀。另外，建模時，裂紋應該怎么處理呀，難道只有畫出一條線嗎？
首先說一下裂紋怎么畫，其實裂紋很簡單啊。只要畫出裂紋的上下表面（線）就可以了，即使是兩個面（線）重合也一定要是兩個面（線）；如果考慮道對稱模型就更好辦了，裂紋尖點左面用一個面（線），右邊用另外一個面（線），加上對稱邊界約束。
再說一下裂尖點附近網格的劃分。ansys提供了一個kscon的命令，主要是使得crack
tip的第一層單元變成奇異單元，用來模擬斷裂奇異性(singularity)。當然這個步驟不是必須的，有的人說起用ansys算強度因子的時候就一定要用奇異單元，其實是誤區（原因下面解釋）
好了，回到強度因子的計算。其實只要學過一些斷裂力學都知道，K的求法很多。就拿Mode
I的KI來說吧，Ansys自己提供了一個辦法（displacement extrapolation）
，中文可能翻譯作“位移外推”法，其實就是根據解析解的位移公式來對計算數據進行fitting的。分3步走，如果你已經算完了：
第一步，先定義一個crack-tip的局部坐標系，這是ansys幫助文件中說的，其實如果你的裂紋尖端就是整體坐標原點的話，而且你的x-axis就順著裂紋，就沒有什么必要了。
第二步，定義一個始于crack-tip的path,什么什么？path怎么定義？？看看幫助吧，在索引里面查找fracture
mechanics，找到怎么計算斷裂強度因子。（my god,我這3步全是在copy幫助中的東東啊）。
第三步，Nodal
Calcs>Stress Int Factr ，別忘了，這是在后處理postproc中啊。
展開 
計算三維問題的應力強度因子的程序
/COM,ANSYS MEDIA REL. 60 (090601) REF. VERIF. MANUAL: REL. 60
/VERIFY,VM143
*CREATE,FRACT,MAC
/NOPR
NSEL,ALL
*GET,N,NODE,,NUM,MAX ! CURRENT MAXIMUM NODE NUMBER
CMSEL,S,CRACKTIP ! SELECT THE TIP NODES
ESLN ! ANY ELEMENTS ATTACHED
*GET,ELMAX,ELEM,,NUM,MAX ! CURRENT MAXIMUM ELEMENT NUMBER
*DO,IEL,1,ELMAX ! LOOP ON MAX ELEMENT
ELMI=IEL
*IF,ELMI,LE,0,EXIT ! NO MORE SELECTED
*GET,ELTYPE,ELEM,ELMI,ATTR,TYPE ! GET ELEMENT TYPE
*IF,ELTYPE,NE,ARG1,CYCLE ! CHECK FOR SELECTED ELEMENT
N3 = NELEM(ELMI,3) ! GET NODE 3 (K)
*IF,NSEL(N3),LE,0,CYCLE ! IT MUST BE SELECTED
N7 = NELEM(ELMI,7) ! GET NODE 7 (L)
*IF,NSEL(N7),LE,0,CYCLE ! IT MUST ALSO BE SELECTED
N1 = NELEM(ELMI,1) ! GET NODE 1 (I)
N2 = NELEM(ELMI,2) ! GET NODE 2 (J)
N5 = NELEM(ELMI,5) ! GET NODE
展開 分享一個通過ABAQUS計算二維穿透裂紋應力強度因子的實例 ￥10
通過ABAQUS計算二維穿透裂紋應力強度因子。理論值為396，有限元計算值為407，誤差小于3%。
Python二次開發-雙孔雙裂紋板計算應力強度因子
通過編寫py腳本文件，實現自動建立雙孔雙裂紋模型，劃分網格，定義seam和crack，控制輸出。
py腳本準備就緒后，在abaqus的PDE中點擊運行，效果如視頻所示。可以實現自動建模，劃分網格，控制輸出。
在腳本中，兩條裂紋的長度，載荷的大小均為變量，可以通過控制變量的輸入，實現參數化分析。
ABAQUS計算三維孔邊角裂紋應力強度因子的實例模型 ￥15
通過ABAQUS，可以計算三維角裂紋的應力強度因子。本實例中對平板孔邊的三維角裂紋進行了模擬。
seam及crack定義如下圖：
網格如下圖：
計算后的位移云圖如下：
對裂尖進行放大觀察：
本實例的難點在于孔邊三維角裂紋的模型的建立，需要經過一系列的布爾操作（merge/cut）得到。
相應的應力強度因子可以在提交job計算完成后，到dat文件中找到。
詳細的模型可參考附件。
斷裂／裂縫／裂紋／應力強度因子和Ｊ積分計算實例（原創，如轉載，請注明出處）
材料：復合材料
分析類型：斷裂力學
技術難點：斷裂 應力強度因子 Ｊ積分
完成人：技術鄰ANSYS專家
網址：http://www.yqgqt.org.cn/content/other/402981
模擬過程：
斷裂力學計算應力強度因子和Ｊ積分
圖１ 裂縫模型
圖２ Ｋ１應力強度因子
圖３ Ｋ２應力強度因子
圖４ Ｋ３應力強度因子
圖５ Ｊ積分
圖６ 裂縫前沿集中應力
基于單元應力的應力強度因子外推法
表1 基于應力的外推法數據表
圖5 基于應力的外推法計算應力強度因子
以上就是基于應力外推裂尖應力強度因子的入門介紹，后期會更新更多斷裂相關內容，一起加油吧~
【11月23-26日 南京】結構斷裂力學與裂紋擴展壽命計算專題培訓
結構斷裂力學與裂紋擴展壽命計算專題培訓
14個實例模型課程中人手一機操作指導
案例01、高壓輸氣管道內壓作用下的表面橢圓裂紋應力強度因子計算
案例02、壓力容器表面任意形狀裂紋的J積分計算
案例03、含V形缺口三點彎曲能量釋放率計算
案例04、三點彎曲PMMA半圓板T應力計算
案例05、軸對稱結構表面裂紋的材料構型力計算
案例06、中心裂紋拉伸作用蠕變過程中的C*積分計算
案例07、基于VCCT技術的粘膠脫粘過程數值計算
案例08、涂層-基體界面裂紋擴展計算
案例09、產品的跌落沖擊連接界面失效計算方法
案例10、三點彎曲平面模型裂紋擴展計算
案例11、三維平板邊裂紋的裂紋擴展計算
案例12、周期對稱循環載荷作用下的裂紋擴展壽命計算
案例13、T型焊接件焊縫表面裂紋的熱-結構耦合應力強度因子計算
案例14、復合材料的分層裂紋擴展
課程差異化
1、專注CAE仿真計算，17年大量的工程案例積累
2、6000多學員反饋、提煉的精選內容與實例，形成的版權課程體系
3、有自己的超算中心，有豐富的項目案例庫
專家團隊
團隊12年專注CAE技術工程應用方法，為客戶提供系統的產品質量提升和優化的技術方案，具備上百例的工程問題解決經驗，熟悉CAE技術應用過程中的難點與關鍵點，團隊提供有價值的CAE技術服務。
課程大綱
時間地點
2018年11月23日-2018年11月26日 南京（第一天報道，上課三天）
培訓費用
3980元/人,住宿可統一安排，費用自理。
展開 
xfem輸出應力強度因子
使用ABAQUS中的xfem裂紋，在歷程輸出中選擇輸出應力強度因子，但是顯示出錯，沒找到裂紋前端
應力強度因子手冊
應力強度因子手冊3.rar
應力強度因子手冊1.rar
應力強度因子手冊2.rar
基于節點位移的應力強度因子外推法
今天木木給大家分享的是基于節點位移求解應力強度因子，相比于上一期出的基于單元應力求解應力強度因子得出的結果更加接近解析解。這一期包括以下內容：（1）簡要講述INP文件（2）運用最小二乘法進行線性擬合（3）對裂尖數據進行特殊處理。
計算平板裂紋強度因子
結果分析與結論：
改變平板長寬比，理論計算與有限元分析匯總如下，α代表應力強度因子，a、b分別為平板長寬，可得出規律：平板長寬比越不平衡，其強度因子越小。
