拉氏变换——边界元法中的数值反演和离散化.pdf

上 海 交 通 大 学 学 报 第 24卷 第 5、6期 JOURNALOFSHANGHAI][AOTONGUNIVERSITYVo1．24№5，6 1990 拉氏变换——边界元法中的数值 反 演 和 离 散 化 张永 元 张向晖 (上海交通大学 ) 摘 要 研究了边界元——扛氏变换法在三堆弹性动力学问遇应用中的若干技术细 节，如 Durbin数值反演技巧、单元离散方案选择厦频率对计算精度的移响等，获 得 了十分有益的结论． 关键词 边界元，拉 氏变换，弹性动力学，瞬杰响应，稳态响应 O 引 言 用边界元法解弹性动力学问题通常有两种途径t一种是时域方法，即在时域中用时间步 长选代法直接求得瞬态解J另一种是变换方法，即用拉氏或富氏变换后，求频域 中的 数值 解，再反演到耐域中求时间响虚或稳态解 ．至4目前为止，运 用拉 氏变换法求解三维弹性动力 学问题，还为数甚少．文献[1]最早用拉 氏变换法研究瞬态弹性动力学 问题， 但仅给 出了二 维问题算例 ，文献[2]研究了含球形空穴的无限体受谐和位移矢量波激励后 的应力集 中 系数 响应 |而文献[3]则研究了有刚性正方形地基的半无限体受单位谐和垂直位移作用后 的柔度 响应，其中，所受的激励是谐和时就无需进行拉氏变换的数值反演，在拉氏变换 空 间 算 出 的结果就是真实解，文献[4]是将拉氏变换法用于三维弹性动力学和三维动态裂纹问题 的初 步尝试，结出了有贯穿裂纹和表面裂纹厚板的稳杏解 |文献(5]对边界元法的几种途 径 在平 面问题中的直用进行了比较研究后认为t在相同精度的前提下，拉氏变换法最经济，而时域 法最耗机时． 本文通过受水平激励的柱体，着重研究拉 氏变换法在三维弹性动力学问题中的某些技术 细节的处理，拉氏变换数值反演精度与单元离散方案对计算精度的影响等问题． ■ ●I t1 1 边界元——拉 氏变换法的基本原理 1．1 基本方程 对于各向同性的均匀弹性体，在拉氏变换域中的运动方程为 本文 1660年4月蚰 日收曼．『 ·国家自然科学基盘资助项 目 ● I ， ● 张水元等：拉氏变换——边界无法中的数值反演和离散化 15f (c}一ci)，【…，+ci，【，-一s，【一·一吉。， (1) 本构方程 为 ∑．，一pE(c}一2ci)【，^，6 +cI(U…+ 【，，，)] 2‘) 边界条件为 U。一U。(p，s)，P，===芝 ．，=P (p，s) P∈厂 (3) 其 l}l：0。(p，)一p[ ．( ，s)+ 。(p，0)+ 【，，( ，0)]JU．(p，s)-Sc (u．(p，f))； ．，(p，)= )} )= (6 ∞)}尸I )一 l( )，c一尘 jcl一 是 拉氏变换参数 ． 利用加权残值法，可{拉氏变换域中任意一点的位移表达式为 u一一jru,PrdF+j，PU·~dF+j。0，，【-dO (4) 将式(4)代八式(2)，搿拉 氏变换域 的应力表达式为 -， 一 jrSk,tU~dF J 一P·d·厂+JD0dO (5) 对于边界厂上一点的位移表达式，通过与弹性静力学问题相 同的处理，搿 c-一，【一一U，Pdr+JrP·，【．J+J口0U7。dO (6) 对于光滑边界，c。 1／2；而 【，7和P ，是弹性