混凝土重力坝深层滑动稳定性评估软件模块研发.pdf

第 3 0卷 第4期 2 0 1 3年 4月 长 江科 学 院 院报 J o u r n a l o f Ya n g t z e Ri v e r S c i e n ti fi c Re s e a r c h I n s t i t u t e Vo l _ 30 No 4 Ap r 2 0 1 3 DO I ： 1 0 3 9 6 9 j i s s n 1 0 0 1 5 4 8 5 2 0 1 3 0 4 0 1 8 2 0 1 3 ， 3 0 ( 0 4 ) ： 8 5 9 0 混凝土重力坝深层滑动稳定性评估软件模块研发 崔书生 ， 苏怀智 ， 张帅 ， 顾浩钦 ( 河海大学 a 水文水资源与水利工程科学国家重点实验室； b 水资源高效利用与工程安全国家工程研究中心； c 水利水电学院， 南京2 1 0 0 9 8 ) 摘要： 基于混凝土重力坝深层滑动面上应力分布情况的数值计算成果， 结合分项系数极限状态设计方法 ， 研究了一 种评估混凝土重力坝深层组合滑动稳定性的改进分项系数法。在此基础上， 利用 D e l p h i 编程工具 ， 结合 A u t o C A D 二次开发技术对该方法实现了程序化 ， 得到了评估混凝土重力坝深层滑动稳定性的可视化软件模块。该方法和软 件模块， 充分利用了有限元数值模拟方法在计算应力场和考虑坝基及坝体变形影响等方面的优势， 又结合了可靠 度设计方法， 得出的稳定安全系数能够较好地反映重力坝的抗滑稳定特性， 且沿用了设计规范中的 K 表示方法， 易 于计算和普及 。 关键词： 重力坝； 抗滑稳定 ； D e l p h i ； A u t o C A D； 二次开发； A e t i v e X 中图分类号： T V 6 4 2 3 ； T P 3 1 1 5 6 文献标志码 ： A 文章编号： 1 0 0 1 5 4 8 5 ( 2 0 1 3 ) o 4 0 0 8 5 0 6 1 重 力坝深层抗滑稳定分析方法 建基面及坝基深层抗滑稳定分析 ， 是重力坝设 计 中一项非常重要的内容 ； 并且在保证施工、 运行的 安全性方面起着举足轻重的作用 ； 同时 ， 准确分析坝 体 的抗滑稳定特性也是重力坝开展其他工作的基本 前提。 目前 重力 坝坝基深层抗 滑稳定分析方法 大致分为 3种 ： 刚体极限平衡法 、 有限单元法和地质 力学模型法 。 随着计算力学和计算机软件及硬件的发展 ， 利 用有限单元法可以比较准确地计算出坝体及其坝基 的应变 、 应力和塑性区。依靠计算机 ， 有限单元法能 够精确建立模型 ， 并考虑各种材料的非均质 、 各 向异 性等特性 ， 因而能较准确地计算接触面 以及深层组 合滑动面上的阻滑力和滑动力 。J 。然而仅仅计算 出阻滑力和滑动力还不能判断大 坝的抗滑稳定性 状 ， 需要找到一种计算方法 ， 以较为准确地反映大坝 的抗滑稳定特性。本文首先利用有限元法准确计算 出重力坝深层滑动面上的阻滑力和滑动力 ， 然后结 合分项系数极限状态设计方法 ， 对每一个剖分单元 运用抗剪断强度公式 ， 定义一种新 的抗滑稳定安全 系数 的计算公式 ， 得到了一种计算混凝土重力坝 深层组合滑 动 的改进分 项系数法 。并采用 D e l p h i 语言， 结合 A u t o C A D二次开发技术对该方法实现了 程序化。该方法充分利用 了有限元能够计算出精确 应力场和考虑坝基及坝体的变形影响的优 势， 又考 虑了可靠度设计方法 ， 所以得出的 K 能够较好地反 映重力坝 的抗滑稳定特性 ， 并且本方法沿用 了设计 规范 中的 K 表示方法 ， 易于计算和普及。 2 混凝 土 重 力坝 深 层 组 合 滑 动 的 改 进分项 系数法 本方法 的基本思想 是 ： 引入 广义力 的概念 ， 将传统的作用效应函数 S ( ) 和结构及构件抗力函 数 R ( ) 广义化 ， 结合混凝土重力坝的实际地质情 况 ， 分析深层滑移模式 ， 从而决定这种广义力以及该 广义力对应的作用效应函数和抗力函数。 下面以平面状 2滑面和3滑面( 如图 1 所示) 为 例 ， 简述混凝土重力坝深层组合滑动的改进分项 系 数法 。根据混凝土承载能力极 限状态 的基 本组合 ， 设计表达式为 1 1 ， 、 Y o s ( G G k ， Q Q K ， k ) I ， 0 k ) 。 ( 1 ) ， d l 、 y m 式中： 。为结构 重要 性 系数； 为设计 状况 系数 ； 收稿 日期 ： 2 0 1 20 8 2 7 ： 修回日期 ： 2 0 1 21 01 7 基金项 目： 国家 自然科学基金 ( 5 1 1 7 9 0 6 6 ， 5 1 1 3 9 0 0 1 ) ； 新世 纪优 秀人 才支持 计划 ( N C E T一1 00 3 5 9 ) ； 国家重 点实验 室专项 经费 资助项 目 ( 2 0 0 9 5 8 6 9 1 2 ， 2 0 0 9 5 8 6 0 1 2 ) ； 中央高校基本科研业务费专项资金资助( 2 0 1 2 B 0 6 6 1 4 ) ； 江苏省第四期“ 3 3 3 工程” 培养资金资助项目 ( B R A 2 0 1 l 1 7 9 ) ； 江苏高校优势学科建设工程资助项 目( 水利工程) ( Y S 1 1 0 0 1 ) 作者简介： 崔书生( 1 9 8 7一) ， 男 ， 安徽六安人 ， 硕士研究生 ， 研究方 向为水 工结构工程安全监控理论和方法 ， ( 电话 ) 1 3 7 7 0 5 1 6 7 1 2( 电子信箱 ) s h u s h e n g j i a y o u 1 6 3 t o m。 长 江科 学院院报 S () 为作用效应函数 ； ( ) 为结构及构件抗力函 数 ； G 为永 久作用标 准值 ； 。为永久 作用 分项 系 数 ； Q 为可 变作用标准值 ； 。为可变作用 分项 系 数； 口 为几何参数标准值； 为承载能力极限状态 基本组合的结构系数 ； 为材料性能标准值 ； 为 材料性能分项系数。 ( a ) 2 滑面( b ) 3 滑面 图 1 平 面状滑面 图 F i g 1 F l a t s l i p s u r f a c e s 对于图 1中所示 的平面状 2滑动面和 3滑动 面， 由于各个单元及滑动面上的滑动力和阻滑力不 可以直接 加 减 ， 所 以必须 将 式 ( 1 )中的 S() ， R( ) 进行广义化。事实上 ， 当大坝及坝基沿着组 合滑动面出现抗滑失稳的时候 ， 滑移体 ( 由坝体及 滑动面以上 的岩体组成 ) 受到边 界的约束 ， 将绕着 滑移体瞬时滑移 中心做转动。在此 ， 将阻滑力和抗 滑力广义为阻滑力偶和抗滑力偶。 具体计算方法及步骤 ： 对于其中的第 1 个滑动 面， 首先根据有限单元法 的应力计算成果 ， 利用式 ( 2 ) 和式 ( 3 ) 将滑动面相关单元 的应力合成为法 向 力 和切向力 Q。 ， 同理第 2个滑动面为法向力 2 和切向力 Q ； 除了确定力的大小外 ， 还需确定力 的 作用点。将滑动面上的各单元法向力看成是作用在 该滑动面上的平行力 系， 利用求解平行力系中心 的 方法可以求得作用点 O ， 和 O 的位置。详细求解方 法此处不再赘述， 见 理论力学 J 。 过作用点 O ， O ：画力 ， 2 的延长线交 于点 0。 利用结构计算成果应力场， 得到下面的法 向力和切向力的计算公式 ： Q= A ； ( 2 ) i = 1 n ， 、 = A ( L+ )。 ( 3 ) i ： l 、 m Tm 式 ( 2 ) 和式( 3 ) 中： n为组合滑动面上划分 的单元个 数 ； A 为第 i 个单元上 的滑动面的面积 ( 二维情况下 为长度 ) ； 为组合滑动面上单元 i 沿滑动面法 向 的正应力 ； 为组合滑动面上单元 i 沿滑动面切 向 的切应力 为组合滑动面上材料之间的摩擦系数； C 为组合滑动 面上材料之问的黏聚力 ； 其余符号意 义同前 。选取点 O作为瞬时转动中心 ， 对 中心取矩 后 ， 得到广义的 5 () ， R ( ) 计算公式 ： S ( )=Q Z l+Q 2 f 2 ； ( 4 ) ( ) = f l + A 1 z 。 + f + A ， 1 f 。 ( 5 ) 式( 4 ) 、 式( 5 ) 中： A ， A 分别 为组合滑动面中第 1 、 第 2滑动面的面积 ； f ， f ： 分别为组合滑动面中第 1 、 第 2滑动面到交点 O的距离； 其余符号意义同上。 如图 1 ( b ) 所示 ， 对于 3滑面上的每一个滑动 面 ， 按照 2滑面的原理得到每个面上 的法向力 和切 向力及其作用点分别为： N ， Q ， O ； N 2 ， Q ， O ； N 3 ， Q ， ， O ， 。和 2滑面不同的是 ， 经过 。 ， 2 和 的 3条直线正常不会交于一点 ， 所以不能确定瞬心。 但是 ， 根据 理论力学中的知识 ， 可 以将任意 2条 直线的交点作为名义瞬心 ， 将另一边 的法向力 向该 点转化为 1个力矩和 1 个力。然后分别对名义瞬心 取矩 ， 方法 同2滑面。具体的广义 s ( ) ， R( ) 计 算公式为 ： S ( ) =Q1 f l +Q 2 f 2+Q 3 f 3一N 2 d 2 ； ( 6 ) R ( ) = f N i + A )f f。 ( 7 ) 式( 6 ) 、 式( 7 ) 中： d 为瞬时转动中心 O点到第 2个滑 动面上法 向合力 2 作用线的距离， 当 2 对 O点的 矩为逆时针时取正， 反之取负。其余符合意义同上。 因此 ， 根据以上对 S ( ) ， R( ) 的计算 ， 定义 混凝土重力坝组合滑面抗滑稳定安全系数 ， 具体 定义为 K s= S 。 ( 8 ) y 0 d 1( ) 对于 3滑面的情况 ， 因为其 中可 以算出 3个不 同的 K ， 在这里取 3个当中最小的作为 3滑面的抗滑稳 定安全系数。 比较式 ( 1 ) 可得混凝土重力坝深层抗滑稳定分 析中的组合滑动面的失稳判据为 、 r 1 稳定 临 ) 除此以外 ， 对于圆弧状单滑面， 平面状单滑面和 多滑面等组合滑动面， 具体计算方法和平面状 2滑 面或 3滑面类似 ， 这里不再赘述 。 8 8 长江科 学院院报 2 0 1 3丘 a c a d ： =C r e a t e O l e O b j e c t ( A u t o C A D A p p l i c a t i o n 1 7 )a s A c a d A p p l i c a t i o n； l a y e r ： a c a dAc t i v e Do c u me n t La y e r s Ad d ( g l i d e p l a n e ) ； 创建名字为g l id e p l a n e 的图层 S h o w M e s s a g e ( 下 面将在 A u t o C A D 的 g l i d e p l a n e图层 中绘出组合滑动面 ! ) ； a c a d Vi s i bl e：=Tr u e； a c a d A c t i v e D o c u me n t A c t i v e L a y e r l a y e r ； f f 设置当前图层为 l a y e r a c a d mo s pa c e： a c a d Ac t i v e Do c u me n t Mo d e l S p a c e ； 提取当前活动文档为模型空间 从 u n i t 表中读取单元号和其 中的节点 ， 然后根 据节点在 c o o r d i n a t e表中查询相应节点号的坐标 ， 作 为绘制直线的坐标 ， 从而逐个绘出有限元剖分单元。 第 3步： 这一步是整个软件的核心， 即找出有限 元网格中和组合滑动面相交的网格 ， 本文称之为相 关单元。当找出所有和组合滑动面相交的相关单元 之后 ， 就可以通过这些单元信息 ， 从有限元计算软件 的所计算的应力场导出和组合滑动面相关的单元应 力信息 ， 最后读取这些 应力 ， 通过弹性力学知识 ， 换 算成滑动面上的正应力和切应力。然后和单元面积 相乘 ， 转化为我们想要知道的滑动力和抗滑力。最 后通过改进的组合系数法 ， 求出抗滑稳定安全系数 ， 判断重力坝深层抗滑是否稳定。 寻找相关单元的核心思想是 ： 通过 A u t o C A D给 出的对象模 型树 中提供 的 I n t e r s e c t Wi t h方法 ， 遍历 每个单元的每条边 ， 求出组合滑动面和单元 的边的 交点 ， 确定个数 ， 如果个数不为零， 则确定该单元就 是相关单元 ， 并且记下滑动面和该单元的交点， 求 出 2交点之间的距离 d和滑动面的方向余弦 f ， m， 为之 后求解抗滑力和下滑力做准备 ， 具体代码这里不再 赘述 。 在找出相关单元并且把相关单元的信息 ， 包括 单元号 、 直线号 、 滑动面和相关单元 的交点读入 r e - l a t e u n i t 数据表之后 ， 通过 d e l p h i 的 a d o q u e r y 控件根 据单元号在数据 u n i t s t r e s s中查询该单元的应力 ， 即 o r ， Or ， 这里所用单元的应力为单元形心处 的 应力 。 然后根据弹性力学中的知识求出组合滑动面上 的应力 ， 即： o r=1 2 or +m or +2 m l r ； ( 1 0 ) 下： ( Z +， n r ) +( mor +Z 丁 ) 一or 。 ( 1 1 ) 式 ( 1 0 ) 、 式( 1 1 ) 中： Z ， m， or， 分别表示滑动面的方 向余弦和斜面上 的正应力和切应力 ； 然后代人上面 式子算出 K ， 再利用式( 9 ) 判断重力坝的稳定情况。 最后 ， 需要强调的是 ： 本文抗滑稳定安全系数的 计算精度与有限元应力结果密切相关 ， 而有 限元的 应力往往是存在误差的 ， 所 以滑动面上 的单元应细 分 ， 以保证应力求解精度。 5 算例分析 某混 凝 土重力 坝 ， 最大 坝高 9 7 m， 坝顶 全 长 3 8 6 5 m， 共 由2 1 个坝段组成。其中， 5 坝段偏于右 岸 ， 是实体重力坝 ， 坝高为 6 0 m， 顶宽和底宽分别为 6 m和 4 5 m， 下游坝坡 为 1 ： 0 7 。坝基 内有倾 向下 游的与水平面呈 1 8 。 夹角的软弱结构 面， 下游抗力 体内有倾 向上游角度为 1 O 。 的缓倾 角软弱结构面。 上述软弱结构面共 同形成连续滑动面( 见图 4) ， 必 须验算坝基深层抗滑稳定。 图 4有限元剖分网格及滑动 面示意图 F i g 4 Fi n i t e e l e me n t me s h e s a n d t h e s l i p s u r f a c e 应用本文所述 的改进分项系数法验算该重力坝 深层抗滑稳定。首先建立该坝段的有限元分析数值 模型。为了简化实例计算 ， 将该坝段视为平面应变 问题 ， 利用平面三角形和 四边形等参单元模拟。考 虑该坝段及其基础条件 的实际情况 ， 选定其有 限元 计算范围为 ： 上游约 1倍坝高， 取为 6 0 m； 下游约 2 倍坝高 ， 取为 1 2 0 m； 坝基 以下取为 6 0 m。根据上述 分析结果 ， 5 坝段有限元 网格剖分结果如图 4所示 ， 该图为利用有限元软件导 出的节点坐标信息 ， 利用 本系统在 C A D中绘制的网格示意图。该模 型共有 1 2 5 2 个 单元 ， 1 2 0 2 个 节点。模 型 中， 上游水 位为 5 8 m， 计算结果见图 5 。 完成有 限元计算之后 ， 将本程序需要的数据按 照一定的格式导出， 通过相关程序处理之后导入到 数据库相应的数据表 中( 库中各表信息在第 4节中 已做详细介绍 ) ， 作为抗 滑稳定 的已知信息 。所需 要 的数据包括 ： 各节点的坐标信息 ， 各单元的信息以 及各单元沿 方向和 Y方向的应力 。 第4期 崔书生 等 混凝土重力坝深层滑动稳定性评估软件模块研发 单位： p a 图5 5 坝段最小主应力计算结果 F i g 5 Ca l c u l a t i o n r e s u l t o f mi n i mu m p r i n c i p a l s t r e s s f o r da m s e g m e n t 5 运行本程序， 完善各计算参数。点击 “ 开始计 算” 按钮 ， 按系统提示 进行操 作， 最终计 算得 K s = 1 2 4 1 ， 根据公式( 9 ) 判断可得坝体处于稳定状态； 而使用传统的刚体极限平衡法可得该坝段深层抗滑 稳定安全系数 =3 1 ， 同样满足规 范规定 的相应 工况下的安全系数要求 。虽然上述采用不同的分析 方法， 但是所得的结论是一样的。因此， 就本算例而 言 ， 使用本文提出的改进 的分项系数法对混凝土重 力坝的深层抗滑稳定性进行分析是可行 的， 其结果 是可信的。但是 ， 该方法仍然需要经过大量 的工程 实例的验证 ， 以确定其可行性和结果可靠性 。 6 结 语 ( 1 )随着有限元软件在水 利工程 中的广 泛应 用， 大量的后续处理工作需要各种软件相互结合 ， 从 而要求各软件之间的数据通信 。运用 A c t i v e X技术 就可 以解决数据通信问题。本文利用该技术对 A u t o C A D进行二次开发 ， 与重力坝深层抗滑稳定分析 相结合 ， 收到了很好的效果。结合有 限元 的计算成 果和 C A D的强大的绘图功能， 从而使水工计算程序 化 、 快速化 、 自动化 、 可视化。 ( 2 )基于 A c t i v e X技术 ， 可 以对支持该技术 的 软件， 如 w o r d ， e x c e l 等， 利用结构化面向对象的编程 语言 ， 并结合具体专业 问题 ， 进行软件 的二次开发 ， 从而对有限元计算结果的后续分析提供一种新 的、 快捷的、 自动化的、 实时的处理方法 。 ( 3 )结合本专业 的具体 问题对 A u t o C A D等支 持 A c t i v e X技术的软件进行二次 开发 的方法 ， 将 会 给水利工程的设计 、 分析 、 监测带来极大的方便 ， 在 这方面还有待进一步开发与研究。 参考文献： 1 李旭 重力坝抗滑稳定问题分析 J 四川水利， 2 0 0 7 ， ( 6 ) ： 3 O一3 3 ( L I X u A n t i s l i d i n g 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 7 ， ( 6 ) ： 3 0 3 3 ( i n C h i n e s e ) ) 2 苏怀智 ， 胡江， 温志萍， 等 基于时变风险率的大坝 使用寿命评估模型 J 中国科学 E辑： 技术科学， 2 0 0 9 ， 5 2( 7 ) ： 1 9 6 61 9 7 3 ( S U H u a i z h i ，HU J i a n g ， W E N Z h i p i n g，e t a 1 E v a l u a t i o n Mo d e l f o r S e r v i c e L i f e o f Da m B a s e d o n T i me Va r y i n g Ri s k P r o b abi l i t y l J 1 S c i e n c e i n C h i n a S e rie s E：T e c h n o l o g i c al S c i e n c e s ，2 0 0 9， 5 2 ( 7 )： 1 9 6 61 9 7 3 ( i n C h i n e s e ) ) 3 苏怀智 ， 温志萍， 吴中如 智能推理机在大坝健康早期 预警中的应用研究 J 水资源管理 ， 2 0 1 1 ， 2 5 ( 6 ) ： 1 5 4 51 5 6 3 ( S U Hu a i - z h i ， WE N Z h i p i n g ， wu Z h o n g - 1 3 1 S t u d y o n A n I n t e l l i g e n t I n f e r e n c e E n gin e i n E a r l y - Wa r n i n g S y s t e m o f Da m He alth I J J Wa t e r Re s o u r c e s M a n a g e m e n t ，2 0 1 1 ，2 5( 6) ：1 5 4 5 1 5 6 3 ( i n C h i - n e s e ) ) 4 段庆伟 ， 陈祖煜， 王杰 ， 等 重力坝抗滑稳定的强度 折减法探讨及应用 J 岩石力学与工程学报， 2 0 0 7 ， 2 6 ( 增 2 ) ： 4 5 1 04 5 1 6 ( D U A N Q i n g - w e i ，C H E N Z u y u，WANG J i e ，e t a 1 Di s c u s s i o n o n S t r e n g t h Re d u c t i o n Me t h o d for G r a v i t y Dam S t ab i l i t y An a l y s i s Ag a i n s t S l i d i n g a n d I t s A p p l i c a t i o n l J 1 C h i n e s e J o u r n al o f R o c k Me c h a n i c s a n d E n g i n e e ri n g ，2 0 0 7 ， 2 6( S u p 2) ： 4 5 1 0 4 5 1 6 ( i n C h i n e s e ) ) 5 张军涛， 姜礼斌， 唐春风基于有限元法的混凝土重力 坝抗滑稳定分析 J 四川建筑科学研究， 2 0 1 0 ， 3 6 ( 6 ) ： 1 3 21 3 4 ( Z H A N G J u n t a o ， J I A N G L i b i n ， T A N G C h u n f e n g Us i n g F E M An aly s i s o n S l i d e o f C o n c r e t e G r a v i t y D a m l J 1 S i c h u a n B u i l d i n g S c i e n c e ，2 0 1 0 ， 3 6 ( 6 ) ： 1 3 21 3 4 ( i n C h i n e s e ) ) 6 张彬洪 ， 姚激， 李泽， 等 基于有限元法的混凝土 重力坝抗滑稳定分析方法 J 水利与建筑工程学报， 2 0 1 0 ， 8 ( 1 ) ： 6 26 4 ( Z H A N G B i n h o n g ，Y A O J i ，L I Z e，e t a 1 An aly s i s o n S t ab i l i t y A g a i n s t S l i d i n g o f C o n - c r e t e G r a v i t y D am B a s e d o n F i n i t e E l e m e n t Me t h o d J J o u r n al o f Wa t e r R e s o u r c e s and A r c h i t e c t u r al E n gi n e e r - i n g ， 2 0 1 0 ， 8 ( 1 ) ： 6 2 6 4 ( i n C h i n e s e ) ) 7 蒋春艳 ， 常晓林， 周伟 用于重力坝抗滑稳定分析的 分项系数有限元方法 J 水力发 电学报， 2 0 0 6 ， 2 5 ( 2) ： 1 62 0 ( J I A N G C h u n y a n ，C HA N G X i a o l i n ， Z HOU W e i P a r t i al C o e ff i c i e n t F E M for Gr a v i t y Da m An t i s l i d i n g S t abi l i t y An aly s i s J J o u r n al o f Hy d r o e l e c t ri c E n g i n e e ri n g ， 2 0 0 6， 2 5 ( 2 ) ： 1 6 2 0 ( i n C h i n e s e ) ) 8 汪亚超 混凝土重力坝抗滑稳定多维时空分析方法研 究 D 南 京： 河 海 大 学， 2 0 1 1 ( WA N G Y a c h a o Mu l t i d i me n s i o n al S p a t i o T e mp o r a l An aly s i s Me t h o d for t h e A n t i s l i d i n g S t a b i l i t y o f C o n c r e t e G r a v i t y D a m D N a n j i n g ： H o h a i U n i v e r s i t y ， 2 0 1 1 ( i n C h i n e s e ) ) 9 陈世民理论力学简明教程 M 北京 ： 高等教育出 版社 ， 2 0 0 1 ( C HE N S h i m i n B ri e f T u t o ri al o f T h e o r e t i c al M e c h a n i c s M B e i j i n g ：Hi g h e r E d u c a t i o n P r e s s 长江科 学院院报 2 0 1 3芷 2 0 0 1 ( i n C h i n e s e ) 戴春来， 陈功 ， 陆信 A u t o C A D A u t o ma t i o n技术的 应用 J 工程设计 C A D与智能建筑， 1 9 9 9 ， ( 2 ) ： 5 6 5 8 ( D A I C h u n l a i ，C H E N G o n g ，L U X i n A p p l i c a t i o n o f A u t o C A D A u t o m a t i o n J C o m p u t e r A i d e d D e s i g n a n d I n t e l l i g e n t B u i l d i n g ， 1 9 9 9 ， ( 2 ) ： 5 6 5 8 ( i n C h i n e s e ) ) 周强 基于 A c t i v e X对象模型树的参数化绘 图技术 的研究与应用 J 南京化工大学学报， 2 0 0 1 ， 2 3 ( 5 ) ： 8 38 6 f Z HO U Q i a n g R e s e a r c h a n d A p p l i c a t i o n o f Gr a ph Pa r a me t r i z i n g Te c h n i q u e Ba s e d o n Ac t i v e X Mo de l T r e e l J I J o u r n a l o f N a n j i n g U n i v e r s i t y o f C h e m i c a l T e c h n o l o g y ， 2 0 0 1 ， 2 3 ( 5 ) ： 8 38 6 ( i n C h i n e s e ) ) 王玉用 V B A开发 A u t o C A D 2 0 0 0应用程序 M 北京 ： 人 民 邮 电 出版 社 ， 1 9 9 9 ( WA N G Y u D e v e l o p me n t o f Au t o C AD 2 0 0 0 Ap p l i c a t i o n s b y Us i n g VB A f M 1 B e i j i n g ： P o s t sT e l e c o m P r e s s ，1 9 9 9 ( i n C h i n e s e ) ) 1 3 赵峰， 毕苏萍， 周振红D e l p h i 通过双接 口控制 A u t o C A D绘 图 J 郑 州 大 学 学 报 ( 工 学 版 ) ， 2 0 0 7 ， 2 8 ( 4 ) ： 6 66 9 ( Z H A O F e n g ，B I S u p i n g ，Z H O U Z h e n h o n g De l p h i Dr i v i ng Au t o CAD t o Dr a w b y Wa y o f Du a l I n t e r f a c e s I J I J o u r n a l o f Z h e n g z h o u U n i v e r s i t y( E n g i n e e r i n g S c i e n c e ) ， 2 0 0 7， 2 8 ( 4 ) ： 6 66 9 ( i n C h i n e s e ) ) 1 4 李兵 ， 赵锡英D e l p h i 在 C A D系统中的应用 J 兰 州工业高等专科学校学报， 2 0 0 4 ， 1 1 ( 4) ： 1 61 9 ( I J I Bi ng，ZHAO Xi y i n g Th e Ap pl i c a t i o n o f De l p h i i n CAD S y s t e m J J o u r n a l o f L a n z h o u P o l y t e c h n i c C o l l e g e ， 2 0 0 4 ， 1 l ( 4 ) ： 1 61 9 ( i n C h i n e s e ) ) ( 编辑 ： 姜小兰 ) S o f t wa r e M o d u l e t o As s e s s t h e De e p S l i d i n g S t a b i l i t y o f Co n c r e t e Gr a v i t y Da m CUI S h u s h e n g ，S U Hu a i z h i ，ZHANG S h u a i ，GU Ha o q i ng ( 1 S t a t e K e y L a b o r a t o r y o f H y d r o l o g y - Wa t e r R e s o u r c e s a n d H y d r a u l i c E n g i n e e ri n g ，H o h a i U n i v e r s i t y ， N a n j i n g 21 0 09 8，Chi n a； 2 Na t i o n a l En g i n e e r i n g Re s e a r c h Ce n t e r o f W a t e r Re s o u r c e s Effi c i e n t Ut i l i z a t i o n a nd E n g i n e e r i n g S a f e t y ， H o h a i U n i v e r s i t y ， N a n j i n g 2 1 0 0 9 8 ， C h i n a ； 3 C o l l e g e o f Wa t e r C o n s e r v a n c y a n d H y d r o p o w e r ， H o h a i U n i v e r s i t y ， N a n j i n g 2 1 0 0 9 8 ， C h i n a ) Abs t r a c t ： To p r e c i s e l y r e fle c t t h e a nt i - s l i di n g s t a bi l i t y o f d a m ，a v i s u a l i z e d s o ftwa r e mo d u l e wa s r e s e a r c h e d a n d d e v e l o p e d t o a s s e s s t he d e e p s e a t e d s l i d i n g s t a b i l i t y o f c o n c r e t e g r a v i t y d a m Fi r s t l y，t h e s t r e s s d i s t rib u t i o n o n t h e d e e p s e a t e d s l i di n g p l a n e o f t h e d a m wa s n u me r i c a l l y c a l c u l a t e dOn t h i s ba s i s ，b y u s i n g t he d e s i g n me t ho d o f t he l i mi t s t a t e o f pa r t i a l c o e ffi c i e n t ，t h e f o r mu l a o f s h e a r s t r e n g t h wa s a p p l i e d t o e a c h F EM e l e me n t t o d e fin e a n e w f o r mu l a f o r t h e s a f e t y c o e ffic i e n t ，t h us a n i mp r o v e d me t h o d o f p a r t i a l c o e f fic i e n t f o r t h e s t a b i l i t y a s s e s s me n t wa s d e v e l o p e dI n s u b s e q u e n c e，t h e me t h o d wa s p r o g r a mme d b y p r o g r a mmi n g t o o l De l p h i i n a s s o c i a t i o n wi t h Au t o CAD S i n c e t h e a d v a n t a g e o f FEM i n c a l c u l a t i n g s t r e s s fi l e d a n d c o n s i d e rin g d a m f o un d a t i o n a n d b o d y d e f o r ma t i o n i s r e fl e e t e d i n t h e me t ho d a n d t h e s o f t wa r e mo d ul e，a n d o wi n g t o t h e r e l i a b i l i t y d e g r e e t a k e n i n t o a c c o u n t a s we l l ，t h e o b t a i ne d s a f e t y c o e ffi c i e n t c o u l d we l l r e