混凝土重力坝动力特性仿真计算方法研究王艳.pdfVIP

第 32 卷第 3 期 　　　　　　　　　　　　　　　　 Vol. 32 , No. 3 　　人 　民　黄 　河 　　　　　　　　　　　　　　　　　 　 　　20 10 年 3 月 　　　　　　　　　　　　　　　　　　　YELLOW 　R IV ER 　　　　　　　　　　　　　　　　　　 　M ar. , 20 10　　 【水利水电工程 】 混凝土重力坝动力特性仿真计算方法研究 王 　艳 ,辛全才 ,刘 　艳 ,崔 　娜 (西北农林科技大学 水利与建筑工程学院 ,陕西 杨凌 712 100) 摘 　要 :在 8度设防烈度下 ,应用大型有限元软件 AN SYS建立某混凝土重力坝的三维有限元模型 ,采用峰值调整后的天 津地震波进行了动力响应时程分析 ,分别采用固定连接和接触算法考虑坝体与地基的相互作用关系 。计算结果表明:在 坝体与地基两种相互作用计算模型中 ,数值仿真结果的基本规律是一致的 ,但接触计算模型中大坝的动力响应显著降 低 ,实际计算时应将两者结合使用 。 关 　键 　词 : 混凝土重力坝 ; 地震载荷 ; 动力分析 ; 有限元仿真 中图分类号 : TV 642　　　文献标识码 : A 　　　do i: 10. 3969 / j. issn. 1000 - 1379. 20 10. 03. 064 　　混凝土重力坝是水利工程中常见的一种大体积挡水建筑 坡面垂直 ,下游坡面系数为 0. 75。按照规范 ,混凝土动态强度 物 ,在其工作周期内常受到地震荷载的作用 。从国内外已有的 和动态弹性模量的标准值可较其静态标准值提高 30% ,混凝土 经受过强震作用的混凝土重力式高坝的震害实例来看 ,大坝虽 动态抗拉强度的标准值可取动态抗压强度标准值的 8% 。坝体 未出现整体失稳或毁坏的情况 ,但却在上部出现了上下游贯穿 弹性模量为 36. 4 kN /mm2 ,泊松比为 0. 167 ,容重为 24 kN /m3 ; 性的裂缝 [ 1 ] ,说明高重力坝存在抗震薄弱部位 。 岩基的弹性模量为 29. 0 kN /mm2 ,泊松比为 0. 3。 计算机技术的迅猛发展 ,使数值仿真计算在大坝抗震安全 2. 2　计算模型及荷载 设计 [2 ] 中起到了越来越重要的作用 ,选择合理的计算模型与仿 地基沿上下游及深度方向各取 1. 5倍坝高范围进行计算 , 真方法 ,对提高有限元计算结果的有效性至关重要 。而在大坝 坝体及地基采用空间 8节点六面体单元剖分 ,地基周边取法向 承受地震荷载的动态响应计算中 ,坝体与周围岩石及地基的连 约束 。按最不利原则组合静态 、坝体 自重 、上下游静水压力及 接方式对数值计算结果有很大的影响 。笔者利用大型有限元 坝基扬压力和动态反应得到综合反应 。地震选用峰值调整后 通用软件 AN SYS建立仿真模型 ,对某大坝与基础分别采用固 的天津地震波 , 由地基底部均匀输入 ,只考虑地基的弹性作用 , 定连接和接触算法进行了计算 ,并比较和分析了计算结果 。 消除坝基对地震的放大作用 ,取垂直方向和顺河向的记录 ,计 算中按照 0. 1 s的间隔取值 ,记录时长 19. 11 s。有限元计算模 1　固定连接与接触的定义 型见图 1。 在有限元中 ,固定连接与接触是用来描述相邻两个结构相 互关系的两种方法 。固定连接是指两个相邻结构的单元之间 共用节点 ,相互之间没有相对位移 ,应力和变形均完整地传递 , 相邻单元的应变是连续的。接触是指相邻结构的单元不共用 节点 ,两者之间有相对位移 , 比如错动 、摩擦等 。在接触关系 中 ,应力和变形不完全传