CN112660410B - Estimation method for weight of high-aspect-ratio wing - Google Patents
Estimation method for weight of high-aspect-ratio wing Download PDFInfo
- Publication number
- CN112660410B CN112660410B CN202011610036.3A CN202011610036A CN112660410B CN 112660410 B CN112660410 B CN 112660410B CN 202011610036 A CN202011610036 A CN 202011610036A CN 112660410 B CN112660410 B CN 112660410B
- Authority
- CN
- China
- Prior art keywords
- wing
- bending moment
- cantilever beam
- weight
- length
- 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
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
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
The invention belongs to the field of optimization design of airplane structures, and relates to a method for estimating the weight of a high-aspect-ratio wing. The method comprises the following steps: equating a high aspect ratio target wing of which the weight needs to be estimated to an equal-section I-shaped cantilever beam fixedly supported at one end, and determining the length and the average height of the equal-section I-shaped cantilever beam; parameters of all the sections in the spanwise direction of the equal-section I-shaped cantilever beam are completely the same; the preset reference wing is also equivalent to a reference equal-section I-shaped cantilever beam fixedly supported at one end; determining the length and average height of a reference constant-profile I-shaped cantilever beam; obtaining equivalent total bending moment load of a target wing with a high aspect ratio and equivalent total bending moment load of a reference wing; the weight of the target wing with high aspect ratio is estimated.
Description
Technical Field
The invention belongs to the field of aircraft structure optimization design, and relates to a method for estimating the weight of a high aspect ratio wing.
Background
In the design stage of the airplane scheme, the comprehensive optimization design of aerodynamic appearance, structural arrangement and strength and rigidity of the airplane wings is required, and the wing weight is an extremely important index for multi-scheme evaluation comparison. The wing weights for different scenarios are estimated using finite element modeling and the conversion coefficients between the model and the actual structure. However, the design parameters of the wing are numerous, more possible schemes need to be compared, a finite element method is adopted for all structural schemes in the scheme stage, high estimation accuracy can be obtained, the period is long, the economy is poor, and the continuous change rule of the weight of the wing and the sensitive parameters cannot be given.
Disclosure of Invention
The purpose of the invention is as follows: for the high-aspect-ratio wing, in order to realize rapid evaluation and comparison among different schemes in the scheme stage and provide a continuous change rule between the wing weight and sensitive parameters, a method for rapidly estimating the wing weight based on the engineering beam theory is provided. The weight of the target wing can be estimated by using a simple estimation formula according to the relative relation of geometric parameters, loads and the like between the target wing and the reference wing.
The technical scheme is as follows:
a method of estimating high aspect ratio wing weight, comprising:
equating a target wing with a large aspect ratio and weight to be estimated to an I-shaped cantilever beam with a fixed and supported equal section at one end, and determining the length and the average height of the I-shaped cantilever beam with the equal section; parameters of all sections in the spanwise direction of the I-shaped cantilever beam with the equal sections are completely the same;
the preset reference wing is also equivalent to a reference equal-section I-shaped cantilever beam fixedly supported at one end; determining the length and average height of the reference equal-section I-shaped cantilever beam;
obtaining equivalent total bending moment load of a target wing with a high aspect ratio and equivalent total bending moment load of a reference wing;
and estimating the weight of the high aspect ratio target wing according to the length and the average height of the equal-section I-shaped cantilever beam, the length and the average height of the reference equal-section I-shaped cantilever beam, the equivalent total bending moment load of the high aspect ratio target wing, the equivalent total bending moment load of the reference wing, the material density and the material allowable stress level of the high aspect ratio target wing, the material density and the material allowable stress level of the reference wing and the weight of the reference wing.
Estimating the weight of the high-aspect-ratio target wing according to the length and the height of the equal-section I-shaped cantilever beam, the length and the height of the reference equal-section I-shaped cantilever beam, the equivalent total bending moment load of the high-aspect-ratio target wing, the equivalent total bending moment load of the reference wing, the material density and the allowable material stress level of the high-aspect-ratio target wing, the material density and the allowable material stress level of the reference wing and the weight of the reference wing, wherein the estimating comprises the following steps:
W＝αWRwherein
W is the weight of the target wing with high aspect ratio to be estimated, WRFor reference wing weight, α is the weight proportionality coefficient; m is the equivalent total bending moment load of the target wing with high aspect ratio, MRFor reference to wing equivalent total moment loading, αMThe equivalent total bending moment load proportionality coefficient; l is the length of the target wing with high aspect ratio, LRFor reference to the length of the wing, αLIs a length scale factor; h is the average height of the target wing with high aspect ratio, HRFor reference to the average height of the wing, alphaHIs an averaged height scale factor; rho is the density of the target wing material with high aspect ratio, rhoRFor reference to wing material density, αρIs the material density proportionality coefficient; sigma is the allowable stress level of the material of the target wing with high aspect ratio, sigmaRFor reference to permissible wing material stress level, ασAllowing for a stress proportionality coefficient for the material.
The length of the cantilever beam is the connecting line length of 50 percent chord length position of each streamlining section of the wing.
The average height of the wing is the height between the central planes of the upper and lower edge strips of the cantilever beam, namely the ratio of the volume of a space surrounded by the shape of the wing to the projection area of the wing on a chord plane.
Obtaining an equivalent total moment load of the wing, comprising:
averaging the existing aerodynamic and inertial total bending moment load envelope lines of the wings, namely drawing a curve of the total bending moment along with the relative span length, and then obtaining the area enclosed by the curve and a transverse shaft as the equivalent total bending moment of the wings;
or if the aerodynamic and inertial total bending moment load envelope of the wing does not exist, the maximum bending moment of the wing root section obtained through estimation is used as the equivalent total bending moment of the wing.
In the change curve, the Y axis is the total bending moment, and the X axis is the relative extension.
The estimated reference wings need to select the wings of one reference airplane with similar overall layout and approximately equivalent takeoff weight.
A computer-readable storage medium, on which a computer program is stored which, when being executed by a processor, carries out the above-mentioned method.
Has the advantages that: the method of the invention adopts a simple estimation formula, can quickly estimate the structural weight of the wing with a large aspect ratio, shortens the scheme evaluation period, saves the design cost, and can provide the continuous change rule between the weight of the wing and the sensitive parameter.
Drawings
FIG. 1 is a schematic view of a high aspect ratio wing;
FIG. 2 is a schematic view of an iso-section I-shaped cantilever;
FIG. 3 is a graph of total bending moment as a function of relative span length.
Detailed Description
The invention provides a method for estimating the weight of a wing with a high aspect ratio, which comprises the following steps:
(1) The I-shaped section cantilever beam is simplified. The target wing with high aspect ratio and needing to estimate the weight is shown in figure 1, and the equivalent constant-section I-shaped cantilever beam with one end fixedly supported is shown in figure 2 (all section parameters in the spanwise direction are identical). Assuming equal areas of the upper and lower edge strips, the area of the web is negligible; the length of the cantilever beam is the connecting line length L of 50 percent chord length position of each streamlining section of the wing; the height between the respective central planes of the upper and lower edge strips of the cantilever beam is the average height H of the wing, namely the ratio of the volume of a space surrounded by the shape of the wing to the projection area of the wing on a chord plane;
(2) Equivalent total moment loading of the wing. The aerodynamic and inertial total bending moment load envelope of the wing is averaged, that is, a curve of the total bending moment along with the relative span length is drawn (the total bending moment is the Y axis, and the relative span length is the X axis), and then the area enclosed by the curve and the X axis is calculated and obtained as M shown in fig. 3. When there is no accurate relevant data, the maximum bending moment of the wing root profile obtained by estimation can be used as M for subsequent calculation and analysis.
(3) And (4) material performance. The density of the target wing structural material to be estimated is ρ and the allowable stress level is σ.
(4) And selecting a reference wing. Taking the wing of a reference airplane with similar overall layout and approximately equivalent takeoff weight as an estimated reference wing, and knowing the structural weight and relevant geometric physical parameters (wing length L) of the reference wingRAnd average height HRThe area of the curve of the homogenized total bending moment or the maximum bending moment M of the rootRMaterial density ρRAllowable stress level of material sigmaR)。
(5) The weight of the target wing is estimated. And (3) taking the equivalent I-shaped cantilever beam with the equal section as an object, neglecting the structural bending rigidity and the structural weight of the web plate, and assuming that the gradients of the normal stresses of the upper edge strip and the lower edge strip of the I-shaped cantilever beam in the height direction are zero and the stress levels are equivalent, estimating the weight of the wing.
W＝αWRWherein α = αMαLαρ/(αHασ)，αM＝M/MR，αL＝L/LR，αH＝H/HR，αρ＝ρ/ρR，ασ＝σ/σR。
Claims (3)
1. A method of estimating the weight of a high aspect ratio wing, comprising:
equating a target wing with a large aspect ratio and weight to be estimated to an I-shaped cantilever beam with a fixed and supported equal section at one end, and determining the length and the average height of the I-shaped cantilever beam with the equal section; parameters of all the sections in the spanwise direction of the equal-section I-shaped cantilever beam are completely the same;
the preset reference wing is also equivalent to a reference equal-section I-shaped cantilever beam fixedly supported at one end;
determining the length and average height of a reference constant-profile I-shaped cantilever beam;
obtaining equivalent total bending moment load of a target wing with a high aspect ratio and equivalent total bending moment load of a reference wing; estimating the weight of the high aspect ratio target wing according to the length and the average height of the equal-section I-shaped cantilever beam, the length and the average height of the reference equal-section I-shaped cantilever beam, the equivalent total bending moment load of the high aspect ratio target wing, the equivalent total bending moment load of the reference wing, the material density and the material allowable stress level of the high aspect ratio target wing, the material density and the material allowable stress level of the reference wing and the weight of the reference wing;
the method comprises the following steps of estimating the weight of a high-aspect-ratio target wing according to the length and the height of an equal-section I-shaped cantilever beam, the length and the height of a reference equal-section I-shaped cantilever beam, the equivalent total bending moment load of a high-aspect-ratio target wing, the equivalent total bending moment load of a reference wing, the material density and the allowable material stress level of the high-aspect-ratio target wing, the material density and the allowable material stress level of the reference wing and the weight of the reference wing, and comprises the following steps:
W＝αWRwherein
Wherein W is the weight of the target wing with high aspect ratio to be estimated, WRFor reference wing weight, α is the weight proportionality coefficient; m is the equivalent total bending moment load of the target wing with high aspect ratio, MRFor reference of wing equivalent total moment loads, alphaMThe equivalent total bending moment load proportionality coefficient; l is the length of the target wing with high aspect ratio, LRFor reference to the length of the wing, αLIs a length scale factor; h is the average height of the target wing with high aspect ratio, HRFor reference to the average height of the wing, alphaHIs an averaged height scale factor; ρ is the high aspect ratio target wing material density, ρRFor reference to wing material density, αρIs the material density proportionality coefficient; sigma is the allowable stress level of the material of the target wing with high aspect ratio, sigmaRFor reference to permissible stress level of wing material, alphaσAllowing a stress proportionality coefficient for the material;
the length of the cantilever beam is the connecting length of 50 percent chord length position of each downwind section of the wing;
the average height of the wing is the height between the respective central planes of the upper and lower edge strips of the cantilever beam, namely the ratio of the volume of a space surrounded by the shape of the wing to the projection area of the wing on a chord plane;
wherein, obtain wing equivalent total moment load, include: averaging the existing aerodynamic and inertial total bending moment load envelope of the wing, namely drawing a curve of the total bending moment along with the relative span length, and then obtaining the area enclosed by the curve and a transverse shaft as the equivalent total bending moment of the wing; or if the aerodynamic and inertial total bending moment load envelope line of the wing does not exist, the maximum bending moment of the wing root section obtained through estimation is used as the equivalent total bending moment of the wing;
the estimated reference wings need to select the wings of one reference airplane with similar overall layout and approximately equivalent takeoff weight.
2. The method of claim 1, wherein the variation curve has a total bending moment on the Y-axis and a relative span on the X-axis.
3. A computer-readable storage medium, on which a computer program is stored, which, when being executed by a processor, carries out the method according to claim 1 or 2.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011610036.3A CN112660410B (en) | 2020-12-29 | 2020-12-29 | Estimation method for weight of high-aspect-ratio wing |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011610036.3A CN112660410B (en) | 2020-12-29 | 2020-12-29 | Estimation method for weight of high-aspect-ratio wing |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112660410A CN112660410A (en) | 2021-04-16 |
| CN112660410B true CN112660410B (en) | 2022-11-01 |
Family
ID=75411063
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011610036.3A Active CN112660410B (en) | 2020-12-29 | 2020-12-29 | Estimation method for weight of high-aspect-ratio wing |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112660410B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113761667B (en) * | 2021-09-26 | 2022-04-26 | 中国空气动力研究与发展中心高速空气动力研究所 | Empty weight estimation method for supersonic passenger plane |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| FR2883967A1 (en) * | 2005-04-04 | 2006-10-06 | Messier Bugatti Sa | Aircraft`s e.g. freight type wide-bodied airliner, weight and position of center of gravity determining method, involves measuring electrical or magnetic characteristics of unit to evaluate characteristic quantity and weight of aircraft |
| CN104029826A (en) * | 2014-06-24 | 2014-09-10 | 中国飞机强度研究所 | Taring computing method for airplane structural strength test |
| CN108100221A (en) * | 2017-11-22 | 2018-06-01 | 中国航空工业集团公司西安飞机设计研究所 | A kind of aircraft wing nonlinear analysis method |
| CN110334427A (en) * | 2019-06-27 | 2019-10-15 | 南京航空航天大学 | An Analysis Method of Aircraft Wing Structure Based on Engineering Beam Theory for TBW Layout |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8577657B2 (en) * | 2008-12-29 | 2013-11-05 | Airbus Operations S.L. | Methods for estimating the effect of variations of design variables on the weight of an aircraft component |
| CN110990957A (en) * | 2019-12-19 | 2020-04-10 | 中国航空工业集团公司沈阳飞机设计研究所 | Wing structure weight analysis method in aircraft overload performance design |
-
2020
- 2020-12-29 CN CN202011610036.3A patent/CN112660410B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| FR2883967A1 (en) * | 2005-04-04 | 2006-10-06 | Messier Bugatti Sa | Aircraft`s e.g. freight type wide-bodied airliner, weight and position of center of gravity determining method, involves measuring electrical or magnetic characteristics of unit to evaluate characteristic quantity and weight of aircraft |
| CN104029826A (en) * | 2014-06-24 | 2014-09-10 | 中国飞机强度研究所 | Taring computing method for airplane structural strength test |
| CN108100221A (en) * | 2017-11-22 | 2018-06-01 | 中国航空工业集团公司西安飞机设计研究所 | A kind of aircraft wing nonlinear analysis method |
| CN110334427A (en) * | 2019-06-27 | 2019-10-15 | 南京航空航天大学 | An Analysis Method of Aircraft Wing Structure Based on Engineering Beam Theory for TBW Layout |
Non-Patent Citations (2)
| Title |
|---|
| 机翼结构重量预测的多学科分析优化方法;余雄庆等;《航空学报》;20150818;第37卷(第01期);第235-243页 * |
| 高空长航时无人机结构重量估算方法研究;张毅等;《航空计算技术》;20070731;第37卷(第04期);第37-40页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112660410A (en) | 2021-04-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Bae et al. | Aerodynamic and static aeroelastic characteristics of a variable-span morphing wing | |
| CN105975706B (en) | A kind of conceptual level wing method for parameter estimation | |
| CN110162822B (en) | Time domain fast unsteady aerodynamic force calculation method of coupling structure mode | |
| CN111274648B (en) | Distributed flight load design method for civil aircraft leading edge flap | |
| CN114858407A (en) | Static aeroelasticity wind tunnel test wing model and composite material skin | |
| CN112660410B (en) | Estimation method for weight of high-aspect-ratio wing | |
| Xie et al. | Linearization method of nonlinear aeroelastic stability for complete aircraft with high-aspect-ratio wings | |
| US9026406B2 (en) | Method, apparatus and computer program product for determining the strain induced at a selected point in a stiffened panel structure in response to a load, taking into account one or more out of plane (OOP) effects | |
| Suleman et al. | Non-linear aeroelastic analysis in the time domain of high-aspect-ratio wings: Effect of chord and taper-ratio variation | |
| CN115374543B (en) | An aerodynamic/structural multidisciplinary design optimization method for a Lambda wing | |
| Zhao et al. | Nonlinear aeroelastic analysis of SUGAR truss-braced wing (TBW) wind-tunnel model (WTM) under in-plane loads | |
| CN105550383B (en) | A kind of design method of unsteady aerodynamic force measurement pilot system | |
| Gil et al. | Prediction of the non-linear aeroelastic behavior of a cantilever flat plate and equivalent 2D model | |
| CN116894353B (en) | Estimation method for rapidly obtaining wake vortex parameters of aircraft | |
| Sitaraman et al. | Coupled CFD/CSD prediction of rotor aerodynamic and structural dynamic loads for three critical flight conditions | |
| CN118088373A (en) | Wind turbine blade integrated design method and system based on isogeometric analysis method | |
| CN116842781A (en) | Quality inspection method for finite element model of horizontal tail structure | |
| CN114638129A (en) | Climbing capacity data obtaining method based on engine state change | |
| Previtali et al. | Performance evaluation of a 3D morphing wing and comparison with a conventional wing | |
| Cardoso et al. | Performance evaluation of a morphing joined wing aircraft configuration | |
| CN114048543B (en) | Wing box section longitudinal structure parameter definition method for associating weight and strength | |
| Diederich et al. | Charts and approximate formulas for the estimation of aeroelastic effects on the loading of swept and unswept wings | |
| CN119026432B (en) | Structural characteristic rapid equivalent analysis method for integrated coupling wing configuration design | |
| CN118532301B (en) | A wind turbine blade scaling criterion construction and testing method based on mass relationship | |
| CN114117950B (en) | Flutter judgment method for shuttle aircraft based on acting principle |
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 |