接触刚度的计算

1、istep(time,0,0d,0.68,-12000d)+step(time,0.68,0d,1.77,0d)+step(time,1.77,0d,2.45,12000d)3.2.3定义齿轮啮合的接触碰撞力为了保证仿真分析的真实性，齿轮之间的啮合运动关系没有被定义成理想化的几何约束关系，而是被定义为基于接触碰撞的力约束关系，即齿轮之间只能通过接触碰撞力(法向)和摩擦力(切向)相互约束，而不存在其他的约束关系。在ADAMS中有两种接触碰撞的计算模型，一种是基于Hertz理论的Impact函数模型，一中是基于恢复系数(Coefficientofrestitution)的泊松(POISSON)模型

2、。两种力模型都来自于法向接触约束的惩罚函数。ADAMS/C+Solver使用惩罚因子来转换所有的接触约束。采用Impact函数来计算各啮合齿轮轮齿之间的接触碰撞力。Impact函数模型将实际中物体的碰撞过程等效为基于穿透深度的非线性弹簧阻尼模型，其计算表达式为：xx1其中K接触刚度系数；x1位移开关量，用于确定单侧碰撞是否起作用;x接触物体之间的实测位移变量；d阻尼达到最大时两接触物体的穿透深度；Cmax最大接触阻尼；x穿透速度；n非线性弹簧力指数。当xx时，两物体不发生接触，接触力为0,当xx时，两物体接触，接11触力大小与接触刚度系数、非线性指数、阻尼系数以及两物体距离的改变量即穿透量有关

3、。由以上公式可知，Impact接触力包括两个部分：弹性分量K(x-x)n，相当于一个非线性弹簧；阻尼分量step(x,x-d,C,x,0)x，其方向与运动方向相反，为了避1max1免阻尼分量突变而使得函数变得不连续，采用了阶跃函数step()来定义阻尼，step()函数是利用三次多项式逼近海赛(Heacisde)阶跃函数，具有连续的一阶导数，但在起始点处二阶导数不连续。在ADAMS中的表达形式为：step(x,x,h,x,h)0011方0=h+(h一h)l(xx)/(x一x)1一2(x一x)/(x01001001hI1xxxxx13-2)其中，x为自变量，当x小于x时，因变量的值为初始值h,当

4、x大于x时，001因变量的值为终止值h；当x在初始值和终止值之间变化时，因变量根据一定规1律光滑过渡，避免出现数值过渡突变、微分值不连续。由此定义的碰撞阻尼系数，如图3-5所示。（1）弹性分量计算物体接触刚度系数与物体的材料属性和接触表面的几何形状有关，在此根据文献1提供的接触刚度计算式来计算各齿轮啮合的齿廓面接触刚度，计算式为43兀(h+h)12RR11/21_R+R123-3)其中，R，R相啮合的两齿廓面在啮合点处的曲率半径。对于渐开线12齿轮，其工作过程中啮合点在齿廓面上的位置是不断变化的，其曲率半径也是不断变化的，主动轮齿廓面上啮合点处的曲率半径由小变大，从动轮尺阔面上啮合点处的曲率半

