CN112630044B - Creep life prediction method of nickel-based single crystal alloy based on crystal orientation - Google Patents
Creep life prediction method of nickel-based single crystal alloy based on crystal orientation Download PDFInfo
- Publication number
- CN112630044B CN112630044B CN202011303059.XA CN202011303059A CN112630044B CN 112630044 B CN112630044 B CN 112630044B CN 202011303059 A CN202011303059 A CN 202011303059A CN 112630044 B CN112630044 B CN 112630044B
- Authority
- CN
- China
- Prior art keywords
- creep
- test piece
- stress
- nickel
- single crystal
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N3/08—Investigating strength properties of solid materials by application of mechanical stress by applying steady tensile or compressive forces
- G01N3/18—Performing tests at high or low temperatures
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/22—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by measuring secondary emission from the material
- G01N23/225—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by measuring secondary emission from the material using electron or ion
- G01N23/2251—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by measuring secondary emission from the material using electron or ion using incident electron beams, e.g. scanning electron microscopy [SEM]
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N3/08—Investigating strength properties of solid materials by application of mechanical stress by applying steady tensile or compressive forces
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N3/24—Investigating strength properties of solid materials by application of mechanical stress by applying steady shearing forces
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/0014—Type of force applied
- G01N2203/0025—Shearing
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/0058—Kind of property studied
- G01N2203/006—Crack, flaws, fracture or rupture
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/0058—Kind of property studied
- G01N2203/0069—Fatigue, creep, strain-stress relations or elastic constants
- G01N2203/0071—Creep
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/0058—Kind of property studied
- G01N2203/0069—Fatigue, creep, strain-stress relations or elastic constants
- G01N2203/0075—Strain-stress relations or elastic constants
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/02—Details not specific for a particular testing method
- G01N2203/022—Environment of the test
- G01N2203/0222—Temperature
- G01N2203/0226—High temperature; Heating means
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/10—Different kinds of radiation or particles
- G01N2223/102—Different kinds of radiation or particles beta or electrons
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/60—Specific applications or type of materials
- G01N2223/632—Specific applications or type of materials residual life, life expectancy
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
- Testing Resistance To Weather, Investigating Materials By Mechanical Methods (AREA)
Abstract
The disclosure relates to the technical field of nickel-based alloys, in particular to a creep life prediction method of a nickel-based single crystal alloy based on crystal orientation. The creep life prediction method comprises the following steps: performing a creep test on a plurality of test pieces of the nickel-based single crystal alloy with different crystal orientations to obtain creep curves of the test pieces; observing creep test processes of a plurality of test pieces to obtain a slip system start rule of the nickel-based single crystal alloy; based on a creep curve and a slip system start law, constructing a creep constitutive model and a creep damage model of the test piece; and constructing a life prediction model of the test piece based on the creep constitutive model and the creep damage model. The creep life prediction method can predict the creep life of the nickel-base single crystal alloy, and can reduce adverse effects caused by crystal orientation, thereby improving the creep performance of the nickel-base single crystal alloy.
Description
Technical Field
The disclosure relates to the technical field of nickel-based alloys, in particular to a creep life prediction method of a nickel-based single crystal alloy based on crystal orientation.
Background
The excellent high-temperature mechanical property of the nickel-based single crystal alloy mainly comes from a gamma/gamma 'two-phase microstructure, and the microstructure is formed by uniformly and coherent gamma' precipitate phases of a face-centered cubic structure with high volume fraction (about 65 percent) in a gamma matrix phase.
Since the high temperature properties of the nickel-based single crystal alloy have remarkable anisotropy, it is required that the crystal orientation with the minimum elastic modulus and the maximum load direction are consistent in the preparation of the single crystal blade to reduce the thermal cycle stress, but in practice, single crystals with the crystal orientation being severely deviated from the load direction or non-preferential growth orientation are often obtained, resulting in the creep life of the nickel-based single crystal alloy being affected.
Currently, there is no method for predicting the creep life of a nickel-base single crystal alloy based on the crystal orientation.
The above information disclosed in the background section is only for enhancement of understanding of the background of the disclosure and therefore it may include information that does not form the prior art that is already known to a person of ordinary skill in the art.
Disclosure of Invention
The present disclosure is directed to a method for predicting creep life of a nickel-base single crystal alloy based on crystal orientation, which can predict the creep life of the nickel-base single crystal alloy and reduce adverse effects due to crystal orientation, thereby improving creep performance of the nickel-base single crystal alloy.
In order to achieve the above purpose, the present disclosure adopts the following technical scheme:
according to an aspect of the present disclosure, there is provided a creep life prediction method of a nickel-based single crystal alloy based on crystal orientation, the creep life prediction method including:
performing a creep test on a plurality of test pieces of nickel-base single crystal alloys with different crystal orientations to obtain creep curves of the test pieces;
Observing creep test processes of a plurality of test pieces to obtain a slip system start rule of the nickel-based single crystal alloy;
Constructing a creep constitutive model and a creep damage model of the test piece based on the creep curve and the slip system start law;
and constructing a life prediction model of the test piece based on the creep constitutive model and the creep damage model.
In an exemplary embodiment of the present disclosure, constructing a creep constitutive model and a creep damage model of the test piece based on the creep curve and the slip train start law includes:
Based on the creep curve, acquiring creep curve parameters of the test piece;
Acquiring a sliding system parameter of the test piece based on the sliding system start rule;
And constructing the creep constitutive model and the creep damage model based on the creep curve parameters and the sliding system parameters.
In one exemplary embodiment of the present disclosure, the curve parameters include a temperature creep parameter, an initial damage rate, a critical split stress, a bergs vector mode, and a material constant; the slippage system parameters comprise slippage direction, slippage surface unit normal vector and material raft rate; the constructing the creep constitutive model and the creep damage model based on the creep curve parameters and the slip system parameters comprises the following steps:
Constructing the creep constitutive model based on the temperature creep parameter, the sliding direction and the sliding surface unit normal vector;
and combining the creep constitutive model, and constructing the creep damage model based on the initial damage rate, the critical slitting stress, the Bergers vector mode, the material constant and the material raft rate.
In one exemplary embodiment of the present disclosure, the creep constitutive model satisfies the following first relational expression:
In the method, in the process of the invention, A creep shear strain rate for the test piece; alpha represents different slip trains of the test piece; a is a first temperature creep parameter of the test piece; n is a second temperature creep parameter of the test piece; τ (α) is the shear stress of the slip system, and
Wherein sigma is the stress tensor of the test piece under the crystal axis system; p (α) is an orientation factor; m (α) is the slip direction of the slip system; n (α) is the unit normal vector of the slip plane.
In one exemplary embodiment of the present disclosure, the creep damage model satisfies the following second relation:
wherein ω (α) is the material damage of the test piece; The damage rate of the test piece is determined; /(I) For the initial injury rate;
τ (α) is the creep splitting stress of the slip system; τ or is the creep threshold stress of the test piece; Creep retarding stress for the test piece; τ c is the critical splitting stress;
beta is a constant; χ is a third temperature creep parameter of the test piece; phi is a fourth temperature creep parameter of the test piece;
An initial creep shear strain rate for the test piece; /(I) Is the steady state creep shear strain rate of the test piece.
In one exemplary embodiment of the present disclosure, the creep threshold stress τ or satisfies the following third relationship:
Wherein G is the shear modulus of the test piece; b is the Berger vector mode; lambda is the material constant; kappa is the current width of the matrix channel of the test piece, and
Wherein κ 0 is the initial matrix channel width of the test piece; c 1 is the material raft rate, and the value range of the material raft rate c 1 is 0.01-0.02 mm/s; t is creep time.
In an exemplary embodiment of the present disclosure, the barrier stressThe following fourth relational expression is satisfied:
where c 2 is a constant value, Is the dislocation density of a gamma matrix in the nickel-based single crystal alloy, and/>The evolution law is:
Wherein k 1 is a material constant characterizing dislocation stress hardening; k 2 is the material constant that characterizes dislocation stress recovery.
In an exemplary embodiment of the present disclosure, the critical slitting stress τ c satisfies the following fifth relationship:
τc＝Sf×σ0.2
Wherein S f is a Schmitt factor; σ 0.2 is the yield stress of the nickel-base single crystal alloy.
In an exemplary embodiment of the present disclosure, the initial creep shear strain rateThe following sixth relation is satisfied:
Wherein R is a gas constant; t is absolute temperature; q is the activation energy.
In one exemplary embodiment of the present disclosure, the life prediction model satisfies the following seventh relational expression:
Wherein t f is the breaking time of the test piece; n is the number of the slip plane starting directions under the slip system, N corresponding to the hexahedral slip system is 6, and N corresponding to the octahedral slip system and the dodecahedral slip system is 12.
In the creep life prediction method of the nickel-base single crystal alloy based on the crystal orientation of the disclosed embodiment, in the operation process, first, a plurality of test pieces of the nickel-base single crystal alloy with different crystal orientations are subjected to creep tests to obtain creep curves of the test pieces; secondly, observing creep test processes of a plurality of test pieces to obtain a slip system start rule sum of the nickel-based single crystal alloy; then, based on the creep curve, the starting law of the sliding system and the sum, constructing a creep constitutive model and a creep damage model of the test piece; and finally, constructing a life prediction model of the test piece based on the creep constitutive model and the creep damage model.
Therefore, on one hand, the service life of the nickel-based single crystal alloy can be calculated according to the life prediction model, so that reference is provided for the actual use of the nickel-based single crystal alloy; on the other hand, the crystal orientation can be optimized based on the service life of the nickel-based single crystal alloy, so that adverse effects caused by the crystal orientation are effectively reduced, and the creep property of the nickel-based single crystal alloy is further improved.
Drawings
The accompanying drawings, which are incorporated in and constitute a part of this specification, illustrate embodiments consistent with the disclosure and together with the description, serve to explain the principles of the disclosure. It will be apparent to those of ordinary skill in the art that the drawings in the following description are merely examples of the disclosure and that other drawings may be derived from them without undue effort.
FIG. 1 is a flow chart of a method for predicting creep life of a crystal orientation-based nickel-base single crystal alloy according to an embodiment of the present disclosure.
Fig. 2 is a schematic diagram of a crystallographic orientation of an embodiment of the present disclosure.
Fig. 3 is a schematic diagram of dislocation bowing out of the matrix phase past the strengthening phase in an embodiment of the present disclosure.
Detailed Description
Example embodiments will now be described more fully with reference to the accompanying drawings. However, the exemplary embodiments may be embodied in many different forms and should not be construed as limited to the examples set forth herein; rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the concept of the example embodiments to those skilled in the art. The described features, structures, or characteristics may be combined in any suitable manner in one or more embodiments. In the following description, numerous specific details are provided to give a thorough understanding of embodiments of the disclosure.
The described features, structures, or characteristics may be combined in any suitable manner in one or more embodiments. In the following description, numerous specific details are provided to give a thorough understanding of embodiments of the disclosure. One skilled in the relevant art will recognize, however, that the disclosed aspects may be practiced without one or more of the specific details, or with other methods, components, materials, etc. In other instances, well-known structures, materials, or operations are not shown or described in detail to avoid obscuring the main technical ideas of the present disclosure.
Although relative terms such as "upper" and "lower" are used in this specification to describe the relative relationship of one component of an icon to another component, these terms are used in this specification for convenience only, such as in terms of the orientation of the examples described in the figures. It will be appreciated that if the device of the icon is flipped upside down, the "up" component will become the "down" component. Other relative terms such as "high," "low," "top," "bottom," "left," "right," and the like are also intended to have similar meanings.
When a structure is "on" another structure, it may mean that the structure is integrally formed with the other structure, or that the structure is "directly" disposed on the other structure, or that the structure is "indirectly" disposed on the other structure through another structure. The terms "a," "an," "the" are used to indicate the presence of one or more elements/components/etc.; the terms "comprising" and "having" are intended to be inclusive and mean that there may be additional elements/components/etc. in addition to the listed elements/components/etc. The terms "first" and "second" and the like are used merely as labels, and are not intended to limit the number of their objects.
Embodiments of the present disclosure provide a method for predicting creep life of a nickel-base single crystal alloy based on crystal orientation, which can predict the creep life of the nickel-base single crystal alloy. As shown in fig. 1, the creep life prediction method of a crystal orientation-based nickel-base single crystal alloy may include the steps of:
Step S110, performing a creep test on a plurality of test pieces of nickel-based single crystal alloys with different crystal orientations to obtain creep curves of the test pieces;
step S120, observing creep test processes of a plurality of test pieces to obtain a slip system start rule of the nickel-based single crystal alloy;
Step S130, a creep constitutive model and a creep damage model of the test piece are constructed based on a creep curve and a slip system start law;
And step S140, constructing a life prediction model of the test piece based on the creep constitutive model and the creep damage model.
According to the creep life prediction method of the nickel-based single crystal alloy based on the crystal orientation, on one hand, the service life of the nickel-based single crystal alloy can be calculated according to the life prediction model, so that reference is provided for practical use of the nickel-based single crystal alloy; on the other hand, the crystal orientation can be optimized based on the service life of the nickel-based single crystal alloy, so that adverse effects caused by the crystal orientation are effectively reduced, and the creep property of the nickel-based single crystal alloy is further improved.
The following describes in detail a creep life prediction method of a crystal orientation-based nickel-base single crystal alloy provided in an embodiment of the present disclosure:
And step S110, performing creep test on a plurality of test pieces of nickel-based single crystal alloys with different crystal orientations to obtain creep curves of the test pieces.
It is easy to understand that the following operations are also performed before step S110:
First, a plurality of nickel-based single crystal alloys having different crystal orientations are cast, for example, the nickel-based single crystal alloys may be three kinds, as shown in fig. 2, the respective crystal orientations may be respectively denoted as [001], [010] and [100], and the [001], [010] and [100] are perpendicular to each other; secondly, processing three nickel-based single crystal alloys to obtain three test pieces with different crystal orientations, wherein the crystal orientations of the three test pieces are respectively [001], [010] and [100], and the three test pieces are not described in detail herein; finally, three test pieces are tested and original microstructure information of each test piece is obtained, and specific testing processes are not described in detail herein.
In step S110, a plurality of test pieces of nickel-based single crystal alloys having different crystal orientations may be divided into three groups, and the three groups of test pieces may be subjected to creep test at a high temperature (1000 ℃ or higher), a medium temperature (800 ℃ to 1000 ℃) and a low temperature (800 ℃ or lower), respectively. As described above, the number of test pieces may be three, and accordingly, the three test pieces are subjected to creep tests at high temperature, medium temperature and low temperature, respectively, to thereby obtain creep curves of the three test pieces.
And S120, observing creep test processes of a plurality of test pieces to obtain a slip system start rule of the nickel-based single crystal alloy.
Specifically, a Scanning Electron Microscope (SEM) is adopted to observe the microstructure evolution process in the creep process, and information such as coarsening, microcrack formation, microcrack expansion, fracture modes and the like of the strengthening phase is obtained, so that the nickel-based single crystal alloy is obtained. Meanwhile, a Transmission Electron Microscope (TEM) is adopted to observe dislocation morphology of a test piece in three orientations of [001], [010] and [100] in the early creep stage, microstructure evolution information such as the type, generation and increment of surface layer dislocation, crystal slip characteristics and the like is obtained, and thus a slip system start rule of the nickel-based single crystal alloy is obtained.
And step S130, constructing a creep constitutive model and a creep damage model of the test piece based on the creep curve and the slip system start law. To describe in detail, step S130 may include the steps of:
Step S1301, based on the creep curve, obtaining the creep curve parameters of the test piece. Specifically, the creep curve parameters may include temperature creep parameters, initial damage rate, critical splitting stress, bergs vector mode, material constants, etc., which are not described herein;
Step S1302, acquiring the sliding system parameters of the test piece based on the sliding system start law. Specifically, the sliding system parameters include a sliding direction, a sliding surface unit normal vector, a material raft rate and the like, which are not described one by one;
step S1303, a creep constitutive model and a creep damage model are constructed based on the creep curve parameters and the slip system parameters.
Specifically, step S1303 includes the following two steps:
Step S13031, constructing a creep constitutive model based on temperature creep parameters, a slip direction and a slip plane unit normal vector; step S13032, combining the creep constitutive model, and constructing the creep damage model based on the initial damage rate, the critical cutting stress, the Berger vector mode, the material constant and the material raft rate.
The creep constitutive model can meet the following first relation:
In the method, in the process of the invention, The creep shear strain rate of the test piece; alpha represents different sliding systems of the test piece, and the sliding systems of the nickel-based single crystal alloy comprise hexahedral sliding systems, octahedral sliding systems and dodecahedral sliding systems; a is a first temperature creep parameter of the test piece, n is a second temperature creep parameter of the test piece, and both the first temperature creep parameter and the second temperature creep parameter can be obtained by a creep curve of the test piece and are changed along with the change of a sliding system; τ (α) is the shear stress of the slip system, and
Wherein sigma is the stress tensor of the test piece under the crystal axis system; p (α) is an orientation factor; m (α) is the slip direction of the slip system, n (α) is the unit normal vector of the slip plane, both of which can be derived from the slip system start law and the slip system start law of the nickel-based single crystal alloy, and will not be described in detail herein.
It should be noted that the constitutive model of the test piece can also be characterized by a creep deformation strain rate, and the creep deformation strain rate can be noted asAnd satisfies the following relation:
In the method, in the process of the invention, Is macroscopic strain rate, and/>And/>Is the strain rate of the elastic part,/>Is the strain rate of the inelastic portion, and
Then decomposing creep strain:
Wherein, Corresponding to the creep strain of hexahedral sliding system,/>Corresponding to octahedral glide system/>Corresponding to the dodecahedron sliding system. If one of the slips is not activated, the corresponding value is zero.
In addition, C e is an anisotropic elastic tensor, an
Wherein C 11、C12 and C 44 are both elastic constants, anE is the elastic modulus, μ is the Poisson's ratio, and G is the shear modulus.
In order to obtain the elastic modulus E, the Poisson's ratio mu and the shear modulus G, a tensile test is required to be carried out on a test piece, a tensile curve of the test piece is obtained, and then the obtained tensile curve is fitted, so that the elastic modulus E, the Poisson's ratio mu and the shear modulus G are solved.
Since the nickel-based single crystal alloy is an anisotropic material, the anisotropic elastic tensor C e is only suitable for the [001] crystal axis, so when the used coordinate system is different from the [001] crystal axis, the C e matrix needs to be subjected to coordinate transformation, and can be obtained by matrix operation:
CXYZ＝[T][C][T]T
Wherein,
Where l, m and n are the directional cosine of the model coordinate axes O-X-Y-Z in the crystal axes O-X-Y-Z, the crystal axes are the three crystal orientations [001], [010] and [100], and will not be described in detail here.
Meanwhile, the creep damage model may satisfy the following second relational expression:
the following describes each parameter in the second relation in detail:
①ω(α) Material damage for the test piece; The damage rate of the test piece; /(I) Is the initial damage rate; ②τ(α) The creep rupture stress of the sliding system is described in detail above and will not be described here again; ③τor In order to measure the creep threshold stress of the test piece, the bypass mechanism of the dislocation can prevent the dislocation from moving in the matrix, and the width of the matrix channel is widened along with the increase of time in the creep process, so that the bypass mechanism of the dislocation under the external force is controlled by the width of the matrix channel. As shown in fig. 3, if the stress required for dislocations to bypass the strengthening phase is less than the stress of the shear strengthening phase, the material plastic deformation results primarily from dislocation motion in the matrix phase.
Thus, the creep threshold stress τ or satisfies the following third relationship:
Wherein G is the shear modulus of the test piece; b is a Bergers vector mode; lambda is the material constant; kappa is the current width of the matrix channel of the test piece, and
Wherein, kappa 0 is the initial matrix channel width of the test piece; c 1 is the material raft rate, and the value range of c 1 is 0.01-0.02 mm/s; t is creep time.
④Stress is hindered for creep of the test piece, and stress/>, is hinderedThe following fourth relational expression is satisfied:
where c 2 is a constant value, Is the dislocation density of a gamma matrix in the nickel-based single crystal alloy, and/>The evolution law is:
Wherein k 1 is a material constant characterizing dislocation stress hardening; k 2 is the material constant that characterizes dislocation stress recovery.
⑤τc For critical splitting stress, based on the theory of crystal plasticity, the relationship between splitting stress τ and macroscopic stress σ resolved onto each slip plane can be expressed as follows:
τ＝Sfσ
wherein S f is Schmid (Schmidt) factor, and the specific values are shown in Table 1.
TABLE 1S of three oriented slip systems f
As described above, the tensile curve of the test piece can be obtained by performing the tensile test on the test piece, and accordingly, the yield stress σ 0.2 can be obtained according to the tensile curve, and the critical splitting stress τ c＝Sf×σ0.2 can be obtained by combining S f, which will not be described in detail herein.
⑥The steady-state creep shear strain rate of the test piece; /(I)Is the initial creep shear strain rate of the test piece, and the initial creep shear strain rate/>The following sixth relation may be satisfied:
Wherein R is a gas constant; t is absolute temperature; q is the activation energy, the Q value of the hexahedral slip system is 7.3X10 -19 J/atom, the Q value of the octahedral slip system is 6.97X10 -19 J/atom, and the Q value of the dodecahedral slip system is required to be determined through experiments, and will not be described in detail herein.
⑦ Chi is a third temperature creep parameter of the test piece, phi is a fourth temperature creep parameter of the test piece, and both can be obtained by a creep curve of the test piece and changed along with the change of a sliding system; beta is a constant and will not be described in detail here.
Step S140, a life prediction model of the test piece is constructed based on the creep constitutive model and the creep damage model, and the life prediction model can meet the following seventh relational expression:
Wherein t f is the breaking time of the test piece; τ c is the critical splitting stress; n is the number of slip plane starting directions under the slip system, specifically, N corresponding to the hexahedral slip system is 6, and N corresponding to the octahedral slip system and the dodecahedral slip system is 12.
It should be noted that before predicting the breaking time (lifetime) of the test piece, the failure mode of the test piece needs to be determined, specifically, firstly, R (α) is used as a creep (long-term) failure criterion and a tensile (short-term) failure criterion, and
Second, whenWhen R (α) is less than or equal to 0, the matrix of the test piece is free from sliding deformation, and the test piece is not damaged at the moment;
When (when) I.e./>When the test piece is in a creep deformation state, the corresponding creep damage model meets the second relation, and the corresponding life prediction model meets the sixth relation, which are not repeated here;
When (when) I.e., τ (α)≥τc, the test piece will experience an instantaneous tensile failure with an external load greater than the material yield strength.
It is to be understood that the disclosure is not limited in its application to the details of construction and the arrangement of components set forth in the disclosure. The disclosure is capable of other embodiments and of being practiced and carried out in various ways. The foregoing variations and modifications are within the scope of the present disclosure. It should be understood that the present disclosure disclosed and defined herein extends to all alternative combinations of two or more of the individual features mentioned or evident from the text and/or drawings. All of these different combinations constitute various alternative aspects of the present disclosure. The embodiments described herein explain the best modes known for practicing the disclosure and will enable others skilled in the art to utilize the disclosure.
Claims (9)
1. A method for predicting creep life of a nickel-base single crystal alloy based on crystal orientation, the method comprising:
performing a creep test on a plurality of test pieces of nickel-base single crystal alloys with different crystal orientations to obtain creep curves of the test pieces;
Observing creep test processes of a plurality of test pieces to obtain a slip system start rule of the nickel-based single crystal alloy;
Constructing a creep constitutive model and a creep damage model of the test piece based on the creep curve and the slip system start law;
constructing a life prediction model of the test piece based on the creep constitutive model and the creep damage model;
Wherein the life prediction model satisfies the following seventh relational expression:
Wherein t f is the breaking time of the test piece; n is the number of the slip plane starting directions under the slip system, N corresponding to the hexahedral slip system is 6, N corresponding to the octahedral slip system and the dodecahedral slip system is 12, Is the initial damage rate; τ c is the critical splitting stress; beta is a constant; phi is a fourth temperature creep parameter of the real blade sample; χ is a third temperature creep parameter of the test piece; r (α) is creep failure criterion and tensile failure criterion; alpha represents different slip trains of the test piece; n is a second temperature creep parameter of the test piece.
2. The creep life prediction method according to claim 1, wherein constructing a creep constitutive model and a creep damage model of the test piece based on the creep curve and the slip system start law comprises:
Based on the creep curve, acquiring creep curve parameters of the test piece;
Acquiring a sliding system parameter of the test piece based on the sliding system start rule;
And constructing the creep constitutive model and the creep damage model based on the creep curve parameters and the sliding system parameters.
3. The creep life prediction method according to claim 2, wherein the curve parameters include a temperature creep parameter, an initial damage rate, a critical splitting stress, a berges vector mode, and a material constant; the slippage system parameters comprise slippage direction, slippage surface unit normal vector and material raft rate;
Based on the creep curve parameters and the slip system parameters, constructing the creep constitutive model and the creep damage model comprises the following steps:
Constructing the creep constitutive model based on the temperature creep parameter, the sliding direction and the sliding surface unit normal vector;
and combining the creep constitutive model, and constructing the creep damage model based on the initial damage rate, the critical slitting stress, the Bergers vector mode, the material constant and the material raft rate.
4. The creep life prediction method according to claim 3, wherein the creep constitutive model satisfies the following first relational expression:
In the method, in the process of the invention, A creep shear strain rate for the test piece; alpha represents different slip trains of the test piece; a is a first temperature creep parameter of the test piece; n is a second temperature creep parameter of the test piece; τ (α) is the shear stress of the slip system, and
Wherein sigma is the stress tensor of the test piece under the crystal axis system; p (α) is an orientation factor; m (α) is the slip direction of the slip system; n (α) is the unit normal vector of the slip plane.
5. The creep life prediction method according to claim 4, wherein the creep damage model satisfies the following second relational expression:
wherein ω (α) is the material damage of the test piece; The damage rate of the test piece is determined; /(I) For the initial injury rate;
τ (α) is the creep splitting stress of the slip system; τ or is the creep threshold stress of the test piece; Creep retarding stress for the test piece; τ c is the critical splitting stress;
beta is a constant; χ is a third temperature creep parameter of the test piece; phi is a fourth temperature creep parameter of the test piece;
An initial creep shear strain rate for the test piece; /(I) Is the steady state creep shear strain rate of the test piece.
6. The method of claim 5, wherein the creep threshold stress τ or satisfies the following third relation:
Wherein G is the shear modulus of the test piece; b is the Berger vector mode; lambda is the material constant; kappa is the current width of the matrix channel of the test piece, and
Wherein κ 0 is the initial matrix channel width of the test piece; c 1 is the material raft rate, and the value range of the material raft rate c 1 is 0.01-0.02 mm/s; t is creep time.
7. The creep life prediction method according to claim 6, wherein the barrier stressThe following fourth relational expression is satisfied:
wherein c 2 is a constant; Is the dislocation density of a gamma matrix in the nickel-based single crystal alloy, and/> The evolution law is:
Wherein k 1 is a material constant characterizing dislocation stress hardening; k 2 is the material constant that characterizes dislocation stress recovery.
8. The creep life prediction method according to claim 7, wherein the critical split stress τ c satisfies the following fifth relation:
τc＝Sf×σ0.2
Wherein S f is a Schmitt factor; σ 0.2 is the yield stress of the nickel-base single crystal alloy.
9. The method of creep life prediction according to claim 8, wherein the initial creep shear strain rateThe following sixth relation is satisfied:
Wherein R is a gas constant; t is absolute temperature; q is the activation energy.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011303059.XA CN112630044B (en) | 2020-11-19 | 2020-11-19 | Creep life prediction method of nickel-based single crystal alloy based on crystal orientation |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011303059.XA CN112630044B (en) | 2020-11-19 | 2020-11-19 | Creep life prediction method of nickel-based single crystal alloy based on crystal orientation |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112630044A CN112630044A (en) | 2021-04-09 |
| CN112630044B true CN112630044B (en) | 2024-06-18 |
Family
ID=75303448
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011303059.XA Active CN112630044B (en) | 2020-11-19 | 2020-11-19 | Creep life prediction method of nickel-based single crystal alloy based on crystal orientation |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112630044B (en) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113420502B (en) * | 2021-06-11 | 2023-08-11 | 北京科技大学 | Re-Ru-containing nickel-based single crystal superalloy creep life prediction method |
| CN114091315B (en) * | 2022-01-24 | 2022-04-15 | 西北工业大学 | High-temperature alloy stress and damage evolution method and device, storage medium and electronic equipment |
| CN114818481A (en) * | 2022-04-13 | 2022-07-29 | 北航（四川）西部国际创新港科技有限公司 | Creep rupture life prediction method of single crystal high temperature alloy |
| CN114813329B (en) * | 2022-04-26 | 2023-07-28 | 北京航空航天大学杭州创新研究院 | A Remaining Life Prediction Method for Nickel-based Single Crystal Superalloys Under Random Creep Loads |
Family Cites Families (24)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5042295A (en) * | 1985-06-21 | 1991-08-27 | General Electric Company | Method for determining remaining useful life of turbine components |
| US4947341A (en) * | 1988-03-15 | 1990-08-07 | Digital Equipment Corporation | Method of predicting fatigue lifetimes of metallic structures |
| JP4186691B2 (en) * | 2003-04-28 | 2008-11-26 | マツダ株式会社 | Engine part failure prediction system and control program therefor |
| JP4618061B2 (en) * | 2005-09-06 | 2011-01-26 | 株式会社Ｉｈｉ | Creep rupture life evaluation method |
| JP5090142B2 (en) * | 2007-11-28 | 2012-12-05 | ポリプラスチックス株式会社 | Prediction method for brittle creep rupture remaining life of molded parts |
| CN102108555B (en) * | 2009-12-23 | 2012-08-29 | 中国科学院金属研究所 | High-temperature perfect oxidation resistance type nickel-based single-crystal alloy and preparation method thereof |
| KR20120105693A (en) * | 2011-03-16 | 2012-09-26 | 한국기계연구원 | Ni base single crystal superalloy with enhanced creep property |
| CN102162049B (en) * | 2011-04-07 | 2012-12-19 | 上海大学 | Nickel-base alloy material for super-supercritical steam turbine and preparation method thereof |
| JP5943829B2 (en) * | 2012-12-25 | 2016-07-05 | 三菱重工業株式会社 | Life prediction system and life prediction method for reinforced high temperature aging pipes |
| CN103344511B (en) * | 2013-07-01 | 2015-09-30 | 北京理工大学 | A kind of Creep-fatigue Life Prediction method of high-power diesel engine all-aluminium piston |
| CN105158084B (en) * | 2015-09-15 | 2018-02-09 | 华东理工大学 | A kind of Creep-Fatigue Life Prediction method of material |
| CN105784508B (en) * | 2016-04-11 | 2018-08-10 | 沈阳工业大学 | A method of characterization monocrystalline Ni based alloy croop properties |
| CN106568655B (en) * | 2016-10-28 | 2019-04-12 | 沈阳工业大学 | A method of prediction heat-resisting alloy creep life |
| CN106557630A (en) * | 2016-11-21 | 2017-04-05 | 中国石油大学(华东) | A kind of creep impairment life-span prediction method of material under multi-axis stress state |
| GB2565063B (en) * | 2017-07-28 | 2020-05-27 | Oxmet Tech Limited | A nickel-based alloy |
| CN107273649B (en) * | 2017-08-16 | 2018-05-04 | 中国石油大学(华东) | A kind of Forecasting Methodology of fragile material failure probability under high-temerature creep state |
| CN109142083B (en) * | 2018-09-05 | 2020-04-24 | 南京航空航天大学 | Creep damage calculation method under variable load process |
| CN110175394A (en) * | 2019-05-23 | 2019-08-27 | 电子科技大学 | A kind of turbo blade spleen tissue extracts damage coupling Probabilistic Life Prediction calculation method |
| CN110688788B (en) * | 2019-08-28 | 2021-06-22 | 南京航空航天大学 | A Creep Deformation and Life Prediction Method for High Temperature Materials |
| CN111027198B (en) * | 2019-12-03 | 2022-07-19 | 西北工业大学 | A creep-life prediction method for nickel-based single crystal alloys considering topologically close-packed phase evolution |
| CN111008495B (en) * | 2019-12-03 | 2022-12-06 | 西北工业大学 | Method for predicting creep residual life of nickel-based single crystal turbine blade |
| CN110940572A (en) * | 2019-12-11 | 2020-03-31 | 北京科技大学 | A creep life prediction method of high Cr ferritic heat-resistant steel |
| CN111534720A (en) * | 2020-05-12 | 2020-08-14 | 山东大学 | Twin crystal strengthened nickel-based high-temperature alloy and preparation method and application thereof |
| CN111500896B (en) * | 2020-05-19 | 2021-09-17 | 西北工业大学 | Gamma' phase reinforced third generation nickel base single crystal high temperature alloy and preparation method thereof |
-
2020
- 2020-11-19 CN CN202011303059.XA patent/CN112630044B/en active Active
Non-Patent Citations (1)
| Title |
|---|
| 镍基单晶合金超温超载蠕变性能研究;梁建伟;《中国博士学位论文全文数据库 工程科技II辑》;摘要，第2章第12-27页，附图2.6-2.9，公式2-11-2-12 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112630044A (en) | 2021-04-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112630044B (en) | Creep life prediction method of nickel-based single crystal alloy based on crystal orientation | |
| CN112630045B (en) | Creep life prediction method of nickel-based single crystal alloy based on real blade sample | |
| Tanaka et al. | Creep deformation of single crystals of new Co–Al–W-based alloys with fcc/L12 two-phase microstructures | |
| Korte-Kerzel | Microcompression of brittle and anisotropic crystals: recent advances and current challenges in studying plasticity in hard materials | |
| Priest et al. | A combined deformation map-ductility exhaustion approach to creep-fatigue analysis | |
| Segersäll et al. | Influence of crystal orientation on the thermomechanical fatigue behaviour in a single-crystal superalloy | |
| Sidey et al. | Low-cycle fatigue damage mechanisms at high temperature | |
| Gao et al. | Microstructural evolution of martensite during deformation in Zr50Cu50 shape memory alloy | |
| Shi et al. | Tensile behavior of the second generation single crystal superalloy DD6 | |
| Dong et al. | Life modeling of anisotropic fatigue behavior for a single crystal nickel-base superalloy | |
| CN105302987B (en) | A kind of method of equivalent prediction Thermomechanical Fatigue Life | |
| Ma et al. | Temperature effect on low-cycle fatigue behavior of nickel-based single crystalline superalloy | |
| Zhang et al. | Creep deformation mechanisms in some modern single-crystal superalloys | |
| Wang et al. | Crystal plasticity theory coupled with meso-damage to predict the ratchetting behavior of nickel-based single crystal superalloy | |
| Bauer et al. | Thermomechanical fatigue behaviour of a third generation γ-TiAl intermetallic alloy | |
| Tian et al. | Microstructural adjusting crack evolution of polycrystalline NCM particle during charge/discharge cycle | |
| Liang et al. | Analysis and prediction of non-isothermal creep behavior in Ni-based single crystal superalloy | |
| Dong et al. | Fatigue life modeling of a single crystal superalloy and its thin plate with a hole at elevated temperature | |
| CN110232223B (en) | Prediction method for microstructure rafting type of nickel-based single crystal superalloys | |
| Deguchi et al. | Damage propagation mechanism in low-cycle creep fatigue of Cu–Cr–Zr alloy | |
| CN109117604A (en) | A kind of acquisition methods of Walker cast material constant of the Crystal Nickel-based Superalloy based on slip system | |
| Zhou et al. | Investigations on the thermo-mechanical fatigue of two Ni-based single-crystal superalloys | |
| Liu et al. | Evolution of defects and deformation mechanisms in different tensile directions of solidified lamellar Ti–Al alloy | |
| Lv et al. | Thickness-dependent creep rate caused by climbing control of a< 100> superdislocation in a Ni-based single crystal superalloy | |
| CN113109190B (en) | A short-crack-based life prediction method under multiaxial thermomechanical loading |
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 |