轮子转动和摩擦问题求助，谢谢

yinhb77

求助，希望大家跑跑命令流，给指点一下

希望转动两根轴，带动四个轮子，轮子通过摩擦使放在上面的管转动。
可是结果不转动，并且管子飞走了（已经加重力）

谢谢

/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 andmesh 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   

!!=============buildand 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

!!!!!===========definesome 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 radialstiffness 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 paircreation-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

yinhb77

anybody help?
Plz!