lsdyna中如何定义弹性体绕质心的转动实例.docx
《lsdyna中如何定义弹性体绕质心的转动实例.docx》由会员分享，可在线阅读，更多相关《lsdyna中如何定义弹性体绕质心的转动实例.docx（12页珍藏版）》请在冰豆网上搜索。
lsdyna中如何定义弹性体绕质心的转动实例
ls-dyna中如何定义弹性体绕质心的转动实例
关于在ls－dyna中如何定义弹性体绕质心的转动实例
在我们学习ls－dyna的过程中，可能会遇到很多情况下不适宜把物体设置为刚体，而设置为弹塑性体是比较符合实际情况的，而在ls－dyna中关于转动的设置对于刚性体来说是比较容易的，而对于弹性体来说就不是很容易了。
在这里，我通过一个实例来说明一下我的一些体会。
假如一个空心圆柱实体，它只有绕轴心线转动的自由度，而别的方向的自由度要固定：
如在轴向不能移动，则可以通过设置边在轴向固定就是了；同时要使圆柱体在径向也没有移动，就是圆柱体除了绕轴线转动外，不能再有其他方向的运动，这是我们就不是很好定义了。
因为在ls－dyna中只有笛卡儿坐标系，其他的形式都是靠它的旋转来得到，没有向ansys中一样的圆柱坐标可以来控制径向的约束。
在这种情况下，我参考了一些资料，就认为可以把变形体与刚体结合起来应用，既是刚体与变形体通过公共节点结合起来，然后设置刚体的自由度来约束变形体的运动。
如我举例中的空心圆柱体，我们就可以在圆柱体的圆面加一个圆盘刚体，让它与圆柱体同心，这样我们就是把刚体的自由度定义为只有绕轴线的转动，而没有别的自由度了。
这样既是受到力，变形体也不会在径向移动了。
要是你感觉这样对变形体的变形不均匀，则可是在两边分别一个刚体盘，保持对称性就更理想了。
但是不容置疑的刚体与变形体连接的地方肯定有变形及相互作用，可能还很大，要是这样对分析有影响，确实不好办，不知道大家还有没有更好的办法来做这个。
请一起相互讨论，共同进步吧
下面是实例（关于轧辊变形体与一个同心刚体联合起来控制轧辊的运动约束的，使轧辊只能绕轴线旋转），在ansys中的命令流，大家可以参考下,希望对部分学习者有用。
!
设置标题
/TITLE，PROROLLINGANSYS
!
进入前处理器
/PREP7
!
定义单元类型
ET,1,SOLID164
!
定义轧件材料模型
MP,DENS,1,785000
MP,EX,1,117e9
MP,NUXY,1,0.362
TB,BISO,1,,,,
TBDAT,1,134e6
TBDAT,2,13.4e6
!
定义工作辊辊身材料模型
MP,DENS,2,785000
MP,EX,2,210e9
MP,NUXY,2,0.3
!
定义支撑辊辊身材料模型
*CSET,1,2,2,3,
MPCOPY,,2,3
TBCOPY,ALL,2,3
!
定义工作辊辊径材料模型
EDMP,RIGI,4,7,4
MP,DENS,4,785000
MP,EX,4,210e9
MP,NUXY,4,0.3
!
定义支撑辊辊径材料模型
*CSET,1,2,4,5,
MPCOPY,,4,5
TBCOPY,ALL,4,5
!
建立轧件材料模型
BLOCK,-0.10,-0.60,,0.019,,0.6,
!
创建关键点
K,,,0.39376,-0.001,
K,,,0.39376,0.79,
K,,,0.77376,0.79,
K,,,0.77376,-0.001,
K,,,1.53376,-0.001,
K,,,1.53376,0.79,
K,,,2.29376,0.79,
K,,,2.29376,-0
.001,
K,,,0.39376,0.89,
K,,,0.62376,0.89,
K,,,0.62376,0.79,
K,,,1.53376,0.89,
K,,,1.99376,0.89,
K,,,1.99376,0.79,
!
创建面
A,16,15,14,13
A,22,21,20,14
A,12,11,10,9
A,19,18,17,10
!
面旋转生成体
VROTAT,7,8,,,,,13,20,360,,
VROTAT,9,10,,,,,9,17,360,,
!
选择体（工作辊）
VSEL,S,,,10,17
!
在选择的体下体编辑
ALLSEL,BELOW,VOLU
!
对体进行GLUE
VGLUE,10,11,12,13,14,15,16,17
!
对面进行ADD
AADD,55,74
AADD,63,76
AADD,70,78
AADD,47,71
!
对线进行ADD
LCOMB,70,98,0
LCOMB,80,99,0
LCOMB,90,100,0
LCOMB,28,97,0
!
选择所有
ALLSEL,ALL
!
选择体（支撑辊）
VSEL,S,,,2,9,
!
在选择的体下体编辑
ALLSEL,BELOW,VOLU
!
对体进行GLUE
VGLUE,2,3,4,5,6,7,8,9
!
对面进行ADD
AADD,25,66
AADD,33,74
AADD,40,78
AADD,17,52
!
对线进行ADD
LCOMB,36,102,0
LCOMB,46,103,0
LCOMB,56,104,0
LCOMB,20,101,0
!
选择所有
ALLSEL,ALL
!
给体定义单元属性（轧件）
CM,_Y,VOLU
VSEL,,,,1
CM,_Y1,VOLU
CMSEL,S,_Y
CMSEL,S,_Y1
VATT,1,,1,0
CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
!
给体定义单元属性（工作辊辊身）
FLST,5,4,6,ORDE,2
FITEM,5,18
FITEM,5,-21
CM,_Y,VOLU
VSEL,,,,P51X
CM,_Y1,VOLU
CMSEL,S,_Y
CMSEL,S,_Y1
VATT,2,,1,0
CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
!
给体定义单元属性（支撑辊辊身）
FLST,5,4,6,ORDE,4
FITEM,5,10
FITEM,5,12
FITEM,5,14
FITEM,5,16
CM,_Y,VOLU
VSEL,,,,P51X
CM,_Y1,VOLU
CMSEL,S,_Y
CMSEL,S,_Y1
VATT,3,,1,0
CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
!
给体定义单元属性（工作辊辊径）
FLST,5,4,6,ORDE,4
FITEM,5,11
FITEM,5,13
FITEM,5,15
FITEM,5,17
CM,_Y,VOLU
VSEL,,,,P51X
CM,_Y1,VOLU
CMSEL,S,_Y
CMSEL,S,_Y1
VATT,4,,1,0
CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
!
给体定义单元属性（支撑辊辊径）
FLST,5,4,6,ORDE,4
FITEM,5,3
FITEM,5,5
FITEM,5,7
FITEM,5,9
CM,_Y,VOLU
VSEL,,,,P51X
CM,_Y1,VOLU
CMSEL,S,_Y
CMSEL,S,_Y1
VATT,5,,1,0
CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
!
定义单元尺寸
LESIZE,9,,,20,,,,,1
LESIZE,10,,,20,,,,,1
LESIZE,11,,,20,,,,,1
LESIZE,12,,,20,,,,,1
LESIZE,13,,,20,,,,,1
LESIZE,15,,,20,,,,,1
LESIZE,21,,,20,,,,,1
LESIZE,23,,,20,,,,,1
LESIZE,29,,,20,,,,,1
LESIZE,39,,,20,,,,,1
LESIZE,49,,,20,,,,,1
LESIZE,63,,,20,,,,,1
LESIZE,73,,,20,,,,,1
LESIZE,83,,,20,,,,,1
LESIZE,17,,,1,,,,,1
LESIZE,19,,,1,,,,,1
LESIZE,25,,,1,,,,,1
LESIZE,27,,,1,,,,,1
LESIZE,34,,,1,,,,,1
LESIZE,44,,,1,,,,,1
LESIZE,54,,,1,,,,,1
LESIZE,68,,,1,,,,,1
LESIZE,78,,,1,,,,,1
LESIZE,88,,,1,,,,,1
LESIZE,2,0.01,,,,,,,1
LESIZE,1,,,4,,,,,1
!
用MESH
划分轧件
VMESH,1
!
定义单元尺寸
LESIZE,24,,,10,,,,,1
LESIZE,65,,,10,,,,,1
LESIZE,75,,,10,,,,,1
LESIZE,85,,,10,,,,,1
LESIZE,66,0.019,,,,,,,1
LESIZE,76,0.019,,,,,,,1
LESIZE,86,0.019,,,,,,,1
LESIZE,93,0.019,,,,,,,1
!
用VSWEEP划分工作辊辊身
VSWEEP,18
VSWEEP,19
VSWEEP,20
VSWEEP,21
!
选择智能控制
SMRT,6
!
用VSWEEP划分工作辊辊径
VSWEEP,11
VSWEEP,13
VSWEEP,15
VSWEEP,17
!
关闭智能控制
SMRT,OFF
!
定义单元尺寸
LESIZE,32,0.019*1.9,,,,,,,1
LESIZE,42,0.019*1.9,,,,,,,1
LESIZE,52,0.019*1.9,,,,,,,1
LESIZE,59,0.019*1.9,,,,,,,1
LESIZE,16,,,15,,,,,1
LESIZE,31,,,15,,,,,1
LESIZE,41,,,15,,,,,1
LESIZE,51,,,15,,,,,1
!
用VSWEEP划分支撑辊辊身
VSWEEP,10
VSWEEP,12
VSWEEP,14
VSWEEP,16
!
选择智能控制
SMRT,6
!
用SWEEP划分支撑辊辊径
VSWEEP,3
VSWEEP,5
VSWEEP,7
VSWEEP,9
!
生成PART
EDPART,CREATE
!
定义接触
EDCGEN,ASTS,1,2,0.35,0.25,0,0,0,,,,,0,10000000,0,0
EDCGEN,ASTS,2,4,0.35,0.25,0,0,0,,,,,0,10000000,0,0
!
选择体（轧件）
VSEL,S,,,1
!
选择体中的节点
NSLV,S,1
!
创建Component
CM,SHEET,NODE
!
选择Y＝0的节点
NSEL,R,LOC,Y,0
!
定义约束
D,ALL,,0,,,,UY,,,,,
!
选择全部
ALLSEL,ALL
!
选择Z＝0的节点
NSEL,R,LOC,Z,0
!
定义约束
D,ALL,UZ,0
!
选择全部
ALLSEL,ALL
!
选择Z＝-0.001的节点
NSEL,R,LOC,Z,-0.001
!
定义约束
D,ALL,UZ,0
!
选择全部
ALLSEL,ALL
!
定义初速度
EDVE,VELO,SHEET,1.5,0,0,0,0,0,,,,,,
!
定义时间变量
*DIM,TIME,ARRAY,2,1,1,,,
!
定义速度变量
*DIM,VELOCITY,ARRAY,2,1,1,,,
!
赋值时间变量
*SET,TIME（2,1,1）,10
!
赋值速度变量
*SET,VELOCITY（1,1,1）,4.21
!
赋值速度变量
*SET,VELOCITY（2,1,1）,4.21
!
施加荷载
EDLOAD,ADD,RBOZ,0,3,TIME,VELOCITY,0,,,,,
!
进入SOLUTION
FINISH
/SOL
!
定义输出内存
EDSTART,0,85000000,7,
!
定义能量模式
EDENERGY,1,1,1,1
!
定义结束时间
TIME,0.4,
!
定义输出格式
EDOPT,ADD,blank,BOTH
!
定义.RST文件输出步数
EDRST,20,
!
定义.HST文件输出步数
EDHTIME,100,
!
输出ASCII文件
EDOUT,RCFORC
!
保存
SAVE
!
求解
SOLVE
!
设置显示比例
/DSCALE,1,1.0
!
进入POST1
FINISH
/POST1
!
设置观察步数
SET,LAST
!
显示应变分布情况
PLNSOL,EPTO,EQV,0,1.0
!
显示动画过程
PLNS,S,EQV
ANDATA,0.5,,0,0,0,1,0,1
!
退出ANSYS
/REPLOT
FINISH
/EXIT,NOSAV
升级会员 