CN113935114A - 一种用于水上飞机着水数值仿真的规则波造波方法 - Google Patents

一种用于水上飞机着水数值仿真的规则波造波方法 Download PDF

Info

Publication number
CN113935114A
CN113935114A CN202111242287.5A CN202111242287A CN113935114A CN 113935114 A CN113935114 A CN 113935114A CN 202111242287 A CN202111242287 A CN 202111242287A CN 113935114 A CN113935114 A CN 113935114A
Authority
CN
China
Prior art keywords
seaplane
water
model
regular wave
simulation
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.)
Pending
Application number
CN202111242287.5A
Other languages
English (en)
Inventor
杨仕福
吕继航
曾毅
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
South China Aircraft Industry Co Ltd of China Aviation Industry General Aircraft Co Ltd
Original Assignee
South China Aircraft Industry Co Ltd of China Aviation Industry General Aircraft Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by South China Aircraft Industry Co Ltd of China Aviation Industry General Aircraft Co Ltd filed Critical South China Aircraft Industry Co Ltd of China Aviation Industry General Aircraft Co Ltd
Priority to CN202111242287.5A priority Critical patent/CN113935114A/zh
Publication of CN113935114A publication Critical patent/CN113935114A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/15—Vehicle, aircraft or watercraft design
    • 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/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • 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
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling
    • 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)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Computational Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Algebra (AREA)
  • Computing Systems (AREA)
  • Fluid Mechanics (AREA)
  • Mathematical Physics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种用于水上飞机着水数值仿真的规则波造波方法。包括下述步骤：a.构建水上飞机静水着水数值仿真模型；b.在所述的静水着水仿真模型上设计初始规则波波浪形态；c.对已设计初始规则波波浪形态的静水着水仿真模型内的质点进行动能赋值，形成水上飞机规则波着水仿真模型。本发明能够大幅提高水上飞机在规则波上的着水数值仿真计算效率和精度。

Description

