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

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

Info

Publication number
CN111811468A
CN111811468A CN202010922824.XA CN202010922824A CN111811468A CN 111811468 A CN111811468 A CN 111811468A CN 202010922824 A CN202010922824 A CN 202010922824A CN 111811468 A CN111811468 A CN 111811468A
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.)
Granted
Application number
CN202010922824.XA
Other languages
English (en)
Other versions
CN111811468B (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 423367DEST_PATH_IMAGE001
其中，H i 为第i个变形监测点到坝表面的垂直距离，单位为m；z i 为第i个变形监测点到参考面的距离，单位为m；n为坝体材料的邓肯-张模型参数；B为坝体变形系数；
S4、根据坝基覆盖层深度、坝基材料的邓肯-张模型参数及变形监测点到参考面的距离，计算施工期坝体压缩产生的牵连变形S li ：
Figure 419136DEST_PATH_IMAGE003
其中，H l 为坝基覆盖层深度，单位为m；n l 为坝基材料的邓肯-张模型参数；B l 为坝基变形系数；
S5、根据材料的蠕变参数，计算施工期的蠕变变形
Figure 402135DEST_PATH_IMAGE004
：
Figure 607DEST_PATH_IMAGE005
式中，K p 为蠕变变形系数；α为坝体堆石料蠕变参数；t j 为坝体开始填筑至变形监测点获得初始测值的时间，单位为天；t c 为坝体开始填筑至进行变形预测的时间，单位为天；t p 为0至t j 中间过程中选定的时间节点，单位为天；p为时间节点变量，1≤p≤j；j为选取的时间节点个数，视精度和预测时间长短而定；e为自然数；
S6、累加沉降变形S di 、牵连变形S li 和蠕变变形
Figure 18241DEST_PATH_IMAGE006
得到施工填筑期坝体任意监测点处的总沉降变形
Figure 919332DEST_PATH_IMAGE007
，S z0 为常数；
S7、判断变形监测点i到参考面的距离是否小于变形监测点i处对应的浸润线高程H wi ，若是进入步骤S8，否则进入步骤S9；
S8、根据变形监测点i到参考面的距离和变形预测点i到坝体表面的垂直距离，计算运行期坝体任意监测点的总沉降变形S zi ：
Figure 389628DEST_PATH_IMAGE008
其中，
Figure 791790DEST_PATH_IMAGE009
、
Figure 663932DEST_PATH_IMAGE010
均为与材料湿化变形系数；n w 为材料湿化变形参数；
S9、根据变形预测点i处对应的浸润线高程和变形预测点i到坝体表面的垂直距离，计算运行期坝体任意监测点处的总沉降变形S zi ：
Figure 735924DEST_PATH_IMAGE011
。
本发明的有益效果为：本方案在进行沉降变形预测时，首先获取预测点的相对位置关系，并基于相对位置关系可以得到监测点所在位置的应力应变状态，最终通过应力应变关系得到每个预测点的沉降变形；由于本方案在预测时充分考虑高心墙堆石坝上述的物理成因，相对现有技术中数学模拟模型而言，综合考虑了预测点的位置关系和应力应变状态，使得本方案的沉降变形预测模型能较好地模拟和反映测点变形趋势，显著提升变形预测精度。
附图说明
图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 427936DEST_PATH_IMAGE012
其中，H i 为第i个变形监测点到坝表面的垂直距离，单位为m；z i 为第i个变形监测点到参考面的距离，单位为m；n为坝体材料的邓肯-张模型参数；B为坝体变形系数；
在本发明的一个实施例中，所述施工期坝体压缩产生的变形S di 的构建方法包括：
S31、构建坝体内任意一点的压缩变形的初始模型：
Figure 368211DEST_PATH_IMAGE013
（1-1）
式中，α d 为坝体压应力修正系数；γ为填筑料的容重；E s 为填筑料的压缩模量；z 1，z 2均为积分变量；
S32、采用邓肯-张模型进行坝体变形计算填筑料的压缩模量E s ：
Figure 94858DEST_PATH_IMAGE014
（1-2）
其中，
Figure 739684DEST_PATH_IMAGE015
分别为坝体内的大、小主应力；
S33、坝体内的大、小主应力
Figure 325517DEST_PATH_IMAGE016
的计算公式分别为：
Figure 803903DEST_PATH_IMAGE017
S34、将公式（1-3）和公式（1-4）带入公式（1-2），可得：
Figure 650636DEST_PATH_IMAGE019
（1-5）
S35、令
Figure 923485DEST_PATH_IMAGE020
，并将公式（1-5）带入公式（1-1）中可得：
Figure 731036DEST_PATH_IMAGE021
。
在步骤S4中，根据坝基覆盖层深度、坝基材料的邓肯-张模型参数及变形监测点到参考面的距离，计算施工期坝体压缩产生的牵连变形S li ：
Figure 13112DEST_PATH_IMAGE022
其中，H l 为坝基覆盖层深度，单位为m；n l 为坝基材料的邓肯-张模型参数；B l 为坝基变形系数；
在本发明的一个实施例中，所述施工期坝体压缩产生的牵连变形S li 的构建方法包括：
S41、当不考虑附加应力在坝基内的应力扩散衰减时，有覆盖层变形引起的坝体牵连变形为：
Figure 714352DEST_PATH_IMAGE023
（2-1）
式中，E ls 为坝基压缩模量；a l 为坝基内压应力修正系数；采用邓肯-张模型进行坝体变形计算填筑料的压缩模量E ls ：
Figure 689261DEST_PATH_IMAGE024
（2-2）
S42、坝基内的大、小主应力
Figure 718528DEST_PATH_IMAGE025
的计算公式分别为：
Figure 538717DEST_PATH_IMAGE026
S43、将公式（2-3）和公式（2-4）带入公式（2-2），可得：
Figure 625621DEST_PATH_IMAGE027
（2-5）
S44、令
Figure 646798DEST_PATH_IMAGE028
，并将公式（2-5）带入公式（2-1）中可得：
Figure 287995DEST_PATH_IMAGE029
实施时，本方案优选所述坝基应力相关系数B l 的计算公式为：
Figure 646295DEST_PATH_IMAGE030
其中，a l 为坝基压应力修正系数；β l 、η l 分别为坝基内的大主应力及小主应力修正系数；γ为填筑料的容重；k el 、R fl 分别为坝基材料邓肯-张模型参数；pa为大气压强。
在步骤S5中，根据材料的蠕变参数，计算施工期的蠕变变形
Figure 322127DEST_PATH_IMAGE031
：
Figure 638839DEST_PATH_IMAGE032
式中，K p 为蠕变变形系数；α为坝体堆石料蠕变参数；t j 为坝体开始填筑至变形监测点获得初始测值的时间，单位为天；t c 为坝体开始填筑至进行变形预测的时间，单位为天；t p 为0至t j 中间过程中选定的时间节点，单位为天；p为时间节点变量，1≤p≤j；j为选取的时间节点个数，视精度和预测时间长短而定；e为自然数。
在本发明的一个实施例中，施工期的蠕变变形
Figure 642698DEST_PATH_IMAGE033
的构建方法包括：
在施工期，土石坝沉降变形由瞬时压缩变形和蠕变变形部分组成，t时刻蠕变量
Figure 804689DEST_PATH_IMAGE034
可利用经验蠕变模型进行计算：
Figure 69449DEST_PATH_IMAGE035
（5-1）
在三轴蠕变试验中，轴向蠕变变形与体积蠕变变形和剪切蠕变变形之间满足如下关系：
Figure 557062DEST_PATH_IMAGE036
（5-2）
其中，
Figure 907272DEST_PATH_IMAGE037
为最终蠕变量；
Figure 138533DEST_PATH_IMAGE038
为轴向蠕变量；
Figure 398744DEST_PATH_IMAGE039
为轴向蠕变量；
Figure 791679DEST_PATH_IMAGE040
为剪切蠕变量；
轴向蠕变计算公式为：
Figure 629185DEST_PATH_IMAGE041
（5-3）
其中，
Figure 386839DEST_PATH_IMAGE042
为坝体材料蠕变参数；q为偏应力；SL为应力水平；
假设中主应力可近似表示为：
Figure 891770DEST_PATH_IMAGE043
（5-4）
令
Figure 330973DEST_PATH_IMAGE044
Figure 921354DEST_PATH_IMAGE045
那么，式（5-3）中一点的最终轴向蠕变可近似表示如式（5-5）所示的形式：
Figure 494418DEST_PATH_IMAGE046
（5-5）
式中，z为空间点上覆土层厚度；
Figure 322697DEST_PATH_IMAGE047
为材料的蠕变参数；λ为中主应力系数。
参考图5，其示出了大坝填筑历时曲线，那么，施工期任意高程点
Figure 323014DEST_PATH_IMAGE048
的蠕变变形就等于该点以下土层在时间
Figure 10478DEST_PATH_IMAGE049
范围内所产生的蠕变变形之和，如式（5-6）所示。
Figure 387233DEST_PATH_IMAGE050
（5-6）
若用
Figure 601176DEST_PATH_IMAGE051
表示
Figure 382182DEST_PATH_IMAGE052
点以下任意位置z i 因填筑产生的最终蠕变变形增量，那么在
Figure 415997DEST_PATH_IMAGE053
范围内任意位置z i 发生的蠕变增量为：
Figure 330863DEST_PATH_IMAGE054
（5-7）
其中，
Figure 274680DEST_PATH_IMAGE055
为蠕变增量；
将坝体划分为N层，每层对应的坝体高度分别为
Figure 616799DEST_PATH_IMAGE056
，对应填筑时间为
Figure 403490DEST_PATH_IMAGE057
；以k表示p点以上的填筑层编号（见图5），那么第k层填筑对第p层产生的蠕变增量可用式（5-8）表示：
Figure 856468DEST_PATH_IMAGE058
（5-8）
其中，
Figure 779424DEST_PATH_IMAGE059
为第k层填筑对第p层产生的蠕变增量，
Figure 167812DEST_PATH_IMAGE060
、
Figure 707377DEST_PATH_IMAGE061
；如此，则最终蠕变量可用式（5-9）进行计算，式中
Figure 698467DEST_PATH_IMAGE062
为填筑到高程
Figure 741509DEST_PATH_IMAGE063
所对应的时间。
Figure 300798DEST_PATH_IMAGE064
（5-9）
其中，
Figure 796501DEST_PATH_IMAGE065
为第p层高度；
根据式（5-9），合并同类项，则任意层的沉降蠕变变形量可如式（5-10）计算：
Figure 591282DEST_PATH_IMAGE066
（5-10）
其中，
Figure 766129DEST_PATH_IMAGE067
为任意层的沉降蠕变变形量；
Figure 496318DEST_PATH_IMAGE068
分别为第p+1层，p+2层…p+n层对第p层蠕变变形的影响系数；
进一步合并式（5-10）中的同类项后可写成如式（5-11）所示的一般形式。
Figure 948160DEST_PATH_IMAGE069
（5-11）
其中，
Figure 953156DEST_PATH_IMAGE070
为第k层对第p层蠕变变形的影响系数；
Figure 705211DEST_PATH_IMAGE071
（5-12）
施工期的蠕变变形
Figure 137461DEST_PATH_IMAGE072
可进一步写成如下形式：
Figure 342177DEST_PATH_IMAGE073
式中，
Figure 744339DEST_PATH_IMAGE074
为蠕变变形系数；
Figure 350901DEST_PATH_IMAGE075
为坝体堆石料蠕变参数；t j 为坝体开始填筑至变形监测点获得初始测值的时间，单位为天；t c 为坝体开始填筑至进行变形预测的时间，单位为天；t p 为0至t j 中间过程中选定的时间节点，单位为天；p为时间节点变量，
Figure 813107DEST_PATH_IMAGE076
；
Figure 646065DEST_PATH_IMAGE077
为选取的时间节点个数，视精度和预测时间长短而定；e为自然数；
在步骤S6中，累加沉降变形S di 、牵连变形S li 和蠕变变形
Figure 586339DEST_PATH_IMAGE078
得到施工填筑期坝体任意监测点处的总沉降变形
Figure 312986DEST_PATH_IMAGE079
，S z0 为常数；
在步骤S7中，判断变形监测点i到参考面的距离是否小于变形监测点i处对应的浸润线高程H wi ，若是进入步骤S8，否则进入步骤S9；
在步骤S8中，根据变形监测点i到参考面的距离和变形预测点i到坝体表面的垂直距离，计算运行期坝体任意监测点的总沉降变形S zi ：
Figure 946093DEST_PATH_IMAGE080
其中，
Figure 266347DEST_PATH_IMAGE081
、
Figure 744733DEST_PATH_IMAGE082
均为与材料湿化变形系数；n w 为材料湿化变形参数；
在步骤S9中，根据变形预测点i处对应的浸润线高程和变形预测点i到坝体表面的垂直距离，计算运行期坝体任意监测点处的总沉降变形S zi ：
Figure 591466DEST_PATH_IMAGE083
。
在本发明的一个实施例中，与材料湿化变形参数相关的回归系数
Figure 129895DEST_PATH_IMAGE084
、
Figure 937445DEST_PATH_IMAGE085
的计算公式分别为：
Figure 219522DEST_PATH_IMAGE086
其中，
Figure 920762DEST_PATH_IMAGE087
、
Figure 630092DEST_PATH_IMAGE088
和
Figure 924938DEST_PATH_IMAGE089
均为材料湿化变形参数；β和η分别为大主应力及小主应力系数；γ sat 为土体饱和容重；φ为坝体材料内摩擦角；Pa为大气压强。
为验证上述近似公式所反映坝体变形规律的合理性，对沉降变形S di 中的坝高幂函数
Figure 745126DEST_PATH_IMAGE090
和牵连变形S li 中覆盖层深度幂函数
Figure 289154DEST_PATH_IMAGE091
设置如下算例进行分析：
理想均质土石坝坝高100m，坝基覆盖层深度50m，由坝体压缩和覆盖层压缩引起的坝体随坝高的变化规律如图3所示（
Figure 434964DEST_PATH_IMAGE092
），图中横坐标为沉降变形的坝高幂函数和覆盖层深度幂函数计算值。
由图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 810582DEST_PATH_IMAGE094
表2 VE1-18参数表（变化参数
Figure 575407DEST_PATH_IMAGE096
Figure 985660DEST_PATH_IMAGE098
表3 VE2-37参数表（变化参数）
Figure 36792DEST_PATH_IMAGE100
Figure 165285DEST_PATH_IMAGE102
基于变形监测点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 592855DEST_PATH_IMAGE104
表5 VE1-25参数表（变化参数）
Figure 998560DEST_PATH_IMAGE106
Figure 486173DEST_PATH_IMAGE108
Figure 836383DEST_PATH_IMAGE110
Figure DEST_PATH_IMAGE112
Figure 5327DEST_PATH_IMAGE114
Figure 265539DEST_PATH_IMAGE116
表6 LD59参数表（变化参数）
Figure 189632DEST_PATH_IMAGE118
Figure 27138DEST_PATH_IMAGE120
Figure 796511DEST_PATH_IMAGE122
基于获取的数据采用本方案的预测方法进行预测，变形监测点VE1-25和LD59和预测数据参考图8和图9。
通过图8和图9可以得出，采用本方案的预测方法得到的预测值与实测值的相对误差在10%以内，精度较高，满足工程应用需求。

Claims (6)

1.高心墙堆石坝沉降变形的预测方法，其特征在于，包括：
S1、选取高心墙堆石坝的建基面作为参考面，根据坝体内已布设的沉降变形监测点及其位置确定各监测点到参考面的距离；
S2、判断高心墙堆石坝为施工期或运行期，若为施工期，则进入步骤S3，若为运行期，则进入步骤S7；
S3、根据坝体材料的邓肯-张模型参数及监测点到参考面的距离，计算施工期坝体压缩产生的沉降变形S di ：
Figure 993176DEST_PATH_IMAGE001
其中，H i 为第i个变形监测点到坝表面的垂直距离，单位为m；z i 为第i个变形监测点到参考面的距离，单位为m；n为坝体材料的邓肯-张模型参数；B为坝体变形系数；
S4、根据坝基覆盖层深度、坝基材料的邓肯-张模型参数及变形监测点到参考面的距离，计算施工期坝体压缩产生的牵连变形S li ：
Figure 254524DEST_PATH_IMAGE003
其中，H l 为坝基覆盖层深度，单位为m；n l 为坝基材料的邓肯-张模型参数；B l 为坝基变形系数；
S5、根据材料的蠕变参数，计算施工期的蠕变变形
Figure 96578DEST_PATH_IMAGE004
：
Figure 835995DEST_PATH_IMAGE006
式中，K p 为蠕变变形系数；α为坝体堆石料蠕变参数；t j 为坝体开始填筑至变形监测点获得初始测值的时间，单位为天；t c 为坝体开始填筑至进行变形预测的时间，单位为天；t p 为0至t j 中间过程中选定的时间节点，单位为天；p为时间节点变量，1≤p≤j；j为选取的时间节点个数；e为自然数；
S6、累加沉降变形S di 、牵连变形S li 和蠕变变形
Figure 260154DEST_PATH_IMAGE004
得到施工填筑期坝体任意监测点处的总沉降变形
Figure 879354DEST_PATH_IMAGE008
，S z0 为常数；
S7、判断变形监测点i到参考面的距离是否小于变形监测点i处对应的浸润线高程H wi ，若是进入步骤S8，否则进入步骤S9；
S8、根据变形监测点i到参考面的距离和变形预测点i到坝体表面的垂直距离，计算运行期坝体任意监测点的总沉降变形S zi ：
Figure 756175DEST_PATH_IMAGE010
其中，
Figure 299283DEST_PATH_IMAGE011
、
Figure 561637DEST_PATH_IMAGE012
均为与材料湿化变形系数；n w 为材料湿化变形参数；
S9、根据变形预测点i处对应的浸润线高程和变形预测点i到坝体表面的垂直距离，计算运行期坝体任意监测点处的总沉降变形S zi ：
Figure 899208DEST_PATH_IMAGE013
。
2.根据权利要求1所述的高心墙堆石坝沉降变形的预测方法，其特征在于，所述坝体应力相关系数B的计算公式分别为：
Figure 201008DEST_PATH_IMAGE015
其中，α d 为坝体压应力修正系数；k e 、R f 分别为坝体材料邓肯-张模型参数；γ为填筑料的容重；β、η分别为坝体内的大主应力及小主应力修正系数；φ为坝体材料内摩擦角。
3.根据权利要求1所述的高心墙堆石坝沉降变形的预测方法，其特征在于，所述坝基应力相关系数B l 的计算公式为：
Figure 531495DEST_PATH_IMAGE017
其中，a l 为坝基压应力修正系数；β l 、η l 分别为坝基内的大主应力及小主应力修正系数；γ为填筑料的容重；k el 、R fl 分别为坝基材料邓肯-张模型参数;
Figure 664667DEST_PATH_IMAGE018
为大气压强；φ l 为坝基材料内摩擦角。
4.根据权利要求1所述的高心墙堆石坝沉降变形的预测方法，其特征在于，与材料湿化变形参数相关的回归系数
Figure 438719DEST_PATH_IMAGE020
、
Figure 477082DEST_PATH_IMAGE022
的计算公式分别为：
Figure 373712DEST_PATH_IMAGE024
其中，
Figure 564653DEST_PATH_IMAGE025
、
Figure 758874DEST_PATH_IMAGE026
和
Figure 35265DEST_PATH_IMAGE027
均为材料湿化变形参数；β和η分别为大主应力及小主应力系数；γ sat 为土体饱和容重；φ为坝体材料内摩擦角；Pa为大气压强。
5.根据权利要求1-4任一所述的高心墙堆石坝沉降变形的预测方法，其特征在于，所述施工期坝体压缩产生的变形S di 的构建方法包括：
构建坝体内任意一点的压缩变形的初始模型：
Figure 458287DEST_PATH_IMAGE029
（1-1）
式中，α d 为坝体压应力修正系数；γ为填筑料的容重；E s 为填筑料的压缩模量；z 1，z 2均为积分变量；
采用邓肯-张模型进行坝体变形计算填筑料的压缩模量E s ：
Figure 815319DEST_PATH_IMAGE030
（1-2）
其中，
Figure 931174DEST_PATH_IMAGE032
分别为坝体内的大、小主应力；
坝体内的大、小主应力
Figure 429283DEST_PATH_IMAGE033
的计算公式分别为：
Figure 170843DEST_PATH_IMAGE035
将公式（1-3）和公式（1-4）带入公式（1-2），可得：
Figure 70797DEST_PATH_IMAGE037
（1-5）
令
Figure DEST_PATH_IMAGE039
，并将公式（1-5）带入公式（1-1）中可得：
Figure DEST_PATH_IMAGE041
。
6.根据权利要求1-4任一所述的高心墙堆石坝沉降变形的预测方法，其特征在于，所述施工期坝体压缩产生的牵连变形S li 的构建方法包括：
当不考虑附加应力在坝基内的应力扩散衰减时，有覆盖层变形引起的坝体牵连变形为：
Figure DEST_PATH_IMAGE043
（2-1）
式中，E ls 为坝基压缩模量；a l 为坝基内压应力修正系数；采用邓肯-张模型进行坝体变形计算填筑料的压缩模量E ls ：
Figure DEST_PATH_IMAGE045
（2-2）
坝基内的大、小主应力
Figure DEST_PATH_IMAGE046
的计算公式分别为：
Figure DEST_PATH_IMAGE048
将公式（2-3）和公式（2-4）带入公式（2-2），可得：
Figure DEST_PATH_IMAGE050
（2-5）
令
Figure DEST_PATH_IMAGE052
，并将公式（2-5）带入公式（2-1）中可得：
Figure DEST_PATH_IMAGE054
。
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 true CN111811468A (zh) 2020-10-23
CN111811468B 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)

