有没有大佬帮忙看看这个程序怎么修改,施加载荷的时候位移载荷不会出错,但是压力载荷却总是出错了?
FINISH
/CLEAR
/PREP7 !进入前处理器
ET, 1, PLANE182 !选择单元类型
ET, 2, TARGE169
ET, 3, CONTA171
KEYOPT, 1, 3, 0 !设置单元关键字(1为第一个单元,3为K3,1为所选项(Plane stress=0; Axisymmetric=1; Plane strain=3))
MP, EX, 1, 2.06E5 !输入材料弹性模量(定义材料属性)
MP, PRXY, 1, 0.28 !输入材料泊松比
TB, BKIN, 1 !输入材料应力应变特性
TBDATE, 1, 450 !输入材料屈服强度
TBDATE, 2, 0.8E5 !输入材料剪切模量
MP, MU, 1, 0.15 !输入摩擦系数(润滑摩擦条件下:钢-钢=0.15,钢-橡胶=0.30)
MP, EX, 2, 2.06E5 !输入材料弹性模量(定义材料属性)
MP, PRXY, 2, 0.28 !输入材料泊松比
TB, BKIN, 2 !输入材料应力应变特性
TBDATE, 1, 835 !输入材料屈服强度
TBDATE, 2, 0.8E5 !输入材料剪切模量
MP, MU, 2, 0.15 !输入摩擦系数(润滑摩擦条件下:钢-钢=0.15,钢-橡胶=0.30)
K, 1, 5, 83, 0 !创建关键点
K, 2, -80, 83, 0
K, 3, -80, 67.5, 0
K, 4, 5, 67.5, 0
K, 5, 0, 52, 0
K, 6, 0, 12, 0
K, 7, -72.26, 12, 0
K, 8, -72.26, 28, 0
K, 9, -6, 52, 0
K, 41, -72.26, 65.64, 0
K, 42, -72.26, 28, 0
K, 43, -6, 52, 0
K, 44, -6, 65.64, 0
K, 45, -9.5, 67.64, 0
K, 46, -9.5, 65.64, 0
K, 47, -13, 67.64, 0
K, 48, -13, 65.64, 0
K, 49, -16.5, 67.64, 0
K, 50, -16.5, 65.64, 0
K, 51, -20, 67.64, 0
K, 52, -20, 65.64, 0
K, 53, -23.5, 67.64, 0
K, 54, -23.5, 65.64, 0
K, 55, -27, 67.64, 0
K, 56, -27, 65.64, 0
K, 57, -30.5, 67.64, 0
K, 58, -30.5, 65.64, 0
K, 59, -34, 67.64, 0
K, 60, -34, 65.64, 0
K, 61, -37.5, 67.64, 0
K, 62, -37.5, 65.64, 0
K, 63, -41, 67.64, 0
K, 64, -41, 65.64, 0
K, 65, -44.5, 67.64, 0
K, 66, -44.5, 65.64, 0
K, 67, -48, 67.64, 0
K, 68, -48, 65.64, 0
K, 69, -51.5, 67.64, 0
K, 70, -51.5, 65.64, 0
K, 71, -55, 67.64, 0
K, 72, -55, 65.64, 0
K, 73, -58.5, 67.64, 0
K, 74, -58.5, 65.64, 0
K, 75, -62, 67.64, 0
K, 76, -62, 65.64, 0
K, 77, -65.5, 67.64, 0
K, 78, -65.5, 65.64, 0
K, 79, -69, 67.64, 0
K, 80, -69, 65.64, 0
* do, i, 41, 79
a = i
b = a + 1
lstr, a, b
* enddo
L, 80, 41
* do, j, 4, 39
c = j
d = j + 1
LFILLT, c, d, 0.5 !生成倒角线
* enddo
LCOMB, 4, 41
* do, k, 5, 39
e = k
f = k + 37
LCOMB, e, f !合并线段
LCOMB, 4, e
* enddo
LCOMB, 4, 40
!LSEL, S, , , 1, 40, 1
!AL, ALL !由所有线生成面-锁定滑块
L, 1, 2 !由关键点生成线段-42
L, 2, 3 !43
L, 3, 4 !44
L, 4, 1 !45
L, 5, 6 !46
L, 6, 7 !47
L, 7, 8 !48
L, 8, 9 !49
L, 9, 5 !50
AL, 1, 2, 3, 4 !由线生成面
AL, 5, 6, 7, 8
AL, 9, 10, 11, 12, 13
/PNUM, KP, 1 !显示关键点编号
/PNUM, LINE, 1 !显示线段编号
/PNUM, AREA, 1 !显示面编号
APLOT !显示面
ALLSEL !选择所有实体
MAT, 2 !锁定滑块
REAL, 3 !实常数
TYPE, 1 !单元类型
MSHKEY, 0 !自由网格化分
SMRTSIZE, 2 !智能网格划分控制
AESIZE, 2 !对面上线段设置网格划分控制
AMESH, 1 !面上生成面单元 - 对面进行网格划分
MAT, 1 !管壁
REAL, 3
TYPE, 1
MSHKEY, 0
SMRTSIZE, 1
AESIZE, 1
AMESH, 2
MAT, 2 !承压头
REAL, 3
TYPE, 1
MSHKEY, 0
SMRTSIZE, 1
AESIZE, 1
AMESH, 3
EPLOT !显示单元
ALLSEL !选择所有实体
!开始定义接触单元
/COM, CONTACT PAIR CREATION - START
CM,_NODECM,NODE
CM,_ELEMCM,ELEM
CM,_KPCM,KP
CM,_LINECM,LINE
CM,_AREACM,AREA
CM,_VOLUCM,VOLU
/GSAV,cwz,gsav,,temp
MP, MU, 1, 0.15
R, 5
REAL, 5
ET, 6, 169
ET, 7, 172
!KEYOPT, 7, 9, 0
!KEYOPT, 7, 10, 2
KEYOPT, 7, 12, 4
R, 5, , , 0.1, 0.1
! Generate the target surface
LSEL,S,,, 2
CM,_TARGET,LINE
TYPE,6
MAT, 2
NSLL,S,1
ESLN,S,0
ESURF
CMSEL,S,_ELEMCM
! Generate the contact surface
LSEL,S,,, 12
CM,_CONTACT,LINE
TYPE,7
MAT, 2
NSLL,S,1
ESLN,S,0
ESURF
ALLSEL
ESEL,ALL
ESEL,S,TYPE,,6
ESEL,A,TYPE,,7
ESEL,R,REAL,,5
/PSYMB,ESYS,1
/PNUM,TYPE,1
/NUM,1
EPLOT
ESEL,ALL
ESEL,S,TYPE,,6
ESEL,A,TYPE,,7
ESEL,R,REAL,,5
CMSEL,A,_NODECM
CMDEL,_NODECM
CMSEL,A,_ELEMCM
CMDEL,_ELEMCM
CMSEL,S,_KPCM
CMDEL,_KPCM
CMSEL,S,_LINECM
CMDEL,_LINECM
CMSEL,S,_AREACM
CMDEL,_AREACM
CMSEL,S,_VOLUCM
CMDEL,_VOLUCM
/GRES,cwz,gsav
CMDEL,_TARGET
CMDEL,_CONTACT
/COM, CONTACT PAIR CREATION - END
/COM, CONTACT PAIR CREATION - START
CM,_NODECM,NODE
CM,_ELEMCM,ELEM
CM,_KPCM,KP
CM,_LINECM,LINE
CM,_AREACM,AREA
CM,_VOLUCM,VOLU
/GSAV,cwz,gsav,,temp
MP,MU, 1, 0.15
R, 6
REAL, 6
ET, 8, 169
ET, 9, 172
!KEYOPT, 9, 9, 0
!KEYOPT, 9, 10, 2
R,6, , , 0.1, 0.1
RMORE,
RMORE, , 0
RMORE, 0
! Generate the target surface
LSEL,S,,,4
CM,_TARGET,LINE
TYPE,8
MAT, 1
NSLL,S,1
ESLN,S,0
ESURF
CMSEL,S,_ELEMCM
! Generate the contact surface
LSEL,S,,, 7
CM,_CONTACT, LINE
TYPE,9
MAT, 2
NSLL,S,1
ESLN,S,0
ESURF
ALLSEL
ESEL, ALL
ESEL, S, TYPE,, 8
ESEL, A, TYPE,, 9
ESEL, R, REAL,, 6
/PSYMB, ESYS, 1
/PNUM, TYPE, 1
/NUM,1
EPLOT
ESEL,ALL
ESEL,S,TYPE,,8
ESEL,A,TYPE,,9
ESEL,R,REAL,,6
CMSEL,A,_NODECM
CMDEL,_NODECM
CMSEL,A,_ELEMCM
CMDEL,_ELEMCM
CMSEL,S,_KPCM
CMDEL,_KPCM
CMSEL,S,_LINECM
CMDEL,_LINECM
CMSEL,S,_AREACM
CMDEL,_AREACM
CMSEL,S,_VOLUCM
CMDEL,_VOLUCM
/GRES,cwz,gsav
CMDEL,_TARGET
CMDEL,_CONTACT
/COM, CONTACT PAIR CREATION - END
/SOLU !进入求解器
ANTYPE, STATIC !选择求解类型
NLGEOM, ON
AUTOT, ON
TIME, 1 !指定计算终止时间
NSUBST, 25000, 100000, 2500 !指定迭代子步数
LNSRCH, ON !激活线性搜索
PRED, 1 !激活自由度求解预测
EPLOT !显示单元
ALLSEL !选择所有实体
!LSEL, S, , , 5 !选择线段(1,9中间间隔8)
!LSEL, A, , , 10 !添加到上面线段选择集
!NSLL, S, 1 !选择线段上所有节点
!D, ALL, UY !施加位移约束
DL,5,, ALL
DL,10,, ALL
!LSEL, S, , , 8
!LSEL, A, , , 9
!NSLL, S, 1
!D, ALL, UX
DL,8,, ALL
DL,9,, ALL
DL,1,, UZ
EPLOT !显示单元
ALLSEL !选择所有实体
!LSEL, S, , , 1
!NSLL, S, 1
!D, ALL, UX, 1
SFL, 1, PRES, 22.4
!LSEL, S, , , 1
!NSLL, S, 1
!SF, ALL, PRES, 22.4
ALLSEL
SOLVE !开始求解计算
FINISH