专利名称：基于三参量控制的振动台软启动控制方法
技术领域：
本发明涉及一种用于三参量控制振动台软启动的控制方法，属于试验技术领域。
背景技术：
三参量控制算法于20世纪70年代提出，主要为了补偿地震模拟振动台控制中只有位移信号控制时，频带不宽，不能实现加速度、速度控制的缺陷。三参量控制中，低频段由 位移控制、中频段由速度控制、高频段由加速度控制，从而达到扩宽控制频带的目的。由于 其控制频带宽、控制方式简便等优点，三参量现已逐渐应用于机械、采矿等领域的振动台控 制。三参量主要通过图1、2方式实现，输入信号通过图1生成位移、速度、加速度三个信号， 而后通过图2将三个信号合成系统驱动指令信号。由图1可得生成信号(加速度、速度、位 移)的传递函数如下式
<formula>formula see original document page 4</formula>
<formula>formula see original document page 4</formula>式中<formula>formula see original document page 4</formula>d_，vfflax, afflax分别为振动台最大设计位移、速度、加速度，根据实际情况由最大功
max
能曲线决定。通过调整Dw、nw、Ra、Rv、Rd等参数可以分别实现相当于不同输入量的控制方 式(1)当nw< ，DW<1时，输入信号与加速度参量相当，引入系统控制可以实现加速度 控制；(2)当取 =V^，DW > 1，RV = Ra时输入信号与速度相当，引入系统控制可实现速 度控制；(3)当nw > (0，Dw < 1，Rd = Ra时输入信号与位移相当，引入系统控制可实现位移 控制。在运用上述三参量实现加速度控制时，如果在启动阶段有高频大幅值信号，会出 现位移大幅度超调现象，例如中小型地震模拟振动台在进行5Hz正弦加速度输入控制时， 振动台启动阶段位移超调量达到600%，一般幅值输入都会大大超出振动台位移最大限值， 从而出现振动台台体与侧壁的碰撞现象，严重时会使实验试件在实验启动阶段毁坏。目前 解决超调的常用方法为将输入信号乘以一个如式x(t) = l-e_at的加窗函数，以减小开始 段激励信号量，对整个系统而言，相当于加入了缓冲。这种方法存在的不足在于(1)在进 入稳定工作段之前的1 2s内，激励信号被过大削弱。三参量较多的用于试验控制系统，可 很多试验往往在几秒内就完成了，比如缩尺比例较大的地震模拟振动台试验一般只有3 6s，峰值加速度往往就在1 2s内，这样很难准确满足试验要求；(2) —旦选定针对系统的 加窗函数，无论输入波产生多大超调量，同一时刻的调整比例都一样，无法根据实际输入波进行具体调整；(3)这种方式在整个试验在线过程中都要进行指数函数的求值及信号乘法 运算，一定程度上增加了在线计算量。因此，研究高效、高精度的修正三参量超调的算法对 提高振动台实验精度具有重大意义。
发明内容
本发明提出了一种离线消除三参量超调的超调修正算法，特别是一种基于三参量 控制的振动台软启动控制方法，该算法将三参量传递函数转换为微分方程，通过微分方程 求解得到三参量控制算法在三角函数(正弦或余弦)输入下的超调量时域解析解，而后结 合傅里叶变换得到任意输入下的超调量近似时域解析解，从而反向修正三参量控制算法位 移超调。具有在线计算量小、超调修正精度高、多输入信号本身无影响等特点，在振动台启 动阶段，能很好的保证信号的保真度，能大大提高振动台实验，特别是大缩尺实验的控制精 度。为了实现上述目的，本发明采用以下技术手段实现通过求解三参量控制算法的微分方程，得到任意输入下三参量产生的位移超调量 (近似解析解)，在三参量生成的位移驱动信号ed基础上减去位移超调量，从而消除三参量 控制算法中产生的超调，达到软启动的目的，包括以下步骤步骤1、根据振动台三参量设计参数求取加速度控制下三参量传递函数；步骤2、将三参量传递函数经拉普拉斯反变换得到三参量控制微分方程；步骤3、求解三角函数输入时三参量控制微分方程，得到此时的位移超调量；步骤4、将任意函数形式输入通过傅里叶变换，将任意输入变换为有限个三角函数 输入之和；步骤5、以步骤3为基础，得到步骤4中每个三角函数输入下三参量控制算法产生 的位移超调，并求和得到任意输入下的位移超调；步骤6、在输入信号经三参量控制算法产生的位移驱动信号中减去步骤5求得的 位移超调，从而消除三参量产生的位移超调。所述步骤1具体为根据振动台设计标准参数最大设计位移d_、速度Vmax、加速度 a_，及相关三参量设计参数R、Rd、Rv、Ra求得三参量控制下对应位移、速度、加速度信号传递 函数，其公式如下
R 2
ea ---G s u0
K "
r n尺 ^ 1ev=Y "TcSU° >
R r l
e, =--Cr -un
“d w D p r1"
^12 R 1 n R \ \ n 1式中Gw ^ -’ nw =—,0 >
S +2Dwnws + nwRd RXR2CXC2Rv R,CX 2nw o)f2
Ri，叫为振动台最大功能曲线对应频率下限、上限值，CpQ
"maxmax
为模拟电路中采用电容值，Dw， 分别为以上求得三参量二阶传函的频率和阻尼比，R为三
5参量生成电路放大电阻阻值，Rd、Rv、民分别为三参量产生位移、速度、加速度对应放大电阻 阻值，u0为输入加速度对应位移，s为拉普拉斯变化因子。所述步骤2具体为将所述传递函数经拉普拉斯反变换得到三参量控制算法对应 微分方程<formula>formula see original document page 6</formula>式中，wn = nw，ζ=Dw，f(t)为参考输入信号，nw、Dw分别 为以上求得三参量二阶传函的频率和阻尼比。所述步骤3具体为设输入为f(t) =Pcos(cot+￥)的三角函数时(t为时间，￥、 P为分别三角函数相位和幅值)，通过求解步骤2中微分方程得到三参量控制算法产生的位 移超调u。(t)，其公式如下<formula>formula see original document page 6</formula>
式中<formula>formula see original document page 6</formula>
所述步骤4具体为傅里叶变换的基本表达式如下
<formula>formula see original document page 6</formula>
式中…为第j点对应下的傅里叶幅值，t为时间，i为虚数单位； 在2^1周期内，以Ji为基点求和后，上式与下式等效<formula>formula see original document page 6</formula>式中wi，＃i， Pi分别为输入记录傅里叶变换后第j个三角函数对应的频率、相位 及幅值。即傅里叶变换将长度为N的任意输入展开成n个余弦函数之和，次数N为偶数n =N/2+1，N 为奇数:n= (N-l)/2+l。所述步骤5具体为即当输入f(t)为任意函数的加速度信号时，通过傅里叶变 换将其变换成有限个三角函数输入之和，而后求解得到三参量控制算法产生的位移超调 ura(t)，其公式如下<formula>formula see original document page 6</formula>式中，
<formula>formula see original document page 6</formula><formula>formula see original document page 7</formula>USti, = pl/W2n，wi， Vi，Pi分别为输入记录傅里叶变换后第i个三角函数对应的频
率、相位及幅值。与现有技术相比，本发明的优点如下1)针对三角函数输入和任意输入信号，直接求出了超调信号的解析解，修正精度 更高；2)本发明安全采用离线计算，在线只有简单的信号叠加，大大减小了在线计算量， 提高了控制器处理效率。
图1三参量信号发生示意图；图2三参量信号合成示意图；图3本发明方法实现流程图；图4使用本发明前后三参量控制算法产生的位移信号。
具体实施例方式本发明的技术方案参见图3所示，结合相关附图，下面详细介绍本发明的实施步 骤(1)、根据振动台设计标准参数最大设计位移d_、速度v_、加速度a_，及附图1 中相关三参量设计参数R、Rd、Rv、Ra求得三参量控制下对应位移、速度、加速度信号传递函 数，其公式如下
<formula>formula see original document page 7</formula>式中<formula>formula see original document page 7</formula>w1，w2为振动台最大功能曲线对应频率下限、上限值，Ci; C2为模拟电路中采用电容值，Dw， 分别为以上求得三参量二阶传函的频率和阻尼比，R为三 参量生成电路放大电阻阻值，Rd、Rv、民分别为三参量产生位移、速度、加速度对应放大电阻 阻值，u0为输入加速度对应位移，s为拉普拉斯变化因子。(2)将三参量传递函数经拉普拉斯反变换得到三参量控制微分方程将以上传递函数经拉普拉斯反变换得到三参量控制算法对应微分方程ii + 2^conu + co2nu = fit)，式中，con =nw, I = Dw，f(t)为参考输入信号； (3)求解三角函数(正弦或余弦)输入时三参量控制微分方程，得到此时的位移超 调量设输入为形如f(t) =Pcos(cot+￥)的三角函数(t为时间，V、P为分别三角函数相
位和幅值)，可通过求解以上微分方程得到 如下参量控制算法产生的位移超调u。(t)，其公式<formula>formula see original document page 8</formula>式中<formula>formula see original document page 8</formula><formula>formula see original document page 8</formula>
<formula>formula see original document page 8</formula>
(4)将任意函数形式输入通过傅里叶变换，将任意输入变换为有限个三角函数输 入之和<formula>formula see original document page 8</formula>式中co” Vi，Pi分别为输入记录傅里叶变换后第j个三角函数对应的频率、相位 及幅值。即傅里叶变换将长度为N(地震记录缩尺后N应扩大相应倍数)的任意输入展开 成n (N为偶数n = N/2+1 ；N为奇数n = (N-l) /+1)个余弦函数之和；(5)以步骤(3)为基础，得到步骤(4)中每个三角函数输入下三参量控制算法产生 的位移超调，并求和得到任意输入下的位移超调即当输入f(t)为任意函数的加速度信号 时，通过傅里叶变换将其变换成有限个三角函数输入之和，而后求解得到三参量控制算法 产生的位移超调u (t)，其公式如下<formula>formula see original document page 8</formula>式中，<formula>formula see original document page 8</formula><formula>formula see original document page 8</formula>
UStt = pjco]^,, ￥ijPi分别为输入记录傅里叶变换后第i个三角函数对 应的频率、相位及幅值；(6)在输入信号经三参量控制算法产生的位移驱动信号ed中减去步骤(5)求得的 位移超调u。,，从而消除三参量产生的位移超调；即由附图2将步骤(1)中产生的三个信号 叠加，并减去u 后可得振动台前馈驱动信号U(l，其公式如下u0 = Ad (ed-ucr) +Avev+Aaea式中Ad、Av、Aa分别为三参量控制下位移、速度、加速度前馈增益。本发明的效果图如图4，可以看出采用该方法后三参量控制算法下位移超调得到 了明显改善，在保证振动台安全运行的同时提高了控制精度。最后应说明的是以上发明内容仅用以说明本发明，而并非限制本发明所描述的 技术方案；本发明并非只适用于振动台的控制，三参量控制算法涉及的相关超调问题的解 决均应涵盖在本发明的权利要求范围当中。
权利要求
一种基于三参量控制的振动台软启动控制方法，其特征在于包括以下步骤步骤1、根据振动台三参量设计参数求取加速度控制下三参量传递函数；步骤2、将三参量传递函数经拉普拉斯反变换得到三参量控制微分方程；步骤3、求解三角函数输入时三参量控制微分方程，得到此时的位移超调量；步骤4、将任意函数形式输入通过傅里叶变换，将任意输入变换为有限个三角函数输入之和；步骤5、以步骤3为基础，得到步骤4中每个三角函数输入下三参量控制算法产生的位移超调，并求和得到任意输入下的位移超调；步骤6、在输入信号经三参量控制算法产生的位移驱动信号中减去步骤5求得的位移超调，从而消除三参量产生的位移超调。
2.根据权利要求1所述的基于三参量控制的振动台软启动控制方法，其特征在于所 述步骤1具体为根据振动台设计标准参数最大设计位移d_、速度v_、加速度amax，及相关 三参量设计参数R、Rd、Rv、Ra求得三参量控制下对应位移、速度、加速度信号传递函数，其公 式如下<formula>formula see original document page 2</formula>式中 <formula>formula see original document page 2</formula>w1,w2，为振动台最大功能曲线对应频率下限、上限值，Cl，C2为模拟电maxmax路中采用电容值，Dw， 分别为以上求得三参量二阶传函的频率和阻尼比，R为三参量生成 电路放大电阻阻值，Rd、Rv、Ra分别为三参量产生位移、速度、加速度对应放大电阻阻值，％为 输入加速度对应位移，s为拉普拉斯变化因子。
3.根据权利要求2所述的基于三参量控制的振动台软启动控制方法，其特征在于所 述步骤2具体为将所述传递函数经拉普拉斯反变换得到三参量控制算法对应微分方程<formula>formula see original document page 2</formula>式中，wn = nw, ζ =Dw,f(t)为参考输入信号，nw、Dw分别为以上求 得三参量二阶传函的频率和阻尼比。
4.根据权利要求3所述的基于三参量控制的振动台软启动控制方法，其特征在于所 述步骤3具体为设输入为<formula>formula see original document page 2</formula>的三角函数时，其中t为时间，V、P为分 别三角函数相位和幅值，通过求解步骤2中微分方程得到三参量控制算法产生的位移超调 u。(t)，其公式如下<formula>formula see original document page 2</formula>式中<formula>formula see original document page 2</formula><formula>formula see original document page 3</formula>
5.根据权利要求4所述的基于三参量控制的振动台软启动控制方法，其特征在于所 述步骤4具体为傅里叶变换的基本表达式如下<formula>formula see original document page 3</formula>式中…为第j点对应下的傅里叶幅值，t为时间，i为虚数单位； 在2JI周期内，以为基点求和后，上式与下式等效<formula>formula see original document page 3</formula>式中Wi，U/yPi分别为输入记录傅里叶变换后第j个三角函数对应的频率、相位及幅 值。即傅里叶变换将长度为N的任意输入展开成n个余弦函数之和，次数N为偶数m = N/2+1，N 为奇数n = (N-l)/2+l。
6.根据权利要求5所述的基于三参量控制的振动台软启动控制方法，其特征在于所 述步骤5具体为即当输入f(t)为任意函数的加速度信号时，通过傅里叶变换将其变换成 有限个三角函数输入之和，而后求解得到三参量控制算法产生的位移超调(t)，其公式如 下<formula>formula see original document page 3</formula>式中，<formula>formula see original document page 3</formula>,=乃/仍 2, q，％，A分别为输入记录傅里叶变换后第i个三角函数对应的频率、相位及幅值。
全文摘要
本发明涉及一种用于消除三参量控制下振动台启动过程中产生位移超调的控制方法。该方法在三参量控制的基础上提出了离线时域超调修正算法通过对时域内三参量控制微分方程求解得出三参量控制算法产生的超调量理论解，进而从直接消除暂态分量出发，以三角函数、傅里叶变换为工具，提出了任意输入下超调量的修正方法。和传统的加窗软启动相比，在降低超调量的同时，大大降低了在线计算量，从而在消除振动台位移超调的同时提高了控制器运行效率。
文档编号G01M7/02GK101832849SQ20101014574
公开日2010年9月15日 申请日期2010年4月9日 优先权日2010年4月9日
发明者唐贞云, 李振宝, 李晓亮, 纪金豹 申请人:北京工业大学