层状盐岩体三维Cosserat介质扩展本构模型的程序实现

摘要

对层状盐岩体内的油(气)储存硐室进行稳定性计算时,如何选用恰当的本构模型来反映层状盐岩体的特点是至关重要的。针对我国大多数盐矿为多层盐岩的地质构造特征,基于宏观平均意义下的考虑细观弯曲效应的三维Cosserat介质扩展本构模型,结合FLAC3D软件的主要计算思路,使用VC++平台开发了该模型的FLAC3D接口程序,并编译成动态链接库DLL文件。经算例验证,该程序的计算结果可靠,不仅可以考虑不同力学特性的岩层先后破坏引起复合体单元的渐进损伤破坏,而且引入了夹层体积含量等统计参数后可大大节省单元网格数量,为大型层状岩体内的地下硐室稳定性数值模拟计算奠定了基础。