支护阻力作用下围岩应力计算自然边界元方法.pdfVIP

支护阻力作用下围岩应力计算的 自然边界元方法1) 2) 赵慧明 董正筑 陈家瑞 臧彤 （中国矿业大学，徐州 221008 ） 摘要 本文首先研究了圆形巷道在支护阻力作用下围岩应力计算的自然边界元方法，然后研究了方 形巷道在支护阻力作用下围岩应力计算的自然边界元与有限元的耦合法，最后研究了自然边界元与 有限元的耦合法在 ANSYS 中的实现。计算实例表明自然边界元方法计算量小，精度高，收敛速度 快，尤其是自然边界元方与有限元的耦合法既应用了有限元方法适应较复杂边界的优点，又充分利 用了自然边界元方法适于求解无界区域且具有较高精度的优点。相同的网格划分下，耦合法能更快 地收敛到准确解。经 ANSYS 软件进一步验证，有限元方法要达到同等精度，需要更大的求解区域、 更精细的网格划分。自然边界元与有限元的耦合法在 ANSYS 中的实现又可使计算结果通过ANSYS 的后处理模块得到更直观的展示。 关键词 支护阻力，围岩应力，自然边界元，有限元，耦合法 0 引言 自然边界元方法是由我国学者冯康、余德浩等人提出并发展起来的一种新型边界元方法[1,2]。所 谓的自然边界归化，就是根据Green函数和Green公式，将微分方程边值问题归化为边界上强奇异积 分方程（或称之为超奇异积分方程，在本方法中称为自然积分方程），然后再化成相应的变分形式在 边界上离散化求解的一类数值计算方法。韩厚德[3]与余德浩几乎同时, 以及稍晚些时候许多外国学 者也分别把这一思想应用于其它椭圆边值问题。J. B. Keller和D. Givoli 在八十年代末也相对独立的 提出了Dirichlet-Neumann映射的表达式[4-5],即自然边界归化的思想方法。 自然边界元方法与传统的边界元方法相比，除具有传统的边界元方法的所有优点外，还具有许 多独特的优点：（1）由Poisson 积分公式可容易的得到自然积分公式。（2）在自然边界归化下，能量 泛函保持不变。（3）用自然边界归化得到的刚度矩阵是对称正定的，尤其边界是圆域时，刚度矩阵 还具有某种循环性。因此只要计算大约半行系数就行了。这与传统的边界元方法相比，计算量确实 大为减少。（4）自然积分方程由原边值问题唯一确定, 它正是此边值问题的解的互补的微分边值之 间的本质的关系, 也即狄里赫莱-诺依曼映射。这些优点使自然边界元方法具有许多理论上和数值计 算上的优越性。 自然边界元方法与有限元方法的耦合也由于二者基于同样的变分原理而具有独特的优越性。耦 合法的应用不仅避免了用有限元方法直接处理无界区域，而且解决了传统的边界元归化方式不能保 持问题原有的性质而产生的诸如精度不高等这样那样的问题，使自然边界元与有限元的耦合法可以 充分利用二者的有点，而避免二者的缺点[2,6-10]。 本文即研究自然边界元方法对支护阻力作用下圆形巷道围岩应力的计算，自然边界元与有限 元的耦合法对支护阻力作用下方形巷道围岩应力的计算，以及自然边界元与有限元的耦合法在 ANSYS 中的实现。 1) 基金项目：国家自然科学基金项目，中国矿业大学校科技基金 (0k060157) 2) 第一作者简介：赵慧明（1971－），男，山东曲阜人，副教授，博士，主要从事计算力学方面的研究。E-mail: lx_zhm@163.com，Tel ：0516 1 支护阻力作用下圆形巷道围岩应力计算的自然边界元方法 考察有光滑边界Γ的平面区域Ω上的平面弹性方程组的第一边值问题 r r μΔu +(λ+μ)grad divu 0， =Ω内 ⎧