黄土高填方场地工后沉降预测新模型

于永堂; 郑建国

doi:10.3969/j.issn.0258-2724.20200872

黄土高填方场地工后沉降预测新模型

doi: 10.3969/j.issn.0258-2724.20200872

于永堂^{1, 2, 3,},
郑建国^{1, 2}

1.
机械工业勘察设计研究院有限公司，陕西西安 710043
2.
西安建筑科技大学土木工程学院，陕西西安 710055
3.
中联西北工程设计研究院有限公司，陕西西安 710077

基金项目: 国家自然科学基金（41790442）；陕西省“三秦学者”创新团队支持计划（2013KCT-13）；陕西省技术创新引导专项(基金)计划（2020CGHJ-002）；陕西省秦创原“科学家 + 工程师”队伍建设项目（2022KXJ-086）

详细信息

作者简介:
于永堂（1983—），男，高级工程师，研究方向为岩土工程监测与测试技术、湿陷性土地基处理技术，E-mail：yuyongtang@126.com

中图分类号: TU444
计量
- 文章访问数: 359
- HTML全文浏览量: 245
- PDF下载量: 43
- 被引次数: 10
出版历程
- 收稿日期: 2021-01-12
- 修回日期: 2021-04-01
- 网络出版日期: 2022-10-09
- 刊出日期: 2021-05-26

New Prediction Model for Post-Construction Settlement of Loess High Fill Site

YU Yongtang^{1, 2, 3
,},
ZHENG Jianguo^{1, 2}

1.
China Jikan Research Institute of Engineering Investigation and Design Co., Ltd., Xi’an 710043, China
2.
College of Civil Engineering, Xi’an University of Architecture and Technology, Xi’an 710055, China
3.
China United Northwest Institute for Engineering Design and Reasearch Co., Ltd., Xi’an 710077, China

摘要

摘要:
黄土高填方场地的工后沉降预测结果是确定地面工程建设时序和空间布局的重要依据. 为准确预测黄土高填方场地的工后沉降量，在分析典型黄土高填方场地沉降数据特点、曲线特征和发展演化规律的基础上，提出了收敛型和发散型两种用于工后沉降预测的新模型；介绍了新模型的基本性质与参数求解方法，并检验了新模型在典型黄土高填方场地工后沉降预测中的应用效果. 结果表明：新模型的内拟合误差和外推预测误差均较小，适合黄土高填方场地的工后沉降预测，其中发散型模型更适合“S”形沉降曲线的外推预测，平均绝对百分比误差（MAPE）为4.6%，相较于传统预测模型的外推预测误差降低了78.7% ~ 95.8%；收敛型模型更适合“J”形沉降曲线的外推预测，平均绝对百分比误差（MAPE）为1.9%，相较于传统预测模型的外推预测误差降低了68.3% ~ 84.4%；新模型具有较好的适应性、通用性和稳定性，可为今后黄土高填方场地的工后沉降预测和评估提供更多的选择和参考.
- 黄土 /
- 高填方场地 /
- 工后沉降 /
- 预测模型
Abstract:
The post-construction settlement prediction of loess high fill site is an important basis for determining the time sequence of construction and the spatial layout of ground engineering. In order to accurately predict the post-construction settlement of loess high fill site, two post-construction settlement predictive models, one is a convergent model and the other is a divergent one, were proposed based on the characteristics of settlement data and the evolution laws of settlement curves of two typical loess high fill sites. The basic properties of the models, the solution to the parameters of the models, and the applicability of the proposed models in predicting the post-construction settlement of typical loess high fill site were all presented in detail.The results show that the proposed models are suitable for post-construction settlement prediction of loess high fill site due to its smaller fitting error and prediction error. Moreover, it is found that the divergent model is more suitable to the grounds whose settlement curves have S-shape, the mean absolute percentage error (MAPE) is 4.6%, which is 78.7%−95.8% lower than the prediction error of the traditional prediction models, while the convergence model is more applicable to the grounds whose settlement curves have J-shape, the mean absolute percentage error (MAPE) is 1.9%, which is 68.3%−84.4% lower than the prediction error of the traditional prediction models. Due to their good adaptability, generality and stability, the proposed models can provide more choices and references for the prediction and evaluation of post-construction settlement of loess high fill site in the future.
- loess /
- high fill site /
- post-construction settlement /
- prediction model

HTML全文

单箱双室箱梁在城市桥梁和公路桥梁中应用广泛. 采用该截面形式桥梁的箱体和桥面宽度较大，在竖向荷载作用有偏心时，其畸变效应不同于单箱单室箱梁. 箱梁的畸变效应研究文献较多，一般采用能量变分法或板元分析法都可得到一个四阶控制微分方程，方程未知量为畸变角或箱梁梁肋挠度w，都能很好地揭示单箱单室箱梁的畸变效应.

张元海等^[1]采用能量变分原理对单箱单室箱梁的畸变效应进行了研究，分析了双层悬臂板对畸变的影响；徐勋等^[2-3]采用广义坐标法研究了单箱单室箱梁的畸变；王兆南等^[4]采用板元分析法研究了单箱单室箱梁的畸变效应；Pezeshky等^[5]研究了宽翼缘钢梁的畸变. 单箱单室简支直线箱梁畸变的研究文献较多，成果较多. 曲线箱梁因扭转、弯曲和畸变相互耦合，畸变研究较为复杂^[6-8]. 畸变的试验研究可对畸变理论分析结果进行验证：刘保东等^[9]对连续刚构桥的扭转和畸变进行了试验研究；狄瑾等^[10]对波形钢腹板组合箱梁的扭转和畸变进行了试验研究，分析了扭转、畸变和弯曲正应力的占比. 波形钢腹板组合箱梁由于结构自重轻，适合大跨度桥梁的设计，其畸变也比混凝土箱梁明显^[11]，畸变研究也较深入. Li等^[12-14]研究了变截面箱梁的畸变，不同于求解常截面箱梁畸变控制微分方程常采用的初参解法，而采用纽马克法求解变截面箱梁的畸变控制微分方程. 双室箱梁桥的畸变不同于单室箱梁桥，蔺鹏臻等^[15]研究了单箱双室箱梁对称弯曲时的局部扭转效应；Chidolue等^[16]研究了多室箱梁的扭转和畸变，多数文献在研究双室箱梁的畸变效应时，只考虑反对称畸变，而忽略了正对称畸变.

单箱双室箱梁相比单箱单室箱梁增加了一道腹板，这使得在偏心竖向荷载作用下，箱梁畸变不仅有反对称畸变，还存在正对称畸变. 若忽略正对称畸变，在理论上将不完善，计算结果误差也较大，且双室箱梁在得出畸变翘曲惯性矩，框架横向抗弯惯性矩时也不同于单箱单室箱梁. 由于有两种畸变模式，因此无法采用一个微分方程来描述畸变效应，需采用两个未知量，对应两个微分方程才能完整描述单箱双室箱梁的畸变效应. 本文从单箱双室箱梁的畸变荷载分析出发，以箱梁的畸变角γ_d为未知量，定义正、反对称畸变对应的畸变角，建立正、反对称畸变的四阶控制微分方程，以揭示单箱双室箱梁的畸变效应，为单箱双室箱梁桥的设计计算提供参考.

