CN114970402A - 飞机地面振动试验中气动载荷模拟计算系统及方法 - Google Patents
飞机地面振动试验中气动载荷模拟计算系统及方法 Download PDFInfo
- Publication number
- CN114970402A CN114970402A CN202210770406.2A CN202210770406A CN114970402A CN 114970402 A CN114970402 A CN 114970402A CN 202210770406 A CN202210770406 A CN 202210770406A CN 114970402 A CN114970402 A CN 114970402A
- Authority
- CN
- China
- Prior art keywords
- function
- time
- test
- test piece
- aerodynamic force
- 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
Links
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/15—Vehicle, aircraft or watercraft design
-
- 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)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Fluid Mechanics (AREA)
- Mathematical Physics (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Automation & Control Theory (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computational Mathematics (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
本发明提供了飞机地面振动试验中气动载荷模拟计算系统及方法,属于飞机振动试验技术领域。系统包括:计算子系统、地面振动试验子系统、半实物仿真子系统。方法包括以下步骤:S1、得到适用于飞机地面振动试验整个时变历程的非定常气动力快速计算模型;S2、对非定常气动力快速计算模型进行试验验证;S3、采用非定常气动力快速计算模型进行气动载荷模拟计算。本发明解决了当前飞机地面振动试验模拟只能应用于定常系统的问题,具有有效完善当前飞机地面振动试验模拟技术,推动飞机地面振动试验模拟工程应用进展的优点。
Description
技术领域
本发明涉及飞机振动试验技术领域,具体是涉及飞机地面振动试验中气动载荷模拟计算系统及方法。
背景技术
飞机地面振动试验模拟技术作为近年来发展的一种以真实结构作为验证对象的半物理仿真颤振验证技术,具有风险小、周期短、精度高的优势。该技术利用传感器测量结构的振动响应信号,将其代入预先计算得到的非定常气动力模型获得结构的等效气动力,再采用激振器将等效气动力实时加载在试验结构上,形成闭环系统在地面实现对结构颤振现象的模拟。通过对试验结构施加初始扰动并观察其响应,即可开展针对真实结构颤振边界测试验证。
在上述飞机地面振动试验模拟的步骤中,考虑到常规结构的颤振发生频率,试验系统完成一次循环计算的周期应不高于0.002秒。为了满足飞机地面振动试验模拟系统对计算效率的需求,飞机地面振动试验模拟设计时依据结构动力学特性对频域气动力影响系数矩阵进行降阶缩聚处理,进而采用有理函数近似以及拉氏变换得到由数个低阶矩阵组成的时域非定常气动力模型。
上述建模方法能够成立的一个大前提是结构的动力学特性不发生改变,而在实际飞行器结构中无论是舵面偏转影响或气动热效应导致的结构热模态时变特性,均会使现有非定常气动力计算方法不再完全适用。
本发明提出一种基于Kriging代理模型的飞机地面振动试验中气动载荷模拟计算系统及方法。本发明依据时间冻结假设建立多个离散时间点的非定常气动力模型后,采用Kriging代理模型对离散时间点的非定常气动力模型进行拟合建模,使得新得到的气动力模型在保证系统运行效率的同时,针对动力学特性时变的结构具有更高的计算精度,有效扩充了地面颤振模拟试验的应用范围。
发明内容
针对当前地面颤振模拟试验只能应用于定常系统的缺陷,本发明给出了可应用于动力学特性时变结构的非定常气动力快速计算方法,有效完善当前地面颤振模拟试验技术,推动其工程应用进展。
为解决上述问题,本发明的技术方案如下:
本发明提供了飞机地面振动试验中气动载荷模拟计算系统,包括:
用于进行飞机地面振动试验气动载荷模拟的计算子系统,计算子系统包括:用于通过控制信号控制地面振动试验运行且搭载有非定常气动力快速计算模型的加载控制器,所述加载控制器还搭载有根据测量得到的试验件响应信号变化趋势自动调整所述非定常气动力快速计算模型中的风速,以实现对颤振临界速度的实时辨识与记录的颤振边界追踪模块,
用于对所述非定常气动力快速计算模型进行验证的地面振动试验子系统,所述地面振动试验子系统包括:与所述加载控制器电性连接的:位于试验件底部并为试验件提供颤振力的激振器,布设在试验件上试验件控制点处的激光位移传感器和加速度传感器,布设在试验件上方的加热器,用于将控制信号放大的功率放大器,
用于将所述非定常气动力快速计算模型计算的等效气动力转化为激振器控制信号的半实物仿真子系统,半实物仿真子系统通过NI VeriStand仿真配置软件对接Matlab.simulink。
本发明还提供飞机地面振动试验中气动载荷模拟计算方法,基于上述飞机地面振动试验中气动载荷模拟计算系统,包括以下步骤:
S1、基于时间冻结假设,通过将多个离散时间点的非定常气动力快速计算模型进行Kriging代理模型拟合,得到适用于飞机地面振动试验整个时变历程的非定常气动力快速计算模型;
S2、通过实际飞机地面振动试验对步骤S1建立的非定常气动力快速计算模型进行试验验证,具体包括以下步骤:
S2-1、对飞机地面振动试验颤振计算中的气动平面和结构间的插值点进行缩聚,通过优化方法选取最优插值点,将最优插值点作为试验件控制点,
S2-2、选择数个采样时间点,将数个采样时间点数据导入步骤S1,得到适用于试验件控制点的非定常气动力快速计算模型,
S2-3、在实验室内搭建地面振动试验子系统,在试验件上方布设加热器,在试验件底部固定激振器以及与激振器电性连接的功率放大器,在试验件上试验件控制点布置激光位移传感器和加速度传感器,功率放大器、激光位移传感器、加速度传感器与加载控制器电性连接,
S2-4、通过激光位移传感器和加速度传感器将获取的试验件控制点响应信号传入加载控制器,加载控制器将试验件控制点响应信号代入步骤S2-1建立的非定常气动力快速计算模型,非定常气动力快速计算模型计算飞机地面振动试验时变历程中试验件的等效气动力,并将等效气动力传入半实物仿真子系统,半实物仿真子系统通过NI VeriStand仿真配置软件对接Matlab.simulink,将试验件的等效气动力转化为激振器控制信号,加载控制器通过激振器控制信号控制激振器将等效气动力实时加载在试验件上,加载完成后回到步骤S2-2,直至飞机振动试验等效模拟过程结束,在此过程中,加载控制器记录飞机振动试验等效模拟过程中的颤振临界速度,
S2-5、通过地面振动试验子系统开展实际飞机地面振动试验,并采用颤振边界追踪模块记录试验件在飞机地面颤振试验时变历程中的颤振临界速度,与步骤S2-4中非定常气动力快速计算模型计算得到的飞机地面振动试验时变历程中的颤振临界速度进行比对,以判定非定常气动力快速计算模型的准确性,若非定常气动力快速计算模型的准确性满足要求,得到试验验证后的非定常气动力快速计算模型,若不满足要求,加密采样时间点回到步骤S2-2,直至非定常气动力快速计算模型的准确性满足要求;
S3、采用步骤S2试验验证后的非定常气动力快速计算模型进行飞机地面振动试验中气动载荷模拟计算。
进一步地,步骤S2-1中最优插值点为:激振点和拾振点。
进一步地,步骤S1具体包括以下步骤:
S1-1、建立试验件控制点承受非定常气动力的计算公式,
S1-2、建立离散时间点的试验件控制点应承受非定常气动力计算模型,
S1-3、结合Kriging代理模型建立对应随机函数
S1-4、建立静态随机过程函数在任意两个采样时间点间的协方差计算模型
在飞机地面振动试验定义的时间范围内,得到的静态随机过程函数在任意两个采样时间点间的协方差计算模型如下:
上式中,为方差,Z(t1)为第一采样时间点的静态随机过程函数,Z(t2)为第二采样时间点的静态随机过程函数,t1为第一采样时间点,t2为第二采样时间点,R(t1,t2)为相关函数,相关函数的形式为高斯函数或三次样条函数,
S1-5、建立非定常气动力快速计算模型,即Kriging代理模型在任意时刻t的预测值计算模型
通过公式(4)构建形如式(3)的样本点数据列向量定义公式如下:
[Aij(t)]s=[Aij(t1)Aij(t2)…Aij(tn)]T=βG+z (5)
上式中,[Aij(t)]s为样本点数据列向量,β为基函数系数,[Aij(t1)Aij(t2)…Aij(tn)]为n个采样时间点的样本点数据,G为基函数列向量,G=[g(t1)g(t2)…g(tn)]T,z为随机函数列向量,z=[Z(t1),Z(t2)…Z(tn)]T,
为满足Kriging代理模型无偏估计E(Z(t))=0的要求,根据公式(5)与公式(8)可得基函数与基函数列向量关系的计算公式:
根据公式(3)、公式(8)与公式(9)得到Kriging代理模型的预测值误差计算公式如下:
上式中,为采样时间点外的时间节点函数值,为Aij(t)通过Kriging代理模型得到的对应随机函数,Aij(t)为任一时变气动力影响系数矩阵中的元素,Z为随机函数列向量,Z(t)是均值为0、方差为的静态随机过程函数,为加权行向量,
根据Kriging代理模型预测值满足无偏估计且均方误差最小的要求,即为在GTct=g(t)的前提下,寻找满足均方误差最小的加权列向量ct,为此引入拉格朗日乘数λ构建拉格朗日函数H(ct,λ)如下:
分别对公式(12)中加权列向量ct和拉格朗日乘数λ求偏导数,得到满足要求的加权列向量ct计算公式为:
将公式(13)带入公式(8)得到非定常气动力快速计算模型,即Kriging代理模型在任意时刻t的非定常气动力预测值计算模型如下:
由公式(5)与公式(13)可知,公式(14)中除ct表达式中的rt外其他参数只与采样时间点有关,因此可以提前计算得到,而rt需要将当前时刻t代入相关函数公式(7)求解,模型的输入参数为时间,因此基于公式(14)可以获得任意时刻下与结构特性相匹配的气动力。
进一步地,步骤S1-1具体包括以下内容:
根据飞机地面振动试验中试验件的结构振动响应计算试验件控制点承受的非定常气动力,试验件控制点承受非定常气动力的计算公式如下:
fs=q∞A(k)zs (1)
上式中,fs为试验件控制点应承受的非定常气动力,q∞为来流动压,A(k)为降阶后的气动力影响系数矩阵,A(k)是通过与马赫数相匹配的非定常气动力理论计算后并依据结构动力学特性进行降阶处理得到的,Zs表示试验件的结构振动位移,k为减缩频率。
更进一步地,减缩频率k的计算公式如下:
上式中,k为减缩频率,ω为简谐运动的圆频率,b为参考长度,对于飞机的翼面类结构,参考长度取根弦长度的一半,V为来流速度。
优选地,步骤S1-2具体包括以下内容:
由于实际飞机地面振动试验需要试验件控制点应承受非定常气动力的时域关系式,将公式(1)通过最小状态法有理函数近似延拓至拉氏域,再通过拉氏反变换得到离散时间点的试验件控制点应承受非定常气动力计算模型如下:
上式中,fs(t)为与时间相关的试验件控制点应承受非定常气动力,A0为位移气动力影响系数矩阵,A1为速度气动力影响系数矩阵,A2为加速度气动力影响系数矩阵,B1为第一辅助气动力影响系数矩阵,B2为第二辅助气动力影响系数矩阵,B3为第三辅助气动力影响系数矩阵,ZS为控制点位移,为控制点速度,为控制点加速度,q∞为来流动压,V为来流速度,b为参考长度,对于飞机的翼面类结构,参考长度取根弦长度的一半。
由于公式(2)的建模方法只能适用于时不变系统,而在部分实际飞机地面振动试验各个时间节点的气动力影响系数矩阵数值均会发生变化,因此,需在公式(2)的基础上引入Kriging代理模型,将公式(1)和公式(2)的求解过程替换为简单矩阵运算,实现针对时变动力学特性结构的非定常气动力快速计算。
优选地,步骤S1-5中:
相关向量rt的计算公式如下:
rt=[R(t1,t)R(t2,t)…R(tn,t)]T (7)
上式中,t为n个时间采样点中任意一个样本点,t∈[t1 t2 … tn],rt为相关向量,R(t1,t)为n个时间采样点中任意一个样本点t与第1个采样时间点的相关函数,R(t2,t)为n个时间采样点中任意一个样本点t与第2个采样时间点的相关函数,R(tn,t)为n个时间采样点中任意一个样本点t与第n个采样时间点的相关函数。
进一步优选地,步骤S1-5还包括:
根据Kriging代理模型预测值满足无偏估计且均方误差最小的要求,通过公式(10)得到Kriging代理模型的均方误差计算公式如下:
上式中,为Kriging代理模型的预均方误差,为采样时间点外的时间节点函数,为Aij(t)通过Kriging代理模型得到的对应随机函数,Aij(t)为任一时变气动力影响系数矩阵中的元素,为加权行向量,E为数学期望,为方差,ct为加权列向量,rt为相关向量,为n个时间采样点之间相关性的相关矩阵。
本发明的有益效果是:
(1)本发明将Kriging代理模型引入地面颤振模拟试验,依据时间冻结假设获取采样时间点数据并提出了根据采样时间点数据的完整建模方法,获得了可应用于动力学特性时变结构的非定常气动力快速计算方法,增强了此类试验系统对于时变颤振系统的处理能力;
(2)本发明所需工作模块只是在现有结构时不变系统所采用非定常气动力模块的基础上增加数个低阶矩阵运算,对地面颤振模拟试验计算子系统运算效率影响极小,能够满足地面颤振模拟试验对于气动力运算实时性的要求,此外,由于Kriging模型在应对非线性系统时具有明显优势,因此在计算非定常气动力时拟合精度较高,能够保证在试验件时变历程中都保持较高的非定常气动力计算精度。
附图说明
图1是实施例1飞机地面振动试验中气动载荷模拟计算系统的系统结构图;
图2是实施例2飞机地面振动试验中气动载荷模拟计算方法流程图;
图3是实施例2中步骤S1的流程图;
图4是实施例2中离散时间点的试验件控制点应承受非定常气动力计算模型结构图;
图5是实施例2中最优插值点位置示意图。
具体实施方式
为了使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明作进一步地详细描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例,都属于本发明保护的范围。
在本发明实施例中使用的术语是仅仅出于描述特定实施例的目的,而非旨在限制本发明。在本发明实施例和所附权利要求书中所使用的单数形式的“一种”、“所述”和“该”也旨在包括多数形式,除非上下文清楚地表示其他含义,“多种”一般包含至少两种。
实施例1
本实施例为飞机地面振动试验中气动载荷模拟计算系统,如图1所示,包括:
用于进行飞机地面振动试验气动载荷模拟的计算子系统,计算子系统包括:计算子系统包括:用于通过控制信号控制地面振动试验运行且搭载有非定常气动力快速计算模型的加载控制器,所述加载控制器还搭载有根据测量得到的试验件响应信号变化趋势自动调整所述非定常气动力快速计算模型中的风速,以实现对颤振临界速度的实时辨识与记录的颤振边界追踪模块,
非定常气动力快速计算模型的模型计算公式,即Kriging代理模型在任意时刻t的非定常气动力预测值计算模型如下:
用于对非定常气动力快速计算模型进行验证的地面振动试验子系统,地面振动试验子系统包括:与加载控制器电性连接的:位于试验件底部并为试验件提供颤振力的激振器,布设在试验件上试验件控制点处的激光位移传感器和加速度传感器,布设在试验件上方的加热器,用于将控制信号放大的功率放大器,
用于将所述非定常气动力快速计算模型计算的等效气动力转化为激振器控制信号的半实物仿真子系统,半实物仿真子系统通过NI VeriStand仿真配置软件对接Matlab.simulink。
实施例2
本实施例为飞机地面振动试验中气动载荷模拟计算方法,基于实施例1飞机地面振动试验中气动载荷模拟计算系统,如图2所示,包括以下步骤:
S1、基于时间冻结假设,通过将多个离散时间点的非定常气动力快速计算模型进行Kriging代理模型拟合,得到适用于飞机地面振动试验整个时变历程的非定常气动力快速计算模型,如图3所示,具体包括以下步骤:
S1-1、建立试验件控制点承受非定常气动力的计算公式
根据飞机地面振动试验中试验件的结构振动响应计算试验件控制点承受的非定常气动力,试验件控制点承受非定常气动力的计算公式如下:
fs=q∞A(k)zs (1)
上式中,fs为试验件控制点应承受的非定常气动力,q∞为来流动压,A(k)为降阶后的气动力影响系数矩阵,A(k)是通过与马赫数相匹配的非定常气动力理论计算后并依据结构动力学特性进行降阶处理得到的,ZS表示试验件的结构振动位移,k为减缩频率,减缩频率k的计算公式如下:
上式中,k为减缩频率,ω为简谐运动的圆频率,b为参考长度,对于飞机的翼面类结构,参考长度取根弦长度的一半,V为来流速度,
S1-2、建立离散时间点的试验件控制点应承受非定常气动力计算模型
由于实际飞机地面振动试验需要试验件控制点应承受非定常气动力的时域关系式,将公式(1)通过最小状态法有理函数近似延拓至拉氏域,再通过拉氏反变换得到离散时间点的试验件控制点应承受非定常气动力计算模型如下:
上式中,fs(t)为与时间相关的试验件控制点应承受非定常气动力,A0为位移气动力影响系数矩阵,A1为速度气动力影响系数矩阵,A2为加速度气动力影响系数矩阵,B1为第一辅助气动力影响系数矩阵,B2为第二辅助气动力影响系数矩阵,B3为第三辅助气动力影响系数矩阵,Zs为控制点位移,为控制点速度,为控制点加速度,q∞为来流动压,V为来流速度,b为参考长度,对于飞机的翼面类结构,参考长度取根弦长度的一半,离散时间点的试验件控制点应承受非定常气动力计算模型结构如图4所示,
S1-3、结合Kriging代理模型建立对应随机函数
S1-4、建立静态随机过程函数在任意两个采样时间点间的协方差计算模型
在飞机地面振动试验定义的时间范围内,得到的静态随机过程函数在任意两个采样时间点间的协方差计算模型如下:
上式中,为方差,Z(t1)为第一采样时间点的静态随机过程函数,Z(t2)为第二采样时间点的静态随机过程函数,t1为第一采样时间点,t2为第二采样时间点,R(t1,t2)为相关函数,相关函数的形式为高斯函数或三次样条函数,
S1-5、建立非定常气动力快速计算模型,即Kriging代理模型在任意时刻t的预测值计算模型
通过公式(4)构建形如式(3)的样本点数据列向量定义公式如下:
[Aij(t)]s=[Aij(t1)Aij(t2)…Aij(tn)]T=βG+z (5)
上式中,[Aij(t)]s为样本点数据列向量,β为基函数系数,[Aij(t1)Aij(t2)…Aij(tn)]为n个采样时间点的样本点数据,G为基函数列向量,G=[g(t1)g(t2)…g(tn)]T,z为随机函数列向量,z=[Z(t1)Z(t2)…Z(tn)]T,
相关向量rt的计算公式如下:
rt=[R(t1,t)R(t2,t)…R(tn,t)]T (7)上式中,t为n个时间采样点中任意一个样本点,t∈[t1 t2…tn],rt为相关向量,R(t1,t)为n个时间采样点中任意一个样本点t与第1个采样时间点的相关函数,R(t2,t)为n个时间采样点中任意一个样本点t与第2个采样时间点的相关函数,R(tn,t)为n个时间采样点中任意一个样本点t与第n个采样时间点的相关函数,
为满足Kriging代理模型无偏估计E(Z(t))=0的要求,根据公式(5)与公式(8)可得基函数与基函数列向量关系的计算公式:
根据公式(3)、公式(8)与公式(9)得到Kriging代理模型的预测值误差计算公式如下:
上式中,为采样时间点外的时间节点函数值,为Aij(t)通过Kriging代理模型得到的对应随机函数,Aij(t)为任一时变气动力影响系数矩阵中的元素,z为随机函数列向量,Z(t)是均值为0、方差为的静态随机过程函数,为加权行向量,
根据Kriging代理模型预测值满足无偏估计且均方误差最小的要求,通过公式(10)得到Kriging代理模型的均方误差计算公式如下:
上式中,为Kriging代理模型的预均方误差,为采样时间点外的时间节点函数,为Aij(t)通过Kriging代理模型得到的对应随机函数,Aij(t)为任一时变气动力影响系数矩阵中的元素,为加权行向量,E为数学期望,为方差,ct为加权列向量,rt为相关向量,为n个时间采样点之间相关性的相关矩阵,
根据Kriging代理模型预测值满足无偏估计且均方误差最小的要求,即为在GTct=g(t)的前提下,寻找满足均方误差最小的加权列向量ct,为此引入拉格朗日乘数λ构建拉格朗日函数H(ct,λ)如下:
分别对公式(12)中加权列向量ct和拉格朗日乘数λ求偏导数,得到满足要求的加权列向量ct计算公式为:
将公式(13)带入公式(8)得到非定常气动力快速计算模型,即Kriging代理模型在任意时刻t的预测值计算模型如下:
S2、通过实际飞机地面振动试验对步骤S1建立的非定常气动力快速计算模型进行试验验证,具体包括以下步骤:
S2-1、如图5所示,对飞机地面振动试验颤振计算中的气动平面和结构间的插值点进行缩聚,通过优化方法选取最优插值点,最优插值点为激振点和拾振点,将最优插值点作为试验件控制点,
S2-2、选择15个采样时间点,将15个采样时间点数据导入步骤S1,得到适用于试验件控制点的非定常气动力快速计算模型,
S2-3、在实验室内搭建地面振动试验子系统,在试验件上方布设加热器,在试验件底部固定激振器以及与激振器电性连接的功率放大器,在试验件上试验件控制点布置激光位移传感器和加速度传感器,功率放大器、激光位移传感器、加速度传感器与加载控制器电性连接,
S2-4、通过激光位移传感器和加速度传感器将获取的试验件控制点响应信号传入加载控制器,加载控制器将试验件控制点响应信号代入步骤S2-1建立的非定常气动力快速计算模型,非定常气动力快速计算模型计算飞机地面振动试验时变历程中试验件的等效气动力,并将等效气动力传入半实物仿真子系统,半实物仿真子系统通过NI VeriStand仿真配置软件对接Matlab.simulink,将试验件的等效气动力转化为激振器控制信号,加载控制器通过激振器控制信号控制激振器将等效气动力实时加载在试验件上,加载完成后回到步骤S2-2,直至飞机振动试验等效模拟过程结束,在此过程中,加载控制器记录飞机振动试验等效模拟过程中的颤振临界速度,
S2-5、通过地面振动试验子系统开展实际飞机地面振动试验,并采用颤振边界追踪模块记录试验件在飞机地面颤振试验时变历程中的颤振临界速度,与步骤S2-4中非定常气动力快速计算模型计算得到的飞机地面振动试验时变历程中的颤振临界速度进行比对,以判定非定常气动力快速计算模型的准确性,若非定常气动力快速计算模型的准确性满足要求,得到试验验证后的非定常气动力快速计算模型,若不满足要求,加密采样时间点回到步骤S2-2,直至非定常气动力快速计算模型的准确性满足要求;
S3、采用步骤S2试验验证后的非定常气动力快速计算模型进行飞机地面振动试验中气动载荷模拟计算。
Claims (8)
1.飞机地面振动试验中气动载荷模拟计算系统,其特征在于,包括:
用于进行飞机地面振动试验气动载荷模拟的计算子系统,所述计算子系统包括:用于通过控制信号控制地面振动试验运行且搭载有非定常气动力快速计算模型的加载控制器,所述加载控制器还搭载有根据测量得到的试验件响应信号变化趋势自动调整所述非定常气动力快速计算模型中的风速,以实现对颤振临界速度的实时辨识与记录的颤振边界追踪模块,
用于对所述非定常气动力快速计算模型进行验证的地面振动试验子系统,所述地面振动试验子系统包括:与所述加载控制器电性连接的:位于试验件底部并为试验件提供颤振力的激振器,布设在试验件上试验件控制点处的激光位移传感器和加速度传感器,布设在试验件上方的加热器,用于将控制信号放大的功率放大器,
用于将所述非定常气动力快速计算模型计算的等效气动力转化为激振器控制信号的半实物仿真子系统,所述半实物仿真子系统通过NI VeriStand仿真配置软件对接Matlab.simulink。
2.飞机地面振动试验中气动载荷模拟计算方法,基于权利要求1所述的飞机地面振动试验中气动载荷模拟计算系统,其特征在于,包括以下步骤:
S1、基于时间冻结假设,通过将多个离散时间点的非定常气动力快速计算模型进行Kriging代理模型拟合,得到适用于飞机地面振动试验整个时变历程的非定常气动力快速计算模型;
S2、通过实际飞机地面振动试验对步骤S1建立的非定常气动力快速计算模型进行试验验证,具体包括以下步骤:
S2-1、对飞机地面振动试验颤振计算中的气动平面和结构间的插值点进行缩聚,通过优化方法选取最优插值点,将最优插值点作为试验件控制点,
S2-2、选择数个采样时间点,将数个采样时间点数据导入步骤S1,得到适用于试验件控制点的非定常气动力快速计算模型,
S2-3、在实验室内搭建地面振动试验子系统,在试验件上方布设加热器,在试验件底部固定激振器以及与激振器电性连接的功率放大器,在试验件上试验件控制点布置激光位移传感器和加速度传感器,功率放大器、激光位移传感器、加速度传感器与加载控制器电性连接,
S2-4、通过激光位移传感器和加速度传感器将获取的试验件控制点响应信号传入加载控制器,加载控制器将试验件控制点响应信号代入步骤S2-1建立的非定常气动力快速计算模型,非定常气动力快速计算模型计算飞机地面振动试验时变历程中试验件的等效气动力,并将等效气动力传入半实物仿真子系统,半实物仿真子系统通过NI VeriStand仿真配置软件对接Matlab.simulink,将试验件的等效气动力转化为激振器控制信号,加载控制器通过激振器控制信号控制激振器将等效气动力实时加载在试验件上,加载完成后回到步骤S2-2,直至飞机振动试验等效模拟过程结束,在此过程中,加载控制器记录飞机振动试验等效模拟过程中的颤振临界速度,
S2-5、通过地面振动试验子系统开展实际飞机地面振动试验,并采用颤振边界追踪模块记录试验件在飞机地面颤振试验时变历程中的颤振临界速度,与步骤S2-4中非定常气动力快速计算模型计算得到的飞机地面振动试验时变历程中的颤振临界速度进行比对,以判定非定常气动力快速计算模型的准确性,若非定常气动力快速计算模型的准确性满足要求,得到试验验证后的非定常气动力快速计算模型,若不满足要求,加密采样时间点回到步骤S2-2,直至非定常气动力快速计算模型的准确性满足要求;
S3、采用步骤S2试验验证后的非定常气动力快速计算模型进行飞机地面振动试验中气动载荷模拟计算。
3.如权利要求2所述的飞机地面振动试验中气动载荷模拟计算方法,其特征在于,所述步骤S2-1中最优插值点为:激振点和拾振点。
4.如权利要求2所述的飞机地面振动试验中气动载荷模拟计算方法,其特征在于,所述步骤S1具体包括以下步骤:
S1-1、建立试验件控制点承受非定常气动力的计算公式,
S1-2、建立离散时间点的试验件控制点应承受非定常气动力计算模型,
S1-3、结合Kriging代理模型建立对应随机函数
S1-4、建立静态随机过程函数在任意两个采样时间点间的协方差计算模型
在飞机地面振动试验定义的时间范围内,得到的静态随机过程函数在任意两个采样时间点间的协方差计算模型如下:
上式中,为方差,Z(t1)为第一采样时间点的静态随机过程函数,Z(t2)为第二采样时间点的静态随机过程函数,t1为第一采样时间点,t2为第二采样时间点,R(t1,t2)为相关函数,相关函数的形式为高斯函数或三次样条函数,
S1-5、建立非定常气动力快速计算模型,即Kriging代理模型在任意时刻t的预测值计算模型
通过公式(4)构建形如式(3)的样本点数据列向量定义公式如下:
[Aij(t)]s=[Aij(t1) Aij(t2)…Aij(tn)]T=βG+z (5)
上式中,[Aij(t)]s为样本点数据列向量,β为基函数系数,[Aij(t1) Aij(t2)…Aij(tn)]为n个采样时间点的样本点数据,G为基函数列向量,G=[g(t1) g(t2)…g(tn)]T,z为随机函数列向量,z=[Z(t1) Z(t2)…Z(tn)]T,
为满足Kriging代理模型无偏估计E(Z(t))=0的要求,根据公式(5)与公式(8)可得基函数与基函数列向量关系的计算公式:
根据公式(3)、公式(8)与公式(9)得到Kriging代理模型的预测值误差计算公式如下:
上式中,为采样时间点外的时间节点函数值,为Aij(t)通过Kriging代理模型得到的对应随机函数,Aij(t)为任一时变气动力影响系数矩阵中的元素,Z为随机函数列向量,Z(t)是均值为0、方差为的静态随机过程函数,为加权行向量,
根据Kriging代理模型预测值满足无偏估计且均方误差最小的要求,即为在GTct=g(t)的前提下,寻找满足均方误差最小的加权列向量ct,为此引入拉格朗日乘数λ构建拉格朗日函数H(ct,λ)如下:
分别对公式(12)中加权列向量ct和拉格朗日乘数λ求偏导数,得到满足要求的加权列向量ct计算公式为:
将公式(13)带入公式(8)得到非定常气动力快速计算模型,即Kriging代理模型在任意时刻t的非定常气动力预测值计算模型如下:
5.如权利要求4所述的飞机地面振动试验中气动载荷模拟计算方法,其特征在于,所述步骤S1-1具体包括以下内容:
根据飞机地面振动试验中试验件的结构振动响应计算试验件控制点承受的非定常气动力,试验件控制点承受非定常气动力的计算公式如下:
fs=q∞A(k)zs (1)
上式中,fs为试验件控制点应承受的非定常气动力,q∞为来流动压,A(k)为降阶后的气动力影响系数矩阵,A(k)是通过与马赫数相匹配的非定常气动力理论计算后并依据结构动力学特性进行降阶处理得到的,Zs表示试验件的结构振动位移,k为减缩频率。
7.如权利要求5所述的飞机地面振动试验中气动载荷模拟计算方法,其特征在于,所述步骤S1-2具体包括以下内容:
由于实际飞机地面振动试验需要试验件控制点应承受非定常气动力的时域关系式,将公式(1)通过最小状态法有理函数近似延拓至拉氏域,再通过拉氏反变换得到离散时间点的试验件控制点应承受非定常气动力计算模型如下:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210770406.2A CN114970402A (zh) | 2022-06-30 | 2022-06-30 | 飞机地面振动试验中气动载荷模拟计算系统及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210770406.2A CN114970402A (zh) | 2022-06-30 | 2022-06-30 | 飞机地面振动试验中气动载荷模拟计算系统及方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114970402A true CN114970402A (zh) | 2022-08-30 |
Family
ID=82967399
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210770406.2A Pending CN114970402A (zh) | 2022-06-30 | 2022-06-30 | 飞机地面振动试验中气动载荷模拟计算系统及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114970402A (zh) |
-
2022
- 2022-06-30 CN CN202210770406.2A patent/CN114970402A/zh active Pending
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Silva et al. | Development of reduced-order models for aeroelastic analysis and flutter prediction using the CFL3Dv6. 0 code | |
Houtzager et al. | Recursive predictor-based subspace identification with application to the real-time closed-loop tracking of flutter | |
CN108121856B (zh) | 一种全飞行域飞行器动稳定性分析方法 | |
CN111368466B (zh) | 一种基于频响函数参数修正的机械振动预测方法 | |
CN112528411B (zh) | 一种基于模态减缩的几何非线性结构噪声振动响应计算方法 | |
CN115630558B (zh) | 一种复合材料构件装配变形预测方法 | |
CN113111547A (zh) | 基于缩减基的频域有限元模型修正方法 | |
Yang et al. | An improved nonlinear reduced-order modeling for transonic aeroelastic systems | |
CN112182954A (zh) | 一种基于lstm的流体仿真数据预测模型 | |
Li et al. | An efficient implementation of aeroelastic tailoring based on efficient computational fluid dynamics-based reduced order model | |
CN114583767B (zh) | 一种数据驱动的风电场调频响应特性建模方法及系统 | |
CN114818550B (zh) | 一种飞机振动试验中时变气动载荷地面等效模拟方法 | |
CN110455477B (zh) | 一种固体火箭舱段结构振动载荷谱的获取方法 | |
CN113722860B (zh) | 基于降阶模型的瞬态热力状态在线评估方法、设备及介质 | |
CN110765560B (zh) | 一种基于时变阻尼机械机构振动预测方法 | |
Vendl et al. | Projection-based model order reduction for steady aerodynamics | |
CN112580855A (zh) | 基于自适应变异pso-bp神经网络的电缆群稳态温升预测方法 | |
CN114970402A (zh) | 飞机地面振动试验中气动载荷模拟计算系统及方法 | |
CN111539138A (zh) | 基于阶跃函数的结构动力学峰值时域响应灵敏度求解方法 | |
CN113221064B (zh) | 模型处理方法、装置、计算机设备和存储介质 | |
Song et al. | Nonlinear reduced roms: Formulation and applications | |
CN111639463B (zh) | 一种基于XGBoost算法的电力系统扰动后频率特征预测方法 | |
CN114491730A (zh) | 一种高速铁路路基结构动力安定分析迭代方法及装置 | |
CN104038132A (zh) | 一种具有时变测量延迟输出和噪声的伺服电机的状态观测方法 | |
Yu | An Accurate three-dimensional deformation measurement method in wind turbine Blade Static Loading Test |
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 |