加权函数法消除长期项中阻带宽度的讨论.pdfVIP

四川建筑科学研究 第34卷 第 1期 60 SichuanBuildingScience 2008年 2月 加权函数法消除长期项中阻带宽度的讨论 范君艳 ，王凤武 ，范建蓓 ，徐志扬 (1．浙江机电职业技术学院，浙江杭州 310053； 2．中交第三航务工程勘察设计院，上海 200032) 摘 要：应用摄动有限元法解决随机结构的动力问题会伴随长期项的产生。采用加权函数法消除长期项，关键是阻带宽度的 确定。通过数值方法讨论了阻带宽度的确定问题，取得了理想的数值分析结果。 关键词：阻带宽度 ；an权函数法；长期项消除；摄动随机有限元法 中图分类号：O174 文献标识码 ：B 文章编号：1008—1933(2008)01—0060—03 互(t) Ⅱ速度激励； 0 引 言 {E}——惯性指示 向量。 在摄动随机有限元法解决随机结构的动力问题 假定结构的某一参数 是随机扰动的，对该参 时，由于方程中摄动格式的使用，将导致一阶和二阶 数建立随机场模型后，可表示为如下形式： 灵敏度方程中出现共振因素并造成动力方程摄动解 = 三(1+R ) (2) 的精度变坏，在一些情况下甚至出现不收敛的结果。 式中 ；一 的均值 ； 这就是所谓的 “长期项 (secularterms)”问题。长期 R：— 的变异系数； 项是随机结构动力方程 由于摄动格式的使用而产生 — — 是均值为零的随机场，它反映了参数 的“副产品”，会使计算结果精度随时间的增长很快 的随机性，若随机场 是正态随机场， 变差 J。因此，为保证摄动随机有限元的计算结果 随机场离散后，可化为服从标准正态 有效和精度，必须消除长期项的影响。目前 ，人们提 [0，1]分布的随机向量 {}。 出很多方法在数值上消除长期项 ，其中，uu等 提 根据摄动随机有限元法，方程 (1)可以表示为 出的用加权函数消除长期项的方法，对于在线性随 如下的递推方程： 机结构动力响应中消去长期项，是 目前所认可的一 [眠]{ }+[C。]{ }+ 种较好的数值方法。 [ ]{，，0}=_[眠]{E (3) 加权函数法是通过给离散化的激励函数 Fou— [眠]掣 +[c。]掣 +[]掣 ： rier谱系数赋予一个加权函数，从而将产生长期项 的某个邻域内的频率项系数消去，以达到消除长期 一 ( + 。EC]lYl+ o 项的目的。按原文献[2]给出的邻域范围消除长期 项的效果有限，本文通过数值方法对邻域范围，即阻 {Yo}) (4) 带宽度进行了讨论，给出了阻带宽度的确定方法。 +[c。]a2 ] = 1 随机结构动力反应控制方程 具有 1／,个 自由度的随机参数结构动力反应方程