Microsoft Word - numero_49_art_38_2453 M. Bannikov et alii, Frattura ed Integrità Strutturale, 49 (2019) 383-395; DOI: 10.3221/IGF-ESIS.49.38 383 Focused on Russian mechanics contributions for Structural Integrity Damage evolution in the AlMg6 alloy during high and very high cycle fatigue Mikhail Bannikov, Dmitry Bilalov, Vladimir Oborin, Oleg Naimark Institute of Continuous Media Mechanics UB RAS, Russia mbannikov@icmm.ru, http://orcid.org/0000-0002-5737-1422 ledon@icmm.ru, http://orcid.org/ 0000-0003-1541-2246 oborin@icmm.ru, https://orcid.org/0000-0003-2836-2073 naimark@icmm.ru, https://orcid.org/0000-0001-6537-1177 ABSTRACT. Paper presents the “in situ” method for determining of irreversible fatigue damage accumulation, based on the analysis of nonlinear manifestations of the feedback signal in a closed system of an ultrasonic fatigue system. During very high cycle (gigacycle) fatigue, the anomalies of the elastic properties of the material are appear, which leads to a nonlinearity effect in the amplitude of oscillations. This effect increases with the initiation and growth of fatigue cracks. The technology was applied to samples of AlMG-6 alloy with preliminary dynamic deformation to determine the moment of initiation and growth of the fatigue crack in very high cycle fatigue regime. This method is applicable for the early detection of fatigue cracks both on the surface and inside the material under cyclic loading in the ultrasonic mode. On the basis of wide-range defining relations for a deformable solid body with mesoscopic defects, a mathematical model has been proposed that can adequately describe behavior of the material during fatigue failure. The results of mathematical modeling are in good agreement with the experimental data KEYWORDS. Damage accumulation; Scaling; High- and very high cycle fatigue; Mathematical modeling; Destruction, mesoscopic defects. Citation: Bannikov, M., Bilalov, D., Oborin, V., Naimark, O., Damage evolution in the AlMg6 alloy during high and very high cycle fatigue, Frattura ed Integrità Strutturale, 49 (2019) 383-395. Received: 01.04.2019 Accepted: 14.06.2019 Published: 01.07.2019 Copyright: © 2019 This is an open access article under the terms of the CC-BY 4.0, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. INTRODUCTION he determination of fatigue cracks origin nature is one of the most important fundamental problem for various application areas, especially speaking about ultrahigh-cycle (gigacycle) [1-3] fatigue, when crack forms inside the material. The series of catastrophes caused by the fatigue destruction of gas turbine engines [4-6], combined with the high cost of resource assessment and the potential cost of developing new structures, stimulated promising concepts of national programs in the field of high and very high cycle fatigue, based on the use of new fundamental results in T http://www.gruppofrattura.it/VA/49/2453.mp4 M. Bannikov et alii, Frattura ed Integrità Strutturale, 49 (2019) 383-395; DOI: 10.3221/IGF-ESIS.49.38 384 assessing fatigue strength. Special attention should be paid to the foreign object damage (FOD) research section, which deals with the destruction of engine parts after preliminary dynamic loading [7-11]. At present, after several decades of research, it is becoming clear that fatigue damage is the result of accumulation of microplastic deformations, which lead to irreversible microstructural and topological changes [12] that determine critical states of damage and the transition to destruction. General knowledge of the mechanisms of microstructural changes leading to fatigue failure allows us to associate these critical states with the development of an ensemble of microcracks originating in localized shifts. The interest to fundamental problems of fatigue damage has sharply increased in connection with the possibility of attaining a fatigue resource corresponding to the so-called very high (gigacycle) fatigue [1-10]. The results of these studies raised, in particular, the question of the existence of a fatigue limit in the transition from high cycle to gigacycle loading conditions. The formation of fatigue damage is traditionally associated with microplastic deformations [1-2, 12], which are occur under cyclic loading conditions, initiating various microstructural mechanisms that control durability and depend on the nature and initial structure of the material. For ductile metals, a special type of fatigue (cyclic) deformation localization (persistent slip bands - PSB) is observed, the development of which traditionally leads to the initiation of cracks in the near-surface zone. The second common case, for example, in the fatigue fracture of high-strength steels, is the onset of damage (microcracks), in the vicinity of inclusions at low load amplitudes corresponding to very high cycle fatigue. A characteristic feature of destruction under gigacycle fatigue conditions is the decisive effect of the initiation stage of a fatigue crack on the fatigue life. At the same time, a qualitative difference is the formation of a fatigue crack in the bulk of the material, which decisively changes the formulation of the problem of assessing the fatigue life, methods for studying the stages of fracture development. In contrast to the established traditions in the field of high-cycle fatigue, where central attention is paid to the crack propagation stage, a fundamental problem arises about formation of the fatigue crack during multiscale damage development processes associated with defects of various nature (inclusions, localized plastic shear bands, microcracks, pores). The staging of destruction is characterized by the effects of “irreversibility”, initiated by the formation of localized shifts that play a key role in the initiation of a fatigue crack [13-15], which can manifest itself in signs of nonlinearity of the elastic behavior of materials, “anomalies of elastic compliance” of fatigue samples. The role of the initiation stage is especially important for gigacyclic loading regimes, which are characterized by "fish-eye" crack formed in the bulk of the material [13-14]. This article presents a method for monitoring the kinetics of accumulation of defects by analyzing the amplitude of the second harmonic of oscillations of the free end of the sample. With the appearance of internal defects, nonlinearity effects arise, which significantly increase with their correlated interaction and the formation of a fatigue crack. MATERIAL AND EXPERIMENT CONDITION n this research were investigated samples of aluminum alloy AMG-6, the chemical composition and mechanical characteristics of which are presented in Tabs. 1 and 2. The material was subjected to dynamic deformation performed by split Hopkinson pressure bar set-up at the strain rates of ~103 s-1 [16]. The material was pre-deformed in order to specify a certain degree of initial damage to assess further of the defects development kinetics. Samples with a geometry corresponding to Fig. 1. Al Cu Mg Mn Si Fe Zn Be Ti 91.1-93.68 0.10 5.8-6.8 0.5-0.8 0.4 0.4 0.20 0.0002-0.005 0.02-0.1 Table 1: The chemical composition of AMg6 (as a percentage) Elastic modulus (GPa) Yield strength (MPa) Tensile strength (MPa) Maximum elongation (%) 71 180 355 25 Table 2: Quasi-static tensile properties of AMg6 I M. Bannikov et alii, Frattura ed Integrità Strutturale, 49 (2019) 383-395; DOI: 10.3221/IGF-ESIS.49.38 385 Figure1: Geometry of the samples (scales in mm) The fatigue loading was carried out on an ultrasonic Shimadzu USF-2000 resonant-type testing machine (Fig. 2) at stress levels of 90-162 MPa and a symmetrical cycle R = -1. The principle of operation testing machine is based on formation of standing wave in the generator-horn-sample system [2]. In this case, the maximum amplitude of displacements occurs at the ends of the sample, and the maximum amplitude of stresses at its center. On the side of the free end of the specimen, an inductive sensor is fixed, which makes possible to determine the displacement of the end with an accuracy of tenths of a micron and a frequency of up to 1 MHz. The machine allows testing materials on the basis of 109-1010 cycles with amplitude from the 1st to several tens of microns with a frequency of 20 kHz, which reduces the test time to several days, in contrast to the classical fatigue installations, where this number of cycles is achieved during the years of testing. Figure 2: The experimental setup: 1 - Horn, 2 - sample, 3 - displacement sensor, 4 - cooling system. The formation of cracks in cylindrical samples undergoing high cycle loading in the range 106-107 cycles begins from the surface. For the dynamic preloading samples AlMg6 subject to cyclic loads in the range exceeding 108 cycles the crack formation is started in the bulk of the sample. The fracture surface in this case exhibits a fatigue zone known as a “fish- eye” which is a distinguishing feature of such fatigue regimes. The central part of this region comprises a fracture nucleus surrounded by the region of refined (submicrocrystalline) structure so-called FGA [13] (Fig. 3a). M. Bannikov et alii, Frattura ed Integrità Strutturale, 49 (2019) 383-395; DOI: 10.3221/IGF-ESIS.49.38 386 a) (b) Figure 3: (a) Characteristic surface relief of a gigacycle fatigue fracture zone of AlMg6, (b) scaled-up fragment of the «fish-eye» zone obtained by the scanning electron microscope Hitachi S3400 (σ=138 MPa Nf=7.51108), the region 1 is the center of crack initiation, the region 2 is FGA zone. NONLINEAR ACOUSTICS METHOD FOR NON-DESTRUCTIVE TESTING OF FATIGUE FAILURE he method is based on the initiation of a longitudinal finite-amplitude perturbation A0 with a frequency ω0 on one side of the sample, while its other end remains free [17-19]. Oscillations of the free end of the sample will contain a number of harmonic components: components with amplitude A1 at the fundamental frequency ω0, amplitude A2 of the second harmonic with frequency 2ω0, and so on. The nonlinearity parameter βe is determined experimentally by measuring the absolute amplitudes of the signals of the first A1 and second A2 harmonics corresponding to the nonlinear law of elasticity: 2 21 1 ... ...2 3 22 2 u u u ue e e eA A A a a a a                                              (1) where σ is the stress, u is the displacement, a is the spatial coordinate, 2 eA and 3 eA - are the second and third order elastic coefficients, respectively. By introducing the nonlinearity coefficient: 3 2( / ) e e eA A   and the wave equation can be represented as: 2 2 2 1 e u u u c a at                     (2) where u is the component of the displacement vector in the direction, c - is the longitudinal speed of sound, t - is time. Its solution, given that the end oscillation u =u1cos(ωt), will be: 0 1 2cos( ) sin 2( ) ...u u u t u t ka      (3) That the amplitude of the second harmonic 2 22 1(1/ 8) eu k u a , where 0 /k   – wave number, we can derive: 2 22 18 / e u k u a  (4) In the study of nonlinear phenomena in the gigaclic fatigue regime, the relative parameter is determined by measuring the amplitudes of the main and second harmonics: T M. Bannikov et alii, Frattura ed Integrità Strutturale, 49 (2019) 383-395; DOI: 10.3221/IGF-ESIS.49.38 387 relative 0/   (5) where β0 refers to intact material. MEASUREMENT OF NONLINEARITY COEFFICIENT o determine the nonlinearity coefficient, we prologue both parts of the Eqn. (4): 2 1log( ) log( ) log( ) 2 log( )k u u    (6) The amplitude of the oscillations is measured in decibels at scales  1 20 1A log u and (2 20 2)A log u , then the ratio takes the form: 2 120 log( ) 20 log( ) 20 log( ) 40 log( ),k u u    (7) 2 120 log( ) 2K A A    (8) Thus, the relative nonlinearity parameter β can be found from the formula: 0 2 1 2 1 020 log( / ) ( 2 ) ( 2 )A A A A      (10) The authors developed a software and hardware complex which analyzes the signal in real time with a variable duration and performs a fast Fourier transform to search for the amplitudes of the fundamental and second harmonics. When a trend is detected to increase the amplitude of the second harmonic, the software stops the experiment and reports the possible occurrence of a fatigue crack. a) b) Figure 4: The change in the non-linearity factor βrelative (dB) during testing of the sample with an internal crack (a) σ = 138 MPa, Nf = 7.51 • 108 (Fish-eye) and a sample with a surface crack (b) σ = 140 MPa, Nf = 2.72 • 108. The dependences of the non-linearity coefficient βrelative on the loading amplitude during fatigue life for aluminum alloy AMg-6 are given. During the experiment, βrelative increases slightly at the initial stage, which is usually attributed to the effect of heating the sample and, as a result, changing its linear dimensions [17-19], then most of the time of the experiment, the amplitude of the second harmonic remains constant and increases significantly with the formation of a fatigue crack and its growth (Fig. 4, a). The amplitude of the second harmonic is associated with the deformation due to defects. For given loading parameters, the system is macroscopically in the elastic region, and all deformations are T M. Bannikov et alii, Frattura ed Integrità Strutturale, 49 (2019) 383-395; DOI: 10.3221/IGF-ESIS.49.38 388 reversible. The defects weakly correlate with each other, and the amplitude of the second harmonic is small and varies little. However, in the process of cyclic loading, irreversible microdamages are appeared, which begin to interact with each other, which subsequently leads to the formation of a fatigue crack. In one of the samples, which were destroyed with the formation of a surface crack, a monotonic decrease in the amplitude of the second harmonic uncharacteristic for others was found during the testing process (Fig. 4.b). This is supposed to be due to the fact that the initial structure of the material was saturated with defects, since the samples were pre-deformed. During cyclic loading, a small amplitude in the material resulted in relaxation of internal stresses [20], which reduced the nonlinearity of the signal. When a stress relaxation resource is depleted, an avalanche-like growth of the second harmonic amplitude is observed, which indicates a high correlation of defects and their interaction on a large scale compared with the internal nucleation of a fatigue crack. In the fig. 5 the growth of the amplitude of the second harmonic is shown. At the time point of 5.415 * 10 ^ 4, the software stopped the experiment on the basis of a significant increase in the amplitude of the second harmonic (6 times), which indicated the formation of a fatigue crack. Later, the experiment was started again, which allowed us to observe how the amplitude of the second harmonic increases with increasing crack length. Thus, the βrelative parameter can be used not only as an indicator of the appearance of a fatigue crack, but also to determine its size, as it was done in [21]. Figure 5. The amplitude of the second harmonic of prestrained sample AMG-6. The difference in the number of cycles between the moment of detection of the nucleation of a fatigue crack and its exit to macroscopic destruction (the output of the system from resonance) was more than 2 • 106 cycles. MATHEMATICAL MODEL he wide-range constitutive equations [22-24], describing the relationship between the kinetics of damage caused by defects and the relaxation properties of materials, were used to construct a mathematical model that is capable of describing the deformation behavior of metals and alloys under fatigue loading. To describe the evolution of the structure of the material, two parameters are introduced: tensor: ( ) 2 s  s lb bl and scalar: 3 0 ( ) R δ r  where: s - tensor of microshears [23]; s - shear intensity; l - unit normal vector to the shear plane; b - unit vector in the direction of shear; δ - scale-invariant, structural parameter; R - the distance between defects, 0r - the characteristic size of the nuclei of defects. T M. Bannikov et alii, Frattura ed Integrità Strutturale, 49 (2019) 383-395; DOI: 10.3221/IGF-ESIS.49.38 389 Averaging s over an elementary volume gives a tensor of microshears density [23]: ( , , )W s dV p s l b (10) where: V - volume, ( , , )W s l b - distribution function of orientation and intensity of microshears. In its physical meaning p is a strain due to defects. The total strain rate ( ε ) consists of three components: plastic ( pε ), elastic ( eε ) and due to defects ( p ): e p  ε ε ε p    From the second law of thermodynamics, it follows that the energy dissipation can be represented as: : : 0p F F TS δ δ         σ ε p p    , (11) where: T - temperature; S - rate of change of entropy; σ - stress tensor; - no equilibrium free energy. According to the Onsager principle, еhe following relations are obtained from (11): 1 2 pl l σ ε p  (12) 2 3 pF l l       ε p p   (13) 4 F l δ δ      (14) where: 1l , 2l , 3l , 4l - kinetic coefficients, in the general case, depending on state parameters, satisfying the constraint: 2 1 3 2 0l l l  . Relations (12)-(14) are complemented by Hooke's law in the rate formulation and approximation for non-equilibrium free energy. Thus, the complete system of constitutive equations looks as follows: ( : ) 2 ( )pλ G   σ D E E ε ε p    (15) 1 2 p d F      ε σ p  (16) 2 3d F      p σ p  (17) 4 F δ δ      (18) 2 2 2 1 2 3 4 : ln( ) 2 2 2 d m p pF c p c c c p p F δ G        σ p (19) d s σ σ σ , 1 ( : ) 3 s σ σ E E M. Bannikov et alii, Frattura ed Integrità Strutturale, 49 (2019) 383-395; DOI: 10.3221/IGF-ESIS.49.38 390 where: 31 2 1 3 2 l l l l    , 22 2 1 3 2 l l l l    , 13 2 1 3 2 l l l l    , 4 4 1 l   ; E - the unit tensor, dσ and sσ - the deviator and spherical parts of the stress tensor; mF , 1c - 4c – potential approximation constants, G – shear modulus; λ – the first parameter Lame; :p  p p . In general, kinetic coefficients depend on state parameters and are representable as follows: 1 i iG   where i - the characteristic relaxation times, which in general can be represented as dependencies: 0 ( , , , , , , , , )exp( ) p p i i i U δ T kT    σ ε ε ε ε p p   , where: k is the Boltzmann constant, ( )iU  is the characteristic activation energy. The hypothesis is accepted that for the studied processes the relaxation times have the following dependencies: 0 0 0exp( ) in i i i U ε U kT       , 1, 2, 3i  (20) 0 4 0 0 4 4 exp( ) nU p U kT        (21) where: 0 2 : 3 c ε ε  ε ε    , 1cε  s-1 is the dimensionless factor; 0U kT , iU , 1...3i  are the constants that have the meaning of the characteristic energy of overcoming barriers, after which new relaxation mechanisms are activated; 0 : c p ε  p p   ; 1 2 3, , ,n n n n - constants, characterizing the rate sensitivity of the material. The meaning of this hypothesis is that the physicomechanical characteristics of the material depend on the strain rate, and the evolution of damage under cyclic loading with small stress amplitudes depends only on the parameter, which in this case can be interpreted as the growth rate of a fatigue crack, related to its characteristic scale. Assuming a weak strain rate sensitivity of the relaxation times, we can expand the expression (20) into a Taylor series with an accuracy of the second term: 0 0 0 0 0 0 0exp( ) 1 1 i i i in n n n i i i iU ε U U ε U U ε U U ε kT kT kT kT kT             then we get: 0 0 0 0 0 1 i i i U i i i i n n n i i kT U ε U ε ε kT          Similarly, for the expression (21): M. Bannikov et alii, Frattura ed Integrità Strutturale, 49 (2019) 383-395; DOI: 10.3221/IGF-ESIS.49.38 391 0 0 0 04 0 0 4 4 4 4 4 4 4 0 0 4 0 0 4 0 0 1 1 exp( ) exp( ) 1 n U n n n n U p U kT kT U p U U p U U p p kT kT                           . Then the kinetic coefficients will have the following form: 0 in i U i ε G    , 1...3i  , 04 4 n U p G      . We introduce the notation: 1 1 0 nε    , 22 0 n p ε    , 33 0 n pε    , 4 0 np     , 1 2n n n   , 2 3pn n n  . The final determining equation is the fracture criterion. In the framework of the proposed model, we can enter a criterion based on the structural scaling parameter: 0δ  (22) The meaning of this criterion is that the percentage ratio of the volume of material occupied by defects tends to 100%, which can be achieved by tending to zero the distance between defects or tending to infinity of the size of defects. The second option is impossible, and the first has a clear physical meaning. In the calculation, the condition (22) is not realizable due to the instability, as it is, therefore, we can replace it with a softer one: fδ δ (23) where 0.4fδ  is the critical value [22], after which the avalanche-like growth of defects begins and the material is considered destroyed. This value is universal for all plastic materials for the proposed model. Critical value fδ δ corresponds to critical value cp . Model Parameter Identification The identification of the parameters of the constructed model (15)-(19) was carried out in the uniaxial case: 2 ( )pσ λ G p         (24) 0 ( ) np p F ε σ p            (25) 0 ( ) pn p p F p ε σ p          (26) 0 n Fδ p δ          (27) 2 2 2 1 2 3 4ln( ) 2 2 2m p p σpF c p c c c p p F δ G        (28) M. Bannikov et alii, Frattura ed Integrità Strutturale, 49 (2019) 383-395; DOI: 10.3221/IGF-ESIS.49.38 392 where:  , σ , pl , p are the corresponding components of the tensors ε , σ , pε , p . At the first stage of identification, the problem of minimizing the discrepancy between the experimental and calculated stress-strain diagrams is solved. This determines the parameters  , p , p ; dimensionless strain rate 0 1ε  (static loading), 1.15δ const  . At the second stage, using the known values  , p , p the problem of minimizing the discrepancy between the experimental and calculated stress-strain diagrams at different strain rates is solved. This determines the parameters n and pn , 1.15δ const  . An illustration of optimization problem solving is presented in Tab. 3. Strain rate, s-1 Yield strength, MPa (the experiment, [25, 26]) Yield strength, MPa (the calculation) 0.0001 165 166 520 175 176 1210 210 213 Table 3: The yield strength in the experiment [25, 26] and the calculation at different strain rates. The final step in identifying the parameters of the model is to determine the values of  and n that are found using experimental fatigue loading data. Two characteristic points on the Wöhler curve were selected: (300,000 cycles, 190 MPa) and (300,000,000 cycles, 156 MPa). Parameters  and n were chosen in such a way that at the loading amplitude of 190 MPa the material was destroyed after 300,000 cycles, and at 156 MPa - after 300,000,000 cycles, respectively. The frequency of loading in this case was equal to 20 000 Hz, as in the experiment. Thus, the complete set of constants for the AlMg6 alloy is as follows: Constants known from literature [24]: 2670ρ  kg / m3, 41  GPa, 27G  GPa. Constants that were defined: 529  (Pa·s)-1, 38.5p  (Pa·s)-1, 2.9p  (Pa·s)-1, 3.46  (Pa·s)-1, 0.967pn n   , 1.985n  . Numerical simulation results A numerical study of the behavior of a real structure and even a laboratory specimen in the very high cycle fatigue mode is not possible due to the requirement of huge computational and time resources. Therefore, instead of solving the boundary value problem, the problem of loading a representative material volume (24)-(28) with loading conditions in the form of applied cyclic stresses (29) with amplitude A and frequency 20000  Hz, as well as initial conditions (21) is considered. ( ) sin( 2 )At t   (29) 00 00 0 0, 0, 0, 1.15p t tt t p δ δ          (30) The loading condition (19) can be written in terms of strain: ( ) ( )sin( 2 )Aε t ε σ t (31) where ( )ε σ is the function that converts stresses into strains for AlMg6 alloy. Differentiating (31) by time, we get: ( ) ( )2 cos( 2 )Aε t ε σ t  (32) The numerically constructed Wöhler curve and its comparison with experimental data of authors and experiment [27] with similar experiment conditions (AlMg-6, R=-1, frequency 10 kHz) are presented in Fig. 6. M. Bannikov et alii, Frattura ed Integrità Strutturale, 49 (2019) 383-395; DOI: 10.3221/IGF-ESIS.49.38 393 Figure 6: Wöhler curve for AlMg6 alloy. Solid line - numerical calculation, points - experimental data of the authors, squares - experiment [27] with conditions: frequency 10 kHz, R=-1. Due to the amplitudes of the applied stresses below the yield strength, the material is macroscopically in the elastic region, and the applied deformations can be considered small. Then expression (24) can be rewritten in the form: 2 ( )pσ λ G p      (33) Since macroplastic strain is not observed, then in the very high cycle fatigue mode, the parameter p is responsible for the accumulation of damage and destruction. The growth of defects leads to a gradual change in the amplitude of the applied voltages, and consequently to a change in the amplitude of the applied displacements, since for some displacement amplitude AU in the loading machine is always assigned Aσ : ( ) 2 ( ( ) )A A Aσ λε σ G ε σ p   (34) After reaching the parameter p of some critical value pc (which corresponds to 0, 4fδ  ), an avalanche-like growth of defects occurs, leading to the destruction of the material, which in the experiment corresponds to an avalanche-like change in the amplitude of the second harmonic of displacements (Fig. 7.). Figure 7: Defect kinetics (dashed line) in the calculation and change in the amplitude of the second harmonic of displacements in the experiment (solid line) . In fig. 7 shows the normalized curves of variation of the parameter p, referred to pc and the amplitude of the second harmonic of displacements, referred to its maximum value. Figure 7 shows the experimental data corresponding to the preloaded specimen. According to the proposed model, the change in structure due to preloading is described in terms of the parameter δ. Namely, a smaller value of the initial condition δ0 corresponds to the preloaded specimen. To obtain the dotted curve (Fig. 7), the Eqns. (24-29, 30, 32) was solved with the initial condition δ0 = 1.02. Thus, the proposed model is able to adequately describe the behavior of the material in the modes of high- and very high cycle fatigue. Based on the results of numerical simulation, we can conclude that the change in the amplitude of the M. Bannikov et alii, Frattura ed Integrità Strutturale, 49 (2019) 383-395; DOI: 10.3221/IGF-ESIS.49.38 394 second harmonic of displacements in the experiment gives adequate information about the accumulation of damage in the material. CONCLUSION n analysis of the free end oscillations during gigacyclic tests showed a significant increase in the second harmonic amplitude and nonlinearity coefficient during the formation of fatigue cracks. In conjunction with the simulation results, we can conclude that the deviation from the linear-elastic law (“anomaly of elastic compliance”) is due to the accumulation of defects in the material, and the avalanche-like growth of the second harmonic displacement corresponds to the avalanche-like growth of defects and the formation of a macroscopic crack. Thus, it is possible to judge the accumulation of damage, including - in the amount of material based on the acoustic signal. Using the proposed equations, it can to predict the fatigue life of the material and identify internal defects and cracks before macroscopic destruction occurs. ACKNOWLEDGMENTS uthors thank Prof. Valery Matveyenko for the proposal to present the paper in the issue. This work was supported by the Russian Science Foundation, grant number 18-72-00142. REFERENCES [1] Botvina, L.R. (2004) Gigacycle fatigue is a new problem in physics and mechanics of destruction // Factory Laboratory. Diagnostics of materials 70(4), pp. 41. (In Russian). [2] Bathias, C., Paris, P.C. (2005), Gigacycle Fatigue in Mechanical Practice, Marcel Dekker Publisher Co., pp. 328. [3] Oborin, V., Bannikov, M., Naimark, O., Palin-Luc, T. (2010). Scaling invariance of fatigue crack growth in gigacycle loading regime, Technical Physics Letters, 36 (11), pp. 1061-1063. DOI: 10.1134/S106378501011026X. [4] Cowles, B.A. (1996). High cycle fatigue in aircraft gas turbines - an industry perspective, International Journal of Fracture 80, pp.147-163. DOI: 10.1007/BF00012667. [5] Shanyavsky, A.A. (2007) Simulation of fatigue damage of metals. Synergetic in aviation, Ufa: Monograph LLC. (In Russian). [6] Nicholas, T. (2006). High Cycle Fatigue. A Mechanics of Material Perspective, Elsevier., pp. 641. [7] Peters, J.O., Ritchie, R.O. (2000). Influence of foreign object damage on crack initiation and early crack growth during high-cycle fatigue of Ti-6Al-4V, Eng. Fract. Mech. 67, pp. 193-207. DOI: 10.1016/S0013-7944(00)00045-X. [8] Spanrad, S., Tong, J. (2011). Characterisation of foreign object damage (FOD) and early fatigue crack growth in laser shock peened Ti–6Al–4V aerofoil specimens, Materials Science and Engineering A, 528, pp. 2128–2136. DOI: 10.1016/j.msea.2010.11.045. [9] Oakley, S.Y., Nowell, D. (2007) Prediction of the combined high- and low-cycle fatigue performance of gas turbine blades after foreign object damage, International Journal of Fatigue, 29, pp. 69–80. DOI: 10.1016/j.ijfatigue.2006.02.042. [10] Chen, Xi (2005). Foreign object damage on the leading edge of a thin blade, Mechanics of Materials, 37, pp. 447–457. DOI: https://doi.org/10.1016/j.mechmat.2004.03.005. [11] Nowell, D., Duó, P., Stewart, I.F. (2003). Prediction of fatigue performance in gas turbine blades after foreign object damage, International Journal of Fatigue, 25, pp. 963-969. DOI: 10.1016/S0142-1123(03)00160-9. [12] Mughrabi, H. (2015). Microstructural mechanisms of cyclic deformation, fatigue crack initiation and early crack growth, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences, 373(2038), 20140132. DOI: 10.1098/rsta.2014.0132. [13] Zhang, Li-Li, et al. (2016). On the formation mechanisms of fine granular area (FGA) on the fracture surface for high strength steels in the VHCF regime, International Journal of Fatigue, 82, pp. 402-410. DOI: 10.1016/j.ijfatigue.2015.08.021. A A M. Bannikov et alii, Frattura ed Integrità Strutturale, 49 (2019) 383-395; DOI: 10.3221/IGF-ESIS.49.38 395 [14] Oborin, V., Sokovikov, M., Bilalov, D., & Naimark, O. (2016). Multiscale study of morphology of the fracture surface aluminum-magnesium alloy with consecutive dynamic and gigacycle loading, Procedia Structural Integrity, 2, pp. 1063- 1070. DOI: 10.1016/j.prostr.2016.06.136. [15] Froustey, C., Naimark, O., Bannikov, M., Oborin, V. (2010). Microstructure scaling properties and fatigue resistance of pre-strained aluminium alloys (part 1: AlCu alloy), European Journal of Mechanics A/Solids, 29, pp. 1008-1014. DOI: 10.1016/j.euromechsol.2010.07.005. [16] Oborin, V.A., Bayandin, Yu. V., Bilalov, D. A., Sokovikov, M. A., Chudinov, V. V, Naimark, O. B. (2018). Self-similar laws of damage development and evaluation of the reliability of alloys D16T and AMg6 under combined dynamic and gigacycle loading, Phys. Mezomekh, 21(6), pp. 135-145 DOI: 10.1134/S1029959919020048. [17] John H. Cantrell, William T. Yost (2001) Nonlinear ultrasonic characterization of fatigue microstructures, Int. J. of Fatigue, 23, pp.487–490. DOI: 10.1016/S0142-1123(01)00162-1. [18] Kumar, A., Torbet, J.C. Pollock, M.T., Jones, W.J. (2010). In situ characterization of fatigue damage evolution in a cast Al alloy via nonlinear ultrasonic measurements, ActaMaterialia, 58(6), pp. 2143-2154. DOI: 10.1016/j.actamat.2009.11.055 [19] Kumar A. et al. (2011) In situ damage assessment in a cast magnesium alloy during very high cycle fatigue. ScriptaMaterialia, 64(1), pp.65-68. DOI: 10.1016/j.scriptamat.2010.09.008 [20] Nazarov A.A., (2018). Nonequilibrium grain boundaries in bulk nanostructured metals and their recovery under the influences of heating and cyclic deformation. Review, Letters on Materials 8 (3), 2018 pp. 372-381 (in Russian). DOI: 10.22226/2410-3535-2018-3-372-381 [21] W. Li, H. Cui, W. Wen, X. Su, C. C. Engler-Pinto Jr.: (2015). In situ Nonlinear Ultrasonic for Very High Cycle Fatigue Damage Characterization of a Cast Aluminum Alloy. Materials Science and Engineering A, 645, pp. 248-254. DOI: 10.1016/j.msea.2015.08.029. [22] Bilalov, D.A., Bayandin, Yu.V., Naimark, O.B. (2018). Mathematical modeling of failure process of AlMg2.5 alloy during high- and very high cycle fatigue, Computational continuum mechanics. 11(3), pp. 323-334. (In Russian). DOI: 10.7242/1999-6691/2018.11.3.24. [23] Naimark, O.B. (2003) Collective properties of defects ensembles and some nonlinear problems of plasticity and fracture, Physical mesomechanics, 6(4), pp. 39-63. [24] Froustey C., Naimark O.B., Panteleev I.A., Bilalov D.A., Petrova A.N., Lyapunova E.A. (2017) Multiscale structural relaxation and adiabatic shear failure mechanisms, Physical Mesomechanics, 20(1), pp. 31-42. DOI: 10.1134/S1029959917010039. [25] Glushak, B.L., Ignatova, O.N., Pushkov, V.A., Novikov, S.A., Girin, A.S., Sinitsyn, V.A. (2000). Dynamic Deformation of Aluminum Alloy AMg-6 at Normal and Higher Temperatures, Journal of Applied Mechanics and Technical Physics., 41(6), pp. 1083-1086. DOI: 10.1023/A:1026662824249. [26] Frolov, K.V. (2001). Mechanical Engineering. Encyclopedia. Volume II-3: Non-ferrous metals and alloys. Composite metallic materials. Moscow, MechanicalEngineering, 880 p. (In Russian). [27] Yakovleva, T. Yu., Matokhnyuk, L.E. (2004). Prediction of fatigue metal resistance characteristics at different loading frequencies, Strength problems, 4, pp. 145-155. (In Russian). << /ASCII85EncodePages false /AllowTransparency false /AutoPositionEPSFiles true /AutoRotatePages /None /Binding /Left /CalGrayProfile (Dot Gain 20%) /CalRGBProfile (sRGB IEC61966-2.1) /CalCMYKProfile (U.S. Web Coated \050SWOP\051 v2) /sRGBProfile (sRGB IEC61966-2.1) /CannotEmbedFontPolicy /Error /CompatibilityLevel 1.4 /CompressObjects /Tags /CompressPages true /ConvertImagesToIndexed true /PassThroughJPEGImages true /CreateJobTicket false /DefaultRenderingIntent /Default /DetectBlends true /DetectCurves 0.0000 /ColorConversionStrategy /CMYK /DoThumbnails false /EmbedAllFonts true /EmbedOpenType false /ParseICCProfilesInComments true /EmbedJobOptions true /DSCReportingLevel 0 /EmitDSCWarnings false /EndPage -1 /ImageMemory 1048576 /LockDistillerParams false /MaxSubsetPct 100 /Optimize true /OPM 1 /ParseDSCComments true /ParseDSCCommentsForDocInfo true /PreserveCopyPage true /PreserveDICMYKValues true /PreserveEPSInfo true /PreserveFlatness true /PreserveHalftoneInfo false /PreserveOPIComments true /PreserveOverprintSettings true /StartPage 1 /SubsetFonts true /TransferFunctionInfo /Apply /UCRandBGInfo /Preserve /UsePrologue false /ColorSettingsFile () /AlwaysEmbed [ true ] /NeverEmbed [ true ] /AntiAliasColorImages false /CropColorImages true /ColorImageMinResolution 300 /ColorImageMinResolutionPolicy /OK /DownsampleColorImages true /ColorImageDownsampleType /Bicubic /ColorImageResolution 300 /ColorImageDepth -1 /ColorImageMinDownsampleDepth 1 /ColorImageDownsampleThreshold 1.50000 /EncodeColorImages true /ColorImageFilter /DCTEncode /AutoFilterColorImages true /ColorImageAutoFilterStrategy /JPEG /ColorACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /ColorImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000ColorACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000ColorImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasGrayImages false /CropGrayImages true /GrayImageMinResolution 300 /GrayImageMinResolutionPolicy /OK /DownsampleGrayImages true /GrayImageDownsampleType /Bicubic /GrayImageResolution 300 /GrayImageDepth -1 /GrayImageMinDownsampleDepth 2 /GrayImageDownsampleThreshold 1.50000 /EncodeGrayImages true /GrayImageFilter /DCTEncode /AutoFilterGrayImages true /GrayImageAutoFilterStrategy /JPEG /GrayACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /GrayImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000GrayACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000GrayImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasMonoImages false /CropMonoImages true /MonoImageMinResolution 1200 /MonoImageMinResolutionPolicy /OK /DownsampleMonoImages true /MonoImageDownsampleType /Bicubic /MonoImageResolution 1200 /MonoImageDepth -1 /MonoImageDownsampleThreshold 1.50000 /EncodeMonoImages true /MonoImageFilter /CCITTFaxEncode /MonoImageDict << /K -1 >> /AllowPSXObjects false /CheckCompliance [ /None ] /PDFX1aCheck false /PDFX3Check false /PDFXCompliantPDFOnly false /PDFXNoTrimBoxError true /PDFXTrimBoxToMediaBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXSetBleedBoxToMediaBox true /PDFXBleedBoxToTrimBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXOutputIntentProfile () /PDFXOutputConditionIdentifier () /PDFXOutputCondition () /PDFXRegistryName () /PDFXTrapped /False /CreateJDFFile false /Description << /ARA /BGR /CHS /CHT /CZE /DAN /DEU /ESP /ETI /FRA /GRE /HEB /HRV (Za stvaranje Adobe PDF dokumenata najpogodnijih za visokokvalitetni ispis prije tiskanja koristite ove postavke. Stvoreni PDF dokumenti mogu se otvoriti Acrobat i Adobe Reader 5.0 i kasnijim verzijama.) /HUN /ITA /JPN /KOR /LTH /LVI /NLD (Gebruik deze instellingen om Adobe PDF-documenten te maken die zijn geoptimaliseerd voor prepress-afdrukken van hoge kwaliteit. De gemaakte PDF-documenten kunnen worden geopend met Acrobat en Adobe Reader 5.0 en hoger.) /NOR /POL /PTB /RUM /RUS /SKY /SLV /SUO /SVE /TUR /UKR /ENU (Use these settings to create Adobe PDF documents best suited for high-quality prepress printing. Created PDF documents can be opened with Acrobat and Adobe Reader 5.0 and later.) >> /Namespace [ (Adobe) (Common) (1.0) ] /OtherNamespaces [ << /AsReaderSpreads false /CropImagesToFrames true /ErrorControl /WarnAndContinue /FlattenerIgnoreSpreadOverrides false /IncludeGuidesGrids false /IncludeNonPrinting false /IncludeSlug false /Namespace [ (Adobe) (InDesign) (4.0) ] /OmitPlacedBitmaps false /OmitPlacedEPS false /OmitPlacedPDF false /SimulateOverprint /Legacy >> << /AddBleedMarks false /AddColorBars false /AddCropMarks false /AddPageInfo false /AddRegMarks false /ConvertColors /ConvertToCMYK /DestinationProfileName () /DestinationProfileSelector /DocumentCMYK /Downsample16BitImages true /FlattenerPreset << /PresetSelector /MediumResolution >> /FormElements false /GenerateStructure false /IncludeBookmarks false /IncludeHyperlinks false /IncludeInteractive false /IncludeLayers false /IncludeProfiles false /MultimediaHandling /UseObjectSettings /Namespace [ (Adobe) (CreativeSuite) (2.0) ] /PDFXOutputIntentProfileSelector /DocumentCMYK /PreserveEditing true /UntaggedCMYKHandling /LeaveUntagged /UntaggedRGBHandling /UseDocumentProfile /UseDocumentBleed false >> ] >> setdistillerparams << /HWResolution [2400 2400] /PageSize [612.000 792.000] >> setpagedevice