一种基于复系数法的磁悬浮扁平转子系统章动频率解析方法技术方案

本发明专利技术涉及一种基于复系数法的磁悬浮扁平转子系统章动频率解析方法，根据欧拉运动学方程，通过对磁悬浮扁平转子系统转动动力学建模分析，建立控制系统参数及转速与转动运动间的关系，利用复系数法，将二自由度转动动力学方程表示为单自由度复变量特征方程，通过变量重构建立阻尼振动频率特征方程，根据章动模态与转子转速同向的特点，获得章动频率与控制系统参数和转速间的解析关系。该方法计算简单，容易实现，可较精确的解算分散PID控制下章动模态频率，提高了章动模态的控制针对性和效率。

A nutation frequency analysis method for magnetic levitation flat rotor system based on complex coefficient method

The invention relates to a magnetic flat rotor system nutation frequency coefficient method based on complex analytic method, based on the Euler kinematic equation, through the analysis of rotational dynamics modeling of the magnetic suspension flat rotor system, establish the control parameters and the relationship between speed and rotational motion, using the complex coefficient method, the two degree of freedom rotational dynamics equation expressed as a single freedom the degree of complex variable characteristic equation, establish characteristic damping vibration frequency equation by variable reconstruction, according to the characteristics of nutation modes and rotor speed in the same direction, obtain the analytical relationship between nutation frequency and control system parameters and the speed of the. The method is simple and easy to implement, and can accurately calculate nutation modal frequencies under decentralized PID control, and improve the pertinence and efficiency of nutation mode control.

