试探施加给定位移得到参考计算结果;然后采用应力函数法在有限元模型的范围内进行回归。回归计算中,顾及上述参考结果以及实测测点应力的最大、最小和平均值,在计算模型的边界上设定若干点的应力分量的值,并使这些边界条件同时参与回归计算;进而得到在建模范围内地应力的第1步近似结果。
第2步:利用第1步的结果计算出在建模边界上所有单元结点的应力,采用类似边界荷载调整法的做法把它们作为远场外荷载施加到有限元模型的边界上。此时由于边界荷载采用的是满足弹性力学中变形协调方程和平衡方程的应力场所求得的,边界上的正应力和剪应力可以共同保证构成平衡力系,并且解决了求剪应力场所遇到的困难。利用上述远场边界荷载,借助有限元程序计算建模范围内的地应力,然后将此结果在实测地应力的测点上与实测地应力进行比较,并且考虑在洞室群范围内地应力的分布情况结合工程概念考察回归计算结果的合理性。根据所得到的趋势,调整在第1步计算中施加在边界上的有限各点的力学边界条件,再次回归计算内部各点的地应力,并由此获得全部边界上由远场位移所引起的所有应力分量。
重复以上两步的分析计算,直到得到满意的应力场为止。
3、围岩稳定计算结果分析及相关问题
3.1 破坏准则与本构模型
在本次计算中,采用双参数准则中的Drucker-Prager准则来判断复杂应力状态下岩体的破坏。该准则是Mises准则的推广,因此又被称为广义Mises准则[2]。该准则认为,处于弹性区的应力满足式(1)时,则认为岩体处于破坏的临界状态。 上一页 [1] [2]
|