功能梯度材料界面开裂分析.docVIP

/FILNAME,Functionally Graded Materials,0 /PREP7 L=0.35 !模型的长度 H1=0.016 !功能梯度材料1的宽度 H2=0.016 !功能梯度材料2的宽度 A=0.03 !功能梯度材料界面裂纹长度 E0=2E11 !功能梯度材料界面处弹性模量 BEITA1=12 !功能梯度材料1的材料非均匀参数 BEITA2=-12 !功能梯度材料2的材料非均匀参数 N=7 !功能梯度材料划分的层数 ET,1,182 !定义功能梯度材料单元 KEYOPT,1,3,2 !定义单元模拟平面应变 ET,2,202 !定义与结构单元相对应的界面单元 KEYOPT,2,3,2 !定义单元模拟平面应变 *DO,I,1,N,1 MP,EX,I,E0*EXP(BEITA1*(H2/N*(2*I-1))) MP,PRXY,I,0.3 *ENDDO *DO,I,1,N,1 MP,EX,I+N,E0*EXP(BEITA2*(-H1/N*(2*I-1))) MP,PRXY,I+N,0.3 *ENDDO TB,CZM,2*N+1,,,EXPO !* 粘结区界面材料模型 TBDATA,1,95E6,0.006e-3,1 !定义关键点 K,1,0,-H1,0 K,2,L,-H1,0 K,3,L,0,0 K,4,A,0,0 K,5,0,0,0 K,6,0,0,0 K,7,0,H2,0 K,8,L,H2,0 !通过关键点定义面 A,1,2,3,4,5 A,3,4,6,7,8 wprot,0,-90,0 !转动工作平面，绕Y轴转动-90度 !通过循环语句施加网格划分 *DO,I,1,N-1,1 wpoff,0,0,H1/N ASBW,ALL *ENDDO wpoff,0,0,-H1/N*(N-1) *DO,I,1,N-1,1 wpoff,0,0,-H2/N ASBW,ALL *ENDDO LESIZE,4,,,6 !设置网格划分份数位6 LESIZE,6,,,6 ESIZE,0.005 !设置总体单元尺寸 *DO,I,1,N,1 ASEL,S,LOC,Y,H1/N*(I-1),H1/N*I MAT,I AMESH,ALL *ENDDO ALLSEL,ALL !选择所有实体 *DO,I,1,N,1 ASEL,S,LOC,Y,-H1/N*I,-H1/N*(I-1) MAT,I+N AMESH,ALL *ENDDO ALLSEL,ALL NSEL,S,LOC,X,A,L TYPE,2 MAT,2*N+1 CZMESH,,,0,Y,0, !划分界面网格 ALLSEL,ALL /SOL NLGEOM,1 !考虑大变形 NROPT,FULL, , !采用完全牛顿法 !施加关键点7和1的位移约束 DK,7, ,0.002, ,0,UY, , , , , , DK,1, ,-0.002, ,0,UY, , , , , , NSEL,S,LOC,X,L !选择在X方向等于L处的节点 D,ALL, , , , , ,ALL, , , , ,!约束所有自由度 ALLSEL,ALL OUTRES,ALL,ALL, TIME,1 NSUBST,200 KBC,0 TSRES,ERASE SOLVE