一种天然河道三维网络生成方法技术

本发明专利技术公开了一种天然河道三维网络生成方法，其特征在于，将天然河道整体模型分解成在水平方向上连续且互不重叠的柱状的子域，然后将各子域规则化后采用网格化软件对各子域进行网格剖分，再合并剖分后的各子域并构成天然河道三维网格组合体。本发明专利技术能够更好地实现天然河道的立体空间剖分，使其能够更好地直接完成天然河道三维水流数值模拟，提高天然河道三维水流数值模拟的准确性、可靠性和稳定性。可靠性和稳定性。可靠性和稳定性。

【技术实现步骤摘要】
一种天然河道三维网络生成方法


[0001]本专利技术涉及天然河道计算流体力学前处理
，具体涉及一种对宽浅不规则的天然河道三维网络生成方法。

技术介绍

[0002]航道整治、河道治理、水电开发等工程建设中常常需要开展天然河道三维水流数值模拟，河道三维网格剖分是数值模拟中最基础和最重要的工作，然而由于天然河道床面凹凸不平、宽浅不一、岸线弯曲，河道整体实体形状极不规则，使用软件自动划分三维网格十分困难，甚至不可能，因此极大限制了天然河道三维水流数模拟的应用。
[0003]CN201110372235.X公开了一种用于河道水文数值模拟的复杂河道网格化方法，首先将带有支流出入汇、江心洲等各种复杂河道区域自由分割为简单区域集合，根据河道数值模拟需要设置区域控制曲线网格间距，生成二维网格，归并简单区域网格得到复杂河道剖分后的初始网格。这样对于河道岸线特定区域，可根据需求灵活增加区域网点密度实现该区域内网格加密划分。另外，CN201810929678.6基于地形特征界线的无汊河道二维结构网格剖分方法，包括(1)获取基础数据，(2)目标河段二维结构网格初步剖分和(3)对初步剖分结果进行横向加密这三个步骤。该方法通过将河道边界线、深泓线、滩槽分界线及水边线等特征地形界线纳入河道区域的二维结构网格剖分过程，这能够使所剖分的网格能够较好地适应河道边界变化。
[0004]但上述方法均只能实现河道二维网格剖分，无法用于河道三维水流数模拟计算的应用。目前天然河道三维水流数值模拟网格处理方法有：（1）平面正交网线变换结合垂向无量纲变换，该方法需对控制方程进行坐标变换，形式十分复杂，对复杂边界条件的天然河道适应性差；（2）结构化矩性网格，该方法网格结点无法准确布置在天然河道复杂的边界上，为增加计算结果的准确性，需加密网格，故网格单元数量十分庞大，计算成本较高；（3）平面非结构网格结合垂向sigam坐标变换，该方法在一定程度上能适应复杂的边界条件，目前应用最广，但变换后控制方程形式极为复杂，计算结果需加以转换，计算难度增大，还可能引入新的误差，某些情况下容易发散，在地形变化复杂处误较大。
[0005]
技术实现思路

