麻烦王老师给看一下,多谢了
finish
/cle
/NOPR
/PMETH,OFF,0
KEYW,PR_SET,1
KEYW,PR_STRUC,1
KEYW,LSDYNA,1
!******************************
/UNITS,SI
!******************************
/plopts,info,0
!/view,,1,1,1
/VIEW, 1 ,,-1
dw=0.001 !1 millimetre毫米 =0.03937 inch 英寸
C_PI=3.1415926 !PI的值
RHO=2800 !铝合金的密度
VTOT=0 !活塞的体积初始化
WT=0 !活塞的质量初始化
SURF=1E7 !压力
C_R=85*dw/2 !C_R为底面半径
C_F=SURF*C_PI*C_R*C_R !C_F压力
ACAZ=0 !Z方向的加速度
/prep7
ET,1,SHELL163 !铝合金和钢的二维有限元模型的et
ET,2,solid164 !铝合金的solid164
ET,3,solid164 !钢的solid164 刚体
MP,DENS,1,2800 !aluminum
MP,EX,1,72e9
MP,NUXY,1,0.33
TB,BKIN,1,,,,
TBDAT,1, 480e6
TBDAT,2, 27e9
MP,DENS,2,7800
MP,EX,2,200e9
MP,NUXY,2,0.3
EDMP,RIGI,2,7,7
wpro,,90.000000,
!建模活塞
RECTNG,0,85/2000,20/1000+dw,31/1000+dw
RECTNG,83/2000,85/2000,0/1000+dw,20/1000+dw
ADEL,ALL
LANG,1,5,-(90-5.5)
LFILLT,1,10,5/1000
LDELE,7,9,1
LDELE,2,6,4
L,3,6
LANG,3,10,90
LANG,2,11,90
LDIV,11
L,3,14
AL,4,1,7,6
AL,7,11,13,3
AL,12,9,2,13
AL,10,5,8,9
CM,HUOSAI,AREA
!分格活塞二维模型
TYPE, 1
mat,1
MSHAPE, 0, 2D
MSHKEY, 1
ESIZE, ,5
AMESH, HUOSAI
!建模分格活塞二维模型finish
K,1000, 0
K,1001, 0, 0,5
TYPE, 2
mat,1
ESIZE, ,20
VROTAT, ALL, , , , , , 1000, 1001, 360,4
CMSEL, S,HUOSAI
ACLEAR, ALL
ALLSEL
EPLOT
!/eof
!1
!***********************
*DO,I,101,8100
*GET,VCAI,ELEM,I,VOLU
VTOT=VTOT+VCAI !总的体积
*enddo
WT=RHO*VTOT !质量WT
!**********************
ACAZ=C_F/WT !Z方向的加速度
!/eof
!2
!建模制动面
wpro,,-90.000000,
k,1011,85.1*dw/2,0,0*dw
k,1012,108*dw/2,0,0*dw
k,1013,80*dw/2,0,-25*dw
k,1014,82*dw/2,0,-25*dw
k,1015,108*dw/2,0,-25*dw
k,1016,80*dw/2,0,-47*dw
k,1017,82*dw/2,0,-47*dw
k,1018,108*dw/2,0,-47*dw
/pnum,kp,1
/pnum,line,1
/pnum,area,1
gplot
l,1011,1012
l,1012,1015
l,1015,1014
l,1014,1011
l,1013,1014
l,1014,1017
l,1017,1016
l,1016,1013
l,1017,1018
l,1018,1015
!/eof
!3
al,82,83,84,85 !area65
al,84,91,90,87 !area66
al,86,87,88,89 !area67
asel,s,,,65
asel,a,,,66
asel,a,,,67
CM,ZHIDONG,AREA
!/eof
!4
TYPE, 1
mat,2
ESYS, 0
MSHAPE, 0, 2D
MSHKEY, 1
ESIZE, ,10
AMESH, ZHIDONG
TYPE, 3
mat,2
ESIZE, ,20
VROTAT, ALL, , , , , , 1000, 1001, 360,4
CMSEL, S,ZHIDONG
ACLEAR, ALL
ALLSEL
/type,1,8
EPLOT
/eof
ESEL,S,MAT,,1
NSLE,S
CM,BOX,NODE
ALLSEL,ALL
ESEL,S,MAT,,2
NSLE,S
CM,TABLE,NODE
ALLSEL,ALL
!EDPART,CREATE
!生成接触ASSC
EDCGEN,ASSC
!EDCONTACT,1.0,,,,,,2,,0
edclist,all
finish
/solu
!修改了速度将公制改为英制 速度方向从y方向转为z方向
EDVE,VELO,BOX,0,0,-200,0,0,0, , , , , ,
*DIM,time,ARRAY,2,1,1, , ,
*DIM,accc,ARRAY,2,1,1, , ,
*SET,TIME(1,1,1) , 0
*SET,TIME(2,1,1) , 3.0e-4
*SET,ACCC(1,1,1) , ACAZ
*SET,ACCC(2,1,1) , ACAZ
EDLOAD,ADD,ACLZ,0,BOX,TIME,ACCC, 0, , , , ,
TIME,3.0e-4, !终止时间
EDRST,50,
EDHTIME,50,
EDDUMP,1,
EDENERGY,1,1,1,1
edhist,box
edout,all
nsel,all
allsel,all
/replot
finish
!/eof
!
/SOLU
/STAT,SOLU
SOLVE
FINISH
/EXIT,ALL |