一种基于非概率可靠性优化的结构振动极点配置控制系统设计方法

文档序号:10723701阅读:203来源:国知局
一种基于非概率可靠性优化的结构振动极点配置控制系统设计方法
【专利摘要】本发明涉及一种基于非概率可靠性优化的结构振动极点配置控制系统设计方法,该方法首先建立结构振动极点配置的主动控制闭环系统,然后建立了闭环输出反馈主动控制系统非概率可靠度计算方法,基于含区间参数的闭环输出反馈主动控制系统特征值的区间估计方法和非概率可靠度计算方法对输出反馈控制器进优化,最终得到满足可靠度要求的基于极点配置的输出反馈控制系统。本发明从可靠性的角度确定输出反馈控制器的参数,从而解决了传统结构振动极点配置输出反馈控制器设计过程中面临的不确定性问题,使得基于极点配置得到的闭环控制系统在参数扰动的情况下也能够满足设计要求。
【专利说明】
一种基于非概率可靠性优化的结构振动极点配置控制系统设 计方法
技术领域
[0001] 本发明涉及结构振动主动控制的技术领域,具体涉及一种基于非概率可靠性优化 的结构振动极点配置控制系统设计方法。
【背景技术】
[0002] 振动与噪声的主动控制是一个快速发展的新领域。振动与噪声的主动控制技术涉 及多学科交叉耦合,具有广泛的应用价值和学术价值,是航空航天结构走向智能化的开始。 振动与噪声主动控制就是通过引入次级干扰(振动或者噪声)来消除初级干扰形成的振动 或噪声。次级干扰是通过驱动器产生,次级干扰的具体形式由控制器经过计算得出。振动与 噪声主动控制技术的发展不是一蹴而就的,也经过了一定时期的积累。现有的控制器设计 方法大都是根据控制理论提出的方法进行的,主要包括极点配置、最优控制、PID控制和鲁 棒控制等理论。当然,随着现代控制理论的发展,结构振动与噪声主动控制系统控制器的设 计也是在前进和发展的。
[0003] 极点配置方法是一种非常简单的反馈控制器设计方法。这种设计方法的目的就是 保证闭环控制系统矩阵的特征值与设计者要求的特征值接近。通常来说,设计者可以把闭 环系统矩阵的特征值设计到任何位置,这在理论上可行的。但是,在考虑各种其他物理因素 的情况下,实际工程中并不这样做。例如,可以通过设计控制器使得闭环系统的特征值更加 靠近复平面的左边,这样就可以提高系统稳定的鲁棒性。Manning等利用系统识别和极点配 置控制方法实现了柔性智能梁第一阶模态的振动主动控制,结果表明极点配置方法可以媲 美速度反馈控制器。Sethi等针对一个3.5m长的复合材料工字型梁,利用压电材料作为传感 器和驱动器设计了极点配置控制器,结果表明基于极点配置设计的控制器充当了阻尼器的 作用。不仅保持了原来系统的稳定性,同时还增加了系统对外部干扰的鲁棒性。Kumar等对 智能结构的振动主动控制系统的极点配置可行性进行了研究,考虑了结构参数的不确定 性,提出了鲁棒极点配置方法,该方法不仅可以保证系统的鲁棒稳定性,而且还能保证系统 的鲁棒性能。Hanagan和Murray等研究了基于极点配置的结构振动抑制方法,通过施加速度 反馈达到对闭环系统极点进行调节的目的,从而实现结构振动响应的控制。JUNKINS等基于 最小灵敏度设计方法得到了极点配置的输出反馈控制器,实现了振动响应的主动控制。 Slater等人利用约束状态和输出反馈控制设计了基于特征值配置的主动控制系统,并将其 应用于柔性结构的振动控制。Bittanti等利用简单的极点配置方法提出了一种新的适用于 直升机的叶片振动控制的方法,该方法具有两个优点:一、可以不影响整个直升机的飞行动 力学;二、控制器的设计和计算简单有效。Ram研究了多输入系统的极点配置方法,单输入系 统不能解决闭环控制系统特征值相关性的问题,没法保证在其他特征值不变的情况下实现 要求的特征值变化,Ram通过多输入控制方法对特征值相关的极点配置问题进行了研究。 Abdelaziz等利用状态微分反馈控制器研究了极点配置在单输入单输出线性系统中的应 用,其求解过程是与Ackermann公式相似的,同时该方法还可以应用于结构振动控制时变系 统。路小波等利用试验对系统进行辨识,设计了基于极点配置的柔性智能结构主动振动控 制器,并对仿真结果进行了试验。李书等针对含不确定参数的结构系统,提出了应用 Householder变换的方法设计具有鲁棒性的极点配置控制器。张家凡针对具有非对称阻尼 和非对称刚度矩阵的一般动力学系统,研究了多输入极点配置问题,给出了控制增益矩阵 的显式解。周星德等提出反向极点配置,该控制策略能够实现部分极点的重新配置,并且保 证了不会出现控制输出力较大的情况,具有很高的工程实用价值。
[0004] 虽然基于极点配置的闭环控制系统的设计方法已被广泛用于振动控制领域,但是 当系统参数存在不确定性时,利用名义系统得到的闭环控制系统有时候并不能满足设计的 可靠度要求。本发明就是从非概率可靠性的角度出发,基于非概率可靠性优化提出的一种 可以考虑模型不确定性的极点配置振动主动控制方法,使得到的闭环控制系统在满足可靠 性要求的前提下,达到要求的振动控制效果。现有专利文献和非专利文献均无相关技术的 报道。

