Page 226 - 《振动工程学报》2026年第5期
P. 226

1430                               振     动     工     程     学     报                     第 39 卷

              行分析,采用系统可控性的概念来描述主控式弹支                            2.2.2    振幅约束
              干摩擦阻尼器对转子振动的控制能力,并由此建立                                由于航空发动机转子增减速过程较为迅速,当
              参数优化数学模型。                                         转子转速靠近临界转速区间时,需提前对转子模态
                  由于阻尼减振效果与弹性支承紧密关联,弹支                          振动进行控制,以保证控制过程有充足时间和位移
              相对应变能越大,阻尼器发挥阻尼效能的条件越有                            幅值裕度。因此,针对转子系统振动控制目标,设置
              利 [30] 。为了定量评估阻尼器在各阶模态下的适用                        转子振动预警幅值,即正常工况允许的转子振动最
              性,本文引入弹性支承总应变能在系统应变能中的                            大幅值为     A max 。
              占比作为评价指标,用以衡量阻尼器发挥效能的有                                                                     (11)
                                                                                 f(A n ) = A max − A n
              利程度,如下式所示:                                        式中,   A n 为转子振幅。当      f(A n ) > 0时,表示该转速的
                               E  eb                            振幅满足约束条件;当           f(A n ) < 0时,表示不满足振幅
                                Li
                                         eb
                       eb
                     P =              ,P ∈ [0,1)       (7)
                                         Li
                       Li
                            eb
                                rb
                          E + E + E  LP
                           Li   Li  Li                          约束条件。控制过程中要求所有模态都满足这一要
              式中,   P 为低压转子第       i 阶模态的弹性支承总应变
                     eb
                     Li                                         求,可以统一表示为:
                       eb
                                                    rb
              能占比;    E 为第   i 阶模态弹支总应变能;          E 为刚性
                                                                                     [
                                                                                          ]
                       Li                           Li                           min f (A n ) > 0        (12)

              支承的总应变能;        E 为低压转轴上的应变能。相关
                               LP
                               Li                               2.2.3    稳定性约束
              参数表达式如下        [31] :
                                                                    为保证转子增减速的稳定性,防止过约束或欠
                            
                                 W ∑
                            
                              eb      eb 2                     约束现象造成的震荡,以转子不平衡响应过程中的
                            
                            E =   0.5S r
                            
                              Li      a  Li,eb,a
                            
                            
                                a=1                            峰值个数代表稳定性,对幅值波动的变换做出相应
                            
                            
                            
                                 Z ∑                   (8)
                              rb      rb 2
                            E =   0.5S r                       要求。针对转子特性,设置最大峰值个数为                    P max 。
                            
                            
                              Li      b  Li,rb,b
                            
                            
                                b=1
                            
                            
                              LP    T  LP                                                               (13)
                                                                                f(P n ) = P max − P n
                             E = 0.5Φ S Φ Li
                              Li     Li
              式中,   S 为第   a  个弹性支承的刚度;        S 为第   b  个刚     式中,   P n 为峰值个数。当      f(P n ) > 0时,表示这一控制
                                                rb
                     eb
                     a                          b
              性支承的刚度;       r Li,rb,a 表示低压转子第   i 阶模态振型         参数下的振幅稳定性满足约束条件;当                   f(P n ) < 0时,
              中第   a  个弹性支承处的归一化位移;            r Li,rb,b 表示低压   表示振幅稳定性不满足约束条件。控制过程中要求
              转子第    i 阶模态振型中第       b  个刚性支承处的归一化              所有模态都满足这一要求,可以统一表示为:
                                                                                     [
                                                                                          ]
              位移;W   和  Z  分别为弹性支承和刚性支承的个数;               S LP                   min f (P n ) > 0        (14)
              为低压转轴的刚度矩阵;           Φ Li 为转子第  i阶振型。
                                                                2.3    优化目标函数
                  弹性支承总应变能占比是为               ESDFD  发挥阻尼
              减振作用提供保障条件的量化指标,其值越大,阻尼                               与单目标优化不同,多目标优化通常无法使每
              器 减 振 效 果 越 突 出, 模 态 可 控 度 越 高 。 实 际 操 作          个子目标同时达到最优,而是获得一个能够综合权
              时,可通过布置加速度传感器等方法,在特定工况下                           衡各目标的折衷解。因此,在发动机增减速过程的
              估计各个模态上的振动能量。本文为表征控制方法                            控制器参数选取中,难以直接确定一个绝对的最优
              特点和效果,根据不同模态分别对控制器参数进行                            点。因此,需要通过优化算法从可以选择的范围中
              优化。                                               选择出一组最优的控制参数。在发动机增减速的过

                                                                程中,期望振动幅值与参考幅值达到高度吻合,同时
              2.2    约束性指标
                                                                也保证转子振动幅值变化尽可能小,所需压力尽可

              2.2.1    作动器限制                                    能降低。因此结合减振率、稳定性和耗费代价(压
                  ESDFD  主要依靠提供的正压力产生摩擦力进行                      力)的考虑,设置目标函数为三者的线性加权,以优
              减振,因此,针对阻尼器的工作特点,应设置正压力                           化振动控制系统的性能表现。同时,需要对发动机
              限制,控制器输出的正压力不能超过压电陶瓷最大                            增减速过程中的控制器参数进行上下限的约束以表
              输出力(或者包含预紧力):                                     示可选范围。优化目标函数的表达式如下:
                                                       (9)                                               (15)
                               f(N n ) = N max − N n                  Fit = (c 1 f 1 +c 2 f 2 +c 3 f 3 +c 4 f 4 )+ f penalty
              式 中,  N n 为 实 际 输 出 力 ;  N max 为 最 大 正 压 力 。 当                   
                                                                               
                                                                                f 1 = λ 11 A 1 +λ 12 A 2
                                                                               
                                                                               
              f(N n ) > 0时,表示这一转速下的正压力满足约束条                                    
                                                                               
                                                                               
                                                                                f 2 = λ 21 N 1 +λ 22 N 2
                                                                                                        (16)
              件;当   f(N n ) < 0时,表示不满足正压力约束条件。控                                
                                                                               
                                                                                f 3 = λ 31 F 1 +λ 32 F 2
                                                                               
                                                                               
                                                                               
              制过程中要求所有模态都满足这一要求,可以统一                                             f 4 = λ 41 P 1 +λ 42 P 2
                                                                               
                                                                               
              表示为:
                                                                式中,   c i (i = 1,2,3,4)为相关参数权重值,具体值根据
                                  [     ]
                               min f (N n ) > 0        (10)     优化侧重点选取;        λ ij 为各子目标函数     f i (i = 1,2,3,4)的
   221   222   223   224   225   226   227   228   229   230   231