刚体圆柱带动外部圆环自转问题，为何外部圆环无法转动？
建立的模型如图，内部为刚体，施加自转角速度，带动外部物体绕中心轴转动，但是无法计算，一直出现变形过大，查看部分计算结果是外部圆环被约束住了，无法转动。请各位大神老师指导一下是什么原因？
/REPLOT,RESIZE
!*
/NOPR
KEYW,PR_SET,1
KEYW,PR_STRUC,1
KEYW,PR_THERM,0
KEYW,PR_FLUID,0
KEYW,PR_MULTI,0
KEYW,LSDYNA,1
KEYW,PR_DYNA,1
/GO
!*
/COM,
/COM,Preferences for GUI filtering have been set to display:
/COM, Structural with LS-DYNA Explicit
!*
/PREP7
!*
ET,1,SOLID164
ET,2,SOLID164
!*
MP,DENS,1,7.8
MP,EX,1,2.1
MP,NUXY,1,0.3
EDMP,RIGI,2,0,0
MP,DENS,2,7.8
MP,EX,2,2.1
MP,NUXY,2,0.3
CYLIND,10,8,0,15,0,360,
CYLIND,8,0,0,15,0,360,
KWPAVE, 9
wpro,,90.000000,
FLST,2,2,6,ORDE,2
FITEM,2,1
FITEM,2,-2
VSBW,P51X
CSYS,0
WPAVE,0,0,0
CSYS,0
!*
wpro,,,90.000000
FLST,2,4,6,ORDE,2
FITEM,2,3
FITEM,2,-6
VSBW,P51X
FLST,2,8,6,ORDE,4
FITEM,2,1
FITEM,2,-2
FITEM,2,7
FITEM,2,-12
VGLUE,P51X
FLST,5,8,4,ORDE,4
FITEM,5,1
FITEM,5,-4
FITEM,5,9
FITEM,5,-12
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1, , ,10, , , , ,1
!*
FLST,5,16,4,ORDE,6
FITEM,5,21
FITEM,5,-28
FITEM,5,37
FITEM,5,-40
FITEM,5,43
FITEM,5,-46
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1, , ,8, , , , ,1
!*
FLST,5,8,4,ORDE,4
FITEM,5,31
FITEM,5,-32
FITEM,5,56
FITEM,5,-61
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1, , ,2, , , , ,1
!*
FLST,5,9,4,ORDE,9
FITEM,5,17
FITEM,5,-18
FITEM,5,29
FITEM,5,-30
FITEM,5,41
FITEM,5,-42
FITEM,5,47
FITEM,5,49
FITEM,5,53
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1, , ,10, , , , ,1
!*
FLST,2,8,6,ORDE,2
FITEM,2,1
FITEM,2,-8
VGLUE,P51X
TYPE, 1
MAT, 1
REAL,
ESYS, 0
SECNUM,
!*
MSHAPE,0,3D
MSHKEY,1
!*
FLST,5,4,6,ORDE,2
FITEM,5,3
FITEM,5,-6
CM,_Y,VOLU
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
!*
VMESH,_Y1
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
!*
TYPE, 2
MAT, 2
REAL,
ESYS, 0
SECNUM,
!*
FLST,5,4,6,ORDE,4
FITEM,5,1
FITEM,5,-2
FITEM,5,7
FITEM,5,-8
CM,_Y,VOLU
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
!*
VMESH,_Y1
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
!*
EDPART,CREATE
PARTSEL,'PLOT'
*DIM,TIME,ARRAY,2,1,1, , ,
*SET,TIME(2,1,1) , 10
*DIM,DISP,ARRAY,2,1,1, , ,
!*
*SET,DISP(2,1,1) , 10
*SET,DISP(1,1,1) , 10
EDLOAD,ADD,RBOZ,0, 2,TIME,DISP, 0, , , , ,
FINISH
/SOL
TIME,10,
EDRST,30,
EDHTIME,50,
EDDUMP,1,
EDOPT,ADD,blank,LSDYNA
EDOUT,GLSTAT
EDENERGY,1,1,1,1
全部回答 (8)
我给你幸福 2017年5月22日 独孤神宇 2017年4月26日
邹正刚 2017年4月24日试了一下你的命令流，发现以下问题：
1 材料 2 虽然定义为 rigid，但是没有给出具体约束了哪些位移；
2 两种材料的性能数据太离谱。由于弹性模量太小 (数值为 21，不知是什么单位？什么材料？)，稍一受力就会产生很大的变形，导致单元高度扭曲，造成负体积，而不收敛。
3 没有施加任何约束，结构产生刚体运动，计算时软件会出现警告信息。
最终以大量单元负体积而停止计算。
对以上问题进行修改：对材料 2 设置约束除 ROTz 外的其它自由度；两种材料的性能数据使用钢的数据：弹性模量 2.1E11；质量密度 7800。在中心线两端添加约束，使其只能有 ROTz 的转动。
终止时间设置为 100 s，计算到 90 s 之后才出现负体积，而终止计算。
一些结果如下：
其中，三个图片分别是不同时间的周向位移，可以看到 10s 时外层也是在转动的；90s 时，外层已经开始发生单元高度扭曲的现象。
txt 文件为修改后的命令流。
奈文摩尔 2017年4月21日 邹正刚 2017年4月18日 邹正刚 2017年4月18日 邹正刚 2017年4月17日
技术工 2017年4月17日
