CN104484503B - 考虑作用点位置的基坑柔性支护主动土压力计算方法 - Google Patents
考虑作用点位置的基坑柔性支护主动土压力计算方法 Download PDFInfo
- Publication number
- CN104484503B CN104484503B CN201410672267.5A CN201410672267A CN104484503B CN 104484503 B CN104484503 B CN 104484503B CN 201410672267 A CN201410672267 A CN 201410672267A CN 104484503 B CN104484503 B CN 104484503B
- Authority
- CN
- China
- Prior art keywords
- msub
- mrow
- theta
- msup
- slip
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Abstract
本发明涉及一种考虑作用点位置的基坑柔性支护主动土压力计算方法,该方法主要是针对库仑土压力和朗肯土压力计算方法的不足,提供了一种考虑作用点位置的主动土压力精确求解方法,包括以下步骤:确定精度控制量eps、基坑几何要素和土体物理力学参数;确定主动土压力合力作用点位置系数na;土体滑裂面用曲线y=s(x)表示;确定函数Φ(x0,s0)的最小值Φmin,判断Φmin是否小于等于eps;最后计算得出滑裂面方程和主动土压力合力。本发明所提供的一种符合实际的基坑柔性支护结构的主动土压力计算方法,可靠性高,有利于基坑工程的合理设计和科学管理。
Description
技术领域
本发明属于岩土工程技术领域,涉及一种基坑支护结构主动土压力的计算方法,特别涉及一种考虑作用点位置的基坑柔性支护结构主动土压力计算方法。
背景技术
支护结构在基坑开挖和使用过程中的受力变形是一个极为复杂的问题,影响因素很多,正确计算基坑工程中侧向土压力的分布和大小是合理设计基坑工程的前提。从广义上讲,土压力是土作用在工程结构上的或作用在被土体所包围的结构物表面上的压力或那些压力的合力。在基坑工程的设计与施工中,土压力是作用在支护结构上的主要荷载,在进行支护结构强度计算和土体整体稳定性验算以及变形量计算时,首先要确定的是作用在结构体上的总土压力分布,土压力是土与挡土结构之间相互作用的结果,它主要取决于支护桩(墙)体的变位方式、方向、大小。基坑工程往往处于城市密集区,周围环境要求极高,如果不能正确计算土压力,就很难较精确地预测墙土体的变形,从而无法进行环境变形控制设计。合理地进行基坑工程设计,首先必须提出简单实用而尽可能合理的土压力计算模式。目前仍然广泛沿用极限平衡理论为基础的土压力理论,如朗肯土压力理论和库仑土压力理论。在无超载情况下,上述两种理论沿墙背的土压力分布自上至下由数值零开始线性分布,合力作用点位于墙体的下三分点处。考虑到柔性支护结构和土体接触面的共同协调变形,导致土压力的分布不一定是三角形分布,甚至由于支护结构刚度的不均匀导致土压力的非线性分布。比如当柔性支护结构顶部的位移较大时,导土压力的分布图形大致为矩形,因而,土压力的作用点大致在中点,与其相对应的墙后的不稳定土体的滑裂面也不是直线。因此,上述的主动土压力理论并不适用于柔性支护结构领域。在基坑工程中,柔性支护结构可以通过调节支护方式、构件位置来调整土压力作用位置,也可能因为施工方法、步骤和施工位置的不同,也会使作用点位置发生变化。因此,研究一种符合实际的基坑柔性支护结构的主动土压力计算方法,是基坑工程合理设计并科学管理的基础。
发明内容
有鉴于此,本发明的目的在于提供考虑作用点位置的基坑柔性支护主动土压力计算方法,该计算方法是考虑作用点位置的主动土压力精确求解方法,适用于交通、水利、市政、建筑等部门中基坑支护方案的设计。
为达到上述目的,本发明提供如下技术方案:
考虑作用点位置的基坑柔性支护主动土压力计算方法,该方法包括以下步骤:
S1:确定精度控制量eps、基坑几何要素和土体物理力学参数;
所述eps是一个接近于0的数值;所述基坑几何要素包括基坑深度H;所述土体物理力学参数包括基坑土体的重度γ,内摩擦角φ,内聚力c,通过取样和实验手段确定它们的数值;
S2:确定主动土压力合力作用点位置系数na;
所述na是作用点位置至坑底的距离与基坑深度的比值,通过进行调查研究,根据柔性支护形式和实际工程工况来合理确定na;
S3:初定滑裂面中心坐标(x0,s0),由此土体滑裂面可用对数螺旋线y=s(x)表示,Pa为作用在支护结构的主动土压力,σ(x)、τ(x)分别为作用在滑裂面上的法向、切向应力;
S4:求出函数Φ(x0,s0)的表达式,搜索滑裂面真实的中心坐标(x0,s0),寻找函数Φ(x0,s0)的最小值Φmin,若Φmin≤eps,则跳转至S5,否则终止计算,表明在此作用点处,支护结构达不到平衡;
S5:确定滑裂面方程;
根据S4搜索的滑裂面真实中心坐标(x0,s0),即可得到滑裂面方程:
S6:计算主动土压力合力;主动土压力合力Pa的计算公式如下;
进一步,步骤S4中所述函数Φ(x0,s0)值的求解步骤如下:
S41:根据初定的滑裂面中心坐标(x0,s0),计算滑裂面O点处角坐标θ0:
S42:求解以下方程,计算滑裂面A点处角坐标θA
其中,
S43:计算滑裂面处的正应力
其中:
S44:求函数Φ(x0,s0);
其中,
由此Φ(x0,s0)的值可以通过数值积分的方法求得;
S45:搜索滑裂面真实中心坐标,寻找函数Φ(x0,s0)的最小值Φmin。
进一步,所述eps取值为1×10-5。
进一步,所述初定滑裂面中心坐标(x0,s0),x0=-2H,s0=2H。
进一步,所述Φ(x0,s0)的最小值Φmin通过非线性函数优化方法求解,采用matlab软件中的 fminsearch来进行搜索求解。
本发明的有益效果在于:本发明所提供的考虑作用点位置的基坑柔性支护主动土压力计算方法,可适用于柔性支护结构领域,在基坑工程中,柔性支护结构可以通过调节支护方式、构件位置来调整土压力作用位置,也可能因为施工方法、步骤和施工位置的不同,也会使作用点位置发生变化。因此,本发明提供的一种考虑作用点位置的基坑柔性支护结构主动土压力计算方法,该方法可靠性高,有利于基坑工程的合理设计和科学管理。
附图说明
为了使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明作进一步的详细描述,其中:
图1为本发明所述方法的流程图;
图2为基坑柔性支护主动土压力计算模型;
图3为实施例1中函数Φ随作用点位置系数na的变化曲线;
图4为实施例1中主动土压力随作用点位置系数na的变化曲线;
图5为实施例1中滑裂面曲线随作用点位置系数na的变化。
具体实施方式
下面将结合附图,对本发明的优选实施例进行详细的描述。
本发明所提供的考虑作用点位置的基坑柔性支护主动土压力计算方法,该方法主要是针对库仑土压力和朗肯土压力计算方法的不足,设计的一种考虑作用点位置的主动土压力精确求解方法,包括以下步骤:
(1)确定精度控制量eps,基坑几何要素和土体物理力学参数。
(2)确定主动土压力合力作用点位置系数na。
(3)初定滑裂面中心坐标(x0,s0),土体滑裂面用对数螺旋线y=s(x)表示,Pa为作用在支护结构的主动土压力,σ(x)、τ(x)分别为作用在滑裂面上的法向、切向应力。
(4)求出函数Φ(x0,s0)的表达式,搜索滑裂面真实的中心坐标,寻找函数Φ(x0,s0)的最小值Φmin,若Φmin≤eps,则跳转至S5,否则终止计算,表明在此作用点处,支护结构达不到平衡;
(5)确定滑裂面方程。
(6)计算主动土压力合力。
步骤(1)中eps是一个接近于0的数值,可取1×10-5;基坑几何要素包括基坑深度H,通过取样和实验手段确定基坑土体的重度γ,内摩擦角内聚力c。
步骤(2)中主动土压力合力作用点位置系数na的值是作用点位置至坑底的距离与基坑深度的比值。应进行调查研究,根据柔性支护形式和实际工程工况来合理确定。
步骤(3)中极限平衡状态下土体滑裂面为对数螺旋线,用曲线y=s(x)表示,Pa为作用在支护结构的主动土压力,σ(x)、τ(x)分别为作用在滑裂面上的法向、切向应力,见图2。
步骤(4)中通过非线性函数优化方法求解函数Φ(x0,s0)最小值Φmin,这里采用matlab软件中的fminsearch来进行搜索求解。
函数Φ(x0,s0)值的求解步骤如下:
1)初定滑裂面中心坐标(x0,s0),如x0=-2H,s0=2H;
2)计算滑裂面O点处角坐标θ0:
3)求解以下方程,计算滑裂面A点处角坐标θA
4)计算滑裂面处的正应力
其中:
5)求函数Φ(x0,s0)的值;
其中
由此Φ(x0,s0)的值可以通过数值积分的方法求得。
步骤(5)中根据步骤(4)得到的滑裂面中心坐标(x0,s0),可以得到滑裂面方程:
步骤(6)中主动土压力合力Pa的计算公式如下;
步骤(4)、(5)、(6)的理论推导如下:
取支护结构后处于极限平衡状态的滑楔体OAB作为研究对象,根据力的平衡方程,
由∑X=0得:
由∑Y=0得:
由∑MO=0得:
上三式中s′=ds/dx,xA为点A的X坐标,na为主动土压力作用点位置系数,其值为作用点位置至坑底的距离与基坑深度的比值。设在滑面上法向应力σ(x)和切向应力τ(x)服从Mohr-Coulomb破坏准则,即:
其中,
考察(1)~(4)式,主动土压力Pa为自变量函数σ(x)和滑面s(x)的泛函极值问题。具体如下:
由式(3)得泛函:
其中
由式(1)得约束条件:
其中
由式(2)得约束条件:
其中
上述为边界待定的条件变分极值中的等周问题。滑裂面曲线的起始点为O点,坐标为 xo=0,yo=0,终点为土体表面上A点坐标(xA,H)待定。
在约束条件下的变分法,用拉格朗日乘子法构造如下的泛函J*,使上述条件极值问题转化为无约束的泛函极值问题:
其中
F为辅助函数,λ1、λ2为拉格朗日乘子。依据等周问题极值存在的必要条件,滑面方程y=s(x)及沿滑面分布的法向应力σ(x)必须满足欧拉微分方程、边界条件及可动边界处的横截条件:
助函数F的Euler微分方程:
(1)积分约束方程:同式(6)、式(7)
(2)边界条件:
固定边界条件:
可动边界条件:
(3)可动边界处的横截条件:
由式(9),可以得到:
引入坐标变换:
设新坐标原点在原坐标中为(x0,s0)。则
令
则微分方程式(14)可化为齐次方程:
分离变量得通解为:
对变换后的坐标换成极坐标:
式(17)变为:
c0、c1为任意积分常数。
由固定边界条件:s(0)=0,原坐标原点O在新坐标中为u0=x0,v0=-s0,O点在新坐标的极坐标为(ro,θo),则滑裂面方程为:
其中,
滑裂面为对数螺旋面。
由式(10)得:
将上式写成新坐标下的极坐标形式:
上式微分方程通解为:
z为积分常数,θ1为任意角度,可取θ1=0
由可动边界处的横截条件式可得滑面A点处正应力:
代入(23)
上述变分问题只包含两个未知常数x0,s0,可以由两个约束条件方程式(6)和式(7)来求出,等同于求下式函数Φ的零值问题:
上式的解可通过求解函数Φ的极小值且极小值为0得到,而当Φ的极小值不为零时,说明所给定的土压力作用点不能使土体平衡,是不合适的。
实施例1:
设一垂直基坑深10米,基坑周围土γ=18kN/m3,内聚力c=10kPa,内摩擦角φ=20°,图3 为函数Φ随作用点位置系数na的变化曲线。计算表明,当作用点位置系数存在上下界限值(下限值nad=0.2677,上限值nau=0.4906),在此范围内,Φ=0值,也即是说当支护结构所提供的支护反力的合力作用点位于离坑底2.677~4.906范围内时,基坑可以达到平衡稳定,否则无论提供多大的支护反力,当土体处于主动状态时,基坑必定失稳。主动土压力及滑裂面随作用点位置系数na的变化曲线见图4、图5所示,在作用点位置系数下限处,主动土压力最小。滑裂面为平面;随着作用点位置的上移,主动土压力呈非线性增长,相应滑裂面为对数螺旋面。表1为不同作用点位置的各主动压力值。
表1 不同作用点位置下的主动土压力值
实施例2:
采用一砂土垂直基坑作为算例,基坑深10米,基坑周围土γ=18kN/m3,内聚力c=0,内摩擦角φ=20°,表2为砂土内摩擦角不同情况下,作用点位置系数上下限处的计算结果。从表中可以看出,对于砂土,在不同内摩擦角下,作用点位置系数下限系数均为nad=1/3,所计算压力与按库仑理论计算的结果完全一致。表2中Δ为作用点位置系数上限值处土压力比下限值处的增大比例,可见作用点位置系数上限值随内摩擦角的增大而增大,其相应的土压力值也随之增加,且摩擦角越大,上限处压力增大比例也越大。
表2 砂土基坑主动土压力算例
最后说明的是,以上优选实施例仅用以说明本发明的技术方案而非限制,尽管通过上述优选实施例已经对本发明进行了详细的描述,但本领域技术人员应当理解,可以在形式上和细节上对其作出各种各样的改变,而不偏离本发明权利要求书所限定的范围。
Claims (4)
1.考虑作用点位置的基坑柔性支护主动土压力计算方法,其特征在于:该方法包括以下步骤:
S1:确定精度控制量eps、基坑几何要素和土体物理力学参数;
所述eps是一个接近于0的数值;所述基坑几何要素包括基坑深度H;所述土体物理力学参数包括基坑土体的重度γ,内摩擦角φ,内聚力c;
S2:确定主动土压力合力作用点位置系数na,所述na是作用点位置至坑底的距离与基坑深度的比值;
S3:初定滑裂面中心坐标(x0,s0),土体滑裂面用对数螺旋线y=s(x)表示,Pa为作用在支护结构的主动土压力,σ(x)、τ(x)分别为作用在滑裂面上的法向、切向应力;
S4:求出函数Φ(x0,s0)的表达式,搜索滑裂面真实的中心坐标,寻找函数Φ(x0,s0)的最小值Φmin,若Φmin≤eps,则跳转至S5,否则终止计算,表明在此作用点处,支护结构达不到平衡;所述函数Φ(x0,s0)的值求解步骤如下:
S41:根据初定的滑裂面中心坐标(x0,s0),计算滑裂面O点处角坐标θ0:
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&theta;</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<mi>arctan</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>0</mn>
</msub>
<mo>/</mo>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
<mo>&le;</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&theta;</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<mo>-</mo>
<mi>&pi;</mi>
<mo>-</mo>
<mi>arctan</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>0</mn>
</msub>
<mo>/</mo>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
<mo>></mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
S42:求解以下方程,计算滑裂面A点处角坐标θA
<mrow>
<mi>H</mi>
<mo>-</mo>
<msub>
<mi>s</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<msub>
<mi>r</mi>
<mn>0</mn>
</msub>
<msup>
<mi>e</mi>
<mrow>
<msub>
<mi>n</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
</mrow>
</msup>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&theta;</mi>
<mo>=</mo>
<mn>0</mn>
<mo>,</mo>
</mrow>
其中,n1=tanφ;
S43:计算滑裂面处的正应力
<mrow>
<mi>&sigma;</mi>
<mo>=</mo>
<msup>
<mi>ze</mi>
<mrow>
<mn>2</mn>
<msub>
<mi>n</mi>
<mn>1</mn>
</msub>
<mi>&theta;</mi>
</mrow>
</msup>
<mo>-</mo>
<mfrac>
<mrow>
<msub>
<mi>&gamma;r</mi>
<mn>0</mn>
</msub>
<msup>
<mi>e</mi>
<mrow>
<msub>
<mi>n</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&theta;</mi>
<mo>-</mo>
<mn>3</mn>
<msub>
<mi>n</mi>
<mn>1</mn>
</msub>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mn>9</mn>
<msubsup>
<mi>n</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
</mrow>
</mfrac>
<mo>-</mo>
<mfrac>
<mi>c</mi>
<msub>
<mi>n</mi>
<mn>1</mn>
</msub>
</mfrac>
</mrow>
其中:
S44:求函数Φ(x0,s0);
<mrow>
<mi>&Phi;</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
<mo>,</mo>
<msub>
<mi>s</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mo>&Integral;</mo>
<msub>
<mi>&theta;</mi>
<mn>0</mn>
</msub>
<msub>
<mi>&theta;</mi>
<mi>A</mi>
</msub>
</msubsup>
<mo>(</mo>
<mrow>
<mfrac>
<msub>
<mi>F</mi>
<mn>0</mn>
</msub>
<mrow>
<msub>
<mi>n</mi>
<mi>a</mi>
</msub>
<mi>H</mi>
</mrow>
</mfrac>
<mo>+</mo>
<msub>
<mi>F</mi>
<mn>1</mn>
</msub>
</mrow>
<mo>)</mo>
<mi>d</mi>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mo>&Integral;</mo>
<msub>
<mi>&theta;</mi>
<mn>0</mn>
</msub>
<msub>
<mi>&theta;</mi>
<mi>A</mi>
</msub>
</msubsup>
<msub>
<mi>F</mi>
<mn>2</mn>
</msub>
<mi>d</mi>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
其中,
<mrow>
<msub>
<mi>F</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>n</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msup>
<mi>s</mi>
<mo>&prime;</mo>
</msup>
<mo>)</mo>
</mrow>
<mi>&sigma;</mi>
<mo>+</mo>
<mi>c</mi>
<mo>;</mo>
<msub>
<mi>F</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>n</mi>
<mn>1</mn>
</msub>
<msup>
<mi>s</mi>
<mo>&prime;</mo>
</msup>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mi>&sigma;</mi>
<mo>+</mo>
<msup>
<mi>cs</mi>
<mo>&prime;</mo>
</msup>
<mo>+</mo>
<mi>&gamma;</mi>
<mi>s</mi>
<mo>-</mo>
<mi>&gamma;</mi>
<mi>H</mi>
<mo>;</mo>
</mrow>
<mrow>
<mi>x</mi>
<mo>=</mo>
<msub>
<mi>r</mi>
<mn>0</mn>
</msub>
<msup>
<mi>e</mi>
<mrow>
<msub>
<mi>n</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
</mrow>
</msup>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&theta;</mi>
<mo>+</mo>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
<mo>;</mo>
<mi>s</mi>
<mo>=</mo>
<msub>
<mi>r</mi>
<mn>0</mn>
</msub>
<msup>
<mi>e</mi>
<mrow>
<msub>
<mi>n</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
</mrow>
</msup>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&theta;</mi>
<mo>+</mo>
<msub>
<mi>s</mi>
<mn>0</mn>
</msub>
<mo>;</mo>
</mrow>
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msup>
<mi>s</mi>
<mo>&prime;</mo>
</msup>
<mo>=</mo>
<mo>-</mo>
<mi>cot</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>&theta;</mi>
<mo>+</mo>
<mi>&phi;</mi>
</mrow>
<mo>)</mo>
</mrow>
<mo>;</mo>
<mi>d</mi>
<mi>x</mi>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>r</mi>
<mn>0</mn>
</msub>
<msup>
<mi>e</mi>
<mrow>
<msub>
<mi>n</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>&theta;</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<mi>&theta;</mi>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>n</mi>
<mn>1</mn>
</msub>
<mi>cos</mi>
<mi>&theta;</mi>
<mo>+</mo>
<mi>sin</mi>
<mi>&theta;</mi>
</mrow>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>&theta;</mi>
<mo>;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>r</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<msqrt>
<mrow>
<msubsup>
<mi>x</mi>
<mn>0</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>s</mi>
<mn>0</mn>
<mn>2</mn>
</msubsup>
</mrow>
</msqrt>
<mo>;</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
Φ(x0,s0)的值可以通过数值积分的方法求得;
S45:搜索滑裂面真实中心坐标,寻找函数Φ(x0,s0)的最小值Φmin;
S5:确定滑裂面方程;
根据S4搜索得到的滑裂面中心坐标(x0,s0),即可得到滑裂面方程:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>x</mi>
<mo>=</mo>
<msub>
<mi>r</mi>
<mn>0</mn>
</msub>
<msup>
<mi>e</mi>
<mrow>
<msub>
<mi>n</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
</mrow>
</msup>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&theta;</mi>
<mo>+</mo>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>s</mi>
<mo>=</mo>
<msub>
<mi>r</mi>
<mn>0</mn>
</msub>
<msup>
<mi>e</mi>
<mrow>
<msub>
<mi>n</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>&theta;</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<mi>&theta;</mi>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</msup>
<mi>sin</mi>
<mi>&theta;</mi>
<mo>+</mo>
<msub>
<mi>s</mi>
<mn>0</mn>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mn>0</mn>
</msub>
<mo>&le;</mo>
<mi>&theta;</mi>
<mo>&le;</mo>
<msub>
<mi>&theta;</mi>
<mi>A</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
S6:计算主动土压力合力;主动土压力合力Pa的计算公式如下;
<mrow>
<msub>
<mi>P</mi>
<mi>a</mi>
</msub>
<mo>=</mo>
<msubsup>
<mo>&Integral;</mo>
<msub>
<mi>&theta;</mi>
<mn>0</mn>
</msub>
<msub>
<mi>&theta;</mi>
<mi>A</mi>
</msub>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>n</mi>
<mn>1</mn>
</msub>
<mi>&sigma;</mi>
<mo>-</mo>
<msup>
<mi>s</mi>
<mo>&prime;</mo>
</msup>
<mi>&sigma;</mi>
<mo>+</mo>
<mi>c</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>x</mi>
<mo>.</mo>
</mrow>
2.根据权利要求1所述的考虑作用点位置的基坑柔性支护主动土压力计算方法,其特征在于:所述eps取值为1×10-5。
3.根据权利要求2所述的考虑作用点位置的基坑柔性支护主动土压力计算方法,其特征在于:所述初定滑裂面中心坐标(x0,s0),x0=-2H,s0=2H。
4.根据权利要求2所述的考虑作用点位置的基坑柔性支护主动土压力计算方法,其特征在于:所述Φ(x0,s0)的最小值Φmin通过非线性函数优化方法求解,采用matlab软件中的fminsearch来进行搜索求解。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410672267.5A CN104484503B (zh) | 2014-11-20 | 2014-11-20 | 考虑作用点位置的基坑柔性支护主动土压力计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410672267.5A CN104484503B (zh) | 2014-11-20 | 2014-11-20 | 考虑作用点位置的基坑柔性支护主动土压力计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104484503A CN104484503A (zh) | 2015-04-01 |
CN104484503B true CN104484503B (zh) | 2017-09-26 |
Family
ID=52759044
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410672267.5A Expired - Fee Related CN104484503B (zh) | 2014-11-20 | 2014-11-20 | 考虑作用点位置的基坑柔性支护主动土压力计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104484503B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108197363A (zh) * | 2017-12-23 | 2018-06-22 | 中铁二院工程集团有限责任公司 | 膨胀土边坡刚性支挡结构设计中墙背侧向土压力分布计算方法 |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104794369B (zh) * | 2015-05-18 | 2017-07-25 | 重庆大学 | 基于土拱效应的抗滑桩桩间距及桩板土压力确定方法 |
CN107871026B (zh) * | 2016-09-26 | 2019-11-15 | 天津大学 | 一种适用于砂土中筒型基础正压下沉阻力计算方法 |
CN109635489B (zh) * | 2018-12-26 | 2022-12-16 | 天津大学 | 一种考虑桩靴脚印卸荷作用对邻近桩基础承载力影响的有限元计算方法 |
CN109918790B (zh) * | 2019-03-08 | 2023-05-12 | 何进 | 一种适用于坡顶荷载与降雨作用下边坡稳定性的判断方法 |
CN110824144B (zh) * | 2019-11-19 | 2022-08-16 | 中国科学院、水利部成都山地灾害与环境研究所 | 松散堆积体内摩擦角的确定方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS62146325A (ja) * | 1985-12-18 | 1987-06-30 | Taiyo Gijutsu Kaihatsu Kk | 擁壁の施工方法 |
CN101514553A (zh) * | 2009-04-03 | 2009-08-26 | 重庆交通大学 | 基于极限平衡理论和应力分析的土坡稳定分析方法 |
-
2014
- 2014-11-20 CN CN201410672267.5A patent/CN104484503B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS62146325A (ja) * | 1985-12-18 | 1987-06-30 | Taiyo Gijutsu Kaihatsu Kk | 擁壁の施工方法 |
CN101514553A (zh) * | 2009-04-03 | 2009-08-26 | 重庆交通大学 | 基于极限平衡理论和应力分析的土坡稳定分析方法 |
Non-Patent Citations (5)
Title |
---|
土压力合力作用点分析与计算;曹振民等;《公路》;19940630;第1994年卷(第6期);全文 * |
基于薄层单元法主动土压力计算的复合遗传算法;陈昌富等;《岩土力学》;20060331;第27卷(第3期);第2页左栏第3-4段,第2页右栏第5-9段,第3页右栏第9段,第4页左栏第1-2段 * |
对库伦主动土压力的误差分析与修正意见;沈雨生等;《科协论坛》;20120229;第2012年卷(第2期);全文 * |
挡土墙后曲面滑裂面下黏性主动土压力计算;王奎华等;《西南交通大学学报》;20111031;第46卷(第5期);全文 * |
计算柔性支挡结构主动土压力的简化方法;陈祖煜等;《岩土工程学报》;20100731;第32卷(第S1期);第2页左栏第4-6段,第2页右栏第5段,第3页左栏第1段,附图4 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108197363A (zh) * | 2017-12-23 | 2018-06-22 | 中铁二院工程集团有限责任公司 | 膨胀土边坡刚性支挡结构设计中墙背侧向土压力分布计算方法 |
CN108197363B (zh) * | 2017-12-23 | 2020-04-17 | 中铁二院工程集团有限责任公司 | 膨胀土边坡刚性支挡结构设计中墙背侧向土压力分布计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104484503A (zh) | 2015-04-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104484503B (zh) | 考虑作用点位置的基坑柔性支护主动土压力计算方法 | |
Adams et al. | The uplift capacity of shallow foundations | |
Braun et al. | Sediment transport mechanisms on soil-mantled hillslopes | |
Kubo | Experimental study of the behavior of laterally loaded piles | |
Seo et al. | Study on the applicability of a retaining wall using batter piles in clay | |
CN104794369B (zh) | 基于土拱效应的抗滑桩桩间距及桩板土压力确定方法 | |
CN104005777B (zh) | 一种大型地下洞室群布置设计方法 | |
CN108108334A (zh) | 改进的基坑整体滑动稳定安全系数计算方法 | |
Du et al. | Analysis of earth pressure for shallow square tunnels in anisotropic and non-homogeneous soils | |
Choi et al. | Rapid forecasting of tsunami runup heights from 2-D numerical simulations | |
Veiskarami et al. | Stability of sheet-pile walls subjected to seepage flow by slip lines and finite elements | |
CN104195956B (zh) | 膺架式水域连续梁现浇支架及其搭建方法 | |
CN105887941A (zh) | 一种桩基础桩长测算方法 | |
CN104573319B (zh) | 基于槽壁整体稳定性的富水砂卵地层地下连续墙成槽施工泥浆参数设计方法 | |
CN105421335B (zh) | 基于场地超孔隙水压比的水泥搅拌桩复合地基抗液化方法 | |
Jafari et al. | Three-dimensional levee and floodwall underseepage | |
Haimson | Design of underground powerhouses and the importance of preexcavation stress measurements | |
Youngs | Effect of the capillary fringe on steady-state water tables in drained lands | |
Vatankhah | Direct solutions for normal and critical depths in standard city-gate sections | |
Chen et al. | Using Boundary Transform Method to Solve Geotechnical Problems with Mixed-Type Boundary Conditions | |
Kacimov et al. | Steady-state water-table depressions caused by evaporation in lands overlying a water-bearing substratum | |
CN109724570A (zh) | 地下跌水的跌水量、跌水宽度、坎上水层厚度的计算方法 | |
CN206616513U (zh) | 兼顾床面稳定和快速排水的河工模型排水系统 | |
Veiskarami et al. | Stability of supported vertical cuts in granular matters in presence of the seepage flow by a semi-analytical approach | |
CN111046468A (zh) | 地下综合管廊抗震设计方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170926 Termination date: 20181120 |