基于大坝安全监测资料的物理力学参数反演.pdfVIP

水 利 学 报 2004 年8 月 SHUILI XUEBAO 第8 期 文章编号: 0559-9350 (2004) 08-0098-05 基于大坝安全监测资料的物理力学参数反演 向衍 苏怀智 吴中如 (河海大学 水利水电学院 江苏 南京 210098) 摘要 本文介绍了采用实数编码遗传算法反演物理力学参数的基本步骤 然后结合通用的商业有限元软件MSC.Mar 并利用其与FORTRAN语言的接口 编写了基于MSC.Marc的物理力学参数的遗传算法反演程序 文中结合某重力坝溢 流坝段的监测资料 利用该方法反演了大坝坝体弹性模量和坝基变形模量 算例表明 采用遗传算法 MSC.Marc 反演大坝坝体及坝基的物理力学参数是可行的 关键词 遗传算法 大坝 反演 力学参数 中图分类号 TV698.1 文献标识码 A 仿效系统识别理论 将正分析成果作为依据 通过相应的理论分析 反求水工建筑物和地基的物理力 学参数等 即为物理力学参数的反演分析[1 2] 反演分析方法目前主要有优化法和逆解法 其中优化法应 用较普遍 随着人工神经网络 遗传算法等前沿学科不断地被引入反分析领域 进一步推动了反分析的发 展 与传统的参数反演算法相比 遗传算法是一种直接的随机寻优方法 其模拟了自然界生物进化过程中 的优胜劣汰和适者生存的法则 将选择 交叉和变异等概念引入到算法中 克服了传统优化方法的连续性 和可导性的限制以及易陷入局部最优解的缺点 搜索具有隐含并行性 可较快地搜索到全局最优解 并且 对目标函数的形态没有具体要求 与此同时 随着计算力学的发展 已开发了一批专业的商用有限元软件 如MSC.Marc[3 4] 自20世纪70年代初诞生至今 一直紧跟有限元方法的理论和计算机软 硬件的最新发展 已发展成为功能强大的有限元软件系统 它拥有丰富和完善的单元库 材料模型库和高效的求解器 接口 的开放性使用户很方便的扩展程序功能 本文利用MSC.Marc软件系统提供的与FORTRAN语言的接口 编写 了基于MSC.Marc的物理力学参数的遗传算法反演程序 1 遗传算法反演 美国的Holland教授及其同事D.E.Goldberg等提出了遗传算法的基本思想 该方法借鉴了物种进化的 思想 从代表问题可能潜在解集的一个种群(Population)开始 且每个种群由经过基因编码的一定数目的 个体组成 产生初始种群之后 按照适者生存和优胜劣汰的原理 逐代演化产生出越来越好的近似解 在 每代中 根据问题域中个体的适应度大小挑选个体 并借助于自然遗传学的遗传算子进行组合交叉和变异 产生出代表新的解集的种群 该过程将导致种群像自然进化一样的后生代种群比前代更加适应于环境 末 代种群中的最优个体经过解码 可以作为问题的近似最优解[5 6] 以待反演的大坝坝体及坝基部位的物理 力学参数作为反演变量 共n个 记作a=(a ,a , a ) 利用遗传算法的基本操作(即选择 交叉和变异) 1 2 n 模拟自然选择和自然遗传过程中的繁殖 交配和变异现象 从待反演参数组成的解种群中逐代产生新的群 体 比较个体 如此循环 最终搜索到最优个体 即为反演的参数 收稿日期: 2003-04-08 基金项目: 国家自然科学基金重点项目 国家973项目(2002CB412707) 作者简介: 向衍(1977-),男,土家族,湖北巴东人,博士生,主要从事大坝安全监控理论和方法研究. 1 水 利 学 报 2004 年8 月 SHUILI XUEBAO 第8 期 1.1 参数遗传反演的基本方法[7 8