平面正交各向异性体材料参数识别算法及软件设计 ALGORITHM AND SOFTWARE DESIGN OF MATERIAL PARAMETER IDENTIFICATION FOR PLANE ORTHOTROPIC MATERIAL BODIES.pdfVIP

2010牟第5期 玻璃钢／复合材料 35 平面正交各向异性体材料参数识别算法及软件设计 黄立新1’2，刘玉印1，郭相武1 (1．广西大学土木建筑工程学院，南宁 530004； 2．广两大学工程防灾与结构安全教育部重点实验室，南宁 530004) 摘要：基于边界元法，本文提出了平面正交各向异性体材料参数识别算法。通过构造以测量位移与边界元计算的相应位 函数的问题。同时在识别过程中考虑了测量位移含正态分布的误差，研究测量误差对参数识别结果的影响。算例表明本文提 出的参数识别方法是有效的。 关键词：平面正交各向异性体；边界元法；参数识别 中图分类号：TP311文献标识码：A 文章编号：1003一0999《2010)05—0035—05 l 引 言 避免地含有误差。采用线性同余算法和中心极限定 正交各向异性材料(特别是先进复合材料)在 理，设计了一个正态分布的随机数发生器，用来模拟 现代工程结构中得到了广泛的应用，其准确的材料 测量误差，研究测量误差对参数识别结果的影响。 参数对工程设计与评价具有很重要的作用。融合测 算例表明本文提出的参数识别方法具有精度高、稳 定性好、抗噪音性能好、效率高的特点。 量技术、数值分析方法和优化技术的参数识别技术 是获取这些材料参数的有效途径…。材料参数识 2平面正交各向异性体参数识别的数 别属于反问题范畴，而正问题是反问题的基础。正 学模型 问题必须正确解决，才能提反问题。正交各向异性 0|： 平面正交各向异性体边界元法控制方程为【l 复合材料的正问题分析比较复杂，一般可借助于有 H(s)M(s)=G(s)f(s)(1) 限元法等数值方法加以解决口J。Kavanagh和 其中，s为参数向量，即s=[slls22 s12 Clough[31首次结合静态测量值与有限元技术，识别 s酯r；u(s)和t(s)分别为边界上的位移和面力； 正交各向异性复合材料及其结构的材料参数。除了 日(s)和G(s)影响矩阵。 测量静态值之外，许多研究人员又研究了测量动态 平面正交各向异性体在荷载作用下，可以根据 值的方法进行材料参数的识别，他们采用不同的数 方程(1)求解产生的位移Ⅱ和应力矿，这是边界元法 值分析方法作为正问题求解方法，Ayofinde【41等采 求解正问题过程。另一方面，根据测量位移Ⅱ反求 用瑞利(Rayleigh)方法；Deobald、Gibson口1、Moussu 材料参数s，即材料参数识别问题，这是基于边界元 和Nivoit【61等采用瑞利-里兹(Rayleigh-Ritz)方法； 法求解反问题过程。材料参数识别问题可以转化成 Cunha、Piranda…和Rikards’等采用有限元法。由 目标函数极小化问题。 于参数识别问题是非线性的，需要反复迭代求解，导 如果选取结构的位移为目标变量，则目标函数 致计算工作量很大。边界元法具有降维、精度高等 可以表示为： 独到的长处，并且域内点的位移和应力可按需计