1、第 18 卷增刊 2地 下 空 间 与 工 程 学 报Vol.182022 年 12 月Chinese Journal of Underground Space and EngineeringDec.2022高切坡挡墙整体稳定性和工作性能的数值分析唐晓松1,2,郑颖人3,陈佳4,董世勇4(1.重庆地质矿产研究院 重庆市地质灾害自动化监测工程技术研究中心,重庆 400042;2.重庆公共运输职业学院 铁道与建筑学院,重庆 402247;3.重庆市地质灾害防治工程研究中心,重庆 400041;4.重庆市三峡库区工作服务中心,重庆 401147)摘要:高切坡支挡工程的健康运营对三峡库区移民迁建区人民
2、的安居乐业和社会的和谐稳定至关重要,目前关于服役期混凝土材料高切坡挡墙运营安全的分析缺少基于力学机制的定量分析方法,为此文章基于数值极限分析方法开展了服役期高切坡挡墙整体稳定性和工作性能的定量分析。首先,采用有限元强度折减法进行了混凝土材料强度衰减过程中高切坡挡墙工程整体稳定性的动态评价。然后,基于数值极限应变法建立的混凝土材料破坏标准,通过数值计算结果直观显示服役期内挡墙结构的易损关键部位和损伤特征。服役期高切坡挡墙整体稳定性和工作性能定量分析方法的提出,可以掌握安全管理的主动权,为其安全评价和监测预警,以及养护管理决策等提供了新思路和新方法。关键词:高切坡挡墙;服役期;整体稳定性;工作性能
3、;数值分析中图分类号:TU432文献标识码:A文章编号:1673-0836(2022)增 2-0910-06Numerical Analysis on Overall Stability and Working Performance of High Cut Slope Retaining WallTang Xiaosong1,2,Zheng Yingren3,Chen Jia4,Dong Shiyong4(1.Chongqing Institute of Geology and Mineral Resources,Chongqing Engineering Research Center o
4、f Automatic Monitoring for Geological Hazards,Chongqing 400042,P.R.China;2.Chongqing Vocational College of Public Transportation,College of Railway and Architecture,Chongqing 402247,P.R.China;3.Engineering and Technology Research Center of Geological Hazard Prevention and Treatment,Chongqing 400041,
5、P.R.China;4.Chongqing Water Resources Bureau,Chongqing Three Gorges Reservoir Area Work Service Center,Chongqing 401147,P.R.China)Abstract:The healthy operation of the high cut slope retaining project is crucial to the peoples living and working in the resettlement and construction area of the Three
6、 Gorges Reservoir area,as well as the harmony and stability of the society.At present,the analysis of the operation safety of the high cut slope retaining wall of concrete materials during service period lacks quantitative analysis method based on mechanical mechanism.Therefore,based on the numerica
7、l limit analysis method,quantitative analysis of overall stability and working performance of high cut slope retaining wall in service period is carried out in the paper.First,the finite element strength reduction method is used to dynamically evaluate the overall stability of high cut slope retaini
8、ng wall in the process of concrete material strength attenuation.Then,the failure standard of concrete material is established based on the numerical limit strain method.The key vulnerable parts and damage characteristics of retaining wall structure in service period are displayed directly through t
9、he numerical calculation results.The quantitative analysis method of 收稿日期:2022-05-03(修改稿)作者简介:唐晓松(1979),男,江苏海门人,博士,正高级工程师,主要从事岩土工程稳定分析及其数值模拟方面的工作。E-mail:hnm97 基金项目:重 庆 市 教 委 科 学 技 术 研 究 项 目(KJZD-K201905801);重 庆 市 技 术 创 新 与 应 用 发 展 专 项 重 点 项 目(cstc2019jscx-tjsb0019)overall stability and working perfo
10、rmance of high cut slope retaining wall in service period can master the initiative of safety management,and provide a new idea and method for safety evaluation,monitoring and early warning,and maintenance management decision.Keywords:high cut slope retaining wall;service period;overall stability;wo
11、rking performance;numerical analysis0引言重庆三峡库区高切坡是三峡库区移民迁建中形成的,广泛分布于三峡库区移民迁建城镇中,并大都紧邻企业、机关、学校、医院、住宅以及生命线工程等,高切坡支挡工程的健康运营对移民迁建区人民的安居乐业和社会的和谐稳定至关重要。在三峡重庆库区高切坡防护工程中,混凝土材料挡墙应用十分广泛,且结构形式多样,有重力式、悬臂式和衡重式等多种类型1-2。由于混凝土材料在荷载和环境等因素的共同作用下,其力学性能会随着服役年限的增长逐渐衰减3-5。尽管通过监测预警可以有效防范潜在的安全隐患6-7,但是仍亟需开展服役期高切坡挡墙整体稳定性分析方法的
12、研究,掌握安全管理的主动权。目前,针对服役期高切坡挡墙整体稳定性的分析大都采用层次分析法、神经网络或模糊逻辑等定性或半定性的方法8-10,以建立结构安全等级评价指标体系和分级标准为主要目的,缺少基于力学机制的定量分析方法。为此,本文应用数值极限分析方法进行了服役期混凝土材料高切坡挡墙的整体稳定性和工作性能分析,其中整体稳定性分析采用有限元强度折减法,工作性能分析则是依据数值极限应变法得到的混凝土材料的破坏标准,这两种方法均属于数值极限分析方法。基于数值极限分析方法开展服役期高切坡挡墙的整体稳定性和工作性能的定量分析,为其安全评价和监测预警,以及养护管理决策等提供了新思路和新方法。1数值极限分析
13、方法1.1有限元强度折减法有限元强度折减法在不断降低岩土类材料力学强度参数(黏聚力 c 和内摩擦角)的同时,进行弹塑性数值计算,根据数值计算收敛与否判断防护工程是否达到破坏状态。当数值计算不收敛时,最后一步计算收敛的荷载增量步所对应的强度折减系数,即为强度储备安全系数。同时,根据计算结果可自动生成破坏面。有限元强度折减法既具有数值方法适应性广的优点,又具有极限分析法贴近岩土工程设计且实用的优点,因此在边(滑)坡工程稳定性分析中被广泛应用11-14。1.2有限元极限应变法阿比尔的等15、李世贵等16提出弹塑性材料从弹性阶段进入屈服阶段,然后经过塑性发展至破坏状态,此时材料剪应变达到极限剪应变。通
14、过单轴压缩试验的数值模拟,可以得到极限荷载下材料极限剪应变的数值解,并进行了室内试验验证。这种采用数值分析方法求解岩土类材料极限剪应变,并据此进行破坏状态判断的方法即为数值极限应变法。数值分析时不同有限元软件采用的剪应变假设各自不同,如 ANSYS 采用等效塑性应变,FLAC采用应变偏张量第二不变量,PLAXIS 则采用总偏应变,但这并不影响数值极限应变法的使用,因为在高切坡挡墙整体稳定性和工作性能分析时,某一类材料的极限应变是在同一力学参数条件下采用同一软件计算得到的。2计算参数的确定 2.1混凝土材料抗剪强度参数的确定根据混凝土轴心抗压强度 fc和轴心抗拉强度ft可以得到两个莫尔应力圆,对
15、两个莫尔圆作公切线,基于莫尔 库仑强度理论可以得到混凝土材料抗剪强度参数黏聚力 c 和内摩擦角 为:c=fcft2(1)tan=fc-ft2fcft(2)采用直线代替混凝土材料实际的破坏包络线尽管存在一定误差17-18,但混凝土材料抗剪强度参数的确定不是本文的主要研究内容,实际工程中可通过小直径芯样的剪切试验等方法准确获取。1192022 年增刊 2唐晓松,等:高切坡挡墙整体稳定性和工作性能的数值分析文中数值分析时采用的不同强度等级混凝土材料的抗剪强度参数如表 1 所示。表 1不同强度等级混凝土材料抗剪强度参数Table 1Shear strength parameters of concre
16、te materials with different strength grades抗剪强度参数混凝土强度等级C15C20C25C30C35C40轴心抗压强度 fc/MPa7.29.611.914.316.721.1轴心抗拉强度 ft/MPa0.911.101.271.431.571.71黏聚力 c/MPa1.281.621.942.262.563.00内摩擦角/()50.86 52.60 53.82 54.90 55.91 58.222.2混凝土材料极限应变的确定采用 PLAXIS 有限元数值程序,参考阿比尔的等15的方法计算混凝土材料的极限剪应变,数值模型的尺寸为 150 mm150 mm150 mm 的立方体,底面施加约束,顶面施加竖直向下的均布荷载,不考虑顶面和底面的摩擦。立方体试块有限元网格如图 1 所示。采用表 1 中的混凝土材料抗剪强度参数。图 1立方体试块有限元网格示意图Fig.1Schematic diagram of FEM grid of cubic specimen以 C35 混凝土材料极限剪应变的数值计算为例,取图 1 中顶面四个角点中任一点处的一个四面体单