弹性动力学的相似边界元法研究.pdfVIP

6期 重庆 建筑 大学 学报 V(l.2三 2x00202 Dec N2o00.06 李番 JournalofChonggingJianzhuUniversity 文章编号:1006一73`29(2000)06一0001一04 弹性动力学的相似边界元法’ 程玉民’， 彭妙娟 “ (1.上海大学 上海市应用数学和力学研究所.上海 200072;z，上海大学 土木工程系，上海 200072) 摘要:讨论了弹性动力学边界元法中边界单元相似时单元之间的一些矩阵关来 建立了相 似边界元法的公式。在一组相似单元中。尺要求得一个单元的相应矩阵，通过比例关系即 可求得其它单元的相应知阵、然后通过迭加建立代数方程纽系数拒阵。与通常的每个单元 都各 自进行积分计算相比，本文万法可大幅度减少计葬t。 关 钮 词:弹性动力学;边界积分方程;相似单元;相似边界元法 中侄W t . 文献标识码:A 在边界元法中，要建立最终可求解的代数方程组，需要在所有单元上进行大量的积分运算。当 单元数目较多时，单元上积分的计算量将大大增加计算时间。 在有限元法的研究中，文 1〔〕讨论了单元相似时，相似单元之间单元刚度矩阵的关系。在无限元 法的研究中，文幻〔也讨论了单元相似时的情况。在边界元法中，若能建立相似单元之间相应矩阵的 关系，则可以减少大量的积分运算。 为此，本文对弹性动力学问题的边界元法讨论了边界单元相似时单元之间的一些矩阵关系，建 立了弹性动力学的相似边界元法。在一组相似单元中，只要求得一个单元的相应矩阵，通过比例关 系即可求得其它单元的相应矩阵，然后通过迭加建立线性代数方程组的系数矩阵.与通常的每个单 元都各自进行积分计算相比，本文方法可大幅度减少计算量。 1弹性动力学边界元法 对于均匀、各向同性的线弹性材料，在小变形条件下，运动微分方程为 (c,一cz)u，。(二，t)+Ciu,(二，t)+.f;(二，，〕=u;(z,t)二C-。 (1) 其中c和。2分别为弹性体内膨胀波和畸变波的传播速度;九为体力分量;口为弹性体所在的域，其 边界为r;2=(T7,z,z8)(对三维问题)或z=(二】,x2)(对二维问题)。 应力和位移满足如下边界条件 t;(z,t)二几f;;=户(z,t) 任r. (2) u;(z,t)一4;(z,t) 〔i.. 和初始条件 绮(了,0十)=ub(z) z任 月 (3) 风(z,0`)一vn;(z) 本构关系为 r;;(z,t)=pCc}一2ci)u-(二，1)u,(z,t)St;+C,(w(x,t)+u;,;(X,t))〕二E。 ta〕 其中r.和r。分别为已知面力和位移的边界，r,Ur.=r;s,;为KronekerDelta, .收稿 日期:2000-as一of 基金项目:国家自然科学基金资助项目(ssosoo;; 作者简介程玉民a(sss-〕男，山西人教卜授.博士导师博士.主要从事计算力学、大跨空间结构研究 2 重庆建筑大学学报 一 2-9盗 对式((1)一(4)作用Laplace或Fourier变换，即得相应的变换域中的方程，再应用加权残数法， 得到弹性动力学向题在变换域中的边界积分方程 CkQ()sci。一Jr:u4