CN101813552A

CN101813552A - 用于补偿试件台面相互作用的地震模拟振动台控制方法

Info

Publication number: CN101813552A
Application number: CN 201010145743
Authority: CN
Inventors: 李振宝; 唐贞云; 纪金豹; 李晓亮
Original assignee: Beijing University of Technology
Current assignee: Beijing University of Technology
Priority date: 2010-04-09
Filing date: 2010-04-09
Publication date: 2010-08-25
Anticipated expiration: 2030-04-09
Also published as: CN101813552B

Abstract

本发明涉及一种用于补偿试验试件与台面相互作用、提高地震记录再现精度的地震模拟振动台控制方法。该方法在将试件按刚性负载进行振动台控制参数设计的基础上提出了所有负载形式均按零负载进行振动台控制参数设计，此外利用试件产生反力实时补偿试件-台面相互作用的控制方法。与传统控制方法相比：该方法提高了系统控制的稳定性，使得地震模拟振动台更具可控性；扩宽了系统可控频带，使得振动台可进行的试验范围更宽；提高了地震记录再现精度，使得地震模拟振动台试验更具科学性。

Description

用于补偿试件台面相互作用的地震模拟振动台控制方法

技术领域

本发明涉及一种用于补偿试件台面相互作用的地震模拟振动台控制方法，属于建筑结构抗震试验技术领域。

背景技术

地震模拟振动台控制系统算法设计时通常将台面与试件看作整体以单自由度系统方式处理：图1为地震模拟振动台组要组成架构，以图中伺服阀为研究对象，可建立伺服阀在零工作点附近的线性压力-流量方程，以作动缸为研究对象，可建立作动缸需求流量关系，；以活塞头为研究对象可建立作动器出力平衡方程。以上三个关系式统称为液压系统三连续方程，其表达式如下：

\begin{matrix} Q_{L} = k_{q} E - K_{c} p_{L} \\ Q_{L} = A_{p} {\overset{\cdot}{x}}_{t} + \frac{V}{4 β} {\overset{\cdot}{p}}_{L} + C_{c} p_{L} \\ (M_{t} + M_{p}) {\overset{\cdot \cdot}{x}}_{t} = A_{p} p_{L} \end{matrix}\}

式中k_q为流量增益，K_c为流量-压力系数，E为阀驱动信号，A_p为活塞有效面积，x_t为活塞、台面位移，V为油缸等效容积，β为油体积弹性模量，p_L为油缸内P_A与P_B间压力差，C_c为活塞总泄系数，M_t为振动台台面总质量，M_p为试验试件总质量。而后利用三参量前馈和三参量反馈作为振动台控制的主要控制算法，主要控制流程如图2。现有控制算法将试件看作台面质量的一部分，按刚性处理，以此求取图2中位移、速度加速度前馈增益(A_d、A_v、A_a)，位移、速度加速度反馈增益(A′_d、A′_v、A′_a)，以此进行地震模拟振动台系统控制，模拟试验试件在实际地震作用中的性能。

这种设计方法使得负载与空台的唯一区别只在于台面质量的增加，可试件与振动台相比二者特性差别较大，常规中小型振动台共振频率在20Hz左右，试验试件一阶频率大都在2～8Hz。因此，将试件与台面统一为一体存在较大偏差，从1988年开始国内外学者相继从理论及试验两方面对试件与台面相互作用进行了研究，结果表明负载在振动台系统中的参与质量比例与频率相关，只有输入波频率小于负载频率时，负载质量才完全作为刚性质量与台面成为一体，因此，负载频率越高参与质量比例也就越高，在大于负载频率3Hz以后的频带范围，负载参与质量不到负载质量的10％。换言之，当输入波频率处在大于负载频率3Hz以后的频带范围时，振动台系统中将不再有负载质量参与，相当于控制空载系统。而且其参与比例随其自振频率的减小而降低，中高频段负载质量很少参与到振动台物理系统中，此时其特性与空载更接近，而常规中小型振动台要求的频率都在50Hz以上，采用刚性负载参数时相当于用一大质量空台参数控制小质量空台，因而当试件质量较大、频率较低时还有可能使得振动台系统不稳定，即不可控状态；此外，试件产生的反力反过来会影响振动台的控制性能，由于共振的影响，试件频率附件的波形产生更大的反力，使得试件频率成为整个振动台系统的第二共振频率，对结构自振频率周围的时程跟踪性能有很大影响，该频率附近的输入波分量再现精度相比其它频率分量低很多；因此研究考虑试件-台面相互作用的振动台控制算法对提高振动台再现地震记录精度具有重要意义。

发明内容

本发明提出了一种能考虑试验试件与台面间相互作用的控制方法，特别是一种用于补偿试件台面相互作用的地震模拟振动台控制方法。该方法利用空载条件下振动台三参量控制参数控制所有试验试件工况；利用实测试验值获取试件对台面的反力，实时补偿试件与台面间相互作用，运用该算法对振动台进行控制具有系统稳定性好、控制精度高、适于试验工况范围广等特点，对基频较高的试验试件同样可以取得较好控制效果，特别是对缩尺比较大的试件在压缩比较大的地震激励作用下的控制效果能有很大幅度提高。

