CN110987662A - 一种考虑冻融循环损伤的节理剪切强度的确定方法 - Google Patents
一种考虑冻融循环损伤的节理剪切强度的确定方法 Download PDFInfo
- Publication number
- CN110987662A CN110987662A CN201911167205.8A CN201911167205A CN110987662A CN 110987662 A CN110987662 A CN 110987662A CN 201911167205 A CN201911167205 A CN 201911167205A CN 110987662 A CN110987662 A CN 110987662A
- Authority
- CN
- China
- Prior art keywords
- damage
- freeze
- joint
- shear strength
- thaw
- 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
- 230000006378 damage Effects 0.000 title claims abstract description 168
- 238000000034 method Methods 0.000 title claims abstract description 49
- 239000011435 rock Substances 0.000 claims abstract description 27
- 238000004364 calculation method Methods 0.000 claims abstract description 21
- 238000004891 communication Methods 0.000 claims abstract description 18
- 230000015556 catabolic process Effects 0.000 claims abstract description 5
- 238000006731 degradation reaction Methods 0.000 claims abstract description 5
- 238000010257 thawing Methods 0.000 claims description 39
- 238000009826 distribution Methods 0.000 claims description 23
- 230000008569 process Effects 0.000 claims description 19
- 230000008859 change Effects 0.000 claims description 17
- 230000008014 freezing Effects 0.000 claims description 17
- 238000007710 freezing Methods 0.000 claims description 17
- 230000009471 action Effects 0.000 claims description 16
- 238000011160 research Methods 0.000 claims description 12
- 230000003902 lesion Effects 0.000 claims description 11
- 208000027418 Wounds and injury Diseases 0.000 claims description 6
- 208000014674 injury Diseases 0.000 claims description 6
- 238000006243 chemical reaction Methods 0.000 claims description 5
- 230000003750 conditioning effect Effects 0.000 claims description 4
- 125000004122 cyclic group Chemical group 0.000 claims description 4
- 239000000463 material Substances 0.000 claims description 4
- 238000012544 monitoring process Methods 0.000 claims description 4
- 238000012424 Freeze-thaw process Methods 0.000 claims description 3
- 230000001351 cycling effect Effects 0.000 claims description 3
- 238000005315 distribution function Methods 0.000 claims description 3
- 229910052500 inorganic mineral Inorganic materials 0.000 claims description 3
- 239000011707 mineral Substances 0.000 claims description 3
- 239000005445 natural material Substances 0.000 claims description 3
- 230000003068 static effect Effects 0.000 claims description 3
- 238000006467 substitution reaction Methods 0.000 claims description 3
- 150000001875 compounds Chemical class 0.000 claims description 2
- 238000012360 testing method Methods 0.000 abstract description 21
- 238000009825 accumulation Methods 0.000 abstract description 2
- 230000001808 coupling effect Effects 0.000 abstract description 2
- 230000006872 improvement Effects 0.000 abstract description 2
- 210000001503 joint Anatomy 0.000 description 18
- 238000010008 shearing Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 230000007246 mechanism Effects 0.000 description 3
- 230000006866 deterioration Effects 0.000 description 2
- 239000004033 plastic Substances 0.000 description 2
- 229920003023 plastic Polymers 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000011083 cement mortar Substances 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000005336 cracking Methods 0.000 description 1
- 238000002425 crystallisation Methods 0.000 description 1
- 230000008025 crystallization Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 229910000514 dolomite Inorganic materials 0.000 description 1
- 239000010459 dolomite Substances 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000001939 inductive effect Effects 0.000 description 1
- 238000013101 initial test Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 150000003839 salts Chemical class 0.000 description 1
- 229920006395 saturated elastomer Polymers 0.000 description 1
- 230000035939 shock Effects 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N3/24—Investigating strength properties of solid materials by application of mechanical stress by applying steady shearing forces
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N3/02—Details
- G01N3/06—Special adaptations of indicating or recording means
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/02—Details not specific for a particular testing method
- G01N2203/0202—Control of the test
- G01N2203/0212—Theories, calculations
- G01N2203/0218—Calculations based on experimental data
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/02—Details not specific for a particular testing method
- G01N2203/026—Specifications of the specimen
- G01N2203/0298—Manufacturing or preparing specimens
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/02—Details not specific for a particular testing method
- G01N2203/06—Indicating or recording means; Sensing means
- G01N2203/067—Parameter measured for estimating the property
- G01N2203/0676—Force, weight, load, energy, speed or acceleration
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开一种考虑冻融循环损伤的节理剪切强度的确定方法,针对寒区岩体节理的复杂性,本发明方法深入分析了节理连通率与剪切强度的关系,基于室内试验提出了剪切强度参数随节理连通率变化的计算公式,便于精确计算复杂节理的剪切强度。本发明基于节理试样的冻融试验和直剪试验,建立的剪切强度参数损伤演化模型,进一步推导出改进的剪切强度计算公式,此计算公式同时考虑了冻融循环的损伤累积和节理连通率的影响,可深入分析冻融循环与节理连通率耦合作用对剪切强度参数的劣化特性,根据后续试验结果,得出改进计算公式的计算结果与实验数据较高的吻合程度,表明了改进准则的合理性。
Description
技术领域
本发明属于工程技术领域,具体涉及一种考虑冻融循环损伤的节理剪切强度的确定方法。
背景技术
天然岩体中赋存的节理在冻融循环作用下引起岩石微观损伤、冻胀碎裂以及破坏失稳等劣化现象,被认为是诱发寒区地质灾害的主要原因之一。在温度交替变化作用下,节理内蕴藏的水分会经历冻融过程,冻融产生的周期性冻胀力引起节理裂隙不断扩展,最终导致岩体失稳。目前,许多学者对此展开了大量的理论探索及试验研究。一方面是考虑冻融作用对岩石宏观力学性质影响,如:M.H.Ghobadi等在相关的研究中讨论了冻融循环作用岩石各项物理力学参数特性的影响(Experimental Studies on the Effects of CyclicFreezing–Thawing,Salt Crystallization,and Thermal Shock on the Physical andMechanical Characteristics of Selected Sandstones,Rock Mechanics&RockEngineering);吴安杰等在相关的研究中分别分析了岩石的各项力学参数与冻融损伤的关系(冻融循环作用下泥质白云岩力学特性及损伤演化规律研究.岩土力学)。另一方面是从微观及理论出发,将传统的损伤力学引入冻融损伤领域,研究岩石冻融损伤机制及其表征方法,如:贾海梁等根据不同的数学模型建立了岩石冻融损伤模型,并通过后期试验验证了所提出模型的的适用性(冻融循环作用下饱和砂岩损伤扩展模型研究.岩石力学与工程学报);Lai等将损伤理论方法应用于冻融损伤场中,分析了冻融损伤与岩石试样应变软化的相关规律(Yield criterion and elasto-plastic damage constitutive model forfrozen sandy soil,International Journal of Plasticity)。刘泉声等研究了冻融条件下岩体力学性能的演化并得出了相关的损伤模型(岩体冻融疲劳损伤模型与评价指标研究.岩石力学与工程学报)。
虽然目前有关于岩石冻融损伤机制的研究,但是关于不同连通率节理在冻融环境下力学特性的劣化,尤其是节理抗剪力学特性损伤演化机制的研究成果较少,而节理抗剪力学性能对岩体工程的稳定性往往起至关重要的作用。
发明内容
本发明的目的是提供一种考虑冻融循环损伤的节理剪切强度的确定方法,解决现有寒区岩土工程多节理且存在冻融循环损伤的条件,岩体剪切强度难以预测的技术难点。
本发明这种考虑冻融循环损伤的节理剪切强度的确定方法,包括以下步骤:
S1、根据节理试样的连通率变化与剪切强度参数之间的关系,构建计算公式;
S2、根据Lemaitre提出的经典损伤公式,定义黏聚力和内摩擦角损伤变量Dc和并计算得出不同连通率节理试样经历冻融后的黏聚力和内摩擦角损伤变量Dc和然后以冻融循环前强度为初始强度状态,以冻融循环后为损伤状态,结合损伤力学理论和Mohr-Coulomb准则,得出冻融循环后节理试样剪切强度τn方程;
S4、根据实验监测数据,拟合得出节理内冻胀损伤力F(n)与冻融循环次数的关系表达式,然后根据F(n)的表达式,转化S3步骤中黏聚力和内摩擦角损伤变量Dc和的表达式,再然后基于S2步骤中计算得出不同连通率节理试样经历冻融后的剪切强度参数损伤变量Dc和确定转化后表达式中模型参数a、和的计算表达式;
所述步骤S1中,构建计算公式的具体方法,包括以下步骤:分析节理试样未冻融时的剪切强度参数分布特征,可得剪切强度参数随连通率变化情况,进行拟合得到剪切强度参数与连通率的计算公式;
所述步骤S2的具体方法为:
S2.1分析节理试样经历不同冻融循环后抗剪强度劣化的劣化程度,得出抗剪强度参数的损伤变量随节理冻融循环累积的损伤规律;根据Lemaitre提出的经典损伤公式,黏聚力和内摩擦角损伤变量可以定义为:
S2.3再基于损伤力学理论和Mohr-Coulomb准则,并以冻融循环前强度为初始强度状态,以冻融循环后为损伤状态,冻融循环后节理试样剪切强度方程为:
cn=c0(1-Dc)(4a)
S3.1将岩石试样视为一种天然材料,其矿物组成和胶结作用各不相同,则可假定节理试样由不同材料强度的微单元组成,总微元数为N;Weibull分布函数适用于微元损伤,因此假定节理试样微元强度服从Weibull分布,其损伤概率密度函数可用损伤力F表示为:
S3.2在冻融循环作用下,节理试样的损伤是由局部微元体不均匀变形破坏引起的,假设在任意损伤力区间[F,F+dF]内产生的微元数为NP(x)dx,当损伤力达到某一水平F时,则已经损伤破坏的微元总数为:
式中,Nf为破坏的微单元数量,N为微单元的总数;
S3.3由于节理试样在冻融循环作用下,节理剪切强度会出现不同程度的损伤劣化,而剪切强度由剪切强度参数表征,剪切强度参数也会呈现损伤劣化,故则粘聚力和内摩擦角的损伤变量Dc和可定义为Nf与N之比,因此,联立式(5)和(6),则Dc和的方程式可分别表示为:
由此得出的方程表示损伤变量Dc和与损伤力F之间的本构关系,但仅适用于不考虑多次循环,冻融损伤力F在整个过程中不变时,损伤本构方程才是静态有效的;对于循环冻融损伤的整个变化动态过程,损伤状态变量Dc和应与多次冻融循环的完整变化过程相结合。
所述S4中步骤:得出冻胀损伤力F(n)与冻融循环次数的关系表达式的具体步骤为:
根据多位学者的试验研究发现,冻融循环条件下,节理试样主要在冻胀力作用下出现宏观裂纹并逐渐扩展,同时试样剪切强度呈现明显劣化,故冻融循环过程中冻胀力是导致的节理强度损伤的主要因素;因此,可将冻融循环过程中节理内的损伤力F和冻胀力F(n)视作近似等效;故采用进口薄膜压力传感器对整个冻融循环过程中冻胀力进行监测,分别记录下冻融循环第1次、第10次、第20次、第30次、第40次节理内最大冻胀力;
传感器记录下的各冻融循环次数下最大冻胀力,根据实验数据显示冻胀力随循环次数增加呈指数函数下降趋势,通过拟合可得出节理内冻胀力F(n)的表达式为:
式中:a为冻胀力随循环次数增加的变化系数,n为冻融循环次数,n≥1。
对上述分布参数进行拟合,(a/Fc)、mc和随节理连通率增长均满足指数型函数变化规律,采用Exponential函数对分布参数与节理连通率的关系进行归一化处理,得出(a/Fc)、和mc、与节理连通率的关系式为:
(a/Fc)=1.897+1.623*exp(-λ/0.181) (10a)
mc=0.210+0.568*exp(-λ/0.266) (10b)
所述S5步骤中,确定冻融剪切强度关系的具体方法为:
根据S4步骤拟合得出的(a/Fc)、mc和公式,并结合式S2步骤中的式(3)和S4步骤中式(9),可推算出节理冻融后剪切强度参数,并代入Mohr-Coulomb准则中,可得出冻融循环条件下的不同节理剪切破坏强度:
本发明的有益效果:1、针对寒区岩体节理的复杂性,本发明方法深入分析了节理连通率与剪切强度的关系,基于室内试验提出了剪切强度参数随节理连通率变化的计算公式,便于精确计算复杂节理的剪切强度。2、本发明基于节理试样的冻融试验和直剪试验,建立的剪切强度参数损伤演化模型,进一步推导出改进的剪切强度计算公式,此计算公式同时考虑了冻融循环的损伤累积和节理连通率的影响,可深入分析冻融循环与节理连通率耦合作用对剪切强度参数的劣化特性,根据后续试验结果,得出改进计算公式的计算结果与实验数据较高的吻合程度,表明了改进准则的合理性。
附图说明
图1为实施例S1步骤中节理试样模型及剪切示意图;其中,(a)节理立体图,(b)节理平面图。
图2为实施例S1步骤中节理试样剪切强度参数随连通率增加的变化规律;其中,(a)粘聚力随连通率变化关系,(b)内摩擦角随连通率变化关系。
图3为实施例S4步骤中λ=40%试样冻融循环过程中冻胀力监测示意图。
图4为实施例S4步骤中不同连通率节理中最大冻胀力随冻融循环次数增加的变化曲线。
图5为实施例S4步骤中节理粘聚力和内摩擦角冻融损伤演化模型拟合曲线;其中,(a)粘聚力损伤变量变化曲线,(b)内摩擦角损伤变量变化曲线,
图7为实施例S5步骤λ=50%节理试样冻融前后剪切强度实验值与解析值对比效果。
具体实施方式
下面对本发明的具体实施方式进行描述,以便于本技术领域的技术人员理解本发明,但应该清楚,本发明不限于具体实施方式的范围,对本技术领域的普通技术人员来讲,只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内,这些变化是显而易见的,一切利用本发明构思的发明创造均在保护之列。
实施例
根据本申请的一个实施例,根据冻融试验和直剪试验结果,推导得出考虑冻融循环损伤的节理剪切强度准则,包括:
S1、根据节理试样的连通率变化与剪切强度参数之间的关系,构建计算公式;
S2、根据Lemaitre提出的经典损伤公式,定义黏聚力和内摩擦角损伤变量Dc和并计算得出不同连通率节理试样经历冻融后的黏聚力和内摩擦角损伤变量Dc和然后以冻融循环前强度为初始强度状态,以冻融循环后为损伤状态,结合损伤力学理论和Mohr-Coulomb准则,得出冻融循环后节理试样剪切强度τn方程;
S4、根据实验监测数据,拟合得出节理内冻胀损伤力F(n)与冻融循环次数的关系表达式,然后根据F(n)的表达式,转化S3步骤中黏聚力和内摩擦角损伤变量Dc和的表达式,再然后基于S2步骤中计算得出不同连通率节理试样经历冻融后的剪切强度参数损伤变量Dc和确定转化后表达式中模型参数a、和的计算表达式;
以下对上述步骤进行详细说明:
S1、选用类岩石材料-水泥砂浆进行制样,试样长、宽、高分别为100mm、100mm、100mm,分别设置节理连通率λ为:0(完整试样)、20%、30%、40%,即节理尺寸为70mm×20mm×3mm、70mm×30mm×3mm、70mm×40mm×3mm,试样示例及剪切原理见图1。选择冻融循环次数n和结构面连通率λ为试验变量,进行2个因素的全面试验。冻融循环次数n分别设置为0次、10次、20次、30次和40次,待完成所有节理试样的冻融循环试验后,再进行直剪试验,记录其剪应力-剪切位移曲线,利用最小二乘法拟合得到峰值剪切强度与正应力关系,结合Mohr-Coulomb准则计算节理试样抗剪强度参数,则不同连通率节理在冻融条件下的粘聚力c和内摩擦角结果如表1。
分析节理试样未冻融时的剪切强度参数分布特征,可得剪切强度参数随连通率增大而减小,且二者之间的关系呈现显著的线性特征,如图2。进行线性拟合得到剪切强度参数与连通率的关系式为:
cλ=c0+cλ (1a)
φλ=φ0+dλ (1b)
表1节理试样冻融前后剪切强度参数
S2、分析节理试样经历不同冻融循环后抗剪强度劣化的劣化程度,得出抗剪强度参数的损伤变量随节理冻融循环累积的损伤规律。根据Lemaitre提出的经典损伤公式,黏聚力和内摩擦角损伤变量可以定义为:
根据式(2)可计算得出不同连通率节理试样经历冻融后的抗剪强度参数损伤变量Dc和结果如表2所示。再基于损伤力学理论和Mohr-Coulomb准则,并以冻融循环前强度为初始强度状态,以冻融循环后为损伤状态,冻融循环后节理试样剪切强度方程为:
cn=c0(1-Dc) (4a)
表2节理试样冻融前后剪切强度参数损伤变量
S3、将岩石试样视为一种天然材料,其矿物组成和胶结作用各不相同,则可假定节理试样由不同材料强度的微单元组成,总微元数为N。Weibull分布函数适用于微元损伤,因此假定节理试样微元强度服从Weibull分布,其损伤概率密度函数可用损伤力F表示为:
在冻融循环作用下,节理试样的损伤是由局部微元体不均匀变形破坏引起的,假设在任意损伤力区间[F,F+dF]内产生的微元数为NP(x)dx,当损伤力达到某一水平F时,则已经损伤破坏的微元总数为:
式中,Nf为破坏的微单元数量,N为微单元的总数。
由于节理试样在冻融循环作用下,节理剪切强度会出现不同程度的损伤劣化,而剪切强度由剪切强度参数表征,剪切强度参数也会呈现损伤劣化,故则粘聚力和内摩擦角的损伤变量Dc和可定义为Nf与N之比,因此,联立式(5)和(6),则Dc和的方程式可分别表示为:
由此得出的方程表示损伤变量Dc和与损伤力F之间的本构关系,但仅适用于不考虑多次循环,冻融损伤力F在整个过程中不变时,损伤本构方程才是静态有效的。对于循环冻融损伤的整个变化动态过程,损伤状态变量Dc和应与多次冻融循环的完整变化过程相结合。
S4、根据多位学者的试验研究发现,冻融循环条件下,节理试样主要在冻胀力作用下出现宏观裂纹并逐渐扩展,同时试样剪切强度呈现明显劣化,故冻融循环过程中冻胀力是导致的节理强度损伤的主要因素。因此,可将冻融循环过程中节理内的损伤力F和冻胀力F(n)视作近似等效。故采用进口薄膜压力传感器Flexiforce A201对整个冻融循环过程中冻胀力进行监测,分别记录下冻融循环第1次、第10次、第20次、第30次、第40次节理内最大冻胀力,试验过程示意图如图3。
传感器记录下的各冻融循环次数下最大冻胀力,试验结果如图4所示,结果显示冻胀力随循环次数增加呈指数函数下降趋势,通过拟合可得出节理内冻胀力F(n)的表达式为:
式中,a为冻胀力随循环次数增加的变化系数,n为冻融循环次数,n≥1。
因此,当节理分布参数a、和一定时,则可计算出一般节理试样在不同冻融循环次数下的损伤粘聚力和内摩擦角。根据表2中的剪切强度参数损伤变量,式(9)中节理试样的参数a、和可分别计算得到,据此可得出(a/Fc)、mc和的值,结果如表3所示。损伤变量随冻融循环次数变化的拟合曲线如图5所示,拟合的相关系数R2均高于0.976,表明推导得出的粘聚力和内摩擦角冻融损伤演化模型是可靠的。
对上述分布参数进行拟合,(a/Fc)、mc和随节理连通率增长均满足指数型函数变化规律,采用Exponential函数对分布参数与节理连通率的关系进行归一化处理,拟合函数曲线见图6,(a/Fc)、和mc、与节理连通率的关系式为:
(a/Fc)=1.897+1.623*exp(-λ/0.181) (10a)
mc=0.210+0.568*exp(-λ/0.266) (10b)
式中,为连通率λ节理冻融循环n次后剪切强度。
同时为验证节理冻融循环损伤节理剪切强度准则,采用与初始试验相同的材料配比和制作方法,制λ=50%节理试样进行冻融试验,分别进行0次、10次、20次、30次及40次冻融循环,并按照前期试验中相同的法向应力水平进行直剪试验,采用Mohr-Coulomb准则计算抗剪强度参数,得出连通率λ=50%的试样经历不同冻融循环次数的黏聚力和内摩擦角。
然后根据式(1)、(2)、(9)、(10)、(11),联立计算得出冻融后节理剪切强度参数的解析解,并与试验值进行对比,结果如图7。由此得出,采用节理冻融损伤下节理剪切强度准则得到的剪切强度参数数值解与实测结果基本上一致。
本发明考虑冻融循环损伤的节理剪切强度准则相对于现有的技术方案,能较为准确的计算出不同节理试样冻融后的剪切强度,准则形式简单便于应用、参数物理意义明确且较易从试验数据中获取。准则预测结果与试验结果具有较高的吻合度,证明所建立的准则合理,能够有效地解决现有冻融环境下节理剪切强度的计算问题。
虽然结合附图对发明的具体实施方式进行了详细地描述,但不应理解为对本专利的保护范围的限定。在权利要求书所描述的范围内,本领域技术人员不经创造性劳动即可做出的各种修改和变形仍属本专利的保护范围。
Claims (8)
1.一种考虑冻融循环损伤的节理剪切强度的确定方法,包括以下步骤:
S1、根据节理试样的连通率变化与剪切强度参数之间的关系,构建计算公式;
S2、根据Lemaitre提出的经典损伤公式,定义黏聚力和内摩擦角损伤变量Dc和并计算得出不同连通率节理试样经历冻融后的黏聚力和内摩擦角损伤变量Dc和然后以冻融循环前强度为初始强度状态,以冻融循环后为损伤状态,结合损伤力学理论和Mohr-Coulomb准则,得出冻融循环后节理试样剪切强度τn方程;
S4、根据实验监测数据,拟合得出节理内冻胀损伤力F(n)与冻融循环次数的关系表达式,然后根据F(n)的表达式,转化S3步骤中黏聚力和内摩擦角损伤变量Dc和的表达式,再然后基于S2步骤中计算得出不同连通率节理试样经历冻融后的黏聚力和内摩擦角损伤变量Dc和确定转化后表达式中模型参数a、和的计算表达式;
2.根据权利要求1所述的考虑冻融循环损伤的节理剪切强度的确定方法,其特征在于,所述步骤S1中,构建计算公式的具体方法,包括以下步骤:分析节理试样未冻融时的剪切强度参数分布特征,可得剪切强度参数随连通率变化情况,进行拟合得到剪切强度参数与连通率的计算公式。
3.根据权利要求1所述的考虑冻融循环损伤的节理剪切强度的确定方法,其特征在于,所述步骤S2的具体方法为:
S2.1分析节理试样经历不同冻融循环后抗剪强度劣化的劣化程度,得出抗剪强度参数的损伤变量随节理冻融循环累积的损伤规律;根据Lemaitre提出的经典损伤公式,黏聚力和内摩擦角损伤变量可以定义为:
S2.3再基于损伤力学理论和Mohr-Coulomb准则,并以冻融循环前强度为初始强度状态,以冻融循环后为损伤状态,冻融循环后节理试样剪切强度方程为:
cn=c0(1-Dc) (4a)
S3.1将岩石试样视为一种天然材料,其矿物组成和胶结作用各不相同,则可假定节理试样由不同材料强度的微单元组成,总微元数为N;Weibull分布函数适用于微元损伤,因此假定节理试样微元强度服从Weibull分布,其损伤概率密度函数可用损伤力F表示为:
S3.2在冻融循环作用下,节理试样的损伤是由局部微元体不均匀变形破坏引起的,假设在任意损伤力区间[F,F+dF]内产生的微元数为NP(x)dx,当损伤力达到某一水平F时,则已经损伤破坏的微元总数为:
式中,Nf为破坏的微单元数量,N为微单元的总数;
S3.3由于节理试样在冻融循环作用下,节理剪切强度会出现不同程度的损伤劣化,而剪切强度由剪切强度参数表征,剪切强度参数也会呈现损伤劣化,故则粘聚力和内摩擦角的损伤变量Dc和可定义为Nf与N之比,因此,联立式(5)和(6),则Dc和的方程式可分别表示为:
5.根据权利要求1所述的考虑冻融循环损伤的节理剪切强度的确定方法,其特征在于,所述S4中步骤:得出冻胀损伤力F(n)与冻融循环次数的关系表达式的具体步骤为:
根据多位学者的试验研究发现,冻融循环条件下,节理试样主要在冻胀力作用下出现宏观裂纹并逐渐扩展,同时试样剪切强度呈现明显劣化,故冻融循环过程中冻胀力是导致的节理强度损伤的主要因素;因此,可将冻融循环过程中节理内的损伤力F和冻胀力F(n)视作近似等效;故采用进口薄膜压力传感器Flexiforce A201对整个冻融循环过程中冻胀力进行监测,分别记录下冻融循环第1次、第10次、第20次、第30次、第40次节理内最大冻胀力;
传感器记录下的各冻融循环次数下最大冻胀力,根据实验数据显示冻胀力随循环次数增加呈指数函数下降趋势,通过拟合可得出节理内冻胀力F(n)的表达式为:
式中:a为冻胀力随循环次数增加的变化系数,n为冻融循环次数,n≥1。
7.根据权利要求4所述的考虑冻融循环损伤的节理剪切强度的确定方法,其特征在于,所述S4步骤中,根据S2步骤中计算得出不同连通率节理试样经历冻融后的粘聚力和内摩擦角损伤变量Dc和确定转换公式中参数的方法为:
对上述分布参数进行拟合,(a/Fc)、mc和随节理连通率增长均满足指数型函数变化规律,采用Exponential函数对分布参数与节理连通率的关系进行归一化处理,得出(a/Fc)、和mc、与节理连通率的关系式为:
(a/Fc)=1.897+1.623*exp(-λ/0.181) (10a)
mc=0.210+0.568*exp(-λ/0.266) (10b)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911167205.8A CN110987662B (zh) | 2019-11-25 | 2019-11-25 | 一种考虑冻融循环损伤的节理剪切强度的确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911167205.8A CN110987662B (zh) | 2019-11-25 | 2019-11-25 | 一种考虑冻融循环损伤的节理剪切强度的确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110987662A true CN110987662A (zh) | 2020-04-10 |
CN110987662B CN110987662B (zh) | 2022-04-29 |
Family
ID=70086516
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911167205.8A Active CN110987662B (zh) | 2019-11-25 | 2019-11-25 | 一种考虑冻融循环损伤的节理剪切强度的确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110987662B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111398019A (zh) * | 2020-04-22 | 2020-07-10 | 中国科学院武汉岩土力学研究所 | 一种快速判断不同应变率荷载下岩石损伤相对大小的方法 |
CN113063675A (zh) * | 2021-03-26 | 2021-07-02 | 长安大学 | 一种软-硬节理峰值剪切强度的计算方法 |
CN113063810A (zh) * | 2021-03-22 | 2021-07-02 | 西安科技大学 | 获得砂岩冻融作用下宏细观损伤演化规律的方法 |
CN113155699A (zh) * | 2021-04-09 | 2021-07-23 | 大连海事大学 | 一种热-水-力共同作用的岩石统计损伤计算方法及其应用 |
Citations (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU920451A1 (ru) * | 1979-04-23 | 1982-04-15 | Московский Ордена Трудового Красного Знамени Инженерно-Строительный Институт Им.Куйбышева | Способ определени сопротивлени грунта сдвигу |
CN101739716A (zh) * | 2009-12-04 | 2010-06-16 | 东北大学 | 工程岩体三维空间结构建模与关键块识别方法 |
CN101936008A (zh) * | 2010-09-30 | 2011-01-05 | 东北大学 | 岩体边坡三维模型及块体滑落分析方法 |
CN101986130A (zh) * | 2010-09-06 | 2011-03-16 | 水利部交通运输部国家能源局南京水利科学研究院 | 薄弱层面剪切强度特性的三轴测试方法 |
JP2013144921A (ja) * | 2011-12-14 | 2013-07-25 | National Agriculture & Food Research Organization | 地盤の現位置剪断強度測定方法および地盤の現位置剪断強度測定装置 |
US20130275099A1 (en) * | 2012-04-17 | 2013-10-17 | Schlumberger Technology Corporation | Determining A Limit Of Failure In A Wellbore Wall |
CN104361211A (zh) * | 2014-10-24 | 2015-02-18 | 中冶长天国际工程有限责任公司 | 岩石统计损伤本构模型的构建和应用方法 |
CN104977321A (zh) * | 2015-07-14 | 2015-10-14 | 中国矿业大学 | 冻融损伤梯度与各向异性测量方法 |
CN107515291A (zh) * | 2017-08-23 | 2017-12-26 | 西安科技大学 | 一种等围压作用下岩石冻融损伤本构模型的构建方法 |
CN108829916A (zh) * | 2018-04-25 | 2018-11-16 | 中铁二院工程集团有限责任公司 | 硬岩冻融损伤长期变形模型的构建方法 |
CN108982328A (zh) * | 2018-08-14 | 2018-12-11 | 中南大学 | 一种计算冻融作用下岩石孔隙体积变形的方法 |
CN109283086A (zh) * | 2018-11-26 | 2019-01-29 | 中国地质大学(北京) | 精准测算反复冻融损伤后寒区隧道岩石弹性模量的方法 |
CN109283215A (zh) * | 2018-11-26 | 2019-01-29 | 中国地质大学(北京) | 反复冻融损伤后寒区隧道围岩冻胀力的测算方法 |
CN109885980A (zh) * | 2019-03-29 | 2019-06-14 | 中南大学 | 基于应力差确定屈服点的节理剪切全过程损伤本构模型 |
CN110441174A (zh) * | 2019-07-09 | 2019-11-12 | 郑州大学 | 一种研究循环动荷载下应变硬化土疲劳损伤判定的方法 |
-
2019
- 2019-11-25 CN CN201911167205.8A patent/CN110987662B/zh active Active
Patent Citations (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU920451A1 (ru) * | 1979-04-23 | 1982-04-15 | Московский Ордена Трудового Красного Знамени Инженерно-Строительный Институт Им.Куйбышева | Способ определени сопротивлени грунта сдвигу |
CN101739716A (zh) * | 2009-12-04 | 2010-06-16 | 东北大学 | 工程岩体三维空间结构建模与关键块识别方法 |
CN101986130A (zh) * | 2010-09-06 | 2011-03-16 | 水利部交通运输部国家能源局南京水利科学研究院 | 薄弱层面剪切强度特性的三轴测试方法 |
CN101936008A (zh) * | 2010-09-30 | 2011-01-05 | 东北大学 | 岩体边坡三维模型及块体滑落分析方法 |
JP2013144921A (ja) * | 2011-12-14 | 2013-07-25 | National Agriculture & Food Research Organization | 地盤の現位置剪断強度測定方法および地盤の現位置剪断強度測定装置 |
US20130275099A1 (en) * | 2012-04-17 | 2013-10-17 | Schlumberger Technology Corporation | Determining A Limit Of Failure In A Wellbore Wall |
CN104361211A (zh) * | 2014-10-24 | 2015-02-18 | 中冶长天国际工程有限责任公司 | 岩石统计损伤本构模型的构建和应用方法 |
CN104977321A (zh) * | 2015-07-14 | 2015-10-14 | 中国矿业大学 | 冻融损伤梯度与各向异性测量方法 |
CN107515291A (zh) * | 2017-08-23 | 2017-12-26 | 西安科技大学 | 一种等围压作用下岩石冻融损伤本构模型的构建方法 |
CN108829916A (zh) * | 2018-04-25 | 2018-11-16 | 中铁二院工程集团有限责任公司 | 硬岩冻融损伤长期变形模型的构建方法 |
CN108982328A (zh) * | 2018-08-14 | 2018-12-11 | 中南大学 | 一种计算冻融作用下岩石孔隙体积变形的方法 |
CN109283086A (zh) * | 2018-11-26 | 2019-01-29 | 中国地质大学(北京) | 精准测算反复冻融损伤后寒区隧道岩石弹性模量的方法 |
CN109283215A (zh) * | 2018-11-26 | 2019-01-29 | 中国地质大学(北京) | 反复冻融损伤后寒区隧道围岩冻胀力的测算方法 |
CN109885980A (zh) * | 2019-03-29 | 2019-06-14 | 中南大学 | 基于应力差确定屈服点的节理剪切全过程损伤本构模型 |
CN110441174A (zh) * | 2019-07-09 | 2019-11-12 | 郑州大学 | 一种研究循环动荷载下应变硬化土疲劳损伤判定的方法 |
Non-Patent Citations (10)
Title |
---|
DAXING LEI等: "《Effect of Cyclic Freezing-Thawing on the Shear Mechanical Characteristics of Nonpersistent Joints》", 《ADVANCES IN MATERIALS SCIENCE AND ENGINEERING》 * |
JIAN-QIAO MU等: "《Degradation characteristics of shear strength of joints in three rock types due to cyclic freezing and thawing》", 《COLD REGIONS SCIENCE AND TECHNOLOGY》 * |
PHILIPP MAMOT等: "《A temperature-and stress-controlled failure criterionfor ice-filled permafrost rock joints》", 《THE CRYOSPHERE》 * |
YANLONG CHEN等: "《Effects of Freezing and Thawing Cycle on Mechanical Properties and Stability of Soft Rock Slope》", 《ADVANCES IN MATERIALS SCIENCE AND ENGINEERING》 * |
刘泉声等: "《岩体冻融疲劳损伤模型与评价指标研究》", 《岩石力学与工程学报》 * |
刘红岩等: "《循环冻融条件下节理岩体损伤破坏试验研究》", 《岩土力学》 * |
李泽栋等: "《节理试件峰后剪切变形特性研究》", 《湖北工业大学学报》 * |
杨更社等: "《冻融环境下岩体损伤力学特性多尺度研究及进展》", 《岩石力学与工程学报》 * |
林杭等: "《强度折减法在Hoek-Brown准则中的应用》", 《中南大学学报(自然科学版)》 * |
王乐华等: "《节理岩体冻融力学特性试验研究》", 《水利水电技术》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111398019A (zh) * | 2020-04-22 | 2020-07-10 | 中国科学院武汉岩土力学研究所 | 一种快速判断不同应变率荷载下岩石损伤相对大小的方法 |
CN113063810A (zh) * | 2021-03-22 | 2021-07-02 | 西安科技大学 | 获得砂岩冻融作用下宏细观损伤演化规律的方法 |
CN113063675A (zh) * | 2021-03-26 | 2021-07-02 | 长安大学 | 一种软-硬节理峰值剪切强度的计算方法 |
CN113063675B (zh) * | 2021-03-26 | 2023-06-23 | 长安大学 | 一种软-硬节理峰值剪切强度的计算方法 |
CN113155699A (zh) * | 2021-04-09 | 2021-07-23 | 大连海事大学 | 一种热-水-力共同作用的岩石统计损伤计算方法及其应用 |
CN113155699B (zh) * | 2021-04-09 | 2024-04-16 | 大连海事大学 | 一种热-水-力共同作用的岩石统计损伤计算方法及应用 |
Also Published As
Publication number | Publication date |
---|---|
CN110987662B (zh) | 2022-04-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110987662B (zh) | 一种考虑冻融循环损伤的节理剪切强度的确定方法 | |
Liu et al. | A prediction model for uniaxial compressive strength of deteriorated rocks due to freeze–thaw | |
Hughes et al. | Numerical analysis of high strain rate splitting-tensile tests | |
CN111476404B (zh) | 一种冻融损伤软岩长期变形的预测方法 | |
Wu et al. | Damage degree evaluation of masonry using optimized SVM-based acoustic emission monitoring and rate process theory | |
Hariri-Ardebili et al. | Estimation of probable damages in arch dams subjected to strong ground motions using endurance time acceleration functions | |
CN105046051A (zh) | 一种实际环境冻融作用次数预估方法 | |
CN105891029A (zh) | 一种基于断裂力学理论的沥青疲劳性能评价方法 | |
Yahaghi et al. | Experimental, theoretical and numerical modelling of the deterioration and failure process of sandstones subject to freeze–thaw cycles | |
Birkeland et al. | The stuffblock snow stability test: comparability with the rutschblock, usefulness in different snow climates, and repeatability between observers | |
CN111159794A (zh) | 多裂隙类岩石试样力学性质的几何损伤流变分析方法 | |
CN109541192A (zh) | 一种道路沥青自愈行为特性的评价方法 | |
Secor et al. | Viscoelastic response of asphalt paving slabs under creep loading | |
CN115235884A (zh) | 一种饱和砂岩的损伤状况判别及剩余强度预估方法 | |
CN111400926A (zh) | 一种基于flac3d数值软件的巷道围岩热物理参数反演方法 | |
CN110987661A (zh) | 一种改进Harris分布的结构面剪切损伤本构模型的方法 | |
Park et al. | Temperature correction of backcalculated moduli and deflections using linear viscoelasticity and time-temperature superposition | |
US8880383B1 (en) | Determining conductivity of aged fracture proppants for simulating flow in a fractured reservoir | |
Kulhawy et al. | Uncertainty in basic properties of geomaterials | |
Sherafati et al. | An investigation into the time dependency of shear strength of clay brick walls; an approximate approach | |
CN117805351B (zh) | 一种沥青混合料抗疲劳性能综合评价方法 | |
Pavements | Prediction of low-temperature creep and thermal strain in asphalt concrete pavements | |
CN116935983B (zh) | 一种冻融后岩石物理力学性质衰减程度的预测方法 | |
Casini et al. | Coupled phenomena induced by freezing in a granular material | |
Poulos et al. | Method of calculating long-term creep settlements |
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 | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20240112 Address after: 230000 floor 1, building 2, phase I, e-commerce Park, Jinggang Road, Shushan Economic Development Zone, Hefei City, Anhui Province Patentee after: Dragon totem Technology (Hefei) Co.,Ltd. Address before: Yuelu District City, Hunan province 410083 Changsha Lushan Road No. 932 Patentee before: CENTRAL SOUTH University |