用有限元强度折减法求边坡稳定安全系数

合集下载

用有限元强度折减法计算边坡安全系数

用有限元强度折减法计算边坡安全系数
有着广泛 的适用性 和 良好的应用前景 。 用的屈服准则主要有摩尔一 库仑准则 ( 简称 M C准则 ) — 和德鲁克一

致的 。基本公式为 : 折减后 的抗剪强度指标分别为 :
c =
音t =∞。 ,p a m
2 路面翻浆治理 。路 面翻浆一 般发生在 春天融 消时期 , ) 对大 水 , 防止水分蒸 发 返盐 , 行 土壤保 湿 ; 物成 活后 进 行抚 育 管 进 植 面积路面翻浆 , 采用更换基层 、 垫层粘性 盐渍土 的办法 , 因地 制宜 理 ; 利用植物每年 自然返青 的规律达到永久性边坡防护的 目的。 填筑明砂或砾石 材料 。小 面积 路面 翻浆 , 未形 成翻 浆之 前 , 在 采 参考文献 : 用在路面打砂桩 的办 法 。这种 办法 简单 实用 , 本低 , 理 速度 [ ] 周 华 , 成 治 1 李全胜 , 朱瑞 成. 疆库 尔勒地 区盐 溃土的 工程 地 新 快, 不影 响正常交通 , 治理 效果 好。 刷沥青等防腐材 料。对产 生孔 隙 、 裂缝 、 落部 位及 时采 用树 脂 剥 质特征 [] 西部探矿 工程 ,o o 3 :3 - 3 J. 2 o ( ) 121 . 3 探 讨[ ] 公路 工程地 质 ,90 8 4 :07 . J. 19 , ( ) 7 -2 3 桥涵防腐养护维修措施 。对 桥涵易腐 蚀部位清 除集盐 , ) 涂 [ ] 黄立度 , 2 席元伟. 中盐分和积聚迁移规 律 以及 盐胀机 制的 土
第3 8卷 第 4期 20 12年 2月
山 西 建 筑
SHANXI ARCHnE CTURE
Vo . 8 No 4 13 .
F b 2 1 e . 02
・9 ・ 7
文章编号 :0 9 6 2 (0 2 0 —0 70 10 — 8 5 2 1 )4 0 9 —2

强度折减法计算安全系数实例

强度折减法计算安全系数实例

强度折减法计算安全系数实例【原创版】目录一、引言二、强度折减法计算安全系数的原理三、实例分析四、结果对比与讨论五、结论正文一、引言在工程领域,边坡稳定性分析是一项重要的研究内容。

为了确保边坡的安全稳定,工程技术人员需要对其进行安全系数计算。

安全系数是指边坡在实际工况下的承载能力与实际荷载之间的比值,该值越大,边坡越稳定。

目前,常用的计算方法包括有限元强度折减法、极限平衡法等。

本文将以强度折减法为例,介绍计算安全系数的实例。

二、强度折减法计算安全系数的原理强度折减法是一种常用的边坡稳定性分析方法,其核心思想是按照一定的折减系数对边坡岩土体的强度进行折减,以考虑工程荷载作用下边坡稳定性的影响。

具体来说,首先需要建立边坡岩土体的有限元模型,然后对模型进行强度折减,最后计算出边坡的安全系数。

三、实例分析假设有一个边坡工程,边坡高度为 100 米,边坡底部宽度为 100 米,边坡顶部宽度为 50 米。

为了计算边坡的安全系数,首先需要建立有限元模型,包括以下几个部分:1.建立边坡岩土体的几何模型;2.划分有限元网格;3.设定材料参数,包括弹性模量、泊松比、密度等;4.应用强度折减法,对模型进行强度折减;5.计算边坡的安全系数。

在计算过程中,需要选用合适的材料模型和参数,以保证计算结果的准确性。

同时,需要注意考虑边坡的实际情况,如边坡的倾斜角度、边坡底部的支撑条件等。

四、结果对比与讨论通过上述实例计算,可以得到边坡的安全系数。

为了验证计算结果的准确性,可以将其与极限平衡法等其他方法进行对比。

在实际工程中,由于地质条件、工程荷载等因素的复杂性,不同方法计算出的安全系数可能存在一定的差异。

整体来说,有限元数值方法的计算结果会更加准确。

然而,在有限元方法中,由于模型的建立、材料参数的选择等因素的影响,计算结果可能存在一定程度的误差。

为了提高计算精度,可以采用多种方法,如选用合适的材料模型、合理设定材料参数、考虑边坡的实际情况等。

求解安全系数的有限元法

求解安全系数的有限元法

求解安全系数的有限元法
在边坡稳定性分析中,有限元法(Finite Element Method, FEM)被广泛用于求解土坡的安全系数。

安全系数是衡量边坡稳定性的指标,它代表了边坡实际的抗滑力与潜在滑动力之间的比值。

传统的极限平衡法通过确定可能的滑动面并计算作用于该面上的剪切强度和力矩平衡来估算安全系数。

然而,在有限元框架下,求解安全系数通常采用以下两种方法:
1. **有限元强度折减法 (Finite Element Strength Reduction Method, FSRM)**:
- 此方法基于逐步减少土体材料的抗剪强度参数(如内摩擦角或粘聚力),模拟土体逐渐趋向破坏的过程。

- 在每个折减步长上,重新进行有限元分析以获得新的位移场和应力状态。

- 当土体出现明显的塑性流动或达到预设的位移增量时,停止折减过程,并根据最后一次非线性迭代的结果计算出相应的安全系数。

- 这种方法得到的安全系数往往偏高,因为它考虑了整个土体的非线性响应,而非仅限于单一滑动面。

2. **结点位移法**:
- 结点位移法也是强度折减法的一种形式,通过监测随着抗剪强度降低,某些关键节点(如可能的滑裂带上的节点)的位移变化情况。

- 当位移突然增大时,表示潜在的滑动面已接近失稳状态,此时的抗剪强度折减比例可以用来反推安全系数。

有限元迭代解法也可以应用于边坡稳定分析中的复杂问题,例如当滑动面不明确或者滑动模式非常复杂时。

这种方法要求更为精细的网格划分和更为严谨的收敛条件控制,确保计算结果的准确性和可靠性。

有限差分强度折减法求解边坡安全系数

有限差分强度折减法求解边坡安全系数

