本帖最后由 黑洞菩提 于 2012-7-20 11:37 编辑
ls-dyna程序如下:
修改前,辊筒为刚性体程序能完美运行为:
/PREP7
ET,1,SOLID164
UIMP,1,DENS, , ,78000000,
UIMP,1,EX, , ,1.2e11,
UIMP,1,NUXY, , ,0.32,
TB,BISO,1
TBDAT,1,1.09e8,
TBDAT,2, ,
UIMP,2,DENS, , ,79000000,
UIMP,2,EX, , ,3e11,
UIMP,2,NUXY, , ,0.29,
EDMP,RIGID,2,7,5
K, ,-0.02,-0.085,0,
K, ,0.02,-0.085,0,
K, ,-0.02,-0.080,0,
K, ,0.02,-0.080,0,
K, ,0.02,0,0,
K, ,-0.02,0,0,
GPLOT
LSTR, 1, 2
LSTR, 2, 4
LSTR, 4, 3
LSTR, 3, 1
AL,all
FLST,2,1,5,ORDE,1
FITEM,2,1 !area1
FLST,8,2,3
FITEM,8,5 !keypoint(5,6)
FITEM,8,6
VROTAT,P51X, , , , , ,P51X, ,360, ,
BLOCK,-0.01,0.01,-0.092,-0.082,0.06,0.14,
vatt,1,1,1
LESIZE,44, , ,40, , , , ,1 !the number of line75 is 40
LESIZE,39, , ,15, , , , ,1
LESIZE,38, , ,8, , , , ,1
MSHAPE,0,3D
MSHKEY,1
vmesh,5 !volume5 is meshed
vatt,2,1,1
LESIZE,17, , ,20, , , , ,1 !the number of line75 is 40
LESIZE,19, , ,15, , , , ,1
LESIZE,14, , ,3, , , , ,1
LESIZE,21, , ,20, , , , ,1
MSHAPE,0,3D
MSHKEY,1
vmesh,1,4 !meshing volume(1,2,3,4)
VSEL,S, , , 5 ! select volume5
NSLV,S,1 !the nodes is attached to volume
NPLOT
CM,slab,NODE !define component(slab)
ALLSEL,ALL
APLOT
edpart,create !creatw part ID
EDCGEN,ASTS,SLAB,2,0.35,0.18,0,0,0, , , , ,0,10000000 !define contact
ASEL,S, , , 23 !select area39
NSLA,S,1 !nodes attached to area
NPLOT
D,ALL, ,0, , , ,UY, , , , , !apply load to nodes(uy=0)
ALLSEL,ALL
APLOT
*DIM,time,ARRAY,2,1,1, , ,
*set,TIME(2,1,1),1
*DIM,rots,ARRAY,2,1,1, , ,
*set,ROTS(1,1,1),6.28
*set,ROTS(2,1,1),6.28
EDVEL,VELO,SLAB,0,0,-0.53,0,0,0,0,0,0,0,
EDLOAD,ADD,RBOX,0,2,TIME,ROTS,0,,,,,
EDENERGY,1,1,1,1
TIME,0.8,
edopt,add,blank,both
EDRST,100,
EDHTIME,1000,
把刚性体辊筒定义为柔体后有问题的程序为:
/PREP7
ET,1,SOLID164
UIMP,1,DENS, , ,10000000,
UIMP,1,EX, , ,6.1e6,
UIMP,1,NUXY, , ,0.49,
TB,BISO,1
TBDAT,1,9.24e6,
TBDAT,2,2.9e6 ,
UIMP,2,DENS, , ,79000000,
UIMP,2,EX, , ,3e11,
UIMP,2,NUXY, , ,0.29,
K, ,-0.02,-0.085,0,
K, ,0.02,-0.085,0,
K, ,-0.02,-0.080,0,
K, ,0.02,-0.080,0,
K, ,0.02,0,0,
K, ,-0.02,0,0,
GPLOT
LSTR, 1, 2
LSTR, 2, 4
LSTR, 4, 3
LSTR, 3, 1
AL,all
FLST,2,1,5,ORDE,1
FITEM,2,1 !area1
FLST,8,2,3
FITEM,8,5 !keypoint(5,6)
FITEM,8,6
VROTAT,P51X, , , , , ,P51X, ,360, ,
BLOCK,-0.01,0.01,-0.092,-0.082,0.06,0.14,
vatt,1,1,1
LESIZE,44, , ,40, , , , ,1 !the number of line75 is 40
LESIZE,39, , ,15, , , , ,1
LESIZE,38, , ,8, , , , ,1
MSHAPE,0,3D
MSHKEY,1
vmesh,5 !volume5 is meshed
vatt,2,1,1
LESIZE,17, , ,20, , , , ,1 !the number of line75 is 40
LESIZE,19, , ,15, , , , ,1
LESIZE,14, , ,3, , , , ,1
LESIZE,21, , ,20, , , , ,1
MSHAPE,0,3D
MSHKEY,1
vmesh,1,4 !meshing volume(1,2,3,4)
VSEL,S, , , 5 ! select volume5
NSLV,S,1 !the nodes is attached to volume
NPLOT
CM,slab,NODE !define component(slab)
ALLSEL,ALL
APLOT
edpart,create !creatw part ID
EDCGEN,ASTS,1,2,0.35,0.18,0,0,0, , , , ,0,10000000 !define contact
ASEL,S, , , 23 !select area39
NSLA,S,1 !nodes attached to area
NPLOT
D,ALL, ,0, , , ,UY, , , , , !apply load to nodes(uy=0)
vsel,s,mat,,2
nslv,s,1
cm,gunti,node
nsel,r,loc,x,0.02
d,all,,0,,,,ux,,,,,
allsel,below,volu
nsel,r,loc,x,-0.02
d,all,,0,,,,ux,,,,,
ALLSEL,ALL
APLOT
*DIM,time,ARRAY,2,1,1, , ,
*set,TIME(2,1,1),1
*DIM,rots,ARRAY,2,1,1, , ,
*set,ROTS(1,1,1),6.28
*set,ROTS(2,1,1),6.28
EDVEL,VELO,SLAB,0,0,-0.53,0,0,0,0,0,0,0,
EDLOAD,ADD,omgx,0,gunti,TIME,ROTS,0,,,,,
EDENERGY,1,1,1,1
TIME,0.55,
edopt,add,blank,both
EDRST,100,
EDHTIME,1000,
定义辊筒为刚性辊的时候程序能够完美的运行,当我把辊筒改为柔体把edload,add,rbox,0,2,time,rots,0,,,,,改为edload,add,omgx,0,gunti,time,rots,0,,,,,后,就运行出错了,哪位高手能帮我看一下么,可能就是一个小小的问题而已
还想知道像这样的例子能用workbench中的explicit dynamic或者transient制作吗,如果可以的话还请指点一下,我还是比较喜欢workbench的,不过是不知道这种类型的能否制作而已。
希望斑竹不要把我的帖子封了,我也是没有办法才在这个板块发ls-dyna的帖子的,请手下留情
|