工程车辆传动系统扭转振动特性研究与分析

1工程车辆扭转振动动力学模型的建立

工程车辆传动系统一端通过离合器与发动机相连，输出端通过轮胎与工程车辆平动质量相连，组成了一个多质量的弹性扭转振动系统。在计算整个系统的固有频率和振型时，通常可忽略系统的阻尼，将整个传动系统看成是由多个刚性圆盘通过弹性轴连接的无阻尼振动系统。现在某型装备四缸柴油机的中型装载机传动系统为例，其扭转振动力学模型如图1-1所示。

1.1 当量转动惯量的计算

当量转动惯量J是指将传动系统中与发动机曲轴不同转速旋转的零部件的转动惯量换算成与曲轴同转速旋转下的转动惯量，这种换算方法的原理是能量守恒。设传动轴的转动惯量为J,实际转速为?曲轴转速为?0，则将传动轴换算成曲轴转速?0的当量转动惯量为

112J?2?Jd?022??Jd?J????0式中，ig为变速器的传动比。 1.2当量扭转刚度的计算

设两圆盘之间弹性轴的当量扭转刚度为Kd，则可以根据弹性变形量守恒的原理将系统中的时间扭转刚度K换算过来。现以后桥半轴为例，相应的当量扭转刚度为

?1??????J??i???g?22

?1?? Kd?K??ii??g0?2式中，i0为主减速器的传动比。

2传动系统扭转动力学方程

根据图1-1所示的简化的传动系统模型，可建立系统动力学方程组为

J1?1?K（?T11?1-?2）J2?2-K（?K（?T21?1-?2）2?2-?3）J3?3-K2(?2??3)?K3(?3??4)?T3??J10?10-K（（?10-?11）?09?9??10)?K10J11?11?K10（?10-?11）?0方程组（1）中，?1-?11分别为对应质量的扭转角位移；T1-T4分别为发动机1-4缸的有效输出转矩。

为了简单起见，可以将（1）改为矩阵形式的动力学方程一般式，即

。。。。。。。。。。 (1)

J??C??K??T

。。。0?J1?J2?式中，当量转动惯量矩阵J????0J10???阻尼矩阵C=[0]；刚度矩阵； 圆盘的角位移矩阵????1???? ??J11???2?3?40??11?T。

T??T1T2T3T4T0?0?

一般以发动机振动激励为系统输入矩阵，则

2.1扭转系统固有特性的分析

这里的固有特性是指固有频率和主振型，多自由度系统的固有频率和主振型可以根据系统的无阻尼自由振动方程得到，即

J??K??0 （2）

假设方程的解为

。。??Aeti?n （3）

式中，A为系统自由振动时的振幅列向量，A??Am1Am2Am3?Am11?。

T将式（3）及其二阶导数代入方程式（2），并消去

2n??Aeti?n，得到主振型方程为

?K??J?A?0

2令H?K??nJ，则有

HA=0

式中，H为系统的特征矩阵。 根据线性代数知识可知，主振型方程存在非零解的条件是系统的特征矩阵的行列式必须等于0，即

2H?K??nJ?0 (4)

式（4）称为系统的特征方程。当系统微分方程组较多时，解上述方程的工作量是相当大的，

调用MATLAB中的eig（）函数来求解系统特征方程，可以非常方便地计算出系统的固有频率，其对应的特征矢量就是该固有频率所对应的振型。 此外，还可以根据计算出的振型画出振型图。下面直接给出该型装载机传动系统各部分的转动惯量和各轴段的扭转刚度参数值，如表1-1和表1-2所示。

