2016院报1期（4校）3月8日

Size: px

Start display at page:

Download "2016院报1期（4校）3月8日"

换徐
7 years ago
Views:

1 第 14 卷第 1 期中国水利水电科学研究院学报 Vol.14 No 年 1 月 Journal o China Institute o Water Resources and Hydropower Research February,2016 文章编号 (2016) 组合网格法在无压稳定渗流中的应用 1,2 1,2 张顺福, 丁留谦 (1. 中国水利水电科学研究院防汛抗旱减灾研究所, 北京 ;2. 水利部防洪抗旱减灾中心, 北京 ) 摘要针对无压稳定渗流分析过程中常遇到的线形结构问题, 提出将组合网格法应用于无压稳定渗流计算采用两套独立的网格, 对整体区域采用尺寸较大的粗网格进行模拟, 不考虑线形结构的影响, 在线形结构附近则采用尺寸比较小的细网格进行模拟, 考虑线形结构的影响计算在粗细两套网格之间迭代调整, 直到达到收敛精度组合网格法大大简化离散工作量, 且适应于非规则网格数值计算结果表明, 组合网格法合理可行, 计算程序可靠, 这为线形结构的无压稳定渗流分析提供了一种思路和方法关键词无压稳定渗流 ; 组合网格 ; 线形结构 ; 改进初流量法 ; 有限单元法中图分类号 TV 文献标识码 A doi /j.cnki.jiwhr 研究背景在饱和无压稳定渗流有限元数值模拟中, 常常会遇到这样一些结构, 它们在某一方向的尺寸与整体计算区域相比非常小, 比如岩体中大量存在的断层蚀变带软弱夹层等结构面 ; 土石坝中常用的土工膜防渗面板防渗心墙防渗帷幕等局部构造等等 ( 为方便描述, 下文中统称这些结构为 [1] 线形结构或奇异区域 ) 当这些结构中的水流运动可以忽略不计的时候, 可以采用结点双编号或内 [2] 部隔水边界对这些结构进行处理, 否则对这些结构需要进行精细离散, 但这些结构在某一方向的尺寸可能达到厘米甚至毫米级别, 与整体模拟区域的尺寸相差太远, 根据有限元误差理论, 为保持计算精度其它区域必然也需要精细的离散, 从而导致离散所需结点数目和单元数目将会大大增加, 且考虑到结构的空间特性, 整体计算区域的离散难度将会大幅度提升, 计算所需时间也大大增加, 甚至到达一般计算机难以承受的程度对于这类局部特性强烈的线形结构, 一般是采用非规则网格局部加密的方法对这些结构及其附近区域进行局部的离散加密, 但局部加密网格法的要求比较高, 不仅要求整个计算区域的网格质量很好, 而且还要求网格的尺寸过渡非常光滑, 这些条件在二维情况下尚且能够很好的满足, 但在三维情况下, 对复杂区域的非规则网格划分常常遇到下述问题直接划分失败, 划分成功但单元质量差, 或者划分成功但单元尺寸过渡不光滑等为避免局部加密所遇到的问题, 一种可行的方案是采用非拟一致网格 (Non-conorming Grid) 技术, 非拟一致网格技术的主要思想是将整体区域分解为相互之间只在区域表面有交集的不重叠小区域, 区域采用独立的网格离散, 区域之间不需要满足网格协调的要求, 区域之间的连续性条件通过多点约束技术 (Multiple-Point Constraints) 实现, 通过运用这种局部区域协调整体区域非协调的网格离散方式, 达到不同的局部区域上采用不同的网格离散或者不同的数值计算方法的目的, 也就是说 [3] 可以在某个区域采用有限元方法进行计算, 而在其它区域采用有限体积或者有限差分进行计算非拟一致网格技术很好的将局部特性控制在小范围内, 使局部区域的特性不至于导致整体计算的急剧增长但有效的适应局部特性目标往往与计算求解过程相冲突, 包括 (1) 由于离散尺寸的差收稿日期基金项目中国科技部国际合作项目 (2010DFA74520) 作者简介张顺福 (1982-), 男, 广西人, 博士, 工程师, 主要从事渗流与控制研究 zhangshunu@live.cn 16