[0006]针对上述现有技术的不足，本专利技术所要解决的技术问题是：怎样提供一种能够更好地实现河道的立体空间剖分的天然河道三维网络生成方法，使其能够更好地直接实现天然河道三维水流数值模拟，无需对水流控制方法进行坐标变换，提高天然河道三维水流数值模拟的准确性和稳定性。
[0007]为了解决上述技术问题，本专利技术采用了如下的技术方案：一种天然河道三维网络生成方法，其特征在于，将天然河道整体模型分解成在水平方向上连续且互不重叠的柱状的子域，然后将各子域规则化后采用网格化软件对各子域
进行网格剖分，再合并剖分后的各子域并构成天然河道三维网格组合体。
[0008]本专利技术基于“离散”的思想，将天然河道整体分解成在水平方向上连续且互不重叠的柱状的子域，这样方便以最少的误差代价将各子域规则化，然后再对规则化后的子域可以很方便地采用软件进行自动网格剖分，最终再合并各子域构成天然河道三维网格组合体。可方便快捷地实现了复杂天然河道的三维网络剖分，具有效率高，生成网格质量好，表征天然河道三维复杂边界特性精确度高的优点。
[0009]进一步地，所述将各子域规则化是指将各子域下底面的河床曲面转化为平面。
[0010]这样，转化为规则体方便后续网格剖分。
[0011]进一步地，采用的网格化软件为Gambit软件。
[0012]Gambit软件是一种可以帮助分析者和设计者建立并网格化计算流体力学（CFD）模型的现有软件，能够直接实现对规则立体区域的网格剖分，功能强大，计算快捷。
[0013]进一步地，本专利技术具体包括以下的步骤：a获取需要计算三维水流数值模拟区域的河道的四条边界线（包括两条边岸线和两条边岸线对应端点之间的连线），构建一个四条边界线所围成的水平面，给定每条边界线的节点数，采用能够沿水平面拼接延伸的多边形，以相邻两个边界节点为多边形的一个边，对四条边界线之间的区域进行二维网格剖分，使其剖分为若干个多边形，导入河床高程，插值得到每个多边形单元中各个节点（即多边形角所在位置）对应河床位置的地形高程，进而得到由若干多边形单元拼接构成的河床底面模型；b以河床底面模型中各多边形单元为下底面，对多边形单元各节点赋予一定的相同的顶高程，形成上顶面，并得到一个柱状的三维立体的子域，逐个将各子域的结点坐标和高程编制成数据文件，每个子域均按照先下底面，后上顶面的方式，上、下面的节点按逆序编排；c将各子域的结点数据文件导入Gambit软件，使用Gambit编程语言编制Journal文件，逐个自动生成线、面和三维的多面体，对各多面体的每条边指定剖分节点数，对多面体自动进行网格剖分,待全部子域网格生成后，使用face connect命令对共面进行合并，最终生成天然河道三维网格组合体。
[0014]这样，本方法能够方便快捷地实现对天然河道的三维网格剖分，生成天然河道三维网格组合体，以满足天然河道三维水流数值模拟的需求，提高了后续三维水流数值模拟的精确性、可靠性和稳定性，同时方法利用了已有的简单模型三维剖分软件实现了对复杂空间的三维剖分，方法自身具有步骤简单，可操作性强，自动化程度高等优点。
[0015]进一步地，a步骤中所述多边形为三角形或者四边形，c步骤中的多面体为五面体或六面体。
[0016]这样采用三角形或者四边形，多边形节点可方便布置在边界上，对平面复杂边界适应性好，建模更加简单，且方便实现沿水平面的拼接延伸，方便实现二维网格剖分。
[0017]进一步地，a步骤中，获取需要计算三维数值区域的天然河道的四条边界线，是根据现有河道地形图获取，保证其数据真实性。
[0018]进一步地，a步骤中，进行二维网格剖分时，多边形的边长范围从数厘米到数米范围内取值。具体地说在不影响计算机处理效率的基础上可取尽量小的范围，取值越小越有利于提高后续河道三维水流数值模拟计算精度。
[0019]进一步地，a步骤中，在二维网格剖分时，对河道内重点关心区域位置的多边形边长范围取值更小，使其多边形分布密度更大。实现平面网格加密效果。
[0020]这样，可以提高重点关心区域的河道三维数值模拟计算精度，而不关心区域的多边形边长可取值较大，总体减少三维单元数量，以提高计算效率。
[0021]进一步地，所述河道内重点关心区域为河道内水利设施建设或河道整治涉及位置区域。
[0022]这样，可以提高工程相关区域的河道三维数值模拟计算精度，使其利于河道工程建设。
[0023]进一步更好的选择是，c步骤中，在垂直方向上边线的节点布置采用非平均分布，在后续计算关联高度位置上分布更密，非关联高度位置上分布更疏。
[0024]这样实现了垂向网格加密，相比于平均分布，能够更好地提高后续对应的应用情况计算的精度。例如在后续涉及河床底部泥沙冲刷情况模拟计算需求下，可以采用上端疏，下端密；涉及河床本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种天然河道三维网络生成方法，其特征在于，将天然河道整体模型分解成在水平方向上连续且互不重叠的柱状的子域，然后将各子域规则化后采用网格化软件对各子域进行网格剖分，再合并剖分后的各子域并构成天然河道三维网格组合体。2.如权利要求1所述的天然河道三维网络生成方法，其特征在于，所述将各子域规则化是指将各子域下底面的河床曲面转化为平面。3.如权利要求1所述的天然河道三维网络生成方法，其特征在于，采用的网格化软件为Gambit软件。4.如权利要求1所述的天然河道三维网络生成方法，其特征在于，本方法具体包括以下的步骤：a获取需要计算三维水流数值区域的天然河道的四条边界线，构建一个四条边界线所围成的水平面，给定每条边界线的节点数，采用能够沿水平面拼接延伸的多边形，以相邻两个边界节点为多边形的一个边，对四条边界线之间的区域进行二维网格剖分，使其剖分为若干个多边形，导入河床高程，插值得到每个多边形单元中各个节点对应河床位置的地形高程，进而得到由若干多边形单元拼接构成的河床底面模型；b以河床底面模型中各多边形单元为下底面，对多边形单元各节点赋予一定的相同的顶高程，形成上顶面，并得到一个柱状的三维立体的子域，逐个将各子域的结点坐标和高程编制成数据文件，每个子域均按照先下底面，后上顶面的方式，上、下面的节点按逆序编排；c将各子域的结点数据文件导入Gambit软件，使用Gambit编程语言编制Journal文件，逐个自动生成线、面和三维的多面体，对各多面体的每条边指定剖分节点数，对多面体自动进行网格剖分,待全部子域网格生成后，使用face connect命令对共面进行合并，最终生成天然河道三维网格组合体。5.如权利要求4所述的天然河道三维网络生成方法，其特征在于，a步骤...

【专利技术属性】
技术研发人员：杨忠超，杨睿，胡雪梅，
申请(专利权)人：重庆交通大学，
类型：发明
国别省市：