危岩崩塌落石运动计算excel
危岩破坏后的运动计算浅析
危岩破坏后的运动计算浅析孙智慧;唐勇;刘涛【期刊名称】《内江科技》【年(卷),期】2017(038)012【总页数】2页(P58,85)【作者】孙智慧;唐勇;刘涛【作者单位】内江职业技术学院土木工程系;内江职业技术学院土木工程系;内江职业技术学院土木工程系【正文语种】中文危岩体发生破坏，其运动方式受下部斜坡的物质组成、边坡坡角等的影响，运动距离各不相同。
本文主要对危岩破坏后的运动计算进行了浅析，为危岩破坏后的运动计算总结了相应的计算公式与理论依据。
根据R·M·Spang（1978）的研究成果，崩落体只有坡度角小于一定临界值（约27°）时，才停积于崖脚，随坡度角增大，可分别表现为滑动、滚动、跳跃和自由崩落等方式，大部分或全部堆积于坡脚。
勘查区内危岩崩落、运动的斜坡坡度区间值为35°～80°，因此岩体在产生变形破坏后，其运动方式表现为滚动、跳跃或自由崩落的方式向坡脚运动，最后堆积于坡脚缓坡地带，直接影响坡下乡村公路和居民住房及人身的安全。
石块在斜坡上的运动形式是比较复杂的，既有滑动、滚动还有跳跃运动，甚至在整个运动过程中三者兼而有之。
根据能量守恒定律，在物体下落过程中动能的增加等于势能的减少，机械能的总量保持不变。
即：根据地形剖面可计算出斜坡坡度β和碰撞时的切向速度Vt与法向速度Vn，即：落石与斜坡松散层坡面的法向碰撞可认为是塑性碰撞，所以Vn=0。
切向碰撞参考Hungr等人的研究，切向损失率采用10%，即落石第一次在斜坡上碰撞后维持其继续运动的动能为1/2 m（0．9Vt）2 。
块石在斜坡上的继续运动是以滚动和滑动为主的综合形式运动，其摩擦角称为综合摩擦角。
根据功能原理，落石的势能变化等于动能变化和克服摩擦所做的功：式中：Vi为落石在斜坡面上任意位置处所具有的速度；di为各直线段斜坡的平均坡度；△hi为各直线段斜坡的铅直高度；Φ为落石与坡面的综合摩擦角；Li为各直线段斜坡的长度。
边坡落石轨迹计算例题
边坡落石轨迹计算例题根据本项目地质勘查报告，边坡稳定性系数为 1.067~1.124，处于欠稳定~基本稳定状态。
在暴雨及地震情况下，边坡稳定性系数较小，可能失稳。
边坡失稳后，其崩落体将沿坡面滑落，并堆积于坡脚，下面采用rockfall软件计算其崩落体滑移过程对挡土堤的影响。
崩落体在坡面的运动模式有两种，一种是滚动或滑动模式，另一种为弹跳模式，如1所示。
图1 落石（土）坡面运动模式图采用rock fall软件，可以对边坡崩落体进行速度、弹起高度和落点等运动特征进行模拟分析计算，并自动计算其统计学规律，从而确定拦截系统的位置和高度。
该模拟计算分析，基于以下假定：1、由于失稳前的位移长期积累，崩落体的初速度一般很小，以初速度为零作自由落体运动，模拟计算时不考虑随机因素的影响。
2、采用偏保守的方法模拟时未考虑崩落体可能摔碎、破裂的情况。
3、假设崩落体重量对落点的路径和抛射距离无影响。
坡面为植被稀少或裸露的硬质土体，计算状态参数如表6-1所示。
崩落土体模拟的初始状态参数为：初始水平线速度为0m/s，初始垂直线速度为0m/s，初始角速度为0rad/s；落石块为球形、弹性体；重量以5000kg估算。
表1 坡面状况参数表DMSZDZSJ-04）为例，模拟分析结果如下图6-2～6-4。
图2 落点路径计算模型图图3 落点弹射轨迹及弹射高度图图4 落点弹射轨迹及总运动能量图由上述图中可以看到，从坡体顶部落下的崩落体运动轨迹都经过坡体中部平台（长度约7m）向下运动，在平台上的弹跳最大高度不超过1.5m，运动至浆砌块石挡土堤处高度远远小于0.5m；落点在平台坡脚处总运动能量最大（约700KJ），运动过程中被平台上填土逐渐吸收，至浆砌块石挡土堤处能量远远小于50KJ。
故本次设计浆砌块石挡土堤距离坡脚平均距离约7.0～7.5m，墙高约1.5m，墙后填土厚度约1.0m，墙体上余0.5m高度用于拦截崩落的土体。
根据珠海区域地质灾害的发育规律及发生机制，本市地质灾害一般发生在下雨过程中或之后的若干小时，并且以浅表层崩滑为主，方量一般不大。
危岩稳定性分析及崩塌落石计算
危岩稳定性分析及崩塌落石计算许可【摘要】随着我国经济的快速发展，公路、铁路等工程的兴建，崩塌灾害日益显著。
危岩稳定性分析及崩塌落石计算的准确性，对防治工程设计起着至关重要的作用。
根据工程实例，对危岩进行稳定性分析并对崩塌落石进行计算，为防治工程设计提供依据。
%Collapse rockfall and unstable rock mass threaten the security of communication installation. This paper deals with stability analysis of unstable rock mass and collapse rockfalls calculation by the example of concrete engineering.【期刊名称】《四川地质学报》【年(卷),期】2016(036)004【总页数】4页(P609-612)【关键词】崩塌;落石;稳定性;计算【作者】许可【作者单位】中冶成都勘察研究总院有限公司，成都 610023【正文语种】中文【中图分类】P642.21崩塌是较陡斜坡上的岩体在重力作用下突然脱离母体崩落、滚动、堆积于坡脚（或沟谷）的地质现象。
危岩是发生崩塌的初始物质条件，是由多组结构面切割并位于陡崖或陡坡上稳定性差的岩石块体及其组合。
危石指已脱离母岩，残留于坡面并有滚落可能性的块石。
随着我国经济的快速发展，公路、铁路等工程的兴建，崩塌灾害日益显著。
危岩稳定性分析及崩塌落石计算的准确性，对防治工程设计起着至关重要的作用。
本文根据工程实例，对危岩进行稳定性分析并对崩塌落石进行计算。
某工程于坡脚铺设管道，北侧山体基岩裸露，且坡体残留大量危石，该区域呈带状分布于山体中上部。
其变形迹象主要表现在岩体受构造及风化作用发育多条裂隙，在前缘临空条件下，产生卸荷裂隙，岩体经结构面切割成块，在外部营力作用下失稳并发生崩塌。
对三类危岩崩塌后影响斜坡稳定性的定量计算 (1)
对三类危岩崩塌后影响斜坡稳定性的定量计算摘要：三类不同运动轨迹的危岩与斜坡撞击后对斜坡稳定性的影响不同。
本文通过刚体运动学的理论知识，将三类危岩的崩塌体与斜坡作为一个系统进行研究，应用质心定理，能量守恒定律以及动量定理分别对三类危岩崩塌体与斜坡构成的系统的稳定性作了定量计算，最后给出了每个系统最终滑移距离的计算公式。
关键词：危岩斜坡定量计算1 前言危岩是指位于岩质陡坡或陡坡的崩塌源被结构面切割且稳定性较差的岩块体。
外力的作用，如地震作用，人工爆破和分化作用等使得危岩体后部主控结构面失稳断裂和贯通，大块岩体或岩石群突然从陡坡坠落。
危岩体失稳破坏的这个过程也称之为崩塌。
危岩崩塌是山岭地区最主要的一种地质灾害现象。
大量的危岩崩塌体突然从陡坡坠落，崩塌体在向下的运动过程中，垂直运行的距离远远大于水平运行距离，大块的危岩体或群体在重力作用下，获得了巨大的能量。
当不稳定斜坡受到危岩崩塌体的冲击后，危岩崩塌体的动力作用就成为了斜坡失稳的起搏器，诱使其形成崩塌滑坡。
滑坡的滑移距离能否危及该地区人民的生命安全是我们最为关心的问题。
鉴于此，具体定量的分析各类危岩崩塌体对斜坡的稳定性的影响就显得非常重要。
根据陈洪凯，唐红梅等人对危岩具体研究，可将危岩体化分为以下类：（1）坠落式危岩（2）倾倒式危岩（3）滑塌式危岩，据实地调查，陈洪凯，唐红梅等人对危岩类型的划分符合实际情况。
具体分析这三类危岩运动轨迹后发现，危岩与斜坡撞击后运动轨迹受到斜坡地形地貌和崩塌体自身形状等因素影响，很难准确地予以宏观测定以及类比分析其运动轨迹。
故采用多刚体运动学把崩塌体与斜坡作为一个系统进行研究，对这三类危岩崩塌体应用质心定理，功能转化原理以及动量定理进行定量分析计算，并认为这一细化的定量分析方法基本可信，可以为防灾治灾工作提供计算依据。
2计算过程分析2.1 坠落式危岩—斜坡系统联合运动分析计算坠落式危岩—斜坡系统：高悬于陡崖上端和岩岩腔顶部的岩体受裂隙切割脱离母岩，下部受结构面切割脱离母岩，上部及后部母岩尚未脱离，在重力作用下基本不受阻力便失稳崩塌冲击陡崖下的不稳定斜坡后联合运动。
矿山拦石墙缓冲厚度、危岩崩塌支撑柱(墙)反力、棚洞设计、水文与水力、采空塌陷防治范围计算
坏面为矩形时可取 1/6； e ——地震作用水平系数； Q ——地震力（kN/m）；
H ——后缘裂隙上端到未贯通段下端的垂直距离（m）； h ——后缘裂隙深度（m）； a0 ——崩塌体重心到潜在破坏面的水平距离（m）； b0 ——崩塌体重心到潜在破坏面形心的铅垂距离（m）； l ——柱撑体距离主控裂隙面在危岩底部出露点的水平距离（m）； k ——崩塌体抗拉强度标准值，单位为千帕（kPa），根据岩石抗拉強度标准值乘以 0.20
R1 Fst G Q tan cH 0
（I.5）
R2
Fst (Ga0
Qb0 ) l
k
H
2 0
（I.6）
R Max(R1，R2 )
（I.7）
式中： H0 ——崩塌体后缘潜在破坏面高度（m）；
k ——崩塌体抗拉强度标准值（kPa），根据岩石抗拉强度标准值乘以 0.30 折减系数确定；
其他符号意义同前。
1 n
mn
( n) cos (1 n)sin (m n)
对左侧（回填土石坡面向下倾斜）侧压力系数计算为：
1 n
mn
( n) cos (1 n)sin (m n)
当回填土石坡面水平时，即 m ，且侧压力作用方向水平，即 0 时：
186
（J.8） （J.9）
DB 32/×××.3—2020
0.15
极坚硬：易被拇指指甲划伤
0.20
坚固：难于被拇指指甲划伤
0.20
极软岩：可被拇指指甲划伤
0.15
较软岩：地质锤尖击打可破碎，易被小刀切削
0.15
kt
0.50 0.55 0.65 0.70 0.75 0.80～0.85 0.75 0.75
危岩体稳定性分析
———————————————————————————————— 作者:
———————————————————————————————— 日期：
ﻩ
附件２ 危岩体稳定性分析
１、WY-０１危岩体稳定性定量评价
1计算模型
从工程防治的角度按照危岩失稳类型进行分类，可将危岩概化分为滑移式危岩、倾倒式危岩和坠落式危岩3类。WY-01危岩体为滑移式危岩；其软弱结构面倾向山外，上覆盖体后缘裂隙与软弱结构面贯通,在动水压力、地震和自重力作用下,缓慢向前滑移变形,形成滑移式危岩，其模式见图（图３-1)。
图3-３ 危岩崩塌破坏运动图示
根据落石的运动情况，可以分为两种状态:启动阶段、运动阶段。
1启动阶段
滑移（错断）式危岩体附着于母岩上，以一定角度的裂隙面相接,在危岩体自重和地表水渗入裂隙等因素的作用下,裂隙面锁固部位被贯通,危岩体沿母岩(或基岩)发生剪切滑移破坏。如图３-４所示。
图3－４滑移式破坏初始运动状态
WＹ-０1
滑移式
１.6５
1．37
1.3６
1.13
未贯通
1．3９
1.１4
1.18
０.9４
后缘切割面贯通4０%，暴雨时完全充水
1.3３
１.09
1.１３
0.பைடு நூலகம்0
后缘切割面贯通５0%，暴雨时完全充水
1．38
1.1２
1.１７
0.９3
后缘切割面贯通６０%,暴雨时完全充水
1.２1
0．9８
1.0３
0.81
后缘切割面贯通70％，暴雨时完全充水
1．15
0.９３
0.98
０.77
后缘切割面贯通8０%,暴雨时完全充水
excel计算崩塌公式
excel计算崩塌公式
计算崩塌公式需要更具体的信息，因为不同的崩塌模型和情况会有不同的计算公式。
以下是一个常见的崩塌计算公式示例：
1. 找到所需的参数：根据崩塌模型，确定需要的参数，如土壤重度（γ）、坡度（θ）、抗剪强度（c）等。
2. 计算切线力：使用公式 T = γ乘 h乘 sin(θ)，其中 T 表示切线力，γ表示土壤重度，h 表示坡高，θ表示坡度。
3. 计算抗剪强度：根据具体情况使用相应的抗剪强度计算公式，如考虑土壤内聚力和摩擦力的情况下可以使用τ = c + σ乘 tan(φ)，其中τ表示抗剪强度，c 表示土壤的内聚力，σ表示正应力，φ表示土壤的内摩擦角。
4. 判断稳定性：比较计算得到的切线力和抗剪强度，如果切线力大于抗剪强度，则表示存在崩塌的可能性。
请注意，以上仅为一个示例，实际的崩塌计算公式会根据具体情况有所不同。
在进行崩塌计算时，建议参考专业的地质工程文献或咨询相关领域的专业人士以获取准确和可靠的计算方法。
危岩稳定性计算表格-滑移式-倾倒式-坠落式-完整版
后缘裂隙深度（h）（m）
裂隙水高度（裂隙1/3）(hw)（m） 0.00 危岩 的破 后缘裂隙上端到未贯通段下端的垂直 坏由 距离(H)(m) 底部 危岩体重心到倾覆点的水平距离(a)(m) 危岩体与基座接触面倾角(α )（° ） 岩体 危岩体重心到倾覆点的垂直距离(h0)(m) 抗拉 水容重（kN/m） 9.8 强度 岩石质量(W)（kN·m） 0.0 控制 地震水平系数（ζ e） 0.05 地震力（Q)（kN·m） 0.00 危岩抗弯力矩计算系数(ζ ) 后缘 后缘裂隙上端到未贯通段下端的垂直距离 有陡 (H)(m) 倾裂 重心到潜在破坏面的水平距离(a0)（m） 隙的 悬挑 式危 坠 岩 落
后缘 有陡 倾裂 隙的 重心到过潜在破坏面形心的铅垂距离(b0)(m) 悬挑 地震水平系数（ζ e） 0.05 式危 地震力（Q)（kN·m） 0 坠 岩 稳定性系数(K) 落 式 后缘 危岩抗弯力矩计算系数(ζ ) 有陡 危岩体后缘潜在破坏面高度（H0）（m） 倾裂 重心到潜在破坏面的水平距离(a0)（m） 隙的重心到过潜在破坏面形心的铅垂距离(b0)(m) 悬挑 地震水平系数（ζ e） 0.05 式危 地震力（Q)（kN·m） 0 岩 稳定性系数(K)
0 9.8
#DIV/0!
0 9.8
#DIV/0! 9.8
抗拉强度标准值(flk)(kPa)
重力加速度(m/s)
后缘裂隙未贯通段下端到倾覆点 之间的水平距离(b)(m) 危岩体与基座接触面倾角(α )（° ） 后缘裂隙倾角(β )（° ）
岩石容重（kN/m） 岩石体积（m³/m） 裂隙水压力（V）（kN·m） 稳定性系数(K)
0 #DIV/0! 9.8
Байду номын сангаас
危岩抗拉强度标准值(flk)(kPa)
危岩稳定性计算表格-滑移式-倾倒式-坠落式-完整版
后缘裂隙深度（h）（m）
裂隙水高度（裂隙1/3）(hw)（m） 危岩 后缘裂隙上端到未贯通段下端的垂直 体重 距离(H)(m) 心在 危岩体重心到倾覆点的水平距离(a)(m) 倾覆 危岩体重心到倾覆点的垂直距离(h0)(m) 点之 水容重（kN/m） 外 岩石质量(W)（kN·m） 地震水平系数（ζ e） 地震力（Q)（kN·m） 0
#VALUE!
#DIV/0!
其他格子无需改动，黄色格子为稳定性系数 裂隙水高度（裂隙1/3）（m） 岩石体积（m³/m） 裂隙水压力（V）（kN·m） 重力加速度(m/s) 结构面倾角（α ）(°） 后缘结构面摩擦角（φ ）(°） 稳定性系数（K） 裂隙水高度（裂隙1/3）（m） 岩石体积（m³/m） 裂隙水压力（V）（kN·m） 重力加速度(m/s) 结构面倾角（α ）(°） 后缘结构面摩擦角（φ ）(°） 稳定性系数(K)
后缘裂隙深度（h）（m）
裂隙水高度（裂隙1/3）(hw)（m） 0.00 危岩 的破 后缘裂隙上端到未贯通段下端的垂直 坏由 距离(H)(m) 底部 危岩体重心到倾覆点的水平距离(a)(m) 危岩体与基座接触面倾角(α )（° ） 岩体 危岩体重心到倾覆点的垂直距离(h0)(m) 抗拉 水容重（kN/m） 9.8 强度 岩石质量(W)（kN·m） 0.0 控制 地震水平系数（ζ e） 0.05 地震力（Q)（kN·m） 0.00 危岩抗弯力矩计算系数(ζ ) 后缘 后缘裂隙上端到未贯通段下端的垂直距离 有陡 (H)(m) 倾裂 重心到潜在破坏面的水平距离(a0)（m） 隙的 悬挑 式危 坠 岩 落
后缘 有陡 倾裂 隙的 重心到过潜在破坏面形心的铅垂距离(b0)(m) 悬挑 地震水平系数（ζ e） 0.05 式危 地震力（Q)（kN·m） 0 坠 岩 稳定性系数(K) 落 式 后缘 危岩抗弯力矩计算系数(ζ ) 有陡 危岩体后缘潜在破坏面高度（H0）（m） 倾裂 重心到潜在破坏面的水平距离(a0)（m） 隙的重心到过潜在破坏面形心的铅垂距离(b0)(m) 悬挑 地震水平系数（ζ e） 0.05 式危 地震力（Q)（kN·m） 0 岩 稳定性系数(K)
云台山景区公路边坡危岩体稳定性计算及落石运动轨迹研究
云台山景区公路边坡危岩体稳定性计算及落石运动轨迹研究作者：郭龙龙耿国建丛颖来源：《西部资源》2017年第04期摘要：介绍云台山景区内公路边坡危岩体的稳定性计算，并根据计算结果应用Rock Fall 软件模拟计算危岩体失稳后落石的运动轨迹，为危岩体的防治提供科学依据。
以编号为TW23的危岩体为例，采用静力计算的方法分析计算其在不同工况下的稳定系数，采用数值模拟软件研究危岩体失稳形成落石后的运动距离、速度、弹跳高度、冲量等运动特征，根据数值模拟结果，选取合适的防治措施，将危岩体失稳后可能造成的损害降到最低，同时为景区内其他部位危岩体的治理提供依据。
关键词：云台山；危岩体；落石；稳定性；运动轨迹引言危岩体是指发育在边坡斜体上的随时可能在各种因素作用下发生失稳破坏的岩体，具有突发性、速度快、冲击力大的特点，是山区常见的一种地质灾害。
云台山景区位于河南省焦作市境内，西北方向与山西省晋城市接壤，是我国首批5A级世界地质公园之一。
近年来景区内公路上发育的危岩体给景区的正常运行带来了安全隐患，通过现场调查，选取编号为JW23危岩体为研究对象，计算其稳定性并研究失稳后所形成的落石运动轨迹，最终为危岩体的防治提供建议。
JW23位于前往叠彩洞景点的盘山公路边坡的顶部，大体呈柱状，长2.1m，高4.6m，宽1.5m，体积14.49m3，质量38.7t，倾角65°，岩性为灰岩，危岩后部被裂隙切割，基本贯穿整个危岩体，坡面角度30°，植被生长较茂盛。
根据重庆市地方标准《地质灾害防治工程设计规范》（DB50/5029-2004）对该危岩体进行分类，属于滑移型危岩体。
1.稳定性评价标准根据陈洪凯等（2011）在《地质灾害理论与控制》一书中有关危岩体稳定性评价标准内容，可依据危岩体稳定性系数将危岩体稳定性分为：不稳定、基本稳定、稳定，对应的具体稳定性评价标准见表1。
2.稳定性计算2.1基本假设根据前人的研究成果，滑移式危岩体稳定性计算依据以下假设条件：危岩体变形发展过程中，尤其是在其破坏失稳运动以前，将危岩体视为刚体；把复杂的空间运动问题简化成平面问题；危岩体与稳定坡体之间无摩擦力。
边坡落石运动轨迹计算新方法
边坡落石运动轨迹计算新方法边坡落石是指悬崖或陡坡上个别岩块在自重或外力作用下突然脱离母体而急剧下落的现象。
落石是一种突发的地质灾害，虽然规模不如滑坡巨大，但落石发生往往很突然，同样会造成巨大损失。
落石运动是一种复杂的运动，影响边坡坡面的地质与性质、落石大小、形状等。
落石本身及地面条件是多样的，落石的运动形式也是多种多样，要实现对落石的运动轨迹精确计算有很大的困难。
现实中落石与边坡碰撞的受力情况有很大的不确定性。
落石运动可分为5种形式：自由落体、斜抛运动、碰撞、滑动和滚动等运动方式的一种或多种组合沿着坡面向下快速运动。
最后在较平缓的地带或有障碍物附近静止下来的一个动力学过程。
1.自由落体当岩块脱离母体，只在重力作用下自由坠落时，或者在只有一个初始角速度时，危石突然倾倒破坏后坠落，忽略空气阻力作用，角速度在坠落过程中保持不变，则落石的速度为v。
式中：1H为自1由落体高度，g为重力加速度。
2.斜抛运动阶段当落石在碰撞后法向恢复系数和切向恢复系数＞0时，落石将离地做斜抛运动，落石做斜抛运动的危害极大。
工程设计时常采用阻挡拦截结构物。
因此；需要确定落石做斜抛运动的弹跳最大高度和最大水平运动距离。
落石做斜抛运动时没有和地面接触，忽略空气阻力影响，可以根据能量守恒得到弹跳的最大高度为：23r 1t211h =v cos -v sin 2gθθ（） 式中3h 为弹跳最高点与坡面碰撞点的竖直距离。
最大水平运动距离为→2r 1v vt gL ⎡=⎢⎣ 式中1L 为起跳点与下一个碰撞点的水平距离。
3. 滑动阶段当落石与坡面碰撞都得法向速度为0时而切向速度不为0，且r d =0时，落石会沿坡面滑动，下滑速度为4vt 式中：f 为滑动摩擦系数；2l 为滑动的距离；3vt 为碰撞后的切向初速度；4vt 为滑动后的切向速度。
如果已知滑动后坡面的坡度不变且坡面较长时，可根据上式计算出落石滑动停止时所经过的切向距离。
4. 滚动阶段（椭圆落石）当落石与坡面碰撞后法向速度为0，而r d ﹤0，切向速度反向时，落石会向后滚动。
新疆库尔勒铁门关崩塌危岩体稳定性研究
新疆库尔勒铁门关崩塌危岩体稳定性研究作者：***来源：《新疆地质》2023年第04期關键词：崩塌；危岩；GIS；RocFall崩塌是一种多发于山区的地质灾害，具点多面广、快速突然和随机分布特点。
在地质作用及外界作用综合影响下，崩塌运动方式包括一块或多块岩石离开原来位置，随后以坠落、滑动、滚动、弹跳等方式中的一种或多种运动形式组合[1-2]。
新疆地区崩塌地质灾害较普遍，前人进行了大量研究工作，主要集中在崩塌地质灾害发育特征及成因机制研究、崩塌稳定性分析、崩塌运动学特征研究及崩塌灾害危险性评估等方面[3-15]。
新疆库尔勒铁门关景区为国家“AAA”级风景旅游区，在旅游旺季时，景区道路车流量可达2～4辆/分钟，人流量10～20人/分钟。
景区道路位于霍拉山低中山孔雀河峡谷中，两侧山坡陡峭，岩体结构面发育，易发生崩塌落石地质灾害，威胁来往车辆和行人生命财产安全。
因此，亟需对新疆库尔勒铁门关崩塌危岩体进行稳定性分析。
1 工程概况1.1 地理位置研究区位于新疆库尔勒市天山街道双拥社区铁门关生活片区北侧，铁门关景区公主墓西侧区域，中心地理坐标东经86°11′28.83″，北纬41°49′6.44″。
研究区距库尔勒市区约9 km，有柏油道路相连接，交通十分便利。
1.2 气象水文研究区属暖温带大陆性荒漠气候，冷热差异悬殊，温度的年、月变化大，最热月与最冷月平均气温差约36 ℃，日气温变化差平均17 ℃～18 ℃。
区内年平均降水量61.26 mm，年平均蒸发量1 674.56 mm。
研究区西部为孔雀河，发源于博斯腾湖，流经焉耆盆地西南，由东向西南经铁门关峡谷，全长785 km，年平均径流量11.872 4×108 m3。
1.3 工程地质条件研究区地处霍拉山低中山区，微地貌为孔雀河下切冲蚀霍拉山形成的峡谷地貌，地貌类型分为两类，分别为构造侵蚀低中山地貌和侵蚀峡谷地貌。
区内地层主要为上元古界震旦系和第四系。
危岩体稳定性分析汇总
附件2危岩体稳定性分析1、WY-01危岩体稳定性定量评价1计算模型从工程防治的角度按照危岩失稳类型进行分类，可将危岩概化分为滑移式危 岩、倾倒式危岩和坠落式危岩 3类。
WY-01危岩体为滑移式危岩；其软弱结构 面倾向山外，上覆盖体后缘裂隙与软弱结构面贯通， 在动水压力、地震和自重力 作用下，缓慢向前滑移变形，形成滑移式危岩，其模式见图（图3-1）。
危岩后缘2计算公式①后缘有陡倾裂隙、滑面缓倾时，滑移式危岩稳定性按下式计算:K(Wcos Qsin V sin V)tg c l K 二后缘裂隙图3-1滑移式危岩示意图ne地下水位危岩体 扬压力U危岩前缘9Wcos e : 7W轨弱结构面静水压力V=h图3-2滑移式危岩稳定性计算示意图（后缘有陡倾裂隙）Wsin 日+ Q COST+ V cos6式中：V——裂隙水压力(kN/m),,hw ――裂隙充水高度(m)，取裂隙深度的1/3。
yw——取10kN/m。
Q——地震力(kN/m)，按公式Q = e W确定，式中地震水平作用系数七级烈度地区e取0.075;K——危岩稳定性系数；c――后缘裂隙粘聚力标准值(kPa)；当裂隙未贯通时，取贯通段和未贯通段粘聚力标准值按长度加权和加权平均值，未贯通段粘聚力标准值取岩石粘聚力标准值的0.4倍；――后缘裂隙内摩擦角标准值(kPa)；当裂隙未贯通时，取贯通段和未贯通段内摩擦角标准值按长度加权和加权平均值，未贯通段内摩擦角标准值取岩石内摩擦角标准值的0.95倍；二一一软弱结构面倾角(°，外倾取正，内倾取负；W——危岩体自重(kN/m3)。
3危岩稳定性计算结果根据危岩结构特征和形态特征，②区危岩破坏模式主要为滑移式。
(1)计算参数：崩塌区出露地层为第四系崩坡积物和石炭系太原组，根据附近工程岩体参数及工程类比得出物理力学参数见表：注：由于坡表白云岩、灰岩多为强〜弱风化强卸荷岩体，其参数均参考类比相似强〜弱风化强卸荷岩体参数。
危岩崩塌落石稳定性运动计算总表
F-F
Ⅰ-Ⅰ 1-2危岩
1-1危岩
Ⅱ-Ⅱ 2-1危岩
坡度α
82 34 67 13
48 40 31
57 29
58
38 48 35 24
53 33
53 57 51 41 28 36 31
48 32 45 37
22.7 212.5 113.85
57 36.7
0.57345 0.627168 0.696762 0.567072 0.5175
1.429089737 1.111442854 0.870048201 1.483523463 2.145745998
2-2危岩 Ⅲ-Ⅲ
Ⅳ-Ⅳ Ⅴ-Ⅴ
A-A B-B C-C D-D E-E Ⅰ-Ⅰ Ⅱ-Ⅱ
810 34.44512917 65.9148 810 87.39405048 424.3178 810 106.8517638 634.2944 810 110.0617933 672.9777
810 96.79035455 520.4652 810 100.8051537 564.5377
cotβ
2160 77.83511833 126.2147 2160 92.94850368 179.9880
13.66542288 0.2 24.06965571 0.2 36.46992413 0.2 39.73597972 0.2 39.74587404 0.2 42.63369652 0.2 43.49288715 0.2
2700 2700 2700 2700
2160 2160 2160 2160
122.2805984 172.3342859 189.8534329 183.1415298
危岩崩塌落石稳定性运动计算总表(秦皇岛资源环境勘察院)
落石腾跃高度计算使用说明：
初始速度：输入你所要求解断面中坠 落后的速度。例如：A-A断面，AB段落 在B点后的速度。最大拦截高度在于初 始速度大小。用于设置挡墙位置与高
度。
反射速度V0 8.876515308
石块第一 次弹跳的 最大高度
Hmax
相应的最大 水平距离
Lmax
0.107179 7.35650433
30-60° K=0.543-0.0048α+0.000162α² 60-90° K=1.05-0.0125α+0.0000025α
落石运动计算使用说明：
坡度处输入各剖面分段坡度角。 垂直距离处输入分段内垂直高度。 崩塌摩擦系数K按右侧表取值。 落石体积自己输入 红色是需要自己输入的，蓝色是计算得到的。
48.15250372 25 2.924018 75.99893 0.466038
AB
坡脚弹跳 CD
崩塌 落石 弹跳 计算
λ
ρ
φ
tanφ cosφ
0.4
0.3
13 0.230747 0.974396
落石冲击力及缓冲填土层厚度的计算
四号危岩段设置当地碎块石 土作为缓冲材料，容重γ =21.75
最大块径=2*1.4*1
13.58 0.743658 0.75429 10.81043 18.74 0.561018 1.54085 12.90951
21.7 41.5 91.9 44.1 21.5 60.4 44.8
0.743658 0.75429 13.66542 0.795738 0.650125 24.06966 0.719562 0.810531 36.46992 0.618522 1.151212 39.73598 0.5304 1.881851 39.74587 0.580152 1.377304 42.6337 0.549882 1.665314 43.49289
落石撞击能量计算1
落石撞击能量计算一、不同边坡落石速度及坠落点的计算根据前苏联尼•米•罗依尼会维里教授提出的落石运动速度的计算方法，用电脑可能计算出 各种不同边坡形式的落石速度，进而计算出撞击能量及落石距边坡的距离。
为了便于分析，我们把边坡简化为以下四种形式进行分析。
1边坡分类1.1 I--单一坡度边坡：包括山坡是台阶式的，但各个台阶的高度小于 5m ，及山坡为折线,但其各段长度小于10m 或相邻坡度差在5°以内者。
1.2 II--缓折线形山坡：其中缓山坡的坡度角aV 30°,陡坡段坡度角aW 60°,坡段长超过10m ，相邻坡段的坡度角相差5°以上。
IV--直立边坡：边坡坡度角70 VaV 90，仅有一个台阶或没有。
2 落石运动速度计算2.1 I 型边坡计算公式：坠落石块沿单一山坡运动时的计算速度，可用任意形状物体滚动、 滑动、跳跃运动的公式，即__V=2gH=H(1-1)1.3 III--陡折线形山坡：上部坡段为极陡坡aV 其高度超过10m ， F 部坡段坡度较缓。
V 1V 2（o ）/1.4 H 1V RH 2aV t （0）” a口=V 1-Kctg a& =W2g (1-2)式中H---石块坠落高度(m ) ; g---重力加速度(m/s 2)) a —山坡坡度角(度);K---石块沿山坡运动所受一切有关因素综合影响的阻力特性系数。
卩、&值见表1。
2.2II型边坡计算公式：最高一个坡段坡脚的速度按公式(1-1) (1-2)计算，其余坡段终端的速度为：V j(i)= VV0(i)+2gH i(1-K i ctga i) =VV0(i) 2+ £i2H i (1-3)式中V o(i)---石块运动所考虑坡段的起点的初速度，可按下列不同情况考虑，若a (i-1 ) >a i时，则V o (i) =V j(i-1)C0S( a (i-1 ) - a i)；若％(i-1 )<a i 时，则V o (i) =V j(i-1)。
单体危岩崩塌定量计算方法与治理设计中的应用
单体危岩崩塌是地质灾害中一种常见的灾害类型，具有较大的危害性和不确定性。
在治理设计中，需要采用合适的方法对单体危岩崩塌进行定量计算，并据此设计治理方案。
本文将就单体危岩崩塌的定量计算方法和在治理设计中的应用进行探讨。
一、单体危岩崩塌的定量计算方法1. 单体危岩崩塌的形成机理单体危岩崩塌是指由于地质体内部存在一定的构造缺陷或岩层间的滑动面、节理面等内在缺陷，受外力作用下所引发的大规模岩石和土石体的破坏现象。
单体危岩崩塌的形成机理是多方面的，主要包括地质构造、地下水、外力作用等因素的综合影响。
2. 单体危岩崩塌的定量计算方法（1）地质勘察调查：对危岩体进行详细的地质勘察调查，确定构造缺陷、节理面等内在缺陷的分布情况。
（2）力学参数的测试：对危岩体进行力学参数的测试，包括岩石的强度、变形特性、破裂模式等参数的测定。
（3）数值模拟分析：采用数值模拟方法，对单体危岩崩塌的稳定性进行分析，包括有限元法、离散单元法等方法的应用。
（4）定量评价：根据地质勘察调查和力学参数测试结果，采用定量评价方法对单体危岩崩塌的潜在危险性进行评估。
二、单体危岩崩塌治理设计中的应用1. 危岩崩塌的治理原则（1）综合治理：危岩崩塌的治理应该采取综合治理的原则，包括加固、排除、隔离等多种手段的综合应用。
（2）科学设计：危岩崩塌的治理设计应该科学合理，充分考虑地质环境、工程条件等因素，设计出合适的治理方案。
2. 危岩崩塌的治理措施（1）围护结构：采用围护结构对危岩体进行加固，包括挡土墙、引导坡、挡石墙等结构的设计和施工。
（2）爆破除爆：对危岩体进行爆破除爆处理，消除潜在的危险因素。
（3）排水处理：对危岩体周围的地下水进行排水处理，降低地下水对危岩体稳定性的影响。
（4）监测预警：对治理后的危岩体进行监测预警，及时发现异常情况并采取相应的应对措施。
三、结语单体危岩崩塌的定量计算方法和治理设计在地质灾害治理中具有重要的意义。
通过对单体危岩崩塌的形成机理的研究，采用合适的定量计算方法，可以科学地评估危岩崩塌的潜在危险性，为治理设计提供科学依据。
危岩体计算模板
表3.6-1 W1危岩体特征、稳定性评价及整治方案建议说明表位置剖面1计算受力分析图aβαbvGOHhhwh危岩照片岩性变质砂岩基座岩性变质砂岩岩层产状273°∠28°分布高程（m）861.90～857.42 危岩形态不规则块体块体规模（m）（长×高×厚）3.3×3.75×1.5 危岩体积（m3）18.56主崩方向（°）76°崩塌方式滑移式最大垂直落差(m)136.90m最大水平落距(m)95.67m控制结构面及其描述①组裂隙，产状为87°∠55°，延伸3～4m，间距0.50～2.00m，张开度20～30mm，裂面粗糙，局部粘土或碎石充填；②组裂隙，产状为200°∠76°，延伸3～5m，间距3.00～4.00m，张开度5～20cm，裂面粗糙，局部粘土或碎石充填；③岩层产状：273°∠28°④边坡坡向：75°∠70°剖面、立面图比例尺：1:10075°51'结构面赤平投影分析图NESW1234①岩层产状：273°∠28°②裂隙LX1产状：87°∠55°③裂隙LX2产状：200°∠76°④边坡坡向：75°∠70°稳定性定性分析基本稳定稳定性定量分析计算参数γ天=23.6kN/m3,γ饱和=23.8kN/m3 ,S=29.56㎡ Q=34.88 kN/m，α=70°,L=3.35m，h w=2.12m。
计算结果及稳定性分析工况1 F=1.324工况2 F=1.190工况3 F=1.268稳定性综合评价工况1 稳定工况2 欠稳定工况3 基本稳定治理措施建议清除表3.6-2 W2危岩体特征、稳定性评价及整治方案建议说明表位置剖面2计算受力分析图aβαbvGOHhhwh危岩照片岩性变质砂岩基座岩性变质砂岩岩层产状273°∠28°分布高程（m）851.12～848.29 危岩形态不规则柱状块体规模（m）（长×高×厚）3.14×4.31×1.4 危岩体积（m3）8.44主崩方向（°）123°崩塌方式倾倒式(重心在内) 最大垂直落差(m)125.32m最大水平落距(m)60.50m控制结构面及其描述①组裂隙，产状为130°∠40°，延伸3～8m，间距1.00～5.00m，张开度10～20mm，裂面粗糙，无充填；②组裂隙，产状为200°∠78°，延伸3～5m，间距3.00～4.00m，张开度5～30cm，裂面粗糙，无充填；③岩层产状：273°∠28°④边坡坡向：123°∠55°剖面、立面图比例尺：1:100123°37'结构面赤平投影分析图NESW1234①岩层产状：273°∠28°②裂隙LX1产状：130°∠40°③裂隙LX2产状：200°∠78°④边坡坡向：123°∠55°稳定性定性分析稳定稳定性定量分析计算参数γ天=23.6kN/m3,γ饱和=23.8kN/m3 ,S=2.77㎡ Q=3.27 kN/m，α=66°，β=80°,f lk=108kPa， H=3.2m，h=3.19m，a=0.58m，b=0.79m,h0=1.61m。
RocFall软件在公路边坡危岩崩塌运动路径分析中的应用
由总动 能分 布 图可 以看 出： 在 落石 到 达 Ｘ一 １ ０ ｍ 处时 ， 总动 能 达到最 大 ， 约 １ ． ５ ×１ Ｏ ｋ Ｊ ， 碰 撞之 后
总 动能 急 剧 减 小 。总 动 能 与 平 移 速 度 图形 基 本 相 似， 几 乎 同时到 达最 值 、 同时增 大 、 同时 减小 ， 这 与其 计 算原 理密 切相 关 ， 同 时可 以推 理 出其 旋 转 速 度 也 与 总动 能 图相似 ， 在 此不列 举 。
拦 石墙 的位 置与 高度 。 Ｒ ｏ ｃ Ｆ ａ ｌ ｌ 软 件在 岩 土 工 程 问题 的处 理 中具 有 方
各 个 水平 位置 的 弹跳 高度分 布 。
便、 适用 的辅 助 价值 ， 运用 地质 勘查 资料 能方 便快 捷 地 体现 岩体 的崩 落路 径 及 运 动 模 式 ， 对 于危 岩 体 防
Ｅ ５ ３ 罗 田． 岩质边坡危岩落石 运动特 征和 防护研究 Ｅ Ｄ］ ． 成 都： 西南交通大学 ， ２ ０ １ ３ ． 收 稿 日期 ： ２ ０ １ ４ —１ １ —１ ９
研究 Ｅ Ｄ ３ ． 成都 ： 成都 理工大学 ， ２ ０ ０ ８ ．
［ ４ ］ 陈洪 凯 ， 唐红梅 ， 叶四桥 ， 等． 三 峡 库 区 危 岩 发 育 链 式 机 理及失稳 运动路径研 究 ［ Ａ３ ． 第 八 次 全 国 岩 石 力 学 与
工 程 学 术 大会 论 文 集 Ｅ ｃ ３ ． ２ ０ ０ ４ ．
图 ４ 障碍 设 置轨 迹 图
由 图 ５可 以看 出 ： 在 落 石 到 达 Ｘ 一８ ． ５ ｍ 处时，
弹跳 高度 最 大 ， 约 １ ９ ． ５ ｄ ｍ。该 点 之 后 ， 经过碰撞 、 滚 动 摩擦 等 ， 弹跳 高度 逐渐 变小 ， 直 至停 止 。表 ２为
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
1 2 0 0 2.93 3.37 1.51 0 0
0.1 0.3 0.1 0.3 0.1 0.1 0.3 0.1 0.3
18.06846 16.96015 11 8 18.43806 25.67645 20.09675 13.1
10剖面3#崩塌
第1段 第2段 第3段
跳跃 跳跃 跳跃 跳跃 跳跃 运动模式 0
0.5 0.5 0.5 0.5 0.5 0.5 0.5
0.3 0.3 0.3 0.3 0.3 0.3 0.3
0.463503 0.176555 0.0636 0.02273 0.008115
26.57021 10.12097 3.645881 1.302973 0.465181
0.434078 24.88346 0.163986 9.400499
0.5 0.5
0.3 0.3
0.26747 15.33268 0.09751 5.589726
13剖面
第4-1段 第4段-2
跳跃 跳跃
41 17.23557
0.5 0.5
0.3 0.3
0.555331 31.83426 0.217957 12.49434
9剖面2#崩塌
9剖面2#崩塌
第4-1段 第4段-2
跳跃 跳跃 跳跃 跳跃 跳跃 跳跃 跳跃 跳跃
41.3 44.6 36.5 0
57.6 42.4 80.1 0
2.36 2.56 2 0
0.1 0.1 0.3 0.3
17.91114 23.17041 24.04678 15.3
上部治理后 跳跃 0 0 0 0.3 11
落石撞击斜坡后最大偏离计算 运动模式 19剖面 2 计算坡段 坡段初始 开3次方 反射角β 坡度α 速度Vj （Vj） 61 34 55.9 6.2 22 16 tanα
跳跃 跳跃
36.7 14.89713
0.5 0.5
0.3 0.3
0.488959 28.02948 0.187682 10.75886
10剖面3#崩塌
第4-1段 第4段-2
跳跃 跳跃
31.8 12.47838
0.5 0.5
0.3 0.3
0.416681 23.88618 0.156693 8.982388
第2段 第2段
跳跃 跳跃 跳跃
1.837091 85.25558 1.801754 2.802039 77.30877 0.674071 2.519842 68.62317 1.475421
3剖面WYD2 第2段
第三段坡肩处 第四段 跳跃
72 38
15.1 27.4
2.5 3
45.44 3.071025 70.60741 0.780744
坡脚外12米 坡脚外16米
运动模式 0 9000 9000 9000 9000 9000 4597 1180.884 1324.928 2154.53 3549.196 1360.529 605.7666 0 直径1.2米 18剖面下部9# 第3段 第4段 0 2660 2660 2660 413.2746 508.1419 480.3301 0 第1段 直径1.2米 17剖面下部 第2段 第3段 跳跃 跳跃 跳跃 跳跃 跳跃 第1段 第2段 跳跃 跳跃
直径1.5米
2剖面
坡脚外32米 坡脚外36米
第6段-1 第6段-2
跳跃 跳跃
0 0
9000 9000 9000 9000 9000 9000 9000 9000
1468.547 1152.542 1519.379 717.1734 539.7648 325.1245 0 544.5 220.5
直径1.5米
0.5 0.5 0.5 0.5 0.5
0.3 0.3 0.3 0.3 0.3
0.851504 0.150598 0.740157 0.571202 0.794833
48.81234 8.633002 42.42938 32.74405 45.56367
9剖面
第3-1段 第3-2段
跳跃 跳跃
21 7.80198
12剖面4#崩塌
第4-1段 第4段-2
跳跃 跳跃
30 11.64396
0.5 0.5
0.3 0.3
0.390928 22.40989 0.146068 8.373351
13剖面5#崩塌
第8-1段 第8段-2
跳跃 跳跃
36.5 14.79368
0.5 0.5
0.3 0.3
0.485942 27.85654 0.186351 10.68252
坡度
40 48
35 0 41 52 33
第4段diyitiao 跳跃 落石的质 滚石动能 量kg E（KJ） 912.4 912.4 912.4 912.4 912.4 116.9289 900.863 710.9975 380.5656 200.2022 0 1944.299 884.6403 397.3681 0 1672.807 845.6769 0 1800 832.505 3009.39 1489.697 0 1095.12 544.5 1971.216 913.7668 0 103.68 832.505 1744.543 7070.97 8371.232 0 3353.805 1624.861 0 9000 2177.366 9剖面2#崩塌 第1段 跳跃
57.7 38 37.5 29 21.8
第3段 15剖面第二道 第4段 第5段 0 第1段 第2段 3剖面WYD2
跳跃 跳跃 跳跃
37.5 29 21.8
190 87 36.5
2.11 1 1
0.4 0.4 0.4
29.08434 19.61829 13.14844
跳跃 跳跃 跳跃 跳跃 0
55.9 28.2 0 0 58.7 30 0 0 58 26.1 0 0 47.7 55 37.6 0 0 52.7
19.28043 13.70869 20 13.6 25.86026 18.19461 15.6 11 20.9296 14.24989 4.8 13.6 19.68949 39.63996 43.13089 27.3 19 21.9968
第1段 第2段 5剖面WYD3 3-1段
跳跃 跳跃 跳跃 跳跃 0
第1段 第2段 7剖面WYD4
跳跃 跳跃 跳跃 跳跃 0
8剖面WYD5
第1段 第2段 第3段
跳跃 跳跃 跳跃 跳跃 跳跃 0
第1段 9剖面2#崩塌
跳跃
9剖面2#崩塌
第2段 第3段
跳跃 跳跃 跳跃 跳跃 0
25.6 36.7 0 0 52.4 68.2 31.8 0 0 坡度
15.8 33.4 0 0 39.6 35.6 85.3 0 0
落石撞击落石平台后最大偏离计算
运动模式
山坡入射 恢复系数 瞬间摩擦 反射角γ 角ψ ρ 系数λ (数字)
反射角γ （度数）
6剖面
第6-1段 第6段-2
跳跃 跳跃
28
0.5
0.3
0.362767 20.79558
2剖面
第2段 第滚动
58 12 52 42 55
13剖面5#崩塌
第8-1段 第8段-2
跳跃 跳跃
36.5 14.79368
0.5 0.5
0.3 0.3
0.485942 27.85654 0.186351 10.68252
每延米系统 每根锚杆设 锚杆设计锚 设计推力（KN/m） 计轴向拉力 锚杆面积 固力 （KN/m） （KN/m）
锚固长度（砂浆） 锚固长度（土体）
运动模式
山坡入射 恢复系数 瞬间摩擦 反射角γ 角ψ ρ 系数λ (数字)
反射角γ （度数）
第3段 第4段 18剖面9#崩塌 第5段 第6段 第4段-2 17剖面8#崩塌 第4-1段 第4段-2
滚动 跳跃 滚动 滚动 跳跃 跳跃 跳跃
35 14.03135 5.097647 1.823861 0.65124 33 13.04976
80 72 7 42 29 55
8977.5 2660 落石的质 滚石动能 量kg E（KJ） 0 11376 11376 11376 11376 11376 11376 9900 9000 9002 1958.297 1961.226 662.3374 · 341.482 321.9558 · 178.2232 99.98066 3042 1458.324 第1段 第2段 第3段 第4段 第5段 跳跃 跳跃 跳跃 跳跃 跳跃 38 54 42 53 36
2-2剖面
6-6剖面
2-2剖面
0 第1段 第2段 第3段 9-9 poumian 第4段 第5段 跳跃 跳跃 49 37 38 51 21 34.7 40 52 50 12.4 3 2.75 2.07 2.15 2.9 1 0.1 0.3 0.1 0.1 0.3 16.19934 17.15891 21.88114 28.08399 17.38792 16.23418
坡度
坡段垂直 崩落摩擦 瞬间摩擦 距离 系数ε 系数λ 26.5 409.5 190.4 87.2 36.5 3.11 2.15 2.11 1 1 0.3 0.4 0.4 0.3 0.3
坡段末速 度V （m/s） 16.0097 44.43772 39.47812 28.88266 20.94868
9.1 坡段末速 坡段垂直 崩落摩擦 瞬间摩擦 度V 距离 系数ε 系数λ （m/s） 13.5 120.6 58.5 0 0 35.1 167.5 60 20 57.6 42.4 80.1 0 0 1.99 3.12 1.11 0 0 3.19 2.7 1 1.99 2.36 2.56 2 0 0 0.1 0.1 0.3 0.3 0.4 0.1 0.1 0.1 0.1 0.1 0.1 0.3 0.1 0.3 7.311727 34.79852 23.06941 15.1 10.5 18.89924 38.64552 33.97516 31.62786 33.52854 34.43008 29.82852 19.1 13.2