关键设置:
冲击下的接触单元/目标单元设置
et,2,TARGE170 ! 刚性墙为目标单元170
et,3,CONTA174 ! 圆杆接触面为接触单元174
! 采用增广拉格朗日公式
! 接触检测点方向为目标面节点方向,在刚性冲击问题下需要进行该设置
!以提高收敛性,尤其在接触面几何不规则,该设置尤为重要,这是因为
!接触面几何不规则,会产生不对称的接触力,破坏收敛性。
keyopt,3,7,4 ! 冲击约束下的时间载荷步控制,该设置在冲击问题下尤为重要,严重影响
!结果的收敛性
keyopt,3,10,2 ! 每一载荷步计算完成后,接触刚度更新
弹塑性材料设置
mp,ex,1,117e9 ! Young's Modulus for copper bar (Pa)
mp,nuxy,1,0.35 ! Poisson's ration
mp,dens,1,8930 ! Density of copper
tb,biso,1 ! Bilinear isotropic definition
tbdata,1,400e6 ! 屈服应力指定
tbdata,2,100e6 ! 切线模量指定
初始速度加载
ic,all,uz,,-227 !节点上指定初始条件,可指定初始位移、速度
瞬态分析控制
关于求解器的计算方法和积分算法读者可自行根据问题所需进行设置,对比计算结果。
/solu
antype,trans ! Perform a transient analysis
nlgeom,on ! 大变形
trnopt,full, , , , ,HHT ! 指定全分析计算方法和HHT时间积分算法
tintp,0.1 ! 积分阻尼指定,默认为0.005
time,80e-6 ! 计算时间
nsub,100,10000,100 ! 计算步设置
outres,all,all ! 结果输出
solve
finish