三维裂纹J积分求解

一、建part

草图尺寸100，建草图为矩形（-10,50）、（50，-10）

拉伸厚度10，ok，得到六面体板

二、材料属性，具体参数自己设置吧，一步步都应该清楚。

3、装配（creat instance）

装配要用第二个（独立的），装配之后要切出一条缝，为后面定义seam，用下图右边的那个

这个是选择分割草图的一条边这里应该选择红色面的右边， 这个是分割需要的草图，黄色的线就是分割线，然后，分割tool-partition

选择cell-extrude/sweep edges

选择红色的圈

这个选择

选中的这个然后选择z轴，看好方向如果方向正确就选ok，如果不正确就选flip然后再选

按住shfit选择涂黄的两条线（就是裂纹线）

（1）先定义seam，special-creat-assign选择红色的区域，就定义了seam

（2）定义crake，special-creat-crake，用围线积分继续，这里选择

的时候容易选择的只是一个面

所以先隐藏那个小圆面，选择的那个就是隐藏的命令，隐藏小圆面后就可以选择crake front了。

Crake Front 亮红色的就是选择的裂纹线

学裂纹尖端有奇异性，所以这儿要按下图设置，不懂的看理论吧。此时crake定义好了

5、荷载步

按静力学分析，继续默认设置就行，谈后在历史输出步中要编辑

Domain这个要选择之前定义的crake-1下面的围线数我一般用的是8，这个我感觉8左右就行，但是不能就2,3个，也不必要太多（路径无关性）。

6、load

这个我简单加载两端均布拉力，做好加一个约束边界条件，这里我加了一个右边的UZ

7、划分网格，

感觉断裂最麻烦的就是划分网格了，划分网格之前必须对实体分割能比较规则的，而且在裂纹尖端需要时放射状，其实尖端是四边形也能算，但是所有的例子都是放射状的。

因为划分网格比较麻烦我就补一步步的说了那个之前做的例子来示意一下

最好是按照下图进行分割，

矩形里面最好再分一个圆这样可以有效保证网格效果比较好。

分割之后就是布种子，布种子时可以约束选择第三个。并且在厚度方向的线也要布种子，最好都布，而且数目相同，要不然会划分出来的网格不是一层一层的，求J积分的时候就会提示错误不能求J积分，但是不影响力学分析。

然后就是mesh contral黄色区域设置扫略，绿色区域设置结构话，再选择单元类型，这里用二阶单元，得到的网格就是下图这样。当然了如果会用hypermesh的可以用hypermesh，本人还不会hypermesh所以只能这么笨的方法划分网格了。然后就可以提交计算了。

辛苦了一中午为了保护自己的权益，写个自己的名字吧！

2015年5月13日星期三于大连理工大学

张鹏