一种用于水上飞机着水数值仿真的规则波造波方法
技术领域
本发明属于结构物入水流固耦合数值仿真中液体波浪动力特性仿真计算技术领域，具体涉及一种用于水上飞机着水数值仿真的规则波造波方法。
背景技术
水上飞机着水数值仿真分析，属于结构物入水流固耦合仿真分析的一种。在分析技术上，经历多年的发展，由开始的水上飞机静水着水数值仿真分析，发展到规则波着水数值仿真分析、不规则波着水数值仿真分析，模型建模越来越复杂，网格数量越来越大，以至于计算量急剧增加，数值仿真的实用性受到很大制约。
水上飞机规则波/不规则波着水数值仿真分析，目前的建模方法主要参考摇板造波、推板造波的理论方法，利用摇波板或推板的往返运动，推动液体的运动进而形成波浪，然后进行水上飞机着水数值仿真分析。该造波方法的缺点明显：造波耗时巨大；得到的波浪参数误差较大；造波计算过程中累计误差影响容易中断计算。
水上飞机着水性能的评估，主要参考水上飞机在规则波上的着水性能进行评估。因此，水上飞机规则波着水数值仿真分析，对促进水上飞机的发展起到重大作用。
发明内容
本发明的目的是：提供了一种用于水上飞机着水数值仿真的规则波造波方法。本发明能够大幅提高水上飞机在规则波上的着水数值仿真计算效率和精度。
本发明的技术方案是：一种用于水上飞机着水数值仿真的规则波造波方法，包括下述步骤：
a.构建水上飞机静水着水数值仿真模型；
b.在所述的静水着水仿真模型上设计初始规则波波浪形态；
c.对已设计初始规则波波浪形态的静水着水仿真模型内的质点进行动能赋值，形成水上飞机规则波着水仿真模型。
前述的水上飞机着水数值仿真的规则波造波方法的步骤a中，水上飞机静水着水数值仿真模型的构建原则为：考虑飞机在波浪水面着水后的滑行距离，建立“固-液-气”耦合水上飞机静水着水数值仿真模型。
前述的水上飞机着水数值仿真的规则波造波方法的步骤a中，水上飞机静水着水数值仿真模型的构建方法为：
依据飞机总体气动外形，划分飞机结构有限元模型，不考虑水上飞机机体结构弹性特性，采用刚性体材料属性；
依据水上飞机重量特性，给飞机赋质量特性数据，定义“飞机-水体-空气”的耦合关系，定义流体边界条件；
定义水上飞机初始着水速度和输出项，采用ALE法或EULER法，完成水上飞机静水着水数值仿真计算模型的构建。
前述的水上飞机着水数值仿真的规则波造波方法中，所述的水上飞机静水着水数值仿真模型包括飞机结构模型、水流体模型和空气模型；所述的飞机结构模型的网格单元为拉格朗日单元；水流体模型和空气模型的网格单元均为ALE单元或EULER单元。
前述的水上飞机着水数值仿真的规则波造波方法的步骤b中，初始规则波波浪形态的设计方法为：
根据水上飞机静水着水数值仿真分析模型生成的仿真计算输入文件和规则波自由液面方程(1)，使模型中网格单元的节点坐标发生移动，形成初始规则波波浪形态；
η(x,t)＝a cos(kx-ωt) (1)
式中，η为距离平衡面的位移；a表示微幅波振幅；k表示波数；ω为波动圆频率；t为时间；x为站位。
前述的水上飞机着水数值仿真的规则波造波方法的步骤c中，质点的动能赋值方法如下：
基于初始规则波波浪形态，得到的新的网格单元位置信息，结合式(2)和式(3)，计算新的网格单元的节点速度值；
z为水深，g为重力加速度，λ为波长，则有：
k＝2π/λ＝ω2/g (2)
u为水质点水平方向速度，v为水质点垂直方向速度，规则波的水质点速度为：
Figure BDA0003319630740000031
前述的水上飞机着水数值仿真的规则波造波方法的步骤a中，仿真模型构建时，通过优化建模参数，使飞机的静水着水仿真计算精度满足工程需求。
前述的水上飞机着水数值仿真的规则波造波方法中，仿真模型构建时，优化建模参数的方法为：将模型中处于规则波波浪液面附近及水上飞机着水位置的网格单元的尺寸细化。
本发明的优点是：本发明基于数值造波在水上飞机规则波着水数值仿真中的重要作用，设计了一种用于水上飞机规则波着水数值仿真计算中的规则波造波方法，该方法相对于常用的摇板造波、推板造波数值仿真方法，取得了如下有益效果：
1)流体网格数量减少30％～50％，没有了摇板造波、推板造波模型的造波区和流场末端的削波区。
2)计算效率得到大幅度的提高，本发明水上飞机规则波着水数值仿真计算时间主要花在数值造波上，水上飞机着水数值仿真计算只有1～3秒左右时间，而摇板造波、推板造波从开始造波到形成稳定的波浪需要1分钟多，计算时间减少约80％。
3)计算精度大幅提高，利用本发明生成的规则波，波形和波内能直接对应理论模型，而摇板造波和推板造波理论基于微幅波假设，当遇到波高较大的规则波时，其造出来的波浪与实际存在一定误差，需要修正。
4)减少大量的模型调试时间，即减少了摇板造波或推板造波模型调试时间，在常用的水上飞机规则波着水数值仿真计算中，模型调试时间可占到总计算时间的20％以上。
综上所述，本发明能够大幅提高水上飞机在规则波上的着水数值仿真计算效率和精度。
附图说明
图1为规则波造波流程示意图，描述了水上飞机规则波着水数值仿真计算流程。
图2为水上飞机静水着水数值仿真模型示意图，其中，1—空气网格，2—水面，3—着水区域网格加密，4—水上飞机。
图3为程序计算处理后得到的局部放大规则波形状图。
图4为程序计算处理后得到的水上飞机规则波着水数值仿真模型示意图。
具体实施方式
下面结合附图和实施例对本发明作进一步的说明，但并不作为对本发明限制的依据。
实施例1。一种用于水上飞机着水数值仿真的规则波造波方法，如图1-4所示，具体步骤如下：
步骤一、依据水上飞机总体外形、着水速度和所需计算的水上滑行距离，建立“固-液-气”耦合数值仿真模型，水上飞机网格单元为拉格朗日单元，水和空气流体网格单元为ALE单元或EULER单元，优化着水区域网格尺寸；
步骤二、依据水上飞机重量特性，给水上飞机赋质量特性数据，定义“飞机-水体-空气”的耦合关系，定义流体边界条件；
步骤三、定义水上飞机初始着水速度和输出项(主要为加速度、位移，耦合压力输出)，采用ALE法或EULER法，建立水上飞机静水着水数值仿真计算模型；
步骤四、水上飞机规则波水面着水参数分析，给出描述规则波的数学方程，包括自由液面方程、质点运动速度方程；
步骤五、利用规则波液面方程和水上飞机静水着水数值仿真模型，编程计算并调整仿真模型节点位置，使形成初始规则波水面形态；
步骤六、利用规则波质点运动方程步骤五得到的水上飞机着水数值仿真模型，编程计算并给水质点赋予初始速度，使水流体具有初始动能；
步骤七、在步骤五和步骤六前提下，综合所有对水上飞机静水着水数值仿真模型的更改，得到水上飞机规则波着水数值仿真模型，省去了利用摇板或推板进行造波的仿真计算；
步骤八、提交商用有限元软件计算，得到水上飞机规则波着水数值仿真计算结果。
实施例2。一种用于水上飞机规则波着水数值仿真计算中的规则波造波方法，参见图1-4，结合水上飞机静水着水数值仿真建模方法特点和规则波数学描述方法，具体步骤如下：
步骤一：根据需求建立水上飞机静水着水数值仿真模型
考虑水上飞机在波浪水面着水后的滑行距离，建立“固-液-气”耦合水上飞机静水着水数值仿真模型，仿真模型包括水上飞机结构模型、水流体模型和空气模型。通过优化建模参数，使水上飞机在静水着水数值仿真计算精度满足工程需求。
步骤二：水上飞机规则波着水数值仿真模型中初始规则波波浪形态计算
根据水上飞机静水着水数值仿真分析模型生成的仿真计算输入文件，和规则波自由液面方程式(1)，利用编写的程序，使流体模型中网格节点坐标发生移动，形成规则波浪的水面形态。
规则波自由表面方程为
η(x,t)＝a cos(kx-ωt) (1)
式中，η为距离平衡面的位移，a表示微幅波振幅，k表示波数，ω为波动圆频率，t为时间，x为站位。
步骤三：水上飞机着水数值仿真模型规则波内部质点动能赋值
跟据步骤2得到的新的流体模型网格点位置信息和式(2)和式(3)，利用编写的程序，计算波浪内部流体网格节点速度值，生成水体模型的初始速度。z为水深，g为重力加速度，λ为波长，则有：
k＝2π/λ＝ω2/g (2)
u为水质点水平方向速度，v为水质点垂直方向速度，规则波的水质点速度为：
Figure BDA0003319630740000061
步骤四：生成水上飞机规则波着水数值仿真模型
综合步骤三和步骤四对水上飞机静水着水数值仿真模型的修改，利用编写的程序生成水上飞机规则波着水数值仿真模型，该模型中水流体已具有初始规则波面形态和初始动能，模拟了规则波的初始运动特性。
实施例3。一种用于水上飞机着水数值仿真的规则波造波方法，参见图1-图4，包括下述步骤：
a.构建水上飞机静水着水数值仿真模型；
b.在所述的静水着水仿真模型上设计初始规则波波浪形态；
c.对已设计初始规则波波浪形态的静水着水仿真模型内的质点进行动能赋值，形成水上飞机规则波着水仿真模型。
前述的步骤a中，水上飞机静水着水数值仿真模型的构建原则为：考虑飞机在波浪水面着水后的滑行距离，建立“固-液-气”耦合水上飞机静水着水数值仿真模型。
前述的步骤a中，水上飞机静水着水数值仿真模型的构建方法为：
依据飞机总体气动外形，划分飞机结构有限元模型，不考虑水上飞机机体结构弹性特性，采用刚性体材料属性；
依据水上飞机重量特性，给飞机赋质量特性数据，定义“飞机-水体-空气”的耦合关系，定义流体边界条件；
定义水上飞机初始着水速度和输出项(主要为加速度、位移，耦合压力输出)，采用ALE法或EULER法，完成水上飞机静水着水数值仿真计算模型的构建。
前述的水上飞机静水着水数值仿真模型包括飞机结构模型、水流体模型和空气模型；所述的飞机结构模型的网格单元为拉格朗日单元；水流体模型和空气模型的网格单元均为ALE单元或EULER单元。
前述的步骤b中，初始规则波波浪形态的设计方法为：
根据水上飞机静水着水数值仿真分析模型生成的仿真计算输入文件和规则波自由液面方程(1)，使模型中网格单元的节点坐标发生移动，形成初始规则波波浪形态；
η(x,t)＝a cos(kx-ωt) (1)
式中，η为距离平衡面的位移；a表示微幅波振幅；k表示波数；ω为波动圆频率；t为时间；x为站位。
前述的步骤c中，质点的动能赋值方法如下：
基于初始规则波波浪形态，得到的新的网格单元位置信息，结合式(2)和式(3)，计算新的网格单元的节点速度值；
z为水深，g为重力加速度，λ为波长，则有：
k＝2π/λ＝ω2/g (2)
u为水质点水平方向速度，v为水质点垂直方向速度，规则波的水质点速度为：
Figure BDA0003319630740000081
前述的步骤a中，仿真模型构建时，通过优化建模参数，使飞机的静水着水仿真计算精度满足工程需求。
仿真模型构建时，优化建模参数的方法为：将模型中处于规则波波浪液面附近及水上飞机着水位置的网格单元的尺寸细化。
以上所述实施仅表达了本发明的实施方式，其描述较为具体和详细，但并不能因此而理解为对本发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进，这些都属于本发明的保护范围。

