CN108182328A - A kind of big angle of attack Nonlinear Aerodynamic reduced-order model suitable for stall flutter - Google Patents

A kind of big angle of attack Nonlinear Aerodynamic reduced-order model suitable for stall flutter Download PDF

Info

Publication number
CN108182328A
CN108182328A CN201810011636.4A CN201810011636A CN108182328A CN 108182328 A CN108182328 A CN 108182328A CN 201810011636 A CN201810011636 A CN 201810011636A CN 108182328 A CN108182328 A CN 108182328A
Authority
CN
China
Prior art keywords
model
attack
aerodynamic
nonlinear
nonlinear aerodynamic
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
CN201810011636.4A
Other languages
Chinese (zh)
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.)
Beihang University
Original Assignee
Beihang University
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 Beihang University filed Critical Beihang University
Priority to CN201810011636.4A priority Critical patent/CN108182328A/en
Publication of CN108182328A publication Critical patent/CN108182328A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/15Vehicle, aircraft or watercraft design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/04Architecture, e.g. interconnection topology
    • G06N3/045Combinations of networks
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/08Learning methods
    • G06N3/084Backpropagation, e.g. using gradient descent
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling 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)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Biophysics (AREA)
  • Molecular Biology (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Biomedical Technology (AREA)
  • Computer Hardware Design (AREA)
  • Computational Linguistics (AREA)
  • Data Mining & Analysis (AREA)
  • General Health & Medical Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Computing Systems (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • Automation & Control Theory (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

A kind of big angle of attack Nonlinear Aerodynamic reduced-order model suitable for stall flutter of the invention designs sinusoidal multilevel signal as Nonlinear Aerodynamic identification signal first, and the amplitude of modelled signal need to meet substantially movement needs.Nonlinear Aerodynamic identification is carried out secondly based on the Recognition with Recurrent Neural Network model of deep learning, training output signal is calculated using CFD approach, dynamic Time-Domain Nonlinear Aerodynamic Model when substantially being vibrated with the Recognition with Recurrent Neural Network model foundation wing based on deep learning, and the system calculated with CFD under test signal responds, compared with the identification result of network model, the performance of model is verified.The present invention saves Flight Vehicle Design cost, improves stall flutter design efficiency, has developed the reduced-order model of Nonlinear Aerodynamic, so as to fast prediction stall flutter.

Description

一种适用于失速颤振的大攻角非线性气动力降阶模型A Nonlinear Aerodynamic Reduced-Order Model for Large Angle of Attack Applicable to Stall Flutter

技术领域technical field

本发明属于航空飞行器设计与系统辨识领域,是一种针对快速预测飞行器失速颤振特性,提高计算效率而发展出的非线性降阶模型。The invention belongs to the field of aircraft design and system identification, and is a non-linear reduced-order model developed for rapidly predicting aircraft stall flutter characteristics and improving calculation efficiency.

背景技术Background technique

失速颤振是飞行器翼面或舵面处于大攻角时,气流分离引起的非线性气动力和弹性结构耦合所发生的自激振动,失速颤振体现出强非线性特性。飞行器大攻角飞行(比如歼击机或者导弹具有高机动性和敏捷性,需要大迎角范围内飞行)或遭遇阵风时,可能会使升力面处于失速状态,当迎角达到某一临界值时容易导致气动失速,进而可能发生失速颤振,影响飞行器的结构安全;尤其是螺旋桨叶和直升机叶片,桨尖较易发生失速及失速颤振的现象。Stall flutter is a self-excited vibration caused by nonlinear aerodynamic force and elastic structure coupling caused by airflow separation when the aircraft wing or rudder surface is at a large angle of attack. Stall flutter exhibits strong nonlinear characteristics. When an aircraft flies at a high angle of attack (for example, a fighter or missile has high maneuverability and agility and needs to fly within a range of a large angle of attack) or encounters a gust of wind, it may cause the lifting surface to stall. When the angle of attack reaches a certain critical value, it is easy to This will lead to aerodynamic stall, which may lead to stall flutter, which will affect the structural safety of the aircraft; especially for propeller blades and helicopter blades, the tip of the blade is more prone to stall and stall flutter.

目前对机翼失速颤振的研究方法有试验研究、经验模型和CFD-CSD流固耦合仿真。试验方法耗时耗力,成本较高。经验模型计算依赖模型和试验数据,结果精度不高。CFD-CSD计算能获得高精度的结果,但是计算量巨大,耗时很长,增加计算成本和飞机研制周期。At present, the research methods of wing stall flutter include experimental research, empirical model and CFD-CSD fluid-structure interaction simulation. The test method is time-consuming and labor-intensive, and the cost is high. The calculation of the empirical model depends on the model and test data, and the accuracy of the result is not high. CFD-CSD calculation can obtain high-precision results, but the calculation is huge and time-consuming, which increases the calculation cost and aircraft development cycle.

要实现飞行器失速颤振的预测,首先要解决的关键技术问题是对大攻角非线性气动力的准确辨识,主要存在以下问题:To realize the prediction of aircraft stall flutter, the key technical problem to be solved is the accurate identification of nonlinear aerodynamic forces at large angles of attack. The main problems are as follows:

1.传统的信号在非设计点的非线性气动力预测效果较差。要准确捕捉大攻角极限环振动的特性,需合理设计合理的辨识输入信号,设计的辨识输入信号需准确刻画极限环振动的规律、大运动幅值和频率要求。1. Traditional signals are less effective in predicting non-linear aerodynamics at non-design points. In order to accurately capture the vibration characteristics of the limit cycle with large angle of attack, it is necessary to design a reasonable identification input signal. The designed identification input signal must accurately describe the law of limit cycle vibration, large motion amplitude and frequency requirements.

2.对于设计好的辨识输入信号,需要发展鲁棒的泛化的降阶模型和算法来准确辨识动态的大攻角分离的非线性气动力的数学模型。2. For a well-designed identification input signal, it is necessary to develop a robust generalized reduced-order model and an algorithm to accurately identify the dynamic mathematical model of nonlinear aerodynamic forces separated by large angles of attack.

基于上述情况,需要在输入信号设计和辨识算法上提出新的模型。Based on the above situation, it is necessary to propose a new model in the input signal design and identification algorithm.

发明内容Contents of the invention

针对上述问题,为了节省飞行器设计成本,提高失速颤振设计效率,提出了一种适用于失速颤振的大攻角非线性气动力降阶模型,从而快速预测失速颤振。In view of the above problems, in order to save the cost of aircraft design and improve the efficiency of stall flutter design, a nonlinear aerodynamic reduction model with large angle of attack is proposed for stall flutter, so as to quickly predict stall flutter.

本发明适用于失速颤振的大攻角非线性气动力降阶模型,通过下述步骤得到:The present invention is applicable to the large angle of attack nonlinear aerodynamic reduction model of stall flutter, obtained through the following steps:

步骤1:根据机翼极限环振动时的幅值、频率和振动规律,设计左右对称的多级正弦信号,作为深度学习模型系统的输入数据。Step 1: According to the amplitude, frequency and vibration law of the wing limit cycle vibration, design a left-right symmetrical multi-level sinusoidal signal as the input data of the deep learning model system.

步骤2:将步骤1种设计好的正弦信号输入到CFD软件中,获得不同速度下的机翼在该信号激励下的气动力系数,作为深度学习模型系统的输出数据。Step 2: Input the sinusoidal signal designed in step 1 into the CFD software, and obtain the aerodynamic coefficients of the wing at different speeds under the excitation of the signal, as the output data of the deep learning model system.

步骤3:基于深度学习的循环神经网络模型进行非线性气动力辨识,得到大攻角非线性气动力降阶模型。Step 3: Perform nonlinear aerodynamic identification based on the cyclic neural network model of deep learning, and obtain the nonlinear aerodynamic reduction model for large angle of attack.

将上述的大攻角非线性气动力降阶模型带入到结构和气动力耦合计算方程中,利用4阶龙格-库塔方法进行时域推进计算,预测各模态随时间的响应过程,从而达到预测失速颤振的目的。The above-mentioned nonlinear aerodynamic reduction model at large angle of attack is brought into the structural and aerodynamic coupling calculation equation, and the 4th-order Runge-Kutta method is used to perform time-domain propulsion calculations to predict the response process of each mode over time, so that To achieve the purpose of predicting stall flutter.

本发明的优点在于:The advantages of the present invention are:

1、本发明适用于失速颤振的大攻角非线性气动力降阶模型,通过对失速颤振振动特点的分析,设计了多级正弦辨识信号,这种信号覆盖了振动的幅值、频率范围,能够很好的预测振动过程的非线性气动力。1. The present invention is applicable to the large angle of attack nonlinear aerodynamic reduction model of stall flutter. Through the analysis of the characteristics of stall flutter vibration, a multi-level sine identification signal is designed. This signal covers the amplitude and frequency of vibration range, it can well predict the nonlinear aerodynamic forces of the vibration process.

2、本发明适用于失速颤振的大攻角非线性气动力降阶模型,是基于深度学习的循环神经网络模型的非线性气动力辨识模型,在动态非线性流场建模方面具有巨大的优势和应用前景。2. The present invention is applicable to a large angle of attack nonlinear aerodynamic reduction model for stall flutter, and is a nonlinear aerodynamic identification model based on a deep-learning cyclic neural network model, which has huge advantages in dynamic nonlinear flow field modeling Advantages and application prospects.

3、本发明适用于失速颤振的大攻角非线性气动力降阶模型,适应强,辨识精度高,操作方法比较简单,利用MATLAB编程即可实现。3. The present invention is suitable for a large angle of attack nonlinear aerodynamic reduction model of stall flutter, and has strong adaptability, high identification accuracy, relatively simple operation method, and can be realized by using MATLAB programming.

4、本发明适用于失速颤振的大攻角非线性气动力降阶模型,随着CFD技术的发展,省去了复杂的实验研究,同时计算效率又可比CFD计算提高一到两个数量级,辨识方法亦可用在其他非线性气动力辨识系统中,具有一定的通用性。4. The present invention is applicable to the large angle of attack nonlinear aerodynamic reduction model of stall flutter. With the development of CFD technology, complicated experimental research is omitted, and the calculation efficiency can be increased by one to two orders of magnitude compared with CFD calculation. The identification method can also be used in other nonlinear aerodynamic identification systems, and has certain versatility.

附图说明Description of drawings

图1为本发明适用于失速颤振的大攻角非线性气动力降阶模型设计流程图;Fig. 1 is the flow chart of the design of the large angle of attack nonlinear aerodynamic order reduction model applicable to stall flutter of the present invention;

图2为本发明适用于失速颤振的大攻角非线性气动力降阶模型的辨识信号示意图;Fig. 2 is a schematic diagram of the identification signal of the large angle of attack nonlinear aerodynamic reduction model applicable to stall flutter in the present invention;

图3为本发明适用于失速颤振的大攻角非线性气动力降阶模型中采用的基于深度学习的循环神经网络模型结构图;Fig. 3 is the structural diagram of the cyclic neural network model based on deep learning adopted in the large angle of attack nonlinear aerodynamic reduction model applicable to stall flutter in the present invention;

图4为基于深度学习的循环神经网络模型中神经元单元展开计算结构图。Fig. 4 is a structural diagram of neuron unit expansion calculation in the deep learning-based cyclic neural network model.

具体实施方式Detailed ways

下面结合附图对本发明作进一步详细说明。The present invention will be described in further detail below in conjunction with the accompanying drawings.

本发明适用于失速颤振的大攻角非线性气动力降阶模型,如图1所示,通过下述步骤得到:The present invention is applicable to the large angle of attack nonlinear aerodynamic reduction model of stall flutter, as shown in Figure 1, obtained through the following steps:

步骤1:设计输入信号。Step 1: Design the input signal.

输入信号是为预测机翼失速颤振而设计的,机翼发生大攻角极限环振动时,每一个周期的振动形式为正弦振动。本发明设计方法为,针对一种机翼模型,根据其极限环振动时的幅值、频率和振动规律,设计左右对称的多级正弦信号,这些正弦包含了不同的幅值和频率信息,且覆盖了所需的幅值和频率信息;其中,频率一般为低频振动,因此信号频率范围覆盖失速颤振振动频率的范围即可。如图2所示。在针对具体的模型时,需要根据振动幅值和频率的需要来设计。The input signal is designed to predict the stall flutter of the wing. When the wing has a large angle of attack limit cycle vibration, the vibration form of each cycle is sinusoidal. The design method of the present invention is, for a wing model, according to the amplitude, frequency and vibration law of its limit cycle vibration, design left and right symmetrical multi-level sinusoidal signals, these sinusoids include different amplitude and frequency information, and The required amplitude and frequency information are covered; the frequency is generally low-frequency vibration, so the signal frequency range only needs to cover the range of the stall flutter vibration frequency. as shown in picture 2. When targeting a specific model, it needs to be designed according to the needs of vibration amplitude and frequency.

步骤2:开发CFD求解器,计算气动力系数。Step 2: Develop a CFD solver to calculate the aerodynamic coefficients.

将步骤1种设计好的正弦信号输入到CFD软件中,考虑到大攻角条件下,不同的来流速度对计算出的气动力参数会有影响,因此通过CFD软件获得不同速度下的机翼在该信号激励下的气动力系数,包括升力系数、力矩系数、阻力系数等。上述气动力系数为深度学习模型系统的输出数据,而设计的正弦信号为深度学习模型系统的输入数据。Input the sinusoidal signal designed in step 1 into the CFD software. Considering that under the condition of large angle of attack, different incoming flow velocities will have an impact on the calculated aerodynamic parameters, so the wings at different speeds are obtained through the CFD software The aerodynamic coefficient under the signal excitation includes lift coefficient, moment coefficient, drag coefficient, etc. The above-mentioned aerodynamic coefficient is the output data of the deep learning model system, and the designed sinusoidal signal is the input data of the deep learning model system.

步骤3:非线性气动力辨识。Step 3: Nonlinear aerodynamic force identification.

本发明中,非线性气动力辨识基于深度学习的循环神经网络模型,即步骤2中的深度学习模型,如图3所示。由于非线性流场的复杂性,传统模型对经验要求很高,很难给出一个合适的显式表达式,而基于深度学习的神经网络(Neural Network,NN)模型具有不需要给出辨识系统输入/输出之间的显示数学表达式的优点。该神经网络模型通过学习训练型号中输入/输出的特征,以及输入影响输出的方式,从而获取带辨识系统的特征,并用学习的“经验”预测新的输入下的输出。In the present invention, the nonlinear aerodynamic identification is based on the deep learning cycle neural network model, that is, the deep learning model in step 2, as shown in FIG. 3 . Due to the complexity of the nonlinear flow field, the traditional model has high requirements for experience, and it is difficult to give a suitable explicit expression, while the neural network (Neural Network, NN) model based on deep learning does not need to give an identification system Advantages of displaying math expressions between input/output. The neural network model obtains the characteristics of the identification system by learning the characteristics of the input/output in the training model and the way the input affects the output, and uses the learned "experience" to predict the output under the new input.

基于深度学习的循环神经网络模型分为输入层、隐含层和输出层。其中,隐含层数大于4层,具体层数可以根据实际模型进行调整。从输入层到输出层,每层神经元单元和前后层的神经元单元连接并传输数据。隐含层的每一个神经元单元会反馈回来作为自身的输入。The recurrent neural network model based on deep learning is divided into input layer, hidden layer and output layer. Among them, the number of hidden layers is greater than 4 layers, and the specific number of layers can be adjusted according to the actual model. From the input layer to the output layer, the neuron units of each layer and the neuron units of the previous and subsequent layers are connected and transmit data. Each neuron unit in the hidden layer feeds back as its own input.

上述每个的神经元单元的内部结构展开图如图4所示。An expanded view of the internal structure of each of the above neuron units is shown in FIG. 4 .

神经元单元是按时间序列展开,对该展开图做以下说明:The neuron unit is expanded in time series, and the expansion diagram is described as follows:

1)x(t)为t时刻训练样本的输入,同样的,x(t-1)和x(t+1)分别为t-1时刻和t+1时刻训练样本的输入。1) x (t) is the input of training samples at time t, and similarly, x (t-1) and x (t+1) are the inputs of training samples at time t-1 and time t+1 respectively.

2)h(t)为t时刻模型的隐藏状态,h(t)由x(t)和h(t-1)共同决定。2) h (t) is the hidden state of the model at time t, and h (t) is jointly determined by x (t) and h (t-1) .

