利用C＋＋对FLAC中的双线性的莫尔-库仑及非线性的虎克-布朗屈服函数进行了复合， 复合后的线性及非线性屈服函数交界处的临界应力对计算结果的影响是一个值得研究的问题。针对圆形巷道围岩的塑性区、应力、位移及剪切应变增量进行了数值计算。计算结果表明， 随着临界应力的值的增加， 隧洞围岩的塑性区尺寸、被监测单元或节点的剪切应变增量及位移减小；塑性区内的被监测单元的应力增加。这是由于随着临界应力的值的增加， 本文提出的非线性本构模型中线性剪切屈服部分的屈服极限提高所致。在高地应力（莫尔-库仑本构模型不再适用）及侧