1. 基本假定和畸变荷载

在薄壁箱梁畸变基本理论假定^[1]的基础上，结合单箱双室箱梁在偏心荷载作用下的位移变形和实际设计需考虑的因素等，补充矩形截面单箱双室箱梁畸变分析的假定：单箱双室箱梁的边腹板厚度相同，箱梁横向变形时角点的转角位移可忽略.

单箱双室箱梁一般都会存在一个对称轴y，箱梁截面尺寸如图1所示. 图1中：d、a₁、a₂分别为悬臂板、腹板、顶板（底板）的宽度；t₁、t₂、t₃、t₄、t_h分别为腹板AD、底板、腹板BC、顶板、腹板KF的厚度；e为荷载偏心距；P为单箱双室箱梁顶板上线性分布的偏心竖向荷载 . 采用右手坐标系，O为原点. 箱梁截面各角点分别采用A、B、C、D、K、F表示.

图 1 单箱双室箱梁截面形式

Figure 1. Cross section of single box double-cell box girder

下载: 全尺寸图片幻灯片

如图1所示，沿梁纵向取单位长框架后，P可等效成一力矩Pe和一作用于K点的竖向力P. 力矩可分解出使箱梁发生反对称畸变的荷载 $P_{}e/{a_2}$ ，竖向力P可分解出使箱梁发生正对称畸变的荷载 $P/3$ 和 $2P/3$ ，如图2所示. 畸变正应力可对应叠加，继续分解的畸变荷载如式（1）所示.

图 2 正对称和反对称畸变荷载

Figure 2. Positive-symmetrical and antisymmetrical distortion load

下载: 全尺寸图片幻灯片