为了实现上述目的，本发明采取了如下技术方案：

一种用于补偿试件台面相互作用的地震模拟振动台控制方法，其为了保证系统稳定性，采用空载系统设计参数控制所有试验试件工况；为了消除试件与台面之间相互作用，通过实测分析得到试件对振动台台面反力，而后经伺服阀逆传函求解得到试件反力等效驱动信号，与原驱动信号叠加，消除二者相互作用，达到提高振动台控制精度的目的。

主要包括以下步骤：

步骤1、根据地震模拟振动台基本设计参数，设计其在没有试验试件条件下的系统控制参数；

步骤2、在空载设计参数控制条件下，通过试验实测、分析得到试验试件对振动台系统的反作用力；

步骤3、求取振动台电液伺服系统逆传递函数；

步骤4、结合实测反力及电液伺服系统逆传函计算试件产生反力的等效驱动信号；

步骤5、将试件反力产生的等效驱动信号与原有驱动信号叠加，得到新的驱动信号，对振动台系统进行控制。

所述步骤1中，在进行地震模拟振动台控制参数设计时，试验试件与振动台台体组成的被控系统同作动器的出力之间的关系由下面的力平衡公式表示：

(M_{t} + M_{p}) {\overset{\cdot \cdot}{x}}_{t} = A_{p} p_{L}

式中M_t为振动台台面总质量，M_p为试验试件总质量，x_t为台面位移，A_p为作动器活塞有效面积，p_L为伺服阀产生压差；

令M_p≡0，根据振动台已有或设计硬件参数，运用三参量控制参数常规求取方法，求出位移、速度、加速度输入增益(A_d、A_v、A_a)，位移、速度、加速度反馈增益(A′_d、A′_v、A′_a)作为任何负载条件振动台的控制参数。

所述步骤2具体为：当振动台作动器反力不可测时，测量试验试件各自由度对应加速度，运用公式计算试验试件产生的反力，式中m_n、

分别为试验试件第n个自由度对应质量和绝对加速度；当振动台作动器反力可测时，直接测量作动器反力，运用公式

计算试验试件产生的反力，式中F_act为作动器反力，M_t、

分别为振动台台面质量及加速度。

所述步骤3具体为：计算振动台伺服阀系统逆传函的公式如下：

H_{f}^{- 1} (s) = \frac{M_{s} {Vs}^{2} + 4 β M_{t} (K_{c} + C_{c}) s + 4 β A_{p}^{2}}{4 {βA}_{p} s},

式中M_t为振动台台面总质量，V为油缸等效容积，β为油体积弹性模量，K_c为流量-压力系数，C_c为活塞总泄系数，A_p为活塞有效面积。

所述步骤4具体为：结合步骤2中的反力值和步骤3中的传递函数综合计算出反力等效驱动信号，其计算公式为：

所述步骤5具体为：新驱动信号计算公式为：

u₀＝A_de_d+A_ve_v+A_ae_a-A′_de′_d-A′_ve′_v-A′_ae′_a-Δu

式中A_d、A_v、A_a分别为位移、速度、加速度前馈增益，A′_d、A′_v、A′_a分别为位移、速度、加速度反馈增益，e_d、e_v、e_a分别为三参量算法产生的前馈位移、速度、加速度信号，e′_d、e′_v、e′_a分别为三参量算法中位移、速度、加速度反馈信号。

与现有控制算法相比，本发明的优点如下：

(1)所有负载形式均按零负载进行振动台控制参数设计提高了振动台系统控制稳定性，避免了按刚性负载条件设计的三参量控制参数可能造成的系统失控，甚至使得试验试件在短时间内被毁坏的现象发生，同时提高了控制精度；

(2)利用试件反力实时补偿大大减小了试件与台面相互作用，大大减小了在试件各阶频率附件出现的波峰和波谷，扩宽了系统可控频带，使得振动台可进行的试件缩尺范围、波形压缩范围更宽；提高了地震记录再现精度，使得地震模拟振动台试验更具科学性。

附图说明

图1地震模拟振动台基本硬件组成；

图2地震模拟振动台控制系统组成；

图3本发明方法实现流程图；

图4实时反力补偿控制算法实现架构；

图5使用本发明前后振动台控制性能Bode图。

附图标记：1-进油口；2-阀芯位移；3-出油口；4-阀芯；5-阀孔；6-伺服阀体；7-作动器；8-作动缸腔；9-振动台台体；10-试验试件

具体实施方式

本发明的技术方案参见图3所示，结合相关附图，下面详细介绍本发明的实施步骤：

(1)根据地震模拟振动台基本设计参数，设计其在没有试验试件条件下的系统控制参数：在进行地震模拟振动台控制参数设计时，试验试件与振动台台体组成的被控系统同作动器的出力之间的关系由下面的力平衡公式表示：

