一个简单的刚性磙子(刚体)在平面上滚压的例子<BR><BR><BR>程序:<BR>/PREP7 <BR><BR>L_plate=0.02 !定义平板的长度(20mm)<BR>H_plate=0.002 !定义平板的厚度( 2mm)<BR>W_plate=0.01<BR><BR>R_wheel=0.005 !滚轮的半径( 5mm )<BR>W_wheel=0.002<BR><BR>dis1=0.00015 !单次向下进给量(0.15mm)<BR><BR>rot_wheel=-15 !滚轮的旋转角速度(旋转中心为z轴)<BR>v_plate=R_wheel*rot_wheel !平板的初始进入速度(x轴负方向;单位:m/s)<BR><BR>!定义单元类型<BR>ET,1,solid164 !定义单元类型<BR><BR>!定义材料类型<BR>MP,DENS,1,7800 !定义材料的密度 :kg/m3<BR>MP,EX,1,0.7e11 !定义材料的弹性模量:pa<BR>MP,NUXY,1,0.334 !定义材料的泊松比 <BR>TB,bkin,1,,,, !定义材料为双线性随动强化模型<BR>TBDAT,1,3.44e7 !屈服强度 a<BR>TBDAT,2,1.38e8 !剪切弹性模量 a <BR><BR>EDMP,RIGI,2,7,4 !定义第二类材料模型(刚体,约束所有的平动自由度和沿X、Y轴的转动自由度) <BR>MP,DENS,2,7900 !定义材料的密度 :kg/m3<BR>MP,EX,2,2E+011 !定义材料的弹性模量:pa<BR>MP,NUXY,2,0.292 !定义材料的泊松比 <BR><BR>!建立模型<BR>BLOCK,0.0013,L_plate+0.0013,-R_wheel+dis1,-(R_wheel+H_plate)+dis1,-W_plate/2,W_plate/2<BR>WPOFFS,,,-W_wheel/2<BR>CYL4,,,R_wheel,,,360,W_wheel !CYL4, XCENTER, YCENTER, RAD1, THETA1, RAD2, THETA2, DEPTH <BR><BR>!平板划分网格<BR>TYPE, 1 <BR>MAT, 1<BR>REAL, <BR>ESYS, 0 <BR>SECNUM, <BR>!* <BR>FLST,5,3,4,ORDE,3 <BR>FITEM,5,5 <BR>FITEM,5,-6 <BR>FITEM,5,10 <BR>CM,_Y,LINE <BR>LSEL, , , ,P51X <BR>CM,_Y1,LINE <BR>CMSEL,,_Y <BR>!* <BR>LESIZE,_Y1,0.0005, , , , , , ,1 <BR>!* <BR>MSHAPE,0,3D <BR>MSHKEY,1<BR>!* <BR>CM,_Y,VOLU <BR>VSEL, , , , 1 <BR>CM,_Y1,VOLU <BR>CHKMSH,'VOLU' <BR>CMSEL,S,_Y <BR>!* <BR>VMESH,_Y1 <BR>!* <BR>CMDELE,_Y <BR>CMDELE,_Y1 <BR>CMDELE,_Y2 <BR><BR>!滚轮划分网格<BR>TYPE, 1 <BR>MAT, 2<BR>REAL, <BR>ESYS, 0 <BR>SECNUM, <BR>!* <BR>ESIZE,0.0005,0, <BR>FLST,2,1,5,ORDE,1 <BR>FITEM,2,8 <BR>AESIZE,P51X,0.001, <BR>CM,_Y,VOLU <BR>VSEL, , , , 2 <BR>CM,_Y1,VOLU <BR>CHKMSH,'VOLU' <BR>CMSEL,S,_Y <BR>!* <BR>VSWEEP,_Y1 <BR>!* <BR>CMDELE,_Y <BR>CMDELE,_Y1 <BR>CMDELE,_Y2 <BR><BR>!生成part<BR>EDPART,CREATE <BR><BR>!*定义接触 <BR><BR>EDCGEN,FSTS, 2, 1,0.3,0.1,0,0,0, , , , ,0,10000000,0,0 <BR>!施加约束<BR>FLST,2,1,5,ORDE,1 <BR>FITEM,2,3 <BR>!* <BR>/GO <BR>DA,P51X,UY, <BR><BR>!施加初速度<BR>EDPV,VGEN, 1,v_plate,0,0,0, , ,0,0,0,0,0,0, <BR>!施加载荷<BR>*DIM,TIME,ARRAY,2,1,1, , , ! <BR>*DIM,DIS,ARRAY,2,1,1, , , ! <BR>*SET,TIME(2,1,1) , 0.2 !给时间数组的第2行第1 <BR>*SET,DIS(1,1,1) , rot_wheel !给速度数组的第1行第1 <BR>*SET,DIS(2,1,1) , rot_wheel !给速度数组的第2行第1 <BR>EDLOAD,ADD,RBOZ,0, 2,TIME,DIS, 0, , , , , <BR>FINISH <BR><BR>/SOL<BR>TIME,0.2 ! <BR>EDRST,10 !定义写.rst<BR>EDHTIME,10 ! <BR>EDDUMP,1 ! <BR>!* <BR>ALLSEL,ALL ! <BR>SAVE ! <BR>solve <BR>FINISH <BR><BR>/POST1 <BR><BR>SET,1,LAST,1, <BR>!* <BR>PLNS,S,EQV <BR>ANDATA,0.5, ,0,0,0,1,0,1 <BR> |