CN102176120B - 随机地震激励***磁流变阻尼最优控制的方法 - Google Patents
随机地震激励***磁流变阻尼最优控制的方法 Download PDFInfo
- Publication number
- CN102176120B CN102176120B CN 201110044934 CN201110044934A CN102176120B CN 102176120 B CN102176120 B CN 102176120B CN 201110044934 CN201110044934 CN 201110044934 CN 201110044934 A CN201110044934 A CN 201110044934A CN 102176120 B CN102176120 B CN 102176120B
- Authority
- CN
- China
- Prior art keywords
- control
- active
- optimal control
- random
- design
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Vibration Prevention Devices (AREA)
- Fluid-Damping Devices (AREA)
Abstract
本发明公开了一种随机地震激励***磁流变阻尼最优控制的方法,其根据随机地震激励结构控制***,采用具有两步优化特征的物理随机最优控制策略进行主动最优控制设计;采用***二阶统计量评价准则进行主动最优控制增益设计,得到最优的控制律参数,由此确定用于磁流变阻尼最优控制***设计的目标主动最优控制力;以跟踪实现目标主动最优控制力为控制准则,采用限界Hrovat半主动随机控制策略进行磁流变阻尼控制设计。本发明基于概率密度演化理论框架,突破了以Ito随机微分方程描述动力***的经典随机最优控制的藩篱,为半主动随机最优控制开辟了新的道路。
Description
技术领域
本发明属于地震工程领域,涉及一种随机地震激励***磁流变阻尼最优控制的方法。
背景技术
作用于工程结构的地震荷载在时间、空间和大小上往往具有明显的随机性,结构材料力学特性一般亦具有随机性,因此,工程结构在地震作用下的反应性态表现出不同程度的随机性。传统的结构性态确定性控制与设计方法,已不能满足结构工程发展的要求,结构性态随机控制越来越受到人们的关注。近年来出现了关于磁流变阻尼随机最优控制策略和方法的若干尝试性探讨,如,Dyke等人提出了结构磁流变阻尼控制的限幅线性二次高斯(LinearQuadratic Gaussian,LQG)控制策略;Ying等人基于拟Hamilton***和随机平均法,提出了电/磁流变阻尼***的随机半主动控制策略。
仔细分析可知,上述工作均是经典随机最优控制理论框架下的展开,即随机激励或量测噪声假定为白噪声或过滤白噪声。事实上,(过滤)白噪声过程假定与地震动过程相去甚远,而基于相空间的经典随机最优控制理论恰是依赖于这一假定的。
发明内容
本发明的目的在于提供一种随机地震激励***磁流变阻尼最优控制的方法,根据物理随机最优控制的基本框架,结合限界Hrovat半主动控制算法,以跟踪实现目标主动最优控制力为准则,适用于一般随机地震激励结构控制***磁流变液阻尼器最优参数的设计方法。
为达到以上目的,本发明所采用的解决方案是:
一种随机地震激励***磁流变阻尼最优控制的方法,其包括以下步骤:
1)根据随机地震激励结构控制***,采用具有两步优化特征的物理随机最优控制策略进行主动最优控制设计;
2)采用***二阶统计量评价准则进行主动最优控制增益设计,得到最优的控制律参数,由此确定用于磁流变阻尼最优控制***设计的目标主动最优控制力;
3)以跟踪实现目标主动最优控制力为控制准则,采用限界Hrovat半主动随机控制策略进行磁流变阻尼控制设计。
所述两步优化:一是对于随机参数Θ的每一个实现样本θ,通过性能泛函的最小化,建立控制律参数集合与控制增益集合之间的映射关系;二是根据目标性态,优化控制增益,确定最优的控制律参数。
所述性能泛函为随机过程:
其中,U(t)为依赖于状态量;Z(t)的控制力向量;Q为半正定的状态权矩阵;R为正定的控制力权矩阵;P(t)为Riccati矩阵。
所述性能泛函最小化极值,根据Pontryagin极大值原理构造Euler-Lagrange随机微分方程组,或者根据最优性原理推导Hamilton-Jacobi-Bellman方程求解。
所述***二阶统计量评价准则为:
即在给定约束条件下,
所述限界Hrovat半主动随机控制策略为:根据限界Hrovat控制算法,建立磁流变阻尼装置控制力的分量表达式,根据携带概率信息的样本关于位移、加速度、黏滞阻尼系数、库伦力与主动控制力的关系,结合最小二乘拟合曲线,获得磁流变液阻尼器的最优设计参数。
由于采用了上述方案,本发明具有以下特点:本发明根据物理随机最优控制的基本框架,结合限界Hrovat半主动控制算法,以跟踪实现目标主动最优控制力为准则,提出了适用于一般随机地震激励结构控制***磁流变液阻尼器最优参数的设计方法,突破了以Ito随机微分方程描述动力***的经典随机最优控制的藩篱,为半主动随机最优控制开辟了新的道路。
附图说明
图1为物理随机最优控制、线性二次高斯控制、确定性控制响应的位移与权矩阵系数比的关系图。
图2为物理随机最优控制、线性二次高斯控制、确定性控制响应的控制力与权矩阵系数比的关系图。
图3为半主动最优控制力在均方根意义上追踪主动最优控制力的示意图。
图4为磁流变液阻尼器的阻尼力曲线及模型比较的位移与阻尼力关系图。
图5为磁流变液阻尼器的阻尼力曲线及模型比较的速度与阻尼力关系图。
图6为随机地震激励***磁流变阻尼最优控制设计方法框图。
具体实施方式
以下结合附图所示实施例对本发明作进一步的说明。
1.物理随机最优控制的基本框架
考察随机动力***:
对于关心的物理量
上述随机动力***的最优控制涉及特定二次性能泛函的最大化或最小化,在随机参数向量Θ的背景下,二次性能泛函定义为:
其中,U(t)为依赖于状态量Z(t)的控制力向量,显然其与随机参数向量Θ构成的***也满足广义密度演化方程(3);Q为半正定的状态权矩阵;R为正定的控制力权矩阵;P(t)为Riccati矩阵。
在经典LQG控制中,性能泛函定义为对式(4)取期望的形式,因此性能泛函为确定性的时变函数,而本发明中定义的性能泛函是一随机过程,由此发展的物理随机最优控制方法涉及两步优化:一是对于随机参数Θ的每一个实现样本θ,通过性能泛函的最小化,建立控制律参数集合与控制增益集合之间的映射关系;二是根据目标性态,优化控制增益,确定最优的控制律参数。
2.***二阶统计量评价准则
依赖于性能泛函形式的***控制,无论是经典的确定性泛函、还是随机泛函式(4),控制效果都取决于设定的控制律参数。事实上,结构***性态精细化控制的核心是控制律及其参数的设计与优化,而其准则恰恰是依赖于结构反应性态的。本发明提出了基于***二阶统计量评价的控制律设计方法。即在给定约束条件下,***二阶统计量评价准则为:
显然,***二阶统计量评价准则可以综合评价***状态和控制力,寻求各物理量之中出现最小均值或者最小标准差时的控制律参数,以达到柔性意义上的最佳均衡。采用这一准则,物理随机最优控制可以客观、合理地控制***的性态,而采用名义白噪声输入的经典线性二次高斯控制,不能合理设计土木工程结构控制***。以随机地震动作用下某剪切型框架结构主动锚索控制为例进行比较分析,如图1所示:物理随机最优控制具有与确定性控制相似的位移-控制力关系,表明了物理随机最优控制的样本轨道描述特性和概率密度控制特征;而线性二次高斯控制在较小权矩阵系数比时会低估所需要的控制力,在较大权矩阵系数比时会高估所需要的控制力。
3.限界Hrovat半主动随机控制策略
根据限界Hrovat控制算法,建立磁流变阻尼装置控制力的分量表达式:
式(6)中,Ua(Θ,t)是主动最优控制力;是磁流变液阻尼器在t时刻可能实现的最大阻尼力;Cd为黏滞阻尼系数,Udc,max,Udc,min分别是磁流变液阻尼器的最大、最小库伦阻尼力。其中,Udc,max,Udc,min,Cd均为磁流变液阻尼器参数。不难看出,随机动力***的最优控制力分量Us(Θ,t),Ua(Θ,t),速度分量均满足广义密度演化方程式(3),因此它们的概率密度函数和样本轨迹表现均可以方便地得到。
对于前述的半主动控制***,为了获得与主动随机最优控制相近的控制效果,设计磁流变液阻尼器,使包括被动黏滞阻尼力在内的最大阻尼力等于最大(极值)主动最优控制力。假定磁流变阻尼控制与主动最优控制的效果一致(设计最大控制力时层间速度相同),即
假定最小库伦阻尼力Udc,min=0,阻尼力的可调倍数s。则由式(7)有
于是,被动黏滞阻尼系数
磁流变液阻尼器的最大库伦力由式(7)和式(8)得到
上述表明,对于考虑随机地震动输入的结构控制***,***状态量及其相关的最优控制力为随机过程,在样本空间***控制律参数亦表现出不确定性,如Udc,max,Cd的客观表现形式是依赖于Θ的。然而,***控制律参数是预先设定的确定模式。即:无论针对什么样的随机输入样本,***控制的基本设计参数是不变的。基于这一理念,无论是开环控制(如被动控制***)、还是闭环控制(如输入反馈或状态反馈的主动、半主动控制***),控制律参数设计是满足***目标性态的一组确定性设计值。
根据携带概率信息的样本关于位移、加速度、黏滞阻尼系数、库伦力与主动控制力的关系,并结合最小二乘拟合曲线,可以得到磁流变液阻尼器的最优设计参数。
4.磁流变液阻尼器的工作性态
以前述随机地震动作用的结构磁流变阻尼控制为例,进行半主动最优控制力与主动最优控制力的比较,以及磁流变液阻尼器的阻尼力曲线与试验模型的比较,如图2和图3所示。
图2表明,在二阶统计特征意义上,半主动控制器可以达到与主动控制器几乎相同的控制效果,其实施控制力的细微差异在于:当主动控制力与层间速度同向,或者主动控制力与层间速度反向、主动控制力大于磁流变液阻尼器可能实现的最大阻尼力时,半主动控制实施的阻尼力仅与最大主动控制力相关,而不是实时跟踪主动最优控制力。
图3给出了代表性地震动作用下磁流变液阻尼器的阻尼力曲线(位移与阻尼力关系、速度与阻尼力关系)。可以看出:速度与阻尼力几乎只出现在第1、III象限,即半主动控制力的方向基本上在所有的时刻都与速度方向相反,这证明了采用磁流变阻尼力装置可以实时调整阻尼力、跟踪实现主动最优控制力的合理性。同时,图中示意了Dyke等人关于磁流变液阻尼器性能的仿真预测与试验比例轮廓曲线(他们以LORD公司生产的200kN磁流变液阻尼器为考察对象,采用扩展Bouc-Wen模型建模,得到了窄带高斯激励下的阻尼力曲线)。对比表明,限界Hrovat控制策略能充分发挥磁流变液阻尼器的动态阻尼性能,恢复力曲线表现出了似Bouc-Wen模型的强度退化、刚度退化以及捏拢效应。
随机地震激励***磁流变阻尼最优控制设计方法流程如图4所示:
●根据随机地震激励结构控制***,采用具有两步优化特征的物理随机最优控制策略进行主动最优控制设计。关于性能泛函(4)的条件极值问题,既可以根据Pontryagin极大值原理构造Euler-Lagrange随机微分方程组,也可以根据最优性原理推导Hamilton-Jacobi-Bellman方程(HJB方程)。实施例中,采用了离散动态规划算法求解HJB方程,并采用与离散动态规划算法相同的一阶向前差分格式进行广义密度演化方程中速度量求解的确定性动力反应分析。
●采用***二阶统计量评价准则进行主动最优控制增益设计,得到最优的控制律参数,由此确定目标主动最优控制力。
●采用限界Hrovat半主动随机控制策略进行磁流变阻尼控制设计。根据携带概率信息的样本关于位移、加速度、黏滞阻尼系数、库伦力与主动控制力的关系,结合最小二乘拟合曲线,获得磁流变液阻尼器的最优设计参数。
上述的对实施例的描述是为便于该技术领域的普通技术人员能理解和应用本发明。熟悉本领域技术的人员显然可以容易地对这些实施例做出各种修改,并把在此说明的一般原理应用到其他实施例中而不必经过创造性的劳动。因此,本发明不限于这里的实施例,本领域技术人员根据本发明的揭示,对于本发明做出的改进和修改都应该在本发明的保护范围之内。
Claims (1)
1.一种随机地震激励***磁流变阻尼最优控制的方法,其特征在于:其包括以下步骤:
1)根据随机地震激励结构控制***,采用具有两步优化特征的物理随机最优控制策略进行主动最优控制设计;
2)采用***二阶统计量评价准则进行主动最优控制增益设计,得到最优的控制律参数,由此确定用于磁流变阻尼最优控制***设计的目标主动最优控制力;
3)以跟踪实现目标主动最优控制力为控制准则,采用限界Hrovat半主动随机控制策略进行磁流变阻尼控制设计,
所述两步优化:一是对于随机参数Θ的每一个实现样本θ,通过性能泛函的最小化,建立控制律参数集合与控制增益集合之间的映射关系;二是根据目标性态,优化控制增益,确定最优的控制律参数;
所述限界Hrovat半主动随机控制策略为:根据限界Hrovat控制算法,建立磁流变阻尼装置控制力的分量表达式,根据携带概率信息的样本关于位移、加速度、黏滞阻尼系数、库伦力与主动控制力的关系,结合最小二乘拟合曲线,获得磁流变液阻尼器的最优设计参数;
所述性能泛函为随机过程:
所述性能泛函最小化极值,根据Pontryagin极大值原理构造Euler-Lagrange随机微分方程组,或者根据最优性原理推导Hamilton-Jacobi-Bellman方程求解;
所述***二阶统计量评价准则为:
即在给定约束条件下,
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110044934 CN102176120B (zh) | 2011-02-24 | 2011-02-24 | 随机地震激励***磁流变阻尼最优控制的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110044934 CN102176120B (zh) | 2011-02-24 | 2011-02-24 | 随机地震激励***磁流变阻尼最优控制的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102176120A CN102176120A (zh) | 2011-09-07 |
CN102176120B true CN102176120B (zh) | 2013-08-14 |
Family
ID=44519311
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201110044934 Expired - Fee Related CN102176120B (zh) | 2011-02-24 | 2011-02-24 | 随机地震激励***磁流变阻尼最优控制的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102176120B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102866631B (zh) * | 2011-12-31 | 2015-09-09 | 上海绿地建设(集团)有限公司 | 基于在线支持向量机的建筑结构地震反应智能控制设计方法 |
JP6523854B2 (ja) * | 2015-07-29 | 2019-06-05 | 株式会社東芝 | 最適制御装置、最適制御方法、コンピュータプログラム及び最適制御システム |
CN106776471A (zh) * | 2016-11-25 | 2017-05-31 | 上海电机学院 | 一种基于随机平均法的概率密度追踪控制方法 |
CN109507891B (zh) * | 2019-01-21 | 2021-07-27 | 闽江学院 | 一种半主动模糊控制方法 |
CN112904718B (zh) * | 2021-01-14 | 2022-10-21 | 吉林大学 | 一种基于Hammerstein模型的磁流变阻尼器控制***及方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101016758A (zh) * | 2007-01-22 | 2007-08-15 | 哈尔滨工业大学 | 结构振动控制的新型混合质量驱动变阻尼控制装置 |
CN101592919A (zh) * | 2009-07-08 | 2009-12-02 | 江南大学 | 一类不确定中立***时滞依赖鲁棒镇定策略 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4197371B2 (ja) * | 1999-03-03 | 2008-12-17 | 大成建設株式会社 | 制振方法 |
US7216018B2 (en) * | 2004-05-14 | 2007-05-08 | Massachusetts Institute Of Technology | Active control vibration isolation using dynamic manifold |
-
2011
- 2011-02-24 CN CN 201110044934 patent/CN102176120B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101016758A (zh) * | 2007-01-22 | 2007-08-15 | 哈尔滨工业大学 | 结构振动控制的新型混合质量驱动变阻尼控制装置 |
CN101592919A (zh) * | 2009-07-08 | 2009-12-02 | 江南大学 | 一类不确定中立***时滞依赖鲁棒镇定策略 |
Non-Patent Citations (1)
Title |
---|
JP特开2000-250638A 2000.09.14 |
Also Published As
Publication number | Publication date |
---|---|
CN102176120A (zh) | 2011-09-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102176120B (zh) | 随机地震激励***磁流变阻尼最优控制的方法 | |
Zhou et al. | Adaptive fuzzy control for nonlinear building-magnetorheological damper system | |
CN103809458B (zh) | 一种基于改进蛙跳算法的磁流变阻尼控制方法 | |
Tang et al. | Regenerative semi-active control of tall building vibration with series TMDs | |
CN105487384A (zh) | 一种基于事件触发机制的整车悬架控制***及其设计方法 | |
CN106843306B (zh) | 浮置板轨道及其振动控制方法和装置 | |
Mitchell et al. | Active control of highway bridges subject to a variety of earthquake loads | |
CN102866631B (zh) | 基于在线支持向量机的建筑结构地震反应智能控制设计方法 | |
CN105155717A (zh) | 一种基于殖民竞争算法的磁流变阻尼模糊控制器设计方法 | |
CN105160100A (zh) | 安装弹簧-质量***的tmd最优设计方法 | |
Esmaeili et al. | Load-dependent LPV/H2 output-feedback control of semi-active suspension systems equipped with MR damper | |
Zhan et al. | Joint optimization approach to building vibration control via multiple active tuned mass dampers | |
CN103197541A (zh) | 一种基于混沌***的模糊控制方法 | |
Xu et al. | Neural network predictive control of vehicle suspension | |
Koumboulis et al. | A three term controller for ride comfort improvement | |
CN103942394A (zh) | 一种基于磁流变阻尼器限幅最优的方法 | |
Halperin et al. | Using constrained bilinear quadratic regulator for the optimal semi-active control problem | |
Metwally et al. | Reducing vibration effects on buildings due to earthquake using magneto-rheological dampers | |
Oliveira et al. | Predictive control for earthquake response mitigation of buildings using semiactive fluid dampers | |
Bitaraf et al. | Structural performance improvement using MR dampers with adaptive control method | |
Chae et al. | Evaluation of structural control strategies for improving seismic performance of buildings with MR Dampers using real-time large-scale hybrid simulation | |
Villamizar et al. | Experimental verification of a Backstepping controller for magnetorheological MR dampers in structural control | |
Pnevmatikos et al. | Actively and semi-actively controlled structures under seismic actions: modeling and analysis | |
Guo et al. | Semi-active seismic isolation control based on MR damper and neural networks | |
Kumar et al. | Development of An Adaptive Linear Quadratic Gaussian (Lqg) Controller For Structural Control Using Particle Swarm Optimization |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20130814 Termination date: 20160224 |