资源描述：
《基于吸收边界条件的有限元方法和离散切比雪夫展开的二维缺陷识别.pdf》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。
1、第22卷第6期计算力学学报Vol.22,No.62005年12月ChinesejournalofComputationalmechanicsDecember2005文章编号:1007-4708(2005)06-0762-05基于吸收边界条件的有限元方法和离散切比雪夫展开的二维缺陷识别常建梅1,冯怀平1,冯文杰#2(1.石家庄铁道学院土木工程分院0500432.石家庄铁道学院力学与科学系050043)摘要:结合考虑了吸收边界条件的有限元正演方法,通过离散切比雪夫多项式在特定离散点的展开逼近待识别函数,实现了无限区域内二维缺陷的识别O正问
2、题计算中引入的吸收边界条件虽然不能严格地模拟无限域对近场波动的影响,但能够以满足实际需要的精度模拟这一影响,且具有解耦特性,极大地简化了数值计算,在反演计算中,离散点取在切比雪夫多项式的零点处,使逼近时最大偏差最小,数值算例表明此方法的有效性O关键词:吸收边界条件有限元离散切比雪夫多项式缺陷识别中图分类号:0347文献标识码:Al引言2基本原理~散射场的计算和反演方法当今,缺陷识别问题是应用力学研究领域极为2.l基本原理活跃的前沿课题之一,无论在理论研究还是在工程如图1所示,对于平面应变的动力问题,在各应用中,都受到人们的很大关注,并
3、且已广泛地应向同性非均匀弹性介质中传播的波动方程为用于许多领域O如在材料和医疗等方面,利用超声8z(/8kuk)8(8uz8zu)]Fz=0uz波对材料及人体进行无损伤探测及医学成像在岩式中uz为位移分量,Fz为体力分量,/和是材料的拉梅系数,0是密度,0表示背景介质的密度,土工程中,利用弹性波反射信号做地质超前预报,02砼参数检测,桩基质量检测等在军事上,通过雷达888uz=,uz=2,假定=0,当体力为0时,8z8Iz8t及声纳技术,利用电磁波和声波的反射信号确定目设E8-wtkuk=-$6,则上式通过变分变换可表示为标的位置和形
4、状等1-5]OBorn直接反演方法基于线1$k22性原理,过程简单,运算速度快,但仅适用弱散射情-0n$=0(1)0(1)况O本文利用有限元方法计算散射波场,正问题的-0式中02-(1)=,E=06,n(1)=1/2计算中引入一种新的吸收边界条件,此吸收边界条00~60件虽然不能严格地模拟无限域对近场波动的影响,但能够以满足实际需要的精度模拟这一影响,且具有解耦特性,极大地简化了数值计算,在反演计算中通过离散切比雪夫多项式逼近待识别参数,离散点取在切比雪夫多项式的零点处,这样的逼近使最大偏差最小,从理论分析上可知与一般的多项式展开方法
5、相比,获得同样的精度所需的多项式项数要少得多O最终用最小二乘法建立目标函数,通过迭代实现二维缺陷的识别O收稿日期:2003-11-10修改稿收到日期:2004-01-20-作者简介:常建梅(1976-),女,讲师,硕士-冯文杰#(1967-),男,教授,博士-图1散射示意图Fig.1Sketchmap!or"cattering第6期常建梅,等:基于吸收边界条件的有限元方法和离散切比雪夫展开的二维缺陷识别763O当波在无界区域内传播时,应满足Sommer-k0=,假设c08Z11feld辐射条件lim~1+Jk0Z=0,运用有限6(1)
6、=,(1)=(2)1-O81]20(1)0c元方法,用有界区域代替无界区域时,在人为的边已知背景介质的参数为6B(1)和B(1),未知介界上须附加一定的边界条件[7],此条件应使得波的质参数为6(1)和(1),参考介质参数为6和将00,传播过程中不在边界上发生反射,在这方面有许多已知介质参数与未知介质参数间的差别通过参考专家做了研究,本文通过引入二阶BGT算子B近2介质参数6和正规化为下式[6]00:似Sommerfeld辐射条件[8],6(1)6B(1)66(1)=-=6(1)-6B(1)606085851B2Zs=+Jk0++Jk
7、0+Zs=O9/2(1)B(1)(8121)(8121)(1)6(1)=-=(1)-B(1)(3)00(7)从式中可以看出,只有当坐标位于未知的介质中通过变换,上式得时,66(1)和6(1)不为零值,将式(2)和式(3)代入式(1),得-8Zs8Zs=-=-28n81V(6VZ)+k0Z=0(4)12537k0假设总波场由入射波场和散射波场组成:22k0+2+J21k0+Zs+2(1+1k0){[41(41)]Z(1)=Zzn(1)+Zs(1)2-1k08Zs由Green理论,可得式(4)的解2+J2(8)(11)86}2Zs(1)=
8、k0[B(1)-1]Z(s)Gdxdy+将式(8)代入式(6),可将方程组写成矩阵的形式:s22V{[6B(1)-1]VZ(s)}Gdxdy+(S-k0T+A+JB){Zs}=(Sz-k0Tz){Zzn}(9)s式中66
