04_面状裂纹的虚拟裂纹闭合法.pptVIP

Fracture Mechanics 华中科技大学船海学院 袁锐 Fracture Mechanics 华中科技大学船海学院 袁锐 * HUST-XD-B-0111-2014 面状裂纹的虚拟裂纹闭合法 * * 4.1 面状裂纹的基本计算公式 实际裂纹前端沿虚拟裂纹扩展线上的应力分量 虚拟裂纹尖端后部沿虚拟裂纹扩展线上的张开位移分量 裂纹前沿 裂纹面 图4.1 面状裂纹的裂纹前沿及局部坐标系 * 4.1 面状裂纹的基本计算公式 1 ΔA 3 4 裂纹前沿 * 4.2 面状裂纹的哑节点断裂单元 假设裂纹位于（X,Y）平面内，且其局部的 轴与总体坐标的Z轴一致，如图4.3所示。这种断裂单元需要9个节点。同针对线状裂纹的断裂单元一样，只有位于裂纹尖端的节点1和节点2对单元的刚度矩阵产生影响。节点3和节点4用来计算裂纹张开位移。其他节点（节点5至节点9）用来确定局部坐标系和计算虚拟裂纹扩展面积。在ABAQUS的定义单元子程序中，该9节点的完整位移矢量为： 1 2 3 4 5 6 7 8 9 裂纹前沿 1-2影响单元刚度矩阵 3-4计算裂纹张开位移 5-9确定局部坐标系和虚拟扩展面积 为了计算裂纹尖端处的节点力，在节点1和节点2之间放置强弹簧（刚度值很大）。因此，参与刚度矩阵的位移矢量为： 总体坐标轴方向上的节点力分量为： 1 2 * 单元的刚度矩阵为： 1 2 3 4 8 节点3和节点4是哑节点，它们可以用来计算裂纹张开位移，如下： 节点5至节点9也是哑节点，他们可以用来计算虚拟裂纹扩展面积?A。可以利用它们的节点坐标和高斯积分法求出面积。 到目前为止，所求的节点力和张开位移都是总体坐标下面的值，为了将他们投射到局部坐标系中，必须先构建好局部坐标系。利用节点1和节点8来计算倾斜角，即： 4.2 面状裂纹的哑节点断裂单元 * 局部坐标系中的节点力： 局部坐标系中的张开位移： 最终计算能量释放率： 注意：此处的能量释放率均是在局部坐标下计算的！ 4.2 面状裂纹的哑节点断裂单元 * SUBROUTINE UEL(RHS,AMATRX,SVARS,ENERGY,NDOFEL,NRHS,NSVARS,PROPS, 1 NPROPS,COORDS,MCRD,NNODE,U,DU,V,A,JTYPE,TIME,DTIME,KSTEP,KINC, 2 JELEM,PARAMS,NDLOAD,JDLTYP,ADLMAG,PREDEF,NPREDF,LFLAGS, 3 MLVARX,DDLMAG,MDLOAD,PNEWDT,JPROPS,NJPROP,PERIOD) C INCLUDE aba_param.inc C DIMENSION RHS(MLVARX,*),AMATRX(NDOFEL,NDOFEL),SVARS(NSVARS), 1 ENERGY(8),PROPS(*),COORDS(MCRD,NNODE), 2 U(NDOFEL),DU(MLVARX,*),V(NDOFEL),A(NDOFEL),TIME(2), 3 PARAMS(3),JDLTYP(MDLOAD,*),ADLMAG(MDLOAD,*), 4 DDLMAG(MDLOAD,*),PREDEF(2,NPREDF,NNODE),LFLAGS(*),JPROPS(*) C DIMENSION FORCE(8), XE(4), YE(4), GPX(2), GPY(2) C COMMON/A1/ XE,YE,GG,HH C DATA GPX/-0.577350269,0.577350269/ DATA GPY/-0.577350269,0.577350269/ C DO 105 K1=1,6 RHS(K1,NRHS)=0.D0 DO 105 K2=1,6 AMATRX(K2,K1)=0.D0 105 CONTINUE C STX=PROPS(1) STY=PROPS(2) STZ=PROPS(3) C AMATRX(1,1)= STX AMATRX(1,4)=-STX AMATRX(2,2)= STY AMATRX(2,5)=-STY AMATRX(3,3)= STZ AMATRX(3,6)=-STZ AMATRX(4,1)=-STX AMATRX(4,4)= STX AMAT