您能给我看一下命令吗,我想得到一个贯通的塑性区,初学阶段才疏学浅,始终办不到,拜托了
finish
/clear
/PREP7
! KeyPoints
K,2,0.03,0
K,259,3.70E-02,-1.9E-02
K,260,3.68E-02,-1.73736000000001E-02
K,261,3.66999999998281E-02,-1.64644499999999E-02
K,262,3.66999999998281E-02,-1.56550500000001E-02
K,263,3.68E-02,-.0147573
K,264,.0367,-1.17477000000001E-02
K,265,.0366,-1.05079499999999E-02
K,266,.0367,-8.38755000000015E-03
K,267,3.7E-02,-7.24755000000005E-03
K,268,3.73E-02,-.0056259
K,269,3.79E-02,-2.73029999999994E-03
K,270,3.82E-02,-8.32200000000057E-04
K,272,.1,0
K,273,.1,-1.9E-02
K,274,0.03,-1.9E-02
K,282,3.82E-02,0
K,159,3.70E-02,-1.9E-02
K,160,3.68E-02,-1.73736000000001E-02
K,161,3.66999999998281E-02,-1.64644499999999E-02
K,162,3.66999999998281E-02,-1.56550500000001E-02
K,163,3.68E-02,-.0147573
K,164,.0367,-1.17477000000001E-02
K,165,.0366,-1.05079499999999E-02
K,166,.0367,-8.38755000000015E-03
K,167,3.7E-02,-7.24755000000005E-03
K,168,3.73E-02,-.0056259
K,169,3.79E-02,-2.73029999999994E-03
K,170,3.82E-02,-8.32200000000057E-04
K,182,3.82E-02,0
LSTR, 260, 259
LSTR, 261, 260
LSTR, 262, 261
LSTR, 263, 262
LSTR, 264, 263
LSTR, 265, 264
LSTR, 266, 265
LSTR, 267, 266
LSTR, 268, 267
LSTR, 269, 268
LSTR, 270, 269
LSTR, 282, 270
LSTR, 282,2
LSTR, 274,2
LSTR, 273,272
LSTR, 274,259
LSTR, 160, 159
LSTR, 161, 160
LSTR, 162, 161
LSTR, 163, 162
LSTR, 164, 163
LSTR, 165, 164
LSTR, 166, 165
LSTR, 167, 166
LSTR, 168, 167
LSTR, 169, 168
LSTR, 170, 169
LSTR, 182, 170
LSTR, 182, 272
LSTR, 273,159
allsel,all,all
LSEL,S,LINE,,1,12
LCOMB,ALL
ALLSEL,ALL
LSEL,S,LINE,,17,28
LCOMB,ALL
ALLSEL,ALL
numcmp,all
ALLSEL,ALL
k,9,3.70E-02,-0.6000000
k,10,3.82E-02,0.300000
k,11,0.6,-0.6000000
k,12,0.6,0.300000
LSTR, 9,2
LSTR, 10,3
LSTR, 10,12
LSTR, 12,11
LSTR, 9,11
allsel,all,all
a,2,3,5,6
a,7,1,8,4
lsel,S,line,,7,13
lsel,a,line,,4
cm,YUTI,line
al,YUTI
!土体参数
E_1=10E6
U_1=0.35
DENS_1=1900
C_1=30E3
FI_1=12
DFI_1=6
MP,EX,1,E_1
MP,PRXY,1,U_1
MP,DENS,1,DENS_1
TB,DP,1
TBDATA,1,C_1,FI_1,DFI_1
MP,MU,1,0.2
!锚固体参数
E_2=30E10
U_2=0.2
DENS_2=2500
MP,EX,2,E_2
MP,PRXY,2,U_2
MP,DENS,2,DENS_2
ET,1,82
KEYOPT,1,3,1
ET,3,targe169
ET,4,conta172
KEYOPT,4,3,1
KEYOPT,4,10,2
KEYOPT,4,8,2
keyopt,4,12,2
!界面层接触面实常数
R,6
lsel,s,line,,1,6,5
LESIZE,ALL,,,100
lsel,s,line,,2,5
LESIZE,ALL,,,15
lsel,s,line,,7
LESIZE,ALL,,,60,2
lsel,s,line,,8
LESIZE,ALL,,,60,0.5
lsel,s,line,,9
LESIZE,ALL,,,100,0.05
lsel,s,line,,10
LESIZE,ALL,,,80,0.05
lsel,s,line,,11,13
LESIZE,ALL,,,50
ASEL,S,AREA,,1
AATT,1,,1
mshkey,0
mshape,1,2D
AMESH,ALL
ASEL,S,AREA,,2
AATT,2,,1
mshkey,0
mshape,1,2D
AMESH,ALL
ASEL,S,AREA,,3
AATT,1,,1
mshkey,0
mshape,1,2D
AMESH,ALL
allsel,all,all
!创建目标面(锚固体)
lsel,s,line,,1
nsll,s,1
esln,s,0
type,3
real,6
mat,2
esurf
allsel,all,all
!创建接触面(界面层)
lsel,s,line,,6
nsll,s,1
esln,s,0
type,4
real,6
mat,1
esurf
allsel,all,all
LSEL,S,LINE,,11,13,2
NSLL,S
D,ALL,Uy
LSEL,S,LINE,,12
NSLL,S
D,ALL,Ux
LSEL,s,,,2
NSLL,S
SF,all,PRES,-4e4
LSEL,s,,,3
NSLL,S
SF,all,PRES,1e4
allsel,all,all
/solu
antype,static
nlgeom,on
autots,on
time,1
NSUBST,200
OUTRES,ALL,ALL
NROPT,full,,off
NROPT,unsym,,off
solve
allsel,all,all
ASEL,S,AREA,,1,2
ESLA,S |