基于模拟退火的重力坝最危险滑动面全局搜索.pdfVIP

第45卷 第 19期 人 民 长 江 V01．45．No．19 2014年 10月 Yangtze River Oct．， 2014 文章编号 ：1001—41792【014 19—0060一o3 基于模拟退火的重力坝最危险滑动面全局搜索 李 静 ，陈 健 云 ，徐 强 ，王 舒2 1．大连理工大学 海岸和近海工程国家重点实验室，辽宁 大连 116024； 2．大连理工大学 建设工程学部，辽 宁 大连 116024 摘要 ：为了弄清地震荷载作用下重力坝坝头部位的失稳规律，采用有限元非线性时程分析方法，得到重力坝在 地震作用过程 中的应力场 ，进而应用模拟退火算法进行重力坝临界滑动面的全局搜索。该方法可 以在所有增 量步的应力场 中搜索 出坝头区域的最脆弱滑动面，并得到与之对应的抗滑稳定系数及 出现时刻。分析结果表 明，重力坝在地震荷载作用下更容易向上游方 向失稳 。在一条特定的地震波作用下，即使调整了其幅值，坝头 部位最危险软弱滑动面 出现的时刻几乎仍是相 同的，且 同一时刻下滑动面形式非常相近 。 关 键 词 ：非线性有限元 ；模拟退火算法 ；滑动面；抗滑稳定系数 ；重力坝 中图法分类号 ：TV642 文献标志码 ：A 目前 ，有关学者主要是从强度和稳定两方面来评 1 非线性有限元时程分析 价混凝土重力坝 的安全性 ¨。当重力坝遭受地震这 一 特殊荷载时，整个坝体坝头部位的反应通常会很强 选用阿海水电站某一非溢流坝段作为工程 背景。 烈 ，因此，有必要校核该情况下坝头部位的稳定性。传 该坝段高90m，顶宽 12m，上游坝面为竖直的平面，下 统的稳定分析方法系采用刚体极限平衡理论，以安全 游坝面坡度为 1：0．75。模型采用平面有 限元剖分 。 系数作为工程安全度的评价标准。有限元法 的引入 ， 根据许多研究结果，重力坝在地震荷载作用下，塑性区 不仅可 以考虑结构变形 、应力重分布带来的影响，还能 主要分布在坝踵和坝头部位，而坝头部位的塑性区主 计人材料非线性等各种非线性 的影响，比刚体极 限平 要位于下游折坡点附近 ，故本文选择的搜索 区域为下 衡法前进了一步 。但受认识水平、工程经验等多方面 游折坡点附近的一片区域。 的限制 ，目前还缺乏与该计算方法相配套 的稳定度判 有限元分析 中考虑的荷载主要有坝体 自重 ，上游 别标准 。 静、动水压力及 、y方向的双 向地震荷载。根据设计 随着人工智能算法在水工结构领域中的应用 ，王 资料 ，上游正常水位为84m。地震波用的是 Koyna地 成华、夏绪勇等曾将非线性有限元分析和遗传搜索算 震波，地震持续时间 10s。计算模型采用刚性地基模 法相结合 ，研究了基于应力场的任意形状临界滑动面 型，材料非线性模 型采用的是塑性损伤模型 ，弹模 E 的遗传算法搜索技术 。李守巨、刘迎曦基于改进 的 3．1×10 Pa，泊松 比 0．2，结构阻尼采用 Ray． 模拟退火算法建立了边坡稳定分析最小安全系数全局 1eigh阻尼 。 搜索方法 。本文将非线性有 限元时程分析与模拟 根据规范要求 ，动力计算时应考虑结构和水体 退火算法结合起来 ，基于应力场来定义抗滑稳定安全 的动力相互作用 。地震过程 中库水动水压力作用采用 系数，进而应用模拟退火算法搜索重力坝在地震作用 附加质量来模拟，附加质量的计算采用 Westergaard公 下最脆弱滑动面及相应 的安全系数。 式 ] 收稿 日期 ：2014—05—17 基金项 目：国家重点基础研究发展 “973”计划课题 “极端条件下高碾压混凝土坝动力响应机制 ” 2013CB035905 ；国家 自然科 学基金资助项 目 2013CB035905 作者简介：李 静，女，副教授，主要从事结构抗震安全评估和风险分析研究。E—mail：lijing@dlut．edu．cn 第 19期 李 静，等：基于模拟退火 的重力坝最危险滑动面全局搜索 61 ÷ f1+6f2 ／2 1 22·目标函数 目标函数用于在模拟退火算法 中判断是否接受新 式 中，m 为坝体上游面第 i个节点上的附加质量 ；^为 解 ，结合安全系数 的定义 ，采用分段线性来模拟任意形 水库中水的深度 ；y为节点 i处到水面的距离；p 为水 状的滑动面 ，即用直线段连接 个离散点近似表示 的密度 ；b；和 b 则分别为节点 i两侧单元 的边长。 滑动面曲线。在应力场 已经确定 的条件下，优化 目标 2 模拟退火算法的基本原理 函数 即抗滑稳定安全系数 可以表示为 N 一1 模拟退火 simulatedannealing，SA 算法是一种通 ’ orfl,f +Iic k L ——一 3 用的启发式算法，是对局部搜索