表1-1 某型装载机传动系统转动惯量参数值 零部件名称 发动机第一曲轴连杆机构 发动机第二曲轴连杆机构 发动机第三曲轴连杆机构 发动机第四曲轴连杆机构 发动机飞轮及离合器总成 变速器第一轴、中间轴总成及倒档齿轮 变速器第二轴总成、取力器总成 主减速器主动锥齿轮及差速器总成 主减速器从动锥齿轮及差速器总成 两个半轴及两个驱动车轮 装载机平动质量的当量转动惯量 零部件名称 发动机第一、第二杆轴颈间轴段 发动机第二、第三杆轴颈间轴段 发动机第三、第四杆轴颈间轴段 发动机第四轴颈杆中部与飞轮之间轴段 变速器第一轴 变速器第二轴及中间轴 传动轴总成轴段 主减速器主动锥齿轮 后桥两个半轴 转动惯量/（kg.m2） J1=0.001964 J2=0.00191 J3=0.00191 J4=0.0019241 J5=0.1055 J6=0.0030074 J7=0.031771 J8=0.006066 J9=0.13724 J10=6.2947 J11=159.47 扭转刚度/（N.m /rad） K1=8.113e5 K2=8.113 e5 K3=8.113 e5 K4=8.2443 e5 K5=1.1962 e6 K6=1.8904 e5 K7=1.7945 e5 K8=1.056 e7 K9=1.1529 e6 表1-2 装载机传动系统轴段扭转刚度参数值 两个驱动车轮 K10=2.0735 e5 对于n个自由度的系统，存在n个固有频率和n个相应的主振型，第i阶固有频率对应第i阶主振型Ai。它们总是成对出现，并描述系统的一个独立特性。需要注意的是，将特征值代回至系统微分方程组只能得到n个未知量Ai之间的比值。如果Ai是方程的解，则Ai乘以任何非零的常数都是方程的解。因此，可以说，主振型的形态是确定的，但是振幅是不确定的；换句话说，主振型只是确定系统按照某一阶固有频率自由振动时的各个坐位移的比值，而振幅的数量则可以是任意的。 因此，在计算机主振型时，无需求出具体初始条件下系统某阶主振动时各个坐标幅值得具体绝对数值，而只是一般性地描述系统某阶主振型的形式，任意规定其中某一坐标的幅值。 2.2 仿真计算

经过以上的理论分析，下面进行程序的编制，可自行编写函数eigfre-niuzheng（）来求解系统特征方程，以得到特征值和特征向量。该函数为求解各种振动系统特征方程的通用函数，本书还将利用本函数进行悬架双质量系统的固有特性分析。 具体程序如下：

function[w,A]=eigfre_niuzheng(J,K) w=[ ];

H=inv(J)*K; [A,d]=eig(H); n=size(H); for i =1:n w(i)=d(i,i);

A(:,i)=A(:,i)/A(1,i); end

for j =1:n-1; for i=j+1:n; if w(j)>w(i); t=w(i); w(j)=w(i); w(j)=t; q=A(:,i); A(:,i)=A(:,j); A(:,j)=q; end end end m=5; for i=1:m subplot(m,1,i) plot(A(:,i)) end

对于当量转动惯量和当量扭转刚度的计算，特别要注意的是：在离合器完全结合的情况下，通常可以认为离合器、机械变速器第一轴转速等于发动机曲轴的转速?0；而变速器第二轴、

传动轴以及主减速器主动锥齿轮为同一转速?0ig；主减速器从动锥齿轮、差速器总成、两根

驱动半轴以及两个驱动车轮为同一转速?0/(igi0)。为了分段进行当量转动惯量和当量扭转刚度的计算，现自行编制MATLAB函数文件niuzheng_para。本函数输入参数为变速器的档位gi，输出参数为传动系统的当量转动惯量矩阵j和当量扭转刚度矩阵k。 具体程序如下：

function[j,k]=niuzheng_para(gi) ig=[5.78 4.57 2.91 1.72]; i0=7.28;

jj=[0.001964 0.00191 0.00191 0.0019241 0.1055 0.0030074 0.031771 0.006066 0.13724 6.2947 159.47];

kk=[8.113e5 8.113e5 8.113e5 8.2443e5 1.19962e6 1.8904e5 1.7945e5 1.056e7 1.1529e6 2.0735e5] for i =1:11; if i>4&i<9

j(i)=jj(i)./(ig(gi))^2; elseif i>8

j(i)=jj(i)./(ig(gi)*i0)^2; else

j(i)=jj(i); end end

for i=1:10; if i>3&i<8

k(i)=kk(i)./(ig(gi))^2; elseif i>7

k(i)=kk(i)./(ig(gi)*i0)^2; else

k(i)=kk(i); end end

以某中型装载车辆为例，该工程车辆所用发动机是四缸四冲程柴油机，根据四缸柴油机的特定分析，其扭转激励主谐次为k=1,2,3…该发动机正常怠速约为700r/min，根据工况特性考虑发动机常用工作转速ne=1200-2500r/min。 工作转速范围会有相对应的基本频率范围，即

fe??nek30

?nekfe??(251.327?523.599)Hz,k?230工程上通常认为低频振动对动力传动系统部件影响最大，因此在这里只考虑低频特性。综合这些因素可以得出，要想避开共振区域，动力传动系统的固有频率就不能落在这些频率范围内或其附近。下面分析的主要目的就是校验动力传动系统的固有频率是否落在这些频率范围内或其附近。下面分析的主要目的就是校验动力传动系统的固有频率是否落在这些范围内。

?(125.664?261.799)Hz,k?1

百度搜索“77cn”或“免费范文网”即可找到本站免费阅读全部范文。收藏本站方便下次阅读，免费范文网，提供经典小说教育文库工程车辆传动系统扭转振动特性研究与分析在线全文阅读。