冲击力仿真计算与试验研究

本资料包含pdf文件1个，下载需要1积分

抗 冲击性能是材料 的重要属性之- ，而材料的抗冲击性能要依靠冲击试验测得。冲击试验应用的设备是材料冲击试验机，通过摆锤冲击试样后得出的冲击吸收功和冲击力来衡量材料的抗冲击性能。

冲击过程是个非常复杂的过程，本文根据接触力学理论，应用虚拟平台，对材料冲击试验机冲击过程进行了仿真模拟和计算，得出了冲击力的具体曲线，并和实际的冲击试验数据进行 了对 比，分析总结 了两者的差别 ，为更深人 的研究打下基础 。

2 碰撞力模型根据 Hertz contact theory，采用非线性等效弹簧阻尼模型作为接触力的计算模型。当计算两个构件之间的接触力时，接触力由两个部分组成：-个是由于两个构件之间的相互切入而产生的弹性力；另- 个是由相对速度产生的阻尼力。其广义形式可以表示为：F - K C其中， 为法向接触力，单位为 N；K(Stif～ness)为 Hertz接触刚度，表示接触表面的刚度，单位为 N/mm； (Penetration Depth)为接触点的法向穿透深度，单位为 mm；e(Force Exponent)为碰撞指数，剐度项的贡献因子，反映了材料的非线性程度，其推荐值：金属与金属材料为 1.5，橡胶材料为2；C(Damping)为阻尼系数，单位为 N·sec/mm，通常取刚度值的 0.1～1%；V 为接触点的法向相对速度； 是 的导数，单位为 mm/sec。

接触 刚度的表达式为 ：K- E[收藕日期] 2O12-12-14[作者简介] 蒋东霖(1979-)，男，硕士，高级工程师，主要从事机械系统设计和理论研究工作。

工程与试验其中，R- 为两碰撞物体在碰撞 时间单位为 s纵坐标为冲击力单位为 kN。

处的曲率半径； -百1-u 1-u，u 、u。分别是两物体的泊松比，E 、E 分别是两物体的弹性模量 。

3 仿真与试验根据碰撞力的模型，试验中采用冲击刀半径为2.Smm，采用 10×1O×55的标准试样。本文采用虚拟仿真平台，因为其具有 良好的界面和强大的计算功能。虚拟样机 中常用的积分器有 3种：GS-TIFF、WSTIFF与BDF。不同的刚性积分器，计算的效率和稳定性也有所不同。

GSTIFF是最常用的积分器，也是系统默认的积分器 。虚拟样机 中有 3种积分格式 ：SI3、SI2和SI1。SI3格式求解速度快，但精度太差；SI2格式可以避免 Jacobian矩阵的病态，而且考虑了约束方程，虽然求解时间变长，但求解精度高；SI1求解精度比SI2还高，但它计算量太大，-般不采用。积分误差决定了在求解动力学方程的过程中，某-步的预测值与校正值之间所能接受的差值。积分误差过大，计算容易进行，但最终结果会产生过大误差；积分误差过小，求解时间太长。本例模拟计算积分误差设为 0.001。

图 1为模拟计算的结果图，即冲击碰撞力随冲击时间的变化曲线，图中横坐标为时间，单位为 S，纵坐标为冲击力，单位为 N。

图 1 模拟计算结果 图图2为通过试验得出的结果图，即实际测得的冲击碰撞力随冲击时间的变化曲线，图中横坐标为· 72 ·图 2 试验结 果图如图 1、图 2所 示，图 l的最 大 冲击力 大约 为25kN，整个冲击过程的时间区间在0.78s0.788s之间。图 2中最大 冲击力大约为 23kN，整个冲击过程时间区间要比模拟计算的大，但是冲击力起始的时间很接近。模拟计算和试验数据的整个曲线趋势是-致的，冲击力产生的前期阶段，冲击力变化率都比较大，而冲击的后期阶段，冲击力变化率都趋于缓慢。

模拟计算与试验测得数据的差距产生的原因主要是因为，模拟是理想状态，认为试样材料是均匀的，试样材料本身的物理属性是理想的，排除了试样夹具座各部件、冲击刀、以及整个冲击摆的变形、振动。此外，模拟计算时采取的各个参数也并不能和实际的情况是完全-致的，但是整个模拟的结果与实际之间是 比较接近的。

4 结束语本文根据 Hertz接触理论，利用虚拟平台，模拟了冲击试验机摆锤冲击试样的整个过程，绘制了冲击力随时间变化的曲线，并和实际试验数据做了比较，分析了造成两者数据差别的原因，最终得出，模拟计算的结果和实际差别相差不大 ，可 以作为摆锤冲击试验机冲击力计算的参考依据，为设计摆锤式冲击试验机打下基矗