CN106094518B - 一种基于非概率可靠性优化的结构振动极点配置控制方法 - Google Patents
一种基于非概率可靠性优化的结构振动极点配置控制方法 Download PDFInfo
- Publication number
- CN106094518B CN106094518B CN201610457443.2A CN201610457443A CN106094518B CN 106094518 B CN106094518 B CN 106094518B CN 201610457443 A CN201610457443 A CN 201610457443A CN 106094518 B CN106094518 B CN 106094518B
- Authority
- CN
- China
- Prior art keywords
- control system
- pole
- closed
- reliability
- structural vibration
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 66
- 238000005457 optimization Methods 0.000 title claims abstract description 43
- 238000013461 design Methods 0.000 claims abstract description 29
- 238000004364 calculation method Methods 0.000 claims abstract description 14
- 230000004044 response Effects 0.000 claims description 13
- 238000010586 diagram Methods 0.000 claims description 9
- 238000013076 uncertainty analysis Methods 0.000 claims description 5
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 claims description 4
- 238000005259 measurement Methods 0.000 claims description 4
- 230000008569 process Effects 0.000 claims description 4
- 238000012938 design process Methods 0.000 abstract description 2
- 239000011159 matrix material Substances 0.000 description 41
- 238000011160 research Methods 0.000 description 5
- 241000764238 Isis Species 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000001629 suppression Effects 0.000 description 2
- 206010063385 Intellectualisation Diseases 0.000 description 1
- 241000053208 Porcellio laevis Species 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 238000011217 control strategy Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000006880 cross-coupling reaction Methods 0.000 description 1
- 238000013016 damping Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000011426 transformation method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
- G05B13/042—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Feedback Control In General (AREA)
Abstract
本发明涉及一种基于非概率可靠性优化的结构振动极点配置控制方法,该方法首先建立结构振动极点配置的主动控制闭环系统,然后建立了闭环输出反馈主动控制系统非概率可靠度计算方法,基于含区间参数的闭环输出反馈主动控制系统特征值的区间估计方法和非概率可靠度计算方法对输出反馈控制器进优化,最终得到满足可靠度要求的基于极点配置的输出反馈控制系统。本发明从可靠性的角度确定输出反馈控制器的参数,从而解决了传统结构振动极点配置输出反馈控制器设计过程中面临的不确定性问题,使得基于极点配置得到的闭环控制系统在参数扰动的情况下也能够满足设计要求。
Description
技术领域
本发明涉及结构振动主动控制的技术领域,具体涉及一种基于非概率可靠性优化的结构振动极点配置控制方法。
背景技术
振动与噪声的主动控制是一个快速发展的新领域。振动与噪声的主动控制技术涉及多学科交叉耦合,具有广泛的应用价值和学术价值,是航空航天结构走向智能化的开始。振动与噪声主动控制就是通过引入次级干扰(振动或者噪声)来消除初级干扰形成的振动或噪声。次级干扰是通过驱动器产生,次级干扰的具体形式由控制器经过计算得出。振动与噪声主动控制技术的发展不是一蹴而就的,也经过了一定时期的积累。现有的控制器设计方法大都是根据控制理论提出的方法进行的,主要包括极点配置、最优控制、PID控制和鲁棒控制等理论。当然,随着现代控制理论的发展,结构振动与噪声主动控制系统控制器的设计也是在前进和发展的。
极点配置方法是一种非常简单的反馈控制器设计方法。这种设计方法的目的就是保证闭环控制系统矩阵的特征值与设计者要求的特征值接近。通常来说,设计者可以把闭环系统矩阵的特征值设计到任何位置,这在理论上可行的。但是,在考虑各种其他物理因素的情况下,实际工程中并不这样做。例如,可以通过设计控制器使得闭环系统的特征值更加靠近复平面的左边,这样就可以提高系统稳定的鲁棒性。Manning等利用系统识别和极点配置控制方法实现了柔性智能梁第一阶模态的振动主动控制,结果表明极点配置方法可以媲美速度反馈控制器。Sethi等针对一个3.5m长的复合材料工字型梁,利用压电材料作为传感器和驱动器设计了极点配置控制器,结果表明基于极点配置设计的控制器充当了阻尼器的作用。不仅保持了原来系统的稳定性,同时还增加了系统对外部干扰的鲁棒性。Kumar等对智能结构的振动主动控制系统的极点配置可行性进行了研究,考虑了结构参数的不确定性,提出了鲁棒极点配置方法,该方法不仅可以保证系统的鲁棒稳定性,而且还能保证系统的鲁棒性能。Hanagan 和Murray等研究了基于极点配置的结构振动抑制方法,通过施加速度反馈达到对闭环系统极点进行调节的目的,从而实现结构振动响应的控制。JUNKINS等基于最小灵敏度设计方法得到了极点配置的输出反馈控制器,实现了振动响应的主动控制。Slater等人利用约束状态和输出反馈控制设计了基于特征值配置的主动控制系统,并将其应用于柔性结构的振动控制。Bittanti等利用简单的极点配置方法提出了一种新的适用于直升机的叶片振动控制的方法,该方法具有两个优点:一、可以不影响整个直升机的飞行动力学;二、控制器的设计和计算简单有效。Ram研究了多输入系统的极点配置方法,单输入系统不能解决闭环控制系统特征值相关性的问题,没法保证在其他特征值不变的情况下实现要求的特征值变化,Ram通过多输入控制方法对特征值相关的极点配置问题进行了研究。Abdelaziz等利用状态微分反馈控制器研究了极点配置在单输入单输出线性系统中的应用,其求解过程是与Ackermann 公式相似的,同时该方法还可以应用于结构振动控制时变系统。路小波等利用试验对系统进行辨识,设计了基于极点配置的柔性智能结构主动振动控制器,并对仿真结果进行了试验。李书等针对含不确定参数的结构系统,提出了应用Householder变换的方法设计具有鲁棒性的极点配置控制器。张家凡针对具有非对称阻尼和非对称刚度矩阵的一般动力学系统,研究了多输入极点配置问题,给出了控制增益矩阵的显式解。周星德等提出反向极点配置,该控制策略能够实现部分极点的重新配置,并且保证了不会出现控制输出力较大的情况,具有很高的工程实用价值。
虽然基于极点配置的闭环控制系统的设计方法已被广泛用于振动控制领域,但是当系统参数存在不确定性时,利用名义系统得到的闭环控制系统有时候并不能满足设计的可靠度要求。本发明就是从非概率可靠性的角度出发,基于非概率可靠性优化提出的一种可以考虑模型不确定性的极点配置振动主动控制方法,使得到的闭环控制系统在满足可靠性要求的前提下,达到要求的振动控制效果。现有专利文献和非专利文献均无相关技术的报道。
发明内容
本发明要解决的技术问题是:克服现有极点配置方法无法处理模型存在不确定性的缺点,提供一种基于非概率可靠性优化的结构振动极点配置主动振动控制方法,从而提高极点配置主动控制系统的可靠性。
本发明解决上述技术问题采用的技术方案为:一种基于非概率可靠性优化的结构振动极点配置控制方法,其步骤如下:
第一步:根据结构振动有限元方程和期望得到的闭环控制系统的特征值,建立结构振动极点配置控制系统;
第二步:在第一步的基础上进行结构振动极点配置控制系统的特征值分析,利用区间不确定性分析方法得到结构振动极点配置控制系统特征值的上下界,即其中,λ为闭环控制系统的特征值或者称为极点,λ,分别为特征值的上下界;
第三步:结构振动极点配置控制系统非概率可靠度计算,经过第一步的结构振动极点配置控制系统的建立和第二步的结构振动极点配置控制系统的特征值计算,得到基于极点配置方法的结构振动极点配置控制系统特征值的区间,利用非概率可靠性度量指标的计算方法,对基于极点配置的主动输出反馈控制系统进行可靠度计算,得到基于极点配置的主动控制系统非概率可靠度Pos(sys);
第四步:对给定的期望特征值进行非概率可靠性优化,在满足非概率可靠度Pos(sys)的基础上,得到优化后的闭环控制系统极点特征值;优化目标为使得控制器的输出u(t)的最大值max(u(t))最小;
第五步:利用非概率优化后的极点特征值进行结构振动极点配置控制系统的设计。
其中,对步骤四进行控制器输出的最大值进行限定,使得基于极点配置的闭环控制系统在满足可靠度约束下,所需控制力最小。
其中,所述第三步中,非概率可靠性度量指标的计算方法:
设计要求为闭环控制系统响应不超过某一特定值Ycri时,利用如下的计算公式进行非概率可靠度计算:
其中Pos(sys)为基于极点配置的闭环控制系统的非概率可靠度,λcri为闭环控制系统的临界值,即期望设计值,λ为闭环控制系统特征值的下界,为闭环控制系统特征值的上界。
其中,所述第四步中对期望的闭环控制系统极点进行非概率可靠性优化时采用的非概率可靠性优化模型如下:
findλ
min||u(t)||∞
s.t.Pos(sys)≥Rcri
其中:λ是极点配置的闭环控制系统的特征值,即设计优化变量;||u(t)||∞为控制器输出力H∞范数,Pos(sys)为基于极点配置闭环控制系统的非概率可靠度;Rcri为设计人员要求的可靠度,为给定值。
其中,为了满足控制器输出最小,给定的设计可靠度Rcri等于1。
其中,所述第五步,利用经过非概率可靠性优化后的基于极点配置的控制器设计得到结构振动主动控制系统过程如下:
(1)利用第一步给出的极点的期望值设计初始的基于极点配置的结构振动主动控制系统,在Matlab/Simulink中建立相应的框图;
(2)利用Matlab中的优化控制箱对基于极点配置的闭环控制系统的特征值进行优化;
(3)得到基于极点配置的闭环控制系统控制器参数后设计最终反馈控制系统。
本发明与现有技术相比的优点在于:
(1)本发明将区间不确定性分析引入到结构振动极点配置设计中,通过非概率可靠性优化得到满足可靠度设计要求的期望极点,利用该极点进行闭环主动控制器设计,使得结构振动主动控制系统在不确定条件下能够满足可靠性要求。
(2)本发明提出了主动输出反馈控制系统的非概率可靠性指标计算方法。该方法可以得到主动输出反馈控制系统在各种不确定性情况下的可靠度指标。
附图说明
图1为控制前系统极点分布示意图;
图2为控制前状态变量x1的响应曲线示意图;
图3为控制后系统极点分布示意图;
图4为控制后状态变量x1的响应曲线示意图;
图5为目标函数的迭代曲线示意图;
图6为控制后状态变量x1的区间响应曲线示意图;
图7为本发明的实现流程图。
具体实施方式
下面将结合附图对本发明作进一步的详细说明本发明的实施方式。
本发明技术解决方案:一种基于非概率可靠性优化的结构振动极点配置控制方法,首先,建立结构振动极点配置的主动控制闭环系统。然后,针对建立的闭环输出反馈主动控制系统,提出了含区间参数的闭环输出反馈主动控制系统特征值的区间估计方法。在此基础上,建立了闭环输出反馈主动控制系统非概率可靠度计算方法。基于含区间参数的闭环输出反馈主动控制系统特征值的区间估计方法和非概率可靠度计算方法对输出反馈控制器进优化,最终得到满足可靠度要求的基于极点配置的输出反馈控制系统。本发明从可靠性的角度确定输出反馈控制器的参数,从而解决了传统结构振动极点配置输出反馈控制器设计过程中面临的不确定性问题,使得基于极点配置得到的闭环控制系统在参数扰动的情况下也能够满足设计要求。如图7所示,其实现步骤如下:
第一步:根据结构振动有限元方程和期望得到的闭环控制系统的特征值,建立基于极点配置的结构振动主动控制系统
假定期望得到的结构的刚度矩阵为K0,质量矩阵为M0,则由动力学方法可以得到期望的动力系统为
其中,为结构结点的加速度向量,x为结构结点的位移向量。
已知现有的结构动力系统为
其中,M,K分别为开环系统的质量矩阵和刚度矩阵。Bs为控制输出力的定位矩阵,fc为控制输入向量(控制力),Cp为输出定位向量。y为系统的输出。以输出反馈为例,
fc=-Gy(t) (3)
最终,可以求得
其中,G为闭环系统输出增益矩阵,其他变量如上所示。
第二步:在第一步的基础上进行结构振动输出反馈控制系统的特征值分析,利用区间不确定性分析方法得到输出反馈控制系统特征值的上下界,即其中,λ为闭环控制系统的特征值或者称为极点,λ,分别为特征值的上下界;
含有不确定性参数的声振耦合系统状态空间表示可以写成如下的形式:
其中,Ac=A(bc),Bc=B(bc),和e=[-1,1]。输出矩阵C和直接转换矩阵D是与不确定参数无关的,在此不做讨论。为闭环系统的状态对时间的导数,X为闭环系统的状态,Ac、Bc分别为系统矩阵的名义值。e=[-1,1]为单位区间向量。U为主动控制力,p为结构不确定性参数的个数。
不失一般性的认为整个系统模型完全可控和完全可观,并采用如下负状态反馈控制器,
U(t)=-G·X(t) (6)
其中,G为增益矩阵,X为闭环系统的状态,U为主动控制力。将式(6)代入到式(5)中,可以得到:
式(7)中的第一个式子表征了系统的稳定性,第二个式只是输出矩阵,因此,考虑系统的稳定性,只需关心:
考虑如下的闭环控制系统
其中,为一区间矩阵。闭环控制系统特征值的基本表达式为:
其中,λ为系统的特征值,ν为特征向量。正如上面提到的,状态矩阵为一区间矩阵,是不确定变量b=[b1,b2,...bp]T的函数,即:
要解决的问题是寻找满足式(11)的所有特征值:
通常情况下式所示的边界非常复杂,同时要找到所有的特征值也是没有意义的,但是可以寻找一组区间界限使得下式成立:
其中,λI为闭环系统特征值区间,n为结构的自由度。为了得到特征值的区间边界,采用一阶泰勒摄动方法。
考虑如下的标准特征值问题:
其中,为未摄动项,参数的确定性部分。是的第i阶右特征向量,是的第i阶左特征向量,为闭环系统第i阶特征值的名义值,n为结构的自由度。通常情况下,右特征向量可以正交化:
其中,n为结构的自由度;
现在,给定矩阵一个较小的摄动量为未知量,可以得到如下的摄动特征值问题:
其中,λi是系统的第i阶摄动特征值。通常情况下,和为已知的量。则对应的特征值λi的摄动量为δλi。将代入到方程(16)可以得到:
考虑到式(17)中的第一项,可以将式(17)简化为:
方程(18)两边同时左乘可以得到:
即:
由于主动控制系统增加了控制器,破坏了原系统矩阵的对称性,因此矩阵的特征值通常情况下为复数,公式(20)可以将特征值的实部和虚部分开写:
其中,λir和λiy分别是特征值的实部和虚部,和分别是系统特征值实部和虚部的名义值。
δλir和δλiy分别是特征值的实部和虚部的摄动量。考虑到正交性可以得到:
考虑特征向量的实部和虚部,可以得到:
展开方程(23),分别考虑实部和虚部可以得到:
首先,针对特征值的实部利用区间扩张原理可以得到特征值实部的区间为:
其中,
为区间矩阵的半径。利用区间算法,并令Δλir为特征值λi实部的半径,则:
其中,同样,可以得到特征值λi,i=1,2,...,2n的虚部的区间:
其中,可以得到区间矩阵的区间特征值:
其中,
其中,和λir分别为第i个特征值实部的上下界,和λiy分别为第i个特征值虚部的上下界。和分别为第i个特征值实部和虚部的名义值。和分别为第i个特征向量实部和虚部的名义值。
第三步:结构振动输出反馈控制系统非概率可靠度计算,经过第一步的结构振动输出反馈控制系统的建立和第二步的输出反馈主动控制系统的特征值计算,得到基于极点配置方法的闭环输出反馈控制系统特征值的区间,利用非概率可靠性度量指标的计算方法,对基于极点配置的主动输出反馈控制系统进行可靠度计算,得到基于极点配置的主动控制系统非概率可靠度Pos(sys):
其中,Pos(sys)为基于极点配置的闭环控制系统的非概率可靠度,λcri为闭环控制系统的临界值,即期望设计值。λ为闭环控制系统特征值的下界,为闭环控制系统特征值的上界。
第四步:所述第四步中对期望的闭环控制系统极点进行非概率可靠性优化时采用的非概率可靠性优化模型如下:
findλ
min||u(t)||∞
s.t.Pos(sys)≥Rcri
其中,λ是极点配置的闭环控制系统的特征值,即设计优化变量;||u(t)||∞为控制器输出力H∞范数,Pos(sys)为基于极点配置闭环控制系统的非概率可靠度;Rcri为设计人员要求的可靠度,为给定值。find为寻找设计变量,min为最小优化目标,s.t.为要满足的约束条件。
第五步:利用非概率优化后的极点特征值进行结构振动极点配置控制系统的设计。
(1)利用第一步给出的极点的期望值设计初始的基于极点配置的结构振动主动控制系统,在Matlab/Simulink中建立相应的框图;
(2)利用Matlab中的优化控制箱对基于极点配置的闭环控制系统的特征值进行优化;
(3)得到基于极点配置的闭环控制系统控制器参数后设计最终反馈控制系统。
本发明具体实施例介绍如下:
考虑具有以下状态空间形式的多输入多输出动力系统:
x=Ax+Bu
y=Cx+Du
其中,x为系统状态变量,y为系统输出变量,u为控制输入;
系统的不确定量为ΔA=abs(0.01×A),假定控制输入矩阵B不存在不确定性。给定初始条件为Condition=[1 0 0 0 0 0]T,利用极点配置方法实现对动力响应输出的抑制。由于系统为多输入多输出系统,本算例只关心第一个自由度的变化,首先可以计算出系统的特征值为:3,0.8105,0.5±1.3229j,-0.9053±1.2837j,如图1所示。可以发现系统的特征值存在实部大于零的情况,因此,该系统为不稳定系统。在给定的初始条件下,系统输出响应随着时间的推进会变为无穷大,如图2所示。为了达到抑制系统响应的目的,期望系统的极点能够出现在以下几个位置:-1,-2,-3,-4,-1±j。
第一步:根据结构振动有限元方程和期望得到的闭环控制系统的特征值,建立基于极点配置的结构振动主动控制系统。
假定期望得到的结构的刚度矩阵为K0,质量矩阵为M0,则由动力学方法可以得到期望的动力系统为:
已知现有的结构动力系统为:
其中,Bs为控制输出力的定位矩阵,fc为控制输入向量(控制力)。以输出反馈为例:
fc=-Gy(t) (3)
最终,可以求得:
根据极点配置方法可以得到状态反馈矩阵G如下,并且满足以上的极点要求如图3所示。利用得到的反馈矩阵G构建闭环控制系统,得到的状态变量x1的响应曲线如图4所示。
第二步:在第一步的基础上进行结构振动输出反馈控制系统的特征值分析,利用区间不确定性分析方法得到输出反馈控制系统特征值的上下界,即
含有不确定性参数的声振耦合系统状态空间表示可以写成如下的形式:
其中,Ac=A(bc),Bc=B(bc),和e=[-1,1]。输出矩阵C和直接转换矩阵D是与不确定参数无关的,在此不做讨论。不失一般性的认为整个系统模型完全可控和完全可观,并采用如下负状态反馈控制器:
U(t)=-G·X(t) (6)
其中,G为增益矩阵,将式(6)代入到式(5)中,可以得到:
式(7)中的第一个式子表征了系统的稳定性,第二个式只是输出矩阵,因此,考虑系统的稳定性,只需关心:
考虑如下的闭环控制系统:
其中,为一区间矩阵。闭环控制系统特征值的基本表达式为:
其中,λ为系统的特征值,ν为特征向量。正如上面提到的,状态矩阵为一区间矩阵,是不确定变量b=[b1,b2,...bp]T的函数,即:
要解决的问题是寻找满足式(11)的所有特征值:
通常情况下式所示的边界非常复杂,同时要找到所有的特征值也是没有意义的,但是可以寻找一组区间界限使得下式成立:
其中,为了得到特征值的区间边界,采用一阶泰勒摄动方法。
考虑如下的标准特征值问题:
其中,为未摄动项,参数的确定性部分。是的第i阶右特征向量,是的第i阶左特征向量。通常情况下,右特征向量可以正交化:
现在,给定矩阵一个较小的摄动量为未知量,可以得到如下的摄动特征值问题:
其中,λi是系统的第i阶摄动特征值。通常情况下,和为已知的量。则对应的特征值λi的摄动量为δλi。将代入到方程(16)可以得到:
考虑到式(17)中的第一项,可以将式(17)简化为:
方程(18)两边同时左乘可以得到:
即:
由于主动控制系统增加了控制器,破坏了原系统矩阵的对称性,因此矩阵的特征值通常情况下为复数,公式(20)可以将特征值的实部和虚部分开写:
其中,λir和λiy分别是特征值的实部和虚部。δλir和δλiy分别是特征值的实部和虚部的摄动量。考虑到正交性可以得到:
考虑特征向量的实部和虚部,可以得到:
展开方程(23),分别考虑实部和虚部可以得到:
首先,针对特征值的实部利用区间扩张原理可以得到特征值实部的区间为:
其中,
为区间矩阵的半径。利用区间算法,并令Δλir为特征值λi实部的半径,则:
其中,同样,可以得到特征值λi,i=1,2,...,2n的虚部的区间:
其中,可以得到区间矩阵的区间特征值:
其中,
最终得到闭环控制系统的极点如下表1所示
表1传统设计方法的闭环控制系统极点及系统可靠度
第三步:结构振动输出反馈控制系统非概率可靠度计算,经过第一步的结构振动输出反馈控制系统的建立和第二步的输出反馈主动控制系统的特征值计算,得到基于极点配置方法的闭环输出反馈控制系统特征值的区间,利用非概率可靠性度量指标的计算方法,对基于极点配置的主动输出反馈控制系统进行可靠度计算,得到基于极点配置的主动控制系统非概率可靠度Pos(sys):
其中,Pos(sys)为基于极点配置的闭环控制系统的非概率可靠度,λcri为闭环控制系统的临界值,即期望设计值。λ为闭环控制系统特征值的下界,为闭环控制系统特征值的上界。此时,有传统方法得到的闭环系统的可靠度以及在表1中给出。
第四步:所述第四步中对期望的闭环控制系统极点进行非概率可靠性优化时采用的非概率可靠性优化模型如下:
findλ
min||u(t)||∞
s.t.Pos(sys)≥Rcri
其中,λ是极点配置的闭环控制系统的特征值,即设计优化变量;||u(t)||∞为控制器输出力H∞范数,Pos(sys)为基于极点配置闭环控制系统的非概率可靠度;Rcri为设计人员要求的可靠度,为一给定的值。
从表1中可以看出,由于不确定性的存在,闭环控制系统极点的实部已经有大于期望值的可能性,这样利用本发明提出的非概率可靠极点配置方法,将闭环控制系统的极点实部最大值配置到期望值处,即将系统的可靠度设计为1。经过迭代以后,图5给出了相应的优化迭代曲线,得到的控制器如下:
表2本发明设计的闭环控制系统极点与可靠度
第五步:利用非概率优化后的极点特征值进行结构振动极点配置控制系统的设计。
(1)利用第一步给出的极点的期望值设计初始的基于极点配置的结构振动主动控制系统,在Matlab/Simulink中建立相应的框图;
(2)利用Matlab中的优化控制箱对基于极点配置的闭环控制系统的特征值进行优化;
(3)得到基于极点配置的闭环控制系统控制器参数后设计最终反馈控制系统。
图6给出了基于非概率可靠优化的结构振动极点配置控制系统的最终响应输出。表3给出了传统极点配置方法和本发明所提出的非概率可靠极点配置方法的响应输出的比较,可以看出,本发明所提出的方法满足可靠度要求。
表3两种方法控制效果比较
Claims (6)
1.一种基于非概率可靠性优化的结构振动极点配置控制方法,其特征在于步骤如下:
第一步:根据结构振动有限元方程和期望得到的闭环控制系统的特征值,建立结构振动极点配置控制系统;
第二步:在第一步的基础上进行结构振动极点配置控制系统的特征值分析,利用区间不确定性分析方法得到结构振动极点配置控制系统特征值的上下界,即其中,λ为闭环控制系统的特征值或者称为极点,分别为特征值的上下界;
第三步:结构振动极点配置控制系统非概率可靠度计算,经过第一步的结构振动极点配置控制系统的建立和第二步的结构振动极点配置控制系统的特征值计算,得到基于极点配置方法的结构振动极点配置控制系统特征值的区间,利用非概率可靠性度量指标的计算方法,对基于极点配置的主动输出反馈控制系统进行可靠度计算,得到基于极点配置的主动控制系统非概率可靠度Pos(sys);
第四步:对给定的期望特征值进行非概率可靠性优化,在满足非概率可靠度Pos(sys)的基础上,得到优化后的闭环控制系统极点特征值;优化目标为使得控制器的输出u(t)的最大值max(u(t))最小;
第五步:利用非概率优化后的极点特征值进行结构振动极点配置控制系统的设计。
2.根据权利要求1所述的基于非概率可靠性优化的结构振动极点配置控制方法,其特征在于:对步骤四进行控制器输出的最大值进行限定,使得基于极点配置的闭环控制系统在满足可靠度约束下,所需控制力最小。
3.根据权利要求1所述的基于非概率可靠性优化的结构振动极点配置控制方法,其特征在于:所述第三步中,非概率可靠性度量指标的计算方法:
设计要求为闭环控制系统响应不超过某一特定值Ycri时,利用如下的计算公式进行非概率可靠度计算:
其中Pos(sys)为基于极点配置的闭环控制系统的非概率可靠度,λcri为闭环控制系统的临界值,即期望设计值,λ为闭环控制系统特征值的下界,为闭环控制系统特征值的上界。
4.根据权利要求1所述的基于非概率可靠性优化的结构振动极点配置控制方法,其特征在于:所述第四步中对期望的闭环控制系统极点进行非概率可靠性优化时采用的非概率可靠性优化模型如下:
findλ
min||u(t)||∞
s.t.Pos(sys)≥Rcri
其中:λ是极点配置的闭环控制系统的特征值,即设计优化变量;||u(t)||∞为控制器输出力H∞范数,Pos(sys)为基于极点配置闭环控制系统的非概率可靠度;Rcri为设计人员要求的可靠度,为给定值。
5.根据权利要求4所述的基于非概率可靠性优化的结构振动极点配置控制方法,其特征在于:为了满足控制器输出最小,给定的设计可靠度Rcri等于1。
6.根据权利要求1所述的基于非概率可靠性优化的结构振动极点配置控制方法,其特征在于:所述第五步,利用经过非概率可靠性优化后的基于极点配置的控制器设计得到结构振动主动控制系统过程如下:
(1)利用第一步给出的极点的期望值设计初始的基于极点配置的结构振动主动控制系统,在Matlab/Simulink中建立相应的框图;
(2)利用Matlab中的优化控制箱对基于极点配置的闭环控制系统的特征值进行优化;
(3)得到基于极点配置的闭环控制系统控制器参数后设计最终反馈控制系统。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610457443.2A CN106094518B (zh) | 2016-06-22 | 2016-06-22 | 一种基于非概率可靠性优化的结构振动极点配置控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610457443.2A CN106094518B (zh) | 2016-06-22 | 2016-06-22 | 一种基于非概率可靠性优化的结构振动极点配置控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106094518A CN106094518A (zh) | 2016-11-09 |
CN106094518B true CN106094518B (zh) | 2018-12-21 |
Family
ID=57237816
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610457443.2A Active CN106094518B (zh) | 2016-06-22 | 2016-06-22 | 一种基于非概率可靠性优化的结构振动极点配置控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106094518B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6731666B2 (ja) * | 2017-12-19 | 2020-07-29 | 国立大学法人神戸大学 | 設計支援装置、設計支援方法、及び設計支援プログラム |
CN108732933B (zh) * | 2018-06-01 | 2021-03-16 | 广东工业大学 | 基于lmi的周期分段振动协同系统的时变连续控制器设计 |
CN109766522B (zh) * | 2019-03-18 | 2022-12-09 | 西安科技大学 | 一种刮板输送机链轮的非概率可靠性灵敏度分析方法 |
CN113110071B (zh) * | 2021-05-31 | 2023-04-14 | 郑州轻工业大学 | 一种基于模态约束的奇异振动结构鲁棒镇定方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH10323071A (ja) * | 1997-05-15 | 1998-12-04 | Toyo Electric Mfg Co Ltd | 2慣性共振系速度制御の2次直列補償器の係数決定方 法 |
CN102142830A (zh) * | 2011-01-31 | 2011-08-03 | 上海大学 | 参考信号自提取的压电智能结构振动主动控制方法 |
CN102323754A (zh) * | 2011-07-29 | 2012-01-18 | 北京理工大学 | 一种基于区域极点配置的鲁棒增益调度控制方法 |
CN103019091A (zh) * | 2012-12-20 | 2013-04-03 | 北京航空航天大学 | 一种基于线性扩张状态观测器的柔性航天器容错姿态控制方法 |
CN103269196A (zh) * | 2013-06-14 | 2013-08-28 | 中国船舶重工集团公司第七一二研究所 | 一种抑制五相电机输出电流不平衡的方法 |
CN103412488A (zh) * | 2013-08-12 | 2013-11-27 | 北京航空航天大学 | 一种基于自适应神经网络的小型无人旋翼机高精度控制方法 |
CN104730922A (zh) * | 2015-01-21 | 2015-06-24 | 浙江工业大学 | 基于扩张状态观测器的伺服系统线性反馈控制和极点配置确定参数方法 |
CN105676903A (zh) * | 2016-04-07 | 2016-06-15 | 北京航空航天大学 | 一种基于非概率可靠性优化的振动最优控制系统设计方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7917256B2 (en) * | 2007-05-02 | 2011-03-29 | Lockheed Martin Corporation | Active vibration damping (AVD) system for precision pointing spacecraft |
-
2016
- 2016-06-22 CN CN201610457443.2A patent/CN106094518B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH10323071A (ja) * | 1997-05-15 | 1998-12-04 | Toyo Electric Mfg Co Ltd | 2慣性共振系速度制御の2次直列補償器の係数決定方 法 |
CN102142830A (zh) * | 2011-01-31 | 2011-08-03 | 上海大学 | 参考信号自提取的压电智能结构振动主动控制方法 |
CN102323754A (zh) * | 2011-07-29 | 2012-01-18 | 北京理工大学 | 一种基于区域极点配置的鲁棒增益调度控制方法 |
CN103019091A (zh) * | 2012-12-20 | 2013-04-03 | 北京航空航天大学 | 一种基于线性扩张状态观测器的柔性航天器容错姿态控制方法 |
CN103269196A (zh) * | 2013-06-14 | 2013-08-28 | 中国船舶重工集团公司第七一二研究所 | 一种抑制五相电机输出电流不平衡的方法 |
CN103412488A (zh) * | 2013-08-12 | 2013-11-27 | 北京航空航天大学 | 一种基于自适应神经网络的小型无人旋翼机高精度控制方法 |
CN104730922A (zh) * | 2015-01-21 | 2015-06-24 | 浙江工业大学 | 基于扩张状态观测器的伺服系统线性反馈控制和极点配置确定参数方法 |
CN105676903A (zh) * | 2016-04-07 | 2016-06-15 | 北京航空航天大学 | 一种基于非概率可靠性优化的振动最优控制系统设计方法 |
Non-Patent Citations (5)
Title |
---|
Active Control of Vibrations in Helicopters via Pole Assignment Techniques;Sergio Bittanti 等;《IEEE TRANSACTIONS ON CONTROL SYSTEMS TECHNOLOGY》;19941231;第2卷(第4期);第343-350页 * |
Pole assignment of friction-induced vibration for stabilisation through state-feedback control;Huajiang 等;《Journal of Sound and Vibration》;20100524;第339卷(第11期);第1985-1991页 * |
基于区间分析的结构非概率可靠性优化设计;祁武超 等;《中国科学》;20131231;第43卷(第1期);第85-91页 * |
模型不确定压电柔性结构的鲁棒振动控制;徐亚兰 等;《机械强度》;20061231;第28卷(第2期);第185-189页 * |
结构非概率鲁棒可靠性准则;邱志平 等;《计算力学学报》;20040228;第21卷(第1期);第1-5页 * |
Also Published As
Publication number | Publication date |
---|---|
CN106094518A (zh) | 2016-11-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106094518B (zh) | 一种基于非概率可靠性优化的结构振动极点配置控制方法 | |
Hiramoto et al. | Optimal sensor/actuator placement for active vibration control using explicit solution of algebraic Riccati equation | |
Zeng et al. | Adaptive feedforward control for gust load alleviation | |
Liao et al. | Subspace identification for fractional order Hammerstein systems based on instrumental variables | |
Khot et al. | Modeling and response analysis of dynamic systems by using ANSYS© and MATLAB© | |
Singla et al. | A new optimal sliding mode controller design using scalar sign function | |
US20060025973A1 (en) | Methods and systems for model reduction and system identification of dynamic systems with multiple inputs | |
Theis et al. | Modal matching for LPV model reduction of aeroservoelastic vehicles | |
Enokida et al. | Stability analysis and comparative experimentation for two substructuring schemes, with a pure time delay in the actuation system | |
Xing et al. | Delay-independent sliding mode control for a class of quasi-linear parabolic distributed parameter systems with time-varying delay | |
Frihauf et al. | Nash equilibrium seeking for games with non-quadratic payoffs | |
CN105892284A (zh) | 一种基于非概率可靠性优化的结构振动pid控制系统设计方法 | |
Abdelaziz | Robust pole assignment for linear time-invariant systems using state-derivative feedback | |
Moulin et al. | Aeroservoelastic structural and control optimization using robust design schemes | |
Knöös et al. | Nonlinear dynamic inversion and block backstepping: a comparison | |
Chang | Dynamic output feedback sliding mode control for uncertain mechanical systems without velocity measurements | |
Li et al. | Optimization of the time-delay margin of L1 adaptive controller via the design of the underlying filter | |
Yao et al. | Revisiting the separation principle for improved transition control | |
Menon et al. | Computer-aided synthesis of nonlinear autopilots for missiles | |
Liu et al. | Decentralized state feedback control of linear time-invariant system with free-interface substructures | |
Zhang et al. | Dynamic output feedback H∞ attitude control for hypersonic gliding vehicles | |
Yuan et al. | A grey-box system identification procedure for scale model helicopters based on frequency-domain estimation methods | |
Su et al. | An unsymmetric Lanczos algorithm for damped structural dynamics systems | |
Zhu et al. | A Modified L 1 Adaptive Output Feedback Control and its Application in Control of Wind Tunnel | |
Jafarnejadsani et al. | Convex Multi-Objective Filter Optimization for Output Feedback L1-Adaptive Controller |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |