平面正交各向异性体材料参数识别算法及软件设计.pdfVIP

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