板壳结构动力响应的随机有限元分析.pdfVIP

‘工程力学)增刊 1999年 ·473· 板壳结构动力响应的随机有限元分析 雷震宇 陈虬 (西敕逼大学应用力学与工程景，虚蕾6l∞31) 提要 本文利用随机有限元法对板壳结构进行弹性动力响应分析．分析中计及材料性能和几何尺寸 的随机性，随机性的表述是和用等参局部平均模型．研制了三角形板壳单元的计算程序．数值计算结果显 示本文的方法和程序均是有效的． 关键词扳壳结构，动力分析，随机有限元法 l、结构随机性的描述 结构工程的诸多方面往往包含着许多不确定因素。其中，结构自身属性的不确定性是非 常重要的一类，譬如：结构的材料性能和几何尺寸都有可能是具有随机扰动的参数，这类结 构通常被称为“随机结构”。这些因素往往被表述为空间的随机场．平板结构中材料性能(弹 性模量E、泊松系数u)和几何尺寸(厚度t)的随机性，是可以用二维局部平均模型来近似表 述的。 设叫毛_’，)是一个霉均值的二维连续平稳随机场．将该随机场用矩形的局部平均单元来 划分。设第i个单元的边长分别为TⅪ和k，定义该随机场的局部平均为： qG，力=去麟瞄G，脚 m 对于三维的壳体结构。一般需要将随机场的局部平均模型推广到三维。而本文采用一种 简化的方法。将壳体结构的随机场近似视为二维的来处理。具体做法是：将壳体结构的投影 面积(如球壳)或展开面积(如柱壳)近似作为随机场的有效面积。即不计及血率对随机性 的影响。在进{亍编程计算时，可单独对随机场进行划分，也就是说，对随机场建立一个与有 霄囊宇．女．t970．12出生。工掌博士 田塞自然基盘蠢助瑗目(j％7瓤，9) ·474- ‘工程力学，增刊 1999年 限单元独立的二维坐标系，这样就可近似处理随机场了。 为了实际应用的方便，可采用等参局部平均模型Il】，将离散随机场单元推广到任意直边 四边形．此时，定义单元e的局部平均为z 卟之孓抽由 ‰y啦 其中Ae为单元e的面积，Q。为单元e所占有的区域，它是任意的直边四边形．任意两个 局部平均随机场的协方差为： 叫％，吩)=i暑LL小一工’，y—y’岫。咖’ ④ 这里，a2是随机场a(x，y)的方差，p代y)是d(x，y)的相关函数．引入等参变换后的计算过 程可参阅文【l】． 为了节省计算量，还可以采用特征正交化技术，将协方差矩阵对角化为方差矩阵． 2摄动随机有限元动力分析列式 线性结构的运动方程为： 【M】(s(t))+【c《s(t斗+【K佟(t)}=妒(t)) (4) 如果材料性能和戚几何特性有随机扰动，质量矩阵【M、阻尼矩阵【c】以及刚度矩 阵【司就可能成为不确定的了。另外，载荷向量{F}也可能是随机扰动的，因此节点位移 也应具有随机扰动的性质闭。 假设结构的某一参数z是随机扰动的，对该参数建立随机场模型后，其扰动量可以用 一个随机小参数口来表示，即将z表示为确定部分和随机部分之和。 z=；(1+功 (5) 式中．-为z的均值；口是均值为零的随机场．它反映了参数z的随机性。对随机场离散 后，羽以化为随机向量f耐．将质量矩阵【M、阻尼矩阵【c】、刚度矩阵【明和载荷 列阵{F)在{口)的均值处按寨勒级数展开，并略去二阶以上项，经整理就得到了关于随 机变量f口}的 !三墨查堂!苎型 !!竺兰 ：!!!：． 零阶(均值)式： 【埘l。槲+【c】。¨+陋1。斜。={F)。 ∞’ L J0 ～JO 一阶变异式： 【乩料．+【c】。料。州荆．巾}』l l (7) {F