應力強度因子 ansys的案例
在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中啊。
展開 基于單元應力的應力強度因子外推法
表1 基于應力的外推法數據表
圖5 基于應力的外推法計算應力強度因子
以上就是基于應力外推裂尖應力強度因子的入門介紹，后期會更新更多斷裂相關內容，一起加油吧~
xfem輸出應力強度因子
使用ABAQUS中的xfem裂紋，在歷程輸出中選擇輸出應力強度因子，但是顯示出錯，沒找到裂紋前端
應力強度因子手冊
應力強度因子手冊3.rar
應力強度因子手冊1.rar
應力強度因子手冊2.rar
ABAQUS xfem裂紋擴展計算應力強度因子
Delete or suppressthe interaction, or change the crack to allow growth.輸入文件未生成，作業未提交分析
想請教一下，這個應該怎么改，未勾選允許裂紋生長，想要輸出應力強度因子，但是一值提交不上
應用三維有限單元法計算應力強度因子
應用三維有限單元法計算應力強度因子
來源：中國機械工程 作者：林曉斌
摘要 描述了兩種基于有限單元計算面形裂紋應力強度因子的方法，建議了一種創造三維有限單元網格的途徑。計算方法的精度通過和其它解析解或數值解的比較得到了說明。
關鍵詞 應力強度因子 有限元 損傷容限設計 斷裂評定
無論在損傷容限設計還是在缺陷評定階段，工程師們需要知道正在分析的構件中裂紋的應力強度因子，因為判斷含裂紋構件的斷裂，或者計算剩余疲勞壽命大多依賴于這一參量。因此，在斷裂力學發展中，如何求取應力強度因子一直是一個重要的課題。當前已有許多方法可用來計算應力強度因子，較為典型的有解析法、邊界配位法、有限單元法、邊界元素法、體力法、權函數法和線彈簧模型。利用這些方法，大量的應力強度因子解已經獲得，已出版的應力強度手冊［1］中收編了許多典型的解。盡管如此，工程師們仍然會感到自己所需要的應力強度因子解很難找到，這是因為要解決的工程問題往往是一些受復雜載荷的構件，包含的裂紋也往往是一些不規則裂紋。
本文簡單介紹了兩種基于三維有限單元法計算面形裂紋應力強度因子的方法。有限單元法已經成為工程設計分析領域中一個強有力的計算工具，它能模擬非常復雜的構件。基于有限元的應力強度因子計算方法，自然也將具有卓越的工程能力。除了計算方法的介紹以外，還將簡單描述一種簡化網格的生成方法。最后提供了一些所得到的典型應力強度因子解，并和大家熟知的解進行了比較，以說明本文所描述的方法的可靠性。
展開 計算三維問題的應力強度因子的程序
/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) !
展開 基于節點位移的應力強度因子外推法
今天木木給大家分享的是基于節點位移求解應力強度因子，相比于上一期出的基于單元應力求解應力強度因子得出的結果更加接近解析解。這一期包括以下內容：（1）簡要講述INP文件（2）運用最小二乘法進行線性擬合（3）對裂尖數據進行特殊處理。
怎么在abaqus提取xfem下，動態三維應力強度因子？
求解答
斷裂力學—有限寬板含雙邊裂紋的應力強度因子計算 ￥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計算二維穿透裂紋應力強度因子的實例 ￥10
通過ABAQUS計算二維穿透裂紋應力強度因子。理論值為396，有限元計算值為407，誤差小于3%。
一個斷裂分析實例（J積分和應力強度因子）
只有計算應力強度因子才需要。這里只需要保證全局坐標系的X方向與裂紋平行就是了。
csys,0
！這里應該有一個定義path，這里沒有寫出。
Python二次開發-雙孔雙裂紋板計算應力強度因子
通過編寫py腳本文件，實現自動建立雙孔雙裂紋模型，劃分網格，定義seam和crack，控制輸出。
py腳本準備就緒后，在abaqus的PDE中點擊運行，效果如視頻所示。可以實現自動建模，劃分網格，控制輸出。
在腳本中，兩條裂紋的長度，載荷的大小均為變量，可以通過控制變量的輸入，實現參數化分析。
ABAQUS計算三維孔邊角裂紋應力強度因子的實例模型 ￥15
通過ABAQUS，可以計算三維角裂紋的應力強度因子。本實例中對平板孔邊的三維角裂紋進行了模擬。
seam及crack定義如下圖：
網格如下圖：
計算后的位移云圖如下：
對裂尖進行放大觀察：
本實例的難點在于孔邊三維角裂紋的模型的建立，需要經過一系列的布爾操作（merge/cut）得到。
相應的應力強度因子可以在提交job計算完成后，到dat文件中找到。
詳細的模型可參考附件。
分享一個斷裂分析實例（J積分和應力強度因子）
另外，ANSYS里面給了一個路徑項求導的操作：general postproc->path operation ->differentiate 。（differentiate是不是求導，請指教） 那這個東西用來求偏導數不行嗎？ 為什么help 里面要那么來求偏導數。
3、這個實例的建模
這個實例的建模，我是建的1/4模型。1/8模型也可以建出來，但是我對于1/8模型還能不能用對稱邊界條件有懷疑。1/4模型用對稱邊界條件是絕對正確的。另外，對稱邊界條件得到的約束條件在載荷步里面查看到，約束是發生在環向的（柱坐標系下），環向位移約束為零。加約束的時候直接加環向位移為零，也是可以的。
四、 命令流（log文件另附）
log.rar
1、 建模和求解部分（這里的建模網格劃分比較密，可能不是很實用,這里的網格劃分不好，在裂紋尖端第一行單元沒有奇異性，最好還是用kscon 來做）
2、 J積分部分（對于半邊裂紋，如果你已經定義了路徑的話，直接把這部分命令流輸入進去就可以了）
/post1
!local,11,0,0.034,0,0 !這里不應該建立局部坐標系。只有計算應力強度因子才需要。這里只需要保證全局坐標系的X方向與裂紋平行就是了。
csys,0
！這里應該有一個定義path，這里沒有寫出。
3、應力強度因子
（1）方法一、先建立局部坐標系：原點在裂紋尖端，x方向與裂紋平行，Y與裂紋垂直，笛卡爾坐標系。定義路徑，直接點一下菜單路徑就出來了，或者用kcalc就可以了。
（2）方法二、線彈性情況下。先算出J積分然后根據J積分與應力強度因子的關系來求應力強度因子。對于平面應變模型，J積分＝應力強度因子的平方×（1－泊松比×泊松比）/彈性模量
全尺寸裂紋模型前面所建立的模型都只有裂紋的半邊。
展開 