有沒(méi)有大佬幫忙看看這個(gè)程序怎么修改，施加載荷的時(shí)候位移載荷不會(huì)出錯(cuò)，但是壓力載荷卻總是出錯(cuò)了？
FINISH
/CLEAR
/PREP7 !進(jìn)入前處理器
ET, 1, PLANE182 !選擇單元類型
ET, 2, TARGE169
ET, 3, CONTA171
KEYOPT, 1, 3, 0 !設(shè)置單元關(guān)鍵字（1為第一個(gè)單元，3為K3，1為所選項(xiàng)（Plane stress=0; Axisymmetric=1; Plane strain=3））
MP, EX, 1, 2.06E5 !輸入材料彈性模量（定義材料屬性）
MP, PRXY, 1, 0.28 !輸入材料泊松比
TB, BKIN, 1 !輸入材料應(yīng)力應(yīng)變特性
TBDATE, 1, 450 !輸入材料屈服強(qiáng)度
TBDATE, 2, 0.8E5 !輸入材料剪切模量
MP, MU, 1, 0.15 !輸入摩擦系數(shù)（潤(rùn)滑摩擦條件下：鋼-鋼=0.15，鋼-橡膠=0.30）
MP, EX, 2, 2.06E5 !輸入材料彈性模量（定義材料屬性）
MP, PRXY, 2, 0.28 !輸入材料泊松比
TB, BKIN, 2 !輸入材料應(yīng)力應(yīng)變特性
TBDATE, 1, 835 !輸入材料屈服強(qiáng)度
TBDATE, 2, 0.8E5 !輸入材料剪切模量
MP, MU, 2, 0.15 !輸入摩擦系數(shù)（潤(rùn)滑摩擦條件下：鋼-鋼=0.15，鋼-橡膠=0.30）
K, 1, 5, 83, 0 !創(chuàng)建關(guān)鍵點(diǎn)
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 !由關(guān)鍵點(diǎn)生成線段-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 !顯示關(guān)鍵點(diǎn)編號(hào)
/PNUM, LINE, 1 !顯示線段編號(hào)
/PNUM, AREA, 1 !顯示面編號(hào)
APLOT !顯示面
ALLSEL !選擇所有實(shí)體
MAT, 2 !鎖定滑塊
REAL, 3 !實(shí)常數(shù)
TYPE, 1 !單元類型
MSHKEY, 0 !自由網(wǎng)格化分
SMRTSIZE, 2 !智能網(wǎng)格劃分控制
AESIZE, 2 !對(duì)面上線段設(shè)置網(wǎng)格劃分控制
AMESH, 1 !面上生成面單元 - 對(duì)面進(jìn)行網(wǎng)格劃分
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 !選擇所有實(shí)體
!開(kāi)始定義接觸單元
/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 !進(jìn)入求解器
ANTYPE, STATIC !選擇求解類型
NLGEOM, ON
AUTOT, ON
TIME, 1 !指定計(jì)算終止時(shí)間
NSUBST, 25000, 100000, 2500 !指定迭代子步數(shù)
LNSRCH, ON !激活線性搜索
PRED, 1 !激活自由度求解預(yù)測(cè)
EPLOT !顯示單元
ALLSEL !選擇所有實(shí)體
!LSEL, S, , , 5 !選擇線段（1，9中間間隔8）
!LSEL, A, , , 10 !添加到上面線段選擇集
!NSLL, S, 1 !選擇線段上所有節(jié)點(diǎn)
!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 !選擇所有實(shí)體
!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 !開(kāi)始求解計(jì)算
FINISH
