附加质量法在土石混合料密度检测中的适用性研究

徐思远; 陈生水; 傅中志; 刘斯宏; 次仁云旦

doi:10.11779/CJGE2023S10022

附加质量法在土石混合料密度检测中的适用性研究 English Version

徐思远^{1, 2,},
陈生水^{1, 2},
傅中志^{1, 2, ,},
刘斯宏³,
次仁云旦⁴

1.
南京水利科学研究院岩土工程研究所, 江苏南京 210024
2.
水利部土石坝破坏机理与防控技术重点实验室, 江苏南京 210024
3.
河海大学水利水电学院, 江苏南京 210024
4.
西藏自治区水利电力规划勘测设计研究院, 西藏拉萨 850004

基金项目:

国家重点研发计划资助项目 2021YFC3090101

国家自然科学基金面上项目 52279135

江苏省卓越博士后计划项目 2022084

水力发电工程学会青年人才托举工程项目水电学秘字〔2023〕15号

南京水利科学研究院青年基金项目 Y323004

详细信息

作者简介:
徐思远(1993—)，男，工程师，主要从事土石坝渗漏病害探测、岩土原位测试等方面的研究工作。E-mail: syxu@nhri.cn

通讯作者:
傅中志, E-mail: fu_zhongzhi@yahoo.com

中图分类号: TU43
计量
- 文章访问数: 269
- HTML全文浏览量: 37
- PDF下载量: 111
出版历程
- 收稿日期: 2023-07-04
- 网络出版日期: 2023-11-23
- 刊出日期: 2023-10-31

Applicability of additive mass method in density tests on soil-mixture materials

XU Siyuan^{1, 2,},
CHEN Shengshui^{1, 2},
FU Zhongzhi^{1, 2, ,},
LIU Sihong³,
CIREN Yundan⁴

1.
Geotechnical Engineering Department, Nanjing Hydraulic Research Institute, Nanjing 210024, China
2.
Key Laboratory of Failure Mechanism and Safety Control Techniques of Earth-Rock Dam of the Ministry of Water Resources, Nanjing 210024, China
3.
College of Water Conservancy and Hydropower Engineering, Hohai University, Nanjing 210024, China
4.
Tibet Hydraulic Waterpower Planning and Designing Institution, Lasha 850000, China

摘要

摘要: 附加质量法已经在多座土石坝工程的不同筑坝料密度检测中得到了应用，随着近年来水利水电工程的发展，宽级配高细粒含量的土石混合料也逐渐被应用在坝体填筑中，而附加质量法在土石混合料中的适用研究开展的仍较少。依托某抽水蓄能电站工程，对不同风化玢岩含量的混合填筑料开展了附加质量法试验和坑测密度试验，并对附加质量法密度求解的准确性进行了讨论和分析。经过研究发现：附加质量法在玢岩含量较低的试料中表现良好，但随着混合料玢岩比例的提高，试料的实测干密度以及附加质量法测得的频率和刚度均表现出一定的随机性，部分测点结果的合理性较差，需要开展 ${\omega ^{ - 2}}{\text{ - }}\Delta m$ 相关性分析来验证测点结果合理性才能取得准确的密度结果。
- 土石混合料 /
- 附加质量法 /
- 压实密度检测 /
- 刚度相关法 /
- 压实质量
Abstract: The additive mass method (AMM) has been applied in the density tests on different materials of many earth-rockfill dams. With the development of hydropower projects in recent years, the soil-rock mixtures (S-RM) with wide gradation and high fines content have been used for dam construction. While the researches on the applicatbility of the AMM in the soil-rock mixtures are insufficient. In this study, the AMM is adopted to evaluate the compaction densities of the S-RM composing of fresh dolomite and strongly weathered diorite porphyrite in a pumped storage power station in China. The compaction densities evaluated by the AMM are compared with those measured by the pit tests. The results demonstrate that for the S-RM containing low contents of diorite porphyrite, the compaction densities can be reasonably evaluated by the AMM. For the S-RM with high contents of diorite porphyrite, the analysis of the ${\omega ^{ - 2}}{\text{ - }}\Delta m$ correlation obtained by the AMM in advance is required to conduct reasonable compaction density evaluation.
- soil-rock mixture /
- additive mass method /
- compaction density evaluation /
- stiffness-density correlation /
- compaction quality

HTML全文

0. 引言

城市排水管网肩负着城市污水、雨水排放的重任,是现代城市重要的基础设施。但中国大部分市政排水管道铺设于20世纪80年代以前,随着使用年限的临近,管道老化严重,管道下部脱空、沉降、管道腐蚀、渗漏、断裂等灾害不断发生,导致管道无法正常使用^[1],更严重的是因此产生的地面塌陷,地下水污染等事故频发。排水管道的传统修复方式主要是现场整体开挖进行局部修复或管道整体更换,该方式成本昂贵且对城市道路及周围居民生活造成严重影响。高聚物注浆修复技术是一种高效的非开挖修复技术,通过向管道结构脱空渗漏处注射高聚物注浆材料,利用非水反应类高聚物材料本身的膨胀性及固化后的高黏结性,达到填充脱空、封堵渗漏和抬升沉降管道的目的（如图1,2所示）,目前该技术已成功应用于多种地下管道修复项目中^[2-3]。

图 1 高聚物注浆填充管道脱空

Figure 1. Polymer grouting fills the void outside pipe

下载: 全尺寸图片幻灯片

图 2 高聚物注浆抬升管道沉降

Figure 2. Polymer grouting lifts the settlement pipe

下载: 全尺寸图片幻灯片

