粗粒土burgers模型参数确定方法及考虑湿化作用模型构建方法技术

本发明专利技术为粗粒土burgers模型参数确定方法及考虑湿化作用模型构建方法，首先通过粗粒土直剪蠕变试验验获得不同湿化程度下土体蠕变变形数据。然后试验数据可简单确定burgers模型所需本构参数，建立不同湿化程度下土体参数公式，使参数可直接用于数值模拟。对实际工程建立数值模型，利用上述构建方法获得本构参数已知的考虑湿化作用的burgers模型，能获得合理的本构和外部条件，用来预测长期变形。试验装置中将直剪盒设置为外部直剪盒、内部直剪盒、井字支架和储水袋组装形式，代替由上部、下部两个等大的剪切盒构成的传统直剪盒，能够实现浸泡、干湿循环两种工况的长期剪切试验，即实现了长期浸泡、干湿循环工况的长期加载。干湿循环工况的长期加载。干湿循环工况的长期加载。

【技术实现步骤摘要】
粗粒土burgers模型参数确定方法及考虑湿化作用模型构建方法


[0001]本专利技术属于岩土工程
，具体来说涉及一种适用于粗粒土体的burgers模型参数确定方法及考虑湿化作用模型构建方法。

技术介绍

[0002]由粗粒土构成的土石坝是我国目前大坝的主要坝型，然而在土石坝的实际运行中其坝体会发生比较明显的变形，且这种变形通常持续时间较久。引起坝体的长期变形的重要原因是其主要建筑材料
‑
粗粒土体。目前对于粗粒土体的长期变形预测是通过构建本构模型来进行数值模计算，但大多数模型基本只考虑了荷载因素对土体长期变形预测的影响，并未考虑湿化作用。根据现有观测资料表明：坝体在降雨前后的变形的速率明显超过非降雨时间段。如何考虑湿化作用对粗粒土体的长期变形影响是目前研究领域的空白。
[0003]FLAC3D程序为模拟材料的长期变形提供了一种非线性变形模型——burgers模型，但是其所需要的本构参数一般并不能由试验直接确定，本构参数的确定往往存在很大的人为因素与不确定性。为了将湿化作用考虑其中，提供了一种新的burgers模型的构建方法，实现模型的客观便捷使用。

技术实现思路