3)o(t)代表t时刻模型的输出,o(t)只由模型当前的隐藏状态h(t)决定。3) o (t) represents the output of the model at time t, and o (t) is only determined by the current hidden state h (t) of the model.

4)L(t)代表t时刻模型的损失函数。4) L (t) represents the loss function of the model at time t.

5)y(t)代表t时刻训练样本的真实输出。5) y (t) represents the real output of the training sample at time t.

6)U,W,V这三个矩阵是模型的线性关系参数,它在整个循环网络中是共享的,体现了循环网络模型的“循环反馈”的思想。6) The three matrices U, W, and V are the linear relationship parameters of the model, which are shared in the entire recurrent network, reflecting the idea of "circular feedback" in the recurrent network model.

应用上述基于深度学习的循环神经网络模型,可以逼近任意的动态非线性映射系统,可很好的辨识动态非线性非定常气动力模型。整个循环神经网络的前向传播算法和普通前馈神经网络相同,其中隐含层激活函数选用连续双曲正切函数,输出层激活函数选用softmax函数。循环神经网络参数训练时选用误差反向传播算法中的梯度下降法,通过反复调整循环神经网络内部的连接权值和阀值,最终使网络输出逼近实际的输出。Applying the above-mentioned deep learning-based cyclic neural network model can approximate any dynamic nonlinear mapping system, and can well identify dynamic nonlinear unsteady aerodynamic models. The forward propagation algorithm of the entire cyclic neural network is the same as that of the ordinary feed-forward neural network. The activation function of the hidden layer is selected from the continuous hyperbolic tangent function, and the activation function of the output layer is selected from the softmax function. The gradient descent method in the error backpropagation algorithm is used in the training of the cyclic neural network parameters, and the network output is finally approached to the actual output by repeatedly adjusting the connection weights and thresholds inside the cyclic neural network.