【发明内容】

[0005] 本发明要解决的技术问题是:克服现有极点配置方法无法处理模型存在不确定性 的缺点,提供一种基于非概率可靠性优化的结构振动极点配置主动振动控制系统设计方 法,从而提高极点配置主动控制系统的可靠性。
[0006] 本发明解决上述技术问题采用的技术方案为:一种基于非概率可靠性优化的结构 振动极点配置控制系统设计方法,其步骤如下:
[0007] 第一步:根据结构振动有限元方程和期望得到的闭环控制系统的特征值,建立结 构振动极点配置控制系统;
[0008] 第二步:在第一步的基础上进行结构振动极点配置控制系统的特征值分析,利用 区间不确定性分析方法得到结构振动极点配置控制系统特征值的上下界,即1 ,其 中,λ为闭环控制系统的特征值或者称为极点,立,f分别为特征值的上下界;
[0009] 第三步:结构振动极点配置控制系统非概率可靠度计算,经过第一步的结构振动 极点配置控制系统的建立和第二步的结构振动极点配置控制系统的特征值计算,得到基于 极点配置方法的结构振动极点配置控制系统特征值的区间,利用非概率可靠性度量指标的 计算方法,对基于极点配置的主动输出反馈控制系统进行可靠度计算,得到基于极点配置 的主动控制系统非概率可靠度Pos(sys);
[0010] 第四步:对给定的期望特征值进行非概率可靠性优化,在满足非概率可靠度P〇S (sys)的基础上,得到优化后的闭环控制系统极点特征值;优化目标为使得控制器的输出u (t)的最大值max(u(t))最小;
[0011] 第五步:利用非概率优化后的极点期望值进行结构振动极点配置控制系统的设 计。
[0012] 其中,对步骤四进行控制器输出的最大值进行限定,使得基于极点配置的闭环控 制系统在满足可靠度约束下,所需控制力最小。
[0013] 其中,所述第三步中,非概率可靠性度量指标的计算方法:
[0014] 设计要求为闭环控制系统响应不超过某一特定值时Yai,利用如下的计算公式进 行非概率可靠度计算:
[0016] 其中P〇S(SyS)为基于极点配置的闭环控制系统的非概率可靠度,Acri为闭环控制 系统的临界值,即期望设计值d为闭环控制系统特征值的下界,无为闭环控制系统特征值 的上界。
[0017] 其中,所述第四步中对期望的闭环控制系统极点进行非概率可靠性优化时采用的 非概率可靠性优化模型如下:
[0018] f?η?λ
[0019] min I I u(t) I |〇〇
[0020] s.t.Pos(sys)^Rcri
[0021] 其中:λ是极点配置的闭环控制系统的特征值,即设计优化变量;I |u(t) | μ为控制 器输出力范数,P〇S(Sys)为基于极点配置闭环控制系统的非概率可靠度;R cri为设计人员 要求的可靠度,为给定值。
[0022] 其中,为了满足控制器输出最小,给定的设计可靠度RCT1等于1。
[0023] 其中,所述第五步,利用经过非概率可靠性优化后的基于极点配置的控制器设计 得到结构振动主动控制系统过程如下:
[0024] (1)利用第一步给出的极点的期望值设计初始的基于极点配置的结构振动主动控 制系统,在Matlab/Simulink中建立相应的框图;
[0025] (2)利用Matlab中的优化控制箱对基于极点配置的闭环控制系统的特征值进行优 化;
[0026] (3)得到基于极点配置的闭环控制系统控制器参数后设计最终反馈控制系统。
[0027] 本发明与现有技术相比的优点在于:
[0028] (1)本发明将区间不确定性分析引入到结构振动极点配置设计中,通过非概率可 靠性优化得到满足可靠度设计要求的期望极点,利用该极点进行闭环主动控制器设计,使 得结构振动主动控制系统在不确定条件下能够满足可靠性要求。
[0029] (2)本发明提出了主动输出反馈控制系统的非概率可靠性指标计算方法。该方法 可以得到主动输出反馈控制系统在各种不确定性情况下的可靠度指标。
【附图说明】
[0030] 图1为控制前系统极点分布示意图;
[0031] 图2为控制前状态变量X1的响应曲线示意图;
[0032] 图3为控制后系统极点分布示意图;
[0033]图4为控制后状态变量:χι的响应曲线不意图;
[0034]图5为目标函数的迭代曲线示意图;
[0035]图6为控制后状态变量:χι的区间响应曲线不意图;
[0036] 图7为本发明的实现流程图。
【具体实施方式】
[0037] 下面将结合附图对本发明作进一步的详细说明本发明的实施方式。
[0038] 本发明技术解决方案:一种基于非概率可靠性优化的结构振动极点配置控制系统 设计方法,首先,建立结构振动极点配置的主动控制闭环系统。然后,针对建立的闭环输出 反馈主动控制系统,提出了含区间参数的闭环输出反馈主动控制系统特征值的区间估计方 法。在此基础上,建立了闭环输出反馈主动控制系统非概率可靠度计算方法。基于含区间参 数的闭环输出反馈主动控制系统特征值的区间估计方法和非概率可靠度计算方法对输出 反馈控制器进优化,最终得到满足可靠度要求的基于极点配置的输出反馈控制系统。本发 明从可靠性的角度确定输出反馈控制器的参数,从而解决了传统结构振动极点配置输出反 馈控制器设计过程中面临的不确定性问题,使得基于极点配置得到的闭环控制系统在参数 扰动的情况下也能够满足设计要求。如图7所示,其实现步骤如下:
[0039] 第一步:根据结构振动有限元方程和期望得到的闭环控制系统的特征值,建立基 于极点配置的结构振动主动控制系统
[0040] 假定期望得到的结构的刚度矩阵为Ko,质量矩阵为Mo,则由动力学方法可以得到期 望的动力系统为
[0041 ] M0X + K0X = 0: (1)
[0042]其中,S为结构结点的加速度向量,X为结构结点的位移向量。
[0043]已知现有的结构动力系统为
[0045] 其中,M,K分别为开环系统的质量矩阵和刚度矩阵。Bs为控制输出力的定位矩阵,f。 为控制输入向量(控制力),C P为输出定位向量。y为系统的输出。以输出反馈为例,
[0046] fc = -Gy(t) (3)
[0047] 最终,可以求得
[0048] G = Β^Μ(Μ^Κ, - Μ^κ) Cj ' (句
[0049] 其中,G为闭环系统输出增益矩阵,其他变量如上所示。
[0050] 第二步:在第一步的基础上进行结构振动输出反馈控制系统的特征值分析,利用 区间不确定性分析方法得到输出反馈控制系统特征值的上下界,即Ae[丕I],其中,λ为闭 环控制系统的特征值或者称为极点,分别为特征值的上下界;
[0051] 含有不确定性参数的声振耦合系统状态空间表示可以写成如下的形式: X = ? Ac + ΔΑ -e)X+{BC + ΔΒ ·
[0052] 、 , v > 〇) Y=CX+DU
[0053]
和e =
[_1,1]。输出矩阵C和直接转换矩阵D是与不确定参数无关的,在此不做讨论。X为闭环系统 的状态对时间的导数,X为闭环系统的状态,分别为系统矩阵的名义值。e=[-l,l]为单 位区间向量。U为主动控制力,p为结构不确定性参数的个数。
[0054]不失一般性的认为整个系统模型完全可控和完全可观,并采用如下负状态反馈控 制器,
[0055] U(t)=-G · X(t) (6)
[0056] 其中,G为增益矩阵,X为闭环系统的状态,U为主动控制力。将式(6)代入到式(5) 中,可以得到:
[0058] 式(7)中的第一个式子表征了系统的稳定性,第二个式只是输出矩阵,因此,考虑 系统的稳定性,只需关心:
[0059] Χ(/) = (Α1' +ΔΑ·?->-Β1' -G -ΔΒ-C -ejXO) (:8)
[0060] 考虑如下的闭环控制系统
[0061 ] X(/) = A-X(/) (9)
[0062] 其中,X = (Ae + ΔΑ · e-Be G - ΔΒ ·G · e)为一区间矩阵。闭环控制系统特征值的基 本表达式为:
[0063] Kv=之v (10)
[0064] 其中,λ为系统的特征值,ν为特征向量。正如上面提到的,状态矩阵5为一区间矩 阵,是不确定变量b = [bi,b2,. . .bP]T的函数,即:
[0065] A = A(b) (11)
[0066] 要解决的问题是寻找满足式(11)的所有特征值:
[0067] Γ = ·|:/1:Λ e:C,A.(b)v二 1\',b eb1} (12)
[0068] 通常情况下式所示的边界非常复杂,同时要找到所有的特征值也是没有意义的, 但是可以寻找一组区间界限使得下式成立:
[0069] Γ c; '=[立,>1].=(為!),.為r =[心為],f .= .1,2.,…,(13)
[0070]
41为闭环系统特征值区间,η为结构的自由度。为了 得到特征值的区间边界,采用一阶泰勒摄动方法。
[0071 ]考虑如下的标准特征值问题:
[0073] 其中,1为未摄动项,参数的确定性部分。< 是的第i阶右特征向量,μ丨是Α的第 i阶左特征向量,#为闭环系统第i阶特征值的名义值,η为结构的自由度。通常情况下,右特 征向量(i = l,2,. . .,2η可以正交化:
[0074] (V:.) V:. =1. I = 1,2,.,.,2//- (15)
[0075] 其中,η为结构的自由度;
[0076] 现在,给定矩阵七一个较小的摄动量,Μ为未知量,可以得到如下的摄动特 征值问题:
[0077] + (>'入)V〉= Λν;. / 二 1,2,…,2" (16)
[0078] 其中,Ai是系统的第i阶摄动特征值。通常情况下,义和为已知的量。则对应的 特征值h的摄动量为δλ,。将為+M,.代入到方程(16)可以得到:
[0079] Acv^' + δ Αν": - l-cvf + ^Λν^' i = 1,2"." 2/? (17)
[0080] 考虑到式(17)中的第一项,可以将式(17)简化为:
[0081 ] cjAv1; = i - 1,2,..., 2η Π 8)
[0082] 方程(18)两边同时左乘>可以得到:
[0083] =(ν;)! ?^Αν; (19)
[0084] 艮p :
[0085] -^ +(ν;)Τ δλν^ (20)
[0086] 由于主动控制系统增加了控制器,破坏了原系统矩阵的对称性,因此矩阵Α的特 征值通常情况下为复数,公式(20)可以将特征值的实部和虚部分开写:
[0088] 其中,Air和λ?γ分别是特征值的实部和虚部,2丨和為〖.分别是系统特征值实部和虚部 的名义值。别是特征值的实部和虚部的摄动量。考虑到正交性可以得到:
[0089] (v;)1 t)Av; =Λ<; +ν-??^, (22)
[0090] 考虑特征向量的实部和虚部,可以得到:
[0092] 展开方程(23),分别考虑实部和虚部可以得到:
[0093] V , ^ ;τ ^ (24) <Η>=(ν·.) <5Αν;+(ν;) ?5Αν;
[0094] 首先,针对特征值的实部利用区间扩张原理可以得到特征值实部的区间为:
[0095] AJ.=Ay.+<SAl (25)
[0096] 其中,
[0098] Δ1为区间矩阵孟的半径。利用区间算法,并令Λλ"为特征值\实部的半径,则:
[0099] = (27)
[0100] 其中,4. = < -Δ4Λ. = 4 +ΔΛ ' Δ4. = |ν;.|τ么五卜';;.| + |<.「δΧ|^|。同样,可以 得到特征值Μ,i = 1,2,. . .,2η的虚部的区间:
[0101] ,ζ'.二 1-,,2"巧.2样 (28)
[0102] 其中,4 = 4 & = g + M,.,= |<「"K.I+1<.「"I。可以得到区 间矩阵S的区间特征值:
[0106] 其中,4和^ir分别为第i个特征值实部的上下界,疋.和&V分别为第i个特征值虚部 的上下界。<和2,丨分别为第i个特征值实部和虚部的名义值。%和 <.分别为第i个特征向量 实部和虚部的名义值。
[0107] 第三步:结构振动输出反馈控制系统非概率可靠度计算,经过第一步的结构振动 输出反馈控制系统的建立和第二步的输出反馈主动控制系统的特征值计算,得到基于极点 配置方法的闭环输出反馈控制系统特征值的区间,利用非概率可靠性度量指标的计算方 法,对基于极点配置的主动输出反馈控制系统进行可靠度计算,得到基于极点配置的主动 控制系统非概率可靠度Pos(sys):
[0109] 其中,P〇S(SyS)为基于极点配置的闭环控制系统的非概率可靠度,A cri为闭环控制 系统的临界值,即期望设计值为闭环控制系统特征值的下界,I为闭环控制系统特征值 的上界。
[0110] 第四步:所述第四步中对期望的闭环控制系统极点进行非概率可靠性优化时采用 的非概率可靠性优化模型如下:
[0111] findA
[0112] min | | u(t) | |〇〇
[0113] s. t.Pos(sys)^Rcri
[0114]其中,λ是极点配置的闭环控制系统的特征值,即设计优化变量;I |u(t)I μ为控制 器输出力范数,P〇S(SyS)为基于极点配置闭环控制系统的非概率可靠度;R cri为设计人员 要求的可靠度,为给定值。find为寻找设计变量,min为最小优化目标,s. t.为要满足的约束 条件。
[0115] 第五步:利用非概率优化后的极点期望值进行结构振动极点配置控制系统的设 计。
[0116] (1)利用第一步给出的极点的期望值设计初始的基于极点配置的结构振动主动控 制系统,在Matlab/Simulink中建立相应的框图;
[0117] (2)利用Matlab中的优化控制箱对基于极点配置的闭环控制系统的特征值进行优 化;
[0118] (3)得到基于极点配置的闭环控制系统控制器参数后设计最终反馈控制系统。
[0119] 本发明具体实施例介绍如下:
[0120]考虑具有以下状态空间形式的多输入多输出动力系统:
[0121] x=Ax+Bu
[0122] y = Cx+Du
[0123] 其中,x为系统状态变量,y为系统输出变量,u为控制输入;
[0125] 系统的不确定量为AA = abs(0.01XA),假定控制输入矩阵B不存在不确定性。给 定初始条件为Condition=[l 0 0 0 0 0]τ,利用极点配置方法实现对动力响应输出的抑 制。由于系统为多输入多输出系统,本算例只关心第一个自由度的变化,首先可以计算出系 统的特征值为:3,0.8105,0.5± 1.3229 j,-0.9053 ± 1.2837 j,如图1所示。可以发现系统的 特征值存在实部大于零的情况,因此,该系统为不稳定系统。在给定的初始条件下,系统输 出响应随着时间的推进会变为无穷大,如图2所示。为了达到抑制系统响应的目的,期望系 统的极点能够出现在以下几个位置:-1,-2,_3,-4,-1 土 j。
[0126] 第一步:根据结构振动有限元方程和期望得到的闭环控制系统的特征值,建立基 于极点配置的结构振动主动控制系统。
[0127] 假定期望得到的结构的刚度矩阵为Ko,质量矩阵为Mo,则由动力学方法可以得到期 望的动力系统为:
[0128] Μ0Χ + Κ0Χ = 0 (1)
[0129] 已知现有的结构动力系统为:
[0131]其中,Bs为控制输出力的定位矩阵,f。为控制输入向量(控制力)。以输出反馈为例:
[0132] fc = -Gy(t) (3)
[0133] 最终,可以求得:
[0135]根据极点配置方法可以得到状态反馈矩阵G如下,并且满足以上的极点要求如图3 所示。利用得到的反馈矩阵G构建闭环控制系统,得到的状态变量^的响应曲线如图4所示。

