基于重整化群的水工混凝土结构整体破坏概率研究.pdf

第 4期 2 0 1 0年 1 2月 水利水运工程学报 HYDRO． SCI ENCE AND ENGI NEERI NG N0 ． 4 De c ． 2 01 0 基于重整化群的水工混凝土结构整体破坏概率研究 顾 培英 ，黄勤红 ，邓 昌 ，汤 雷 1 ．南京水利科学研究 院 水 利部水科学与水工程重点实 验室 ， 江苏 南京2 1 0 0 2 9 ； 2 ．张家港市长江 防洪工程 管理处 ，江苏 张家港2 1 5 6 2 5 摘要 在借鉴重整化群理论和局部损伤概率的基础上， 形成水工混凝土结构整体破坏概率研究方法． 针对单元 损伤概率不等的正方体结构 ， 推导了 1 级原胞损 伤概率 ， 再依次 确定 级原胞的损伤概率 ， 从而得到结构整体 的 破坏概率 ． 传统 的岩石等概率破碎模 型只是针对单元损伤概率 相等情况 下 的一个 特例 ， 而该数 学模 型是 一个通 式． 以某大体积混凝土结构为例 ， 应用所建立 的数学模 型进 行了实例 分析． 计算结 果表 明， 当单 元损伤概 率不超 过 0 ． 7时 ， 结构较为安全 ， 否则结 构整体安全性会迅速 降低 ， 外界环境一旦恶化 ， 结构极有可能发 生垮 塌． 关 键 词 重整化群理论； 损伤概率； 原胞 ； 安全评价；整体破坏 中图分类号 T V 3 3 1 文献标志码 A 文章编号 1 10 9 6 4 0 X 2 0 1 0 0 4 0 0 0 1 0 5 随着科技 的飞速发展 ， 许多水工混凝土结构承担着生命线工程 的重任 ， 其失事所造成 的直接或间接经济 损 失将 不可估 量 ． 此外 ， 我 国很 多水 工混凝 土结构 已进 入 中老 年期 ， 由于地 质复 杂 、 勘 测设 计 不周 、 施工 缺 陷 、 管理不善、 结构老化等原因， 往往存在各种损伤与病害， 可能导致结构失稳或强度破坏 ， 严重影响工程正常运 行 ． 此外 ， 我 国 因地 震 活动频 繁 、 强度 大 、 震 源 浅 、 分 布 广 ， 是 一 个 震灾 严 重 的 国家 ． 四川 汶 川 大地 震 对 水 利工 程 破坏严 重 ， 灾后 重建 不仅 是恢 复 到地震 前 的状态 ， 而应 有更 长 远 的考 虑 ， 震 损后 结 构 要 保 持 足够 的安 全储 备 ， 避免 将来 遭遇 大地 震 时发生 垮塌 ． 所 以 ， 有 必要 开 展水 工 混凝 土 结 构整 体 安 全 性研 究 ． 目前 ， 混 凝 土 结构 损 伤状 况评价 是通 过外 观破 损调 查 ， 并借 助 于多种 无 损或局 部 有损 专项 检测 技术 进行 的． 如何 根据 局部 检测 结果 正 确评价 整体 安全 性 ， 目前 还是 个 难题 ． 老 化病 害 常 采 用 层 次 评 估 法 ⋯ ， 该 方 法 属 半 经 验半 理论 方 法 ． 根据 技术 发展 趋势 ， 老 化结 构评 估方 法将 会走 向基 于 可靠度 理论 的统计 分析 方法 ， 目前该 方法 已取 得初 步成 果 ， 但应用还主要停 留在构件可靠度水平． 国内外对水工混凝土结构抗震安全性进行 了较 多的研究 。 ， 但对损伤后 的整体安全性还缺乏充分 的 认识 ． 陈厚 群 院士提 出对 西部 所 有重 大高 坝工 程 的抗震 安全 性 重新进 行研 究 与评 估 ， 特别 要大 力开 展在最 大 可信地震作用下 ， 不发生溃坝严重灾变后果 的深入研究 ， 力求给出突破传统理念和方法 的更可靠和有共识 的 定量 判别 准则 ． 还有 不少 学者 开 展 了基 于 大 坝监 控 系 统 的安 全 评 价 ， 朱 伯 芳 院 士 等 提 出 在人 工 巡 视 和 仪器 监控之 外 ， 增加数 字 监控 ， 该方 法需 借 助部 分监 测仪 器 ． 本文 提 出基 于重 整化群 的混 凝 土结 构 整体 破 坏概率研究． 该方法可通过破损结构局部试验或统计数据 ， 得到特定荷载下 的局部损伤概率 ， 借鉴重整化群 方法从 局部 损伤 概率 出发 研究 整体 安全 性 ． 重整 化 群方 法是 通 过 改变 物 体 的粗 视 化 程度 来 研 究 观 测物 理 量 变化规 律 的方法 ， 可处 理各 种相 变 和临界 点 问题 ． 重整 化 群方 法成 功地 解 释 了逾 渗模 型 的突 变现 象 ， 包 括 岩石 破坏 、 岩石 多孔 介质 中的渗 流 问题 等 ， 由于混凝 土 与岩 石均 属脆 性材 料 ， 性 能较为 相似 ， 所 以将 重整化 群 方法 应用 于混凝 土 结构 ． 首 先介 绍传 统 的岩石 等概 率破 碎模 型 ， 由于传 统 岩石 等概 率破 碎模 型 目的是确定 破 碎 临界点 ， 对应 临界 点 的破 碎概 率也较 大 ， 并考 虑原 沿用 习惯 ， 所 以仍 称之 为 破碎 概 率 ． 而本 文 提 出的方 法 是 建 立水 工 混 凝 土结 收 稿 日期 2 0 0 1 00 31 6 基金项 目国家 自然科学基金资助项 目 5 0 7 7 9 0 3 6 ； “ 十一五 ” 国家科 技支撑计划项 目课题 2 0 0 9 B A K 5 6 B 0 4 作者简介 顾培英 1 9 6 8一 ， 女 ， 江苏南通人 ， 教授级高级工程师 ， 博士 ， 主要从 事结构健康诊断 、 安全评估 、 抗震及振 动控 制 等 的 研 究 ．E ma i l p y g u n h r i ． c n 学兔兔 w w w .x u e t u t u .c o m 2 水 利 水 运 工 程 学 报 2 0 1 0年 1 2月 构不等概率局部损伤与整体安全性之间的关系 ， 目的是得到结构整体安全性从渐变降低到突变的全过程 ， 为 了区别起见称之为局部损伤概率 ； 在推导公式时， 为了区分单元和原胞的损伤概率 ， 所以， 又分别称为单元损 伤概率和原胞损伤概率． 1 岩石等概 率破碎模 型 取 1级单元为最小单位 ， 按重整化群方法组成越来越大的原胞 ， 最低一级 由 8个单元组成 1级原胞． 按 照 A l l e g r e等人的做法 ， 假定原胞 中每个单元要么破碎 ， 要么完整． 根据每个单元破碎概率 P 。 确定 1级原 胞的破碎概率 P ． A l l e g r e假设如果原胞两面之间有一不破碎 的“ 支柱” 相连， 则认为原胞是完整的． 知道了单 元破 碎概 率后 ， 需求 出原 胞 的破碎 概率 ． 对 每个原 胞来说 ， 它包含 的破 碎单 元数 可从 0变 到 8个 ， 因此 可 能 的 情 况共有 2 。 2 5 6种 ． 去除 重复情 况 ， 从拓 扑学 角度共 有 2 2种不 同情 况 ， 岩石破 碎情 况示 意 图见 图 1 ． 图中 的 数字为破碎单元个数 ， a ， b ， ⋯⋯ ， f 为某破碎单元数不变情况下的不 同组合 ， 括号中的数字表示 2 5 6种情况中 出现 该情 况 的次数 ． 破碎 单元 用 角上 的黑点 表示 ． 4 f ， 5 c ， 6 b ， 6 c ， 7和 8的情况 是破 碎 的． 4 a 6 4 b 8 4 c 2 4 4 d 2 4 4 e 6 4 f 2 5 a 2 4 5 b 2 4 5 c 8 6 a 0 2 1 6 b 1 2 6 c 4 7 8 8 1 图 1 岩石破碎情况 Fi g ． 1 Ro c k b r e a k i ng o f pr i mi t i v e c e l l l e v e l 1 该模型假定 8 个单元的破碎概率相等， 用p 。 表示， 则 8 个单元均破碎的概率为P ， 7个单元破碎 1 个单 元 完整 的概 率是 P 1一 P 。 ， 依 次类 推 ． 考 虑 到每 种 情 况 出现 的次 数 ， 将 出现 的次数 乘 以各 自破碎 情 况 的概 率 ， 即为该情 况 的 1 级 原胞 破碎 概率 ． 最后 将 6种情 况下 的 概率 求 和 ， 1 级 原 胞破 碎 概 率 P 与 1级单 元 破 损 概率 P 。 之间存在以下关系 P P 8 p 1 一 P 0 1 6 p 1一 P 0 8 p 1 一P 0 2 p 1一 P 0 P 3 p 一8 p 4 p 2 1 重整化 后 ， 同样 的方 法可 用 于更 高 级数 ， 将式 1 中的下 标 1和 0替换成 n1和 ／ 7 , ， 有 P P 3 p 一8 p 4 p 2 2 如果 1级原胞特征长度为 2 ， 则第 n级原胞的特征长度为 2 ． 2 整体破坏概率数 学模 型的建立 岩石破碎模型是针对各个单元损伤概率均相等的情况 ， 对于实际结构而言， 单元损伤概率不等 ， 此 时式 2 不成立． 下面推导单元损伤概率不等时的结构整体破坏概率模型． 图2为三维实体模型． 图 2 a 为 3级原胞 ； 图 2 b 为 1 级原胞 ； 数字 1 8为 1级原胞中的 8个单元编 学兔兔 w w w .x u e t u t u .c o m 第 4期 顾 培英 ， 等 基于重整化群的水工混凝土结构整体破坏概率研究 3 号， 假定各单元损伤概率依次为P P ⋯⋯ ， P 本节就是根据每个单元损伤概率 ， 确定 1级原胞 的损伤概 率 p ． 参 照 图 1 ， 原假 设依 然成 立 ， 即如果 原胞 两 面之 间有一 不破 碎 的 “ 支 柱 ” 相连 ， 则认 为原 胞是 完整 的 a 3级原胞 b 1级原胞 图 2三维实体模型 Fi g． 2 Thr e e d i me n s i o n a l mo d e l 图 1中4 f ， 5 c ， 6 b ， 6 c ， 7和 8的情况是破碎的， 下面首先推导各情况下 的 1级原胞损伤概率 ， 分别用 p ， p 。， p ， p6 ，P ， P 表示． 推导得 以下公式 p 4 f P o zP 0 4P o P o7 1 一 P 。 1一 P 0 ] 1 一 P 0 6 1 一 P o P o ,P o 3P o6 P 。8 1 一 P 。 1 一 P o 4 1 一 P 05 1一 P o 7 3 p 。P o 3 P 0 4 P 0 l P 0 6 P 0 8 1一 P o 2 1一P o 5 1一 P o 7 P o , P o 2 P o 3 P o 6 P o 8 1一P 0 4 1一P 0 5 1一P o 7 P o 4 P 0 1 P 0 2 P o 5 P o 7 1一P 0 3 1一P 0 6 1一P 0 8 P o z P 0 3 P 0 4 P 0 5 P 0 7 1一P 0 1 1～P 0 6 1一P 0 8 P o v P 0 8 P 0 5 P o 2 P 0 4 1一P o 6 1一P o 1 1一P 0 3 p 0 5 p 0 6 P 0 7 P 0 2 P 0 4 1一P 0 8 1一P 0 1 1一P 0 3 P o s P 0 5 P 0 6 P o l P o 3 1一P o 7 1一P 0 2 1一P 0 4 p 0 6 P 0 7 P o 8 P 0 1 P 0 3 1一P o 5 1 ‘ 一 P 0 2 1一P 0 4 4 p6 1 P o 3 P 0 4 P 0 8 P 0 7 P o 1 P o 6 1一P 0 2 1一P 0 5 P 0 3 P 0 4 P 0 8 P o 7 P 0 2 P o 5 1一P 0 1 1一P 0 6 P 。 1 P 。 5 P 0 8 P 0 4 P 0 2 P 0 7 1一P 0 3 1一P 0 6 P 0 1 P 0 5 P 。 8 P 0 4 P o 3 P 0 6 1一P o 2 1一P 0 7 P 0 1 P 0 5 P 0 6 P 0 2 P o 4 p o 7 1一P 0 3 1一P 0 8 P o 1 P 0 5 P 0 6 p 0 2 P 0 3 P 0 8 1一P o 4 1一p 0 7 P o 2 P 0 6 P 0 7 P o 3 P 0 1 P o 8 1一P 0 4 1一P 0 5 P o 2 P o 6 P 0 7 P 0 3 P 0 4 P 0 5 1一P 0 1 1一P 0 8 P 0 l P 。 2 P 0 3 P o 4 P 0 5 P 0 7 1一P 0 6 1一P o 8 P o 1 P o 2 P 0 3 P 0 4 P o 6 P 0 8 1一P o 5 1一P 0 7 P 0 5 P 0 6 P 0 7 P 0 8 P 0 l P o 3 1一P o 2 1一P 0 4 P o 5 P o 6 P o 7 P 0 8 P 0 2 P o 4 1一P o 1 1一P 0 3 5 p P o 3 P o 4 P 0 1 P 0 5 P 0 6 P o 7 1一P o 2 1一P 0 8 P o l P 0 2 P 0 3 P o 7 P 0 8 P 0 5 1一P 0 4 1一P 0 6 P o 4 P 0 l P 0 2 P 0 6 P 0 7 P 0 8 1一 P 0 3 1一P 0 5 P o 2 P 0 3 P 0 4 P 0 8 P o 5 P o 6 1一P 0 1 1一P o 7 P P o 1 P 0 2 P o 3 P 。 4 P 0 5 P o 6 P 0 7 1一P 0 8 P 0 l P 0 2 P o 3 P 0 4 P 0 6 P 0 5 P 0 8 1一P 0 7 P 0 1 P 0 2 P 0 3 P 0 4 P 0 5 P 0 8 P 0 7 1一P 0 6 P 0 1 P 0 2 P 0 3 P 0 4 P 0 6 P 0 7 P 0 8 1一P 0 5 P o l P o 2 P o 3 P o 5 P 0 6 P 0 7 P o 8 1一P 0 4 P o 4 P o l P 0 2 P 0 5 P 0 6 P 0 7 P 0 s 1～P o 3 P 0 3 P 0 4 P 0 I P 0 5 P 0 6 P 0 7 P o 8 1一P 0 2 P 0 2 P o 3 P o 4 P 0 5 P o 6 P 0 7 P 0 8 1～p 0 1 8 pI p0 1 p0 2 po 3 p0 4 p0 5 po 6 p0 7 po 根 据式 3 ～ 8 ， 求 和 得到 1级原 胞 损伤 概率 P 为 P。 p4 。 p 。 p p P P 6 7 8 9 重 整化后 ， 同样 的方法 可用 于更 高级数 P 。p p i p 十 p 十 p p 1 0 若 最 大级原 胞对 应 于某 一有 限 的混凝土 结 构体 ， 则 可通过 式 1 0 递 推得 到该 结构 的整 体破 坏概 率 P ． 由 于实际结构各单元损伤概率极有可能不等 ， 并非是一个常数 ， 所以， 式 1 o 是一个通式 ； 而传统的岩石 等概 率破碎模型式 2 只是针对单元损伤概率相等情况下的一个特例． 3 整体 安全 性评价 算例 根据以上算法 ， 利用 MA T L A B自编 了计算程序． 以某一正方体大体积混凝土结构为例 ， 对边长进行 8等 学兔兔 w w w .x u e t u t u .c o m 4 水 利 水 运 工 程 学 报 2 0 1 0年 1 2月 份 ， 划分成 8 即 5 1 2 个小正方体 1级单元 ， 计算了 8 1 组工况． 结构底面以上 8层 1 级单元依次称为 1 层 、 2层、 ⋯⋯、 8层 见 图 3 ． 对 于 实际 结构 而 言 ， 结 构 底 部 的损 伤对 抗 渗 、 受 力 不 利 ， 极有 可能使 损伤 进一 步加 剧 ， 对 于遭 受过 地震 的结 构 ， 上部 结 构 地震响应较大 ， 损伤概率也较大． 假设 1～ 3层和 6～8层 、 4和 5层的 1级单元损伤概率分别相等， 概率变化范围为 0 ． 1 ～0 ． 9 ． 计算得到结构整体破坏概率 P曲线如图 4所示 ， 图中横坐 标为 1～ 3层和 6～ 8层的 1级单元损伤概率值 用 表示 ， 并以 4层和 5层的单元损伤概率 P P 为系列． 由于 P P ≤ 0 ． 5时 ， 整体破坏概率 P趋近于 0 ， 所以， 图中主要列出了 0 ． 5 ≤ P P ≤0 ． 9具有代表性的曲线． 根据图 4 ， 从结构整体破坏概率 曲线趋势看 ， 整体破坏概率 P随单元损伤概率 的增大而增大． 当 P P ≤0 ． 5时 ， 即使其他 单元损伤概率较大 ， P值仍很小， 趋近于 0 ； 当 P P 0 ． 6 ， 一旦 P o I0 ． 8时 ， P开 始 明显 增 大 ； 当 P P ≥0 ． 7 ， P o ≥0 ． 8时 ， P随 P 的增大而急剧增大． 所以， 总体上讲 ， 当单元损伤概率不超过 0 ． 7时， 整体破坏概率很小 ， 结构较为安全， 不会发生垮塌突发 事件． 但单元损伤概率一旦超过 0 ． 7 ， 整体安全性会 迅速 降低． 此时 ， 若外界环境恶化 ， 极有可能发生垮塌事件． 这类似于传统 重整化群方法的临界值问题． 4 结 语 图 3 计算 网格 F i g ． 3 Co mpu t i n g g r i ds ％ 图 4 结构整体破坏概率曲线 Fi g ． 4 Pr o ba bi l i t y c u r v e s o f ma s s i v e s t r uc t u r e da ma g e 针对单元损伤概率不等的正方体结构 ， 推导了 1 级原胞损伤概率 ， 再依次确定 ／ 7 , 级原胞 的损伤概率 ， 从 而得到结构整体的破坏概率． 传统的岩石等概率破碎模型只是针对单元损伤概率相等情况下的一个特例 ， 而 该数学模型是一个通式． 以某大体积混凝土结构为例 ， 应用所建立 的数学模型进行了实例分析． 从破坏概率 曲线趋势看 ， 整体破坏概率 P随着单元损伤概率的增大而增大． 当单元损伤概率不超过 0 ． 7时， 结构较为安全 ， 单元损伤概率一旦超过 0 ． 7 ， 结构整体安全性会迅速降低． 若外界环境发生恶化， 结构极有可能发生垮塌． 值得一提的是 ， 本文所提出的方法虽然也涉及单元问题 ， 但与有限元有本质的区别． 有限元是 以结点位移 为基本未知量 ， 通过结点力与结点位移间的力学关系， 求出结点位移， 进而求 出应力． 而本文是为整体安全性 评 价提 出 的一种新 思 路 ， 建 立 了局部 破损 概率 与整 体安 全性 之 间的关 系 ， 从 定量 的角 度评 价 了正方 体结 构整 体破坏概率 ； 但仍有这几方面需要进一步研究 1 非正方体结构整体破坏概率模型 的建立． 本模型只适用 于正方体结构 ， 并划分成正方体单元 ， 而一般混凝土结构不具备此条件 ， 经常包含斜面或曲面， 所以需分别对 结构体和单元体进行处理 ， 建立符合实际混凝土结构的整体破坏概率模型． 2 混凝土结构损伤不可避免会 受到多种不利因素的影响， 即使局部损伤概率相同 ， 但损伤位置 、 工作环境不同 ， 局部损伤对结构整体安全性 的影响也不 同， 所 以还需考虑权值问题下的整体破坏概率． 3 进一步开展外力及其他环境因素对结构整体 破坏概率的影响规律研究． 目前已对以上部分问题进行 了初步探讨 ， 将陆续在其他相关专题里进一步阐述． 参考 文 献 [ 1 ]杨华舒 ， 施 延华 ， 杨 志刚．用与承载能力相关 的指标诊断水工混凝 土结构 的老化[ J ] ．水 电能源科学 ， 2 0 0 2 ， 2 0 2 1 7 - 1 9 ， 7 5 ． Y A N G Hu a ． s h u ，S HI Y a n h u a ，Y A N G Z h i - g a n g ．D i a g n o s i n g a g i n g o f c o n c r e t e c o n s t r u c t i o n o f c i v i l e n g i n e e r i n g b y c a r r y i n g c a p a c i t y d e s c e n d i n g [ J ] ．H y d r o e l e c t r i c E n e r g y ， 2 0 0 2 ， 2 0 2 1 7 1 9 ， 7 5 ． i n C h i n e s e [ 2 ]林皋．混凝土大坝抗震 技术 的发 展现状 与展 望 I [ J ] ．水科 学与 工程 技术 ，2 0 0 4 6 1 3 ． L I N G a o ．D e v e l o p m e n t s t a t u s a n d p r o s p e c t o f s e i s m i c t e c h n i q u e a b o u t c o n c r e t e d a m I [ J ] ．He b e i Wa t e r R e s o u r c e s a n d H y d r o p o w e r E n g i n e e r i n g ， 2 0 0 4 6 1 ． 3 ． i n C h i n e s e 学兔兔 w w w .x u e t u t u .c o m 第 4期 顾培英 ，等 基于重整化群 的水工混凝土结构整体破坏概率研究 5 [ 3 ] [ 4 ] 林皋．混凝 土大 坝抗 震技 术 的发展 现状 与展 望 1 I [ J ] ．水科 学 与工程 技 术 ， 2 0 0 5 1 1 3 ． L I N G a o ．D e v e l o p m e n t s t a t u s a n d p r o s p e c t o f s e i s mi c t e c h n i q u e a b o u t c o n c r e t e d a m Ⅱ [ J ] ．He b e i Wa t e r R e s o u r c e s a n d H y d r o p o w e r E n g i n e e r i n g ， 2 0 0 5 1 1 3 ． i n Ch i n e s e L I N G a o ， H U Z h i - q i a n g ， S H I Y u n ． x i a o ， e t a 1 ．S o me p r o b l e m s o ff t h e s e i s m i c d e s i g n o f l a r g e c o n c r e t e d a m s [ c ] ff P r o c e e d i n g o f 1 3 t h Wo r l d C o n f e r e n c e o n E a r t h q u a k e E n g i n e e r i n g WC EE，Va n c o u v e r ，B C C a n a d a ，2 0 0 41 0 8 5 ． [ 5 ]林皋．汶川 大地 震 中大 坝震 害 与 大坝 抗震 安 全性 分 析 [ J ] ．大连 理 工大 学学 报 ，2 0 0 9 ， 4 9 5 6 5 7 6 6 6 ． L I N G a o ． P e r f o r m a n c e o f d a m s s u f f e r e d i n We n c h u a n e a r t h q u a k e a n d s e i s mi c s a f e t y a n a l y s i s o f d a m s [ J ] ． J o u r n a l o f D a l i a n U n i v e r s i t y o f T e c h n o l o g y ，2 0 0 9，4 9 56 5 7 6 6 6 ． i n Ch i n e s e [ 6 ]涂劲 ， 陈厚群 ，张伯 艳．高拱坝体 系整体抗 震 安全 评价 方法 研究 [ J ] ．世 界地 震工 程 ， 2 0 0 7 ， 2 3 1 3 1 - 3 7 ． T u J i n ， C H E N H 0 u q u n ， Z H A N G B o ． v a n ．T h e m e t h o d o f h o l i s t i c e a r t h q u a k e s a f e t y e v a l u a t i o n o f a h i g h a r c h d a m s y s t e m[ J ] ． Wo r l d Ea r t h q u a k e En g i n e e r i n g， 2 0 0 7，2 3 1 3 1 3 7 ． i n C h i n e s e [ 7 ]李周顺 ， 何江达 ， 苏 向震 ， 等．某重力坝坝体动力结构 特性 及坝基动力稳定性分析 [ J ] ．四川水利 ， 2 0 0 8 ， 2 8 1 1 3 1 7 ． L I Z h o u s h u n ，H E J i a n g d a ， S U X i a n g z h e n ，e t a 1 ．D y n a mi c c h a r a c t e r i s t i c s a n d d y n a m i c s t a b i l i t y o f g r a v i t y d a m[ J ] ．S i c h u a n Wa t e r R e s o u r c e s ， 2 0 0 8 ， 2 8 1 1 3 1 7 ． i n C h i n e s e [ 8 ] 韦凤 年． 分 析总结汶川地震经验教训加 紧组织水 工抗 震设计规范修订一访 中国工程 院院士陈厚群 [ J ] ．中国水 利 ， 2 0 0 8 ， 6 0 5 1 1 3 _ 4 ． WE I F e n g ． n i a n ．D r a w i n g l e s s o n s f r o m We n c h u a n e a h q u a k e a n d r e v i s i o n o f c o d e f o r s e i s mi c d e s i g n o f h y d r a u l i c s t r u c t u r e s i n t e r v i e w wi t h a c a d e mi c i a n C HEN Ho u q u n[ J ] ．C h i n a Wa t e r Re s o u r c e s ，2 0 0 8，6 0 5 1 1 3 - 4 ． i n C h i n e s e [ 9 ]苏怀智 ，吴 中如 ， 戴会超 ．初探大坝安全智 能融合监控体系 [ J ] ．水力发 电学报 ， 2 0 0 5 ， 2 4 1 1 2 2 1 2 6 ， 5 2 ． S U H u a i z h i ．wu Z h o n g r u ，D A I H u i c h a o ．S y s t e m o f t h e o r i e s a n d m e t h o d s m o n i t o r i n g i n t e l l i g e n t l y d a m s a f e t y[ J ] ．J o u r n a l o f Hy d r o e l e c t r i c En g i n e e r i n g ，2 0 0 5，2 4 1 1 2 2 1 2 6，5 2 ． i n C h i n e s e [ 1 O ]朱伯芳 ， 张 国新 ， 贾金 生 ， 等．混凝土坝的数字监控一提 高大坝监控水平 的新途径 [ J ] ．水 力发 电学 报 ， 2 0 0 9 ， 2 8 1 1 3 0 1 3 6 ． Z HU B o f a n g， Z HANG Gu o x i n ，J I A J i n s h e n g ，e t a 1 ．Nu me r i c a l mo n i t o r i n g o f c o n c r e t e d a ms a n e w w a y f o r i mp r o v i n g t h e s a f e t y c o n t r o l o f c o n c r e t e d a ms [ J ] ．J o u r n a l o f H y d r o e l e c t r i c E n g i n e e r i n g ， 2 0 0 9 ， 2 8 1 1 3 0 1 3 6 ． i n C h i n e s e [ 1 1 ]特科特 D L ．分形与混沌 在地质 学和地球物理学 中的应用 [ M] ．陈颞 ， 郑捷 ， 季颖 ， 译 ． 北京 地震出版社 ， 1 9 9 3 2 1 2 0 1 ． T U R C O T T E D L ．F r a c t a l s a n d c h a o s a p p l y i n g i n g e o l o g y a n d g e o p h y s i c s [ M] ．C H E N Y o n g ， Z H E N G J i e ，J I Y i n g ， T r a n s l a t i n g ．B e i j i n g E a r t h q u a k e P u b l i s h C o mp a n y ，1 9 9 3 2 1 2 0 1 ． i n C h i n e s e [ 1 2 ]L I U B i n ，S HI F e n g ，G A O Y u j i n ．D y n a m i c l o a d b a l a n c i n g b a s e d o n r e s t r i c t e d m u l t i c a s t t r e e i n h o m o g e n e o u s mu h i p r o c e s s o r s y s t e m s [ J ] ．J o u r n a l o f B e i j i n g I n s t i t u t e o f T e c h n o l o g y ， 2 0 0 8 ，1 7 2 1 8 4 1 8 8 ． [ 1 3 ]马志涛 ， 谭 云亮．岩石破坏演化 细观非均 质物理元胞 自动机模拟研究 [ J ] ．岩石力学 与工 程学报 ， 2 0 0 5 ， 2 4 1 5 2 7 0 4 2 7 0 8 ． MA Z h i t a o ，T A N Y u n l i a n g ．S i m u l a t i o n s t u d y o f r o c k f a i l u r e b a s e d o n M H P C A mo d e l [ J ] ．C h i n e s e J o u r n a l o f R o c k Me c h a n i c s a n d E n g i n e e r i n g ，2 0 0 5，2 4 1 5 2 7 0 4 2 7 0 8 ． i n Ch i n e s e [ 1 4 ]J E N A B K， D H I L L O N B S ．A s s e s s m e n t o f r e v e r s i b l e m u l t i s t a t e k - o u t o f - n G ／ F ／ L o a d S h a r i n g s y s t e m s w i t h fl o w g r a p h mo d e l s [ J ] ．R e l i a b i l i t y E n g i n e e r i n g a n d S y s t e m S a f e t y ， 2 0 0 6 ， 9 1 7 7 6 5 7 7 1 ． Da m a g e p r o b a bi l i t y o f h y dr a u l i c c o n c r e t e s t r u c t u r e b a s e d o n r e no r ma l i z a t i o n g r o u p t he o r y G U P e i ． y i n g ，HU A N G Q i n g ． h o n g ，D E N G C h a n g 。 ，T A N G L e i 1 ．K e y L a b o r a t o r y o f W a t e r S c i e n c e a n d E n g i n e e r i n g ， M i n i s t r y o f W a t e r R e s o u r c e s ， N a n j i n g H y d r a u l i c R e s e a r c h I n s t i t u t e ， Na n j i n g 2 1 0 0 2 9 ，C h