Signorini传输问题的有限元—边界元耦合.pdfVIP

维普资讯 中 国 柳 荸 (A 辑) 第 3l卷 第 2期 SCIENCEINCHINA (SeriesA) 2001年2月 Signorini传输 问题 的有 限元．边界元耦合 胡齐芽 余德浩 (中国科学院数学与系统科学研究院计算教学与科学工程计算研究所，北京 l00咖 ) 摘要 对一类带 Signorini接触条件的非线性传输 问题 ，提出了新 的有限元．边界元 耦合框槊．为求解所得到的耦合变分不等式，设计 了一种 区域分解型迭代方法 ，并 对其做 了完全的收敛性分析． 关键词 sil啊 接触 有 限元一边界元耦台 变分不等式 迭代方法 收敛率 设n是 中的有界区域，r是其Lipschitz边界．假定r= 砸 ，其中 和 是r中 的非空 、不相交开集 ．以 表 r的从n指向n ：= 、历 的单位法 向量．假定 ，∈L2(n)和 0∈H言(r)及 z0∈ (r)是已知函数． 在 n 的内部 ，考虑非线性 的偏微分方程 div(p(f u1)·Vu)+f=0， (0．1) 其中P：[0， )一 [0， )是一连续 函数 ，且 ·p(f)关于 单调增． 在 n 的外部 ，考虑 Laplace方程 Au =0， (0．2) u满足无穷远条件 u()= 。+0(1) ( 一 *)， (0．3) 其 中 。是一给定 的常数 ． 在 上考虑传输条件 u I =uI u。I ，p( u )3 i=8n +tOIp, (0．4) n 在 上考虑 Signorini条件 u ≤‰ I， p(1Vu1) i=8n 『 I≤。， 0：p(Iu1)弗l·(+u。一u)l． (0．5) 这类问题有很强的实际背景，如可看做一个线弹性的无界介质与一个具非线性应力应变关系 的可变形体之间的两体接触模型 (见文献 [1～4])、与此类似且更简单 的界面问题见文献 [5～ 7]．有限元与边界元相耦合是处理这类 问题的最有效方法，它能充分利用有 限元法 (处理 2000~04—17收藉 *国家 自然科学基盘 (批准号 ：198O1030)、国家博士后基叠、中国科学院王宽诚奖励基盘、国家 自然科学重点基金 (批准 号 ：l933l023)和 国家重点基础研究专项经费 (批准号 ：资助项 目 维普资讯 第 2期 胡齐芽等 ：Signorini传输问题 的有 限元 ·边界元耦合 复杂方程)和边界元法 (处理无界区域)各 自的优点(见文献 [1，2，5—7j)． 文献 [2]对 问题 (O．1)一(0．5)介绍 了一类耦合 的有限元．边界元变分不等式 ，并得到了相 应离散 问题 的渐近误差估计 ．对该耦合框架 ，定义在 r上 的边界积分算子不能显式地表示 ． 因此 ，为得到离散系统 ，需求解许多边界积分方程 ，从而导致昂贵的计算 ．此外 ，由于该边界积 分算子将产生稠密的刚度矩阵，它会破坏有限元部分的稀疏带状结构 ，故对离散系统整体地求 解会相当困难 ． 本文提 出另一耦合框架 ，它可避免在离散化过程 中求解边界积分方程 ，而只计算一些奇 异积分 ，且相应 的逼近解 同样具有最优 的误差估计 ．为求解得到 的耦合变分不等式 ，本文设计 了一种全新 的迭