刚性磙子（刚体）在平面上滚压问题

sweetboy9063

我做的一个简单的刚性磙子（刚体）在平面上滚压的例子，欲算出平面部分的应力状态以及磙子所受的支反力。程序和示意图如下，为什么显示不出等效应力的分布情况呢。请高手指点，多谢了！

程序：
/PREP7

L_plate=0.02 !定义平板的长度（20mm）
H_plate=0.002 !定义平板的厚度（ 2mm）
W_plate=0.01

R_wheel=0.005 !滚轮的半径（ 5mm ）
W_wheel=0.002

dis1=0.00015 !单次向下进给量（0.15mm）

rot_wheel=-15 !滚轮的旋转角速度（旋转中心为z轴）
v_plate=R_wheel*rot_wheel !平板的初始进入速度(x轴负方向；单位：m/s)

!定义单元类型
ET,1,solid164 !定义单元类型

!定义材料类型
MP,DENS,1,7800 !定义材料的密度 :kg/m3
MP,EX,1,0.7e11 !定义材料的弹性模量:pa
MP,NUXY,1,0.334 !定义材料的泊松比
TB,bkin,1,,,, !定义材料为双线性随动强化模型
TBDAT,1,3.44e7 !屈服强度 a
TBDAT,2,1.38e8 !剪切弹性模量 a

EDMP,RIGI,2,7,4 !定义第二类材料模型(刚体,约束所有的平动自由度和沿X、Y轴的转动自由度)
MP,DENS,2,7900 !定义材料的密度 :kg/m3
MP,EX,2,2E+011 !定义材料的弹性模量:pa
MP,NUXY,2,0.292 !定义材料的泊松比

!建立模型
BLOCK,0.0013,L_plate+0.0013,-R_wheel+dis1,-(R_wheel+H_plate)+dis1,-W_plate/2,W_plate/2
WPOFFS,,,-W_wheel/2
CYL4,,,R_wheel,,,360,W_wheel !CYL4, XCENTER, YCENTER, RAD1, THETA1, RAD2, THETA2, DEPTH

！平板划分网格
TYPE, 1
MAT, 1
REAL,
ESYS, 0
SECNUM,
!*
FLST,5,3,4,ORDE,3
FITEM,5,5
FITEM,5,-6
FITEM,5,10
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1,0.0005, , , , , , ,1
!*
MSHAPE,0,3D
MSHKEY,1
!*
CM,_Y,VOLU
VSEL, , , , 1
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
!*
VMESH,_Y1
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2

！滚轮划分网格
TYPE, 1
MAT, 2
REAL,
ESYS, 0
SECNUM,
!*
ESIZE,0.0005,0,
FLST,2,1,5,ORDE,1
FITEM,2,8
AESIZE,P51X,0.001,
CM,_Y,VOLU
VSEL, , , , 2
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
!*
VSWEEP,_Y1
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2

！生成part
EDPART,CREATE

!*定义接触

EDCGEN,FSTS, 2, 1,0.3,0.1,0,0,0, , , , ,0,10000000,0,0
!施加约束
FLST,2,1,5,ORDE,1
FITEM,2,3
!*
/GO
DA,P51X,UY,

!施加初速度
EDPV,VGEN, 1,v_plate,0,0,0, , ,0,0,0,0,0,0,
!施加载荷
*DIM,TIME,ARRAY,2,1,1, , , !
*DIM,DIS,ARRAY,2,1,1, , , !
*SET,TIME(2,1,1) , 0.2 !给时间数组的第2行第1
*SET,DIS(1,1,1) , rot_wheel !给速度数组的第1行第1
*SET,DIS(2,1,1) , rot_wheel !给速度数组的第2行第1
EDLOAD,ADD,RBOZ,0, 2,TIME,DIS, 0, , , , ,
FINISH

/SOL
TIME,0.2 !
EDRST,10 !定义写.rst
EDHTIME,10 !
EDDUMP,1 !
!*
ALLSEL,ALL !
SAVE !
solve
FINISH

/POST1

SET,1,LAST,1,
!*
PLNS,S,EQV
ANDATA,0.5, ,0,0,0,1,0,1

sweetboy9063

大家来看看！帮帮我

hyfhyf88

这是以前做的，看对你是否有帮助？

sweetboy9063

不错不错。请问如何得到你那个辊子在整个辊压过程中所受到的平板的向上的反作用力和时间的关系曲线。

sweetboy9063

to:hyfhyf88

你好！请问你做的这个东西，你对做的结果的可信度能说出多少，你用的材料模型是什么样的，摩擦系数选用多少啊！