```
!----------------------生成接触对
MP,MU,1,0.15 !定义摩擦系数
MAT,2 !激活材料属性2
R,3 !定义实常数3
REAL,3 !激活实常数3
ET,2,170 !定义单元类型2 contact单元
ET,3,174 !定义单元类型3
KEYOPT,3,9,0 !设定单元类型4的关键项9
KEYOPT,3,10,2 !设定单元类型3的关键项10
R,3,,,0.1,0.1,0, !设定法向接触刚度为0.1
RMORE,,,1.0E20,0.0,1.0,
RMORE,0.0,0,1.0,,1.0,0.5
RMORE,0,1.0,1.0,0.0,,1.0
NSEL,S,,,OUT_NODE !生成目标面
TYPE,2 !激活单元类型2
ESLN,S,0 !选择节点依附的单元
ESURF !在当前选择的单元上覆盖生成单元
ESEL,ALL !选择所有的单元
NSEL,S,,,in_NODE !生成接触面
TYPE,3 !激活单元类型3
ESLN,S,0 !选择节点依附的单元
ESURF !在当前选择的单元上覆盖生成单元
ESEL,ALL !选择所有的单元
ALLSEL !选择所有的图元
WPCSYS,-1,0
WPSTYLE,,,,,,,,0
FINISH
!********************边界条件
/SOLU
ASEL,S,LOC,Y,0
DA,ALL,SYMM
ALLSEL
NSEL,S,,,CM_sym
CM,_sym,NODE
dsym,symm,Z !在节点上施加对称边界条件
NSEL,S,,,CM_Fix
CM,_Fix,NODE
d,all,all !在节点上施加全部自由度约束
NSEL,S,,,CM_Fixx
CM,_Fixx,NODE
d,all,UX
NSEL,S,,,CM_FixY
CM,_FixY,NODE
d,all,UX
*set,F,1479000
SFL,13,PRES,F/(d0*b)
!************************************************!非线性求解设置
lnsrch,on !打开线性搜索
pred,on !打开预测矫正
autot,on !打开自动时间步
nsubst,8,25,2 !设定子步数
ALLSEL !选择全部图元```这是apdl代码,请加注释,并指出问题