基于Zienkiewicz-Pande屈服准则的弹塑性本构模型在FLAC3D中的二次开发及应用

Zienkiewicz-Pande屈服准则是Mohr-Coulomb准则的改进,在p-q子午面和π平面上都是光滑曲线,不存在尖点,不仅在数值迭代计算过程中易于处理,而且在一定程度上考虑了屈服曲线与静水压力的关系以及中主应力σ2对屈服的影响。该准则已应用于多个水电站地下洞室开挖支护计算中,但并未被大型数值计算软件所采用。根据弹塑性力学原理,详细推导了基于Zienkiewicz-Pande屈服准则的弹塑性本构模型在FLAC3D中的增量迭代计算格式;在FLAC3D提供的二次开发程序接口基础上,利用VC++程序语言实现了计算过程,并编译成动态链接库文件进行加载和调用。将自编的本构模型对一理想圆形隧洞进行模拟,计算结果与Mohr-Coulomb准则一致,验证了计算格式的正确性。将此本构模型应用于一大型地下洞室群的开挖围岩稳定分析中,得到的塑性区、位移场、应力场分布与Mohr-Coulomb准则计算结果分布规律相符,进一步验证了模型的适用性和可靠性。     

地下洞室开挖爆破围岩松动圈的数值分析计算

根据数值计算提出了在地下洞室开挖爆破中,确定围岩松动圈的方法。给出了地下洞室开挖爆破的三维弹塑性损伤有限元计算方法。该方法论述了爆破荷载、重力荷载和开挖荷载的计算、施加以及迭代的方法,可以合理反应爆破对围岩稳定的影响和锚固支护施加的效果;结合岩石松动圈的实测方法,推导了判定围岩开始松动的损伤系数阈值公式,可以结合数值计算得到的洞周围岩损伤系数分布,确定在开挖爆破作用下围岩松动圈的范围。对实际工程的分析计算表明,根据数值计算判定的围岩松动范围与实测松动范围的规律基本一致,量值吻合较好,其结果为地下洞室开挖爆破的围岩松动圈的确定提供了有效的方法。     

大型地下洞室复杂地质断层数值模拟分析方法

根据大型地下洞室中复杂断层对围岩稳定的影响,提出了隐含断层单元的数值模拟计算方法。该方法将复杂断层单元隐含在岩体单元中,使得岩体单元的剖分不受断层切割的影响。使用这种隐含的断层单元,不但将受复杂断层切割的大型复杂地下洞室群的有限单元离散简单化,加快分析计算周期,而且能够有效地反映断层对复杂地下洞室围岩稳定的影响。通过对三峡地下厂房复杂断层结构的分析计算,证明该方法为复杂断层结构的分析提供了一种十分有效的计算方法和思路。