边坡稳定分析的一些进展.pdf

文章编号 10012831X20010420262210 边坡稳定分析的一些进展 Ξ 郑颖人,赵尚毅,时卫民,林丽 后勤工程学院,重庆 400041 摘 要本文介绍近年在用极限平衡法与有限单元法,分析边坡稳定性方面的一 些进展,在基于极限平衡的解析法上,我们导出了多阶斜坡的稳定安全系数与滑裂面 角度的计算公式。这对岩质边坡的设计有很高的实用意义。导出了目前采用的各种 条分法的统一计算公式。 对于非严格条分法,用一个平衡方程并假设条间力的作用方 向,即能求得安全系数;对严格条分法,用二个平衡方程,并假设条间力的作用方向或 条间力的作用点位置,就能求出安全系数。 统一式是一简单的迭代式,因而计算简便, 并有很高精度。 提出了两种用有限元法求边坡稳定安全系数的方法一种基于极限平 衡法,对土质边坡采用圆弧搜索法,对岩质边坡采用在滑移面上布置节理单元的方 法。另一种采用有限元强度折减法,便于采用大型软件,是一种很有前途的求边坡稳 定安全系数的新方法。 关键词边坡;稳定性;极限平衡法;有限元法;边坡稳定安全系数 中图分类号 TU 457 文献标识码 A 1 前言 随着我国的改革开放,国内基础设施建设蓬勃兴起,尤其是西部大开发,为边坡工程的发 展创造了良好条件。然而,边坡的稳定分析,参数选用,至今仍是一大技术难题,特别是岩质边 坡尚没有实用的分析方法。近年来,我们在用极限平衡法与有限单元法,分析边坡稳定性方面 作了一些工作,本文就是介绍这方面的进展。 2 极限平衡法的进展 211 边坡稳定安全系数定义 极限平衡法是边坡稳定分析中最常用的方法。 它是通过分析在临近破坏状况下,土体外力 与内部强度所提供抗力之间的平衡,计算土体在自身和外荷作用下的土坡稳定性程度,通常以 边坡稳定系数表示 第21卷 第4期 地 下 空 间 Vol . 21 No. 4 2001年12月 UNDERGROUND SPACE Dec . 2001 Ξ收稿日期 2001206215 作者简介郑颖人19342 , 浙江宁波人,教授,博士生导师,从事岩土本构关系理论及数值分析方面教 学与研究工作。 1994-2007 China Academic Journal Electronic Publishing House. All rights reserved. F s Σ c c′ tgΥ tgΥ ′ 式中 S抗剪强度; Σ实际剪应力; c、tgΥ 土体实际的抗剪强度参数; c′、tgΥ 土体达到极限状态时的抗剪强度参数。 边坡中最小的稳定系数称为边坡稳定安全系数,它表示了该边坡的稳定程度。 212 解析法 土质边坡的滑动面一般为曲线,但在某些情况下可假设为直线,如计算库仓土压力时。岩 质边坡的滑动面一般为直线或折线。 这种情况下常可用解析法求边坡稳定系数或安全系数,因 而解析法是一种有效的实用方法。对于常见的单阶竖直边坡,已有许多计算公式。本文推导单 阶斜坡与多阶斜坡的解析式。 21211 单阶斜坡的稳定性分析 如图1所示边坡,坡角为 Α1,坡高h,滑裂角 Η,土体重度为 Χ,内摩擦角 Υ,粘聚力c,则稳定 系数为 F cl tgΥ ∫ 1 0Ρd l ∫ 1 cΣ d l ch sinΗ GtgΥcosΗ GsinΗ ch Gsin2Η tgΥ tgΗ 1 2c Χh sinΑ sin Α1-Η sin Η tgΥ tgΗ 直接用式1求F最小值有困难,将三角函数用坡高h和水平距离x表示,由图1几何关 系得 sinΗ h h2x 2, tgΗ h x , G Χh 2 x-hctgΑ12 代入式1得 F 2ch2x 2 Χh2x-hctgΑ xtgΥ h 3 下面求稳定系数最小值,令其导数为零得 2 cΧhtgΥx 2 - 4chıctgΑ1-2Χh2tgΥ ıctgΑ1x - 2ch2Χh3tgΥıctg2Α1 0 4 解方程得 xm inhctgΑ1h 2c1 ctg2Α1 2cΧhtgΥ 5 上式就是边坡稳定系数最小时的x值,将其代入3式就可得边坡最小的稳定系数,即边 坡稳定安全系数。利用图1几何关系可得滑裂角 Ηm in Ηm inar ctg h xm in arctg 1 ctg Α1 2c1 ctg2a1 2cΧhtgΥ 6 arctg sinΑ1 cos Α1 2c 2cΧhtgΥ 3622001年第4期 郑颖人等边坡稳定分析的一些进展 1994-2007 China Academic Journal Electronic Publishing House. All rights reserved. 应当注意,式6中 Ηm in是边坡稳定系数最小的滑裂角,而不是极限状态下的滑裂角,只有 当式6中h为边坡的临界高度hcr时,才是极限状态下滑裂角。临界高度可按文献[2]求得 hcr 4c Χ sinΑ1cosΥ 1 - cos Α1-Υ 7 由式6和7可求得如下几个特殊情况下的结果 1当c 0时,Ηm inΑ1,由3式得F tgΥ tgΑ1; 2当a1 Π 2 ,Υ 0时,Ηm inΠ 4; 3当hhcr时,Ηm in ΑΥ 2 ,F 1; 当hhcr时,Ηm in ΑΥ 2 ,F 1. 15,所以该边坡是稳定的。 213 条分法 21311 各种条分法的统一计算公式 条分法有多种计算方法,因为条分法都需要做假设,由于假设不同而形成不同方法。非严 格条分法中要假定条间力的作用方向,即Β角才能求得稳定系数。 严格条分法既要保证整体平 衡,又要保证条块的力矩平衡,因而多了一个平衡方程,但同时也多了一个未知数,即条间力作 用点的位置,所以仍要作出假设。假设条间力作用点位置的方法就是J anbu法;假设条间力的 5622001年第4期 郑颖人等边坡稳定分析的一些进展 1994-2007 China Academic Journal Electronic Publishing House. All rights reserved. 作用方向是M orgenstern2P rice法。假设 Β角为常数是S pencer法,它有足够的精度,又计算简 便,下面推导条分法统一公式时,设 Β为常数或为 Αi的函数。 如图4所示,当不计外力时,沿着土条底部和垂直于土条底部的合力分别为 SiwisinΑi-∃Pi cos Β-Αi13 NiwicosΑi-∃Pi sin Β-Αi14 由莫尔-库仑准则可得 Si ci∃xisecΑi F Ni tgΥ F 15 将15式代入13式,得出 ∃1i wisinΑi-ci∃xisecΑiNitgΥF cos Β-Αi 16 使土条满足整体力平衡的条件是6∃Pi 0,由此得 F 6 ci∃xisecΑiNitgΥ 6 wisinΑ 17 联立14 16 17式 Ni wicosΒci∃xisecΑi sin Β-ΑiF cos Β-Αi - tanΥsin Β-ΑiF 18 F 1 6 wisinΑi 6 c∃xisecΑi cos Β-Αi wicosΒtanΥ cos Β-Αi - tanΥ sin ΒΑiF 19 式19既可用于严格条分法,也可用于非严格条分法。对于非严格条分法只要假定 Β值, 问题就变成静定。当 ΒΑi时是瑞典园弧法其值与Fellenius法相同 , 也是我国的不平衡推力 法。当 Β 0时为简化B ishop法,当 ΒΒi时是陆军师团法;当 Β Β1Αi 2 时为L o w e2K araf iath 法。对于严格条件法,多了一个土条力矩平衡的方程,但也多了一个力作用点位置hi的未知 力,因而平衡方程中有二个未知力。如图4,对土条度中点取矩。 PisinΒ ∃xi 2 P1∃Pi sin Β ∃xi 2 -PicosΒh1 tanΑi ∃xi 2 20 Pi∃Pi cos Βhi 1-tanΑi ∃xi 2 0 化简得 hi 1 Pi Pi∃1ih i ∃xi 2 tan Β-tanΑi1 Pi Pi∃Pi 21 因此,将14代入16中,消去Ni,求解 ∃1i,并化简得 ∃1i -ci∃xisecΑiFwisinΑi-wicosΑitanΥ cos Β-Αi [ F- tan Β-Αi tan Υ] 22 将22式代入21中,可求出hi 1得 hi 1 Pi Pi -ci∃xisecΑiFwisinΑ-cosΑitanΥ cos Β-Αi [ F- tan Β-Αi tan Υ] hi- 662地 下 空 间 第21卷 1994-2007 China Academic Journal Electronic Publishing House. All rights reserved. ∃xi 2 tan Β-tanΑi 1 Pi Pi -ci∃xisecΑiFwisinΑi-wicosΑitanΥ cos Β-Αi [ F- tan Β-Αi tan Υ] 23 由以上可知,hi 1是关于hi,Pi,Β,F的函数,我们可假设初设已知条件hi,P1的值为零,这 样,我们就可以通过假设不同的 Β 、F值来求hi 1,直至满足其最后边界值为零为止。通过这种 方法求出的解即满足了力与力矩平衡又符合边界条件,另外它还避免了非严格条分法中人为 的假设条间力方向 Β值所引起的误差。 通过假设条间力作用点位置hi,求解F所用的方法是J anbu法,而通过假设条间力方向 Β 求解F所用的方法有M orgenstern2P rice法和S pencer法。用严格条分法求解,不仅能解出F, 还能确定条间力方向 Β,据 Β还可反算出滑移面位置。 21312 算例 1验算用统一公式中非严格条分法计算安全系数的正确性 某边坡,其容重 Χ 25kNm 3, 粘结力c 42kPa,内摩擦角 Υ 17,坡高H 20m ,坡角 Β 30,求圆心座标xc 67. 139m ,yc 11. 240m ,圆弧深度D3 30. 530m ,土条分成20块时的圆 弧滑移面分别用统一公式及土质边坡稳定分析程序S TAB95使用手册计算各种方法的安全 系数值见表1。 计算结果表1 安全系数值F 方 法 用统一公式解 用边坡稳定分析 程序S TAB95求解 瑞典法ΒΑi1. 3411. 343 Fellenius法ΒΑi1. 3411. 343 简化B isho法Β 01. 3971. 400 陆军师团法Β 301. 4001. 395 L ow e2K araf iatp法Β 30Αi21. 3751. 380 不平衡推力法ΒΑi1. 3411. 343 2验算用统一公式中严格条分法计算安全系数的正确性 同上例,计算指定任意形状滑移面的安全系数F,分别用统一公式及土质边坡稳定分析程 序S TAB95使用手册计算,假设推力线倾角值为常数。计算结果见表2。 计算结果表2 用统一公式计算S TAB95计算 F 14202, Β 19. 275F 1. 4196 21313 对各种条分法的评述 由式17可见,当Ni大,F也大。而Ni是 Β与 Αi的函数,对式18求偏导得 5Ni 5Β -wi sin Αi ciLiFwitanΥcosΑiF cos Β-Αi - tanΥ sin Β-Αi F 2 24 7622001年第4期 郑颖人等边坡稳定分析的一些进展 1994-2007 China Academic Journal Electronic Publishing House. All rights reserved. 公式24为零,其解为极小解。当分母为零时,即有 F ciLiwitgΥcosΑi wisinΑi 25 式25为瑞典法的安全系数公式,可见瑞典法是极小解。令 ΒΑi,式24分母为1,当导 数为零时,也可得式25。 可见Fellenius法,与我国的不平衡推力法也是极小解,它们都具有同 样的最小安全系数。这在前人的算例中也可得到证明。 简化B ishop法设 Β 0,即考虑条间的水平力,由于水平力起主要作用,因而简化B ishop 法有较高计算精度,且计算方便。由于F与 Βi没有太多直接联系,因而陆军师团法、L o w e2 K araf iath法没有多少计算依据。 严格条分法完全满足平衡条件,因而精度高于非严格条分法。S pencer法设 Β为常数,计算 简便,又有足够精度,是一种值得推荐的方法。按作者对50个计算数据的统计,瑞典法安全系 数最小,简化B ishop法的安全系数平均高出瑞典法6◊ ～7◊,而S pencer法高出简化B ishop 法平均2◊ ～3◊ 。 4 有限元法 411 有限元圆弧搜索法 通过有限元分析可以得到每个节点的应力张量。 假定一个滑动面,即可用有限元数据给出 滑动面一点的法向应力 Ρn和剪应力 Σ 。根据莫尔2库仑准则可以求得该点的抗滑力。由此即能 求得滑动面上每点节的下滑力与抗滑力,再对滑动面上的下滑力与抗滑力进行积分,就可求得 每一滑动面的稳定系数S F. S F ∫ 1 0 cΡntgΥds ∫ 1 0 Ρ x-Ρy 2 sin2 Σxycos2Αds 假设滑动面为圆弧型,并不断变换圆心与半径进行搜索,就可得到最小的安全系数,即边 坡稳定安全系数。图5中出示了算例,算得边坡安全系数为1135,按照简化B ishop法为1134, 两者是很接近的。 图5 图7 412 用有限单元法求岩质边坡稳定安全系数 岩体是弱面体,岩体滑动一般是在结构面上,因而不能采用圆弧搜索法求滑动面,需要地 862地 下 空 间 第21卷 1994-2007 China Academic Journal Electronic Publishing House. All rights reserved. 图6 质分析与力学分析方法结合确定滑面。 对于一些简 单情况,如坡高不大的边坡可按外倾结构面确定滑 面;对于复杂岩体如何确定滑动面,还需要深入研 究。岩体强度是很难确定的。当岩体无裂缝时,强度 接近岩块强度;当岩体有裂隙时,应采用平均强度 即岩体强度 ; 但当有外倾结构面时,在结构面上应 采用结构面强度。 结构面强度通常低于岩体强度,并 起控制作用。当前边坡数值计算中不考虑结构面强 度,有可能导致边坡不安全。边坡稳定计算主要是控制强度,因而网格划分时只需在某些可能 滑动的面上设置节理单元,而不必把所有结构面都模拟出来。图6是岩质边坡实例,算得的安 全系数为1104,采用极限平衡解析法,其安全系数为1108,误差为318◊ 。 413 有限元法强度折减法 在有限元法中,通过强度折减,直至计算到不收敛为止,其折减的倍数即为稳定系数。 这种 方法在国外八十年代就采用,但是由于力学概念上的不十分明确,而且要受到计算程度及计算 精度的影响,因而这种方法至今没有在国内流行。 随着各类通用大型软件的出现和计算水平的提高,我们认为这种方法已有条件应用于边 坡稳定安全系数的计算。 通过实例计算,我们发现只要严格控制计算精度,并采用莫尔2库仑准 则,其计算结果就与传统方法十分接近。 由于采用莫尔2库仑条件求解十分困难,许多大型软件 如AN S YS、M A RC、PA TRAN、N A S TRAN等只给出了广义米赛斯准则 FΑI1J2k26 式中I1,J2分别为应力张量的第一不变量和应力偏张量的第二不变量。 这是一个通用表达式,通过变换 Α 、k的表达式就可以在有限元中实现不同的屈服准则。Α, k是与岩土材料内摩擦角 Υ和内聚力c有关的常数。对于不同的圆图7,有不同的 Α,k。 1当 Α,k满足下列表达式时, Α 2sinΑ 3 3 -sinΥ k 2ccosΥ 33 -sinΥ 27 屈服面在 Π平面上为不等角度的六边形的外接圆。 2当 Α 、k满足下列表达式时, Α sinΥ 33 sin2Υ k 3ccosΥ 3 sin2Υ 28 屈服面在 Π平面上为不等角度的六边形的内切圆,在国内特指此圆为D2P准则。 3当 Α,k满足下列表达式时, Α 3sinΥ 3 3 cosΗΡ-sinΗΡsinΥ , k 3 cosΥ 3 cosΗΡ-sinΗΡsinΥ 29 ΗΡ sin- 1 - 2 3A sinΥ 4 9A 2sin2Υ- 4 sin2Υ 3 A 2 3 -1 1 2 2 sin2Υ 3 1 9622001年第4期 郑颖人等边坡稳定分析的一些进展 1994-2007 China Academic Journal Electronic Publishing House. All rights reserved. A Π 9 - sin2Υ 63 这是徐干成、 郑颖人1990提出的莫尔2库仑等面积圆屈服准则[7],它的面积等于不等角 六边形莫尔2库仑屈服准则的图形面积,它比当前采用的逼近不等角六边形的近似屈服曲线有 更高的计算精度。 稳定系数用 Ξ表示,它满足稳定系数的定义,此时屈服准则写成如下形式 F Α Ξ 1J2 k Ξ 30 计算到不收敛时的 Ξ,即为边坡稳定安全系数。令F1Α1I1 J2k1,为实际采用的准 则,F2Α2I1J2k2。 为等面积圆莫尔2库仑准则,获得二准则之间的比值 Γ F1 F2 Α1 Α2 k1 k2 fΥ , 由此即能将按实际采用准则求得的安全系数转换成按莫尔2库仑准则的系数。 414 算例比较 均质边坡,坡高H 20m,土容重 Χ 25kNm 3, 粘结力c 42kPa,内摩擦角 Υ 17,求坡 角 Β 30、35 、40、45、50 时边坡的安全系数。计算结果见表3。 各种方法计算得的安全系数表3 坡角度 方法 3035404550 有限元法外接圆屈服准则11781162114811361129 有限元法莫尔2库仑等面积圆屈服准则11471134112211121106 简化B ishop法1139411259111531106201992 S pencer法11463113181121211115110384 从表中计算结果可以看出,采用外接圆屈服准则计算的安全系数比传统的方法大许多,采 用莫尔2库仑等面积圆屈服准则计算的结果与传统方法计算的结果十分接近,说明采用莫尔2 库仑等面积圆屈服准则来代替莫尔2库尔不等角六边形屈服准则是可行的,这样使计算大为方 便。而采用外接圆屈服准则计算的安全系数要比莫尔2库仑等面积圆屈服准则计算的结果大 Γ 倍。 415 有限元强度折减法的优点 如果使有限元法保持足够的计算精度,那么有限元法较传统的方法具有如下优点 1能够对具有复杂地貌、 地质的边坡进行计算; 2考虑了土体的本构关系,以及变形对应力的影响; 3能够模拟土坡的边坡过程及其滑移面形状如图8211,滑移面大致在水平位移突变 的地方及塑性变形发展严重的部位; 4能够模拟土体与支护的共同作用,图12为无锚杆时边坡稳定安全系数为111,图13 为有锚杆支护时安全系数为115,且塑性区后移。 5求解安全系数时,可以不需要假定滑移面的形状,也无需进行条分。 072地 下 空 间 第21卷 1994-2007 China Academic Journal Electronic Publishing House. All rights reserved. 图8 变形前的网格图 图9 变形后的网格图 图10 水平方向位移等值云图 图11 变形后的塑性区 图12 不加锚杆 Ξ θ 1. 1 时的塑性区 图13 加锚杆 Ξ θ 1. 5 时的塑性区 参考文献 [1] 时卫民,叶晓明,郑颖人,王敬林 1 阶梯形边坡稳定性分析 1 三峡库区城镇移民建设区高边坡支挡结构 的合理选型及其综合利用课题中间报告[R ]12001171 [2] 时卫民,郑颖人 1 滑移面为直线假设下的斜坡稳定分析[J ]1 公路交通技术, 2001,No. 3. [3] 林丽 1 路堤边坡的极限平衡法稳定性分析[D ].重庆后勤工程学院硕士论文 1 [4] 赵尚毅,郑颖人,等 1 用有限单元法求边坡稳定安全系数 1 三峡库区城镇移民建设区高边坡支挡结构的 合理选型及其综合利用课题中间报告[R ]12001171 [5] D. V. Griffiths and P. A. lane.Slope stability analysis by finite elements[J ].Geotechnique, 493 387- 403, 1999 [6] 郑颖人,龚晓南 1 岩土塑性力学基础[M ]1 北京中国建筑工业出版社, 19891 [7] 徐干成,郑颖人 1 岩土工程中屈服准则应用的研究[J ].岩土工程学报, 1990,第2期 1 1722001年第4期 郑颖人等边坡稳定分析的一些进展 1994-2007 China Academic Journal Electronic Publishing House. All rights reserved. UNDERGROUND SPACE Vol . 21 No. 4SerialNo. 87Q uarterly Dec. 2001 CON TEN TS AND ABSTRACTS Papers Investigation on M etro Construction in RussiaQ I AN Q ihu241⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯ [Abstract]Based on investigation into Russian metro and acquaintancew ith the metro in large cities of other countries, the change and development ofmetro design and advance in construction in Russiawere intro2 duced.Introduced in particular is the principle of metro lines selection, the distance between metro stations in large cities of different countries, design of metro lines, protection from noise and vibration and design of metro buildings and structures .The versions and classification of metro stationswere analyzed aswell . [Keywords]Russia; metro design; transition station; separated running; parallel running D igital Panoram icBorehole Camera Technique and D igitalBoreholeGE Xiu2run et al254⋯ [Abstract]D igital Panoram icBorehole Camera TechniqueDPBCTrepresents a direction of development for borehole camera technology.Itsmain characteristic ispanoram ic technique and digital technique . The digital in2 ation gained in borehole is the basis of setting D igitalBoreholeDB . In this paper, the ation, basic theory, key technique and application of the DPBCT are discussed systematically.The basic conception of the DB. and the relationship between DPBCT and DB is analyzed. [Key words]panoram ic borehole camera, digital borehole, expanded i mage, virtual coregraph Progress in A nalysis of Slope StabilityZHEN G Yingren et al . 262⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯ [Abstract]In thispaper the progress in analysison slope stability by li m it equilibrium and finite elementmeth2 odswas introduced. Based on li m it equilibrium, an equation for calculating the safety factorsof stability ofmulti2 step slope and the angle of sliding surface has been derived, which is i mportant for design of rock slope. A united equation for present various slice s has been derived. For non2strict slice the safety factor can be obtained by one equilibrium equation given the orientation of force action between the slices .For strict slice the safety factor can be obtained by two equilibrium equations given the orientation or the point of force action between the slices . The united equation is a si mple iterative one and the calculation can be pered si m2 ply w ith high accuracy. Twos for obtaining safety factorsof slope stability by FEA have been presented. 1994-2007 China Academic Journal Electronic Publishing House. All rights reserved. One ofwhich is based on the li m it equilibrium w ith arch search for soil slope and arrangement of joint el2 ement on sliding surface for rock slope. The another one is a strength reduction of finite elements for application of large software, which is a new for safety factor of slope stability. [Keywords]slope; stability; li m it equilibrium; FEA; safety factor for slope stability A nalyses on Strength ofM arble underU niaxialCompression and the Influences of Strain Rates XU Songlin et al . 272 ⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯ [Abstract]U sually the shear fractureoccurs in the rupture surface near the top and the bottom ofmarble spec2 i mens under uniaxial compression. Henceforth, the strength of speci mens could be determ ined on account of its complex distribution of stress .The result consistswellw ith the strength of triaxial compression. W ith the in2 crease of strain rates, the C value increased whileΥdecreased. [Keywords]marble; uniaxial compression; strength; influence of strain rates A nalysis on Seepage through Dam Foundation of Fractured Rock M assCHA IJunrui276⋯ [Abstract]It is very i mportant to investigate the seepage through dam foundation of fractured rock mass, which is one of the bases for design of dam construction.The principle of seepage through fracture network for discontinuousmedia is utilized to analyze the seepage through dam foundation w ith curtain grouting and drain holes in different combination.It can be illustrated from analysis results that the curtain grouting controls seep2 age discharge and the drain holes have obvious effect on drainage and pressure decreasing. [Keywords]fractured rock mass; dam foundation; seepage; fracture network Experi mentalA nalysis andM athematicalSi mulation on the Strength CharacteristicsofL i mestone GUO Zhonghua et al . 280 ⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯ [Abstract]The strength characteristicsof li mestone are studied in this paper . The non2linear relationship be2 tween uniaxial compressive strength and the ratio of highness to diameter, and the non2linear relationship be2 tween uniaxial compressive strength and cross section in natural and saturation condition are presented separate2 ly. The relationship ulas are very convenient to provide consultation and reference for uating the strength of the resemble rock in the related engineering. The rock under in2situ conditions is al ways in the state of three2di mensional compression, so the laboratory triaxial compression test is the relatively true reflection of the stress2state of the rock. A i m ing at the resultsof triaxialcompression test, the curvesof stress strain relation2 ship of li mestone are analyzed and compared and then fitted aswell .The results indicate that the procedure of li mestone deation can be divided into four different stages .The stress2strain relationship is affected a great deal by the lateral stress and the stress2stra