整体元能量矩阵病态特性的改善.pdfVIP

第三届全喝青年岩土力学与工程会议论文集，南豪，1999年4月 整体元能量矩阵病态特性的改善 唐秋生 (上海铁道大学上海200331) 摘要：本文分析整体元解的结果的不稳定性的原因，提出了改善整体元能量矩阵病态特性 的方法，从而提高了它的计算精度并扩大了它的适用范围。 关键词：整体元病态方程能量矩阵 l问题的提出 整体元方法是将整个研究对象作为一个单元，利用最小势陡原理及条件设值原理得到 垒场近似解析解。从理沦上说，它能求解酃些具有任 意几何形状及任意边界条件的问题”1。然而，实际上 要达到整体元理沦上能达到的适用范围，尚有相当一 部分阃题需要处理。以下的i十算就很能说明问冠，设 有一压力隧洞(如图1)，其内径为l，外径为2，受 径向压力q=10，隧洞的弹性模量E=2，泊抡比jFO， 设其位移函数为： “=∑∑如x“Y (1) I=0 2=0 盘 v=∑∑民x2，1 (2】 出，当册较小时，Ⅲ的增大确实能捷整体元解与精确值之间的误差减小．然丽当Ⅲ≥6时， 整体元解出现了不稳定性，通过大繁算例可以发现这种现象具有普遍性，它捷得整体元的 实用性大为降低，它限定了整体元精度的提高，同时也限制了它所能解决的阔题的复杂程 度．可以这样说，这个问题在没有得到彻底解决之前，整津元可能仅是橱窗垦的一种摆设。 卒文想就上述问题出现的原因及对策提出一些初步设想。 lB6 一 。r y m=2 忉=3 ，=4 小；5 卅。6 n1=7 m=8精确解 1 00 000—2 595 一】0 55 -476 -078 ～5 46 “51 —1712000 1 20 000 一1418—719 t56 一l 23 -467 -283 一132l _407 仉 40 000 一1240 q82 -265 ．4 19 -579 一l 98 —1082 -553 1 60 O00 一1062 -353 叫61 -819 —8 23 -385 -804 —8 13 1 80 O00 -885 -333 ～8 03 一ll70 一ll26 —8 37 -207 —9 2I 00 -0 2 000 —707 o21 -125t 一1307 一1374 ～1506 90 一1000 2整体元能量矩阵的特性分析 整体元算法最终归结勾解一线性方程组，这一点与其他计算力学方法类似．然而整体 元昕形应的线性方程组的系数矩阵有其自己独有的特点，由于它是由弹性体的能量表达式 对各待定系数求偏导而得，为了叙述方便，暂称其为能量矩阵，用符号￡，表示．本节将对 能量矩阵U的特性作一讨沧．相信随着讨论的深人，能够找到本文第～节所提问题的原因 及解决此问题的相应对策． 设一平面问题的位移函数取为： 品 “=乙4妒，(z，Y) (3) }：j 丑 r=2-A，仍(r，y) (4) …1 设上述位移函数已经满芷该。平面问题的已知位移条件(如果所设位移函数不能满足 已知位移条件，则可将已知位移条件代人所设函数通过对愿函数重新组合，使其既满足已