连续介质力学行为分析的静态拉格朗日方法-中国烟草学报.PDF

第 26 卷 第 7 期 岩石力学与工程学报 Vol.26 No.7 2007 年 7 月 Chinese Journal of Rock Mechanics and Engineering July ，2007 连续介质力学行为分析的静态拉格朗日方法 刘建华 ( 山东大学 土建与水利学院，山东 济南 250061) 摘要：通过对 FLAC 的研究，得出一个结论：拉格朗日方法无论从计算收敛的角度还是从物理意义的角度，基于 运动方程求解都是不必要的。提出了一种新的求解思想，建立了所谓静态拉格朗日方法：(1) 由节点不平衡力和 节点刚度系数确定节点位移增量；(2) 由节点位移增量确定单元应变和应力增量；(3) 由单元应力增量确定单元应 力状态，进而确定新的节点不平衡力。步骤(1)～(3)的过程重复进行，体系将达到平衡状态。给出了静态拉格朗日 方法的基本公式和基本处理方法，包括节点不平衡力、节点刚度系数、单元应变和应力公式以及介质离散、单元 应变和应力求解、介质开挖与支护、大变形问题等的处理方法。基于静态拉格朗日方法编制了计算软件，能模拟 多荷载步、介质开挖与支护等复杂问题以及介质线弹性、弹塑性、流变等多种力学行为。软件可视化的后处理功 能，可给出介质位移矢量图、塑性区分布图及节点位移、单元应力求解过程线。给出了静态拉格朗日方法算例， 并与 FLAC 计算结果作了比较，通过算例说明了静态拉格朗日法的特点。将静态拉格朗日方法应用于实际工程， 对黄河小浪底水利枢纽地下厂房大型洞室群围岩进行了力学分析。基于现场压缩试验，将岩体作为黏弹性介质， 采用广义 Kelvin 模型，并确定了岩体力学参数；阐述了黏弹性模型静态拉格朗日方法的计算过程；根据地下洞室 群开挖施工过程确定了数值计算的荷载步；建立了三维网格；基于实测计算了原始地应力场；数值计算给出了考 虑流变效应的岩体位移和应力状态，给出了岩体流变稳定时间，为评价岩体的稳定性和设计支护方案提供了 依据。概述连续介质静态拉格朗日分析方法的特点，该方法属于一种松弛或迭代方法，不需要通过节点运动 方程求解，也不需要节点质量、阻尼力、惯性力、节点运动速度、加速度、时间步长等这些物理力学量，使 求解方法简化。较 FLAC 方法的优点是，求解过程基本上是静态的，介质没有振动，或只有很小的振动，计算精 度提高，克服了 FLAC 求解过程中介质振动、计算结果物理意义不明确这一缺点；算例的结果表明，对于较 小的几何模型(单元和节点数较少) ，静态拉格朗日方法具有更快的计算速度，而对于较大的几何模型(单元和节 点数较多) ，FLAC 方法计算速度更快，在计算速度方面静态拉格朗日方法是否存在固有缺陷，尚需要进一步 研究。 关键词：数值计算；连续介质；力学行为分析；静态拉格朗日方法；弹塑性；流变 中图分类号：O 241 ；O 33 文献标识码：A 文章编号：1000 –6915(2007)07 –1438 –10 STATIC LAGRANGIAN METHOD FOR ANALYSIS OF CONTINUUM MECHANICAL BEHAVIORS LIU Jianhua (School of Civil and Hydraulic Engineering ，Shandong University ，Jinan ，Shandong 250061 ，China) Abstract ：Basic characteristics and solving processes of FLAC are summarized ；the advantages and disadvantages of FLAC are pointed out. Based on FLAC，a conclusion is obtained that it is not necessary f