ANSYS计算受压闭合裂纹的应力强度因子实例APDL命令流（原创）.pdfVIP

! ! 师访，中国矿业大学力学与建筑工程学院，江苏徐州 ! SHI Fang, China University of Mining Technology ! Website: , QQ: 1549221758 ! Email: fshi@ /126.com ! Date: March 28, 2014 ! 倾斜闭合裂纹应力强度因子计算实例 ! /thread-1102153-1-1.html ! 这里给出一个利用ANSYS 计算压剪裂纹应力强度因子的APDL 命令流，全参数 化控制。 注意事项： （1）定义裂纹面接触关系； （2 ）利用KSCON 定义裂尖奇异点:KSCON,关键点号,半径,1,10,0.5； （3 ）在裂尖建立局部坐标系，命令是:LOCAL,11,0,x,y,z,转角,转角,转角,1,1，建 立的局部坐标系的x 轴必须与裂纹面平行； （4 ）建立路径，命令：PATH 和PPATH ，路径上布置5 个节点； （5 ）在后处理/POST1 中激活局部坐标系，命令RSYS,11； （6 ）KCALC 命令计算应力强度因子。 y 方向应力： ANSYS 输出的应力强度因子： 命令流： finish /clear /Titie,SHI Fang,CUMT /PREP7 /DSCALE,ALL,1.0 /eshape,on !控制参数 w = 2.05 h = 4.05 a = w/5*sqrt(2) !裂纹半长 angle = 45 !裂纹与 σ1方向的夹角（度） f = 100e6 !σ1 COF = 0.3 !摩擦系数 ! theta = angle*3.1415926/180 !裂纹与 σ1方向的夹角（弧度） /DSCALE,ALL,1 /PNUM,KP,1 /PNUM,LINE,1 ! et,1,183 KEYOPT,1,3,2 ! mp,ex,1,65e9 mp,prxy,1,0.25 MP,MU,1,COF ! k,1,0,0 k,2,w,0 k,3,w,h k,4,0,h k,5,0,h/2-w/2/tan(theta) k,6,w,h/2+w/2/tan(theta) k,7,w/2+a*sin(theta),h/2+a*cos(theta) k,8,w/2-a*sin(theta),h/2-a*cos(theta) k,9,w/2,h/2 k,19,w/2,h/2 k,21,w/2,h k,22,w/2,0 !111 和222 关键点用于定义路径，以便计算应力强度因子 k,111,w/2+a/2*sin(theta),h/2+a/2*cos(theta) k,222,w/2+a/2*sin(theta),h/2+a/2*cos(theta) ! a,1,22,19,8,5,1 a,22,2,6,7,222,19,22 a,9,111,7,6,3,21,9 a,5,8,9,21,4,5 ! dk,1,all,0 DL,1, ,uy, DL,6, ,uy, ! allsel,all SFL,14,PRES,f, SFL,17,PRES,f, ! lesize,all,w/60 KSCON,7,a/10,1,10,0.5 !定义裂纹尖端 KSCON,8,a/10,1,10,0.5 !定义裂纹尖端 MSHKEY,0 AMESH,all !对面进行网格划分 ! allsel,all LSEL,S,,,3,9,6 LSEL,a,,,10,, NSLL,S,1 CM,lowPair,NODE allsel,all LSEL,S,,,11