混凝土结构孔隙侵蚀数值模拟.pdf
《混凝土结构孔隙侵蚀数值模拟.pdf》由会员分享,可在线阅读,更多相关《混凝土结构孔隙侵蚀数值模拟.pdf(4页珍藏版)》请在咨信网上搜索。
1、第 1 2 卷第 5 期 2 0 1 4年 1 O月 南水北调与 水利科技 S o u t h - t o - No r t h Wa t e r Tr a n s f e r s a n d Wa t e r S c i e n c e& Te c h n o l o g y Vol _ 1 2 No 5 Oc t 2 O 1 4 D OI : 1 0 1 3 4 7 6 j c n k i n s b d q k 2 0 1 4 0 5 0 O 8 混凝土结构 孔隙侵蚀数值模 拟 葛金金 , 詹美礼 , 吴金萍 , 盛金 昌 , 张 羽。 ( 1 河海大学 水利水电学院, 南京 2 1 0
2、 0 9 8 ; 2 长江水利委员会水文局 长江下游水文水资源勘测局, 南京 2 1 0 0 1 1 ; 3 广东省水利电力勘测设计研究院, 广州 5 1 0 6 3 5 ) 摘要: 混凝土材料所具有的多孔隙结构是由其本身的物理、 化学性质所决定的, 而孔隙结构又决定了混凝土整体的 强度和渗透性 。将混凝土复杂的孔隙结构简化为单根毛细管 , 根据 F i c k定律和质量守恒原理推导出混凝土孔隙溶 解控制方程, 并运用 B e a r 理论中关于多孔介质毛细孔渗透率的研究理论及饱和溶质运移扩散方程建立侵蚀模型。 模拟结果表明, 该模型能够有效地反映混凝土孔隙中固相钙离子溶解、 浓度变化和孔隙边
3、壁侵蚀厚度变化过程。 关键词: 混凝土; 孔隙结构; 侵蚀模型; 数值模拟 中图分类号 : T U3 7 5 文献标 志码 : A 文章 编号 : 1 6 7 2 1 6 8 3 ( 2 0 1 4 ) 0 5 0 0 3 2 0 4 Nu m e r i c a l s i m u l a t i o n o f po r e e r o s i o n i n c o n c r e t e s t r u c t ur e G E J i n - j i n , Z HAN Me i l i , WU J i n - p i n g , S HE NG J i n - c h a n g
4、 , Z HANG Yu 。 ( 1 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 Hy d r o p o we r E n g i n e e r i n g, Ho h a i U n i v e r s i t y, Na i n g 2 1 0 0 9 8 , C h i n a ; 2 L o we r Re a c h e s o J t h e C h a n g J i a n g Hy d r o l o g y a n d Wa t e r R e s o u r c e s S u r v e y B u
5、r e a u , C h a n g j i a n g Wa t e r Re s o u r c e s C o mmi s s i o n Bu r e au , Na i n g 2 1 0 0 1 1 , Ch i n a; 3 C man g d o n g Hy d r o p o we r Pl a n n i n g a n d De s i g n I n s t i t u t e, Gu a n g z h o u 5 1 0 6 3 5, C h i n a) Ab s t r a c t : Th e c o n c r e t e ma t e r i a l
6、h a s a p o r o u s s t r u c t u r e d u e t O i t s o wn p h y s i c a l a n d c h e mi c a l p r o p e r t i e s , a n d t h e p o r e s t r u c t u r e d e t e r mi n e s t h e s t r e n g t h a n d p e r me a b i l i t y o f c o n c r e t e I n t h i s p a p e r , t h e c o mp l e x p o r o u s
7、 s t r u c t u r e o f c o n c r e t e wa s s i mp l i f i e d t O t h e s i n g l e c a p i l l a r y Ac c o r d i n g t o t h e Fi c k s l a w a n d ma s s c o n s e r v a t i o n p r i n c i p l e , t h e d i s s o l u t i o n e q u a t i o n o f c o n c r e t e p o r e s wa s o b t a i n e d , a
8、 n d t h e e r o s i o n mo d e l wa s e s t a b l i s h e d b a s e d o n t h e r e s e a r c h o f c a p i l l a r y p e r me a b i l i t y o f p o r o u s me d i u m i n Be a r t h e o r y a n d s a t u r a t e d s o l u t e t r a n s p o r t d i f f u s i o n e q u a t i o n s Th e s i mu l a t
9、i o n r e s u l t s s h o we d t h a t t h e mo d e l c a n r e f l e c t t h e v a r i a t i o n s o f d i s s o l u t i o n a n d c o n c e n t r a t i o n o f s o l i d p h a s e c a l c i u m i o n s a n d e r o s i o n t h i c k n e s s o f s i d e wa l l i n p o r e s Ke y wo r d s : c o n c r
10、 e t e ; p o r o u s s t r u c t u r e ; e r o s i o n mo d e l ; n u me r i c a l s i mu l a t i o n 混凝土结构的侵蚀是指在流动的软水或具有侵蚀性的 溶液作用下, 混凝土孔隙边壁固相水化物钙离子 自动剥落或 与侵蚀性水溶液中的其它离子发生交换、 扩散、 迁移, 从而造 成了钙离子从混凝土结构中缓慢流失的现象, 这是坝体、 坝 基破坏 的重要原 因之一 l 1 。 。因此研 究混 凝土 侵蚀 机理 , 建 立混凝土孔隙侵蚀溶解模型, 预测混凝土侵蚀程度, 具有重 要的理论和实际意义。 目前 ,
11、模拟混凝 土侵蚀 过程 的数值 模 型 比较 少 。其 中 , G6 r a r d 模型L 4 是由热力学理论发展而来的, 其基本原理如 下 : 混凝土内部孔隙固相与液相之间存在的浓度梯度, 在液 相石灰浓度未能达到各自水化物的极限浓度的情况下 , 该浓 度梯度促使混凝土固相水化物( 氢氧化钙、 水化硅酸钙等) 中 钙离子溶解并向外扩散, 因此根据溶出液中钙离子的浓度变 化就可以判断混凝土受侵蚀的程度。G 6 r a r d模型的缺点是 未能考虑混凝土孔隙液相中钙离子浓度 的动态平衡。D e t l dE 等人针对混凝土水泥基质是不均匀混合体的性质, 发 展了一种关于骨架孔隙溶液中的钙离子浓
12、度平滑函数 , 用于 模拟对孔隙框架内基质扩散率、 孔隙率 以及表面扩散率等控 制性参数, 其局限性是没有考虑到对流项的传递, 而是只考 虑孔隙 中的扩散问题 。 本文在分析混凝土渗流侵蚀机理的基础上, 避开以往对 混凝土表面侵蚀的研究方式 6 , 分别建立孔隙边壁溶解方 程、 溶质运移扩散方程和渗流方程, 最后形成混凝土侵蚀溶 解、 扩散控制模型, 模拟并分析混凝土孔隙中固相钙离子溶 解、 浓度变化和孔隙边壁侵蚀厚度变化规律。 1 混凝土渗透侵蚀模型 混凝土作为多孔介质, 虽然具有较复杂的孔隙结构 , 但 收稿 日期 : 2 0 1 4 0 3 0 4 修回 日期 : 2 0 1 4 0 7
13、 2 8 网络 出版时 间: 2 0 1 4 0 8 2 7 网络 出版地址 : h t t p : ww w c n k i n e t k c ms d o i 1 0 1 3 4 7 6 j c n k i n s b d q k 2 0 1 4 0 5 0 0 8 h t ml 基金项 目: 国家 自然科学基金( 5 1 0 7 9 0 3 9 ; 5 1 0 0 9 0 5 3 ) 作者简介: 葛金金( 1 9 8 8 一 ) , 女, 安徽合肥人, 主要从事岩土渗透破坏与边坡稳定方面的研究。E - ma i l : 9 2 4 0 8 5 6 7 0 q q c o rn 娩 |
14、试 验 豫 灸 葛金金等 混凝土结构孔隙侵蚀数值模拟 对其内部结构被溶液侵蚀的机理可以简化为一个单根毛细 孔隙的侵蚀过程_ 8 。由于毛细管壁在化学渗透水流的作用 下发生了侵蚀溶解, 同时侵蚀析出的钙离子在对流、 扩散的作 用下运移, 使得毛细管中钙离子浓度发生变化, 毛细管中不同 位置钙离子的浓度会影 响该处水化物进一步侵蚀的速度 。 1 1 孔 隙边 壁溶 解 方程 混凝土结构 中的毛细 孔主要 是 由于混 凝 土 中的水 泥基 质经过水解作用引起的。采用 B e a r 理论所提到的直平行毛 细管结构对毛细孔结构进行模拟, 并做出以下假设 。 ( 1 ) 毛细管道中的渗流为层流。 ( 2
15、 ) 管壁水化物为均质, 溶解特性为各向同性。 ( 3 ) 毛细管道中溶解后所产生的不可溶杂质不会沉淀堆 积在管道中堵塞管道。 ( 4 ) 固相水化物的溶解过程相对于扩散过程时问很短, 忽 略钙离子的溶解时 间。 ( 5 ) 不可溶质对离子扩散运移的影响忽略不计 。 ( 6 ) 忽略水泥水化温度的变化。 取一根毛细管( 图 1 ) , 沿毛细管水流方向上取坐标 X, 管 道半径为 r , 横截面管壁的侵蚀变化量为 d r , 根据费克( F i c k ) 扩散定律 1 0 1 1 有 J= -D, 8 C d 式中 : 为溶质 的扩散通量 ; D 为 固相溶解 扩散 系数 ; C为液 相浓度
16、。 图 1 毛 细管壁侵 蚀溶解示意 图 Fig 1 S c h e ma t i c d i a g r a m o f d i s l u t io n a n d e r o s i o n i n c a p i l l a r y wa l l s 根据物质质量平衡原理, 对于微原体 r d r的溶 解有 : 移带走, 从而使得 处的液相钙离子浓度不能达到石灰的极 限浓度。钙离子在孔隙溶液中的扩散、 运移符合多孔介质中 的对流扩散方程 : 一 d i v ( D VC ) 一VVC+f ( 4 ) 式中: C为孔隙中液相钙离子浓度; D为水动力弥散系数, 包 括分子扩散系数 和机械弥
17、散 系数; V为孔隙中流速矢量; f 为源项 。 1 2 1 弥散系数确定 多孔介质中的分子扩散系数 取决于溶液中的分子扩 散系数 D 和多孔基质的弯曲率 T , 即: 一Dd T 。 机械弥散系数 D 取决于宏观渗透速度 。在各向同性 的介质 中机械弥散系数是二 阶对称 张量 , 对 于 x o y平面 内的 二维流为 D 一 ( a +口 ) D 一 ( n T +口 L ) D 一口 T V ( 5 ) D 一吉( L n T ) V s D 一0 D 一O ( 6 ) 式中: a 、 a 分别为纵向和横向弥散度; 为分子扩散系数; 为平均水流速度 。 所以, 水动力弥散系数为 。 一
18、c 1 2 2源项的确定 设 F为固相表面的水化物浓度; f ( c , F) 为单位时间从 单位体积固相中进入液相的水化物质量, 则单位时间内溶解 NNNNN: _厂 ( c , F) , 由 d r D M ( GC =f ) ,得 出 时间 内毛细管道 所溶解的摩尔数为 m一 一2 rD ( c c)C L cd f( 8 ) 单位体积溶液中, 单位时间溶入的溶质质量为 , d r。 r d 一 d r r d r ( 2 , : 式中: C L r 为沿孔隙管道的长度 ; 为水泥水化物 的平均密 度 ; M为水泥水化物的摩尔质量 ; d r为溶解边壁的厚度。其 余符号 含义与式 ( 1
19、 ) 相 同。 联立式( 1 ) 、 式( 2 ) , 并引进薄膜理论得单根孔 隙侵蚀模 型 的溶 解方 程为 d r Df ( C C, ) ( 3 ) u _ D 式 中 : C 为孔 隙液相 中钙离子 的极 限浓度 ; C _f 为 t时刻管 道 处的液相钙离子浓度。其 中, 固相水化物的扩散系数 D 采用 G 6 r a r d实验 中测定 的与孔 隙溶液浓度 有关 的参 数 ( 2 0 5 6 0 1 0 4 C ) 1 0 。 m2 s计算。其余符号含义与式 ( 2 ) 相 同。 1 2 溶质 运 移扩散 方 程 毛细管道边壁水化物溶解后, 液相钙离子浓度会逐渐增 加。但是, 在有
20、渗透水流的作用下 , 钙离子会被水流不断运 : 垒 二 : 堂 : : 二Q m 0d x 出 4 9 ,( _ C - - 一C)( 9 ) 1 3 渗 流方 程 根据 B e a r 理论 , 假设 混凝 土 多孔介 质 中的孔 隙为 许多 平行直毛细管组成, 渗透系数为 Kgn c G ( 1 0 ) 式中: K为多孔介质的渗透系数 ; g为重力加速度; 为水的 运动黏质 系数 , 取 1 0 1 0 m 2 s ; d T为毛细 管道 z处 的直 径; 为孔隙率。 渗流连续性方程的微分形式为 + 警一 一 d iv (1D ( 1 1 ) 式中: p 、 是宏观量, 分别为流体密度和流
21、速。 孑 L 隙介质骨架不可压, 但孔隙介质结构因侵蚀而变化, 且流体为等温均质不可压时, 与时间 有关, 即 O n o 、lD g| 蠢 3 3 第 1 2卷 总第 7 4期 南水北调与水利科技 2 0 1 4年第 5期 C时, 则式( 1 1 ) 可变为 姿 丝+ + 一 一 一O n ( 1 2 ) 。 O y。a Ot 、 在以水头为因变量的微分形式中, 取坐标轴方向与介质 各 向异性 的主方 向一致 , 达西定律的表达式为 f 一一K l l , ,aH 一 一 K O y ( 1 3 ) l 一 一 K 当介质为各向同性 K 一 一K K时, 将式( 1 3 ) 带 入式 ( 1
22、 2 ) 得 暑 ( K ) + ( K ) + 未 (K ) 一 ( 1 4 ) 对于x o y平面内的二维渗流问题, 将式( i 0 ) 代入上式得 ( D H ) -H O( OH )O n 3 2 v 0 y 3 2 v O y O t ) 、 细管 内壁 溶解 后, 内径增 加 到 一 + d r , 又 ”一 NO r 8 2 4 ) L L=N 4 , 式 ( 1 5 ) 变为 ( 3 2 v ) +0 y( 3 2 v 0 y ) 一 0 t ( 1 6 ) 、 、 联立方程( 3 ) 、 ( 4 ) 和( 1 6 ) 得到混凝土毛细管道边壁侵蚀 溶解 、 扩散控制方程如下:
23、一 M ( G _C) ac d i v ( D c ) 一 vc + : : 4 D, ( C - C) a r O H、 +旦 r 、32 v av、32 v0y af 一 S _ _ d r 2 混凝土侵蚀数值模 拟 2 1 模拟参数与条件 采用多物理场耦合有限元软件 C OMS OI , 对混凝土在 渗透溶解以及浓度扩散、 运移过程中的钙离子浓度变化和毛 细孔隙壁的溶解进行模拟。根据实验测定的初始渗透系数 , 利用 Ko z e n y - C a r ma n t ” 妇公式计算出孔隙直径大小, 并以此 孔径作为初始条件, 具体参数详见表 1 , 表 1 孔隙侵蚀模 拟参数 Ta b
24、 1 Po r e e r o s i o n s i m u l a t i o n p a r a me t e r s 参数 值 参数 值 孔径 d m 7 3 e - 7 R T m 0 0 0 2 孔 长 m 0 2 固 相 沓 化 物 质 量 7 4 ( gm o l ) 水 化物密度 。 ( k g m。 ) 。 u u 黏滞系数 1 0 1 0 a L m 0 0 0 5 I ) d ( m S ) 1 5 e - l 1 分别给定边界条件和初始条件如下。 ( 1 ) 边界条件 。 H 一2 0 0 m, HI S 一0 ; 孔隙边壁为不透水边界。 3 4 试验研究 l S 1
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 混凝土结构 孔隙 侵蚀 数值 模拟
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【k****o】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【k****o】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。