CN110990924A - 一种注浆工程加固路基浆液扩散范围公式推导方法 - Google Patents

一种注浆工程加固路基浆液扩散范围公式推导方法 Download PDF

Info

Publication number
CN110990924A
CN110990924A CN201911194674.9A CN201911194674A CN110990924A CN 110990924 A CN110990924 A CN 110990924A CN 201911194674 A CN201911194674 A CN 201911194674A CN 110990924 A CN110990924 A CN 110990924A
Authority
CN
China
Prior art keywords
slurry
grouting
diffusion
curve
range
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.)
Granted
Application number
CN201911194674.9A
Other languages
English (en)
Other versions
CN110990924B (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.)
Southwest Jiaotong University
Original Assignee
Southwest Jiaotong University
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 Southwest Jiaotong University filed Critical Southwest Jiaotong University
Priority to CN201911194674.9A priority Critical patent/CN110990924B/zh
Publication of CN110990924A publication Critical patent/CN110990924A/zh
Application granted granted Critical
Publication of CN110990924B publication Critical patent/CN110990924B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/11Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G06F17/12Simultaneous equations, e.g. systems of linear equations

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • Operations Research (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Consolidation Of Soil By Introduction Of Solidifying Substances Into Soil (AREA)

Abstract

本发明涉及一种注浆工程加固路基浆液扩散范围公式推导方法,包括以下步骤:1)利用CFD‑DEM耦合计算方法,模拟计算出浆液‑土体双向耦合微观层面上的破坏范围;2)选取一种工况,将破坏范围导入相关软件中并将其描绘出包络图;3)将包络图以数学表达式的形式表达出来;4)将其他注浆压力下的浆液扩散曲线数学表达式做出假设;5)建立多个方程组对表达式的系数进行求解,得出曲线方程。本发明对注浆浆液扩散范围的计算方法进行了推导,得出了浆液扩散曲线的公式,填补了国内研究空白,并为未来注浆工程的研究提出了可行的方向。

Description

