abaqus应力强度因子计算例子.pdfVIP

如何使用ABAQUS 计算应力强度因子 Simwefanhj （fanhjhj@ ） 2011.9.9 问题描述：以无限大平板含有一贯穿裂纹为例，裂纹长度为 10mm （2a ）， 2 在远场受双向均布拉应力 ＝100N/mm 。按解析解，此 I 型裂纹计算出的应力 强度因子K I  a ＝396.23 （N.mm-3/2 ） 以下为使用ABAQUS6.10 的计算该问题的过程。 第一步：进入part 模块 ①建立平板part （2D Planar;Deformation;shell ），平板的尺寸相对于裂纹足够大， 本例的尺寸为100×50 （mm ）。 ②使用Partation Face:sketch 工具，将part 分隔成如图 1 形式。 图1 第二步：进入property 模块 ①建立弹性材料； ②截面选择平面问题的solid ，homogeneous ； ③赋予截面。 第三步：进入Assembly 模块 不详述。需注意的是：实体的类型（instance type ）选择independent 。 第四步：进入mesh 模块 除小圈内使用CPS6 单元外，其它位置使用CPS8 单元离散 （图2 ）。裂纹尖端的 奇异在interaction 模块中 （图4 ）考虑。 四边形单元CPS8 三角单元CPS6 图2 第五步：进入interaction 模块 ①指定裂纹special/creak/assign seam ，选中示意图3 中的黄色线，done ！ ②生成裂纹crack 1，special/crack/create,name ：crack 1，type: contour integral. 当提示选择裂纹前端时，选则示意图的红圈区域，当提示裂纹尖端区域时选择红 圈的圆心，用向量q 表示裂纹扩展方向 （示意图3 绿色箭头）。用同样的方法建 立crack 2 （示意图3 中的蓝色区域）。 special/crack/edit ，对两个裂纹进行应力奇异的设置，如图4 所示。 Crack 1 Crack 2 Assign seam 裂纹方向矢量 裂纹方向矢量 图3 图4 第六步：进入step 模块 ①在initial 步之后建立static ，general 步。 ②在output/history output requests/create/ 中创建2 个输出变量，如图5 所示。（分 别定义crack 1 和crack2 ）。其中Number of contours 的意义是： SIF 是通过J 积分算出的，J 积分本来是线积分，在有限元计算中不方便，后来有人研究可以 用环形区域积分代替，Abaqus 里J 积分的计算就是这样做的，积分围数 1 就是 指计算J 积分时取裂尖处奇异单元外第一圈单元作为环形区域算J 积分，不同积 分围数算出的应力强度因子应该相差不多，如果相差太多，就说明算的有问题。 图5 第七步：进入load 模块 定义位移和荷载边界，如图6 所示。 图6 第八步：进入job 模块 提交计算。Mises 应力分布见图7 ，在.dat 文件中 （图8 ）查看应力强度因子。 图7