Cited By (1)

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

Citations (7)

* 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 中国电建集团西北勘测设计研究院有限公司 一种面板堆石坝强制预沉降施工方法
JP2017198017A (ja) * 2016-04-28 2017-11-02 株式会社大林組 ロックフィルダムにおけるロックゾーンの品質管理方法
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 水发规划设计有限公司 一种水库土石坝渗透变形和稳定性的数值模拟方法及装置

Patent Citations (7)

* 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 中国电建集团西北勘测设计研究院有限公司 一种面板堆石坝强制预沉降施工方法
JP2017198017A (ja) * 2016-04-28 2017-11-02 株式会社大林組 ロックフィルダムにおけるロックゾーンの品質管理方法
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 (7)

* Cited by examiner, † Cited by third party
Title
HAO-FENG XING ET AL.: ""Construction of concrete-faced rockfill dams with weak rocks"", 《JOURNAL OF GEOTECHNICAL AND GEOENVIRONMENTAL ENGINEERING》 *
XIONGXIONG ZHOU ET AL.: ""Study on wetting deformation characteristics of coarse granular materials and its simulation in core-wall rockfill dams"", 《INT.J.NUMER.ANAL.METHODS GEOMECH.》 *
XIONGXIONG ZHOU ET AL.: ""Wetting deformation of core-wall rockfill dams"", 《INT.J.GEOMECH》 *
刘东升 等: ""高心墙堆石坝蓄水期沉降统计预测模型及应用"", 《人民长江》 *
卢祥 等: ""高心墙堆石坝施工期变形时空预测模型研究"", 《工程科学与技术》 *
方国宝: ""面板堆石坝堆石体湿化变形分析方法研究"", 《中国优秀硕士学位论文全文数据库 工程科技II辑》 *
陆阳洋 等: ""湿化变形对斜心墙堆石坝受力变形特性的影响"", 《扬州大学学报（自然科学版）》 *