2 异变化而导致方程求解器会降低效率甚至求解失败 ;(2) 由于考虑了非结构网格, 数据结构变得繁琐 ;(3) 计算机体系结构可能会不能够高效的处理非规则网格的影响 ( 比如, 向量化被抑制 );(4) 即使是在离散过程中, 这个目标也会带来重重困难, 对于有限差分格式而言, 很难为非规则网格构造出相应的高精度差分格式 ; 对于有限元计算而言, 一组质量好的非结构化网格的自动生成意味着高额的花销 [4-6] 为避免非拟一致网格的困境,Mc Cormick [6] 提出了一种有叠加的区域分解法快速自适应组合网格 (FAC) 方法 FAC 方法首先形成一套对应于整体区域的规则粗网格, 然后对规则粗网格上相应于关注的局部区域形成一套或者多套与整体粗网格嵌套的细网格, 计算的时候先对粗网格求解, 然后将粗网格的计算结果通过延拓算子延伸到细网格上, 进而获得细网格上的误差, 将这些误差再通过 [7-8] 限定算子返回粗网格上, 粗细网格之间循环迭代, 反复对误差进行校正直至收敛由于规则网的处理 ( 离散与求解 ) 相对不规则网而言要容易得多, 从而使得快速自适应组合网格 (FAC) 方法成为非常有效的方法, 但是快速自适应组合网格 (FAC) 方法高效性有 2 个基本要求 (1) 嵌套和结构化或规则的网格又往往成为其应用的限制条件, 由于工程实际问题中区域的复杂性使得规则网格和嵌套的规则网格往往很难甚至无法得到 ;(2)FAC 方法相关的数值计算格式以嵌套为基础, 如果不满足嵌套条件, 那么几乎所有相关的计算都无从下手这些均使得 FAC 的应用收到了限制为此, 本文引入一种改进的快速自适应组合网格法组合网格法 (Composite Grid Method, CGM) CGM 继承了 FAC 反复在粗细两套网格进行误差校正的思想, 并且可以使用非规则网格, 细 [9] 网格和粗网格相互独立, 不要求网格嵌套组合网络法, 该方法已广泛应用于水工结构油田开采数 [10-11] [12] [13] [14-16] 值模拟电磁场电磁机械耦合焊接数值模拟等领域 2 无压稳定渗流基本理论 [17] 根据微分体质量守恒定律, 在渗流区域 D, 三维地下水稳定渗流微分方程为 - 3 v i + W = 0 (,2,3 ) (1) x i 式中 x i 为笛卡尔直角坐标系的坐标轴 ;v i 为沿坐标轴 x i 方向的渗流速度,m/s;W 为源或汇项,1/s, 以获得水流为正同时满足 Darcy 定律 3 v i = - k ij (,2,3 ) (2) x j 式中 k ij 为渗透系数 ;H 为总水头,m,H = x 3 + p γ ;p 为孔隙水压强,Pa;γ 为水的容重,N/m 3 定解条件如下 ( 如图 1 所示 ) (1) 定水头边界 Γ 1 H Γ1 = H 0 ( x 1,x 2,x 3 ) (3) (2) 流量边界 Γ 2 ( 以流入为正 ) æ ç k ij ø (3) 自由面边界 Γ 3 æ H Γ3 = x 3, 且 ç x 1,x 2,x 3 (4) x j n i Γ 2 = q ( ) k ij n x i = 0 (5) j ø Γ 3 图 1 渗流计算边界条件 17

种有叠加的区域分解法快速自适应组合网格 (FAC) 方法 FAC 方法首先形成一套对应于整体区域的规则粗网格, 然后对规则粗网格上相应于关注的局部区域形成一套或者多套与整体粗网格嵌套的细网格, 计算的时候先对粗网格求解,