对于高聚物修复管道的力学性能,目前多利用有限元软件进行模拟或进行试验探究,对修复后管道计算理论和方法的研究相对不足。如徐建国等^[4]采用ABAQUS软件计算了地下管道高聚物注浆修复后的受力情况,结果表明高聚物注浆后,管道内最大拉压应力和最大竖向位移均明显减小；李斌等^[5]对修复后管道在车辆荷载下的力学性能进行了有限元分析,结果显示管道经高聚物修复后,其受力变形恢复到了正常管道水平；张倩倩^[6]利用有限元软件分析了高聚物浆液膨胀扩散对管道应力应变的影响；Wang等^[7]进行了高聚物注浆修复钢筋混凝土排水管道在冲击荷载与车辆荷载下的足尺试验,结果表明高聚物注浆可以有效地修复脱空管道。

通常对埋地管道进行计算分析时,需要考虑管与周围介质的相互作用。国内外学者对此进行了大量研究,如Kjartanson等^[8]通过将管道模化为一系列梁单元,管周的填充土模化为一个非线性材料,从而提出了不开槽埋设管道的管-土相互作用有限元分析法；Lee等^[9]利用有限元模型来反映土壤特性、回填和原位条件,提出了考虑管-土相互作用的地下管道水平挠度方程；Jung等^[10]利用Mohr-Coulomb屈服准则建立了粒状土壤中管-土相互作用的有限元模型；可文海等^[11]利用Pasternak模型模拟管-土相互作用,计算了盾构隧道下穿管道施工引起的管道处土体竖向位移；王小龙等^[12]基于Winkler地基梁理论,建立了埋地钢管在局部悬空时管-土相互作用的力学模型；刘全林等^[13]基于地埋管道受力特征的实测结果,建立了地埋管道与土的相互作用分析组合模型,并给出了模型参数的确定方法。上述文献研究多集中于管-土相互作用,而对于高聚物注浆修复管道则需要考虑管-土-高聚物之间的相互作用。因此本文基于弹塑性理论建立考虑管-土-高聚物相互作用的Winkler模型和Vlazov模型,采用传递矩阵法,编制了MATLAB计算程序,经对地下管道高聚物注浆修复前后的受力变形特性的对比计算,并通过参照现场试验结果,探讨了Winkler模型和Vlazov模型的准确性和合理性。最后采用Vlazov模型分析计算了多工况下管道脱空高聚物注浆修复后的受力变形特性。

1. 管-土-高聚物相互作用分析模型

1.1 管-土-高聚物相互作用Winkler模型

Winker所提出的土介质理想化弾性模型是假设土介质表面任一点处的位移 $ω$ 与作用在该点的应力 $q$ 成正比而与作用在其他各点的应力无关^[14],即

$q (x, y) = k ω (x, y)$ ,

(1)

式中, $q (x, y)$ 为作用在土介质表面的压力, $ω (x, y)$ 为土介质表面的沉降, $k$ 为基底反力系数。

方程（1）是Winker模型的响应函数。该模型将土介质理想化为一系列各自独立的弹簧单元,其弹簧常数为 $k$ 。将土与高聚物都简化为弹性材料即可得出管-土-高聚物相互作用Winkler模型,如图3所示。

图 3 管-土-高聚物相互作用Winkler模型

Figure 3. Winkler model for pipe-soil-polymer interaction

下载: 全尺寸图片幻灯片

对于其中的模型参数 $k_{1}$ ,可按地基承载板试验结果确定,或根据加州承载比试验结果确定。而模型参数 $k_{2}$ 可通过求解无限弹性介质中圆形边界上均匀径向位移,从而求得基底反力系数 $k_{2}$ 。

对于平面应变的轴对称问题,取应力函数^[15]为 $Φ = A \ln r + B r^{2} \ln r + C r^{2} + D$ ,当不计体力时则有

$\begin{array}{l} σ_{r} = \frac{1}{r} \frac{\partial Φ}{\partial r} = \frac{A}{r^{2}} + B (1 + 2 \ln r) + 2 C, \\ σ_{θ} = \frac{\partial^{2} Φ}{\partial r^{2}} = - \frac{A}{r^{2}} + B (3 + 2 \ln r) + 2 C, \\ τ_{r θ} = 0 。 \end{array}}$

(2)

式中, $Φ$ 为应力系数, $A, B, C, D$ 为系数, $r$ 为半径, $σ_{r}$ , $σ_{θ}$ 分别为径向应力与切向应力。

由位移单值性可得 $B = 0$ ,其次在距离原点无穷远处,由圣维南原理可知 $σ_{r} = σ_{θ} \overset{}{\to} 0$ 从而可得 $C = 0$ 。将 $B = 0$ , $C = 0$ 代入式（2）化简可得

$\begin{array}{l} σ_{r} = \frac{1}{r} \frac{\partial Φ}{\partial r} = \frac{A}{r^{2}}, \\ σ_{θ} = \frac{\partial^{2} Φ}{\partial r^{2}} = - \frac{A}{r^{2}} 。 \end{array}}$

(3)

由平面应变的胡克定律,并将式（3）代入可得

$\begin{array}{l} ε_{r 1} = \frac{1}{E_{g}} (σ_{r} - ν_{g} σ_{θ}) = \frac{A (1 + ν_{g})}{E_{g} r^{2}}, \\ ε_{r 2} = \frac{1}{E_{s}} (σ_{r} - ν_{s} σ_{θ}) = \frac{A (1 + ν_{s})}{E_{s} r^{2}} 。 \end{array}}$

(4)

式中, $ε_{r 1}$ , $ε_{r 2}$ 分别为高聚物层与土层的应变, $E_{g}$ , $ν_{g}$ 分别为高聚物的弹性模量与泊松比, $E_{s}$ , $ν_{s}$ 分别为土的弹性模量与泊松比。

圆形边界上均匀径向位移可为

$u = \int_{R}^{R + D} ε_{r 1} d r + \int_{R + D}^{\infty} ε_{r 2} d r$ ,

