CN117521379A - Equivalent calculation method for system installation strength under amphibious aircraft vibration environment - Google Patents
Equivalent calculation method for system installation strength under amphibious aircraft vibration environment Download PDFInfo
- Publication number
- CN117521379A CN117521379A CN202311521080.0A CN202311521080A CN117521379A CN 117521379 A CN117521379 A CN 117521379A CN 202311521080 A CN202311521080 A CN 202311521080A CN 117521379 A CN117521379 A CN 117521379A
- Authority
- CN
- China
- Prior art keywords
- acceleration
- fatigue
- mean square
- root mean
- system installation
- 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
- 238000009434 installation Methods 0.000 title claims abstract description 38
- 238000004364 calculation method Methods 0.000 title claims abstract description 26
- 230000001133 acceleration Effects 0.000 claims abstract description 50
- 238000001228 spectrum Methods 0.000 claims abstract description 22
- 238000000034 method Methods 0.000 claims abstract description 20
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 10
- 230000009194 climbing Effects 0.000 claims description 5
- 230000007774 longterm Effects 0.000 claims description 3
- 238000012549 training Methods 0.000 claims description 3
- 238000013461 design Methods 0.000 abstract description 7
- 238000013016 damping Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000033001 locomotion Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 239000000446 fuel Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
Classifications
-
- 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
- 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
- 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
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Automation & Control Theory (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
The invention belongs to the technical field of aircraft strength calculation, and relates to a system installation strength equivalent calculation method under an amphibious aircraft vibration environment. Dividing the vibration environment spectrum into independent phase spectrums at different rotating speeds to calculate the response of system installation, and obtaining the duration and the number of frames occupied by the rotating speeds of the propellers at each task stage of ground shutdown, ground/water surface taxi, take-off & climb, cruising and the like by statistics, and equivalently converting the duration and the number of frames into the duration related to the rotating speeds; and converting the fatigue S-N curve into parameter equations such as root mean square acceleration, fatigue stress factors and the like, calculating the root mean square acceleration in the total flight time according to an equivalent fatigue damage criterion, determining the limiting acceleration of system installation based on Gaussian distribution, comparing with maneuver overload, and selecting a value with severe acceleration for calculation. The safety and reliability of the aircraft in the running process are greatly improved; design level and efficiency can be effectively improved.
Description
Technical Field
The invention belongs to the technical field of aircraft strength calculation, and relates to a system installation strength equivalent calculation method under an amphibious aircraft vibration environment.
Background
The amphibious aircraft power plant adopts a turbo-propeller engine and a 6-blade propeller. In the normal operation process of the aircraft, the vibration of the machine body structure is mainly caused by the unbalance of the mass of the propeller blade, the unbalance of aerodynamic force caused by the manufacturing error of the blade or the installation inclination, and the vibration is transmitted to the aircraft from the rotating shaft through the engine bracket, wherein the transmission path is as follows: shaft-engine-wing-fuselage. In turn, movement of the body structure excites vibration of the fuel, hydraulic, pneumatic and other system devices coupled thereto. Therefore, it is necessary to ensure that the installation of the system devices is able to withstand the vibration environment encountered during service. Currently, the calculation of system installation only considers motor overload in the running process of the aircraft, and the system installation strength in the vibration environment is not fully considered due to the lack of a proper method.
Disclosure of Invention
Object of the Invention
The invention aims to provide a system installation strength equivalent calculation method under an amphibious aircraft vibration environment, which is used for guiding the intensity calculation of the system installed under the vibration environment.
Technical proposal
A system installation intensity equivalent calculation method under amphibious aircraft vibration environment divides vibration environment spectrum into independent phase spectrum under different rotation speeds to calculate system installation response, obtains the time and the number of frames occupied by the rotation speed of a propeller in each task phase of ground shutdown, ground/water surface coasting, take-off & climbing, cruising and the like by statistics, and converts the time and the number of frames into time related to the rotation speed in an equivalent way; and converting the fatigue S-N curve into parameter equations such as root mean square acceleration, fatigue stress factors and the like, calculating the root mean square acceleration in the total flight time according to an equivalent fatigue damage criterion, determining the limiting acceleration of system installation based on Gaussian distribution, comparing with maneuver overload, and selecting a value with severe acceleration for calculation.
Further, the equivalent calculation method for the system installation strength under the amphibious aircraft vibration environment comprises the following steps:
step 1, dividing a full envelope spectrum of a vibration environment according to the rotating speed of an airplane propeller, and calculating the response of system installation;
step 2, obtaining the duration and the number of times of the rotating speed of the propeller in each task stage of ground shutdown, take-off, climb, cruising and the like by using statistics, and equivalently converting the duration and the number of times into duration related to the rotating speed;
step 3, converting the fatigue S-N curve into a related equation of parameters such as root mean square acceleration, fatigue stress factors and the like;
step 4, calculating root mean square acceleration in the total flight time according to an equivalent fatigue damage criterion;
step 5, determining the limit acceleration of system installation based on Gaussian distribution;
and 6, comparing the maneuver overload and selecting a value with a severe acceleration for calculation.
Further, in the step 1, dividing the full envelope spectrum of the vibration environment into independent spectrums at different rotating speed stages according to the rotating speed of the aircraft propeller;
further, in the step 2, the time length and the number of frames of each stage of long-term flight, short-term flight, training flight and task flight are counted according to the fatigue typical task;
further, in the step 2, the duration of each stage is calculated according to ground shutdown, ground sliding, water surface sliding, take-off & climbing, high-level cruising, low-level cruising, descending & approaching, water surface sliding, ground sliding and ground shutdown;
further, in step 3, the fatigue S-N curve equation is:
S m ·N=C
wherein S is a fatigue stress peak; n is the load cycle number before failure; c is the level of damage required for failure to occur; m is fatigue stress index;
further, in step 3, it is assumed that the number of failure cycles is proportional to the failure duration, and the fatigue stress peak is proportional to the root mean square acceleration;
further, in step 4, the root mean square acceleration is calculated as:
in G_RMS flight Root mean square acceleration for total time of flight; G_RMS alt The root mean square acceleration is the cruising stage; G_RMS take_off For take-off&Root mean square acceleration at climbing stage; t (T) alt Is equivalent cruising duration; t (T) take_off Is equivalent take-off duration; t (T) flight Is the total flight time; m is the fatigue stress index.
The beneficial effects of this application lie in:
the invention provides a system installation strength equivalent calculation method under an amphibious aircraft vibration environment, aiming at random vibration environments suffered by amphibious aircraft system installation and related systems in the ground, water surface and air operation processes, the fatigue equivalent method is adopted to convert the vibration of each stage of aircraft operation into root mean square acceleration under the total flight time and is used for strength check, the situation that the conventional aircraft design system installation only considers static load and does not consider dynamic load such as vibration environment is broken through, and the safety and reliability in the aircraft operation process are greatly improved; the method solves the problem that the intensity calculation method of the system installation design is deficient in a vibration environment, the calculation process is simple and easy to understand, the flow is clear, the method can be used for guiding the design and calculation of the installation intensity of the aircraft system, and the design level and efficiency can be effectively improved.
Drawings
FIG. 1 is a schematic diagram of a system installation;
wherein: 1-system support, 2-hydraulic hose and 3-hydraulic hard tube
FIG. 2 is a typical vibration environment spectrum in the present invention.
FIG. 3 is a Gaussian graph of the present invention.
Fig. 4 is a schematic diagram of a system installation calculation result.
Detailed Description
The invention is further described below with reference to examples. The following description is of some, but not all embodiments of the invention. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
Fig. 1 shows a typical installation of an aircraft system, with hydraulic hoses and hydraulic hard tubes being fastened to the aircraft structure by means of system brackets.
The invention relates to a system installation strength equivalent calculation method under an amphibious aircraft vibration environment, which comprises the following specific implementation processes:
step one, preparing parameters including the rotating speed of the propeller, the full envelope spectrum of the vibration environment, damping coefficients, fatigue stress factors and the like.
And step two, taking the rotating speeds of different propellers as central frequencies, taking a narrow frequency band as a range of +/-5% of the central frequencies, and keeping a spectrum peak value unchanged to form a vibration environment spectrum at each rotating speed.
Step three, calculating the natural frequency f of the system installation structure n 。
Step four, calculating a displacement transfer function Q:
in the formula (1), f n The natural frequency of the system installation structure is in Hz; f (f) 0 Is random motion and forcing frequency, in Hz; alpha is the damping coefficient.
Step five, taking the vibration environment spectrum in the step two as an acceleration spectrum density input, solving to obtain an acceleration spectrum density response, wherein a calculation formula is as follows:
ASD OUT =Q 2 ·ASD IN ..................................................(2)
in the formula (2), ASD OUT In g, the unit of the acceleration spectral density response 2 /Hz;ASD IN In g as a vibration environmental spectrum 2 /Hz。
Step six, calculating step five to obtain square root acceleration G_RMS of the area under the acceleration spectral density response curve:
and step seven, counting the time length of each stage of long-term flight, short-term flight, training flight and mission flight and the number of frames according to the fatigue typical mission.
And step eight, converting the statistical result in the step seven into the duration occupied by flight stages such as ground stop, ground sliding, water surface sliding, take-off and climbing, high-level cruising, low-level cruising, descending and approaching, water surface sliding, ground stop and the like.
Step nine, accumulating the duration under the same rotation speed in the step eight to obtain the equivalent task duration T under each rotation speed i 。
Step ten, based on the assumption that the number of failure cycles N is proportional to the failure duration, and that the fatigue stress peak S is proportional to the root mean square acceleration g_rms, the S-N curve can be written as:
G_RMS m ·T=K..............................................(4)
in the formula (4), m is fatigue stress index; t is the task duration; k is the proportion of damage level required for failure to occur.
Step eleven, calculating root mean square acceleration in the total flight time according to an equivalent fatigue damage criterion:
G_RMS flight m ·T flight =∑G_RMS i m ·T i ...................................(5)
in formula (5), G_RMS flight Root mean square acceleration as total time of flight; t (T) flight Is the total flight time; G_RMS i Is the square root acceleration at each rotational speed.
Step twelve, according to the Gaussian distribution curve of the attached figure 3, the ratio of the transient acceleration to the root mean square acceleration is 68.3% of time within +/-1, the curve also shows that 99.73% of time is within +/-3, and the limiting vibration acceleration is 3 times of the root mean square acceleration, so that the system installation acceleration can be basically covered within 100% of time.
Step thirteen, comparing with the motor overload, and selecting an acceleration severity value to calculate aiming at the system bracket in the figure 1.
According to the equivalent calculation method for the system installation strength under the amphibious aircraft vibration environment, the vibration environment spectrum is divided into independent phase spectrums under different rotation speeds to calculate the response of system installation, the time and the number of frames occupied by the rotation speeds of the propellers in each task phase of ground shutdown, take-off, climb, cruise and the like are obtained through statistics, and the equivalent conversion is carried out to the time related to the rotation speeds; and converting the fatigue S-N curve into a related equation of parameters such as root mean square acceleration, fatigue stress factors and the like, calculating the root mean square acceleration in the total flight time according to an equivalent fatigue damage criterion, determining the limit acceleration of system installation based on Gaussian distribution, comparing with maneuver overload, and selecting a harsh value of the acceleration for strength calculation. The method solves the problem that the intensity calculation method of the system installation design in the vibration environment is not available, the flow is clear, and the design level and the efficiency can be effectively improved.
It will be understood by those skilled in the art that, unless otherwise defined, all terms (including technical and scientific terms) used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs. It will be further understood that terms, such as those defined in commonly used dictionaries, should be interpreted as having a meaning that is consistent with their meaning in the context of the prior art and will not be interpreted in an idealized or overly formal sense unless expressly so defined herein. While the foregoing is directed to embodiments of the present invention, other and further details of the invention may be had by the present invention, it should be understood that the foregoing description is merely illustrative of the present invention and that no limitations are intended to the scope of the invention, except insofar as modifications, equivalents, improvements or modifications are within the spirit and principles of the invention.
Claims (8)
1. A method for equivalently calculating the system installation strength under the vibration environment of an amphibious aircraft is characterized in that the vibration environment spectrum is divided into independent phase spectrums under different rotation speeds to calculate the system installation response, and the time and the number of frames occupied by the rotation speeds of propellers in each task stage of ground shutdown, ground/water surface taxi, take-off & climb, cruising and the like are obtained by statistics and are equivalently converted into time related to the rotation speeds; and converting the fatigue S-N curve into parameter equations such as root mean square acceleration, fatigue stress factors and the like, calculating the root mean square acceleration in the total flight time according to an equivalent fatigue damage criterion, determining the limiting acceleration of system installation based on Gaussian distribution, comparing with maneuver overload, and selecting a value with severe acceleration for calculation.
2. The method of claim 1, characterized by the steps of:
step 1, dividing a full envelope spectrum of a vibration environment according to the rotating speed of an airplane propeller, and calculating the response of system installation;
step 2, obtaining the duration and the number of times of the rotating speed of the propeller in each task stage of ground shutdown, take-off, climb, cruising and the like by using statistics, and equivalently converting the duration and the number of times into duration related to the rotating speed;
step 3, converting the fatigue S-N curve into a related equation of parameters such as root mean square acceleration, fatigue stress factors and the like;
step 4, calculating root mean square acceleration in the total flight time according to an equivalent fatigue damage criterion;
step 5, determining the limit acceleration of system installation based on Gaussian distribution;
and 6, comparing the maneuver overload and selecting a value with a severe acceleration for calculation.
3. The method according to claim 2, characterized in that in step 1, the vibration environment full envelope spectrum is divided into separate spectrums at different rotational speed phases depending on the rotational speed of the aircraft propeller.
4. A method according to claim 3, wherein in step 2, the long-term flight, short-term flight, training flight, the duration of each stage of the mission and the number of frames are counted according to the fatigue-typical mission.
5. The method of claim 4, wherein in step 2, the phase durations are calculated in terms of ground stop, ground taxi, water taxi, take-off & climb, high level cruise, low level cruise, descent & approach, water taxi, ground stop.
6. The method of claim 5, wherein in step 3, the fatigue S-N curve equation is:
S m ·N=C
wherein S is a fatigue stress peak; n is the load cycle number before failure; c is the level of damage required for failure to occur; m is the fatigue stress index.
7. The method of claim 6, wherein in step 3, the number of failure cycles is assumed to be proportional to the length of the failure time and the fatigue stress peak is assumed to be proportional to the root mean square acceleration.
8. The method of claim 7, wherein in step 4, the root mean square acceleration is calculated as:
in G_RMS flight Root mean square acceleration for total time of flight; G_RMS alt The root mean square acceleration is the cruising stage; G_RMS take_off For take-off&Root mean square acceleration at climbing stage; t (T) alt Is equivalent cruising duration; t (T) take_off Is equivalent take-off duration; t (T) flight Is the total flight time; m is the fatigue stress index.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311521080.0A CN117521379A (en) | 2023-11-15 | 2023-11-15 | Equivalent calculation method for system installation strength under amphibious aircraft vibration environment |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311521080.0A CN117521379A (en) | 2023-11-15 | 2023-11-15 | Equivalent calculation method for system installation strength under amphibious aircraft vibration environment |
Publications (1)
Publication Number | Publication Date |
---|---|
CN117521379A true CN117521379A (en) | 2024-02-06 |
Family
ID=89741450
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311521080.0A Pending CN117521379A (en) | 2023-11-15 | 2023-11-15 | Equivalent calculation method for system installation strength under amphibious aircraft vibration environment |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117521379A (en) |
-
2023
- 2023-11-15 CN CN202311521080.0A patent/CN117521379A/en active Pending
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Igra | Research and development for shrouded wind turbines | |
Rauch | Design study of an air pump and integral lift engine ALF-504 using the Lycoming 502 core | |
Hütter | Optimum wind energy conversion systems | |
US20210108572A1 (en) | Advance ratio for single unducted rotor engine | |
US4350898A (en) | Lighter than air wind energy conversion system utilizing an external radial disk diffuser | |
Abrego et al. | Performance study of a ducted fan system | |
US4350899A (en) | Lighter than air wind energy conversion system utilizing a rearwardly mounted internal radial disk diffuser | |
US4486146A (en) | Aircraft propulsion means | |
Igra | The shrouded aerogenerator | |
Kentfield | Fundamentals/wind-driven water | |
WO2016010450A1 (en) | Dual rotor wind power assembly (variants) | |
CN103661925B (en) | A kind of automatic Variable Pitch screw based on composite determines method | |
Mort et al. | Aerodynamic characteristics of a 4-foot-diameter ducted fan mounted on the tip of a semispan wing | |
CN117521379A (en) | Equivalent calculation method for system installation strength under amphibious aircraft vibration environment | |
Platt Jr | Static tests of a shrouded and an unshrouded propeller | |
Webster et al. | The effect of rotor blade sweep on the performance of the Wells turbine | |
CN111056036A (en) | Rapid iteration generation method for high-altitude propeller | |
McCloud et al. | Full Scale Wind-tunnel Tests of Blowing Boundary-layer Control Applied to a Helicopter Rotor | |
Vogeley | Sound-level measurements of a light airplane modified to reduce noise reaching the ground | |
CN111271179B (en) | Power performance test method for ram air turbine | |
Balaji et al. | Effect of angle of attack on pressure distribution of NACA5520 airfoil blade | |
RU2792994C1 (en) | Propeller with external annular frame and tension blades | |
Dunham et al. | Low-speed wind-tunnel tests of single-and counter-rotation propellers | |
Mort | Summary of large-scale test of ducted fans | |
Patterson Jr | A wind tunnel investigation of jet-wake effect of a high-bypass engine on wing-nacelle interference drag of a subsonic transport |
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 |