碾压混凝土重力坝整体三维仿真分析.pdf

科研 与管理 水利 规划 与设计 2 0 1 1年 第 3期 碾压混凝土重力坝整体三维仿真分析 吕 文 丽 ( 广东 省水利水 电科 学研 究院 广 州 5 1 0 6 1 0 ) 【 摘要】 本文运用有限单元法对光照碾压混凝土重力坝进行施工 、蓄水及运行 全过程的温度 、应力应 变综合 仿真分析 ，得出了温度场 、位移场 、应力 场的分布规律 。文 中还对坝体 坝基系统进行 了全程综合仿真 分析与传 统不仿真分析 ，并将两种方案结果进行对 比，为工程设计和施工提供依据 。同时也为在狭窄河谷 地区建造 2 0 0 m 级 高的碾压混凝土坝提供可靠 的理论和实践经验。 【 关键词】 碾压混凝土重力坝温度应力应变全程综合仿 真分 析 【 中图分类号】T V2 2 2 1 【 文献标识 码】 B 【 文章编号】 1 6 7 2 2 4 6 9( 2 0 1 1 )0 3 一O 0 4 O O 5 【 D0I 编 码 】 1 0 3 9 6 9 j i s s n 1 6 7 22 4 6 9 2 0 1 1 0 3 0 1 3 1 引言 1 计算原理与分析方法 据 对 国内外 已建碾压 混凝土 坝 的运 行观测 ，总 体来说碾压混凝土筑坝技术是成功的，但己运行的 一 部 分坝 中出 现较 多 裂缝 ，如 加拿 大 的 R e v e l s t o k e 、美 国的 D wo r s h a k以及 中国 的水 东 等大 坝 ， 相继 出现裂缝 和漏水 。这 不仅影 响 了大 坝 的安全运 行，而且使人们对碾压混凝土坝的筑坝技术产生疑 虑，因而引起 国际坝工界 的高度重视 1 。 碾压混凝土坝产生裂缝的主要原因是在受约束 条件下，温度拉应力超过了混凝土的允许拉伸变形 能力从 而导 致开裂 。为 了寻找裂 缝 的形 成 、发展 和 扩展 的规律 ，在设 计 阶段 必须进 行温度 场及应 力 场 仿真分析。由于施工过程中的温度变化对碾压混凝 土坝 的应力状 态有 重要 影 响 ，随着碾 压 混凝 土坝 高 的不断增 加 ，这 种 影 响 会 越 来 越 大 。且 现 在 在 建 、 拟建的碾压混凝土坝绝大多数是碾压混凝土重力坝， 因此 ，碾压混凝 土重 力坝 成层 浇筑过 程 中 的温度 场 及应力场的仿真计算就成为一个重要研究课题。 本文 运用有 限单元 法对光 照碾压 混凝 土重力 坝 进行 施工 、蓄水 及运行 全过程 的温度 、应力 应变 综 合仿真分析，得 出了温度场、位移场、应力场的分 布规律。还对坝体 坝基系统进行 了全程综合仿 真 分析与传统不仿真分析 ，并将两种方案结果进行对 比，为工程设计和施工提供依据。 40 2 1 温 度场计 算 的基本原 理 根据热传导理论，三维非稳定温度场 应 满足 下列偏微分方程 ： 苦 等 + 等 + 筹 + 一 十 十 十 式 中 ：T 温度 ( ) ； a 导温系数，口一 印， 单位( m h ) ； Q 由于水化热作用，单位时间内单位体 积 中发 出 的热量 ，( k J m。 h ) ； C 混 凝土 比热 ，( k J k g) ； p 密度 ，( k g m。 ) ； r 时间 ，( h ) 。 2 2 徐变 应力计算 的基本 公 式 4 由于弹模 E和徐变度都 随时间而变化，不能 采用 常 规 的 方 法 ，可 以采 用 增 量 法计 算 ，把 时 间 划 分 为一 系列 时 段 ( 一 1 ，2 )。在 时段 内产生 的应 变增量 为 ： 一 e + e ) + ) + +We ) 式 中 ： e ) 、 Ae 、 e ) 、 、 ) 分 别为弹性、徐变、温度、自生体积、干缩应变增量； 应 力增 量和应 变增量 关 系为 ： 作者简介 ：吕文丽 ( 1 9 8 2年一) ，女，汉族 ，河南人 ，硕士研 究 生 。 科研 与管理 水利规 划与设计 2 0 I 1 年 第 3期 一 D ( e 一 刁 ) 十 + + A ) 计算 的平 衡方程 为 ： K ) 一 P + P ) + P + P + AP 其 中： K为 整 体 刚 度 矩 阵； P ) 、 P 、 P 、 P 、 P 分别 为外 载 荷 、徐变、温度、自身体积变形和干缩引起的结点 载 荷增量 。 在计算应力时，将求得 的节点温度作为载荷 ， 同时考虑徐 变的影 响 ，根据混凝 土实 际的浇筑情 况 设 置荷载步 ， 读取 相应 的温 度荷 载 ，徐 变度 的拟 合 公 式为 ： C ( t ，r ) 一( z 1 + ,2 7 2 Z - x 3 ) 1 一 e - X 4 + ( 5 + z 6 r却) 1 一 e - X 8 其中： z ，5 17 ，5 1 7 是 8个待定的参数 ，由 于参数比较多，如何根据试验资料来确定这些参 数 ，是混凝土徐 变 问题 中 的一个 难点 。这 样就涉 及 到 多参数 的曲线拟合 问题 ，求解 方法有 复形法 、罚 函数法 、序列线 性规划法 、可行 方 向法等 等 。本 文 采用的是复形法对这 8 个参数进行拟合。 3 计算模型与计算方案 3 I 工程简介及计算模型 光照水 电站 位于贵州 省关岭 县与晴隆 县交界 的 北 盘江 中游 ，是 北盘江干 流梯级 的龙头 电站 ，电站 以发 电为 主 ，其 次航运 ，兼顾灌 溉 、供水 及其他 综 合效益。光照水电站碾压混凝土重力坝最大坝高达 2 0 0 5 0 m，最 大底 宽 1 5 9 0 5 m，顶 宽 1 2 m，坝 顶 长 4 1 0 m。大坝 由左 、右岸 挡水 坝段 及 溢流 坝段 组 成 ， 长度 分 别 为 1 6 3 m、1 5 6 m 和 9 1 m。正 常 蓄 水 位 为 7 4 5 m。对大坝整体进行分析，有限元计算模型基本 采用六面体八结 点单元 。基础 向底 部延伸 3 0 0 m，向 上、延伸 2 5 0 m，下游 3 5 0 m。为了确保模 拟精 度， 计算模型取 3 m一个浇筑层。网格共划分 1 0 3 3 7 5 5个 单元，1 0 1 2 0 4 2 个结点，其 中坝体部分 6 4 9 3 5 3个单 元，6 8 9 0 9 8 个结点。计算网格图见图 1 。 图 I 整体模型网格 图 图 2 整体模型网格 ( 坝体放大) 3 2 施 工过程和 蓄水过程 本方案将趾槽部分( 即 5 5 0 m 高程到 5 5 5 m高 程)作为一个浇筑层，开始施工时间为 2 0 0 5年 1 2 月 2 5日，从 2 0 0 6年 2月 6日开始 5 5 5 m高程以上 的浇筑施工 。碾压混凝土坝一般采用通仓薄层浇 筑 ，坝体浇筑浇筑过程按下述原则简化：以溢流坝 段 为基 准 ，同时考 虑陡坡 、底孔两 个坝段 ，根据施 工进度表 ，原则 上按 3 m 一个 浇筑层 。 根据关键点的温度监测值，反演影响温度场的 主要参数 ，即浇筑温度、通水温度 、通水时间等。 一 期通水冷却和二期通水冷却按分段函数模拟。一 期冷却中，采用等效法计算，把冷却水管看成内部 热源 ，在平均意义上考虑水管的冷却效果。其中在 处 理有热 源 的水管 冷却 问题时 ，混 凝土 的绝 热温升 公式采用双曲线型。进行坝体二期冷却时 ，水泥水 化热 已基 本 散 发 完 毕 ，可 以看 成 一 个 初 温均 匀 分 布 、无热源的温度场进行分析，只考虑水管冷却作 用计算 的混凝土平 均温度 作为绝热 温升 。 蓄水 计划严格 按照设 计院 给定 的 水 库蓄水进 展 统计表 ，水位 每 天都 有 变化 ，但 变化 较 小 。本 方案将水位每变化 1 0 m作为 1 次蓄水过程。水位 从 6 0 2 3 5 m 升 高到 7 4 5 0 m，整个 蓄水分 为 1 6次 。 3 3 计算 方案 为 了 比较考虑 温度徐 变应力 和不考虑温 度徐变 的区别 ，分别采用 了两种计 算方案 。 ( 1 ) 方案一 ( 仿真 方 案 ) ：考 虑温 度 和 徐变 作 用 ，同时 考虑坝体 实 际 浇筑 过程 和 分期 蓄水 过 程 。 考虑的荷载有坝体 自重、温度荷载、渗透荷载及作 用在坝体上的静水压力。计算时按照边浇注边蓄水 的实 际施 工 过 程 ，坝 体 分 层 浇 注 ，每 层 3 5 m， 2 0 0 7年 1 2月 3 1日开始 蓄水，蓄水过程分 为 1 6 次 ，每次都是在坝体分层浇注到一定高程后 ，再一 次性将水位升到相应高程。蓄水计戈 4 与温度场计算 中的蓄水计 划一致 。 ( 2 )方案二 ( 不仿真方案) ：不考虑温度和徐 变作用，不考虑坝体浇筑过程和分期蓄水过程。考 虑的荷载有坝体 自重 、渗透荷载及作用在坝体上的 静水 压力 。 3 4仿真分 析流程 图 在温度场和应力场的仿真计算中，由于要考虑 外界环境、材料参数、施工过程和冷却措施等因素 在时间和空间中的不断变化 ，使得即使采用现有软 件和程序，仍必须将软件和程序的有关功能加以整 合 ，灵 活运用 。为 了实 现仿 真 分 析这 一 复杂 过程 ， 41 科研与管理 水利规划与设计 2 0 1 1 年第 3期 编制仿真温度场和应力场分析控制程序 ，使建立模 型和分析过程实现参数化 。 仿真分析过程流程图见图 3 。 图 3 仿真计算流程图 4 计算成果及分析 4 1 温 度场仿 真分析 在计算 的过程中按照实际的施工进度模拟出了 任何一个时段的温度分布规律，限于篇幅，文中选 出溢流坝段典型断面几个 有代 表性 的结 果进行分 析 。 ( 1 )图 4为 2 0 0 8 0 7 0 5年典型坝段温度等 值线图，图 5是 2 0 0 8年坝体浇筑完工 3个月后 的 温度等值线 图。从图 4可以看出，施工期间最高温 度出现在高程 6 7 0 7 0 0 m 之间。这是 因为，这部 分混凝土浇筑 时间为夏季高温季节，入仓温度较 高。从温度等值线图可以看出，坝体表面的温度基 本等于外界气温值，内部由于水化热温升，温度值 较大 ，在坝体边界处温度等值线较密，温度变化幅 度较大。从图 5可以看出坝体 内部的温度逐步降 低，上下游水温气温影响的范围向坝体 内部深入。 ( 2 )从图 6中可以明显看出这个高程有一期通 水和二期通水两个过程，A点整个曲线刚开始随气 温变化 ，在粘土墙铺筑后温度升高，待通水一段时 间后 ，温度开始降低 ，通水结束，温度继续升高 ， 直到二期通水 ，温度再次下降，最后趋 于平稳。B 点变化滞后于 A 点。C点处于坝体中心 ，开始升 温，待通水一段时 间后 ，温度开始 降低 ，通水结 束 ，温度 继续升高 ，直到二期通水 ，温 度再次下 42 降，最后趋于平稳。图 7 所显示的高程均在 6 1 5以 上，在初期 阶段，A，B，C三条直线 的温度变 化 规 律均 是开始 升温 ，这是 由于水 化热温 升远 大于冷 却水管的降温作用 ，待通水一段时间后 ，温度开始 降低，直至通水结束 ，然后温度继续升高，A 曲线 随着气温的变化而变化 ，到 0 8年 3月后 ，A 曲线 呈现梯状变化 ，出现 明显的温度陡升和陡降 ，这是 由于蓄水时水位先升后降造成的，最后水温稳定在 库底水 温 1 3 左右 。 图 4 ( 左 ) 2 0 O 8 0 7 一O 5 典 型坝段温度等值线图 图 5 ( 右 ) 2 0 0 8 1 2 -1 0 典 型坝段温度等值线 图 1 I L ， 啪 26 1 I 图 6( 左) 5 8 5高程三个特征点温度历程线图 图 7 ( 右) 6 4 0高程三个特征点温度历程线 ( 圈 6中 A点位 于上游裹面 B点距 上游面 2 5 m；C点距上 游面 3 0 0 m；圈 7中 A点位于上游襄面 B点距上游面 1 O m；C点 距上游面 2 0 0 m) 4 2 应 力场仿 真分析 4 2 1位移 分析 在这 2种计算工况中位移分布有较大的不 同， 方案 1 计算工况下坝体 向下游倾倒，顺河向位移指 向下游 ，铅直向位 移沿高程从上到下位移逐渐增 大 ，方向向下 ，坝顶位移最小，由于其考虑分层浇 筑 ，每层浇筑块在浇筑前位移为 0 ，其位移是其在 浇筑过程中的实际位移 ，在静水压力的作用下，会 沿着水压力作用点发生扭转 。而方案 2 计算中，顺 河向位移指向下游 ，铅直向位移向下 ，坝顶位移最 大 。2种计算工况下特征点 的位移对 比见 表 1 ，2 种 工 况 下 最 大 的 不 同 之 处 在 坝 顶 ，差 别 可 达 3 9 mm。在坝踵和坝趾二者相差很小 ，仅为 2 mm。 科研 与管理 水 利规划 与设计 2 0 1 1 年 第 3期 图 8 ( 左) 方案 1蓄水后典型断面位移矢量图 图 9( 右) 方案 2蓄水后典型断面位移矢量 图 表 1 两种 方 案 下 断面 特 征 点 位移 比较( 单 位 ：mm) 工况 特征点位置 坝顶 坝踵 坝趾 铅直向 一3 3 5 4 3 4 9 2 4 4 O 仿真分析 顺 河 向 l 9 2 l 1 4 9 8 1 6 3 7 铅 直向 一4 2 1 3 。 一3 9 6 3 2 2 6 3 不仿 真分析 顺 河 向 2 8 7 1 7 2 1 1 6 2 8 4 2 2 应 力 分 析 图 1 0 图 1 7为两种计 算工 况下 典 型断 面铅 直 向应力等值线 图 、 顺 河 向应 力 等值 线 图、应 力 等 值 线 图及 应 力等 值 线 图。从 图 中可 以看 出两 种 计 算 工况 F 力的分 布规 律 总体 趋 势是 相似 。只是 在 方案 】中出现许多分布较为密集的层状等值线封闭 区 。 使得方案 1应力等值 线分 布不 均匀 ， 出现层状应 力 集 中。这是 因为在 坝体浇筑 和蓄水过程 中温度 作 用 一直存在 ， 只是随 着高程 的增 加 ， 坝体 上部浇筑 块 J 。 F游方 向厚 度逐 渐减 小 温 升较 之坝 体 下部 基 础 强 约束 内要小 ， 温度作 用也有所 减小 ， 由于温度应 力 造成 的层状 应力 集 中 区也渐 渐减 小 ； 但 是在 施 工 过程中坝体下游表面 ， 由于一直直接与空气接触 ， 受 气温影响， 表面出现一定范围的受拉区， 坝体表面与 寒冷的空气相接触， 而内部又由于水化热作用不断 放热 ， 使得下游 溢流堰表 面 出现较 大范 围的受拉 区 。 从 整个 浇筑和 蓄水 过程来 看 ， 温 度应 力 的 确会 造成 坝体 表面及其 内部 出现 拉应 力 集 中区 ， 特 别在 每次 浇筑块交接面上 ， 应力集中较为明显， 但由于温度作 用造成的坝体内的拉应力较小。 表 3 列出了特征点位置 的应力值 ， 两种计算 工况 下， 应力值方向相同， 量值的差别均在 1 7 MP a以内。 图 1 0 ( 左) 方案 1 典型断面铅直 向应力等值线 图 1 1 ( 右) 方案 2 典型断面铅直向应力等值线 图 1 2 ( 左) 方案 1 典型断面顺 河向应力等值线 图 图 1 3 ( 右) 方案 2 典 型断面顺河向应力等值线 ： 一 一 - 一 、一 、 I _l1 。 - 、一 a 、 0 0 一、 、 一 一 一 、 一 ， m l 、t ， s- ? ，1 ： 图 1 4 ( 左) 方案 1典型断 面第一主应力等值线 图 1 5 ( 右) 方案 2 典 型断面第一主应力等值线 攀 so,w) 一s=-is l；=-a 图 1 6 ( 左) 方案一典型坝段第三主应力等值线 图 1 7 ( 右) 方案二典型断 面第三主应力等值线 表 2 两种方案下断面特征点应力比较( 单位 ： MP a ) 工况 位置 铅直 向应力 顺河向应力 第一主应力 第三主应 坝踵 一3 1 2 1 4 6 1 2 6 3 2 3 仿真分析 坝 趾 一8 3 8 2 0 8 1 4 5 8 9 9 坝踵 1 4 4 0 3 3 0 8 8 1 6 2 不仿真分析 坝 趾 一8 4 0 2 1 0 1 O 8 1 0 1 8 5 结 论 ( 1 ) 只控制浇筑温度 ， 由于会出现气温倒灌 ， 而 同 时采用通水冷却措施， 则不仅能有效的降低整个施工 过程中的最高温度值 ， 也同时控制了最高温升值。 ( 2 ) 2种计算工况 中， 重力坝坝体 坝基系统的 位移规律有明显的不同， 方案 2计算工况下坝顶铅 直 向位移最 大 ， 方案 1 工 况下 由于分层浇筑 ， 分次蓄 水 的影 响 ， 坝 顶铅直 向位 移最小 。 ( 3 )在 这2 种 工 况 下 ， 应 力 值 相 差 约 为 1 7 MP a 。在方案 1 计算中坝体下部基 础强约束区 和坝体的上下游边缘的温度应力较大 ， 由于温度影 响 ， 层 间应力呈 现锯齿状 。 ( 4 ) 由于变形 、 应力在 高碾压混 凝土重 力坝设计 中都是重要指标， 所 以综合仿真分析具有意义 。同 时也 为在狭窄河 谷地 区建 造 2 0 0 m( 下转第 7 0页 ) 43 设计 与施 工 水利 规划 与设计 2 0 1 1 年 第 3 期 4 2 9 m， 按 锚索 内锚 段 长 8 9 m 考 虑 ， 最 长 锚 索 的 理论 长度不应 小于 5 2 m。 根据计算分析结果， 结合推测的边坡地质状况， 为确保 锚索 内锚段 位 于较 好 的 围岩 段 ， 锚 固设 计 长 度为 5 0 m， 现场实际施工过程 中按 4 0 5 0 m之间。 4 4敏感性 分析 结 构面物 理力学 参数是 影响边 坡稳定 的关键 参 数， 水荷载是主要外荷。在边坡楔形体稳定计算 中 对边坡结构 面连通率与暴雨情 况进行 了敏感性分 析， 边坡连通率取不同百分数， 采用加权平均方法计 算综合力学参数， 暴雨作用 以不同的地下水折减系 数进行计算 ， 主要结论如下 。 ( 1 ) 结构面连通率每降低 l 0 ， 边坡安全系数 值大概提高 0 1 ； 当第和第组这两组结构面的 力 学参数 取地质 提 供 的最 高 值 时 ， 两 组 结 构 面 的连 通率在 6 0 以下时就可以满足规范稳定要求 ； 当两 组结构面的力学参数取地质提供 的最低值时 ， 两组 结 构面 的连通 率 在 5 O 以下 时就 可 以满 足 规 范 稳 定要 求 。 ( 2 ) 加 固力 与水压力 折减 系数成 线性关 系 ， 折 减 系数 降低 0 1 ， 所需加 固力减少约 7 0 0 0 t ； 结构面连 通率对加固力影响也较大 ， 一般结构面连通率降低 1 O ， 所需 加 固力 可减少 6 0 0 0 8 0 0 0 t 。 4 5 边坡 支护参 数 的调 整 根据 开挖揭 露 的地质 情况 和现 场施工 的实 际状 况 ， 对支护参数做了以下调整。 ( 1 ) 土 体 锚 杆 的调 整 ： 9 9 6高 程 以 上 部 分 长 锚 杆、 锚索施工过程 中， 在深度 1 5 3 0 m 范围内仍有 塌孔 、 卡钻 等情况 较严重 ， 同时锚 杆孔 进行 固壁注浆 时 ， 在压力很小的情况下 ， 部分孔注浆量相当大。 业主组织专家进行专题讨论， 将原方案中土体 锚杆角度由上倾 2 O 。 调整为上倾 5 。 1 0 3 0 m高程以 上未施 工锚 杆人 岩 长度 由 2 0 m 调 整 为 9 m， 现 场 已 造孔 的达到或 超过上 述长度 的锚 杆立 即完成 后序施 工 1 0 3 0 m 高程以上第一期在厂房后边坡位置布 置 1 6根 I O 0 0 KN拉压分散型锚索。 ( 2 ) 2 0 0 8年 6月 份 ， 厂 房 后 边 坡 E L 1 0 1 0至 E L 1 0 0 0坡 面 ， 厂 左 0 +0 1 4 2 5 m厂 右 0 +0 2 0 7 5 m 范围内出现小范围卸荷坍塌， 坍塌深度 2 3 m。在 坍塌部 位边坡 基岩 面布设 1 6钢筋 网， 间距 2 0 2 0 c m， 边 坡 先 期 浇 筑 5 0 c m 厚 C 2 0贴 坡 混 凝 土 ； 补 打 2 5砂浆 锚杆 。 ( 3 ) 取 消 E L 9 9 0 E L 9 7 1 m 边坡范 围的网格 梁 ， E L I O I O E L 1 0 3 0 m 边坡 范 围适 当增 加 锚 索 ， 局 部 设 网格梁 。 5 边 坡 监 测 为 了在施 工 期 及运 行 期对 边 坡 进行 观测 ， 了解 边坡的变形情况 ， 在厂房后边坡共 布置 2个典 型观 测断面， 断面布置锚索测力计 、 多点位移计 、 锚杆应 力计 等 ， 目前各项 观测 数据波 动很 小 ， 呈 稳定收 敛趋 势 。 6 结 语 天花板水 电站 厂房后边 坡 自然坡高约 2 4 3 m， 永久 边坡 开挖高 度 8 5 m， 厂址 出露 地层 以第 四 系松 散堆 积物 为主 ， 基 础 主要 为 中厚层 至厚 层 状岩 屑石 英砂岩， 发育三组主要裂隙 ， 边坡的稳定计算分析与 支护 是厂房 设计 的关键 问题 。 边坡稳 定分 析与支 护设计 以地 质专业 提供 的岩 体力学参数为基础， 对原始边坡的稳定性进行 了分 析 ， 根 据边坡 设计 规范 的要求 ， 计算 分析 了各种 工况 下开挖边坡的浅层滑动稳定、 深层 滑动稳定及楔形 体稳定进行计算 ， 对关键支护措施的参数进行研究， 最后综合分析比较， 选取合理的开挖坡 比与支 护参 数， 对类似的工程有一定参考与借鉴。 ( 上接第 4 3页) 级高的碾压混凝土坝提供可靠的理论 和实践经验。 ( 5 ) 本 文尚未考虑 应力和渗流的耦合作用 5 。 有证据表明这种耦合作用对坝基岩体和成层碾压混 凝土的影响较大 ， 下一步研究将给予考虑。 ( 6 ) 本 文 的 结 果 主 要 适 用 于 重 力 坝 ， 对 拱 坝 结 构 ， 由于其 复杂性 ， 尚需进 行专题 研究 。 参考 文献 1 沈崇刚 中国碾压混凝土坝的发展成就与前景( 上) J 贵 70 州水力发 电， 2 0 0 2 ， 6 ( 1 6 ) 2 李春敏等 浅忻 我 国碾压混 凝土拱 坝技术 E J 水 利水 电技术 ， 2 0 0 1 ， 1 1 ( 3 2 ) ： 7 8 3 朱伯芳 有限单元法原理 与应用I- M 2版 ， 北京 ：中国 水利水 电出版社 ， 1 9 9 8 4 朱伯芳 大 体积 混凝 土温 度应 力 与温 度控 制 M 北 京 ： 中国电力 出版社 ，1 9 9 8 5 吉小明 ， 杨春和 ，白世伟 岩体结 构与岩体 水力耦合计 算模型I- J - 1 岩土力学 ， 2 0 0 6 ， 2 7 ( 5 ) ： 7 6 3 7 6 8