《Abaqus裂纹模拟心得(Contour Integral不是XFEM)》由会员分享，可在线阅读，更多相关《Abaqus裂纹模拟心得(Contour Integral不是XFEM)（2页珍藏版）》请在金锄头文库上搜索。
1、Abaqus 裂纹模拟心得（Contour Integral 不是 XFEM）最近由于项目需要，做 了一些裂纹相关的模拟，在此把一些心得体会贴到论坛上与大家分享，如有不当之处，欢迎大家指正！本帖主要侧重于介绍裂纹定义过程中各个选项的意义，具体的操作过程论坛里已经有高手做了很好的教程，至于断裂力学理论推荐大家看一下沈成康写的断裂力学一书。裂纹的定义和输出需要用到 interaction 模块和 step 模块：一、Interaction 模块1.1 预制裂纹（步骤：菜单/special/crack/assign seam）注意：并不是作裂纹分析都要定义 seam，如果你的裂纹不是一条缝，而是一个
2、缺口，则不需要 assign seam，直接走下一步（定义裂纹）就行。1.2 创建裂纹（步骤：菜单/special/crack/create，type：contour integral）crack front：crack front 是用来定义第一围线积分的区域，2D 下我们可以选择包围裂尖点的面，3D 则选择包围裂尖线的面；另外还有一种定义 crack front 的方法，就是直接选择裂尖点（2D ）或裂尖线 3D），用这个方法定义 crack front 不需要再定义下一步的 crack tip/line，比较简便，两种方法算出的结果没有明显的差别，其实只是影响积分路线的问题，但是 J 积
3、分值是路径无关的，看个人喜好吧crack tip/line：这个比较好理解就是裂尖点（2D ）或线（3D ），如果我们在上一步中用方法二定义 crack front，这一步就直接跳过了crack extension direction（定义裂纹扩展方向）：这里定义的其实是一个虚拟的裂纹扩展方向，定义了这个参考方向后，我们才能通过输出的角度判断裂纹扩展方向，可以通过两种方法：o q vector：输入一个方向，用来作为计算裂纹的扩展方向的参考方向；o normal to crack plane：crack plane 表示裂纹的对称面（当裂纹在一个平面内时，可能需要分开定义多个裂纹），这种方法下
4、我们只需定义裂纹面的法线方向，通过 （t 表示裂纹尖端的切线）, 会在每个节点得出一个 q 方向（如下图）；o 注意：q 的方向对输出的应力强度因子，J 积分等都会有影响，一般情况下， q 最好在裂纹平面内，且垂直于裂尖线的切线，否则算出的应力强度因子，J 积分值等等在不同围线积分中会差别较大。二、step 模块定义好了裂纹相关参数后，我们需要返回 step 模块定义输出变量：步骤：菜单/output/history output requests/create，domain ：crack，可以输出的值包括：J-integral，Ct-integral，stress intensity fac
5、tor，T-stressJ-integral ：用于应变率无关材料的准静态分析过程，包括线弹性，非线性弹性，弹塑性材料（单调加载工况）的静态分析。J-integral 的优点是和积分路径无关，从而可以避开尖端塑性区的影响。Ct-integral：用于蠕变分析（一般较少用到）应力强度因子：o 只能用于分析线弹性材料，表示裂纹尖端的应力场强度；o 有三个应力强度因子 K1，K2，K3，分别对应于张开型，滑开型和撕开型裂纹的应力强度因子o 在输出应力强度因子时也会输出一个 J-integral 值，因为算法不同，这个值和直接输出的 J-integral 会略有差异；o 方向判断准则： Maximum tangential stress（在 dat 文件中输出的 MTS 值就是通过这个准则算出的裂纹扩展方向），Maximum energy release rate（dat 中用 MERR 表示），K2=0（dat 中的 K20） T-stress（表示裂纹尖端平行于裂纹面方向的应力）