Claims (8)

1.一种用于水上飞机着水数值仿真的规则波造波方法,其特征在于，包括下述步骤：
a.构建水上飞机静水着水数值仿真模型；
b.在所述的静水着水仿真模型上设计初始规则波波浪形态；
c.对已设计初始规则波波浪形态的静水着水仿真模型内的质点进行动能赋值，形成水上飞机规则波着水仿真模型。
2.根据权利要求1所述的水上飞机着水数值仿真的规则波造波方法，其特征在于，步骤a中，水上飞机静水着水数值仿真模型的构建原则为：考虑飞机在波浪水面着水后的滑行距离，建立“固-液-气”耦合水上飞机静水着水数值仿真模型。
3.根据权利要求1所述的水上飞机着水数值仿真的规则波造波方法，其特征在于，步骤a中，水上飞机静水着水数值仿真模型的构建方法为：
依据飞机总体气动外形，划分飞机结构有限元模型，不考虑水上飞机机体结构弹性特性，采用刚性体材料属性；
依据水上飞机重量特性，给飞机赋质量特性数据，定义“飞机-水体-空气”的耦合关系，定义流体边界条件；
定义水上飞机初始着水速度和输出项，采用ALE法或EULER法，完成水上飞机静水着水数值仿真计算模型的构建。
4.根据权利要求1所述的水上飞机着水数值仿真的规则波造波方法，其特征在于，所述的水上飞机静水着水数值仿真模型包括飞机结构模型、水流体模型和空气模型；所述的飞机结构模型的网格单元为拉格朗日单元；水流体模型和空气模型的网格单元均为ALE单元或EULER单元。
5.根据权利要求1所述的水上飞机着水数值仿真的规则波造波方法，其特征在于，步骤b中，初始规则波波浪形态的设计方法为：
根据水上飞机静水着水数值仿真分析模型生成的仿真计算输入文件和规则波自由液面方程(1)，使模型中网格单元的节点坐标发生移动，形成初始规则波波浪形态；
η(x,t)＝acos(kx-ωt) (1)
式中，η为距离平衡面的位移；a表示微幅波振幅；k表示波数；ω为波动圆频率；t为时间；x为站位。
6.根据权利要求1所述的水上飞机着水数值仿真的规则波造波方法，其特征在于，步骤c中，质点的动能赋值方法如下：
基于初始规则波波浪形态，得到的新的网格单元位置信息，结合式(2)和式(3)，计算新的网格单元的节点速度值；
z为水深，g为重力加速度，λ为波长，则有：
k＝2π/λ＝ω2/g (2)
u为水质点水平方向速度，v为水质点垂直方向速度，规则波的水质点速度为：
Figure FDA0003319630730000021
7.根据权利要求1所述的水上飞机着水数值仿真的规则波造波方法，其特征在于，步骤a中，仿真模型构建时，通过优化建模参数，使飞机的静水着水仿真计算精度满足工程需求。
8.根据权利要求7所述的水上飞机着水数值仿真的规则波造波方法，其特征在于，仿真模型构建时，优化建模参数的方法为：将模型中处于规则波波浪液面附近及水上飞机着水位置的网格单元的尺寸细化。
CN202111242287.5A 2021-10-25 2021-10-25 一种用于水上飞机着水数值仿真的规则波造波方法 Pending CN113935114A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111242287.5A CN113935114A (zh) 2021-10-25 2021-10-25 一种用于水上飞机着水数值仿真的规则波造波方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111242287.5A CN113935114A (zh) 2021-10-25 2021-10-25 一种用于水上飞机着水数值仿真的规则波造波方法

