弹性动力学方程外边值问题的边界元区域分解算法研究.pdfVIP

408-415 力李与工租一 杜庆华院士八十寿辰庆贺文集.北京:清华大学出版社.1999.4 弹性动力学方程外边值问题的边界元区域分解算法 祝家麟 李开海 (重庆建筑大学，400045.重庆) 摘要 本文把边界元法引入区城分解算法来求解弹性动力学外边值问题。通过Fourier变 换，把瞬态弹性动力学方程转换为依核于振动频率的德态椭团型偏橄分方程后，采用t全 型区城分解界法求解 在子区城把问题转化为边界积分方程后用边界元离散.本文讨论了 边界积分方程的几种形式和葬法的几何收教性. 关健词 弹性动力学方祖;边界元;区域分解法 一、引 言 机械振动、地震动作用下结构的弹性动力学性质在工程技术上有重要意义。随着 计算机和有限元等数值计算方法的发展，使得定义在任意形状区域的弹性动力学方程 的数值求解成为可能。对于无界区域问题，使用边界元方法求解在无穷远处可以不设 人工边界，较之其它数值方法具有独特优点。对于有复杂形状区域的大型问题，把边 界元方法和区域分解算法结合起来求解，则又可以提高计算效率。它可以将大问题化 为若干个小问题求解，小问题又归结为求解边界积分方程，整个问题的计算可以使用 并行算法。区域分解算法允许在不同的子域选用不同的计算方法，把边界元方法引人 区域分解法中也有多种模式，可以选择不同类型的边界积分方程，还可以选用子区域 重叠型和非重叠型等交替算法。对于一个具体问题求解，则应根据问题本身的特点来 作出具体的选择。在区域分解算法的框架下，使用边界元仍然是希望在提高精度的同 时提高计算效率，为并行计算创造条件。但是用边界积分方程来计算动力学问题也要 注意特征值问题。对于弹性动力学问题的外边值问题，不存在特征值。但有界区域上 的内边值问题却存在特征值问题，用边界积分方程求解时，相应的边界积分方程也存 在特征值。事实上，稳态弹性动力学方程原始边值问题的特征值和对应的边界积分方 程的特征值之间有明确的对应关系，一方面我们可以通过求解边界积分方程来确定原 始边值问题的特征值。另一方面也必须注意到这样一个事实:仅管外边值问题存在唯 一解，但当通过某一类型的边界积分方程求解时，如果振动频率正好对应于某类内边 值问题的特征值时，这种边界积分方程可能无解或者解不唯一。这时，可以选择其它 形式的边界积分方程或者采用修正基本解等实用的数值处理方法[[31 本文不具体讨论边界元方法实现的细节，着重从区域分解法的角度，讨论求解平 面弹性动力学方程的第一类外边值问题 (Dirichlet外问题)的重叠型边界元区域分解算 一 法及收敛性。 之 产 - / 祝家麟等:弹性动力学方程外边值问题的边界元区城分解算法 二、关于平面弹性动力学方程 弹性动力学方程有多种表现形式，考虑方程 pAu+(A+p)graddivu+pj=pa (1) 在稳定振动状态下，位移场u(x,t)可用u(x,w)e的实部表示 ux(,t)=Re(ux(,ro)e“) (2) 将((2)代人((l)消去因子e0就得到依赖于振动频率w的稳态弹性动力学Navier方程 p4u(x,rv)+(A+p)graddivu(x,oc)+Pf+P)ou(x,w)=o (3) 其中A,,u是Lam常‘数，P为质量密度，.f是单位质童的体积力，这种形式