混凝土重力坝深层滑动稳定性评估软件模块研发.pdf
《混凝土重力坝深层滑动稳定性评估软件模块研发.pdf》由会员分享,可在线阅读,更多相关《混凝土重力坝深层滑动稳定性评估软件模块研发.pdf(6页珍藏版)》请在咨信网上搜索。
1、第 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
2、 水文水资源与水利工程科学国家重点实验室; b 水资源高效利用与工程安全国家工程研究中心; c 水利水电学院, 南京2 1 0 0 9 8 ) 摘要: 基于混凝土重力坝深层滑动面上应力分布情况的数值计算成果, 结合分项系数极限状态设计方法 , 研究了一 种评估混凝土重力坝深层组合滑动稳定性的改进分项系数法。在此基础上, 利用 D e l p h i 编程工具 , 结合 A u t o C A D 二次开发技术对该方法实现了程序化 , 得到了评估混凝土重力坝深层滑动稳定性的可视化软件模块。该方法和软 件模块, 充分利用了有限元数值模拟方法在计算应力场和考虑坝基及坝体变形影响等方面的优势, 又结合
3、了可靠 度设计方法, 得出的稳定安全系数能够较好地反映重力坝的抗滑稳定特性, 且沿用了设计规范中的 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 重 力坝深层抗滑稳定分析方法 建基面及坝基深层抗滑稳定分析 , 是重力坝设 计 中一项非常重要的内容 ; 并且在保证施工、 运行的 安全
4、性方面起着举足轻重的作用 ; 同时 , 准确分析坝 体 的抗滑稳定特性也是重力坝开展其他工作的基本 前提。 目前 重力 坝坝基深层抗 滑稳定分析方法 大致分为 3种 : 刚体极限平衡法 、 有限单元法和地质 力学模型法 。 随着计算力学和计算机软件及硬件的发展 , 利 用有限单元法可以比较准确地计算出坝体及其坝基 的应变 、 应力和塑性区。依靠计算机 , 有限单元法能 够精确建立模型 , 并考虑各种材料的非均质 、 各 向异 性等特性 , 因而能较准确地计算接触面 以及深层组 合滑动面上的阻滑力和滑动力 。J 。然而仅仅计算 出阻滑力和滑动力还不能判断大 坝的抗滑稳定性 状 , 需要找到一种计
5、算方法 , 以较为准确地反映大坝 的抗滑稳定特性。本文首先利用有限元法准确计算 出重力坝深层滑动面上的阻滑力和滑动力 , 然后结 合分项系数极限状态设计方法 , 对每一个剖分单元 运用抗剪断强度公式 , 定义一种新 的抗滑稳定安全 系数 的计算公式 , 得到了一种计算混凝土重力坝 深层组合滑 动 的改进分 项系数法 。并采用 D e l p h i 语言, 结合 A u t o C A D二次开发技术对该方法实现了 程序化。该方法充分利用 了有限元能够计算出精确 应力场和考虑坝基及坝体的变形影响的优 势, 又考 虑了可靠度设计方法 , 所以得出的 K 能够较好地反 映重力坝 的抗滑稳定特性 ,
6、 并且本方法沿用 了设计 规范 中的 K 表示方法 , 易于计算和普及。 2 混凝 土 重 力坝 深 层 组 合 滑 动 的 改 进分项 系数法 本方法 的基本思想 是 : 引入 广义力 的概念 , 将传统的作用效应函数 S ( ) 和结构及构件抗力函 数 R ( ) 广义化 , 结合混凝土重力坝的实际地质情 况 , 分析深层滑移模式 , 从而决定这种广义力以及该 广义力对应的作用效应函数和抗力函数。 下面以平面状 2滑面和3滑面( 如图 1 所示) 为 例 , 简述混凝土重力坝深层组合滑动的改进分项 系 数法 。根据混凝土承载能力极 限状态 的基 本组合 , 设计表达式为 1 1 , 、 Y
7、 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 ) ;
8、中央高校基本科研业务费专项资金资助( 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 ()
9、为作用效应函数 ; ( ) 为结构及构件抗力函 数 ; 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()
10、 , R( ) 进行广义化。事实上 , 当大坝及坝基沿着组 合滑动面出现抗滑失稳的时候 , 滑移体 ( 由坝体及 滑动面以上 的岩体组成 ) 受到边 界的约束 , 将绕着 滑移体瞬时滑移 中心做转动。在此 , 将阻滑力和抗 滑力广义为阻滑力偶和抗滑力偶。 具体计算方法及步骤 : 对于其中的第 1 个滑动 面, 首先根据有限单元法 的应力计算成果 , 利用式 ( 2 ) 和式 ( 3 ) 将滑动面相关单元 的应力合成为法 向 力 和切向力 Q。 , 同理第 2个滑动面为法向力 2 和切向力 Q ; 除了确定力的大小外 , 还需确定力 的 作用点。将滑动面上的各单元法向力看成是作用在 该滑动面上的
11、平行力 系, 利用求解平行力系中心 的 方法可以求得作用点 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 沿滑动面法 向 的正应力 ; 为组合滑动面上单元
12、 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
13、) 所示 , 对于 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+
14、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滑面的抗滑稳 定安全系数。 比较式
15、 ( 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 m
16、e 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:
17、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步: 这一步是整个软件的核心, 即找出有限 元网格中和组合滑动面相交的网格 , 本文称之为相 关单元。当找出所有和组合滑动面相交的相关单元 之后 , 就可以通过这些单元信息 , 从有限元计算软件 的所计算的应力场导出和组合滑动面相关的单元应 力信息 , 最后读取
18、这些 应力 , 通过弹性力学知识 , 换 算成滑动面上的正应力和切应力。然后和单元面积 相乘 , 转化为我们想要知道的滑动力和抗滑力。最 后通过改进的组合系数法 , 求出抗滑稳定安全系数 , 判断重力坝深层抗滑是否稳定。 寻找相关单元的核心思想是 : 通过 A u t o C A D给 出的对象模 型树 中提供 的 I n t e r s e c t Wi t h方法 , 遍历 每个单元的每条边 , 求出组合滑动面和单元 的边的 交点 , 确定个数 , 如果个数不为零, 则确定该单元就 是相关单元 , 并且记下滑动面和该单元的交点, 求 出 2交点之间的距离 d和滑动面的方向余弦 f , m,
19、 为之 后求解抗滑力和下滑力做准备 , 具体代码这里不再 赘述 。 在找出相关单元并且把相关单元的信息 , 包括 单元号 、 直线号 、 滑动面和相关单元 的交点读入 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 +
20、, 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 个坝段组成。其中,
21、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 应用本文所述 的改进分项系数法验算该重力坝 深
22、层抗滑稳定。首先建立该坝段的有限元分析数值 模型。为了简化实例计算 , 将该坝段视为平面应变 问题 , 利用平面三角形和 四边形等参单元模拟。考 虑该坝段及其基础条件 的实际情况 , 选定其有 限元 计算范围为 : 上游约 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
23、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
24、 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 , 同样满足规 范规定 的相应 工况下的安全系数要求 。虽然上述采用不同的分析 方法, 但是所得的结论是一样的。因此, 就本算例而 言 , 使用本文提出的改进 的分项系数法对混凝土重 力坝的深层抗滑稳定性进行分析是可行 的, 其结果 是可信的。但是 , 该方法仍然需要经过大
25、量 的工程 实例的验证 , 以确定其可行性和结果可靠性 。 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技术 , 可 以对支持该技术 的 软件, 如
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 混凝土 重力坝 深层 滑动 稳定性 评估 软件 模块 研发
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【fir****x8】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【fir****x8】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。