METHOD OF PREDICTING LOW-CYCLE FATIGUE CRACK INITIATION AND PROPAGATION BEHAVIORS UNDER MULTI-SCALE FRAMEWORK

Abstract

A method of predicting low-cycle fatigue crack initiation and propagation behaviors under a multi-scale framework includes the following steps: S1, providing a calculation method for low-cycle fatigue crack initiation and propagation damages under a multi-scale framework; S2, determining a slip system where a maximum damage is located and an accumulated damage of all slip systems by calculation using the calculation method in S1; S3, a crack initiating and propagating in a direction towards the slip system where the maximum damage is located when the accumulated damage reaches a critical value; and S4, conducting calculation repeatedly until a predicted crack length reaches a fracture length of a low-cycle fatigue specimen under test conditions.

Claims

1. A method of predicting low-cycle fatigue crack initiation and propagation behaviors under a multi-scale framework, characterized by comprising the following steps: S1, providing a calculation method for low-cycle fatigue crack initiation and propagation damages under a multi-scale framework: ${\stackrel{.}{d}}_{\mathrm{initial}}={.Math.}_m\frac{2{G\left({\gamma }_{n,m,e}\right)}^2}{\pi \left(1-v\right).Math.d^3.Math.w_{m,\mathrm{critical}}};$ ${\stackrel{.}{d}}_{\mathrm{growt}▯}={.Math.}_m{\mathrm{\xi \lambda }}_d\frac{\int {\tau }_{n,m,e}d{\gamma }_{n,m,e}}{w_{m,\mathrm{critical}}};$ wherein {dot over (d)}.sub.initial is a damage rate of low-cycle fatigue crack initiation, {dot over (d)}.sub.growth is a damage rate of low-cycle fatigue crack propagation, m is a number of slip systems, G is a shear modulus, v is a Poisson's ratio, d is an average grain diameter, γ.sub.n,m,e is an effective shear strain on the slip systems m, τ.sub.n,m,e is an effective shear stress on the slip systems m, w.sub.m,critical is a fracture energy corresponding to each slip system, ξ is a material parameter, and λ.sub.d is a mean free path of dislocations; S2, determining a slip system where a maximum damage is located and an accumulated damage of all the slip systems m by calculation using the calculation method in S1; S3, a crack initiating and propagating in a direction towards the slip system where the maximum damage is located when the accumulated damage reaches a critical value; and S4, conducting calculation repeatedly until a predicted crack length reaches a fracture length of a low-cycle fatigue specimen under test conditions.

2. The method of predicting the low-cycle fatigue crack initiation and propagation behaviors under the multi-scale framework according to claim 1, further comprising: building a low-cycle fatigue finite element model considering a microstructure; and calculating the effective shear stress and effective shear strain on different slip systems of the slip systems m in each grain by using an orientation of each grain and a macroscopic mechanical response of a material.

3. The method of predicting the low-cycle fatigue crack initiation and propagation behaviors under the multi-scale framework according to claim 1, wherein formulas for calculating the effective shear stress and effective shear strain are described below:
τ.sub.n,m,e=(n.sub.m).sup.Tσ.sub.p(n.sub.m),
γ.sub.n,m,e=(n.sub.m).sup.Tε.sub.p(n.sub.m); wherein n.sub.m is a normal vector of the slip systems m, γ.sub.n,m is a shear strain on the slip systems m, ε.sub.p is a macroscopic plastic strain of a material, and σ.sub.p is a macroscopic plastic stress of the material.

4. The method of predicting the low-cycle fatigue crack initiation and propagation behaviors under the multi-scale framework according to claim 1, further comprising: building a fracture energy calculation model under a molecular dynamics system, wherein an XZ plane is defined as a slip plane, a Y direction is defined as a slip direction, and a tensile load at a constant rate is applied in the Y direction.

5. The method of predicting the low-cycle fatigue crack initiation and propagation behaviors under the multi-scale framework according to claim 1, wherein a method of calculating the fracture energy of the slip system is as follows:
w.sub.m,critical=∫.sub.s.sub.1.sup.s.sup.2τ.sub.n,mds.sub.n,m; wherein s.sub.1 is a displacement corresponding to a peak stress; s.sub.2 is a corresponding displacement when a fracture occurs; s.sub.n,m is a tensile displacement on the slip systems m; and τ.sub.n,m is a tensile stress on the slip systems m.

