本帖最后由 gdyu_yu 于 2010-8-11 13:52 编辑
我做的一个简单的刚性磙子(刚体)在平面上滚压的例子,欲算出平面部分的应力状态以及磙子所受的支反力。程序和示意图如下,为什么显示不出等效应力的分布情况呢。请高手指点,多谢了!
程序:
/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 |