一种基于首次穿越理论的振动主动控制系统非概率动力可靠性分析方法
【专利摘要】本发明公开了一种基于首次穿越理论的振动主动控制系统非概率动力可靠性分析方法。该方法综合考虑材料的区间不确定性对主动控制系统适定性的影响。基于一阶状态空间方程格式和状态反馈控制原理，推导了主动控制系统的动力响应表达式；基于区间数学思想和时变不确定性传播分析方法，构建了受控结构动力响应的非概率区间过程模型；进而利用首次穿越理论和结构的控制率失效准则，定义了振动主动控制系统非概率动力可靠性指标，并探索了高效稳健的求解方法。本发明在进行可靠度计算时体现了材料分散性对受控结构动力响应的综合影响，可为控制器的精细化设计提供必要且可行的理论参考。
【专利说明】
-种基于首次穿越理论的振动主动控制系统非概率动力可靠 性分析方法
技术领域
[0001] 本发明设及振动主动控制系统的安全性评估技术领域，特别设及考虑材料分散性 作用下结构安全态势的合理辨识和有效度量，为进一步规划大型结构可靠控制理念的工程 化应用W及面向主动控制系统优化设计方案的制定提供可借鉴的理论依据。
【背景技术】
[0002] 在工程结构问题中，振动现象无处不在，振动效应对结构造成的影响是难W避免 的。多数情况下，振动效应对结构造成的影响是有害的。如何减小振动效应带来的危害，必 须要考虑对振动如何加W控制的问题。结构振动控制技术是振动力学里的一个很重要的方 面。现有的振动控制包含W下两方面:主动控制和被动控制。在众多的工程应用中，被动控 制具有结构简单，不依赖外界能源而且减震性能良好等优点。但随着科学技术的发展，实际 工程结构变得越来越复杂，由此导致振动控制问题也越来越复杂，传统的被动控制在处理 复杂的工程问题时凸显出了如控制效果差、控制效率低等问题。面对实际工程领域的广泛 需求，振动主动控制技术得W发展，并受到了广泛的关注。
[0003] 然而，振动主动控制技术在发展过程中暴露了一个明显问题，就是对于结构自身 存在的不确定性因素，特别是材料分散性极为敏感。当系统中出现的不确定因素对系统影 响较大时，毫无疑问，会降低控制系统的可控性和稳定性。关于振动主动控制系统的可靠性 问题也引起了学术界和工程界的高度重视。
[0004] 当前，国内外学者与工程技术人员对含主动控制结构的不确定性分析方法研究主 要集中在两个方面：（1)基于鲁棒控制理论的结构不确定性传播影响预测技术；（2)基于准 静态假设的概率可靠度计算方法研究。上述工作一定程度上确保了振动主动控制系统的适 定性，但是忽略了鲁棒约束的过冗余性W及时间累积效应下结构失效事件的相关性，大大 限制了主动控制理论的工程实用化进程。
[0005] 由于实际工程中贫信息、少数据的情况无法避免，建立基于非概率集合理论的不 确定性分析技术、动力可靠度建模与求解技术对主动控制系统具有显著的现实意义。本发 明将重点探究材料分散性作用对受控结构安全性的综合影响，为其控制器的优化设计提供 技术保障。
【发明内容】
[0006] 本发明要解决的技术问题是：克服现有技术的不足，提供一种针对振动主动控制 系统的安全性评价方法，充分考虑实际工程问题中普遍存在的材料分散性效应，W提出的 非概率动力可靠性度量作为评判受控结构安全与否的量化标准，所得到的校核结果更加符 合真实情况，工程适用性更强。
[0007] 本发明采用的技术方案为:一种基于首次穿越理论的振动主动控制系统非概率动 力可靠性分析方法，该方法实现步骤如下：
[0008] 第一步:建立n自由度系统在外部载荷和控制载荷共同作用下的振动微分方程：
[0009]
[0010] 其中，M，P，K分别为总体质量矩阵、总体阻尼矩阵和总体刚度矩阵，撕)，拂)，x(t) 分别表示受控结构的加速度、速度和位移向量，x(t〇)和i?蛇)分别表示初始时刻to对应的位 移和速度向量，^〇和而分别是初始位移和初始速度的给定输入条件，f(t)为外部载荷向量， fa(t)代表控制载荷向量，并可表示为：
[0011] fa(t)=Ball(t)
[0012] 其中，Ba为nXr阶权重矩阵，表示驱动器的作用位置，u(t)为r Xl阶驱动载荷向 量。
[0013] 第二步:基于经典主动控制理论，建立第一步中振动微分方程的一阶状态空间方 程：
[0014]
[0015] y(t)=Cz(t)
[0016] 其中
I示化Xl阶的状态向量，u(t)和y(t)分别对应为驱动 载荷向量和输出向量，&灼表示状态向量对时间的导数，命名为状态响应向量，C为输出矩 阵，A，B，E为状态空间离散的特征矩阵，具体表示为：
[0017]
[001引其中，功n X n阶单位矩阵，0为零矩阵，M-嗦示总体质量矩阵M的逆矩阵。
[0019] 假定结构状态是可测的，应用状态反馈控制原理，可知：
[0020] u(t) =-Gz(t)
[0021 ]其中，G为控制器。代入到前式可得到状态响应向量的具体表达如下：
[0022]
[0023] 第S步:综合考虑存在于第一步和第二步中材料参数的分散性，引入区间向量b = [bl，b2, ...，bm]，得到有限样本条件下振动控制系统中不确定性参数的数学表达：
[0024]
[002引其中，b嗦示向量b的区间集，巧日i分别为向量b的下界和上界，b。和A b分别为向 量b的均值和半径，C定义为所有元素包含在[-1，1 ]内的m阶向量集合。运里，向量b均值和半 径的具体表达式为：
[0026]
[0027] 其中，琴和Ab康示向量b中第j个元素b说均值和半径，j为计数指标。
[0028] 第四步:将第=步中表征的材料分散性参数代入到第二步的一阶状态空间方程和 状态响应向量的表达式中，可得：
[0029]
[0030] y(t) =Cz(b,t)
[0031] 和
[0032]
[0033] 上式中，由于材料分散性所带来的不确定性效应会导致矩阵M，P，K的取值发生变 化，状态空间离散的特征矩阵实际上也是区间向量b的函数，故表示为4化），8化），6化）；此 夕h由于结构响应计算采用迭代方法求解，材料分散性带来的时间累积效应不可忽略，于是 z(t)^z(b，t)和雖）->游，〇<;综上，从不确定性传播分析的角度出发，受控结构的动力响 应可由一个域集r给予包络表达，即：
[00341
[0035] 其中，x(b，t)表示考虑材料分散性后主动控制系统的位移向量，〇[?]反映了位 移向量X(b，t)和状态向量z(b，t)之间的映射关系。借助区间传播分析方法，位移向量x(b， t)的变化前亲元九W下反间化古，
[0036]
[0037] 其中，xi(b，t)为位移向量x(b，t)的区间集，互(b，t)和-V-0M)分别为x(b，t)的下界 和上界。
[0038] 第五步:根据时变不确定性分析方法，构建出第四步中位移向量x(b，t)中各个元 素又1化，〇4 = 1，2，...，11的非概率区间过程模型，并分别定义该区间过程模型的特征量函 数，具体包括:均值过程函数怒(我巧、半径过程函数.<(&，〇 W及任意不同时刻ti和t2的相关 系数函数其中，均值过程函数為(々，〇和半径过程函数為(6,〇的具体表达式如下：
[0039]
[0040] 任意不同时刻*1和*2的相关系数函数八,也^^)可由下式求解：
[0041]
[0042] 其中，表示协方差函数。
[0043] 第六步:结合控制系统响应的许用值向量)T，构建与第五步获得的位移向量x(b， t)区间历程间的应力-强度区间过程干设模型，并建立时变极限状态函数如下：
[0044]
[0045] 基于首次穿越理论，通过时同离散化时段，定义如下穿越事件Ei化)发生的可能度 指标：
[0046] F*os{Ei(^k)}=F*os{gi(b,kAt)〉〇ngi(b,(k+l)At)<0}
[0047] 其中，Posl ? }表示事件发生的可能度，gi(b，kAt)〉0表示受控结构在kAt时刻安 全(位移响应小于许用值），gi(b，化+1) AtKO表示受控结构在(k+1) At时刻失效(位移响 应大于许用值），符号"n"表示事件的交运算，k为计数指标，At表示时间增量。运里，事件 Ei化)表示为受控结构比A t，化+1) A t]在时间段内发生了一次穿越，A t通常为一微小量。
[0048] 第屯步：遍历所有时间段内的可能度指标Pos {Ei化）}，建立主动控制系统非概率 动力可靠度指标：
[0049]
[0050] 其中，Rs(T)为动力可靠度，T为完整响应历程，Pos(O)表示受控结构初始失效的可 能度，求解上式即可实现主动控制系统动力安全态势的有效评估。
[0051] 本发明与现有技术相比的优点在于:本发明提供了考虑材料分散性效应下振动主 动控制系统安全性评价的新思路，弥补和完善了传统基于鲁棒控制理论和静态概率可靠性 理论方法的局限性。所构建的非概率动力可靠性度量模型，一方面弥补了鲁棒控制所带来 的设计冗余度，另一方面较静态可靠性方法，基于首次穿越理论的动力可靠度计算方法，更 加合理地计及了动力响应的时间相关性，确保控制器全时间域内的有效性，为控制器精细 化设计理论的完善作出了积极贡献。
【附图说明】
[0052] 图1是本发明针对振动主动控制系统的非概率动力可靠性评估流程图；
[0053] 图2是本发明基于状态反馈原理的状态空间离散策略示意图；
[0054] 图3是本发明定义的位移响应Xi (b，ti)和Xi (b，t2)的几何可行域示意图；
[005引图4是本发明定义的位移响应xi(b，ti)和Xi(b，t2)的相关性示意图；
[0056] 图5是本发明提出的微小时间段内穿越失效可能度计算方法示意图；
[0057] 图6是本发明实施例中含主动控制十杆巧架结构的模型示意图；
[0058] 图7是本发明实施例中含主动控制十杆巧架结构的位移响应历程示意图；
[0059] 图8是本发明实施例中含主动控制十杆巧架结构的可靠度计算结果示意图，其中， 图8(a)为受控结构节点1处水平方向的动力可靠度计算结果，图8(b)为受控结构节点1处竖 直方向的动力可靠度计算结果。
【具体实施方式】
[0060] 下面结合附图W及【具体实施方式】进一步说明本发明。
[0061] 如图1所示，本发明提出了一种基于首次穿越理论的振动主动控制系统非概率动 力可靠性分析方法，包括W下步骤：
[0062] (1)建立n自由度系统在外部载荷和控制载荷共同作用下的振动微分方程：
[0063]
[0064] 其中，M，P，K分别为总体质量矩阵、总体阻尼矩阵和总体刚度矩阵，W〇，i：(/)，x(t) 分别表示受控结构的加速度、速度和位移向量，X (to)和却4)分别表示初始时刻to对应的位 移和速度向量，XO和為分别是初始位移和初始速度的给定输入条件，f(t)为外部载荷向量， fa(t)代表控制载荷向量，并可表示为：
[00化]
[0066] 其中，Ba为nXr阶权重矩阵，表示驱动器的作用位置，u(t)为r X I阶驱动载荷向 量。
[0067] (2)其干经化丰动梓制巧换.建立第一步中振动微分方程的一阶状态空间方程： [006引
[0069] y(t) = Cz(t)
[0070] 其中
表示化Xl阶的状态向量，u(t)和y(t)分别对应为驱动 载荷向量和输出向量，的表不状态向量对时间的导数，命名为状态响应向量，C为输出矩 阵，A，B，E为状态空间离散的特征矩阵，具体表示为：
[0071]
[00巧其中，I为nXn阶单位矩阵，0为零矩阵，M-嗦示总体质量矩阵M的逆矩阵。
[0073] 假定结构状态是可测的，应用状态反馈控制原理，可知：
[0074] u(t) =-Gz(t)
[0075] 其中，G为控制器。代入到前式可得到状态响应向量的具体表达如下：
[0076]
[0077] 于是，上述基于状态反馈的振动主动控制系统状态空间离散策略如图2所示。
[0078] (3)综合考虑存在于第一步和第二步中材料参数的分散性，引入区间向量b=[bi， b2, ...，bm]，得到有限样本条件下振动控制系统中不确定性参数的数学表达：
[0079]
[0080] 具甲，巧不同重b的K间集，叫U6分別刃同量b的下界和上界，b嘴Ab分别为向 量b的均值和半径，C定义为所有元素包含在[-1，1 ]内的m阶向量集合。运里，向量b均值和半 径的具体表达式为：
[0081]
[0082] 其中，6巧日Abj表示向量b中第j个元素 b说均值和半径，j为计数指标。进而，区间 向量b满足如下关系式：
[0083] b = b"+5, |5|^ Ab
[0084] 其中，向量S在几何上被界定于A b所定义的超长方体域内。
[0085] (4)将第=步中表征的材料分散性参数代入到第二步的一阶状态空间方程和状态 响应向量的表试式中，可得：
[0086]
[0087]
[008引
[0089]
[0090] 上式中，由于材料分散性所带来的不确定性效应会导致矩阵M，P，K的取值发生变 化，状态空间离散的特征矩阵实际上也是区间向量b的函数，故表示为4化），8化），6化）；此 夕h由于结构响应计算采用迭代方法求解，材料分散性带来的时间累积效应不可忽略，于是 z(t)^z(b，t)和ミ的^ミ成巧。综上，从不确定性传播分析的角度出发，受控结构的动力响 应可由一个域集r给予包络表达，即：
[00911
[0092] 其中，x(b，t)表示考虑材料分散性后主动控制系统的位移向量，(6[-]反映了位 移向量x(b，t)和状态向量z(b，t)之间的映射关系。借助区间传播分析方法，位移向量x(b, t)的变化可表示为W下区间形式：
[0093]
[0094] 其中，xi(b，t)为位移向量x(b，t)的区间集，2[(b，t)和S(ZM)分别为x(b，t)的下界 和上界。基于一阶泰勒级数展开法，受控结构不确定动力响应的上下界可W显式求解。当 然，首先应计算特征矩阵A(b)和B(b)、状态向量Z (b，t)的近似表达，即：
[0095]
[0096]
[0097]
[009引 C-A(y) rA(d') cz{h\t]
[0099] 其中，一:^，一:^，'分别表示特征矩阵A(b)和B(b)及状态向量z(b，t) ch. V I) oh, J J J 在均值bK处相对于变量扁导数。将上面表达式代入到状态响应向量的公式中，可得：
[0100]
[0101]
[0102]
[0103]
[0104]
[0105] 借助于精确时间积分方法，上式可被快速求解。结合自然区间扩张原理，可知：
[0106] zi(b,t) = z(bc，t)+A zi(t)
[0107] 其中，因此，位移向量x(b，t)的下界和上界可W最终表示 为：
[010 引
[0109] (5)根据时变不确定性分析方法，构建出第四步中位移向量x(b，t)中各个元素Xi 化，〇4 = 1，2，...，11的非概率区间过程模型，并分别定义该区间过程模型的特征量函数， 具体包括:均值过程函数為(々坤、半径过程函数抑伊句W及任意不同时刻ti和t2的相关系 数函数A;(如2)。其中，均值过程函数.和半径过程函数看(6,〇的具体表达式如下：
[0110]
[om]任意不同时刻tl和t2的相关系数函数A(/|，M可由下式求解：
[0112]
[0113] 其中，表示协方差函数。为了计算出Cw, 的表达式，我们需要借 助正则化手段，于是有：
[0114]
[0115]
[0116]
[0117]其中，Cl和C2分别为标准区间变量。运样，一个C坐标系下的标准方形域可W获得W 等效量化位移响应Xi(b，ti)和Xi(b，t2)的可行范围（如图3所示）。显然，受到相关性条件的 约束，存在多个不同形状的偏转矩形域包络于标准方形域内，并与相关性指标一一对应(如 图4所示）。因此，Cov, (/|，/:)可定义为：
[011 引
[0119] 其中，d表示如图4所示矩形域边长的一半。受到随机过程理论的启发，可进一步计 算：
[0120]
[01別]于是，任意不同时刻ti和t2的相关系数函数A,也O最终表示为：
[0122]
[012；3] 运里，如是一个无量纲量，其大小代表了 xi(b，ti)和xi(b，t2)的线性相关度。 [0124]综上，本发明实现了对受控结构不确定性位移响应的定量表征，为后续非概率动 力可靠性建模及求解提供了必要的理论依据。
[0125] (6)结合控制系统响应的许用值向量)T，构建与第五步获得的位移向量x(b，t)区 间历程间的应力-强度区间过程干设模型，并建立时变极限状态函数如下：
[0126]
[0127]基于首次穿越理论，通过时间离散化时段，定义如下穿越事件Ei化)发生的可能度 指标：
[01 巧]F*os{EiA)}=F*os{gi(b,kAt)〉〇ngi(b,A+l)At)<0}
[0129] 其中，Posl ? }表示事件发生的可能度，gi(b，kAt)〉0表示受控结构在kAt时刻安 全(位移响应小于许用值），gi(b，化+1) AtKO表示受控结构在(k+1) At时刻失效(位移响 应大于许用值），符号"n"表示事件的交运算，k为计数指标，At表示时间增量。运里，事件 Ei化)表示为受控结构比A t，化+1) A t]在时间段内发生了一次穿越，A t通常为一微小量。 运里，引入面积比思想（如图5所示），化s{Ei化）}可定义为穿越几何条件与极限状态可行域 的干设面积与过程中总可行域(偏转矩形)面积之比，即：
[0130]
[0131] 4^ 影的 计算通常是一个分段函数，需要结合几何边界与可行域的相交条件分类讨论。
[0132] (7)遍历所有时间段内的可能度指标化s{Ei化）}，建立主动控制系统非概率动力 可靠度指标：
[0133]
[0134] 其中，Rs(T)为动力可靠度，T为完整响应历程，Pos(O)表示受控结构初始失效的可 能度，求解上式即可实现主动控制系统动力安全态势的有效评估。结合上一步的计算结果， 最终Rs (T)可表示为：
[0135]
[0136] 实施例；
[0137] 为了更充分地了解该发明的特点及其对工程实际的适用性，本发明针对如图6所 示的十杆巧架结构进行主动控制后非概率动力可靠性分析。水平杆和竖直杆的长度为 50Omm，杆的横截面积为40Omm2，结构的整体阻尼矩阵为P = aM+邸，其中，a = 1.85 59，0 = 6.2516X10-6。本例中，存在两组外部载荷fl(t)和f2(t)，并且有：
[013 引 fi(t) = 1000e(-3t)sin(1003it)kN 和 f2(t) = 800e(-2t)cos(7 加 t)kN，tG[0,2]年
[0139]此外，两组驱动器定位于节点2和4处，并分别施加主动控制fal(t)和fa2(t)。于是， 定义节点1在水平和竖直方向的位移被控制在许用值x" = 7.5mm内。考虑材料的分散性效 应，假定弹性模量E和密度P分别为区间变量，具体参数信息见表1，其中，变化系数K反映了 不同的区间跨度。
[0140] 表1
[0141]
[0142] 结合本发明提出的受控结构响应计算方法和区间过程模型，得到节点I的位移响 应历程(如图7所示）。基于首次穿越理论，该主动控制系统的非概率动力可靠性分析结果如 图8和表2所示。从结果中可W看出，相较于传统的准静态(时不变)处理方法，本发明提出的 非概率动力可靠性分析方法，由于计及了时间相关性效应，结果上更为保守且贴合真实情 况;此外，随着变化系数K的增大，主动控制系统的可靠度降低明显，所带来的安全性隐患加 剧。
[0143] 综上所述，本发明提出了一种基于首次穿越理论的振动主动控制系统非概率动力 可靠性分析方法。该方法通过区间变量表征材料的分散性，进而讨论区间不确定性对振动 主动控制系统安全性的综合影响。首先，基于状态反馈控制原理和区间过程理论，推导了主 动控制系统动力响应的时变不确定性模型;进而利用首次穿越理论，定义了振动主动控制 系统非概率动力可靠性指标;基于时间离散策略和面积比思想，通过遍历所有时间段内的 穿越可能度，叠加运算后最终给出了本发明所提出的非概率动力可靠度显式求解策略。
[0144] 表2
[0145]
[0146] W上仅是本发明的具体步骤，对本发明的保护范围不构成任何限制；其可扩展应 用于振动主动控制系统的可靠设计领域，凡采用等同变换或者等效替换而形成的技术方 案，均落在本发明权利保护范围之内。
[0147]本发明未详细阐述部分属于本领域技术人员的公知技术。
【主权项】
1. 一种基于首次穿越理论的振动主动控制系统非概率动力可靠性分析方法，其特征在 于实现步骤如下： 第一步:建立η自由度系统在外部载荷和控制载荷共同作用下的振动微分方程：其中，Μ，Ρ，Κ分别为总体质量矩阵、总体阻尼矩阵和总体刚度矩阵，对〇, ±(/)，x(t)分别 表示受控结构的加速度、速度和位移向量，X (to)和)分别表示初始时刻to对应的位移和 速度向量，XdPie分别是初始位移和初始速度的给定输入条件，f (t)为外部载荷向量，fa (t) 代表控制载荷向量； 第二步:基于经典主动控制理论，建立第一步中振动微分方程的一阶状态空间方程：其中，ζ(〇 = V (?)表示2n X 1阶的状态向量，u(t)和y (t)分别对应为驱动载荷 向量和输出向量，汾)表示状态向量对时间的导数，命名为状态响应向量，A，B，E为状态空间 离散的特征矩阵，C为输出矩阵，假定结构状态是可测的，应用状态反馈控制原理，可知： u (t) = -Gz (t) 其中，G为控制器，代入到前式可得到状态响应向量的具体表达如下：第三步:综合考虑存在于第一步和第二步中材料参数的分散性，引入区间向量b=[h， b2,...，bm]，得到有限样本条件下振动控制系统中不确定性参数的数学表达：其中，b1表不向量b的区间集，h和石分别为向量b的下界和上界，b。和Δ b分别为向量b的 均值和半径，ξ定义为所有元素包含在[_1，1]内的m阶向量集合； 第四步:将第三步中表征的材料分散性参数代入到第二步的一阶状态空间方程和状态 响应向量的表达式中，可得：上式中，由于材料分散性所带来的不确定性效应会导致矩阵M，P，K的取值发生变化，状 态空间离散的特征矩阵实际上也是区间向量b的函数，故表示为A(b)，B(b)，E(b);此外，由 于结构响应计算采用迭代方法求解，材料分散性带来的时间累积效应不可忽略，于是z(t) -z(b，t)和综上，从不确定性传播分析的角度出发，受控结构的动力响应可 由一个域集F给予包络表达，即：其中，x(b，t)表示考虑材料分散性后主动控制系统的位移向量，φ[ ·]反映了位移向 量x(b，t)和状态向量z(b，t)之间的映射关系； 第五步:根据时变不确定性分析方法，构建出第四步中位移向量x(b，t)中各个元素 X1 (b，t)，i = l，2, . . .，n的非概率区间过程模型，并分别定义该区间过程模型的特征量函数， 具体包括：均值过程函数.半径过程函数以及任意不同时刻乜和乜的相关系 数函数Α?): 第六步:结合控制系统响应的许用值向量构建与第五步获得的位移向量x(b，t)区 间历程间的应力-强度区间过程干涉模型，并建立时变极限状态函数如下：基于首次穿越理论，通过时间离散化时段，定义如下穿越事件E1(Ic)发生的可能度指标： Pos{Ei(k)} =Pos{gi(b，kAt)>0 Hgi(b，（k+1) At)〈0} 其中，Pos{ · }表示事件发生的可能度，gi(b，kAt)>0表示受控结构在kAt时刻安全，此 时位移响应小于许用值，81(13，仏+1)八〇〈0表示受控结构在(1^+1)八七时刻失效，此时位移 响应大于许用值，符号"η"表示事件的交运算，k为计数指标，At表示时间增量； 第七步:遍历所有时间段内的可能度指标PosiEdk)}，建立主动控制系统非概率动力可 靠度指标：其中，Rs(T)为动力可靠度，T为完整响应历程，Pos(O)表示受控结构初始失效的可能度， 求解上式即可实现主动控制系统动力安全态势的有效评估。2. 根据权利要求1所述的一种基于首次穿越理论的振动主动控制系统非概率动力可靠 性分析方法，其特征在于:所述第一步中控制载荷向量f a(t)可表示为： fa(t) =BaU(t) 其中，BaSnXr阶权重矩阵，表示驱动器的作用位置，u(t)为r X 1阶驱动载荷向量。3. 根据权利要求1所述的一种基于首次穿越理论的振动主动控制系统非概率动力可靠 性分析方法，其特征在于:所述第二步中状态空间离散的特征矩阵分别表示为：其中，I为η X η阶单位矩阵，0为零矩阵，ITi表示总体质量矩阵M的逆矩阵。4. 根据权利要求1所述的一种基于首次穿越理论的振动主动控制系统非概率动力可靠 性分析方法，其特征在于:所述第三步中向量b均值和半径的具体表达式为：其中，和Ab请示向量b中第j个元素 h的均值和半径，j为计数指标。5. 根据权利要求1所述的一种基于首次穿越理论的振动主动控制系统非概率动力可靠 性分析方法，其特征在于:所述第四步中域集F通常无法显式获得，但借助区间传播分析方 法，位移向量x(b，t)的变化可表示为以下区间形式：其中，X1(I^t)为位移向量x(b，t)的区间集，^(b，t)和分别为x(b，t)的下界和上 界。6. 根据权利要求1所述的一种基于首次穿越理论的振动主动控制系统非概率动力可靠 性分析方法，其特征在于:所述第五步中均值过程函数蛘和半径过程函数的具 体表达式如下： 其中，Cov_, 表示协方差函数。7. 根据权利要求1所述的一种基于首次穿越理论的振动主动控制系统非概率动力可靠 性分析方法，其特征在于:所述第六步中事件Ei(k)表示为受控结构[kAt，（k+1) At]在时 间段内发生了一次穿越，At通常为一微小量。
【文档编号】G06F19/00GK106021909SQ201610325648
【公开日】2016年10月12日
【申请日】2016年5月17日
【发明人】王磊, 王晓军, 管闯闯, 樊维超, 苏欢
【申请人】北京航空航天大学