一种注浆工程加固路基浆液扩散范围公式推导方法
技术领域
本发明涉及于路基注浆的研究技术领域,具体为一种注浆工程加固路基浆液扩散范围公式推导方法。
背景技术
本发明涉及路基注浆的研究技术领域,属于注浆过程中浆液扩散范围的理论计算方法。
随着我国高铁网络建设的日益完善,现场建设技术日趋成熟,但是随着高速铁路的布局完善,新的工程问题也出现。到2018年底,中国高速铁路营业里程达2.9万多公里。中国幅员辽阔,在这么长的高速运营里程中,难免不会出现各种各样的地质情况,从而导致路基沉降,对高速铁路的运营造成重大的安全隐患,特别是在软土地基上的情况。而注浆目前作为操作方便,止沉效果卓越的一种沉降治理方法,受到各界研究人员的重视。然而,由于土体的隐蔽性和土体性质的复杂性,从而给注浆工程的研究带来了困难。
发明内容
发明目的:目前关于软土地基中注浆治理沉降的方法基本以注浆止沉为主,但是关于注浆工程的研究工作尚不完善,为填补国内研究的空白,本发明基于CFD-DEM双向耦合计算注浆全过程,观察土体微观层面上接触键的破坏规律,并利用多种软件组合应用,解出了浆液在土体中的扩散计算公式。为未来的注浆工程的研究提出了可行的方向。
为达到上述目的,本发明的技术方案为:
一种注浆工程加固路基浆液扩散范围公式推导方法,包括以下步骤:
一种注浆工程加固路基浆液扩散范围公式推导方法,其特征在于,包括以下步骤:
步骤1:利用CFD-DEM耦合计算方法,模拟计算出浆液-土体双向耦合微观层面上的破坏范围,得到破坏范围的半径Rp;
步骤2:选取一种工程情况,将步骤1中的破坏范围Rp与选取的工程情况结合,导入相关软件中并将其描绘出包络图;
步骤3:将包络图曲线坐标导入EXCEL中进行再处理,处理后将包络图以数学表达式的形式表达出来;
步骤4:对其他注浆压力下的浆液扩散曲线形状进行观察,其形状与步骤2的包络图曲线基本一致的,可以假设其数学表达式与步骤3的数学表达式基本一致,则可依据步骤3的数学表达式列出方程;
步骤5:建立多个方程组对表达式,结合对应的工程情况,对方程组的系数进行求解,得出曲线方程。
优选的,在步骤1中,运用CFD-DEM流固耦合计算方法对注浆工程的全过程进行模拟。
优选的,步骤2中,相关软件为AUTOCAD。
相对于现有技术,本发明的有益效果为:
本发明通过前期对CFD-DEM耦合模拟分析了不同注浆压力下,浆液在土体中的运动方式,并从不同工况的图中发现出浆液的扩散规律,并推导到出了浆液在土体中的扩散的计算方法。本发明方法从前期所做的数值模拟工作出发,将多种软件进行综合应用分析,对浆液在土体中的扩散规律进行具体的数学表达式的推导,为注浆工程的进一步研究提供新的思路,并且可以将此公式应用到实际工程中,指导实际现场工作。
附图说明
图1为导入到AUTOCAD绘制的包络图。
图2为EXCEL中的包络曲线。
图3为包络曲线的前半部分。
图4为在曲线上选取的两个点。
图5为注浆工程加固路基浆液扩散范围公式推导方法流程图。
图6为其他注浆压力下的扩散曲线
具体实施方式
下面结合附图1~5和具体实施方式对本发明技术方案做进一步详细描述:
1、利用CFD-DEM耦合计算方法,模拟计算出浆液-土体双向耦合微观层面上的破坏范围,如图1所示。
根据软土地基注浆止沉的机理,本发明选用CFD-DEM耦合计算方法模拟分析不同注浆压力作用下,浆液与土体的双向作用效果,与其他模拟方式相比,此方法更加贴合注浆过程的实际情况,更能准确,模拟注浆的全过程,共计算7种不同注浆压力作用下的微观层面上,浆液与土体的相互影响情况。
2、选取注浆压力为0.4MPa时,将破坏范围导入相关软件中并将其描绘出包络图如图1所示。
选取注浆压力为p=0.4MPa时的计算结果导入到AUTOCAD中,并利用AUTOCAD的作图优势,将其浆液的影响范围包络图描绘出来,为后面的计算做出铺垫。
3、将包络图以数学表达式的形式表达出来。
将做出的包络图曲线坐标导入到EXCEL中进行再处理,如图2所示,观察到曲线可分为两部分,前半部分为对数曲线,如图3所示,后半部分为圆曲线。所以,对曲线坐标进行拟合得出此种工况下,浆液扩散曲线的拟合方程为
Figure BDA0002294393860000041
4、对其他注浆压力下的浆液扩散曲线数学表达式做出假设。
通过观察得出,其他注浆压力下的浆液扩散曲线大致相同,则对其他注浆压力下的浆液扩散曲线假设为:
Figure BDA0002294393860000042
5、建立多个方程组对表达式的系数进行求解,得出曲线方程。
注浆管分别取端点和连接处的点,假设其坐标为(x1,y1)、(x2,y2),如附图4所示。由于在前期工作中,得到关于扩散范围与注浆压力的关系式为R=64.53p1.331+b,并由曲线关系得到以下方程组:
Figure BDA0002294393860000051
分别将a、c、d、r解得:
Figure BDA0002294393860000052
由于在实际工程中,x1、y1是可以直接得出,x1即可由注浆管的方程计算得出,其计算方程为x2+y2=b2,y1即为注浆管上射浆孔的半径。x2、y2为未知量,故将x2、y2以其他已知量表示。计算得出:
Figure BDA0002294393860000053
其中
Figure BDA0002294393860000054
N=2[2(R-b)+(3R+b)(B-A)],A=ln x1,
Figure BDA0002294393860000055
故其得到的注浆影响范围曲线方程为:
Figure BDA0002294393860000061
其中x1为注浆管半径,y1为射浆孔半径,x2、y2为两曲线连接处的坐标,b为注浆管半径。R为土体微观层面接触键破坏范围半径。
参数验证:当p=0.3MPa时,b=4cm,R=16.996,y1=1,
Figure BDA0002294393860000062
A=ln x1=ln3.873=1.354,
Figure BDA0002294393860000063
B=ln x2=ln13.47=2.621
将以上结果带入到上述方程中,得到M=915.850,N=191.324,y2=2.187
Figure BDA0002294393860000064
Figure BDA0002294393860000065
则前半部分的数学表达式为y=0.905lnx-0.269
Figure BDA0002294393860000066
r=R-d=16.996-14.635=2.361
所以后半部分的数学表达式为(x-14.635)2+y2=2.3612,x>0,y>0。
将其他注浆压力下的计算参数同样带入可得如图6的曲线。
与数值模拟出的注浆扩散曲线基本一致。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何不经过创造性劳动想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求书所限定的保护范围为准。

