CN114476123B - 一种用于软管式空中加油对接模拟试验装置的设计方法 - Google Patents
一种用于软管式空中加油对接模拟试验装置的设计方法 Download PDFInfo
- Publication number
- CN114476123B CN114476123B CN202210277521.6A CN202210277521A CN114476123B CN 114476123 B CN114476123 B CN 114476123B CN 202210277521 A CN202210277521 A CN 202210277521A CN 114476123 B CN114476123 B CN 114476123B
- Authority
- CN
- China
- Prior art keywords
- hose
- refueling
- aerodynamic
- coefficient
- section
- 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 44
- 238000012360 testing method Methods 0.000 title claims abstract description 41
- 238000000034 method Methods 0.000 title claims abstract description 40
- 238000013461 design Methods 0.000 title claims abstract description 27
- 210000001503 joint Anatomy 0.000 title abstract description 11
- 238000003032 molecular docking Methods 0.000 claims abstract description 36
- 230000008569 process Effects 0.000 claims abstract description 10
- 230000004044 response Effects 0.000 claims abstract description 7
- 239000000446 fuel Substances 0.000 claims description 20
- 239000000463 material Substances 0.000 claims description 14
- 239000011159 matrix material Substances 0.000 claims description 11
- 238000009826 distribution Methods 0.000 claims description 9
- 239000012530 fluid Substances 0.000 claims description 9
- 230000001133 acceleration Effects 0.000 claims description 6
- 230000005484 gravity Effects 0.000 claims description 6
- 239000000945 filler Substances 0.000 claims description 4
- 230000000630 rising effect Effects 0.000 claims description 4
- 238000005452 bending Methods 0.000 claims description 2
- 230000017105 transposition Effects 0.000 claims description 2
- 239000000126 substance Substances 0.000 claims 2
- 230000006872 improvement Effects 0.000 description 8
- 238000004364 calculation method Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 5
- 230000007246 mechanism Effects 0.000 description 5
- 230000008901 benefit Effects 0.000 description 4
- 238000005457 optimization Methods 0.000 description 4
- 230000007547 defect Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000001808 coupling effect Effects 0.000 description 2
- 238000013016 damping Methods 0.000 description 2
- 238000009795 derivation Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- NCGICGYLBXGBGN-UHFFFAOYSA-N 3-morpholin-4-yl-1-oxa-3-azonia-2-azanidacyclopent-3-en-5-imine;hydrochloride Chemical compound Cl.[N-]1OC(=N)C=[N+]1N1CCOCC1 NCGICGYLBXGBGN-UHFFFAOYSA-N 0.000 description 1
- RZVHIXYEVGDQDX-UHFFFAOYSA-N 9,10-anthraquinone Chemical compound C1=CC=C2C(=O)C3=CC=CC=C3C(=O)C2=C1 RZVHIXYEVGDQDX-UHFFFAOYSA-N 0.000 description 1
- 206010034719 Personality change Diseases 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000007664 blowing Methods 0.000 description 1
- 230000001066 destructive effect Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000005764 inhibitory process Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 238000012806 monitoring device Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
Images
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64F—GROUND OR AIRCRAFT-CARRIER-DECK INSTALLATIONS SPECIALLY ADAPTED FOR USE IN CONNECTION WITH AIRCRAFT; DESIGNING, MANUFACTURING, ASSEMBLING, CLEANING, MAINTAINING OR REPAIRING AIRCRAFT, NOT OTHERWISE PROVIDED FOR; HANDLING, TRANSPORTING, TESTING OR INSPECTING AIRCRAFT COMPONENTS, NOT OTHERWISE PROVIDED FOR
- B64F5/00—Designing, manufacturing, assembling, cleaning, maintaining or repairing aircraft, not otherwise provided for; Handling, transporting, testing or inspecting aircraft components, not otherwise provided for
- B64F5/60—Testing or inspecting aircraft components or systems
-
- 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)
- Manufacturing & Machinery (AREA)
- Transportation (AREA)
- Aviation & Aerospace Engineering (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种用于软管式空中加油对接模拟试验装置的设计方法,所述方法包括:步骤1)建立软管式空中加油对接模拟装置的相似准则,所述空中加油对接模拟装置包括加油机、加油软管、伞锥和受油机;步骤2)基于相似准则,结合柔性多体动力学,建立空中加油对接模拟装置的动力学模型;步骤3)根据模拟试验要求,基于动力学模型,根据相似准则实现软管式空中加油对接模拟装置的设计,利用多体动力学求解器计算受油机对接伞锥后软管的动力学响应,进而实现加油对接过程中加油软管和伞锥的动态现象及动态特性的模拟。
Description
技术领域
本发明属于航空技术领域,尤其涉及一种用于软管式空中加油对接模拟试验装置的设计方法。
背景技术
空中加油,是当今航空领域,特别是军用飞行器在不着陆条件下,实现燃油快速补给、提升战机滞空时间、及时投入战场作战的重要技术手段;它对于增加战机航程与挂载能力,扩大作战半径具有重要意义,是飞机作战效能的倍增器。
从世界范围看,根据加油管路的不同设计,空中加油可以分为软管加油和硬管加油两种方式。硬管加油系统主要由伸缩管、压力加油机构、控制舵和监控装置等组成。该系统优点是耐压性好,可采用增压设备提高输油速度;缺点是系统复杂,对飞行控制、加油系统制造技术要求比较高,成本也较高。我国的空中加油技术,普遍采用“软管+锥套”的软管式加油系统。
软管加油系统的优点是结构简单、成本较低,一架加油机可同时安装多套系统,实现同时为多架飞机加油。软管加油系统因其自身具有一定柔性,故对加油机和受油机的相对位置、飞行稳定性要求相对较低;但其缺点是对于气流扰动比较敏感,软管柔性材料形变、燃油输出压力、大气紊流/阵风扰动、加油机尾流(如图1所示)、机体振动、受油机头波等内外干扰因素的影响。同时,在对接过程中对于受油机的操控(如:对接冲击速度)要求较高,操作不当轻则导致加油操作的失败,重则致使加受油机受损、危及飞行安全。这其中,容易出现一种极具破坏力的现象,加油软管由于受油机与锥套对接冲击力作用而过度松弛,从而诱发剧烈甩动,即“鞭甩”现象。
在软管式空中加油过程中出现的一系列动态、非定常、非稳态现象,是空气动力学、多体动力学、机械/材料学、自动控制等多学科耦合作用的结果。这些现象,制约了空中加油任务的成功率,并对飞行安全造成严重影响。研究并揭示这些现象形成的内在机理与成因,并建立地面模拟试验技术,提出有效的抑制/解决措施,优化加受油系统设计,对于提高空中加油的安全性、可靠性具有十分重要的意义。
发明内容
本发明的目的在于克服现有技术缺陷,提出了一种用于软管式空中加油对接模拟试验装置的设计方法。
为了实现上述目的,本发明提出了一种用于软管式空中加油对接模拟试验装置的设计方法,所述方法包括:
步骤1)建立软管式空中加油对接模拟试验装置的相似准则,所述空中加油对接模拟装置包括加油机、加油软管、伞锥和受油机;
步骤2)基于相似准则,结合柔性多体动力学,建立空中加油对接模拟试验装置的动力学模型;
步骤3)根据模拟试验要求,基于动力学模型,根据相似准则实现软管式空中加油对接模拟装置的设计,利用多体动力学求解器计算受油机对接伞锥后软管的动力学响应,进而实现加油对接过程中加油软管和伞锥的动态现象及动态特性的模拟。
作为上述方法的一种改进,所述步骤1)具体包括:
建立流体动力学相似准则和结构动力学相似准则,其中,
所述流体动力学相似准则包括:几何外形相似、速度场相似、动压相似、流场压缩性相似和重力影响相似;
所述结构动力学相似准则包括:质量分布相似、固有频率相似、振型相似和刚度分布相似。
作为上述方法的一种改进,所述流体动力学相似准则满足下式:
λV=1
其中,λV为运动速度系数;
λα=1
其中,λα为角度系数
其中,λg为重力加速度系数,λD为结构尺寸缩比系数;
其中,Dm为加油软管-伞锥结构地面模型结构尺寸,Dp为加油软管-伞锥结构风洞模型结构尺寸;
λP=λρ
其中,λP为空气动压力系数,λρ为风洞介质密度。
作为上述方法的一种改进,所述结构动力学相似准则满足下式:
其中,λm为伞锥质量系数;
其中,λEJ为加油软管抗弯刚度;
其中,λI为伞锥转动惯量;
其中,λf为减缩频率系数;
λt=λD
其中,λt为时间系数;
作为上述方法的一种改进,所述步骤2)具体包括:
步骤2-1)基于相似准则,将加油机和受油机均简化为刚体,采用飞行动力学的气动力模型分别描述加油机和受油机的气动力和气动力矩;
步骤2-2)利用基于任意欧拉-拉格朗日描述的绝对节点坐标索单元对包括加油软管和伞锥的加油管路进行建模,用于反映加油软管的变形特性以及变长度特性;
步骤2-3)计算加油软管的气动力载荷;
步骤2-4)计算伞锥的气动力载荷。
作为上述方法的一种改进,所述步骤2-1)具体包括:
加油机与受油机的气动力均表示为:
其中,FX、FY和FZ分别表示在x方向、y方向和z方向的气动力,q∞为来流动压,S为特征面积,V为来流速度,b为翼展,c为平均气动弦长,和分别为质心侧向速度和质心升沉速度,为质心升沉加速度,和分别为俯仰角速度、滚转角速度和偏航角速度,CL0、CD0和Cm0分别为零攻角对应的升力系数、阻力系数和俯仰力矩系数,CL、CD、CY分别为升力系数、阻力系数和侧向力系数,α、β和δe分别为迎角、侧滑角和升降舵偏角;上式中
加油机与受油机的气动力矩均表示为:
其中,MX、MY和MZ分别表示在x方向、y方向和z方向的气动力矩,Cl、Cm和Cn分别为滚转力矩系数、俯仰力矩系数和偏航力矩系数,上式中
作为上述方法的一种改进,所述步骤2-2)具体包括:
设置两结点梁索单元的广义坐标q为:
其中,r1,r2分别表示梁索单元的第一结点全局坐标和第二结点全局坐标,T表示转置,p1,p2分别表示梁索单元的第一结点物质坐标和第二结点物质坐标;梁索单元内部任意物质点的位置r满足下式:
r=Neqe
软管与加油机的收放机构相互作用满足基于任意欧拉-拉格朗日描述ALE的梁索单元结点上的约束:
rp-rp0=0
其中,rp0为加油机端软管连接处在全局坐标系下位置,rp为加油机处软管ALE结点p的位置。
作为上述方法的一种改进,所述步骤2-3)具体包括:
将加油软管均分为N段,靠近伞锥为第N段;
根据下式计算加油管路第K段软管的气动力载荷QK:
QK=mKg+(Dl,K+1+Dl,K)/2
其中,0<K<N,mK为第K段软管的质量,g为重力加速度,Dl,K+1为第K+1段软管受到的气动力,Dl,K为第K段软管受到的气动力,满足下式:
Dl,K=Dt,K+Dn,K,
其中,Dt,K为第K段软管受到的切向气动力,Dn,K为第K段软管受到的法向气动力,Dt,K满足下式:
其中,vt,K=[(vK-ωK)·nK1]nK1,为第K段软管相对于其所处流场风速的速度沿软管的切向分量,dK为加油软管的直径,lK为第K段软管的长度,Ct,K是第K段软管的切向气动阻力系数,vK为第K段软管的速度向量,ωK为当地风速向量,为加油机尾涡诱导速度与大气紊流速度的合速度向量,nK1为第K段软管轴线单位矢量,ρ∞为空气密度;
Dn,K满足下式:
其中,vn,K=vK-ωK-vt,K,为第K段软管相对于其所处流场风速的速度沿软管的法向分量,||·||表示取模。
作为上述方法的一种改进,所述步骤2-4)具体包括:
软管与伞锥连接的质点处受到的除拉力的合外力QN为:
QN=(mN/2+mdro)g+Dl,N/2+Ddro
式中,mN为第N段软管的质量,mdro为伞锥的质量,Dl,N为第N段软管受到的气动力,Ddro为伞锥受到的气动力,满足下式:
式中,qdro为伞锥处来流动压,αdro为伞锥的迎角,αdro0为伞锥的平衡迎角,βdro为伞锥的侧滑角,Cdro0为伞锥的3×1基本气动力系数矩阵,Cdroα为伞锥的3×1迎角系数矩阵,Cdroβ为伞锥的3×1侧滑角系数矩阵。
与现有技术相比,本发明的优势在于:
1、本发明基于Π定理首次提出了针对软管式空中加油对接模拟试验装置的相似准则,同时考虑了流体动力学和空气动力学相似;
2、本发明基于相似准则,结合柔性多体动力学,建立空中加油对接模拟试验装置的动力学模型,利用多体动力学求解器计算模拟试验装置的动力学响应,在试验前判断设计的试验装置是否可以复现在加油对接过程中加油管线系统的动态现象及动态特性;
3、基于本发明的设计方法建立的仿真模型,主要用于前期的系统设计与优化,可以较为快速的对设计方案进行多轮优化、改进,使设计不至于盲目,节省前期设计阶段在时间、人力、物力、财力等各方面成本。
附图说明
图1是编队飞行、空中加油中多机干扰特性示意图;
图2是软管式空中加油对接模拟装置设计原理图;
图3是加油系统柔性多体动力学建模示意图;
图4是加油软管-伞锥结构缩比模型相似率示意图。
具体实施方式
本发明提出了一种用于软管式空中加油对接模拟试验装置的设计方法,所述方法包括:
步骤1)建立软管式空中加油对接模拟试验装置的相似准则,所述空中加油对接模拟装置包括加油机、加油软管、伞锥和受油机;
步骤2)基于相似准则,结合柔性多体动力学,建立空中加油对接模拟试验装置的动力学模型;
步骤3)根据模拟试验要求,基于动力学模型,根据相似准则实现软管式空中加油对接模拟装置的设计,利用多体动力学求解器计算受油机对接伞锥后软管的动力学响应,进而实现加油对接过程中加油软管和伞锥的动态现象及动态特性的模拟。
下面结合附图和实施例对本发明的技术方案进行详细的说明。
实施例
本发明的实施例提出了一种用于软管式空中加油对接模拟试验装置的设计方法,以下具体分析。
1、相似律与相似准则
相似律推导与相似准则的建立,是软管式空中加油对接模拟试验装置的重要前提和理论依据。
受限于风洞尺寸及其它技术约束,在风洞中进行试验时需要对模型等相关装置进行必要缩比,以满足试验需求。对于高速风洞定常试验而言,其相似律主要关注来流速度相似(马赫数Ma)、模型几何相似,而忽略了雷诺数等参数相似。对于空中加油这一问题的高速风洞试验而言,涵盖空气动力学、刚体/柔性体动力学、材料学等多个学科,涉及定常、非定常流动现象。因此,须针对这一问题推导相似律、建立相似准则;在此基础上,科学的、有选择的进行相似参数的取舍。使所建立的相似准则,既能准确描述空中加油这一实际过程,又便于风洞试验的工程实现与实际操作。
加油对接过程中,“软管-伞锥”组合体柔性结构所展现出的运动形态,受加油对接操纵、软管柔性材料形变、大气紊流、加油机尾流、姿态变化、机体振动、受油机头波等内外部干扰因素的影响,是典型的非线性柔性结构与空气动力相互耦合作用的复杂问题。在“软管-伞锥”式空中加油风洞试验缩比模型设计时,需要同时考虑流体动力学相似和结构动力学相似。量纲分析法和缩比实验法在工程中应用十分广泛,这种方法的理论基础是Π定理。
在进行风洞试验设计时,常遇到无法同时满足所有物理量相似,则需要根据实际状况及对物理现象的认识,对各物理量相似律进行取舍。图2中对本问题相似律的取舍进行了说明。总起来讲,就是通过质量分布相似和刚度分布相似,保证结构动力学相似,放弃了材料特性相似(图2中椭圆外的虚线框表示放弃的相似关系),包括密度相似和杨氏模量相似,系统的结构阻尼不易控制,也不得不放弃。而流体动力学相似则通过几何外形、速度场、Strouhal数、Ma数、欧拉数来保证,放弃了重力影响、流体阻尼影响的相似。表1备注反映了各物理量相似律的取舍,并给出取舍原因。依据这一相似准则,则可设计本问题的模型、装置及其它相关部件。
表1“软管-伞锥”结构风洞试验缩比模型相似律
2、“软管式空中加油”动力学建模与仿真计算
主要基于柔性多体动力学,建立“加油机-加油软管-伞锥-受油机系统”动力学模型。图3是软管式空中加油对接模拟装置设计原理图。
a)由于主要关注加油机与受油机的整体运动以及“加油软管-伞套”系统的动力学响应,因此飞机模型简化为刚体;加油机与受油机受到的气动力和气动力矩采用飞行力学中的气动力模型给出;利用基于任意欧拉-拉格朗日描述(ALE)的绝对节点坐标(ANCF)索/梁单元对加油管路进行建模,该方法能够反映加油管路的大幅度变形特性以及变长度特性。
通过软管张力控制系统,选取合适的控制律调整软管长度以抑制软管的剧烈摆动及“鞭甩”现象。与传统拉格朗日描述不同,ALE变长度索单元的结点广义坐标不仅包括结点的全局坐标r,还包含物质坐标p,即起始点到结点的绳长,如图4所示。
两结点梁索单元的广义坐标为:
其中:下标1、2表示结点编号。为了描述单元内部任意物质点的位置r,引入形函数Ne:
r=Neqe (2)
其中,
ξ=(2p-p1-p2)/(p1-p2)
与Lagrange索单元不同的是,由于边界物质流动导致ALE索单元形函数的自然坐标ξ是随时间变化的,因此ALE索单元的形函数也是随时间变化的。对式r=Neqe求导
(3)式第一项为局部导数,第二项为由于物质输运引起的迁移导数,此处体现了ALE索单元与Lagrange索单元的差异。软管与加油机收放机构相互作用可以处理成ALE索单元结点上的约束:
rp-rp0=0 (5)
rp0为加油机端软管连接处在全局坐标系下位置,rp是加油机处软管ALE结点p的位置,不约束物质坐标,软管可以沿加油机收放机构滑动。
b)“加油软管-伞套”的气动力建模简化为切向分布气动力和法向分布气动力。切向与法向分布气动力系数通过圆柱绕流、Hoerner切向摩擦经验公式、CFD计算和风洞吹风试验相结合得到。加油机与受油机气动力与气动力矩表示为
其中,分别为质心航向、侧向、升沉速度,分别为俯仰角速度、滚转角速度、偏航角速度,δe是升降舵偏角,q为来流动压,S为特征面积,V为来流速度,b为翼展,c为平均气动弦长,α、β分别为迎角和侧滑角,CL、CD、CY分别为升力系数、阻力系数、侧向力系数,Cl、Cm、Cn分别为滚转力矩系数、俯仰力矩系数、偏航力矩系数,CL0、CD0、Cm0分别为零攻角升力系数、阻力系数、俯仰力矩系数。式中
c)“加油机-加油软管-锥套的气动力载荷计算总体分为软管段的气动力载荷计算和锥套的气动力载荷计算两大部分,而软管段的气动力载荷还可以细分为软管段切向气动力与法向气动力计算。
软管切向气动力的表达式为
QK=mKg+(Dl,K+1+Dl,K)/2 (8)
式中,mK为第K段软管质量,g为重力加速度,Dl,K+1为第K+1段软管受到的法向气动力,Dl,K为第K段软管受到的气动力,Dl,K=Dt,K+Dn,K,其中,Dn,K为第K段软管受到的法向气动力,Dt,K为第K段软管受到的切向气动力,满足下式:
式中,vt,K={(vK-ωK)·nK1}nK1,为第K段软管相对于其所处流场风速的速度沿软管的切向分量。dK为加油软管的直径,lK为第K段软管的长度,Ct,K是第K段软管的切向气动阻力系数,其大小与当地气流雷诺数有关,根据Hoerner关于圆柱体的理论计算作用于软管的切向气动力,获得相应雷诺数下的气动阻力系数:
切向气动力系数Ct,K由当地的气流雷诺数(Re=||vk-ωk||L/γ)决定。ωK为当地风速向量,为加油机尾涡诱导速度与大气紊流速度的合速度向量。γ是空气的粘度系数。特征长度L=πdK/2sinα,其中α是软管和相对气流之间的夹角,nK1为第K段软管轴线单位矢量,ρ∞为空气密度。
软管法向气动力计算式14中
Dl,K=Dt,K+Dn,K (11)
其中Dn,K为第K段软管受到的法向气动力,同样根据流体力学,作用于物体前端的压力大于其后端,这样会使得一种压差阻力出现在物体表面即为法向气动力或法向摩擦力,其表示为
式中,vn,K=vK-ωK-vt,K,为第K段软管相对于其所处流场风速的速度沿软管的法向分量。dK为加油软管的直径,lK为第K段软管的长度,Cn,K分别是第K段软管的法向气动阻力系数。同样,在计算法向气动阻力的时候,不同的雷诺数对应的压力阻力系数如下:
其中Re=||vn,k||dk/γ。ωK同样为当地风速,即为加油机尾涡诱导速度与大气紊流速度的合速度。γ是空气的粘度系数。
d)在飞行过程中,锥套同样会受到外力的影响,最后一段软管与锥套连接的质点处受到的除拉力的合外力可表示为
QN=(mN+mdro)g+Dl,N/2+Ddro (14)
式中,Dl,N为第N段即最后一段软管受到的气动力,Ddro为可控锥套受到的气动力,mN为第N段即最后一段软管质量的一半,mdro为伞锥套的质量,其受力表达式为
式中,q为动压,α为锥套的迎角,其中α=tan-1(vwdz/vwdx),α0表示锥套的平衡迎角,β为锥套的侧滑角,其中β=sin-1(vwdz/||vwd||),vwd为锥套相对于气流的速度向量在锥套体坐标系Sd下的表示vwdx、vwdy、vwdz为vwd在锥套体坐标系下的三个相互垂直的分量。C0=[CD0CY0 CL0]T为锥套的3×1基本气动力系数矩阵。Cα=[CDα CYα CLα]T为锥套的3×1迎角系数矩阵,Cβ=[CDβ CYβ CLβ]T为锥套的3×1侧滑角系数矩阵。
e)在“加油机-加油软管-伞锥-受油机”系统动力学模型基础上,利用多体动力学求解器计算受油机对接伞锥后软管的动力学响应,进而复现在加油对接过程中加油管线系统的动态现象及动态特性。
基于上述动力学理论建立的仿真模型,主要用于前期的系统设计与优化,可以较为快速的对设计方案进行多轮优化、改进,使设计不至于盲目,节省前期设计阶段在时间、人力、物力、财力等各方面成本。
最后所应说明的是,以上实施例仅用以说明本发明的技术方案而非限制。尽管参照实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,对本发明的技术方案进行修改或者等同替换,都不脱离本发明技术方案的精神和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (6)
1.一种用于软管式空中加油对接模拟试验装置的设计方法,所述方法包括:
步骤1)建立软管式空中加油对接模拟试验装置的相似准则,所述空中加油对接模拟装置包括加油机、加油软管、伞锥和受油机;
步骤2)基于相似准则,结合柔性多体动力学,建立空中加油对接模拟试验装置的动力学模型;
步骤3)根据模拟试验要求,基于动力学模型,根据相似准则实现软管式空中加油对接模拟装置的设计,利用多体动力学求解器计算受油机对接伞锥后软管的动力学响应,进而实现加油对接过程中加油软管和伞锥的动态现象及动态特性的模拟;
所述步骤1)具体包括:
建立流体动力学相似准则和结构动力学相似准则,其中,
所述流体动力学相似准则包括:几何外形相似、速度场相似、动压相似、流场压缩性相似和重力影响相似;
所述结构动力学相似准则包括:质量分布相似、固有频率相似、振型相似和刚度分布相似;
所述流体动力学相似准则满足下式:
λV=1
其中,λV为运动速度系数;
λα=1
其中,λα为角度系数
其中,λg为重力加速度系数,λD为结构尺寸缩比系数;
其中,Dm为加油软管和伞锥结构的地面模型结构尺寸,Dp为加油软管和伞锥结构的风洞模型结构尺寸;
λP=λρ
其中,λP为空气动压力系数,λρ为风洞介质密度;
所述结构动力学相似准则满足下式:
其中,λm为伞锥质量系数;
其中,λEJ为加油软管抗弯刚度;
其中,λI为伞锥转动惯量;
其中,λf为减缩频率系数;
λt=λD
其中,λt为时间系数;
2.根据权利要求1所述的用于软管式空中加油对接模拟试验装置的设计方法,其特征在于,所述步骤2)具体包括:
步骤2-1)基于相似准则,将加油机和受油机均简化为刚体,采用飞行动力学的气动力模型分别描述加油机和受油机的气动力和气动力矩;
步骤2-2)利用基于任意欧拉-拉格朗日描述的绝对节点坐标索单元对包括加油软管和伞锥的加油管路进行建模,用于反映加油软管的变形特性以及变长度特性;
步骤2-3)计算加油软管的气动力载荷;
步骤2-4)计算伞锥的气动力载荷。
3.根据权利要求2所述的用于软管式空中加油对接模拟试验装置的设计方法,其特征在于,所述步骤2-1)具体包括:
加油机与受油机的气动力均表示为:
其中,FX、FY和FZ分别表示在x方向、y方向和z方向的气动力,q∞为来流动压,S为特征面积,V为来流速度,b为翼展,c为平均气动弦长,和分别为质心侧向速度和质心升沉速度,为质心升沉加速度,和分别为俯仰角速度、滚转角速度和偏航角速度,CL0、CD0和Cm0分别为零攻角对应的升力系数、阻力系数和俯仰力矩系数,CL、CD和CY分别为升力系数、阻力系数和侧向力系数,α、β和δe分别为迎角、侧滑角和升降舵偏角;上式中
加油机与受油机的气动力矩均表示为:
其中,MX、MY和MZ分别表示在x方向、y方向和z方向的气动力矩,Cl、Cm和Cn分别为滚转力矩系数、俯仰力矩系数和偏航力矩系数,上式中
4.根据权利要求2所述的用于软管式空中加油对接模拟试验装置的设计方法,其特征在于,所述步骤2-2)具体包括:
设置两结点梁索单元的广义坐标q为:
其中,r1,r2分别表示梁索单元的第一结点全局坐标和第二结点全局坐标,T表示转置,p1,p2分别表示梁索单元的第一结点物质坐标和第二结点物质坐标;梁索单元内部任意物质点的位置r满足下式:
r=Neqe
软管与加油机的收放机构相互作用满足基于任意欧拉-拉格朗日描述ALE的梁索单元结点上的约束:
rp-rp0=0
其中,rp0为加油机端软管连接处在全局坐标系下位置,rp为加油机处软管ALE结点p的位置。
5.根据权利要求2所述的用于软管式空中加油对接模拟试验装置的设计方法,其特征在于,所述步骤2-3)具体包括:
将加油软管均分为N段,靠近伞锥为第N段;
根据下式计算加油管路第K段软管的气动力载荷QK:
QK=mKg+(Dl,K+1+Dl,K)/2
其中,0<K<N,mK为第K段软管的质量,g为重力加速度,Dl,K+1为第K+1段软管受到的气动力,Dl,K为第K段软管受到的气动力,满足下式:
Dl,K=Dt,K+Dn,K,
其中,Dt,K为第K段软管受到的切向气动力,Dn,K为第K段软管受到的法向气动力,Dt,K满足下式:
其中,vt,K=[(vK-ωK)·nK1]nK1,为第K段软管相对于其所处流场风速的速度沿软管的切向分量,dK为加油软管的直径,lK为第K段软管的长度,Ct,K是第K段软管的切向气动阻力系数,vK为第K段软管的速度向量,ωK为当地风速向量,为加油机尾涡诱导速度与大气紊流速度的合速度向量,nK1为第K段软管轴线单位矢量,ρ∞为空气密度;
Dn,K满足下式:
其中,vn,K=vK-ωK-vt,K,为第K段软管相对于其所处流场风速的速度沿软管的法向分量,||·||表示取模,Cn,K是第K段软管的法向气动阻力系数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210277521.6A CN114476123B (zh) | 2022-03-21 | 2022-03-21 | 一种用于软管式空中加油对接模拟试验装置的设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210277521.6A CN114476123B (zh) | 2022-03-21 | 2022-03-21 | 一种用于软管式空中加油对接模拟试验装置的设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114476123A CN114476123A (zh) | 2022-05-13 |
CN114476123B true CN114476123B (zh) | 2023-04-14 |
Family
ID=81488098
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210277521.6A Active CN114476123B (zh) | 2022-03-21 | 2022-03-21 | 一种用于软管式空中加油对接模拟试验装置的设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114476123B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115597819B (zh) * | 2022-11-30 | 2023-03-21 | 中国空气动力研究与发展中心低速空气动力研究所 | 一种风洞模拟空中加油软管-锥套收放装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104180808B (zh) * | 2014-08-05 | 2017-02-15 | 南京航空航天大学 | 一种用于自主空中加油的圆形锥套视觉位姿解算方法 |
CN110110404A (zh) * | 2019-04-24 | 2019-08-09 | 中国船舶重工集团公司第七一九研究所 | 一种空间拖曳系统的介质阻力获取方法及装置 |
CN113335559A (zh) * | 2021-06-01 | 2021-09-03 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种空中加油软管锥套动态的模拟方法 |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110287643B (zh) * | 2019-07-09 | 2023-04-07 | 中国空气动力研究与发展中心高速空气动力研究所 | 空腔流致振动与流致噪声耦合特性的风洞实验模拟方法 |
CN111103905B (zh) * | 2019-12-26 | 2021-08-10 | 南京航空航天大学 | 一种基于电机积分鲁棒控制驱动的加油软管振动抑制方法 |
CN111159920B (zh) * | 2020-01-13 | 2022-03-15 | 西北工业大学 | 一种自主空中加油软管锥套动力学建模方法 |
CN111258231B (zh) * | 2020-01-20 | 2021-03-30 | 北京航空航天大学 | 仿猛禽视觉导航的自主空中加油对接半物理系统及其方法 |
CN111964862A (zh) * | 2020-07-20 | 2020-11-20 | 中国航天空气动力技术研究院 | 一种机弹分离动力学相似风洞试验方法 |
CN112069589B (zh) * | 2020-08-05 | 2022-10-11 | 中国航空工业集团公司哈尔滨空气动力研究所 | 一种低速风洞空中加油软管-锥套动态特性试验方法 |
CN111649906B (zh) * | 2020-08-06 | 2021-08-10 | 中国空气动力研究与发展中心低速空气动力研究所 | 一种风洞试验的受油机模型支撑装置的滑膜跟随控制方法 |
CN112947522B (zh) * | 2021-02-10 | 2022-08-02 | 西北工业大学 | 一种基于有限时间观测器的硬式空中加油姿态控制方法 |
CN113806871B (zh) * | 2021-11-18 | 2022-02-08 | 北京航空航天大学 | 一种考虑结构非线性的柔性飞行动力学建模方法 |
-
2022
- 2022-03-21 CN CN202210277521.6A patent/CN114476123B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104180808B (zh) * | 2014-08-05 | 2017-02-15 | 南京航空航天大学 | 一种用于自主空中加油的圆形锥套视觉位姿解算方法 |
CN110110404A (zh) * | 2019-04-24 | 2019-08-09 | 中国船舶重工集团公司第七一九研究所 | 一种空间拖曳系统的介质阻力获取方法及装置 |
CN113335559A (zh) * | 2021-06-01 | 2021-09-03 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种空中加油软管锥套动态的模拟方法 |
Also Published As
Publication number | Publication date |
---|---|
CN114476123A (zh) | 2022-05-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103970957B (zh) | 一种弹性乘波体高超声速飞行器仿真方法 | |
Isogai | On the transonic-dip mechanism of flutter of a sweptback wing | |
Li et al. | Dynamics modeling and simulation of flexible airships | |
Williamson et al. | Controllable drogue for automated aerial refueling | |
Kumar et al. | Mechanism of side force generation and its alleviation over a slender body | |
CN114476123B (zh) | 一种用于软管式空中加油对接模拟试验装置的设计方法 | |
CN110162929A (zh) | 一种自主空中加油软管锥套平衡位置计算方法 | |
Singh et al. | Dynamic Modelling and control of an aerodynamically controlled capturing device for ‘in-air-capturing’of a reusable launch vehicle | |
Vicroy | Blended-wing-body low-speed flight dynamics: summary of ground tests and sample results | |
Mahdi et al. | Effects of coarse riblets on air flow structures over a slender delta wing using particle image velocimetry | |
Peake et al. | Three-dimensional flow separations on aircraft and missiles | |
Nielsen | Nonlinearities in missile aerodynamics | |
Keener et al. | Side forces on Forebodies at high angles of attack and Mach numbers from 0.1 to 0.7: two tangent ogives, paraboloid and cone | |
Amoozgar et al. | The Effect of Thrust Vectoring on Aeroelastic Stability of Electric Aircraft | |
Fujino et al. | Flutter characteristics of an over-the-wing engine mount business-jet configuration | |
Khaware et al. | Numerical simulation of store separation trajectories for EGLIN test case using overset mesh | |
Thompson | The unsteady aerodynamics of a delta wing undergoing large-amplitude pitching motions | |
Stainback et al. | The Langley Research Center's Low-Turbulence Pressure Tunnel | |
Ericsson et al. | Effect of nose slenderness on forebody flow control | |
Ding et al. | Studying the Relationship Between Aircraft Stability and Shape Based on Hydrodynamics | |
Roy et al. | Effect of Anhedral–Dihedral Angles on 76-40° Double-Delta Wing at Low Speeds | |
Chu et al. | Overall Layout Optimization of an Unmanned Hovercraft Considering Multiple Aerodynamic Interferences | |
Su et al. | Moving platform lectotype and aerodynamic characteristics of Flying Cable-Driven Parallel Robot | |
Rui et al. | Determination of take-off constraints for the aerospace vehicle | |
Wangyi et al. | Aerodynamic Modeling and Analysis of Receiver Under the Disturbance of Tanker’s Wake Flow |
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 |