应甩微分算子理论建立基本解.pdfVIP

清 华 大 学 学 报 (自然 科 学 版) 第 3l卷 第 2期 JOURNAL OF TS1NGHUA UN1VERS1TY 应甩微分算子理论建立基本解 工程力学系 常亮明 杜庆华 文 摘：阐述了应用微分算于理论建立基本解的具体步骤．并应用它和求解常撒分方程 组的若干技巧建立了二维和三维弹性静力学，二维弹性动力学 ，Reissagr板 以及正交各 向异性薄板和平面向题的基本解．验证了所提出方法的有救性。 关键词；基本解．边界元．微分算于理论．弹性力学 、 分娄号 ：O175．3．O348．5 0 引 言 基本解是应用边界元法的出发点，但它的获得通常并非易事。建立基本解现行的常 用方法是积分变换法，即首先将偏微分 方 程 组 经付立叶变换化为常微分方程组，然后 求解此常微分方程组，最后对其解进行积分反演。这三个基本步骤中的积分反演常常相 当繁杂和困难，即使对最常用的弹性力学 Kelvin 解推导起来也较费力。 本文建议采用微分算子理论建立各种问题的基本解，并用此方法具体地建立了二维 和三维弹性静力学、Reissner板 以及各向异性的 Kirehhoff板和二维弹性静力学等 问题的基本解 ，验证 了方法 的有效性。 1 建立基本解的步骤 对于偏微分方程组，采用微分算子的张量形式，可 简记为 (告) )=一F ) (1) 对二变量问题有 ：(1， 2)， ，，=1，2j对三变量 问题 ， = l， 3)，i， ： 1，2，3。△为微分算子矩阵， ，为基本未知 向量 ，Ft为载荷 向量 ，重复下标表 示 求 和，逗号表示对 自变量的求导，下同。 基本解 U，J(，y)是 (1)式对于单位载荷的解，也即满足下列方程 的两 点 函 数： Air()u x，y)=一 一)“ (2) 本文于 1990年 6月 2Ft收刊 清 华 大 学 学 报 第 3l卷 式中：d(y— )为 DiracDelta 函数，y， 分别代表集 中力作用点 (源点)和域 内 该集 中力产生 的点 (场 点)，d 为 Kronecker记号。 按照微分算子理论 “，记 l^ 为 △ 的伴随矩阵，则有 △IJ△7}= (det△)d． (3) 定义 u ，)= ( )咖，) (4) 则由 (2) (3)知两点标量函数 e(x，y)是方程 (5)的解： (deta)P( ，Y)= 一d( — ) (5) 由于基本解是无限大平面或空间的特解，考虑到对称性，有： 对二维 (轴对称)， d(一 =d(r)／2ur 对三维 (球对称)， d(— )=d(r)／4 r (6) 式 中 r= 【—Y【表示场点与源点间的距离。因此通常更方便的是定义基本解 ， 、 uJ(，y)=△7( )P(r) (7) 、 … 相 应地 P(r)是下式的解； (deta)e(r)： 一d( — ) (8) 这样，应用微分