步骤4:利用训练好的网络模型,预测机翼的失速颤振。将网络模型带入到结构和气动力耦合计算方程中,利用4阶龙格-库塔方法进行时域推进计算,预测各模态随时间的响应过程,从而达到预测失速颤振的目的。Step 4: Use the trained network model to predict the stall flutter of the wing. The network model is brought into the structural and aerodynamic coupling calculation equation, and the 4th-order Runge-Kutta method is used for time-domain propulsion calculation to predict the response process of each mode over time, so as to achieve the purpose of predicting stall flutter.

Claims (5)

1.一种适用于失速颤振的大攻角非线性气动力降阶模型,其特征在于:通过下述步骤得到:1. A large angle of attack nonlinear aerodynamic reduction model applicable to stall flutter, characterized in that: obtain by the following steps: 步骤1:根据机翼极限环振动时的幅值、频率和振动规律,设计左右对称的多级正弦信号,作为深度学习模型系统的输入数据;Step 1: According to the amplitude, frequency and vibration law of the wing limit cycle vibration, design a left-right symmetrical multi-level sinusoidal signal as the input data of the deep learning model system; 步骤2:将步骤1种设计好的正弦信号输入到CFD软件中,获得不同速度下的机翼在该信号激励下的气动力系数,作为深度学习模型系统的输出数据;Step 2: Input the sinusoidal signal designed in step 1 into the CFD software, and obtain the aerodynamic coefficient of the wing at different speeds under the excitation of the signal, as the output data of the deep learning model system; 步骤3:基于深度学习的循环神经网络模型进行非线性气动力辨识,得到大攻角非线性气动力降阶模型。Step 3: Perform nonlinear aerodynamic identification based on the cyclic neural network model of deep learning, and obtain the nonlinear aerodynamic reduction model for large angle of attack. 2.如权利要求1所述一种适用于失速颤振的大攻角非线性气动力降阶模型,其特征在于:步骤1里多级正弦信号包含不同的幅值和频率信息。2. A large angle-of-attack nonlinear aerodynamic reduced-order model suitable for stall flutter as claimed in claim 1, characterized in that: in step 1, the multi-level sinusoidal signals contain different amplitude and frequency information. 3.如权利要求1所述一种适用于失速颤振的大攻角非线性气动力降阶模型,其特征在于:基于深度学习的循环神经网络模型中,隐含层激活函数选用连续双曲正切函数,输出层激活函数选用softmax函数。3. A kind of large angle of attack nonlinear aerodynamic reduction model applicable to stall flutter as claimed in claim 1, characterized in that: in the cyclic neural network model based on deep learning, the hidden layer activation function selects continuous hyperbolic The tangent function, the output layer activation function selects the softmax function. 4.如权利要求1所述一种适用于失速颤振的大攻角非线性气动力降阶模型,其特征在于:基于深度学习的循环神经网络模型进行参数训练时选用误差反向传播算法中的梯度下降法。4. A kind of large angle-of-attack nonlinear aerodynamic reduction model applicable to stall flutter as claimed in claim 1, is characterized in that: when the cyclic neural network model based on deep learning is used for parameter training, the error backpropagation algorithm is selected The gradient descent method. 5.如权利要求1所述一种适用于失速颤振的大攻角非线性气动力降阶模型,其特征在于:将大攻角非线性气动力降阶模型带入到结构和气动力耦合计算方程中,利用4阶龙格-库塔方法进行时域推进计算,预测各模态随时间的响应过程,从而达到预测失速颤振的目的。5. A large angle of attack nonlinear aerodynamic reduced order model suitable for stall flutter as claimed in claim 1, characterized in that: the large angle of attack nonlinear aerodynamic reduced order model is brought into the structure and aerodynamic coupling calculation In the equation, the 4th-order Runge-Kutta method is used to carry out time-domain advance calculations to predict the response process of each mode over time, so as to achieve the purpose of predicting stall flutter.
CN201810011636.4A 2018-01-05 2018-01-05 A kind of big angle of attack Nonlinear Aerodynamic reduced-order model suitable for stall flutter Pending CN108182328A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810011636.4A CN108182328A (en) 2018-01-05 2018-01-05 A kind of big angle of attack Nonlinear Aerodynamic reduced-order model suitable for stall flutter

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810011636.4A CN108182328A (en) 2018-01-05 2018-01-05 A kind of big angle of attack Nonlinear Aerodynamic reduced-order model suitable for stall flutter