$\left\{ \begin{gathered} {P_{41}} = {P_{42}} = {P_{21}} = {P_{22}} = P_{}e/(8{a_1}), \\ {P_{11}} = {P_{12}} = {P_{31}} = {P_{32}} = P_{}e/(4{a_2}), \\ {R_{41}} = {R_{42}} = {R_{21}} = {R_{22}} = P_{}{a_2}/(6{a_1}), \\ {R_{11}} = {R_{12}} = {R_{31}} = {R_{32}} = P_{}/3, \\ \end{gathered} \right.$

(1)

式中： ${P_{41}}$ 、 ${P_{42}}$ 、 ${P_{12}}$ 、 ${P_{22}}$ 、 ${P_{21}}$ 、 ${P_{11}}$ 、 ${P_{31}}$ 和 ${P_{32}}$ 为反对称畸变荷载； ${R_{41}}$ 、 ${R_{42}}$ 、 ${R_{12}}$ 、 ${R_{22}}$ 、 ${R_{21}}$ 、 ${R_{11}}$ 、 ${R_{31}}$ 和 ${R_{32}}$ 为正对称畸变荷载.

2. 单箱双室箱梁反对称畸变

采用板元分析法研究单箱双室箱梁的反对称畸变，可得出箱梁的畸变翘曲刚度、框架横向抗弯刚度和畸变荷载之间的关系. 箱梁发生反对称畸变时，σ_dA和σ_dD为箱梁角点A、D的畸变正应力，令β₁= σ_dA/σ_dD，畸变正应力在各板件上的分布如图2所示. β₁可由各板件上的畸变正应力合成的对坐标轴力矩之和为0的条件求得.

2.1 各板元面内力系分析

沿梁纵向取一微段，离散各板件，如图3所示. 图中：q_xB、q_xA、q_xK为箱梁腹板对顶板的横向约束反力；T_odz为微段上腹板对顶板的纵向约束反力；Q_o、M_o分别为顶板上产生的面内剪力和力矩；q_xC、q_xD、q_xF为腹板对底板的横向约束反力；T_udz为腹板对底板的纵向约束反力；Q_u、M_u分别为底板上产生的面内剪力和力矩；q_yA、q_yD分别为箱梁顶板和底板对左腹板的横向约束反力；Q_c、M_c分别为左腹板面内剪力和力矩.

图 3 箱梁各板件的面内受力

Figure 3. In-plane force of box girder plates

下载: 全尺寸图片幻灯片

由顶板、底板和腹板的面内力矩平衡，顶板和底板面内沿x轴，腹板面内沿y轴方向列取力系平衡方程. 考虑到各板件之间在畸变荷载作用下维持力系平衡状态，略去高阶微量后得出：

$\begin{split} & \frac{{{{\rm{d}}^2}{M_{\text{c}}}}}{{{\rm{d}}{{\textit{z}}^2}}} + \frac{{{a_1}}}{{2{a_2}}}\left(\frac{{{{\rm{d}}^2}{M_{\rm{o}}}}}{{{\rm{d}}{{\textit{z}}^2}}} + \frac{{{{\rm{d}}^2}{M_{\text{u}}}}}{{{\rm{d}}{{\textit{z}}^2}}}\right) + \frac{{{a_1}({P_{41}} + {P_{21}})}}{{{a_2}}} \;+\\ &\quad {P_{11}} - \left[\frac{{{a_1}({q_{{\rm{x}}1}} + {q_{{\rm{x}}2}})}}{{2{a_2}}} + {q_{\rm{y}}}\right] = 0 , \end{split}$

(2)

式中： $q_{\mathrm{x} 1}=q_{\mathrm{x}A} + q_{\mathrm{x}B} + q_{\mathrm{x}K}$ ， $q_{\mathrm{x}K}$ 作用在板件AB的跨中，方向和 $q_{\mathrm{x}A}$ 相同； $q_{\mathrm{x} 2}=q_{\mathrm{x}C} + q_{\mathrm{x}D} + q_{\mathrm{x}F}$ ， $q_{\mathrm{x}F}$ 作用在板件CD的跨中，方向和 $q_{\mathrm{x}C}$ 相同； $q_{\mathrm{y}}=q_{\mathrm{y}A} + q_{\mathrm{y}D}$ .

箱梁在反对称畸变时各角点的位移如图4所示，图中：Δ_HA、Δ_HB、Δ_HC、Δ_HD为对应各角点水平位移；Δ_VA、Δ_VB、Δ_VC、Δ_VD为对应各角点竖向位移；α₁、α₂、α₃为对应各板件畸变后的角度改变.

图 4 箱梁各角点的位移

Figure 4. Corner displacements of box girder

下载: 全尺寸图片幻灯片

定义单箱双室箱梁反对称畸变的畸变角γ_d1为

${\gamma _{{\rm{d1}}}} = \frac{{{\alpha _1} + {\alpha _2}}}{2} + {\alpha _3}{\text{ = }}\frac{{{\varDelta _{{\rm{V}}A}} + {\varDelta _{{\rm{V}}D}}}}{{{a_2}}} + \frac{{{\varDelta _{{\rm{H}}D}} + {\varDelta _{{\rm{H}}A}}}}{{{a_1}}}.$

(3)

各板件上的力矩M_o、M_u、M_c之间存在对应关系^[4]. 对γ_d1关于z轴求二次微分后，考虑到箱梁各角点位移和面内力矩的关系^[1]，得到

${\gamma ''_{{\rm{d1}}}} = - \frac{{4{M_{\text{c}}}}}{{{a_2}E{J_1}}},$

(4)

式中：E为材料弹性模量；J ₁为腹板 AD （腹板 BC ）的面内惯性矩.

对式（4）再求二次微分后，可将式（2）化简为

$\begin{split} & - E\frac{{a_2^2{J_1}(1 + {\beta _1}) + a_1^2({\beta _1}{J_4} + {J_2})}}{{4{a_2}(1 + {\beta _1})}}{\gamma^{(4)}_{{\rm{d}}1}} +\\ &\quad \frac{{{a_1}}}{{{a_2}}}({P_{41}} + {P_{21}}) + {P_{11}} - \left[\frac{{{a_1}({q_{{\rm{x}}1}} + {q_{{\rm{x}}2}})}}{{2{a_2}}} + {q_{\rm{y}}}\right] = 0, \end{split}$

(5)

式中：J_n=t_n $a_n^3$ /12 (n=1,2)，J₄=t₄ (a₂ + 2d)³/12，J₂、J₄分别为底板、顶板的面内惯性矩.

2.2 各板元面外力系分析

采用箱梁各板件的面外力系分析化简式（5）中的 ${a_1}({q_{{\rm{x}}1}} + {q_{{\rm{x}}2}})/(2{a_2}) + {q_{\rm{y}}}$ 项. 在分析过程中考虑箱梁腹板KF对框架横向抗弯惯性矩的贡献，单箱双室箱梁在反对称荷载作用下发生畸变后，离散箱梁各板件，各板件的面外受力如图5所示，图中变量m_AK、m_AD等为各板件板端的畸变横向弯矩.

图 5 反对称畸变时各板件的面外受力

Figure 5. Out-of-plane force with antisymmetrical distortion

下载: 全尺寸图片幻灯片

反对称畸变时，取出板件AK和DF分析面外受力可得： ${q_{{\rm{y}}A}}{a_2} = 2({m_{AK}} + {m_{KA}})$ ， ${q_{{\rm{y}}D}}{a_2} = 2({m_{DF}} + {m_{FD}})$ . 取出板件AD和BC分析面外受力可得： ${q_{{\rm{x}}A}}{a_1} = {m_{AD}} + {m_{DA}}$ ， ${q_{{\rm{x}}B}}{a_1} = {m_{BC}} + {m_{CB}}$ ， ${q_{{\rm{x}}A}} = {q_{{\rm{x}}D}}$ ， ${q_{{\rm{x}}B}} = {q_{{\rm{x}}C}}$ . 取出板件KF分析面外受力可得： ${q_{{\rm{x}}K}}{a_1} = {m_{KF}} + {m_{FK}} - {m_{KA}} - {m_{KB}} - {m_{FC}} - {m_{FD}}$ .

板件KB和FC的面外力系分析和板件AK等相同. 为得出各板件的板端畸变横向弯矩和板端位移的关系，可采用力矩分配法进行分析，见式（6）.

$\left\{ \begin{gathered} {m_{{{AK}}}} = \frac{{24E{I_4}{i_1}}}{{a_2^2({i_1} + {i_4})}}{\varDelta _{\rm{V}}} + \frac{{12E{I_1}{i_4}}}{{a_1^2({i_1} + {i_4})}}{\varDelta _{\rm{h}}}, \\ {m_{{{AD}}}} = - \frac{{12E{I_1}{i_4}}}{{a_1^2({i_1} + {i_4})}}{\varDelta _{\rm{h}}} - \frac{{24E{I_4}{i_1}}}{{a_2^2({i_1} + {i_4})}}{\varDelta _{\rm{V}}}, \\ {m_{{{DF}}}} = \frac{{24E{I_2}{i_1}}}{{a_2^2({i_1} + {i_2})}}{\varDelta _{\rm{V}}} + \frac{{12E{I_1}{i_2}}}{{a_1^2({i_1} + {i_2})}}{\varDelta _{\rm{h}}}, \\ {m_{{{DA}}}} = - \frac{{12E{I_1}{i_2}}}{{a_1^2({i_1} + {i_2})}}{\varDelta _{\rm{h}}} - \frac{{24E{I_2}{i_1}}}{{a_2^2({i_1} + {i_2})}}{\varDelta _{\rm{V}}}, \\ {m_{{{KA}}}} = \frac{{24E{I_4}{i_{\rm{h}}}}}{{a_2^2(2{i_4} + {i_{\rm{h}}})}}{\varDelta _{\rm{V}}} + \frac{{12E{I_{\rm{h}}}{i_4}}}{{a_1^2(2{i_4} + {i_{\rm{h}}})}}{\varDelta _{\rm{h}}}, \\ {m_{{{FD}}}} = \frac{{24E{I_2}{i_{\rm{h}}}}}{{a_2^2(2{i_2} + {i_{\rm{h}}})}}{\varDelta _{\rm{V}}} + \frac{{12E{I_{\rm{h}}}{i_2}}}{{a_1^2(2{i_2} + {i_{\rm{h}}})}}{\varDelta _{\rm{h}}}, \\ \end{gathered} \right.$

(6)

式中：i₁、i₂、i₄、i_h分别为箱梁边腹板、底板、顶板和中腹板的线刚度；I₁、I₂、I₄、I_h分别为箱梁边腹板、底板、顶板和中腹板的面外惯性矩；Δ_V、Δ_h分别为板件板端竖向位移和水平位移.

在 ${a_1}({q_{{\rm{x}}1}} + {q_{x2}})/(2{a_2}) + {q_{\rm{y}}}$ 中代入q_x1、 q_x2、 q_y的表达式，结合式（6），考虑到式（3）中有Δ_V = Δ_VA =Δ_VD，Δ_h = Δ_HD = Δ_HA，得到

$\begin{split} & \frac{a_{1}\left(q_{\mathrm{x} 1}+q_{\mathrm{x} 2}\right)}{2 a_{2}}+q_{\mathrm{y}}=\frac{48 E}{a_{2}}\left(\frac{I_{\mathrm{h}} I_{4}}{4 I_{4} a_{1}+I_{\mathrm{h}} a_{2}}\;+\right. \\ &\quad \left.\frac{I_{\mathrm{h}} I_{2}}{I_{\mathrm{h}} a_{2}+4 I_{2} a_{1}}+\frac{I_{1} I_{4}}{I_{1} a_{2}+2 I_{4} a_{1}}+\frac{I_{1} I_{2}}{I_{1} a_{2}+2 I_{2} a_{1}}\right) \gamma_{\mathrm{d} 1}. \end{split}$

(7)

将 ${P_{41}}$ 、 ${P_{21}}$ 、 ${P_{11}}$ 代入式（5）中的 $\dfrac{{{a_1}}}{{{a_2}}}({P_{41}} + {P_{21}}) \;+ {P_{11}}$ 项，化简得到 $P_{}e/(2{a_2})$ .

2.3 反对称畸变控制微分方程

根据以上分析可将式（5）化简为单箱双室箱梁在偏心竖向荷载作用下，以反对称畸变角γ_d1为未知量的四阶畸变控制微分方程，如式（8）.

$E{I_{{{{\text{ω}} {\rm{d}}1}}}}{\gamma ^{(4)}_{{\rm{d}}1}} + E{I_{{\rm{R}}1}}{\gamma _{{\rm{d}}1}} = P_{}e/2 ,$

(8)

式中：I_ωd1为单箱双室箱梁反对称畸变的翘曲惯性矩，单位m⁶，如式（9）；I_R1为单箱双室箱梁反对称畸变的横向抗弯惯性矩，单位m²，如式（10）.

${I_{{\text{ω}} {\rm{d}}1}}{\rm{ = }}\frac{{a_2^2(1 + {\beta _1}){J_1} + a_1^2({\beta _1}{J_4} + {J_2})}}{{4(1 + {\beta _1})}},$

(9)

$\begin{split} &I_{\mathrm{R} 1}=48\left(\frac{I_{\mathrm{h}} I_{4}}{4 I_{4} a_{1}+I_{\mathrm{h}} a_{2}}+\frac{I_{\mathrm{h}} I_{2}}{I_{\mathrm{h}} a_{2}+4 I_{2} a_{1}}\;+\right. \\ &\quad \left.\frac{I_{1} I_{4}}{I_{1} a_{2}+2 I_{4} a_{1}}+\frac{I_{1} I_{2}}{I_{1} a_{2}+2 I_{2} a_{1}}\right). \end{split}$

(10)

3. 单箱双室箱梁正对称畸变

3.1 正对称畸变时箱梁畸变翘曲惯性矩

如图2所示，单箱双室箱梁在正对称畸变荷载作用下，产生正对称畸变翘曲变形. 按假定，畸变正应力在各板件上呈线性分布，关于y轴正对称. 考虑到变形亦呈正对称，顶板和底板有竖向位移，而水平位移很小几乎可以忽略，因此可采用一半结构AKFD进行分析，如图6所示，图中：h为箱梁高度. 考虑结构AKFD各角点角位移不相等的情况，当箱梁各板件厚度各不相等时，现有单箱单室箱梁畸变微分方程将不再适用，因按照假定和分析过程，其都建立在箱梁截面至少关于一个坐标轴对称的基础上.

图 6 箱梁正对称畸变扇性坐标ω_d

Figure 6. Positive-symmetrical distortion sector coordinates ω_d

下载: 全尺寸图片幻灯片

采用一个参数β₁描述畸变正应力在各板件上分布的模式不再满足箱梁各板件厚度都不同的情形. 现采用两个参数β₁和β₂分别描述畸变正应力在腹板、顶板（底板）上的分布. 对各板件厚度都不同的箱梁结构，令β₂=σ_dA/σ_dK （σ_dK为箱梁角点 K 的畸变正应力），认为由畸变正应力引起的面内力矩分布在顶、底板上的反弯点位置相同，但不再二等分. 在求出β₁和β₂后，可确定各板件厚度都不同的箱梁的畸变翘曲位移模式.

β₁和β₂可由各板件面内畸变正应力对x、y轴形成的力矩之和为0的条件，且联立式（11）、（12）求得.

${\beta _1}{\text{ = }}\frac{{2{t_2}{a^3}(\beta _2^2 + 1 - {\beta _2}) + 3{a^2}h({t_1}\beta _2^2 + {t_{\rm{h}}})}}{{{\varGamma _1} + {\varGamma _2}}},$

(11)

${\beta _2} = \frac{{3{t_4}\beta _1^2({a^2} - {d^2}) + 3{t_2}{a^2} + 2ah{t_{\rm{h}}}(\beta _1^2 + 1 - {\beta _1})}}{{3{t_4}\beta _1^2{{(a + d)}^2} + 3{t_2}{a^2} + 2ah{t_1}(\beta _1^2 + 1 - {\beta _1})}},$

(12)

式中：a为分割后单箱的顶板或底板宽度； ${\varGamma _1} = 2{t_4}(a + d)[(\beta _2^2 - {\beta _2} + 1){a^2} + {d^2}{({\beta _2} + 1)^2}]$ ； ${\varGamma _2} = 2{t_4}(a + d)\times ad(2\beta _2^2 + {\beta _2} - 1) + 3{a^2}h({t_1}\beta _2^2 + {t_{\rm{h}}})$ .

求得β₁、β₂后，根据 ${I_{\text{ω}{\rm{ d}}}} = \displaystyle\int_S {\omega _{\rm{d}}^2{\rm{d}}S}$ 可求得结构AKFD的畸变翘曲惯性矩I_ωd2，如式（13）所示，S为板件截面积.

$\begin{split} & {I_{{\text{ω}} {\rm{d}}2}} = \frac{{{h^2}{t_4}\beta _1^2[{{({\beta _2}a +{\beta _2}d + d)}^3} + {a^3}]}}{{48\beta _2^2(1 + {\beta _2}){{(1 + {\beta _1})}^2}}} + \\ &\quad \frac{{{h^2}{t_2}{a^3}(\beta _2^3 + 1)}}{{48\beta _2^2(1 + {\beta _2}){{(1 + {\beta _1})}^2}}} + \frac{{({t_{\rm{h}}} + \beta _2^2{t_1})(\beta _1^3 + 1){a^2}{h^3}}}{{48\beta _2^2{{(1 + {\beta _1})}^3}}}. \end{split}$

(13)

3.2 正对称畸变时框架横向抗弯惯性矩

结构AKFD在正对称畸变荷载作用下发生如图7所示的变形. 在图7中：设定杆件AK和DF杆端的弯矩、转角以逆时针为正；杆件AK和DF皆沿y轴有相同的竖向位移；M_AK、M_KA分别为AK杆件杆端A、K的横向弯矩；γ_d2 = Δ_V/a，为单箱双室发生正对称畸变时的畸变角.

图 7 结构AKFD的变形和横向弯矩

Figure 7. Deformation and transverse bending moment of structure AKFD

下载: 全尺寸图片幻灯片

根据分析，可忽略杆件AD和KF水平x向位移，采用杆件的转角位移公式，可用畸变角γ_d2将角点D的畸变横向弯矩表示出来，有M_DF = I_Dγ_d2. 同理可得其他板件的板端弯矩：M_KA = I_Kγ_d2、 M_AK = I_Aγ_d2、 M_FD = I_Fγ_d2. I_D、I_K、I_A和I_F见式（14），同时可方便地确定结构的横向弯曲应变能.

$\qquad\left\{ \begin{array}{l} {I_{{D}}} = \dfrac{{[(4{\varTheta _1} + {\varTheta _3}){i_2}{i_{\rm{h}}} - 3{\varTheta _1}i_2^2]({\varTheta _2} + 3{\varTheta _1}){i_1}{i_4}} + {{[(6{\varTheta _1} - {\varTheta _2}){i_2}{i_{\rm{h}}} - 3{\varTheta _1}i_2^2][3{\varTheta _1}i_1^2 - (4{\varTheta _1} + {\varTheta _3}){i_1}{i_4}]}}}{{{\varGamma _3}}} ， \\ {I_{{K}}} = \dfrac{{[(4{\varTheta _3} + {\varTheta _1}){i_1}{i_4} - 3{\varTheta _3}i_4^2]({\varTheta _2} + 3{\varTheta _3}){i_2}{i_{\rm{h}}}} + {{[(6{\varTheta _3} - {\varTheta _2}){i_1}{i_4} - 3{\varTheta _3}i_4^2][3{\varTheta _3}i_{\rm{h}}^2 - (4{\varTheta _3} + {\varTheta _1}){i_2}{i_{\rm{h}}}]}}}{{{\varGamma _4}}} ， \\ {I_{{A}}} = \dfrac{{2({\varTheta _2} + 3{\varTheta _1}){i_1}{i_4}}}{{{\varTheta _3}{i_4} + {\varTheta _1}{i_1}}}， \\ {I_{{F}}} = \dfrac{{2({\varTheta _2} + 3{\varTheta _3}){i_2}{i_{\rm{h}}}}}{{{\varTheta _1}{i_2} + {\varTheta _3}{i_{\rm{h}}}}}， \end{array} \right.$

(14)

式中：参数Θ₁、Θ₂、Θ₃与箱梁各板件线刚度有关； ${\varGamma _3}$ = $[3{\varTheta _1}i_1^2 - (4{\varTheta _1} + {\varTheta _3}){i_1}{i_4}]{\varTheta _1}({i_{\rm{h}}} - {i_2}) +[(4{\varTheta _1} + {\varTheta _3}){i_2}{i_{\rm{h}}} - 3{\varTheta _1}i_2^2] {\varTheta _1} ({i_1} - {i_4})$ ； ${\varGamma _4}$ = $[3{\varTheta _3}i_{\rm{h}}^2 - (4{\varTheta _3} + {\varTheta _1}){i_2}{i_{\rm{h}}}]{\varTheta _3}({i_1} - {i_4}) + [(4{\varTheta _3} + {\varTheta _1}){i_1}{i_4} - 3{\varTheta _3}i_4^2]{\varTheta _3}({i_{\rm{h}}} - {i_2})$ .

3.3 正对称畸变控制微分方程

箱梁的总横向弯曲应变能可表示为 ${u_{\rm{e}}} = \dfrac{1}{2}\displaystyle\int_L {\displaystyle\int_s {\dfrac{{M^2{{(s)}}}}{{EI}}{\rm{d}}s} {\rm{d}}L}$ ，其中：M(s)为箱室各板件周向s处的横向弯矩；I为各板件的面外惯性矩；L为梁长.

沿梁轴z向取单位长框架，可得箱梁各板件的横向弯曲应变能之和，如式（15）所示. 也可写为 $\dfrac{1}{2}E{I_{{\rm{R}}2}}\gamma _{{\rm{d}}2}^2$ ，I_R2为箱室各板件厚度均不相等时的框架横向抗弯惯性矩，如式（16）.

$\begin{split} &u_{{{\rm{e}}}}= \frac{\gamma_{{{\rm{d}}} 2}^{2}}{6 E}\left[\frac{a\left(I_{{A}}^{2}+I_{{K}}^{2}-I_{{A}} I_{{K}}\right)}{I_{4}}+\frac{h\left(I_{{K}}^{2}+I_{{F}}^{2}-I_{{K}} I_{{F}}\right)}{I_{{{\rm{h}}}}}+\right.\\ &\quad \left.\frac{a\left(I_{{{{D}}}}^{2}+I_{{F}}^{2}-I_{{{{D}}}} I_{{F}}\right)}{I_2}+\frac{h\left(I_{{A}}^{2}+I_{{D}}^{2}-I_{{A}} I_{{{{D}}}}\right)}{I_1}\right],\\[-16pt] \end{split}$

(15)

$\begin{split} & I_{\mathrm{R} 2}= \frac{1}{3 E^{2}}\left[\frac{a\left(I_{{A}}^{2}+I_{{K}}^{2}-I_{{A}} I_{{K}}\right)}{I_{4}}+\frac{h\left(I_{{K}}^{2}+I_{{F}}^{2}-I_{{K}} I_{{F}}\right)}{I_{\rm{h}}}+\right.\\ &\quad \left.\frac{a\left(I_{{D}}^{2}+I_{{F}}^{2}-I_{{D}} I_{{F}}\right)}{I_{2}}+\frac{h\left(I_{{A}}^{2}+I_{{D}}^{2}-I_{{A}} I_{{D}}\right)}{I_{1}}\right].\\[-16pt] \end{split}$

(16)

通过以上分析可得，单箱双室箱梁发生正对称畸变时，取一半结构得出的畸变微分方程如式（17）. I_R2单位为m²， ${I_{{\text{ω}}{\rm{d}}2}}$ 单位为m⁶.

$E{I_{{\text{ω}} {\rm{d}}2}}{\gamma ^{(4)}_{{\rm{d}}2}} + E{I_{{\rm{R}}2}}{\gamma _{{\rm{d}}2}} = Pa/3.$

(17)

采用畸变角γ_d1、γ_d2分别描述箱梁的反对称和正对称畸变，得出两个畸变微分方程，方程可采用初参数法或弹性地基梁比拟法进行求解^{[1, 4]}.

4. 数值算例及参数分析

在等高度矩形截面简支箱梁上设置端横隔板，计算跨径L=80 m，截面如图8所示. E=35 GPa，泊松比μ=0. 在箱梁顶板上作用偏心竖向分布力矩1 kN · m/m，箱梁截面关于y轴对称. 采用本文单箱双室箱梁畸变理论进行计算分析，同时去掉中腹板按照单箱单室箱梁计算对比，得到的各畸变数值见表1，表中：B_d为畸变双力矩；ω_dA、ω_dD为角点A、D的畸变扇性坐标. B_d、γ_d的值为弹性地基梁比拟法中单位荷载作用下的计算值，σ_d为实际荷载作用下的值.

图 8 箱梁截面尺寸（单位：cm）

Figure 8. Box girder cross section （unit： cm）

下载: 全尺寸图片幻灯片

表 1 箱梁反对称畸变计算值

Table 1. Antisymmetrical distortional values

方法	I_ωd1/m⁶	I_R1/m²	ω_dD/m²	ω_dA/m²	β₁	γ_d1/（ × 10⁻¹⁰ rad）	B_d/（N·m²）	σ_dA/kPa	σ_dD/kPa
单箱单室箱梁畸变	12.1728	0.0040	2.5426	1.2075	0.4749	3.3961	2.6279	0.5213	1.0978
单箱双室箱梁反对称畸变	12.1728	0.0096	2.5426	1.2075	0.4749	1.7660	2.1110	0.2094	0.4409

下载: 导出CSV

| 显示表格

单箱单室箱梁增加一道腹板成为单箱双室箱梁后，σ_dA由0.5213 kPa变为0.2094 kPa，应力减小率为：（0.5213−0.2094）/ 0.5213 × 100%=59.83%，σ_dD的减小率也为59.83%. 因单室和双室箱梁反对称畸变的假定相同，角点应力比β₁数值不变，各角点的畸变扇性坐标数值不变. 由于中腹板增强了截面的横向抗弯刚度，箱梁的框架横向抗弯惯性矩从单室的0.0040增大到双室的0.0096.

算例箱梁正对称畸变计算值如表2所列. 单箱双室箱梁角点A的正对称畸变正应力为反对称畸变值的28.08%，角点D为26.76%. 若取的一半结构按照传统的单箱单室箱梁畸变方法计算角点畸变正应力，角点D的值为单箱双室箱梁反对称畸变值的1.97倍，为按照本文正对称畸变计算值的7.37倍.

表 2 箱梁正对称畸变计算值

Table 2. Positive-symmetrical distortional values

方法	角点应力比	γ_d/（ × 10⁻¹⁰ rad）	σ_dA/kPa	σ_dD/kPa
单箱双室箱梁正对称畸变	β₁ = 0.5000， β₂ = 0.6241	0.0160	0.05897	0.11795
单箱单室箱梁畸变	0.3229	3.7516	0.28055	0.86898

下载: 导出CSV

| 显示表格

最终单箱双室箱梁角点A、D的畸变正应力为箱梁正、反对称畸变计算结果的叠加，如表3. 并与有限元结果进行对比，有限元计算采用Shell-63单元建模分析. 从角点A、D的畸变正应力误差来看，本文方法解析解和有限元解最大误差的绝对值不超过8.71%，角点A的畸变正应力误差仅为1.73%，有限元解和解析解吻合良好.

表 3 箱梁角点畸变正应力

Table 3. Corner distortion normal stress of box girder

畸变正应力	单箱双室箱梁畸变解析解/kPa			有限元解 ②/kPa	误差（（①−②）/② × 100） /%
畸变正应力	反对称畸变	正对称畸变	最终结果 ①	有限元解 ②/kPa	误差（（①−②）/② × 100） /%
σ_dA	0.20940	0.05897	0.26837	0.26380	1.73
σ_dD	0.44090	0.11795	0.55885	0.61220	−8.71

下载: 导出CSV

| 显示表格

本文方法得到的双室箱梁和按单箱单室箱梁计算的角点A、D的应力值如表3和表1所示. 考虑箱梁的正对称畸变后，箱梁由单箱单室变成单箱双室时，角点A畸变正应力的减小率由59.83%变为48.52%，角点D畸变正应力的减小率由59.83%变为49.09%，如忽略正对称畸变的影响，畸变效应和有限元结果的误差将会增大. 双室箱梁的畸变正应力比单室箱梁的要小. 当箱梁横向宽度较大时，增设中腹板可有效减小箱梁的畸变效应.

为更好地比较两种截面形式箱梁的畸变效应差别. 单箱单室和单箱双室简支箱梁跨中截面的畸变角、角点畸变正应力的对比如图9、10所示. 由图9可知：单箱双室箱梁的反对称畸变角沿梁长的分布在畸变角数值的绝对值上小于单箱单室箱梁；中腹板对箱梁畸变的削弱非常明显，跨中截面畸变角迅速变小；在梁端一定区域内，双室箱梁的畸变角变化比单室的平稳，双室箱梁在此区域的反向变形小于单室箱梁，畸变角较早地衰减趋向于0，偏心竖向荷载引起的畸变较单箱单室箱梁的小，较好地反映了增加的箱室腹板对畸变的限制作用，与增加梁跨内横隔板效果相同.

图 9 畸变角沿梁长变化

Figure 9. Distortional angle changes of box girder

下载: 全尺寸图片幻灯片

图 10 角点A畸变正应力沿梁长变化

Figure 10. Changes of distortional normal stress of point A

下载: 全尺寸图片幻灯片

由图10可知：对箱梁畸变正应力沿梁长的分布进行分析，单室箱梁和双室箱梁角点A的畸变正应力沿梁长的分布有着明显的区别；双室箱梁角点A畸变正应力在数值的绝对值上明显小于单室箱梁，在跨中截面有最大值；双室箱梁角点A的畸变正应力向梁端衰减很快，而单室箱梁则衰减较慢. 从以上分析可以看出，双室箱梁畸变角、畸变正应力计算值都小于单室箱梁，中腹板对箱梁畸变的限制是非常明显的.

取中腹板厚度为0.10、0.15、0.20、0.25、0.30、0.35、0.40、0.45、0.50 m，其他条件皆不变，分析箱梁跨中截面反对称畸变角随中腹板厚度变化的规律. 中腹板厚度变化时，双室箱梁反对称畸变的畸变角沿梁长的变化情况如图11所示. 选取中腹板厚度为0.10、0.30、0.50 m的数据列出，0.10 m厚的板为薄板，0.30 m厚的板为薄板和厚板的界限，0.50 m厚的板为厚板.

图 11 中腹板厚度对箱梁畸变角的影响

Figure 11. Influence of mid-web thickness on distortional angle of box girder

下载: 全尺寸图片幻灯片

从图11可以看出：不同中腹板厚度的双室箱梁反对称畸变角沿梁长的变化趋势在箱梁跨中截面区域内基本相同，在箱梁梁端附近变化趋势有所不同；随着中腹板厚度的增加，在梁端附近，畸变角的变化逐渐平稳，双室箱梁在梁端附近畸变翘曲的程度较小；中腹板厚度为0.10 m的双室箱梁畸变角在梁端附近的反向变形比0.50 m板厚的要突出，薄壁箱梁的畸变翘曲程度比厚壁的要明显；随着板厚的增加，双室箱梁的畸变逐渐减小，仅在跨中一定长度范围内畸变较为明显，梁两端逐渐减弱，梁端变形程度较小的区域逐渐向跨内延伸；0.50 m板厚的双室箱梁梁端畸变较小的区域较0.10 m板厚的区域要长.

5. 结　论

1）在偏心竖向荷载作用下，单箱双室箱梁不仅存在反对称畸变，还存在正对称畸变. 正对称畸变正应力较小，本文算例中，最大只占到反对称畸变正应力的28.08%.

2）考虑了正对称畸变效应后，双室箱梁的畸变分析在理论上将更加完善. 通过数值计算显示，畸变正应力解析解和有限元解的误差减小，误差绝对值最大为8.71%，相互吻合更好.

3）当箱体较宽的单箱单室箱梁增设中间腹板成为单箱双室箱梁后，箱梁角点处的畸变正应力可减少到单室箱梁的49.09%. 为减小宽箱梁在偏心竖向荷载作用下的畸变，设置中腹板效果明显.

4）单箱双室简支箱梁的畸变比单箱单室简支箱梁的要小，变形主要在跨中一定的区域内，梁端附近的区域反向变形不如单箱单室的明显. 中腹板厚度的变化可使双室箱梁的畸变发生较明显的改变，厚度较小时畸变较为明显，厚度增大时畸变逐渐减弱.

5）采用两个参数β₁、β₂描述正对称畸变时，单箱双室箱梁各板件上的畸变正应力分布比一个参数的合理. 采用一个参数的箱梁畸变理论计算的正对称畸变效应，箱梁角点畸变正应力最大可为本文方法的7.37倍.

图 1 典型黄土高填方场地的工后沉降曲线

Figure 1. Duration curves of post-construction settlement of typical loess high fill site

下载: 全尺寸图片幻灯片

图 2 新模型Ⅰ的参数变化对曲线形态的影响

Figure 2. Influences of parameter change of the new model I on the shape of curves

下载: 全尺寸图片幻灯片

图 3 新模型Ⅰ的沉降速率和沉降加速度全程变化曲线

Figure 3. Settlement velocity and acceleration curves of the new model I

下载: 全尺寸图片幻灯片

图 4 新模型Ⅱ的参数变化对曲线形态的影响

Figure 4. Influences of parameter change of the new model Ⅱ on the shape of curves

下载: 全尺寸图片幻灯片

图 5 新模型Ⅱ的沉降速率和沉降加速度全程变化曲线

Figure 5. Settlement velocity and acceleration curves of the new model Ⅱ

下载: 全尺寸图片幻灯片

图 6 采用新模型进行拟合及预测的沉降曲线（工程Ⅰ）

Figure 6. Settlement curves fitted and predicted by the new models （ project Ⅰ）

下载: 全尺寸图片幻灯片

图 7 采用新模型进行拟合及预测的沉降曲线（工程Ⅱ）

Figure 7. Settlement curves fitted and predicted by the new models （ project Ⅱ）

下载: 全尺寸图片幻灯片

图 8 新模型与传统模型的工后沉降拟合与预测曲线

Figure 8. Fitting and prediction curves of post-construction settlement between the new models and some other conventional models

下载: 全尺寸图片幻灯片

表 1 新模型的回归模型参数及预测效果评价结果

Table 1. Regression parameters and evaluation results of the new models

模型类型	监测点	模型参数				拟合精度指标	预测精度指标
模型类型	监测点	a	b	c	d	R²	M_APE/%	M_FE
新模型Ⅰ	J_C6	2.296	0.410	8.121×10⁻²	0.376	1.000	4.2	−15.1
	T₁₃	0.555	0.440	5.409×10⁻²	0.278	1.000	1.2	−2.8
	S₁	0.025	0.486	1.086×10⁻²	0.381	0.998	3.5	−0.4
	S₂	0.024	0.458	9.726×10⁻³	0.376	0.998	4.6	−1.8
新模型Ⅱ	J_C6	1.219	0.830	0.202	2.407×10⁻³	0.999	1.9	6.6
	T₁₃	1.256	0.846	1.058	2.975×10⁻³	0.999	0.7	1.5
	S₁	0.023	1.592	1.340	0.538×10⁻³	0.998	11.2	5.0
	S₂	0.013	1.634	1.012	0.324×10⁻³	0.998	6.2	2.3

下载: 导出CSV

参考文献(19)

[1]	宰金珉,梅国雄. 全过程的沉降量预测方法研究[J]. 岩土力学,2000,21(4): 322-325. doi: 10.3969/j.issn.1000-7598.2000.04.003 ZAI Jinmin, MEI Guoxiong. Forecast method of settlement during the complete process of construction and operation[J]. Rock and Soil Mechanics, 2000, 21(4): 322-325. doi: 10.3969/j.issn.1000-7598.2000.04.003
[2]	潘林有,谢新宇. 用曲线拟合的方法预测软土地基沉降[J]. 岩土力学,2004,25(7): 1053-1058. doi: 10.3969/j.issn.1000-7598.2004.07.010 PAN Linyou, XIE Xinyu. Observational settlement prediction by curve fitting methods[J]. Rock and Soil Mechanics, 2004, 25(7): 1053-1058. doi: 10.3969/j.issn.1000-7598.2004.07.010
[3]	刘射洪,袁聚云,赵昕. 地基沉降预测模型研究综述[J]. 工业建筑,2014,44(增1): 738-741,681. doi: 10.13204/j.gyjz2014.s1.008 LIU Shehong, YUAN Juyun, ZHAO Xin. Review of settlement prediction models of foundation[J]. Industrial Construction, 2014, 44(S1): 738-741,681. doi: 10.13204/j.gyjz2014.s1.008
[4]	杨涛,李国维,杨伟清. 基于双曲线法的分级填筑路堤沉降预测[J]. 岩土力学,2004,25(10): 1551-1554. doi: 10.3969/j.issn.1000-7598.2004.10.008 YANG Tao, LI Guowei, YANG Weiqing. Settlement prediction of stage constructed embankment on soft ground based on the hyperbolic method[J]. Rock and Soil Mechanics, 2004, 25(10): 1551-1554. doi: 10.3969/j.issn.1000-7598.2004.10.008
[5]	陈善雄,王星运,许锡昌,等. 路基沉降预测的三点修正指数曲线法[J]. 岩土力学,2011,32(11): 3355-3360. doi: 10.3969/j.issn.1000-7598.2011.11.025 CHEN Shanxiong, WANG Xingyun, XU Xichang, et al. Three-point modified exponential curve method for predicting subgrade settlements[J]. Rock and Soil Mechanics, 2011, 32(11): 3355-3360. doi: 10.3969/j.issn.1000-7598.2011.11.025
[6]	刘宏,李攀峰,张倬元. 九寨黄龙机场高填方地基工后沉降预测[J]. 岩土工程学报,2005,27(1): 90-93. doi: 10.3321/j.issn:1000-4548.2005.01.015 LIU Hong, LI Panfeng, ZHANG Zhuoyuan. Prediction of the post-construction settlement of the high embankment of Jiuzhai—Huanglong Airport[J]. Chinese Journal of Geotechnical Engineering, 2005, 27(1): 90-93. doi: 10.3321/j.issn:1000-4548.2005.01.015
[7]	王海英,常肖,阮祺,等. 建筑垃圾填埋路基沉降预测的三点-星野法[J]. 铁道科学与工程学报,2017,14(3): 473-479. doi: 10.3969/j.issn.1672-7029.2017.03.006 WANG Haiying, CHANG Xiao, RUAN Qi, et al. Subsidence prediction of subgrade filled by construction waste based on three point-hoshino algorithm[J]. Journal of Railway Science and Engineering, 2017, 14(3): 473-479. doi: 10.3969/j.issn.1672-7029.2017.03.006
[8]	王伟,卢廷浩. 基于Weibull曲线的软基沉降预测模型分析[J]. 岩土力学,2007,28(4): 803-806,811. doi: 10.3969/j.issn.1000-7598.2007.04.033 WANG Wei, LU Tinghao. Study on prediction model of soft foundation settlement based on Weibull curve[J]. Rock and Soil Mechanics, 2007, 28(4): 803-806,811. doi: 10.3969/j.issn.1000-7598.2007.04.033
[9]	徐洪钟,施斌,李雪红. 全过程沉降量预测的Logistic生长模型及其适用性研究[J]. 岩土力学,2005,26(3): 387-391. doi: 10.3969/j.issn.1000-7598.2005.03.010 XU Hongzhong, SHI Bin, LI Xuehong. Logistic growth model and its applicability for predicting settlement during the whole process[J]. Rock and Soil Mechanics, 2005, 26(3): 387-391. doi: 10.3969/j.issn.1000-7598.2005.03.010
[10]	余闯,刘松玉. 路堤沉降预测的Gompertz模型应用研究[J]. 岩土力学,2005,26(1): 82-86. doi: 10.3969/j.issn.1000-7598.2005.01.017 YU Chuang, LIU Songyu. A Study on prediction of embankment settlement with the gompertz model[J]. Rock and Soil Mechanics, 2005, 26(1): 82-86. doi: 10.3969/j.issn.1000-7598.2005.01.017
[11]	邓英尔,谢和平. 全过程沉降预测的新模型与方法[J]. 岩土力学,2005,26(1): 1-4. doi: 10.3969/j.issn.1000-7598.2005.01.001 DENG Yinger, XIE Heping. New model and method of forecasting settlement during complete process of construction and operation[J]. Rock and Soil Mechanics, 2005, 26(1): 1-4. doi: 10.3969/j.issn.1000-7598.2005.01.001
[12]	罗战友,龚晓南,杨晓军. 全过程沉降量的灰色verhulst预测方法[J]. 水利学报,2003,34(3): 29-32,36. doi: 10.3321/j.issn:0559-9350.2003.03.006 LUO Zhanyou, GONG Xiaonan, YANG Xiaojun. Grey Verhulst prediction of settlement during the whole process of construction and operation[J]. Journal of Hydraulic Engineering, 2003, 34(3): 29-32,36. doi: 10.3321/j.issn:0559-9350.2003.03.006
[13]	钟才根,丁文其,王茂和,等. 神经网络模型在高速公路软基沉降预测中的应用[J]. 中国公路学报,2003,16(2): 31-34. doi: 10.3321/j.issn:1001-7372.2003.02.008 ZHONG Caigen, DING Wenqi, WANG Maohe, et al. Application of artificial neural network in settlement prediction of highway soft foundation[J]. China Journal of Highway and Transport, 2003, 16(2): 31-34. doi: 10.3321/j.issn:1001-7372.2003.02.008
[14]	葛苗苗,李宁,张炜,等. 黄土高填方沉降规律分析及工后沉降反演预测[J]. 岩石力学与工程学报,2017,36(3): 745-753. doi: 10.13722/j.cnki.jrme.2016.0014 GE Miaomiao, LI Ning, ZHANG Wei, et al. Settlement behavior and inverse prediction of post-construction settlement of high filled loess embankment[J]. Chinese Journal of Rock Mechanics and Engineering, 2017, 36(3): 745-753. doi: 10.13722/j.cnki.jrme.2016.0014
[15]	魏桦, 许文学, 周立新, 等. 某迁建工程试验段变形与沉降监测报告[R]. 北京: 空军工程设计研究局, 2012.
[16]	杜伟飞,郑建国,刘争宏,等. 黄土高填方地基沉降规律及排气条件影响[J]. 岩土力学,2019,40(1): 325-331. doi: 10.16285/j.rsm.2017.1079 DU Weifei, ZHENG Jianguo, LIU Zhenghong, et al. Settlement behavior of high loess-filled foundation and impact from exhaust conditions[J]. Rock and Soil Mechanics, 2019, 40(1): 325-331. doi: 10.16285/j.rsm.2017.1079
[17]	赵明华,刘煜,曹文贵. 软土路基沉降变权重组合S型曲线预测方法研究[J]. 岩土力学,2005,26(9): 1443-1447. doi: 10.3969/j.issn.1000-7598.2005.09.019 ZHAO Minghua, LIU Yu, CAO Wengui. Study on variable-weight combination forecasting method of S-type curves for soft clay embankment settlement[J]. Rock and Soil Mechanics, 2005, 26(9): 1443-1447. doi: 10.3969/j.issn.1000-7598.2005.09.019
[18]	张鸿燕,耿征. Levenberg-Marquardt算法的一种新解释[J]. 计算机工程与应用,2009,45(19): 5-8. doi: 10.3778/j.issn.1002-8331.2009.19.002 ZHANG Hongyan, GENG Zheng. Novel interpretation for Levenberg-Marquardt algorithm[J]. Computer Engineering and Applications, 2009, 45(19): 5-8. doi: 10.3778/j.issn.1002-8331.2009.19.002
[19]	LEWIS C. Industrial and business forecasting methods: a practical guide to exponential smoothing and curve fitting[M]. London: Butterworths, 1982.

施引文献

期刊类型引用(5)

1.	朱彦鹏，董怀堂，韦杰强，赵骏，王怡轩. 改进的分层总和法在黄土地区高填方地基沉降变形中的应用研究. 建筑科学与工程学报. 2024(03): 159-168 . 百度学术
2.	于永堂，孙茉，曹静远，朱建民，张龙. 黄土高填方场地水分迁移规律. 长江科学院院报. 2024(07): 103-109 . 百度学术
3.	孙茉，于永堂，马勃，韩文斌，曹静远. 深厚土层内部沉降监测新装置的研制. 地下空间与工程学报. 2024(S1): 350-357 . 百度学术
4.	Xiaokun Hou，Shengwen Qi，Yongtang Yu，Jianguo Zheng. Long-Term Settlement Characterization of High-Filling Foundation in the Mountain Excavation and City Construction Area of the Yan ' an New District, China. Journal of Earth Science. 2023(06): 1908-1915 . 必应学术
5.	马云峰，彭森，葛苗苗，朱才辉. 黄土高填方地基上部建筑施建控制标准研究. 防灾减灾工程学报. 2023(06): 1415-1424 . 百度学术