6. The method of predicting the low-cycle fatigue crack initiation and propagation behaviors under the multi-scale framework according to claim 1, wherein the low-cycle fatigue specimen is martensitic heat-resistant steel, and the test conditions comprise: adopting strain loading with a loading waveform being triangular at a target test temperature.

7. The method of predicting the low-cycle fatigue crack initiation and propagation behaviors under the multi-scale framework according to claim 1, further comprising: evaluating reliability of the low-cycle fatigue specimen based on a predicted result obtained by the S4.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0025] The drawings required for use in the description of the embodiments or the prior art will be briefly described hereinafter in order to more clearly explain the embodiments of the present invention or the technical solution in the prior art. It will be apparent that the drawings described herein are only embodiments of the present invention, and other drawings may be obtained from the drawings provided without any creative effort by those of ordinary skill in the art.

[0026] FIG. 1 is a schematic diagram of a research route for a method of predicting low-cycle fatigue crack initiation and propagation behaviors under a multi-scale framework in the present invention.

[0027] FIG. 2 is a fracture energy calculation model under a molecular dynamics system and a tensile stress-displacement relation diagram.

[0028] FIG. 3 is a schematic diagram of a low-cycle fatigue finite element model considering a microstructure.

[0029] FIG. 4 is a schematic diagram of a crack propagation direction criterion in the present invention.

[0030] FIG. 5 is a comparison diagram between predicted results of low-cycle fatigue crack length and morphology proposed by the present invention and test results.

[0031] FIG. 6 is a comparison diagram between predicted results of a low-cycle fatigue crack fracture slip system family proposed by the present invention and test results.

DETAILED DESCRIPTION

[0032] The technical solutions in the embodiments of the present invention will be clearly and completely described hereafter in conjunction with the accompanying drawings in the embodiments of the present invention, and it will be apparent that the described embodiments are only a part of the embodiments of the present invention rather than all embodiments. Based on the embodiments in the present invention, all other embodiments obtained by those of ordinary skill in the art without any creative effort are within the scope of protection of the present invention.

[0033] Referring to FIG. 1, the present invention provides a method of predicting low-cycle fatigue crack initiation and propagation behaviors under a multi-scale framework, comprising: a calculation method for low-cycle fatigue crack initiation and propagation damages under a multi-scale framework, a fracture energy calculation method under a molecular dynamics system, a low-cycle fatigue stress-strain calculation method considering a microstructure, and a crack renewal method based on a low-cycle fatigue condition under the multi-scale framework. The present invention will make further explanation by a low-cycle fatigue test and finite element simulation, wherein G115 martensitic heat-resistant steel is used as a test material, the test temperature is 650° C., the low-cycle fatigue test adopts strain loading, and a loading waveform is triangular.

[0034] In step (1), the calculation method for the low-cycle fatigue crack initiation and propagation damages under the multi-scale framework is provided.

[00002]${\stackrel{.}{d}}_{\mathrm{initial}}=\underset{m}{.Math.}\frac{2{G\left({\gamma }_{n,m,e}\right)}^2}{\pi \left(1-v\right).Math.d^3.Math.w_{m,\mathrm{critical}}}$${\stackrel{.}{d}}_{\mathrm{growth}}=\underset{m}{.Math.}\xi {\lambda }_d\frac{\int {\tau }_{n,m,e}d{\gamma }_{n,m,e}}{w_{m,\mathrm{critical}}}$

[0035] wherein, {dot over (d)}.sup.initial is a damage rate of low-cycle fatigue crack initiation, {dot over (d)}.sub.growth is a damage rate of low-cycle fatigue crack propagation, m is a number of slip systems, G is a shear modulus, v is a Poisson's ratio, d is an average grain diameter, γ.sub.n,m,e is an effective shear strain on slip system m, τ.sub.n,m,e is an effective shear stress on slip system m, w.sub.m,critical is a fracture energy corresponding to each slip system, ξ is a material parameter, and λ.sub.d is a mean free path of dislocations. Firstly, crack initiation damages begin to accumulate under a fatigue load. When a value of an accumulated low-cycle fatigue crack initiation damage is equal to 1, a crack initiates. After that, a low-cycle fatigue crack propagation damage starts to be calculated. When the value of an accumulated low-cycle fatigue crack propagation damage is equal to 1, the crack propagates.

