双盘转子扭转振动分析.docx

双盘转子扭转振动分析《动力学与控制学报》2014年第二期1盘面摩擦力与摩擦力矩考虑到发动机转子的结构非对称性以及固有的制造误差假设两转盘所在的位置将整个转轴分为三个不等长轴段如图1所示各轴段长度分别为l1l2和l3且l=l1+l2+l3;设转子的左、右轴承处等效支承刚度为kb质量为mb1与mb2;转子自转转速为Ω;转轴面二次极矩为Ip剪切模量为G杨氏模量为E转轴截面径向惯性主矩为I;转盘i(i=12)处的等效质量为mri偏心距为ei转盘i的半径为Ri;转动惯量为Ii封严内半径可近似为Rseal;以左端轴承处集中质量为中心建立笛卡尔坐标系o－xyzxy平面平行于转盘平面oz轴与转轴中线重合．由于温度的非均匀性使得高速运转的转子系统产生非均匀温度应力导致转盘产生由高温部分向低温部分的轴向窜动且窜动位移量与转、静件初始轴向间隙往往为同一个量级严重情况下会产生转静子盘面挤压．在计算盘面摩擦力与力矩时不考虑扭转振动、弯曲振动与拉压振动之间的耦合关系．在前述假设下轴向碰摩近似为两平面的接触挤压轴向单位面积接触等效刚度ka为一个与材料属性有关的常量摩擦系数为μ．由文献［12－13］分析可知受温度环境的影响转子主轴在热膨胀作用下将产生窜动位移使得盘面将朝向低温段方向窜动．盘处产生的挤压量δi≥0时发生轴向环带面接触其中阴影部分S={(uφ)|0≤φ≤2πD2≤u≤D1}为接触面或挤压面．在图2所示的转盘随体极坐标架Oruφ下D2与D1为对应的内外边界方程微元ds的速度包含三部分:随转盘质心的刚体平动速度vtran、随转盘自转的速度vroto和盘面扭转速度vtors微元总速度vds在惯性坐标系o－xyz下的形式为:2双盘转子扭转振动动力学建模为研究转轴扭转振动将转盘与轴承处的相对转角变化θ1和θ2取为扭转模型的广义坐标如图3所示系统扭转振动的动能:3双盘转子弯扭耦合振动微分方程航空发动机在通常的运行情况下转轴将会发生扭转与弯曲组合变形．本节将建立转子主轴在该组合变形下双盘转子系统的动力学模型并讨论相关的动力学特性．取(θ1θ2xr1xr2yr1yr2xb1xb2yb1yb2)为广义坐标采用柔度影响系数法与D''''Ale-mber原理建立系统的运动微分方程．转轴在oxz与oyz面内的弯曲如图4所示:考虑oxz面内的弯曲变形(oyz内的情况类同)设转轴上转盘形心在ox方向上的弯曲柔度矩阵为H=(hij)H－1即为转轴弯曲刚度矩阵．应用弯曲变形理论可计算出H的分量为:4面摩擦力矩对扭振的影响较大的温度梯度所引发的转子轴向窜动可能导致转静件盘面接触以及挤压由此产生的盘面摩擦力是改变整个系统动力学行为的重要因素．当转盘i处的挤压量δi≥0时(δi=0仅发生盘面接触摩擦力为零)在该转盘处将发生转静盘面接触挤压与摩擦．同时产生作用在转盘形心并与轴向方向平行的周期激励扭矩引起转轴的扭转振动．为了考察不同工作参数对双盘转子系统扭转振动的影响采用四阶Runge－Kutta积分法对转子动力学无量纲运动微分方程进行数值求解．选取系统参数如表1所示这里针对不同工况分9组(见表2)分别进行计算．定义响应s(t)的幅值为:从转盘处的振动响应(图5、图6)图上可以看出扭转响应的量级是毫秒级而弯曲振动响应的幅值是丝米级相比而言扭转振动响应的影响很小．图7给出了给定挤压量条件下由于轴向窜动而引起的摩擦力矩．从表3－4可以看出:比较数组1与数组2增加转盘2处的挤压量将增加2处摩擦扭矩进而增加了该处的扭振幅值而转盘1处的扭振影响较小;比较数组2、3当转盘1处的挤压量较大时转盘1与2处扭振的幅值与平衡位置均随之增加且对转盘1处的影响较明显;比较数组2、4与6当转盘处的偏心距增加时偏心激振力加将影响双盘转子系统弯曲振动间接地改变扭矩激励但对转盘扭转振动的影响较小．比较数组2、5与7当增加面接触刚度或转静盘面处挤压量值较大时摩擦扭矩随之加大将造成较为严重的扭转振动;从数组8与9可以看出当仅仅一个转盘发生摩擦另一处扭矩为零时扭转振动的幅度较小且无摩擦扭矩的转盘将产生被动扭振扭振幅值较存在摩擦的转盘约小一个量级故其影响可以忽略．因此扭转振动的幅值随周期摩擦激励的增加而增加．比较图8～11可以看出扭转振动的幅值与均值与其相对应的摩擦扭矩具有同步变化规律扭振的平均位置受该处摩擦扭矩的均值制约即当忽略掉轴承处的扭矩而将两端视为扭转固定端时双转盘处受扭矩均值的作用将对轴系产生静扭转量直观上造成双盘的扭转搓动效果．5结论温度非均匀分布是航空发动机转子的典型运行环境．转静部件的热膨胀系数不一致以及非均匀温度场导致转盘轴向位移以及转静件轴向间隙的改变是产生轴向碰摩的重要因素之一