CN102968525B - Method for determining overweight ratio of plane flutter model - Google Patents
Method for determining overweight ratio of plane flutter model Download PDFInfo
- Publication number
- CN102968525B CN102968525B CN201210451712.6A CN201210451712A CN102968525B CN 102968525 B CN102968525 B CN 102968525B CN 201210451712 A CN201210451712 A CN 201210451712A CN 102968525 B CN102968525 B CN 102968525B
- Authority
- CN
- China
- Prior art keywords
- mrow
- mtd
- mfrac
- model
- section
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 25
- 238000004364 calculation method Methods 0.000 claims description 7
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 3
- 210000005069 ears Anatomy 0.000 description 2
- 230000002349 favourable effect Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Feedback Control In General (AREA)
Abstract
The invention belongs to the field of aviation structural mechanics, and relates to a design method for obtaining the overweight ratio of a flutter model in the design of the flutter model. The method comprises the following steps of: calculating the equivalent mass mBE of a reference beam frame according to a reference beam frame model, and calculating the equivalent mass ratio nBE of the reference beam frame in combination with the component mass mA; calculating a design scale coefficient ndss; specifying a section shape coefficient ns; specifying the mass ratio kb of a flutter model beam; and calculating the overweight ratio nG of the flutter model. In According to the method disclosed by the invention, a reverse design method is used for replacing a trial and error method adopted in the conventional model design, so that various factors possibly influencing model weight and mathematical relations thereof are considered at the start of model design, the overweight ratio of the model is further obtained, the uncertainty of model design is reduced, and the model design period is shortened greatly.
Description
Technical Field
The invention belongs to the field of aeronautical structure mechanics, and particularly relates to a method for determining an overweight ratio of an airplane flutter model.
Background
The flutter model can be used to obtain the flutter characteristics of the aircraft and its components. In flutter model design, the interference of overweight problem is often encountered, especially the design of high-speed flutter model usually needs small beam frame mass to provide large wing surface rigidity, and the overweight problem of flutter model design is more emphasized. In addition, for some low-speed flutter model design problems, the situation is just opposite, namely the target mass is overlarge, if the flutter model design meets the target mass, the static deformation of the model is larger, and further the test is not favorable, so that an 'overweight ratio' which is smaller than 1 in numerical value is correspondingly added.
In the past, a basic scale of model design is selected and a scale is derived, a certain structural form is selected according to working experience, a flutter model is directly designed, after a model main body is designed, comparison is carried out according to estimated quality and target quality, and if the estimated quality of a model exceeds the target quality, redesign is often carried out at an additional overweight ratio, or redesign is carried out after the structural form of the model and the design scale are adjusted. This method has the following disadvantages:
firstly, the model design seriously depends on the experience of designers, and the uncertainty of the model design is increased;
secondly, the super-weight ratio is acquired by carrying out multiple rounds of redesign and calculation, the time is long, the efficiency is low, and the design period of the model after severe dragging influences the development progress of the airplane.
Disclosure of Invention
The purpose of the invention is:
the invention aims to obtain a method for determining the overweight ratio of a flutter model, which improves the working efficiency, shortens the model design period and changes the original laggard method which depends heavily on experience and needs to be designed by repeated trial and error in the past.
The technical scheme of the invention is as follows:
a method for determining the overweight ratio of airplane flutter model features that the flutter model with reference beam frame is ensured, which includes the mass m of airplaneABeam frame section characteristic moment of inertia Ix、IyHarmonic polar moment of inertia J, high speed flutter model design basic scale speed-pressure ratio etaqAnd density ratio ηρOr the basic scale speed ratio eta of the low-speed flutter modelVAnd density ratio ηρOn the basis of the above-mentioned parameters, determining overweight ratio n of flutter model designG(ii) a Order toUnifying high and low speed flutter models; the specific determination steps are as follows:
step one, calculating the equivalent mass ratio n of the reference beam frame of the airplaneBEThe calculation method is as follows:
1.1, calculating to obtain the equivalent mass m of the reference beam frameBE:
Moment of inertia I according to beam frame section characteristicsx、IyAnd polar inertia moment J, calculating the corresponding rectangular section size with the lug, and further obtaining a reference beam frame model with the section size; giving a beam frame material density to obtain the equivalent mass m of the reference beam frameBE(ii) a From the section characteristic moment of inertia Ix、IyAnd the polar inertia moment J, the method for obtaining the rectangular section size with the lug is as follows:
1.1.1 for a rectangular cross-section, the iterative formula for the dimensions of the rectangular cross-section is:
wherein,
first, alpha is set to alpha00.96, with formula [1]B/a is calculated and is according to the formula [2 ]]Calculating alpha, substituting into equation [1 ]]Through iteration, alpha is converged to a stable value finally, and then a and b can be obtained;
1.1.2 for rectangular cross-section with ear, let the ear thickness
t=0.2b ……………………………[3]
1.1.3 adjusting J
J1=J-(t/2b)2] ……………………………[4]
1.1.4 according to the target value IxAnd J1And repeating the step 1.1.1 to obtain a1And b1;
1.1.5 according to the target value IyAnd a1、b1And t value, obtaining the half width l of the section with the lug;
finally obtaining the cross section size 2a of the rectangular beam with the lug plates1、2b1T and 2 l;
1.2, calculating the equivalent mass ratio n of the reference beam frameBE:
nBE=mBE/mA ……………………………[6]
Step two, calculating a coefficient n of a design scaledss:
Step three, specifying a section shape coefficient ns:
nsThe designation of (c) is divided into two cases:
if the actually adopted main beam section form is also a rectangular beam section form with lugs
ns=1 ……………………………[8]
If the actually adopted section form of the main beam is a rectangular thin-walled beam with lugs, the section form of the main beam is as follows
0.3≤ns≤1 ……………………………[9]
In this case nsThe specific numerical value can be obtained according to specific conditions, and the specific method is obtained by adjusting the thickness parameter of the section thin wall;
step four, specifying the mass ratio k of the flutter model beamb:
0.25≤kb≤0.4 ……………………………[10]
Wherein,
in the case of the low-speed flutter model,kbis small; in the case of the high-speed flutter model, kb2, partial enlargement;
step five, calculating the overweight ratio n of the flutter modelG:
nG=nsndssnBE/kb ……………………………[11]
Overweight ratio n of flutter modelGThe control is carried out in the following range as much as possible:
05≤nG≤20 ……………………………[12]
if it is out of this range, the formula [7 ] can be adjusted]To the formula [10]Obtaining the overweight ratio n of the flutter model meeting the requirementsG;
If 0.9. ltoreq. nGIf n is less than or equal to 1.1, directly taking nG=1。
The invention has the advantages that:
because the invention uses the reverse design thinking to replace the trial and error method adopted in the past model design, various factors which can influence the model weight and the mathematical relationship of the influencing factors are considered at the beginning of the model design, and the model overweight ratio is further obtained. According to the invention, the weight characteristic of the flutter model beam can be conveniently obtained only by adjusting related design parameters, so that the uncertainty of model design is reduced; the work which needs to be finished for several days in the past can be finished within several minutes at present, so that the model design period is greatly shortened; the invention has clear and smooth form, is very suitable for the programming and operation of various general calculation programs and is easy to be mastered by engineering technicians, thereby improving the design efficiency of the flutter model and bringing convenience to the test work of the flutter model and the development of an airplane.
Drawings
FIG. 1 is a schematic view in rectangular section; where a is the half width of the rectangle calculated in step one of the present invention; b is the half height of the rectangle calculated in step one of the present invention; the origin o of the two-dimensional coordinate system in fig. 1 is the center of the rectangle, the x-axis is parallel to the width direction of the rectangle, the positive direction is towards the right, and the positive direction of the y-axis is towards the upper.
FIG. 2 is a schematic view of a rectangular cross section with tabs; the two-dimensional coordinate system in fig. 2 is the same as in fig. 1. In the figure a1Is the half width of the rectangle calculated in step one of the invention; b1 is the half height of the rectangle calculated in step one of the present invention, t is the thickness of the tab calculated in step one of the present invention; in the figure, l is the calculated half width of the cross section at the first step of the present invention.
Detailed Description
The present invention will be described in further detail with reference to the accompanying drawings, referring to fig. 1 and 2.
A method for determining the overweight ratio of airplane flutter model features that the flutter model with reference beam frame (including the mass m of airplane) is ensuredABeam frame section characteristic moment of inertia Ix、IyPolar moment of inertia J) and high-speed flutter model design basic scale speed-pressure ratio etaqAnd density ratio ηρOr the basic scale speed ratio eta of the low-speed flutter modelVAnd density ratio ηρOn the basis of the above-mentioned parameters, determining overweight ratio n of flutter model designG(ii) a The basic scale is defined in handbook of aeroelasticity mechanics of airplanes (tubular, published by aeronautical industry, 1994); order toTherefore, the high-speed flutter model and the low-speed flutter model can be unified under a set of design system; the specific determination steps are as follows:
step one, calculating the equivalent mass ratio n of the reference beam frame of the airplaneBEThe calculation method is as follows:
1.1, calculating to obtain the equivalent mass m of the reference beam frameBE:
Moment of inertia I according to beam frame section characteristicsx、IyAnd polar inertia moment J, calculating the corresponding rectangular section size with the lug, and further obtaining a reference beam frame model with the section size; giving a beam frame material density to obtain the equivalent mass m of the reference beam frameBE(ii) a From the section characteristic moment of inertia Ix、IyAnd the polar inertia moment J, the method for obtaining the rectangular section size with the lug is as follows:
1.1.1 referring to the rectangular cross-section shown in FIG. 1, the iterative formula for the dimensions of the rectangular cross-section:
wherein,
first, alpha is set to alpha00.96, with formula [1]B/a is calculated and is according to the formula [2 ]]Calculating alpha, substituting into equation [1 ]]Through iteration, α is finally converged to a stable value, and a and b can be obtained.
1.1.2 rectangular cross-section with ears as shown in reference to FIG. 2, making the ears thick
t=0.2b ……………………………[3]
1.1.3 adjusting J
J1=J[1-(t/2b)2] ……………………………[4]
1.1.4 according to the target value IxAnd J1And repeating the step 1.1.1 to obtain a1And b1。
1.1.5 according to the target value IyAnd a1、b1And t, the half width l of the section with the ear can be obtained.
Finally obtaining the cross section size 2a of the rectangular beam with the lug plates1、2b1T and 2 l.
1.2, calculating the equivalent mass ratio n of the reference beam frameBE:
nBE=mBE/mA ……………………………[6]
nBEThe structural efficiency of the aircraft component is reflected to a certain extent;
step two, calculating a coefficient n of a design scaledss:
Step three, specifying a section shape coefficient ns:
nsThe physical meaning of the expression is that the moment of inertia I is the same in section characteristicx、IyAnd the ratio of the area corresponding to the actually adopted main beam section form to the area corresponding to the lug-carrying rectangular beam section form under the polar inertia moment J; there are two cases:
3.1, if the section form of the actually adopted main beam is also the section form of the rectangular beam with the lugs
ns=1 ……………………………[8]
3.2, if the actually adopted section form of the main beam is the section form of the rectangular thin-walled beam with the lugs, then
0.3≤ns≤1 ……………………………[9]
In this case nsThe specific numerical value can be obtained according to specific conditions, and the specific method is obtained by adjusting the thickness parameter of the section thin wall;
step four, specifying the mass ratio k of the flutter model beamb:
0.25≤kb≤0.4 ……………………………[10]
Wherein,
in the case of the low-speed flutter model, kbGenerally, the size is small; in the case of the high-speed flutter model, kbGenerally, it is larger;
kbthe physical meaning is the ratio of the beam mass in the flutter model to the total mass of the model;
step five, calculating the overweight ratio n of the flutter modelG:
nG=nsndssnBE/kb ……………………………[11]
Overweight ratio n of flutter modelGThe control is carried out in the following range as much as possible:
0.5≤nG≤2.0 ……………………………[12]
if it is out of this range, the formula [7 ] can be adjusted]To the formula [10]Obtaining the overweight ratio n of the flutter model meeting the requirementsG;
If 0.9. ltoreq. nGIf n is less than or equal to 1.1, directly taking nG=1。
Examples
Taking a low-speed and high-speed flutter model of a certain airplane wing and vertical fin as an example, calculating and verifying:
(1) calculating the equivalent mass m of the reference beam frame according to the reference beam frame modelBEMass m of the bonded memberACalculating the equivalent mass ratio n of the reference beam frameBE;
(2) Calculating the coefficient n of the design scale according to the basic scale of the flutter model designdss;
(3) Specifying a coefficient of cross-sectional shape ns;
(4) Mass ratio k of flutter mode beamb;
(5) Calculating the overweight ratio n of the flutter modelG。
Table 1 shows the corresponding data results, and the units of the data are dimensionless. Because the method of the invention uses a reverse design method to replace the prior trial and error method, when the model is not designed, the overweight ratio of the flutter model can be obtained only according to the corresponding design parameters, the calculation efficiency is greatly improved, the work which needs to be finished for several days in the past can be finished, and only several minutes are needed at present.
TABLE 1 calculation data and result of overweight ratio of airplane flutter model
Claims (1)
1. A method for determining the overweight ratio of airplane flutter model features that the flutter model with reference beam frame is ensured, which includes the mass m of airplaneABeam frame section characteristic moment of inertia Ix、IyHarmonic polar moment of inertia J, high speed flutter model design basic scale speed-pressure ratio etaqAnd density ratio ηρOr the basic scale speed ratio eta of the low-speed flutter modelVAnd density ratio ηρOn the basis of the above-mentioned parameters, determining overweight ratio n of flutter model designG(ii) a Order toUnifying high and low speed flutter models; the method is characterized by comprising the following specific steps of:
step one, calculating the equivalent mass ratio n of the reference beam frame of the airplaneBEThe calculation method is as follows:
1.1, calculating to obtain the equivalent mass m of the reference beam frameBE:
Moment of inertia I according to beam frame section characteristicsx、IyAnd polar inertia moment J, calculating the corresponding rectangular section size with the lug, and further obtaining a reference beam frame model with the section size; giving a beam frame material density to obtain the equivalent mass m of the reference beam frameBE(ii) a From the section characteristic moment of inertia Ix、IyAnd the polar inertia moment J, the method for obtaining the rectangular section size with the lug is as follows:
1.1.1 for a rectangular cross-section, the iterative formula for the dimensions of the rectangular cross-section is:
wherein,
first, alpha is set to alpha00.96, with formula [1]B/a is calculated and is according to the formula [2 ]]Calculating alpha, substituting into equation [1 ]]Through iteration, alpha is converged to a stable value finally, and then a and b can be obtained; wherein a is the half width of the rectangle calculated in the step; b is the half height of the rectangle calculated in this step;
1.1.2 for rectangular cross-section with ear, let the ear thickness
t=0.2b……………………………[3]
1.1.3 adjusting J
J1=J[1-(t/2b)2]……………………………[4]
1.1.4 according to the target value IxAnd J1And repeating the step 1.1.1 to obtain a1And b1(ii) a Wherein a is1Is the half width of the rectangle calculated in this step; b1Is the half height of the rectangle calculated in this step;
1.1.5 according to the target value IyAnd a1、b1And t value, obtaining the half width l of the section with the lug;
finally obtaining the cross section size 2a of the rectangular beam with the lug plates1、2b1T and 2 l;
1.2, calculating the equivalent mass ratio n of the reference beam frameBE:
nBE=mBE/mA……………………………[6]
Step two, calculating a coefficient n of a design scaledss:
Step three, specifying a section shape coefficient ns:
nsThe designation of (c) is divided into two cases:
if the actually adopted main beam section form is also a rectangular beam section form with lugs
ns=1……………………………[8]
If the actually adopted section form of the main beam is a rectangular thin-walled beam with lugs, the section form of the main beam is as follows
0.3≤ns≤1……………………………[9]
In this case nsThe specific numerical value can be obtained according to specific conditions, and the specific method is obtained by adjusting the thickness parameter of the section thin wall;
step four, specifying the mass ratio k of the flutter model beamb:
0.25≤kb≤0.4……………………………[10]
Wherein,
in the case of the low-speed flutter model, kbIs small; in the case of the high-speed flutter model, kb2, partial enlargement;
step five, calculating the overweight ratio n of the flutter modelG:
nG=nsndssnBE/kb……………………………[11]
Overweight ratio n of flutter modelGThe control is carried out in the following range as much as possible:
0.5≤nG≤2.0……………………………[12]
if it is out of this range, the formula [7 ] can be adjusted]To the formula [10]Obtaining the overweight ratio n of the flutter model meeting the requirementsG;
If it is not0.9≤nGIf n is less than or equal to 1.1, directly taking nG=1。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210451712.6A CN102968525B (en) | 2012-11-12 | 2012-11-12 | Method for determining overweight ratio of plane flutter model |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210451712.6A CN102968525B (en) | 2012-11-12 | 2012-11-12 | Method for determining overweight ratio of plane flutter model |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102968525A CN102968525A (en) | 2013-03-13 |
CN102968525B true CN102968525B (en) | 2015-04-29 |
Family
ID=47798662
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210451712.6A Active CN102968525B (en) | 2012-11-12 | 2012-11-12 | Method for determining overweight ratio of plane flutter model |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102968525B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104143014B (en) * | 2014-04-17 | 2017-06-16 | 中国航空工业集团公司沈阳飞机设计研究所 | One kind is without independent counterweight space tubular article flutter model method for designing |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080296432A1 (en) * | 2007-06-04 | 2008-12-04 | Kernkamp Industries Corporation | Method of extending length of aircraft to increase interior space |
CN101612996A (en) * | 2008-06-26 | 2009-12-30 | 中国航空工业第一集团公司沈阳飞机设计研究所 | A kind of plate muscle construction design method |
CN102622479B (en) * | 2012-03-02 | 2013-09-18 | 浙江大学 | Reverse engineering computer-aided design (CAD) modeling method based on three-dimensional sketch |
CN102661730B (en) * | 2012-05-11 | 2014-11-12 | 中国航空工业集团公司西安飞机设计研究所 | Method for determining hollow girder cross-section size of airplane high-speed vibration girder model |
CN102682168B (en) * | 2012-05-11 | 2014-03-19 | 中国航空工业集团公司西安飞机设计研究所 | Method for determining dimension of rectangular beam section of low-speed flutter main beam model of airplane |
-
2012
- 2012-11-12 CN CN201210451712.6A patent/CN102968525B/en active Active
Also Published As
Publication number | Publication date |
---|---|
CN102968525A (en) | 2013-03-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP2466288B1 (en) | Method of designing natural laminar flow wing for reynolds numbers equivalent to actual supersonic aircraft | |
CN108052772A (en) | A kind of geometrical non-linearity static aeroelastic analysis method based on structure reduced-order model | |
CN105122248A (en) | Method for modelling a part, in particular a turbine blade | |
CN108363843A (en) | A kind of full machine Calculate Ways of geometrical non-linearity aeroelastic effect based on structure reduced-order model | |
CN107391891A (en) | A kind of high aspect ratio wing Optimization Design based on Model Fusion method | |
CN109726437B (en) | Cabin door pneumatic load equivalent node force processing method | |
CN107742050B (en) | Entropy correction method for hybrid grid high-precision resistance prediction | |
CN113848963B (en) | Control law parameter design method of flight control system | |
CN110717216A (en) | Method for forecasting rolling response of helicopter with flexible air bag under irregular wave | |
CN102968525B (en) | Method for determining overweight ratio of plane flutter model | |
Coder et al. | Contributions to the 6th AIAA CFD drag prediction workshop using structured grid methods | |
Chen et al. | Overset Euler/Boundary-Layer solver with panel-based aerodynamic modeling for aeroelastic applications | |
CN103577649B (en) | The defining method of cargo hold floor load when transport class aircraft cargo drops | |
CN117494611A (en) | Supersonic flutter analysis method based on viscous local flow piston theory | |
CN113486512B (en) | Flutter analysis method for functional gradient variable-thickness blade model | |
Ramlan et al. | Comparison between Solidworks and Ansys CFX flow simulation on aerodynamic studies | |
CN114065423B (en) | Method for rapidly evaluating flutter of fan blade of aircraft engine | |
Desai | Relative roles of computational fluid dynamics and wind tunnel testing in the development of aircraft | |
Bozeman et al. | Mixed-Element USM3D Contributions to the 4th AIAA High Lift Prediction Workshop | |
Jones et al. | A qualitative study on the effects of mesh guideline modification for unstructured mesh generation of the NASA high lift common research model (HL-CRM) | |
Nagata et al. | Validation of new CFD tool using Non-orthogonal Octree with Boundary-fitted Layer Unstructured Grid | |
Kashyap et al. | Computational analysis of flap Camber and ground clearance in double-element inverted airfoils | |
CN102661729B (en) | Method for confirming sectional dimension of I-shaped hollow beam of high-speed fluttering model of airplane | |
CN112035947A (en) | Method for calculating wing section load with integral oil tank | |
Zhou et al. | Characteristics analysis and optimization of flying-wing vehicle structure |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |