CN118095019A - Method and device for calculating vibration fatigue crack extension life of engine structure - Google Patents
Method and device for calculating vibration fatigue crack extension life of engine structure Download PDFInfo
- Publication number
- CN118095019A CN118095019A CN202410512958.2A CN202410512958A CN118095019A CN 118095019 A CN118095019 A CN 118095019A CN 202410512958 A CN202410512958 A CN 202410512958A CN 118095019 A CN118095019 A CN 118095019A
- Authority
- CN
- China
- Prior art keywords
- crack
- life
- stress
- model
- engine structure
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 59
- 238000001228 spectrum Methods 0.000 claims abstract description 70
- 230000003938 response to stress Effects 0.000 claims abstract description 62
- 238000004364 calculation method Methods 0.000 claims abstract description 43
- 238000012360 testing method Methods 0.000 claims abstract description 29
- 238000004458 analytical method Methods 0.000 claims abstract description 23
- 238000011056 performance test Methods 0.000 claims abstract description 17
- 238000012545 processing Methods 0.000 claims abstract description 14
- 238000005070 sampling Methods 0.000 claims description 13
- 230000003595 spectral effect Effects 0.000 claims description 10
- 238000005457 optimization Methods 0.000 claims description 8
- 125000004122 cyclic group Chemical group 0.000 claims description 7
- 230000005284 excitation Effects 0.000 claims description 7
- 238000012512 characterization method Methods 0.000 abstract description 2
- 230000004044 response Effects 0.000 description 10
- 239000000463 material Substances 0.000 description 7
- 238000013461 design Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 230000007246 mechanism Effects 0.000 description 4
- 238000012986 modification Methods 0.000 description 4
- 230000004048 modification Effects 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 3
- 229910000831 Steel Inorganic materials 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 230000001419 dependent effect Effects 0.000 description 2
- 230000000977 initiatory effect Effects 0.000 description 2
- 239000007788 liquid Substances 0.000 description 2
- 230000000149 penetrating effect Effects 0.000 description 2
- 230000035515 penetration Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 239000010959 steel Substances 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 230000001627 detrimental effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000013213 extrapolation Methods 0.000 description 1
- 238000010998 test method 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
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
The invention discloses a calculation method and a calculation device for the vibration fatigue crack growth life of an engine structure, relates to the technical field of mechanical property test characterization, aims to obtain reliable fracture mechanical parameters and crack growth rate curves, provides accurate quantized crack growth life, and provides method guidance for repeated use engine life assessment. The calculation method comprises the following steps: acquiring a time domain random stress spectrum of stress response of an assessment section in a finite element model of an engine structure; based on the time domain random stress spectrum, carrying out finite element analysis processing on a crack model of the engine structure, and determining fracture mechanical parameters of the crack-containing structure; performing vibration fatigue performance test on a crack growth rate model established according to a preset crack growth mode, and determining a target parameter value corresponding to the crack growth rate model by utilizing parameter identification; and calculating the crack extension life by combining the parameter value corresponding to the crack extension rate model and a preset critical crack length value to obtain the target crack extension life.
Description
Technical Field
The invention relates to the technical field of mechanical property test characterization, in particular to a method and a device for calculating the expansion life of a vibration fatigue crack of an engine structure.
Background
The reuse of the engine is one of the main development directions of the spaceflight carrying power system, and under the reuse requirement, the fatigue strength design and service life assessment requirements of the engine are very outstanding. The conventional fatigue analysis method performs life prediction based on a stress/strain-life curve (S-N curve), and can give a life prediction result with a certain degree of reliability. However, when the fatigue analysis method is used for reusing a structure, the following exists: the S-N curve life model has fewer parameters, and a fatigue damage process mechanism cannot be completely described; in addition, the load factor in fatigue analysis mainly considers the magnitude and occurrence frequency of cyclic stress, and it is difficult to consider the influence of load sequence on life. This can lead to a large error in life prediction results, which is detrimental to structural design and evaluation.
The damage tolerance method and the crack propagation technology based on fracture mechanics are widely applied in the aviation field, and at present, the method mainly adopts a linear elastic fracture mechanics method, and is suitable for macroscopic long crack modeling, material performance test and quasi-static fatigue life calculation. The vibration load of the liquid rocket engine is a remarkable load characteristic in the repeated use process, and the fatigue life of the liquid rocket engine is mostly in a crack initiation stage under the vibration environment, so that the structural integrity is required to be maintained for ensuring the normal operation of the engine, and the occurrence of macro-cracks or cracks penetrating through the wall thickness is generally not allowed. Therefore, structural response characteristics under vibration load, damage modes of unperforated wall thickness cracks and effective performance data are important to consider when designing engine fatigue and evaluating service life based on fracture mechanics. At present, the linear elastic fracture mechanics have the limitations of uncertain parameters and the like when the linear elastic fracture mechanics are used for unperforated cracks in the crack initiation stage, in addition, the elastoplasticity and other nonlinear conditions cannot be considered by directly introducing the cracks in the traditional vibration response calculation based on a modal method, and the transient dynamics calculation scale and calculation efficiency of the crack problem limit the general application of the linear elastic fracture mechanics to engineering. Therefore, the current fracture mechanics method is difficult to adapt to the requirement of fatigue crack life assessment of an engine in a vibration environment.
Based on the above, it is necessary to develop a fatigue crack growth life prediction method for an engine structure under vibration load, so as to obtain reliable fracture mechanical parameters and crack growth rate curves, give accurate quantized crack growth life, and provide method guidance for repeated use of engine life assessment.
Disclosure of Invention
The invention aims to provide a method and a device for calculating the vibration fatigue crack growth life of an engine structure, so as to obtain reliable fracture mechanical parameters and crack growth rate curves, provide accurate quantized crack growth life and provide method guidance for repeated use of engine life assessment.
In a first aspect, the present invention provides a method for calculating the vibration fatigue crack growth life of an engine structure, the method comprising:
acquiring a time domain random stress spectrum of stress response of an assessment section in a finite element model of an engine structure;
based on the time domain random stress spectrum, carrying out finite element analysis processing on a crack model of the engine structure, and determining fracture mechanical parameters of the crack-containing structure;
Performing vibration fatigue performance test on a crack growth rate model established according to a preset crack growth mode, and determining a target parameter value corresponding to the crack growth rate model by utilizing parameter identification;
And calculating the crack extension life by combining the parameter value corresponding to the crack extension rate model and a preset critical crack length value to obtain the target crack extension life.
Under the condition of adopting the technical scheme, after obtaining the time domain random stress spectrum of the stress response of the checking section in the finite element model of the engine structure, based on the time domain random stress spectrum, finite element analysis processing can be carried out on the crack model of the engine structure, so that fracture mechanical parameters of the crack-containing structure are determined, vibration fatigue performance tests are carried out on the crack expansion rate model established according to the preset crack expansion mode, the parameter identification is utilized to determine target parameter values corresponding to the crack expansion rate model, and finally, crack expansion life calculation can be carried out by combining the parameter values corresponding to the crack expansion rate model and the preset critical crack length values, so that the target crack expansion life is obtained. Based on the method, the method integrates the contents of structural random vibration response calculation, vibration stress response time domain simulation, crack modeling, fracture mechanical parameter calculation, crack expansion rate model parameter, initial crack length determination, crack expansion service life calculation and the like, can obtain structural random vibration steady-state response and time domain stress random spectrum, give accurate quantized fracture mechanical parameters, determine crack expansion rate model parameter and initial crack length based on fatigue performance test data, realize structural crack expansion service life calculation, and overcome the defects of large calculation result deviation caused by unclear mechanism, difficult parameter acquisition and unsuitable model when the traditional fatigue analysis method and the macroscopic crack analysis method of linear elastic fracture mechanics are used for calculating the service life of the non-penetrating thickness crack in a vibration environment.
Therefore, the calculation method of the vibration fatigue crack growth life of the engine structure can obtain reliable fracture mechanical parameters and crack growth rate curves, give accurate quantized crack growth life and provide method guidance for repeated use of engine life assessment.
In a second aspect, the present invention also provides a device for calculating the vibration fatigue crack growth life of an engine structure, for implementing the method for calculating the vibration fatigue crack growth life of an engine structure according to the first aspect, the computing device comprising:
The acquisition module is used for acquiring a time domain random stress spectrum of stress response of the examination section in the finite element model of the engine structure;
the first determining module is used for carrying out finite element analysis processing on a crack model of the engine structure based on the time domain random stress spectrum and determining fracture mechanical parameters of the crack-containing structure;
The second determining module is used for performing vibration fatigue performance test on the crack growth rate model established according to the preset crack growth mode, and determining a target parameter value corresponding to the crack growth rate model by utilizing parameter identification;
The obtaining module is used for carrying out crack extension service life calculation by combining the parameter value corresponding to the crack extension rate model and the preset critical crack length value to obtain the target crack extension service life.
Optionally, the acquiring module includes:
The first building unit is used for building a finite element model of the engine structure;
The first determining unit is used for applying random vibration load excitation to the finite element model of the engine structure and determining the power spectral density of stress response of the checking section;
The first obtaining unit is used for sampling the power spectrum density to obtain a time domain random stress spectrum of the stress response of the examination section.
Optionally, the first obtaining unit includes
The acquisition subunit is used for acquiring a stress amplitude probability density function corresponding to the power spectrum density;
The obtaining subunit is used for sampling the time stress response based on the stress amplitude probability density function and the peak probability density function to obtain a stress response sample;
And the determining subunit is used for carrying out randomization processing on the stress response sample and determining a time domain random stress spectrum of the stress response of the checking section.
Optionally, the first determining module includes:
The second building unit is used for building a crack model of the engine structure;
A second determining unit for determining an external load applied to the crack model based on the time domain random stress spectrum;
And the third determining unit is used for carrying out elastic-plastic finite element analysis processing on the crack model with the external load applied to determine fracture mechanical parameters corresponding to the lengths of the cracks.
Optionally, the fracture mechanics parameter comprises at least a stress intensity factor or a J-integral.
Optionally, the second determining module includes:
A third establishing unit, configured to establish a crack propagation rate model according to a preset crack propagation mode;
The acquisition unit is used for acquiring a vibration fatigue stress-life test curve or a vibration fatigue strain-life test curve of the engine structure;
A fourth determining unit for determining a test crack growth life based on the vibration fatigue stress-life test curve or the vibration fatigue strain-life test curve in combination with initial values of the initial crack length and crack growth rate parameters;
And the fifth determining unit is used for optimizing the initial value of the crack expansion rate parameter by using a parameter optimization method, and determining the parameter value with the smallest error as the parameter value corresponding to the crack expansion rate model.
Optionally, the crack growth rate model includes:
，
Wherein a is crack length, N is cycle number, C, N, p, q are undetermined model parameter values, f is crack opening/closing function, R is stress ratio, ΔK is stress intensity factor amplitude, ΔK th is crack propagation threshold, K max is stress intensity factor peak value, and K C is fracture toughness.
Optionally, the obtaining module includes:
The second obtaining unit is used for calculating the crack extension life by using a cyclic joint cycle method or a block spectrum average life calculation method in combination with the time domain random stress spectrum, the fracture mechanical parameter, the target parameter value corresponding to the crack extension rate model and the preset critical crack length value to obtain the target crack extension life.
The beneficial effects of the device for calculating the vibration fatigue crack growth life of the engine structure provided in the second aspect are the same as those of the method for calculating the vibration fatigue crack growth life of the engine structure described in the implementation manner of the first aspect, and are not described here in detail.
Drawings
The accompanying drawings, which are included to provide a further understanding of the invention and are incorporated in and constitute a part of this specification, illustrate embodiments of the invention and together with the description serve to explain the invention and do not constitute a limitation on the invention. In the drawings:
FIG. 1 is a flow chart of steps of a method for calculating the life of vibration fatigue crack growth of an engine structure according to an embodiment of the present invention;
FIG. 2 is a flow chart of steps of a method for calculating vibration fatigue crack growth life of an engine structure according to another embodiment of the present invention;
FIG. 3 is a graph of the power spectral density of the vibration stress response of a structural assessment site according to an embodiment of the present invention;
FIG. 4 (a) is a typical narrow band distribution curve of stress response;
FIG. 4 (b) is a typical broadband profile of stress response;
FIG. 5 is a schematic diagram of a time domain stress random spectrum in an embodiment of the present invention;
FIG. 6 is a schematic diagram of a crack model of an engine structure in an embodiment of the present invention;
FIG. 7 is a schematic illustration of a surface crack in an embodiment of the present invention;
FIG. 8 is a schematic diagram of crack growth rate model parameters and initial crack length determination in an embodiment of the present invention;
FIG. 9 is a schematic diagram of a device for calculating the vibration fatigue crack growth life of an engine structure according to an embodiment of the present invention.
Detailed Description
In order to clearly describe the technical solution of the embodiments of the present invention, in the embodiments of the present invention, the words "first", "second", etc. are used to distinguish the same item or similar items having substantially the same function and effect. For example, the first threshold and the second threshold are merely for distinguishing between different thresholds, and are not limited in order. It will be appreciated by those of skill in the art that the words "first," "second," and the like do not limit the amount and order of execution, and that the words "first," "second," and the like do not necessarily differ.
In the present invention, the words "exemplary" or "such as" are used to mean serving as an example, instance, or illustration. Any embodiment or design described herein as "exemplary" or "for example" should not be construed as preferred or advantageous over other embodiments or designs. Rather, the use of words such as "exemplary" or "such as" is intended to present related concepts in a concrete fashion.
In the present invention, "at least one" means one or more, and "a plurality" means two or more. "and/or", describes an association relationship of an association object, and indicates that there may be three relationships, for example, a and/or B, and may indicate: a alone, a and B together, and B alone, wherein a, B may be singular or plural. The character "/" generally indicates that the context-dependent object is an "or" relationship. "at least one of" or the like means any combination of these items, including any combination of single item(s) or plural items(s). For example, at least one (one) of a, b or c may represent: a, b, c, a and b, a and c, b and c, or a, b and c, wherein a, b, c can be single or multiple.
As shown in fig. 1, an embodiment of the present invention provides a method for calculating an engine structural vibration fatigue crack growth life, including:
s101: and obtaining a time domain random stress spectrum of stress response of the checking section in the finite element model of the engine structure.
In the application, a finite element method is adopted to carry out dynamic modeling on a structure, and random vibration load excitation is applied to obtain dynamic responses such as power spectral density (Power Spectrum Density, PSD) and Root Mean Square (RMS) of stress response of an examination part.
And analyzing and giving a probability density function of a time domain stress response amplitude and a peak value according to the obtained power spectrum density of the vibration stress response, and then sampling to obtain a stress response time domain random number and form a time domain random stress spectrum.
S102: based on the time domain random stress spectrum, finite element analysis processing is carried out on the crack model of the engine structure, and fracture mechanical parameters of the crack-containing structure are determined.
Specifically, a crack propagation mode is first assumed, a crack form is given, and a crack propagation path is set. Then, crack modeling is carried out on the engine structure, and two points are mainly considered in the structural crack model: firstly, the expansion range from the initial crack to the critical crack length can be enveloped; secondly, the steel comprises main structural characteristics such as fillets, holes, variable thickness and the like which have influence on fracture mechanical parameters of cracks with no penetration thickness. And finally, selecting load points in the time domain random stress spectrum as external load to be applied to the crack model for linear elasticity or elastic-plastic finite element calculation, and obtaining fracture mechanical parameters Ki (stress intensity factors) or Ji (J integral) corresponding to each crack length ai.
S103: and carrying out vibration fatigue performance test on the crack growth rate model established according to the preset crack growth mode, and determining a target parameter value corresponding to the crack growth rate model by utilizing parameter identification.
Specifically, firstly, vibration fatigue performance test data such as a fatigue stress-life curve or a strain-life curve of a material or a structural member under vibration load are obtained, a crack expansion rate model is given, and an initial crack form is assumed; then, fracture mechanics calculation is carried out on the vibration fatigue performance test piece, and fracture mechanics parameters such as stress intensity factors Ktest are obtained; finally, presetting a crack growth rate model parameter initial value and an initial crack length initial value for each stress level of vibration fatigue performance test data (curve), calculating the test crack growth life of each stress level, and obtaining a parameter optimization result with the smallest error with the test life result as a final parameter identification value by a parameter optimization method.
Based on the method, a material fatigue damage model represented by a single mechanical parameter S (stress or strain) and a service life N is converted into a material fatigue damage model represented by an initial crack size and a fracture mechanical property parameter, and a crack propagation rate model parameter value is obtained.
S104: and calculating the crack extension life by combining the parameter value corresponding to the crack extension rate model and a preset critical crack length value to obtain the target crack extension life.
Specifically, substituting a time domain random stress spectrum and fracture mechanics parameters into a crack expansion rate model, and calculating crack expansion life by specifying critical crack length af and initial crack length a0 and adopting a cyclic joint-cyclic method or a block spectrum average life calculation method to obtain final crack expansion life. It should be appreciated that a set of stress random spectra corresponds to a set of crack growth life values.
Compared with the prior art, the method for calculating the vibration fatigue crack extension life of the engine structure provided by the embodiment of the invention can be used for carrying out finite element analysis processing on the crack model of the engine structure based on the time domain random stress spectrum after obtaining the time domain random stress spectrum of the stress response of the checking section in the finite element model of the engine structure, so as to determine the fracture mechanical parameters of the crack-containing structure, carrying out vibration fatigue performance test on the crack extension rate model established according to the preset crack extension mode, determining the target parameter value corresponding to the crack extension rate model by utilizing parameter identification, and finally carrying out crack extension life calculation by combining the parameter value corresponding to the crack extension rate model and the preset critical crack length value to obtain the target crack extension life. Based on the above, the embodiment of the invention integrates the contents of structural random vibration response calculation, vibration stress response time domain simulation, crack modeling, fracture mechanical parameter calculation, crack expansion rate model parameter, initial crack length determination, crack expansion service life calculation and the like, can obtain structural random vibration steady-state response and time domain stress random spectrum, give accurate quantitative fracture mechanical parameter, determine crack expansion rate model parameter and initial crack length based on fatigue performance test data, realize structural crack expansion service life calculation, and overcome the defects of large calculation result deviation caused by unclear mechanism, difficult parameter acquisition and inapplicability of the model when the traditional fatigue analysis method and the macroscopic crack analysis method of linear elastic fracture mechanics are used for service life calculation of the non-penetrated thickness crack in a vibration environment.
Therefore, the calculation method of the vibration fatigue crack growth life of the engine structure provided by the embodiment of the invention can obtain reliable fracture mechanical parameters and crack growth rate curves, give accurate quantized crack growth life and provide method guidance for repeated use engine life assessment.
As shown in fig. 2, the embodiment of the invention further provides another method for calculating the vibration fatigue crack growth life of the engine structure, and the specific steps of the method for calculating the vibration fatigue crack growth life of the engine structure provided by the embodiment of the invention will be described in detail with reference to fig. 2 to 8.
S201: and establishing a finite element model of the engine structure.
S202: and (3) applying random vibration load excitation to the finite element model of the engine structure, and determining the power spectral density of stress response of the examination section.
In the application, a finite element method is adopted to carry out dynamic simulation modeling on a structure, and random vibration load excitation is applied to obtain power spectral density (Power Spectrum Density, PSD) and Root Mean Square (RMS) of stress response of an examined part, as shown in figure 3. FIG. 3 schematically shows a stress power spectral density curve of a severe vibration stress response part of an engine structure, namely an examination section, wherein the abscissa represents frequency in Hz, and the ordinate represents stress power spectral density in MPa 2/Hz.
S203: and sampling the power spectrum density to obtain a time domain random stress spectrum of stress response of the examination section.
Specifically, the probability density function giving the time domain stress response is analyzed for the power spectrum density of the obtained vibration stress response, then the time domain random number of the stress response is obtained by sampling, and the time domain random stress spectrum is formed.
The step S203 includes the following sub-steps:
substep A1: obtaining a stress amplitude probability density function corresponding to the power spectrum density;
substep A2: sampling the time stress response based on the stress amplitude probability density function and the peak probability density function to obtain a stress response sample;
substep A3: and carrying out randomization treatment on the stress response sample, and determining a time domain random stress spectrum of stress response of the checking section.
Analyzing the PSD spectrum of the vibration stress response, FIG. 4 (a) is a typical narrow band distribution curve of the stress response; fig. 4 (b) is a typical broadband profile of stress response. As shown in fig. 4 (b), if the broadband is randomly distributed, a Dirlik model is adopted to obtain a stress amplitude probability density function of the broadband distribution; as shown in fig. 4 (a), if the Narrow-band distribution is adopted, a Narrow-band model is adopted to obtain a stress amplitude probability density function of the Narrow-band distribution; and simultaneously, combining with a hypothesized peak probability density function obeying Gaussian distribution, sampling the time stress response to obtain a group of stress response samples with a sufficient number, and finally obtaining a time domain random stress spectrum by a randomization method, wherein in fig. 5, the abscissa represents time, the ordinate represents stress and the ordinate represents Mpa, and the abscissa represents time.
S204: based on the time domain random stress spectrum, finite element analysis processing is carried out on the crack model of the engine structure, and fracture mechanical parameters of the crack-containing structure are determined.
Wherein the fracture mechanics parameter comprises at least a stress intensity factor or a J integral.
The step S204 includes the following sub-steps:
substep B1: establishing a crack model of the engine structure, as shown in fig. 6;
substep B2: determining an external load applied to the crack model based on the time domain random stress spectrum;
Substep B3: and carrying out elastic-plastic finite element analysis treatment on the crack model with the external load, and determining fracture mechanical parameters corresponding to the lengths of the cracks.
The specific test method comprises the following steps: first, assuming a Crack propagation mode, a Crack form, such as edge angle Crack (EDGE CRACK), surface Crack (Surface Crack), is given as shown in fig. 7. Or deep Crack (Embedded Crack), etc., to set a Crack propagation path. Then, finite element modeling is carried out on a crack model of the engine structure, and two points are mainly considered in the structural crack model: firstly, the expansion range from the initial crack to the critical crack length can be enveloped; secondly, the steel comprises main structural characteristics such as fillets, holes, variable thickness and the like which have influence on fracture mechanical parameters of cracks with no penetration thickness. And finally, applying the time domain random stress spectrum of the stress response of the examination section as an external load to a crack model for carrying out elastic-plastic finite element calculation to obtain fracture mechanical parameters Ki (stress intensity factors) or Ji (J integral) corresponding to each crack length ai.
Based on the method, the three-dimensional crack modeling of the actual complex structure is simplified into the problem of proper-scale crack modeling reflecting the main characteristics and stress response characteristics of the structural crack, meanwhile, the structural dynamic response is used as load input to carry out fracture mechanics calculation, the dynamic problem is converted into static problem, and then the elastoplastic fracture mechanics parameter solution can be obtained, so that the problem that elastoplastic analysis and nonlinear factor simulation cannot be carried out in steady-state vibration response calculation based on a modal method is solved.
S205: and establishing a crack growth rate model according to a preset crack growth mode.
Specifically, given a crack growth rate model, the model considers factors such as stress intensity factor amplitude, stress intensity factor maximum, stress intensity factor threshold, crack closure and the like, for example, a NASGRO crack growth rate formula is adopted:
，
Where a is the crack length, N is the number of cycles, C, N, p, q are the values of the model parameters to be determined, f is the crack opening/closing function, R is the stress ratio, ΔK is the stress intensity factor amplitude, ΔK th is the crack propagation threshold, K max is the stress intensity factor peak, and K C is the fracture toughness.
S206: and obtaining a vibration fatigue stress-life test curve or a vibration fatigue strain-life test curve of the engine structure.
S207: based on the vibration fatigue stress-life test curve or the vibration fatigue strain-life test curve, the initial crack length and the initial value of the crack growth rate parameter are combined to determine the test crack growth life.
S208: and optimizing the initial value of the crack expansion rate parameter by using a parameter optimization method, and determining the parameter value with the minimum error as the parameter value corresponding to the crack expansion rate model.
In the present application, vibration fatigue stress-life or strain-life test data (curve) of a material or a structural member is obtained, an initial value (a 0 ini,Cini,nini,pini,qini) of an initial crack length a0 and a crack growth rate parameter is assumed for each stress/strain level of a test, and a test crack growth life is calculated, and a parameter optimization result having the smallest error with a test life result is obtained as a final target parameter value by a parameter optimization method, as shown in fig. 8.
Based on the method, the material fatigue damage model represented by the single mechanical parameter S (stress or strain) and the service life N is converted into the material fatigue damage model represented by the crack size and the fracture mechanical property parameter thereof, and the crack expansion rate model parameter value is obtained.
S209: and calculating the crack extension life by combining the parameter value corresponding to the crack extension rate model and a preset critical crack length value to obtain the target crack extension life.
In the application, a target crack propagation life is calculated by combining a time domain random stress spectrum, a fracture mechanical parameter, a target parameter value corresponding to a crack propagation rate model and a preset critical crack length value and by using a cyclic joint cyclic method or a block spectrum average life calculation method, so as to obtain the target crack propagation life.
It should be appreciated that a set of stress random spectra corresponds to a set of crack growth life values.
The embodiment of the invention has the beneficial effects that:
(1) Through steady-state vibration response calculation and stress response time domain simulation, not only the frequency domain characteristics of structural stress response are considered, but also the time domain statistical characteristics of stress response are considered, the stress characteristics under vibration load are completely reserved, and the service life calculation is more complete and fine;
(2) The time domain stress response is used as an external load to carry out crack modeling, so that a load basis is provided for converting the integral crack modeling of the structure into the crack modeling of the local structure, the model scale is reduced, and the calculation efficiency is greatly improved. On the other hand, the dynamic calculation is converted into the static calculation, so that more real elastoplasticity and nonlinear solution can be obtained, and the accuracy of fracture mechanical parameters is obviously improved;
(3) The fatigue performance test data are fully utilized to identify the crack propagation rate model parameters and the initial crack length distribution, so that the source reliability of key performance parameters for calculating the crack propagation life is ensured, the modeling mechanism is clearer, the uncertainty of the traditional extrapolation method through a macroscopic long crack model is avoided, and the reasonable and reliable calculation model and calculation result are ensured.
As shown in fig. 9, an embodiment of the present invention further provides a device 300 for calculating the propagation life of a vibration fatigue crack of an engine structure, which is configured to implement the method for calculating the propagation life of a vibration fatigue crack of an engine structure in the above embodiment, where the computing device includes:
The acquisition module 301 is configured to acquire a time domain random stress spectrum of an assessment segment stress response in a finite element model of an engine structure;
the first determining module 302 is configured to perform finite element analysis on a crack model of the engine structure based on the time domain random stress spectrum, and determine fracture mechanical parameters of the crack-containing structure;
a second determining module 303, configured to perform a vibration fatigue performance test on a crack growth rate model established according to a preset crack growth mode, and determine a target parameter value corresponding to the crack growth rate model by using parameter identification;
and the obtaining module 304 is configured to calculate a crack growth lifetime by combining a parameter value corresponding to the crack growth rate model and a preset critical crack length value, so as to obtain a target crack growth lifetime.
Optionally, the acquiring module 301 includes:
The first building unit is used for building a finite element model of the engine structure;
The first determining unit is used for applying random vibration load excitation to the finite element model of the engine structure and determining the power spectral density of stress response of the checking section;
The first obtaining unit is used for sampling the power spectrum density to obtain a time domain random stress spectrum of the stress response of the examination section.
Optionally, the first obtaining unit includes
The acquisition subunit is used for acquiring a stress amplitude probability density function corresponding to the power spectrum density;
The obtaining subunit is used for sampling the time stress response based on the stress amplitude probability density function and the peak probability density function to obtain a stress response sample;
And the determining subunit is used for carrying out randomization processing on the stress response sample and determining a time domain random stress spectrum of the stress response of the checking section.
Optionally, the first determining module 302 includes:
The second building unit is used for building a crack model of the engine structure;
A second determining unit for determining an external load applied to the crack model based on the time domain random stress spectrum;
And the third determining unit is used for carrying out elastic-plastic finite element analysis processing on the crack model with the external load applied to determine fracture mechanical parameters corresponding to the lengths of the cracks.
Optionally, the fracture mechanics parameter comprises at least a stress intensity factor or a J-integral.
Optionally, the second determining module 303 includes:
A third establishing unit, configured to establish a crack propagation rate model according to a preset crack propagation mode;
The acquisition unit is used for acquiring a vibration fatigue stress-life test curve or a vibration fatigue strain-life test curve of the engine structure;
A fourth determining unit for determining a test crack growth life based on the vibration fatigue stress-life test curve or the vibration fatigue strain-life test curve in combination with initial values of the initial crack length and crack growth rate parameters;
And the fifth determining unit is used for optimizing the initial value of the crack expansion rate parameter by using a parameter optimization method, and determining the parameter value with the smallest error as the parameter value corresponding to the crack expansion rate model.
Optionally, the crack growth rate model includes:
，
Wherein a is crack length, N is cycle number, C, N, p, q are undetermined model parameter values, f is crack opening/closing function, R is stress ratio, ΔK is stress intensity factor amplitude, ΔK th is crack propagation threshold, K max is stress intensity factor peak value, and K C is fracture toughness.
Optionally, the obtaining module 304 includes:
The second obtaining unit is used for calculating the crack extension life by using a cyclic joint cycle method or a block spectrum average life calculation method in combination with the time domain random stress spectrum, the fracture mechanical parameter, the target parameter value corresponding to the crack extension rate model and the preset critical crack length value to obtain the target crack extension life.
The beneficial effects of the device for calculating the vibration fatigue crack extension life of the engine structure provided by the embodiment of the invention are the same as those of the method for calculating the vibration fatigue crack extension life of the engine structure described in the above embodiment, and are not repeated here.
Although the invention is described herein in connection with various embodiments, other variations to the disclosed embodiments can be understood and effected by those skilled in the art in practicing the claimed invention, from a study of the drawings, the disclosure, and the appended claims. In the claims, the word "comprising" does not exclude other elements or steps, and the "a" or "an" does not exclude a plurality. A single processor or other unit may fulfill the functions of several items recited in the claims. The mere fact that certain measures are recited in mutually different dependent claims does not indicate that a combination of these measures cannot be used to advantage.
Although the invention has been described in connection with specific features and embodiments thereof, it will be apparent that various modifications and combinations can be made without departing from the spirit and scope of the invention. Accordingly, the specification and drawings are merely exemplary illustrations of the present invention as defined in the appended claims and are considered to cover any and all modifications, variations, combinations, or equivalents that fall within the scope of the invention. It will be apparent to those skilled in the art that various modifications and variations can be made to the present invention without departing from the spirit or scope of the invention. Thus, it is intended that the present invention also include such modifications and alterations insofar as they come within the scope of the appended claims or the equivalents thereof.
Claims (10)
1. A method for calculating the propagation life of a vibration fatigue crack of an engine structure, the method comprising:
acquiring a time domain random stress spectrum of stress response of an assessment section in a finite element model of an engine structure;
based on the time domain random stress spectrum, performing finite element analysis processing on a crack model of the engine structure, and determining fracture mechanical parameters of the crack-containing structure;
performing vibration fatigue performance test on a crack growth rate model established according to a preset crack growth mode, and determining a target parameter value corresponding to the crack growth rate model by utilizing parameter identification;
And calculating the crack extension life by combining the parameter value corresponding to the crack extension rate model and a preset critical crack length value to obtain the target crack extension life.
2. The method for calculating the vibration fatigue crack growth life of the engine structure according to claim 1, wherein the obtaining the time-domain stochastic stress spectrum of the stress response of the examination section in the finite element model of the engine structure comprises:
Establishing a finite element model of the engine structure;
applying random vibration load excitation to the finite element model of the engine structure, and determining the power spectral density of stress response of the assessment section;
and sampling the power spectrum density to obtain a time domain random stress spectrum of the stress response of the examination section.
3. The method for calculating the vibration fatigue crack growth life of the engine structure according to claim 2, wherein the sampling the power spectrum density to obtain the time domain random stress spectrum of the stress response of the test section comprises:
acquiring a stress amplitude probability density function corresponding to the power spectrum density;
Sampling the time stress response based on the stress amplitude probability density function and the peak probability density function to obtain a stress response sample;
and carrying out randomization treatment on the stress response sample, and determining a time domain random stress spectrum of the stress response of the examination section.
4. The method for calculating the vibration fatigue crack growth life of an engine structure according to claim 1, wherein the performing finite element analysis processing on the crack model of the engine structure based on the time domain random stress spectrum to determine fracture mechanics parameters of the crack-containing structure comprises:
Establishing a crack model of the engine structure;
Determining an external load applied to the crack model based on the time domain random stress spectrum;
And carrying out elastic-plastic finite element analysis treatment on the crack model applied with the external load, and determining fracture mechanical parameters corresponding to the lengths of the cracks.
5. The method of claim 4, wherein the fracture mechanics parameter comprises at least a stress intensity factor or a J-integral.
6. The method for calculating the vibration fatigue crack growth life of an engine structure according to claim 1, wherein the performing the vibration fatigue performance test on the crack growth rate model established according to the preset crack growth mode, determining the target parameter value corresponding to the crack growth rate model by using parameter identification, comprises:
Establishing the crack growth rate model according to the preset crack growth mode;
acquiring a vibration fatigue stress-life test curve or a vibration fatigue strain-life test curve of the engine structure;
determining a test crack growth life based on the vibration fatigue stress-life test curve or the vibration fatigue strain-life test curve in combination with initial values of initial crack length and crack growth rate parameters;
And optimizing the initial value of the crack expansion rate parameter by using a parameter optimization method, and determining the parameter value with the minimum error as the parameter value corresponding to the crack expansion rate model.
7. The method of calculating a vibratory fatigue crack growth life for an engine structure of claim 6, wherein the crack growth rate model comprises:
，
Where a is the crack length, N is the number of cycles, C, N, p, q are the values of the model parameters to be determined, f is the crack opening/closing function, R is the stress ratio, ΔK is the stress intensity factor amplitude, ΔK th is the crack propagation threshold, K max is the stress intensity factor peak, and K C is the fracture toughness.
8. The method for calculating the crack growth life of the vibration fatigue of the engine structure according to claim 1, wherein the calculating the crack growth life by combining the parameter value corresponding to the crack growth rate model and the preset critical crack length value to obtain the target crack growth life comprises:
And carrying out crack extension life calculation by using a cyclic joint circulation method or a block spectrum average life calculation method in combination with the time domain random stress spectrum, the fracture mechanical parameter, a target parameter value corresponding to the crack extension rate model and a preset critical crack length value to obtain a target crack extension life.
9. A computing device for engine structural vibration fatigue crack growth life, the computing device comprising:
The acquisition module is used for acquiring a time domain random stress spectrum of stress response of the examination section in the finite element model of the engine structure;
the first determining module is used for carrying out finite element analysis processing on the crack model of the engine structure based on the time domain random stress spectrum to determine fracture mechanical parameters of the crack-containing structure;
the second determining module is used for performing vibration fatigue performance test on a crack expansion rate model established according to a preset crack expansion mode, and determining a target parameter value corresponding to the crack expansion rate model by utilizing parameter identification;
The obtaining module is used for carrying out crack extension service life calculation by combining the parameter value corresponding to the crack extension rate model and a preset critical crack length value to obtain a target crack extension service life.
10. The engine structural vibration fatigue crack growth life calculation device of claim 9, wherein the acquisition module comprises:
A first building unit for building a finite element model of the engine structure;
the first determining unit is used for applying random vibration load excitation to the finite element model of the engine structure and determining the power spectral density of stress response of the checking section;
The first obtaining unit is used for sampling the power spectrum density to obtain a time domain random stress spectrum of the stress response of the examination section.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410512958.2A CN118095019B (en) | 2024-04-26 | 2024-04-26 | Method and device for calculating vibration fatigue crack extension life of engine structure |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410512958.2A CN118095019B (en) | 2024-04-26 | 2024-04-26 | Method and device for calculating vibration fatigue crack extension life of engine structure |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN118095019A true CN118095019A (en) | 2024-05-28 |
| CN118095019B CN118095019B (en) | 2024-09-17 |
Family
ID=91155275
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202410512958.2A Active CN118095019B (en) | 2024-04-26 | 2024-04-26 | Method and device for calculating vibration fatigue crack extension life of engine structure |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN118095019B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN118332883A (en) * | 2024-06-17 | 2024-07-12 | 西安航天动力研究所 | Simulation calculation method and device for detail stress of metal hose |
| CN118332883B (en) * | 2024-06-17 | 2024-10-29 | 西安航天动力研究所 | Simulation calculation method and device for detail stress of metal hose |
Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0452573A1 (en) * | 1988-10-14 | 1991-10-23 | Prescott, Rhona Margaret Helen | Fatigue monitoring |
| CN102129512A (en) * | 2011-02-24 | 2011-07-20 | 西北工业大学 | Fatigue life analyzing method based on Paris formula |
| CN102567567A (en) * | 2011-11-15 | 2012-07-11 | 北京宇航系统工程研究所 | Finite element analysis based pipeline random-vibration fatigue life analyzing method |
| CN103616179A (en) * | 2013-12-05 | 2014-03-05 | 广西大学 | Transmission gear fatigue life assessment method based on defect modeling |
| US20140107948A1 (en) * | 2012-10-16 | 2014-04-17 | Christian Amann | Method and system for probabilistic fatigue crack life estimation |
| CN105760577A (en) * | 2016-01-28 | 2016-07-13 | 北京航空航天大学 | Estimation method for sound vibration fatigue life containing uncertain metal structure |
| CN111735695A (en) * | 2020-06-09 | 2020-10-02 | 西北工业大学 | Fatigue life prediction method based on CT test piece |
| CN111737901A (en) * | 2020-06-23 | 2020-10-02 | 石家庄铁道大学 | Cutter fatigue life prediction method and application thereof |
| CN112347668A (en) * | 2020-09-29 | 2021-02-09 | 华东交通大学 | Steel bridge deck fatigue reliability assessment method based on probabilistic fracture mechanics |
| US20210132599A1 (en) * | 2019-11-05 | 2021-05-06 | Sankar Narayanan | System, apparatus, and method for estimating life of components |
| WO2022121203A1 (en) * | 2020-12-08 | 2022-06-16 | 江苏科技大学 | Method for calculating spherical shell surface three-dimensional crack propagation fatigue life |
| CN117556652A (en) * | 2023-10-08 | 2024-02-13 | 上海交通大学 | Titanium alloy fatigue crack two-stage life prediction method based on finite element analysis |
-
2024
- 2024-04-26 CN CN202410512958.2A patent/CN118095019B/en active Active
Patent Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0452573A1 (en) * | 1988-10-14 | 1991-10-23 | Prescott, Rhona Margaret Helen | Fatigue monitoring |
| CN102129512A (en) * | 2011-02-24 | 2011-07-20 | 西北工业大学 | Fatigue life analyzing method based on Paris formula |
| CN102567567A (en) * | 2011-11-15 | 2012-07-11 | 北京宇航系统工程研究所 | Finite element analysis based pipeline random-vibration fatigue life analyzing method |
| US20140107948A1 (en) * | 2012-10-16 | 2014-04-17 | Christian Amann | Method and system for probabilistic fatigue crack life estimation |
| CN103616179A (en) * | 2013-12-05 | 2014-03-05 | 广西大学 | Transmission gear fatigue life assessment method based on defect modeling |
| CN105760577A (en) * | 2016-01-28 | 2016-07-13 | 北京航空航天大学 | Estimation method for sound vibration fatigue life containing uncertain metal structure |
| US20210132599A1 (en) * | 2019-11-05 | 2021-05-06 | Sankar Narayanan | System, apparatus, and method for estimating life of components |
| CN111735695A (en) * | 2020-06-09 | 2020-10-02 | 西北工业大学 | Fatigue life prediction method based on CT test piece |
| CN111737901A (en) * | 2020-06-23 | 2020-10-02 | 石家庄铁道大学 | Cutter fatigue life prediction method and application thereof |
| CN112347668A (en) * | 2020-09-29 | 2021-02-09 | 华东交通大学 | Steel bridge deck fatigue reliability assessment method based on probabilistic fracture mechanics |
| WO2022121203A1 (en) * | 2020-12-08 | 2022-06-16 | 江苏科技大学 | Method for calculating spherical shell surface three-dimensional crack propagation fatigue life |
| CN117556652A (en) * | 2023-10-08 | 2024-02-13 | 上海交通大学 | Titanium alloy fatigue crack two-stage life prediction method based on finite element analysis |
Non-Patent Citations (1)
| Title |
|---|
| 李小锋;李奕方;李建新;: "基于断裂力学的钢拱桥剩余寿命评估", 中国市政工程, no. 05, 25 October 2008 (2008-10-25) * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN118332883A (en) * | 2024-06-17 | 2024-07-12 | 西安航天动力研究所 | Simulation calculation method and device for detail stress of metal hose |
| CN118332883B (en) * | 2024-06-17 | 2024-10-29 | 西安航天动力研究所 | Simulation calculation method and device for detail stress of metal hose |
Also Published As
| Publication number | Publication date |
|---|---|
| CN118095019B (en) | 2024-09-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Pugno et al. | Evaluation of the non-linear dynamic response to harmonic excitation of a beam with several breathing cracks | |
| Quigley et al. | Review and assessment of frequency-based fatigue damage models | |
| Andreaus et al. | Fatigue crack growth, free vibrations, and breathing crack detection of aluminium alloy and steel beams | |
| Wang et al. | An Experimental‐Numerical Combined Method to Determine the True Constitutive Relation of Tensile Specimens after Necking | |
| Jiang et al. | Experimental design and validation of an accelerated random vibration fatigue testing methodology | |
| Karolczuk et al. | Modelling of stress gradient effect on fatigue life using Weibull based distribution function | |
| Rastogi et al. | Fatigue crack growth prediction in nuclear piping using Markov chain Monte Carlo simulation | |
| Zhang | A fuzzy residual strength based fatigue life prediction method | |
| Caiza et al. | Consideration of the runouts and their subsequent retests into SN curves modelling based on a three-parameter Weibull distribution | |
| WO1996041158A1 (en) | Method and apparatus for predicting structural integrity by estimating modal damping factor | |
| Marques et al. | Damage detection and fatigue life estimation under random loads: A new structural health monitoring methodology in the frequency domain | |
| CN114295491A (en) | Prediction method for creep damage and time evolution behavior of deformation | |
| CN118095019B (en) | Method and device for calculating vibration fatigue crack extension life of engine structure | |
| Turker et al. | Vibration based damage identification of concrete arch dams by finite element model updating | |
| Tao et al. | Notch fatigue life prediction considering nonproportionality of local loading path under multiaxial cyclic loading | |
| Habtour et al. | Modeling approach for predicting the rate of frequency change of notched beam exposed to Gaussian random excitation | |
| Mordini et al. | Damage detection on stay cables using an open source-based framework for finite element model updating | |
| Mahmud et al. | On the need to adopt strain-based probabilistic approach in predicting fatigue life | |
| Strauss et al. | Monitoring information and probabilistic-based prediction models for the performance assessment of concrete structures | |
| Rodríguez‐Reyes et al. | Comparative evaluation of fatigue life estimation under variable amplitude loading through damage models based on entropy | |
| CN114692331A (en) | Quantitative analysis method for response characteristics of aging state of rubber O-shaped ring sealing structure | |
| Gaier et al. | A Multiaxial Fatigue Analysis Method for Ductile, Semi-Ductile and Brittle Materials | |
| Yu et al. | Alternative determination of cable forces using flexural theory of axially loaded member | |
| Abdelkrim et al. | Fatigue life estimation of components with use a non-linear energy model coupled a finite element method | |
| Xiong et al. | A practical randomization approach of deterministic equation to determine probabilistic fatigue and fracture behaviours based on small experimental data sets |
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 |