CN113155612A - 微纤维混合硅溶胶固化钙质砂的变形预报方法 - Google Patents
微纤维混合硅溶胶固化钙质砂的变形预报方法 Download PDFInfo
- Publication number
- CN113155612A CN113155612A CN202110413264.XA CN202110413264A CN113155612A CN 113155612 A CN113155612 A CN 113155612A CN 202110413264 A CN202110413264 A CN 202110413264A CN 113155612 A CN113155612 A CN 113155612A
- Authority
- CN
- China
- Prior art keywords
- epsilon
- strain
- stress
- time
- silica sol
- 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
- 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/08—Investigating strength properties of solid materials by application of mechanical stress by applying steady tensile or compressive 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/24—Investigating strength properties of solid materials by application of mechanical stress by applying steady shearing forces
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)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明涉及微纤维混合硅溶胶固化钙质砂的变形预报方法,对微纤维混合硅溶胶渗流固化钙质砂后形成的硅凝胶‑微纤维‑砂复合体,更为准确的预报其受力过程中的体积变形和体积变形率。
Description
技术领域
本发明属于岩土工程研究领域,尤其涉及一种微纤维混合硅溶胶固化钙质砂的变形预报方法。
背景技术
硅溶胶为纳米级的二氧化硅颗粒悬浮在水中,且硅溶胶具有近似水的粘度,硅溶胶快速渗流经过砂土后,二氧化硅颗粒聚集成三维网状结构的硅凝胶,从而胶结固化土体。但是硅凝胶中存在大量微裂纹,因此在硅溶胶中混合微纤维,然后将此硅溶胶渗流经过土体,形成硅凝胶-微纤维-砂复合体,微纤维可以抑制硅凝胶微裂纹扩展,从而提高硅溶胶渗流固化砂土的强度。
但是微纤维混合硅溶胶固化钙质砂后,固化土体受力变形时,同时受钙质砂破碎影响、硅凝胶胶结失效影响和微纤维增强影响,而用常规方法预报硅凝胶-微纤维-砂复合体的变形存在较大误差,因此缺乏一种方法,对微纤维混合硅溶胶渗流固化钙质砂后形成的硅凝胶-微纤维-砂复合体,可以准确预报其受力过程中的体积变形和体积变形率。
发明内容
本发明为了对微纤维混合硅溶胶渗流固化钙质砂后形成的硅凝胶-微纤维-砂复合体,更为准确的预报其受力过程中的体积变形和体积变形率,本发明提供了一种微纤维混合硅溶胶固化钙质砂的变形预报方法。
本发明涉及一些简写和符号,以下为注解:
σ1′:颗粒集合体受到的竖向应力;
σ′2和σ3′:颗粒集合体受到的水平应力,σ′2和σ3′的方向垂直;
ε1、ε2和ε3:应变,且分别与应力σ1、σ2和σ3方向相同;
ηcritical:临界应力比η;
ηpeak:为峰值应力比η;
εv:体应变,εv=ε1+ε2+ε3;
t0,t1,t2,…,ti,…,tn:加载过程中记录的起始时刻为t0,之后记录的时刻从小到大依次为t1,t2,…,ti,…,tn,这里1≤i≤n,n+1为记录的时刻点数;
(εv)i:第ti时刻对应的体应变εv;
(εs)i:第ti时刻对应的剪应变εs;
(△εv)i:体应变增量,(△εv)i=(εv)i-(εv)i-1;
(△εs)i:剪应变增量,(△εs)i=(εs)i-(εs)i-1;
ηi:ti时刻(1≤i≤n)应力比ηi;
f(εs):f(εs)是剪应变εs的函数;
α:材料系数;
K:材料系数;
e0:初始孔隙比;
ec:临界孔隙比;
e:ti时刻(1≤i≤n)孔隙比;
κ:回弹系数;
M:材料参数且等于临界应力比ηcritical;
m:材料参数。
本发明的技术方案:一种微纤维混合硅溶胶固化钙质砂的变形预报方法,具体包含如下步骤:
步骤1:对于微纤维混合硅溶胶渗流固化钙质砂后形成的硅凝胶-微纤维-砂复合体,设设其竖向受应力为σ1′,水平面上受应力分别为σ′2和σ3′,其中σ′2和σ3′的方向垂直,颗粒集合体的应变为ε1、ε2和ε3,其中应变ε1、ε2和ε3的方向分别与应力σ1′、σ′2和σ3′方向相同;定义平均有效应力p′、剪应力q、应力比η、体应变εv和剪应变εs:
εv=ε1+ε2+ε3 (4)
步骤2:设加载过程中记录的起始时刻为t0,之后记录的时刻从小到大依次为t1,t2,…,ti,…,tn,这里1≤i≤n,n+1为记录的时刻点数,设第ti时刻对应的体应变εv为(εv)i,设第ti时刻对应的剪应变εs为(εs)i,设相邻两个时刻差产生的剪应变增量(△εs)i相等,定义体应变增量(△εv)i和剪应变增量(△εs)i:
(△εv)i=(εv)i-(εv)i-1 (6)
(△εs)i=(εs)i-(εs)i-1 (7)
这里每个时刻剪应变(εs)i和剪应变增量(△εs)i为已知设定值;
步骤3:计算ti时刻(1≤i≤n)体应变增量(△εv)i:
(△εv)i=f(εs)·(△εs)i (8)
步骤4:计算ti时刻(1≤i≤n)体应变(εv)i和孔隙比e:
(εv)i=(εv)i-1+(△εv)i (9)
e=e0-(1+e0)(εv)i (10)
步骤5:计算ti时刻(1≤i≤n)应力比ηi:
应力比η的计算公式为:
上式中,α和K为材料系数,ηpeak为峰值应力比η;式(11)中,η取ti时刻(1≤i≤n)应力比ηi,εs取ti时刻剪应变(εs)i;
式(12)中,取其中ec为临界孔隙比,e为ti时刻(1≤i≤n)孔隙比,M为材料参数且等于临界应力比ηcritical,m为材料参数,这里的临界状态指应力比η随剪应变变化曲线中过峰值后的稳定状态、且体应变随剪应变变化曲线中过峰值后的稳定状态,如果体应变随剪应变变化曲线中过峰值后一直没有稳定,则取加载最后时刻记录临界应力比ηcritical和临界孔隙比ec;取后,式(12)写为:
步骤6:重复步骤3-步骤5,直至计算得到每个ti时刻(1≤i≤n)对应的体应变(εv)i和应力比ηi。
上述计算步骤1-6中,计算得到应力比η,当水平向有效应力σ′2和σ3′恒定且相等时,即可得到竖向有效应力σ1′、以及平均有效应力和剪应力:
上述计算步骤1-6中,当水平向应变ε2和ε3相等时,即可得到竖向应变ε1和水平向应变ε3:
附图说明
图1为微纤维混合硅溶胶渗流固化钙质砂后形成的硅凝胶-微纤维-砂复合体,其受竖向和水平向应力示意图;
图2为应力比随剪应变变化曲线,体应变随剪应变变化曲线;
图3为体应变增量与剪应变增量比值随应力比变化曲线;
1.微纤维混合硅溶胶渗流固化钙质砂后形成的硅凝胶-微纤维-砂复合体。
具体实施方式
为了使本发明实现的技术手段、创新特征、达成目的与功效易于明白了解,下面结合具体图示,进一步阐述本发明。
本发明涉及一些简写和符号,以下为注解:
σ1′:颗粒集合体受到的竖向应力
σ′2和σ3′:颗粒集合体受到的水平应力,σ′2和σ3′的方向垂直ε1、ε2和ε3:应变,且分别与应力σ1、σ2和σ3方向相同
ηcritical:临界应力比η
ηpeak:为峰值应力比η
εv:体应变,εv=ε1+ε2+ε3
(εv)i:第ti时刻对应的体应变εv
(εs)i:第ti时刻对应的剪应变εs
(△εv)i:体应变增量,(△εv)i=(εv)i-(εv)i-1
(△εs)i:剪应变增量,(△εs)i=(εs)i-(εs)i-1
ηi:ti时刻(1≤i≤n)应力比ηi
f(εs):f(εs)是剪应变εs的函数
α:材料系数
K:材料系数
e0:初始孔隙比
ec:临界孔隙比
e:ti时刻(1≤i≤n)孔隙比
κ:回弹系数
M:材料参数且等于临界应力比ηcritical
m:材料参数
本发明的技术方案:一种微纤维混合硅溶胶固化钙质砂的变形预报方法,具体包含如下步骤:
步骤1:如图1所示,硅溶胶和碳化硅纳米线混合后渗流经过钙质砂,形成硅凝胶-微纤维-砂复合体1,设设其竖向受应力为σ1′,水平面上受应力分别为σ′2和σ3′,其中σ′2和σ3′的方向垂直,颗粒集合体的应变为ε1、ε2和ε3,其中应变ε1、ε2和ε3的方向分别与应力σ1′、σ′2和σ3′方向相同;定义平均有效应力p′、剪应力q、应力比η、体应变εv和剪应变εs:
εv=ε1+ε2+ε3 (4)
在这一系列试验中,竖向加载方式为应变加载,即水平向位移为常数,试验为排水固结三轴试验,水平向有效应力σ′2和σ3′恒定且相等,且σ3′为围压,围压分别设置为100kPa,200kPa和600kPa,水平向应变ε2和ε3相等,有如下3个关系式:
步骤2:设加载过程中记录的起始时刻为t0,之后记录的时刻从小到大依次为t1,t2,…,ti,…,tn,这里1≤i≤n,n+1为记录的时刻点数,设第ti时刻对应的体应变εv为(εv)i,设第ti时刻对应的剪应变εs为(εs)i,设相邻两个时刻差产生的剪应变增量(△εs)i相等,定义体应变增量(△εv)i和剪应变增量(△εs)i:
(△εv)i=(εv)i-(εv)i-1 (6)
(△εs)i=(εs)i-(εs)i-1 (7)
这里每个时刻剪应变(εs)i和剪应变增量(△εs)i为已知设定值;
步骤3:计算ti时刻(1≤i≤n)体应变增量(△εv)i:
(△εv)i=f(εs)·(△εs)i (8)
步骤4:计算ti时刻(1≤i≤n)体应变(εv)i和孔隙比e:
(εv)i=(εv)i-1+(△εv)i (9)
e=e0-(1+e0)(εv)i (10)
步骤5:计算ti时刻(1≤i≤n)应力比ηi:
应力比η的计算公式为:
上式中,α和K为材料系数,ηpeak为峰值应力比η;式(11)中,η取ti时刻(1≤i≤n)应力比ηi,εs取ti时刻剪应变(εs)i;
式(12)中,取其中ec为临界孔隙比,e为ti时刻(1≤i≤n)孔隙比,M为材料参数且等于临界应力比ηcritical,m为材料参数,这里的临界状态指应力比η随剪应变变化曲线中过峰值后的稳定状态、且体应变随剪应变变化曲线中过峰值后的稳定状态,如果体应变随剪应变变化曲线中过峰值后一直没有稳定,则取加载最后时刻记录临界应力比ηcritical和临界孔隙比ec;取后,式(12)写为:
步骤6:重复步骤3-步骤5,直至计算得到每个ti时刻(1≤i≤n)对应的体应变(εv)i和应力比ηi。
作为对照的试验样品为选取硅溶胶浓度为20%,微纤维为碳化硅纳米线,其浓度为0.01%,硅溶胶和碳化硅纳米线混合后渗流经过钙质砂,形成硅凝胶-微纤维-砂复合体1,在这一系列试验中,竖向加载方式为应变加载,即水平向位移为常数,试验为排水固结三轴试验,水平向有效应力σ′2和σ3′恒定且相等,且σ3′为围压,围压分别设置为100kPa,200kPa和600kPa。
如图2所示,其上部a部分显示用这里给出的方法预报出围压为100kPa,200kPa和600kPa时应力比随剪应变变化曲线、体应变随剪应变变化曲线,可以看出预报曲线和实际试验曲线吻合良好;而其下部b部分显示用传统的体应变增量与剪应变增量比值公式(△εv)i=(M-η)·(△εs)i得到的曲线,其他部分的计算与这里提出的方法相同,可以看到,b部分显示的预报体应变曲线误差更大,特别是剪应变从0-5%时其误差非常明显。
Claims (5)
1.一种微纤维混合硅溶胶固化钙质砂的变形预报方法,其特征为:
其包括如下步骤:
步骤1:对于微纤维混合硅溶胶渗流固化钙质砂后形成的硅凝胶-微纤维-砂复合体,设设其竖向受应力为σ′1,水平面上受应力分别为σ′2和σ′3,其中σ′2和σ′3的方向垂直,颗粒集合体的应变为ε1、ε2和ε3,其中应变ε1、ε2和ε3的方向分别与应力σ′1、σ′2和σ′3方向相同;定义平均有效应力p′、剪应力q、应力比η、体应变εv和剪应变εs:
εv=ε1+ε2+ε3 (4)
步骤2:设加载过程中记录的起始时刻为t0,之后记录的时刻从小到大依次为t1,t2,…,ti,…,tn,这里1≤i≤n,n+1为记录的时刻点数,设第ti时刻对应的体应变εv为(εv)i,设第ti时刻对应的剪应变εs为(εs)i,设相邻两个时刻差产生的剪应变增量(△εs)i相等,定义体应变增量(△εv)i和剪应变增量(△εs)i:
(△εv)i=(εv)i-(εv)i-1 (6)
(△εs)i=(εs)i-(εs)i-1 (7)
这里每个时刻剪应变(εs)i和剪应变增量(△εs)i为已知设定值;
步骤3:计算ti时刻体应变增量(△εv)i:
(△εv)i=f(εs)·(△εs)i (8)
步骤4:计算ti时刻体应变(εv)i和孔隙比e:
(εv)i=(εv)i-1+(△εv)i (9)
e=e0-(1+e0)(εv)i (10)
步骤5:计算ti时刻应力比ηi:
应力比η的计算公式为:
上式中,α和K为材料系数,ηpeak为峰值应力比η;式(11)中,η取ti时刻应力比ηi,εs取ti时刻剪应变(εs)i;
步骤6:重复步骤3-步骤5,直至计算得到每个ti时刻对应的体应变(εv)i和应力比ηi。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110413264.XA CN113155612B (zh) | 2021-04-16 | 2021-04-16 | 微纤维混合硅溶胶固化钙质砂的变形预报方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110413264.XA CN113155612B (zh) | 2021-04-16 | 2021-04-16 | 微纤维混合硅溶胶固化钙质砂的变形预报方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113155612A true CN113155612A (zh) | 2021-07-23 |
CN113155612B CN113155612B (zh) | 2022-09-30 |
Family
ID=76868424
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110413264.XA Active CN113155612B (zh) | 2021-04-16 | 2021-04-16 | 微纤维混合硅溶胶固化钙质砂的变形预报方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113155612B (zh) |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100071309A1 (en) * | 2008-05-01 | 2010-03-25 | Brangan Charles D | Soil stabilization system, stabilized soil comprising same, and a method of stabilizing soil |
CN101779204A (zh) * | 2007-08-17 | 2010-07-14 | 宝洁公司 | 广义本构建模方法和*** |
CN103528897A (zh) * | 2013-10-17 | 2014-01-22 | 三峡大学 | 高围压下粗粒土大三轴试验数据临界状态参数确定方法 |
US20150057988A1 (en) * | 2013-08-26 | 2015-02-26 | The University Of Akron | Method for Identifying Stochastic Information of Heterogeneous Materials |
CN106814017A (zh) * | 2017-01-18 | 2017-06-09 | 浙江科技学院(浙江中德科技促进中心) | 纳米溶胶加固砂土的渗透性与动模量关系测试装置及方法 |
CN109115167A (zh) * | 2018-09-11 | 2019-01-01 | 天津城建大学 | 基于直角四面体角应变的三维应变状态测试装置及方法 |
CN109469039A (zh) * | 2018-11-12 | 2019-03-15 | 浙江科技学院 | 微纤维增强作用下纳米硅溶胶快速渗透加固砂土的方法 |
CN111046608A (zh) * | 2019-12-19 | 2020-04-21 | 东南大学 | 一种玄武岩纤维沥青胶浆有限元模型的建模方法 |
CN111062162A (zh) * | 2019-12-12 | 2020-04-24 | 王靖涛 | 一种岩土材料精确本构模型的数值建模与应用方法 |
CN111368481A (zh) * | 2020-03-29 | 2020-07-03 | 中国石油天然气集团有限公司 | 一种定向穿越孔眼稳定性分析方法 |
CN111783332A (zh) * | 2020-06-12 | 2020-10-16 | 广西大学 | 压硬性非线性变化和剪缩突变特性材料振动累积变形的有限元仿真方法 |
-
2021
- 2021-04-16 CN CN202110413264.XA patent/CN113155612B/zh active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101779204A (zh) * | 2007-08-17 | 2010-07-14 | 宝洁公司 | 广义本构建模方法和*** |
US20100071309A1 (en) * | 2008-05-01 | 2010-03-25 | Brangan Charles D | Soil stabilization system, stabilized soil comprising same, and a method of stabilizing soil |
US20150057988A1 (en) * | 2013-08-26 | 2015-02-26 | The University Of Akron | Method for Identifying Stochastic Information of Heterogeneous Materials |
CN103528897A (zh) * | 2013-10-17 | 2014-01-22 | 三峡大学 | 高围压下粗粒土大三轴试验数据临界状态参数确定方法 |
CN106814017A (zh) * | 2017-01-18 | 2017-06-09 | 浙江科技学院(浙江中德科技促进中心) | 纳米溶胶加固砂土的渗透性与动模量关系测试装置及方法 |
CN109115167A (zh) * | 2018-09-11 | 2019-01-01 | 天津城建大学 | 基于直角四面体角应变的三维应变状态测试装置及方法 |
CN109469039A (zh) * | 2018-11-12 | 2019-03-15 | 浙江科技学院 | 微纤维增强作用下纳米硅溶胶快速渗透加固砂土的方法 |
CN111062162A (zh) * | 2019-12-12 | 2020-04-24 | 王靖涛 | 一种岩土材料精确本构模型的数值建模与应用方法 |
CN111046608A (zh) * | 2019-12-19 | 2020-04-21 | 东南大学 | 一种玄武岩纤维沥青胶浆有限元模型的建模方法 |
CN111368481A (zh) * | 2020-03-29 | 2020-07-03 | 中国石油天然气集团有限公司 | 一种定向穿越孔眼稳定性分析方法 |
CN111783332A (zh) * | 2020-06-12 | 2020-10-16 | 广西大学 | 压硬性非线性变化和剪缩突变特性材料振动累积变形的有限元仿真方法 |
Non-Patent Citations (18)
Also Published As
Publication number | Publication date |
---|---|
CN113155612B (zh) | 2022-09-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105868484B (zh) | 膨胀土地基胀缩作用下无砟轨道低矮路堤顶面升降量计算方法 | |
JP5503653B2 (ja) | 衝撃破砕機のための複合衝撃材 | |
CN109255191B (zh) | 一种定量计算铁路路基堆载预压后所产生沉降量的方法 | |
JP5254978B2 (ja) | カーボンナノチューブを活用した複合焼結材料及びその製造方法。 | |
JP5484468B2 (ja) | 階層複合材料 | |
JP6197571B2 (ja) | ソイルセメントの圧縮強さ推定方法 | |
CN105178126B (zh) | 一种多年冻土区路桥过渡带组合结构及施工方法 | |
Granger et al. | Effect of composition on basic creep of concrete and cement paste | |
CN112203786B (zh) | 具有截头圆锥形***件的复合齿 | |
CN111622260B (zh) | 一种模块式加筋土挡墙极限状态墙面水平位移的确定方法 | |
CN113155612B (zh) | 微纤维混合硅溶胶固化钙质砂的变形预报方法 | |
CN105917057A (zh) | 用于制造混凝土构件的方法、混凝土构件的预制的结构元件以及混凝土构件 | |
Jeong et al. | Densification and compressive strength of in-situ processed Ti/TiB composites by powder metallurgy | |
DE68916515T2 (de) | Verfahren zur Herstellung von Verbundmaterial mit einer Metallmatrix mit kontrollierter Verstärkungsphase. | |
CN109553256A (zh) | 一种快速固化淤泥的方法 | |
Ishihara et al. | Characterization of undrained behaviour of soils in the reclaimed area of Kobe | |
CN113125263B (zh) | 硅溶胶固化非可破碎砂土受力变形的预报方法 | |
CN115828605A (zh) | 一种顶端与沉陷道路刚性连接钢管桩竖向承载力计算方法 | |
CN113125262A (zh) | 可破碎钙质砂受载荷过程中变形的快速预报方法 | |
Habibi et al. | Monotonic behavior of geotextile reinforced soils under discrete rotation of principal stresses | |
Deng et al. | Ageing and viscous effects on the deformation of clay in 1D compression | |
Li | Basic concepts on the compaction of soil | |
CN113008676A (zh) | 硅溶胶加固砂土的检测方法 | |
CN213062106U (zh) | 一种砖砌闸井井壁底部抗滑结构 | |
JP3479802B2 (ja) | 締固め砂杭による砂質地盤の改良工法 |
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 |