冲击响应计算.docx

冲击响应计算
一:模型
由PRO/E模型导入到ASNSYWORKBENCH中进行简化与网格划分,再导入到ANSYS中进行冲击响应分析。导入到ANSYS后的模型如下图。

二:冲击激励
利用ANSYS的瞬态冲击分析功能进行冲击响应分析，采用国军标冲击试验环境中常用的半正弦脉冲作为激励，其相关参数如图a。在分析过程中将半正弦加速度激励转换成位移激励（如图5.19b所示），沿Z轴正方向施加在分析对象的基础上。

a半正弦加速度激励b积分得到的位移激励
三:应力和应变的时域响应结果
从模型中选取节点82，绘制节点82的应力与应变随时间变化曲线图，如下图。

节点82的Z方向应变曲线图

节点82的Z方向应力曲线图
从节点82的应力应变曲线图可以看出，应力和应变以相同的频率振荡，并且两者的最大值的时刻几乎相同。随着仿真时间的增长，等效应力和等效应变的响应也应该衰减，但是由于系统的阻尼较小，两者响应的衰减作用并不明显。
四：应力和应变的分布
为了了解在激励作用下分析对象中的等效应力和等效应变的分布情况，提取了在时刻（0.033）时模型某部分Z方向应力和Z方向应变的云图结果。如下图。

Z方向应力云图

Z方向应变云图
五：整体应力和应变的分布
为了了解在激励作用下模型整体的等效应力和等效应变的分布情况，提取了在时刻（0.033）时模型Z方向应力和Z方向应变的云图结果。如下图。

整体的Z方向应力云图

整体的Z方向应变云图
六:总结
1,冲击载荷加载在模型各四个角上的一个节点上.
2,对于此模型,由于在进行冲击分析时ANSYS求解每一个子步需要约30分钟,所以在离散半正弦脉冲时只是离散了33个节点进行求解,所需时间大约为13个小时求解.
3,求解结果应变比较小,最大值大约为1mm.应力也比较小,最大值不超过1Mpa