(5)

式中,D为高聚物层厚度。

将式（4）代入式（5）从而得到径向位移：

$u = \frac{A}{R + D} (\frac{1 + ν_{s}}{E_{s}} - \frac{1 + ν_{g}}{E_{g}}) + \frac{A (1 + ν_{g})}{E_{g} R}$ 。

(6)

按Winker模型中弹性参数 $k$ 的定义,将式（3）,（6）代入,则得到高聚物层的基底反力系数为

$k_{2} = \frac{σ_{r}}{u} = \frac{1}{\frac{R^{2}}{R + D} (\frac{1 + ν_{s}}{E_{s}} - \frac{1 + ν_{g}}{E_{g}}) + \frac{R (1 + ν_{g})}{E_{g}}}$ 。

(7)

1.2 管-土-高聚物相互作用Vlazov模型

Winkler模型由于在描述自然土体的连续性态中所固有的缺陷以及Winkler模型只适用于地基土（压缩层）很薄的情况,而当压缩层较厚时就会产生较大误差。为了改进Winkler模型的上述缺陷,通过在弹性连续介质模型中引入约束或简化位移分布与应力的某些假设,从而得到一种新的模型,称为双参数模型^[14]。Vlazov模型是一种双参数模型,它所列出的模型公式是以变分法的应用为基础,通过引进一些能简化各向同性线弹性连续介质基本方程的位移约束而导出的。其响应函数为

$q (x, y) = K ω (x, y) - T \nabla^{2} ω (x, y)$ ,

(8)

式中, $q (x, y)$ 为作用在土介质表面的压力, $ω (x, y)$ 为土介质表面的沉降, $\nabla^{2}$ 为拉普拉斯算子, $K$ 为压应力作用下对土介质变形的度量,T为对邻近土单元体传递性的度量。

该模型也适用于非均匀的弹性层,将土与高聚物都简化为弹性材料即可得出管-土-高聚物相互作用分析Vlazov模型,如图4所示。

图 4 管-土-高聚物相互作用Vlazov模型

Figure 4. Vlazov model for pipe-soil-polymer interaction

下载: 全尺寸图片幻灯片

可通过求解双参数模型的地基控制方程来求得模型参数。弹性层在平面应变状态下,在 $x$ - $z$ 平面内假设土层中任一点M的位移满足

$\bar{u} (x, z) = 0, \bar{ω} (x, z) = ω (x) φ (z), [φ (0) = 1, φ (H) = 0]$ ,

(9)

式中, $φ (z)$ 为描述 $z$ 方向位移 $ω (x, z)$ 的变化函数。

利用平面应变条件下的应力与应变关系可得

$\begin{array}{l} σ_{z 1} = \frac{E_{1}}{1 - ν_{1}^{2}} ω (x) φ^{'} (z), \\ τ_{x z 1} = \frac{E_{1}}{2 (1 + ν_{1})} ω^{'} (x) φ (z), \\ σ_{z 2} = \frac{E_{2}}{1 - ν_{2}^{2}} ω (x) φ^{'} (z), \\ τ_{x z 2} = \frac{E_{2}}{2 (1 + ν_{2})} ω^{'} (x) φ (z) 。 \end{array}}$

(10)

式中 $E_{1} = E_{g} / (1 - ν_{g}^{2})$ , $ν_{1} = ν_{g} / (1 - ν_{g})$ ； $E_{2} = E_{s} /$ $(1 - ν_{s}^{2})$ , $ν_{2} = ν_{s} / (1 - ν_{s})$ ； $E_{g}$ , $ν_{g}$ 为高聚物的弹性模量与泊松比, $E_{s}$ , $ν_{s}$ 为土的弹性模量与泊松比。

利用拉格朗日虚功原理,取虚位移 $δ \bar{u} = 0$ , $δ \bar{ω} =$ $φ (z) δ ω$ ,得到外力所作的虚功为

$\begin{matrix} U_{e} = b q (x) φ (0) δ ω d x + b \int_{0}^{H_{1}} \frac{\partial τ_{x z 1}}{\partial x} φ (z) δ ω d x d x + \\ b \int_{H_{1}}^{H} \frac{\partial τ_{x z 2}}{\partial x} φ (z) δ ω d x d z \end{matrix}$ 。

(11)

内力所作的虚功为

$U_{i} = - b \int_{0}^{H_{1}} σ_{z 1} φ^{'} (z) δ ω d x d z - b \int_{H_{1}}^{H} σ_{z 2} φ^{'} (z) δ ω d x d z$ 。

(12)

由 $U_{e} + U_{i} = 0$ ,代入式（10）并消去 $δ ω d x$ 可得

$q (x) = K ω (x) - 2 T \frac{d^{2} ω (x)}{d x^{2}}$ 。

(13)

式中, $K = \frac{E_{1}}{1 - ν_{1}^{2}} \int_{0}^{H_{1}} {[φ^{'} (z)]}^{2} d z + \frac{E_{2}}{1 - ν_{2}^{2}} \int_{H_{1}}^{H} {[φ^{'} (z)]}^{2} d z$ , $T =$ $\frac{E_{1}}{4 (1 + ν_{1} ）} \int_{0}^{H_{1}} φ^{2} （ z) d z + \frac{E_{2}}{4 (1 + ν_{2} ）} \int_{H_{1}}^{H} φ^{2} （ z) d z$ 。

对于函数 $φ (z)$ ,Vlazov提出了线性与非线性表达式分别为

$φ_{1} (z) = 1 - z / H, φ_{2} （ z) = \frac{sh [γ (H - z) / L]}{sh [γ h / L]} 。$

式中 H为计算弹性层厚度,L为地埋管道计算长度, $γ$ 为反映基础梁特性的常数,通常取1.5。本文中对于高聚物层,取线性表达式,对于土层,采用非线性表达式。将其代入K,T计算表达式可得模型参数如下：