Cited By (1)

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

Also Published As

Publication number Publication date
CN111811468B (zh) 2020-12-08

Similar Documents

Publication Publication Date Title
CN107016185B (zh) 一种冰碛湖溃决洪水峰值流量预测的计算方法
CN108487904B (zh) 一种基于图版的消除末端效应的相渗曲线校正方法
Göğüş et al. Broad-crested weirs with rectangular compound cross sections
CN107515025B (zh) 古城墙病害等级评估方法
Bian et al. Evaluating the effect of soil structure on the ground response during shield tunnelling in Shanghai soft clay
CN111811468B (zh) 高心墙堆石坝沉降变形的预测方法
CN113806852B (zh) 一种深水薄壁钢圆筒稳定性预测方法
CN105155425A (zh) 满堂支架流体仿真模拟预压施工方法
Shi et al. A statistical model of deformation during the construction of a concrete face rockfill dam
CN110969370B (zh) 建筑结构构件的质量风险分析方法
CN104111137B (zh) 基于三次多项式模型的锚杆承载力预测方法
CN113987663A (zh) 基于实测数据考虑计算参数变化的软土路基沉降预测方法
CN111950065A (zh) 高心墙堆石坝施工期水平变形的预测方法
CN110287634A (zh) 一种基于体积力施加的坝肩变形模拟方法
CN110197015A (zh) 一种坝基预应力锚索有效拉应力测定方法
CN114662341A (zh) 岩体临界滑动面极限分析方法
Urzua et al. New approach to analysis of consolidation data at early times
CN115221595A (zh) 一种基于复合地基的碎石粒料桩承载力计算方法
CN109933861B (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)
ERGUN et al. COMPARISON OF GEOMETRICAL IMPERFECTION DEFINITIONS ON ENCASED COMPOSITE COLUMNS
Mozer et al. Time-dependent behavior of concrete beams
RU2457455C2 (ru) Способ определения герметичности запорных арматур нефтепровода
Wessels et al. Flow-gauging structures in South African rivers part 2: calibration
CN110020792A (zh) 基于洪水洪峰-洪量联合的高堆石坝施工度汛风险预测方法

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