薄体各向异性位势边界条件识别反问题正则化边界元方法.pdfVIP

第29卷专辑 固体力学学报 V01．29嚣。№明 2IX)M 2008年12月 CH丁卜J-ESEJOURNALOFSOLIDMECHANICSDecember 薄体各向异性位势边界条件识别反问题 正则化边界元方法木 周焕林” 王剑 牛忠荣 (合肥工业大学土木与水利工程学院工程力学系，合肥，230009) 摘要针对二维薄体各向异性位势Cauchy边界条件识别反问题，提出了解析积分和截断奇异值分解联合正 则化算法。解析积分用于该薄体位势问题边界元法中几乎奇异积分的正则化。截断奇异值分解技术用来求解系统 方程。提出了波动曲线法来截取得到适用的奇异值，该法基于位势梯度计算结果．数值算例表明该算法的有效性 和精确性。 关键词反问题，边界元法，几乎奇异积分，解析积分，截断奇异值分解，波动曲线法 O 引言 界f上，位势和法向位势梯度边界条件都已知，但 几乎奇异积分正则化边界元方法求解薄体问题 在边界‘上，所有边界条件都未知，则为Cauchy 具有自由度少、精确度高的优越性。对二、三维弹 问题。 、 性力学和位势薄体问题，导出了几乎奇异积分的半 区域内点J，的边界积分方程可写成如下形式： 解析或解析积分算法Il圳。 反问题是边界元方法研究的重要内容之一【4】。 “(y)=*qdF—lrq*ud厂(2) 如果部分边界上的边界条件可以全部获取，而剩余 位势基本解为： 边界上的边界条件无法测量，可通过数值方法反求 解．此即Cauchy边界条件识别反问题。 ln吾， V7(3) 兀1胀“l 小孺12 数值反问题通常是不稳定的，测量数据的误差 会影响数值计算结果的精确性。正则化技术是反问 其中。 题研究的关键【51。截断奇异值分解算、法16’7J可用于弹 性力学Cauchy边界条件识别反问题的边界元法【8J。 陌l砘k磁=隐kkl21(4) 本文采用解析积分算法处理了二维薄体各向异 令： 性位势问题边界元法中的几乎奇异积分。系统方程 用截断奇异值分解技术求解。提出了波动曲线法截 …～=bkl2l七2t 1-1_／SII．312322]J 取奇异值。 1 各向异性位势问题边界积分方程 且设源点坐标为◇1，弛)，场点坐标为@l，jc2)，则r可 在二维位势问题中，当材料具有各向异性时， 表示为： 其控制方程为： 毛。筹％豪峨：象=。㈣对应胁法向矢鲂向位势黻基本舰 砀是各向异性材料特性系数．设区域几何边界 g’2碉1·虹掣竽业 为厂=正u正，f，‘≠囝，正n正=o．如果在边其中，捍。，刀：为边界外法向矢量方向余弦。 1) ·安徽高校省级自然科学研究重点项目(KJ2008A04 ·t通讯作者． E-mail：huanlin_zhou@163．COnL ·46· 固体力学学报 2008年第29卷 、_源点他‘j-边界厂上，则得位势问题边界积分 力_Jl!：