CN111811468B - 高心墙堆石坝沉降变形的预测方法 - Google Patents

高心墙堆石坝沉降变形的预测方法 Download PDF

Info

Publication number
CN111811468B
CN111811468B CN202010922824.XA CN202010922824A CN111811468B CN 111811468 B CN111811468 B CN 111811468B CN 202010922824 A CN202010922824 A CN 202010922824A CN 111811468 B CN111811468 B CN 111811468B
Authority
CN
China
Prior art keywords
deformation
dam
dam body
high core
fill
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.)
Active
Application number
CN202010922824.XA
Other languages
English (en)
Other versions
CN111811468A (zh
Inventor
裴亮
陈建康
周正军
李艳玲
吴震宇
张瀚
陈辰
周靖人
李俊儒
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Sichuan University
Original Assignee
Sichuan University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Sichuan University filed Critical Sichuan University
Priority to CN202010922824.XA priority Critical patent/CN111811468B/zh
Publication of CN111811468A publication Critical patent/CN111811468A/zh
Application granted granted Critical
Publication of CN111811468B publication Critical patent/CN111811468B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C5/00—Measuring height; Measuring distances transverse to line of sight; Levelling between separated points; Surveyors' levels
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B21/00—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
    • G01B21/32—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring the deformation in a solid

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种高心墙堆石坝沉降变形的预测方法，其包括：以高心墙堆石坝的建基面作为参考面，依据坝体内已布设的沉降变形监测点及其位置确定各监测点到参考面的距离z i ；判断拟进行变形预测的高心墙堆石坝处于施工期还是运行期；当高心墙堆石坝处于施工期时，依据沉降变形物理成因、组成以及距离z i ，分别计算坝体沉降变形、牵连变形和蠕变变形，并累加得到施工期坝体沉降变形；当高心墙堆石坝处于运行期时，根据监测点与浸润线相对位置和距离z i ，采用对应的总沉降变形计算方法进行总沉降变形的计算。本发明所述沉降变形预测方法更加符合工程实际，能够为高心墙堆石坝的变形调控和安全管控提供更加准确的数据支撑。

Description

高心墙堆石坝沉降变形的预测方法
技术领域
本发明涉及大坝安全监测技术领域，具体涉及一种高心墙堆石坝沉降变形的预测方法。
背景技术
高心墙堆石坝由于具有对地基条件适应性好、能就地取材和充分利用建筑物开挖料、投资省、施工方法简单以及抗震性能好等优点，在坝工建设中被广泛采用。目前，在我国西部水电开发中，将有一批高心墙堆石坝工程投入建设和运行，如长河坝（坝高240m）、两河口（坝高295m）、古水（坝高305m）、双江口（坝高314m）等，这些高心墙堆石坝变形机理复杂，变形预测难度大，因不协调变形产生的变形破坏风险高，因此，有必要从沉降变形机理与物理成因入手，研究沉降变形的高精度预测方法，为高堆石坝变形调控和安全管控提供技术支撑。
目前心墙堆石坝变形预测一般采用传统的数学模型法，这种方法是通过建立描述效应量与因变量之间的统计学模型，确定它们之间的定量关系。传统的统计回归模型进行测点变形预测，仅仅是采用多元线性回归模型进行系数拟合，未考虑测点变形的物理成因，极易出现拟合精度低、鲁棒性差、适用性低等情况，而且对于时效分量的求解，按照上述形式构建回归方程很有可能造成回归模型起始点偏差较大的现象。
发明内容
针对现有技术中的上述不足，本发明提供的高心墙堆石坝沉降变形的预测方法解决了现有技术未考虑变形监测点变形的物理成因致使预测精度差的问题。
为了达到上述发明目的，本发明采用的技术方案为：
提供一种高心墙堆石坝沉降变形的预测方法，其包括：
S1、选取高心墙堆石坝的建基面作为参考面，根据坝体内已布设的沉降变形监测点及其位置确定各监测点到参考面的距离；
S2、判断高心墙堆石坝为施工期或运行期，若为施工期，则进入步骤S3，若为运行期，则进入步骤S7；
S3、根据坝体材料的邓肯-张模型参数及监测点到参考面的距离，计算施工期坝体压缩产生的沉降变形S di ：
Figure DEST_PATH_IMAGE001
其中，H i 为第i个变形监测点到坝表面的垂直距离，单位为m；z i 为第i个变形监测点到参考面的距离，单位为m；n为坝体材料的邓肯-张模型参数；B为坝体变形系数；
S4、根据坝基覆盖层深度、坝基材料的邓肯-张模型参数及变形监测点到参考面的距离，计算施工期坝体压缩产生的牵连变形S li ：
Figure 849998DEST_PATH_IMAGE002
其中，H l 为坝基覆盖层深度，单位为m；n l 为坝基材料的邓肯-张模型参数；B l 为坝基变形系数；
S5、根据材料的蠕变参数，计算施工期的蠕变变形
Figure DEST_PATH_IMAGE003
：
Figure 581807DEST_PATH_IMAGE004
式中，K p 为蠕变变形系数；α为坝体堆石料蠕变参数；t j 为坝体开始填筑至变形监测点获得初始测值的时间，单位为天；t c 为坝体开始填筑至进行变形预测的时间，单位为天；t p 为0至t j 中间过程中选定的时间节点，单位为天；p为时间节点变量，1≤p≤j；j为选取的时间节点个数，视精度和预测时间长短而定；e为自然数；
S6、累加沉降变形S di 、牵连变形S li 和蠕变变形
Figure 158019DEST_PATH_IMAGE003
得到施工填筑期坝体任意监测点处的总沉降变形
Figure DEST_PATH_IMAGE005
，S z0 为常数；
S7、判断变形监测点i到参考面的距离是否小于变形监测点i处对应的浸润线高程H wi ，若是进入步骤S8，否则进入步骤S9；
S8、根据变形监测点i到参考面的距离和变形预测点i到坝体表面的垂直距离，计算运行期坝体任意监测点的总沉降变形S zi ：
Figure 266265DEST_PATH_IMAGE006
其中，
Figure DEST_PATH_IMAGE007
、
Figure 713165DEST_PATH_IMAGE008
均为与材料湿化变形系数；n w 为材料湿化变形参数；
S9、根据变形预测点i处对应的浸润线高程和变形预测点i到坝体表面的垂直距离，计算运行期坝体任意监测点处的总沉降变形S zi ：
Figure DEST_PATH_IMAGE009
。
本发明的有益效果为：本方案在进行沉降变形预测时，首先获取预测点的相对位置关系，并基于相对位置关系可以得到监测点所在位置的应力应变状态，最终通过应力应变关系得到每个预测点的沉降变形；由于本方案在预测时充分考虑高心墙堆石坝上述的物理成因，相对现有技术中数学模拟模型而言，综合考虑了预测点的位置关系和应力应变状态，使得本方案的沉降变形预测模型能较好地模拟和反映测点变形趋势，显著提升变形预测精度。
附图说明
图1为高心墙堆石坝沉降变形的预测方法的流程图。
图2为坝内任意点与坝高、覆盖层深度的位置关系图。
图3为坝体任意点沉降变形受坝体压缩和坝基压缩的影响。
图4为5-3FLAC3d数值计算算例中沉降变形随坝高的变化。
图5为坝填筑高度与时间的对应关系。
图6为施工期VE1-18测点沉降测值和计算值。
图7为施工期VE2-37测点沉降测值和计算值。
图8为运行期VE1-25 实测值与计算值对比。
图9为LD59 实测值与计算值对比。
具体实施方式
下面对本发明的具体实施方式进行描述，以便于本技术领域的技术人员理解本发明，但应该清楚，本发明不限于具体实施方式的范围，对本技术领域的普通技术人员来讲，只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内，这些变化是显而易见的，一切利用本发明构思的发明创造均在保护之列。
参考图1，图1示出了高心墙堆石坝沉降变形的预测方法的流程图；如图1所示，该方法S包括S1至S9。
在步骤S1中，选取高心墙堆石坝的建基面作为参考面，根据坝体内已布设的沉降变形监测点及其位置确定各监测点到参考面的距离；高心墙堆石坝的结构示意参考图2。
在步骤S2中，判断高心墙堆石坝为施工期或运行期，若为施工期，则进入步骤S3，若为运行期，则进入步骤S7；
在步骤S3中，根据坝体材料的邓肯-张模型参数及监测点到参考面的距离，计算施工期坝体压缩产生的沉降变形S di ：
Figure 108153DEST_PATH_IMAGE010
其中，H i 为第i个变形监测点到坝表面的垂直距离，单位为m；z i 为第i个变形监测点到参考面的距离，单位为m；n为坝体材料的邓肯-张模型参数；B为坝体变形系数；
在本发明的一个实施例中，所述施工期坝体压缩产生的变形S di 的构建方法包括：
S31、构建坝体内任意一点的压缩变形的初始模型：
Figure DEST_PATH_IMAGE011
（1-1）
式中，α d 为坝体压应力修正系数；γ为填筑料的容重；E s 为填筑料的压缩模量；z 1，z 2均为积分变量；
S32、采用邓肯-张模型进行坝体变形计算填筑料的压缩模量E s ：
Figure 139038DEST_PATH_IMAGE012
（1-2）
其中，
Figure DEST_PATH_IMAGE013
分别为坝体内的大、小主应力；
S33、坝体内的大、小主应力
Figure 752773DEST_PATH_IMAGE014
的计算公式分别为：
Figure DEST_PATH_IMAGE015
S34、将公式（1-3）和公式（1-4）带入公式（1-2），可得：
Figure 818293DEST_PATH_IMAGE016
（1-5）
S35、令
Figure DEST_PATH_IMAGE017
，并将公式（1-5）带入公式（1-1）中可得：
Figure 712079DEST_PATH_IMAGE018
。
在步骤S4中，根据坝基覆盖层深度、坝基材料的邓肯-张模型参数及变形监测点到参考面的距离，计算施工期坝体压缩产生的牵连变形S li ：
Figure DEST_PATH_IMAGE019
其中，H l 为坝基覆盖层深度，单位为m；n l 为坝基材料的邓肯-张模型参数；B l 为坝基变形系数；
在本发明的一个实施例中，所述施工期坝体压缩产生的牵连变形S li 的构建方法包括：
S41、当不考虑附加应力在坝基内的应力扩散衰减时，有覆盖层变形引起的坝体牵连变形为：
Figure 433523DEST_PATH_IMAGE020
（2-1）
式中，E ls 为坝基压缩模量；a l 为坝基内压应力修正系数；采用邓肯-张模型进行坝体变形计算填筑料的压缩模量E ls ：
Figure DEST_PATH_IMAGE021
（2-2）
S42、坝基内的大、小主应力
Figure 821254DEST_PATH_IMAGE022
的计算公式分别为：
Figure DEST_PATH_IMAGE023
S43、将公式（2-3）和公式（2-4）带入公式（2-2），可得：
Figure 6860DEST_PATH_IMAGE024
（2-5）
S44、令
Figure DEST_PATH_IMAGE025
，并将公式（2-5）带入公式（2-1）中可得：
Figure 274810DEST_PATH_IMAGE026
实施时，本方案优选所述坝基应力相关系数B l 的计算公式为：
Figure DEST_PATH_IMAGE027
其中，a l 为坝基压应力修正系数；β l 、η l 分别为坝基内的大主应力及小主应力修正系数；γ为填筑料的容重；k el 、R fl 分别为坝基材料邓肯-张模型参数；pa为大气压强。
在步骤S5中，根据材料的蠕变参数，计算施工期的蠕变变形
Figure 716506DEST_PATH_IMAGE028
：
Figure DEST_PATH_IMAGE029
式中，K p 为蠕变变形系数；α为坝体堆石料蠕变参数；t j 为坝体开始填筑至变形监测点获得初始测值的时间，单位为天；t c 为坝体开始填筑至进行变形预测的时间，单位为天；t p 为0至t j 中间过程中选定的时间节点，单位为天；p为时间节点变量，1≤p≤j；j为选取的时间节点个数，视精度和预测时间长短而定；e为自然数。
在本发明的一个实施例中，施工期的蠕变变形
Figure 376770DEST_PATH_IMAGE030
的构建方法包括：
在施工期，土石坝沉降变形由瞬时压缩变形和蠕变变形部分组成，t时刻蠕变量
Figure DEST_PATH_IMAGE031
可利用经验蠕变模型进行计算：
Figure 387189DEST_PATH_IMAGE032
（5-1）
在三轴蠕变试验中，轴向蠕变变形与体积蠕变变形和剪切蠕变变形之间满足如下关系：
Figure DEST_PATH_IMAGE033
（5-2）
其中，
Figure 91619DEST_PATH_IMAGE034
为最终蠕变量；
Figure 554699DEST_PATH_IMAGE035
为轴向蠕变量；
Figure 113594DEST_PATH_IMAGE036
为轴向蠕变量；
Figure 40837DEST_PATH_IMAGE037
为剪切蠕变量；
轴向蠕变计算公式为：
Figure 190933DEST_PATH_IMAGE038
（5-3）
其中，
Figure 734784DEST_PATH_IMAGE039
为坝体材料蠕变参数；q为偏应力；SL为应力水平；
假设中主应力可近似表示为：
Figure 362950DEST_PATH_IMAGE040
（5-4）
令
Figure 879119DEST_PATH_IMAGE041
Figure 456510DEST_PATH_IMAGE042
那么，式（5-3）中一点的最终轴向蠕变可近似表示如式（5-5）所示的形式：
Figure 159762DEST_PATH_IMAGE043
（5-5）
式中，z为空间点上覆土层厚度；
Figure 60460DEST_PATH_IMAGE044
为材料的蠕变参数；λ为中主应力系数。
参考图5，其示出了大坝填筑历时曲线，那么，施工期任意高程点
Figure 431136DEST_PATH_IMAGE045
的蠕变变形就等于该点以下土层在时间
Figure 923035DEST_PATH_IMAGE046
范围内所产生的蠕变变形之和，如式（5-6）所示。
Figure 379162DEST_PATH_IMAGE047
（5-6）
若用
Figure 614709DEST_PATH_IMAGE048
表示
Figure 371050DEST_PATH_IMAGE049
点以下任意位置z i 因填筑产生的最终蠕变变形增量，那么在
Figure 269736DEST_PATH_IMAGE050
范围内任意位置z i 发生的蠕变增量为：
Figure 502175DEST_PATH_IMAGE051
（5-7）
其中，
Figure 275834DEST_PATH_IMAGE052
为蠕变增量；
将坝体划分为N层，每层对应的坝体高度分别为
Figure 886682DEST_PATH_IMAGE053
，对应填筑时间为
Figure 720383DEST_PATH_IMAGE054
；以k表示p点以上的填筑层编号（见图5），那么第k层填筑对第p层产生的蠕变增量可用式（5-8）表示：
Figure 151102DEST_PATH_IMAGE055
（5-8）
其中，
Figure 728452DEST_PATH_IMAGE056
为第k层填筑对第p层产生的蠕变增量，
Figure 724965DEST_PATH_IMAGE057
、
Figure 198409DEST_PATH_IMAGE058
；如此，则最终蠕变量可用式（5-9）进行计算，式中
Figure 952362DEST_PATH_IMAGE059
为填筑到高程
Figure 35200DEST_PATH_IMAGE060
所对应的时间。
Figure 355060DEST_PATH_IMAGE061
（5-9）
其中，
Figure 264985DEST_PATH_IMAGE062
为第p层高度；
根据式（5-9），合并同类项，则任意层的沉降蠕变变形量可如式（5-10）计算：
Figure 404717DEST_PATH_IMAGE063
（5-10）
其中，
Figure 58290DEST_PATH_IMAGE064
为任意层的沉降蠕变变形量；
Figure 498237DEST_PATH_IMAGE065
分别为第p+1层，p+2层…p+n层对第p层蠕变变形的影响系数；
进一步合并式（5-10）中的同类项后可写成如式（5-11）所示的一般形式。
Figure 579063DEST_PATH_IMAGE066
（5-11）
其中，
Figure 760687DEST_PATH_IMAGE067
为第k层对第p层蠕变变形的影响系数；
Figure 686792DEST_PATH_IMAGE068
（5-12）
施工期的蠕变变形
Figure 748289DEST_PATH_IMAGE069
可进一步写成如下形式：
Figure 265596DEST_PATH_IMAGE070
式中，
Figure 379920DEST_PATH_IMAGE071
为蠕变变形系数；
Figure 172034DEST_PATH_IMAGE072
为坝体堆石料蠕变参数；t j 为坝体开始填筑至变形监测点获得初始测值的时间，单位为天；t c 为坝体开始填筑至进行变形预测的时间，单位为天；t p 为0至t j 中间过程中选定的时间节点，单位为天；p为时间节点变量，
Figure 852151DEST_PATH_IMAGE073
；
Figure 540359DEST_PATH_IMAGE075
为选取的时间节点个数，视精度和预测时间长短而定；e为自然数；
在步骤S6中，累加沉降变形S di 、牵连变形S li 和蠕变变形
Figure 407559DEST_PATH_IMAGE076
得到施工填筑期坝体任意监测点处的总沉降变形
Figure 292380DEST_PATH_IMAGE077
，S z0 为常数；
在步骤S7中，判断变形监测点i到参考面的距离是否小于变形监测点i处对应的浸润线高程H wi ，若是进入步骤S8，否则进入步骤S9；
在步骤S8中，根据变形监测点i到参考面的距离和变形预测点i到坝体表面的垂直距离，计算运行期坝体任意监测点的总沉降变形S zi ：
Figure 295846DEST_PATH_IMAGE078
其中，
Figure 420534DEST_PATH_IMAGE079
、
Figure 509451DEST_PATH_IMAGE080
均为与材料湿化变形系数；n w 为材料湿化变形参数；
在步骤S9中，根据变形预测点i处对应的浸润线高程和变形预测点i到坝体表面的垂直距离，计算运行期坝体任意监测点处的总沉降变形S zi ：
。
在本发明的一个实施例中，与材料湿化变形参数相关的回归系数
Figure 766918DEST_PATH_IMAGE080
、
Figure 796929DEST_PATH_IMAGE079
的计算公式分别为：
Figure 638721DEST_PATH_IMAGE082
其中，
Figure 865344DEST_PATH_IMAGE083
、
Figure 577823DEST_PATH_IMAGE084
和
Figure 513156DEST_PATH_IMAGE085
均为材料湿化变形参数；β和η分别为大主应力及小主应力系数；γ sat 为土体饱和容重；φ为坝体材料内摩擦角；Pa为大气压强。
为验证上述近似公式所反映坝体变形规律的合理性，对沉降变形S di 中的坝高幂函数
Figure 311085DEST_PATH_IMAGE086
和牵连变形S li 中覆盖层深度幂函数
Figure 521224DEST_PATH_IMAGE087
设置如下算例进行分析：
理想均质土石坝坝高100m，坝基覆盖层深度50m，由坝体压缩和覆盖层压缩引起的坝体随坝高的变化规律如图3所示（
Figure 117903DEST_PATH_IMAGE088
），图中横坐标为沉降变形的坝高幂函数和覆盖层深度幂函数计算值。
由图3显示的变形规律可知，因坝体压缩产生的沉降变形随计算点所在的坝高位置表现为先增后减的规律，而坝基压缩对坝体变形的影响表现为越远离大坝建基面，则因坝基覆盖层压缩产生的坝体牵连变形越小。为进一步验证计算成果的合理性，对上述算例采用FLAC3d进行了对比分析。结果表明，FLAC3d计算的坝体沉降变形随坝高的分布规律（见图4）与图3解析计算的分布规律一致，表明本方案在施工期引入坝体沉降变形物理力学成因，能够提高预测精度。
下面结合结合瀑布沟高心墙堆石坝，对本方案的预测方法的效果进行说明：
选取瀑布沟高心墙堆石坝的坝底作为参考面，选取瀑布沟高心墙堆石坝0+240m监测断面进行验证，在施工期阶段，选取变形监测点VE2-37在大坝施工期2008年02月23日~2009年01月05日期间及变形监测点VE1-18在大坝施工期2009年02月20日~2009年08月25日期间的相对位置、大坝的材料参数及这两个变形监测点的实测数据参考表1至表3：
表1 VE1-18、VE2-37参数表（固定参数）
Figure 692979DEST_PATH_IMAGE089
表2 VE1-18参数表（变化参数
Figure 978204DEST_PATH_IMAGE090
Figure 484313DEST_PATH_IMAGE091
表3 VE2-37参数表（变化参数）
Figure 968122DEST_PATH_IMAGE092
Figure 510837DEST_PATH_IMAGE093
基于变形监测点VE1-18和VE2-37的相对位置及大坝的材料参数，采用本方案的预测方法进行预测，变形监测点VE1-18和VE2-37的预测数据和实测数据分别参见图6和图7。
通过图6和图7可以得出，变形监测点计算值与实测值的平均相对误差均在10%以内，精度较高，可以采用本方案的方法得到的施工期预测数据与实测值历时过程对比示例，本方案从监测点变形序列的发展过程来看，本方案的预测方法能较好地描述测点变形的历时过程，预测效果较好。
在运行期间，选取变形监测点VE1-25在大坝施工期2009年10月11日~2012年10月30日期间及变形监测点LD59在大坝施工期2009年11月06日~2010年10月05日期间的水位高及材料湿化变形参数，变形监测点到坝表面的高度及两个变形监测点的实测数据表4至表6：
表4 VE1-25、LD59参数表（固定参数）表中无单位
Figure 548938DEST_PATH_IMAGE094
表5 VE1-25参数表（变化参数）
Figure 304142DEST_PATH_IMAGE095
Figure 111299DEST_PATH_IMAGE096
Figure 795221DEST_PATH_IMAGE097
Figure 586197DEST_PATH_IMAGE098
Figure 410671DEST_PATH_IMAGE099
Figure 158089DEST_PATH_IMAGE100
表6 LD59参数表（变化参数）
Figure 42606DEST_PATH_IMAGE101
Figure 789720DEST_PATH_IMAGE102
Figure 417885DEST_PATH_IMAGE103
基于获取的数据采用本方案的预测方法进行预测，变形监测点VE1-25和LD59和预测数据参考图8和图9。
通过图8和图9可以得出，采用本方案的预测方法得到的预测值与实测值的相对误差在10%以内，精度较高，满足工程应用需求。

Claims (5)

1.高心墙堆石坝沉降变形的预测方法，其特征在于，包括：
S1、选取高心墙堆石坝的建基面作为参考面，根据坝体内已布设的沉降变形监测点及其位置确定各监测点到参考面的距离；
S2、判断高心墙堆石坝为施工期或运行期，若为施工期，则进入步骤S3，若为运行期，则进入步骤S7；
S3、根据坝体材料的邓肯-张模型参数及监测点到参考面的距离，计算施工期坝体压缩产生的沉降变形S di ：
Figure 784935DEST_PATH_IMAGE001
其中，H i 为第i个变形监测点到坝表面的垂直距离，单位为m；z i 为第i个变形监测点到参考面的距离，单位为m；n为坝体材料的邓肯-张模型参数；B为坝体变形系数；
S4、根据坝基覆盖层深度、坝基材料的邓肯-张模型参数及变形监测点到参考面的距离，计算施工期坝体压缩产生的牵连变形S li ：
Figure 131734DEST_PATH_IMAGE002
其中，H l 为坝基覆盖层深度，单位为m；n l 为坝基材料的邓肯-张模型参数；B l 为坝基变形系数；
S5、根据材料的蠕变参数，计算施工期的蠕变变形
Figure 560179DEST_PATH_IMAGE003
：
Figure 128039DEST_PATH_IMAGE004
式中，K p 为蠕变变形系数；α为坝体堆石料蠕变参数；t j 为坝体开始填筑至变形监测点获得初始测值的时间，单位为天；t c 为坝体开始填筑至进行变形预测的时间，单位为天；t p 为0至t j 中间过程中选定的时间节点，单位为天；p为时间节点变量，1≤p≤j；j为选取的时间节点个数；e为自然数；
S6、累加沉降变形S di 、牵连变形S li 和蠕变变形
Figure 212407DEST_PATH_IMAGE003
得到施工填筑期坝体任意监测点处的总沉降变形
Figure 628476DEST_PATH_IMAGE005
，S z0 为常数；
S7、判断变形监测点i到参考面的距离是否小于变形监测点i处对应的浸润线高程H wi ，若是进入步骤S8，否则进入步骤S9；
S8、根据变形监测点i到参考面的距离和变形预测点i到坝体表面的垂直距离，计算运行期坝体任意监测点的总沉降变形S zi ：
Figure 911428DEST_PATH_IMAGE006
其中，
Figure 384610DEST_PATH_IMAGE007
、
Figure 753012DEST_PATH_IMAGE008
均为与材料湿化变形系数；n w 为材料湿化变形参数；
S9、根据变形预测点i处对应的浸润线高程和变形预测点i到坝体表面的垂直距离，计算运行期坝体任意监测点处的总沉降变形S zi ：
Figure 238351DEST_PATH_IMAGE009
；
与材料湿化变形参数相关的回归系数
Figure 110230DEST_PATH_IMAGE010
、
Figure 512172DEST_PATH_IMAGE011
的计算公式分别为：
Figure 367870DEST_PATH_IMAGE012
其中，
Figure 860162DEST_PATH_IMAGE013
、
Figure 852127DEST_PATH_IMAGE014
和
Figure 434156DEST_PATH_IMAGE015
均为材料湿化变形参数；β和η分别为大主应力及小主应力系数；γ sat 为土体饱和容重；φ为坝体材料内摩擦角；Pa为大气压强。
2.根据权利要求1所述的高心墙堆石坝沉降变形的预测方法，其特征在于，所述坝体应力相关系数B的计算公式分别为：
Figure 275685DEST_PATH_IMAGE016
其中，α d 为坝体压应力修正系数；k e 、R f 分别为坝体材料邓肯-张模型参数；γ为填筑料的容重；β、η分别为坝体内的大主应力及小主应力修正系数；φ为坝体材料内摩擦角。
3.根据权利要求1所述的高心墙堆石坝沉降变形的预测方法，其特征在于，所述坝基应力相关系数B l 的计算公式为：
Figure 335783DEST_PATH_IMAGE017
其中，a l 为坝基压应力修正系数；β l 、η l 分别为坝基内的大主应力及小主应力修正系数；γ为填筑料的容重；k el 、R fl 分别为坝基材料邓肯-张模型参数;
Figure 949298DEST_PATH_IMAGE018
为大气压强；φ l 为坝基材料内摩擦角。
4.根据权利要求1-3任一所述的高心墙堆石坝沉降变形的预测方法，其特征在于，所述施工期坝体压缩产生的变形S di 的构建方法包括：
构建坝体内任意一点的压缩变形的初始模型：
Figure 702228DEST_PATH_IMAGE019
（1-1）
式中，α d 为坝体压应力修正系数；γ为填筑料的容重；E s 为填筑料的压缩模量；z 1，z 2均为积分变量；
采用邓肯-张模型进行坝体变形计算填筑料的压缩模量E s ：
Figure 266939DEST_PATH_IMAGE020
（1-2）
其中，
Figure 918280DEST_PATH_IMAGE021
分别为坝体内的大、小主应力；
坝体内的大、小主应力
Figure 415995DEST_PATH_IMAGE021
的计算公式分别为：
Figure 575712DEST_PATH_IMAGE022
将公式（1-3）和公式（1-4）带入公式（1-2），可得：
Figure 627719DEST_PATH_IMAGE023
（1-5）
令
Figure 262575DEST_PATH_IMAGE024
，并将公式（1-5）带入公式（1-1）中可得：
Figure 83639DEST_PATH_IMAGE025
。
5.根据权利要求1-3任一所述的高心墙堆石坝沉降变形的预测方法，其特征在于，所述施工期坝体压缩产生的牵连变形S li 的构建方法包括：
当不考虑附加应力在坝基内的应力扩散衰减时，有覆盖层变形引起的坝体牵连变形为：
Figure 414257DEST_PATH_IMAGE026
（2-1）
式中，E ls 为坝基压缩模量；a l 为坝基内压应力修正系数；采用邓肯-张模型进行坝体变形计算填筑料的压缩模量E ls ：
Figure 953560DEST_PATH_IMAGE027
（2-2）
坝基内的大、小主应力
Figure 159152DEST_PATH_IMAGE028
的计算公式分别为：
Figure 333256DEST_PATH_IMAGE029
将公式（2-3）和公式（2-4）带入公式（2-2），可得：
Figure 622328DEST_PATH_IMAGE030
（2-5）
令
Figure 150392DEST_PATH_IMAGE031
，并将公式（2-5）带入公式（2-1）中可得：
Figure 894095DEST_PATH_IMAGE032
。
CN202010922824.XA 2020-09-04 2020-09-04 高心墙堆石坝沉降变形的预测方法 Active CN111811468B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010922824.XA CN111811468B (zh) 2020-09-04 2020-09-04 高心墙堆石坝沉降变形的预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010922824.XA CN111811468B (zh) 2020-09-04 2020-09-04 高心墙堆石坝沉降变形的预测方法

Publications (2)

Publication Number Publication Date
CN111811468A CN111811468A (zh) 2020-10-23
CN111811468B true CN111811468B (zh) 2020-12-08

Family

ID=72860638

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010922824.XA Active CN111811468B (zh) 2020-09-04 2020-09-04 高心墙堆石坝沉降变形的预测方法

Country Status (1)

Country Link
CN (1) CN111811468B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113449437B (zh) * 2021-07-26 2022-04-12 青岛理工大学 一种黏土斜墙土石坝渗漏预警方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109917450A (zh) * 2019-04-10 2019-06-21 长江水利委员会长江科学院 一种全断面网络化布设的土石坝位移监测系统及方法
CN110987618A (zh) * 2019-12-13 2020-04-10 中国电建集团中南勘测设计研究院有限公司 一种堆石料试验装置及其方法
CN111210519A (zh) * 2019-12-25 2020-05-29 四川大学 一种高心墙堆石坝变形参数反演方法
CN111475981A (zh) * 2020-04-22 2020-07-31 水发规划设计有限公司 一种水库土石坝渗透变形和稳定性的数值模拟方法及装置

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN203249668U (zh) * 2013-05-07 2013-10-23 山东大学 一种坝体浸润线监测模型试验装置
CN104196010A (zh) * 2014-08-14 2014-12-10 中国电建集团西北勘测设计研究院有限公司 一种面板堆石坝强制预沉降施工方法
JP6696290B2 (ja) * 2016-04-28 2020-05-20 株式会社大林組 ロックフィルダムにおけるロックゾーンの品質管理方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109917450A (zh) * 2019-04-10 2019-06-21 长江水利委员会长江科学院 一种全断面网络化布设的土石坝位移监测系统及方法
CN110987618A (zh) * 2019-12-13 2020-04-10 中国电建集团中南勘测设计研究院有限公司 一种堆石料试验装置及其方法
CN111210519A (zh) * 2019-12-25 2020-05-29 四川大学 一种高心墙堆石坝变形参数反演方法
CN111475981A (zh) * 2020-04-22 2020-07-31 水发规划设计有限公司 一种水库土石坝渗透变形和稳定性的数值模拟方法及装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
"面板堆石坝堆石体湿化变形分析方法研究";方国宝;《中国优秀硕士学位论文全文数据库 工程科技II辑》;20071215(第06期);正文第3.3.3小节，第3.4节 *
"高心墙堆石坝施工期变形时空预测模型研究";卢祥 等;《工程科学与技术》;20170731;第49卷(第4期);正文第1.2.1小节、图1,2 *

Also Published As

Publication number Publication date
CN111811468A (zh) 2020-10-23

Similar Documents

Publication Publication Date Title
CN107016185B (zh) 一种冰碛湖溃决洪水峰值流量预测的计算方法
CN108487904B (zh) 一种基于图版的消除末端效应的相渗曲线校正方法
CN107515025B (zh) 古城墙病害等级评估方法
Göğüş et al. Broad-crested weirs with rectangular compound cross sections
CN111811468B (zh) 高心墙堆石坝沉降变形的预测方法
Bian et al. Evaluating the effect of soil structure on the ground response during shield tunnelling in Shanghai soft clay
CN105155425A (zh) 满堂支架流体仿真模拟预压施工方法
Shi et al. A statistical model of deformation during the construction of a concrete face rockfill dam
CN104111137B (zh) 基于三次多项式模型的锚杆承载力预测方法
CN113987663A (zh) 基于实测数据考虑计算参数变化的软土路基沉降预测方法
CN114493702A (zh) 一种基于有限元分析及地理信息的农村房产损失评估方法
CN110969370B (zh) 建筑结构构件的质量风险分析方法
CN111950065A (zh) 高心墙堆石坝施工期水平变形的预测方法
Wang et al. Dynamic tensile test of mass concrete with Shapai Dam cores
De Vos et al. The current draft of the application rules of Eurocode 7 in Belgium for the design of pile foundations
Jostad et al. Cyclic capacity of monopiles in sand under Partially drained conditions: a numerical approach
CN110197015A (zh) 一种坝基预应力锚索有效拉应力测定方法
CN114662341A (zh) 岩体临界滑动面极限分析方法
Abedian et al. Settlement Evaluation of a Concrete Face Rock-Fill Dam (CFRD) Using a Back-Analysis Method Based on Measurement Results (A Case Study of Siah-Bisheh Dam)
CN111898176B (zh) 一种成层地基逐渐加荷条件下固结度的计算方法
ERGUN et al. COMPARISON OF GEOMETRICAL IMPERFECTION DEFINITIONS ON ENCASED COMPOSITE COLUMNS
RU2457455C2 (ru) Способ определения герметичности запорных арматур нефтепровода
Mozer et al. Time-dependent behavior of concrete beams
CN110020792A (zh) 基于洪水洪峰-洪量联合的高堆石坝施工度汛风险预测方法
Wessels et al. Flow-gauging structures in South African rivers part 2: calibration

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