abaqus裂紋j積分的案例
ABAQUS直裂紋、斜裂紋圍道積分計算裂紋尖端J積分
之前算過一個關于裂紋擴展的問題，當時創建裂紋選擇的是contour intergral，后來又有人咨詢我裂紋尖端J積分的計算問題。我才恍然大悟，其實圍道積分方法還是適用于計算裂紋尖端在某時刻的J積分，至于動態擴展問題，還是交給XFEM吧（雖然也不太好）。
計算了幾種情況下的裂紋尖端J積分，包括直裂紋、斜裂紋以及裂紋尖端傾斜等三種情況。
部分試件的應力分布及J積分結果如圖所示：
展開 3D-Crack的一個簡單的彈性例子，算J積分和裂紋擴展的
我是必葫蘆畫瓢作的，不當之處，請多多指教
拿出來大家一起學習吧
呵呵
所用軟件：MSC.Marc2005
以下為menten源代碼：
*define h 4
*define b 8
*define l 15
*define a1 1.5
*define a2 1
*define a2 1.5
*define a1 1
*define force -5e5
*define moment -1e5
*define nexpand 22
*define da 0.5
*add_solids
0 0 0 b h -l
*set_solid_type cylinder
*add_solids
0 0 0 0 0 -1 1.5 1.5
*move_reset
*set_move_scale_factors
a2/a1 1 1
*move_solids
2
# | End of List
*move_reset
*set_move_translations
b h 0
*move_solids
2
# | End of List
*fill_view
*solids_subtract
1 2
# | End of List
*reset_view
*fill_view
*rot_model_cspace_x_for
*rot_model_cspace_x_for
*rot_model_cspace_y_rev
*rot_model_cspace_y_rev
*rot_model_cspace_z_rev
*fill_view
*convert_solid_faces_surfaces
all_existing
*remove_solids
all_existing
*set_curve_div_type_fix_avgl
*set_curve_div_type_fix_avgl
展開 