薄钢板冲压分析

/FILNAME,Analysis of Sheet Metal,0
/PREP7
L=0.1 !钢板的宽度
H=0.004 !钢板的厚度
L1=0.08 !凹模的半长
H1=0.02 !凹模的深度
L2=L1+H !凸模的半长
H2=H1 !凸模的深度
L3=0.04 !凹模的底半长
R1=0.009 !模具的圆角
R2=0.009 !模具的圆角
R3=0.009 !模具的圆角
R4=0.009 !模具的圆角


ET,1,PLANE42 !定义单元
KEYOPT,1,3,2 !设置单元关键字，指定单元模拟平面应变行为
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,2E11 !定义材料弹性模量
MPDATA,PRXY,1,,0.3 !定义材料泊松比
MPDATA,DENS,1,,7580 !定义密度
TB,BISO,1,1,2, !设置材料为双线性等向强度
TBTEMP,0
TBDATA,,2E8,2E9,,,, !定义材料的屈服应力和切线模量


!定义关键点
K,1,0,0,0
K,2,L1,0,0
K,3,L1,-H1,0
K,4,L1+L3,-H1,0
K,5,0,H2+H,0
K,6,L2,H2+H,0
K,7,L2,H,0
K,8,L1+L3,H,0
!通过关键点定义线
L,1,2
L,2,3
L,3,4
L,5,6
L,6,7
L,7,8

RECTNG,L1+L3-L,L1+L3,0,H, !定义矩形
!定义模具的圆角
LFILLT,1,2,R1, ,
LFILLT,2,3,R2, ,
LFILLT,4,5,R3, ,
LFILLT,5,6,R4, ,

ESIZE,0.005/5 !设置总体单元尺寸
MSHKEY,1 !设置网格划分方法为映射网格
AMESH,1 !划分面的网格

!使用接触向导定义带控制节点的接触对
CM,_NODECM,NODE
CM,_ELEMCM,ELEM
CM,_KPCM,KP
CM,_LINECM,LINE
CM,_AREACM,AREA
CM,_VOLUCM,VOLU
/GSAV,cwz,gsav,,temp
MP,MU,1,0.3
MAT,1
R,3
REAL,3
ET,2,169
ET,3,172
KEYOPT,3,9,0
KEYOPT,3,10,2
R,3,
RMORE,
RMORE,,0
RMORE,0
! Generate the target surface
LSEL,S,,,4
LSEL,A,,,5
LSEL,A,,,6
LSEL,A,,,13
LSEL,A,,,14
CM,_TARGET,LINE
TYPE,2
LATT,-1,3,2,-1
TYPE,2
LMESH,ALL
! Create a pilot node
KSEL,S,,,5
KATT,-1,3,2,-1
KMESH,5
! Generate the contact surface
LSEL,S,,,9
CM,_CONTACT,LINE
TYPE,3
NSLL,S,1
ESLN,S,0
ESURF
*SET,_REALID,3
ALLSEL
ESEL,ALL
ESEL,S,TYPE,,2
ESEL,A,TYPE,,3
ESEL,R,REAL,,3
LSEL,S,REAL,,3
/PSYMB,ESYS,1
/PNUM,TYPE,1
/NUM,1
EPLOT
ESEL,ALL
ESEL,S,TYPE,,2
ESEL,A,TYPE,,3
ESEL,R,REAL,,3
LSEL,S,REAL,,3
CMSEL,A,_NODECM
CMDEL,_NODECM
CMSEL,A,_ELEMCM
CMDEL,_ELEMCM
CMSEL,S,_KPCM
CMDEL,_KPCM
CMSEL,S,_LINECM
CMDEL,_LINECM
CMSEL,S,_AREACM
CMDEL,_AREACM
CMSEL,S,_VOLUCM
CMDEL,_VOLUCM
/GRES,cwz,gsav
CMDEL,_TARGET
CMDEL,_CONTACT
!*
!*
!*
CM,_NODECM,NODE
CM,_ELEMCM,ELEM
CM,_KPCM,KP
CM,_LINECM,LINE
CM,_AREACM,AREA
CM,_VOLUCM,VOLU
/GSAV,cwz,gsav,,temp
MP,MU,1,0.3
MAT,1
R,4
REAL,4
ET,4,169
ET,5,172
KEYOPT,5,9,0
KEYOPT,5,10,2
R,4,
RMORE,
RMORE,,0
RMORE,0
! Generate the target surface
LSEL,S,,,1
LSEL,A,,,2
LSEL,A,,,3
LSEL,A,,,11
LSEL,A,,,12
CM,_TARGET,LINE
TYPE,4
LATT,-1,4,4,-1
TYPE,4
LMESH,ALL
! Create a pilot node
KSEL,S,,,1
KATT,-1,4,4,-1
KMESH,1
! Generate the contact surface
LSEL,S,,,7
CM,_CONTACT,LINE
T

