浙江大学学报(工学版)

浙江大学学报(工学版), 2023, 57(10): 2116-2125 doi: 10.3785/j.issn.1008-973X.2023.10.020

土木工程

跨尺度网格识别方法及在防渗墙分析中的应用

余翔,, 赖远平, 王钰轲,, 屈永倩, 郑浩然

1. 中国地球物理学会工程物探检测重点实验室，湖北 武汉 430000

2. 郑州大学 水利与土木工程学院，河南 郑州 450001

3. 大连理工大学 水利工程学院，辽宁 大连 116024

Identification method of cross-scale mesh and application in analysis of cutoff wall

YU Xiang,, LAI Yuan-ping, WANG Yu-ke,, QU Yong-qian, ZHENG Hao-ran

1. Key Laboratory of Engineering Geophysical Prospecting and Detection of Chinese Geophysical Society, Wuhan 430000, China

2. School of Water Conservancy and Civil Engineering, Zhengzhou University, Zhengzhou 450001, China

3. School of Hydraulic Engineering, Dalian University of Technology, Dalian 116024, China

通讯作者: 王钰轲，男，教授. orcid.org/0000-0003-1849-4857. E-mail： ykewang@163.com

收稿日期: 2022-11-30  

基金资助: 中国地球物理学会工程物探检测重点实验室开放研究基金资助项目（CJ2021D05）；国家自然科学基金资助项目（52192670, 52109151, 51809034）；中国博士后科学基金资助项目（2021M692938）；河南省省重点研发与推广专项资助项目（222102320098）；云南省重大科技专项计划资助项目（202102AF08002）

Received: 2022-11-30  

Fund supported: 中国地球物理学会工程物探检测重点实验室开放研究基金资助项目（CJ2021D05）；国家自然科学基金资助项目（52192670,52109151,51809034）；中国博士后科学基金资助项目（2021M692938）；河南省省重点研发与推广专项资助项目（222102320098）；云南省重大科技专项计划资助项目（202102AF08002）

作者简介 About authors

余翔（1988—），男，副教授，从事数值仿真方法开发与应用研究.orcid.org/0000-0002-7488-3550.E-mail：xiangyu@zzu.edu.cn , E-mail：xiangyu@zzu.edu.cn

摘要

为了高效地生成二维比例边界有限元-有限元（SBFEM-FEM）跨尺度耦合网格，提出跨尺度网格识别方法. 该方法结合CAD人为可控的特点，对图形线段进行识别、裁剪、整理，生成有效节点与线段；根据节点与线段的拓扑关系以及封闭域的构造，生成合理的多边形比例边界单元或有限单元；最后组装节点与单元信息，输出可以用于数值分析的二维SBFEM-FEM跨尺度耦合网格. 为了验证生成SBFEM-FEM跨尺度耦合网格的有效性及计算精度，依托某堤坝工程，对比分析了防渗墙不同网格剖分方法获得的墙体应力变形分布规律及峰值. 结果表明，采用常规大尺度有限元网格模拟获得的主应力误差可以超过48%，基于所提跨尺度精细数值网格获得的结果误差不能超过5%，且其网格单元量大幅度降低. 跨尺度网格识别方法可以为堤坝工程防渗结构提供有力的支持.

关键词： SBFEM-FEM耦合 ; 跨尺度网格生成 ; 精细模拟 ; 堤坝防渗墙 ; 可控性

Abstract

A cross-scale mesh identification method was presented to efficiently create two-dimensional scaled boundary finite element method-finite element method (SBFEM-FEM) cross-scale coupled meshes. The method combined the artificially controllable characteristics of CAD to identify, cut, and organize graphic line segments generating effective nodes and line segments. A reasonable polygonal-scaled boundary element or finite element was generated based on the topological relationship between the nodes and the line segments and the construction of the closed domain. The nodes and element information were assembled. The two-dimensional SBFEM-FEM cross-scale coupled mesh suitable for numerical analysis was produced. The stress and deformation distribution law and peak value of the wall obtained by different mesh generating methods were compared and studied to verify the validity and calculation accuracy of the generated SBFEM-FEM cross-scale coupling mesh based on a dam project. Results showed that the error of the principal stress obtained by the conventional large-scale finite element mesh simulation could exceed 48%, the error of the results obtained based on the proposed cross-scale fine numerical mesh could not exceed 5%, and the number of mesh elements was greatly reduced. The cross-scale mesh identification method can provide strong support for the anti-seepage structures in dam engineering.

Keywords： coupled scaled boundary finite element method-finite element method ; generation of cross-scale mesh ; refined simulation ; cutoff wall of embankment ; controllability