Claims (3)

1.一种注浆工程加固路基浆液扩散范围公式推导方法,其特征在于,包括以下步骤:
步骤1:利用CFD-DEM耦合计算方法,模拟计算出浆液-土体双向耦合微观层面上的破坏范围,得到破坏范围的半径Rp;
步骤2:选取一种工程情况,将步骤1中的破坏范围Rp与选取的工程情况结合,导入制图软件中并将其描绘出包络图;
步骤3:将包络图曲线坐标导入EXCEL中进行再处理,处理后将包络图以数学表达式的形式表达出来;
步骤4:对其他注浆压力下的浆液扩散曲线形状进行观察,其形状与步骤2的包络图曲线基本一致的,可以假设其数学表达式与步骤3的数学表达式基本一致,则可依据步骤3的数学表达式列出方程;
步骤5:建立多个方程组,结合对应的工程情况,对方程组的系数进行求解,得出曲线方程。
2.根据权利要求1所述一种注浆工程加固路基浆液扩散范围公式推导方法,其特征在于,在步骤1中,运用CFD-DEM流固耦合计算方法对注浆工程的全过程进行模拟。
3.根据权利要求1所述一种注浆工程加固路基浆液扩散范围公式推导方法,其特征在于,步骤2中,制图软件为AUTOCAD。
CN201911194674.9A 2019-11-28 2019-11-28 一种注浆工程加固路基浆液扩散范围的确定方法 Expired - Fee Related CN110990924B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911194674.9A CN110990924B (zh) 2019-11-28 2019-11-28 一种注浆工程加固路基浆液扩散范围的确定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911194674.9A CN110990924B (zh) 2019-11-28 2019-11-28 一种注浆工程加固路基浆液扩散范围的确定方法

Publications (2)

Publication Number Publication Date
CN110990924A true CN110990924A (zh) 2020-04-10
CN110990924B CN110990924B (zh) 2021-02-26

Family

ID=70087925

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911194674.9A Expired - Fee Related CN110990924B (zh) 2019-11-28 2019-11-28 一种注浆工程加固路基浆液扩散范围的确定方法

Country Status (1)