[0036] In step (2), a fracture energy calculation model under a molecular dynamics system is built, wherein an XZ plane is defined as a slip plane, a Y direction is defined as a slip direction, and a tensile load at a constant rate is applied in the Y direction, as is illustrated in FIG. 2(a). It should be ensured in the fracture energy calculation model that the lengths in X and Y directions are greater than 20 n.sub.m, and the length in a Z direction is greater than 4 times of a lattice constant. Calculation methods for fracture energies of different slip systems are put forward according to the law of a stress-displacement curve, as is illustrated in FIG. 2(b).

w.sub.m,critical=∫.sub.s.sub.1.sup.s.sup.2τ.sub.n,mds.sub.n,m

[0037] wherein s.sub.1 is a displacement corresponding to a peak stress, s.sub.2 is a corresponding displacement when fracture occurs, s.sub.n,m is a tensile displacement on slip system m, and τ.sub.n,m is a tensile stress on slip system m.

[0038] In step (3), a low-cycle fatigue finite element model considering a microstructure is built, as is illustrated in FIG. 3. The model has a length consistent with a diameter of a specimen, and a loading direction thereof is consistent with that in the test. The effective shear stress and effective shear strain on different slip systems in each grain are calculated by using an orientation of each grain and a macroscopic mechanical response of a material.

τ.sub.n,m,e=(n.sub.m).sup.Tσ.sub.p(n.sub.m) γ.sub.n,m,e=(n.sub.m).sup.Tε.sub.p(n.sub.m)

[0039] wherein, n.sub.m is a normal vector of slip system m, γ.sub.n,m is a shear strain on slip system m, ε.sub.p is a macroscopic plastic strain of the material, and σ.sub.p is a macroscopic plastic stress of the material. Effective shear stress/strain is a projection of a plastic stage in macroscopic stress/strain on slip system m.

[0040] In step (4), an accumulated damage of grains under a cyclic load is calculated using the formulas in steps (1) to (3). When a value of the accumulated crack initiation/propagation damage reaches 1, the slip system where a maximum damage is located is determined by using the formula in step (1). A crack usually breaks along the slip system under the low-cycle fatigue condition. Therefore, a predicted crack is renewed along the slip system where the maximum damage is located. FIG. 4 is a schematic diagram of a crack propagation criterion in the method.

[0041] In step (5), steps (1) to (4) are executed repeatedly until a predicted crack length reaches a fracture length of a low-cycle fatigue specimen under test conditions. It is found by observing the fracture of G115 steel in a low-cycle fatigue test at 650° C. that the test stops when the crack propagates to 3 mm. In order to verify the effect of the method of predicting the low-cycle fatigue crack initiation and propagation behaviors based on the multi-scale framework provided by the present invention, low-cycle fatigue crack prediction results, obtained by the method, of the G115 steel at 650° C. are compared with test results, as is illustrated in FIG. 5 and FIG. 6. The prediction results and the test results show high accuracy in characterizing the crack length and a proportion of a fracture slip system family. Therefore, the method of predicting the low-cycle fatigue crack initiation and propagation behaviors under the multi-scale framework provided by the present invention can well predict a low-cycle fatigue fracture mechanism.

[0042] The embodiments in the description are described in a progressive manner and each embodiment focuses on differences from the other ones. Identical and similar parts of the embodiments can be referred to by one another. Apparatuses disclosed in the embodiments are described relatively simple since they correspond to the methods disclosed in the embodiments, and the relevant parts can refer to the description of the methods.

[0043] The above description of the disclosed embodiments enables those skilled in the art to practice or use the present invention. Various modifications to these embodiments will be apparent to those skilled in the art and the general principles defined herein may be implemented in other embodiments without departing from the spirit or scope of the present invention. Therefore, the present invention will not be limited to these embodiments shown herein but shall accord with the widest scope consistent with the principles and novel features disclosed herein.