CN113218615A - Equivalent method of distributed aerodynamic force and limited excitation point excitation load - Google Patents
Equivalent method of distributed aerodynamic force and limited excitation point excitation load Download PDFInfo
- Publication number
- CN113218615A CN113218615A CN202110620988.1A CN202110620988A CN113218615A CN 113218615 A CN113218615 A CN 113218615A CN 202110620988 A CN202110620988 A CN 202110620988A CN 113218615 A CN113218615 A CN 113218615A
- Authority
- CN
- China
- Prior art keywords
- excitation
- aerodynamic
- aerodynamic force
- point
- modal
- 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
- 230000005284 excitation Effects 0.000 title claims abstract description 31
- 238000000034 method Methods 0.000 title claims abstract description 22
- 238000006073 displacement reaction Methods 0.000 claims description 12
- 239000011159 matrix material Substances 0.000 claims description 7
- 239000013598 vector Substances 0.000 claims description 4
- 238000012360 testing method Methods 0.000 abstract description 18
- 230000008569 process Effects 0.000 abstract description 4
- 230000009471 action Effects 0.000 abstract description 2
- 238000005516 engineering process Methods 0.000 abstract description 2
- 238000004364 calculation method Methods 0.000 description 4
- 238000002474 experimental method Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000010998 test method Methods 0.000 description 3
- 238000009434 installation Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000001687 destabilization Effects 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M9/00—Aerodynamic testing; Arrangements in or on wind tunnels
- G01M9/06—Measuring arrangements specially adapted for aerodynamic testing
-
- 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
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- 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
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Automation & Control Theory (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Aviation & Aerospace Engineering (AREA)
- Fluid Mechanics (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
An equivalent method of distributed aerodynamic force and limited excitation point excitation load belongs to the field of aircraft ground aerodynamic tests, and specifically relates to an aerodynamic force equivalent technology in the aerodynamic test. The invention aims to solve the problem that the pneumatic action of the distributed aerodynamic force on the structure is equivalently simulated by a limited excitation point in the structure pneumatic flutter test. In the process of the structure pneumatic flutter test, the aerodynamic force is integrated along the surface of the structure, the excitation load of the equivalent finite excitation point is made to work equally to the structure by the aerodynamic force based on the energy equivalence principle, and then the equivalent excitation load of each excitation point is obtained by solving.
Description
Technical Field
The invention belongs to the field of aircraft ground pneumatic tests, and particularly relates to a structural pneumatic test method, in particular to a pneumatic equivalent technology in a pneumatic test.
Background
Aerodynamic flutter is a aerodynamic phenomenon that seriously threatens the stability and safety of an aircraft. The existing test research method is mainly a wind tunnel method, but the wind tunnel has the outstanding problems of complex structure, high test cost, high experiment technical difficulty and the like, so that the finding of a simple and efficient aircraft fluid-solid coupling ground wind tunnel test method has important practical significance.
The existing ground wind tunnel experiment method mainly tests the vibration displacement and the speed of the aircraft structure through an eddy current sensor, a laser displacement sensor and the like, obtains the aerodynamic force value of a measuring point based on empirical formulas such as a piston theory and the like, and applies corresponding aerodynamic force to the structure by using a vibration exciter so as to analyze the vibration response and the aerodynamic stability of the aircraft structure. Some of the existing test methods are based on finite element methods, and have low calculation speed and no real-time property; there are some methods based on equivalence, but the equivalence process lacks a strict theoretical basis.
The aerodynamic load borne by the aircraft structure belongs to distributed load, and the ground test is difficult to completely and truly simulate. The distributed load can be reduced to a plurality of concentrated loads according to an equivalent algorithm and then applied to the structure by the actuators. The equivalent algorithm needs to synthesize the test results of a plurality of sensors and calculate the result of the equivalent aerodynamic force in real time, so that high requirements on the accuracy and the real-time performance of the algorithm are provided, and the method is a key and difficult point for realizing the ground aerodynamic test of the aircraft.
Disclosure of Invention
The invention aims to solve the problem that the pneumatic action of distributed aerodynamic force on a structure is simulated equivalently by a limited excitation point in a structural pneumatic flutter test, provides an equivalent method of the distributed aerodynamic force and the excitation load of the limited excitation point, and realizes a simple, efficient and accurate structural pneumatic flutter ground test.
In order to achieve the purpose, the technical scheme adopted by the invention is as follows:
an equivalent method of distributed aerodynamic force and limited excitation point excitation load comprises the following steps:
arranging m piezoelectric actuators on the upper surface of the wallboard to apply excitation to simulate aerodynamic force, wherein the laying positions of the m piezoelectric actuators are aj(xj,yj) The force applied is fjWherein j is 1 … m;
n piezoelectric sensors are arranged on the upper surface or the lower surface of the wall plate for measuring displacement and speed, and the laying positions of the n piezoelectric sensors are sk(xk,yk) Measuring the displacement wkAnd velocity wkWherein k is 1 … n;
for a four-sided simply supported plate, the lateral displacement at any point on the plate at time t is expressed as:
wherein i is a modal order, l is a truncated modal order, and q is a linear orderi(t) is the modal coordinates at time t,to satisfy the lattice function of the boundary condition of the wall plate, as shown in fig. 3, a and b are the geometric dimensions of the flat plate, i.e., the length and width;
wherein M is a matrix formed by first-order modal vectors of the structure, and each-order modal vector comprises the modal of n measuring points; the expression of M is as follows:
obtaining q by the above formulaiAndwith respect to wk,Calculating and solving the displacement and the speed of any point on the wall plate;
work P of aerodynamic force in unit time of whole wall plateΔPIs equal to the work done by m piezoelectric actuators per unit time of excitationNamely:
wherein:
wherein, Δ p is the acting force of the airflow at each point on the surface of the wall plate, and can be obtained by the aeroelasticity theory; delta is a dirac function, and when the independent variable is not 0, the function value is constantly zero;
the formula (5) and the formula (6) may be taken into the formula (4):
namely:
wherein:
solving the m piezoelectric actuators according to the formulas (1) to (7) to obtain the force f required to be applied to the m piezoelectric actuatorsj,ΔpiWork, K, done for aerodynamic forcesijDenotes a symbol introduced for simplifying the expression of the formula (7), and has no particular physical meaning, v∞Infinite incoming flow velocity and gamma is the adiabatic coefficient. To represent
Further, regarding the matrix M, it is necessary to intercept the modes that are equal to the number of stations and are uncorrelated, so as to ensure that M is an invertible square matrix, i.e. n ═ l.
Compared with the prior art, the invention has the beneficial effects that: the equivalent method provided by the invention based on the energy equivalent principle has the outstanding advantages of simple calculation process, high processing efficiency and the like. In the process of the structure pneumatic flutter test, the aerodynamic force is integrated along the surface of the structure, the excitation load of the equivalent finite excitation point is made to work equally to the structure by the aerodynamic force based on the energy equivalence principle, and then the equivalent excitation load of each excitation point is obtained by solving.
Drawings
FIG. 1 is a hardware composition diagram of a wall flutter test in an embodiment;
FIG. 2 is a schematic view of the wall panel boundary installation of the embodiment;
FIG. 3 is a schematic diagram of the sensor laying position in the embodiment;
FIG. 4 shows dynamic pressure λ 580 < λcrA time history chart of transverse vibration of the time wall plate;
FIG. 5 shows dynamic pressure λ 620 > λcrTime history chart of transverse vibration of the time wall plate.
Detailed Description
The technical solution of the present invention will be clearly and completely described below with reference to the accompanying drawings and embodiments. It is obvious that the described embodiments are only a part of the embodiments of the invention, not all embodiments, and all other embodiments obtained by those skilled in the art without any inventive work are within the scope of the present invention.
Example 1:
pneumatic flutter ground equivalent experiment of plate structure
The plate-shaped structure is very common, and the aerodynamic characteristics of the plate-shaped structure relate to the safe operation of engineering structures such as wings and blades. In this embodiment, a plate structure flutter experiment is taken as an example, a composition diagram of experimental hardware is shown in fig. 1, and a structure installation is shown in fig. 2. The aerodynamic loads are equivalent to the excitation loads of 8 excitation points as shown in fig. 3, according to the method described in the patent content. The specific parameters of the plate structure are as follows, the geometric parameters are 0.8m 0.6m 0.002m in length and width, and the density is 2689kg/m3The modulus of elasticity was 68.3 GPa.
In the embodiment, a rectangular plate structure is taken as an example, exciting force is applied at 8 exciting points, the test result and the numerical calculation result show that the pneumatic instability critical speed of the plate is 717.5m/s and has better consistency, and the calculation result is efficient, accurate and reliable.
As shown in fig. 4 and 5, it can be seen that when the dynamic pressure is taken 580, the structure is in a stable state, and the structural vibration caused by external excitation is gradually attenuated and approaches to zero; when dynamic pressure is taken 620, the structure is in a flutter destabilization state, and external disturbance causes the structure to enter a self-excited vibration state of periodic motion. Therefore, the method can well simulate the aerodynamic flutter phenomenon of the aircraft structure through ground tests, and also proves the feasibility, reliability and convenience of the equivalent method of distributed aerodynamic force and limited excitation point excitation load.
Claims (2)
1. An equivalent method of distributed aerodynamic force and limited excitation point excitation load is characterized in that: the method comprises the following steps:
arranging m piezoelectric actuators on the upper surface of the wallboard to apply excitation to simulate aerodynamic force, wherein the laying positions of the m piezoelectric actuators are aj(xj,yj) The force applied is fjWherein j is 1 … m;
n piezoelectric sensors are arranged on the upper surface or the lower surface of the wall plate for measuring displacement and speed, and the laying positions of the n piezoelectric sensors are sk(xk,yk) Measuring the displacement wkAnd velocityWherein k is 1 … n;
for a four-sided simply supported plate, the lateral displacement at any point on the plate at time t is expressed as:
wherein i is a modal order, l is a truncated modal order, and q is a linear orderi(t) is the modal coordinates at time t,in order to satisfy the matrix function of the boundary condition of the wall plate, a and b are the geometric dimensions of the flat plate, namely the length and the width;
wherein M is a matrix formed by first-order modal vectors of the structure, and each-order modal vector comprises the modal of n measuring points; the expression of M is as follows:
obtaining q by the above formulaiAndwith respect to wk,Calculating and solving the displacement and the speed of any point on the wall plate;
work P of aerodynamic force in unit time of whole wall plateΔPIs equal to the work done by m piezoelectric actuators per unit time of excitationNamely:
wherein:
wherein, Δ p is the acting force of the airflow at each point on the surface of the wall plate, and can be obtained by the aeroelasticity theory; delta is a dirac function, and when the independent variable is not 0, the function value is constantly zero;
the formula (5) and the formula (6) may be taken into the formula (4):
namely:
wherein:
solving the m piezoelectric actuators according to the formulas (1) to (7) to obtain the force f required to be applied to the m piezoelectric actuatorsj,ΔpiWork, K, done for aerodynamic forcesijDenotes a symbol introduced for simplifying the expression of the formula (7), and has no particular physical meaning, v∞Infinite incoming flow velocity and gamma is the adiabatic coefficient.
2. The equivalent method of distributed aerodynamic and finite excitation point excitation loads according to claim 1, characterized in that: regarding the matrix M, it is necessary to intercept the modes that are equal to the number of stations and are uncorrelated to ensure that M is an invertible square matrix, i.e. n ═ l.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110620988.1A CN113218615A (en) | 2021-06-03 | 2021-06-03 | Equivalent method of distributed aerodynamic force and limited excitation point excitation load |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110620988.1A CN113218615A (en) | 2021-06-03 | 2021-06-03 | Equivalent method of distributed aerodynamic force and limited excitation point excitation load |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113218615A true CN113218615A (en) | 2021-08-06 |
Family
ID=77082636
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110620988.1A Pending CN113218615A (en) | 2021-06-03 | 2021-06-03 | Equivalent method of distributed aerodynamic force and limited excitation point excitation load |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113218615A (en) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102866637A (en) * | 2012-10-07 | 2013-01-09 | 西北工业大学 | Quadratic order-reduction based method for simulating unsteady aerodynamic force of aerofoil with operation surface |
CN104133933A (en) * | 2014-05-29 | 2014-11-05 | 温州职业技术学院 | Pneumatic elastic mechanical characteristic analytical method of hypersonic speed aircraft in thermal environment |
CN109889094A (en) * | 2019-03-12 | 2019-06-14 | 哈尔滨工业大学 | A kind of double oscillator piezoelectric harvesters of tunable aerofoil profile flutter Exciting-simulator system |
CN109933876A (en) * | 2019-03-03 | 2019-06-25 | 西北工业大学 | A kind of unsteady aerodynamic force order reducing method based on broad sense aerodynamic force |
CN111324991A (en) * | 2019-12-10 | 2020-06-23 | 中国飞机强度研究所 | Reconstruction method of aerodynamic model in ground flutter test |
-
2021
- 2021-06-03 CN CN202110620988.1A patent/CN113218615A/en active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102866637A (en) * | 2012-10-07 | 2013-01-09 | 西北工业大学 | Quadratic order-reduction based method for simulating unsteady aerodynamic force of aerofoil with operation surface |
CN104133933A (en) * | 2014-05-29 | 2014-11-05 | 温州职业技术学院 | Pneumatic elastic mechanical characteristic analytical method of hypersonic speed aircraft in thermal environment |
CN109933876A (en) * | 2019-03-03 | 2019-06-25 | 西北工业大学 | A kind of unsteady aerodynamic force order reducing method based on broad sense aerodynamic force |
CN109889094A (en) * | 2019-03-12 | 2019-06-14 | 哈尔滨工业大学 | A kind of double oscillator piezoelectric harvesters of tunable aerofoil profile flutter Exciting-simulator system |
CN111324991A (en) * | 2019-12-10 | 2020-06-23 | 中国飞机强度研究所 | Reconstruction method of aerodynamic model in ground flutter test |
Non-Patent Citations (1)
Title |
---|
邵崇晖: ""超声速流中壁板颤振的抑制和地面试验研究"", 《中国优秀博士学位论文全文数据库工程科技Ⅱ辑》 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Arrieta et al. | Modelling and configuration control of wing-shaped bi-stable piezoelectric composites under aerodynamic loads | |
EP3770575B1 (en) | Multi-dimensional vibration control method for sting-mounted aircraft model | |
Tschanz et al. | The base balance technique for the determination of dynamic wind loads | |
Bryant et al. | Development of an aeroelastic vibration power harvester | |
Epureanu et al. | Reduced-order models of unsteady transonic viscous flows in turbomachinery | |
Yao et al. | Chaotic motion of a composite laminated plate with geometric nonlinearity in subsonic flow | |
Zeng et al. | GVT-based ground flutter test without wind tunnel | |
Barzegari et al. | Aeroelastic characteristics of cantilever wing with embedded shape memory alloys | |
CN103942381B (en) | State near field dynamics method used for predicting airplane aluminum alloy structure performance | |
CN109840349A (en) | A kind of fixed wing aircraft gust response modeling and analysis methods | |
Maheri et al. | WTAB, a computer program for predicting the performance of horizontal axis wind turbines with adaptive blades | |
Fazelzadeh et al. | Active control law design for flutter suppression and gust alleviation of a panel with piezoelectric actuators | |
CN113218615A (en) | Equivalent method of distributed aerodynamic force and limited excitation point excitation load | |
Kim et al. | An optimal reduced-order aeroelasltic modeling based on a response-based modal analysis of unsteady CFD models | |
CN111651920B (en) | Transmission analysis method for wind-induced mechanical energy of long-span power transmission tower under tower-line system | |
CN109946029B (en) | Multi-frequency fatigue test method for turbocharger | |
Li et al. | Improving the aeroelastic flutter characteristics of supersonic beams using piezoelectric material | |
Zhang et al. | A nonlinear flow-induced energy harvester by considering effects of fictitious springs | |
Taylor et al. | A structure-coupled CFD method for time-marching flutter analysis | |
Nandi et al. | Prediction and analysis of the nonsteady transitional boundary layer dynamics for flow over an oscillating wind turbine airfoil using the γ-reθ transition model | |
Leng et al. | Structural dynamic analysis of a tidal current turbine using MBDyn | |
An et al. | Large Deformation Prediction of Wing Model and Geometric Nonlinear Aeroelastic Analysis | |
Arena et al. | Dynamic morphing of elastic plates via principal parametric resonance | |
Yu et al. | Modelling of flutter characteristics of aircraft wing with pylon-mounted engine nacelle | |
CN113358311A (en) | Plate/beam structure transverse displacement determination method based on limited measuring points and vibration modes |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20210806 |
|
RJ01 | Rejection of invention patent application after publication |