【技术实现步骤摘要】
一种基于复系数法的磁悬浮扁平转子系统章动频率解析方法
本专利技术涉及磁悬浮扁平转子系统控制的
，具体涉及一种基于复系数法的磁悬浮扁平转子系统章动频率解析方法，可用于分散PID控制条件下磁悬浮扁平转子系统的章动频率解析计算。
技术介绍
磁悬浮轴承具有无接触、无摩擦、无需润滑、使用寿命长等诸多优点，随着磁悬浮控制技术的日趋成熟，其广泛应用于磁悬浮分子泵、磁悬浮储能飞轮、磁悬浮控制力矩陀螺、磁悬浮高能量密度电机等机构。磁悬浮轴承控制系统通常采用分散PID控制，但是随着磁悬浮转子转速的升高，积分环节和滤波环节的引入使系统的频率特性产生相位滞后，由于强陀螺耦合效应，涡动(进动和章动)模态得不到足够的阻尼而容易造成高速转子的进动或章动失稳，对于数字PID控制，系统的采样频率和运算速度有限造成的控制延时，也是高速旋转时陀螺效应导致系统失稳的主要原因。因此，对转子涡动模态的控制是磁悬浮转子系统高速稳定运行的必要条件。研究发现，磁悬浮高速转子的涡动具有一定的频率特性。随着转速的升高，转子系统的两个涡动模态出现频率分叉，章动频率随转速同步上升，逐渐趋向于转动惯量比(极转动惯量Jz与赤道转动惯量Je之比)与转速的乘积，进动频率随转速上升而下降，逐渐趋向于零。高速时，章动模态成为影响系统稳定性的主导因素，章动模态的稳定控制决定了磁悬浮转子系统所能达到的最高转速。目前，国内外研究主要是针对章动模态的定性分析，而对于具有强陀螺效应的高速转子系统在闭环控制条件下的章动频率定量分析缺乏深入的研究。
技术实现思路
本专利技术要解决的技术问题是：克服现有技术的不足，利用复系数法，提供一种基于复系数法的磁悬浮扁平转子系统章动频率解算方法，该方法可建立章动频率与控制系统参数和转速间的关系，降低了高速磁悬浮转子由于章动模态突发造成失稳的风险，提高了章动模态抑制的针对性和效率。本专利技术解决上述技术问题采用的技术方案是：一种基于复系数法的磁悬浮扁平转子系统章动频率解析方法，包括以下步骤：(1)根据欧拉运动学方程，建立高速磁悬浮扁平转子转动动力学模型：式中Jx为x轴方向转动惯量，Jy为y轴方向转动惯量，Jz为z轴方向转动惯量，Ω为转速，α为转子绕x轴的偏转角度，β为转子绕y轴的偏转角度，lm为磁轴承力到质心的距离，fax、fbx、fay、fby为径向四通道ax、ay、bx、by方向的磁轴承力。(2)对于分散PID磁轴承控制系统，各控制通道等效为控制参数相同的单自由度控制系统，根据磁悬浮转子系统前向控制通道，可将控制电流i与转子位移qs间的频域关系表示为：i(s)＝-φi(s)qs式中qs为传感器位置处检测到的转子径向位移，φi(s)是由控制器Gk、功放Ga、传感器放大倍数ks等环节构成的控制增益，其频域关系可表示为φi(s)＝ksGk(s)Ga(s)。由于磁轴承力可在参考工作点近似线性化为f＝kii+khqh，因此，可将径向四通道ax、ay、bx、by方向的磁轴承力由控制系统参数表示如下：式中ki为电流刚度，kh为位移刚度，ls为传感器探头到转子质心的距离，表示卷积算子。(3)回转体结构转子，其x、y轴转动惯量相等，即Jx＝Jy＝Je，令Je为赤道转动惯量，可得系统时域磁悬浮转子转动动力学方程为：磁悬浮转子转动动力学方程可表示为仅与控制系统参数、转速和转子转动角度相关的函数。(4)当转速不为零时，转子绕x、y轴的转动运动相互耦合。采用复系数法，对转动变量重构，将二自由度转动动力学方程表示为单自由度复变量形式，原MIMO系统转变为SISO系统，将原转动动力学方程化简为：(5)对变量重构后的转动动力学方程进行拉普拉斯变换，零初始条件下令s＝jω，则可得求解阻尼振动频率的特征方程为：式中ω为磁悬浮转子涡动角速度，|φi(jΩ)|为控制算子。(6)根据阻尼振动频率的特征方程，解得分散PID控制下磁悬浮扁平转子正反两个涡动模态的阻尼振动频率解析解：其中ω+＞0，ω_＜0。根据章动与转速方向相同的特性，可知章动频率ωn＝ω+。因此，磁悬浮扁平转子系统章动频率解析解可由控制系统参数及转速表示为：本专利技术的原理是：磁悬浮轴承控制系统中积分环节和滤波环节的引入使系统的频率特性产生相位滞后，由于强陀螺耦合效应，涡动模态得不到足够的阻尼而容易造成高速转子涡动失稳，在高速状态下，章动模态成为影响磁悬浮转子系统稳定性的主导因素。对于章动模态频率精确解析计算的作用，就是对于章动模态特性有预见性的分析，从而有效的控制章动模态。通过对磁悬浮转子系统转动动力学建模分析，建立控制系统参数及转速与转动运动间的关系，利用复系数法，通过变量重构建立阻尼振动频率特征方程，根据章动模态与转速同向的特点，获得控制系统参数与转速和章动频率间的解析关系。本专利技术与现有技术相比的优点在于：(1)本专利技术首次提出了一种基于复系数法的磁悬浮扁平转子系统章动频率解算方法，该方法计算简单，易于实现，可适用于大多数磁悬浮高速转子的章动模态分析控制场合，提高了章动模态控制方法的针对性及效率。(2)本专利技术所采用的方法无需借助任何设备，利用复系数法，结合磁悬浮转子控制系统参数及转速信号就可实现对章动频率的精确解析计算。附图说明图1为磁悬浮转子结构控制示意图，其中，1为传感器放大倍数，2为控制器，3为功放，4为扁平转子系统A端传感器探头，5为扁平转子系统A端磁轴承，6位磁悬浮扁平转子，7为扁平转子系统B端磁轴承，8为扁平转子系统B端传感器探头；图2为本专利技术的一种基于复系数法的磁悬浮扁平转子系统章动频率解析方法流程图；图3为转子轴心轨迹图，其中，9代表章动，10代表转速方向。具体实施方式下面结合附图以及具体实施方式进一步说明本专利技术。如图1-3所示，本专利技术一种基于复系数法的磁悬浮扁平转子系统章动频率解析方法的具体步骤如下：(1)如图1-2所示，本专利技术针对磁悬浮扁平转子系统在闭环控制条件下，分析转子绕x轴、y轴偏转运动时所受到的由各通道磁轴承控制力产生的控制力矩，根据欧拉运动学方程，建立高速磁悬浮扁平转子转动动力学模型：式中Jx为x轴方向转动惯量，Jy为y轴方向转动惯量，Jz为z轴方向转动惯量，Ω为转速，α为转子绕x轴的偏转角度，β为转子绕y轴的偏转角度，lm为磁轴承力到质心的距离，fax、fbx、fay、fby为径向四通道ax、ay、bx、by方向的磁轴承力。(2)如图1所示，其中，1为传感器放大倍数，2为控制器，3为功放，4为扁平转子系统A端传感器探头，5为扁平转子系统A端磁轴承，6位磁悬浮扁平转子，7为扁平转子系统B端磁轴承，8为扁平转子系统B端传感器探头；对于分散PID磁轴承控制系统，各控制通道等效为控制参数相同的单自由度控制系统，控制电流由传感器反馈的转子位移信号引起，根据磁悬浮转子系统前向控制通道，可将控制电流i与转子位移qs间的频域关系表示为：i(s)＝-φi(s)qs(2)式中qs为传感器位置处检测到的转子径向位移，φi(s)是由控制器Gk、功放Ga、传感器放大倍数ks等环节构成的控制增益，其频域关系可表示为φi(s)＝ksGk(s)Ga(s)。选择Gk为工业控制的典型PID算法：式中kp为比例环节系数，ki为积分环节系数，kd为微分环节系数。功放环节Ga可等效为典型一阶惯性环节，具体表达形式如下：式中kw功本文档来自技高网...