显式有限差分方法在数值计算 中不需要像有
限单元法那 样 将 每个 单 元 组装 成 庞大 的总 刚度 矩 阵, 其突 出特 点是计算 速度 快 , 储量 小 , 存 不存在 计 算收敛性 问题 . 目前 , 在连 续 性介 质 力 学 问题 的范 畴 内, 显式有 限差分法 已经被 证 明具有 与有 限单 元
计算 处理从 而 获得边 坡 稳定 性 安 全 系数 和有 关 变 量 的等值 阴影 图.
4 边 坡 的 失稳 判 据
目前 , 断边 坡 是 否 失 稳 还 没有 一 个 统 一 的 判 标准 , 即安 全 系数求解 过 程 的终 止 条 件不 统 一 . 我
F 为折减系数 ) 然后对边坡稳定性 系数进行数 , 值分 析 , 不断 地 增 大 折 减 系 数 , 复 计 算 , 至达 反 直
影 响较 大 , 判据 c 法 单 独地 判 断边 坡是 否 失稳 , 无 因此 在下 面 的算 例 分 析 中 , 联 合 采 用 塑 性 区 从 将
坡 脚 贯通 至坡顶 和速度 矢 量 图来 判 断 边坡 是 否 失
稳. 因为 当分 析 的结果 中 出现贯 通 的 塑性 区时 , 仅
仅 能 表示该 边 坡 的某 些 区域 已 出现 了塑 性 破 坏 ,
关键词 : 边坡安全系数 ; 稳定性 ; 强度折减法 ; 三维快速拉格朗 日法 ; 数值模拟
中 图分 类号 :U 5 T 47 文 献 标识 码 : A d i1 .9 9 ji n 17 -8 9 2 1. 40 5 o:0 36 /.s . 642 6 .0 20 .0 s
0 引 言
通过 工程 实际经 验发 现任 何 边坡 的失 稳破 坏 都 伴 随着 一定 的位 移 变化 , 边 坡 内 的某 特 征 点 位 移 当 发生 突变 时 , 们 就 可 以认 为 该 边 坡 处 于 临 界 失 我 稳状 态 . 是 , 但 该准则 受 特征 点 所选 位 置 以及 位 移 容许 限值 的影 响 比较 大 , 以运 用起来 不 方便 . 所 c 以塑 性 区从 坡脚 至坡 顶贯 通 为标 准 . .

有限元强度折减法研究岩质边坡稳定性

有限元强度折减法研究岩质边坡稳定性
袭 l 糟瑾力学’披计算取值
√ 一一应力张量的第一不变量;
√ 一一 力 张 第 不 量。 应 偏 量的 二 变
这是个通用表达式, 通过变换 、k的
表 达式 就可 以在有 限 元中实现 不 同的屈服 准 则, 图l 如
, =

.+ tql ( f rn) 4 ( y c a ̄ d )
给 定的强度折减参数 C 口 , 土体形成的的 广义剪应 变自 角底 部下方 向坡 顶贯通 , 坡 则认
将 等式 两边 同。 时积分 :
为对 于稍 高于 c 和 的土体强度参数 c 和 . 使得土体处于临界状 态 , 而与 q和 矗 对应的强 。 为安全系数。 但莫 尔 库仑准 则在 平面上 的图形 为不规 度折减 系数 l 则的六 角形截 面, 在尖顶 和棱角 , 存 给数值 计 算 带来 困难 。 广义 米赛 斯 准 则在c 力空 但 式 主应 间的屈服 面为一 圆锥面 , 7 平面上为圆形, 在 t " = 不 存在 尖顶 处 的数 值计 算 问题 , 目前 国际上 流行的ANS 以 及美国MS YS C公司的MARC, N S RA A T N等均采用广义米赛 斯准则。 2 屈服准 则的选用 本文计算采用的是理想弹 塑性 模型。 用 使 有限元软件ANS 进行分析, YS 采用Dr c e uk r P a e准 则。 rg r () 3 在常用的 极限平 衡方法 中, 安全 系统 定义 F cl 4 2 k =d + J = () 1 为沿滑动面抗 剪强度与滑动面 上实 际剪力的比 式 中: 值。
— — — — — — 一
I 2 00 l 5
】2 5 O 01 , 8 】 1 4
: ● n 】 j c I 3
l 4 ¨¨ J ?

有限元强度折减法在边坡稳定性分析中的应用

有限元强度折减法在边坡稳定性分析中的应用

55 6 50 6 55 5
5o 5
55 4 50 4 55 3
50 3
55 2 50 2
0 1 2 3 4 5 6 7】 8 9 1 0 1 2 0 0 0 0 0 0 【 0 0 0 1 0 1 0
阕 1 滑坡稳定性计算剖面
. 陶系 中 、 统 (2+3 龟 裂 纹 状 灰 岩 。 上 部 以 青 灰 色 厚 层 龟 裂 纹 2 3 计 算参数 的选取 上 0 ) 滑动面抗 剪强度参 数 的选取 关系 到滑坡稳 定性 预测 和滑坡 状灰岩为主 , 该层 泥 质 含量 质 抗 岩
水, 大气降水经溶 沟 ( )溶 蚀裂 隙下 渗 , 洞 、 储存 于溶 沟 ( ) 溶蚀 参数包 括峰值 强 度 和残 余 强度 两 种 。其峰 值 强度 平 均值 = 洞 、 1 .。C =4 . P , 3 9,。 3 2k a残余强度平 均值 =1 . 。C =3 . P 。 2 3 ,, 9 0k a 裂隙中 , 水量分 布受 岩溶 发育程 度及 气候 的影 响较 大 , 由于该处 并 选 龟 裂纹 状 灰 岩 岩 溶 发 育 一 般 , 多 以 浅 层 发 育 为 主 , 下 水 位 较 结合滑动带残余剪 切试 验指 标 , 结合 反算 的结果 , 取滑 动带 大 地 计 算 参 数 为 =1。C 0k a 6 , =1 P 。 高 , 调 查 , 段 边 坡 上 下 水 位 离地 面 一 般 为 1 . 据 该 00m~1 . 5 0m。
有 限元 强度 折减 法 在 边 坡 稳 定 性 分 析 中 的应 用
周 攀峰

韩 利 光
要: 应用弹塑性 强度折减 有限元计算边坡稳定安全 系数 , 通过算例对极 限平衡法和 强度折减 法计算 的安全 系数进行

