岩土数值极限分析方法的发展与应用_郑颖人
青树---郑颖人院士
岩土领域的常青树---郑颖人院士（全文转载）记著名岩土工程专家、中国工程院院士、解放军后勤工程学院郑颖人教授内容摘要: 郑颖人，我国著名岩土工程专家、中国工程院院士、解放军后勤工程学院教授。
2004年，郑颖人和他的学生首创性地应用有限元强度折减法革新了隧洞设计计算方法，有力地推动隧洞设计水平再上新台阶。
以郑颖人为带头人的“岩土本构关系和军事地下工程稳定性”学术团队郑颖人院士在兰成渝输油管道重庆末站现场指导地基处理（左三为郑颖人）指导学生做实验人物小传郑颖人，浙江镇海人，1933年1月出生，1951年入伍，岩土工程与地下工程专家，现为解放军后勤工程学院专业技术1级教授、博士生导师，中国工程院院士，重庆市地质灾害防治工程技术研究中心主任，兼任空军工程技术顾问，中国岩石力学与工程学会、中国土木工程学会隧道与地下工程分会及防护工程分会、中国力学学会岩土力学专业委员会等顾问，重庆市科协副主席、重庆市土木工程学会名誉理事长等多种学术职务。
主编与参编国标、军标与地方标准10部，发表论文400余篇，出版专著10部。
先后获国家科技进步二、三等奖各1项，国家科技大会奖1项，军队和部委级科技进步一、二等奖9项。
培养博士、硕士研究生近百名。
3次被评为全军优秀教师，获总后勤部“一代名师”“重庆直辖十年建设功臣”“新中国成立60周年重庆杰出贡献英模”等荣誉称号，并获国土资源部全国地质灾害防治科技进步特别贡献奖。
郑颖人，我国著名岩土工程专家、中国工程院院士、解放军后勤工程学院教授。
曾经一度有些专家断言：“岩土是门经验科学，没有多少理论。
”但郑颖人不信邪。
他一生孜孜不倦、执着追求、大胆探索、严谨治学，为发展和完善岩土力学理论，辛勤耕耘半个多世纪，在岩土塑性力学基础理论、岩土数值极限分析方法与地下工程弹塑性与粘弹塑性理论方面硕果累累，为军内外解决了军事与城市地下工程、地质灾害防治工程、特殊土及其地基处理等军民两用工程一些重大技术难题！如今，年近八旬的郑老依然在岩土工程领域奋斗不息，他对记者说：“兴趣支撑着我走过这么多年，人老了，志趣不减，很多年轻人还叫我…长青树‟。
郑颖人 院士 演讲PPT茅以升与国际会议 有限元极限分析法在隧道稳定分析中的应用
1975, Zienkiewicz put forward FEM strength reduction. Through reducing c and of the rock and soil until failure happens, calculations show that the safety factor is the multiple of reducing the strength.
Pressure arch of deep-buried tunnel （普氏压力拱）
Shallow-buried theory by Terzaghi（太沙基）
Modern failure mechanism of tunnels
（现代的隧洞破坏机理）
► Based on elastic-plastic theories ► Rabcewicz crack wedge
（深埋隧洞破坏的现场监测）
（1）The anchor on the arch top bear small force, while that in the two sides bear large force. 顶部锚杆受力小 侧面锚杆受力大 侧向位移大 垂直位移小
（2）The average horizontal displacement of steel frame is 19mm .The average difference of settlement between arch top and arch foot is 9.8mm.
Buried depth:12m safety factor:0.7
When the buried depth is between 10m and 18m, the safety factor keeps constant basically and two crack surfaces appear respectively on the arch top and side wall.
岩土数值极限分析方法的发展与应用_郑颖人
DEVELOPMENT AND APPLICATION OF NUMERICAL LIMIT ANALYSIS FOR GEOLOGICAL MATERIALS
ZHENG Yingren1
，2
(1. Department of Civil Engineering， Logistical Engineering University， Chongqing 400041， China； 2. Chongqing Engineering and Technology Research Center of Geological Hazard Prevention and Treatment，Chongqing 400041，China)
1
引
言
[1-2]
土稳定性问题包括边坡稳定、地基承载力、土压力 等，其理论基础是极限分析理论，土体的极限分析 已有百年以上的历史， 在岩土 法起始于 1773 年的库仑定律，20 世纪 20 年代建立
极限分析方法
收稿日期：2012–03–13；修回日期：2012–06–05 基金项目：国家重点基础研究发展计划(973)项目(2011CB710606)；重庆市自然科学基金项目(2010BC8002) 作者简介：郑颖人(1933–)，男，1956 年毕业于石油学院储运专业，现任中国工程院院士、教授、博士生导师，主要从事岩土力学、岩土工程与地下 工程方面的教学与研究工作。E-mail：cqdzzx@
三剪情况 公式 J2 C
q ( 3 cos sin sin ) 3 1 3 tan sin 3 3 tan 2 4 3 tan sin
1 3 k
q ( 3 cos 3
莫尔–库仑
岩土工程：数值分析在岩体力学中的应用和发展.doc
岩土工程：数值分析在岩体力学中的应用和发展（一）数值分析方法的分类在岩石力学有关领域的数值分析方法应用中，主要使用的方法为有限元法，边界单元，离散单元法，拉格朗日单元法及块体理论等（二）有限元法原理及其应用要点原理：通过变分原理（或加权余量法）和分区插值的离散化处理把基本支配方程转化为线性代数方程，把待解域内的连续函数转化为求解有限个离散点（节点）处的场函数值。
应用要点：1.正确划分计算范围与边界条件2.正确输入岩体参数及初始地应力场3.采用特殊单元来考虑岩体的非连续性和边界效应（三）岩石力学问题的其他数值分析方法1.边界单元法有限元法是对问题的微分近似表达式给出了精确解，它实质上属于微分法。
与微分法相对应的是积分法，积分法所涉及的边界可包围整个问题域，而数值分析的离散化仅在边界上近似。
下图表示了在外部问题模拟时微分法与积分法之间的区别。
2.离散单元法离散单元法完全强调岩体的非连续性。
它认为，岩体中的各离散单元，在初始应力作用下各块体保持平衡。
岩体被表面或内部开挖以后，一部分岩体就存在不平衡力，离散单元法对计算域内的每个块体所受的四周作用力及自重进行不平衡计算，并采用牛顿运动定律确定该岩块内不平衡力引起的速度和位移。
反复逐个岩块进行类似计算，最终确定岩体在已知荷载作用下是否将破坏或计算出最终稳定体系的累计位移。
3.块体理论块体理论就是针对个性各异的岩体中具有结构面这一共性，根据集合论柘朴学原理，运用矢量分析和全空间赤平投影图形方法，构造出可能有的一切块体类型，进而将这些块体和开挖面的关系分成可移动块体和不可移动块体，对几何可移动块体在按力学条件分为稳定块体、潜在关键块体、关键块体。
此外，在计算方法上，还有半解析法、加权残余法以及松弛法中的经松弛法以及上述方法的耦合应用。
有限元强度折减法在土坡与岩坡中的应用_郑颖人
第23卷第19期岩石力学与工程学报23(19)：3381～3388 2004年10月 Chinese Journal of Rock Mechanics and Engineering Oct.，2004有限元强度折减法在土坡与岩坡中的应用郑颖人赵尚毅(后勤工程学院土木工程系重庆 400041)摘要通过有限元强度折减，使边坡达到破坏状态时，滑动面上的位移将产生突变，产生很大的且无限制的塑性流动，有限元程序无法从有限元方程组中找到一个既能满足静力平衡又能满足应力-应变关系和强度准则的解，此时不管是从力的收敛标准，还是从位移的收敛标准来判断有限元计算都不收敛，因此采用力和位移的收敛标准作为边坡破坏的判据是合理的。
对有限元强度折减法的计算精度和影响因素进行了详细分析，包括屈服准则、流动法则、有限元模型本身以及计算参数对安全系数计算精度的影响，并给出了提高计算精度的具体措施。
研究表明：采用徐干成、郑颖人(1990年)提出的摩尔-库仑等面积圆屈服准则求得的稳定安全系数与传统Spencer法的误差在5%左右，证实了其实用于工程的可行性。
在平面应变条件下则可采用摩尔匹配DP准则。
该文还将此法应用于岩质边坡的稳定分析，得到了岩质边坡的滑动面和安全系数，开创了求节理岩质边坡滑动面与稳定安全系数的先例。
关键词边坡稳定分析，有限元强度折减法，摩尔-库仑等面积圆屈服准则，精度分析，节理岩质边坡分类号 O 319.56 文献标识码 A 文章编号1000-6915(2004)19-3381-08APPLICATION OF STRENGTH REDUCTION FEMIN SOIL AND ROCK SLOPEZheng Yingren，Zhao Shangyi(Department of Civil Engineering，Logistical Engineering University， Chongqing 400041 China)Abstract With the c-tanϕ reduction，the FEM model of slope reaches instability，and the value of the nodal displacement just after slope failure has a big jump compared with that before failure. This actually means that no stress distribution can be achieved to satisfy both the yield criterion and global equilibrium. Slope failure and numerical non-convergence take place at the same time. So non-convergence in finite element program can be taken as a suitable evaluation criterion of slope failure. The influence on safety factor precision of different yield criterions，flow rule，FEM itself is thoroughly analyzed. At the same time some measures to improve the precision are put forward. A cone characterizes the Mohr-Coulomb surface in three-dimensional stress space with the vertices in deviatoric cross section. It brings difficulty to numerical analysis. For convenience this surface can be replaced by a smooth surface yield criterion，Mohr-Coulomb equivalent area circle DP yield criterion，which was proposed by professor Xu Gancheng and Zheng Yingren in 1990. The results show that the average error of safety factors obtained by FEM with Mohr-Coulomb equivalent area circle DP yield criterion and by Spencer method is about 5%. The average error of safety factor obtained by FEM with the plane strain Mohr-Coulomb matching DP yield criterion and by Spencer method is about 2%. The strength reduction FEM can also be used in the jointed rock slope. Through a series of case studies，the applicability of the proposed method is clearly exhibited.Key words slope stability analysis，strength reduction FEM，Mohr-Coulomb equivalent area circle yield criterions，error analysis，jointed rock slope2004年3月27日收到来稿，2004年5月7日收到修改稿。
郑颖人院士讲座学习体会--岩土类材料的空间莫尔圆理论及其屈服准则
线不是条直线，求无围压时的c、φ 值，
试验方法：做无法向力的直剪试验，得 到c值，按单轴抗压强度做莫尔圆，
从c点出发做莫尔圆
切线得到φ 角，
以C25砼为例得 c=3.2MPa, φ=61.3°
用数值极限分析法，验证抗剪强度的准确性。
混凝土物理力学参数
围岩类别
C25 C30 C35 C40
国内主要有姚仰平提出的基于 SMP准则的统一强度准则，俞茂宏[56]提出的基于双剪应力的统一强度准 则，以及高红-郑颖人提出的，基于传 统空间Mohr应力圆三剪能量强度准则。
σc =σcc 时，公式（5）、（6）、（7）给出了 以拉为正的等强度三剪能量屈服准则：
psinq(
3
3cosσ sinσsin)
3.89
C40 41.68 40.00
4.20
表明围压为零时c、φ值准确
三、传统的空间莫尔应力圆理论
软钢
σc= σcc
单向受拉(压)与双向受拉(压)强度相等 Tresca不考虑中主应力,Mises考虑中主应力
τ
σ 2
M ises
σ 3 O
σ 2
σ 1 σ
internal Tresca
σ 3
σ 1
空间Mohr应力圆 内接Tresca六边形
大(/)max的作用面（＝45+/2）
岩土材料破坏面 子午面上 能量屈服曲线
偏平面上 能量屈服曲线
岩土类材料一般以压为正
psinq(
3
3cosσ sinσsin)
2ccos 12 1 3
1 3tanσsin 33tan2σ 4 3tanσsin
，
c 256
有限元极限分析法发展及其在岩土工程中的应用研究
有限元极限分析法发展及其在岩土工程中的应用研究【摘要】有限元极限分析法适用于岩土工程的设计与分析。
笔者在本文中，主要介绍了岩土工程安全系数、方法和失稳判据等，以及有限元极限分析法在土坡、土基扩大以及基岩边坡基岩的应用，实现革新设计方法的目标。
【关键词】有限元；极限分析法；岩土工程；应用研究在岩土工程中，极限分析法得到了良好的应用，但是由于这一方法需要做假设，而且求解的范围有限，所以方法的应用受到了很大的限制。
但是有限元数值方法，具有很强的适应性，但是由于无法计算出稳定安全系数F，所以其应用也受到一定的限制。
在本文中，笔者探讨了有限元极限分析法的发展，以及其在岩土工程中的应用。
1 有限元极限分析法的发展20世纪70年代中期，英国科学家Zienkiewicz首先提出了有限元极限分析法，并且在岩土工程极限荷载与安全系数的计算中进行了应用。
在随后的1980年代和90年代，这种方法在边坡及地基稳定性分析中也有了良好的应用。
不过，由于当时的技术条件有限，缺乏可靠、强大的大型有限元程序、强度准则等，致使计算精度不够，在岩土工程中没有得到广泛的应用。
20世纪末，关于有限元极限分析法，国际上又出现了多种相关的研究文章，研究的方向主要集中在有限元强度折减法求解均质土坡安全稳定系数F方面。
但是由于计算结果与之前的研究结果比较相似，所以逐步为主流学术界所接受。
一些学者认为，这标志着有限元强度折减法分析边坡的稳定性，进入了一个崭新的时期。
1999年，美国的D. V. Griffith等人用该方法分析了边坡的稳定性，创新点在孔隙水压力与模拟水位两方面，同时也对库水下降情况下的边坡稳定性做了分析。
而我国有限元极限分析法在20世纪末才开始，主要是在土坡分析中的应用。
21世纪初期，国内的一些学者在边坡稳定性的分析中，采用了有限元强度折减法。
这是国内比较早的研究有限元强度折减法的文章，研究的方向集中在基本理论及计算精度两方面。
随着计算精度的不断提高，逐渐被设计单位和岩土工程部门所重视。
边坡稳定分析的一些进展_郑颖人
第21卷 第4期 地 下 空 间 V o l.21 N o.4 2001年12月 U N DERG RO U N D SP ACE Dec.2001 文章编号:1001-831X(2001)04-0262-10边坡稳定分析的一些进展郑颖人,赵尚毅,时卫民,林丽(后勤工程学院,重庆 400041)摘 要:本文介绍近年在用极限平衡法与有限单元法,分析边坡稳定性方面的一些进展,在基于极限平衡的解析法上,我们导出了多阶斜坡的稳定安全系数与滑裂面角度的计算公式。
这对岩质边坡的设计有很高的实用意义。
导出了目前采用的各种条分法的统一计算公式。
对于非严格条分法,用一个平衡方程并假设条间力的作用方向,即能求得安全系数;对严格条分法,用二个平衡方程,并假设条间力的作用方向或条间力的作用点位置,就能求出安全系数。
统一式是一简单的迭代式,因而计算简便,并有很高精度。
提出了两种用有限元法求边坡稳定安全系数的方法:一种基于极限平衡法,对土质边坡采用圆弧搜索法,对岩质边坡采用在滑移面上布置节理单元的方法。
另一种采用有限元强度折减法,便于采用大型软件,是一种很有前途的求边坡稳定安全系数的新方法。
关键词:边坡;稳定性;极限平衡法;有限元法;边坡稳定安全系数中图分类号:TU457 文献标识码:A1 前言随着我国的改革开放,国内基础设施建设蓬勃兴起,尤其是西部大开发,为边坡工程的发展创造了良好条件。
然而,边坡的稳定分析,参数选用,至今仍是一大技术难题,特别是岩质边坡尚没有实用的分析方法。
近年来,我们在用极限平衡法与有限单元法,分析边坡稳定性方面作了一些工作,本文就是介绍这方面的进展。
2 极限平衡法的进展2.1 边坡稳定安全系数定义极限平衡法是边坡稳定分析中最常用的方法。
它是通过分析在临近破坏状况下,土体外力与内部强度所提供抗力之间的平衡,计算土体在自身和外荷作用下的土坡稳定性程度,通常以边坡稳定系数表示:收稿日期:2001-06-15作者简介:郑颖人(1934-),浙江宁波人,教授,博士生导师,从事岩土本构关系理论及数值分析方面教学与研究工作。
隧道稳定性分析与设计方法讲座之三_隧道设计理念与方法_郑颖人
对深埋隧道，按照散体力学压力拱理论，给出隧道 围岩的松散压力，并以结构力学计算方法计算衬砌的 安全系数，作为设计依据。但这种算法与深埋隧道实 际受力情况相差较远，也与当前隧道及时支护的施工 方法不相适应。拱形深埋隧道并不存在压力拱，其破 坏主要来 自 隧 道 两 侧，围 岩 主 要 承 受 形 变 压 力。对 于 施 工 中 ，隧 道 围 岩 块 体 脱 落 和 局 部 失 稳 ，也 应 按 块 体理论计算松散压力。浅埋隧道一般按松散压力计 算 ，采 用 浅 埋 隧 道 相 应 的 计 算 公 式 ，通 常 实 测 压 力 会 远小于计 算 压 力，计 算 偏 于 保 守。但 浅 埋 隧 道 还 应 考 虑 雨 水 、环 境 等 对 围 岩 强 度 降 低 的 影 响 ，尤 其 是 土 体隧道和具 有 不 良 裂 隙 的 岩 体 隧 道，浅 埋 隧 道 计 算 公式还有待改进。
现行的隧道深、浅埋分界是基于散体理论，以隧道 顶部上方能否形成压力拱为原则，按压力拱原理，深埋 隧道能形成压力拱可以自稳，因而深埋隧道一定是稳 定的；浅埋隧道按理是不稳定的，但按浅埋围岩压力公 式计算，围岩可能稳定，也可能不稳定。至于隧道深、 浅埋分界线的确定原则上是基于压力拱高度，而目前 并没有从理论上导出压力拱高度，因而一般采用经验 值。如普氏压力拱拱高按岩石坚硬系数确定，而国内 规范采用经验统计值，所以至今仍是人为经验值。按 这种分法深、浅埋分界高度不仅与洞形、洞跨有关，还 与围岩的强度有关。
Tunnel Design Idea and Tunnel Design Method
有限元极限分析法发展及其在岩土工程中的应用
有限元极限分析法发展及其在岩土工程中的应用摘要：有限元极限分析法实际应用于岩土工程中，能够对岩土工程的安全系统、失稳数据等做出判断，但是在应用的过程中，需要做出假设，并且求解范围相对有限，在应用上有一定的限制。
关键词：有限元极限分析法；发展；岩土工程；应用；在实际应用过程中，是需要做出假设并求解的，而且应用的范围有一定的局限性，这是有限元极限分析法应该创新的地方，在科技进步之下，对方法进行完善，让其适用的范围有所扩大，同时也推动在岩土工程中应用的价值。
1有限元极限分析法发展历程1.1有限元极限法最初的提出者是英国科学家，时间在20世纪70年代中期，这也是首次将有限元极限分析法应用于岩土工程中，计算出岩土工程额极限荷载及其安全系数。
在20世纪90年代，该方法又应用于边坡和地基的稳定性分析中，但当时收到技术限制，并没有较强大和可靠的元程序支持，计算的精度也不够，在岩土工程中的推广使用收到了限制。
1.2在20世纪末，国际又对有限元极限分析法做出了新的研究，主要以有限元强度折减法的求解上比较集中，计算结果和之前的结果仍然很相似，慢慢也就被学术界接受到，从此有限元极限分析法也就进入了一个新的发展时期。
直到20世纪末，有限元分析法才在我国开始应用，主要是应用于土坡分析上。
在21世纪初，我国学者分析边坡稳定性上，有效应用了有限元折减法，这也是我国最早对有限元强度折减法的应用，并在基本理论以及计算精度上做出了细致研究。
在这两方面，我国也得到了较好的应用，并向着长远发展目标推进。
1.3在研究方面，有限元强度折减法主要集中在安全系数与滑面系数方面，而有限元增量超载法主要是在地基极限车承载力方面。
这方面的研究文献虽然不多，但是却取得了可观的研究成果。
这两种方法，统称为有限元极限分析法，从根本上来说，均为采用数值分析方法求解的一种极限分析法。
在国际上，有限元极限分析法大都采用编数值分析程序比较多，而该方法的应用范围仅局限于二维平面土基与土坡分析中。
岩土数值极限分析方法的发展与应用
岩土数值极限分析方法的发展与应用一、本文概述随着科学技术的不断进步和工程实践的日益深化，岩土工程的数值极限分析方法在工程安全评估、优化设计以及风险控制等方面发挥着越来越重要的作用。
本文旨在全面概述岩土数值极限分析方法的发展历程、现状以及未来趋势，并深入探讨其在各类岩土工程中的应用。
本文将首先回顾岩土数值极限分析方法的起源与发展，梳理其从早期的简单理论模型到现代复杂数值分析技术的演变过程。
接着，文章将重点介绍当前主流的数值极限分析方法，包括有限元法、有限差分法、离散元法等，并分析它们各自的优缺点和适用范围。
本文还将探讨岩土数值极限分析方法在岩土工程中的应用案例，如边坡稳定性分析、隧道开挖模拟、地下工程安全评估等，以展示其在实际工程中的重要作用。
本文将展望岩土数值极限分析方法的未来发展趋势，包括技术创新、方法优化、多学科交叉融合等方面，以期为相关领域的研究人员和实践工作者提供有益的参考和启示。
通过本文的阐述，希望能够推动岩土数值极限分析方法在岩土工程领域的进一步发展与应用。
二、岩土数值极限分析方法的发展历程岩土数值极限分析方法的发展历程可以追溯到20世纪中期，随着计算机技术的飞速发展和数值计算方法的不断创新，岩土数值极限分析逐渐成为一种重要的研究手段。
其发展过程大致可以分为以下几个阶段：初期探索阶段：在20世纪50至60年代，研究者开始尝试运用数值方法对岩土体的极限状态进行分析。
当时主要采用有限元法等基本的数值计算方法，对岩土体的应力、应变和位移等进行了初步的探索。
这一阶段的研究虽然较为基础，但为后续的发展奠定了坚实的基础。
方法发展阶段：随着计算机技术的不断进步和数值计算方法的日益成熟，岩土数值极限分析方法在20世纪70至80年代得到了快速发展。
研究者开始尝试运用更加复杂和精确的数值方法，如离散元法、边界元法、有限差分法等，对岩土体的力学特性、破坏模式和极限承载能力等进行了深入的研究。
这些方法的出现极大地丰富了岩土数值极限分析的手段，提高了分析的准确性和可靠性。
岩土数值极限分析方法的发展与应用
4、应用实例
以某岩石边坡工程为例，该边坡由不同岩石组成，表现出明显的非线性特征。 采用非线性破坏准则下的有限元法进行极限分析，首先建立边坡的三维模型，考 虑岩石的弹塑性性质和节理面的接触特性。然后，通过加载模拟不同的工况条件， 得到边坡的应力分布和位移变化。根据分析结果，可以评估出边坡在不同工况下 的稳定性。
一、岩土工程极限分析
岩土工程的极限分析主要集中在研究材料或结构的响应以及其所能承受的最 大负载。它提供了理解和预测在特定条件下岩石和土壤行为的框架。这个分析通 常需要考虑许多变量，包括材料类型、压力、湿度、温度以及时间等。
二、有限元法
有限元法是一种数值分析方法，它通过将连续的问题离散化成一组有限个微 小元素（或单元）的组合来解决问题。每个元素都有其自身的物理和/或化学特 性，这些特性通过数学模型进行模拟。通过这种方式，可以模拟出大规模、复杂 的物理现象，从而得到更精确的结果。
岩土数值极限分析方法的发展 与应用
01 引言
03 应用领域
目录
02 发展历程 04 参考内容
引言
岩土工程领域涉及许多复杂的物理现象，如变形、稳定性、渗流等，这些现 象的分析与预测对于工程建设的安全性与可靠性至关重要。为了更好地解决这些 问题，岩土工程师们不断发展并应用新的理论和计算方法，其中岩土数值极限分 析方法就是一种非常重要的工具。本次演示将详细介绍岩土数值极限分析方法的 发展历程、应用领域、主要方法以及未来发展方向。
然而，目前对于非线性破坏准则下的岩土极限分析方法研究仍存在不足之处， 例如对于不同类型岩土体的非线性特性缺乏深入了解、模型的参数确定仍存在困 难等。因此，未来需要进一步开展相关研究，完善非线性破坏准则下的岩土极限 分析方法，解决实际工程中的关键问题。
我国岩石力学的研究现状及其进展
第19卷 增刊西安矿业学院学报V ol.19Suppl. 1999年9月 JO U RNA L OF X I'A N M I NI NG INST I T U T E Sept.1999我国岩石力学的研究现状及其进展杨更社(西安科技学院建筑工程系,西安710054)摘 要:论述了我国岩石力学的研究现状及其进展,回顾了岩石力学在我国的发展历史以及岩石力学专家们一些年来所取得的主要成果;总结了我国岩石力学与工程的发展特色,并对可预期的进展及其前景进行了展望分析。
关键词:岩石力学;研究;进展中图分类号:T U452 文献标识码:A 文章编号:1001-7127(1999)S0-005-07我国的岩石工程有着长时期的发展历史。
在古代,著名的都江堰水利工程和闻名全球、被誉为世界八大奇观之一的万里长城以及由北京直达杭州的古老运河等都是代表性的佳作。
在当时,先辈们凭借丰富的实践经验设计施工,还没有建立岩土力学的概念。
新中国成立以后,各项经济建设事业取得了极大的发展,同时,也遇到了许多与工程地质及岩土力学密切相关的技术难题。
如特殊的区域性构造地质、松散破碎复杂岩基、高地应力作用下的极软岩、大跨洞室围岩的大变形、水工隧洞群之间的相互受力作用、高陡岩坡的持续稳定、岩体内的不稳态渗流,以及“三下”(铁路下、水下和建筑物下)采煤等等工程建设中遇到的十分突出的问题。
交通、能源、水利水电与采矿工业各个经济领域的需要对岩石力学与工程学科在我国的发展起到了有力的促进作用[1],[2]。
从50年代末开始,我国有历史意义的大型水利水电工程设计勘测的大规模展开,为岩石力学的试验和理论研究以及实际的工程应用注入了巨大的活力[3],[4]。
80年代末,中国政府决定正式兴建长江三峡工程,更大量的岩石力学与工程问题摆在中国专家、学者们的面前,如长达6km、坡高最大达170m的永久船闸高边坡岩体开挖,其整体稳定性与变形机制、岩体流变与地下水渗流等等极为复杂多变的岩石力学课题[5]。
（郑颖人）极限分析有限元法讲座_有限元强度折减法中边坡失稳的判据探讨
第26卷第2期 岩 土 力 学 V ol.26 No.2 2005年2月 Rock and Soil Mechanics Feb. 2005收稿日期：2004-08-02 修改稿收到日期：2004-08-02作者简介：赵尚毅：男，1969年生，博士，从事岩土工程稳定性极其数值分析研究。
E-mail:********************文章编号：1000－7598－(2005) 02－0332－05极限分析有限元法讲座——Ⅱ有限元强度折减法中边坡失稳的判据探讨赵尚毅1，郑颖人1 ，张玉芳2（1.后勤工程学院 土木工程系，重庆 400041；2.铁科院深圳铁科岩土工程公司，广东 深圳 518034）摘 要：边坡失稳，滑体滑出，滑体由稳定静止状态变为运动状态，同时产生很大的且无限发展的位移，这就是边坡破坏的特征。
有限元中通过强度折减使边坡达到极限破坏状态，滑动面上的位移和塑性应变将产生突变，且此位移和塑性应变的大小不再是一个定值，有限元程序无法从有限元方程组中找到一个既能满足静力平衡又能满足应力-应变关系和强度准则的解，此时，不管是从力的收敛标准，还是从位移的收敛标准来判断有限元计算都不收敛。
塑性区从坡脚到坡顶贯通并不一定意味着边坡破坏，塑性区贯通是破坏的必要条件，但不是充分条件，还要看是否产生很大的且无限发展的塑性变形和位移，有限元计算中表现为塑性应变和位移产生突变。
在突变前计算收敛，突变之后计算不收敛，表征滑面上土体无限流动，因此可把有限元静力平衡方程组是否有解，有限元计算是否收敛作为边坡破坏的依据。
- 关 键 词：边坡稳定分析；有限元强度折减法; 失稳判据 中图分类号：O 319.56 文献标识码：AStudy on slope failure criterion in strength reduction finite element methodZHAO Shang-yi 1, ZHENG Ying-ren 1, ZHANG Yu-fang 2(1 Department of Civil Engineering, Logistical Engineering University, Chongqing 400041, China ；2 Shenzhen TieKe Geotechnical Engineering Co. Ltd., Shenzhen 518034, China)Abstract: Slope collapse and the slide body come into moving state from stable static state simultaneously, and are accompanied by a dramatic increase in displacement of slide body. Furthermore, the displacement is not a definite value, but an infinite increase. This is the definition of overall collapse of a slope. In finite element model, the slope reaches instability with the strength reduction, value of the nodal displacement just after slope failure has a sudden change compared to the one before failure. This actually means that no stress distribution can be achieved to satisfy both the yield criterion and global equilibrium. Slope failure and numerical non-convergence take place at the same time. An element stress reaching the yield criterion state not always means that infinite “plastic flow” occurred. It is determined by boundary condition. The plastic zone developed from slope toe to top not means the overall collapse occurred. On the other hand, the distribution of plastic zone was influenced by many factors such as Poisson's ratio, flow rule, etc. So non-convergence in finite element program can be taken as a suitable evaluating criterion of slope failure. Through a series of case studies, the applicability of the proposed method was clearly exhibited. Key words: slope stability analysis; strength reduction FEM; criterion of slope failure1 引 言随着计算机软硬件及非线性弹塑性有限元计算技术的发展，采用理论体系更为严密的有限元法分析边坡的稳定性已经成为可能。
岩土工程极限分析有限元法及其应用
岩土工程极限分析有限元法及其应用摘要：通过研究分析发现，将工程结构离散化是极限分析有限元法的核心内容，简单地说实际的工程结构是通过想象进行离散一定数量的规则单元组合体，然后分析这些组合，结果应用于实际的结构中，通过这种实践在一定程度上解决了工程建设过程中的问题。
因此，本文笔者将详细对极限分析有限元法进行分析阐述。
关键字：岩土工程；极限分析有限元法；应用引言自上世纪初，岩土工程的极限分析方法（包括极限平衡法、滑移线场法、上下限分析法）取得了较好进展，在实际工程得到了广泛的应用。
其中一些方法需要一些人工架设，一些方法的解决方案非常有限，这限制了该方法的开发和应用。
其中有限元法数值方法适应力较强且应用广泛，但在工程设计中，不能求出稳定安全系数 F 和极限承载力，从而限制了岩土工程中有限元数值分析方法的运用。
一、经典岩土极限分析法的发展及问题基于力学的极限分析方法，土体处于理想的弹塑性或者刚塑性状态，处于极限平衡状态，即土体滑动面上各点的剪应力与土体的抗剪强度相等或者滑动面上的作用力与抗剪力相等。
极限平衡状态下的土体有两个力学性质：第一是土体处于不稳定的状态，所以它可以作为一个岩土工程破坏失稳的判据；第二是岩土材料强度充分发挥，达到最大经济效益，因此，在岩土工程中常把土体极限平衡作为设计依据。
有两种方法可以将地基或土坡引入极限状态：一是增量加载，如地基的极限承载力；二是强度折减，如土坡的稳定安全系数。
经典极限分析方法普遍应用于均质材料。
极限状态的设计计算仅参考破坏条件及屈服条件，不需要参考岩土复杂的本构关系，从而大大简化了岩土工程的设计计算。
极限状态计算应满足以下条件：（1）屈服条件或者破坏条件。
（2）静力平衡条件和力的边界条件。
（3）应变、位移协调条件和位移边界条件。
目前主要采用以下4种经典极限分析法：上、下限分析法、滑移线场法、变分法与极限平衡法。
每种都具有各自的特点，但还有一些需作假定，如上限法、滑移线场法、极限平衡法等都需对临界滑动面作假定，不适用于非均质材料，特别是岩石工程强度的不均性，从而限制了极限分析法的应用，这正是极限分析法在经典岩土工程的缺陷。
郑颖人有限元极限分析法在边坡与隧道工程中的应用
各类事故，
坍塌占68%
（1）地质条件复杂
（2）勘察不到位
（3）机理不清
（4）设计欠科学
（5）施工不当
(６) 管理不善 矿山法施工事故类型分布图
62
软弱围岩地铁车站坍塌举例：
某地铁车站开挖跨度27m，最大32m，埋 深7~15m，软岩和软弱结构面，拱顶二衬 完成，开挖侧墙下半部时，发生长30m坍 塌，损失1千万。
39
监测结果：
（1）4根桩前后地层位移表面最大， 安装时变动大，目前最大50mm。
位移/mm
6#桩后监测孔监测滑坡水平位移曲线 60 50 40 30 20 10
0 1 4 7 10 13 16 19 22 25 28 深度/m
2008-9-13 20009-8-3
6＃桩后监测孔监测滑坡水平位移曲线 40
（2）桩上位移监测结果：
两根桩桩顶位移，最大3.8mm
（3）桩上推力测试，最大推力 90KPa，呈拱形
（4）桩上抗力测试，最大抗力 80KPa，呈三角形
B12号桩桩身推力分布图
压力/KPa
90
80
70
60
2008-9-13
50
2008-11-11
40 30
顺
层
顺 层倾
边 坡
状层
边状
坡边 屈
坡曲
边
坡
顺层滑动 屈曲破坏
直线型滑面
折线型滑面
11
具有两组平行节理面的岩质边坡
两组方向不同的节理，贯通率100%，第一组 软弱结构面倾角30度，平均间距10m，第二组 软弱结构面倾角75度，平均间距10m.
12
计算结果
计算方法
安全系数
与时俱进,应用先进科技提高岩土工程建设水平
础 理论 、数 值计 算方 法 与 基于 监测 与数 值计 算 结 合 的信 息化 设计 、施 工 方法 ，将 成 为今 后岩
４
第４ ２ 卷
郑 颖人 ： 与 时俱进 ，应 用先进 科技 提 高岩 土工程 建设 水平
第 ９期
与监测方法的进步 ， 也是近年发展的趋势 。 此外 ，
力无关 ， 称 为粘聚力 ； 二是 固、 液、 气 三相耦合 。
上述 基本特点 决定 了岩 土力学难度很 大 。
步 ，但一 些先进物探设 备多数来 自国外 。
１ ． ２ 地 基 与 基 础
岩 土工 程 的力 学 问题基 本有 三类 ：对于岩
我 国土地 广 阔 ，尤其是 特殊 土分布 范 围很
1956年毕业于北京石油学院现为后勤工程学院军事建筑工程系技术一级教授博士生导师重庆市地质灾害防治工程技术研究中心名誉主任兼任空军工程技术顾问中国岩石力学与工程学会中国土木工程学会隧道与地下工程分会及防护工程分会等顾问中国力学学会岩土力学专业委员会副主任重庆市科协副主席等多种学术职务先后在军事地下工程边坡工程地下工程围岩力学岩土塑性力学岩土极限分析方法区域性土等理论研究与应用领域取得了丰硕的成果发表论文四百余篇出版专著10部主编国标军标5部
势 ，但 在数 值计 算 方 法 以及相 应 软件 编制 研究
方 面 ，与 国外相 比有较 大差距 。
于工程 实践 ， 有 待进 一步发展 。 在桩基设计 理论 、
新桩 型 、新工艺方 面得到长 足发展 。
当前 国内室 内与现场 的测 试技 术有 了较大
改 革 开 放 以 来 ，我 国 引 进 和 自主 创 新
强 度可 分 为两 部分 ，一 部分 与颗 粒 间 的法 向应
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
• 1298 •
岩石力学与工程学报
2012 年
了极限平衡法，40 年代，又相继出现了滑移线场法 (特征线法)，50 年代又提出了极限分析的上、下限 法。极限分析法经过百年的发展已逐趋成熟。从工 程实践上看，极限分析法具有很好的应用效果，解 决了岩土工程的一些设计问题，尤其是岩土稳定问 题。但对复杂的层状、非均质岩土材料及不同工程 的复杂情况，这一方法往往无能为力。随着岩土力 学数值方法的发展，逐渐兴起了数值极限分析方法， 它既有很广的适用性，又有很好的实用性。1975 年， O. C. Zienkiewicz 等 提出了有限元强度折减法与 超载法，可以用数值方法求解材料的稳定安全系数 与极限荷载，在国际上的边(滑)坡稳定分析中广为 应用
摘要：固体材料受力后从弹性发展到塑性、再发展到破坏，表明屈服与破坏是不同的。本文简述了岩土材料的屈 服准则，并提出岩土材料的屈服准则体系，提出材料应力场中点破坏和整体面破坏 2 种概念及其定义，并用传统 极限分析方法定义材料的整体面破坏准则。回顾分析传统极限分析方法的求解特点，由此指出传统极限分析法与 数值分析法的不足，两者结合形成近年发展的数值极限分析方法，使极限分析方法的适用范围大幅扩大。就传统 极限分析法与数值极限法的内涵、特点做深入的分析，论证了方法的可靠性，同时指出数值极限分析法的优越性及其 存在的问题，最后举例说明数值极限分析方法在岩土工程中的边(滑)坡、地基以及隧道工程中的广泛适用性。 关键词：数值分析；屈服准则；破坏准则；极限分析法；数值极限分析法；极限荷载；稳定安全系数 中图分类号：O 241 文献标识码：A 文章编号：1000–6915(2012)07–1297–20
sin sin ) c cos
德鲁克–普拉格 岩 土 材 料
I1 J 2 k 0
k 2 sin 3(3 sin ) 6c cos 3(3 sin ) 2 sin 3(3 sin ) 6c cos 3(3 sin )
Abstract： Solid material develops from elastic to plastic then to failure after undertaking some load， which means that yield is different from failure. The system of yield criteria for geomaterial is discussed and two definitions which are point failure and surface failure in stress field are put forward. And the criteria of surface whole failure is defined through traditional limit analysis method. Simultaneously，the deficiencies for traditional limit analysis and numerical analysis are pointed out by reviewing and analyzing the solution characteristics of traditional limit analysis method. Based on the combination of these two methods，a newly developed numerical limit analysis method is built，which has enlarged the application rage of limit analysis. The significances and features of traditional limit analysis method and numerical limit method are studied，as well as the reliabilities. The advantages and disadvantages of numerical limit analysis are listed；and the wide applicability of numerical limit analysis in the slope(landslide) engineering，foundation engineering and tunnel engineering is illustrated. Key words： numerical analysis； yield criteria； failure criteria； limit analysis method； numerical limit analysis method； ultimate load；stability safety factor 工程中被广泛应用，并作为设计的依据。经典的岩
三剪情况 公式 J2 C
q ( 3 cos sin sin ) 3 1 3 tan sin 3 3 tan 2 4 3 tan sin
1 3 k
q ( 3 cos 3
莫尔–库仑
p sin
高 红–郑颖人
2c cos
1
引
言
[1-2]
土稳定性问题包括边坡稳定、地基承载力、土压力 等，其理论基础是极限分析理论，土体的极限分析 已有百年以上的历史， 在岩土 法起始于 1773 年的库仑定律，20 世纪 20 年代建立
极限分析方法
收稿日期：2012–03–13；修回日期：2012–06–05 基金项目：国家重点基础研究发展计划(973)项目(2011CB710606)；重庆市自然科学基金项目(2010BC8002) 作者简介：郑颖人(1933–)，男，1956 年毕业于石油学院储运专业，现任中国工程院院士、教授、博士生导师，主要从事岩土力学、岩土工程与地下 工程方面的教学与研究工作。E-mail：cqdzzx@
DEVELOPMENT AND APPLICATION OF NUMERICAL LIMIT ANALYSIS FOR GEOLOGICAL MATERIALS
ZHENG Yingren1
，2
(1. Department of Civil Engineering， Logistical Engineering University， Chongqing 400041， China； 2. Chongqing Engineering and Technology Research Center of Geological Hazard Prevention and Treatment，Chongqing 400041，China)
第 31 卷
第7期
郑颖人：岩土数值极限分析方法的发展与应用 表1 应力表述的屈服准则体系[20]
• 1299 •
Table 1 System of yield criteria expressed by stress[20]
单剪情况 材料类型 屈服准则 金属材料 屈瑞斯卡 公ห้องสมุดไป่ตู้ 屈服准则 米赛斯
p sin
[4-11] [3]
制作用，不会出现破坏。所以屈服并不等于破坏， 但屈服使材料进入塑性，并造成材料损伤。当塑性 发展到一定程度后，就会在应力集中的地方出现局 部裂隙，可称为材料的点破坏，对此人们还缺乏足 够的研究。继续加载后，材料的局部裂隙就会贯通， 直至材料中破坏面形成，发生整体面破坏失稳。虽 然目前还没有公认的材料整体面破坏准则，但传统 极限分析实质上已经提供了材料的整体破坏条件， 并在工程中应用。因此，可以求出工程的极限荷载 或稳定安全系数，这正是极限分析法的魅力所在。 2.1 材料的屈服准则 屈服是材料达到弹性极限进入塑性，是一个过 程，从初始屈服、后继屈服达到塑性极限。材料从 弹性进入塑性要通过屈服准则来判别，材料的屈服 准则目前已有很多，一般是依据某种理论或者是某 种试验现象而建立的，也有依据经验而建立的。对 金属材料有屈瑞斯卡与米赛斯准则， 岩土中有莫尔– 库仑准则(简称 M-C 准则)和德鲁克–普拉格准则 (简称 D-P 准则)。上述准则都可依据弹性力学理论 推出，因而材料的弹性极限既可用应力表述，也可 用应变表述，因为两者是一一对应的。由于莫尔– 库仑准则没有考虑中间主应力的影响，从而出现了 由真三轴试验获得的三剪应力屈服准则，如 H. Matsuoka 的 SMP 准则、Lade 准则等。近年来，高 红 等[19]依据能量理论与三剪应力矢量，采用 2 种方法 导出了三剪能量屈服准则，它是岩土材料与金属材 料共同的屈服准则，国际上各种著名屈服准则都是 它的特例。当不考虑内摩擦角时，即成为米赛斯准 则；当不考虑内摩擦角同时又不考虑中间主应力时 就成为屈瑞斯卡准则；如果考虑内摩擦角而不考虑 中间主应力，即简化成为单剪状态，就成为莫尔– 库仑准则；如又假定洛德角为常数时，就是德鲁 克–普拉格准则。而且得到的屈服面形状与国内外 真三轴试验结果一致，因而可以由此列出岩土与金 属材料的屈服准则体系[20]，见表 1。 当材料的强度极限曲线为直线时，上述准则都 可以按严格的理论导出，其他由试验拟合得出的准 则与按近似理论导出的准则不包含在内。由试验可 知，土体的强度极限曲线一般为直线，而岩石的强 度极限曲线是非线性的。但目前工程中都将岩石视 作线性处理，此时应将极限曲线分段视作直线处理， 以获得合理的黏聚力 c 和内摩擦角 值。强度极限 曲线为二次曲线的岩石屈服准则，可参见郑颖人和 孔 亮[2]的研究。
2
岩土材料的受力破坏过程分析
固体材料随着受力的增大，一般都是先进入弹
性状态，随后材料中有些部位达到弹性极限，即这 部分材料屈服进入塑性状态，然后由塑性发展直到 塑性极限，进入破坏
[17-18]
。屈服的本质就是材料中