各向异性层合板弯曲刚度系数反演分析 王晓纯 沈新普3 徐.pdfVIP

第 14 卷 第 2 期 复 　 合 　 材 　 料 　 学 　 报 V o l. 14　N o. 2 1997 年　 5 月 A CTA M A T ER IA E COM PO S ITA E S IN ICA M ay　　 1997 各向异性层合板弯曲刚度系数的反演分析 王晓纯　　沈新普 　　徐秉业 ( ) ( ) 北方工业大学建筑学院, 北京 100041 　　 清华大学工程力学系, 北京 100081 摘　要　本文提出了通过测量位移来获取复合材料正交异性板的弯曲刚度系数的反演方法。给出 了最小二乘法的迭代格式。详细讨论了不同的反分析模型对数值结果的影响以及数值解的收敛性 及灵敏度系数计算等问题。应用实例表明, 本文的方法是实用的、有效的。 关键词　复合材料, 弯曲刚度, 反演分析, 最小二乘法 中图分类号　 241. 2 O 　　用有限元法求解固体力学问题时, 最终归结为求解大型线性代数方程组 {F } = [K ]{u } 。 其中 [K ] 与材料性质、离散方式以及位移模式有关。对于金属材料, 由于各向同性的原因, 材 料性质以及材料参数的实验是比较容易的。对于复合材料, 问题要相对复杂一些。因为复合材 料的各种参数不仅与各组分材料有关, 而且与加工方式有关。对于一个完整的复合材料层合 板, 由于不能准确测量其铺层方式、材料间的粘合程度以及各层材料的性质, 因此, 也就无法计 算与各层材料性质有关的各种参数, 必须采用能够直接测试复合材料整体性质的方法[ 1 , 2 ]。已 有的一些方法只能解决单层板或者近似作为单层板的测试问题, 应用范围受限制。这样, 如何 直接测试层合板的各种物性参数方法的研究, 一直受到人们的重视。 [ 7 ] 　　基于上述原因, 本文根据系统辨识原理 , 针对复合材料层合板弯曲刚度系数的测试问 题, 提出用最小二乘递推公式反演弯曲刚度系数的方法, 给出了较为优化的反演模型, 讨论了 灵敏度系数的计算以及测点位置的选择对反演计算中收敛速度的影响等问题。几种典型算例 的结果表明, 文中提出的反演分析方法, 对解决复合材料层合板弯曲刚度系数的测试问题是行 之有效的。 1　反演材料参数的最小二乘递推公式 [ 3 ] 　　这里我们采用控制边界位移误差的方法反演材料参数 。首先设定, 针对一个具体计算模 型, 如果给出材料参数的初植之后, 可以进行有限元的正分析计算, 同时, 可以准确测量该计算 模型上若干个点的位移。 　　变量 {x }k 是第 k 次递推计算时的材料参数向量; {u }k 是与 {x }k 相对应, 由有限元的正分 - k 析计算出的位移向量; 设 {u } 是位移的测量值。反演计算的目标是不断改变材料参数 {x } , 使 　 k - 位移的计算值 {u } 与位移测量值 {u } 的差值的平方和为最小, 数学表达式为 　 　　本文于 1996 年 3 月 26 日收到修改稿, 1996 年 1 月 4 日收到初稿 © 1995-2005 Tsinghua Tongfang Optical Disc Co