发明名称 一种外倾顺层节理岩体被动侧向岩石压力的计算方法
摘要 本发明是一种顺层节理岩体被动侧向岩石压力的计算方法,属于岩土力学中岩土压力分析技术领域。本发明以外倾顺层节理岩体为研究对象,将顺层节理岩体按节理面离散为有限数量的刚性岩块,基于潘家铮最大值原理,以被动侧向岩石压力为目标函数,将节理面的剪力和法向力作为未知量,构建满足破坏岩体的平衡方程、节理面的屈服条件、岩体与支护结构接触面的屈服条件,建立求解支护结构上被动侧向岩石压力的数学规划模型,采用有效集合算法对线性数学规划模型进行求解,求解得到外倾顺层节理岩体被动侧向岩石压力的最大极限值。发明方法概念明确、计算精度高、工程应用简便,可将其应用于节理岩质边坡挡土墙、抗滑桩、桩锚支护等结构的设计中。
申请公布号 CN106547942A 申请公布日期 2017.03.29
申请号 CN201610850797.3 申请日期 2016.09.26
申请人 昆明理工大学 发明人 李泽;周宇;薛龙
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 代理人
主权项 一种外倾顺层节理岩体被动侧向岩石压力的计算方法,其特征在于:以外倾顺层节理岩体为研究对象,将顺层节理岩体按节理面离散为有限数量的刚性岩块,基于潘家铮最大值原理,以被动侧向岩石压力为目标函数,将节理面的剪力和法向力作为未知量,构建满足破坏岩体的平衡方程、节理面的屈服条件、岩体与支护结构接触面的屈服条件,建立求解支护结构上被动侧向岩石压力的数学规划模型,采用有效集合算法对线性数学规划模型进行求解,求解得到外倾顺层节理岩体被动侧向岩石压力的最大极限值;具体步骤如下:(一)拟定外倾顺层节理岩体的计算参数根据外倾顺层节理岩体的实际情况,拟定其计算参数,包括:几何参数、节理地质参数,岩体、荷载参数信息、节理材料的物理力学参数,节理材料的物理力学参数包括容重、凝聚力、摩擦角等;(二)将外倾顺层节理岩体沿节理面离散为刚性岩块沿各层节理面将破坏岩体离散为刚性岩块,并以相邻岩块之间节理面的法向力、剪切力以及支护结构与岩体交界面的法向力、切向力为未知量,同时要求岩块满足力的平衡条件、节理面满足摩尔库伦屈服条件;(三)建立求解极限被动侧向岩石压力的线性数学规划模型(1)岩体被动破坏的受力分析外倾顺层节理岩体发生被动破坏时,ABC为发生被动破坏的岩体;ABC形心作用有自重G;AB为支护结构与岩体交界面,其上作用有水平被动侧向岩石压力E<sub>p</sub>和剪切力V<sub>p</sub>,规定E<sub>p</sub>以受压为正,V<sub>p</sub>以向上为正,两者的合力为P<sub>p</sub>;BC为发生破坏的节理面,其上作用有法向力N<sub>R</sub>和剪切力V<sub>R</sub>,规定N<sub>R</sub>以受压为正,V<sub>R</sub>以指向左下为正;破坏的岩体ABC沿节理面离散成岩块以后,取岩块i作为研究对象;节理面i作用有法向力<img file="FDA0001120329640000011.GIF" wi="58" he="63" />和切向力<img file="FDA0001120329640000012.GIF" wi="75" he="63" />支护结构与岩体交界面i作用有被动侧向岩土压力<img file="FDA0001120329640000013.GIF" wi="57" he="71" />和剪力<img file="FDA0001120329640000014.GIF" wi="75" he="70" />岩块i形心作用有自重G<sub>i</sub>,在这些力的作用下,岩块保持极限状态平衡;(2)岩体被动破坏的极限状态①岩块的平衡方程:<maths num="0001"><math><![CDATA[<mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msubsup><mi>E</mi><mi>p</mi><mi>i</mi></msubsup><mo>+</mo><mrow><mo>(</mo><msubsup><mi>N</mi><mi>R</mi><mi>i</mi></msubsup><mo>-</mo><msubsup><mi>N</mi><mi>R</mi><mrow><mo>&prime;</mo><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msubsup><mo>)</mo></mrow><mi>s</mi><mi>i</mi><mi>n</mi><mi>&theta;</mi><mo>-</mo><mrow><mo>(</mo><msubsup><mi>V</mi><mi>R</mi><mi>i</mi></msubsup><mo>-</mo><msubsup><mi>V</mi><mi>R</mi><mrow><mo>&prime;</mo><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msubsup><mo>)</mo></mrow><mi>c</mi><mi>o</mi><mi>s</mi><mi>&theta;</mi><mo>=</mo><mn>0</mn><mo>,</mo></mrow></mtd><mtd><mrow><mo>(</mo><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mo>...</mo><mo>,</mo><msub><mi>n</mi><mi>b</mi></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mrow><msubsup><mi>V</mi><mi>p</mi><mi>i</mi></msubsup><mo>+</mo><mrow><mo>(</mo><msubsup><mi>N</mi><mi>R</mi><mi>i</mi></msubsup><mo>-</mo><msubsup><mi>N</mi><mi>R</mi><mrow><mo>&prime;</mo><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msubsup><mo>)</mo></mrow><mi>c</mi><mi>o</mi><mi>s</mi><mi>&theta;</mi><mo>-</mo><mrow><mo>(</mo><msubsup><mi>V</mi><mi>R</mi><mi>i</mi></msubsup><mo>-</mo><msubsup><mi>V</mi><mi>R</mi><mrow><mo>&prime;</mo><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msubsup><mo>)</mo></mrow><mi>s</mi><mi>i</mi><mi>n</mi><mi>&theta;</mi><mo>-</mo><msub><mi>G</mi><mi>i</mi></msub><mo>=</mo><mn>0</mn><mo>,</mo></mrow></mtd><mtd><mrow><mo>(</mo><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mo>...</mo><mo>,</mo><msub><mi>n</mi><mi>b</mi></msub><mo>)</mo></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0001120329640000015.GIF" wi="1286" he="167" /></maths>上式中:(i=1,...,n<sub>b</sub>),n<sub>b</sub>为岩块数量,G<sub>i</sub>为岩块i的自重;<img file="FDA0001120329640000016.GIF" wi="62" he="61" />为节理面i的法向力,<img file="FDA0001120329640000017.GIF" wi="49" he="63" />为节理面i的切向力;<img file="FDA0001120329640000018.GIF" wi="93" he="63" />为节理面i‑1的法向力,<img file="FDA0001120329640000019.GIF" wi="86" he="64" />为节理面i‑1的切向力;<img file="FDA00011203296400000110.GIF" wi="61" he="71" />为支护结构与岩体交界面i的被动侧向岩土压力,<img file="FDA00011203296400000111.GIF" wi="48" he="68" />为支护结构与岩体交界面i的切向力,θ为顺层节理的倾角,逆时针方向为正;②节理面屈服条件当支护结构后的节理岩体处于被动破坏的临界状态时,破坏节理面应满足Mohr‑Coulomb屈服条件,节理面的屈服条件为:<img file="FDA0001120329640000021.GIF" wi="789" he="70" />其中:(i=1,...,n<sub>j</sub>),n<sub>j</sub>为节理面的数量,c为节理面的凝聚力,<img file="FDA0001120329640000022.GIF" wi="38" he="46" />为节理面的内摩擦角,l<sub>i</sub>为节理面i的长度,并规定N<sub>R</sub>以受压为正,V<sub>R</sub>以指向左下为正;③支护结构与岩体交界面作用力的约束条件支护结构与岩体交界面作用力的被动侧向岩石压力<img file="FDA0001120329640000023.GIF" wi="59" he="70" />与切向力<img file="FDA0001120329640000024.GIF" wi="51" he="71" />需满足以下关系:<img file="FDA0001120329640000025.GIF" wi="654" he="71" />其中:(i=1,...,n<sub>E</sub>),n<sub>E</sub>为支护结构与岩体交界面的数量,<img file="FDA0001120329640000026.GIF" wi="57" he="70" />为支护结构与岩体交界面i的被动侧向岩土压力,<img file="FDA0001120329640000027.GIF" wi="46" he="70" />为支护结构与岩体交界面i的切向力,δ为支护结构与岩体交界面之间的内摩擦角;④目标函数本发明的目的是求解极限被动侧向岩石压力E<sub>p</sub>,因此设E<sub>p</sub>为目标函数:<img file="FDA0001120329640000028.GIF" wi="733" he="71" />其中:(i=1,...,n<sub>E</sub>),n<sub>E</sub>为支护结构与岩体交界面的数量,<img file="FDA0001120329640000029.GIF" wi="55" he="71" />为支护结构与岩体交界面i的被动侧向岩土压力;因此,求解极限被动侧向岩石压力E<sub>p</sub>最大化的目标函数为:<img file="FDA00011203296400000210.GIF" wi="950" he="71" />(3)求解极限被动侧向岩石压力的线性数学规划模型根据以上推导得到的目标函数和约束条件,求解极限被动侧向岩石压力的数学规划模型具体表达式为:<img file="FDA00011203296400000211.GIF" wi="1525" he="437" />
地址 650093 云南省昆明市五华区学府路253号