自己已经解决了,希望后来的人可以参考下。
/prep7
et,1,plane182
et,2,169
et,3,172
et,4,plane42
mp,ex,1,14.04
mp,nuxy,1,0.499
mp,mu,1,0.2
TB,HYPER,1,,2,MOONEY
tbdata,1,1.87
tbdata,2,0.47
mp,ex,4,2e5
mp,nuxy,4,0.3
R=7 ! RADIUS OF CYLINDER,(mm)
a=50 !外度宽
h0=24 !外部高
b=24 !槽宽
c=12 !槽深
h=6 !上板厚
blc4,-a/2,0,a,h0
wprota,,,-90
wpoffs,,,b/2
asbw,all
wpoffs,,,-b
asbw,all
wprota,,-90
wpoffs,,,h0-c
asbw,all
asel,s,loc,x,-b/2,b/2
asel,r,loc,y,h0-c,h0
adele,all,,,1
allsel,all
wpcsys,-1
wpoffs,b/2-r,h0-c+r
csys,4
cswpla,11,1
cyl4,0,0,r
csys,0
wpcsys,-1
wpoffs,0,h0-c+2*r
blc4,-a/2,0,a,h
allsel,all
asel,s,loc,x,b/2-2*r,b/2
asel,r,loc,y,h0-c,h0-c+2*r
csys,0
wpcsys,-1
wpoffs,b/2-r,h0-c+r
wprota,,90
asbw,all
wprota,,,90
asbw,all
cm,rubber_area,area
aatt,1,2,1
esize,1
mshkey,1
amesh,all
allsel,all
cmsel,u,rubber_area
aatt,4,,4
esize,2
amesh,all
allsel,all
!!!接触对
r,1
type,2
mat,1
real,1
lmesh,12
lmesh,15
lmesh,21
csys,11
nsel,s,loc,x,r
type,3
real,1
esurf
allsel,all
cmsel,u,rubber_area
allsel,below,area
csys,0
nsel,r,loc,y,h0-c+2*r+h
CP,1,UY,ALL
*GET,NCEN,NODE,,NUM,MIN ! GET MINIMUM NODE NUMBER FROM SELECTED SET
allsel,ALL
nsel,s,loc,y,0
d,all,all,0
allsel,all
nsel,s,loc,x,-a/2
nsel,r,loc,y,h0-c+2*r,h0-c+2*r+h
d,all,ux,0
nsel,s,loc,x,a/2
nsel,r,loc,y,h0-c+2*r,h0-c+2*r+h
d,all,ux,0
allsel,all
/solu
nlgeom,on
solc,on
OUTRES,ERASE
OUTRES,all,all
autots,on
nsubst,20,1000,10
time,1
D,NCEN,UY,-1.7
allsel,all
lswri,1
lsel,s,,,4,7,3
SFl,all,PRES,5
allsel,all
time,2
lswri,2
LSSOLVE,1,2,1,
|