ADAMS软件中IMPACT函数参数确定.docxVIP

3.3.3 ADAMS 软件中 IMPACT 函数参数的确定 (1) 碰撞函数中的非线性弹簧参数的确定 在使用 ADAMS 碰撞函数中，确定碰撞力大小 Fn 主要由等效刚度 k 以及幂指 数 q 两个参数确定。为了合理的确定等效刚度以及幂指数，通常使用 Hertz 弹性 碰撞模型来计算出正确的 k 与 q 。Hertz 弹性碰撞模型如下： v A Fn A A B v B B F （ a）碰撞前 （ b）碰撞 （ c）碰撞后 图 3-3 弹性碰撞示意图 如图 3-3 所示，设某一半径为 R的圆球体 A 以速度 v 飞向平面 B，则该过程分为以下三个主要过程： 当碰撞发生时，物体 A 与物体 B 之间有力的作用，其中在碰撞时因为碰撞 产生了弹性力 Fn 同时物体 A 也发生了一定的变形 ， Hertz 模型指出 Fn 与 满足 如下关系： 2 1 3 9Fn (3-14) 16ER2 其中： 1 R 1R 1 R （ R1,R2 分别为碰撞点处两物体的曲率半径） 1 2 1 (1 12 ) (1 22) E E1 E2（E1,E2 分别为两物体的材料弹性模量 , 1 , 分别为两物体的材料泊松比）由式 (3-14)我们可知 Fn 2 * 16RE 2 3 9 故而可以确定 16RE2 k 9 (2) 碰撞函数中的阻尼参数 cmax 的确定 25/7 式(3-12) 可以写成 Fn k n D (3-16) 其中 D 为该式中的阻尼参数， 为两物体接触时的相对移动速度。 Hund 与 Grossley[83] 曾提出确定 D 的方法： D  n (3-17) 在式 (3-17) 中 被称为滞后阻尼因子。 在碰撞非线性弹簧阻尼模型中， 碰撞过程 会导致能量的损失， 这是由碰撞的阻尼项造成的。 图 3-4 表示在碰撞过程中碰撞力的变化： 图 3-4 碰撞力的变化 在该图中 t( ) 表示开始碰撞前的那一时刻， t (m ) 表示碰撞达到最大穿透量的那一时 刻，t ( ) 表示开始碰撞结束后两物体分离的那一时刻。 在式 (3-17) 中阻尼系数 D 或 滞后阻尼因子 需要确定。应用经典的冲量定理与能量守恒原理可以确定以上两 个变量。 根据能量守恒定理可知在碰撞前与碰撞后两物体的动能损失可以用恢复系 数 e 与两物体间的相对速度 ( )来表示 其中 () Vi ( ) Vj( ) ( Vi ( ) 表示物体 i 碰撞接触前的速度， Vj( ) 表示物体 j 碰 撞接触前的速度 ) 26/7 T 1 (eff ) ( )2 2 (3-18) m (1 e ) 2 m(eff ) mi mj (3-19) mj mi 在碰撞过程中的能量损失也可以对碰撞力进行积分来获得 [84] ： m 2 m eff( ) 3 T D d n d 2 n d (3-20) 0 3 K 式 (3-20) 中 K 与 n 的取值与非线性弹簧项中对应的取值相同 . 由式 (3-20) 可知滞后因子 可以算得 : 3K (1 e2 ) (3-21) 4 ( ) 即在式 (3-12) 中 cmax D n 由以上可知阻尼系数与两碰撞物体本身材料属性有关 ( 杨氏模量 , 泊松比 ) 。同时，也与两碰撞物体的外形轮廓有关。要得到恢复系数 e 的大小 , 该值一般需由实验测定。 碰撞函数中 dmax 的确定 在 IMPACT函数中 , 参数 dmax 的作用在于两物体接触后当两物体的穿刺深度 dmax 时，令非线性阻尼系统中的阻尼大小的取值为 D , 当 0 dmax 时, 其非线性弹簧阻尼系统中的阻尼大小由 STEP 函数决定 , 其大小与穿刺深度的关系变化见图 3-5 。 27/7 图 3-5 由碰撞动力学模型可知当两物体一产生接触时就产生阻尼 , 并且阻尼系数的 大小在一次完整的碰撞过程中是不变的 , 因此 dmax 的取值应该越小越好 . 同时考虑到 ADAMS 中的数值收敛性,一般可采用 ADAMS 中推荐的取值 dmax 0.01mm 碰撞函数中摩擦系数的确定 由式 (3-13)可知摩擦力由 Coulomb 模型确定 ,即摩擦力大小取决于作用于物体之上的正压力大小与摩擦系数的乘积 . 正压力的大小可以由 (3-12)式得出； 摩擦系数分为动摩擦系数 d 与静摩擦系数 s ，当两接触物体处于静止状 态时，二者之间的摩擦系数即取为静摩擦系数 ,当两接触物体的相对运动速度较 大处于运动状态时二者之间的摩擦系数即取为动摩擦系数 .一般而言动摩擦系数 和静摩擦系数与两接触物体的材料各属性以及表面的粗糙程度等相关因素有关 , 需要通过实验获得数据，以下表是常用的不同材料相接触的摩擦系数： 28/7 表 3-1