Page 136 - 2025年第56卷第2期

［ 15］ 郭盼盼，龚晓南，汪亦显 . 考虑土与结构非线性接触特性的格形地下连续墙围护结构力学性状研究［J］ 岩土
.
工程学报，2021，43（7）：1201-1209.
［ 16］ 翁永红，张练，徐唐锦，等 . 高水头下大型导流洞新型堵头-围岩相互作用规律与安全评价［J］ 岩土力学，
.
2020，41（1）：242-252.
［ 17］ 何坤，龙军飞 . 毛尔盖水电站导流洞封堵施工组织及温控设计［J］ 水电站设计，2013，29（4）：4-7，18.
.
.
［ 18］ 杨光伟 . 二滩水电站导流洞封堵施工［J］ 水电站设计，2001，17（3）：29-33.
［ 19］ 杨静安，伍鹤皋，苏凯，等 . 大型水电站导流洞堵头结构与稳定分析［J］ 水电能源科学，2007，25（1）：
.
94-97.
.
［ 20］ 朱伯芳 . 大体积混凝土温度应力与温度控制［M］ 2 版 . 北京：中国电力出版社，2012.
［ 21］ 杨静安，郭红彦 . 拉西瓦水电站导流洞堵头结构与稳定分析［J］ 电网与清洁能源，2008，24（11）：71-75.
.

Study on the temperature control simulation and stability safety of diversion tunnel plug
based on contact model

1
1
2
CHENG Jing ，XIAO Wenjun ，ZHOU Wei ，ZHAO Lanhao 1

（1. College of Water Conservancy and Hydropower Engineering，Hohai University，Nanjing 210098，China；

2. School of Water Resources and Hydropower Engineering，Wuhan University，Wuhan 430072，China）

Abstract：To address the potential sliding problem between the concrete plug and surrounding rocks of the diversion
tunnel，a mixed finite element method is introduced to establish a full process thermal stress field simulation model

that takes into account the contact nonlinearity，with the HZY Hydropower Station as an example. Based on the conse⁃

quent calculation results，the sliding-stability safety factor of plug is checked by the finite element limit equilibrium

method. In order to verify the rationality of the proposed approach，the rigid limit equilibrium method and the elasto⁃
plastic finite element method are adopted for calculating the safety factor respectively for comparison. The results

show that，the safety factor of anti-sliding stability of the plug under the three stability safety review methods can all

meet the specification requirements. However，the thermal stress simulation analysis model considering the contact
nonlinearity can better reflect the continuous-discontinuous process between plug and surrounding rocks during con⁃
struction and operation period，thus it can provide a more accurate evaluation for the contact stress state and plug sta⁃

bility under various temperature control schemes and water impoundment schemes. To improve the calculation effi⁃
ciency of thermal stress problems with contact nonlinearity，a simplified method is given for the composition of matrix

equation，in light of the character of contact iteration function. The research here can be referred for similar thermal
stress problems with contact nonlinearity in civil and hydraulic projects.

Keywords：diversion tunnel plug；anti-sliding stability；mixed finite element method；nonlinear contact；tempera⁃

ture control and crack prevention
（责任编辑：韩 昆）

— 276 —