CN115544822A - 一种确定高填方库盆现场压实土石混合料流变参数的方法 - Google Patents
一种确定高填方库盆现场压实土石混合料流变参数的方法 Download PDFInfo
- Publication number
- CN115544822A CN115544822A CN202211015714.0A CN202211015714A CN115544822A CN 115544822 A CN115544822 A CN 115544822A CN 202211015714 A CN202211015714 A CN 202211015714A CN 115544822 A CN115544822 A CN 115544822A
- Authority
- CN
- China
- Prior art keywords
- test
- reservoir basin
- site
- rock mixture
- parameters
- 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.)
- Pending
Links
- 239000011435 rock Substances 0.000 title claims abstract description 42
- 239000000203 mixture Substances 0.000 title claims abstract description 39
- 238000000034 method Methods 0.000 title claims abstract description 18
- 239000000463 material Substances 0.000 title claims description 15
- 238000012360 testing method Methods 0.000 claims abstract description 74
- 238000006073 displacement reaction Methods 0.000 claims abstract description 39
- 238000010276 construction Methods 0.000 claims abstract description 22
- 238000012544 monitoring process Methods 0.000 claims abstract description 13
- 230000003068 static effect Effects 0.000 claims abstract description 9
- 230000002068 genetic effect Effects 0.000 claims description 10
- 238000005192 partition Methods 0.000 claims description 7
- 238000004458 analytical method Methods 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000003860 storage Methods 0.000 abstract description 8
- 238000005096 rolling process Methods 0.000 abstract description 7
- 238000000518 rheometry Methods 0.000 abstract description 5
- 238000005457 optimization Methods 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 5
- 230000007774 longterm Effects 0.000 description 4
- 238000005056 compaction Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005429 filling process Methods 0.000 description 2
- 238000011065 in-situ storage Methods 0.000 description 2
- 230000002411 adverse Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000011798 excavation material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000003892 spreading Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/06—Multi-objective optimisation, e.g. Pareto optimisation using simulated annealing [SA], ant colony algorithms or genetic algorithms [GA]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/20—Hydro energy
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开了一种确定高填方库盆现场压实土石混合料流变参数的方法，涉及技术领域，具体为一种确定现场压实土石混合料流变参数的方法，包括以下步骤：S1、在现场碾压结束后的土石混合料库盆表面开展平板载荷试验和直剪试验，获取加载‑沉降和剪应力‑剪切位移数据；S2、反演确定现场压实土石混合料的静力本构模型参数，并计算出库盆填筑产生的瞬时位移；S3、结合施工期位移监测资料，分离出流变产生的位移，进而反演确定流变参数。本发明避免了利用运行期位移监测资料反演出的流变参数难以在填筑施工阶段应用于库盆变形控制方案优化设计的问题，通过施工期监测资料获取的流变参数更具工程指导意义。
Description
技术领域
本发明涉及土石混合料流变参数确定技术领域，具体为一种确定高填方库盆现场压实土石混合料流变参数的方法。
背景技术
近年来，我国抽水蓄能电站建设进入了高速发展阶段，为抬高上库水头并实现电站内部的土方平衡，不少电站的上水库库盆填筑料均采用下水库的土石混合开挖料。为控制土石混合料高填方的沉降变形，工程中通常采用薄层重碾+激振的碾压方式。尽管如此，由于土石混合料具有较为显著的流变特性，其产生的长期变形仍会对运行期的库盆防渗系统安全运行造成不利影响，因此，若能在施工阶段准确地预测库盆长期变形，对于变形控制措施的提前设置以及竣工后预沉降时间的确定、保障库盆防渗系统安全运行至关重要。
目前，确定填筑料流变参数主要包括室内流变试验和反演分析两种方法。对于室内试验，由于试样尺寸和试验时间上均存在缩尺效应，无法反映现场真实填筑料的流变特性，所得参数难以准确预测库盆长期变形。对于反演分析方法，通常情况下采用运行期的位移监测资料进行流变参数反演，但反演结果只能对竣工后的库盆运行安全性进行评价，无法指导填筑施工阶段的变形控制方案实施。此外，对于施工期监测到的坝体位移，通常包括填筑引起的瞬时位移和流变引起的位移两部分，反演时难以将两者准确分开，这使得反演得到的流变参数不具唯一性，难以合理预测库盆长期变形。
发明内容
针对现有技术的不足，本发明提供了一种确定高填方库盆现场压实土石混合料流变参数的方法，解决了上述背景技术中提出的问题。
为实现以上目的，本发明通过以下技术方案予以实现：一种确定高填方库盆现场压实土石混合料流变参数的方法，包括以下步骤：
S1、当库盆填筑至一定高度时，在不同材料分区的库盆表面开展平板载荷试验，获取加载-沉降（p-s）曲线的数据；
S2、在平板载荷试验位置，同步开展大型直剪试验，获取剪应力-剪切位移曲线的数据，并确定压实土石混合料的内摩擦角和粘聚力；
S3、选用邓肯E-B模型作为压实土石混合料的静力模型，建立平板载荷试验和直剪试验的有限元模型，利用平板载荷试验和直剪试验结果，基于改进加速遗传算法反演压实土石混合料E-B模型参数；
S4、根据库盆的实际填筑施工分级以及材料分区，建立库盆三维有限元模型；
S5、选用南科院五参数模型作为流变模型，利用施工期的库盆位移监测资料，采用改进加速遗传算法反演流变参数，以计算沉降量与实测沉降量的均方根误差最小为目标函数，得到一组最优的计算参数，即为反演结果。
可选的，所述步骤S1、步骤S2中，在开展平板载荷试验和直剪试验时，对应的库盆填筑高度不小于20m。
可选的，所述步骤S3中，同时采用平板载荷试验和直剪试验结果反演E-B模型参数，反演分析时，对平板载荷试验和直剪试验两套有限元网格进行交替计算，目标函数为：
其中，j为平板载荷试验加载级数的编号，n为平板载荷试验总共加载级数，
为第j级荷载下实测沉降，
为对应的计算值；k为直剪试验中剪切力的编号，m为直剪试验剪切力的数量，
为剪切力k对应的实测剪切位移，
为对应的计算值。
可选的，所述步骤S5中的施工期库盆位移测点起测时间至少6个月，库盆填筑产生的瞬时变形由步骤S3确定的静力本构模型参数计算。
本发明提供了一种确定高填方库盆现场压实土石混合料流变参数的方法，具备以下有益效果：
针对现场压实土石混合料，本发明提出基于原位载荷试验、直剪试验以及填筑施工期位移监测资料的现场压实土石混合料流变参数反演方法，首先通过平板载荷试验和直剪试验得到竖向荷载-沉降曲线、剪切应力-水平位移曲线，建立平板载荷试验和直剪试验的有限元模型，采用改进加速遗传算法反演现场压实土石填筑料的静力模型参数，然后结合施工期位移监测资料，建立考虑实际填筑过程的库盆三维有限元模型，利用改进加速遗传算法反演得到压实土石混合料的流变参数，该方法可准确分离出施工期产生的瞬时变形和随时间的流变变形，从而确保反演结果的准确性、可靠性和合理性。
本发明有效克服了利用运行期位移监测资料反演的流变参数无法指导施工期变形控制方案实施的问题；同时，克服了施工期位移中填筑荷载和流变产生的位移难以准确分离的问题，不仅确保了反演结果的准确性和合理性，还使反演结果更具工程指导意义。
附图说明
图1为本发明的载荷试验示意图；
图2为本发明的p-s曲线示意图；
图3为本发明的直剪试验示意图；
图4为本发明的剪应力-剪切位移曲线示意图；
图5为本发明的参数反演分析示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。
实施案例一
请参阅图1至图5，本发明提供一种技术方案：一种确定高填方库盆现场压实土石混合料流变参数的方法，包括以下步骤：
S1、当库盆填筑至一定高度时，在不同材料分区的库盆表面开展平板载荷试验，获取加载-沉降（p-s）曲线的数据；
S2、在平板载荷试验位置，同步开展大型直剪试验，获取剪应力-剪切位移曲线的数据，并确定压实土石混合料的内摩擦角和粘聚力；
其中，步骤S1、步骤S2中，在开展平板载荷试验和直剪试验时，对应的库盆填筑高度不小于20m；
S3、选用邓肯E-B模型作为压实土石混合料的静力模型，建立平板载荷试验和直剪试验的有限元模型，利用平板载荷试验和直剪试验结果，基于改进加速遗传算法反演压实土石混合料E-B模型参数；
采用平板载荷试验和直剪试验结果反演E-B模型参数，反演分析时，对平板载荷试验和直剪试验两套有限元网格进行交替计算，目标函数为：
其中，j为平板载荷试验加载级数的编号，n为平板载荷试验总共加载级数，
为第j级荷载下实测沉降，
为对应的计算值；k为直剪试验中剪切力的编号，m为直剪试验剪切力的数量，
为剪切力k对应的实测剪切位移，
为对应的计算值；
S4、根据库盆的实际施工分级以及材料分区，建立库盆三维有限元模型；
S5、选用南科院五参数模型作为流变模型，利用施工期的库盆位移监测资料，采用改进加速遗传算法反演流变参数，以计算沉降量与实测沉降量的均方根误差最小为目标函数，得到一组最优的计算参数，即为反演结果；
施工期库盆位移测点起测时间至少6个月，库盆填筑产生的瞬时变形由步骤S3确定的静力本构模型参数计算。
实施案例二
如图1-5所示，一种确定现场压实土石混合料流变参数的方法，包括以下步骤：
步骤一：预先埋设4个1.2´1.2m的大型剪切框，然后按照正常库盆填筑流程进行铺料、碾压，本实施案例中采用26吨的大功率振动碾对土石混合料进行碾压，碾压层厚为60cm，行车速度为2~3km/h，碾压遍数为8遍，激振力为410kN；
步骤二：在压实后的土石混合料库盆表面开展平板载荷试验，获得不同荷载下压实土石混合料的沉降位移；
步骤三：在压实后的土石混合料库盆表面开展大型直剪试验，获得不同竖向荷载下压实土石混合料的剪应力-剪切位移曲线，并利用摩尔-库伦强度准则确定现场压实土石混合料的内摩擦角和粘聚力；
步骤四：取邓肯E-B模型作为材料静力本构模型，根据平板载荷试验获得的p-s曲线和直剪试验获得的剪应力-剪切位移曲线，建立目标函数，采用改进加速遗传算法对现场压实土石混合料的E-B模型参数K，Kb，n、m进行反演；
其中，j为平板载荷试验加载级数的编号，n为平板载荷试验总共加载级数，
为第j级荷载下实测沉降，
为对应的计算值；k为直剪试验中剪切力的编号，m为直剪试验剪切力的数量，
为剪切力k对应的实测剪切位移，
为对应的计算值。
步骤五：根据库盆的实际填筑分级和材料分区建立库盆三维有限元模型，利用确定的E-B模型参数计算填筑产生的瞬时位移，结合施工期的位移监测资料，分离出施工期各测点由流变产生的位移；
步骤六：选用南科院五参数流变模型，建立目标函数，采用改进加速遗传算法反演获得现场压实土石混合料的流变参数。
综上所述，本发明的基本思路是先利用平板载荷试验和现场直剪试验得到的p-s曲线以及剪应力-剪切位移曲线作为邓肯E-B模型参数反演的监测资料，反演得到现场压实土石填筑料的邓肯E-B模型参数，再利用施工期的位移监测资料，分离出流变产生的位移，进而反演得到现场压实土石混合料的流变参数。
以上，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，根据本发明的技术方案及其发明构思加以等同替换或改变，都应涵盖在本发明的保护范围之内。
Claims (4)
1.一种确定高填方库盆现场压实土石混合料流变参数的方法，其特征在于，包括以下步骤：
S1、当库盆填筑至一定高度时，在不同材料分区的库盆表面开展平板载荷试验，获取加载-沉降曲线的数据；
S2、在平板载荷试验位置，同步开展大型直剪试验，获取剪应力-剪切位移曲线的数据，并确定压实土石混合料的内摩擦角和粘聚力；
S3、选用邓肯E-B模型作为压实土石混合料的静力模型，建立平板载荷试验和直剪试验的有限元模型，利用平板载荷试验和直剪试验结果，基于改进加速遗传算法反演压实土石混合料E-B模型参数；
S4、根据库盆的实际填筑施工分级以及材料分区，建立库盆三维有限元模型；
S5、选用南科院五参数模型作为流变模型，利用施工期的库盆位移监测资料，采用改进加速遗传算法反演流变参数，以计算沉降量与实测沉降量的均方根误差最小为目标函数，得到一组最优的计算参数，即为反演结果。
2.根据权利要求1所述的一种确定高填方库盆现场压实土石混合料流变参数的方法，其特征在于，所述步骤S1、步骤S2中，在开展平板载荷试验和直剪试验时，对应的库盆填筑高度不小于20m。
3.根据权利要求1所述的一种确定高填方库盆现场压实土石混合料流变参数的方法，其特征在于，所述步骤S3中，采用平板载荷试验和直剪试验结果反演E-B模型参数，反演分析时，对平板载荷试验和直剪试验两套有限元网格进行交替计算，目标函数为：
其中，j为平板载荷试验加载级数的编号，n为平板载荷试验总共加载级数，
为第j级荷载下实测沉降，sj为对应的计算值；k为直剪试验中剪切力的编号，m为直剪试验剪切力的数量，
为剪切力k对应的实测剪切位移，dk为对应的计算值。
4.根据权利要求1所述的一种确定高填方库盆现场压实土石混合料流变参数的方法，其特征在于，所述步骤S5中的施工期库盆位移测点起测时间至少6个月，填筑产生的库盆瞬时变形由步骤S3确定的静力本构模型参数计算。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211015714.0A CN115544822A (zh) | 2022-08-23 | 2022-08-23 | 一种确定高填方库盆现场压实土石混合料流变参数的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211015714.0A CN115544822A (zh) | 2022-08-23 | 2022-08-23 | 一种确定高填方库盆现场压实土石混合料流变参数的方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN115544822A true CN115544822A (zh) | 2022-12-30 |
Family
ID=84725215
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211015714.0A Pending CN115544822A (zh) | 2022-08-23 | 2022-08-23 | 一种确定高填方库盆现场压实土石混合料流变参数的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115544822A (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117910285A (zh) * | 2024-03-20 | 2024-04-19 | 中国电建集团昆明勘测设计研究院有限公司 | 库盆设计曲面模型的创建方法、装置、设备及存储介质 |
-
2022
- 2022-08-23 CN CN202211015714.0A patent/CN115544822A/zh active Pending
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117910285A (zh) * | 2024-03-20 | 2024-04-19 | 中国电建集团昆明勘测设计研究院有限公司 | 库盆设计曲面模型的创建方法、装置、设备及存储介质 |
| CN117910285B (zh) * | 2024-03-20 | 2024-06-04 | 中国电建集团昆明勘测设计研究院有限公司 | 库盆设计曲面模型的创建方法、装置、设备及存储介质 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107748111B (zh) | 一种岩体结构面长期抗剪强度的确定方法 | |
| CN107543789B (zh) | 一种研究爆破动力失稳机理的边坡模型试验方法及系统 | |
| CN105259035B (zh) | 一种岩石材料时效-弹塑性力学本构模型的构建方法 | |
| CN107679293B (zh) | 堆石坝室内力学参数和碾压施工参数相互优选的方法 | |
| CN108444813A (zh) | 多尺度土石混合体-基岩界面剪切特性测试装置和方法 | |
| CN103915030B (zh) | 一种具扩缩功能的地下工程掘进模型试验装置 | |
| CN204694560U (zh) | 一种室内冻土模型桩静力加载实验装置 | |
| Li et al. | Effects of the soil water content and relative roughness on the shear strength of silt and steel plate interface | |
| CN106294984B (zh) | 一种基于细观裂纹增长率确定岩体损伤门槛的方法 | |
| Zhang et al. | Creep behavior of rocks with rough surfaces | |
| CN104328776A (zh) | 一种预测动力强夯对土体及周边环境影响的方法 | |
| CN110243681B (zh) | 一种粗粒料原位载荷试验确定颗粒破碎率的方法及应用 | |
| CN115544822A (zh) | 一种确定高填方库盆现场压实土石混合料流变参数的方法 | |
| Wang et al. | Study on acoustic emission characteristics of sandstone under different fracture modes | |
| Xie et al. | Numerical study on cracking behavior and fracture failure mechanism of fractured rocks under shear loading | |
| CN107391803A (zh) | 一种土石混填路基土的cbr数值试验方法 | |
| CN111896202A (zh) | 一种用于地震动力边坡模型试验相似材料配比的确定方法 | |
| Liu et al. | Movement and embedding characteristics of interlayer aggregates during roller-compacted concrete compaction process using discrete element simulation | |
| CN113008683B (zh) | 一种模拟采动非均布加载的实验装置及方法 | |
| CN206736984U (zh) | 现场地基沉降模拟装置 | |
| CN206521775U (zh) | 一种下压静载作用下的桩基础承载特性模型试验装置 | |
| CN115964901A (zh) | 基于离散单元法的水致岩石强度劣化的模拟方法及系统 | |
| Zhao et al. | The Load-Bearing Mechanism of Rock-Socketed Piles Considering Rock Fragmentation | |
| CN114112685B (zh) | 一种现场压实土石混合料前期固结应力确定方法 | |
| CN117405586B (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 | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20240401 Address after: No.1, No.2, Baiguang Road, Xicheng District, Beijing Applicant after: STATE GRID XINYUAN Co.,Ltd. Country or region after: China Applicant after: JIANGSU JURONG PUMPED STORAGE CO.,LTD. Applicant after: HOHAI University Address before: 605, Floor 6, Building 18, Luomashi Street, Xicheng District, Beijing Applicant before: STATE GRID XINYUAN Co.,Ltd. Country or region before: China Applicant before: JIANGSU JURONG PUMPED STORAGE CO.,LTD. Applicant before: HOHAI University Applicant before: TIBET AUTONOMOUS REGION WATER CONSERVANCY AND ELECTRIC POWER PLANNING & EXPLORATION DESIGN INSTITUTE |
|
| TA01 | Transfer of patent application right |