本帖最后由 yinhb77 于 2015-4-8 22:58 编辑
/prep7
et,1,45
ex,1,200.e9
nuxy,1,0.3
dens,1,8700.
mu,1,0.2
ex,2,300e9
nuxy,2,0.3
dens,2,7850
!!===========build and mesh the pipe==========
cyl4,,,0.124,0,0.136,180,0.65
cyl4,,,0.124,180,0.136,360,0.65
vglue,all
numcmp,all
lsel,s,,,9,12,1
lesize,all,,,24 !along pipe length
lsel,s,,,2,4,2
lsel,a,,,13,14,1
lesize,all,,,18 ! along pipe circumference
lsel,s,,,1,3,2
lesize,all,,,1 ! along pipe thickness
mshape,0,3d
mshkey,1
type,1
mat,1
r,1
real,1
vmesh,all
!!=============build and mesh all the wheels===
cyl4,,,0.15,0,0.178,180,0.05
cyl4,,,0.15,180,0.178,360,0.05
vglue,3,4
numcmp,all
lsel,s,,,25,28,1
lesize,all,,,2 ! along length
lsel,s,,,18,20,2
lsel,a,,,29,30,1
lesize,all,,,14
lsel,s,,,22,24,2
lesize,all,,,1
mshape,0,3d
mshkey,1
type,1
mat,1
r,2
real,2
vmesh,3,4
vgen,1,3,4,,-0.2,-0.2421,0.05,,,1
vgen,2,3,4,,+0.4,-0.0000,0.0,,,0
vgen,2,3,6,,+0.0,-0.0000,0.45,,,0
!!!!!===========define some node for pipe and wheels
alls
*get,nmaxtotal,node,,num,max
nd_rb=nmaxtotal+1
n,nd_rb,-0.2,-0.2421,-0.1
nd_nrb=nmaxtotal+2
n,nd_nrb,-0.2,-0.2421,0.075
nd_frb=nmaxtotal+3
n,nd_frb,-0.2,-0.2421,0.525
nd_lb=nmaxtotal+4
n,nd_lb,+0.2,-0.2421,-0.1
nd_nlb=nmaxtotal+5
n,nd_nlb,+0.2,-0.2421,0.075
nd_flb=nmaxtotal+6
n,nd_flb,+0.2,-0.2421,0.525
!** Shaft
et,2,188
keyopt,2,3,3
sectype,1,BEAM,CSOLID,Shaft, 0
secoffset,CENT
secdata,0.05
type,2
mat,2
r,3
real,3
secnum,1
e,nd_rb,nd_nrb
e,nd_nrb,nd_frb
et,3,188
keyopt,3,3,3
sectype,2,BEAM,CSOLID,Shaft, 0
secoffset,CENT
secdata,0.05
type,3
mat,2
r,4
real,4
secnum,2
e,nd_lb,nd_nlb
e,nd_nlb,nd_flb
!====================
!** bearings radial stiffness modeling
et,4,184
keyopt,4,1,1
r,5
wpcsys, -1 , 0
NWPAVE,nd_rb
CSWPLA,11,1
CSYS,11
alls
nsel,s,loc,x,0.148,0.152 ! for near right bottom wheel
nsel,r,loc,z,0.173,0.177
*get,nmaxnrb,node,,num,max
*get,nminnrb,node,,num,min
type,4
real,5
mat,2
*do,k,nminnrb,nmaxnrb
*if,nsel(k),eq,1,then
e,nd_nrb,k
*endif
*enddo
alls
nsel,s,loc,x,0.148,0.152 ! for far right bottom wheel
nsel,r,loc,z,0.623,0.627
*get,nmaxfrb,node,,num,max
*get,nminfrb,node,,num,min
type,4
real,5
mat,2
*do,k,nminfrb,nmaxfrb
*if,nsel(k),eq,1,then
e,nd_frb,k
*endif
*enddo
wpcsys, -1 , 0
NWPAVE,nd_lb
CSWPLA,12,1
CSYS,12
alls
nsel,s,loc,x,0.148,0.152 ! for near left bottom wheel
nsel,r,loc,z,0.173,0.177
*get,nmaxnlb,node,,num,max
*get,nminnlb,node,,num,min
type,4
real,5
mat,2
*do,k,nminnlb,nmaxnlb
*if,nsel(k),eq,1,then
e,nd_nlb,k
*endif
*enddo
alls
nsel,s,loc,x,0.148,0.152 ! for far left bottom wheel
nsel,r,loc,z,0.623,0.627
*get,nmaxflb,node,,num,max
*get,nminflb,node,,num,min
type,4
real,5
mat,2
*do,k,nminflb,nmaxflb
*if,nsel(k),eq,1,then
e,nd_flb,k
*endif
*enddo
!==========================================
et,5,TARGE170 !STIFFER-LOWER ORDER-COARSER-FLAT/CONCAVE-LARGER SURFACE
et,6,CONTA174
allsel,all
wpcsys, -1 , 0
csys,1
/com,contact pair creation-start
r,5
keyopt,6,5,1
keyopt,6,9,0
keyopt,6,10,1
keyopt,6,12,1
asel,s,,,13,43,10
asel,a,,,19,49,10
cm,area_target,area ! close to origin
mat,1
type,5
real,5
nsla,s,1
esln,s
esurf
esel,all
asel,s,,,3,9,6
cm,area_contact,area
mat,1
type,6
real,5
nsla,s,1
esln,s
esurf
allsel
alls
vsel,s,,,1,2,1
eslv,s
cm,epipe,elem
!**components
alls
esel,s,type, ,2 ! for right bottom shaft
cm,beam_r,elem
esel,s,type, ,3 ! for left bottom shaft
cm,beam_l,elem
!==========================================
alls
save
fini
/sol
antype,trans
/nerr,,1000000
csys,0
!acel,,-9.8
cmacel,epipe,,-9.8
pi=acos(-1)
vr_rad=2.0*pi/6
cmomega,beam_r,,,vr_rad,-0.2,-0.2421,-0.1,-0.2,-0.2421,0.525
cmomega,beam_l,,,vr_rad,+0.2,-0.2421,-0.1,+0.2,-0.2421,0.525
d,nd_nrb,UX,0
d,nd_nrb,UY,0
d,nd_nrb,UZ,0
d,nd_frb,UX,0
d,nd_frb,UY,0
d,nd_frb,UZ,0
d,nd_nlb,UX,0
d,nd_nlb,UY,0
d,nd_nlb,UZ,0
d,nd_flb,UX,0
d,nd_flb,UY,0
d,nd_flb,UZ,0
alls
NLGEOM,ON
autots,on
time,1
deltim,0.05,0.005,0.5
outres,basic,last
allsel,all
solve
|
IP卡
狗仔卡
发表于 2015-4-8 22:40:16
提升卡
沉默卡
喧嚣卡
变色卡
抢沙发
千斤顶
显身卡