Publications (1)

Publication Number Publication Date
CN108182328A true CN108182328A (en) 2018-06-19

Family

ID=62550012

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810011636.4A Pending CN108182328A (en) 2018-01-05 2018-01-05 A kind of big angle of attack Nonlinear Aerodynamic reduced-order model suitable for stall flutter

Country Status (1)

Country Link
CN (1) CN108182328A (en)

Cited By (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109063290A (en) * 2018-07-20 2018-12-21 中国航空工业集团公司沈阳飞机设计研究所 A kind of flutter prediction technique based on nerual network technique
CN109086501A (en) * 2018-07-20 2018-12-25 中国航空工业集团公司沈阳飞机设计研究所 A kind of flutter prediction technique
CN110826600A (en) * 2019-10-18 2020-02-21 北京航空航天大学 Engine surge prediction method based on adaptive resonance network online incremental learning
CN111898327A (en) * 2020-06-30 2020-11-06 西北工业大学 An abnormal data expansion method for flutter signal of aeroelastic system
CN112380794A (en) * 2020-12-08 2021-02-19 中北大学 Multi-disciplinary parallel cooperation optimization design method for aviation turbine engine blade
CN112800543A (en) * 2021-01-27 2021-05-14 中国空气动力研究与发展中心计算空气动力研究所 Nonlinear unsteady aerodynamic modeling method based on improved Goman model
CN112948973A (en) * 2021-03-04 2021-06-11 北京航空航天大学 Wing stall flutter closed-loop control method for continuously variable camber trailing edge
CN113673031A (en) * 2021-08-11 2021-11-19 中国科学院力学研究所 A flexible airship service angle-of-attack identification method based on the fusion of strain response and deep learning
CN115034152A (en) * 2022-05-17 2022-09-09 浙江大学 A data-driven nonlinear reduced-order prediction method and device for fluid-structure coupled systems
CN115343012A (en) * 2022-07-07 2022-11-15 中国航空工业集团公司哈尔滨空气动力研究所 Unsteady-state large-amplitude oscillation test method
CN115408931A (en) * 2022-08-16 2022-11-29 哈尔滨工业大学 Vortex vibration response prediction method based on deep learning
CN115422497A (en) * 2022-08-16 2022-12-02 哈尔滨工业大学 Ordinary Differential Equation Recognition Method Based on Convolutional Differential Operator and Symbolic Network
CN117992761A (en) * 2024-04-07 2024-05-07 西北工业大学 Intelligent prediction method for aerodynamic force of dynamic stall of wind turbine blade
CN118133429A (en) * 2024-04-15 2024-06-04 北京航空航天大学 Flutter suppression design and verification method based on reinforcement learning and variable camber trailing edge
CN118862654A (en) * 2024-07-03 2024-10-29 北京航空航天大学 An efficient and high-precision prediction method for stall flutter based on energy graph and deep neural operator network

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6606612B1 (en) * 1998-08-13 2003-08-12 The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration Method for constructing composite response surfaces by combining neural networks with other interpolation or estimation techniques
CN107525647A (en) * 2017-07-26 2017-12-29 北京航空航天大学 A kind of dynamical bifurcation generating means of aerodynamic stalling

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6606612B1 (en) * 1998-08-13 2003-08-12 The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration Method for constructing composite response surfaces by combining neural networks with other interpolation or estimation techniques
CN107525647A (en) * 2017-07-26 2017-12-29 北京航空航天大学 A kind of dynamical bifurcation generating means of aerodynamic stalling

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
IAN GOODFELLOW等: "《深度学习》", 31 August 2017, 人民邮电出版社 *
THI-THU-HUONG LE等: "Analyzing Effective of Activation Functions on Recurrent Neural", 《JOURNAL OF MULTIMEDIA AND INFORMATION SYSTEM》 *
XIANG ZHENGPING等: "Nonlinear Aerodynamic Identification and Stall Flutter Analysis", 《INTERNATIONAL FORUM ON AEROELASTICITY AND STRUCTURAL DYNAMICS》 *
刘晓燕等: "适用于气动弹性的气动力自适应辨识方法", 《工程力学》 *
叶正寅等: "《流固耦合力学基础及其应用》", 30 June 2010, 哈尔滨工业大学出版社 *
杨超等: "考虑机动载荷不确定性的机翼气动弹性优化设计", 《中国科学:技术科学》 *

Cited By (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109063290A (en) * 2018-07-20 2018-12-21 中国航空工业集团公司沈阳飞机设计研究所 A kind of flutter prediction technique based on nerual network technique
CN109086501A (en) * 2018-07-20 2018-12-25 中国航空工业集团公司沈阳飞机设计研究所 A kind of flutter prediction technique
CN110826600A (en) * 2019-10-18 2020-02-21 北京航空航天大学 Engine surge prediction method based on adaptive resonance network online incremental learning
CN111898327B (en) * 2020-06-30 2022-06-07 西北工业大学 Flutter signal abnormal data expansion method for aeroelastic system
CN111898327A (en) * 2020-06-30 2020-11-06 西北工业大学 An abnormal data expansion method for flutter signal of aeroelastic system
CN112380794A (en) * 2020-12-08 2021-02-19 中北大学 Multi-disciplinary parallel cooperation optimization design method for aviation turbine engine blade
CN112380794B (en) * 2020-12-08 2022-11-08 中北大学 A Multidisciplinary Parallel Collaborative Optimal Design Method for Aero Turbine Engine Blades
CN112800543A (en) * 2021-01-27 2021-05-14 中国空气动力研究与发展中心计算空气动力研究所 Nonlinear unsteady aerodynamic modeling method based on improved Goman model
CN112800543B (en) * 2021-01-27 2022-09-13 中国空气动力研究与发展中心计算空气动力研究所 Nonlinear unsteady aerodynamic modeling method based on improved Goman model
CN112948973A (en) * 2021-03-04 2021-06-11 北京航空航天大学 Wing stall flutter closed-loop control method for continuously variable camber trailing edge
CN112948973B (en) * 2021-03-04 2022-05-03 北京航空航天大学 Wing stall flutter closed-loop control method for continuously variable camber trailing edge
CN113673031A (en) * 2021-08-11 2021-11-19 中国科学院力学研究所 A flexible airship service angle-of-attack identification method based on the fusion of strain response and deep learning
CN113673031B (en) * 2021-08-11 2024-04-12 中国科学院力学研究所 Flexible airship service attack angle identification method integrating strain response and deep learning
CN115034152A (en) * 2022-05-17 2022-09-09 浙江大学 A data-driven nonlinear reduced-order prediction method and device for fluid-structure coupled systems
CN115343012B (en) * 2022-07-07 2023-04-07 中国航空工业集团公司哈尔滨空气动力研究所 Unsteady-state large-amplitude oscillation test method
CN115343012A (en) * 2022-07-07 2022-11-15 中国航空工业集团公司哈尔滨空气动力研究所 Unsteady-state large-amplitude oscillation test method
CN115408931A (en) * 2022-08-16 2022-11-29 哈尔滨工业大学 Vortex vibration response prediction method based on deep learning
CN115422497A (en) * 2022-08-16 2022-12-02 哈尔滨工业大学 Ordinary Differential Equation Recognition Method Based on Convolutional Differential Operator and Symbolic Network
CN117992761A (en) * 2024-04-07 2024-05-07 西北工业大学 Intelligent prediction method for aerodynamic force of dynamic stall of wind turbine blade
CN117992761B (en) * 2024-04-07 2024-06-11 西北工业大学 An intelligent prediction method for dynamic stall aerodynamics of wind turbine blades
CN118133429A (en) * 2024-04-15 2024-06-04 北京航空航天大学 Flutter suppression design and verification method based on reinforcement learning and variable camber trailing edge
CN118862654A (en) * 2024-07-03 2024-10-29 北京航空航天大学 An efficient and high-precision prediction method for stall flutter based on energy graph and deep neural operator network
CN118862654B (en) * 2024-07-03 2025-01-24 北京航空航天大学 An efficient and high-precision prediction method for stall flutter based on energy graph and deep neural operator network

Similar Documents

Publication Publication Date Title
CN108182328A (en) A kind of big angle of attack Nonlinear Aerodynamic reduced-order model suitable for stall flutter
Xia et al. Unsteady aerodynamics and vortex-sheet formation of a two-dimensional airfoil
Brunton et al. Reduced-order unsteady aerodynamic models at low Reynolds numbers
Antoniadis et al. Assessment of CFD methods against experimental flow measurements for helicopter flows
CN104298805A (en) CFD pneumatic modeling method for hypersonic aerocraft
CN113127976B (en) A kind of boundary layer separation induced transition prediction method, device, equipment and medium
Glaz et al. A surrogate-based approach to reduced-order dynamic stall modeling
CN113221237B (en) A high angle of attack flutter analysis method based on reduced-order modeling
CN114004023A (en) A method for identifying aerodynamic parameters of aircraft based on cyclic neural network
Liu et al. Efficient reduced-order aerodynamic modeling in low-Reynolds-number incompressible flows
Cayon et al. Fast Aeroelastic Model of a Leading-Edge Inflatable Kite
Reade et al. Deformable blade element and unsteady vortex lattice fluid-structure interaction modeling of a 2d flapping wing
CN117556727A (en) Flapping time-varying aerodynamic force prediction method and device based on deep learning
Manar et al. Transient response of a single degree-of-freedom wing at high angle-of-attack
CN112948973B (en) Wing stall flutter closed-loop control method for continuously variable camber trailing edge
Şugar-Gabor et al. Fast and accurate quasi-3D aerodynamic methods for aircraft conceptual design studies
Miao et al. The Aerodynamic Characteristics of a Diamond Joined‐Wing Morphing Aircraft
Paturle et al. Dynamic detection of flow separation using integral formulation of unsteady boundary layer equations
Gardner et al. Rotating wing dynamic stall: State of the art and future directions
Wu et al. Influence of area distribution along the span direction on flapping wing aerodynamics in hover based on numerical modeling analysis
Shen et al. Aerodynamic characteristics of multiple flapping wing configurations
Abdol-Hamid et al. DPW-6 Results Using FUN3D With Focus on k-kL-MEAH2015 Turbulence Model
Zhu et al. The aerodynamic performance of flexible wing in plunge
Lankford Experimental and coupled CFD/CSD investigation of flexible MAV-scale flapping wings in hover
Melville et al. An efficient model for aeroelastic tailoring of aircraft wings under gust loads

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