基于ANSYS的有限元强度折减法确定边坡安全系数

基于ANSYS的有限元强度折减法确定边坡安全系数

基于 ANSYS的有限元强度折减法确定边坡安全系数摘要:本文基于ANSYS,采用D-P外角点外接圆屈服准则对国内某矿区边坡进行稳定性计算分析,通过不断对边坡强度参数黏聚力和内摩擦角进行折减,直到软件计算不收敛为止,其折减的倍数即为边坡稳定安全系数。

计算结果显示,利用ANSYS自带D-P本构模型计算得到的边坡安全系数远大于极限平衡法计算得到的边坡安全系数。

最后应用不同屈服准则安全系数的转换关系得到该边坡平面应变下与M-C匹配的D-P准则的安全系数,并与极限平衡法结果对比,吻合较好。

据此得出结论:在估算边坡安全系数方面,采用有限元强度折减法是一种值得信赖的方法,但计算中采用理想弹塑性材料模型时,屈服准则的选择会对边坡安全系数的计算产生较大影响。

关键词:有限元;强度折减;屈服准则;边坡稳定;安全系数1 引言目前,边坡稳定性分析发方法较多,主要有定性分析法(图解法、类比法)、定量分析法(极限平衡法、数值分析法)、非确定分析法(模糊分析评判法、可靠性分析法)。

而对于边坡安全系数,许多学者大多用定量分析法[1-2]。

传统的极限平衡法首先要确定一个潜在的滑动面,基于一系列简化假定后,由力系平衡或能量守恒求得滑动面的安全系数,用它作为评价边坡安全性的指标。

这些方法有瑞典条分法、简化毕肖普法、简布法、不平衡系数传递法等。

这些方法的基本出发点是一样的,即刚塑性假定,不同之处在于对条间力所作的假定不同。

由于这些假定的物理意义不一样,因此它们所能满足的平衡条件也不相同,计算步骤有繁有简,为了检验所列的各方法和其他边坡稳定性分析方法的精确性,许多学者在过去几十年里从不同角度做了大量研究并进行了系统总结[3-7]。

传统的极限平衡法由于没有考虑土体内部应力与应变的关系,故无法模拟分析土体发生变形甚至破坏的过程。

随着计算机技术的发展,数值计算方法在边坡稳定分析中得到了广泛的应用。

其最大的优点是求解安全系数时,不需要假定滑移面的形状和位置,也无需进行条分,可以分析任何形状的几何体,不但能进行线性分析还可进行非线性分析。

边坡稳定分析中的有限元强度折减法研究

边坡稳定分析中的有限元强度折减法研究

边坡稳定分析中的有限元强度折减法研究【摘要】在边坡稳定性分析中采用有限元强度折减法,主要包括破坏判据的确定,屈服准则的影响和选用。

本文详细介绍有关有限元强度折减法的内容。

【关键词】有限元强度折减法;边坡稳定分析一、有限元强度折减安全系数对于摩尔一库仑材料,强度折减安全系数定义与边坡稳定分利用有限单元法,通过强度折减来求边坡稳定安全系数。

通过强度折减,使系统达到不稳定状态时,有限元计算将不收敛,此时的折减系数就是安全系数。

计算时,首先选取初始折减系数,折减土体强度参,将折减后的参数作为输入,进行有限元计算,若程序收敛,则土体仍处于稳定状态,然后再增加折减系数,直到不收敛为止,此时的折减系数即为边坡的稳定安全系数,此时的滑移面即为实际滑移面,这种方法称为土体强度折减系数法。

对于摩尔一库仑材料,强度折减安全系数定义与边坡稳定分析的极限平衡条分法安全系数定义形式是一致的。

有限元强度折减系数法的基本原理是将坡体强度参数:粘聚力c和内摩擦角值同时除以一个折减系数fs,得到一组新的c′、′值,然后作为新的资料参数输入,再进行试算,当计算不收敛时,对应的fs被称为坡体的最小稳定安全系数,此时坡体达到极限状态,发生剪切破坏,同时可得到坡体的破坏滑动面。

二、有限元中边坡破坏的判据有限元强度折减法分析边坡稳定性的一个关键问题是如何根据有限元计算结果来判别边坡是否达到极限破坏状态。

土体破坏的标准有如下几种:1.以有限元静力平衡计算不收敛作为边坡整体失稳的标志。

2.以塑性区(或者等效塑性应变)从坡脚到坡顶贯通作为边坡整体失稳的标志。

3.土体破坏标志应当是滑动土体无限移动,此时土体滑移面上应变和位移发生突变且无限发展。

三、屈服准则的影响和选用研究表明,采用该准则与传统摩尔一库仑屈服准则的计算结果有较大误差,不管是评价边坡稳定性,还是地基极限承载力等等,在实际工程中如果采用该准则是偏于不安全的。

四、有限元法进行边坡稳定分析的优点如果使有限元法保持足够的计算精度,那么有限元法较传统的方法具有如下优点:1.能够对具有复杂地貌、地质的边坡进行计算。

基于非线性动力的有限元强度折减法边坡稳定性分析

基于非线性动力的有限元强度折减法边坡稳定性分析

