车辆智能革速新式齿轮程序的情态化摹拟

　　系统的状态方程状态方程是借以表征系统内部状态随时间变化的物理变量，根据系统的键合图模型，以惯性元的广义动量p和容性元的广义位移q为状态变量，各单元的势和流分别用e和f表示，对各键进行统一的编号120，其中变量p和q的下标代表部件编号，变量e和f的下标代表所在键的编号。系统各储能元件都具有积分因果关系，设整个系统的状态变量为：X=[q0，p4，qH4，p1，q1，p1c，qH2，p2]T(1)根据系统因果关系和功率流方向得到该系统键合图仿真状态空间方程：q0=-1I4 Sfp4=1C0q0-1CH4qH4qH4=1I4(2)式中：m1=11 p;m2=p1 p;m3=11 p行星齿轮特性参数：p=z2z1=10848，pc=z3z1c=84364系统的动态仿真分析根据行星齿轮传动系统参数和材料力学特性等，由材料力学[7]推荐公式来确定状态方程中惯性元I和容性元C等参数。系统的基本参数见。

　　系统的基本参数类别太阳轮Z1行星轮Z6内齿圈Z2太阳轮Z1c行星轮Z5内齿圈Z3齿数4830108362484齿宽/mm5040模数/mm3压力角/(b)20由于各部件相应的参数差别较大，所以对状态方程中系数作相应无量纲处理，采用四阶龙格-库塔法对系统状态方程式(2)进行微分方程求解。为得到系统的真实输出(T，X)，对得到的广义动量p除以相应I元件的惯量就能得到转速，对得到的广义位移q除以相应C元件的柔度就能得到转矩T[6].得到系统行星架的变形曲线和输出仿真结果。

　　结束语由仿真曲线可以看出：①由于行星架2的刚度比行星架4的小，变形较大，但是变形过程较4缓慢，而行星架4与输入轴作为一体，受到输入力矩、内齿圈和太阳轮的同时作用，受力复杂，变形产生波动。

　　②由公式X=p/I知，系统是由2输出，其转速随2的广义动量增加而增加。可中输出转速曲线能正确反映实际情况。而输出转矩受行星架4变形的影响，而产生波动。③键合图作为一种分析方法，简便易行、准确可靠，在系统动力学分析中具有广阔的应用前景。