3 (4) 渗流逸出边界 Γ 4 æ H Γ = x 3, 且 4 ç k ij n x i < 0 (6) j ø Γ 4 式中 H 0 为已知水头函数,m;q 为已知流量,m 3 /s;n i 为边界外法线方向余弦求解 [18] 无压稳定渗流控制方程属于椭圆型偏微分方程, 本文采用基于子域积分改进的初流量法进行 3 组合网格法原理 3.1 组合网格法基本思想组合网格法是一种有叠加的区域分解算法, 其基本思想是采用粗网格和细网格两套网格进行分析计算, 粗网格对应整个区域, 不考虑线形结构的影响, 细网格则为包含线形结构在内的局部区域两套网格之间反复迭代直至收敛, 网格之间的信息交互通过网格空间位置插值实现, 网格之间无需嵌套, 因此适应性 [19] 强, 规则网格和非规则网格均可使用组合网格法是文献 [20-21] 所提出的方法的改进从数学上证明了用两套网格所求得的数值解可以在粗网格区域达到粗网格的精度, 在细网格区域达到细网格的精度, 从理论上证明了组合网格法的有效性 3.2 算法描述对于图 2 所示的整体区域 = c ( 下标 c L c 图 2 L 计算区域 [18,22-24] c 分别代表细网格区域及细网格以外的区域 ), 考虑采用初流量法进行求解的三维无压稳定渗流问题 ì L c v = c c í (7) îl u = v Γ = u Γ 式中 u v 分别为细网格区域未知量及细网格以外区域的水头值 ;L p = 3 v 0 æ x ç k ij 为相应区域 i x j ø i 的偏微分算子 ; p = -W x p,v 0 [22] i 为初始流速,W p 为相应区域的源或汇项 ;Γ = c 为两个区 i 域的交界的边界采用尺寸比较大的粗网格 T H 对整体区域进行网格剖分, 得到有限元离散空间 S H ; 采用尺寸比较小的细网格 T h 对区域也进行网格剖分, 得到有限元离散空间 S h 根据虚位移原理, 式 (7) 对应的虚功方程为式中,( a,b ) = a T bd ì í î v,δv ) =( c c,δv ) c ( L u,δu ) = ( ),δu ;v Γ = u Γ (8) 将式 (8) 中两式相加并注意到 ( v ) = ( u c ) ; ( δv ) = ( δu ), 整理可得 ì í î v - c,δv ) = u c - c,δv ) - ( ) ( L u,δu ) = ( ) L u -,δv (9),δu ;v Γ = u Γ 18

1 组合网格法基本思想组合网格法是一种有叠加的区域分解算法, 其基本思想是采用粗网格和细网格两套网格进行分析计算, 粗网格对应整个区域, 不考虑线形结构的影响, 细网格则为包含线形结构在内的局部区域两套网格之间反复

