刚体圆柱带动外部圆环自转问题,为何外部圆环无法转动?

浏览:939 回答:14

微信图片_20170415232514.png

建立的模型如图,内部为刚体,施加自转角速度,带动外部物体绕中心轴转动,但是无法计算,一直出现变形过大,查看部分计算结果是外部圆环被约束住了,无法转动。请各位大神老师指导一下是什么原因?

微信图片_20170415232709.png

微信图片_20170417175036.png



/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)
默认 最新
我给你幸福
您好,能加下楼主的qq吗 我也是做空心圆柱的仿真,不会施加约束 想请假你一下
2017年5月22日
评论 1 点赞
独孤神宇
EDMP,RIGI,2,0,0 改为 EDMP,RIGI,2,7,4 , 我也做这个刚体带动柔体的 仿真? 可以一起交流一下 QQ 523074926
2017年4月26日
评论 点赞
邹正刚

  试了一下你的命令流,发现以下问题:

  1  材料 2 虽然定义为 rigid,但是没有给出具体约束了哪些位移;

  2  两种材料的性能数据太离谱。由于弹性模量太小 (数值为 21,不知是什么单位?什么材料?),稍一受力就会产生很大的变形,导致单元高度扭曲,造成负体积,而不收敛。

  3  没有施加任何约束,结构产生刚体运动,计算时软件会出现警告信息。

  最终以大量单元负体积而停止计算。

  对以上问题进行修改:对材料 2 设置约束除 ROTz 外的其它自由度;两种材料的性能数据使用钢的数据:弹性模量 2.1E11;质量密度 7800。在中心线两端添加约束,使其只能有 ROTz 的转动。

  终止时间设置为 100 s,计算到 90 s 之后才出现负体积,而终止计算。

  一些结果如下:

out_is_fixed000.jpg

out_is_fixed002.jpg

out_is_fixed003.jpg

out_is_fixed.txt

 其中,三个图片分别是不同时间的周向位移,可以看到 10s 时外层也是在转动的;90s 时,外层已经开始发生单元高度扭曲的现象。

  txt 文件为修改后的命令流。

2017年4月24日
评论 1 点赞
奈文摩尔
既然提示了这个错误,你得看看你是否有约束外面节点自由度的操作
2017年4月21日
评论 点赞
邹正刚
1 材料 2 定义为 rigid,但是没有约束任何自由度。 2 整个模型也没有施加任何约束;只有一个转动载荷 ROTz。 3 计算时出现 160 个错误,主要是负体积问题,因此计算被终止,没有结果。
2017年4月18日
已采纳 评论 3 点赞
邹正刚
你的计算中出现了大量负体积的错误,能有计算结果吗?
2017年4月18日
评论 点赞
邹正刚
命令流中: *SET,DISP(2,1,1) , 10 *SET,DISP(1,1,1) , 10 是否有误?改为: *SET,DISP(2,1,1) , 10 *SET,DISP(1,1,1) , 0 试试。
2017年4月17日
评论 1 点赞
技术工
@dyna_focus @ 蓝牙@奈文摩尔
2017年4月17日
评论 点赞

没解决?试试专家一对一服务

换一批