一种改进的基于代理模型的重频结构振动特征值的随机摄动方法
【专利摘要】本发明公开了一种改进的基于代理模型的重频结构振动特征值的随机摄动方法，该方法首先对重复频率结构的特征值进行了摄动分析，获得了在重复频率结构的刚度和质量等参数发生扰动和变化后，关于原重复频率结构自由振动特征值的一阶摄动量的矩阵方程。然后基于多项式混沌展开方法，构建关于重频特征值一阶摄动量的代理模型。结合摄动方法和基于多项式混沌展开的代理模型技术，提出了改进的重复频率结构自由振动特征值的近似计算方法，并基于该近似计算方法，进一步得到了重频结构特征值在参数扰动情况下的均值和方差的表达式。本发明解决了传统摄动方法中，重频结构特征值的一阶摄动量无法被显式表达，因而无法进一步研究其统计特征的困难。
【专利说明】
一种改进的基于代理模型的重频结构振动特征值的随机摄动 方法
技术领域
[0001] 本发明适用于重复频率结构系统的特征值分析，用以求解重复频率结构系统在经 受各种参数扰动的情况下，其特征值的统计学特性和变化范围，本发明可为重复频率结构 系统的特征值分析技术、结构分析、模型更新和设计优化提供指导。
【背景技术】
[0002] 在工程实际中常会出现重频系统或者密频系统，如平流层飞艇、飞机、大型火箭 等。特别是对于空间柔性结构，其自由振动的主要特点之一就是固有频率重复或者密集出 现。重频是实际工程结构中固有频率出现重叠的现象，在工程中，若不考虑重频的影响，将 重频模态处理为孤立模态，便会引起模态的遗漏。具有重特征值的系统称为退化系统，退化 系统有两个重要的特点：一是结构参数变化后，原来的一组重复特征值可能分离为一组非 重复特征值;二是结构参数变化后，特征向量可能产生跳跃现象，即相同编号的新旧两个特 征向量之差不是一个小量。产生这种现象的原因，主要是在结构具有重复频率的情况下，原 系统的特征向量的选取具有一定的随意性。
[0003] Kaminski M和Solecka M将结构特征值和特征向量用混纯多项式(PCE)方法进行 展开，研究了线性随机系统的受迫振动响应分析。X.W.YANG和S.H.CHEN将Pade approximation应用到矩阵摄动理论中，求出了特征向量和特征值变化量的表达式。关于结 构实模态特征值的统计特性，Qiu Z.P和Qiu H.C提出了直接方差分析方法(DVA方法），无需 已知或假定结构参数的相关系数矩阵，通过矩阵摄动理论和概率理论便能直接计算得到实 模态结构的随机特征值的方差。Zhao和Xu等人针对重频阻尼振动系统，提出了基于松弛组 合近似方法求解重分析问题的快速算法。通过将特征向量表示成为基向量和系数向量的组 合形式，避免了求解大规模方程组，简化了复杂求解运算。Palej和Krowiak利用解析的方法 进行了不同自由度下，重频弹簧质量块系统的模态分析。Vessel和Ram等人给出了两种含不 确定参数的阻尼振动系统重复特征值导数的计算方法。Wang和Zhang研究了将两个相近频 率处理成重复频率带来的误差。Chen和Jiao构造出一个弹簧质量块系统，并且提出能够通 过调整系统的参数，产生任意类型的重复频率振动系统。Xia和Shi等人将最大化结构的第 一阶独立或者重复固有频率为优化目标，研究了基于水平集方法的结构拓扑优化问题。 Reimherr进行了具有重复频率特征值的方程的回归分析。Fan和Xiang等人通过引入子域方 法，将随机摄动方法的处理范围扩展到了参数具有较大范围不确定性的领域以及重复特征 值的领域。
[0004] 矩阵摄动方法作为一种能够进行快速结构重分析和结构快速灵敏度分析的实用 工具，具有计算效率高、计算成本低和易于实现等特点，已经同时在理论基础和工程应用中 受到了广泛的关注并取得了长足的进步。但目前国内外学者对于摄动方法的研究大多集中 在基于泰勒级数展式的摄动方法上，当处理重频特征值时，由于特征值的导数不唯一，其分 析过程遇到困难。因此，目前关于使用摄动方法处理重复频率结构特征值问题的研究工作 还很少。
【发明内容】
[0005] 本发明要解决的技术问题是：克服现有技术的不足，提供一种改进的基于代理模 型的重频特征值随机摄动方法，用于计算重复频率结构在参数发生扰动和变化后，其特征 值的近似表达式，并基于该近似计算方法，进一步得到了重频结构特征值在参数扰动情况 下的均值和方差。本发明解决了传统摄动方法中，重频结构特征值的一阶摄动量无法被显 式表达，因而无法进一步研究其统计特征的困难。同时本发明方法只需少数的配点计算就 能建立比较精确的近似模型，也避免了蒙特卡洛方法大规模高成本的样本点计算问题。因 此本发明极大的方便了重复频率结构的动态分析、模型更新和优化设计等工作。
[0006] 本发明采用的技术方案为:一种改进的基于代理模型的重频结构振动特征值的随 机摄动方法，适用于具有重复自由振动频率的结构系统。首先根据矩阵摄动理论，推导了重 复频率结构系统在刚度矩阵、质量矩阵等参数发生扰动和变化后，关于其自由振动特征值 的一阶摄动量的矩阵方程。然后引入多项式混沌展开技术，建立针对重频结构自由振动特 征值一阶摄动量的代理模型。进一步结合矩阵摄动方法和多项式混沌展开技术，给出了重 复频率结构自由振动特征值的近似计算方法。基于该近似计算方法，得到了重频结构特征 值在参数扰动情况下的均值和方差的表达式，其实现步骤如下：
[0007] 第一步:根据矩阵摄动理论，针对具有重复自由振动频率的结构系统，推导了在其 刚度矩阵、质量矩阵等结构参数发生扰动和变化后，关于该系统自由振动特征值的一阶摄 动量的矩阵方程，这一步是建立重复频率结构自由振动特征值近似算法的基础，后续的算 法构建都是基于此方法实现的；
[0008] 第二步:针对第一步建立的矩阵方程中，方程的根即特征值一阶摄动量无法被结 构参数显式表示的困难，引入多项式混沌展开技术，将上述特征值一阶摄动量表示成为关 于系统中结构参数的混沌多项式展式的形式，建立关于重频特征值一阶摄动量的代理模 型；
[0009] 第三步:结合摄动方法和基于多项式混沌展开的代理模型技术，提出了改进的重 复频率结构自由振动特征值的近似计算方法，基于该近似计算方法，进一步得到了重频结 构特征值在参数扰动情况下的均值和方差的表达式。
[0010]所述第一步具体实现如下：
[0011] 步骤（11 )、给出退化系统具有m重的特征值本，和相应m个两两正交的特征向量^ 的特征方程，其中K〇和Mo分别表示重复频率结构系统的原刚度矩阵和原质量 矩阵；
[0012] 步骤（12)、根据矩阵摄动理论，建立关于重复频率结构自由振动特征值一阶摄动 m 量的矩阵方程]. I!州'-4[械.1财 1-)^/.=外％ (灸=H 讲），其中wVi:!，〗，…， i=i m)为原重复频率系统的特征向量，aj(j = 1,2,…，m)是待定系数，Ki和Mi分别是系统刚度矩 阵和质量矩阵的一阶摄动量，￥和<分别为重频系统自由振动的原特征值和扰动后的一阶 摄动量；
[0013] 所述第二步具体实现如下：
[0014] 步骤(21)、针对结构参数u2,…，U茜足高斯分布的情况，使用Hermite多项式展 式，将心表示为，⑷ ⑷ Zj. -1 ij =1 u =1 Z| -1 i-?-l i,=\ 式中C = ^，q,C^2，…，C;w表示待定的系数矢量，l=(U2，-_, In)为服从高斯分布 的随机变量，見(式，，…,戋):表示n次多维Hermite多项式；
[0015] 步骤(22)、将步骤(21)中的Y(U用有限项截断并用二阶Hermite多项式近似展开， 得到K⑷=c〇.2+J f备+￡ q,_(A 4) +乞，上式中n为随机变量的维数，C0,2， j~l z-1 i-\ j> i (^,2，(^,2，(^,2表示最高阶次为二阶的1161'1]1；^6多项式展开式的待定系数；
[0016]步骤(23)、选择好配点组合，依次将不同的配点代入原系统中，生成相应的系统响 应函数，最后建立如下的方程，本发明由于利用最高阶次为2阶的Hermite多项式，故配点应 当取为〃
的根，即jc = ，0，对于n维随机变量来说，采样点的数目，即随 机变量取值的组合数为3n。
L〇〇18」 b，U，…，b为米样点，N为米样点数目，s为待足糸数的数目，利用最小二乘法对上 式进行回归分析，就能计算出混沌多项式的展开系数(^(^^^，(^^，(^^，从而建立关于重 频特征值一阶摄动量的代理模型；
[0019]所述第三步具体实现如下：
[0020]步骤(31)、结合矩阵摄动方法和基于多项式混沌展开的代理模型技术，建立改进 的重复频率结构自由振动特征值的近似计算方法
[0022] 原重频特征值序列在受到扰动之后，得到一个新的按大小升序排列的特征值序 列，其中上标Y表示排序为Y的特征值；
[0023] 步骤(32)、根据步骤(31)给出的重频特征值的近似计算方法表达式，推导出重频 特征值在结构参数发生扰动和变化后的均值和方差。
[0024] 本发明与现有技术相比的优点在于：
[0025] (1)本发明解决了经典摄动方法中，重频结构特征值的一阶摄动量无法被显式表 达，因而无法直接研究其统计特征的困难；
[0026] (2)本发明方法只需少数的配点计算就能建立比较精确的近似模型，避免了蒙特 卡洛方法大规模高成本的样本点计算问题；
[0027] (3)本发明避免了基于泰勒级数展式的摄动方法处理重复频率结构的自由振动特 征值问题时，由于特征值的导数非唯一而造成的分析困难，因此本发明的应用范围更加广 泛。
【附图说明】
[0028] 图1为本发明一种改进的基于代理模型的重频结构振动特征值的随机摄动方法的 实现流程图；
[0029] 图2为本发明一种改进的基于代理模型的重频结构振动特征值的随机摄动方法的 实施例示意图，其中图2(a)为穹顶结构侧视图，图2(b)为穹顶结构俯视图。
【具体实施方式】
[0030] 下面结合附图以及【具体实施方式】进一步说明本发明。
[0031] 本发明提出了一种改进的基于代理模型的重频结构振动特征值的随机摄动方法， 其具体实施步骤是：
[0032] 第一步:根据矩阵摄动理论，针对具有重复自由振动频率的结构系统，推导了在其 刚度矩阵、质量矩阵等结构参数发生扰动和变化后，关于该系统自由振动特征值的一阶摄 动量的矩阵方程，这一步是建立重复频率结构自由振动特征值近似算法的基础，后续的算 法构建都是基于此方法实现的，下面给出具体的过程：
[0033] (1)给出重频结构自由振动的特征值方程
[0034]设重频结构(退化系统）的某个特征值#为m重特征值，相应有m个两两正交的特征 向量#(1 = 1，2,…，m)，并且满足方程：
[0035 ] K0wl =
[0036] 不仅如此，w1的线性组合也是该重特征值馬的特征向量，即： Mq = + (X.r, <j()n W
[0037] 「 i r = [w? ???,?,"] =W a
[0038]其中a是待定的向量。
[0039] (2)根据矩阵摄动理论，建立关于重复频率结构自由振动特征值一阶摄动量的矩 阵方程
[0040] 受到使用环境的变化和材料批次的差异等外界扰动时，退化系统的结构参数发生 变化，相应的有质量和刚度的变化发生，可以表示为： = i^+痛,
[0041] \
[0042] 扰动后的系统不再出现重特征值，即应该有m个不同的特征值^,i2,…，其对应 的特征向量为uUi = l，2,…，m)。类似于孤立的特征值情况，A1和u1仍可以表示成为￡的级 数，
[0043] +s2A!y + >--
[0044] tt ~ + G lty, + * ?
[0045] 根据摄动方法，将以上四式同时代入Ku = AMu，展开之后略去0(e3)项之后，比较e 的同次幂系数可以得到：
[0046] K0u[ +K{u!0 = A^Mqu[ + M{MQu}0
[0047] 将< +a2w2 +…+%>1^代入上式，得到： m m m
[0048] JC0?J +Kl'^ajwJ ~/^M0u{ +AiiMl'^jaJwJ +x1,M0^c(,ywi j-l j'=I j=l
[0049] 用[wk]T左乘上式两边，则有： m m. m
[0050] [wkf Kiw，ai - K\.wkf M0u\ M^w}(Xj Mow'，ai }~\ ；=1
[0051 ] 考察上式，同时根据[wY夏。B卜[:?丨f4aT碼#，戽 以及正交关系式[wHowbSkj，可得： m
[0052] y]([w/']/ K{w' -/!1![w/']/ Myy.t= {k = \,m) j=i .
[0053]不妨令 ak j = [ wk ] TMiwj，bk j = [ wk ] TKiwj，贝 lj 上式又可以改写成： m
[0054] YXbfi~^)&i=^a^ k = \,2,--,m
[0055] 上式的矩阵形式为=沒，其中矩阵D的元素为4%，(k，j = 1,2, ??? ,m) 〇
[0056] 上式的矩阵方程是一个m阶的标准特征值方程，求解此特征值问题可以得到相应 于重特征值4的一阶摄动々(./ = 1，2,…，），以及待定系数⑴。#仍然按照以前的规定从小 到大升序排列，若上式矩阵方程无重特征值，则所求的冬和&是唯一的，简便起见假定上式 矩阵方程无重复特征值。求解之，等价于求解：
[0057] 6tef[i>-^r，xw] = 0
[0058] 可以得到#(_/ = 1,2,…,m)，将其带回-= 可以得到q(j = l，2,…， 11〇，求得戽(_/ = 1，2，一，》)和(^(」=1，2，-.，111)之后，根据4=，.议，可以得到与重特征值4 相应的特征向量（./_ = /，/+丨，…，/ + m-l}。
[0059] 第二步:针对第一步建立的矩阵方程中，方程的根即特征值一阶摄动量无法被结 构参数显式表示的困难，引入多项式混沌展开技术，将上述特征值一阶摄动量表示成为关 于系统中结构参数的混沌多项式展式的形式，建立关于重频特征值一阶摄动量的代理模 型。具体实施步骤如下：
[0060] (1)针对高斯随机场，我们使用齐次Hermite多项式展开，随机响应Y(u可以表示 为：
[0061] n幻=(。+10私)+玄1>,1;卢2(以:)+si;fcv#(tA4) +… ；1 -i zj =1. i-, -1 ；|=1. ^ =1 i-i ~I
[0062]式中C = ，…)是待定的系数矢量，| = |n)为服从高 斯分布的随机变量，丑,,(^,4:,~4")表示11次多维此^1^6多项式。上式可以被截断，并且 用有限项来近似展开，不妨取s项，则上式可以简写为其中$为要求解的 7=0 系数，厂^.⑷=厂；(4，《:，?.、，备,)为j阶广义Wiener-Askey混纯多项式。这里Hn(U2，…，U 和r j(|) 对应，，％2，…，ev2; 3，…和& 对应。
[0063] (2)我们用二阶Hermite多项式展开构建代理模型，表示随机响应YU)，则有：
[0064] [⑷=c0.2 + 艺+ [c".2 (《' -1) + [ i=i i=i j=l j>i
[0065]上式中n为随机变量的维数，(30,2，(^,2，(^,2，叫,2表示最高阶次为二阶的1161'1]1；^多 项式展开式的待定系数。根据上式可以归纳出二阶Hermite随机响应多项式的展开式中待 定系数的数目为
，其中P表示多项式展开式的最高阶次，当p = 2时，s = (n+2)(n+ 1)/2。
[0066] (3)进行配点的选取，从而计算得出多项式混沌展式的各待定系数
[0067] 考虑到本文构造的PCE展开式的最高阶次为p = 2,故配点通常取为（p + 1)阶 Hermite多项式的根，即:r_=±^,CL我们不妨认为系统中的随机变量li(i = 1，2，…，n)均满足 高斯分布#(A,，di = l，2,…，n，那么每一个L配点的取值有三个片-a/Jct,.，片+ ^7;， 对于n维随机变量来说，配点的组合数目为3n。并且每一个配点组可以表示为：
[0068] Cp= (Cpi,Cp2, ,^.Pn) ,0= 1,2, ??? ,3n m
[0069] 将每一个选定配点组合依次代入=次％中，并将 M Mi，Ki表不成为：
[0071 ]然后求解方程式如下： 12345 -WfM，卞"=4% (灸=1，2,…，和） q=\ 2 把所得的值从小到大排列，得到一组解1|1=[41，2'￥1,...，％ 1]，重复以上步骤， 3 m 4 将每一个配点顺次代入方程-々ivy /w>;>,=次《， 则得到解向量： 5 j=i -戽4 4… ^ 4 4 ??? 4 4 U
[0074] XI, ^ 1^]= i：：:；： %私私…仏 ? ? ? * ^ ::. :; '?. ? ，' -巧木:；
[0075] 考察上式中任意阶的觯界=[况4私…?]，丫=1，2,…，m，可以建立以 下方程组： 乂仏）c(4)…/；-#)]「<]「々r n(4) r^) ... r^) cr ?J
[0076] :.. : . :;.：=： n(4) n(4) ??
[0077] 上式由3n个方程组成，要求解s个待定系数。对于PCE方法，通常来说3n>s，所以可 以基于最小二乘法进行回归分析，求出一组4…从而建立关于重频特征值 一阶摄动量的代理模型。
[0078]第三步:结合摄动方法和基于多项式混沌展开的代理模型技术，提出了重复频率 结构自由振动特征值的近似计算方法，基于该近似计算方法，进一步得到了重频结构特征 值在参数扰动情况下的均值和方差的表达式。具体实施步骤如下：
[0079] (1)建立改进的重复频率结构自由振动特征值的近似计算方法
[0080] 经过以上的求解过程，对于重频结构特征值的一阶摄动，其次序为y的值可以表 示为 n n n-\ n
[0081 ]代)=c。, + 艺(f 1) + 乞艺c#级 M: ni： j>r
[0082]其中| = (U2，…，U为服从高斯分布州/4,<)，纟=1，又…，《的随机变量，经过第 -?步的求角牛过程，CQ，q，Ci，q，Cii，q，Cij，q成为已知里。
[0083]这样我们只考虑到一阶摄动项，可以得到：
[0084] A/ =A^ +slf = A^ + s c〇r + ^c.y^. +'^cjjy^r - l) + XScy.;^^y 乂 i=l 1=1 i=l j>i ) 123 为了简便起见，在下文中，将上式的上下脚标y去掉，即： /- x n n n-\ n 2 k - A0 + ￡/>i1= Z0 + e C0 + + Xcs _ 0 + 、 /-I /-I 1-1 j>i j 3 (2)构建重频特征值在结构参数发生扰动和变化后的均值和方差表达式 [0088] 对上式求期望得到： E[X\ = E[k, + S^ = E -l) + ) \. /^1 1 } >i J
[0089] 「 n A v n-l fi ~ =aQ+se c〇+Jc.4+￡c;7 (^2 -1)+Y^LCM, i=l /=1 ：/-.! />/
[0090] 同样的，使用方差算符，并且考虑Hermite多项式的正交性，我们有：
[0092]综上，本发明结合摄动方法和基于多项式混沌展开的代理模型技术，提出了改进 的重复频率结构自由振动特征值的近似计算方法，并基于该近似计算方法，进一步建立了 重频结构自由振动特征值在参数扰动情况下的均值和方差的表达式。本发明解决了经典摄 动方法中，重频结构特征值的一阶摄动量无法被显式表达，因而无法直接研究其统计特征 的困难;避免了基于泰勒级数展式的摄动方法处理重复频率结构的振动特征值问题时，由 于特征值的导数非唯一而造成的分析困难;同时本发明只需少数的配点计算就能建立较精 确的特征值近似模型，避免了大规模高成本的样本点计算问题。
[0093] 实施例：
[0094] 为了更充分的了解该发明的特点及其对工程实际的适用性，本发明以图2的穹顶 结构系统为例进行重复频率结构的特征值分析验证。图2中的穹顶结构由52根杆和21个节 点组成，每一根杆件的截面积均为l〇cm 2。穹顶结构相对于Y轴和X轴对称，其杆件的弹性模 量、材料密度和泊松比的名义值分别为￡ = 2106?&4 = 78501^/1113# = 0.3。重要的具有代表 性的节点的坐标如表1所示。
[0095] 表 1
[0098]在本实施例中，将弹性模量和材料密度处理为随机变量，并且分别满足分布N (210，10.52)和N( 7850，392.52)。应用本发明提出的方法，得出此重复频率结构的自由振动 特征值的均值和方差如表2所示。为了验证本发明所提出的方法，同样采用Monte-Carlo方 法计算了本算例中的重频结构特征值的均值和方差。在随机数取值为1〇5时，由Monte-Carlo方法计算得到的结果和本发明得到的结果对比如表2所示。
[0099]表 2
针对如本实施例所示的穹顶结构，和经_典的蒙特卡洛方法相比，计算的精度是能够令人 满意的..不仅如此，本发明所提出的方法，在计算时只需要在原系统中对各随机变量进行有 限的几次配点计算，无需如统计方法那样对所有样本点进行重复计算。因此在计算效率方面 也有明显的优势。表3给出/本发明方法和蒙特卡洛方法在计算时间消耗方面的比较^
[0101] 由表3可以看出，本发明所提出的方法对于处理大型复杂空间对称结构的重频特 征值问题具有计算精度高、计算耗时少的显著优势。能够有效处理重复频率结构的动态分 析，模型更新，不确定性分析和设计优化等问题。经过二次开发，具有成为成熟商业软件的 潜力。以上实施例验证了本方法针对复模态结构随机特征值求解的可行性和优越性。
[0102] 以上仅是本发明的具体步骤，对本发明的保护范围不构成任何限制。
[0103] 本发明未详细阐述部分属于本领域技术人员的公知技术。
【主权项】
1. 一种改进的基于代理模型的重频结构振动特征值的随机摄动方法，其特征在于：该 方法适用于具有重复的自由振动频率的结构系统，并且考虑其结构参数发生变化和扰动的 情况，包括W下步骤： 第一步:根据矩阵摄动理论，针对具有重复自由振动频率的结构系统，推导了在其刚度 矩阵、质量矩阵结构参数发生扰动和变化后，关于该系统自由振动特征值的一阶摄动量的 矩阵方程； 第二步:针对第一步建立的矩阵方程中，方程的根即特征值一阶摄动量无法被结构参 数显式表示的困难，引入多项式混浊展开技术，将上述特征值一阶摄动量表示成为关于系 统中结构参数的混浊多项式展式的形式，建立关于重频特征值一阶摄动量的代理模型； 第=步：结合摄动方法和基于多项式混浊展开的代理模型技术，提出了改进的重复频 率结构自由振动特征值的近似计算方法，基于该近似计算方法，进一步得到了重频结构特 征值在参数扰动情况下的均值和方差的表达式。2. 根据权利要求1所述的一种改进的基于代理模型的重频结构振动特征值的随机摄动 方法，其特征在于:所述第一步具体实现如下： 步骤(11)、给出退化系统具有m重的特征值為，和相应m个两两正交的特征向量Wi的特征 方程其中Ko和Mo分别表示重复频率结构系统的原刚度矩阵和原质量矩阵； 步骤（12)、根据矩阵摄动理论，建立关于重复频率结构自由振动特征值一阶摄动量的 矩阵方程为 原重复频率系统的特征向量，a^j = l，2,…，m)是待定系数，Ki和Mi分别是系统刚度矩阵和 质量矩阵的一阶摄动量，為和絮分别为重频系统自由振动的原特征值和扰动后的一阶摄动 量。3. 根据权利要求1所述的一种改进的基于代理模型的重频结构振动特征值的随机摄动 方法，其特征在于:所述第二步具体实现如下： 步骤(21)、针对结构参数…，Cn满足高斯分布的情况，使用化rmite多项式展式， 将乂1'表示^巧 中C = (c。，，…，(?，…，，…)表示待定的系数矢量，C = ( Cl，C2，? ? ?，Cn)为服从高斯分布的 随机变量，巧,(写，每，…，写,)表示n次多维化rmite多项式； 步骤(22)、将步骤(21)中的Y(U用有限项截断并用二阶化rmite多项式近似展开，得到，上式中n为随机变量的维数，CO,2，Ci,2， Cii, 2, Cij, 2表示最高阶次为二阶的化rmite多项式展开式的待定系数； 步骤(23)、选择好配点组合，依次将不同的配点代入原系统中，生成相应的系统响应函 数，最后建立如下的方程，由于利用最高阶次为2阶的Hermite多项式，故配点应当取为向根，即;^=二^/^.()，对于n维随机变量来说，采样点的数目，即随机变量 取值的组合数为3。，其中向点，…，Cn为采样点，N为采样点数目，S为待定系数的数目，利用最小二乘法对 上式进行回归分析，就能计算出混浊多项式的展开系数C0,2，Ci,2，Cii,2，C化2,从而建立关于 重频特征值一阶摄动量的代理模型。4.根据权利要求1所述的一种改进的基于代理模型的重频结构振动特征值的随机摄动 方法，其特征在于:所述第=步具体实现如下： 步骤(31)、结合矩阵摄动方法和基于多项式混浊展开的代理模型技术，建立改进的重 复频率结构自由振动特征值的近似计算方法： 、''，. 原重频特征值序列在受到扰动之后，得到一个新的按大小升序排列的特征值序列，其 中上标T表示排序为T的特征值； 步骤(32)、根据步骤(31)给出的重频特征值的近似计算方法表达式，推导出重频特征 值在结构参数发生扰动和变化后的均值和方差。
【文档编号】G06F17/16GK105912508SQ201610278883
【公开日】2016年8月31日
【申请日】2016年4月28日
【发明人】邱志平, 仇翯辰, 王晓军, 何巍, 许孟辉, 李云龙, 胡永明, 郑宇宁, 张泽晟
【申请人】北京航空航天大学, 中国航空工业集团公司, 中国特种飞行器研究所