本发明涉及一种楔形体稳定性分析技术，尤其涉及一种基于赤平极射投影理论和变形分析的楔形体稳定性分析方法。
背景技术：
自上世纪60年代，国内外学者开始进行滑坡的预报研究工作。尽管学者们在滑坡预报方面做出了不懈的努力，也取得了很多成果，但由于岩土体的离散型和复杂性，滑坡预报仍没有得到有效解决。
楔形体破坏是岩质边坡的一个主要失稳模式，在边坡失稳模式中占有重要的位置。在实际工程中，如何根据楔形体的变形来判断其稳定状态是一个重要的问题。
当抗滑力不足以抵抗下滑力时，将会发生楔形体失稳破坏。楔形体失稳破坏是个累积性的破坏过程。坡体在破坏之前会产生一定的变形，当变形达到一定的程度时，即实际应变达到最大剪切应变时，坡体会发生失稳破坏。进行楔形体稳定性分析时，只考虑楔形体的几何方向，很难准确地得到楔形体的稳定性系数。同时，楔形体是由两条或两条以上的结构面对岩体切割而成，滑移方向沿着两个面发生。由于两个滑面的产状、抗剪强度参数等方面存在差异，其对楔形体的抗滑作用也不同。如果从楔形体的位移出发就需考虑两个结构面的相互关系。
因此，建立结构面对楔形体位移的关系，对于实际工程中的监测和预报坡体变形有重要的指导意义。
技术实现要素：
本发明的目的是提供一种基于赤平极射投影和变形分析的楔形体稳定性分析方法，以解决现有技术中存在的问题。
为实现本发明目的而采用的技术方案是这样的，一种基于赤平极射投影和变形分析的楔形体稳定性分析方法，包括以下步骤：
1)运用赤平极射投影确定边坡楔形体的空间几何形态和力学特征；识别潜在不稳定楔形体w。
2)对步骤1)所识别潜在不稳定楔形体w进行受力分析。
3)对楔形体w进行位移分析；利用楔形体结构的受力特点，求解出结构面的权重参数。
4)运用步骤3)求解得到的权重参数，求解楔形体w的最大剪切位移。
5)设定楔形体w安全系数，推导临界安全位移。
6)根据监测点与结构面交线l位置关系，判断监测点位移与最大剪切应变的关系。
进一步，步骤3)中，对楔形体w进行位移分析时，选取只考虑内摩擦角
进一步，步骤3)中，利用hoek极限平衡法求解公式，建立结构面对坡体位移的权重参数。
进一步，步骤6)中，结构面交线l上位移δx1为：
式中，δx为监测点位移值。θ1为结构面交线l与开挖面间的夹角。α为开挖面与自然坡面延长面间的夹角。
进一步，步骤6)中，当楔形体w下滑存在后缘裂缝时，结构面交线l的位移为：
式中，δx’为后缘裂缝宽度，θ1为结构面交线l与开挖面间的夹角。α为开挖面与自然坡面延长面间的夹角。
本发明的技术效果是毋庸置疑的：
a.根据楔形体的受力特点，求解出楔形体的最大剪切位移，避免了传统方法中对楔形体位移求解的不准确性；
b.方法概念明确，工程应用简便。
附图说明
图1为楔形体变形示意图；
图2为楔形体分析要素图；
图3为水压力分布图；
图4为赤平极射投影图；
图5为垂直开挖示意图；
图6为倾斜开挖示意图；
图7为楔形体后缘裂缝示意图；
图8为方法流程图。
图中：楔形体w、结构面交线l、结构面a1、交线ⅰ101、交线ⅲ102、结构面b2、交线ⅱ201、交线ⅳ202、自然坡面3、开挖面4。
具体实施方式
下面结合实施例对本发明作进一步说明，但不应该理解为本发明上述主题范围仅限于下述实施例。在不脱离本发明上述技术思想的情况下，根据本领域普通技术知识和惯用手段，做出各种替换和变更，均应包括在本发明的保护范围内。
实施例1：
赤平极射投影(stereographicprojection)是一种极点投影法,可将空间线、面的三维相对位置关系表达在二维空间。结构面与投影球相交得到结构面的球面投影。通过一个球极向结构面球面投影的圆周各点发出射线,该射线和赤道平面的交点轨迹构成结构面的大圆赤平投影。
参见图8，本实施例公开一种基于赤平极射投影和变形分析的楔形体稳定性分析方法，包括以下步骤：
1)运用赤平极射投影确定边坡楔形体的空间几何形态和力学特征。在多结构面的岩体中寻找潜在滑移的楔形体，并初步判别楔形体w的滑动类型，为进一步的稳定分析提供必要的参数。
参见图1和图2，在本实施例中，岩体边坡受结构面a1和结构面b2的切割,形成一个孤立的潜在不稳定楔形体w。开挖面4与结构面a1的交线为交线ⅰ101，与结构面b2的交线为交线ⅱ201。自然坡面3与结构面a1的交线为交线ⅲ102，与结构面b2的交线为交线ⅳ202。结构面a1和结构面b2的交线为结构面交线l。楔形体w本身不透水。楔形体水压力分布如图3所示。水从楔形体w顶端沿交线ⅲ102和交线ⅳ202进入，并从坡面上的交线ⅰ101和交线ⅱ201处渗出。最大压力出现在结构面交线l上，且在交线ⅰ101、交线ⅱ201、交线ⅲ102和交线ⅳ202上压力为零。
2)对潜在不稳定楔形体w进行受力分析。根据e.hoek图解法，将结构面a1、结构面b2、自然坡面3和开挖面4绘制在赤平极射投影图上。参见图4，4个圆弧有5个交点，分别代表了5条线，各线之间的夹角可在图中测出。在赤平极射投影图上，外圆为赤道大圆，代表赤道平面(既水平面)，其上、下、左、右分别代表北(n)、南(s)、西(w)、东(e)方位，并按360°方位角分度。各结构面投影在赤道大圆上反映为一条圆弧线：弧线与赤道大圆两侧交点的连线为大圆的一条直径，代表了结构面的走向；弧线弯曲方向代表结构面的倾向；弧线中点距大圆圆心的距离代表了结构面的倾角，弧线中点距圆心越近、倾角越大(通过圆心，倾角为90°)，弧线中点距圆心越远、倾角越小(通过大圆，倾角为0°)。图中，na为结构面a1的极点，即结构面a1法线的投影。nb为结构面b2的极点，即结构面b2法线的投影。根据测得的角度，求出楔形体w的几何形状参数：
式中，θ13为交线ⅰ101和交线ⅲ102之间的夹角，°；θ24为交线ⅱ201和交线ⅳ202之间的夹角，°；θ35为结构面交线l和交线ⅲ102之间的夹角，°；θ45为结构面交线l和交线ⅳ202之间的夹角，°；θ1nb为交线ⅰ101与结构面b2法线之间的夹角，°；θ2na为交线ⅱ201与结构面a1法线之间的夹角，°；θna,nb为结构面a1法线与结构面b2法线之间的夹角，°；ψa为结构面a1的倾角，°；ψb为结构面b2的倾角，°；ψ5为结构面交线l的倾角，°；
则楔形体w稳定性系数k为：
k＝hcr/h(5)
式中，h为楔形体w的高度，hcr为楔形体的临界高度。其中，
ψ5＝tanψacos(αa-α5)＝tanψbcos(αb-α5)(8)
式中，αa为结构面a1的倾向，°；αb为结构面b2的倾向，°；α5为结构面交线l的倾向，°；ca为结构面a1的黏聚力，kpa；
3)对楔形体w进行位移分析。利用楔形体结构的受力特点，求解出结构面的权重参数。
在楔形体位移求解中，考虑到结构面a1和结构面b2都对楔形体w位移产生滑移作用，引用hoek提出的极限平衡法求解公式，建立结构面a和结构面b的对坡体位移的权重参数。
当只考虑内摩擦角作用时，楔形体安全系数计算公式为：
根据式(9)建立结构面a1和结构面b2的权重参数为：
式(10)和式(11)中，ma,mb分别为只考虑内摩擦角时结构面a和结构面b的权重参数。
建立楔形体w位移公式为
εmax＝ma·εmax,a+mb·εmax,b(12)
式(12)中，εmax,a为结构面a1的最大剪切位移；εmax,b为结构面b2的最大剪切位移。
当考虑内摩擦角和黏聚力作用时，楔形体安全系数计算公式为：
根据式(13)建立结构面a和结构面b的权重参数为
式(14)和式(15)中，na,nb分别为考虑内摩擦角和黏聚力时结构面a和结构面b的权重参数。
建立楔形体位移公式为
εmax＝na·εmax,a+nb·εmax,b(16)
式(16)中，εmax,a为结构面a的最大剪切位移；εmax,b为结构面b的最大剪切位移。
当考虑内摩擦角、黏聚力和水头作用时，楔形体安全系数计算公式为：
根据式(17)建立结构面a和结构面b的权重参数为
式(18)和式(19)中，pa,pb分别为考虑内摩擦角、黏聚力和地下水作用时结构面a和结构面b的权重参数。
建立楔形体位移公式为
εmax＝na·εmax,a+nb·εmax,b(20)
式(20)中，εmax,a为结构面a的最大剪切位移；εmax,b为结构面b的最大剪切位移。
4)参考试验资料及相关规范，运用步骤3)求解得到的权重参数，求解楔形体w的最大剪切位移。
5)设定楔形体w安全系数，推导临界安全位移。
6)根据监测点与结构面交线l位置关系，判断监测点位移与最大剪切应变的关系。
结构面交线l的倾角为θ2，结构面交线l与开挖面4夹角为θ1，开挖面4与自然坡面3延长面间夹角为α。
参见图5，当坡面垂直开挖时，由监测点位移值δx可得结构面交线上位移：
参见图6，当坡面倾斜开挖时，由监测点位移值δx可得结构面交线上位移：
参见图7，当坡面倾斜开挖时，若楔形体下滑存在后缘裂缝，且测得后缘裂缝宽度为δx’，则可得结构面交线位移为：
7)应用到工程实践。
值得说明的是，本实施例利用赤平极射投影法，对楔形体进行受力分析。由于楔形体结构面抗剪强度对楔形体稳定性影响程度不同，运用其受力特点，推导出楔形体结构面权重参数。分别列出只考虑内摩擦角作用，考虑内摩擦角和黏聚力作用，综合考虑内摩擦角、黏聚力、水头作用时的权重参数。参考试验资料及相关规范，列出单结构面坡体的最大剪切应变，运用求解得到的权重参数求解出楔形体的最大剪切应变。设定楔形体安全系数推导出临界位移值，根据监测点与结构面交线位置关系，判断监测点位移与最大剪切应变的关系。在实际工程中，对于结构面多而复杂的岩质边坡，如果发生楔形体滑移失稳的岩体，需要选取两组结构面，其稳定判断方法同上。