接觸分析一個剛性加載的例子

mark

2006年8月28日 22:34

THICK=0.1
RADIUS=0.5
LENGTH=1.5
/PREP7
ET,1,182,1,,2
ET,2,169
ET,3,171
R,2,0,0,1,0.1
RMORE,0.1,0.01
R,3,0,0,1,0.1
RMORE,0.1,0.01
!* !* Material Properties !*
MP,EX,1,16E6
MP,NUXY,1,0.33
MP,MU,1,0.3
TB,MISO,1,0,9
TBPT,,0.000625,10000
TBPT,,0.0025,15000
TBPT,,0.005,21000
TBPT,,0.01,29000
TBPT,,0.015,32600
TBPT,,0.02,34700
TBPT,,0.04,36250
TBPT,,0.10,39000
TBPT,,0.20,40250
K,1
K,2,LENGTH
KGEN,2,1,2,1,,THICK,,2
A,1,2,4,3
TYPE,1
MAT,1
LESIZE,1,,,70
LESIZE,2,,,5
LESIZE,3,,,70
LESIZE,4,,,5
AMESH,ALL
NUMMRG,ALL
ASEL,NONE
LSEL,NONE
WPAVE,,THICK+RADIUS
CSWPLA,11,1
K,20,RADIUS,180
K,21,RADIUS,-90
K,22,RADIUS,0
K,23,0,0
k,24,radius,90
l,20,21
L,21,22
l,22,24
l,24,20
TYPE,2
REAL,2
LMESH,ALL
ALLS
CSYS,0
NSEL,S,LOC,Y,THICK
TYPE,3
ESURF
nsel,s,loc,y,0
real,3
esurf
ASEL,NONE
LSEL,NONE
WPAVE,,-THICK
CSWPLA,12,1
k,30
k,31,thick,0
k,32,thick,90
k,33,thick,180
l,31,32
l,32,33
type,2
lmesh,all
KMESH,23
KSEL,S,KP,,23
NSLK,S
*GET,N_LOAD,NODE,,NUM,MAX
csys,0
!* !* Boundary Conditions !*
asel,all
esla,s
nsle,s
nsel,r,loc,x,0
d,all,ux,0
nsel,r,loc,y,thick
d,all,uy,0
alls
FINISH
/SOLU
ANTYPE,STATIC
NLGEOM,ON
SOLC,ON
NSUBST,20,500,10
OUTRES,ALL,ALL
D,N_LOAD, ,2, , , ,ROTZ, , , , ,
SOLVE

登錄后免費查看全文

立即登錄