【技术保护点】
一种基于复系数法的磁悬浮扁平转子系统章动频率解析方法，其特征在于：通过如下步骤实现：(1)根据欧拉运动学方程，建立高速磁悬浮扁平转子转动动力学模型：

【技术特征摘要】
1.一种基于复系数法的磁悬浮扁平转子系统章动频率解析方法，其特征在于：通过如下步骤实现：(1)根据欧拉运动学方程，建立高速磁悬浮扁平转子转动动力学模型：式中Jx为x轴方向转动惯量，Jy为y轴方向转动惯量，Jz为z轴方向转动惯量，Ω为转速，α为转子绕x轴的偏转角度，β为转子绕y轴的偏转角度，lm为磁轴承力到质心的距离，fax、fbx、fay、fby为径向四通道ax、ay、bx、by方向的磁轴承力；(2)对于分散PID磁轴承控制系统，各控制通道等效为控制参数相同的单自由度控制系统，根据磁悬浮转子系统前向控制通道，可将控制电流i与转子位移qs间的频域关系表示为：i(s)＝-φi(s)qs式中qs为传感器位置处检测到的转子径向位移，φi(s)是由控制器Gk、功放Ga、传感器放大倍数ks等环节构成的控制增益，其频域关系可表示为φi(s)＝ksGk(s)Ga(s)；由于磁轴承力可在参考工作点近似线性化为f＝kii+khqh，因此，可将径向四通道ax、ay、bx、by方向的磁轴承力由控制系统参数表示如下：式中ki为电流刚度，kh为位移刚度，ls为传感器探头到转子质心的距离，表示卷积算子；(3)回转体结构转子，其x、y轴转动惯量相等，即Jx＝Jy＝Je，令Je为赤道转动惯量，可得系统时域磁悬浮转子转动动力学方程为：

【专利技术属性】
技术研发人员：郑世强，杨景裕，周冲，韩邦成，刘刚，
申请(专利权)人：北京航空航天大学，
类型：发明
国别省市：北京,11