PDF (3206KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

余翔, 赖远平, 王钰轲, 屈永倩, 郑浩然. 跨尺度网格识别方法及在防渗墙分析中的应用. 浙江大学学报(工学版)[J], 2023, 57(10): 2116-2125 doi:10.3785/j.issn.1008-973X.2023.10.020

YU Xiang, LAI Yuan-ping, WANG Yu-ke, QU Yong-qian, ZHENG Hao-ran. Identification method of cross-scale mesh and application in analysis of cutoff wall. Journal of Zhejiang University(Engineering Science)[J], 2023, 57(10): 2116-2125 doi:10.3785/j.issn.1008-973X.2023.10.020

土石堤坝工程是我国重要的基础设施工程，为了保证设计的可靠性以及工程的安全性，常常会对工程结构的工作状态进行分析. 网格剖分是有限元数值模拟的基础，因此网格尺寸大小对数值计算结果的可靠性和精确性影响很大[1]. 防渗墙作为控制堤坝渗漏的关键结构，是堤坝安全的控制部位. 为了准确定位在不同运行条件下，防渗墙的薄弱部位会对防渗墙进行精细的网格剖分. 防渗墙与堤体尺度跨越巨大，若按照均匀尺度的原则对堤坝工程进行网格剖分，防渗墙的小尺度网格尺寸标准将会造成数值模型具有大量自由度，计算量难以承受. 相比于防渗墙，堤体的应力梯度一般较为平缓，常规尺度网格能够满足其精度要求. 数值模拟对防渗墙的核心结构进行精细网格剖分，其它部分采用常规尺度网格剖分，获得了跨尺度网格，并引入数值计算方法，数值分析的计算量和精度均能够得到可靠保证. 数值分析的计算量和精度同时得到保证是目前数值分析进一步推广应用的难题.

根据比例边界有限元法（scaled boundary finite element method, SBFEM）获得的网格单元对复杂几何边界适应能力强，具有实现跨尺度、精细化离散的特点. 邹德高等[2-10]对比例边界有限元法进行发展和改进，取得了一定的成果. 结合比例边界有限元法对构成单元节点没有限制的特点以及有限元法常规单元技术成熟易于求解的特点，SBFEM-FEM耦合求解分析逐渐发展起来. 邹德高等[11-13]将该耦合法应用于结构应力分析，验证了该耦合法的可行性. 殷德胜等[14-17]将耦合法应用于裂缝扩展和损伤破坏. Ye等[18]将该耦合法应用于弹性板结构与多层无界弹性地基相互作用问题的求解. 目前SBFEM-FEM跨尺度分析所采用的数值网格通常由四分树法获得. 四分树法获得的网格单元形状固定、可控性差，无法适应土石堤坝中分层填筑和材料分区复杂的边界条件，难以根据实际情况相应做出合理的网格剖分方式，且可能造成需要精细模拟的部位会被忽略[19].

为了增强网格剖分的可控性，在CAD软件中进行人为可控的网格划分，再识别并生成数值网格. 李华等[20-22]提出二维有限元网格的识别和生成方法. 田林等[23-27]在CAD环境下，对二维网格生成方法进行改善优化，提高了网格的质量和生成效率. 籍冉冉等[28]对莫尔斯（Morse）算法进行优化，提高了所生成的四边形网格质量. 目前均是针对三角形和四边形单元的有限元数值网格，而有关非常规的多边形单元或SBFEM-FEM耦合数值网格识别与生成的研究十分少见.

本研究结合CAD人为可控的特点，对图形线段进行预处理，生成有效节点与线段，根据所生成的节点和线段之间的拓扑关系及封闭域构造，构建可以直接用于数值分析的网格模型，以及二维SBFEM-FEM耦合网格生成方法. 本研究将跨尺度网格识别方法应用于某实际堤坝工程的数值分析研究中，验证了网格生成方法获得的SBFEM-FEM耦合网格的有效性和精度.

1. SBFEM-FEM耦合分析方法

比例边界有限元法是由Wolf等[29-30]提出，该方法在边界处利用Galerkin方法进行数值离散，降低了一个空间的计算维度，可以在边界内部不需要基本了解的情况下利用解析的方法求解. 比例边界有限元法不仅具有较高的精度，而且对单元形状的容许度更高，不需要根据单元的结点数量开发相应单元类型，突破了常规有限元法的限制. 图1为典型的多边形比例边界单元，图中ξ为单元径向坐标，s为环向坐标.

图 1

新窗口打开| 下载原图ZIP| 生成PPT

图 1   比例边界有限元法中的边界离散单元

Fig.1   Boundary discrete element in scaled boundary finite element method


对于任意一个扇形区，通过边界积分点和特征值分解技术可求得单元形函数和应变位移转换矩阵：

$ {\boldsymbol{B}}(\xi ,s) = [{{\boldsymbol{B}}_1}(s){{\boldsymbol{\varPsi }}_{\text{u}}}[ - {{\boldsymbol{S}}_n}]+{{\boldsymbol{B}}_{\text{2}}}(s){{\boldsymbol{\varPsi }}_{\text{u}}}]{\xi ^{ - {{\boldsymbol{S}}_n} - {\boldsymbol{I}}}}{\boldsymbol{\varPsi }}_{\text{u}}^{{{ - 1}}} . $

式中： ${\boldsymbol{\varPsi}}_{\rm{u}}$为位移模式对应的转换矩阵，可以通过特征值分解求得； ${\boldsymbol{S}}_n $为特征值分解获得的特征值矩阵；n为单元边数； ${\boldsymbol{B}}(\cdot) $为应变-位移的转换矩阵；I为单位矩阵. 通过采用有限元的原理以及域内积分点，即可求出单元刚度矩阵：

$ {{\boldsymbol{K}}_{{\text{ep}}}} = \sum\limits_{i = 1}^{3n} {{{\boldsymbol{B}}^i}(\xi ,s){\boldsymbol{D}}_{{\text{ep}}}^i{{\boldsymbol{B}}^i}(\xi ,s){F_i}} \text{，} $

$ {{\boldsymbol{R}}_{{\text{int}}}} = \sum\limits_{i = 1}^{3n} {{{\boldsymbol{B}}^i}{{(\xi ,s)}^{\rm{T}}}{\sigma _i}(\xi ,s){F_i}} \text{，} $

$ {\sigma _i}(\xi ,s) = \sum\limits_{i = 1}^{3n} {{\boldsymbol{D}}_{{\text{ep}}}^i} {\varepsilon _i}(\xi ,s) . $

式中： ${\boldsymbol{K}}_{\rm{ep}} $为弹塑性刚度矩阵， ${\boldsymbol{D}}_{\rm{ep}} $为弹塑性本构矩阵， $F $为每个高斯点所代表的面积， ${\boldsymbol{R}}_{{\rm{int}}} $为内力向量，σ为单元应力场，ε为应变。

通过网格跨尺度剖分，实现对核心结构网格加密，通常需要用多级过渡，如图2（a）所示. 单元尺寸过渡中必然会引起5节点及以上节点数目的单元，而这类单元需要单独构造形函数、积分方法，计算较为麻烦. 不需要根据单元的结点数量开发相应单元类型是比例边界有限元法的一个特点，网格跨尺度通过一级比例边界有限单元就可以过渡完成，如图2（b）所示. 过渡单元两侧采用有限单元模拟，且采用本研究方法进行网格任意加密，有效地解决了复杂几何区域的网格剖分.

图 2

新窗口打开| 下载原图ZIP| 生成PPT

图 2   跨尺度网格过渡方法的示意图

Fig.2   Schematic diagram of cross-scale mesh transition method


2. SBFEM-FEM耦合网格生成方法

2.1. CAD数据处理

本研究提出的SBFEM-FEM耦合网格生成方法主要分3步：处理CAD数据；对节点和线段信息预处理；生成数值网格信息. 以下结合图3所示的结构，说明耦合网格的生成方法.

图 3

新窗口打开| 下载原图ZIP| 生成PPT

图 3   悬臂梁结构的示意图

Fig.3   Schematic diagram of cantilever beam structure


在CAD软件中，根据模拟对象创建节点与线段. 在CAD中绘出结构的二维简化图形，如图4所示. 提取图中6条线段、12个节点的信息，作进一步处理. 首先求解线段的交点，并将线段在交点处裁剪为多个子线段；判断节点的公用次数，若公用次数不大于1，则该节点无效，且该节点所属线段也无效；最后对所有有效节点和线段重新进行排序生成基本网格，并存储节点和线段信息如图5所示.

图 4

新窗口打开| 下载原图ZIP| 生成PPT

图 4   悬臂梁结构的简化示意图

Fig.4   Simplified schematic diagram of cantilever beam structure


图 5

新窗口打开| 下载原图ZIP| 生成PPT

图 5   悬臂梁结构的基本网格图

Fig.5   Basic mesh diagram of cantilever beam structure


2.2. 节点与线段信息预处理

1）创建节点对称矩阵A，如表1所示. 通过节点关系矩阵A快速确定节点之间的连接关系，相应的节点连接关系，如图6所示. 若是 $ {{A}}(i,j) = 1 $，则i节点和j节点可组成ij线段. 矩阵A的上三角矩阵中值为1的个数是整个网格中的线段数量，应与上一步所识别出的有效线段数量相等.

表 1   节点连接关系矩阵A

Tab.1  Node connection relation matrix A

节点编号 1 2 3 4 5 6 7 8
1 — 1 — — — 1 — —
2 1 — 1 1 — — — —
3 — 1 — — 1 — — —
4 — 1 — — 1 — 1 —
5 — — 1 1 — — — 1
6 1 — — — — — 1 —
7 — — — 1 — 1 — 1
8 — — — — 1 — 1 —

新窗口打开| 下载CSV


图 6

新窗口打开| 下载原图ZIP| 生成PPT

图 6   节点连接的关系图谱

Fig.6   Relationship diagram of node connection


2）创建单个节点连接向量VNL，如表2所示. $ {\rm{VNL}}(i) $用于存储矩阵A的上三角矩阵第i行中值为1的个数.

表 2   节点连接的数目

Tab.2  Number of node connections vector

节点编号 1 2 3 4 5 6 7 8
连接数目 2 2 1 2 1 1 1 0

新窗口打开| 下载CSV


3）处理共用某一节点的所有线段. 创建3个向量NLN、NLC与NLA，分别存储共用节点线段的编号、在相应共用线段中的节点排序（值为1或2）以及以节点为基准的线段与向量(1,0)的夹角，并根据NLA由小到大将这些线段进行逆时针排序.

2.3. 数值网格信息生成

2.3.1. 遍历有效节点

1）根据节点编号顺序，依次判别. 当节点编号为 $ a $时，判别 $ {\rm{VNL}}(a) $大小（VNL(a)为矩阵A的上三角矩阵第a行中值为1的个数），若 $ {\rm{VNL}}(a) \geqslant 2 $，则该节点可生成新单元，否则不可生成，且进行下一节点的判别. 以图7为例，判别出节点 $ a $可以生成新单元，则进入第2.3.2节遍历共用当前节点的线段. 2）当所有节点遍历结束，则网格识别完成.

图 7

新窗口打开| 下载原图ZIP| 生成PPT

图 7   单元生成的判别方法及流程

Fig.7   Identification method and process of unit generation


2.3.2. 遍历共用当前节点的线段

1）识别共用节点的线段，并确定过渡点. 每相邻的2条线段( ${\rm{ NLN}}(i+ 1)$与 ${\rm{ NLN}}(i) $)为一组，其中必有 $ {\rm{ NLA}}(i+1) > {\rm{ NLA}}(i) $. 确定过渡点后进入第2.3.3步. 2）所有相邻共用线段遍历结束，则回至2.3.1步中第1）步重新判别下一个节点的VNL值大小.

以共用节点 $ a $的线段为例，若单元有n个节点，则节点 $ {{a}} $为单元的起始点，即第1个节点. 线段 ${\rm{ NLN}}(i+1) $为单元的终结边，该线段另一节点b为单元的终结点，即第n个节点；线段 $ {\rm{ NLN}}(i) $为单元的起始边，该线段另一节点m为单元的过渡点.

2.3.3. 判别过渡点与终结点的连接关系

1）若过渡点与终结点可连成线段，即 $ {{A}}(i+j) = 1 $，则新单元生成，回至第2.3.2节判断共用节点的下一组线段. 2）若 $ {{A}}(i+j) \ne 1 $，则该过渡点为单元中间点，继续寻找下一个过渡点. 3）顺序循环执行本步第1）、2）步.

以上述共用节点a为例，若是 $ {{A}}(b,m) = 1 $，则组成单元a-m-b；若是 $ {{A}}(b,m) \ne 1 $，则m组成新单元的中间点，记录该节点并继续查找组成单元的第n-1个节点. 查询共用节点 $ {{m}} $的所有线段，确定线段 $ {\rm{ NLN}}(i) $在共用节点 $ m $的所有线段中的排序j，基于组成单元的所有节点须逆时针排序的原则，定位线段 $ {\rm{NLN}}(j - 1) $，获得与节点 $ m $共同组成线段 $ {\rm{NLN}} (j - 1) $的另一节点，该节点为此时的过渡节点m，然后继续判别 $ {{A}}({{b}}, {{m}}) $是否等于1. 该例最终生成单元a- m1-m2- m3-b.

用上述方法将图3所示结构划分为用于数值分析的网格，如图8所示. 单元尺寸为在CAD软件中所画线段的实际长度，通过所提方法的识别节点和线段直接确定单元尺寸，材料区域的数量对网格的识别无任何影响. 在网格识别成功后，根据材料分区线、填筑分层线以及界面位置等控制信息自动完成材料区域的划分，接触单元生成计算分析前的关键步骤.

图 8

新窗口打开| 下载原图ZIP| 生成PPT

图 8   最终识别出的悬臂梁结构的跨尺度网格

Fig.8   Cross-scale mesh of cantilever beam structure of final identification


3. 应用案例及分析

某堤坝位于河南省信阳市罗山县，是以防洪、灌溉为主的小（1）型水库. 该水库工程等级为第Ⅳ级，主要建筑物为4级. 该水库主坝为均质土质堤坝，坝基由重粉质壤土、细砂及云母石英片岩组成. 坝基土体渗透系数为 $ 1.6 \times {10^{ - 6}} \sim 8.9 \times {10^{ - 5}}{\text{ cm/s}} $，平均渗透系数为 $ 2.6 \times {10^{ - 5}}{\text{ cm/s}} $，具有微弱的透水性，在堤坝轴线剖位采用混凝土构筑悬挂式防渗墙进行防渗.

3.1. 计算模型

图9为工程的二维模型，X为顺河向，Y为高程方向. 坝高15.00 m，坝顶宽5.60 m，坝底宽89.66 m. 上游坡比1∶2.73和1∶2.58，混凝土防渗墙位于坝体的中部，高17 m，嵌入坝基 $ {\text{2 m}} $，厚度50 cm. 坝基高20 m，基岩高2 m，坝基、基岩总水平长度290 m. 坝体分层填筑和蓄水模拟共分为29个荷载步：

图 9

新窗口打开| 下载原图ZIP| 生成PPT

图 9   堤坝的材料分布及荷载步的示意图

Fig.9   Sketch of material distribution and load steps of embankment


1）分16步（地基作为初始应力）对坝体从坝基面（高程22 m）填筑至坝顶（高程37 m），墙体材料按土体进行计算；2）第17步将墙体位置的土体换成混凝土材料进行计算；3）再分12步从坝基面蓄水至坝顶. 其中，蓄水引起的静水压力以面力的形式施加于防渗墙的上游面与下游面（嵌入坝基部分），如图10所示.

图 10

新窗口打开| 下载原图ZIP| 生成PPT

图 10   堤坝防渗系统受力示意图

Fig.10   Schematic diagram of force on anti-seepage system


3.2. 材料参数

堤坝体与堤坝基的静力分析了采用邓肯-张E-B非线性模型[31]，坝体、坝基材料参数取自实际试验结果[32]. 参数如表3所示. 其中 $ {\gamma _{\text{d}}} $、 $ {\gamma _{\text{f}}} $分别为干重度和浮重度； $ c $、 $ \varphi $和 $ \Delta \varphi $分别为凝聚力、摩擦角和摩擦角随围压增大10倍的减小量，三者均为强度参数； $ k $和 $ \;\beta $分别为弹性模量系数和指数； $ {R_{\text{f}}} $为破坏比； $ {k_{\text{b}}} $和q分别为体积模量系数和指数参数； $ {k_{{\text{ur}}}} $和 $ {n_{{\text{ur}}}} $分别为卸载模量系数和指数.

表 3   坝体和坝基的静力参数

Tab.3  Static parameters of dam body and foundation

材料 $ {\gamma _{\text{d}}} $/(kg·m−3) $ {\gamma _{\text{f}}} $/(kg·m−3) $ k $ $ {k_{{\text{ur}}}} $ $ \;\beta $ $ {n_{{\text{ur}}}} $
坝体 16.3 982 300 360 0.34 0.34
坝基 16.5 982 320 390 0.30 0.30
材料 $ {R_{\text{f}}} $ $ {k_{\text{b}}} $ $ q $ $ c $/kPa $ \varphi $/(°) $ \Delta \varphi $/(°)
坝体 0.95 200 0.30 22.2 11.3 0
坝基 0.95 215 0.30 21.6 11.8 0

新窗口打开| 下载CSV


混凝土防渗墙和基岩材料均采用线弹性本构模型计算，材料参数如表4所示，其中 $ \;\rho $为密度、E为弹性模量、 $ \nu $为泊松比.

表 4   墙、基岩的静力参数

Tab.4  Static parameters of wall and bedrock

材料 $\; \rho$/(kg·m−3) $ E/{\text{MPa}} $ $ \nu $
混凝土防渗墙 2 400 30 000 0.167
基岩 2 400 200 0.350

新窗口打开| 下载CSV


采用Goodman接触单元模拟墙与土之间的接触面，用Clough-Duncan双曲线模型[33]表达其应力-应变关系，该模型认为单元切向应力与切向相对位移之间的关系呈双曲线关系. 则切向刚度为

$ \tau = \frac{\omega }{{u+v\omega }} , $

$ {k_{\text{s}}} = K{\gamma _{\text{w}}}{\left(\dfrac{{{\sigma _{\text{n}}}}}{P}\right)^j}{\left(1 - {R_{\text{f}}}\frac{\tau }{{{\sigma _{\text{n}}}\tan \delta }}\right)^2} . $

式中：τ为切向刚度，ω为应变力，u为实验参数，K 、Rf和j为非线性参数，δ为接触面的界面摩擦角，γw为水的重度，P为大气压，σn为法向应力.

董景刚[34]对接触面单元的法向刚度应取值为相邻材料弹性模量的50~100倍，接触面模型的法向刚度与切向刚度的关系大致为 $ {k_{\text{n}}} = 2.5\;{k_{\text{s}}} $. 法向刚度 $ {k_{\text{n}}} $= $ {\text{1}}{\text{.5}} \times {\text{1}}{{\text{0}}^{12}}{\text{ kN/}}{{\text{m}}} $，是50倍的混凝土弹性模量. 初始切向刚度= $ 6.0 \times {10^{11}}{\text{ kN/}}{{\text{m}}} $，是法向刚度的0.4倍. 接触面材料参数取自实际试验结果[35]参数如表5所示.

表 5   接触面的静力参数

Tab.5  Static parameters of contact surface

位置 $ K $ $ j $ $ \delta /(^\circ ) $ $ {R_{\text{f}}} $ $ c/{\text{kPa}} $
墙与坝体 757 0.8 11.0 0.89 10.5
墙与坝基 757 0.8 11.0 0.89 10.5

新窗口打开| 下载CSV


3.3. 计算网格模型及验证

为了验证本研究数值网格生成方法的有效性以及SBFEM-FEM耦合数值网格的精度，分别根据以下3种情况剖分堤坝模型，并采用本研究建立的方法进行数值网格的生成. 堤坝主体耦合网格剖分图和网格剖分细图如图11所示. 模型1：大尺度常规单元网格（墙体单元竖向长度与土体单元保持一致）. 模型2：耦合的SBFEM-FEM单元网格（只在墙体进行局部加密，土体单元与模型1一致）. 模型3：精细尺寸常规单元网格（以模型2中墙体单元尺寸为基准，均匀剖分堤坝）.

图 11

新窗口打开| 下载原图ZIP| 生成PPT

图 11   墙体与土体局部单元示意图

Fig.11   Diagram of local elements of wall and soil


图12为模型2耦合SBFEM-FEM网格通过本研究方法识别生成的有限元模型局部网格示意图. 其中，常规单元单元类型号为1、多边形单元单元类型号为44、接触单元单元类型号为4表示. 可以看出，图12与图11（b）完全一致，表明了本研究建立的数值网格生成方法的有效性.

图 12

新窗口打开| 下载原图ZIP| 生成PPT

图 12   耦合SBFEM-FEM网格识别后的单元分布

Fig.12   Coupling SBFEM-FEM mesh after mesh recognition


表6和图13为3种网格下堤坝节点、单元信息和墙体单元信息，其中N*为节点总数量，N为模型总单元数量，n*为防渗墙单元数量，M为模型编号，H1和H2分别为墙体竖向和横向长度. 模型2耦合SBFEM-FEM网格单元数量与模型1相差很小，而墙体单元数量是模型1的3.8倍. 模型2墙体单元数量与模型3墙体单元数量相等，而总单元数量是模型3的约1/3.

表 6   节点、单元的信息

Tab.6  Information of node and unit

M N* N n* H1/m H2/m
1 14 906 14 560 144 0.500 0.125
2 15 381 14 968 552 0.125 0.125
3 48 075 47 634 552 0.125 0.125

新窗口打开| 下载CSV


图 13

新窗口打开| 下载原图ZIP| 生成PPT

图 13   3种模型下堤坝单元和墙体的单元数量

Fig.13   Element number of dam and wall for three conditions


本软件生成的计算模型采用大连理工大学抗震研究所自主研发的大型岩土工程静、动力非线性分析软件GEODYNA进行数值计算. 图14为堤坝蓄水完成后，3种模型下堤体的水平、竖向位移. 可以看出网格对堤坝的整体变形结果影响很小，3种模型下的变形等值线基本吻合. 当研究堤坝整体变形分布规律时，常规有限元模拟及分析方法也能获得较高的精度.

图 14

新窗口打开| 下载原图ZIP| 生成PPT

图 14   蓄水后堤体的位移分布图

Fig.14   Distribution of dam displacement after impoundment


图15为蓄水后墙体的水平位移. 图中 $ \Delta S $为水平位移， $ H $为墙体竖向高程. 从图15可以看出，墙体位移分布规律没有变化. 墙体最大位移出现在墙顶，墙体最大位移为7.3 cm，这是蓄水后上游水压力及上游堤体浮托力共同作用的结果. 防渗墙位于堤坝中部，其变形完全受堤体控制，墙体变形对单元尺寸不敏感. 当研究堤坝防渗墙变形规律时，采用单元尺寸较大的大尺度网格模拟.

图 15

新窗口打开| 下载原图ZIP| 生成PPT

图 15   蓄水后墙体的水平位移示意图

Fig.15   Horizontal displacement of wall after impoundment


图16~18为蓄水后墙体上、下游侧和蓄水后堤坝体（不含防渗墙）最大、最小主应力分布曲线. 图中 $ {\sigma _1} $为最大主应力， $ {\sigma _3} $为最小主应力. 从图16~18可以看出，3种网格获得的应力分布规律相同，但大尺度常规网格的计算结果与其它2种差别较大，本研究耦合SBFEM-FEM网格的计算结果与精细常规网格计算结果相差甚微. 大尺度常规网格计算结果低估了墙体峰值，其计算结果不能作为评价墙体安全性能的依据.

图 16

新窗口打开| 下载原图ZIP| 生成PPT

图 16   蓄水后墙体上游面的最大、最小的主应力分布图

Fig.16   Maximum and minimum principal stress distribution of wall in upstream after impoundment


图 17

新窗口打开| 下载原图ZIP| 生成PPT

图 17   蓄水后墙体下游面的最大、最小主应力分布图

Fig.17   Maximum and minimum principal stress distribution of wall in downstream after impoundment


图 18

新窗口打开| 下载原图ZIP| 生成PPT

图 18   蓄水后堤坝的最大、最小主应力分布图

Fig.18   Maximum and minimum principal stress distribution of embankment in downstream after impoundment


表7为蓄水后防渗墙主应力的峰值. 其中F为峰值，以压应力为正，拉应力为负；D为误差. 与精细网格相比，大尺度网格获得的墙体大主应力误差可以达到29.5%，小主应力误差接近50%，而耦合网格的误差均小于5%. 耦合SBFEM-FEM模拟与分析方法可用少量网格获得较高精度的结果.

表 7   3种模型下墙体的主应力

Tab.7  Principal stress of wall for three conditions

模型 上游 $ {\sigma _1} $ 上游 $ {\sigma _3} $ 下游 $ {\sigma _1} $ 下游 $ {\sigma _3} $
F/MPa D/% F/MPa D/% F/MPa D/% F/MPa D/%
1 3.76 19.7 −0.16 48.80 0.26 29.50 −3.37 20.88
2 4.63 1.11 −0.29 4.43 0.36 3.74 −4.20 1.43
3 4.69 — −3.04 — 0.38 — −4.26 —

新窗口打开| 下载CSV


4. 结　论

（1）对于堤坝工程，可以采用常规大尺度有限元网格分析堤体整体变形规律以及堤坝中部防渗墙的变形规律. 当研究防渗墙应力特性时，防渗墙须采用精细网格剖分. 与常规精细化有限元网格的结果相比，常规大尺度有限元网格模拟获得主应力的误差可以超过48%.

（2）基于SBFEM-FEM耦合数值网格的数值模拟，计算效率及精度均较高. SBFEM-FEM耦合数值模拟分析融合了SBFEM的灵活通用性以及FEM求解高效的特点，简便地完成了大尺度堤体有限元网格与精细墙体有限元网格的过渡. 相应跨尺度数值网格获得的墙体应力与常规精细化有限元网格的结果相比，误差不超过5%，而跨尺度数值网格单元量大幅度降低.

参考文献

崔焕平, 崔燕平, 王宗敏

混凝土非线性有限元分析中的网格尺寸效应

[J]. 混凝土, 2007, (6): 27- 29

[本文引用: 1]

CUI Huan-ping, CUI Yan-ping, WANG Zong-min

Mesh size effect in nonlinear finite element analysis of concrete

[J]. Concrete, 2007, (6): 27- 29

[本文引用: 1]

邹德高, 刘锁, 陈楷, 等

基于四叉树网格和多边形比例边界有限元方法的岩土工程非线性静动力分析

[J]. 岩土力学, 2017, 38 (增2): 33- 40

[本文引用: 1]

ZOU De-gao, LIU Suo, CHEN Kai, et al

Nonlinear static and dynamic analysis for geotechnical engineer-ing based on quadtree mesh and polygon scaled bou-ndary finite element method

[J]. Rock and Soil Mechanics, 2017, 38 (增2): 33- 40

[本文引用: 1]

ZOU D, CHEN K, KONG X, et al

An enhanced octree polyhedral scaled boundary finite element meth-od and its applications in structure analysis

[J]. Engineering Analysis with Boundary Elements, 2017, 84: 87- 107

DOI:10.1016/j.enganabound.2017.07.007     

ZOU D, TENG X, CHEN K, et al

An extended polygon scaled boundary finite element method for the nonlinear dynamic analysis of saturated soil

[J]. Engineering Analysis with Boundary Elements, 2018, 91: 150- 161

DOI:10.1016/j.enganabound.2018.03.019     

邹德高, 陈楷, 刘锁, 等

非线性比例边界有限元在面板坝分析中的应用

[J]. 土木与环境工程学报(中英文), 2019, 41 (3): 11- 18

ZOU De-gao, CHEN Kai, LIU Suo, et al

Application of nonlinear scaled boundary polygon element method in analysis of concrete face rockfill dam

[J]. Journal of Civil and Environmental Engineering, 2019, 41 (3): 11- 18

ZOU D, SUI Y, CHEN K

Plastic damage analysis of pile foundation of nuclear power plants under beyond-design basis earthquake excitation

[J]. Soil Dynamics and Earthquake Engineering, 2020, 136: 106179

DOI:10.1016/j.soildyn.2020.106179     

GENES M C, KOCAK S

Dynamic soilstructure interaction analysis of layered unbounded media via a coupled finite element/boundary element/scaled boun-dary finite element model

[J]. International Journal for Numerical Methods in Engineering, 2005, 62 (6): 798- 823

DOI:10.1002/nme.1212     

BEHNKE R, KALISKE M

Thermomechanical mod-eling of crack propagation in dynamically loaded ela-stomer specimens using a scaled boundary finite ele-ment approach

[J]. Pamm, 2015, 15 (1): 121- 122

DOI:10.1002/pamm.201510051     

李上明

基于比例边界有限元法动态刚度矩阵的坝库耦合分析方法

[J]. 工程力学, 2013, 30 (2): 313- 317

LI Shang-ming

Transient analysis method for dam-r-eservoir interaction based on dynamic stiffness of S-BFEM

[J]. Engineering Mechanics, 2013, 30 (2): 313- 317

YU X, KONG X, ZOU D, et al

Linear elastic and plastic-damage analyses of a concrete cut-off wall co-nstructed in deep overburden

[J]. Computers and Geotechnics, 2015, 69: 462- 473

[本文引用: 1]

邹德高, 陈楷, 余翔, 等

基于跨尺度精细方法的面板坝面板损伤演化尺寸效应分析

[J]. 土木与环境工程学报(中英文), 2019, 41 (6): 36- 42

[本文引用: 1]

ZOU De-gao, CHEN Kai, YU Xiang, et al

Size effect analysis of face slab damage evolution for high concrete face dam under earthquakes based on cross- scale fine method

[J]. Journal of Civil and Environmental Engineering, 2019, 41 (6): 36- 42

[本文引用: 1]

邹德高, 陈楷, 张仁怡, 等

基于SBFEM的心墙坝基座跨尺度精细应力分析

[J]. 人民长江, 2019, 50 (9): 168- 174

ZOU De-gao, CHEN Kai, ZHANG Ren-yi, et al

Cross-scale refined stress analysis on base-support of core wall concrete dam based on scaled boundary finite element method

[J]. Yangtze River, 2019, 50 (9): 168- 174

ZOU D, CHEN K, KONG X, et al

An approach integrating BIM, octree and FEM-SBFEM for highly efficient modeling and seismic damage analysis of building structures

[J]. Engineering Analysis with Boundary Elements, 2019, 104: 332- 346

DOI:10.1016/j.enganabound.2019.03.038      [本文引用: 1]

殷德胜, 尹栓, 周宜红

裂缝分析的比例边界有限元与有限元耦合的虚拟结构面模型

[J]. 计算力学学报, 2014, 31 (6): 735- 741

[本文引用: 1]

YIN De-sheng, YIN Shuan, ZHOU Yi-hong

Coupled SBFEM and FEM for crack analysis based on virtual discontinuous surface method

[J]. Chinese Journal of Computational Mechanics, 2014, 31 (6): 735- 741

[本文引用: 1]

YANG Z J, DEEKS A J

Fully-automatic modelling of cohesive crack growth using a finite element–scaled boundary finite element coupled method

[J]. Engineering Fracture Mechanics, 2007, 74 (16): 2547- 2573

DOI:10.1016/j.engfracmech.2006.12.001     

CHEN K, ZOU D, KONG X, et al

Elasto-plastic fine-scale damage failure analysis of metro structures based on coupled SBFEM-FEM

[J]. Computers and Geotechnics, 2019, 108: 280- 294

DOI:10.1016/j.compgeo.2018.12.030     

QU Y, ZOU D, KONG X, et al

Seismic cracking evolution for antiseepage face slabs in concrete faced rockfill dams based on cohesive zone model in explicit SBFEM-FEM frame

[J]. Soil Dynamics and Earthquake Engineering, 2020, 133: 106106

DOI:10.1016/j.soildyn.2020.106106      [本文引用: 1]

YE W, LIU J, FANG H, et al

High-performance an-alysis of the interaction between plate and multi-laye-red elastic foundation using SBFEM-FEM

[J]. Composite Structures, 2019, 214: 1- 11

[本文引用: 1]

陈楷. 基于比例边界有限元的岩土工程精细化分析方法及应用[D]. 大连: 大连理工大学, 2019: 27-38.

[本文引用: 1]

CHEN Kai. Research on technique and application in refined analysis of complicated geotechnical engineering structures based on scaled boundary finite element method [D]. Dalian: Dalian University of technology, 2019: 27-38.

[本文引用: 1]

李华, 刘晓俭, 钟万勰

二维有限元网格图象的识别

[J]. 大连理工大学学报, 1993, (4): 373- 379

[本文引用: 1]

LI Hua, LIU Xiao-jian, ZHONG Wan-xie

Recognition for two-dimension finite element mesh drawings

[J]. Journal of Dalian University of Technology, 1993, (4): 373- 379

[本文引用: 1]

赵大洲, 王延红

一种人机结合式二维有限元网格生成法

[J]. 华北水利水电学院学报, 1999, 2: 22- 24

ZHAO Da-zhou, WANG Yan-hong

An interactive method to generate 2D finite element mesh

[J]. Journalof North China Institute of Water Conservancy and Hydroelectric Power, 1999, 2: 22- 24

徐思浩, 杨建新

判别二维有限元网格图的简便算法

[J]. 计算力学学报, 2001, 18 (1): 61- 63

DOI:10.3969/j.issn.1007-4708.2001.01.011      [本文引用: 1]

XU Si-hao, YANG Jian-xin

A simple algorithm usedto judge the 2D finite element mesh drawing

[J]. Chinese Journal of Computational Mechanics, 2001, 18 (1): 61- 63

DOI:10.3969/j.issn.1007-4708.2001.01.011      [本文引用: 1]

田林, 孙志林

AutoCAD下非结构网格自动生成

[J]. 水动力学研究与进展(A辑), 2002, 1: 124- 128

[本文引用: 1]

TIAN Lin, SUN Zhi-lin

Automatic generation of unstructured grid under AutoCAD

[J]. Chinese Journal of Hydrodynamics, 2002, 1: 124- 128

[本文引用: 1]

杨星, 蔡开玺. 二维无结构三角形网格自适应生成技术[J]. 人民长江, 2008, 39(19): 101-104.

YANG Xing, CAI Kai-xi, Adaptive generation technology of two-dimensional unstructured triangular mesh[J]. Yangtze River, 2008, 39(19): 101-104.

冯占荣, 于建群, 王玉杰, 等

一种二维非结构网格改进算法的研究

[J]. 燕山大学学报, 2010, 34 (5): 405- 410

FENG Zhan-rong, YU Jian-qun, WANG Yu-jie, et al

An improved algorithm for generating two-dimension-al grid

[J]. Journal of Yanshan University, 2010, 34 (5): 405- 410

郭新强. 边界面法四边形网格生成研究与应用[D]. 长沙: 湖南大学, 2011: 10-27.

Guo Xin-qiang. Study on quadrilateral mesh generation for the boundary face method and their impleme-ntation [D]. Changsha: Hunan University, 2011: 10-27.

徐青, 奚鹏飞

四边形网格自动生成方法改进及工程应用

[J]. 中国农村水利水电, 2018, 9: 182- 186

[本文引用: 1]

XU Qing, XI Peng-fei

The method improvement of quadrilateral mesh automatic generation and engineering application

[J]. China Rural Water ang Hydropower, 2018, 9: 182- 186

[本文引用: 1]

籍冉冉, 郑晓朋, 雷娜, 等

用于保特征四边形网格生成的改进Morse算法

[J]. 大连理工大学学报, 2020, 60 (6): 647- 653

[本文引用: 1]

JI Ran-ran, ZHENG Xiao-peng, LEI Na, et al

Improved Morse algorithm for feature-preserving quadrilateral mesh generation

[J]. Journal of Dalian University of Technology, 2020, 60 (6): 647- 653

[本文引用: 1]

SONG C, WOLF J P

The scaled boundary fniteelement methodalias consistent infinitesimal fnite-element cellmethod for elasto-dynamics

[J]. Computer Methods in Applied Mechanics and Engineering, 1997, 147 (3-4): 329- 355

DOI:10.1016/S0045-7825(97)00021-2      [本文引用: 1]

WOLF J P, SONG C

Dynamic-stiffness matrix in time domain of unbounded medium by infinitesimal finite element cell method

[J]. Earthquake Engineering and Structural Dynamics, 1994, 11 (23): 1073- 1086

[本文引用: 1]

DUNCAN J M, CHANG C

Nonlinear analysis of stress and strain in soils

[J]. Journal of the Soil Mechanics and Foundations Division, 1970, 96: 1629- 1653

[本文引用: 1]

LI J, ZHANG J, WANG Y, et al

Seismic response of earth dam with innovative polymer antiseepage wall

[J]. International Journal of Geomechanics, 2020, 20 (7): 4020079

DOI:10.1061/(ASCE)GM.1943-5622.0001664      [本文引用: 1]

WAYNE C G, M. D J. Finite element analyses of retaining wall behavior [J]. Journal of the Soil Mechanics and Foundations Division, 1971, 97(12).

[本文引用: 1]

董景刚. 土与结构接触面力学特性研究[D]. 大连: 大连理工大学, 2011: 55-59.

[本文引用: 1]

DONG Jing-gang. Study on mechanical properties of interface between soil and structure[D]. Dalian: Dalian university of technology, 2011: 55-59.

[本文引用: 1]

傅华, 章为民. 坝基混凝土防渗墙与泥皮接触面试验研究 [C]// 第一届中国水利水电岩土力学与工程学术讨论会. 昆明: 中国水利学会, 2006: 498-500.

[本文引用: 1]

FU Hua, ZHANG Wei-min. Experimental study on c-ontact surface between concrete cutoff wall and mud cake in dam foundation [C]// 1th National Symposium on Geotechnical Engineering. Kunming: CHES, 2006: 458-500.

[本文引用: 1]

/

〈 〉