(M_{t} + M_{p}) {\overset{\cdot \cdot}{x}}_{t} = A_{p} p_{L}

式中M_t为振动台台面总质量，M_p为试验试件总质量，x_t为台面位移，A_p为作动器活塞有效面积，p_L为伺服阀产生压差。

令M_p≡0，根据振动台已有或设计硬件参数，运用三参量控制参数常规求取方法，求出位移、速度加速度输入增益(A_d、A_v、A_a)，位移、速度加速度反馈增益(A′_d、A′_v、A′_a)作为任何负载条件振动台的控制参数；

(2)在空载设计参数控制条件下，通过试验实测、分析得到试验试件对振动台系统的反作用力：当振动台作动器反力不可测时，测量试验试件各自由度对应加速度，运用公式

计算试验试件产生的反力，式中m_n、

分别为试验试件第n个自由度对应质量和绝对加速度；当振动台作动器反力可测时，直接测量作动器反力，，运用公式

计算试验试件产生的反力，式中F_act为作动器反力，M_t、

分别为振动台台面质量及加速度；

(3)求取振动台电液伺服系统逆传递函数：计算振动台伺服阀系统逆传函的公式如下：

H_{f}^{- 1} (s) = \frac{M_{s} {Vs}^{2} + 4 β M_{t} (K_{c} + C_{c}) s + 4 β A_{p}^{2}}{4 {βA}_{p} s},

式中M_t为振动台台面总质量，V为油缸等效容积，β为油体积弹性模量，K_c为流量-压力系数，C_c为活塞总泄系数，A_p为活塞有效面积；

(4)利用试验实测反力及电液伺服系统逆传函，计算试件产生反力的等效驱动信号：结合步骤2中的反力值和步骤3中的传递函数综合计算出反力等效驱动信号，其计算公式为：

(5)将试件反力产生的等效驱动信号与原有驱动信号叠加，得到新的驱动信号，对振动台系统进行控制，振动台控制架构中的实现如图4所示，新驱动信号计算公式为：

u₀＝A_de_d+A_ve_v+A_ae_a-A′_de′_d-A′_ve′_v-A′_ae′_a-Δu

式中A_d、A_v、A_a分别为位移、速度、加速度前馈增益，A′_d、A′_v、A′_a分别为位移、速度、加速度反馈增益，e_d、e_v、e_a分别为三参量算法产生的前馈位移、速度、加速度信号，e′_d、e′_v、e′_a分别为三参量算法中位移、速度、加速度反馈信号

本发明的效果图如图5，可以看出采用该方法后在保证振动台稳定前提下，控制频带得到了明显改善，在保证振动台安全运行的同时提高了控制精度。

最后应说明的是：以上发明内容仅用以说明本发明，而并非限制本发明所描述的技术方案；本发明并非只适用于地震模拟振动台的控制，其他领域振动台涉及的相关相互作用问题的解决均应涵盖在本发明的权利要求范围当中。

Claims

1.一种用于补偿试件台面相互作用的地震模拟振动台控制方法，其特征在于：包括以下步骤：

步骤3、求取振动台电液伺服系统逆传递函数；

2.根据权利要求1所述的用于补偿试件台面相互作用的地震模拟振动台控制方法，其特征在于：所述步骤1具体为，在进行地震模拟振动台控制参数设计时，试验试件与振动台台体组成的被控系统同作动器的出力之间的关系由下面的力平衡公式表示：

(M_{t} + M_{p}) {\overset{\cdot \cdot}{x}}_{t} = A_{p} p_{L}

3.根据权利要求2所述的用于补偿试件台面相互作用的地震模拟振动台控制方法，其特征在于：所述步骤2具体为：当振动台作动器反力不可测时，测量试验试件各自由度对应加速度，运用公式计算试验试件产生的反力，式中m_n、

计算试验试件产生的反力，式中F_act为作动器反力，M_t、

分别为振动台台面质量及加速度。

4.根据权利要求3所述的用于补偿试件台面相互作用的地震模拟振动台控制方法，其特征在于：所述步骤3具体为：计算振动台伺服阀系统逆传函的公式如下：

H_{f}^{- 1} (s) = \frac{M_{t} {Vs}^{2} + 4 β M_{t} (K_{c} + C_{c}) s + 4 β A_{p}^{2}}{4 β A_{p} s},

5.根据权利要求4所述的用于补偿试件台面相互作用的地震模拟振动台控制方法，其特征在于：所述步骤4具体为：结合步骤2中的反力值和步骤3中的传递函数综合计算出反力等效驱动信号，其计算公式为：

Δu = H_{f}^{- 1} (s) R_{tm} .

6.根据权利要求5所述的用于补偿试件台面相互作用的地震模拟振动台控制方法，其特征在于：所述步骤5具体为：新驱动信号计算公式为：

u₀＝A_de_d+A_ve_v+A_ae_a-A′_de′_d-A′_ve′_v-A′_ae′_a-Δu