2024年5月2日发(作者：)

转载一个例子，说明EDLCS和EDLOAD命令。另结合EDMP可以实现刚体绕过质点

的任意轴恒定转动，即自转。

在对刚体施加不过其质心的任意轴恒定转动时候，暂时遇到了困难，还不知道如

何施加。知道的朋友可以留言告诉一声。

如：实现刚体自转和公转时候，利用EDLCS和EDLOAD及其EDMP可以实现其恒定

速度的自转，却无法实现恒定速度的公转。不过可以利用EDPVEL施加初始速度

的公转。

!采用国际单位制

/prep7

!define the geometry parameters of model

r1=0.10 !第一道轧制的轧辊半径

r2=0.09 !第二道次轧制的轧辊半径

L=0.5 !轧件长度，示意性质，只用于足够的模拟计算长度即可

w1=0.05 !第一个辊子的工作长度

w2=0.06 !第二个工作辊的工作长度

w=0.04 !轧件宽度

D=0.22 !第一个轧辊和第二个轧辊之间的距离

H1=0.03 !轧件原始厚度

H2=0.025 !第一道次轧制后轧件的厚度

H3=0.02 !第二道次轧制后轧件的厚度

rv=600 !工作辊旋转速度

T1=1000 !前张力

T0=0 !后张力

rv=6.28 !定义轧辊的转速，以线速度表示

!define the material property parameters for the strip

dens1=7830

ex1=2.1e11

nuxy1=0.32

yieldstr1=1.09e7

tangentM1=0

!define the material property parameters for the rigid workroll

dens2=7930

ex2=3e11

nuxy2=0.3

!define real constants

r,1

R,2

RMODIF,2,1, , ,0.002, , , ,

R,3

RMODIF,3,1, , ,0.002, , , ,

!define element type used

ET,1,SOLID164

ET,2,shell163

ET,3, shell163

!define material property for the strip

mp,dens,1,dens1

mp,ex,1,ex1

mp,nuxy,1,nuxy1

tb,biso,1,,,,

tbdat,1,yieldstr1

tbdat,2,tangentM1

!define material property for rigid workroll

EDMP,RIGI,2,7,4

MP,DENS,2,dens2

MP,EX,2,ex2

MP,NUXY,2,nuxy2

EDMP,RIGI,3,7,4

MP,DENS,3,dens2

MP,EX,3,ex2

MP,NUXY,3,nuxy2

!define the model

!define the model of the first workroll

k,1,L+r1,H2,-w1/2

K,2,L+r1,H2,w1/2

l,1,2

k,3,L+r1,H2+r1,-w1/2

k,4,l+r1,H2+r1,w1/2

arotat,1,,,,,,3,4,360

k,11,0,0,-w/2

k,12,0,0,w/2

k,14,0,H1,-w/2

K,13,0,H1,w/2

l,11,12

l,12,13

l,13,14

l,14,11

al,13,14,15,16

vext,5,,,l,

!define the model of the second workroll

k,30,L+r1+D,H3,-w2/2

k,31,L+r1+D,H3,w2/2

l,30,31

k,32,L+r1+D,H3+r2,-w2/2