含空洞无限平面奇异解与有限元耦合分析地下工程问题.pdf

第22港 第6期 庆建筑大学学报 VOL22 2000年12月 Journaloff 几三c N2o00.容 Chou朋inqJianzhuUniversity 文童编号:1006-7329(2000)06一0058一04 含空洞无限平面奇异解与有限元藕合 分析地下工程问题 ‘ 唐寿高， 曹志远 (同济大学 工怪力学与技术系 固体力学教育部重点实验室，上海 200092) 摘哭:提出一种411单直接的辐合方法，即含有空洞的无限平面奇异解与有限元机合来求解 无限平面问题。地下二程葬例表明，本方法正确反映 了人二边界外无限城对边界内结构的 弹性效应。 关 锐 词:奇异解;有限元;栩合;,s土板伤 中圈分类号扣31介“’ 文献标识码:A 用有限元分析无限域力学问题，常用的方法是取一足够大的人工边界，忽略该人工边界以外区 域的弹性性质，而仅对人工边界内部区域进行有限元网格划分并将人工边界作为固定边界条件来 处理。为了提高精度，可采用边界元与有限元祸合法求解I[ll。本文提出一种更为简单直接的辆合方 法，即将该人工边界以外的区域看成一含空洞的弹性无限平面问题，设该人工边界(即空洞边界)上 分布有，个结点(与内部区域有限单元在该人工边界上的结点相一致)，通过分析该无限域空洞边 界上的奇异解，求得含空洞无限平面关于洞边二个结点的刚度阵，再将此刚度阵的元素叠加到人 工边界内结构有限元总刚度阵中相应位置上去，然后按一般有限元法求解。计算表明，本方法正确 反映了人工边界外无限域对边界内结构的弹性效应，改进了纯有限元法分析此类问题因人工边界 而带来的误差。 1含空洞无限平面奇异解 应用数学弹性力学理论，可得含有单空洞的无限平面奇异解为C[z.a7 9(z,z,)=一(1+v)/(8n)(X+iY)In(z一z,)+gb(z,z,) (1) 0,(z,z,)!1/(8r)C(3一40(X一Y:)ln(z一z,)+(1+二)(X+,Y)z,/(z一、)〕+A(z,z,) 其中Zy为洞边荷载作用点(源点),gs(Z,z,)和OU(z,补)是关于洞边(为满足该边边界条件)的映射 函数，不难求得为: %(z,z,)=一A,In(z,一z)+ALlnz YU(z,z,)=一A,(In(z,一z)+lnz一(3一v)Rz/zz一R}/Cz(z,一Z))} (2) 式中A,=(X+zY)/(870a 2 含空洞无限平面刚度阵 在源点z，作用单位集中，在场点z产生的位移可用上述奇异解表示为 ·收稿日期 2000-05-01 荃金项目。上海市固休力学重点学科研究项目 作者筒介:唐寿高(1952-).男.七海人，同济大学教搜 博士.主要从攀结构工程与岩土工捏分析的计算机方法研究. ，J弋梦肛 号夏下丁 ~，-一r~户尸..，，~ 一 一 一 第6期 唐寿高子:含空洞无限平面奇异解与有限元拐合分析地下工程问题 。二ReC(3一4v)9+(二，:，)一(二，:*)一O,(z,z,))/2(; yr=一IM03一40R(二，^F)一砂 (二，z,,)一p,(二，z,)J/ZC (3) 式中()为‘关于场点的导数;V为材料剪切模量;下标I表示源点力的作用方间 等于I表示义二1 Y=0，等于2表示X=0.Y一1, 若在洞边布置。个边界结点，这些结点与洞内结构有限元相应结点重合 由结构力学中柔度阵 的定义，可由式(3)求得含空洞无限平面结构的总柔度阵各元素，并集成总柔度阵为 、