CN109726440B - Aeroelasticity analysis method considering dynamic characteristics of internal fluid - Google Patents
Aeroelasticity analysis method considering dynamic characteristics of internal fluid Download PDFInfo
- Publication number
- CN109726440B CN109726440B CN201811476560.9A CN201811476560A CN109726440B CN 109726440 B CN109726440 B CN 109726440B CN 201811476560 A CN201811476560 A CN 201811476560A CN 109726440 B CN109726440 B CN 109726440B
- Authority
- CN
- China
- Prior art keywords
- fluid
- internal fluid
- matrix
- considering
- equation
- 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
- 239000012530 fluid Substances 0.000 title claims abstract description 58
- 238000004458 analytical method Methods 0.000 title claims abstract description 24
- 239000011159 matrix material Substances 0.000 claims abstract description 30
- 238000000034 method Methods 0.000 claims abstract description 10
- 239000007787 solid Substances 0.000 claims description 15
- 238000006073 displacement reaction Methods 0.000 claims description 8
- 230000008878 coupling Effects 0.000 claims description 5
- 238000010168 coupling process Methods 0.000 claims description 5
- 238000005859 coupling reaction Methods 0.000 claims description 5
- 230000035515 penetration Effects 0.000 claims description 2
- 238000009795 derivation Methods 0.000 claims 1
- 230000017105 transposition Effects 0.000 claims 1
- 230000009286 beneficial effect Effects 0.000 abstract description 3
- 238000005094 computer simulation Methods 0.000 abstract description 2
- 238000005516 engineering process Methods 0.000 abstract description 2
- 230000000694 effects Effects 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 4
- 239000000446 fuel Substances 0.000 description 2
- 238000010586 diagram Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
Images
Classifications
-
- 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
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
The invention belongs to the aeroelasticity analysis technology of airplanes, and particularly relates to an aeroelasticity analysis method considering the dynamic characteristics of internal fluid; the method comprises the steps of firstly generating an additional mass matrix considering internal fluid, then assembling a motion equation considering the additional mass matrix, and then carrying out aeroelastic analysis considering the influence of the internal fluid. The method provided by the invention can consider the influence of the dynamic characteristics of the internal fluid in the dynamic simulation model of the airplane with a large amount of internal fluid, thereby being beneficial to improving the aeroelasticity analysis precision of the whole airplane and better guiding the aeroelasticity design. The proposed analysis method is simple and convenient to operate and has high engineering application value.
Description
Technical Field
The invention belongs to the technology of aeroelasticity analysis of airplanes, and particularly relates to an aeroelasticity analysis method considering dynamic characteristics of internal fluid.
Background
In order to improve the combat efficiency, modern aircraft, especially large-scale fuel dispensers and fighters, carry a large amount of fuel. For the aeroelastic characteristic evaluation of the aircrafts, the traditional calculation method takes the internal fluid as a concentrated mass, only considers the mass and inertia characteristics, and does not consider the dynamic characteristics of the internal fluid and the influence between the internal fluid and a container. The problem of structural vibration caused by such internal fluid-to-structure coupling, and the effect of fluid elasticity on the aeroelastic properties of the aircraft, is not clear. The influence of the factors is rarely considered in aeroelastic analysis at home at present.
Disclosure of Invention
The purpose of the invention is as follows:
an aeroelastic analysis method is presented that takes into account the effects of internal fluid dynamics.
The technical scheme of the invention is as follows: .
A method of aeroelastic analysis that takes into account the dynamic characteristics of the internal fluid; the method comprises the steps of firstly generating an additional mass matrix considering internal fluid, then assembling a motion equation considering the additional mass matrix, and then carrying out aeroelastic analysis considering the influence of the internal fluid.
The invention has the advantages and beneficial effects that:
the invention analyzes the influence of the dynamic characteristic of the internal fluid on the aeroelasticity characteristic of the airplane and provides an aeroelasticity analysis method considering the influence of the dynamic characteristic of the internal fluid. The method can consider the influence of the dynamic characteristics of the internal fluid in a dynamic simulation model with a large number of internal fluid airplanes, thereby being beneficial to improving the aeroelasticity analysis precision of the whole airplane and better guiding aeroelasticity design. The proposed analysis method is simple and convenient to operate and has high engineering application value.
Drawings
Figure 1 is a drawing of the present invention with internal fluid zone and interface definitions,
FIG. 2 is a comparison curve of the analysis results of the present invention and the calculation results of the fixed mass model.
Detailed Description
The invention provides a method for aeroelastic analysis considering the influence of the dynamic characteristics of internal fluid. By adding an additional mass matrix to the mass matrix of the equation of motion, consideration of the influence of the internal fluid dynamics is achieved.
(1) Generating additional mass matrix considering internal fluid
Assuming that the internal fluid is non-viscous and incompressible, the area occupied by the solid is Ω S The area occupied by the internal fluid is omega F The free surface is gamma, the interface between the fluid and the solid is sigma, and the schematic diagram is shown in figure 1. Euler equation of the internal fluid is
Where p, p F And u F Pressure, density and displacement of the fluid, respectively, and t is time;
the continuity equation of the inviscid incompressible fluid is
Wherein v is F Is the fluid velocity;
the equation (a) and (b) can be used to derive the equation in the region omega F Inner part
Δφ＝0 (c)
Wherein phi is a potential function, satisfies
Assuming that both fluid motion and structural vibrations are simple harmonic vibrations, i.e.
The conservation of momentum at the fluid-solid interface sigma is
Wherein sigma S Indicating the structure in the solid region omega S Stress of p s And u S Respectively representing the structure density and the displacement, n S Normal to the surface of the structure;
the non-penetration condition at the fluid-solid coupling interface is
Wherein n is F Is the normal direction on the fluid side to the solid side on the fluid-solid interface sigma.
The potential function phi satisfies the boundary condition on the free surface gamma
φ＝0 (f)
Deducing according to formulas (a) to (f) to obtain a matrix equation system
Wherein M is S Is a structural mass matrix, K is a structural stiffness matrix, M F λ = ω for the internal fluid mass matrix 2 ω is the circular frequency, T is the matrix of transfer conditions at the fluid-solid coupling interface (equations (d) and (e)), superscript T is the transpose, u is the structural displacement vector, and Φ is the fluid potential function.
The second row of the matrix equation (g) is obtained
Simplified equation (g) of
K 11 -λ(M S +M A ) 11 ＝0 (i)
Wherein
I.e. an additional mass matrix taking into account the internal fluid mass.
(2) Assembling equations of motion after considering additional mass matrices
The motion equation of the system consisting of the internal fluid and the container is derived as
K U -λ(M S +M A )u＝0
Equation (i) where M S Is a structural mass matrix, K is a structural stiffness matrix, and lambda = omega 2 ω is the circular frequency and u is the structural displacement vector.
(3) Performing aeroelastic analysis taking into account internal fluid effects
The updated mass matrix is used to calculate the aeroelastic properties taking into account the internal fluid effects.
The analysis method is characterized in that the mass matrix in the motion equation is added to consider the influence of the dynamic characteristics of the internal fluid.
Table 1 defines the calculation states of the present invention, where state 1 is the analysis result of the model considering the influence of the internal fluid, and state 2 is the calculation result of the fixed mass model.
TABLE 1 inventive computational State definition
| State numbering | State definition |
| 1 | Model of considering internal fluid effects |
| 2 | Lumped mass model |
Claims (2)
1. A method of aeroelastic analysis that takes into account the dynamic behavior of the internal fluid; the method is characterized in that: the method comprises the following steps:
step one, generating an additional mass matrix considering the internal fluid,
step two, assembling the motion equation after considering the additional mass matrix,
thirdly, performing aeroelastic analysis considering the influence of the internal fluid;
the derivation method of the additional quality matrix in the first step is as follows: assuming the internal fluid is non-viscous and incompressible, the area occupied by the solid is Ω S The internal fluid occupies an area of Ω F With free surface of gamma, intersection of fluid and solidThe interface is sigma;
euler equation of the internal fluid is
Where p, p F And u F Pressure, density and displacement of the fluid, respectively, and t is time;
the continuity equation of the inviscid incompressible fluid is
Wherein v is F Is the fluid velocity;
from the equations (a) and (b), the region Ω can be derived F Inner part
Δφ＝0 (c)
Wherein phi is a potential function, satisfies
Assuming that the fluid motion and the structural vibration are both simple harmonic vibrations, i.e. p = ρ F ω 2 Phi, the conservation of momentum at the fluid-solid interface sigma is
σ S n S ＝ρ F ω 2 φn F (d)
Wherein sigma S Shows the structure in the solid region omega S Stress of p S And u S Respectively representing the structure density and the displacement, n S Normal to the surface of the structure;
the non-penetration condition at the fluid-solid coupling interface is
Wherein n is F A normal direction pointing from the fluid side to the solid side on the fluid-solid interface Σ;
the potential function phi satisfies the boundary condition on the free surface gamma
φ＝0 (f)
Deducing according to formulas (a) - (f) to obtain a matrix equation system
Wherein M is y Is a structural mass matrix, K is a structural stiffness matrix, M F Is an internal fluid mass matrix, λ = ω 2 Omega is the circular frequency, T is the transfer condition matrix at the fluid-solid coupling interface, namely the equations (d) and (e), the superscript T is the transposition, u is the structure displacement vector, and phi is the fluid potential function;
the second row of the matrix equation (g) is obtained
Simplifying the equation (g) to
Ku-λ(M S +M A )u＝0 (i)
Wherein
I.e. an additional mass matrix taking into account the internal fluid mass.
2. The aeroelastic analysis method considering internal fluid dynamics according to claim 1, wherein: the equation of motion in the second step is as follows:
Ku-λ(M S +M A )u＝0
wherein M is S Is a structural mass matrix, K is a structural stiffness matrix, and lambda = omega 2 ω is the circular frequency and u is the structure displacement vector.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811476560.9A CN109726440B (en) | 2018-12-04 | 2018-12-04 | Aeroelasticity analysis method considering dynamic characteristics of internal fluid |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811476560.9A CN109726440B (en) | 2018-12-04 | 2018-12-04 | Aeroelasticity analysis method considering dynamic characteristics of internal fluid |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109726440A CN109726440A (en) | 2019-05-07 |
| CN109726440B true CN109726440B (en) | 2023-04-18 |
Family
ID=66294979
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811476560.9A Active CN109726440B (en) | 2018-12-04 | 2018-12-04 | Aeroelasticity analysis method considering dynamic characteristics of internal fluid |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109726440B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110887636B (en) * | 2019-11-21 | 2021-09-14 | 中国特种飞行器研究所 | Hydrodynamic model test device and method for additional mass of airship |
Family Cites Families (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA2537933C (en) * | 2003-08-01 | 2015-10-13 | Cidra Corporation | An apparatus and method for providing a density measurement augmented for entrained gas |
| US8155794B2 (en) * | 2007-05-09 | 2012-04-10 | Tao Of Systems Integration, Inc. | System and method for control of aeroelasticity effects |
| DE102009002392A1 (en) * | 2009-04-15 | 2010-11-04 | Airbus Deutschland Gmbh | System and method for determining local accelerations, dynamic load distributions and aerodynamic data in an aircraft |
| US20160052621A1 (en) * | 2009-07-10 | 2016-02-25 | Peter Ireland | Energy efficiency improvements for turbomachinery |
| US8366057B2 (en) * | 2009-07-28 | 2013-02-05 | University Of Kansas | Method and apparatus for pressure adaptive morphing structure |
| CN102012953B (en) * | 2010-11-04 | 2013-05-08 | 西北工业大学 | CFD (computational fluid dynamics)/CSD (circuit switch data) coupled solving nonlinear aeroelasticity simulation method |
| CN102938003B (en) * | 2012-10-17 | 2014-12-03 | 北京航空航天大学 | Method for predicting aeroelasticity stability numerical value of turbomachinery with error frequency included |
| EP2948369A4 (en) * | 2013-01-25 | 2017-01-18 | Peter S. Ireland | Energy efficiency improvements for turbomachinery |
| CN105843073B (en) * | 2016-03-23 | 2018-08-17 | 北京航空航天大学 | A kind of wing structure aeroelastic stability analysis method not knowing depression of order based on aerodynamic force |
| CN107066653B (en) * | 2016-12-15 | 2020-12-29 | 中国航空工业集团公司西安飞机设计研究所 | Aeroelasticity analysis method considering dynamic characteristics of engine |
-
2018
- 2018-12-04 CN CN201811476560.9A patent/CN109726440B/en active Active
Non-Patent Citations (2)
| Title |
|---|
| 叶正寅,王刚,杨永年,杨炳渊.基于欧拉方程的一种机翼气动弹性计算方法.西北工业大学学报.2002,第20卷(第02期),第98-102页. * |
| 马戎 ; 常兴华 ; 赫新 ; 张来平 ; .流动/运动松耦合与紧耦合计算方法及稳定性分析.气体物理.2016,第1卷(第06期),第36-49页. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109726440A (en) | 2019-05-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Smith et al. | CFD-based analysis of nonlinear aeroelastic behavior of high-aspect ratio wings | |
| Blumenthal et al. | Computational investigation of a boundary-layer-ingestion propulsion system | |
| CN104133933A (en) | Pneumatic elastic mechanical characteristic analytical method of hypersonic speed aircraft in thermal environment | |
| Chiba et al. | Multidisciplinary design optimization and data mining for transonic regional-jet wing | |
| Goura | Time marching analysis of flutter using computational fluid dynamics | |
| Tamaki et al. | Wall modeling for large-eddy simulation on non-body-conforming Cartesian grids | |
| CN109726440B (en) | Aeroelasticity analysis method considering dynamic characteristics of internal fluid | |
| Gupta et al. | Aeroelastic simulation of hypersonic flight vehicles | |
| Blumenthal et al. | Computational investigation of a boundary layer ingestion propulsion system for the common research model | |
| He et al. | A coupled newton-krylov time spectral solver for wing flutter and lco prediction | |
| Iyer et al. | Wall-modeled LES of the three-dimensional speed bump experiment | |
| Houtman et al. | Resolvent analysis of large aircraft wings in edge-of-the-envelope transonic flow | |
| CN108595893A (en) | A kind of three-dimensional mechanical Modal Analysis analogy method based on three layers of pretreatment | |
| Ishikawa et al. | Sonic boom assessment in primary boom carpet of low-boom supersonic airplane (NASA C25D) | |
| Zheng et al. | Flutter analyses of complete aircraft based on hybrid grids and parallel computing | |
| CN107944116B (en) | A kind of efficient predicting method of transient energy response for Time variable structure | |
| Perkins et al. | Coupled Analysis of a Slotted Natural-Laminar-Flow Transonic Truss-Braced Wing Aircraft Configuration | |
| Xu et al. | Modern Computational Aeroelasticity | |
| Lawless et al. | CFD-Based Kriging Surrogate Models Compared to Axisymmetric Missile Concept From Mach 0.60 to 3.95 | |
| Mehar | Transient effect of aircraft propeller blade by using composites | |
| Blades et al. | Aeroelastic stability predictions of a business jet landing gear door using high fidelity fluid-structure interaction tools | |
| Arizono et al. | Flutter simulations of a T‐tail configuration using non‐linear aerodynamics | |
| Hoshi et al. | Numerical simulations of transonic flutter on a three-dimensional wing | |
| Rizzi et al. | Computation of Aeroelastic Effects on F-16XL at Flight Conditions FC70 and FC25 | |
| Jirasek et al. | Aeroelastic Analysis of the CAWAPI F-16XL Configuration at Transonic Speeds |
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 |