CN108166545A - 一种任意荷载作用下多层地基的沉降计算方法 - Google Patents
一种任意荷载作用下多层地基的沉降计算方法 Download PDFInfo
- Publication number
- CN108166545A CN108166545A CN201810033124.8A CN201810033124A CN108166545A CN 108166545 A CN108166545 A CN 108166545A CN 201810033124 A CN201810033124 A CN 201810033124A CN 108166545 A CN108166545 A CN 108166545A
- Authority
- CN
- China
- Prior art keywords
- laplace
- layer
- stress
- load
- foundation
- 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.)
- Pending
Links
Classifications
-
- E—FIXED CONSTRUCTIONS
- E02—HYDRAULIC ENGINEERING; FOUNDATIONS; SOIL SHIFTING
- E02D—FOUNDATIONS; EXCAVATIONS; EMBANKMENTS; UNDERGROUND OR UNDERWATER STRUCTURES
- E02D33/00—Testing foundations or foundation structures
-
- E—FIXED CONSTRUCTIONS
- E02—HYDRAULIC ENGINEERING; FOUNDATIONS; SOIL SHIFTING
- E02D—FOUNDATIONS; EXCAVATIONS; EMBANKMENTS; UNDERGROUND OR UNDERWATER STRUCTURES
- E02D1/00—Investigation of foundation soil in situ
- E02D1/08—Investigation of foundation soil in situ after finishing the foundation structure
Landscapes
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Mining & Mineral Resources (AREA)
- Paleontology (AREA)
- Civil Engineering (AREA)
- General Engineering & Computer Science (AREA)
- Structural Engineering (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Soil Sciences (AREA)
- Consolidation Of Soil By Introduction Of Solidifying Substances Into Soil (AREA)
Abstract
一种任意荷载作用下多层地基的沉降计算方法,包括:(1)将固结方程进行Laplace变换,并对初始边值条件以及多层地基之间的层间连续条件进行Laplace变换;(2)根据第一步变换后得到的公式,推导特定层的应力解系数矩阵;(3)对任意荷载的波形函数进行Laplace变换,并作为变换域中的边界条件;(4)根据应力解系数矩阵及边界条件,采用矩阵传递法,求得各层在Laplace变换域中的应力及变形解答;(5)采用Laplace逆变换方法,即可得到任意点的应力及总沉降。本发明实现了对任意荷载作用下成层地基沉降规律的定量计算和预测分析,对工程中处于变荷载作用下地基的设计、施工提供科学与有效的依据。
Description
技术领域
本发明涉及一种任意荷载作用下多层地基的沉降计算方法,属于岩土工程技术领域。
背景技术
公路路堤、机场跑道、储藏设施(储油罐、堆场)以及工业厂房等建构筑物的地基所承受的荷载,往往是长时间往复施加的变荷载。这类地基在固结过程中受上述往复荷载作用,相比于静载作用下,会产生更大的附加变形及沉降。因此研究任意变化往复荷载作用下地基沉降,对于保护上方建构筑物及附属设施,具有重要的指导价值和应用意义。
Terzaghi之后,有许多学者对固结问题进行了卓有成效的研究。然而,目前针对变荷载下土体一维固结的计算方法,往往局限于特定荷载形式,如Wilson、Baligh等对矩形荷载下固结问题进行了研究,胡亚元分别对梯形与三角形荷载下的固结问题进行了研究,但这些研究不能针对工程中复杂的任意变化荷载形式进行计算。吴世明等推导了以积分形式表达的任意荷载的一维固结方程的通解,然而这是针对单层弹性地基进行的研究,工程中天然地基大多是由性状不一的层状土组成。
综上所述,目前还没有针对层状弹性地基在任意荷载作用下固结沉降的研究成果。
发明内容
本发明针对国内相关方面的理论研究空白,提出一种任意荷载作用下多层地基的沉降计算方法。
针对受任意荷载作用的多层地基模型,假定地基表面荷载随时间任意变化,其余假定同Terzaghi假定,如图2所示。
本发明提供一种任意荷载作用下多层地基的沉降计算方法,包括以下步骤:
步骤1:将固结方程进行Laplace变换,并对初始边值条件以及多层地基之间的层间连续条件进行Laplace变换;
对于第i层地基,固结方程为:
式中:
z为计算点深度,单位符号为m;
t为计算时刻,单位符号为s;
σ'i(z,t)为z深度土体在t时刻的有效应力,单位符号为Pa;
ki为第i层的土体渗透系数,单位为m/s;
Ei为第i层的弹性模量,单位为Pa;
γw为水的重度,单位符号为N/m3;
cvi为第i层的固结系数,单位符号为m2/s。
固结方程进行Laplace变换后可以得到:
s为Laplace变换参数。
对于双面透水地基,对问题的初值边界条件,及各层之间的层间应力连续、流量连续条件,进行Laplace变换,可以得到:
hi为第i层土层厚度,单位符号为m;
ri定义为
A11,A12,An1,An2等为计算参数;
步骤2:由Laplace变换域中的固结方程,根据Laplace变换后的边界条件,可以求解得:
由Laplace变换域中的应力连续、流量连续条件,可以得到第i层应力解系数矩阵与第i+1层应力解系数矩阵之间的关系:
[Di+1]=[Mi][Di];
其中Di为第i层应力解系数矩阵:
其中,Mi为两层系数矩阵之间的关系矩阵;αi表示为:αi=riki/ri+1ki+1;
步骤3:对外荷载进行Laplace变换,常见的荷载波形示意图如图3所示。
图3(a)的矩形波形荷载可以表示为:
p(t+2T)=p(t)
在Laplace变换域中的表达式为:
图3(b)的三角形波形荷载可以表示为:
p(t+2T)=p(t)
在Laplace变换域中的表达式为:
图3(c)的正弦波形荷载可以表示为:
在Laplace变换域中的表达式为:
图3(d)的梯形波形荷载可以表示为:
在Laplace变换域中的表达式为:
以上为常见的随时间变化的循环荷载,而对于任意荷载的波形函数,可以进行傅里叶级数展开成正余弦函数,随后再进行Laplace变换,得到Laplace变换域中荷载的表达式并作为变换域中的边界条件;
步骤4:根据第i层应力解系数矩阵,以及第i层与第i+1层应力解系数矩阵之间的关系,采用矩阵传递法,得到第i层与第1层土体应力解系数矩阵之间的关系:
随后根据作用于第一层土体上任意荷载,可以得到:
随后可以解得第1层土体的应力解系数矩阵:
其中,m11、m12、m21、m22分别表示相应系数;
至此已求得各层在Laplace变换域中的应力及变形解答:
步骤5:采用Laplace变换域中的应力及变形结果,采用Durbin方法进行Laplace逆变换,即可得到任意点的应力及总沉降。逆变换公式如下:
a,T,NUSM为逆变换参数。
附图说明
图1为本发明任意荷载作用下多层地基的沉降计算方法步骤流程图;
图2是本发明任意荷载作用下多层地基模型;
图3是常见的荷载波形示意图;
图3(a)为矩形波形载荷;图3(b)为三角形波形载荷;图3(c)为正弦波形载荷;图3(d)为梯形波形载荷;
图4是本发明应用于某三层地基的沉降-时间曲线图;
图5是某堆场上所作用的任意变化荷载示意图;
图6是本发明应用于某堆场的沉降-时间曲线图。
具体实施方式
下面结合附图,针对两个具体实例,对本发明的具体实施方法进行说明。
本发明的具体实施方式如图1所示。任意荷载作用下多层地基的沉降计算方法步骤流程如下:
(1)将固结方程进行Laplace变换;
(2)根据Laplace变换后的边界条件,推导特定层的应力解系数矩阵;
(3)对任意荷载的波形函数进行Laplace变换;
(4)求各层在Laplace变换域中的应力及变形解答;
(5)得到任意点的应力及总沉降。
实施例1:
本发明实施例以某三层地基作为案例,进行算例分析。
在该三层地基分别上作用矩形、三角形、正弦波形、梯形波载及瞬加荷载,考虑荷载周期T=30天,荷载幅值P0=1.5*105Pa。三层地基厚度分别为1m、1.5m、1m,弹性模量分别为1*107Pa、2*107Pa、1*107Pa,渗透系数分别为3*10-10m/s、2*10-10m/s、3*10-10m/s。
计算得到在多种不同循环荷载作用下,该三层地基的沉降随时间变化曲线,如图4所示,可以看出,随着时间的增长,瞬加荷载作用下沉降随时间增长并趋于稳定,而循环荷载作用下的成层地基沉降变形增长模式与瞬加荷载作用下沉降变形增长模式有本质区别,循环荷载下沉降增长与荷载形式有较大关系。将循环荷载等效为瞬加荷载,将会导致计算沉降有较大的误差。
实施例2:
本发明实施例以某进出口堆场作为计算案例,进行计算。该堆场修建于软土地基上,软土地基经过地基处理后可视为双层地基,第一层为处理层,厚度16.5m,地基处理后刚度为6.85*106Pa,渗透系数为3*10-9m/s。第二层为未经地基处理的软土层,厚度19m,刚度为2.95*106Pa,渗透系数为1.1*10-7m/s。
在堆场投入正式使用之前,有一个预载期,预载期堆载与时间曲线如图5所示。预载期堆载大小P0=75kPa,T=100天,t1=70天。该荷载不同于循环变化的规则荷载,为任意荷载,首先进行傅里叶级数展开成正余弦函数,随后再进行Laplace变换,该任意变化荷载在Laplace变换域中的表达式为:
其中,
图6为计算得到的预载期沉降,以及相应的实测沉降值,计算结果表明经过地基处理的处理层沉降相对较小,沉降主要由下卧层沉降组成。本发明计算得到的总沉降与实测总沉降吻合良好,证明了本发明的适用性。
Claims (1)
1.一种任意荷载作用下多层地基的沉降计算方法,其特征在于,所述方法包括以下步骤:
步骤1:将固结方程进行Laplace变换,并对初始边值条件以及多层地基之间的层间连续条件进行Laplace变换;
假定地基表面荷载随时间任意变化,其余假定同Terzaghi假定,可以得到第i层土体的固结公式:
式中:z为计算点深度,单位符号为m;t为计算时刻,单位符号为s;σi'(z,t)为z深度土体在t时刻的有效应力,单位符号为Pa;ki为第i层的土体渗透系数,单位为m/s;Ei为第i层的弹性模量,单位为Pa;γw为水的重度,单位符号为N/m3;cvi为第i层的固结系数,单位符号为m2/s;
固结方程进行Laplace变换后可以得到:
s为Laplace变换参数;
对于双面透水地基,对问题的初值边界条件,及各层之间的层间应力连续、流量连续条件,进行Laplace变换,可以得到:
hi为第i层土层厚度,单位符号为m;
ri定义为ri 2=s/cvi;
A11,A12,An1,An2等为计算参数;
步骤2:由Laplace变换域中的固结方程,根据Laplace变换后的边界条件,可以求解得:
由Laplace变换域中的应力连续、流量连续条件,可以得到第i层应力解系数矩阵与第i+1层应力解系数矩阵之间的关系:
[Di+1]=[Mi][Di];
其中Di为第i层应力解系数矩阵:
其中,Mi为两层系数矩阵之间的关系矩阵;αi为αi=riki/ri+1ki+1;
步骤3:对任意荷载的波形函数,可以进行傅里叶级数展开成正余弦函数,随后再进行Laplace变换,得到Laplace变换域中荷载的表达式并作为变换域中的边界条件;
步骤4:根据第i层应力解系数矩阵,以及第i层与第i+1层应力解系数矩阵之间的关系,采用矩阵传递法,得到第i层与第1层土体应力解系数矩阵之间的关系:
随后根据作用于第一层土体上任意荷载,可以得到:
其中,
随后可以解得第1层土体的应力解系数矩阵:
其中,m11、m12、m21、m22分别表示相应系数;
至此已求得各层在Laplace变换域中的应力及变形解答:
步骤5:采用Laplace变换域中的应力及变形结果,采用Durbin方法进行Laplace逆变换,即可得到任意点的应力及总沉降;逆变换公式如下:
a,T,NUSM为逆变换参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810033124.8A CN108166545A (zh) | 2018-01-14 | 2018-01-14 | 一种任意荷载作用下多层地基的沉降计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810033124.8A CN108166545A (zh) | 2018-01-14 | 2018-01-14 | 一种任意荷载作用下多层地基的沉降计算方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108166545A true CN108166545A (zh) | 2018-06-15 |
Family
ID=62514248
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810033124.8A Pending CN108166545A (zh) | 2018-01-14 | 2018-01-14 | 一种任意荷载作用下多层地基的沉降计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108166545A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108984969A (zh) * | 2018-08-22 | 2018-12-11 | 华东交通大学 | 一种软土地基盾构隧道运营期沉降计算方法 |
CN109518573A (zh) * | 2018-11-30 | 2019-03-26 | 长沙理工大学 | 车辆荷载作用下路基动附加应力的计算方法 |
CN112881241A (zh) * | 2021-01-19 | 2021-06-01 | 华东交通大学 | 一种确定颗粒材料模量软化和恢复的方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0892990A (ja) * | 1994-09-27 | 1996-04-09 | Mitsubishi Heavy Ind Ltd | 土構造物の沈下発生装置 |
CN103953023A (zh) * | 2014-04-08 | 2014-07-30 | 珠海电力设计院有限公司 | 一种真空联合堆载预压法沉降量的测量方法 |
-
2018
- 2018-01-14 CN CN201810033124.8A patent/CN108166545A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0892990A (ja) * | 1994-09-27 | 1996-04-09 | Mitsubishi Heavy Ind Ltd | 土構造物の沈下発生装置 |
CN103953023A (zh) * | 2014-04-08 | 2014-07-30 | 珠海电力设计院有限公司 | 一种真空联合堆载预压法沉降量的测量方法 |
Non-Patent Citations (1)
Title |
---|
徐长节等: "任意荷载下成层弹性地基的一维固结", 《土木工程学报》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108984969A (zh) * | 2018-08-22 | 2018-12-11 | 华东交通大学 | 一种软土地基盾构隧道运营期沉降计算方法 |
CN109518573A (zh) * | 2018-11-30 | 2019-03-26 | 长沙理工大学 | 车辆荷载作用下路基动附加应力的计算方法 |
CN112881241A (zh) * | 2021-01-19 | 2021-06-01 | 华东交通大学 | 一种确定颗粒材料模量软化和恢复的方法 |
CN112881241B (zh) * | 2021-01-19 | 2022-10-28 | 华东交通大学 | 一种确定颗粒材料模量软化和恢复的方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Cai et al. | Effect of anisotropic consolidation stress paths on the undrained shear behavior of reconstituted Wenzhou clay | |
CN108166545A (zh) | 一种任意荷载作用下多层地基的沉降计算方法 | |
Brocklehurst et al. | Interaction of hydro-elastic waves with a vertical wall | |
Dan et al. | Coupled hydro-mechanical response of saturated asphalt pavement under moving traffic load | |
Xie et al. | Analytical theory for one‐dimensional consolidation of clayey soils exhibiting rheological characteristics under time‐dependent loading | |
Li et al. | Dynamic response of a stratified transversely isotropic half-space with a poroelastic interlayer due to a buried moving source | |
Wang et al. | Perturbation analysis for upheaval buckling of imperfect buried pipelines based on nonlinear pipe-soil interaction | |
Asgari et al. | Numerical evaluation of seismic response of shallow foundation on loose silt and silty sand | |
Liu et al. | A general solution for 1D consolidation induced by depth-and time-dependent changes in stress | |
Pender et al. | Effect of permeability on the cyclic generation and dissipation of pore pressures in saturated gravel layers | |
Li et al. | Settlement prediction of surcharge preloaded low embankment on soft ground subjected to cyclic loading | |
Dammala et al. | Nonlinear seismic ground response analysis in northeastern India considering the comprehensive dynamic soil behavior | |
Qian et al. | Dynamic shakedown analysis of flexible pavement under traffic moving loading | |
Oumeraci et al. | Hydraulic performance, wave loading and response of PBA revetments and their foundations | |
Chen et al. | Effects of local scouring on capacities of monopile-wheel composite foundations in sand under lateral loading | |
Ghalesari et al. | The settlement behavior of piled raft interaction in undrained soil | |
Panique Lazcano et al. | Long-term dynamic bearing capacity of shallow foundations on a contractive cohesive soil | |
Paul et al. | One dimensional consolidation under cyclic loading | |
Ram et al. | 1D, 2D, and 3D Seismic Response Analysis of Shallow and Deep Foundation Resting on Stratified Soil Deposit | |
Ausilio et al. | Seismic response of alluvial valleys to SH waves | |
Chen et al. | Dynamic constitutive model for soils considering asymmetry of skeleton curve | |
Shafiee et al. | Application of fuzzy set theory to evaluate the effect of pore pressure build-up on the seismic stability of Karkheh dam, Iran | |
Majidi et al. | Estimating compaction parameters of marl soils using multi-layer perceptron neural networks | |
Banerjee et al. | Numerical analysis of piles under cyclic lateral load | |
Sawicki et al. | Wave-induced stresses and pore pressures near a mudline |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20180615 |