其中A和TOLER的数值在开头已经定义了分别是25和0.1E-5。语句已经加粗，拜托了，各位大神。
/PREP7
DIS1=0.1
DIS2=14
L=150
H=4
A=25
TOLER=0.1E-5

ET,1,182 !2维4节点结构实体单元
KEYOPT,1,1,2 ! 增强的应变公式
KEYOPT,1,3,2 !设置单元行为平面应变
ET,2,182
KEYOPT,2,1,2
KEYOPT,2,3,2

ET,3,202 !2维4节点粘结区单元
KEYOPT,3,3,2 !设置单元为平面应变

MP,EX,1,1.353E5 !EX = 135.3 GPA
MP,PRXY,1,0.25

G1C=0.58 !临界能量释放率
G2C=0.82
G3C=0.82
TB,CGCR,1,,3,LINEAR !线性断裂准则
TBDATA,1,G1C,G2C,G3C

! 定义有限元模型
RECTNG,0,L,H/2 !定义面
RECTNG,0,L,0,-H/2
LSEL,S,LINE,,2,8,2 !设置线的网格划分数
LESIZE,ALL,H/8
LSEL,INVE
LESIZE,ALL, , ,300
ALLSEL,ALL
TYPE,1 !对面2划分网格
MAT,1
LOCAL,11,0,0,0,0
ESYS,11
AMESH,2
CSYS,0
TYPE,2 ! 对面1划分网格
ESYS,11
AMESH,1
CSYS,0
NSEL,S,LOC,X,A-TOLER,L
NUMMRG,NODES
ESLN
TYPE,3
MAT,5
CZMESH,,,1,Y,0, !划分界面单元
ALLSEL,ALL
NSEL,S,LOC,X,L !施加约束条件
D,ALL,ALL
NSEL,ALL
!
ESEL,S,ENAME,,202 !* 选择裂纹扩展路径上的界面单元
CM,CPATH,ELEM !* 定义单元组件

NSLE
NSEL,S,LOC,X,A
NSEL,R,LOC,Y,0
CM,CRACK1,NODE !* 定义裂纹尖端节点组件
ALLSEL,ALL
FINISH
/SOLU
NSEL,S,LOC,X
NSEL,R,LOC,Y,H/2 !* 在顶端施加位移载荷
D,ALL,UY,DIS1
ALLSEL,ALL
NSEL,S,LOC,X
NSEL,R,LOC,Y,-H/2 !* 在底端施加位移载荷
D,ALL,UY,-DIS1
ALLSEL,ALL
AUTOTS,ON
TIME,1

CINT,NEW,1 !* 定义裂纹尖端号
CINT,TYPE,VCCT !* 声明进行能量释放率计算
CINT,CTNC,CRACK1 !* 裂纹尖端节点组件
CINT,NORM,0,2

! 裂纹扩展模拟设置
CGROW,NEW,1 !* 声明开始CRACK GROWTH SET
CGROW,CID,1 !* 为VCCT设置ID号
CGROW,CPATH,CPATH !* 裂纹扩展路径
CGROW,FCOP,MTAB,1 !* 定义断裂准则
CGROW,DTIME,2E-3
CGROW,DTMIN,2E-3
CGROW,DTMAX,2E-3
CGROW,FCRA,0.98
NSUB,4,4,4
ALLSEL,ALL
OUTRES,ALL,ALL
SOLVE
TIME,2
NSEL,S,LOC,X
NSEL,R,LOC,Y,H/2 !*在顶端施加位移载荷
D,ALL,UY,DIS2
ALLSEL,ALL
NSEL,S,LOC,X
NSEL,R,LOC,Y,-H/2 !*在底端施加位移载荷
D,ALL,UY,-DIS2
ALLSEL,ALL
NSUBST,1000,1000,10
OUTRES,ALL,ALL
SOLVE !*开始求解
FINISH
/POST1
SET,LAST
PRCI,1
FINISH@wuming524 返回小木虫查看更多