用Abaqus 计算 J- Integral的INP 最简结构-本人原创

wrj001

*Node 节点
      1,           0.,         -50.
...
*Element, type=CPS8,elset=SET1 裂纹区单元及其它单元
1,  886,  420,  876, 2213, 2214, 2215
...
*Solid Section, elset=SET1, material=STL 材料定义
,
*Nset, nset=SET2 裂纹区节点
   14,   15,   17,  300,  301,  302,  303,  304,  305,  306,  307,  308,  309,  310,  311,  312
...
*Nset, nset=SET4 裂纹尖
17,
*Nset, nset=SET8 约束节点
    1,    7,    8,  162,  163,  164,  165,  166,  167,  168,  169,  170,  171,  172,  173,  174
......
*Nset, nset=SET9 载荷节点
......
*Material, name=STL 材料定义
*Elastic
2.1e+11, 0.3
*Boundary 边界
SET8, 1, 1
SET8, 2, 2
*Step, name=Step-1 载荷步
*Static
1., 1., 1e-05, 1.
*Cload
SET9, 2, 121.951
*Restart, write, frequency=0
*Print, solve=NO
*Output, field, variable=PRESELECT
*Output, history, frequency=0
*Contour Integral, crack name=Crack-1_CRACK-1, contours=4, crack tip nodes
SET2, SET4, -1., 0., 0.
*End Step
注：
此为本人原创，其它如*Heading 均非必需，已删除。

此处注意，裂纹单元虽然看起来是三角形单元，但使用的是CSP8。这些单元的一个边长为0，即三个点在一个位置。另两边的节点移动到距裂尖1/4边长处。

wrj001

J积分结果见，生成的DAT 文件。

wrj001

在以上分析中，单元可使用CPE8单元。