CN111967082B - 一种波浪对结构物作用的模拟分析方法 - Google Patents
一种波浪对结构物作用的模拟分析方法 Download PDFInfo
- Publication number
- CN111967082B CN111967082B CN202010851096.8A CN202010851096A CN111967082B CN 111967082 B CN111967082 B CN 111967082B CN 202010851096 A CN202010851096 A CN 202010851096A CN 111967082 B CN111967082 B CN 111967082B
- Authority
- CN
- China
- Prior art keywords
- pressure
- chamber
- closed
- wave
- moment
- 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.)
- Active
Links
- 238000004088 simulation Methods 0.000 title claims abstract description 27
- 230000009471 action Effects 0.000 title claims abstract description 15
- 238000004458 analytical method Methods 0.000 title claims abstract description 10
- 238000004364 calculation method Methods 0.000 claims abstract description 18
- 238000000034 method Methods 0.000 claims abstract description 18
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 26
- 230000008569 process Effects 0.000 claims description 8
- 230000003993 interaction Effects 0.000 claims description 7
- 239000002245 particle Substances 0.000 claims description 5
- 230000002459 sustained effect Effects 0.000 claims description 3
- 239000012530 fluid Substances 0.000 abstract description 5
- 238000007789 sealing Methods 0.000 abstract description 5
- 239000011148 porous material Substances 0.000 description 4
- 230000009194 climbing Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 230000021715 photosynthesis, light harvesting Effects 0.000 description 2
- 230000005514 two-phase flow Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 230000000254 damaging effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Mathematical Optimization (AREA)
- General Engineering & Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Computing Systems (AREA)
- Fluid Mechanics (AREA)
- Mathematical Physics (AREA)
- Algebra (AREA)
- Architecture (AREA)
- Civil Engineering (AREA)
- Structural Engineering (AREA)
- Computational Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供了一种波浪对结构物作用的模拟分析方法，属于海洋工程技术领域。波浪对结构物作用时形成的腔室满足封闭条件时，持续封闭次数i＝i+1，判断i是否大于1，若否，计算此封闭时刻t腔室内封闭气体体积Vt+1，并判断是否满足时间终止条件；若是，计算t+i时刻腔室内封闭气体体积为Vt+i，利用理想气体状态方程，求解t+i时刻腔室内封闭气体压力Pt+i；利用压力Pt+i更新封闭腔室内自由表面压力，再利用更新后的压力更新计算域内流体速度，并判断是否满足时间终止条件。本发明可以更加真实的模拟波浪在结构物附近的运动特征和流场特性，能够较大程度提高计算效率，节约计算成本。
Description
技术领域
本发明属于海洋工程技术领域，具体涉及一种波浪对结构物作用的模拟分析方法。
背景技术
海岸结构物常被用来抵御外海波浪的冲击，并通过反射消能和紊动消能等形式削减波浪对掩护区域的破坏作用，保护海岸及沿岸设施与生命安全。结构物水动力性能及结构物附近的流场变化特性分析，一直是研究海岸结构物安全性和稳定性的重要内容。近年来，随着计算流体力学数值模拟技术的不断成熟，精确性和可靠性的不断提高，数值模拟逐渐成为分析波浪与结构物相互作用的有效手段。对一般模拟波浪与结构物的相互作用而言，具有单相流(水)模型和两相流(空气+水)模型。传统单相流模型通常不考虑空气的作用，通过在结构物表面设置边界条件，求解空间离散后模拟过程中各时刻的速度和压力，实现流体与结构物的相互作用模拟。单相流模型能够较好的预测波浪的水动力特性，且计算快速、计算稳定性高。但当波浪冲击结构物时，有时会在结构物周围形成封闭腔室，如波浪在开孔沉箱结构上的爬高，波浪对透空式结构物底板的冲击等。在封闭腔室内，需要考虑空气压缩变化对水体运动的影响。单相流模型因不考虑空气作用，无法有效实现封闭腔室内流体运动的准确模拟。双相流模型虽然考虑了腔室内空气变化对水体运动的影响，但计算时间长、计算成本高。
发明内容
针对现有技术中存在不足，本发明提供了一种波浪对结构物作用的模拟分析方法，能够显著提高计算效率，节约计算成本。
本发明是通过以下技术手段实现上述技术目的的。
一种波浪对结构物作用的模拟分析方法，包括如下步骤：
步骤(1)，建立波浪数值水槽，模拟波浪与结构物相互作用；
步骤(2)，根据波浪对结构物的作用特征，设置压力更新条件，所述设置压力更新条件为判断波浪对结构物作用形成的腔室是否封闭；并标记持续封闭次数i＝0；
步骤(3)，求解得到模拟过程中各时刻波浪水质点的速度和压力；
步骤(4)，判断所形成的腔室是否封闭，满足腔室封闭条件时，持续封闭次数i＝i+1，执行步骤(5)；不满足腔室封闭条件时，i＝0，执行步骤(6)；
步骤(5)，判断i是否大于1；若否，计算此封闭时刻t腔室内封闭气体体积Vt+1，执行步骤(6)；若是，计算t+i时刻腔室内封闭气体体积为Vt+i，利用常温下理想气体状态方程Pt+ iVt+i＝PVt+1，求解t+i时刻腔室内封闭气体压力Pt+i；利用压力Pt+i更新封闭腔室内自由表面压力，再利用更新后的压力更新计算域内水质点速度，执行步骤(6)；其中P为每次初始封闭时刻腔室内封闭气体压力，一般选取1个标准大气压；
步骤(6)，判断是否满足时间终止条件t>tend，若满足，则结束计算；若不满足，t＝t+△t，返回执行步骤(3)。
进一步的技术方案，所述波浪数值水槽的建立是通过处理边界条件实现的。
进一步的技术方案，所述模拟过程中各时刻水质点的速度和压力，是利用所建立的波浪数值水槽和求解Navier-Stokes方程得到。
进一步的技术方案，所述设置压力更新条件还能够判断空气是否被压缩。
进一步的技术方案，所述Vt+i通过积分腔室内空气体积获得。
本发明的有益效果为：本发明以传统单相流模型为基础，根据波浪冲击结构物时形成腔室的特点，考虑波浪在腔室内的运动特征，引入理想气体状态方程和腔室封闭条件判断标准，计算得到腔室封闭时刻的腔室内空气压力，更新腔室内自由表面压力，可以更加真实的模拟波浪在结构物附近的波面特征和流场特性，因计算过程仍然属于单相流模型计算，能够较大程度提高计算效率，节约计算成本。
附图说明
图1为本发明所述波浪对结构物作用的模拟分析方法流程图；
图2为本发明所述开孔沉箱结构示意图；
图3为传统单相流模型模拟结果图；
图4为本发明模拟结果图；
图5为数值模拟结果与试验结果对比图。
具体实施方式
下面结合附图以及具体实施例对本发明作进一步的说明，但本发明的保护范围并不限于此。
如图1所示，一种波浪对结构物作用的模拟分析方法，包括如下步骤：
步骤(1)，处理边界条件，建立波浪数值水槽，在波浪数值水槽中进行波浪与结构物相互作用模拟；
步骤(2)，根据波浪对结构物的作用特征，设置压力更新条件，即判断结构物所形成的腔室是否封闭或者空气是否被压缩，本实施例优选判断结构物所形成的腔室是否封闭；并标记持续封闭次数i＝0；
步骤(3)，利用建立的波浪数值水槽和N-S方程(Navier-Stokes，纳维-斯托克斯方程)，求解得到模拟过程中各时刻波浪水质点的速度和压力；
步骤(4)，判断所形成的腔室是否封闭，满足腔室封闭条件时，持续封闭次数i＝i+1，执行步骤(5)；不满足腔室封闭条件时，i＝0，执行步骤(6)；
步骤(5)，判断i是否大于1，当i≤1时，计算此封闭时刻t腔室内封闭气体体积Vt+1(通过积分腔室内空气体积获得)，执行步骤(6)；当i＞1时，计算t+i时刻腔室内封闭气体体积为Vt+i(通过积分腔室内空气体积获得)，基于常温下理想气体状态方程PV＝C(常数)，本发明中的理想气体状态方程为Pt+iVt+i＝PVt+1，求解t+i时刻腔室内封闭气体压力Pt+i；利用压力Pt+i更新封闭腔室内自由表面压力(即利用压力Pt+i更新水气交界面压力)，再利用更新后的压力更新计算域内水质点速度，执行步骤(6)；其中P为初始封闭时刻腔室内封闭气体压力(一般选取1个标准大气压)；
步骤(6)，判断是否满足时间终止条件t>tend，若满足，则结束计算；若不满足，t＝t+Δt，返回执行步骤(3)。
本实施中结构物以开孔沉箱结构为例，具体描述波浪对开孔沉箱结构作用的模拟分析方法，具体如下：
开孔沉箱结构因消浪室的存在(图2)，波浪通过开孔前墙进入消浪室内，利用消浪室内外波浪相位差和消浪室内水体的紊动耗能，能够有效降低结构反射系数、波浪力和越浪量。波峰作用于开孔沉箱时，波浪沿开孔沉箱开孔前墙爬高，将爬高高度记为hi(i＝1，2，…n)，开孔前墙最高开孔位置记为hm；当h1>hm时，判定消浪室处于封闭状态，记录为每次封闭状态的初始时刻t1＝t，计算t1时刻消浪室内空气体积为Vt+1；当t2＝t+Δt，波浪继续爬高，h2>hm，消浪室仍处于封闭状态，计算此时消浪室内空气体积为Vt+2，利用理想气体状态方程求得压力Pt+2，更新封闭腔室内自由表面压力，再利用自由表面压力更新计算域内流体速度；当ti＝t+(i-1)Δt时，爬高结束，开孔前墙前的水体落回水槽过程中，如消浪室仍然满足封闭条件，即hi>hm，继续执行计算t+i时刻消浪室内的空气体积Vt+i和压力Pt+i，更新封闭腔室内自由表面压力，再利用更新后的压力更新计算域内流体速度；直到hn<hm，不满足消浪室封闭条件，不再执行压力更新。一个周期结束，继续执行下一个周期，重复上述判断过程，直到满足终止条件。
波浪与开孔沉箱相互作用的单相流模型模拟结果和本发明的模拟结果如图3-4所示，模拟结果与试验结果的对比，见图5。
从图3和图4中可以看出：传统单相流模型(图3)，因忽略了消浪室内的空气存在，水体近乎充满整个消浪室。而当本发明执行压力更新后，消浪室内考虑了空气的存在(图4)，自由表面与消浪室顶部存在空腔。图5为单相流模型模拟结果、本发明模拟结果与试验结果对比，可以发现：单相流模型模拟结果与试验结果存在偏差；本发明模拟结果与试验结果接近一致，本发明可以有效的模拟水体在消浪室内的运动。
所述实施例为本发明的优选的实施方式，但本发明并不限于上述实施方式，在不背离本发明的实质内容的情况下，本领域技术人员能够做出的任何显而易见的改进、替换或变型均属于本发明的保护范围。
Claims (1)
1.一种波浪对结构物作用的模拟分析方法，其特征在于，包括如下步骤：
步骤(1)，建立波浪数值水槽，模拟波浪与结构物相互作用；所述波浪数值水槽的建立是通过处理边界条件实现的；
步骤(2)，根据波浪对结构物的作用特征，设置压力更新条件，所述设置压力更新条件为判断波浪对结构物作用形成的腔室是否封闭；并标记持续封闭次数i＝0；所述设置压力更新条件还能够判断空气是否被压缩；
步骤(3)，求解得到模拟过程中各时刻波浪水质点的速度和压力；所述模拟过程中各时刻水质点的速度和压力，是利用所建立的波浪数值水槽和求解Navier-Stokes方程得到；
步骤(4)，判断所形成的腔室是否封闭，满足腔室封闭条件时，持续封闭次数i＝i+1，执行步骤(5)；不满足腔室封闭条件时，i＝0，执行步骤(6)；
步骤(5)，判断i是否大于1；若否，计算此封闭时刻t腔室内封闭气体体积Vt+1，执行步骤(6)；若是，计算t+i时刻腔室内封闭气体体积为Vt+i，利用常温下理想气体状态方程Pt+iVt+i＝PVt+1，求解t+i时刻腔室内封闭气体压力Pt+i；利用压力Pt+i更新封闭腔室内自由表面压力，再利用更新后的压力更新计算域内水质点速度，执行步骤(6)；其中P为每次初始封闭时刻腔室内封闭气体压力，P值选取为1个标准大气压；所述Vt+i通过积分腔室内空气体积获得；
步骤(6)，判断是否满足时间终止条件t>tend，若满足，则结束计算；若不满足，t＝t+△t，返回执行步骤(3)。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010851096.8A CN111967082B (zh) | 2020-08-21 | 2020-08-21 | 一种波浪对结构物作用的模拟分析方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010851096.8A CN111967082B (zh) | 2020-08-21 | 2020-08-21 | 一种波浪对结构物作用的模拟分析方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111967082A CN111967082A (zh) | 2020-11-20 |
| CN111967082B true CN111967082B (zh) | 2023-12-26 |
Family
ID=73389970
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010851096.8A Active CN111967082B (zh) | 2020-08-21 | 2020-08-21 | 一种波浪对结构物作用的模拟分析方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111967082B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113361217B (zh) * | 2021-07-07 | 2022-10-11 | 中国海洋大学 | 一种高效的两相流无网格数值模型实施方法、装置 |
| CN113552026B (zh) * | 2021-07-22 | 2022-12-13 | 河海大学 | 一种非侵入式检测波致底床液化的振荡水槽的使用方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2003331208A (ja) * | 2002-05-09 | 2003-11-21 | Toshiba Corp | 流体シミュレーション装置、流体シミュレーション方法、プログラム |
| CN108053736A (zh) * | 2017-11-15 | 2018-05-18 | 浙江海洋大学 | 一种海洋结构物水动力特性研究用实验装置 |
| CN108168841A (zh) * | 2018-02-09 | 2018-06-15 | 水利部交通运输部国家能源局南京水利科学研究院 | 一种臂式土工离心机中的波浪力模拟装置 |
| CN108287965A (zh) * | 2018-01-25 | 2018-07-17 | 鲁东大学 | 一种计算波浪与不规则结构物相互作用的方法 |
-
2020
- 2020-08-21 CN CN202010851096.8A patent/CN111967082B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2003331208A (ja) * | 2002-05-09 | 2003-11-21 | Toshiba Corp | 流体シミュレーション装置、流体シミュレーション方法、プログラム |
| CN108053736A (zh) * | 2017-11-15 | 2018-05-18 | 浙江海洋大学 | 一种海洋结构物水动力特性研究用实验装置 |
| CN108287965A (zh) * | 2018-01-25 | 2018-07-17 | 鲁东大学 | 一种计算波浪与不规则结构物相互作用的方法 |
| CN108168841A (zh) * | 2018-02-09 | 2018-06-15 | 水利部交通运输部国家能源局南京水利科学研究院 | 一种臂式土工离心机中的波浪力模拟装置 |
Non-Patent Citations (2)
| Title |
|---|
| 海上浮式结构物的水动力仿真研究;韩翔希;邱昂;范少涛;王海璎;;船舶工程(第04期);全文 * |
| 海浪环境对航行体出水特性影响研究;周敬国;权晓波;程少华;;导弹与航天运载技术(第03期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111967082A (zh) | 2020-11-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111967082B (zh) | 一种波浪对结构物作用的模拟分析方法 | |
| Zhang et al. | DualSPHysics: a numerical tool to simulate real breakwaters | |
| Chaplin et al. | Ringing of a vertical cylinder in waves | |
| Shao et al. | SPH computation of plunging waves using a 2‐D sub‐particle scale (SPS) turbulence model | |
| Song et al. | Boundary element study of wave impact on a vertical wall with air entrapment | |
| CN110705188A (zh) | 一维冰水耦合运动的高精度格式模拟方法 | |
| Shih et al. | Experimental study on the optimization of an oscillating water column structure considering real gas compressibility | |
| Horko | CFD optimisation of an oscillating water column wave energy converter | |
| CN112836440B (zh) | 一种非稳态聚焦波浪生成方法 | |
| CN117094247A (zh) | 基于人工造波的数值波浪水槽模型的建立方法及其应用 | |
| Ning et al. | Hydrodynamic investigation of a dual-chamber OWC Wave Energy Converter | |
| CN116362001A (zh) | 基于滑体质心高度势能计算高位滑坡入水速度及涌浪方法 | |
| Rahman et al. | An eddy viscosity model with near‐wall modifications | |
| Fomin et al. | Simulation and analysis of sea floods in the Don River Delta | |
| CN112861451A (zh) | 突涌水优势通道扩展过程的渗流演化cfd计算方法 | |
| CN114492229B (zh) | 一种基于Fluent软件的原油储罐泄漏模拟方法 | |
| Iele et al. | Boundary conditions for the simulation of wave breaking | |
| CN111199069A (zh) | 基于门限回归理论的面板堆石坝坝顶沉降经验预测方法 | |
| Uchiyama | Wetting and drying scheme for POM and its applications to San Francisco Bay | |
| Thampi et al. | Filter approach to ocean structure response prediction | |
| Li et al. | CFD and full scale examination of the pressure distribution around a building corner | |
| Yokoyama et al. | Study on statistical prediction method of storm surges in seto inland sea | |
| Abualtayef et al. | Three-dimensional hydrostatic modeling of a bay coastal area | |
| Pawitan et al. | MULTIPHASE SPH SIMULATION FOR INCHAMBER IMPACT PRESSURE ON VERTICAL BREAKWATER WITH WAVE ABSOPTION CHAMBER | |
| CN115952751A (zh) | 一种swan模型的风速修正方法及系统 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |