基于Z-P准则的煤巷围岩稳定性数值模拟研究.pdfVIP

|mm ~w 学丰 自然科学版 第201343年卷4第月2期 AVporI3230N14o．2 FLAC。。提供 的程序接 口，实现了粘弹性广义Kelvin模型的FIAC。。二次开发 ；姚仰平等”。]以修正剑桥模型 为基础实现了统一硬化模型在 FLAC。。中的开发 ；赵 同彬[4将 自定义 的B—K锚 固本构模型编译到可调用 的 动态链接库中，大大提高了数值软件 的运行效率。在数值模拟模型开发 中，在满足相关岩石强度准则准确性 的同时又满足计算机程序的实现要求 ，已成为岩土工程 中一个非常具有研究价值以及应用价值的新方 向。 本研究基于 Z—P(Zienkiewicz—Pande)准则，通过 FLAC。。所提供的UDM 程序接 口，进行二次开发 。并 采用开发的Z—P准则和 FLAC。。内置的M—C准则对木城涧煤矿二槽巷道进行数值模拟研究。 1 Z-P准则及其 FLAC二次开发 1．1 Z-P准则 Mohr—Coulomb(M—C)屈服面是 比较可靠 的，但存在尖顶和棱角 的缺点 ，致使计算结果往往不够精确 。 Zienkiewicz等 。为此提出了Z—P准则，在 7c平面上抹 圆了M—C准则 的六个尖角 ，如图 1所示。在 P—q子 午面上通过调整系数a来逼近 M-C准则 ，以消除奇异点，如图2所示 。 M—C准则 ／、L ／ Z．P准则 p P 图 1 平面上的屈服准则示意图 图 2 P—q子午平面上 的屈服准则示意 图 Fig．1 Yieldcriterionon plane Fig．2 Yieldcriterionon qplane Z—P屈服准则的破坏包络线[7]为 F — + a1 一 k+ 一 0。 (1) 其中：k一 (f一36a)sin9／3／(3一sin9) ；J8一一12sin。9／(3一sin9) ； 一~／．，2／g(Oo)； 一 (+ 2+ )／a；叽一2In29／(amSm )2；g(0 一 F 二二-一二二__ 五 _二 ；K一 分别为材料的洛德角、粘聚力及内摩擦角 ；a为待定系数 ，取值与材料的c，值和允许抗拉强度有关 ， 。 为主应力分量 ；J 为偏应力第二不变量 。 1．2 FLAC。二次开发 在 FLAC本构模型的二次开发中，主要工作是由当前时步的应力状态和应变增量计算下一时步的应力 增量和应力状态。材料单元处于弹塑性状态可 由屈服准则和应力状态确定是否屈服。若 F≥O，单元进入屈 服 ；若 Fd0，单元仍处于弹性状态。当材料单元处于弹性状态时，迭代过程 中弹性本构关系保持不变 ；当材 料单元处于塑性状态时，应该用弹塑性屈服流动法则进行修正，拉伸屈服 的单元也要进行应力修正。 FIAC∞二次开发的具体流程如图3所示 ，基于 FLAC∞所提供 的二次开发程序接 口，利用 VC++6．0程序 语言编写 z—P准则的计算程序 ，并编译成动态链接库文件 Z—P．dll，由主程序进行调用和加载，成功实现 Z—P 准则的FLAC。。二次开发。 作为 自定义本构模型内部注册过程的一部分 ，需要一个静态的实例 ，并给出头文件定义的所有私有变量 的初始值 ，一般赋值为 0；Properties()函数给出了本构模型中所有参数名称 的字符 串格式，以0结束 ；Get— Property()函数是返回各模型参数的值 ，返回值是按照序列数 的顺序进行的，从 1开始 ，且这个顺序与 Prop— JournalofShandongUniversityofScienceandTechnology