CN102567594A - 一种近岸岛礁型人工鱼礁群流场仿真建模的方法 - Google Patents
一种近岸岛礁型人工鱼礁群流场仿真建模的方法 Download PDFInfo
- Publication number
- CN102567594A CN102567594A CN2012100380406A CN201210038040A CN102567594A CN 102567594 A CN102567594 A CN 102567594A CN 2012100380406 A CN2012100380406 A CN 2012100380406A CN 201210038040 A CN201210038040 A CN 201210038040A CN 102567594 A CN102567594 A CN 102567594A
- Authority
- CN
- China
- Prior art keywords
- software
- island
- coordinate
- reef
- numerical value
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Artificial Fish Reefs (AREA)
Abstract
本发明提出一种近岸岛礁型人工鱼礁群流场仿真建模的方法,该方法通过平面插值对海底地形进行网格化,通过SolidWorks软件构建复杂海底起伏地形,作为一数值水槽的底面,构建海平面后得到一个三维立体数值水槽,并通过岛屿形状的建模在此水槽中抠除岛屿陆地部分;在GAMBIT软件中导入上述包含岛屿及人工鱼礁区范围的数值水槽,通过计算机辅助设计对人工鱼礁进行三维立体造型,把鱼礁单体或多鱼礁组合体设置于上述数值水槽内,通过平面插值得出鱼礁位置及其所处水深,在数值水槽底部抠除人工鱼礁固体部分,经过合并后在GAMBIT中对上述数值水槽进行面网格剖分,把面网格数据导入Tgrid软件进行体网格划分,体网格数据导入CFD软件即可进行仿真模拟计算。
Description
技术领域
本发明涉及数值模型建模领域,尤其涉及一种近岸岛礁型人工鱼礁群流场仿真建模的方法。
背景技术
人工鱼礁是投放在海洋中,用于改善海域生态环境、增殖渔业资源的人工构造物,可分为沉底式人工鱼礁和浮式人工鱼礁。进入21世纪,我国沿海各地都在建设不同规模的人工鱼礁项目。人工鱼礁投放在海底后,受到海流的作用力,会在迎流面侧产生上升流,在背流面侧产生背涡流。上升流能将海洋底层的营养盐源源不断地输送到真光层中,从而提高海域的初级生产力,进而提高经济鱼类的增殖能力;背涡流能为一些鱼类提供流速相对较缓的栖息场所;鱼礁投放后,在礁体周围形成强弱流交替的三维流场,流速不等的湍流空间能有效诱集不同的经济鱼类。因此,为使人工鱼礁发挥最大的流场调控效能,达到预期设计目的,对鱼礁投放前后的流场仿真模拟十分重要。
小规模流场的测算模型建立在一个数值水槽中,通常整个水槽的水深相同或变化规则。在实际海域中,水深和岸线复杂多变,特别是目前我国的人工鱼礁较多地建设在近岸多岛礁的海域,这就需要建立具有实际水深变化以及岛礁岸线的数值仿真模型。
发明内容
本发明的目的是提供一种一种近岸岛礁型人工鱼礁群流场仿真建模的方法,在近岸岛礁型人工鱼礁群流场仿真建模过程中,实现符合海域实际海底水深起伏以及复杂岛礁海岸线变化的问题。
本发明为解决其技术问题所采用的技术方案是,
一种近岸岛礁型人工鱼礁群流场仿真建模的方法,包括以下步骤:
a、划定需要进行建模的矩形目标海域,目标海域离开边界距离500m以上,矩形的长边平行于当地潮流主流轴方向,矩形的宽边垂直于当地潮流主流轴方向,对矩形进行网格化,长、宽边的水平网格分辨率相同,划分为m×n节点(假设m≥n),以实际海域的水深和岛屿高程为依据,利用MATLAB软件的插值函数内插得到m×n节点处的水深(平均海面以下,为负值)或岛屿高程(平均海面以上,为正值,可将最大值取为当地的最大潮高),按顺序输出n个*.sldcrv文件,文件中含m行数据,每行3列,包括有东西向的x坐标,南北向的y坐标,以及垂向的z坐标,z坐标尺度对应水深或高程值,为便于后续建模中的视觉判定,将水深或高程值扩大10倍;
b、在计算机辅助设计软件Solidworks中按顺序导入a步骤中所生成的n个*.sldcrv文件,得到n条垂向起伏的等间距线段,逐条连接这n条线段,使其产生一个底面,判断地形的起伏及岛屿分布是否与实际海域相同,拉伸这个底面至当地平均海面高程,得到一个地形起伏、包含岛礁岸线、海面水平的大型数值水槽,通过岛屿形状的建模在该数值水槽中抠除岛屿陆地部分,或者直接根据实际高程使岛屿部分不参与计算,并缩短z坐标尺度为a步骤中的1/10,使z坐标尺度与实际海域相符,导出数值水槽模型文件;
c、利用CFD软件的前处理软件导入b步骤中所导出的数值水槽模型文件,所述CFD软件的前处理软件优选Gambit和Tgrid;
d、在Gambit中生成面网格的基础上,利用Tgrid软件生成用于未投礁海域流场仿真计算的四面体网格;
e、在计算机辅助设计软件Solidworks中建立人工鱼礁三维模型,所述人工鱼礁三维模型的构造及尺寸与实际礁体完全相同;
f、在MATLAB软件中用插值函数内插得到鱼礁投放水域的实际水深,输出每个鱼礁单体对应的x,y,z坐标值,利用Gambit在c步骤所建数值水槽基础上导入e步骤中的人工鱼礁三维模型,将鱼礁单体根据上述坐标移动到位;
g、剔除f步骤中需在三维数值水槽模型内设置的人工鱼礁的三维模型所在区域,得到三维数值水槽模型剩余部分,利用CFD软件的前处理软件对剩余部分进行网格划分,得到网格模型文件;在Gambit中重新生成面网格的基础上,利用Tgrid软件生成最后用于投礁后流场仿真计算的四面体网格。
另外,b步骤中在计算机辅助设计软件Solidworks中按顺序导入a步骤所生成的n个*.sldcrv文件,得到n条垂向起伏的等间距线段,逐条连接这n条线段,使其产生一个底面,为判段地形的起伏及岛屿分布是否与实际海域相同,z坐标数据此前扩大10倍,判断所建地形正确后,将所得到的一个地形起伏、包含岛礁岸线、海面水平的大型数值水槽的垂向尺度缩短为此前的1/10;所述人工鱼礁群为复数个人工鱼礁单体的排列组合。
本发明提出的仿真建模方法可以用于近岸岛礁型人工鱼礁群流场仿真的建模,通过计算机辅助设计软件Solidworks对海底地形和岛礁分布以及人工鱼礁进行三维立体造型,人工鱼礁可以是鱼礁单体或多鱼礁单体的排列组合,把鱼礁单体或多鱼礁单体的排列组合设置于包含起伏海底地形及岛礁的大型数值水槽内。
最后需要剔除b步骤中大型数值水槽模型内的人工鱼礁的三维模型所在区域,然后对三维数值水槽模型进行三维网格剖分,为体现鱼礁体的复杂变化,本发明采用四面体网格,然后把得到的网格数据导入CFD软件,CFD软件优选Fluent软件。
综上所述,本发明应用计算机辅助设计技术来建立近岸岛礁型人工鱼礁群流场仿真模型,可以有效的克服现有建模过程中存在的各种问题,例如海底实际地形起伏变化、岛屿岸线走向复杂等,并且能保证建模的精度优于目前常用的按比例缩小的物理模型等建模方法,使用本发明方法能实现准确数值建模,几何精度可达厘米级。
附图说明
图1为solidworks导入实际海域水深曲线;
图2为solidworks通过放样曲面形成海底地形曲面;
图3为solidworks导入画基准面的两条曲线画基准面;
图4为solidworks通过凸台拉伸形成体;
图5为solidworks导入岛屿地形曲线;
图6为solidworks转换实体引用对岛屿进行拉伸切除;
图7为solidworks通过比例缩放,垂向尺度缩短为此前的1/10;
图8为Gambit导入(import)数值水槽模型的*.step文件;
图9是Gambit将人工鱼礁体与海底形成一体;
图10是Gambit通过split分块划分网格。
具体实施方式
为了使本发明实现的技术手段、创作特征、达成目的与功效易于明白了解,下面结合图示与具体实施例,进一步阐述本发明。
本实施例中,具体建模过程如下:
1.划定需要进行建模的矩形海域4000m×4000m,水平空间分辨率相同,划分为m×n节点(假设m≥n),以实际海域的水深或岛屿高程为依据,利用MATLAB软件的插值函数内插得到m×n节点处的水深或岛屿处的高程,按顺序输出n个*.sldcrv文件中,文件中含m行数据,每行3列(包含东西向x坐标,南北向y坐标,垂向z坐标,即水深或高程);
2.在solidworks内导入1中输出的实际海域水深曲线(图1),通过放样曲面形成海底地形曲面(图2),再导入画基准面(平均海面)的两条曲线画基准面(图3),通过凸台拉伸形成体(图4),然后通过导入岛屿地形曲线(图5),转换实体引用对岛屿进行拉伸切除(图6),最后通过比例缩放,垂向尺度缩短为此前的1/10(图7),导出数值水槽模型的*.step文件;
3.在Gambit下导入(import)数值水槽模型的*.step文件(图8),根据实际海域中投放人工鱼礁的位置和排列方式,在MATLAB等软件中插值函数内插得到鱼礁投放水域的实际水深,输出每个鱼礁单体对应的x,y,z坐标值,将鱼礁单体根据上述坐标移位,再通过抠除(subtract)操作将人工鱼礁体与海底形成一体(图9),为了进一步提高网格的精度,并尽量消除面网格的瓶颈,通过split分块划分网格(图10);
4.三维数值水槽模型各个面进行属性设置:
水流进口面(In):速度进口(velocity-inlet),流向为x轴方向,流速大小设为恒定值(即为人工鱼礁区的实测潮周期平均流速);
水流出口面(Out):速度出口(velocity-inlet),大小为进口面积×进口流速/出口面积,取负值;
单礁体模型表面(Reef):固体边界(wall),并根据礁体和表面粗糙度设置固体边界的糙率;
第一侧面和第二侧面(Side1、Side2):选用对称边界条件,在数值水槽宽度足够大的前提下(侧边界处流场已不受礁体影响),可认为在对称平面上,既无质量的交换,也无热量等物理量的交换,因此,在对称边界上,垂直于边界的速度分量为零,任何量的梯度也为零;
顶部(Top):沉底式鱼礁的投放水域一般同时需保证通航需要,礁体顶部到海面具有一定安全水深距离,因此可忽略礁体对水位的影响,因此选用对称边界(moving wall);
底部:选用固体边界(wall),可根据底质粒径设置固体边界的糙率,由于粗糙程度对单礁体周围整体流态的影响有限,因此对于固体边界的糙率设置的误差不会对最终测算结果造成影响;
5.将Gambit画好的网格导入到Tgrid中,首先对面网格进行修正,使其扭曲度尽量保证在0.5以下,然后生成体网格,为了降低体网格的扭曲度,进行网格加密,使其扭曲度降到0.8以下,在进行网格refine的时候,尽量保证体网格数在400万以内。
以上显示和描述了本发明的基本原理、主要特征和本发明的优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等同物界定。
Claims (3)
1.一种近岸岛礁型人工鱼礁群流场仿真建模的方法,其特征在于,包括以下步骤:
a、划定需要进行建模的矩形目标海域,目标海域离开边界距离500m以上,矩形的长边平行于当地潮流主流轴方向,矩形的宽边垂直于当地潮流主流轴方向,对矩形进行网格化,长、宽边的水平网格分辨率相同,划分为m×n节点,m≥n,以实际海域的水深和岛屿高程为依据,利用MATLAB软件的插值函数内插得到m×n节点处的水深或岛屿高程,按顺序输出n个*.sldcrv文件,文件中含m行数据,每行3列,包括有东西向的x坐标,南北向的y坐标,以及垂向的z坐标,z坐标尺度对应水深或高程值,为便于后续建模中的视觉判定,将水深或高程值扩大10倍;
b、在计算机辅助设计软件Solidworks中按顺序导入a步骤中所生成的n个*.sldcrv文件,得到n条垂向起伏的等间距线段,逐条连接这n条线段,使其产生一个底面,判断地形的起伏及岛屿分布是否与实际海域相同,拉伸这个底面至当地平均海面高程,得到一个地形起伏、包含岛礁岸线、海面水平的大型数值水槽,通过岛屿形状的建模在该数值水槽中抠除岛屿陆地部分,或者直接根据实际高程使岛屿部分不参与计算,并缩短z坐标尺度为a步骤中的1/10,使z坐标尺度与实际海域相符,导出数值水槽模型文件;
c、利用CFD软件的前处理软件导入b步骤中所导出的数值水槽模型文件,所述CFD软件的前处理软件优选Gambit和Tgrid;
d、在Gambit中生成面网格的基础上,利用Tgrid软件生成用于未投礁海域流场仿真计算的四面体网格;
e、在计算机辅助设计软件Solidworks中建立人工鱼礁三维模型,所述人工鱼礁三维模型的构造及尺寸与实际礁体完全相同;
f、在MATLAB软件中用插值函数内插得到鱼礁投放水域的实际水深,输出每个鱼礁单体对应的x,y,z坐标值,利用Gambit在c步骤所建数值水槽基础上导入e步骤中的人工鱼礁三维模型,将鱼礁单体根据上述坐标移动到位;
g、剔除f步骤中需在三维数值水槽模型内设置的人工鱼礁的三维模型所在区域,得到三维数值水槽模型剩余部分,利用CFD软件的前处理软件对剩余部分进行网格划分,得到网格模型文件;在Gambit中重新生成面网格的基础上,利用Tgrid软件生成最后用于投礁后流场仿真计算的四面体网格。
2.根据权利要求1所述的近岸岛礁型人工鱼礁群仿真建模的方法,其特征在于,b步骤中在计算机辅助设计软件Solidworks中按顺序导入a步骤所生成的n个*.sldcrv文件,得到n条垂向起伏的等间距线段,逐条连接这n条线段,使其产生一个底面,为判断地形的起伏及岛屿分布是否与实际海域相同,z坐标数据此前扩大10倍,判断所建地形正确后,将所得到的一个地形起伏、包含岛礁岸线、海面水平的大型数值水槽的垂向尺度缩短为此前的1/10。
3.根据权利要求1所述的一种近岸岛礁型人工鱼礁群仿真建模的方法,其特征在于,所述人工鱼礁群为复数个人工鱼礁单体的排列组合。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210038040.6A CN102567594B (zh) | 2012-02-20 | 2012-02-20 | 一种近岸岛礁型人工鱼礁群流场仿真建模的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210038040.6A CN102567594B (zh) | 2012-02-20 | 2012-02-20 | 一种近岸岛礁型人工鱼礁群流场仿真建模的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102567594A true CN102567594A (zh) | 2012-07-11 |
CN102567594B CN102567594B (zh) | 2014-01-15 |
Family
ID=46412986
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210038040.6A Expired - Fee Related CN102567594B (zh) | 2012-02-20 | 2012-02-20 | 一种近岸岛礁型人工鱼礁群流场仿真建模的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102567594B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105787967A (zh) * | 2015-10-15 | 2016-07-20 | 上海海洋大学 | 一种复杂地形岛礁海域海洋牧场区建设面积的测算方法 |
CN112132927A (zh) * | 2020-09-15 | 2020-12-25 | 成都工具研究所有限公司 | 在网页界面上生成二维及三维模型的绘图系统及方法 |
CN113325424A (zh) * | 2021-06-02 | 2021-08-31 | 上海海洋大学 | 根据多波束水深数据估算人工鱼礁空方量的方法 |
CN115994410A (zh) * | 2023-03-22 | 2023-04-21 | 中国人民解放军国防科技大学 | 基于八叉树细化四面体网格的飞行器仿真驱动设计方法 |
CN116702476A (zh) * | 2023-06-09 | 2023-09-05 | 中国人民解放军92493部队试验训练总体研究所 | 一种基于气囊式结构的岛岸靶区仿真模型验证方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH08124A (ja) * | 1994-06-20 | 1996-01-09 | Ishikawajima Harima Heavy Ind Co Ltd | 人工魚礁 |
CN101694678A (zh) * | 2009-10-27 | 2010-04-14 | 上海海洋大学 | 一种测算人工鱼礁流场调控规模的方法 |
CN102169523A (zh) * | 2011-05-18 | 2011-08-31 | 中国水产科学研究院东海水产研究所 | 一种人工鱼礁结构优化设计方法 |
-
2012
- 2012-02-20 CN CN201210038040.6A patent/CN102567594B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH08124A (ja) * | 1994-06-20 | 1996-01-09 | Ishikawajima Harima Heavy Ind Co Ltd | 人工魚礁 |
CN101694678A (zh) * | 2009-10-27 | 2010-04-14 | 上海海洋大学 | 一种测算人工鱼礁流场调控规模的方法 |
CN102169523A (zh) * | 2011-05-18 | 2011-08-31 | 中国水产科学研究院东海水产研究所 | 一种人工鱼礁结构优化设计方法 |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105787967A (zh) * | 2015-10-15 | 2016-07-20 | 上海海洋大学 | 一种复杂地形岛礁海域海洋牧场区建设面积的测算方法 |
CN112132927A (zh) * | 2020-09-15 | 2020-12-25 | 成都工具研究所有限公司 | 在网页界面上生成二维及三维模型的绘图系统及方法 |
CN112132927B (zh) * | 2020-09-15 | 2023-10-20 | 成都工具研究所有限公司 | 在网页界面上生成二维及三维模型的绘图系统及方法 |
CN113325424A (zh) * | 2021-06-02 | 2021-08-31 | 上海海洋大学 | 根据多波束水深数据估算人工鱼礁空方量的方法 |
CN113325424B (zh) * | 2021-06-02 | 2023-07-25 | 上海海洋大学 | 根据多波束水深数据估算人工鱼礁空方量的方法 |
CN115994410A (zh) * | 2023-03-22 | 2023-04-21 | 中国人民解放军国防科技大学 | 基于八叉树细化四面体网格的飞行器仿真驱动设计方法 |
CN116702476A (zh) * | 2023-06-09 | 2023-09-05 | 中国人民解放军92493部队试验训练总体研究所 | 一种基于气囊式结构的岛岸靶区仿真模型验证方法 |
CN116702476B (zh) * | 2023-06-09 | 2023-12-22 | 中国人民解放军92493部队试验训练总体研究所 | 一种基于气囊式结构的岛岸靶区仿真模型验证方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102567594B (zh) | 2014-01-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103514631B (zh) | 基于钻孔数据的海量三维地质模型网格式并行构建方法 | |
CN109271672B (zh) | 一种河-湖-泵站相互影响作用下的河道水面线计算方法 | |
CN102521863B (zh) | 一种基于粒子系统的三维流体标矢量统一动态表示方法 | |
CN102567594B (zh) | 一种近岸岛礁型人工鱼礁群流场仿真建模的方法 | |
CN101694678B (zh) | 一种测算人工鱼礁流场调控规模的方法 | |
Ge et al. | A dike–groyne algorithm in a terrain-following coordinate ocean model (FVCOM): Development, validation and application | |
CN104933268B (zh) | 一种基于一维非恒定流数值模型的洪水分析方法 | |
CN103278115B (zh) | 一种基于dem计算淤地坝淤积量的方法及系统 | |
CN106372290B (zh) | 一种v型河道库容量算方法 | |
CN106682262A (zh) | 一种获取飞行器流场的数值模拟方法 | |
CN104331599A (zh) | 一种非结构化网格嵌套波浪数值模拟方法 | |
CN112784502B (zh) | 一种水文水力学动态双向耦合的洪水预测方法 | |
CN110929443B (zh) | 一种基于高精度地形概化的二维洪水模拟方法 | |
CN104091065A (zh) | 一种求解浅水问题模拟间断水流数值的方法 | |
CN104750972A (zh) | 复杂水网地区楔形水体容积的计算方法 | |
CN109960838A (zh) | 一种体现河流基本特征的河道地形自动生成方法 | |
CN103256914B (zh) | 一种基于dem计算淤地坝淹没面积的方法及系统 | |
CN106844963B (zh) | 模拟开挖至运行全过程的拱坝三维网格模型自动剖分方法 | |
Pan et al. | Impacts of deep waterway project on morphological changes within the north passage of the Changjiang Estuary, China | |
CN103207410B (zh) | 一种针对崎岖海底的混合网格模型建立方法 | |
CN110147646B (zh) | 一种数值模拟框架下线性挡水构筑物的过流处理方法 | |
CN106326588B (zh) | 一种bim辅助大坝仓面设计的方法 | |
CN103412991B (zh) | 一种网箱对水流影响的模拟方法 | |
CN116204956A (zh) | 竖缝式鱼道池室隔板结构优化方法和装置 | |
CN114925624A (zh) | 一种天然河道三维水流数值模拟方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20140115 Termination date: 20150220 |
|
EXPY | Termination of patent right or utility model |