动刚度分析的后处理方法

  在《NVH一本通》中讲解了四种不同的处理方法，分别是平均法、指数法、面积法和三分之一倍频程法等，其详细流程见书籍。本文主要讲解一下面积法，其余方法参阅书籍内容。

其中面积法是通过IPI曲线所包围的面积计算得到，其公式17-24推导如下。IPI （ Input Point Inertance源点导纳分析）是指在一定的频率范围内，通过在加载点施加单位力作为输入激励，同时将该点作为响应点， 测得该点在该频率范围内的加速度作为输出响应， 用于考察该点的局部动刚度。  

源点加速度导纳公式为：

由上式可知，面积法所得到的动刚度是通过响应曲线与横坐标频率围成的面积进行计算得到，如下所示的叉的区域，如图17.37所示。

图17.37 面积法动刚度结果面积图示  

采用面积法进行动刚度计算，需要设定一个目标参考值，即如参考刚度值为1200N/mm，通过面积公式17-24可得需要计算的动刚度公式17-25。

其计算公式为：

Kx={absarea((p4w1c4.x),(p4w1c4.y))*1200./absarea((p4w1c1.x),(p4w1c1.y))}

   1200是目标动刚度值，单位N/mmm;absarea表示求面积

   ((p4w1c4.x),(p4w1c4.y))是目标动刚度曲线的x轴和Y轴；

   ((p4w1c1.x),(p4w1c1.y))是激励点1001X方向的加速度响应曲线的x轴和Y轴；

可得1001点三个方向的动刚度值分别为811.4N/mm、1944N/mm、660.7N/mm。

若要计算50-200Hz频率范围内的动刚度可以按如下操作：

Kx_f50_200= {absarea(50:200,lininterp(p4w1c4.x,p4w1c4.y,50:200))*1200./absarea(50:200,lininterp(p4w1c1.x,p4w1c1.y,50:200))}

通过计算可得1001三个方向的动刚度分别为864.6N/mm、2755.9N/mm、662.1N/mm。

     通过各种方法得到的动刚度有所差异，但在实际工程中，在性能目标定义时应该统一同一种方法即可，包括CAE分析、竞品测试或实车测试，这样的分析优化等才具有一定的实际指导意义。  

参考文献：《汽车NVH一本通，建模、优化及应用》2023.10机械工业出版社