的阻尼 比相 同 ,根 据 《 工建 筑 物 抗 震 设 计 规 范》 水 取边坡介 质材 料 的 阻尼 比为 0 0 。为 了得 到地震 作 .5
用下边坡有 限元模 型微分方程 的阻尼矩 阵 ,对边坡进 行模态分析 ,取前 两阶振 型对 比例 因子计算 。
1 2 有 限元 强度折 减法基本原理 .
{/ },动力荷载 可表示 为 p =一[ 】 五 () /( ) () { £); { () , 五 ) { () 分别 为相 对加 速 度 向量 、速 五 ) { () , M f)
地震诱 发的边坡失稳是主要地质灾 害之一 ,如 中 国海 原大地 震 ( . ,12 ) 中形成 的大规模 黄土 M 85 9 0
的影 响 ,采用有 限元 强度 折减法对边坡稳定性 进行分 析和评价 ,并提 出了边坡 的处治措施。 1 基本理论与计算方法 1 1 非线性 动力响应分析 . 在地震 作用下 ,边坡的动力平衡方程用下式描述
【 { ) +[ 】 ( ) K】 u £ ) ( )( ) M】 ( ) C { ) +[ { ( ) :p £ 1
度向量 、位移 向量 。阻尼矩阵采 用 R y i al g e h阻尼 ,将
阻尼 矩阵简化为质量矩阵和刚度 矩阵的线性组合 ,计 算式 为
[ 】= a[ 】 [ 】 C +b () 2
滑坡群 ,滑坡面积达 3 0 m ;四川炉霍境 内地震 0k 8
( . ,17 )诱 发各 种 规 模 滑 坡 17处 ;美 国 M79 93 3 Lm rt o aPi a地 震 ( . ,1 8 ,触 发 的 滑 坡 达 e M 7 1 9 9) 1 0 处 ,范 围 达 1 0 m 0 3 50 0 k ;中 国 汶 川 大 地 震

基于有限元强度折减法确定某边坡的安全系数

基于有限元强度折减法确定某边坡的安全系数
边坡构成都 比较复杂 , 大多含有 多种介质 , 尤其是岩 质边坡 , 际 安全系数 甜可表示为 : 实 中含有大量 的结构面 。 边坡稳定分析涉及复杂 的地质地形边 界条件 、 材料的应力一 应变 的非线性行 为、 初始地应力 、 水压力 、 地震荷载的耦合分析等
r C t /J tn =( +盯a ) c=C + a n £ () 1
第3 6卷 第 1 0期 2 0 10 年 4 月
山 西 建 筑
S HANXI ARCHI TECTURE
Vl. 6 No. 0 D 3 】 1
A r 2 1 p . 00
・11 ・ 1
文章编号 :0 96 2 (00)00 1 -2 10 —8 52 1 1 —1 10
5 结论和建 议

( )2 . 9 :3
[ ] 陆基孟 . 4 地震勘探 原理[ . M]北京 : 油大学出版社 ,9 3 石 19 .
当前 , 煤炭三维地震勘探 的重点 已转 入中部和西部 的富煤省 份, 这些地 区地表施 工条件 相对 复杂 , 山地 三维地震 就是其 中之

[] 张爱敏 . 区高分辨率 三维地震 勘探 [ . 州 : 5 采 M]徐 中国矿 业
基 于 有 限 元 强度 折减 法确 定非线性有 限元分析软件 MA C, R 采用 强度折减法确定某边坡 的安全系数 , 将塑性 区贯 通作 为边坡破坏 的判
据, 分析结果表 明, 此方法可以有效 克服传统 的缺 陷, 有较强的实用价值 。
关键词 : 边坡 , 强度折减 法, 安全系数
中图分类号 : 143 6 TJ 1 .2 文献标识码 : A
0 引言
下水 、 施工过程对边坡 稳定性 的影 响 , 以考虑 各种支 挡结构 与 可

用有限元强度折减法求算边坡安全系数初探

用有限元强度折减法求算边坡安全系数初探
s a f e t y c o e f i c i e n t o f s l o e p i s c lc a u l a t e d b y he t s t r e n g t h r e d u c t i o n i f n i t e e l e m e n t m e t h o d ② F o r n o n - h o m o g e n e o u s s l o e, p h t e v lu a e s o f C, h a v e m o s t i n l f u e n c e o n t h e l i n k i n g a r e a o f e q u i v a l e n t p l a s t i c a r e a ③ T h e d i s p l a c e m e n t v e c t o r o n he t s l o e p f a c e s h o u l d r e l f e c t he t c u r v e d t a n g e n t d i r e c t i o n o f t h e c e r t i a n p o i n t l o c a t e d o n t h e s l i d i n g f a c e .④ Wh e t h e r t h e s a t -
熊传 治 , 杜 维吾 ,彭 状 ,孙晓然 ,王新 建
( 长沙矿冶研究 院, 湖南 长沙 4 1 0 0 1 2 )

要: 对非均质边坡岩体失稳用有 限元强度折 减法开展 了初步分析研究 , 对 已有文献提出的 3条失稳判 断之 一进行 了修正 , 即坡
பைடு நூலகம்
面位移矢量方 向应反映相应滑面点 的切线方 向。结合可 可托 海露天矿的强度折减 法实例 , 得到 的 4 点初 步结论是 : ① 用有 限元强

有限元强度折减法在边坡稳定计算中的应用

有限元强度折减法在边坡稳定计算中的应用
时考虑到该点 的抗剪强度 ,可用 M0 C u0 b破坏准则表 示 边坡稳定 的影响。 h一 0 l m 合理的做法是采用既不采用关联流动法则 (, I I = 为: T eoa‘ t +-n = t D 则该点土体在这个预定剪切 面上的安全系数为:
T c f +叮 t n P a ‘
维普资讯
民 蕾科 技
丽 而
建筑・ 卜 凇 设计
有限元强度折减法在边坡稳定计算中的应用
王 中发 孙荣鸿
( 湖北水利水电职业技术学院水利建筑工程 系, 湖北 武汉 40 7 ) 3 00
摘 要: 对于边坡稳 定性分析 , 有限元强度折减 系数法正逐渐成为求解边坡稳 定分析研究的新趋势, 其基本原理是将岩土材料
对于边坡稳定性 的分析 , Fl nu 于 12 自 el is 9 7年提出圆弧滑 的假定 : e 其一 为破坏准则和屈服准则 ; 其二为流动法则 。对 于这
动法 以来 , 至今 已先后 出现极限平衡法 、 限分析法 、 极 有限元法 方面假定的不同可以构成不同 的土体本构模型。 等数十种分析方法。目前在分析边坡稳定性时 , 最常采用极 限平 衡条分法和弹塑性有限元数值分析两种方法。 21 . 破坏准则 和屈服准则 破坏准则f 判断土体加卸载过程 中破坏与否的标准 。 , 是 屈服

() 4
为了验 证强度折减法能够应用于 边坡的稳定性计算, 这里
析方法 的研究 ( 作者 : 关立军 ; 导师迟世春 ;0 3 ) 2 o 年 中的一个 算 例 ,边 坡土体参 数如 下 :土 体重度  ̄ 2k /’ / 0 Nm ,黏 聚力 c = = 4 k a 内摩擦角 q 2  ̄弹性模量 E 2 M a 泊松 比 v 03 边坡 0p, o 0, = =0 p , = .,

利用有限元强度折减法分析岩质边坡的稳定性

利用有限元强度折减法分析岩质边坡的稳定性

2 数值 模拟
2 1 算例 .
某岩质路基边坡坡高 为 3 坡 角为 4。坡脚到左端边界 的 0 m, 5,
1 强度 折减 法原理
1 1 安全 系数 的定 义 .
. 坡顶到右端边界的距离为坡高的 2 5倍 , . 且 D na 对边坡岩土体 的剪切强度进行 折减 , u cn 重新定义 了边 距离为坡高的 15倍 , 总高为 2 坡高 。有限元模型的边界条件 是底面为 固定约束 , 倍 J 坡的安全系数 , 即假 设岩 土体重 力加速 度恒 定不 变, 通过 局部 降 l 8 n 低 岩土体 抗剪强度 , C 值 同时除 以折减 系数 , 到新 的一 组 左右边界条件为水平约束 。岩土体单元 采用 Pae2号实体单元 将 , 得 ( 4节点四边形高 阶单元 ) 。流 动法则 采 用相关 联流动 法则 。岩 强度指标 c 进 行有限差分 , , 通过反复 计算 , 到边坡达到 临界 直 B5 2 89 4工程岩体分级标准取值 。 破 坏状态 , 此时得到的岩土体强度指标 与原有 的强度 指标 的比值 体物理 力学参数根据 G 0 1 - 即是边坡的安全系数。
其表达式如下 :
2 2 计 算成 果分析 .
利 用 A S S软件 , 立 二维岩 质边坡 模 型。岩 土体 的本 构 NY 建
c奇 a1 , , r1 = c t 8
以边坡的位移 计算不收敛及 塑性 区贯通作 为边坡失稳判据 , 得到边坡 的安全系数及破坏 滑动面。通 过与成熟 的极 限平 衡法做 比
较 , 明边坡稳定性安全 系数是合理 的, 而也说 明强度折减法在岩质边坡稳定性分析中的优越性 。 证 从 关键词 : 强度折减 , 岩质边坡 , N Y , A S S 安全系数

求解边坡稳定安全系数两种方法的比较

求解边坡稳定安全系数两种方法的比较

求解边坡稳定安全系数两种方法的比较摘要:目前,边坡稳定性分析主要有刚体极限平衡法和有限元强度折减法,本文就理论基础、安全系数的定义及优缺点对以上两种方法进行了简要评述。

基于极限平衡法的发展起来的各种方法物理意义简单,便于计算,但是需要许多假设。

有限元强度折减法不需要假设,可以直接搜索临界滑动面并求出相应的安全系数,同时考虑了岩土体的弹塑性和边坡的破坏失稳过程。

通过对两种方法的认识比较,给岩土边坡工作者设计施工提供一定的参考价值。

关键词:边坡稳定性;极限平衡法;有限元法;安全系数引言边坡稳定分析是一个非常复杂的问题,从20世纪50年代以来,许多专家学者致力于这一研究,因此边坡稳定分析的内容十分丰富。

总体上来说,边坡稳定分析方法可分为两大类:定性分析方法和定量分析方法。

定性分析方法主要是通过工程地质勘探,可以综合考虑影响边坡稳定性的多种因素,对边坡岩土体的性质及演化史、影响边坡稳定性的主要因素、可能的变形破坏方式及失稳的力学机制等进行分析,从而给出边坡稳定性评价的定性说明和解释。

然而,人们更关心的是如何定量表示边坡的稳定性,即边坡稳定性分析的计算方法,定量方法将影响边坡稳定的各种因素都作为确定的量来考虑,通常以计算稳定安全系数为基础。

边坡稳定分析的定量方法有很多种,如条分法、数值分析方法、可靠度方法和模糊数学方法等[1-3]。

目前,边坡稳定分析方法中,人们较为熟知且广泛应用的有条分法和有限元方法。

条分法在边坡稳定分析中最早使用,因其力学模型概念清楚、简单实用,故广泛应用于实际工程中,已经逐渐成为边坡稳定分析的成熟方法。

随着计算机技术的发展,数值分析方法在工程领域应用越来越成熟,有限元方法考虑了土体的非线性应力-应变关系,同时弥补了条分法的不足,近年来有限元方法得到了极大的发展。

[4-6]刚体极限平衡法刚体极限平衡法是人们提出的最早的一类方法,是边坡分析的经典方法,只需要少许力学参数就能提供便于设计应用的稳定性指标即安全系数。

边坡稳定性的有限元强度折减分析法

边坡稳定性的有限元强度折减分析法
】 1 1 66 (00 1 - 09 0 0 - 84 21 )1 07 - 3 0
THE ANALYS S 0l S I F L0PE S TABⅡ rY BY r THE TRE S NGTH REDUCT ON EM I F
L U Je JANG J n to, JA a ・ n Z I i , I u — I NG Xi f g, HAO C n a o e a
2 工 程 实 例 分 析 2 1 工 程 概 况 .
拟建工程位于天津市塘沽 区 , 由于周 围场 地较宽 敞 , 采 用放坡开挖 , 边坡 的几何尺 寸如图 1 示 , 比为 12, 所 坡 : 其平
( rht tr n ier gIs tt Taj nvri Taj 0 0 2 C ia A c i cua E g ei tu e l n n n i e,ini U iesy,ini 3 0 7 , hn ) n t n
A s a tT i pp ri b sdo es eghrd ci E ,o iigata po c,nlz g bt c : hs ae ae nt t n t e ut nF M cmbnn c l r e ta a i r s h r o u j yn
拟出边坡的实 际滑移 面。因 此 , 限元 法近 年来 在工程 中 有 的应用 取得了快速 的发展 。 1 有限元强度折减法原理 强 度折 减法就是将土体 的抗剪强度指 标 C和 妒用 一个 折减系数 按 ( )式进行 折减 , 1 然后用折 减后 的抗剪 强度
I ———— 一 —— —
t e s b l y o o a d t r e d me so a lp . h e ut r o a e i i l e ih p meh d h t i t f w n e — i n in l o e T e r s l a e c mp r d w t smp i d B s o t o a i t h s s h i f t h c h c u a y o e a ay i. ti n i ae h tt e sr n t e u t n F o c e k t e a c r c f t n l ss I s i d c td t a h t gh r d ci EM n lss fr s p h e o a ay i o l e o sa i t s fa i l . t sb t rt o s e e a t a i ai n b fr k h ie o o o r e d — t bl y i e s e I i et c n i r h c u ls u t e oe wema e a c o c ft r h e — i i b e o d t t o w t

用有限元强度折减法分析谈山隧道进口高边坡稳定性

用有限元强度折减法分析谈山隧道进口高边坡稳定性
C uo b准则是一种适 用于摩擦 类材料 的剪 切强 度准 ol m 则 ,是岩土力学研 究 中应用最为广泛 的强度 准则 ,在
数 J—F 的关 系 曲线 。郑颖 人等 研究 表 明 :可将 s s— 曲线上 位移 突变 作 为边坡 破 坏 的依 据 ,即可 F 将 s— 曲线上 拐点位 置 的 F 作 为边坡 的整体 安 F 全系数 。 1 3 有 限元强度折减法 的特点 . 有 限元强度折减法具有如下特点 j 。 ( )可以考 虑多 种 不 同的岩 土 体非 线 形 弹塑性 1
文献标 志码 :A
文章编号 :10 82 (02 0 05 — 2 0 3— 85 2 1 )3— 1 1 0
常采用 M h —C uo b准则 作 为 计 算 基 础 ,其 表 达 o r ol m 式为
=c+o a  ̄ -n t p () 1
边坡稳定性分 析一直 是岩 土工程 界 的重要课 题 。
陈世 刚 :用有 限元 强度 折减 法分析谈 山隧道进 1高边坡稳定性 : 3
・5 1 1・
用 有 限元 强度 折减法 分析谈 山隧道进 口高边坡 稳定性
陈世 刚
( 中铁第 四勘察设计 院集 团有 限公 司 ,武汉 406 ) 30 3

要 :有 限元 强度折减法是近些年 来发展起来 的一种边坡稳定性评价方法 ,采用有限元强度折
的边坡 稳定 性分 析 方法 。本 文采 用 有 限元 强度 折 减
变小 ,边坡 的安全 系数不 断降低 ,必然存 在一个 临界 值f ,C ;当f> 0 >c ,边坡 处 于 稳 定状 态 , o 。 厂 ,C 。时 当f f ,c 。 ,边坡 将发 生失稳 破坏 。强 度折减 <  ̄ o ≤c 时 法就是基 于上述基 本原理 ,将边坡岩土体 的强度参数 c 厂值 同时除 以一个折减 系数 F ,得到一组 强度参 , 数c …

用有限元强度折减法求解坝坡抗滑稳定安全系数

用有限元强度折减法求解坝坡抗滑稳定安全系数

ta h alr ufc su d rdf rn P c tr nweeco e ocasc l to h ttefi es r e n e i ee t u a f D. r e o r ls dt lsia h d,a dte sft a trf m i i me n h aeyfco r D— o
和材料非线性等问题 。 关键词 : 尾矿坝 ; 有限元强度折减法 ; 力稳定 ;N Y 静 A SS 中圈分类号 :96T 4 13 X 3 ;U 4.5 文献标识 码 : A
S u y o n ik dn t bl y s f t a t ro a l p t d n a t i i g sa i t a ey f co fd m so e s i
全系数是对坝体稳定性评价的重要问题之一。坝体
0 引言
尾矿 坝 的安全稳 定 问题对 于 矿 山的安 全生 产 和
稳定性的计算方法有传统 的极限平衡 法、 数值 方法 和可 靠 性 分 析 法 。极 限 平 衡 法 中 主 要 有 Fl n s eei l u
用 有 限元 强度 折 减 法 求解 坝 坡 抗 滑 稳 定 安 全 系数
李 明 , 梁 力 , 王 伟 , 刘冬 霞
( 东北大学资源 与土木工 程学院 , 阳 10 1 ) 沈 18 9

要: 通过 对有 限元 强度折减法的理论研究 , 将此法应 用到尾 矿坝的静力 稳定计算 中。采用 有
限元强度折减法和极 限平衡法对某 尾矿坝 进行 了坝坡抗 滑稳定 安全 系数 的计 算 , 算结 果表 明 , 计
在应用有限元强度折减法时 , 采用不同的 D— P准则 计算 出的滑动面基本相 同且与传统方 法相 近 ,
其 中 D—P 3准则计算的安全系数更加接近瑞典 圆弧法 。因此应用 有限元 强度折 减法求 解尾矿 坝 静力稳 定安全系数是可行 的, 同时说明有限元法可 以较好 的模拟 实际工程 中的复杂 几何断 面形 状

基于ANSYS的强度折减法求解边坡安全系数

基于ANSYS的强度折减法求解边坡安全系数

基于ANSYS的强度折减法求解边坡安全系数摘要:介绍了一种基于有限元的强度折减法的边坡稳定性分析方法,讨论了该方法的基本原理、安全系数的物理意义、屈服准则和流动法则的选用及边坡破坏的判据等。

算例通过不断减小边坡强度参数黏聚力和内摩擦角,得到新的一组黏聚力和内摩擦角,再输入软件计算至不收敛,此时的折减系数就是边坡的安全系数。

计算结果显示,用ANSYS计算边坡的安全系数有一定的实用性和可靠性。

关键词:强度折减;边坡稳定;ANSYS中图分类号:O 319.56 文献标志码:A 文章编号:1674-0696(2011)基金项目:重庆市科学基金重点资助项目(035679);2022年高等学校博士学科点专项科研资助项目(20020183061)0 引言当前,研究边坡稳定性的传统方法主要有:瑞典条分法,简化的毕肖普法等。

这些方法主要是建立在极限平衡理论基础上的,在进行稳定性分析时没有考虑土体内部的应力应变关系,无法分析边坡破坏的发生和发展过程,在求安全系数时通常需要假定滑裂面形状为折线、圆弧、对数螺旋线等。

而有限单元法不但满足力的平衡条件,而且考虑了材料的应力应变关系。

求解安全系数时,不需要假定滑移面的形状和位置,也无需进行条分,这使得计算结果更加合理精确。

有限元法较传统的方法有如下优点:①能够对具有复杂地貌、地质的边坡进行计算;②考虑了土体的非线性弹塑性本构关系,以及变形对应力的影响;③能够模拟土坡的失稳过程及其滑移面形状。

滑移面大致在水平位移突变的地方及塑性变形发展严重的部位,成条带状;性变形发展严重的部位,呈条带状;④能够模拟土体与支护的共同作用。

1.强度折减法在 ANSYS 中的实现有限元法全面满足边坡静力许可、应变相容和应力、应变之间的本构关系,使得有限元分析边坡稳定成为近年来的发展趋势,因为有限元法不仅仅能计算出土体的内应力场分布,还可以了解边坡逐步破坏机理,跟踪边坡内塑性区的发展情况。

1.1 边坡失稳的判别依据边坡滑动面塑性区贯通是边坡发生破坏的必要条件,但不是充分条件,采用塑性区贯通来判别趋于保守。

用有限元强度折减法进行边坡稳定分析

用有限元强度折减法进行边坡稳定分析

坡角(度)
30
35
40
45
50
有限元法(外接圆屈 服准则)
1.78
1.62
1.48
1.36
1.29
有限元法(莫尔-库仑 等面积圆屈服准则)
1.47
1.34
1.22
1.12
1.06
简化Bishop法
1.394 1.259 1.153 1.062 0.99
Spencer法
1.463 1.318 1.212 1.115 1.04
• 随着计算机技术的发展,尤其是岩土材料的非线性弹塑性 有限元计算技术的发展,有限元强度折减法近来在国内外受到 关注,对于均质土坡已经得到了较好的结论,但尚未在工程中 实用,本文采用有限元强度折减法,对均质土坡进行了系统分 析,证实了其实用于工程的可行性,对节理岩质边坡得到了坡 体的危险滑动面和相应的稳定安全系数。该方法可以对贯通和 非贯通的节理岩质边坡进行稳定分析,同时可以考虑地下水、 施工过程对边坡稳定性的影响,可以考虑各种支挡结构与岩土 材料的共同作用,为节理岩质边坡稳定分析开辟了新的途径。
从表中计算结果可以看出,采用外接圆屈服准则计 算的安全系数比传统的方法大许多,而采用莫尔-库 仑等面积圆屈服准则计算的结果与传统极限平衡方 法(Spencer法)计算的结果十分接近,说明采用莫 尔-库仑等面积圆屈服准则来代替莫尔-库仑不等角 六边形屈服准则是可行的。
通过4组计算方案(改变内摩擦角、内聚力、坡角β、
表2 不同流动法则的影响
3.2 有限元法引入的误差
3.2.1 网格的疏密
有限元单元网格划分
表3 网格疏密对计算结果的影响
3.2.2 模型边界范围
表4 边界条件对折减系数的影

有限元强度折减理论求安全系数及存在的问题

有限元强度折减理论求安全系数及存在的问题

加工过 程 中 , 也将使 自己的管理水 平得 到大 幅度提高 。 当前 , 工 项 目成 本 管 理 工 作 没 有 得 到 很 好 地 开 施 展, 其症 结 在 于 对 成 本 核 算 工 作 的模 糊 认 识 和 缺 乏 重 视 。而加强项 目成本 核算 , 将是 建筑 企业 进入 成本 竞 争
1 有 限元法强 度折减 法边坡 稳定 安全 系数 的定义
统储存 , 是成本核算工作得 以高效实施 的保障, 也是企 业成本 战略实 施 的关键 环节 。在施 工项 目管理 机 构 中 , 应要求每位项 目管理人员都具备一专多能的素质 , 既是 工程质 量检查 、 度 监 督人 员 , 进 又是 成 本控 制 和 核 算 人 员。管理人员每天结束工作前应保证 1 h的内部作业时 间, 中成本核 算工作 就是重 要 的 内业 之 一 。通过 项 目 其 管理软 件的开 发和项 目局域 网络 的建 立 , 每位 管理 人 员 的核算结果将按既定核算体系由计算机汇总后 , 将加工 信息提交 项 目经 理 , 为其 制 定 成本 控 制 措 施 的依 据 。 作 项 目经 理及管 理人员在 进行 成本 核算 、 据 总 、 数 [ 整理 、
合论述 和分析 比较 , 并提 出一些建议 。 关键词 : 限元 强度折 减 法; 有 边坡 稳 定安全 系数 ; 临界 失稳 状 态判 据 ; 算精 度 计 中图分 类号 : TU74 文献标 识码 : 文章编 号 :O 4 5 1 (O 7O一O O一 O 1 B 1 0— 7 62O )2 23 4
有限元或者有限差分法, 考虑土体的非线性应力应变关 系, 求得边坡 内部每 一计算 点 的应 力应变 以及 变形 , 这样 能比较真实的反映应力应变情况, 接着再降低岩土体材 料 的抗剪强度参数 , 至边坡 达到 临界破坏状 态 , 直 从而得 到边坡的安全系数。这样的做法与实际的边坡失稳过程 较为吻合 , 即大部分边坡失稳都是 由于岩土体材料的抗 剪强度降低所致。这样不仅可以了解岩土体结构物随抗 剪强度恶化而呈现出的渐进失稳过程, 还可以得到极限 状态下边坡 的失效形 式L。 1 ] 鉴于 以上优点 , 边坡稳 定安 全系 数计算 的有 限元 强 度折减法在 国 内外 引起不 少学 者 的注意 , 相继 有 多篇 相 关 论文 发表 [。, 用方 法 的不 同主要 集 中在 边坡 失稳 2。 - 所 1 极限状态的确定 , 计算结果存在差异 , 主要是因为在计算 过程 中所采用 的屈 服准 则 、 动法 则等 的差异 。就对 以 流 上两方面的研究现状进行分析综述。
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

!
屈服条件的转换
采用不同的屈服条件得到的边坡稳定安全系数是
不同的, 但这些屈服条件可以互相转换。下面提出如 何将按实际所采用的屈服准则求得的安全系数转换成 莫尔 ! 库仑条件下的安全系数。 以外接圆屈服准则为例, 外接圆屈服准则表示为 或 & * " "$ - " ! %* # * $ % * , 其中%* %* # * $ "$ - " % * , -./’ ) " 12. # # 。 , %* " " (0 ! ./’ ) ( 0 0 0 ./’ ! " " # #) 莫尔 ! 库仑等面积圆屈服准则表示为 & - " "$ ./’ # 其中%- " , %" ! %- # * $ % - , (" 0 012. ./’ ! " & & ’ ’ ./’ #) " 0 " 12. " # 。 表达式中的符号含义同前。 012. " & & ’ ! ./’ ’ ./’ # % 因为( " %* " * " %%即%* " ( % * " (% - , &* %- , & ( ! %- # * $ % - ) , 则 * "( &! 3 0 $ ./’ # "( , & #) 0" 0 0 ! ./’ # " !%* # * $ % * " !( %- # * $(% ( ) ! %- # * $ % "( " (。 ! %- # * $ % -
(""# $ "+ $ ", ! 收稿日期:
现求得的安全系数大小与程序采用的屈服准则密切相
万方数据 关, 不同的准则得出不同的安全系数。
3::






-LL- 年
图! !"#$ %
有限元网格划分 !"&"’( ()(*(&’ *(+, 图 % !平面上不同!, $ 的屈服曲线 !"#$ G B,( H"()2 +51/86( 0& ’,( 2(I"8’01"6 9)8&(
./)0 12345678, ./9:; <8456=>4, 1/? @>86A84, @):; B8456C84
( !"#$%&$’() *+#$+,,-$+# .+$/,-%$&0,12"+#3$+# !"""!#, 12$+()
!/&+-#,+: 4+ (+()0%$% "+ %(5,&0 5(’&"- "5 %)"6, &2-"7#2 ! $ ! -,87’&$"+ ()#"-$&29 :0 5$+$&, ,),9,+&% $% 6-,%,+&,8 D ;2,+ &2, %0%&,9 -,(’2,% $+%&(:$)< $&0,&2, +79,-$’() +"+6 ’"+/,-#,+’, "’’7-% %$97)&(+,"7%)0 D =2, %(5,&0 5(’&"- $% &2,+ ":&($+,8 :0 ! $ ! -,87’&$"+ ()#"-$&29 D =2, 5(’&"- $% -,)(&,8 &" &2, 0$,)8 ’-$&,-$"+ D =2$% 6(6,- 6-,%,+&,8 ( ’"96(-$%"+ "5 %,/,-() 0$,)8 ’-$&,-$"+% $+ ’"99"+ 7%, (+8 8,87’,8 &2, %7:%&$&7&$/, -,)(&$"+%2$6 "5 &2,9 D >"- ’"+/,+$,+’, &2, ?"2-61"7)"9: ’-$&,-$"+ $% -,6)(’,8 :0 ?"2-61"7)"9: ,37$/(),+& (-,( ’$-’), ’-$&,-$"+,@2$’2 @(% 6-"6"%,8 :0 6-"5,%%"- A7 B(+’2,+# (+8 C2,+# D$+#-,+ $+ #,," D =2-"7#2 ( %,-$,% "5 ’(%, %&78$,%,&2, %(5,&0 5(’&"- "5 >*? $% 5($-)0 ’)"%, &" &2, -,%7)& "5 &-(8$&$"+() )$9$& ,37$)$:-$79 9,&2"8 D =2, (66)$’(:$)$&0 "5 &2, 6-"6"%,8 9,&2"8 @(% ’),(-)0 ,E2$:$&,8 D 7*% 8(-2&: %)"6, %(5,&0 5(’&"-; 5$+$&, ,),9,+&%; 0$,)8 ’-$&,-$"+%
传统的极限平衡法采用莫尔 $ 库仑准则, 但因莫尔
!


!
目前, 研究边坡稳定性的传统方法主要有: 极限平 衡法, 极限分析法, 滑移线场法等。这些建立在极限平 衡理论基础上的各种稳定性分析方法没有考虑土体内 部的应力应变关系, 无法分析边坡破坏的发生和发展 过程, 没有考虑土体与支挡结构的共同作用及其变形 协调, 在求安全系数时通常需要假定滑裂面形状为折 线、 圆弧、 对数螺旋线等。而有限单元法不但满足力的 平衡条件, 而且考虑了材料的应力应变关系, 使得计算 结果更加精确合理。在有限元法中通过强度折减, 使 系统达到不稳定状态, 有限元计算不收敛, 此时的折减 系数就是安全系数, 这种方法在国外 +" 年代就采用, 但由于力学概念不十分明确, 而且要受到计算程序及 计算精度的影响, 因而这种方法至今没有在国内流行。 随着计算机技术的发展, 尤其是岩土材料的非线性弹 塑性有限元计算技术的发展, 出现了许多适合于岩土 材料的大型通用有限元软件, 其前、 后处理的功能越来 越强大, 为利用有限元法进行边坡稳定分析创造了条
[#] , 通过计算发 件。本文引用有限元强度折减系数法
$ 库仑准则的屈服面为不规则的六角形截面的角锥体 表面, 存在尖顶和棱角, 给数值计算带来困难。为与传 统方法比较, 本文采用了徐干成、 郑颖人 (#,,") 提出的
[(] 莫尔 $ 库仑等面积圆屈服准则代替莫尔 $ 库仑准则 ,
并导出各准则间的换算关系, 由此可将求得的安全系数 折算成莫尔 $ 库仑等面积圆屈服准则下的安全系数。
!E
图# !"#$ 3
水平方向位移等值云图
40&’"&505+ 60&’051+ 0/ ,01"70&’8) 2"+9)86(*(&’
用此准则时, 塑性区最 国内特指此圆为 . D A 准则, 大。 (3) 当!, $ 满足下列表达式时 +"& " (! ) 3 360+ +"& +"& D ! # # $ $ " 3 % 60+ ! " $ E 360+ +"& +"& D ! # # $ $ "
% -
有限元软件 ;<=>=, 以及美国 ?=4 公司的 ?;@4、 A;BC
[ ] 在国外 @;<、 <;=B@;< 均采用了广义米赛斯准则 3 , 被称为德鲁克 D 普拉格准则 (. D A 准则) 。
这是徐干成、 郑颖人 (%KKL) 提出的莫尔 D 库仑等 [-] 它的面积等于不等角六边形莫尔 面积圆屈服准则 , 它比当前采用的逼近不等角的近似 D 库仑屈服准则, 屈服曲线有更高的计算精度。
(%) 当!, $ 满足下列表达式时 -+"& " (3 D +"& ) 3 ! " (-) J % 60+ " $ E 3 3 D +"& !( ") 屈服面在 !平面上为不等角度的六边形的外接圆。
!E
图" !"#$ -
变形后的网格图 .(/01*(2 *(+,
(-) 当!, $ 满足下列表达式时 +"& " 3 3 F +"&-") !( (3) 3 % 60+ ! " $ E !3 F +"& " 屈服面在 !平面上为不等角度的六边形的内切圆, 在
相关文档
最新文档