5、径由大变小，因此以两轮齿在节点处啮合作为计算点，则R二mZisina；i二1,2，m为模数，z为啮合齿轮的齿i2数，a为节圆压力角，对于标准啮合，节圆压力角等于分度圆压力角20o。1一u2h，h材料参数，定义为：h二i；i二1,2，u为泊松比，E为12i兀E弹性模量非线性指数n，根据Hertz理论，一般取1.5较合适2。阻尼分量计算采用文献2给出的非线性阻尼模型来计算轮齿啮合接触阻尼系数，该阻尼模型认为物体表面接触碰撞过程中的能量损失是由接触阻尼引起的，在基于等效能量损失的基础上给出接触阻尼的计算式：3-4)C=3K(1一”8n4U其中，K接触刚度；e弹性恢复系数，一般定义为碰撞前法向速度差值

6、与碰撞后法向速度差值的比值。跟物体的材料、碰撞表面曲率半径、碰撞速度、以及润滑介质的粘度有关，一般通过实验测定。8穿透深度，对应ADAMS取最大阻尼系数时的穿透深度，在此取8-0.1mm；n非线性指数，一般取1.52。U碰撞速度，以相啮合的两个齿轮在节点处的线速度的差值代替。在齿廓面间的动态碰撞力的作用下，相啮合的两个齿轮在节点处的线速度是不等的，且随时间变化，即碰撞速度也不是定值在阻尼分量中，弹性恢复系数和碰撞速度与实际的工况有关，其具体的取值在各仿真工况下确定。摩擦力计算接触体之间的摩擦力采用库仑摩擦模型，考虑静摩擦和动摩擦，有润滑时，取静摩擦系数为0.08，动摩擦系数为0.05，无润滑时

7、，取静摩擦系数为0.15，动摩擦系数为0.1。在ADAMS中摩擦力采用下面的函数表达式计算：F=-N*step(v,-Vs,-1,Vs,1)*step(ABS(v),Vs,Cst,Vtr,Cdy)(3-5)其中：N法向力；v表面相对滑移速度；Vs最大静摩擦对应的相对滑移速度；Cst静摩擦系数；Vtr动摩擦对应的相对滑移速度；Cdy动摩擦系数；摩擦系数与相对滑移速度的关系，如图3-5所示。图3-5摩擦系数与相对滑移速度的关系曲线4发射装置动力学仿真4.1ADAMS动力学分析算法及求解器的选择ADAMS对动力学微分方程，根据机械系统特性，选择不同的积分算法：对刚性系统，采用变系数的BDF(Back

8、wardsDifferentiationFormulation)刚性积分程序，它是自动变阶、变步长的预估校正法(PECE，Predict-Evaluate-Correct-Evaluate),并分别有Index3、SI2、SI1三种积分格式，在积分的每一步采用了修正的Newton-Raphson迭代算法，对应的求解器有GSTIFF、WSTIFF、DSTIFF、CONSTANT_BDF；对高频系统，采用坐标分块法将微分一代数(DAE)方程简化为常微分(ODE)方程，分别利用ABAM(Adams-BashforthandAdams-Moulton)方法和龙格一库塔(RKF45)方法求解，对应的求解

9、器有ABAM求解器和RKF45求解器。4.11BDF刚性积分法步骤预估阶段用Gear预估一校正算法可以有效地求解微分一代数方程。首先，根据当前时刻的系统状态变量，用泰勒级数预估下一时刻系统的状态矢量值：(4-1)Qy丫1d2yy=y+nh+nh2+n+1yQt2!Qt2其中，时间步长h=t-1。这种预估算法得到的新时刻的系统状态矢量值n+1n通常不准确，可以用Geark+1阶积分求解程序(或其他向后差分积分程序)来校正。4-2)y二-hPy+ayn+10n+1in-i+1i=1其中，y为y（t）在t=t时的近似值；P和a为Gear积分程序系统值。n+1n+10i上式经过整理，可表示为：yn+1

10、-1hp0yn+1-ayin-i+1i=14-3)（2）校正阶段求解系统方程G,如G（y,y,t）=0,则方程成立，此时的y为方程的解，否则继续。求解Newton-Raphson线性方程，得到Ay，以更新y,使系统方程G更接近于0。JAy=G（y,y,t），其中J为系统的雅可比矩阵。n+1采用Neeton-Raphson迭代，更新y:yk+1=yk+Ayk重复以上步骤直到Ay足够小。（3）误差控制阶段预估计积分误差并与误差精度比较，如积分误差过大，则舍弃此步。计算优化的步长h和阶数n。如达到仿真结束时间，则停止，否则t=t+At，重新进入第一步。4.1.2坐标缩减的微分方程求解过程步骤（1）坐

11、标分离将坐标的约束方程进行矩阵的满秩分解，可将系统的广义坐标列阵訂分解成独立坐标列阵Qh口非独立坐标列阵V），即q=$iqd丄。（3）预估用Adams-Bashfoh显示公式，根据独立坐标前几个小时步长的值，预估t时刻的独立坐标值ii，P表示预估值。（4）校正用Adams-Moulton隐示公式对上面的预估值，根据给定的收敛误差限进行校正，以得到独立坐标的校正值*，C表示校正值。（5）确定相关坐标确定独立坐标的校正值后，可由相应公式计算出非独立坐标和其他系统状态变量值。（6）积分误差控制与上面预估校正算法积分误差控制过程相同，如果预估值与校正的差值小于给定的积分误差限，接受该解，进行下一时刻的求解。否则减小积分