CN103912315B - Structural dynamics design method of rotor of aerial engine - Google Patents
Structural dynamics design method of rotor of aerial engine Download PDFInfo
- Publication number
- CN103912315B CN103912315B CN201410146849.XA CN201410146849A CN103912315B CN 103912315 B CN103912315 B CN 103912315B CN 201410146849 A CN201410146849 A CN 201410146849A CN 103912315 B CN103912315 B CN 103912315B
- Authority
- CN
- China
- Prior art keywords
- rotor
- support
- foundation system
- omega
- mode
- 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
Landscapes
- Manufacture Of Motors, Generators (AREA)
Abstract
Disclosed is a structural dynamics design method of the rotor of an aerial engine. Based on the characteristics of the working status of a high propulsion ratio aerial engine, the structural dynamics design method of the rotor of the aerial engine presents the concept of thermal modal to adapt to structural dynamics design of the rotor of the aerial engine under variable working conditions. The structural dynamics design method of the rotor of the aerial engine is characterized by optimizing the parameters of the rotor and bearings, enabling the thermal modal of a rotor system to avoid the modal of the rotor under the absolute rigidity of the bearings, meanwhile, giving play to the effects of a damper and enabling the rotor system to meet the requirements of vibration standards in the thermal modal. The structural dynamics design method of the rotor of the aerial engine comprises taking the margin between the modal of the rotor with elastic bearings and the modal of the rotor under the absolute rigidity of the bearings as optimized parameters, meanwhile giving play to the damping effects of the damper to meet the dynamics requirements of the rotor system. The structural dynamics design method of the rotor of the aerial engine changes the traditional design concept of design first and checking calculation later and performs active design of the dynamics characteristics of the rotor to achieve circulating design, thereby improving the working efficiency and having significant engineering practical value.
Description
Technical field
The present invention relates to aeroengine field, specifically a kind of dynamic (dynamical) method of designed engines rotor structure.
Background technique
Reduce vibration, control gap, and reduction supporting structure load is the key content involved by rotor dynamics design criterion.Along with to the lifting of opportunity of combat operational requirements and the diversity of aerial mission, aeroengine variable working condition feature is further outstanding.Therefore, aeroengine rotor Structural Dynamic Design is extremely important.
The design criterion that existing aeroengine rotor Structural Dynamic Design is continued to use always " ensureing that working speed and critical speed of rotation leave enough nargin ", no matter the rotor of i.e. subcritical or overcritical operation, its working speed must keep enough nargin with critical speed of rotation, and such as 15%.Due in overcritical work, therefore damper will must be arranged to be decreased through vibratory response during critical speed of rotation.Height push away than performance motor during operation, rotor frequently crosses that single order is critical, even second order critical speed of rotation, and cannot ensure the margin requirement between working speed and critical speed of rotation, critical speed of rotation position or neighborhood are even dropped in operation point.In other words, critical speed of rotation may become working speed.Therefore, the design criterion " ensureing that working speed and critical speed of rotation leave enough nargin " no longer meets the requirement of this kind of aeroengine rotor dynamics Design.
Content about the rotor dynamics of rotor-support-foundation system in domestic and international Patents document and paper focuses mostly in the design of certain rotor experiment table and the technique study of rotor dynamic Epidemiological Analysis.Its research background depends on current version and rotor.Not by the whole design phase of rotor dynamics design through rotor-support-foundation system, its expansion is subject to a definite limitation.Be disclose one " dynamics characteristic experimental apparatus of heavy type gas turbine pull rod rotor " in the innovation and creation of CN 101329220B at publication number, the invention provides a kind of experimental setup.This experimental setup comprises rod fastening rotor and supporting vibration test system, is suitable for carrying out rotor dynamics experimental research.This invention is a kind of experimental setup of certain rotor, does not obtain the dynamic (dynamical) design method of rotor structure.Be in the innovation and creation of 201410033310.3, disclose one " aircraft engine high pressure rotor Structural Design " at application number, this invention determines the relation between high pressure rotor model design parameter and rotor oscillation characteristic, thus providing design method and criterion for the dynamics Design of high pressure rotor, the design for engine high pressure rotor has important directive significance.Rotor is thought of as stiffness rotor by this invention, is only applicable to the design of aircraft engine high pressure rotor, also needs the design method improving low pressure rotor further.
L.E.Barrett, E.J.Gunter and P.E.Allaire describes a kind of quick approximation method calculating the damping of multi rotor system optimal support in paper " Optimum bearing and supportdamping for unbalance response and stability of rotating machinery " (ISSN:0022-0825), and this optimal damping makes the minimum stability of rotor-support-foundation system unbalance response near first critical speed the highest.The work of this paper belongs to the analysis of rotor dynamics, and the method for proposition makes rotor dynamics analysis be simplified, and saves time and is convenient to application, but does not propose the concrete grammar step of rotor structure dynamics Design.Jiang Yun-fan, Liao Ming-fu, Wang Siji devises a kind of aeroengine of simulating to the supporting means turning two-spool dynamics simulations device and two kinds of high-pressure shaft according to dynamic similarity principle in paper " The Design of Counter-Rotating Dual Rotor Experimental Apparatus " (ISBN:978-7-5612-3441-9).Carry out double rotor with turning, to rotational kinematics and dynamics experiment, have studied birotary engine dynamic balancing technique.The double rotor tester designed in this paper, according to dynamic similarity principle, is a kind of passive method, and object is the dynamics of research birotary engine, instead of the method for active designs rotor.
Summary of the invention
In order to adapt to highly to push away designing requirement more dynamic (dynamical) than performance aeroengine rotor, the present invention proposes a kind of dynamic (dynamical) method of designed engines rotor structure.
Push away than performance aeroengine for height, the present invention proposes a kind of definition of rotor mode, be hot-die state, it is critical that corresponding critical speed of rotation is defined as heat.Why be defined as hot-die state, one is that rotor frequently crosses critical speed of rotation, and critical speed of rotation even may become working speed because motor run duration; Two is because under working state, and material behavior, cooperation rigidity and coupling stiffness generation significant change, cause critical speed of rotation to change in a big way, is difficult to rotating speed nargin when ensureing design; Three is because under hot-die state, the criterion of rotor structure dynamics Design is, the hot threshold response of rotor be controlled under the limits value allowed, and be no longer the rotating speed nargin deliberately ensureing to expect.
Concrete steps of the present invention are:
The first step: set up rotor dynamics model, obtains the vibratory response of rotor-support-foundation system at critical speed of rotation place and the relation of each design parameter.The relation of the vibratory response of described rotor-support-foundation system at critical speed of rotation place and each design parameter comprises the relation of rotor-support-foundation system threshold response and Structural Parameters of its Rotor and the relation of rotor-support-foundation system threshold response and rotor modal parameter.
The equation matrix expression (1) of rotor dynamics model is:
In formula (1):
q is the response of rotor-support-foundation system, and r is the amount of deflection at rotor c.g. place,
for the deflection angle at rotor c.g. place, r
b1, r
b2be respectively the amount of deflection of rotor two supportings place.
Rotor-support-foundation system is at the threshold response of supporting place
The threshold response at rotor c.g. place
Formula (17), formula (18) describe the vibratory response Q of flexibly mounted rotor-support-foundation system at critical speed of rotation place
1cr, Q
2cr, include the modal parameter of rotor-support-foundation system structural dynamic parameter and rotor-support-foundation system.Formula (17), formula (18) describe flexibly mounted rotor-support-foundation system threshold response and rotor critical speed during supporting absolute rigidity
relation.
In formula (17) and (18),
Q
1crfor the threshold response amplitude at rotor c.g. place; Q
2crfor the threshold response amplitude at rotor bearing place.
Ω
crifor rotor i-th rank critical speed of rotation.
The modal matrix that [Φ] is rotor-support-foundation system,
Be respectively the first rank and second-order damping ratios, d
cr1, d
cr2be respectively the first rank and the second-order modal damping coefficient of rotor-support-foundation system, m
cr1, I
cr2be respectively the first rank and second-order modal mass, Ω
cr1, Ω
cr2be respectively the first rank and the second-order critical speed of rotation of rotor.
when being respectively supporting rigidity, the critical speed of rotation of rotor-support-foundation system.
U
1for the amount of unbalance of rotor-support-foundation system.
Ω
crifor the critical speed of rotation of rotor, Q
2crfor rotor is in the threshold response amplitude of supporting place.
Second step: the hot-die state during work of setting rotor-support-foundation system.
According to operating rotational speed range and the performance requirement of motor, the hot-die state of setting rotor-support-foundation system.When setting described hot-die state, must meet the following conditions successively: first step mode critical speed of rotation is positioned at below slow train rotating speed, and allow in slow train position; Second-order mode critical speed of rotation and the 3rd rank mode critical speed of rotation are all positioned at below cruising speed, nargin equal 10%; Critical speed of rotation is not established between cruising speed and maximum (top) speed; When the design parameter of rotor requires to establish critical speed of rotation between the cruising speed and maximum (top) speed of rotor-support-foundation system, then quadravalence mode critical speed of rotation is arranged on maximum (top) speed, nargin 8-10%.
3rd step: the mode of rotor-support-foundation system when determining to support absolute rigidity.
If each support stiffness of rotor-support-foundation system is infinitely great, transfer matrix method or finite element method is adopted to determine the mode of rotor-support-foundation system with this understanding.
4th step: the initial value choosing support stiffness.The initial value of support stiffness is selected according to the rigidity of supporting absolute rigidity rotor-support-foundation system.Each support stiffness of rotor-support-foundation system should be the 20-50% of rigid support rotor rigidity, adopts finite element method or determines the lateral stiffness of rotor at barycenter place by formula (19):
E is the Young's modulus of material, and I is the cross section moments of inertia of rotating shaft.
Get 50% of s as the initial value of rotor-support-foundation system elastic supporting rigidity:
5th step: verify whether the quiet amount of deformation of rotor during the support stiffness initial value determined exceeds standard.
Adopt finite element method or the quiet distortion according to formula (6) calculating rotor-support-foundation system, make the quiet amount of deformation of rotor be not more than the value of this engine design requirement.
In formula, δ is quiet amount of deformation; M is whole rotor quality; G is gravity accleration; M is the quality of dish; S is rotor centroid position lateral stiffness, namely produces the power needed for unit lateral deformation.
6th step: check the mode of rotor during the support stiffness initial value for determining whether to conform to the hot-die state preset.
For selected support stiffness parameter, adopt the mode of transfer matrix method or Finite element arithmetic rotor-support-foundation system.Whether the mode of inspection rotor-support-foundation system conforms to the hot-die state preset in second step.If the mode of rotor-support-foundation system is not inconsistent with the hot-die state preset, then adjust the parameter S of support stiffness
b1and S
b2, the mode of rotor-support-foundation system is conformed to the hot-die state preset in second step.
After the mode of rotor-support-foundation system conforms to the hot-die state preset in second step, check described hot-die state and rotor-support-foundation system just to prop up whether to leave between mode the nargin of 10%.If nargin is inadequate, the adjustment parameter Sb1 of support stiffness and quality, the dimensional parameters of Sb2 or rotor, make hot-die state and rotor-support-foundation system just prop up to leave between mode the nargin of 10%.Need the quiet deformation requirements of satisfied 5th step simultaneously;
7th step: damper is set in supporting place, determine unbalance response and requirement for dynamic balance:
Introduce squeeze film damper in supporting place, calculate the unbalance response of rotor.The damping function that can play according to damper and counterbalance effect, obtain rotor residual unbalance, limitation standard.So far, the Structural Dynamic Design of aeroengine rotor is completed.
The present invention proposes the concept of aeroengine rotor hot-die state, to adapt to the design of aeroengine variable working condition lower rotor part structural dynamics.Described hot-die state refers to engine operation, and the running speed of this motor frequently can pass through the critical speed of rotation of certain or certain several rotor, the mode corresponding to critical speed of rotation of this certain or certain several rotor.Based on this, the present invention establishes method and the step of aeroengine rotor dynamics Design under hot-die state.The core of the present invention to hot-die state lower rotor part dynamics Design optimizes the parameter of rotor and supporting, the mode of rotor when making the hot-die state of rotor-support-foundation system avoid supporting absolute rigidity, play the effect of damper simultaneously, make rotor-support-foundation system under hot-die state, meet the requirement of vibration standard.
The present invention establishes yielding support rotor structure dynamic model, by theory analysis, obtains the vibratory response of rotor-support-foundation system at critical speed of rotation place and the relation of each design parameter.The relation of the vibratory response of described rotor-support-foundation system at critical speed of rotation place and each design parameter comprises the relation of rotor-support-foundation system threshold response and Structural Parameters of its Rotor and the relation of rotor-support-foundation system threshold response and rotor modal parameter.The relation of rotor mode when the flexibly mounted rotor-support-foundation system threshold response of this specification of a model simultaneously and supporting absolute rigidity.The method designed for rotor dynamics and step provide foundation and criterion.
, plural hot-die state in operating range, may be there is in aeroengine rotor complex structure.The present invention propose, utilize flexibly mounted rotor mode and supporting absolute rigidity rotor mode between nargin as Optimal Parameters, physical significance is definitely.
The present invention pushes away the feature than performance aeroengine operation status from height, the concept of hot-die state is proposed, on the basis of rotor dynamics modeling, by the vibratory response characteristic of theory analysis rotor-support-foundation system at critical speed of rotation place and the relation of each design parameter, propose the nargin between flexibly mounted rotor mode and the rotor mode of supporting absolute rigidity as Optimal Parameters, play the effectiveness in vibration suppression of damper, to reach the demanding kinetics of rotor-support-foundation system simultaneously.The rotor dynamics design method that the present invention proposes and step are a set of flow processs through whole rotor design, change traditional first to design the design philosophy checked again.Active designs is carried out, cyclic design to the dynamics of rotor, improves working efficiency, there is important engineering practical value.
Accompanying drawing explanation
Fig. 1 is rotor modeling;
Fig. 2 is the vibration shape of rotor;
Fig. 3 is the vibration shape of rotor;
Fig. 4 is the vibration shape of rotor;
Fig. 5 is the amplitude-versus-frequency curve of rotor;
Fig. 6 is the amplitude-versus-frequency curve of rotor;
Fig. 7 is flow chart of the present invention.In figure: 1. a first order mode; 2. second_mode; 3. three first order modes
Embodiment
In order to adapt to highly to push away designing requirement more dynamic (dynamical) than performance aeroengine rotor, the present embodiment proposes a set of rotor dynamics design method designed through whole engine rotor.
Push away than performance aeroengine for height, the present embodiment proposes a kind of definition of rotor mode, instant heating mode.Described hot-die state refers to engine operation, and the running speed of this motor frequently can pass through the critical speed of rotation of certain or certain several rotor, the mode corresponding to critical speed of rotation of this certain or certain several rotor.
Why be defined as hot-die state, one is that rotor frequently crosses critical speed of rotation, and critical speed of rotation even may become working speed because motor run duration; Two is because under working state, and material behavior, cooperation rigidity and coupling stiffness generation significant change, cause critical speed of rotation to change in a big way, is difficult to rotating speed nargin when ensureing design; Three is because under hot-die state, the criterion of rotor structure dynamics Design is, response limiting when rotor frequently be passed through rotor critical speed under the limits value allowed, and is no longer the rotating speed nargin deliberately ensureing to expect.
The present embodiment is for the dynamic (dynamical) method of designed engines rotor structure.With rotator model as shown in Figure 1 for design object carries out rotor dynamics design.Concrete steps comprise:
The first step: set up rotor dynamics model, obtains the vibratory response of rotor-support-foundation system at critical speed of rotation place and the relation of each design parameter.The relation of the vibratory response of described rotor-support-foundation system at critical speed of rotation place and each design parameter comprises the relation of rotor-support-foundation system threshold response and Structural Parameters of its Rotor and the relation of rotor-support-foundation system threshold response and rotor modal parameter.
The kinetic equations expression matrix form of rotor-support-foundation system is as formula (1):
In formula,
q is the response of rotor-support-foundation system, and r is the amount of deflection at rotor c.g. place,
for the deflection angle at rotor c.g. place, r
b1, r
b2be respectively the amount of deflection of rotor two supportings place.
Order
q
1it is the response at rotor c.g. place;
q
2the response at rotor bearing place, by matrix expression (1) afterwards two row can obtain:
In formula,
Formula (2) is taken back front two row of matrix expression (1), obtains:
In formula,
If rotor at the threshold response of supporting place is
substitute into equation (3) to obtain:
In formula, Ω
crifor the critical speed of rotation of rotor, Q
2crfor rotor is in the threshold response amplitude of supporting place, β
2crfor rotor is at the phase angle of supporting place threshold response.
Rotor speed is critical speed of rotation this moment, and the Section 2 on equation (4) left side is 0, then have
Introduce as down conversion
After substituting into equation (5), obtain
Can be solved by equation (7)
In formula,
When supporting absolute rigidity, if the modal matrix of rotor-support-foundation system is
In formula
for the first first order mode of rotor during rigid support;
for the second-order vibration shape of rotor during rigid support.Meanwhile, modal stiffness and modal mass is introduced:
In formula,
when being respectively rigid support, the first step mode rigidity of rotor and second-order modal stiffness.
In formula,
when being respectively rigid support, the first step mode quality of rotor and second-order modal mass, wherein subscript cr1 and cr2 represents the first rank and second-order respectively.
In formula (8), the molecule of equal sign right-hand member and denominator are with taking advantage of
namely
According to the character of determinant, and substitute into formula (10) and (11)
In formula (12), equal sign right-hand member molecule and denominator are with taking advantage of
and substitute into formula (13), and abbreviation,
In formula,
when being respectively supporting rigidity, the critical speed of rotation of rotor-support-foundation system.
Formula (14) is substituted into formula (6), namely obtains the vibration of rotor in supporting place:
Introduce modal damping matrix
In formula, when [Φ] is for bearing elastic, the modal matrix of rotor-support-foundation system, dcr1, dcr2 are respectively the first rank and the second-order modal damping coefficient of rotor-support-foundation system,
Be respectively the first rank and second-order damping ratios, m
cr1, I
cr2be respectively the first rank and second-order modal mass, Ω
cr1, Ω
cr2when being respectively support elasticity, the first rank of rotor and second-order critical speed of rotation.
Wushu (16) substitutes into formula (15), finally obtains the threshold response of rotor-support-foundation system in supporting place
In generation, returns the threshold response that formula (2) obtains rotor c.g. place
Formula (17), formula (18) describe the vibratory response Q of flexibly mounted rotor-support-foundation system at critical speed of rotation place
1cr, Q
2cr, include the modal parameter of rotor-support-foundation system structural dynamic parameter and rotor-support-foundation system.Formula (17), formula (18) describe flexibly mounted rotor-support-foundation system threshold response and rotor critical speed during supporting absolute rigidity
relation.For the design of follow-up rotor dynamics is laid a good foundation.
Second step: the hot-die state during work of setting rotor-support-foundation system.
According to operating rotational speed range and the performance requirement of motor, the hot-die state of setting rotor-support-foundation system.When setting described hot-die state, must meet the following conditions successively: first step mode critical speed of rotation is positioned at below slow train rotating speed, and allow in slow train position; Second-order mode critical speed of rotation and the 3rd rank mode critical speed of rotation are all positioned at below cruising speed, nargin equal 10%; Critical speed of rotation is not established between cruising speed and maximum (top) speed; When the design parameter of rotor requires to establish critical speed of rotation between the cruising speed and maximum (top) speed of rotor-support-foundation system, then quadravalence mode critical speed of rotation is arranged on maximum (top) speed, nargin 8-10%.
In the present embodiment, if the slow train rotating speed of rotor is 6000 revs/min, cruising speed is 13000 revs/min, maximum (top) speed 15000 revs/min.First three rank critical speed of rotation of the rotor that consideration need design, so the first step mode critical speed of rotation of rotor need be arranged on less than 6000 revs/min, second-order mode critical speed of rotation need be arranged on less than 13000 revs/min, nargin 10%, in 13000 revs/min to 15000 revs/min speed range, critical speed of rotation is not set, 3rd rank mode critical speed of rotation is arranged on more than 14000 revs/min, nargin 8-10%.
3rd step: the mode of rotor-support-foundation system when determining to support absolute rigidity.
If each support stiffness of rotor-support-foundation system is infinitely great, transfer matrix method or finite element method is adopted to determine the mode of rotor-support-foundation system with this understanding.Transfer matrix method is adopted to determine the mode of rotor-support-foundation system with this understanding in the present embodiment.
For the rotor in the present embodiment, design parameter is as follows:
Rotor span L=0.81m, shaft diameter 0.055m, coil from front fulcrum distance a=0.54m, dish quality 20kg, diameter rotary inertia I
d=0.032kgm, polar moment of inertia I
p=0.064kgm, the Young's modulus of axle is 2.06 × 10
11n/m
2., moment of area of inertia I=4.4918 × 10 of shaft section
-7m
4, rotor material density 7850kg/m
3.
Rotor is divided into 18 sections, the two rank critical speed of rotation adopting transfer matrix method to calculate rotor are respectively: 5776 revs/min, 32461 revs/min, its corresponding Mode Shape as shown in Figure 2.
4th step: the initial value choosing support stiffness.The initial value of support stiffness is selected according to the rigidity of supporting absolute rigidity rotor-support-foundation system; Each support stiffness of rotor-support-foundation system should be the 20-50% of rigid support rotor rigidity, adopts finite element method or determines the lateral stiffness of rotor at barycenter place by formula (19).
For the rotor in the present embodiment, rotor centroid place produces the power required for unit deformation, and namely rotor is at the lateral stiffness at barycenter place, is calculated as follows:
E is the Young's modulus of material, and I is the cross section moments of inertia of rotating shaft.
Get 50% of s as the initial value of rotor-support-foundation system elastic supporting rigidity:
5th step: verify whether the quiet amount of deformation of rotor during the support stiffness initial value determined exceeds standard.
Adopt finite element method or the quiet distortion according to formula (21) calculating rotor-support-foundation system, make the quiet amount of deformation of rotor be not more than the value of this engine design requirement.Require in the present embodiment that quiet amount of deformation is not more than 0.1mm.
In formula, δ is quiet amount of deformation; M is whole rotor quality; G is gravity accleration; M is the quality of dish; S is rotor centroid position lateral stiffness, namely produces the power needed for unit lateral deformation.
For the present embodiment, quiet amount of deformation
Meet the demands.
6th step: check the support stiffness initial value for determining, whether the mode of rotor conforms to the hot-die state preset.
For selected support stiffness parameter, adopt the mode of transfer matrix method or Finite element arithmetic rotor-support-foundation system; Whether the mode of inspection rotor-support-foundation system conforms to the hot-die state preset in second step; If the mode of rotor-support-foundation system is not inconsistent with the hot-die state preset, then adjust the parameter s of support stiffness
b1and s
b2, the mode of rotor-support-foundation system is conformed to the hot-die state preset in second step;
After the mode of rotor-support-foundation system conforms to the hot-die state preset in second step, check described hot-die state and rotor-support-foundation system just to prop up whether to leave between mode the nargin of 10%; If nargin is inadequate, the parameter S of adjustment support stiffness
b1and S
b2or the quality of rotor, dimensional parameters, hot-die state and rotor-support-foundation system are just propped up to leave between mode the nargin of 10%; Need the quiet deformation requirements of satisfied 5th step simultaneously;
For the present embodiment, when support stiffness is provided by formula (20), equally rotor-support-foundation system is divided into 18 sections, the three rank critical speed of rotation adopting transfer matrix method to calculate are respectively 3880 revs/min, 12133 revs/min, 29862 revs/min, and the Mode Shape of its correspondence as shown in Figure 3.
Through inspection, the second-order mode of the present embodiment rotor is not inconsistent with the hot-die state preset, and namely second-order critical speed of rotation 12133 revs/min is too high, and the first rank and Section of three mode all well meet the requirement of hot-die state.So retainer spring stiffness to be reduced to 40% of s, now the quiet amount of deformation of rotor is
meet the demands.The three rank critical speed of rotation calculating now rotor-support-foundation system are respectively 3634 revs/min, 10955 revs/min, 28610 revs/min, and through inspection, rotor-support-foundation system mode now conforms to the hot-die state preset.The Mode Shape of its correspondence as shown in Figure 4.
After reaching the requirement of hot-die state, whether inspection hot-die state has just propped up mode with rotor-support-foundation system is left enough nargin, and it is 10% that the present embodiment requires.Be respectively through the hot-die state of test design rotor and a just mode nargin of rotor-support-foundation system: 37%, 66%, meet the demands.
7th step: arrange damper in supporting place, determines unbalance response and requirement for dynamic balance.
Within the engine, in order to meet demanding kinetics, needing to apply damping at supporting place design squeeze film damper to rotor, especially needing the height repeatedly through critical speed of rotation region to push away and compare performance engine.After introducing squeeze film damper, calculate the unbalance response of rotor.The damping function that can play according to damper and counterbalance effect, finally can provide rotor residual unbalance, limitation standard.
For the present embodiment, rotor is that displacement amplitude is not more than 50 μm at the vibration amplitude limitation standard of supporting place, and apply damper at rear support place, damping constant is 500Ns/m, and rotor amount of unbalance before dynamic balancing is 1 × 10
-4kgm, is calculated by transfer matrix method equally, obtain rotor-support-foundation system two bearing places response be respectively: 1 rank: 73.0 μm, 112.2 μm; 2 rank: 86.4 μm, 47.3 μm; 3 rank: 13.63 μm, 23.2 μm, amplitude-versus-frequency curve as shown in Figure 5.Do not meet amplitude standards, need dynamic balancing be done.After dynamic balancing, the amount of unbalance of rotor is reduced to 0.3 × 10
-4kgm, now rotor-support-foundation system two bearing places response be respectively: 1 rank: 21.9 μm, 33.7 μm; 2 rank: 25.9 μm, 14.3 μm; 3 rank: 4.1 μm, 7.0 μm, all meet amplitude standards.Amplitude-versus-frequency curve as shown in Figure 6.
By calculating above, for the present embodiment, the 1 rank amplitude at rotor rear support place is relatively large, when the amount of unbalance of rotor reaches 0.44 × 10
-4during kgm, the 1 rank amplitude at rotor rear support place has reached standard limit 50 μm, therefore, provides the residual unbalance, standard limit of this rotor for being not more than 0.44 × 10
-4kgm.
So far, whole rotor dynamics design process terminates, and by the present embodiment, makes the structural dynamic characteristics of rotor obtain optimization, meets proposed all demanding kinetics.
Claims (1)
1. a design method for aeroengine rotor structural dynamics, is characterized in that, concrete steps are:
The first step: set up rotor dynamics model: obtain the vibratory response of rotor-support-foundation system at critical speed of rotation place and the relation of each design parameter, the relation of the vibratory response of described rotor-support-foundation system at critical speed of rotation place and each design parameter comprises the relation of rotor-support-foundation system threshold response and Structural Parameters of its Rotor and the relation of rotor-support-foundation system threshold response and rotor modal parameter;
The equation matrix expression (1) of rotor dynamics model is:
In formula (1):
q is the response of rotor-support-foundation system, and r is the amount of deflection at rotor c.g. place,
for the deflection angle at rotor c.g. place, r
b1for the stiffness coefficient of rotor-support-foundation system front support, r
b2for the stiffness coefficient of rotor-support-foundation system rear support;
Rotor-support-foundation system is at the threshold response of supporting place
The threshold response at rotor c.g. place
Formula (17), formula (18) describe the vibratory response Q of flexibly mounted rotor-support-foundation system at critical speed of rotation place
1cr, Q
2cr, include the modal parameter of rotor-support-foundation system structural dynamic parameter and rotor-support-foundation system; Formula (17), formula (18) describe flexibly mounted rotor-support-foundation system threshold response and rotor critical speed during supporting absolute rigidity
relation;
In formula (17) and formula (18),
Q
1crfor the threshold response amplitude at rotor c.g. place; Q
2crfor the threshold response amplitude at rotor bearing place;
Ω
crifor rotor i-th rank critical speed of rotation;
The modal matrix that [Φ] is rotor-support-foundation system,
be respectively the first rank and second-order damping ratios, d
cr1, d
cr2be respectively the first rank and the second-order modal damping coefficient of rotor-support-foundation system, m
cr1, I
cr2be respectively the first rank and second-order modal mass, Ω
cr1, Ω
cr2be respectively the first rank and the second-order critical speed of rotation of rotor;
when being respectively supporting rigidity, the critical speed of rotation of rotor-support-foundation system;
U
1for the amount of unbalance of rotor-support-foundation system;
Ω
crifor rotor i-th rank critical speed of rotation, Q
2crfor rotor is in the threshold response amplitude of supporting place;
Second step: the hot-die state during work of setting rotor-support-foundation system:
According to operating rotational speed range and the performance requirement of motor, the hot-die state of setting rotor-support-foundation system; When setting described hot-die state, must meet the following conditions successively: first step mode critical speed of rotation is positioned at below slow train rotating speed, and allow in slow train position; Second-order mode critical speed of rotation and the 3rd rank mode critical speed of rotation are all positioned at below cruising speed, nargin equal 10%; Critical speed of rotation is not established between cruising speed and maximum (top) speed; When the design parameter of rotor requires to establish critical speed of rotation between the cruising speed and maximum (top) speed of rotor-support-foundation system, then quadravalence mode critical speed of rotation is arranged on maximum (top) speed, nargin 8-10%;
3rd step: the mode of rotor-support-foundation system when determining to support absolute rigidity:
If each support stiffness of rotor-support-foundation system is infinitely great, transfer matrix method or finite element method is adopted to determine the mode of rotor-support-foundation system with this understanding;
4th step: the initial value choosing support stiffness: the initial value selecting support stiffness according to the rigidity of supporting absolute rigidity rotor-support-foundation system; Each support stiffness of rotor-support-foundation system should be the 20-50% of rigid support rotor rigidity, adopts finite element method or determines the lateral stiffness of rotor at barycenter place by formula (19):
E is the Young's modulus of material, and I is the cross section moments of inertia of rotating shaft;
Get 50% of s as the initial value of rotor-support-foundation system elastic supporting rigidity:
5th step: verify whether the quiet amount of deformation of rotor during the support stiffness initial value determined exceeds standard:
Adopt finite element method or the quiet distortion according to formula (21) calculating rotor-support-foundation system, make the quiet amount of deformation of rotor be not more than the value of this engine design requirement;
In formula, δ is quiet amount of deformation; M is whole rotor quality; G is gravity accleration; M is the quality of dish; S is rotor centroid position lateral stiffness, namely produces the power needed for unit lateral deformation;
6th step: check the mode of rotor during the support stiffness initial value for determining whether to conform to the hot-die state preset:
For selected support stiffness parameter, adopt the mode of transfer matrix method or Finite element arithmetic rotor-support-foundation system; Whether the mode of inspection rotor-support-foundation system conforms to the hot-die state preset in second step; If the mode of rotor-support-foundation system is not inconsistent with the hot-die state preset, then adjust the parameter s of support stiffness
b1and s
b2, the mode of rotor-support-foundation system is conformed to the hot-die state preset in second step;
After the mode of rotor-support-foundation system conforms to the hot-die state preset in second step, check described hot-die state and rotor-support-foundation system just to prop up whether to leave between mode the nargin of 10%; If nargin is inadequate, the parameter S of adjustment support stiffness
b1and S
b2or the quality of rotor, dimensional parameters, hot-die state and rotor-support-foundation system are just propped up to leave between mode the nargin of 10%;
Need the quiet deformation requirements of satisfied 5th step simultaneously;
7th step: damper is set in supporting place, determine unbalance response and requirement for dynamic balance:
Introduce squeeze film damper in supporting place, calculate the unbalance response of rotor; The damping function that can play according to damper and counterbalance effect, obtain rotor residual unbalance, limitation standard; So far, the Structural Dynamic Design of aeroengine rotor is completed.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410146849.XA CN103912315B (en) | 2014-04-14 | 2014-04-14 | Structural dynamics design method of rotor of aerial engine |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410146849.XA CN103912315B (en) | 2014-04-14 | 2014-04-14 | Structural dynamics design method of rotor of aerial engine |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103912315A CN103912315A (en) | 2014-07-09 |
CN103912315B true CN103912315B (en) | 2015-06-17 |
Family
ID=51038280
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410146849.XA Active CN103912315B (en) | 2014-04-14 | 2014-04-14 | Structural dynamics design method of rotor of aerial engine |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103912315B (en) |
Families Citing this family (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104536475B (en) * | 2014-11-29 | 2017-02-22 | 四川达宇特种车辆制造厂 | Method for adjusting the critical speed of high-speed combined rotor of micro turbojet engine |
CN105224713B (en) * | 2015-07-30 | 2018-08-21 | 西北工业大学 | A kind of derated design method of aero-engine intershaft bearing |
CN105488286B (en) * | 2015-12-04 | 2018-05-01 | 湖南科技大学 | Vertical multi-stage centrifugal pump rotor dynamics modeling method under a kind of band liquid operating mode |
CN107066725B (en) * | 2017-04-11 | 2020-07-03 | 西北工业大学 | Rotor structure dynamics reverse design method based on fault genes |
CN107025348B (en) * | 2017-04-11 | 2020-04-03 | 西北工业大学 | Early accumulated liquid fault identification method for aircraft engine rotor |
CN108287952B (en) * | 2018-01-04 | 2021-01-19 | 西安交通大学 | Low-temperature liquid expander rotor critical rotation speed prediction method considering radial fluid force influence |
CN109063356B (en) * | 2018-08-15 | 2019-10-15 | 东南大学 | A kind of high-speed electric main shaft rotor-bearing-enclosure system dynamic design approach |
CN110877751B (en) * | 2019-12-04 | 2021-04-30 | 哈尔滨工业大学 | Large-scale high-speed rotation equipment rotational inertia stacking method based on vector projection |
CN111219210B (en) * | 2020-01-08 | 2021-08-10 | 西北工业大学 | Dynamics optimization design method based on reverse rotation double-rotor whirling direction reversing |
CN111553033B (en) * | 2020-04-28 | 2022-07-05 | 哈尔滨工业大学 | Equivalent reduction processing method for heavy-duty gas turbine pull rod rotor model based on strain energy method |
CN113486561B (en) * | 2021-07-20 | 2022-06-28 | 中国科学院工程热物理研究所 | Engine rotor dynamic characteristic improving method based on strain energy distribution |
CN113486560B (en) * | 2021-07-20 | 2022-07-26 | 中国科学院工程热物理研究所 | Dynamic design method for engine center pull rod rotor-supporting system |
CN114491915B (en) * | 2021-10-20 | 2023-06-30 | 中国航发四川燃气涡轮研究院 | Bearing performance analysis method considering rotor and spring bearing influence |
CN114330083A (en) * | 2021-12-15 | 2022-04-12 | 重庆邮电大学 | Centrifugal pump rotor dynamics steady prediction method with uncertain bearing parameters |
CN115680779B (en) * | 2022-10-25 | 2024-04-30 | 西北工业大学 | Design method of multi-extrusion oil film damper of double-rotor system |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101329220A (en) * | 2008-07-11 | 2008-12-24 | 西安交通大学 | Dynamics characteristic experimental apparatus of heavy type gas turbine pull rod rotor |
EP2400119A2 (en) * | 2010-06-23 | 2011-12-28 | Honeywell International, Inc. | Gas turbine engine rotor tip clearance and shaft dynamics system and method |
CN103308313A (en) * | 2013-05-15 | 2013-09-18 | 天津大学 | Base-movable aeroengine dual-rotor system model experiment platform |
CN103471854A (en) * | 2013-09-26 | 2013-12-25 | 沈阳黎明航空发动机(集团)有限责任公司 | Aero-engine body vibration characteristic analysis method |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9046001B2 (en) * | 2011-08-29 | 2015-06-02 | Honeywell International Inc. | Annular bearing support dampers, gas turbine engines including the same, and methods for the manufacture thereof |
-
2014
- 2014-04-14 CN CN201410146849.XA patent/CN103912315B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101329220A (en) * | 2008-07-11 | 2008-12-24 | 西安交通大学 | Dynamics characteristic experimental apparatus of heavy type gas turbine pull rod rotor |
EP2400119A2 (en) * | 2010-06-23 | 2011-12-28 | Honeywell International, Inc. | Gas turbine engine rotor tip clearance and shaft dynamics system and method |
CN103308313A (en) * | 2013-05-15 | 2013-09-18 | 天津大学 | Base-movable aeroengine dual-rotor system model experiment platform |
CN103471854A (en) * | 2013-09-26 | 2013-12-25 | 沈阳黎明航空发动机(集团)有限责任公司 | Aero-engine body vibration characteristic analysis method |
Also Published As
Publication number | Publication date |
---|---|
CN103912315A (en) | 2014-07-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103912315B (en) | Structural dynamics design method of rotor of aerial engine | |
CN102803713B (en) | Wind turbine with compensated motor torque | |
Capuzzi et al. | A novel adaptive blade concept for large-scale wind turbines. Part I: Aeroelastic behaviour | |
CN108256210B (en) | Integral coupling analysis method for offshore wind turbine under earthquake action | |
CN110619145B (en) | Self-adaptive modeling method for flexible supporting gear transmission device | |
Ma et al. | Vibration response analysis caused by rubbing between rotating blade and casing | |
Yang et al. | Vibration reduction optimum design of a steam-turbine rotor-bearing system using a hybrid genetic algorithm | |
Liu et al. | Vibration reduction strategy for wind turbine based on individual pitch control and torque damping control | |
Soykasap et al. | Performance enhancement of a composite tilt-rotor using aeroelastic tailoring | |
CN103208795B (en) | Method for optimizing parameters of supplementary sub-synchronous damping controller for high-voltage direct-current transmission | |
Şener et al. | Evaluation of the effect of spar cap fiber angle of bending–torsion coupled blades on the aero-structural performance of wind turbines | |
Li et al. | Aerodynamic response analysis of wind turbines | |
Larwood | Dynamic analysis tool development for advanced geometry wind turbine blades | |
Zhang et al. | Numerical analysis on propulsive efficiency and pre-deformated optimization of a composite marine propeller | |
Kim et al. | Aerodynamic and structural design of 6kW class Vertical-axis wind turbine | |
Liao et al. | Finite element analysis and Lightweight design of hydro generator lower bracket | |
Jensen et al. | Numerical simulation of gyroscopic effects in Ansys | |
Henriksen et al. | HAWCStab2 user manual | |
Tomas-Rodriguez et al. | Floating Offshore Wind Turbines Oscillations Damping. | |
Gantasala et al. | Aeroelastic simulations of wind turbine using 13 DOF rigid beam model | |
Pan et al. | Effects of aerodynamic fairing on full scale blade fatigue test | |
CN107066653A (en) | A kind of aeroelastic analysis method for considering engine dynamics | |
Li et al. | Dynamic response analysis of blades with damping structures of shroud and snubber | |
Shi et al. | Influences of blade stiffness on structural natural frequencies of offshore wind turbine with monopile foundation | |
Miao et al. | Nonlinear response analysis of variable speed rotor system under maneuvering flight |
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 |