$K_{1} = \frac{E_{2}}{H {(1 - ν_{2})}^{2}} \cdot \frac{γ H}{2 L} (\frac{sh (γ H / L) ch (γ H / L) + γ H / L}{{sh}^{2} (γ H / L)}),$

$T_{1} = \frac{E_{2} H}{12 (1 + ν_{2})} \cdot \frac{3 L}{2 γ H} (\frac{sh (γ H / L) ch (γ H / L) - γ H / L}{{sh}^{2} (γ H / L)}),$

$K_{2} = \frac{E_{1}}{H_{1} {(1 - ν_{1})}^{2}} + \frac{E_{2}}{H {(1 - ν_{2})}^{2}} \cdot \frac{γ H}{2 L} \cdot [sh (γ H / L - γ H_{1} / L) \cdot$

$ch (γ H / L - γ H_{1} / L) + γ H / L - γ H_{1} / L] / {sh}^{2} (γ H / L),$

$T_{2} = E_{1} H_{1} / 12 (1 + ν_{1}) + E_{2} H / 12 (1 + ν_{2}) \cdot 3 L / 2 γ H \cdot [sh (γ H / L -$

$γ H_{1} / L) \cdot ch(γ H / L - γ H_{1} / L) - γ H / L + γ H_{1} / L)] / {sh}^{2} (γ H / L) 。$

2. 基于传递矩阵法的管道平面应变问题分析

2.1 传递矩阵法

传递矩阵法是以受力杆件的基本微分方程为基础,采用矩阵相乘来求出杆件各截面的内力和位移值,从而对杆件进行力学分析的一种方法^[16]。其计算过程为

$S_{i} = U_{i} U_{i - 1} \dots U_{2} U_{1} S_{0}$ ,

(14)

式中, $S_{0}$ 为初始状态矢量, $S_{i}$ 为末端状态矢量, $U_{1}$ , $U_{2}$ , $\dots$ , $U_{i}$ , $U_{i - 1}$ 为各段的传递矩阵。

2.2 管道计算模型描述

为计算简便,将管道简化为平面应变问题,忽略管道横截面的不圆度,将其假设为圆环形,并且考虑作用在管道上的荷载为对称荷载,因此只需对管道的一半进行计算分析。如图5所示,取单位长度管道,管道壁厚为t,半径为R,将管道的一半划分为n个单元,同时由于管道的不同位置,管土的变形趋势不同并且不同位置土体的压实程度也会有区别,所以与管道相接触的高聚物和土体也相应的被划分为三部分,分别为：底部支承区、胸腔区、管顶区。

图 5 管道计算模型

Figure 5. Computational model for buried pipe

下载: 全尺寸图片幻灯片

将划分的n个单元都假设为梁单元,并按顺序刚性连接,取第 $i$ 单元作为研究对象,其受力简图如图6所示。

图 6 单元受力简图

Figure 6. Forces acting on elements

下载: 全尺寸图片幻灯片

图6中, $p_{i}$ , $q_{i}$ 为 $i$ 节点处的径向荷载和切向荷载, $p_{i + 1}$ , $q_{i + 1}$ 为 $i + 1$ 节点处的径向荷载和切向荷载, $l$ 为单元 $i$ 的长度, $θ$ 为 $i$ 单元与总体坐标系中 $X$ 轴的夹角。

设在外力作用下,管道 $i$ 单元中 $i$ 节点处的切向位移为 $u_{i}$ ,径向位移为 $v_{i}$ ,转角为 $ω_{i}$ ,弯矩为 $M_{i}$ ,轴力为 $N_{i}$ ,剪力为 $T_{i}$ 。在局部坐标系中,由平衡条件和物理方程可得