[0137] 第二步:在第一步的基础上进行结构振动输出反馈控制系统的特征值分析,利用 区间不确定性分析方法得到输出反馈控制系统特征值的上下界,即2 e[毛I]。
[0138] 含有不确定性参数的声振耦合系统状态空间表示可以写成如下的形式:
[0140] 其中
和e =
[_1,1]。输出矩阵C和直接转换矩阵D是与不确定参数无关的,在此不做讨论。不失一般性的 认为整个系统模型完全可控和完全可观,并采用如下负状态反馈控制器:
[0141] U(t)=-G · X(t) (6)
[0142] 其中,G为增益矩阵,将式(6)代入到式(5)中,可以得到:
[0144] 式(7)中的第一个式子表征了系统的稳定性,第二个式只是输出矩阵,因此,考虑 系统的稳定性,只需关心:
[0145] X(/) = (Α1·' + ΔΑ · e -BL' ·G --ΔΒ ·G · ? X(t) (8)
[0146] 考虑如下的闭环控制系统:
[0147] Χ(〇-Α·Χ(?) (9)
[0148] 其中,i = (Ae + ΔΑ · e - :Τ · G - ΔΒ · G · e)为一区间矩阵。闭环控制系统特征值的基 本表达式为:
[0149] .Αν = (10)
[0150] 其中,λ为系统的特征值,ν为特征向量。正如上面提到的,状态矩阵1为一区间矩 阵,是不确定变量b = [bi,b2,. . .bP]T的函数,即:
[0151] A - A(b) (il)
[0152] 要解决的问题是寻找满足式(11)的所有特征值:
[0153] = = (12)
[0154] 通常情况下式所示的边界非常复杂,同时要找到所有的特征值也是没有意义的, 但是可以寻找一组区间界限使得下式成立:
[0156] 其中:
。为了得到特征值的区间边界,采用一阶泰勒摄动方
法。
[0157] 考虑如下的标准特征值问题:
[0159] 其中,f为未摄动项,参数的确定性部分。v丨是士的第i阶右特征向量,μ〖是1。的第 i阶左特征向量。通常情况下,右特征向量= l,2,. . .,2η可以正交化:
[0160] (V:.) V:. =1. I = 1,2,.,.,2//- (15)
[0161] 现在,给定矩阵玉&一个较小的摄动量为未知量,可以得到如下的摄动特 征值问题:
[0162] (A + ) vf = i = 1,.;2:,,.., 2η (10)
[0163] 其中,Ai是系统的第i阶摄动特征值。通常情况下,为已知的量。则对应的 特征值h的摄动量为δλ,。将2, =< +5為代入到方程(16)可以得到:
[0164] AL V:. + c)'Av;.:=々.V:.十<rH;v:: / -- 1,2,...,2/7 (17)
[0165] 考虑到式(17)中的第一项,可以将式(17)简化为:
[0166] SAv] --?λ ν] i = 1,2,..., 2/? (18)
[0167] 方程(18)两边同时左乘,可以得到:
[0168] al^(vf)' tfAv; {\J)
[0169] 艮p:
[0170] 2. = 2C + δλ = A; + (v,c )T ?'Αν; (20)
[0171] 由于主动控制系统增加了控制器,破坏了原系统矩阵的对称性,因此矩阵A的特 征值通常情况下为复数,公式(20)可以将特征值的实部和虚部分开写:
[0173]其中,Air和Aiy分别是特征值的实部和虚部。δλι4Ρδλ?γ分别是特征值的实部和虚部 的摄动量。考虑到正交性可以得到:
[0175]考虑特征向量的实部和虚部,可以得到:
[0177]展开方程(23),分别考虑实部和虚部可以得到:
[0179]首先,针对特征值的实部利用区间扩张原理可以得到特征值实部的区间为:
[0180] ^=?^ + β^ (25)
[0181] 其中,
[0183] Δ?为区间矩阵5的半径。利用区间算法,并令Λλ"为特征值\实部的半径,则:
[0184] λ? = [ X., i.r ], / = 1,2,..., 2/? (27)
[0185] 其中,l = 2,r - 4 = < + Δ4,Mr = |( I AS |v;rj + |v;I Δ斗I。同样,可以 得到特征值Μ,i = 1,2,. . .,2n的虚部的区间:
[0186] Λν ~ [iy5^ 2,..., 2? (2.8)
[0187] 其中,4 =Λ;+Δ. Δ4 =卜;;.|τδ1|4| + |<|τδΖ|<|。可以得到区 间矩阵X的区间特征值:
[0191 ]最终得到闭环控制系统的极点如下表1所示
[0192] 表1传统设计方法的闭环控制系统极点及系统可靠度
[0193]
[0194] 第三步:结构振动输出反馈控制系统非概率可靠度计算,经过第一步的结构振动 输出反馈控制系统的建立和第二步的输出反馈主动控制系统的特征值计算,得到基于极点 配置方法的闭环输出反馈控制系统特征值的区间,利用非概率可靠性度量指标的计算方 法,对基于极点配置的主动输出反馈控制系统进行可靠度计算,得到基于极点配置的主动 控制系统非概率可靠度Pos(sys):
[0196] 其中,P〇S(SyS)为基于极点配置的闭环控制系统的非概率可靠度,Acri为闭环控制 系统的临界值,即期望设计值为闭环控制系统特征值的下界,X为闭环控制系统特征值 的上界。此时,有传统方法得到的闭环系统的可靠度以及在表1中给出。
[0197] 第四步:所述第四步中对期望的闭环控制系统极点进行非概率可靠性优化时采用 的非概率可靠性优化模型如下:
[0198] f?η?λ
[0199] min I I u(t) I |〇〇
[0200] s.t.Pos(sys)^Rcri
[0201] 其中,λ是极点配置的闭环控制系统的特征值,即设计优化变量;I |u(t)| μ为控制 器输出力范数,P〇S(Sys)为基于极点配置闭环控制系统的非概率可靠度;R cri为设计人员 要求的可靠度,为一给定的值。
[0202] 从表1中可以看出,由于不确定性的存在,闭环控制系统极点的实部已经有大于期 望值的可能性,这样利用本发明提出的非概率可靠极点配置方法,将闭环控制系统的极点 实部最大值配置到期望值处,即将系统的可靠度设计为1。经过迭代以后,图5给出了相应的 优化迭代曲线,得到的控制器如下:
[0204]表2本发明设计的闭环控制系统极点与可靠度
[0205]
[0206] 第五步:利用非概率优化后的极点期望值进行结构振动极点配置控制系统的设 计。
[0207] (1)利用第一步给出的极点的期望值设计初始的基于极点配置的结构振动主动控 制系统,在Matlab/Simulink中建立相应的框图;
[0208] (2)利用Matlab中的优化控制箱对基于极点配置的闭环控制系统的特征值进行优 化;
[0209] (3)得到基于极点配置的闭环控制系统控制器参数后设计最终反馈控制系统。
[0210] 图6给出了基于非概率可靠优化的结构振动极点配置控制系统的最终响应输出。 表3给出了传统极点配置方法和本发明所提出的非概率可靠极点配置方法的响应输出的比 较,可以看出,本发明所提出的方法满足可靠度要求。
[0211]表3两种方法控制效果比较 [0212]
【主权项】
1. 一种基于非概率可靠性优化的结构振动极点配置控制系统设计方法,其特征在于步 骤如下: 第一步:根据结构振动有限元方程和期望得到的闭环控制系统的特征值,建立结构振 动极点配置控制系统; 第二步:在第一步的基础上进行结构振动极点配置控制系统的特征值分析,利用区间 不确定性分析方法得到结构振动极点配置控制系统特征值的上下界,即^ e [么,其中,λ 为闭环控制系统的特征值或者称为极点,全,I分别为特征值的上下界; 第Ξ步:结构振动极点配置控制系统非概率可靠度计算,经过第一步的结构振动极点 配置控制系统的建立和第二步的结构振动极点配置控制系统的特征值计算,得到基于极点 配置方法的结构振动极点配置控制系统特征值的区间,利用非概率可靠性度量指标的计算 方法,对基于极点配置的主动输出反馈控制系统进行可靠度计算,得到基于极点配置的主 动控制系统非概率可靠度化s(sys); 第四步:对给定的期望特征值进行非概率可靠性优化,在满足非概率可靠度Pos(sys) 的基础上,得到优化后的闭环控制系统极点特征值;优化目标为使得控制器的输出u(t)的 最大值max(u(t))最小; 第五步:利用非概率优化后的极点期望值进行结构振动极点配置控制系统的设计。2. 根据权利要求1所述的基于非概率可靠性优化的结构振动极点配置控制系统设计方 法,其特征在于:对步骤四进行控制器输出的最大值进行限定,使得基于极点配置的闭环控 制系统在满足可靠度约束下,所需控制力最小。3. 根据权利要求1所述的基于非概率可靠性优化的结构振动极点配置控制系统设计方 法,其特征在于:所述第Ξ步中,非概率可靠性度量指标的计算方法: 设计要求为闭环控制系统响应不超过某一特定值时Υαι,利用如下的计算公式进行非 概率可靠度计算:其中Pos(sys)为基于极点配置的闭环控制系统的非概率可靠度,λ"ι为闭环控制系统 的临界值,即期望设计值,全为闭环控制系统特征值的下界,J为闭环控制系统特征值的上 界。4. 根据权利要求1所述的基于非概率可靠性优化的结构振动极点配置控制系统设计方 法,其特征在于:所述第四步中对期望的闭环控制系统极点进行非概率可靠性优化时采用 的非概率可靠性优化模型如下: f ?η?λ mini |u(t) I |〇〇 s.t.F*os(sys)>Rcri 其中:λ是极点配置的闭环控制系统的特征值,即设计优化变量;||u(t) I μ为控制器输 出力He?范数,Pos(sys)为基于极点配置闭环控制系统的非概率可靠度;Rcri为设计人员要求 的可靠度,为给定值。5. 根据权利要求4所述的基于非概率可靠性优化的结构振动极点配置控制系统设计方 法,其特征在于:为了满足控制器输出最小,给定的设计可靠度Rai等于1。6. 根据权利要求1所述的基于非概率可靠性优化的结构振动极点配置控制系统设计方 法,其特征在于:所述第五步,利用经过非概率可靠性优化后的基于极点配置的控制器设计 得到结构振动主动控制系统过程如下: (1) 利用第一步给出的极点的期望值设计初始的基于极点配置的结构振动主动控制系 统,在Matlab/Simulink中建立相应的框图; (2) 利用Matlab中的优化控制箱对基于极点配置的闭环控制系统的特征值进行优化; (3) 得到基于极点配置的闭环控制系统控制器参数后设计最终反馈控制系统。
【文档编号】G05B13/04GK106094518SQ201610457443
【公开日】2016年11月9日
【申请日】2016年6月22日
【发明人】王晓军, 李云龙, 邱志平, 许孟辉, 王磊, 仇翯辰, 陈贤佳, 郑宇宁
【申请人】北京航空航天大学
网友询问留言 已有0条留言
  • 还没有人留言评论。精彩留言会获得点赞!
1