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

第 5 期                     孙镜茹,等:调谐质量阻尼惯容器振动控制性能试验研究                                        1443

              滚珠丝杠的结构特点,可得出滚珠丝杠式惯容器的                            2.3    阻尼系数识别
              出力公式    [16]  如下:
                                                                    系统中的阻尼分为两部分,一部分是主结构的
                              (  ) 2
                               2π                               阻尼,主要来自于系统的摩擦,包括滑轨间的摩擦
                           p =     (I 1 + I 2 )(¨u 2 − ¨u 1 )  (20)
                               l d
                                                                等;另一部分是       TMDI 振子的阻尼,主要来自于振子
              式中,p   为惯容力;     ¨ u 1 和  ¨ u 2 分别为惯容两端点的加速       系统内的摩擦,以及滚珠丝杠的阻尼作用。
                         (  ) 2
                          2π
              度;定义    b =     (I 1 + I 2 ) ,则  b  为惯性系数;  I 1 和 分   主结构的阻尼系数可通过自由衰减试验进行测
                                                         I 2
                          l d
                                                                定,衰减位移时序如图           7  所示。依据自由衰减理论,
              别为滚珠丝杠螺母和飞轮的转动惯量。
                                                                通过分析试验过程中主结构的振动衰减特性,计算
                  试验中,选用型号为         1620  的滚珠丝杠作为惯容元
                                                                得出其阻尼比:
              件驱动件,其导程为         20 mm。由于惯性元件飞轮为圆                                        (     )
                                                   1                                  ln x n /x n+ j     (21)
                                                           2
              环状,因此可通过圆环转动惯量公式:                 I 2 = m(r +r )                  ζ main =
                                                        2
                                                   2   1   2                             2πj
              求得飞轮的理论转动惯量,其中              r 1 和 为圆环的内径          式中,   x n 、x n+ j 分别表示第  n  个、第  n+j 个位移峰值的
                                              r 2
              和外径。根据各飞轮的理论转动惯量,可计算各飞                            大小。根据下式可求得主结构阻尼系数                    c 1 ,其中主
              轮的理论惯性系数。                                         结构的刚度      k 1 和质量 m 1 均可直接测量得到:
                  采用自由衰减试验测量不同飞轮的惯性系数,                                                   √               (22)
                                                                                 c 1 = 2ζ main k 1 m 1

              通过对质量块施加一个初值位移,观察结构在没有                                    30                          无控制
              外 力 作 用 下 的 自 由 振 动 过 程, 根 据 振 动 过 程 中 的                  20
                                                                        10
              响应数据来识别飞轮的惯性系数。系统的整体平动                                  主结构位移 / mm  −10 0
                                                                       −20
              质量由质量块和滚珠丝杠的丝杠轴,以及转动的惯                                   −30
                                                                       −40
              性元件产生的等效惯性质量所构成。已知质量块和                                      0       1      2       3       4
                                                                                       时间 / s
              丝杠轴的物理质量后,用整体平动质量减去其余物
                                                                          图 7 主结构自由衰减位移时序图
              理质量,可以得到该工况下系统的等效惯性质量,即
                                                                Fig. 7 Free  decay  displacement  time  history  of  the  main
              试验惯性系数。对比不同飞轮工况,可以得到不同
                                                                      structure
              飞轮的试验惯性系数和表观质量放大系数,如表                        1
                                                                    根据前述的      TMDI 控制系统的理论公式,振子的
              所示。
                                                                质量、刚度和惯性系数可以直接通过测量获得,而
                 表 1 飞轮理论、试验惯性系数和表观质量放大系数                       振子的阻尼参数则无法直接通过理论或试验确定,
              Tab. 1 Theoretical  and  experimental  inertia  coefficients  and
                                                                因此本试验采用基于强迫振动的参数识别方法来识
                    apparent mass amplification factor of the flywheels
                                                                别阻尼系数。强迫振动试验常用于研究结构或系统
                    理论惯性 实际惯性
               飞轮                   误差/% 实际质量/g 放大系数            的振动特性,以及识别系统的动力学参数。在强迫
                     系数/g    系数/g
                                                                振动试验中,可以使用不同类型的激励信号(如正弦
                      368     360    −2.17     8       45
                J 1
                                                                波、方波、随机信号等),并通过传感器采集振动响
                     1444     1410   −2.35    21      67.14
                J 2
                                                                应数据,如位移、速度或加速度等。通过分析激励
                      674     680     0.89    17       40
                J 3
                                                                信号与振动响应之间的关系,可以得出系统的频率
                     2896     2844   −1.80    42      67.71
                J 4
                                                                响应函数、振动模态、阻尼比等信息。对                   TMDI 振子
                     1342     1350    0.60    34      39.70
                J 5
                                                                施加正弦力激励,采集位移和力信号,使用                    Matlab  基
                     2327     2448    5.20    54      45.33
               MJ 1
                                                                于参数优化算法拟合得到振子阻尼系数。
                     4864     5033   −3.47    128     39.32
               MJ 2

                                                                2.4    参数识别方法有效性验证
                  分析试验数据,理论惯性系数与试验测得的数
              值高度吻合,误差在          5%  左右,充分验证了通过自由                   将识别出的惯性系数和阻尼参数应用于该数值
              衰减法测量飞轮惯性系数的准确性和可靠性,为后                            仿真模型中,在定频强迫激励情况下,选取                    2、2.4  和
              续的性能试验提供了可靠的数据基础。需要指出                             3 Hz 的激励进行交叉验证。将采集到的位移数据划
              的是,误差主要来自两方面:一方面是为了方便安                            分为训练集和测试集,其中,2 Hz 工况作为训练集,
              装,飞轮并未设计为严格的圆环形状,这一非理想形                           2.4  和  3 Hz 工况作为测试集。利用         2 Hz 训练集进行
              态在一定程度上影响了惯性系数的测量精度;另一                            参数识别,得到阻尼参数;利用识别得到的阻尼参数
              方面是阻尼的存在会对惯性测量结果产生一定的                             和之前识别得到的惯性系数,对                2.4  和  3 Hz 测试集
              影响。                                               进行拟合,如图       8  所示。通过比较       2.4  和  3 Hz 仿真结
   234   235   236   237   238   239   240   241   242   243   244