YPE,5
NSLL,S,1
ESLN,S,0
ESURF
*SET,_REALID,4
ALLSEL
ESEL,ALL
ESEL,S,TYPE,,4
ESEL,A,TYPE,,5
ESEL,R,REAL,,4
LSEL,S,REAL,,4
/PSYMB,ESYS,1
/PNUM,TYPE,1
/NUM,1
EPLOT
! Reverse target normals
FLST,5,5,4,ORDE,4
FITEM,5,1
FITEM,5,-3
FITEM,5,11
FITEM,5,-12
CM,_Y,LINE
LSEL, , , ,P51X
CM,_YEL,ELEM
CM,_YND,NODE
NSLL,S,1
ESLN,S,1
ESEL,R,REAL,,_REALID
ESURF,,REVERSE
CMSEL,S,_Y
CMSEL,S,_YEL
CMSEL,S,_YND
CMDELE,_Y
CMDELE,_YEL
CMDELE,_YND
/REPLOT
!*
ESEL,ALL
ESEL,S,TYPE,,4
ESEL,A,TYPE,,5
ESEL,R,REAL,,4
LSEL,S,REAL,,4
/PSYMB,ESYS,1
/PNUM,TYPE,1
/NUM,1
EPLOT
ESEL,ALL
ESEL,S,TYPE,,4
ESEL,A,TYPE,,5
ESEL,R,REAL,,4
LSEL,S,REAL,,4
CMSEL,A,_NODECM
CMDEL,_NODECM
CMSEL,A,_ELEMCM
CMDEL,_ELEMCM
CMSEL,S,_KPCM
CMDEL,_KPCM
CMSEL,S,_LINECM
CMDEL,_LINECM
CMSEL,S,_AREACM
CMDEL,_AREACM
CMSEL,S,_VOLUCM
CMDEL,_VOLUCM
/GRES,cwz,gsav
CMDEL,_TARGET
CMDEL,_CONTACT
!设置接触单元单元属性
KEYOPT,3,4,2
KEYOPT,3,10,2

KEYOPT,5,4,2
KEYOPT,5,10,2
!进入求解
/SOL
ANTYPE,0 !设置分析类型为静力学分析
NLGEOM,1 !设置计算考虑大变形
NROPT,FULL, ,ON !设置采用完全牛顿法

DL,8, ,SYMM !设置对称边界条件
DK,1, , , ,0,ALL, , , , , , !约束关键点1的所有自由度
DK,5, , , ,0,UX,ROTZ, , , , ,!约束关键点5的UX和ROTZ
DK,5, ,-H1, ,0,UY, , , , , , !施加关键点5的Y方向位移为-H1
OUTRES,ALL,ALL, !设置输出所有子步
TIME,1 !定义计算时间步
AUTOTS,-1
NSUBST,500 !设置计算子步为500
KBC,0 !设置载荷为斜坡载荷
LSWRITE,1 !写入载荷文件1

DK,5, ,H1*0.2, ,0,UY, , , , , ,!施加关键点5的Y方向位移为0.2*H2
TIME,2 !定义计算时间步
NSUBST,300 !设置计算子步为300
KBC,0
lswrite,2 !写入载荷文件2
LSSOLVE,1,2,1 !按载荷文件数求解。