弹塑性材料杆件撞击刚性墙面——瞬态分析

芷行说 2021年2月3日浏览：2435 收藏：1

本案例为冲击的一个简单案例，主要目的在于启发对冲击仿真有需求的朋友，但又对冲击不太了解，本案例可促进入门。

弹塑性材料杆件撞击刚性墙面——瞬态分析的图1

关键设置：

冲击下的接触单元/目标单元设置

et,2,TARGE170    ! 刚性墙为目标单元170
et,3,CONTA174    ! 圆杆接触面为接触单元174

keyopt,3,2,0     ! 采用增广拉格朗日公式

keyopt,3,4,2     ! 接触检测点方向为目标面节点方向，在刚性冲击问题下需要进行该设置 
!以提高收敛性，尤其在接触面几何不规则，该设置尤为重要，这是因为接触面几何不规则，
!会产生不对称的接触力，破坏收敛性。

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   !节点上指定初始条件,可指定初始位移、速度

瞬态分析控制

关于求解器的计算方法和积分算法读者可自行根据问题所需进行设置，对比计算结果。

/soluantype,trans         ! Perform a transient 
analysisnlgeom,on             ! 大变形

trnopt,full, , , , ,HHT  ! 指定全分析计算方法和HHT时间积分算法
tintp,0.1                ! 积分阻尼指定，默认为0.005
time,80e-6               ! 计算时间             
nsub,100,10000,100       ! 计算步设置

outres,all,all           ! 结果输出  

solve

finish