4 令细网格上的基函数为 { Φ j } n, 粗网格上的基函数为 { Ψ j } n c, 它们之间存在关系 Ψ = DΦ, 其中 æ d 11 d 1L D = ç 为插值矩阵,L 为粗网格单元结点个数,n ç d n c 分别为细网格和粗网格结点总个数 n 1 d n L ø 推导得到边界 Γ 上满足 [19] 组合网格在无压稳定渗流中的应用张顺福丁留谦 n c n v = v j ψ j = V T Ψ;u = u j Φ j = U T Φ (10) U = D T V (11) 3.3 程序编制组合网格算法在三维稳定无压渗流计算中实现的具体流程如图 3 [25] 所示, 其中 a 场代表粗网格区域,b 场代表采用细网格控制参数的细网格区域,c 场代表采用粗网格控制参数的细网格区域,U U c 分别为 b 场和 c 场的水头值据此编制了相应的三维有限元程序 IIFM3DS-CGM [25] 初始化 ; 迭代步 n=0, Uc U V 赋初值 n=n+1 计算 a 场 V n+1 插值获得细网格边界条件 Γ = D T V n + 1 b 场对 a 场的修正 ( L Φ - ( L c 场对 a 场的修正 Φ - ( c 计算修正项 D é ù ( ë ê L c Φ ú - ( c - D é ù ( û ë ê L Φ ú - ( û 图 3 计算流程 4 算例某无限长土坝宽 m 高 30.0 m, 中部有 0.5 m 的垂直心墙, 心墙距离上游 50.0 m, 如图 4 所示心墙和坝壳均为各向同性材料, 上游水位为 30.0 m, 下游水位为 6.0 m 坝长方向取单位长度, 采用组合网格法, 用细网格模拟心墙及其上下游 17.0 m 区域范围, 整体区域则用粗网格进行离散, 其中细网格心墙处的六面体单元尺寸为 m 1 m 0.25 m, 心墙上下游 2 m 范围内单元尺寸为 0.25 m 1 m 0.25 m, 其它区域单元尺寸为 1 m 1 m 0.25 m, 如图 5 粗网格不考虑心墙的影响, 单元尺寸大小 m 1 m 1 m, 如图 6 为验证组合网格法的精度, 文中同时采用较细的网格对整体区域进行离散, 整体细网格在心墙及其上下游 17 m 的范围内采用与图 6 一样的网格剖分, 其它区域用尺寸为 1 m 1 m 0.25 m 的单元进行离散表 1 给出了坝壳渗透系数 k 1 与心墙渗透系数 k 2 不同比值下, 组合网格法计算得到的心墙下游出水 19

3 程序编制组合网格算法在三维稳定无压渗流计算中实现的具体流程如图 3 [25] 所示, 其中 a 场代表粗网格区域,b 场代表采用细网格控制参数的细网格区域,c 场代表采用粗网格控制参数的细网格区域,U U c 分别为 b 场和 c 场的水

5 Z X 图 4 计算模型 Z X 图 5 不考虑心墙的整体区域及相应粗网格剖分 ( 组合网格法 ) Z X 图 6 考虑心墙的局部区域及相应细网格剖分 ( 组合网格法 ) 表 1 心墙下游出水点高程比较 k1/k2 组合网格 /m 整体细网格 /m 相对误差 /% 点高程与整体区域细网格的计算结果比较, 从表 1 可以看到, 组合网格法计算的结果与整体区域细网格的计算结果基本一致, 最大误差为 2.3 % 图 7 给出了不同渗透系数比值条件下, 组合网格法计算得到的心墙附近自由面与整体区域细网格计算结果的比较从图 7 可以看到, 组合网格法的计算结果与整体区域细网格的结果非常接近 5 结语无压稳定渗流有限元分析中, 线形结构增加了计算区域的离散难度本文介绍了组合网格法的 20

4 点高程与整体区域细网格的计算结果比较, 从表 1 可以看到, 组合网格法计算的结果与整体区域细网格的计算结果基本一致, 最大误差为 2.

6 图 7 组合网格与整体细网格的自由面比较基本原理, 并将其引入到稳定无压渗流有限元分析中, 为线形结构的模拟提供了一种新思路由于组合网格的粗细网格可以独立生成, 因此组合网格法对线形结构的网格剖分具有较强的适应性, 能极大的简化无压渗流计算的网格剖分工作组合网格法在粗网格区域达到粗网格的精度, 在细网格区域达到细网格的精度结合改进的初流量法给出了稳定无压渗流组合网格算法格式和计算流程, 并编制了相应的三维有限元程序 IIFM3DS-CGM 采用垂直心墙算例对算法和程序进行了验证, 结果证明算法是合理的, 程序是有效的可信的本文的研究成果还需结合实际工程进行更深入的研究参考文献 [ 1 ] 王浩然, 朱国荣. 淄博市孝妇河源区地下水资源的开发利用研究 [J]. 高校地质学报,2000,4(2) [ 2 ] 薛禹群, 谢春红. 水文地质学的数值法 [M]. 北京煤炭工业出版社,1980. [ 3 ] Bernardi C,Maday Y,Patera A. Domain decomposition by the mortar inite element method[j]. Asymptotic and Numerical Methods or PDEs with Critical Parameters,1993, [ 4 ] Mc Cormick S,Thomas J. The ast adaptive composite grid(fac)method or elliptic equations[j]. Mathematics o Computation,1986, [ 5 ] 吕涛, 石济民, 林振宝. 区域分解算法 [M]. 北京科学出版社,1992. [ 6 ] Mc Cormick S. Fast adaptive composite grid(fac)methodstheory or the variational case[j]. Deect Correction Methods Computing Supplementum,1984,

IIFM3DS-CGM 采用垂直心墙算例对算法和程序进行了验证, 结果证明算法是合理的, 程序是有效的可信的本文的研究成果还需结合实际工程进行更深入的研究参考文献 [ 1 ] 王浩然, 朱国荣.

7 [ 7 ] 江思珉, 朱国荣, 王浩然, 等. FAC 方法在地下水数值模拟中的应用 [J]. 水利学报,2006,37(11) [ 8 ] 陶文铨. 数值传热学 [M]. 第 2 版. 西安西安交通大学出版社,2001. [ 9 ] 郭胜山, 李德玉, 唐菊珍. 水工结构动接触问题的组合网格算法 [J]. 水力发电,2009(5) [ 10 ] 韩修廷, 梁国平, 吴恩成, 等. 油田套管损坏数值模拟专用软件开发及应用 [J]. 数字石油和化工,2006 (4) [ 11 ] 王伯军, 张士诚, 张劲, 等. 考虑注采关系的三维地应力场数值模拟研究 [J]. 西南石油大学学报自然科学版,2006(6) [ 12 ] 甘艳, 阮江军, 张宇, 等. 组合网格法及其在电磁问题中的应用 [J]. 电工技术学报,2008,23(8)8-14. [ 13 ] 张宇, 阮江军, 刘兵, 等. 组合网格法在电磁 - 机械耦合问题中的应用 [J]. 中国电机工程学报,2007,27 (37) [ 14 ] 陈文平. 组合网格法及其在焊接数值模拟中的应用 [D]. 福州福建师范大学,2009. [ 15 ] 陈文平, 马昌凤, 蒋利华. 组合网格法在搅拌摩擦焊接数值模拟中的应用 [J]. 安徽理工大学学报自然科学版,2009(2) [ 16 ] 陈文平, 马昌凤, 蒋利华. 激光焊接问题的组合网格法 [J]. 福建师范大学学报自然科学版,2009(3) [ 17 ] 王洪涛. 多孔介质污染物迁移动力学 [M]. 北京高等教育出版社,2008. [ 18 ] 张顺福, 丁留谦, 刘昌军, 等. 基于子域积分的初流量法改进 [J]. 水电能源科学,2013,31(3) [ 19 ] 王德生. 组合网格法和非结构化网格自动生成 [D]. 北京中国科学院数学研究所,2001. [ 20 ] Xu J. Iterative Methods by Space Decomposition and Subspace Correction[J]. SIAM Review,1992,34(4) [ 21 ] Xu J. The auxiliary space method and optimal multigrid preconditioning techniques or unstructured grids[j]. Computing,1996,56(3) [ 22 ] 张有天, 陈平, 王镭. 有自由面渗流分析的初流量法 [J]. 水利学报,1988(8) [ 23 ] 李新强. 无压渗流有限元分析的改进初流量法 [J]. 水利学报,2007,38(8) [ 24 ] 王媛. 求解有自由面渗流问题的初流量法的改进 [J]. 水利学报,1998(3) [ 25 ] 张顺福. 子域积分法及组合网格法在无压渗流中的应用 [D]. 北京中国水利水电科学研究院,2013. Application o composite grid method in unconined seepage analysis ZHANG Shunu 1,2,DING Liuqian 1,2 (1. China Institute o Water Resources and Hydro Power Research,Beijing ,China; 2. Research Center on Flood & Drought Disaster Reduction o the Ministry o Water Resources,Beijing ,China) AbstractTo solve the line-style structure problem oten encountered in unconined steady seepage analy sis, Composite Grid Method is adopted. Two sets o independent grids are used in numerical simulation. One is the coarse grid o larger size simulating the entire region without consideration o line-style struc ture and the other is the ine grid o relatively small size simulating the domain including the line-style structure. Adjustments are carried out iteratively between the coarse grid and ine grid until desired conver gence precision is achieved. With Composite Grid Method, the amount o discretization is greatly reduced. The Composite Grid Method can be adapted to non-regular grid. Numerical simulation result shows that the theory is reasonable and the program is reliable. Composite Grid Method provides an approach to solving the line-style structure problem in unconined steady seepage analysis. Key words unconined steady seepage; Composite Grid Method; line-style structure; improved Initial Flow method;finite element method ( 责任编辑李琳 ) 22

考虑注采关系的三维地应力场数值模拟研究 [J]. 西南石油大学学报自然科学版,2006(6)33-35. [ 12 ] 甘艳, 阮江军, 张宇, 等. 组合网格法及其在电磁问题中的应用 [J]. 电工技术学报,2008,23(8)8-14.

何秋琳张立春视觉学习研究进展视觉注意视觉感知

何秋琳张立春视觉学习研究进展视觉注意视觉感知第卷第期年月开放教育研究何秋琳张立春华南师范大学未来教育研究中心广东广州随着图像化技术和电子媒体的发展视觉学习也逐步发展为学习科学的一个研究分支得到研究人员和教育工作者的广泛关注基于此作者试图对视觉学习