随机结构重特征值分析的递推随机有限元法.pdf

第19卷第2期 振 动 工 程 学 报 V01．19No．2 2006年6月 ofVibration Jun．2006 Journal Engineering 随机结构重特征值分析的递推随机有限元法 黄 斌1，刘文军2 (1．武汉理工大学土木工程与建筑学院，湖北武汉430070}2．同济大学土木工程学院，上海210092) 摘要：利用递推随机有限元方法研究了具有随机参数结构的重特征值问题。采用随机收敛的非正交多项式展式表 示未知的随机重特征值和随机特征向量，建立了和摄动法类似的一系列确定的递推方程，通过求解这些递推方程， 得到了重特征值的统计值。算例表明，同基于二阶泰勒展开的摄动随机有限元法相比，递推随机有限元法的结果能 在较宽的随机涨落范围内更好地逼近蒙特卡洛模拟结果。 关键词：随机振动；重特征值；有限元方法；非正交多项式展式；摄动法 中图分类号：0324；0241．82文献标识码：A 文章编号：1004—4523(20a6)02—0156005 (K—AM)9=0(1) 引 言 式中K和肘分别为结构的咒×竹维刚度矩阵和质 量矩阵，A和劝结构的特征值和相应的特征向量。 在实际工程结构如对称结构、复杂空间结构的 当结构刚度或质量为随机量时，方程(1)就变成 动力分析中，常会遇到重特征值问题。重特征值问题 了随机特征值方程。如果确定性的结构出现了重特 的研究对于此类结构的动态设计、参数识别以及控 征值，而且当结构参数为随机量时，每个样本结构具 制等方面具有重要意义。许多研究人员对重特征值 有相同的重特征值现象，那么本文就称方程(1)描述 问题特别是重特征值的灵敏度问题进行了研究。例 的问题为随机重特征值问题。 如，王文亮、胡海昌提出了重特征值的小参数法[1]， 刘中生等对胡海昌的小参数法进行了补充研究[2]， 2随机量的表示 费勤发等针对特征灵敏度分析中参数∥的选择问 题做了进一步讨论[3]。 2．1结构参数中的随机量 考虑到工程结构参数的不确定性，朱位秋研究 为表述方便起见，这里仅假设平面梁系或板系 了具有随机参数结构的重特征值问题，在文献[1]基 结构的弹性模量为空间随机场，那么它的刚度矩阵 础上，结合一阶随机摄动法，得到了重特征值的统计 可表示为 值[4]。然而，当结构参数随机量的变异性较大时，该 ^ 方法的有效性将会受到限制。对此，本文采用文献 K≈K。+≥：e，噩 (2) i=1 [5，6]提出的递推随机有限元方法，对具有较大随机 式中甄为均值参数对应的，zX竹维确定性矩阵， 涨落结构的重特征值问题进行了研究，并和基于泰 K为以X咒维矩阵，￡为独立随机变量。事实上，关 勒级数展开的二阶摄动随机有限元法作了比较。两 于质量、几何尺寸乃至边界支撑弹簧刚度等也可以 个算例结果显示，同二阶