利用三维有限元对土石坝施工期位移及应力进行分析.docVIP

利用三维有限元对土石坝施工期位移及应力进行分析 　　摘要：本文主要通过采用三维有限元对大坝坝体进行建模，并模拟坝体施工期填筑过程，根据竣工期及蓄水期的位移及应力成果可知土石坝在竣工期和蓄水期两种工况下的垂直位移等值线、水平位移等值线、和大小主应力等值线，结果合理，进一步验证了程序的可行性。 　　关键词：土石坝，三维有限元，施工期，位移及应力分析 　　1 前言： 　　坝是截断河流或溪谷，用以拦截江河、形成水库或壅高水位的挡水建筑物。挡水坝按筑坝材料和坝体型式可分为土石坝、拱坝、重力坝三大类。人类自从开始农业生产以来，就建造了各种类型的土石坝。土石坝是指由土、石料等当地材料建成的坝，是历史最为悠久的一种坝型。以后，由于硅酸盐水泥的发明，各国建造混凝土坝成为一种时髦，曾使土石坝的建造速度放慢，在坝中占的比例有所下降。20世纪50年代后期，土石坝又发挥其生命力，许多国家又开始注重于修建土石坝。目前，土石坝是世界坝工建设中应用最为广泛和发展最快的一种坝型。截止到1997年底，我国共建各类水库8.5万余座，居世界第一位，其中90%以上水库的坝型是土石坝。这些水库大坝在防洪、灌溉、发电、城镇供水和水产养殖等方面都发挥了巨大效益。 　　土石坝的迅速发展，势必会给人们提出这样的问题：造坝方法是否合理，坝体是否安全。土石坝虽有最悠久的历史，但直到20世纪中叶，工程设计中人们关心和研究的问题仍然只集中在坝体、坝基渗流控制以及坝坡稳定等方面；对坝体的应力和变形则只做些近似估算，例如按单向压缩分层总和法预估沉降等。直到20世纪60年代初，由于大型电子计算机的发明，才借助于有限元技术，对坝体进行系统的分析。土石坝的应力和变形计算有过多种方法，主要有弹性理论方法，极限平衡理论方法，弹塑性分区方法，弹塑性差分方法，碎块体理论，有限单元法等，其中有限单元法是最准确、最有效的方法。 　　2 实例分析 　　黄壁庄水库位于河北省石家庄市西北约30公里处，是海河流域子牙河系两大支流之一──滹沱河中下游重要的、控制性的大（1）型水利枢纽工程。水库于1958年动工修建，1959年拦洪，1960年蓄水，经受了1963年特大洪水后，于1965年进行扩建，至1968年达到现状规模。水库正常蓄水位120.00米，设计洪水位125.84米，校核洪水位128.00米，它与上游28km处的岗南水库一起控制流域面积23400km，总库容12.1亿立方米。水库主要建筑物有：主坝、副坝、重力坝、正常溢洪道、非常溢洪道、灵正渠涵管及电站等。主坝为水中填土均质坝，坝顶长1843m，最大坝高30.7m，坝顶高程128.7m。副坝亦为水中填土均质坝（左右两端部分为碾压式均质坝），坝顶长6907.3m，最大坝高19.2m，坝顶高程129.2m。坝基由第四系各种成因类型堆积物构成，最厚达70米，一般为4层结构，由上到下为壤土、砂层、卵砾石层和基岩。其中，基岩为硅质灰岩、大理岩与千枚岩互层，基岩风化强烈，溶蚀及结构带发育，工程、水文、地质条件非常复杂。 　　考虑土坝施工过程中荷载随坝体上升逐级施加的特点，为反映施工过程中各阶段的应力和变形情况，采用逐级加荷的增量法。施工加荷过程只计坝料自重，认为填土沉陷一次完成，未考虑湿陷。 　　2.1三维有限元的建立和计算条件 　　2.1.1有限元模型的建立 　　1实体模型的建立 　　坝体从上游至下游依次为：贴坡、坝体土、沥青混凝土防渗墙；坝基从上到下依次为：壤土层、砂层、卵石层、基岩。计算区域的选择：坝体顶宽15米，长度为226米。采用笛卡儿直角坐标系。以垂直于坝轴线的水平方向为x坐标轴，以过上游迎水面坡脚的铅垂线为y轴，以平行于坝轴线的方向为z轴，计算坐标原点选取在防渗墙轴线坝体底面的交汇处。 　　2有限元模型的建立 　　在计算中网格划分采用plane82单元，它是2维4节点单元（PLANE42）的高阶版本plane82单元；按照施工顺序，对面进行自网格划分，单元的边长为l。在对面进行自由网格划分之前先对线进行分段，这样总体划分出来单元会比较整齐。整个断面总共划分607个单元，647个节点单元剖分结果如图2-1所示。 　　同时，采用六面体单元对该模型进行划分，坝体共划分5313个节点，4540个单元；坝基共划分9009个节点，7600个单元。 　　图2-1 网格剖分图 图4-2 坝体模型 　　2.1.2 计算材料参数选取 　　本次计算中，共涉及到7种材料，由于材料参数缺乏，其中坝体部分，即坝体土、采用广义邓肯-张本构模型，沥青混凝土心墙、壤土层、砂层、卵石层、基岩采用Drucker-Prager本构模型，材料参数如表2-1所示。 　　表2-1 坝体坝基力学参数指标 　　材料号 　　材料名称 　　E 　　/（KPa