[0004]针对现有技术的不足，本专利技术的目的在于提供一种适用于粗粒土体的burgers模型参数确定方法及考虑湿化作用模型构建方法，该构建方法所获得的考虑湿化作用的burgers模型适用于不同湿化程度的粗粒土体长期变形预测。
[0005]本专利技术的目的是通过下述技术方案予以实现的：
[0006]第一方面，本专利技术提供一种粗粒土burgers模型参数确定方法参数确定方法，所述确定方法包括以下内容：
[0007]将burgers模型的应力应变方程表示为式(2)
[0008]ε＝A+Bt+C(1
‑
e
‑
Dt
)
ꢀꢀꢀ
(2)，
[0009]其中，t为时间；ε为土体的剪切应变；A、B、C、D分别为参数，其中，t为时间；ε为土体的剪切应变；A、B、C、D分别为参数，τ为材料所受剪应力；m
e
为burgers模型中Maxwell体弹性模量、m
η
为burgers模型中Maxwell体黏性模量、k
e
为burgers模型中kelvin体弹性模量、k
η
为burgers模型中kelvin体黏性模量；
[0010]针对不同参数特性采取直剪蠕变试验和参数拟合的方式确定A、B、C、D参数值；
[0011]参数A的确定过程是：设定目标剪应力值，从开始直剪蠕变试验后至剪应力第一次达到设定目标剪应力值的时间段内粗粒土体的剪切应变量值为参数A，再根据及设定目标剪应力值确定本构参数m
e
；
[0012]参数B、C的确定过程是：
[0013]进行直剪蠕变试验，获得待测对象的剪切应变
‑
时间曲线，当t
→
∞时,式(2)为：ε＝A+Bt+C，此时剪切应变
‑
时间曲线近似为斜直线，做其曲线的渐近线，获得渐近线的斜率和截距的值，渐近线与应变轴交于点(0，A+C)；在确定参数A的值后，根据点(0，A+C)求得参数C，再根据确定本构参数k
e
；所述渐近线的表达式可利用土体剪切应变
‑
时间曲线的末端平缓变形段上任取两点算出；
[0014]根据渐近线斜率的值确定参数B，再根据确定本构参数m
η
，此时τ为t
→
∞时材料所受剪应力；
[0015]参数D拟合：
[0016]再确定参数A、B、C的值后，利用oringin软件，输入公式(2)，并将A、B、C值输入，让oringin软件根据直剪蠕变试验确定的土体剪切应变
‑
时间曲线自动计算出参数D,并根据此求得本构参数k
η
。
[0017]第二方面，本专利技术提供一种考虑湿化作用模型构建方法，所述模型为粗粒土的burgers模型，为一种流变模型，所述构建方法包括以下内容：
[0018]1)在自然含水率状态下开展不同法向应力σ下不同剪应力τ的土体直剪蠕变试验，得到自然含水率状态下不同应力下的土体剪切应变
‑
时间曲线；根据自然含水率状态下不同应力下的土体剪切应变
‑
时间曲线分别利用上述的burgers模型参数确定方法确定对应的A、B、C、D参数值和本构参数值，进而获得自然含水率状态下四个本构参数与应力之间的函数关系，即为自然含水率状态下本构参数表达式；
[0019]2)分别在浸泡、干湿循环状态下开展与自然含水率状态下相同法向应力和剪应力下的土体直剪蠕变试验，得到浸泡、干湿循环状态下不同应力下的土体剪切应变
‑
时间曲线；根据浸泡、干湿循环状态下不同应力下的土体应变
‑
时间曲线分别确定对应的A、B、C、D参数值和本构参数值，进而获得浸泡、干湿循环状态下四个本构参数与应力之间的函数关系；
[0020]3)分别将浸泡、干湿循环状态下的本构参数与自然含水率状态下的相应本构参数进行对比，分别确定浸泡、干湿循环状态下的四个本构参数的浸泡湿化系数和干湿循环湿化系数；
[0021]4)以自然含水率状态下本构参数表达式为基准，对自然含水率状态下的本构参数表达式中引入上述的浸泡湿化系数和干湿循环湿化系数来转化成浸泡、干湿循环下本构参数表达式；
[0022]5)在确定应力状态、湿化情况后，带入步骤4)相应的表达式中获得考虑湿化作用的burgers模型的本构参数，进而建立本构参数已知的考虑湿化作用的burgers模型。
[0023]获得不同种类材料粗粒的浸泡湿化系数和干湿循环湿化系数，建立材料粗粒种类与湿化系数的对应关系，在已知材料种类及湿化情况时，能直接确定相应的湿化系数，进而确定对应的考虑湿化作用的burgers模型。
[0024]将考虑湿化作用的burgers模型直接用于数值模拟，用于不同湿化情况下的长期变形预测。
[0025]自然含水率状态下本构参数表达式为：
[0026]m
e
＝120*σ
‑
2500；
[0027]m
η
＝τ/(1*10
‑5)；
[0028]k
e
＝60*τ+100；
[0029]k
η
＝290*τ+500；
[0030]其中，σ为法向应力，τ为剪应力；
[0031]所述浸泡湿化系数为：修正kelvin体的k
e
的浸泡湿化系数a
ke
、修正kelvin体的k
η
的浸泡湿化系数a
kη
、修正Maxwell体的m
e
的浸泡湿化系数a
me
，
[0032]所述干湿循环湿化本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种粗粒土burgers模型参数确定方法，其特征在于，所述确定方法包括以下内容：将burgers模型的应力应变方程表示为式(2)ε＝A+Bt+C(1
‑
e
‑
Dt
)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(2)，其中，t为时间；ε为土体的剪切应变；A、B、C、D分别为参数，其中，t为时间；ε为土体的剪切应变；A、B、C、D分别为参数，τ为材料所受剪应力；m
e
为burgers模型中Maxwell体弹性模量、m
η
为burgers模型中Maxwell体黏性模量、k
e
为burgers模型中kelvin体弹性模量、k
η
为burgers模型中kelvin体黏性模量；针对不同参数特性采取直剪蠕变试验和参数拟合的方式确定A、B、C、D参数值；参数A的确定过程是：设定目标剪应力值，从开始直剪蠕变试验后至剪应力第一次达到设定目标剪应力值的时间段内粗粒土体的剪切应变量值为参数A，再根据及设定目标剪应力值确定本构参数m
e
；参数B、C的确定过程是：进行直剪蠕变试验，获得待测对象的剪切应变
‑
时间曲线，当t
→
∞时,式(2)为：ε＝A+Bt+C，此时剪切应变
‑
时间曲线近似为斜直线，做其曲线的渐近线，获得渐近线的斜率和截距的值，渐近线与应变轴交于点(0，A+C)；在确定参数A的值后，根据点(0，A+C)求得参数C，再根据确定本构参数k
e
；根据渐近线斜率的值确定参数B，再根据确定本构参数m
η
，此时τ为t
→
∞时材料所受剪应力；参数D拟合：再确定参数A、B、C的值后，利用oringin软件，输入公式(2)，并将A、B、C值输入，让oringin软件根据直剪蠕变试验确定的土体剪切应变
‑
时间曲线自动计算出参数D,并根据此求得本构参数k
η
。2.一种考虑湿化作用模型构建方法，其特征在于，所述模型为粗粒土的burgers模型，所述构建方法包括以下内容：1)在自然含水率状态下开展不同法向应力σ下不同剪应力τ的土体直剪蠕变试验，得到自然含水率状态下不同应力下的土体剪切应变
‑
时间曲线；根据自然含水率状态下不同应力下的土体剪切应变
‑
时间曲线分别利用权利要求1所述的burgers模型参数确定方法确定对应的A、B、C、D参数值和本构参数值，进而获得自然含水率状态下四个本构参数与应力之间的函数关系，即为自然含水率状态下本构参数表达式；2)分别在浸泡、干湿循环状态下开展与自然含水率状态下相同法向应力和剪应力下的土体直剪蠕变试验，得到浸泡、干湿循环状态下不同应力下的土体剪切应变
‑
时间曲线；根据浸泡、干湿循环状态下不同应力下的土体应变
‑
时间曲线分别确定对应的A、B、C、D参数值和本构参数值，进而获得浸泡、干湿循环状态下四个本构参数与应力之间的函数关系；3)分别将浸泡、干湿循环状态下的本构参数与自然含水率状态下的相应本构参数进行
对比，分别确定浸泡、干湿循环状态下的四个本构参数的浸泡湿化系数和干湿循环湿化系数；4)以自然含水率状态下本构参数表达式为基准，对自然含水率状态下的本构参数表达式中引入上述的浸泡湿化系数和干湿循环湿化系数来转化成浸泡、干湿循环下本构参数表达式；5)在确定应力状态、湿化情况后，带入步骤4)相应的表达式中获得考虑湿化作用的burgers模型的本构参数，...

【专利技术属性】
技术研发人员：王子寒，刘建泽，景晓昆，肖成志，黄达，
申请(专利权)人：河北工业大学，
类型：发明
国别省市：