CN111476404A - 一种冻融损伤软岩长期变形的预测方法 - Google Patents
一种冻融损伤软岩长期变形的预测方法 Download PDFInfo
- Publication number
- CN111476404A CN111476404A CN202010189911.9A CN202010189911A CN111476404A CN 111476404 A CN111476404 A CN 111476404A CN 202010189911 A CN202010189911 A CN 202010189911A CN 111476404 A CN111476404 A CN 111476404A
- Authority
- CN
- China
- Prior art keywords
- freeze
- thaw
- rock
- deformation
- damage
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B21/00—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
- G01B21/32—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring the deformation in a solid
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Game Theory and Decision Science (AREA)
- Marketing (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Tourism & Hospitality (AREA)
- Entrepreneurship & Innovation (AREA)
- General Business, Economics & Management (AREA)
- Development Economics (AREA)
- Theoretical Computer Science (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
一种冻融损伤软岩长期变形的预测方法,对寒区岩土工程建设提供危岩体变形评判准则,保证基础工程全寿命周期使用的安全性。包括如下步骤:①构建冻融时效损伤变量多元函数;②结合冻融时效损伤变量多元函数,构建冻融时效损伤耦合蠕变力学本构方程;③对现场软岩进行标准规程采样,记录岩石试样的应变随应力及时间的变化值;④应变时间曲线与本构方程进行非线性拟合,既而确定方程中的各力学参数;⑤岩石力学参数转化为岩体力学参数;⑥通过室内外冻融循环等效方程确定等效冻融循环次数;⑦实测现场岩体应力场空间分布规律及位移增量;⑧预测该岩体在未来tf时刻所对应的应变值。
Description
本明涉及工程技术领域,特别涉及一种冻融损伤软岩长期变形的预测方法。
技术领域
一种冻融损伤软性岩体长期变形的预测方法。
背景技术
软岩是指在一定条件下具有明显塑性变形特征的岩体,以泥岩、页岩及粉砂岩等为主,其主要特征为强度低、孔隙度大、颗粒胶结性差、受风化作用显著、节理面发育及内含大量膨胀黏土矿物,在地壳表面以沉积岩、火成岩或变质岩的形式大量赋存。软岩的特殊性质使其对外界作用的敏感性较硬岩更高,就目前我国对西部地区进行的基础设施建设,面临了大量的寒区冻融岩体,其中以软岩为主受冻融作用较为显著,需以此作为载力体,其时效变形机制的研究及全寿命变形量预测就显得尤为重要。在实际工程中,软岩所承受的外界作用不是单一的循环幅值叠加,而是多场耦合条件下的复合结果,故寒区软岩时效劣化特征就应为冻涨、溶解、长期受荷等叠加作用下的直观反映,首先软岩受冻融损伤显著,具体表现为内部微观颗粒(如大量赋存的膨胀黏土矿物)对冻融作用的敏感性较强,连续的冻融疲劳加卸载破坏了岩体的原生结构,或扩展了岩体的微观缺陷,宏观上岩体的力学性能弱化,其次现实中软岩长期处在一定的应力场条件下,在低于屈服强度的条件下产生以弹性为主的变形后趋于稳定,在高于屈服强度条件下产生以塑性为主的变形后破坏,综合考虑在冻融时效损伤耦合作用下软岩的变形特征及量化反映对于寒区工程建设就具有指导性意义。
近年来国内外学者在试验层面及损伤理论对寒区冻融软岩变形劣化机制进行了许多研究,但对于实际寒区工程中运用的准则及方法较少且不全面,且多数仅限于考虑冻融后软岩瞬时力学行为,而现实中软岩在外荷载作用下经历长期变形到破坏的现象更为广泛;同时,对于构建软岩宏观损伤力学本构模型的方法较少,多数也仅为软岩的常规力学模型,如申请号201510387574.3《一种适用于加卸载全过程的软岩蠕变本构模型的建立方法》通过引入Kelvin体与弹簧元件构建了描述软岩流变力学特征的本构模型,该模型可定量分析出软岩衰减蠕变及卸载后含残余变形的本构关系,但由于模型的构建过程中没有考虑岩石流变的时效受荷损伤特性,其参数也为定常性,使得在一定条件下该模型具有局限性,基于此,在考虑流变的时效损伤同时引入外界损伤因子(如冻融损伤),岩石长期变形为冻融时效损伤耦合条件下的综合体现,且将模型参数转换为岩体力学参数使其更能接近于现场损伤岩体的变形特征,以此作为变形预测指标更为可靠。
发明内容
本发明所要解决的技术问题是提供一种冻融损伤软岩长期变形的预测方法,对寒区岩土工程建设提供危岩体变形评判准则,保证基础工程全寿命周期使用的安全性。
本发明解决上述技术问题所采取的技术方案如下:
本发明的一种冻融损伤软性岩体长期变形的预测方法,包括如下步骤:
①构建冻融时效损伤变量多元函数,包含冻融损伤与时效性疲劳损伤,其中冻融损伤以随冻融循环次数变化的非定常性软岩力学参数来量化,而时效性损伤以统计损伤随机分布函数来定义其长期受荷损伤累积效应,具体表达为:
其中,D(T)为冻融损伤变量,P(T)为冻融循环次数为T时的宏观力学参数,P0为冻融循环次数为0时的宏观力学参数,D(t)为时效性损伤变量,w(t)为损伤密度函数,t为时间,θ与λ为长期受荷损伤参数,e为自然常数;
②结合冻融时效损伤变量多元函数,构建冻融时效损伤耦合蠕变力学本构方程,即:
式中,ε为总变形位移;σ为总应力;EM、EK为弹性模量;ηK、ηM为黏滞系数,DE(T)为弹性冻融损伤变量,Dη(T)为黏性冻融损伤变量;
③对现场软岩进行标准规程采样,将完整性较好且声速波动较小的岩样放入室内进行冻融循环试验及不同冻融循环后的三轴蠕变试验,并记录岩石试样的应变随应力及时间的变化值;
④将步骤③试验所得出的应变时间曲线与步骤②中的本构方程进行非线性拟合,既而确定方程中的各力学参数;
⑤将步骤④确定的岩石力学参数转化为岩体力学参数,主要基于岩体的结构面、风化程度与水文地质条件来定义转化函数,确定出来的包含岩体的应力、弹性模量与黏滞系数,具体为:
式中,σm为岩体应力;Eck、cm为岩体弹性模量;ηck、cm为岩体黏滞系数; mi为完整岩石常数;k值为岩体结构面系数;d为扰动系数;f(mi,k,d)、 g(k,d)与h(k)岩体质量综合评价函数,具体为:
⑥通过室内外冻融循环等效方程确定等效冻融循环次数,等效方程为:
其中,Te为实验室等效冻融循环次数;TF为实验室冻融循环破坏寿命; Ti为现场单次气候特征变化的循环加载次数;TFi为现场单次气候特征变化的循环加载对应的疲劳寿命;
⑦实测现场岩体应力场空间分布规律及位移增量,确定岩体初始应力及应变率,通过对步骤2冻融时效损伤耦合本构方程求关于时间的一阶偏导与实测位移增量进行误差拟合,得出现状岩体历史变形时长,即:
⑧预测该岩体在未来tf时刻所对应的应变值,将步骤⑤确定的现场岩体力学参数、步骤⑥确定的等效冻融循环次数及步骤⑦实测的初始应力带入步骤②中的本构方程构建现场岩体目标变形函数ε(t),线性叠加历史变形时长tm及预测变形时长tf代入ε(t),ε(tm+tf)即为现场岩体在未来tf时刻的变形量。
本发明的有益效果是,力学模型的构建考虑了寒区岩体冻融时效损伤耦合效通过一系列等效方程将力学本构运用于计算预测寒区危岩体变形,对寒区岩土工程建设提供危岩体变形评判准则,可有效保证基础工程全寿命周期使用的安全性。
附图说明
本说明书包括如下五幅附图
图1为流程图;
图2为损伤变量拟合图;
图3a、图3b为力学模型拟合图;
图4为现场实测位移数据图;
图5为预测岩体变形时间曲线图。
具体实施方式
下面结合具体工程实例对本发明进行更为清楚且全面的阐述。
以四川省某国道K14+250沿线冻灾区边坡危岩体变形预测为例,选取当地三叠系页岩进行室内试验及力学模型建立,按如下步骤进行:
(1)取样:取完整度较好的岩体按照相关规定制成高100mm直径 50mm的岩石圆柱体,并通过声波仪测定每组试样的声速,选取声速相差不大的岩样进行试验;
(2)冻融循环试验:将试样放入全自动冻融测试仪内进行冻融循环试验,根据当地气象资料确定冻融温差-10℃~30℃,冻结时长6h,溶解时长6h,冻融循环次数以0、15、30及45次为基准;
(3)冻融后三轴流变试验:将经历不同冻融循环次数后的岩样放入三轴流变试验仪内进行三轴蠕变试验,围压σ2=σ3选取0、5及10MPa,轴向应力σ1以该岩石强度的百分七十五进行分级加载,试验过程中记录每岩样应变随应力与时间的变化值;
(4)将构建的冻融时效损伤本构方程与试验实测数据进行最小二乘法拟合,确定模型弹性模量、粘滞系数、冻融损伤特征变量及长期受荷时效损伤参数,拟合图见附图2,拟合参数如下:
(5)根据现场岩体质量等级、表面风化与结构面发育程度等综合工程地质特征指标对岩石力学参数进行折减,使其过渡为岩体力学参数,转化后的参数如下所示:
对上述冻融时效损伤力学模型参数进行非线性拟合(见图3),构建自变量为冻融循环次数的弹性及黏性损伤变量函数,并引入室内外冻融循环等效方程确定等效冻融循环次数T,将T带入损伤变量函数确定DE(T)与 Dη(T);
(6)现场试验工作需确定岩体所处应力场大小及岩体所经历的历史变形时长,进而可得出寒区危岩体的长期变形方程,首先应力场大小根据现场应力传感器进行实测,其次历史变形时长求解需通过现场应变传感器记录岩体微变形量ΔεΔε1、Δε2、Δε3、…..Δεn(见附图4),并对上述所推导的本构方程求关于时间的一阶偏导建立起等效时差函数在收集n个数据的基础上对等效时差函数求关于实测位移量的离散型方差,并求其极小值sm,极值对应的自变量tm即为历史变形时;
将冻融时效损伤力学模型参数包括E、η、λ、θ、D(T)与σ带入冻融时效损伤力学本构方程建立危岩体长期变形方程,函数图见5,具体方程如下:
ε(t)=0.0234+(1-e-1.323t)×0.59812+1.3756×t×e0.521t
(7)量化该路段危岩体在其使用寿命期间的变形演化特征,对此进行风险性评估,取tf为道路服务年限,将tf+tm带入上述方程求解最终变形量,从图5可以看出,该危岩体的变形量在141.3mm,在服务年限后段其失稳概率较高,属于先期缓慢发育损伤累计到后期变形突增型的岩体,需对此进行相应的支护措施。
Claims (5)
1.一种冻融损伤软性岩体长期变形的预测方法,包括如下步骤:
①构建冻融时效损伤变量多元函数,包含冻融损伤与时效性疲劳损伤,其中冻融损伤以随冻融循环次数变化的非定常性软岩力学参数来量化,而时效性损伤以统计损伤随机分布函数来定义其长期受荷损伤累积效应,具体表达为:
其中,D(T)为冻融损伤变量,P(T)为冻融循环次数为T时的宏观力学参数,P0为冻融循环次数为0时的宏观力学参数,D(t)为时效性损伤变量,w(t)为损伤密度函数,t为时间,θ与λ为长期受荷损伤参数,e为自然常数;
②结合冻融时效损伤变量多元函数,构建冻融时效损伤耦合蠕变力学本构方程,即:
式中,ε为总变形位移;σ为总应力;EM、EK为弹性模量;ηK、ηM为黏滞系数,DE(T)为弹性冻融损伤变量,Dη(T)为黏性冻融损伤变量;
③对现场软岩进行标准规程采样,将完整性较好且声速波动较小的岩样放入室内进行冻融循环试验及不同冻融循环后的三轴蠕变试验,并记录岩石试样的应变随应力及时间的变化值;
④将步骤③试验所得出的应变时间曲线与步骤②中的本构方程进行非线性拟合,既而确定方程中的各力学参数;
⑤将步骤④确定的岩石力学参数转化为岩体力学参数,主要基于岩体的结构面、风化程度与水文地质条件来定义转化函数,确定出来的包含岩体的应力、弹性模量与黏滞系数,具体为:
式中,σm为岩体应力;Eck、cm为岩体弹性模量;ηck、cm为岩体黏滞系数;mi为完整岩石常数;k值为岩体结构面系数;d为扰动系数;f(mi,k,d)、g(k,d)与h(k)岩体质量综合评价函数,具体为:
⑥通过室内外冻融循环等效方程确定等效冻融循环次数,等效方程为:
其中,Te为实验室等效冻融循环次数;TF为实验室冻融循环破坏寿命;Ti为现场单次气候特征变化的循环加载次数;TFi为现场单次气候特征变化的循环加载对应的疲劳寿命;
⑦实测现场岩体应力场空间分布规律及位移增量,确定岩体初始应力及应变率,通过对步骤2冻融时效损伤耦合本构方程求关于时间的一阶偏导与实测位移增量进行误差拟合,得出现状岩体历史变形时长,即:
⑧预测该岩体在未来tf时刻所对应的应变值,将步骤⑤确定的现场岩体力学参数、步骤⑥确定的等效冻融循环次数及步骤⑦实测的初始应力带入步骤②中的本构方程构建现场岩体目标变形函数ε(t),线性叠加历史变形时长tm及预测变形时长tf代入ε(t),ε(tm+tf)即为现场岩体在未来tf时刻的变形量。
3.如权利要求1所述的一种冻融损伤软性岩体长期变形的预测方法,其特征是:所述步骤①中冻融循环试验要求目标试样冻融循环温差与时长应与现场气象资料一致,偏差不高于1℃与20min。
4.如权利要求书1所述的一种冻融损伤软性岩体长期变形的预测方法,其特征在于:所述步骤①中冻融循环后软岩三轴蠕变试验围压选取应与现场实测地应力场分布状况一致,偏差不大于0.1MPa。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010189911.9A CN111476404B (zh) | 2020-03-18 | 2020-03-18 | 一种冻融损伤软岩长期变形的预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010189911.9A CN111476404B (zh) | 2020-03-18 | 2020-03-18 | 一种冻融损伤软岩长期变形的预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111476404A true CN111476404A (zh) | 2020-07-31 |
CN111476404B CN111476404B (zh) | 2023-01-31 |
Family
ID=71748402
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010189911.9A Active CN111476404B (zh) | 2020-03-18 | 2020-03-18 | 一种冻融损伤软岩长期变形的预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111476404B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111307874A (zh) * | 2020-04-09 | 2020-06-19 | 河海大学 | 一种气候环境影响的混凝土电阻率预测方法 |
CN112597638A (zh) * | 2020-12-11 | 2021-04-02 | 浙江省工程勘察设计院集团有限公司 | 一种基于钻孔技术的岩芯冻融循环温度衰减模型建立方法 |
CN112765781A (zh) * | 2020-12-30 | 2021-05-07 | 中铁二院工程集团有限责任公司 | 层状岩体破坏准则及确定方法 |
CN113063810A (zh) * | 2021-03-22 | 2021-07-02 | 西安科技大学 | 获得砂岩冻融作用下宏细观损伤演化规律的方法 |
CN113533410A (zh) * | 2021-07-09 | 2021-10-22 | 长沙理工大学 | 一种冻融循环下路基土永久变形的预估方法 |
CN114486532A (zh) * | 2022-03-28 | 2022-05-13 | 中国矿业大学(北京) | 测试含冰软岩样品蠕变中流体运移规律的装置和方法 |
WO2023029793A1 (zh) * | 2021-09-03 | 2023-03-09 | 中铁一局集团第二工程有限公司 | 一种冻融循环作用下岩石蠕变损伤本构模型的构建方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003055204A (ja) * | 2001-08-21 | 2003-02-26 | Sansho Iyaku Kk | ユビデカレノン組成物及び同組成物を内容物とするカプセル |
CN107515291A (zh) * | 2017-08-23 | 2017-12-26 | 西安科技大学 | 一种等围压作用下岩石冻融损伤本构模型的构建方法 |
CN108829916A (zh) * | 2018-04-25 | 2018-11-16 | 中铁二院工程集团有限责任公司 | 硬岩冻融损伤长期变形模型的构建方法 |
CN109283215A (zh) * | 2018-11-26 | 2019-01-29 | 中国地质大学(北京) | 反复冻融损伤后寒区隧道围岩冻胀力的测算方法 |
-
2020
- 2020-03-18 CN CN202010189911.9A patent/CN111476404B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003055204A (ja) * | 2001-08-21 | 2003-02-26 | Sansho Iyaku Kk | ユビデカレノン組成物及び同組成物を内容物とするカプセル |
CN107515291A (zh) * | 2017-08-23 | 2017-12-26 | 西安科技大学 | 一种等围压作用下岩石冻融损伤本构模型的构建方法 |
CN108829916A (zh) * | 2018-04-25 | 2018-11-16 | 中铁二院工程集团有限责任公司 | 硬岩冻融损伤长期变形模型的构建方法 |
CN109283215A (zh) * | 2018-11-26 | 2019-01-29 | 中国地质大学(北京) | 反复冻融损伤后寒区隧道围岩冻胀力的测算方法 |
Non-Patent Citations (3)
Title |
---|
刘新喜等: "冻融循环作用下岩质页岩蠕变模型研究", 《中国公路学报》 * |
陈国庆等: "冻融后石英砂岩三轴蠕变特性实验研究", 《岩土力学》 * |
陈国庆等: "冻融循环作用下砂岩蠕变特性及损伤模型研究", 《工程地质学报》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111307874A (zh) * | 2020-04-09 | 2020-06-19 | 河海大学 | 一种气候环境影响的混凝土电阻率预测方法 |
CN112597638A (zh) * | 2020-12-11 | 2021-04-02 | 浙江省工程勘察设计院集团有限公司 | 一种基于钻孔技术的岩芯冻融循环温度衰减模型建立方法 |
CN112765781A (zh) * | 2020-12-30 | 2021-05-07 | 中铁二院工程集团有限责任公司 | 层状岩体破坏准则及确定方法 |
CN112765781B (zh) * | 2020-12-30 | 2022-07-29 | 中铁二院工程集团有限责任公司 | 层状岩体破坏准则及确定方法 |
CN113063810A (zh) * | 2021-03-22 | 2021-07-02 | 西安科技大学 | 获得砂岩冻融作用下宏细观损伤演化规律的方法 |
CN113533410A (zh) * | 2021-07-09 | 2021-10-22 | 长沙理工大学 | 一种冻融循环下路基土永久变形的预估方法 |
WO2023029793A1 (zh) * | 2021-09-03 | 2023-03-09 | 中铁一局集团第二工程有限公司 | 一种冻融循环作用下岩石蠕变损伤本构模型的构建方法 |
CN114486532A (zh) * | 2022-03-28 | 2022-05-13 | 中国矿业大学(北京) | 测试含冰软岩样品蠕变中流体运移规律的装置和方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111476404B (zh) | 2023-01-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111476404B (zh) | 一种冻融损伤软岩长期变形的预测方法 | |
Liu et al. | Modeling time-dependent deformation behavior of brittle rock using grain-based stress corrosion method | |
CN110532694B (zh) | 一种基于地质力学模型综合法试验的拱坝安全评价方法 | |
CN105259035A (zh) | 一种岩石材料时效-弹塑性力学本构模型的构建方法 | |
CN110346213B (zh) | 一种隧道下穿采空区围岩稳定性评价方法 | |
CN101592652A (zh) | 湿热地区沥青路面疲劳寿命预估方法 | |
CN105046051B (zh) | 一种实际环境冻融作用次数预估方法 | |
CN111209695A (zh) | 一种基于lstm的结构动力响应预测方法 | |
CN101792992A (zh) | 永久性路面设计方法 | |
CN112504851A (zh) | 一种考虑谷幅变形作用的拱坝变形监控方法 | |
Cescatti et al. | Analysis and evaluations of flat jack test on a wide existing masonry buildings sample | |
CN117309060A (zh) | 一种基于云计算的建筑幕墙结构性能监测系统 | |
Ma et al. | Excess pore water pressure generation in saturated sandy soils subjected to various cyclic stress paths | |
Li et al. | Cyclic freezing-thawing induced rock strength degradation, crack evolution, heave and settlement accounted for by a DEM model | |
CN109117593B (zh) | 锚索预应力损失与岩体蠕变耦合分析方法 | |
CN104196024B (zh) | 基于锚索预应力损失与岩土体蠕变耦合的预应力损失计算方法 | |
CN107894437B (zh) | 严寒地区混凝土伸缩缝施工水化历程监测方法 | |
CN110197015B (zh) | 一种坝基预应力锚索有效拉应力测定方法 | |
CN102660967B (zh) | 寒区单桩经验流变预报方程的确定方法 | |
CN112362502A (zh) | 一种隧道锚碇的安全性评估方法 | |
CN113791107B (zh) | 一种尾矿充填体监测方法 | |
Wang et al. | Research on Prediction Method of Foundation Settlement | |
Olek | Some remarks on the pore water pressure dissipation patterns from the one-dimensional consolidation test | |
Feddema et al. | Lateral loading of pile foundations due to embankment construction | |
Lintzén et al. | Study on basic material properties of artificial snow |
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 |