$\begin{array}{l} \frac{d T}{d x} = - K^{'} y + 2 T \frac{d^{2} y (x)}{d x^{2}} - p_{i} (1 - \frac{x}{l}) - p_{i + 1} \frac{x}{l}, \\ \frac{d N}{d x} = q_{i} (1 - \frac{x}{l}) + q_{i + 1} \frac{x}{l}, \\ \begin{matrix} \frac{d M}{d x} = T & , & \frac{d u}{d x} = - \frac{N}{E S}, \end{matrix} \\ \begin{matrix} \frac{d^{2} v}{d x^{2}} = \frac{M}{E I} & , & \frac{d v}{d x} = ϖ 。 \end{matrix} \end{array}}$

(15)

求解式（15）,并满足边界条件：当 $x = 0$ 时, $u = u_{i}$ , $ν = ν_{i}$ , $ω = ω_{i}$ , $M = M_{i}$ , $N = N_{i}$ , $T = T_{i}$ ；当 $x = l$ 时, $u = u_{i + 1}$ , $ν = ν_{i + 1}$ , $ω = ω_{i + 1}$ , $M = M_{i + 1}$ , $N = N_{i + 1}$ , $T = T_{i + 1}$ 。可得

${S_{i + 1}} = [U_{i, i + 1}] {S_{i}}$ ,

(16)

式中, ${S_{i}} = {u_{i}, ν_{i}, ω_{i}, M_{i}, N_{i}, T_{i}, 1}$ 为 $i$ 节点处的内力和位移列阵, $[U_{i, i + 1}]$ 为 $7 \times 7$ 的传递矩阵, ${S_{i + 1}} =$ ${u_{i + 1}, ν_{i + 1}, ω_{i + 1}, M_{i + 1}, N_{i + 1}, T_{i + 1}, 1}$ 为 $i + 1$ 节点处的内力和位移列阵。

通过坐标转换矩阵 $R$ 可将式（16）变为整体坐标系下的传递关系：

${{\bar{S}}_{i + 1}} = [{\bar{U}}_{i, i + 1}] {{\bar{S}}_{i}} 。$

(17)

式中 $[{\bar{U}}_{i, i + 1}] = R [U_{i, i + 1}] R^{- 1}$ ,为 $i$ 单元的传递矩阵, ${{\bar{S}}_{i}}$ 、 ${{\bar{S}}_{i + 1}}$ 分别为整体坐标系下 $i$ 节点与 $i + 1$ 节点处的内力和位移列阵, $R$ 为 $i$ 单元的坐标转换矩阵。已知各单元的传递矩阵后,将其代入式（14）可得

${{\bar{S}}_{n}} = [{\bar{U}}_{n, n + 1}] [{\bar{U}}_{n - 1, n}] \dots [{\bar{U}}_{2, 3}] [{\bar{U}}_{1, 2}] {{\bar{S}}_{1}}$ 。

(18)

由 $u_{1}$ = $ω_{1}$ = $T_{1}$ , $u_{n}$ = $ω_{n}$ = $T_{n}$ ,代入式（18）即可求得 $ν_{1}, M_{1}, N_{1}$ 和 $ν_{n}, M_{n}, N_{n}$ 。

再以 $1$ 节点为已知量,代入

${{\bar{S}}_{i}} = [{\bar{U}}_{i, i + 1}] [{\bar{U}}_{i - 1, i}] \dots [{\bar{U}}_{2, 3}] [{\bar{U}}_{1, 2}] {{\bar{S}}_{1}}$ ,

(19)

即可求得任意节点的内力和位移。求出节点的位移后,可以确定各单元相应的接触压力。结合土体的弹-塑性变形法则和管道本身的变形量,部分接触压力将被限制以满足土的塑性条件,然后经迭代计算使之满足所有的边界条件和变形准则。

2.3 程序实现

根据上述计算过程,利用MATLAB编制了管道应力计算程序,其计算流程图如图7所示。

图 7 管道应力计算程序流程图

Figure 7. Flow chart of stress program for pipe

下载: 全尺寸图片幻灯片

3. 试验概述及数据处理

管道试验在室外现场区域进行,开挖长、宽、深为8 m×2.3 m×2.3 m的试验坑。选择3个完整的管段和两个半管段,有效长度为8 m。管道放置到沟槽中后,开始回填,每层厚约0.2 m,用小夯夯实土壤。然后对同一根管道先后进行正常、脱空、注浆修复等工况的试验。通过在管内外壁黏贴应变片来测量荷载作用下管道的内外应变。试验所用管道、土壤和高聚物性能参数如表1所示,试验管道剖面图及现场试验图如图8所示。

表 1 试验所用管道、土壤和高聚物性能参数表

Table 1. Performance parameters of pipes, soil and polymer

材料	密度/(kg·m^-3)	弹性模量/MPa	泊松比	黏聚力/kPa	内摩擦角/(°)
管道	2400	30000	0.3	—	—
填土	1980	60	0.3	40	30
高聚物	240	20	0.2	—	—

下载: 导出CSV

| 显示表格

图 8 试验管道剖面图及现场试验图

Figure 8. Sectional drawing of test pipes and pictures of field tests

下载: 全尺寸图片幻灯片

为了将试验值与模型计算结果进行更好的对比,需要对试验数据进行一定处理。根据文献[17]可知管道弯矩与内外管壁应变之间的关系如下：

$M = E I (\frac{ε_{in} - ε_{out}}{h}),$

(20)

式中,M为管道弯矩,E为管道弹性模量,I为截面惯性矩, $ε_{in}$ , $ε_{out}$ 分别为内外管壁应变,h为管壁厚度。

在管道设计中,管道自身承受外部荷载的能力通过管道刚度来定义,其与管壁的几何形状和管道材料类型有关。同时在ASTM D 2412—11^[18]中,管道刚度还被定义为：使管道产生一定百分比变形（即5%）的荷载与沿荷载施加方向管道内径变形量的比值,因此已知管道刚度和荷载大小可以根据

$Δ y = \frac{F}{P S}$

(21)

近似求出相应的管道径向变形。式中, $Δ y$ 为管道径向变形, $F$ 为管壁上的径向应力, $P S$ 为管道刚度。

4. 模型计算及与试验结果对比

依据试验情况,取静荷载作用管段为计算管道,采用不同的管-土-高聚物相互作用模型,其计算结果如图9,10所示。图9（a）,（b）分别表示正常、脱空和高聚物修复情况下,采用Winkler模型计算所得的管道径向变形和环向轴力对比图,同样的图10（a）,（b）表示采用Vlazov模型计算所得的结果图。从图9,10可知,修复后管道的径向变形会略大于正常值,而轴力与正常值基本一致,正常与高聚物修复后的径向变形最大差值都出现在管道底部,分别为0.4 mm（采用Winkler模型）与1 mm（采用Vlazov模型）。由此可见,高聚物注浆后管道的变形与力学性能可以恢复到正常水平。在脱空条件下,管道底部的径向变形与管顶的轴力会出现明显的增大现象,与正常条件相比,底部的径向变形增大了3.9倍（采用Winkler模型）与5.3倍（采用Vlazov模型）,同样的,顶部的环向轴力分别增大了1.6倍和2.0倍。说明底部脱空对于管道的力学特性会有很大的影响。

图 9 Winkler模型计算结果

Figure 9. Calculated results of Winkler model

下载: 全尺寸图片幻灯片

图 10 Vlazov模型计算结果

Figure 10. Calculated results of Vlazov model

下载: 全尺寸图片幻灯片

表2为正常与修复条件下管-土-高聚物相互作用分析模型参数的比较。由表2可知,Winkler模型计算出的修复后地基反力系数为正常值的66%,而Vlazov模型达到了正常值的91%,说明Vlazov模型计算出的修复后结果会与正常值更为接近。图11给出了采用不同分析模型得到的径向变形和环向弯矩与试验值的对比。由图11（a）,（b）可以看出,模型计算得出的径向变形和环向弯矩变化趋势图与试验值基本一致,与试验值相比最大径向变形差值都为2 mm,误差较小,说明建立的计算方法是可靠的。并且Vlazov模型计算结果明显比Winkler模型计算结果大,这是由于Vlazov模型增加了土体的侧向传力作用,考虑了土介质的荷载传递性,比Winkler模型合理。由图11（a）可知,采用Vlazov模型计算的径向变形与试验值较接近,除管顶处,计算误差在2%～25%,在45°管腰处,模型计算结果与试验值最接近,Winkler模型计算结果与试验值误差较大,达9%～55%。由图11（b）可知,在环向弯矩计算中,Winkler模型计算结果与试验值的误差在管底、管侧、管顶分别为35%,21%,12%,同理Vlazov模型与试验值的误差分别为18%,16%,11%,因此,采用Vlazov模型计算的环向弯矩与试验值较接近。综上所述,相比于Winkler模型,Vlazov模型与管道实际受力情况符合得更好。

表 2 相互作用分析模型参数

Table 2. Parameters of interaction analysis model

模型参数/(N·m^-3)	底部支承区（修复）	底部支承区（正常）	$\frac{修复}{正常} / %$
Winkler模型K	68945	105230	66
Vlazov模型K	82842	91434	91
Vlazov模型T	66	91	73

下载: 导出CSV

| 显示表格

图 11 模型计算结果与试验对比

Figure 11. Comparison between model and test results

下载: 全尺寸图片幻灯片

5. 不同因素对修复后管道内力和变形的影响

埋地管道的内力和变形受到管道埋深、管径、管道材料、地面荷载及周围填土的性质等诸多因素的影响,同时由于利用Vlazov模型来模拟管-土-高聚物的相互作用更加符合实际情况,因此下文将采用Vlazov模型对管道脱空高聚物注浆修复后不同工况下的管道进行计算。

5.1 不同埋深的影响

图12分别为覆土深度是1.0,1.5,2 m时修复后管道的径向变形和环向轴力对比图。从图中可以看出,不同埋深时,径向变形和轴力的变化趋势是一致的,随着埋深增大,其数值会明显增大。在图12（a）中,由于选用坐标系的不同,虽然管道底部变形为正值,顶部为负值,但二者都是向下的,顶部变形量会大于底部变形量。而对于管道的环向弯矩和轴力,从图12（b）可看出,轴力分布在管顶处呈压缩状态,在管腰下部呈拉伸状态,且轴向最大拉应力大于最大压应力。

图 12 不同埋深下管道的计算值

Figure 12. Calculated values of pipes at different buried depths

下载: 全尺寸图片幻灯片

5.2 不同附加荷载的影响

图13为管顶施加不同大小荷载时,修复后管道的径向变形和环向轴力对比图。随着管道附加荷载的增加,管顶处的径向变形以及管腰处的轴力都出现了明显的增大。从图13（a）可知,荷载为0.1,0.15 MPa时管底的径向变形分别为0.05 MPa时管底径向位移的1.1倍和1.2倍；其管顶的径向变形分别为0.05 MPa时管顶径向位移的1.3倍和1.5倍,由此可见,同样荷载变化下,管顶的变形变化率要大于管底,这是因为相对于管顶,管底的埋深更大,而随着土壤厚度的增加,与土压力相比,附加载荷对管道的影响会变小,同时管土相互作用加强,这也在一定程度上限制了管道的变形。

图 13 不同荷载下管道的计算值

Figure 13. Calculated values of pipes under different loads

下载: 全尺寸图片幻灯片

5.3 不同土壤弹性模量的影响

图14为改变回填土的弹性模量时,修复后管道的径向变形和环向弯矩对比图。从图14（a）可知,在管腰以下,高聚物注浆修复区的径向变形随回填土弹性模量的增大而增大,在非脱空区则呈现相反趋势,图14（b）中的弯矩也有相同的变化趋势。由此可见,非脱空区的管土相互作用随着土体弹性模量的增大而加强,但在高聚物注浆修复区,由于高聚物的加入,与原来的土层一起形成了复合土层,管道与复合土层之间的相互作用则随着土体弹性模量的增大而减弱。从图14（b）可以看出,在非脱空区,管腰下部和管顶处的弯矩呈相反的变化趋势,这是因为,管道内力随回填土弹性模量的变化而变化是两个方面效应共同作用的结果^[19],增大管周土体弹性模量一方面增加了土体对管道的弹性抗力,增强了管土相互作用,使管壁内力分布更加均匀,从而降低了管道内力,但另一方面也增大了管道基础的刚度导致了更大的管顶竖向应力,这会使管道内力随之增大。在管顶区,由于埋深较浅,管顶竖向应力起主导作用,因此图14（b）中管顶区的弯矩随回填土弹性模量的增大而增大,而在管腰以下,埋深相对较深,管土相互作用起主导作用,弯矩也就随之减小。

图 14 不同土壤弹性模量下管道的计算值

Figure 14. Calculated values of pipes under different elastic moduli of soil

下载: 全尺寸图片幻灯片

5.4 不同管壁厚度的影响

图15为不同管壁厚度下,修复后管道的径向变形和环向弯矩对比图。从图中可以看出,管道的环向弯矩随着管壁厚度的增大而增大,径向变形则呈相反趋势。这是由于,同一直径管道,管壁厚度不同,管道的抗弯刚度也就不同,增加管壁厚度实际就是增加了管道的截面刚度,从而导致管道径向变形的减小,但刚度的增大也会使管土之间的变形协调能力降低,使得管顶竖向应力增大,从而增大了管周的应力分布,因此管道的弯矩也随之增大。

图 15 不同管壁厚下管道的计算值

Figure 15. Calculated values of pipes with different pipe thicknesses

下载: 全尺寸图片幻灯片

6. 结论

为了对高聚物注浆修复脱空渗漏的地下管道开展受力变形特性分析,本文基于弹塑性理论分别建立了考虑管-土-高聚物相互作用的Winkler模型和Vlazov模型,采用传递矩阵法,编制了MATLAB计算程序,经对地下管道高聚物注浆修复前后的受力变形特性的对比计算,并通过参照现场试验结果,探讨了Winkler模型和Vlazov模型的准确性和合理性。最后采用Vlazov模型开展了管道脱空高聚物注浆修复后受力变形特性研究,工况包括管道埋深、荷载集度、回填土弹性模量、管道壁厚等因素,得到以下6点结论。

（1）考虑管-土-高聚物相互作用的Winkler模型和Vlazov模型计算得出的径向变形和环向弯矩计算值与试验值基本一致,验证了模型的准确性和合理性,Vlazov模型计算结果与试验值更为接近。

（2）管道底部的脱空对管道的整体受力变形力学特性有很大影响,底部脱空区使管道的径向变形和轴向内力明显增大。但经高聚物注浆修复后,管道的径向变形和轴向内力接近正常管道,说明高聚物注浆具有较好的修复管道功能。

（3）随着管道埋深和附加荷载的增大,管道的径向变形、环向弯矩和环向轴力都有不同程度的增大。

（4）管道的径向变形随回填土弹性模量的增大而减小,随着管道周围填充土弹性模量的增大,管土相互作用也随之加强,更有利于发挥土介质对外荷载的支承作用,从而对管道的整体受力有利。

（5）管道的环向弯矩和轴力随着管壁厚度的增大而增大,径向变形则呈相反趋势,管壁厚度的增加使得管道截面刚度增大,从而导致管道径向变形的减小。

（6）高聚物注浆修复技术是一种高效的非开挖修复技术,通过向管道结构脱空渗漏处注射高聚物注浆材料,利用非水反应类高聚物材料本身的膨胀性及固化后的高黏结性,达到填充脱空、封堵渗漏和抬升沉降管道的目的,目前该技术已成功应用于地下管道的渗漏脱空修复工程。本研究为管道高聚物注浆修复设计施工提供了一定的理论参考。

图 1 附加质量法的假定模型

Figure 1. Assumed model for additive mass method

下载: 全尺寸图片幻灯片

图 2 4种混合料的级配曲线

Figure 2. Grain-size distribution curves of four different mixtures

下载: 全尺寸图片幻灯片

图 3 现场坑测密度的概率分布以及其均值随玢岩含量的变化

Figure 3. Distribution of compaction densities obtained from field tests and their variation ranges

下载: 全尺寸图片幻灯片

图 4 现场刚度结果的概率分布及其均值随玢岩含量的变化

Figure 4. Distribution of stiffness obtained from field tests and their variation ranges

下载: 全尺寸图片幻灯片

压实干密度 ${\rho _{\text{d}}}$ -刚度参数 $K/r$ 关系图

${\rho _{\text{d}}}$ plotted against $K/r$ and fitting lines

下载: 全尺寸图片幻灯片

图 6 压实密度求解值-坑测值

Figure 6. Compaction density: evaluated from AMM vs. measured by pit tests

下载: 全尺寸图片幻灯片

${\omega ^{ - 2}}$ - $\Delta m$ 相关关系决定系数 ${r^2}$

Determination coefficient ${r^2}$ between ${\omega ^{ - 2}}$ and $\Delta m$

下载: 全尺寸图片幻灯片

图 8 压实密度求解值-坑测值（筛后DC-25，DC-33测点）

Determination coefficient ${r^2}$ between ${\omega ^{ - 2}}$ and $\Delta m$ (DC-25 and DC-33 Screened)

下载: 全尺寸图片幻灯片

参考文献(11)

[1]	朱晟. 粗粒筑坝材料现场压实质量的控制标准研究[J]. 水力发电, 2011, 37(12): 22-26, 30. https://www.cnki.com.cn/Article/CJFDTOTAL-SLFD201112008.htm ZHU Sheng. Study on on-site compaction quality control standard of coarse-grained materials[J]. Water Power, 2011, 37(12): 22-26, 30. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SLFD201112008.htm
[2]	李丕武. 地基承载力动测的附加质量法[J]. 地球物理学报, 1993, 36(5): 683-687. https://www.cnki.com.cn/Article/CJFDTOTAL-DQWX199305015.htm LI Peiwu. Additive mass method measured dynamically foundation bearing capacity[J]. Chinese Journal of Geophysics, 1993, 36(5): 683-687. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-DQWX199305015.htm
[3]	谭峰屹, 姜志全, 李仲秋, 等. 附加质量法在昆明新机场填料压实密度检测中的应用研究[J]. 岩土力学, 2010, 31(7): 2214-2218. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201007033.htm TAN Fengyi, JIANG Zhiquan, LI Zhongqiu, et al. Application of additive mass method to testing compacted density of filling material in Kunming new airport[J]. Rock and Soil Mechanics, 2010, 31(7): 2214-2218. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201007033.htm
[4]	薛云峰, 崔琳, 李丕武, 等. 刚度相关法测试堆石体密度[J]. 岩土工程学报, 2010, 32(6): 987-990. http://www.cgejournal.com/cn/article/id/13441 XUE Yunfeng, CUI Lin, LI Peiwu, et al. Density measurement of rockfill by use of rigidity correlation method[J]. Chinese Journal of Geotechnical Engineering, 2010, 32(6): 987-990. (in Chinese) http://www.cgejournal.com/cn/article/id/13441
[5]	张智, 凡友华, 蔡加兴, 等. 测定堆石体密度的附加质量法理论分析研究[J]. 地球物理学报, 2015, 58(1): 257-266. https://www.cnki.com.cn/Article/CJFDTOTAL-DQWX201501023.htm ZHANG Zhi, FAN Youhua, CAI Jiaxing, et al. Theoretical analysis on density measurement of rockfill body by additive mass method[J]. Chinese Journal of Geophysics, 2015, 58(1): 257-266. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-DQWX201501023.htm
[6]	徐文杰, 胡瑞林. 土石混合体概念、分类及意义[J]. 水文地质工程地质, 2009, 36(4): 50-56, 70. https://www.cnki.com.cn/Article/CJFDTOTAL-SWDG200904014.htm XU Wenjie, HU Ruilin. Conception, classification and significations of soil-rock mixture[J]. Hydrogeology & Engineering Geology, 2009, 36(4): 50-56, 70. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SWDG200904014.htm
[7]	李丕武, 冷元宝, 袁江华. 堆石体密度测定的附加质量法[J]. 地球物理学报, 1999, 42(3): 422-427. https://www.cnki.com.cn/Article/CJFDTOTAL-DQWX201501023.htm LI Peiwu, LENG Yuanbao, YUAN Jianghua. The density measurement of rockfill object by additive mass method[J]. Chinese Journal of Geophysics, 1999, 42(3): 422-427. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-DQWX201501023.htm
[8]	SHEN C M, YU J D, LIU S H, et al. A unified fractional breakage model for granular materials inspired by the crushing tests of dyed gypsum particles[J]. Construction and Building Materials, 2021, 270: 121366.
[9]	WANG T, LIU S H, WAUTIER A, et al. Updated skeleton void ratio for gravelly sand mixtures considering effect of grain-size distribution[J]. Canadian Geotechnical Journal, 2022, 59(1): 12-23.
[10]	CHANG C S, YIN Z Y. Micromechanical modeling for behavior of silty sand with influence of fine content[J]. International Journal of Solids and Structures, 2011, 48(19): 2655-2667.
[11]	YIN Z Y, HUANG H W, HICHER P Y. Elastoplastic modeling of sand-silt mixtures[J]. Soils and Foundations, 2016, 56(3): 520-532.

施引文献

资源附件(0)

图(8)

计量

文章访问数: 269
HTML全文浏览量: 37
PDF下载量: 111
被引次数: 0

0. 引言
1. 管-土-高聚物相互作用分析模型
1.1 管-土-高聚物相互作用Winkler模型
1.2 管-土-高聚物相互作用Vlazov模型
2. 基于传递矩阵法的管道平面应变问题分析
2.1 传递矩阵法
2.2 管道计算模型描述
2.3 程序实现
3. 试验概述及数据处理
4. 模型计算及与试验结果对比
5. 不同因素对修复后管道内力和变形的影响
5.1 不同埋深的影响
5.2 不同附加荷载的影响
5.3 不同土壤弹性模量的影响
5.4 不同管壁厚度的影响
6. 结论

0. 引言
1. 管-土-高聚物相互作用分析模型
1.1 管-土-高聚物相互作用Winkler模型
1.2 管-土-高聚物相互作用Vlazov模型
2. 基于传递矩阵法的管道平面应变问题分析
2.1 传递矩阵法
2.2 管道计算模型描述
2.3 程序实现
3. 试验概述及数据处理
4. 模型计算及与试验结果对比
5. 不同因素对修复后管道内力和变形的影响
5.1 不同埋深的影响
5.2 不同附加荷载的影响
5.3 不同土壤弹性模量的影响
5.4 不同管壁厚度的影响
6. 结论

参考文献(11)

施引文献

资源附件(0)

附加质量法在土石混合料密度检测中的适用性研究 English Version

作者简介: 徐思远(1993—)，男，工程师，主要从事土石坝渗漏病害探测、岩土原位测试等方面的研究工作。E-mail: syxu@nhri.cn

通讯作者: 傅中志, E-mail: fu_zhongzhi@yahoo.com

计量

出版历程

Applicability of additive mass method in density tests on soil-mixture materials

0. 引言

1. 管-土-高聚物相互作用分析模型

1.1 管-土-高聚物相互作用Winkler模型

1.2 管-土-高聚物相互作用Vlazov模型

2. 基于传递矩阵法的管道平面应变问题分析

2.1 传递矩阵法

2.2 管道计算模型描述

2.3 程序实现

3. 试验概述及数据处理

4. 模型计算及与试验结果对比

5. 不同因素对修复后管道内力和变形的影响

5.1 不同埋深的影响

5.2 不同附加荷载的影响

5.3 不同土壤弹性模量的影响

5.4 不同管壁厚度的影响

6. 结论

计量

出版历程

目录

0. 引言

1. 管-土-高聚物相互作用分析模型

1.1 管-土-高聚物相互作用Winkler模型

1.2 管-土-高聚物相互作用Vlazov模型

2. 基于传递矩阵法的管道平面应变问题分析

2.1 传递矩阵法

2.2 管道计算模型描述

2.3 程序实现

3. 试验概述及数据处理

4. 模型计算及与试验结果对比

5. 不同因素对修复后管道内力和变形的影响

5.1 不同埋深的影响

5.2 不同附加荷载的影响

5.3 不同土壤弹性模量的影响

5.4 不同管壁厚度的影响

6. 结论

作者简介:
徐思远(1993—)，男，工程师，主要从事土石坝渗漏病害探测、岩土原位测试等方面的研究工作。E-mail: syxu@nhri.cn

通讯作者:
傅中志, E-mail: fu_zhongzhi@yahoo.com