CN104452770B - 由岩土本构折线段确定土体内力的方法 - Google Patents

由岩土本构折线段确定土体内力的方法 Download PDF

Info

Publication number
CN104452770B
CN104452770B CN201410566691.1A CN201410566691A CN104452770B CN 104452770 B CN104452770 B CN 104452770B CN 201410566691 A CN201410566691 A CN 201410566691A CN 104452770 B CN104452770 B CN 104452770B
Authority
CN
China
Prior art keywords
internal force
broken line
line segment
stress
sigma
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
Application number
CN201410566691.1A
Other languages
English (en)
Other versions
CN104452770A (zh
Inventor
白福波
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Harbin Institute Of Technology Architectural Design And Research Institute Co ltd
Original Assignee
Harbin Institute of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Harbin Institute of Technology filed Critical Harbin Institute of Technology
Publication of CN104452770A publication Critical patent/CN104452770A/zh
Application granted granted Critical
Publication of CN104452770B publication Critical patent/CN104452770B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • EFIXED CONSTRUCTIONS
    • E02HYDRAULIC ENGINEERING; FOUNDATIONS; SOIL SHIFTING
    • E02DFOUNDATIONS; EXCAVATIONS; EMBANKMENTS; UNDERGROUND OR UNDERWATER STRUCTURES
    • E02D17/00Excavations; Bordering of excavations; Making embankments
    • E02D17/02Foundation pits
    • EFIXED CONSTRUCTIONS
    • E02HYDRAULIC ENGINEERING; FOUNDATIONS; SOIL SHIFTING
    • E02DFOUNDATIONS; EXCAVATIONS; EMBANKMENTS; UNDERGROUND OR UNDERWATER STRUCTURES
    • E02D1/00Investigation of foundation soil in situ

Landscapes

  • Engineering & Computer Science (AREA)
  • Mining & Mineral Resources (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Paleontology (AREA)
  • Civil Engineering (AREA)
  • General Engineering & Computer Science (AREA)
  • Structural Engineering (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Soil Sciences (AREA)
  • Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
  • Bulkheads Adapted To Foundation Construction (AREA)

Abstract

由岩土本构折线段刚度确定土体内力的方法,涉及基坑侧壁直立开挖深度和土体内力的确定方法,属于地基岩土工程技术领域。为了解决当应用已有的本构模型进行场地土体内力、应力分析时,由于场地土体本构条件状况变异离散,过程运算繁杂,只能粗略简化为概念确定,影响设计的准确性,存在难以满足经济性要求等问题。技术要点:将与应力相关的土工实验曲线,按要求的精确度采用割线或切线模拟方式;在坐标系内进行区间划分。形成一多段折线形状,每段折线的折点坐标或斜率为已知,建立折线段方程可得主应力圆与抗剪强度关系式及土体中内力与变位的折线性设计方法。本发明适用于基坑侧壁直立挖掘。

Description

由岩土本构折线段确定土体内力的方法
本申请是申请日为2013年04月03日、申请号为2013101152721、发明名称为由岩土本构折线段确定基坑挖深的方法的发明申请的分案申请。
技术领域
本发明涉及一种基坑挖深的确定方法,属于地基岩土工程技术领域。
背景技术
岩土是一种由天然介质混合构成的弹塑性材料,本构关系随自然条件而变。一般是在现场勘探取出原状土样试件,在实验室内由土工实验获得物理力学性能指标及相应参数,如:应力与抗剪强度包线、压缩曲线、侧限变形模量等。当前的做法是将变异离散的实验数据(坐标点),用数学方法统计回归后取过坐标原点近于直线那一段,这段线性特征在材料上被视为具有弹性性能,并作为土体中内力与变位的设计依据,也是建立解算线性方程的广义物理条件。这一表述说明:这种线性设计实质是一种简化方法,对一些土工实验中非线性特征明显的土体材料,产生的误差较大,也不适用。目前对这一土工实验中的非线性问题,已构建了许多本构模型的表达式,但用于场地土体内力、应力分析时,则是一极繁琐的行为,难以满足要求。当场地土体本构条件复杂时,只能粗略简化为概念确定。
发明内容
本发明为了解决利用已有的本构模型进行场地土体内力、应力分析(进行地基承载、基坑侧壁土压、岩土结构设计等),当场地土体本构条件变异离散运算繁杂时,只能粗略简化分析,影响设计的准确性,难以满足工程要求等问题,进而提供了一种由岩土本构折线段确定基坑挖深的方法。
本发明为解决上述技术问题采取的技术方案是:
技术方案一:一种由岩土本构折线段确定基坑挖深的方法,所述方法是按照以下步骤实现的:
步骤一、选取折线段:
折线各折点左右线段判定式为:
当σ3,i≤σ3j取左Lj-1,σ3,i3j取右Lj;
σ 3 , i = Σ i = 0 i + 1 Δσ 3 , i ; i=0;△σ3,0=0;
步骤二、根据选取的折线段利用式(1)或式(4)求得△σ1,△σ1表示基坑侧壁的土体单元垂直向压应力增量,即折线段大主压应力增量;
Δσ 1 = Σ i = 1 n Δσ 13 , i - - - ( 1 )
式中: Δσ 13 , i = ( σ 3 , i k j - σ 3 , i - 1 k j - 1 ) + 2 ( C j k j - C j - 1 k j - 1 ) - - - ( 2 )
σ1,i=σ1,0+△σ1 σ 1 , 0 = 2 C 0 k 0
Δσ 1 = σ 3 , i k j + 2 ( C j k j - C 0 k 0 ) - - - ( 4 ) ;
大主压应力增量公式(1)可写成公式(4);式(1)用于电算,式(4)用于手算;σ1,j、σ3,j、τf1,j、τf3,j:折线各折点的应力、剪应力(坐标值);
字符下标i:坐标轴上i点应力。下标j:第j条折线段Lj;
步骤三、求侧压应力下基坑侧壁深度增量△H1:△H1=△σ1
式中:γ表示土重度;
△H1,将用于确定基坑侧壁直立总深度H。
技术方案二:一种由岩土本构折线段确定基坑挖深的方法,所述方法是按照以下步骤实现的:
步骤一、选取折线段:
折线各折点左右线段判定式为:
当σ1,i≤σ1j取左Lj-1,σ1,i1j取右Lj;
σ 1 , i = Σ i = 0 i + 1 Δσ 1 , j ; i=0;△σ1,0=0;
步骤二、根据选取的折线段利用式(5)或式(8)求得△σ3,△σ3表示
基坑侧壁的土体单元水平向压应力增量,即折线段小主压应力增量;
Δσ 3 = Σ i = 1 n Δσ 31 , i - - - ( 5 )
Δσ 31 , i = ( σ 1 , i k j - σ 1 , i - 1 k j - 1 ) - 2 ( C j k j - C j - 1 k j - 1 ) - - - ( 6 )
Δσ 3 = σ 1 , i k j - 2 ( C j k j - C 0 k 0 ) - - - ( 8 )
小主压应力增量公式(5)可写成公式(8);式(5)用于电算,式(8)用于手算;式中:σ1,i、σ3,i:大、小主压应力;
△σ1、△σ3:大、小主压应力增量;
Cj土体内粘聚力、摩擦角;
σ1,j、σ3,j、τf1,j、τf3,j:折线各折点的应力、剪应力(坐标值);
字符下标i:坐标轴上i点应力;下标j:第j条折线段Lj;
步骤三、求侧压应力下基坑侧壁深度减量△H3:△H3=△σ3
式中:γ表示土重度;
由技术方案一得到的△H1及△H3,可求得基坑侧壁直立总深度H。
技术方案三:本技术方案也适于确定基坑挖深,一种由岩土本构折线段刚度选定确定土体内力的方法,所述方法是按照以下步骤实现的:
步骤一、先选择一条实刚度线段,得线刚度方程S=S0+R/K,由已知外荷载按线性关系求解内力R(r)。
步骤二、若内力R(r)或R(r′)落在实刚度线段Li所对应的内力范围内,如:Ri≤R(r)≤Ri+1,则计算值R(r)即是所求内力值;取R(r)值再按材料力学求解基坑土体内:正应力、剪应力、主应力;
步骤三、若计算的内力R(r)或R(r′)未落在选择的Li所对应的内力范围内,那么r或r′点必对应出现在Li线段左或右端点以外延伸线的虚刚度线段上一点a或a′;由a或a′对内力轴引垂线必与此相交有一条实刚度线段b或b′点,可再得实刚度线段方程,然后重复步骤一、二。
本发明的有益效果是:
利用已有的线弹性设计法,将与应力相关的土工实验曲线,按要求的精确度采用割线或切线方式;在坐标系内进行区间划分。形成一多段折线形状,每段折线的折点坐标或斜率为已知,建立直线方程可得主应力圆与抗剪强度关系式及土体中内力与变位的非线性设计方法。它相对准确有效的解决岩土非线性问题,如:用于确定地基承载、基坑侧壁土压、岩土结构设计等。作法适用合理,节能省材。
附图说明
图1是折线性大主压应力图,图2是折线性小主压应力图,图3是选折线刚度设计方法示意图(图3中,①:实刚度线段、②:虚刚度线段);图4、图5为平面坐标应力图解(图4表示:大主压应力图解;图5表示小主压应力图解)。
具体实施方式
由锚杆担梁双桩单元对基坑侧壁施加预压力。采用PC管桩,桩直径0.3m,桩距1.0m,锚杆水平间距2.0m。土重度:19.数值单位:KN,Kpa,m.
土工试验得抗剪强度包线为一曲线,再化分模拟成三段直线呈折线形。
线段L1:C1=17;C2=11;C3=24。三折线段有两个交点(折点)坐标:(σ3,1f3,1)=(σ1,1f1,1)=(30,23);(σ3,2f3,2)=(σ1,2f1,2)=(63,37)。侧压应力P=26(△σ3);q=22(△σ1)。有:σ3,i=3·△σ3=3×26=78;σ1,i=3·△σ1=3×22=66Kpa
实施例1:求侧压应力下基坑侧壁深度增量△H1。
由式(2)与判别式(3),取第二折点坐标值(63,37)左侧线段L2.
σ3j=63-37/cos23°(1-sin23°)=63-37/0.92(1-0.39)=38.47Kpa
σ3,i=78>σ3j=38.47,取右侧折线段L3
k1=tg2(45°+10°/2)=1.1922=1.42;k2=tg2(45°+23°/2)=1.5112=2.283
k3=tg2(45°+12°/2)=1.2352=1.525.
由式(1),(2)或式(4)得:
△σ1=78×1.525+2×(24×√1.525-17×√1.42)=118.95+18.76=137.71Kpa
基坑深度增量:△H1=△σ1/γ=137.71/19=7.25m
实施例2:求竖向压应力作用下,基坑侧壁深度减量△H3。
由式(6)与判别式(7),取第一折点坐标值(30,23)右侧线段L2
σ1j=30+23/cos23°(1+sin23°)=30+23/0.921×(1+0.391)=64.74Kpa
σ1,i=66>σ1j=64.74,取右侧折线段L2.
由式(5),(6)或式(8)得:
△σ3=66/2.283-2×(11/√2.283-17/√1.42)=28.91+13.96=42.87Kpa
基坑深度减量:△H3=△σ3/γ=42.87/19=2.26m
实施例3:求基坑侧壁深度H
基坑初始临界深度:H0c=2×17/19×tg(45°+10°/2)=2.13m;基坑顶放坡深1.5m。基坑侧壁深度:H=H0c+△H1-△H3=2.13+7.25-2.26=7.12m(取两道锚杆)
基坑侧壁开挖深度:h=7.12+1.5=8.62m
△H1/H×P=7.25/7.12×26=26.48;Npi=26.48×1.0×7.12/3=62.85kN
单束锚杆侧向预拉力取2Npi=126kN。(选一根钢铰线)
图4、图5为实施例平面坐标图解应力。
本发明内容的论证:
一、由土工试验抗剪强度曲线(为非线性可含有突变断点),它可由一系列满足精度要求的切线或割线线段包络,得到一模拟接近广义非线性曲线的多段折线。因多段折线斜率一点座标为已知,按截距式方程,对于单直线段:j=1对于广义折线段: Δσ 1 = Σ i = 1 n Δσ 13 , i - - - ( 1 )
式中: Δσ 13 , i = ( σ 3 , i k j - σ 3 , i - 1 k j - 1 ) + 2 ( C j k j - C j - 1 k j - 1 ) - - - ( 2 )
式中:n≥3;n:增量分段。
△σ3,i=σ3,i3,i-1;i=1;σ3,0=0;j=1,2,3,…。图3中,Lj为折线线段。
σ1,i=σ1,0+△σ1 σ 1 , 0 = 2 C 0 k 0
图1为莫尔圆与非线性抗剪强度折线图,也是公式(1)的几何图解。应力圆在向右按增量(步长)移动时,须确定圆与折线的切点位于折点左右那一侧线段上,对折线各折点j左右线段Lj的选取判定,可由几何关系得折线分界点左右线段判定式。由已知折点坐标(σ3,jf3,j),见图1可写得算式:
当σ3,i≤σ3j取左Lj-1,σ3,i3j取右Lj。
σ 3 , i = Σ i = 0 i + 1 Δσ 3 , i ; i=0;△σ3,0=0
图1表明剪切包线在平面坐标系内可以是任意折线段,折线上一点的应力状态与主应力增量步长路径累加无关。可由式(3)判定主应力对应的折线线段,得主应力增量式(4)。
大主压应力增量公式也可写成:
Δσ 1 = σ 3 , i k j + 2 ( C j k j - C 0 k 0 ) - - - ( 4 )
二、按上述方法由图2可得小主压应力公式:
Δσ 3 = Σ i = 1 n Δσ 31 , i - - - ( 5 )
Δσ 31 , i = ( σ 1 , i k j - σ 1 , i - 1 k j - 1 ) - 2 ( C j k j - C j - 1 k j - 1 ) - - - ( 6 )
折线各折点j左右线段Lj的判定式:
当σ1,i≤σ1j取左Lj-1,σ1,i1j取右Lj。
小主压应力增量公式也可写成:
Δσ 3 = σ 1 , i k j - 2 ( C j k j - C 0 k 0 ) - - - ( 8 )
注:图1,2中,圆弧虚线与应力轴交点为分界圆应力点。
三、折线刚度选定法:图3为内力R与位移S曲线的折线性模拟,折线段包络角与根数满足精度要求。在外荷载作用下折线上一点处于平衡状态。通过折线(实刚度线段)上的某一点求解唯一对应内力R。计算步骤是:
1、先选择一条实刚度线段,得线刚度方程:S=S0+R/K由已知外荷载按线性关系求解内力R(r)。
2、若内力R(r)或R(r′)落在实刚度线段Li所对应的内力范围内,如:Ri≤R(r)≤Ri+1,则计算值R(r)是所求内力。再按材料力学求解正应力、剪应力、主应力。
3、若计算的内力R(r)或R(r′)未落在选择的Li所对应的内力范围内,那么r或r′点必对应出现在Li线段左或右端点以外延伸线上的虚刚度线段②一点a或a′。由a或a′对内力轴引垂线必与此相交有一条实刚度线段①上一点b或b′,可再得实刚度线段方程,重复步骤1、2。
选择一条实刚度线段的过程为:先将曲线模拟化分成三条折线段,分别运算必有一条可选择为“实刚度线段”。

Claims (1)

1.由岩土本构折线段确定土体内力的方法,其特征在于:所述方法是按照以下步骤实现的:
步骤一、先选择一条实刚度线段,得线刚度方程S=S0+R/K,由已知外荷载按线性关系求解内力R(r);
步骤二、若内力R(r)或R(r′)落在实刚度线段Li所对应的内力范围内,如:Ri≤R(r)≤Ri+1,则计算值R(r)即是所求内力值;取R(r)值再按材料力学求解基坑土体内:正应力、剪应力、主应力;
步骤三、若计算的内力R(r)或R(r′)未落在选择的Li所对应的内力范围内,那么r或r′点必对应出现在Li线段左或右端点以外延伸线的虚刚度线段上一点a或a′;由a或a′对内力轴引垂线必与此相交有一条实刚度线段b或b′点,可再得实刚度线段方程,然后重复步骤一、二;
选择一条实刚度线段的过程为:先将曲线模拟化分成三条折线段,分别运算必有一条可选择为“实刚度线段”。
CN201410566691.1A 2013-04-03 2013-04-03 由岩土本构折线段确定土体内力的方法 Expired - Fee Related CN104452770B (zh)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310115272.1A CN103195068B (zh) 2013-04-03 2013-04-03 由岩土本构折线段确定基坑挖深的方法

Related Parent Applications (1)

Application Number Title Priority Date Filing Date
CN201310115272.1A Division CN103195068B (zh) 2013-04-03 2013-04-03 由岩土本构折线段确定基坑挖深的方法

Publications (2)

Publication Number Publication Date
CN104452770A CN104452770A (zh) 2015-03-25
CN104452770B true CN104452770B (zh) 2016-06-01

Family

ID=48717871

Family Applications (2)

Application Number Title Priority Date Filing Date
CN201310115272.1A Expired - Fee Related CN103195068B (zh) 2013-04-03 2013-04-03 由岩土本构折线段确定基坑挖深的方法
CN201410566691.1A Expired - Fee Related CN104452770B (zh) 2013-04-03 2013-04-03 由岩土本构折线段确定土体内力的方法

Family Applications Before (1)

Application Number Title Priority Date Filing Date
CN201310115272.1A Expired - Fee Related CN103195068B (zh) 2013-04-03 2013-04-03 由岩土本构折线段确定基坑挖深的方法

Country Status (1)

Country Link
CN (2) CN103195068B (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103938636A (zh) * 2014-04-21 2014-07-23 国家电网公司 一种输电线路工程岩石基础利用液压劈裂机开挖的方法
CN104328787B (zh) * 2014-10-30 2016-06-22 水利部交通运输部国家能源局南京水利科学研究院 一种土的开挖回弹过程模拟方法
CN104878785B (zh) * 2015-05-18 2017-06-30 卢里尔 桩土相互作用及桩端土本构模型和模型参数确定方法
CN104915574B (zh) * 2015-07-03 2017-08-29 长江水利委员会长江科学院 一种适用于加卸载全过程的软岩蠕变本构模型的建立方法
CN107103168B (zh) * 2017-06-23 2020-11-06 哈尔滨工业大学建筑设计研究院 基坑侧壁挡土结构设计方法
CN112557218B (zh) * 2020-12-24 2024-03-19 中国建筑第八工程局有限公司 通过比较截距判断黏土是否受力破坏的方法
CN114032973B (zh) * 2021-11-24 2023-05-12 上海建工集团股份有限公司 一种基坑深层土体变形监测装置及计算方法

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0776840A (ja) * 1993-09-06 1995-03-20 Sabukon Syst:Kk 山留め計算処理支援システム
JP3353273B2 (ja) * 2001-05-09 2002-12-03 パシフィックコンサルタンツ株式会社 掘削工事用シミュレーション装置に用いる背面側の水平地盤反力係数設定方法及び最小土圧設定方法
CN101871776A (zh) * 2009-04-27 2010-10-27 温兆东 一种排桩桩身变形量的计算方法
CN101615212B (zh) * 2009-07-28 2011-10-26 武汉理工大学 一种排桩-内支撑支护体系简化设计计算方法
CN102561293B (zh) * 2010-12-17 2016-11-23 上海市电力公司 软土深基坑水土压力测试方法
CN102251524A (zh) * 2011-04-25 2011-11-23 哈尔滨工业大学 基坑侧壁预应力锚杆担梁排桩围护结构
CN102383433A (zh) * 2011-08-22 2012-03-21 同济大学建筑设计研究院(集团)有限公司 一种超高层建筑深大基坑的设计方法
CN103234473A (zh) * 2013-04-12 2013-08-07 哈尔滨工业大学深圳研究生院 一种可观测土体内部的渗水力模型试验装置和试验方法

Also Published As

Publication number Publication date
CN103195068A (zh) 2013-07-10
CN103195068B (zh) 2015-05-13
CN104452770A (zh) 2015-03-25

Similar Documents

Publication Publication Date Title
CN104452770B (zh) 由岩土本构折线段确定土体内力的方法
Basile Effects of tunnelling on pile foundations
Andrianopoulos et al. Bounding surface plasticity model for the seismic liquefaction analysis of geostructures
Houlsby et al. A plasticity model for the behaviour of footings on sand under combined loading
Nezami et al. Simulation of front end loader bucket–soil interaction using discrete element method
Huang et al. Axial capacity degradation of single piles in soft clay under cyclic loading
Smith et al. Axisymmetric numerical modeling of a unit cell in geosynthetic-reinforced, column-supported embankments
Huang et al. Effect of sea ice on seismic collapse-resistance performance of wind turbine tower based on a simplified calculation model
Duncan Slope stability then and now
Akhmetova et al. The visual representation of numerical solution for a non-stationary deformation in a solid body
CN110263448A (zh) 一种双隧道开挖对邻近地理管线影响的预测方法
CN106202654A (zh) 一种地下厂房多块体围岩变形的数值模拟方法
Yeter et al. Fragility analysis of an ageing monopile offshore wind turbine subjected to simultaneous wind and seismic load
Robert et al. Analysis of soil-pipeline interaction using ABAQUS/explicit
Bergholz An extended bounding surface model for the application to general stress paths in sand
Yang et al. Computation and analysis of high rocky slope safety in a water conservancy project
Sultanov et al. Method for assessing the reliability of earth dams in irrigation systems
CN206601298U (zh) 土体弹性参数测试装置
Masoumi et al. Optimal monitoring instruments selection using innovative decision support system framework
Numada et al. Material point method for run-out analysis of earthquake-induced long-traveling soil flows
Sutubadi et al. Effect of soil properties on stability of soil-steel culverts
Ibrahim et al. The effects of parameter uncertainties on the numerical plastic analysis of non-prismatic reinforced concrete beams
KR100528572B1 (ko) 물량 계산 방법
Abdelrahman et al. 3D numerical analysis of corner mechanically stabilized earth walls
Yousefi Samangani et al. Numerical and Larg‐Scale Laboratory Study of Rock Column Groups in Sandy Soil Behavior Improvement

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
ASS Succession or assignment of patent right

Owner name: CONSTRUCTION DESIGN RESEARCH INST., HARBIN POLYTEC

Free format text: FORMER OWNER: HARBIN INDUSTRY UNIVERSITY

Effective date: 20150521

C41 Transfer of patent application or patent right or utility model
TA01 Transfer of patent application right

Effective date of registration: 20150521

Address after: 150001 No. 73, the Yellow River Road, Nangang District, Heilongjiang, Harbin

Applicant after: BUILDING DESIGN Research Institute HARBIN INSTITUTE OF TECHNOLOGY

Address before: 150001 Harbin, Nangang, West District, large straight street, No. 92

Applicant before: Harbin Institute of Technology

C14 Grant of patent or utility model
GR01 Patent grant
CP01 Change in the name or title of a patent holder
CP01 Change in the name or title of a patent holder

Address after: 150001 No. 73, Huanghe Road, Nangang District, Harbin City, Heilongjiang Province

Patentee after: Harbin Institute of Technology Architectural Design and Research Institute Co.,Ltd.

Address before: 150001 No. 73, Huanghe Road, Nangang District, Harbin City, Heilongjiang Province

Patentee before: BUILDING DESIGN Research Institute HARBIN INSTITUTE OF TECHNOLOGY

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: 20160601