层合结构中SFEM自动形成的技巧.PDFVIP

层合结构中SFEM自动形成的技巧* 宋云连* （内蒙古工业大学建筑工程学院道交系, 010062, 内蒙古呼和浩特市） 摘要 本文简介了随机有限元法（SFEM）理论，并强调了它与传统有限元的不同。还介绍了 复合材料层合板结构的特殊性，在此基础上着重介绍了SFEM在层合板结构运用中，有限元 网格、节点编号、节点坐标，和局部平均单元网格自动形成的技巧，给出了对应的计算公 式，编写了相应的计算机程序。 关键词： SFEM；层合板；局部平均单元 一、引 言 因复合材料的各向异性及其构造特殊性，实际结构上存在很多随机因素，如与材料、 结构几何尺寸，及外载荷有关的随机因素等，这些都影响结构的安全使用。为此，本文 利用随机有限元法（SFEM ）对纤维增强复合材料层合板结构进行研究。 随机有限元法是将摄动法和有限元法相结合的一种数值计算方法，在假定随机变量 是微小扰动的情况下，它将有限元基本方程中各矩阵按随机变量展开成Taylor 级数得到 非线性方程。当假定随机变量均值和相关系数时，可将非线性方程转化成线性的确定递 归方程组，求解它便可得到位移解的各阶摄动系数，进而可得应力场的均值和方差。为 保证计算结果精度，SFEM 要求随机参数扰动微小。 利用SFEM 对结构分析时，须把计算模型离散成有限网格，并将节点、单元，及局 部平均单元信息等输入计算机，人工准备这些将是一项费时、费力又易出错的工作。本 文对这些实现了自动形成，给出了各自的计算公式，并编写了相应的计算机程序。 二、基于局部平均理论的SFEM 简介 1.1 基于正交局部平均理论的 SFEM 基本列式 本文SFEM利用了随机场的二维局部平均理论，其随机场模型的随机变量协方差矩 阵是满阵，采用特征正交化法可将其变为对角阵，下面给出对应的SFEM列式[1] n K V F K− V Ψ ( ) K V F ο e ∑ a a ο li K V F ο ο ο ， i l l ， 2 ο 2 (1) l 1 1 n 1 n n n n n V V Vare ( ) 其中： ， ，K ( ) ∑( ) F −F(K V Ψ)Ψ Var( )e − K V Var e Ψ 2 e e l ∑∑∑ ∑∑ ο 2 l l 2 2 ij ij il jl ο l i e il l l l 1 l 1 i 1 j 1 l 1 i 1 为整体刚阵K 的均值阵，F 、V 类同；n 是随机场 中向量总数；下标a 等表示变量对 a