CN112036095A - 海上救援直升机的飞行实时仿真方法及其仿真系统 - Google Patents
海上救援直升机的飞行实时仿真方法及其仿真系统 Download PDFInfo
- Publication number
- CN112036095A CN112036095A CN202010900472.8A CN202010900472A CN112036095A CN 112036095 A CN112036095 A CN 112036095A CN 202010900472 A CN202010900472 A CN 202010900472A CN 112036095 A CN112036095 A CN 112036095A
- Authority
- CN
- China
- Prior art keywords
- cfd
- rotor
- blade
- flight
- helicopter
- 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
- 238000004088 simulation Methods 0.000 title claims abstract description 95
- 238000000034 method Methods 0.000 title claims abstract description 60
- 238000004364 calculation method Methods 0.000 claims abstract description 67
- 230000015654 memory Effects 0.000 claims abstract description 42
- 238000010168 coupling process Methods 0.000 claims abstract description 41
- 230000008878 coupling Effects 0.000 claims abstract description 37
- 238000005859 coupling reaction Methods 0.000 claims abstract description 37
- 239000012530 fluid Substances 0.000 claims abstract description 25
- 230000008569 process Effects 0.000 claims description 26
- 230000033001 locomotion Effects 0.000 claims description 10
- 230000005540 biological transmission Effects 0.000 claims description 7
- 238000004891 communication Methods 0.000 claims description 7
- 230000004907 flux Effects 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 4
- 239000000284 extract Substances 0.000 claims description 3
- 238000000605 extraction Methods 0.000 claims description 3
- 238000012549 training Methods 0.000 abstract description 5
- 230000002457 bidirectional effect Effects 0.000 abstract description 4
- 238000005516 engineering process Methods 0.000 description 5
- 230000000007 visual effect Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000005553 drilling Methods 0.000 description 2
- 238000005094 computer simulation Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
Images
Classifications
-
- 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
- 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
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- Computer Hardware Design (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Analysis (AREA)
- Evolutionary Computation (AREA)
- Civil Engineering (AREA)
- Computational Mathematics (AREA)
- Structural Engineering (AREA)
- Architecture (AREA)
- Computer Graphics (AREA)
- Software Systems (AREA)
- Algebra (AREA)
- Computing Systems (AREA)
- Fluid Mechanics (AREA)
- Mathematical Physics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供了一种海上救援直升机的飞行实时仿真方法及其仿真系统,该方法包括:构建计算流体动力学CFD计算所需的计算机集群;选择仿真环境并基于仿真环境建立场景建筑物的三维模型,将直升机旋翼抽象为旋翼模型;使用飞行动力学计算单元计算旋翼受到的气动力并存储在共享内存中供CFD单元调用;CFD单元从共享内存中读取旋翼桨毂位置数据和气动力,计算得到耦合流场数据后保存在共享内存中供飞行动力学单元调用;重复计算直至降落。本发明采用双向耦合的方法模拟直升机海上飞行,提高仿真模拟的准确度,并可以通过CFD和视景分布式仿真子系统提高仿真模拟的实时性,以供飞行员进行飞行模拟训练。
Description
技术领域
本发明涉及计算机仿真技术领域,特别涉及一种海上救援直升机的飞行实时仿真方法及其仿真系统。
背景技术
随着我国直升机技术的日益发展,直升机已经成为了灾难事故中实施援救任务不可或缺的工具,海上救援任务同样不可忽视。直升机在海上着舰或者降落钻井平台等建筑物时,由于海平面上高速的风浪,以及旋翼与舰船甲板、机库和阻挡物之间形成复杂的回旋涡,造成直升机降落困难,也容易发生事故。因此,研制相应的驾驶模拟器非常必要。为了能够更精确模拟直升机与建筑物之间的耦合气流,有必要研究一种更加精确且实时模拟的方法,有必要研究一种新的模拟器系统。
现有技术中,仿真时,无法准确模拟直升机和海上建筑物的耦合气流,无法实时计算海上飞行时直升机旋翼、机身等部件的受力和位姿状态,且无法实时显示直升机飞行过程。由于现有技术中的仿真方法或系统与实际直升机飞行过程相差甚远,不利于驾驶员技术的提高。
发明内容
本发明针对上述现有技术中的缺陷,采用双向耦合的方法模拟直升机海上飞行,首先解算直升机气动力,传输气动力到计算流体动力学(CFD)单元,再将将直升机旋翼和海上建筑物的耦合气流利用CFD方法计算出来,再调用耦合气流数据并载入直升机气动力运算程序中,解算直升机气动力程序。本发明使用CFD单元和飞行动力学计算单元进行交互地、实时地仿真计算直升机海上的飞行过程,提高仿真模拟的准确度;所述多个CFD和视景模块计算机采用分布式的布置方式连接,以提高仿真的实时性。本发明关于海上救援直升机的飞行实时仿真方法和系统,实时地计算直升机降落过程中的受力情况及位姿状态,并可以通过CFD分布式仿真子系统和视景分布式仿真子系统实时仿真过程中的直升机的状态显示,以供飞行员进行飞行模拟训练,以尽快提高驾驶技术,解决了现有技术中的一个或多个技术问题。
本发明的具体方案如下:
本发明提供一种海上救援直升机的飞行实时仿真方法,该方法包括以下步骤:
步骤S1,构建计算流体动力学CFD计算所需的计算机集群,所述计算机集群包括多个计算机节点,将各个计算机节点利用TCP/IP通信协议建立网络连接,并从所述多个计算机节点选择一个计算机节点设置为总节点计算机用于控制CFD计算进程,并为CFD单元和飞行动力学计算单元设置共享内存的,用于在CFD单元和飞行动力学计算单元之间实时相互传输数据;
步骤S2,在总节点计算机控制CFD单元和飞行动力学计算单元初始化时,选择仿真环境并基于所述仿真环境建立场景建筑物的三维模型,将直升机旋翼抽象为旋翼模型,并利用展向和周向分布均匀的结构化网格划分旋翼网格,将所述建筑物的三维模型和所述旋翼网格作为CFD计算的模型,其中,所述仿真环境为直升机在舰上降落或在海上平台降落;
步骤S3,基于所述建筑物的三维模型和所述旋翼网格使用飞行动力学计算单元计算旋翼受到的气动力,输出旋翼桨毂位置数据和桨叶微元气动力数据到所述共享内存中供CFD单元调用,并发出待获取数据指令用于获取CFD计算得到的耦合流场数据;
步骤S4,所述CFD单元从所述共享内存中读取所述旋翼桨毂位置数据,然后更新流体计算域中旋翼的相对位置,并利用动量源法将所述桨叶微元气动力数据加载到所述旋翼网格的各网格单元中,所述CFD单元计算得到耦合流场数据后保存在所述共享内存中,所述CFD单元发出待获取数据指令等待所述飞行动力学计算单元传输下一时间的桨叶微元气动力数据;
步骤S5,所述飞行动力学计算单元从所述共享内存中读取所述耦合流场数据并返回步骤S3计算下一时间的直升机旋翼所受的气动力,并更新桨毂位置数据和桨叶微元气动力数据,判断直升机是否完成降落过程,如果否,则返回步骤S4进行计算,否则,结束飞行仿真计算并保存仿真数据。
更进一步地,所述旋翼模型为一层厚度网格的薄圆柱体,或由一层面网格组成,并采用结构网格划分桨叶圆盘为旋翼网格,然后在CFD单元中采用运动嵌套网格仿真旋翼的旋转运动。
更进一步地,所述CFD单元将计算得到的耦合流场数据通过所述总节点计算机发送至所述共享内存中,所述总节点计算机从所述共享内存中读取所述旋翼桨毂位置和桨叶微元气动力数据并发送到CFD单元中;所述耦合流场数据包括旋翼桨叶各微元中心点、机身质心、尾桨桨毂中心及垂尾和平尾气动中心处的流场速度值。
更进一步地,所述利用动量源法将所述桨叶微元气动力数据加载到所述旋翼网格的各网格单元中的操作为:
步骤S41,确定大地坐标系(Xd Yd Zd):选定水平面某一固定点作为坐标系原点,记为O,OXd指向直升机前进方向,OZd垂直向下,OYd方向遵守右手定则;确定直升机机体坐标系(X Y Z),固连于机体质心Ob,ObX轴位于机身对称面内指向机头位置,ObZ轴位于机身对称面内垂直向下,ObY轴方向遵守右手定则,在初始时刻直升机机体坐标系与大地坐标系重合;确定桨轴固定坐标系(Xh Yh Zh),原点位于旋翼桨毂中心,由机体坐标系旋转桨轴前倾角得到;确定桨轴旋转坐标系(Xs Ys Zs),方向为桨轴固定坐标系绕Zh轴转动ψ方位角度,原点位于桨轴中心;确定铰链坐标系(Xj Yj Zj),原点位于沿Ys距离e处,即旋翼挥舞铰和摆振铰中心处,方向与桨轴旋转坐标系重合;确定桨叶剖面坐标系(Xseg Yseg Zseg),旋翼桨叶绕Zj轴转过摆振角δ,然后绕Xj轴转动β挥舞角度,形成桨叶剖面坐标系,原点位于桨叶根部沿展向伸长距离r处;
步骤S42,旋翼桨叶微元剖面来流速度计算,记旋翼桨毂在桨轴固定坐标系下的线速度为(vhx,vhy,vhz),角速度为(whx,why,whz),桨叶绕固定轴Zh旋转角速度为Ω,则由于桨叶运动剖面来流速度(UPS,UTS,URS)可以表示为:
记CFD单元计算得到的桨叶剖面处耦合气流速度为(UPC,UTC,URC),则桨叶剖面总来流速度的分量形式(UP,UT,UR)为:
(UP,UT,UR)T=(UPS,UTS,URS)T-(UPC,UTC,URC)T
其中,UP正方向垂直于桨叶剖面向下;UT为桨叶切向来流,正方向指向翼型后缘;UR为桨叶展向来流,正方向指向桨叶展向外侧;
其中,ρ为气体密度,将桨叶各微元气动力数据存入共享内存中,供CFD单元提取处理;步骤S44,动量源项计算,令桨叶剖面来流角可以表示为:
所述CFD单元从所述共享内存中提取各微元气动力数据并将其转换到流体计算域,得到动量源项为:
步骤S45,CFD耦合流场计算,将动量源项带入CFD的通用计算标准控制方程中计算耦合流场,所述计算标准控制方程为:
更进一步地,在所述飞行动力学计算单元从所述共享内存中读取所述耦合流场数据并返回步骤S3计算下一时间的直升机旋翼所受的气动力的处理方式,与计算桨叶剖面来流速度同理,
所述直升机的机身、尾桨、平尾和垂尾处速度的计算可以统一表示为:
其中,表示机身质心、尾桨桨毂中心、平尾和垂尾气动力中心点的总速度;表示机身质心、尾桨桨毂中心、平尾和垂尾气动力中心点的运动速度;表示由CFD单元计算得到的位于机身质心、尾桨桨毂中心、平尾和垂尾气动力中心点的耦合流场速度;然后将作为速度值计算气动力。
本发明还提出了一种海上救援直升机的飞行实时仿真系统,该系统包括:
CFD分布式仿真子系统,其包括计算流体动力学CFD计算所需的计算机集群,所述计算机集群包括多个计算机节点,将各个计算机节点利用TCP/IP通信协议建立网络连接,并从所述多个计算机节点选择一个计算机节点设置为总节点计算机用于控制CFD计算进程;并为CFD单元和飞行动力学计算单元设置共享内存,用于在CFD单元和飞行动力学计算单元之间实时相互传输数据;在所述CFD分布式仿真子系统中运行上述任一项所述的飞行实时仿真方法;
视景分布式仿真子系统,其包括多个视景模块计算机,所述多个视景模块计算机采用分布式的布置方式连接,每台计算机控制特定仿真部分的显示画面,最后再将所述多个视景模块计算机的画面拼接在一起显示。
本发明的技术效果为:
本发明提出了一种海上救援直升机的飞行实时仿真方法及系统,该方法包括以下步骤:步骤S1,构建计算流体动力学CFD计算所需的计算机集群,所述计算机集群包括多个计算机节点,将各个计算机节点利用TCP/IP通信协议建立网络连接,并从所述多个计算机节点选择一个计算机节点设置为总节点计算机用于控制CFD计算进程,并为CFD单元和飞行动力学计算单元设置共享内存,用于在CFD单元和飞行动力学计算单元之间实时相互传输数据;步骤S2,在总节点计算机控制CFD单元和飞行动力学计算单元初始化时,选择仿真环境并基于所述仿真环境建立场景建筑物的三维模型,将直升机旋翼抽象为旋翼模型,并利用展向和周向分布均匀的结构化网格划分旋翼网格,将所述建筑物的三维模型和所述旋翼网格作为计算流体动力学CFD计算的模型,其中,所述仿真环境为直升机在舰上降落或在海上平台降落;步骤S3,基于所述建筑物的三维模型和所述旋翼网格使用飞行动力学计算单元计算旋翼受到的气动力,输出旋翼桨毂位置数据和桨叶微元气动力数据到所述共享内存中供CFD单元调用,并发出待获取数据指令用于获取CFD计算得到的耦合流场数据;步骤S4,所述CFD单元从所述共享内存中读取所述旋翼桨毂位置数据,然后更新流体计算域中旋翼的相对位置,并利用动量源法将所述桨叶微元气动力数据加载到所述旋翼网格的各网格单元中,所述CFD单元计算得到耦合流场数据后保存在所述共享内存中,所述CFD单元发出待获取数据指令等待所述飞行动力学计算单元传输下一时间的桨叶微元气动力数据;步骤S5,所述飞行动力学计算单元从所述共享内存中读取所述耦合流场数据并返回步骤S3计算下一时间的直升机旋翼所受的气动力,并更新桨毂位置数据和桨叶微元气动力数据,判断直升机是否完成降落过程,如果否,则返回步骤S4进行计算,否则,结束飞行仿真计算并保存仿真数据。本发明采用双向耦合的方法模拟直升机海上飞行,首先解算直升机气动力,传输气动力到计算流体动力学(CFD)程序,再将将直升机旋翼和海上建筑物的耦合气流利用CFD方法计算出来,再调用耦合气流数据并载入直升机气动力运算程序中,解算直升机气动力。本发明使用CFD单元和飞行动力学计算单元进行交互地、实时地仿真计算直升机海上的飞行过程,提高仿真模拟的准确度;所述多个CFD和视景模块计算机采用分布式的布置方式连接,以提高仿真的实时性。本发明关于海上救援直升机的飞行实时仿真方法和系统,实时地计算直升机降落过程中的受力情况及位姿状态,并可以通过CFD和视景分布式仿真子系统实时仿真过程中的直升机的状态显示,以供飞行员进行飞行模拟训练,以尽快提高驾驶技术。
附图说明
通过阅读参照以下附图所作的对非限制性实施例所作的详细描述,本申请的其它特征、目的和优点将会变得更明显。
图1为一种海上救援直升机的飞行实时仿真方法的流程图;
图2为本发明实施例中的为旋翼流体计算域网格及微元气动力点示意图;
图3为本发明实施例中的为桨轴坐标系示意图;
图4本发明海上救援直升机的飞行实时仿真系统的结构图。
具体实施方式
下面结合附图和实施例对本申请作进一步的详细说明。可以理解的是,此处所描述的具体实施例仅仅用于解释相关发明,而非对该发明的限定。另外还需要说明的是,为了便于描述,附图中仅示出了与有关发明相关的部分。
需要说明的是,在不冲突的情况下,本申请中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本申请。
本发明旨在提供一种海上救援直升机的飞行实时仿真方法,该方法包括以下步骤:
步骤S1,构建计算流体动力学CFD计算所需的计算机集群,计算机集群包括多个计算机节点,将各个计算机节点利用TCP/IP通信协议建立网络连接,并从多个计算机节点选择一个计算机节点设置为总节点计算机用于控制CFD计算进程,并为CFD单元和飞行动力学计算单元设置共享内存的,用于在CFD单元和飞行动力学计算单元之间实时相互传输数据;本发明中通过总节点计算机,也称为主节点计算机。通过分布式仿真的控制,并通过共享内存的方式进行数据的传输,提高了程序之间数据交换的速度,进而提高了仿真的性能,这是本发明的一个重要发明点。
步骤S2,在总节点计算机控制CFD单元和飞行动力学计算单元初始化时,选择仿真环境并基于仿真环境建立场景建筑物的三维模型,将直升机旋翼抽象为旋翼模型,并利用展向和周向分布均匀的结构化网格划分旋翼网格,将建筑物的三维模型和旋翼网格作为CFD计算的模型,其中,仿真环境为直升机在舰上降落或在海上平台降落。
在一个实施例中,指定模拟场景(着舰或海上平台降落),建立场景建筑物(舰船或钻井平台)三维模型用于CFD方法计算。将直升机旋翼简化为一微小高度的圆盘,在CFD中利用展向和周向分布均匀的结构化网格划分旋翼网格,将建筑物和简化后的旋翼作为CFD计算的模型;将旋翼各网格中心点的位置确定为旋翼各桨叶微元气动力点的位置,用于计算中加载气动力。如图2所示,将旋翼周向均匀划分为72等份,展向均匀划分为36等份,网格厚度为一层,高度为0.01m的六面体网格,更进一步地,旋翼模型为一层厚度网格的薄圆柱体,或由一层面网格组成,并采用结构网格划分桨叶圆盘为旋翼网格,然后在CFD计算中采用运动嵌套网格仿真旋翼的旋转运动。网格体心坐标作为各微元气动力点位置,用于添加动量源项。
步骤S3,使用飞行动力学计算单元计算旋翼受到的气动力,输出旋翼桨毂位置数据和桨叶微元气动力数据到共享内存中供CFD单元调用,并发出待获取数据指令用于获取CFD计算得到的耦合流场数据;其中,首次气动力计算中不计耦合流场的影响。
步骤S4,CFD单元从共享内存中读取旋翼桨毂位置数据,然后更新流体计算域中旋翼的相对位置,并利用动量源法将桨叶微元气动力数据加载到旋翼网格的各网格单元中,CFD单元计算得到耦合流场数据后保存在共享内存中,CFD单元发出待获取数据指令等待飞行动力学计算单元传输下一时间的桨叶微元气动力数据。
步骤S5,飞行动力学计算单元从共享内存中读取耦合流场数据并返回步骤S3计算下一时间的直升机旋翼所受的气动力,并更新桨毂位置数据和桨叶微元气动力数据,判断直升机是否完成降落过程,如果否,则返回步骤S4进行计算,否则,结束飞行仿真计算并保存仿真数据。在飞机降落仿真过程中,还采用CFD和视景分布式仿真子系统,其包括多个计算机节点,多个CFD计算机节点采用分布式的布置方式连接,以提高仿真计算的速度;每台视景模块计算机控制特定仿真部分的显示画面,最后再将多个视景模块计算机的画面拼接在一起显示,以实时地显示直升机降落的仿真过程,这也是本发明的另一个重要发明点。
本发明采用双向耦合的方法模拟直升机海上飞行,首先解算直升机气动力,传输气动力到计算流体动力学(CFD)单元,再将将直升机旋翼和海上建筑物的耦合气流利用CFD方法计算出来,再调用耦合气流数据并载入直升机气动力运算程序中,解算直升机气动力程序。本发明使用CFD单元和飞行动力学计算单元进行交互地、实时地仿真计算直升机海上的飞行过程,提高仿真模拟的准确度,这是本申请的一个重要发明点。
在一个实施例,CFD单元将计算得到的耦合流场数据通过总节点计算机发送至共享内存中,总节点计算机从共享内存中读取旋翼桨毂位置和桨叶微元气动力数据并发送到CFD单元中;耦合流场数据包括旋翼桨叶各微元中心点、机身质心、尾桨桨毂中心及垂尾和平尾气动中心处的流场速度值。
在一个实施例中,利用动量源法将桨叶微元气动力数据加载到旋翼网格的各网格单元中的操作为:
步骤S41,确定大地坐标系(Xd Yd Zd):选定水平面某一固定点作为坐标系原点,记为O,OXd指向直升机前进方向,OZd垂直向下,OYd方向遵守右手定则;确定直升机机体坐标系(X Y Z),固连于机体质心Ob,ObX轴位于机身对称面内指向机头位置,ObZ轴位于机身对称面内垂直向下,ObY轴方向遵守右手定则,在初始时刻直升机机体坐标系与大地坐标系重合;确定桨轴固定坐标系(Xh Yh Zh),原点位于旋翼桨毂中心,由机体坐标系旋转桨轴前倾角得到;确定桨轴旋转坐标系(Xs Ys Zs),方向为桨轴固定坐标系绕Zh轴转动ψ方位角度,原点位于桨轴中心,如图3所示;确定铰链坐标系(Xj Yj Zj),原点位于沿Ys距离e处,即旋翼挥舞铰和摆振铰中心处,方向与桨轴旋转坐标系重合;确定桨叶剖面坐标系(Xseg Yseg Zseg),旋翼桨叶绕Zj轴转过摆振角δ,然后绕Xj轴转动β挥舞角度,形成桨叶剖面坐标系,原点位于桨叶根部沿展向伸长距离r处;
步骤S42,旋翼桨叶微元剖面来流速度计算,记旋翼桨毂在桨轴固定坐标系下的线速度为(vhx,vhy,vhz),角速度为(whx,why,whz),桨叶绕固定轴Zh旋转角速度为Ω,则由于桨叶运动剖面来流速度(UPS,UTS,URS)可以表示为:
记CFD单元计算得到的桨叶剖面处耦合气流速度为(UPC,UTC,URC),则桨叶剖面总来流速度的分量形式(UP,UT,UR)为:
(UP,UT,UR)T=(UPS,UTS,URS)T-(UPC,UTC,URC)T
其中,UP正方向垂直于桨叶剖面向下;UT为桨叶切向来流,正方向指向翼型后缘;UR为桨叶展向来流,正方向指向桨叶展向外侧;
其中,ρ为气体密度,将桨叶各微元气动力数据存入共享内存中,供CFD单元提取处理;步骤S44,动量源项计算,令桨叶剖面来流角可以表示为:
CFD单元从共享内存中提取各微元气动力数据并将其转换到流体计算域,得到动量源项为:
步骤S45,CFD耦合流场计算,将动量源项带入CFD的通用计算标准控制方程中计算耦合流场,计算标准控制方程为:
如图3所示,一个具体的实施例为:A点位于桨叶方位角ψ=0,距桨叶根部r=0.4m处,且旋翼桨毂在桨轴固定坐标系下的线速度为(0.5,0.5,0.5),角速度为(0.1,0.2,0.3),桨叶绕绕根部旋转角速度为36.55m/s,挥舞摆振铰偏移距离e=0.1,挥舞角β=0.1,摆振角δ=0.01,挥舞角速度摆振角速度则由于桨叶运动剖面来流速度(UPS,UTS,URS)可以表示为:
记CFD计算得到的桨叶剖面处耦合气流速度为(10.0,2.6,0.5),则桨叶剖面总来流速度(UP,UT,UR)为:
(UP,UT,UR)T=(0.6081 18.6335 0.4540)T-(10.0 2.6 0.5)T=(-9.3919 16.0335-0.0460)T
则
取桨叶微元展向宽度为Δy=0.025,升阻力系数分别为cl=0.1003,cd=0.0003,剖面翼型弦长为cy=0.1,密度ρ=1.293,则桨叶微元所受升阻力为:
桨叶剖面来流角可以表示为:
则CFD计算域下动量源项为:
由此,可解出动量源项,然后将其加入CFD控制方程中计算耦合流场。
在一个实施例中,在飞行动力学计算单元从共享内存中读取耦合流场数据并返回步骤S3计算下一时间的直升机旋翼所受的气动力的处理方式,与计算桨叶剖面来流速度同理,直升机的机身、尾桨、平尾和垂尾处速度的计算可以统一表示为:
其中,表示机身质心、尾桨桨毂中心、平尾和垂尾气动力中心点的总速度;表示机身质心、尾桨桨毂中心、平尾和垂尾气动力中心点的运动速度;表示由CFD计算得到的位于机身质心、尾桨桨毂中心、平尾和垂尾气动力中心点的耦合流场速度;然后将作为速度值计算气动力。
本发明通过上述步骤,实时地计算直升机降落过程中的受力情况及位姿状态,并可以通过CFD和视景分布式仿真子系统实时仿真过程中的直升机的状态显示,以供飞行员进行飞行模拟训练,以尽快提高驾驶技术,这是本发明的重要发明点。
如图4所示,一种海上救援直升机的飞行实时仿真系统,该系统包括:CFD分布式仿真子系统,其包括计算流体动力学CFD计算所需的计算机集群,计算机集群包括多个计算机节点,将各个计算机节点利用TCP/IP通信协议建立网络连接,并从多个计算机节点选择一个计算机节点设置为总节点计算机用于控制CFD计算进程,并为CFD单元和飞行动力学计算单元设置共享内存,用于在CFD单元和飞行动力学计算单元之间实时相互传输数据;在CFD分布式仿真子系统中运行上述任一项方法;视景分布式仿真子系统,其包括多个视景模块计算机,多个视景模块计算机采用分布式的布置方式连接,每台计算机控制特定仿真部分的显示画面,最后再将多个视景模块计算机的画面拼接在一起显示。
本发明的海上救援直升机的飞行实时仿真系统,实时地计算直升机降落过程中的受力情况及位姿状态,并可以通过CFD和视景分布式仿真子系统实时仿真过程中的直升机的状态显示,以供飞行员进行飞行模拟训练,以尽快提高驾驶技术,这是本发明的重要发明点。
最后所应说明的是:以上实施例仅以说明而非限制本发明的技术方案,尽管参照上述实施例对本发明进行了详细说明,本领域的普通技术人员应当理解:依然可以对本发明进行修改或者等同替换,而不脱离本发明的精神和范围的任何修改或局部替换,其均应涵盖在本发明的权利要求范围当中。
Claims (6)
1.一种海上救援直升机的飞行实时仿真方法,其特征在于,其包括以下步骤:
步骤S1,构建计算流体动力学CFD计算所需的计算机集群,所述计算机集群包括多个计算机节点,将各个计算机节点利用TCP/IP通信协议建立网络连接,并从所述多个计算机节点选择一个计算机节点设置为总节点计算机用于控制CFD仿真计算进程,并为CFD单元和飞行动力学计算单元设置共享内存,用于在CFD单元和飞行动力学计算单元之间实时相互传输数据;
步骤S2,在总节点计算机控制CFD单元和飞行动力学计算单元初始化时,选择仿真环境并基于所述仿真环境建立场景建筑物的三维模型,将直升机旋翼抽象为旋翼模型,并利用展向和周向分布均匀的结构化网格划分旋翼网格,将所述建筑物的三维模型和所述旋翼网格作为CFD计算的模型,其中,所述仿真环境为直升机在舰上降落或在海上平台降落;
步骤S3,基于所述建筑物的三维模型和所述旋翼网格使用飞行动力学计算单元计算旋翼受到的气动力,输出旋翼桨毂位置数据和桨叶微元气动力数据到所述共享内存中供CFD单元调用,并发出待获取数据指令用于获取CFD计算得到的耦合流场数据;
步骤S4,所述CFD单元从所述共享内存中读取所述旋翼桨毂位置数据,更新流体计算域中旋翼的相对位置,并将所述桨叶微元气动力数据加载到所述旋翼网格的各网格单元中,所述CFD单元计算得到耦合流场数据后保存在所述共享内存中,所述CFD单元发出待获取数据指令等待所述飞行动力学计算单元传输下一时间的桨叶微元气动力数据;
步骤S5,所述飞行动力学计算单元从所述共享内存中读取所述耦合流场数据并返回步骤S3计算下一时间的直升机旋翼所受的气动力,并更新桨毂位置数据和桨叶微元气动力数据,判断直升机是否完成降落过程,如果否,则返回步骤S4进行计算;如果是,则结束飞行仿真计算并保存仿真数据。
2.根据权利要求1所述的海上救援直升机的飞行实时仿真方法,其特征在于,所述旋翼模型为一层厚度网格的薄圆柱体,或由一层面网格组成,并采用结构网格划分桨叶圆盘为旋翼网格,然后在CFD程序中采用运动嵌套网格仿真旋翼的旋转运动。
3.根据权利要求2所述的海上救援直升机的飞行实时仿真方法,其特征在于,所述CFD单元将计算得到的耦合流场数据通过所述总节点计算机发送至所述共享内存中,所述总节点计算机从所述共享内存中读取所述旋翼桨毂位置和桨叶微元气动力数据并发送到CFD单元中;所述耦合流场数据包括旋翼桨叶各微元中心点、机身质心、尾桨桨毂中心及垂尾和平尾气动中心处的流场速度值。
4.根据权利要求3所述的海上救援直升机的飞行实时仿真方法,其特征在于,在所述步骤S4中,将所述桨叶微元气动力数据加载到所述旋翼网格的各网格单元中的具体步骤为:
步骤S41,确定大地坐标系(Xd Yd Zd):选定水平面某一固定点作为坐标系原点,记为O,OXd指向直升机前进方向,OZd垂直向下,OYd方向遵守右手定则;确定直升机机体坐标系(X YZ),固连于机体质心Ob,ObX轴位于机身对称面内指向机头位置,ObZ轴位于机身对称面内垂直向下,ObY轴方向遵守右手定则,在初始时刻直升机机体坐标系与大地坐标系重合;确定桨轴固定坐标系(Xh Yh Zh),原点位于旋翼桨毂中心,由机体坐标系旋转桨轴前倾角得到;确定桨轴旋转坐标系(Xs Ys Zs),方向为桨轴固定坐标系绕Zh轴转动ψ方位角度,原点位于桨轴中心;确定铰链坐标系(Xj Yj Zj),原点位于沿Ys距离e处,即旋翼挥舞铰和摆振铰中心处,方向与桨轴旋转坐标系重合;确定桨叶剖面坐标系(Xseg Yseg Zseg),旋翼桨叶绕Zj轴转过摆振角δ,然后绕Xj轴转动β挥舞角度,形成桨叶剖面坐标系,原点位于桨叶根部沿展向伸长距离r处;
步骤S42,旋翼桨叶微元剖面来流速度计算,记旋翼桨毂在桨轴固定坐标系下的线速度为(vhx,vhy,vhz),角速度为(whx,why,whz),桨叶绕固定轴Zh旋转角速度为Ω,则由于桨叶运动剖面来流速度(UPS,UTS,URS)可以表示为:
记CFD计算得到的桨叶剖面处耦合气流速度为(UPC,UTC,URC),则桨叶剖面总来流速度的分量形式(UP,UT,UR)为:
(UP,UT,UR)T=(UPS,UTS,URS)T-(UPC,UTC,URC)T
其中,UP正方向垂直于桨叶剖面向下;UT为桨叶切向来流,正方向指向翼型后缘;UR为桨叶展向来流,正方向指向桨叶展向外侧;
其中,ρ为气体密度,将桨叶各微元气动力数据存入共享内存中,供CFD单元提取处理;
步骤S44,动量源项计算,令桨叶剖面来流角可以表示为:
所述CFD单元从所述共享内存中提取各微元气动力数据并将其转换到流体计算域,得到动量源项为:
步骤S45,耦合流场计算,将动量源项带入CFD的通用计算标准控制方程中计算耦合流场,所述计算标准控制方程为:
6.一种根据权利要求1-5之一的海上救援直升机的飞行实时仿真方法的海上救援直升机的飞行实时仿真系统,其特征在于,该系统包括:
CFD分布式仿真子系统,其包括计算流体动力学CFD计算所需的计算机集群,所述计算机集群包括多个计算机节点,将各个计算机节点利用TCP/IP通信协议建立网络连接,并从所述多个计算机节点选择一个计算机节点设置为总节点计算机用于控制CFD计算进程;并为CFD单元和飞行动力学计算单元设置共享内存,用于在CFD单元和飞行动力学计算单元之间实时相互传输数据;在所述CFD分布式仿真子系统中运行权利要求1-5任一项所述的飞行实时仿真方法;
视景分布式仿真子系统,其包括多个视景模块计算机,所述多个视景模块计算机采用分布式的布置方式连接,每台计算机控制特定仿真部分的显示画面,最后再将所述多个视景模块计算机的画面拼接在一起显示。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010900472.8A CN112036095B (zh) | 2020-08-31 | 2020-08-31 | 海上救援直升机的飞行实时仿真方法及其仿真系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010900472.8A CN112036095B (zh) | 2020-08-31 | 2020-08-31 | 海上救援直升机的飞行实时仿真方法及其仿真系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112036095A true CN112036095A (zh) | 2020-12-04 |
CN112036095B CN112036095B (zh) | 2023-04-25 |
Family
ID=73587041
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010900472.8A Active CN112036095B (zh) | 2020-08-31 | 2020-08-31 | 海上救援直升机的飞行实时仿真方法及其仿真系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112036095B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113353284A (zh) * | 2021-06-29 | 2021-09-07 | 贵州理工学院 | 一种纵列双旋翼无人植保系统、构建方法、介质、计算机 |
CN115804926A (zh) * | 2022-12-27 | 2023-03-17 | 中国人民解放军陆军装甲兵学院 | 一种无人机投射灭火弹瞄准系统及其方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8452579B1 (en) * | 2010-05-26 | 2013-05-28 | Textron Innovations Inc. | System and method of simulating a drive-train system |
CN103809464A (zh) * | 2014-01-27 | 2014-05-21 | 中国人民解放军空军航空大学军事仿真技术研究所 | 直升机舰面效应影响的仿真方法 |
CN103810332A (zh) * | 2014-01-27 | 2014-05-21 | 中国人民解放军空军航空大学军事仿真技术研究所 | 直升机舰面尾流影响仿真方法 |
CN106650046A (zh) * | 2016-12-02 | 2017-05-10 | 中国船舶工业系统工程研究院 | 一种舰船空气流场的非定常特性获取方法 |
CN108090302A (zh) * | 2018-01-04 | 2018-05-29 | 北京中航智科技有限公司 | 一种直升机飞行力学模拟方法及系统 |
CN108762295A (zh) * | 2018-02-09 | 2018-11-06 | 华南理工大学 | 基于软件总线的一体化无人机控制系统 |
CN108920811A (zh) * | 2018-06-28 | 2018-11-30 | 南京航空航天大学 | 一种用于直升机飞行仿真的模拟方法及系统 |
-
2020
- 2020-08-31 CN CN202010900472.8A patent/CN112036095B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8452579B1 (en) * | 2010-05-26 | 2013-05-28 | Textron Innovations Inc. | System and method of simulating a drive-train system |
CN103809464A (zh) * | 2014-01-27 | 2014-05-21 | 中国人民解放军空军航空大学军事仿真技术研究所 | 直升机舰面效应影响的仿真方法 |
CN103810332A (zh) * | 2014-01-27 | 2014-05-21 | 中国人民解放军空军航空大学军事仿真技术研究所 | 直升机舰面尾流影响仿真方法 |
CN106650046A (zh) * | 2016-12-02 | 2017-05-10 | 中国船舶工业系统工程研究院 | 一种舰船空气流场的非定常特性获取方法 |
CN108090302A (zh) * | 2018-01-04 | 2018-05-29 | 北京中航智科技有限公司 | 一种直升机飞行力学模拟方法及系统 |
CN108762295A (zh) * | 2018-02-09 | 2018-11-06 | 华南理工大学 | 基于软件总线的一体化无人机控制系统 |
CN108920811A (zh) * | 2018-06-28 | 2018-11-30 | 南京航空航天大学 | 一种用于直升机飞行仿真的模拟方法及系统 |
Non-Patent Citations (6)
Title |
---|
吴文韬: "直升机舰面起降飞行特性及起降风限图研究", 《中国优秀博硕士学位论文全文数据库(硕士) 工程科技Ⅱ辑》 * |
李国如: "无人直升机舰面起降仿真技术研究", 《中国优秀博硕士学位论文全文数据库(硕士) 工程科技Ⅱ辑》 * |
王超: "直升机着舰过程的流场数值模拟", 《中国优秀博硕士学位论文全文数据库(硕士) 基础科学辑》 * |
胡彦宇: "复杂系统仿真数据交换技术研究", 《中国优秀博硕士学位论文全文数据库(硕士) 基础科学辑》 * |
贺少华等: "舰载直升机理论风限图飞行模拟获取方法", 《海军航空工程学院学报》 * |
陈华健等: "直升机动态着舰流场下的旋翼气动载荷分析", 《南京航空航天大学学报》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113353284A (zh) * | 2021-06-29 | 2021-09-07 | 贵州理工学院 | 一种纵列双旋翼无人植保系统、构建方法、介质、计算机 |
CN115804926A (zh) * | 2022-12-27 | 2023-03-17 | 中国人民解放军陆军装甲兵学院 | 一种无人机投射灭火弹瞄准系统及其方法 |
CN115804926B (zh) * | 2022-12-27 | 2023-08-22 | 中国人民解放军陆军装甲兵学院 | 一种无人机投射灭火弹瞄准系统及其方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112036095B (zh) | 2023-04-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Ventura Diaz et al. | High-fidelity computational aerodynamics of multi-rotor unmanned aerial vehicles | |
Czyba et al. | Construction prototyping, flight dynamics modeling, and aerodynamic analysis of hybrid VTOL unmanned aircraft | |
CN106155076B (zh) | 一种多旋翼无人飞行器的稳定飞行控制方法 | |
CN106114854B (zh) | 一种无人驾驶航空器 | |
CN106114853A (zh) | 一种无人驾驶航空器 | |
CN107102565A (zh) | 无人机集群软件在环仿真系统 | |
CN112036095B (zh) | 海上救援直升机的飞行实时仿真方法及其仿真系统 | |
WO2022048543A1 (zh) | 一种飞行控制方法、无人机及存储介质 | |
CN108920811A (zh) | 一种用于直升机飞行仿真的模拟方法及系统 | |
CN103810332B (zh) | 直升机舰面尾流影响仿真方法 | |
Kita et al. | Transition between level flight and hovering of a tail-sitter vertical takeoff and landing aerial robot | |
Lee et al. | Development of Unified High‐Fidelity Flight Dynamic Modeling Technique for Unmanned Compound Aircraft | |
Bolander et al. | Static Trim of a Bio-Inspired Rotating Empennage for a Fighter Aircraft | |
CN211685678U (zh) | 一种多旋翼无人机实时尾迹的仿真分析系统 | |
Reed et al. | Coaxial Rotor Wake and Prop Induction Impact on a Horizontal Tail Using HPCMP CREATE™-AV Helios | |
CN116560249A (zh) | 一种机动飞行的高自由度简化建模与轨迹仿真方法 | |
CN110414110A (zh) | 一种用于飞行失速状态下的飞机受力仿真方法 | |
CN103809464B (zh) | 直升机舰面效应影响的仿真方法 | |
CN207881710U (zh) | 一种多旋翼飞行器的飞行姿态分析平台 | |
Park et al. | Hybrid analysis for quadrotor type uav and modified blade element momentum theory considering gust and flight condition | |
Mancinelli et al. | Dual-axis tilting rotor quad-plane design, simulation, flight and performance comparison with a conventional quad-plane design | |
Alam et al. | Robust hover mode control of a tiltrotor using nonlinear control technique | |
Sasongko et al. | Modeling and simulation of rotor dynamics of a tilt-rotor aircraft | |
Uhlig et al. | Stability characteristics of micro air vehicles from experimental measurements | |
US20220335839A1 (en) | Airspace information provision system, airspace information provision method, and airspace information provision program |
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 |