Publications (1)

Publication Number Publication Date
CN113935114A true CN113935114A (zh) 2022-01-14

Family

ID=79284185

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111242287.5A Pending CN113935114A (zh) 2021-10-25 2021-10-25 一种用于水上飞机着水数值仿真的规则波造波方法

Country Status (1)

Country Link
CN (1) CN113935114A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115525978A (zh) * 2022-10-31 2022-12-27 南京航空航天大学 用于飞行器波浪情况水动力分析的速度入口数值造波方法

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115525978A (zh) * 2022-10-31 2022-12-27 南京航空航天大学 用于飞行器波浪情况水动力分析的速度入口数值造波方法

Similar Documents

Publication Publication Date Title
CN104317985B (zh) 一种基于界带有限元和拉格朗日坐标的流体仿真方法
CN107391891A (zh) 一种基于模型融合方法的大展弦比机翼优化设计方法
CN106777482A (zh) 一种基于网格参数化的结构多学科设计优化方法
CN113935114A (zh) 一种用于水上飞机着水数值仿真的规则波造波方法
Thai et al. Cfd validation of small quadrotor performance using createtm-av helios
Park et al. Simulation of unsteady rotor flow field using unstructured adaptive sliding meshes
Wang et al. Transition-based constrained large-eddy simulation method with application to an ultrahigh-lift low-pressure turbine cascade flow
Lu et al. Flow simulation system based on high order space-time extension of flux reconstruction method
CN112199777B (zh) 一种适用于模化仿生前缘流场特征的方法
Sayed et al. Aeroelastic simulations based on high-fidelity CFD and CSD models
CN107369372A (zh) 水电站河床式厂房动力分析数值模型设计方法
CN111256946A (zh) 波浪模拟装置及方法
CN114077770A (zh) 一种用于飞行器水上迫降sph模拟的直接造波方法
Darracq et al. Computation of unsteady turbulent airfoil flows with an aeroelastic AUSM+ implicit solver
Stock A regularized inviscid vortex sheet method for three dimensional flows with density interfaces
Jung Hamiltonian Paths and Strands for Unified Grid Approach for Computing Aerodynamic Flows
CN114048696B (zh) 一种“扑动-划动-扭转”扑翼飞行器气动仿真分析方法
Lee et al. Optimal shape design of the S-shaped subsonic intake using NURBS
Miyaji et al. On Accuracy of Prediction of Flutter Boundaries on Unstructured Grids
CN104699913A (zh) 一种风力发电机组主机架强度计算方法
Xiao et al. Wing Flutter Simulations Using an Aeroelastic Solver Based on the Predictor—Corrector Scheme
Shanmugam et al. Numerical investigation on thrust production and unsteady mechanisms of three-dimensional oscillating wing
CN117829028B (zh) 浮式风机动力响应全耦合数值模拟方法、装置及相关设备
Tanabe et al. A new hybrid method of overlapping structured grids combined with unstructured fuselage grids for rotorcraft analysis
Wang et al. CFD-CSD method for rotor aeroelastic analysis with free wake model

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