Country Link
CN (1) CN110990924B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112347709A (zh) * 2020-10-21 2021-02-09 山东大学 一种基于dem-cfd耦合的渗透注浆过程模拟方法及系统
CN117034803A (zh) * 2023-08-08 2023-11-10 大秦铁路股份有限公司太原高铁工务段 一种黄土路基注浆盲区的计算方法、装置、设备及介质

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103034765A (zh) * 2012-12-14 2013-04-10 天津大学 基于数值模拟的采空区注浆动态全过程仿真方法
CN104895595A (zh) * 2015-04-18 2015-09-09 山东大学 一种确定锚杆注浆时注浆液扩散范围及规律的方法
CN105716998A (zh) * 2016-03-25 2016-06-29 葛洲坝集团试验检测有限公司 一种混凝土裂缝化学灌浆浆液扩散半径的计算方法
CN109085645A (zh) * 2018-09-19 2018-12-25 东北大学 一种基于微震监测的裂隙岩体劈裂注浆扩散范围的识别方法
CN109815573A (zh) * 2019-01-16 2019-05-28 中铁二院工程集团有限责任公司 一种基于流体颗粒耦合模拟高铁路基注浆微观机理的方法
CN110284902A (zh) * 2019-05-09 2019-09-27 中国铁路总公司 掌子面注浆参数设计方法
KR102037313B1 (ko) * 2019-01-30 2019-10-28 박태영 지압형 사방확산 앵커조립체

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103034765A (zh) * 2012-12-14 2013-04-10 天津大学 基于数值模拟的采空区注浆动态全过程仿真方法
CN104895595A (zh) * 2015-04-18 2015-09-09 山东大学 一种确定锚杆注浆时注浆液扩散范围及规律的方法
CN105716998A (zh) * 2016-03-25 2016-06-29 葛洲坝集团试验检测有限公司 一种混凝土裂缝化学灌浆浆液扩散半径的计算方法
CN109085645A (zh) * 2018-09-19 2018-12-25 东北大学 一种基于微震监测的裂隙岩体劈裂注浆扩散范围的识别方法
CN109815573A (zh) * 2019-01-16 2019-05-28 中铁二院工程集团有限责任公司 一种基于流体颗粒耦合模拟高铁路基注浆微观机理的方法
KR102037313B1 (ko) * 2019-01-30 2019-10-28 박태영 지압형 사방확산 앵커조립체
CN110284902A (zh) * 2019-05-09 2019-09-27 中国铁路总公司 掌子面注浆参数设计方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
刘 健,刘人太,张 霄,李术才: ""水泥浆液裂隙注浆扩散规律模型试验与数值模拟"", 《岩石力学与工程学报》 *
崔韬: ""基于离散元的注浆过程模拟研究"", 《山西建筑》 *
雷进生,刘非,彭刚,王乾峰,夏磊: ""考虑参数动态变化和相互关联的浆液扩散范围研究"", 《长江科学院院报》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112347709A (zh) * 2020-10-21 2021-02-09 山东大学 一种基于dem-cfd耦合的渗透注浆过程模拟方法及系统
CN117034803A (zh) * 2023-08-08 2023-11-10 大秦铁路股份有限公司太原高铁工务段 一种黄土路基注浆盲区的计算方法、装置、设备及介质
CN117034803B (zh) * 2023-08-08 2024-05-14 大秦铁路股份有限公司太原高铁工务段 一种黄土路基注浆盲区的计算方法、装置、设备及介质

Also Published As

Publication number Publication date
CN110990924B (zh) 2021-02-26

Similar Documents

Publication Publication Date Title
CN106934251B (zh) 一种类矩形盾构隧道施工中地表沉降计算方法
CN110990924A (zh) 一种注浆工程加固路基浆液扩散范围公式推导方法
CN103866769B (zh) 基于3d地质模型及实时监控大坝灌浆工程分析控制方法
CN1952921B (zh) 一种确定地基预留沉降量的方法
CN105426568B (zh) 一种估算喀斯特地区土壤流失量的方法
CN103425054B (zh) 一种基于数字化的电力隧道施工控制方法
CN103437381B (zh) 群桩拖带叠加沉降量评估处理方法
CN110702881B (zh) 岩土材料参数变异性结果的预测方法及其应用
CN103729880A (zh) 一种地质勘探矿体三维建模系统
CN103399970A (zh) 用离散裂缝线处理进行数模测算油藏流动状况的方法
CN110765677B (zh) 一种高精度、快速三维地质模型有限元模型的建模方法
CN107908901B (zh) 一种基坑施工中承压水降水对周边环境影响的计算方法
CN110210135A (zh) 一种边坡工程全区域质量评价技术
CN103047975B (zh) 准确快速测量放样系统的方法
CN107857517A (zh) 一种使用红砂岩的路基填料及其制作方法
CN116680934B (zh) 基于卸荷率的地层开挖引起下部隧道竖向变形计算方法
CN105005667B (zh) 劈裂注浆复合土体等效弹性模量及泊松比计算方法
CN106638508B (zh) 高地应力条件下大型洞室轴线选择方法
CN103473404B (zh) 一种基于断面测量数据的实测三心拱巷道实体建模方法
CN112861208B (zh) 一种排桩支护体系转化地下连续墙支护结构的分析方法
Guo et al. Explicit form of exact analytical solution for calculating ground displacement and stress induced by shallow tunneling and its application
CN103246778B (zh) 一种曲线拟合的处理方法
CN110909446B (zh) 一种基于事件点与裂缝扩展双场耦合的复杂缝网重构方法
CN103438860B (zh) 电力隧道盾构施工地面最大沉降值预测方法
CN105274921A (zh) 汽车试验场高速环道Mcconnell曲线路基成型方法

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
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20210226

Termination date: 20211128

CF01 Termination of patent right due to non-payment of annual fee