Next Article in Journal
Asymmetric Contrarians in Opinion Dynamics
Next Article in Special Issue
Effective Surface Nano-Crystallization of Ni2FeCoMo0.5V0.2 Medium Entropy Alloy by Rotationally Accelerated Shot Peening (RASP)
Previous Article in Journal
Many-Body Dynamics and Decoherence of the XXZ Central Spin Model in External Magnetic Field
Previous Article in Special Issue
Prediction of Fatigue Crack Growth Rate Based on Entropy Generation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Low Cycle Fatigue Life Prediction Using Unified Mechanics Theory in Ti-6Al-4V Alloys

1
Department of Applied Mechanics, Indian Institute of Technology, Madras 600036, India
2
Civil, Structural and Environmental Engineering, University at Buffalo, State University of New York, New York, NY 10031, USA
*
Author to whom correspondence should be addressed.
Entropy 2020, 22(1), 24; https://doi.org/10.3390/e22010024
Submission received: 22 November 2019 / Revised: 19 December 2019 / Accepted: 22 December 2019 / Published: 23 December 2019

Abstract

:
Fatigue in any material is a result of continuous irreversible degradation process. Traditionally, fatigue life is predicted by extrapolating experimentally curve fitted empirical models. In the current study, unified mechanics theory is used to predict life of Ti-6Al-4V under monotonic tensile, compressive and cyclic load conditions. The unified mechanics theory is used to derive a constitutive model for fatigue life prediction using a three-dimensional computational model. The proposed analytical and computational models have been used to predict the low cycle fatigue life of Ti-6Al-4V alloys. It is shown that the unified mechanics theory can be used to predict fatigue life of Ti-6Al-4V alloys by using simple predictive models that are based on fundamental equation of the material, which is based on thermodynamics associated with degradation of materials.

1. Introduction

Titanium alloys are popular for their superior mechanical properties, such as high yield strength, long fatigue life, toughness, low density, as well as corrosion resistance. About 80% of the global production of titanium alloys are used by aerospace industries [1]. One of the widely used titanium alloys is Ti-6Al-4V [2] which has a dual-phase crystal structure, namely, hexagonal close packed (HCP) and body centered cubic (BCC) structures. In the composition of Ti-6Al-4V alloy, titanium is the matrix material. Aluminium plays the role of stabilizing the HCP structure and vanadium preserve the BCC structure [3]. Many applications of Ti-6Al-4V alloys, such as aero engines, are subjected to cyclic loading [4]. Hence, it is essential to predict the fatigue life of such structural components, when they are subjected to varying amplitudes of cyclic loading during their service period. It is not always feasible to conduct fatigue experiments corresponding to all service conditions. Hence, predictive models based on fundamental physics of materials are helpful in predicting the fatigue life of structures.
A number of studies have been published to investigate the fatigue life of metals. Most of the damage prediction models are based on statistical test data analysis or on experimental curve fit [5,6,7,8,9,10,11]. Low cycle fatigue life prediction in Ti-6Al-4V alloys are generally done, based on stress [12], strain [5,6,13,14,15,16] or hysteresis loss [17]. Most of them are empirical curve-fit models [7,9,13,18,19,20,21,22] or mechanism based phenomenological models [23,24,25] such as fatigue crack initiation models [16]. A detailed review of such models, applied to metals, can be seen in the review article by Santecchia et al. [26]. A model, based on combined Newtonian mechanics and thermodynamics, instead of material-specific and loading-specific, can capture the mechanisms of fatigue damage without the need for curve fitting process.
If the system is less complicated and we want a quick solution we can opt for a one-dimensional model based on certain assumptions. However, validity of the model depends upon the accuracy of the assumptions made while formulation of one-dimensional analytical model. The interpretation of the results using one-dimensional model is also easy as it can be simple in its form and usage. A number of one-dimensional empirical curve-fit fatigue life prediction models can be seen in the literature [5,6,7,8,11,12,14,15,16,17]. Nevertheless, a physics-based one-dimensional model, which can be easily used to predict the fatigue life of Ti-6Al-4V, under appropriate assumptions, is still not found in the literature. If the system is very complicated to arrive at suitable one-dimensional fatigue life prediction model, we look for another appropriate and convenient method. It is known that, a three dimensional computational model can be incorporated with appropriate material nonlinearities (such as plastic flow), to account for the experimental observations [10,22] and to limit the assumptions in developing the model. However, a large number of cyclic loading simulation in a three dimensional numerical model is computationally very expensive [10]. Hence, it is very useful to have an appropriate physics-based procedure, in conjunction with three-dimensional numerical results, to account for all the nonlinearities associated with the computational model, even as we maintain the simplistic predictive capability of a one-dimensional model. Therefore, the present study is focused on both one-dimensional and three-dimensional, thermodynamics-based modeling of the deformation of standard test specimen to predict the fatigue life of Ti-6Al-4V.
Thermodynamics is a field of science that is developed to study change in the state of matter. The historical development of thermodynamics from its classical form to modern-age form has been reviewed by Haddad et al. [27,28]. Between 1872 and 1875, Boltzmann gave a mathematical expression to second law of thermodynamics for quantification of order/disorder in terms of a measure called entropy. In 1998, Basaran and Yan [29] introduced the unified mechanics theory, which unifies Newtonian mechanics with thermodynamics. In unified mechanics theory [29], in addition to nodal displacements, the entropy generation rate is also necessary to relate microstructural changes in the material with spatial and temporal coordinates. This concept [29] has been successfully implemented for a wide range of materials and has been experimentally and mathematically validated and reported in literature [18,19,20,25,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60,61,62,63,64,65]. The entropy generation rate of any material under any external disturbances like mechanical, thermal, electrical, chemical, radiation, and corrosion can be calculated from principles of physics, using the fundamental equation, with no need for curve fitting phenomenological models or polynomials fit to experimental test data.
In the present study, unified mechanics theory is used to estimate the fatigue damage in Ti-6Al-4V, analytically with a one -dimensional (1-D) model as well as numerically with a three-dimensional (3-D) model, and this damage estimation procedure has been used to predict fatigue life under different loading conditions. Fundamental details of the unified mechanics theory-based fatigue life prediction are summarized in Section 2. The principles described in Section 2, are then applied to Ti-6Al-4V, by considering the plasticity as the dominant energy dissipation mechanism.
In order to establish the validity of the proposed model in cyclic loading, comparison of simulation with experimental results, under both the tensile and compressive loading are necessary. In Section 3, the details of implementation and validation of the proposed model, for both compressive and tensile monotonic loading is presented. After the validation of the proposed model, we introduce two different procedures, to estimate the low cycle fatigue life of Ti-6Al-4V alloys in Section 4. Finally, the observations from the presented work are discussed in Section 5, based on the observations made on the principles, procedure and results from the current study for the fatigue life prediction of Ti-6Al-4V alloys.

2. Unified Mechanics Theory-Based Life Prediction Model

2.1. Unified Mechanics Theory

Unified mechanics theory is just unification of Newton’s universal laws of motion and laws of thermodynamics.

2.1.1. Second Law of Unified Mechanics Theory

Initial momentum of a mass, m, subjected to external force, F is defined by Newton’s second universal law of motion. However, Newton’s laws do not account for energy loss after the initial momentum. Energy loss takes place according to the first and second laws of thermodynamics. As a result, a marriage of laws of second law of Newton and laws of thermodynamic is given by:
F = d P d t = d ( m v ) d t ( 1 Φ )
where, P represents the momentum and v represents the velocity. Assuming a constant mass system,
F = m d [ v ( 1 Φ ) ] d t
where, Φ is the Thermodynamic State Index (TSI), which is normalized non-dimensional form of the second law of thermodynamics. TSI ( Φ ) starts at zero and reaches one when the system reaches maximum entropy and minimum entropy generation rate. The value of TSI ( Φ ) is calculated from the fundamental equation of the material, which accounts for all entropy generation mechanisms in the system under the given load towards a pre-defined failure. The fundamental equation must satisfy the conservation of energy, the first law of thermodynamics at every step. Therefore, TSI ( Φ ) just introduces laws of thermodynamics in to the laws of Newton.

2.1.2. Third Law of Unified Mechanics Theory

All forces between two objects exist in equal magnitude and opposite direction (Action–Reaction). However, resulting deformation, according to Hook’s law, in two objects will change over time because of degradation. The resulting equation can be given by:
F 12 = F 21 [ 1 Φ ]
where, the subscripts 12 and 21 represents the action and reaction, respectively. Based on Hooke’s law, the reaction, F 21 can be given by the following:
F 12 = d U 21 d u 21 = [ d [ 1 2   k 21 [ 1 Φ ]   u 21 2 ] ] d u 21
where, U 21 is the strain energy of the reactionary member, k 21 is the stiffness of the reactionary member, u 21 is the displacement in the reactionary member. If we assume that for the increment of displacement, d u 21 derivative of TSI with respect to d u 21 is smaller than derivative of displacement u 12 by an order of magnitude as the differential in displacement d u 21 goes to zero in the limiting case, we can write the following simple relation:
F 12 = k 21 [ 1 Φ ]   u 21
In unified mechanics theory, it has been shown that the degradation of the stiffness follows the laws of thermodynamics [8,18,20,22,27,29,30,31,32,33,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,56,57,58,59,66,67,68,69]. Combining laws of Newton and thermodynamics requires the modification of Newtonian space-time coordinate system. A new thermodynamic axis must be added to be able to define the thermodynamic state of a point. As a result, the motion of any particle can be defined only in a five-dimensional space that has five linearly independent axes. None of these axes can represent the information of other axes. Hence, entropy generation can be mapped onto a non-dimensional coordinate called Thermodynamics State Index (TSI) which is necessary to locate the thermodynamic state of the particle. Coordinates of a point can be defined by Newton’s laws of motion in the space-time coordinate system. However, thermodynamic state coordinate cannot be defined by space-time coordinate system.
Figure 1 shows the coordinate system in unified mechanics theory. Let us assume there is a 5-year-old boy and 100-year-old man. Using the space-time Cartesian coordinate system, their location can be defined by x , y , z coordinates and age on the time axis. However, this does not give any information about their thermodynamic state. Let us assume that a 5-year-old boy has stage 4 cancer is expected to die in a few days and a 100-year-old is expected to die in few days. This information cannot be represented in x , y , z -time- space coordinate system as shown in Figure 1. However, on TSI axis, 5-year-old boy and 100-year-old will have the same thermodynamic state index coordinate at Φ = 0.999 .
Another example can be given for Newton’s second law. If a soccer ball is given an initial acceleration with a force of F , it will move but eventually will come to a stop. Depending on the path it follows, it will come to a stop. Again, the initial acceleration of the ball is governed by the second law of Newton and slowing down process is governed by the laws of thermodynamics, which is represented by ( 1 Φ ) term. Detailed derivation of TSI can be seen in the literature [29]. We provide a simple summary in the following section.

2.1.3. Thermodynamic State Index (TSI) for Damage in Low Cycle Fatigue of Materials

Entropy and Helmholtz free energy are related by the thermodynamic principles [66] as follows:
Ψ = e T s
where Ψ represents the specific Helmholtz free energy, and e, T, s are the specific internal energy, temperature and specific entropy, respectively. Specific entropy is also related to the disorder parameter through Boltzmann’s equation [29,30] as follows:
s = k B ln ( W )
Total entropy for a volume can be given by:
S = N A k B   l n ( W ) m s
where, N A , k B , m s are the Avogadro number, Boltzmann’s constant and molar mass, respectively and W represent the disorder parameter [29,30,38,39,66]. Relation between the number of microstate, probability of microstates and disorder parameter is discussed extensively in the literature [70,71,72]. Using Equation (8), the TSI is given by:
Φ = Φ c ( 1 exp ( Δ s m s R ) )
where, Φ c , is a user defined parameter, representing the predefined failure criterion. R is gas constant. Δ s is a measure of the total change in entropy at a point. Unified mechanics theory states that when a system undergoes thermodynamic change from state A to state B, the remaining useful life can be defined by a factor in each stage of its life, called thermodynamic state index (TSI), Φ [0,1]. The ultimate failure is represented by a value of TSI equal to 1. Since, the value of Δ s is to be evaluated on the basis of mechanisms of dissipation processes involved in a thermodynamic process, the value of Φ c will be governed by a user-defined ultimate failure criterion.

2.2. Analytical Approach for the Prediction of Damage and Fatigue Life

From Equation (9), the TSI is governed by the change in entropy towards a predefined failure. All the dissipation processes that are related to failure lead to increase in entropy. Therefore, an appropriate measure of dissipation is needed to estimate the life of a process. In Ti-6Al-4V alloys, we consider only the mechanical process of dissipation, under monotonic as well as cyclic loading conditions. Hence the plastic dissipation is considered to be the dominant mechanism in the mechanical loading conditions. Entropy generation in plastic dissipation process can be calculated from a mechanical loading experiment in the following way:
Δ s = 1 ρ T t 1 t 2 σ : d ε p
where, ρ , is the mass density of the material, σ and ε p are the stress and plastic strain, respectively. T represents the temperature. Integral limits t 1 and t 2 represents the time bounds of the mechanical loading process, over which we quantify the change in entropy. For one dimensional case, the total plastic strain, ε p ( t ) is calculated as follows:
ε p ( t ) = ε t o t a l ( t ) σ y 0 E
where, ε t o t a l ( t ) is the total strain at the time of loading, t. σ y 0 and E are the yield stress and Young’s modulus, respectively. In the case of monotonic loading, the plastic dissipation is calculated from the engineering stress-strain graphs. In order to accomplish this, the plot is divided into elastic and plastic regime of loading. The area under the plastic region is computed by trapezoidal integration rule, and the cumulative entropy is evaluated in each stage. This accumulated entropy is used to predict the TSI at each and every strain level. A schematic representation of computing the incremental plastic dissipation is given in Figure 2. Accumulated entropy at n-th strain increment is computed from the Equation (10) as follows:
Δ s n = 1 ρ T i = 1 i = n σ i : Δ ε i p
Using Equations (9), (11) and (12), one dimensional approximation of damage measure is calculated under the assumptions that the damage is uniform within the cross section of the dog-bone test sample, and there are no other geometric or boundary effects in the sample. It is also assumed that the heat generation entropy production is small when compared with the entropy generation due to plastic deformation. In case of low cycle fatigue loading, the plastic dissipation is calculated as the area under the stress-strain hysteresis loop. Each cyclic hysteresis loop of engineering stress-strain graph, which represents the incremental dissipation. Hence, the accumulated entropy can be calculated by summation of incremental entropy. For a strain-controlled experiment, the accumulated entropy is a function of stress. Since the stress level at a given stage of cyclic loading is governed by the thermodynamic state index (TSI), Φ of the material, the TSI can be used to calculate the incremental dissipation from any known stage of loading, as follows:
Π i + 1 p = ( 1 Φ i ) Π i p
where, Π i p and Π i + 1 p represents the hysteresis area at i-th and (i+1)-th cyclic loading, respectively and Φ represents the TSI. Hence, the entropy change at any loading stage can be calculated from the initial loading hysteresis area as follows:
Δ s n = 1 ρ T i = 1 i = n Π i p
Φ i + 1 = Φ c [ 1 exp ( Δ s i m s R ) ]
It is important to point out that the entire thermodynamic response of the material point is mapped onto the TSI axes. Under no circumstances, the material point can exist outside the domain of [0,1]. The above approach has limitations that the one-dimensional approximation should be valid when the prediction is compared with experimental observations. To account for all the boundary and geometric effects related to stiffness, instabilities due to buckling, local cracking, stress concentrations, geometric nonlinearities, etc., we have developed a three-dimensional computational model. The detailed derivation is given in Section 2.3 below.

2.3. Computational 3-D Model for the Prediction of Damage

2.3.1. Derivation of the Computational Model

In this section, a three-dimensional model is derived, based on the unified mechanics theory. Entropy balance equation [4,20,29,30], can be written as follows:
d S d t d i v   J q T + 1 T σ : D ρ T d W e d t + ρ 𝓇 T
The following equation, as written in indicial notation, is known as Clausius-Duhem inequality [67,73,74]:
ɤ = 1 T σ i j D i j ρ T d W e d t 1 T 2 J q i T , i + ρ 𝓇 T 0
where, i and j are the indices, representing the spatial coordinates. ɤ is the specific entropy generation rate. σ denotes the stress tensor and T , i represents the spatial derivative of temperature, namely, the gradient of temperature. J q and 𝓇 , represents the heat flux transfer and internal heat generation, respectively. For small strain problems, rate of deformation tensor D is equal to strain rate tensor ε ˙ . According to Hooke’s law, stress is related to the strain through a constitutive tensor as follows:
σ i j = C i j k l ε k l e
where, C i j k l is the fourth order tangential constitutive tensor at the given stage of loading. ε k l e is the elastic part of strain tensor. Based on assumption of the additive decomposition of strain tensor [67], we can write the following equation for small strain problems:
ε i j t o t a l = ε i j e + ε i j p
where, ε i j t o t a l is the component of total strain tensor. For a given material point, based on unified mechanics theory one can write the following modified version of Equation (18), as follows:
σ i j = ( 1 Φ ) C i j k l 0 ε k l e
where the tangential constitutive tensor C i j k l is related to the virgin state of the same, C i j k l 0 (undamaged state) through TSI, Φ . For linear isotropic materials, undamaged constitutive tensor C i j k l 0 can be written as follows:
C i j k l 0 = λ δ i j δ k l + μ ( δ i k δ j l + δ i l δ j k )
where, λ and μ are the Lame’s parameters and δ i j is the identity tensor. The following inverse relations can also be written for a linear elastic isotropic material:
ε i j e = 1 + ν E σ i j ν E σ k k δ i j
where, E and ν are the elastic modulus and Poisson’s ratio, respectively. The rate form of the Equation (19), can be written as follows,
ε ˙ i j t o t a l = ε ˙ i j e + ε ˙ i j p
From incremental theory of plasticity, one can write the evolution equation for the fluxes, using the continuity of dissipation potential function, p (yield surface) [67] as follows:
ε ˙ i j p = Γ ˙ p σ i j
Effective stress at a point can be defined as follows:
σ i j = 1 1 Φ σ i j
where, σ i j is the component of effective stress tensor. Noting that Δ s is the only function that depends on time, the time rate of change of TSI can be obtained by differentiating Equation (9), yielding:
Φ ˙ = m s R Φ c Δ s ˙ ( exp ( Δ s m s R ) )
Assuming that the process is isothermal for each small load increment and all the dissipation mechanisms other than plastic deformation are negligibly small for the strain-controlled monotonic, quasi-static loading and low cycle fatigue loading in Ti-6Al-4V, we can write the entropy evolution as given in Equation (10). Hence the rate form of the entropy evolution from Equation (10) can be written as follows:
Δ ˙ s = 1 ρ T σ i j ε ˙ i j p
With the above assumption in the absence of kinematic hardening, we consider the following additive decomposition form of the Helmholtz free energy function as:
Ψ ( ε e , h ; Φ ) = Ψ E ( ε e ; Φ ) + Ψ I ( h )
where, Ψ E is the elastic strain energy and Ψ I is the free energy from isotropic hardening process. In the Equation (28), the hardening flux parameter h evolves with plastic strain. From the Equations (26) and (27), the plastic strain is a function of TSI.
Using Equations (20), (21) and (28), we get the following form of free energy:
Ψ ( ε e , h ; Φ ) = 1 2 ( 1 Φ ) ( λ ε k k e ε m m e + 2 μ ε i j e ε i j e ) + ( 1 Φ ) 1 r K h r + 1
We have assumed a power law model for isotropic hardening. Here, K and r are the material parameters which are to be found from the succeeding parts of the formulation and experimental data. The conjugate force is derived from Equation (29) as follows [75]:
σ = ρ ε e Ψ
σ i j = ( 1 Φ ) ( λ ε k k e δ i j + 2 μ ε i j e )
The yield function for Ti-6Al-4V can be given by:
p ( σ , H ; Φ ) = σ e q ( σ y o + H )
where, σ e q is the Von-Mises equivalent stress. σ y o represents the initial yield stress and H represents the hardening stress. Von-Mises equivalent stress is given by the following equation:
σ e q = 3 2 S i j S i j
where, the effective deviatoric stress tensor S i j , is given by the following equation,
S i j = σ i j σ k k 3 δ i j
Hence, from Equations (24) and (32), we get the following relation for plastic strain rate tensor, ε ˙ i j p ,:
ε ˙ i j p = Γ ˙ σ e q σ i j
Further simplification can be done on Equation (35) using the Equations (33) and (34). We get the following form for plastic strain rate tensor, based on normality rule of incremental theory of plasticity:
ε ˙ i j p = Γ ˙ [ 1 ( 1 Φ ) 3 2 S i j σ e q ]
where, Γ ˙ is the consistency parameter. By taking the norm of Equation (36) and by doing some algebra, we get the following equation to quantify the parameter, Γ ˙ :
ε ˙ e q p = 2 3 ε ˙ i j p ε ˙ i j p = Γ ˙ 1 ( 1 Φ )
Equation (37) is an important observation that the field variable, h , representing the isotropic hardening process, is related to the plastic deformation. Hence, we get the following form for h ˙ and ε ˙ i j p :
h ˙ = ε ˙ e q p ( 1 Φ )
ε ˙ i j p = ε ˙ e q p [ 3 2 S i j σ e q ]
From Equation (39), it can be observed that the magnitude of plastic strain is given by the equivalent plastic strain, ε e q p , and the direction of plastic loading is given by the term, [ 3 2 S i j σ e q ] .

2.3.2. Algorithm for the Computational Model

In this section, let us consider that all the variables having a superscript, ‘n’ represents values that are updated based on the previous loading and those variables with superscript, ‘n+1’ denotes the values corresponding to the current state of loading. All the variables having subscript, ‘tr’ represents the trial values. For simplicity in representation, indicial representation of the tensorial quantities are avoided.
Total strain at any increment is given by:
ε total n + 1 = ε total n + Δ ε total
Using Equation (19):
ε e t r n + 1 = ε total n + 1 ε p n
Using Equation (20):
σ t r n + 1 = ( 1 Φ n ) C 0 ε e t r n + 1
σ n + 1 = ( 1 Φ n + 1 ) C 0 ( ε total n + 1 ε p n + 1 )
Let:
w = ( 1 Φ )
then:
σ n + 1 = w n + 1 C 0 ( ε total n + 1 ε p n Δ ε p )
Using Equations (25), (36), (40), (41), and (44) in Equation (45), we get the following:
σ n + 1 = 1 w n σ t r n + 1 1 w n + 1 C 0 Δ Γ [ 3 2 S n + 1 σ e q ]
Let:
p = σ k k 3
Therefore, from Equations (34), (46) and (47), we can write the following expression:
S n + 1 + 1 w n p t r   I 1 3 C 0 Δ Γ w n + 1   [ 3 2 S n + 1 σ e q n + 1 ] = 1 w n σ t r n + 1 1 w n + 1 C 0 Δ Γ [ 3 2 S n + 1 σ e q n + 1 ]
Using Equation (21) in (48), we get the simplified form for the iteration equation in indicial notation, as follows:
σ e q n + 1 2 { δ i k δ j l + 1 w n + 1 C i j k l 0 Δ Γ [ 3 2 1 σ e q n + 1 ] } { δ i k δ j l + 1 w n + 1 C i j k l 0 Δ Γ [ 3 2 1 σ e q n + 1 ] } = σ e q t r n + 1 2 w n 2
Algorithmically derived Equation (49) can be solved by an iteration procedure to find the value of Δ Γ , simultaneously with the update of w . A Newton-Raphson iteration scheme is employed in the integration scheme of the present study to solve the yield function given in Equation (32). Successively, the entropy is updated using Equation (27) and the damage is calculated using Equation (15).

3. Validation of the Computational Model for Monotonic Loading

Prior to the simulation of fully reversed cyclic loading, it is important to check the validity of the model under tensile as well as compressive loading. The computation models described in Section 2.3, is implemented in commercial finite element package, ABAQUS. User material subroutine is written to update the stresses according to the strain increments that are supplied to the subroutine as input. In order to validate the model for tensile as well as compressive loading cases in Ti-6Al-4V, we have used the experimental data, reported by Biswas et al. [2] and Carrion et al. [76].

3.1. Validation of the Numerical Model for Monotonic Tensile Loading

The true stress-strain graph reported in the literature [76] for Ti-6Al-4V alloy, is used for the comparison between experimental data and the numerical predictions of monotonic tensile loading. Mill Annealed hot rolled bars were used [76] in the study. The material parameters are taken from the literature [76], so as to match with the material used for the comparison. Details of the model parameters are given in Table 1. Using the common assumption that the gauge section of a dog bone sample experiences uniform strain, we consider 5 mm length in the computational model. Hence, it can reduce the computational cost as well. Diameter of the specimen is kept the same, like that of the experimentally reported sample by Carrion et al. [76], which is 6.35 mm in diameter. In ABAQUS, linear brick elements, C3D8R are used to mesh the numerical model. One end of the sample is defined with zero displacement (fixed) boundary condition and the other end is subjected to controlled displacement loading in the axial direction. After a mesh convergence analysis, an optimum seed size of 0.9 mm is fixed for all the simulations. A schematic representation of the computational geometry is shown in Figure 3.
It can be observed from the Figure 4 that the true stress-strain graph, predicted for monotonic tensile loading in Ti-6Al-4V alloy, matches well with the experimental observations reported by Carrion et al. [76]. A smooth transition can be seen at point A, shown in Figure 4. This transition from elastic to plastic region can be due to the dislocation motion in the microstructure. Further, dislocation multiplication and interaction with each other and inclusions can be the possible reason behind strain hardening of the bulk material. Hence, the validation of the model under tensile loading can be considered as a basis for tensile loading in any kind of geometry or boundary conditions in the numerical investigation.
A comparative plot between the numerical results for damage obtained from three-dimensional model and analytical results based on one dimensional approximation, as described in Section 2.2, is shown in Figure 5. It is observed that the level of matching between computational and experimental results for monotonic tensile loading is closer in the case of prediction of damage, based on the analytical approach and numerical analysis.

3.2. Validation of the 3-D Numerical Model for Monotonic Compressive Loading

Validation of the 3-D numerical model is done under compressive loading as well. Experimental result for a monotonic compression test, reported in the literature [2] is used to validate the proposed numerical model. The computational model parameters are taken from the literature [2], so as to match with the material used for the comparison. Even though the reported experimental results [2,76] are for Ti-6Al-4V alloys, it is noted that the materials are different in terms of their mechanical properties. Details of the model parameters used for the numerical simulation of monotonic compression test are listed in Table 2. We have considered the same dimensions in the numerical model, as that of the experimental samples [2]. Since, true stress-strain data is given in the literature [2], analytical procedure to compute TSI, requires an additional step. This method is adopted from well-known damage rule based on area reduction [77]. In the current study, damage parameter is represented by the TSI. Hence, the current area is related to the original area of undamaged section through the factor, TSI as follows:
A = ( 1 Φ ) A 0
where, A and A 0 represents the current area and initial area. The engineering stress and true stress are related by the principle of static equilibrium as follows,
σ ¯ A = σ ¯ A 0
where, σ ¯ and σ ¯ represents the true stress and engineering stress respectively. Hence, in order to quantify the entropy, we have estimated the engineering yield stress data as follows:
σ ¯ y i + 1 = σ ¯ y i + 1 ( 1 Φ i )
where, σ ¯ y i + 1 and σ ¯ y i + 1 represents the computed engineering stress and true stress at (i+1)-th strain, respectively. Φ i is calculated based on the i-th strain data. Hence, in an incremental way, the TSI is computed using analytical procedure given in Section 2.2. Computation model in ABAQUS is discretized with linear brick finite elements C3D8R. One of the ends of the computational model is constrained from all the translations and the other end is subjected to displacement controlled compressive loading in the axial direction. A mesh convergence analysis is conducted and an optimum seed size of 0.9 mm is adopted in the simulations. A schematic representation of the computational geometry is shown in Figure 6.
Numerical results for monotonic compressive loading in Ti-6Al-4V alloy, shown in Figure 7, are found to be matching well with the reported experimental results [2]. Hence, the proposed model is taken as a basis to simulate compressive loading cases in the succeeding numerical investigations. Using the experimental [2] stress-strain graph, we have analytically calculated the TSI at every incremental plastic strain, based on the procedure stated in Section 2.2. As shown in Figure 8, both the analytical and numerical predictions for TSI matches very well.

4. Model Predictions for Low Cycle Fatigue Life

Carrion et al. [76] tested Ti-6Al-4V samples under tensile loading condition at a strain rate of the order of 10−3 s−1 at room temperature. Similar quasi-static loading condition is established in our numerical loading by controlling the step time of the numerical model in ABAQUS. The material model used in developing the 3-Dimensional numerical model is independent of the strain rate and the temperature and hence the strain rate hardening behavior and temperature effects, including the thermal dissipation are not considered in our study. Unified mechanics theory-based approach for damage calculation, described in Section 2, is used to predict the low cycle fatigue life of Ti-6Al-4V alloys. Details of the one-dimensional analytical model as well as the three-dimensional numerical model to predict fatigue life or Ti-6Al-4V are given in Section 4.1 below.

4.1. Analytical Approach for Fatigue Life Prediction

Experimental results [76] for the stabilized hysteresis loop is assumed to be closer to the first cycle hysteresis loop. Unified mechanics theory is used to evaluate the damage evolution under cyclic loading and the results are plotted in Figure 9.
Low cycle fatigue life of the Ti-6Al-4V sample is predicted by fixing the TSI at failure as 0.98. This is necessary, as to prevent computational instabilities at the verge of failure that are not recorded by experimental results, are to be taken into account when we compare the mathematical model predictions with the experimental results. A MATLAB script is written to compute the fatigue life, from the stabilized hysteresis loop. The results are shown in comparison with the test data [76] and the corresponding numerical predictions at similar amplitudes, as shown in Figure 10.

4.2. Computational Procedure for Fatigue Life Prediction

It is not feasible to conduct a large number of cyclic loading in the numerical model to predict fatigue life, especially when the amplitude of strain is very small because in ABAQUS this process would take weeks. In this section, we propose an alternate way of fatigue life prediction of Ti-6Al-4V alloys at different strain amplitudes, using a combined numerical-experimental procedure. If the hysteresis loop for a given strain amplitude is found out from the experiment, the same test is simulated by using the proposed model. Computational results after the first cycle of loading are used to find the scaling factor for incremental entropy in the computational model. The scaling factor is calculated as the ratio between the experimental hysteresis loop area for the stabilized loop and the numerically computed dissipation for the first cycle of loading. Then the computational model is used to evaluate the dissipation at different strain amplitudes of loading for a single cycle of loading. This hysteresis loop is used to predict the fatigue life at the given amplitude of strain, as per the procedure detailed in Section 2.3.
To compare the numerical predictions for fatigue life with experimental results [76], the same material data, as listed in Table 1, are used. It is assumed that the experimental results are free from any boundary effects or instabilities. Hence, the numerical analysis is done on the sample, with dimensions and boundary conditions as shown in Figure 3. Hysteresis loops at 1.2% strain amplitude for 50 cycles of loading are plotted in Figure 11a. A comparative hysteresis plot for first cycle and 50th cycle of loading is shown in Figure 11b. It can be observed from Figure 11a,b, that the hysteresis loop area decreases with cyclic loading. This reduction in hysteresis loop area is due to the reduction in strength of the material with the evolution of TSI. The fatigue life can be predicted by extrapolating the numerical results on TSI axis vs number of cycles. A comparison plot between test data and simulations for low cycle fatigue life prediction at different strain amplitudes is shown in Figure 10. In Figure 10, the average values of fatigue life test data [76] are plotted for stain amplitudes of 0.8%, 1.0% and 1.2% and compared with the analytical predictions. Fatigue life test data for other amplitudes of strain are not reported in the literature [76]. Results from the numerical approach, for the strain amplitudes 1.0% and 1.2% are also plotted and the model prediction is extended to a strain amplitude of 2.4%. Response at 0.8% strain amplitude was not computed with 3-D model because 0.8% strain amplitude is within the elastic region of loading.
In Figure 10, it is clear that the one-dimensional analytical approach is underestimating the fatigue life by 1600 cycles at a total strain amplitude of 0.8%, while this discrepancy is less scattered in the test data [76]. This discrepancy in fatigue life prediction using one-dimensional model could be due to the unaccounted three-dimensional confinement effects in material response.

5. Conclusions

The work presented here is based on the unified mechanics theory, where the laws of Newtonian mechanics are combined with laws of thermodynamics, directly. The bridging factor in unified mechanics theory is the definition of thermodynamic state index, given in the Equation (9). The definition of damage proposed in the literature [29], is applied in the case of monotonic as well as low cycle fatigue loading conditions in Ti-6Al-4V alloys. Based on the principles of continuum mechanics, we have presented a numerical model, which account for the damage in case of plastic loading in Ti-6Al-4V. It is observed from the three-dimensional numerical and one-dimensional analytical results of the damage model prediction that, they match very well with the experimental observations in the case of monotonic tensile loading, as shown in Figure 4 and Figure 5. In Figure 4, we have considered the stress-strain graph given in the literature [76] for validation and the corresponding damage prediction (value of TSI is around 0.1) is limited to a strain level of about 3%. Linear extrapolation of the damage curves plotted in Figure 5, can lead to wrong prediction of failure strain (to around 20% in the current study). Entropy at each time increment is dependent on the stress level. Hence, the accuracy of life prediction will be dependent on the constitutive model, used to predict the yielding of the material, in a three-dimensional numerical study.
The monotonic compressive stress-strain graph is matching well with the experimental results reported in the literature, as shown in Figure 7. The path traced by the damage prediction from one-dimensional analytical procedure and three-dimensional numerical procedure, as shown in Figure 8, also matches very well. In the case of compressive loading conditions, the results can be affected by the confining effects. The difference in nature of path traced by damage curves in compressive and tensile loading conditions could be due to the difference in confining effects in compressive loading, when compared with tensile loading. Similar observations for alloys can be seen in the literature [78,79]. Current study may be extended in future, for the detailed experimental and numerical investigations on such confining effects, under compressive loading. Since the current focus of the investigation is to introduce an efficient way of predicting the fatigue life of Ti-6Al-4V using computational tools in conjunction with the experiment, we have limited our study to fatigue life prediction.
Thermodynamics of life of any system, as postulated by the unified mechanics theory, is brought in to application level, for the case of low cycle fatigue life prediction in Ti-6Al-4V. From the comparative study on fatigue life prediction, as shown in Figure 10, the proposed procedures, described in Section 4, are found to be very efficient. Only one cycle experimental data is sufficient to predict the low cycle fatigue in Ti-6Al-4V alloys. Hence, the procedure stated in Section 4, will be useful for practical applications.

Author Contributions

Conceptualization, N.B.J.M. and A.K.; Methodology, N.B.J.M.; Software, N.B.J.M.; Validation, N.B.J.M. and A.K.; Formal Analysis, N.B.J.M.; Investigation, N.B.J.M., A.K., C.L.R. and C.B.; Resources, C.L.R. and C.B.; Data Curation, N.B.J.M.; Writing-Original Draft Preparation, N.B.J.M.; Writing-Review & Editing, N.B.J.M., A.K., C.L.R. and C.B.; Supervision, C.L.R. and C.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mouritz, A.P. Introduction to Aerospace Materials; Woodhead Publishing Limited: Cambridge, UK, 2012; ISBN 9781845695323. [Google Scholar]
  2. Biswas, N.; Ding, J.L.; Balla, V.K.; Field, D.P.; Bandyopadhyay, A. Deformation and fracture behavior of laser processed dense and porous Ti6Al4V alloy under static and dynamic loading. Mater. Sci. Eng. A 2012, 549, 213–221. [Google Scholar] [CrossRef]
  3. Banerjee, D.; Williams, J.C. Perspectives on titanium science and technology. Acta Mater. 2013, 61, 844–879. [Google Scholar] [CrossRef]
  4. Altenberger, I.; Nalla, R.K.; Sano, Y.; Wagner, L.; Ritchie, R.O. On the effect of deep-rolling and laser-peening on the stress-controlled low- and high-cycle fatigue behavior of Ti-6Al-4V at elevated temperatures up to 550 °C. Int. J. Fatigue 2012, 44, 292–302. [Google Scholar] [CrossRef]
  5. CHENG, A.S.; LAIRD, C. Fatigue Life Behavior of Copper Single Crystals. Part I: Observations of Crack Nucleation. Fatigue Fract. Eng. Mater. Struct. 1981, 4, 331–341. [Google Scholar] [CrossRef]
  6. Ren, Y.M.; Lin, X.; Guo, P.F.; Yang, H.O.; Tan, H.; Chen, J.; Li, J.; Zhang, Y.Y.; Huang, W.D. Low cycle fatigue properties of Ti-6Al-4V alloy fabricated by high-power laser directed energy deposition: Experimental and prediction. Int. J. Fatigue 2019, 127, 58–73. [Google Scholar] [CrossRef]
  7. Tanaka, K.; Mura, T. Dislocation Model for Fatigue Crack Initiation. Am. Soc. Mech. Eng. 1981, 48, 97–103. [Google Scholar] [CrossRef]
  8. Guo, Q.; Zaïri, F.; Guo, X. An intrinsic dissipation model for high-cycle fatigue life prediction. Int. J. Mech. Sci. 2018, 140, 163–171. [Google Scholar] [CrossRef]
  9. Sosnovskiy, L.A.; Senko, V.I. Tribo-fatigue. In Proceedings of the ASME International Mechanical Engineering Congress and Exposition, Tribology, Orlando, FL, USA, 5–11 November 2005; pp. 141–148. [Google Scholar]
  10. Wang, T.; Samal, S.K.; Lim, S.K.; Shi, Y. Entropy Production-Based Full-Chip Fatigue Analysis: From Theory to Mobile Applications. IEEE Trans. Comput. Des. Integr. Circuits Syst. 2019, 38, 84–95. [Google Scholar] [CrossRef]
  11. Hashin, Z. A Relnterpretation of the Palmgren- miner Rule for Fatigue Life. J. Appl. Mech. 2016, 47, 324–328. [Google Scholar] [CrossRef] [Green Version]
  12. Smith, K.N.; Topper, T.H.; Watson, P. A stress–strain function for the fatigue of metals (stress-strain function for metal fatigue including mean stress effect). J. Mater. 1970, 5, 767–778. [Google Scholar]
  13. Lemaitre, J.; Desmorat, R. Engineering damage mechanics: ductile, creep, fatigue and brittle failures; Springer: Berlin/Heidelberg, Germany; New York, NY, USA, 2005; ISBN 3540215034. [Google Scholar]
  14. Bhattacharya, B.; Ellingwood, B. Continuum damage mechanics analysis of fatigue crack initiation. Int. J. Fatigue 1998, 20, 631–639. [Google Scholar] [CrossRef]
  15. Ontiveros, V.; Amiri, M.; Kahirdeh, A.; Modarres, M. Thermodynamic entropy generation in the course of the fatigue crack initiation. Fatigue Fract. Eng. Mater. Struct. 2017, 40, 423–434. [Google Scholar] [CrossRef]
  16. Harvey, S.E.; Marsh, P.G.; Gerberich, W.W. Atomic force microscopy and modeling of fatigue crack initiation in metals. Acta Metall. Mater. 1994, 42, 3493–3502. [Google Scholar] [CrossRef]
  17. Kumar, J.; Sundara Raman, S.G.; Kumar, V. Analysis and Modeling of Thermal Signatures for Fatigue Damage Characterization in Ti–6Al–4V Titanium Alloy. J. Nondestruct. Eval. 2016, 35, 1–10. [Google Scholar] [CrossRef]
  18. Sosnovskiy, L.A.; Sherbakov, S.S. Mechanothermodynamic entropy and analysis of damage state of complex systems. Entropy 2016, 18, 268. [Google Scholar] [CrossRef] [Green Version]
  19. Kahirdeh, A.; Khonsari, M.M. Energy dissipation in the course of the fatigue degradation: Mathematical derivation and experimental quantification. Int. J. Solids Struct. 2015, 77, 74–85. [Google Scholar] [CrossRef]
  20. Amiri, M.; Naderi, M.; Khonsari, M.M. An experimental approach to evaluate the critical damage. Int. J. Damage Mech. 2011, 20, 89–112. [Google Scholar] [CrossRef]
  21. Sosnovskiy, L.A.; Sherbakov, S.S. Mechanothermodynamical system and its behavior. Contin. Mech. Thermodyn. 2012, 24, 239–256. [Google Scholar] [CrossRef]
  22. Zhang, M.H.; Shen, X.H.; He, L.; Zhang, K.-S. Application of Differential Entropy in Characterizing the Deformation Inhomogeneity and Life Prediction of Low-Cycle Fatigue of Metals. Materials 2018, 11, 1917. [Google Scholar] [CrossRef] [Green Version]
  23. Sosnovskiy, L.A.; Sherbakov, S.S. A Model of Mechanothermodynamic Entropy in Tribology. Entropy 2017, 19, 115. [Google Scholar] [CrossRef]
  24. Sosnovskiy, L.; Sherbakov, S. Mechanothermodynamics; Springer: Berlin/Heidelberg, Germany, 2016; ISBN 978-3-319-24979-7. [Google Scholar]
  25. Young, C.; Subbarayan, G. Maximum Entropy Models for Fatigue Damage in Metals with Application to Low-Cycle Fatigue of Aluminum 2024-T351. Entropy 2019, 21, 967. [Google Scholar] [CrossRef] [Green Version]
  26. Santecchia, E.; Hamouda, A.M.S.; Musharavati, F.; Zalnezhad, E.; Cabibbo, M.; El Mehtedi, M.; Spigarelli, S. A Review on Fatigue Life Prediction Methods for Metals. Adv. Mater. Sci. Eng. 2016, 2016, 9573524. [Google Scholar] [CrossRef] [Green Version]
  27. Haddad, W.M. Thermodynamics: The unique universal science. Entropy 2017, 19, 621. [Google Scholar] [CrossRef] [Green Version]
  28. Haddad, W.M.; Chellabonia, V.; Nersesov, S.G. Thermodynamics: A Dynamical Systems Approach; Princeton University Press: Princeton, NJ, USA, 2005; ISBN 0-691-12327-6. [Google Scholar]
  29. Basaran, C.; Yan, C.Y. A thermodynamic framework for damage mechanics of solder joints. J. Electron. Packag. Trans. ASME 1998, 120, 379–384. [Google Scholar] [CrossRef] [Green Version]
  30. Basaran, C.; Nie, S. An irreversible thermodynamics theory for damage mechanics of solids. Int. J. Damage Mech. 2004, 13, 205–223. [Google Scholar] [CrossRef]
  31. Basaran, C.; Tang, H. Implementation of a thermodynamic framework for damage mechanics of solder interconnect in microelectronic packaging. ASME Int. Mech. Eng. Congr. Expo. Proc. 2002, 11, 61–68. [Google Scholar]
  32. Basaran, C.; Lin, M.; Ye, H. A thermodynamic model for electrical current induced damage. Int. J. Solids Struct. 2003, 40, 1738–1745. [Google Scholar] [CrossRef]
  33. Gomez, J.; Basaran, C. A thermodynamics based damage mechanics constitutive model for low cycle fatigue analysis of microelectronics solder joints incorporating size effects. Int. J. Solids Struct. 2005, 42, 3744–3772. [Google Scholar] [CrossRef]
  34. Gomez, J.; Basaran, C. Damage mechanics constitutive model for Pb/Sn solder joints incorporating nonlinear kinematic hardening and rate dependent effects using a return mapping integration algorithm. Mech. Mater. 2006, 38, 585–598. [Google Scholar] [CrossRef]
  35. Tang, H.; Basaran, C. A damage mechanics-based fatigue life prediction model for solder joints. J. Electron. Packag. Trans. ASME 2003, 125, 120–125. [Google Scholar] [CrossRef]
  36. Temfack, T.; Basaran, C. Experimental verification of thermodynamic fatigue life prediction model using entropy as damage metric. Mater. Sci. Technol. 2015, 31, 1627–1632. [Google Scholar] [CrossRef]
  37. Wang, J.; Yao, Y. An entropy based low-cycle fatigue life prediction model for solder materials. Entropy 2017, 19, 503. [Google Scholar]
  38. Basaran, C.; Zhao, Y.; Tang, H.; Gomez, J. A damage-mechanics-based constitutive model for solder joints. J. Electron. Packag. Trans. ASME 2005, 127, 208–214. [Google Scholar] [CrossRef]
  39. Basaran, C.; Chandaroy, R. Thermomechanical analysis of solder joints under thermal and vibration loading. J. Electron. Packag. Trans. ASME 2002, 124, 60–66. [Google Scholar] [CrossRef]
  40. Basaran, C.; Li, S.; Abdulhamid, M.F. Thermomigration induced degradation in solder alloys. J. Appl. Phys. 2008, 103, 123520-1–123520-9. [Google Scholar] [CrossRef] [Green Version]
  41. Basaran, C.; Lin, M. Damage mechanics of electromigration induced failure. Mech. Mater. 2008, 40, 66–79. [Google Scholar] [CrossRef]
  42. Basaran, C.; Lin, M. Damage mechanics of electromigration in microelectronics copper interconnects. Int. J. Mater. Struct. Integr. 2007, 1, 16–39. [Google Scholar] [CrossRef]
  43. Basaran, C.; Nie, S. Time dependent behavior of a particle filled composite PMMA/ATH at elevated temperatures. J. Compos. Mater. 2008, 42, 2003–2025. [Google Scholar] [CrossRef]
  44. Basaran, C.; Nie, S. A thermodynamics based damage mechanics model for particulate composites. Int. J. Solids Struct. 2007, 44, 1099–1114. [Google Scholar] [CrossRef] [Green Version]
  45. Li, S.; Basaran, C. A computational damage mechanics model for thermomigration. Mech. Mater. 2009, 41, 271–278. [Google Scholar] [CrossRef] [Green Version]
  46. Lin, M.; Basaran, C. Electromigration induced stress analysis using fully coupled mechanical-diffusion equations with nonlinear material properties. Comput. Mater. Sci. 2005, 34, 82–98. [Google Scholar] [CrossRef]
  47. Shidong, L.; Abdulhamid, M.F.; Basaran, C. Simulating Damage Mechanics of Electromigration and Thermomigration. Simulation 2008, 84, 391–401. [Google Scholar] [CrossRef]
  48. Yao, W.; Basaran, C. Electromigration damage mechanics of lead-free solder joints under pulsed DC: A computational model. Comput. Mater. Sci. 2013, 71, 76–88. [Google Scholar] [CrossRef]
  49. Yao, W.; Basaran, C. Computational damage mechanics of electromigration and thermomigration. J. Appl. Phys. 2013, 114, 103708. [Google Scholar] [CrossRef] [Green Version]
  50. Cuadras, A.; Crisóstomo, J.; Ovejas, V.J.; Quilez, M. Irreversible entropy model for damage diagnosis in resistors. J. Appl. Phys. 2015, 118, 165103-1–165103-8. [Google Scholar] [CrossRef] [Green Version]
  51. Cuadras, A.; Romero, R.; Ovejas, V.J. Entropy characterisation of overstressed capacitors for lifetime prediction. J. Power Sources 2016, 336, 272–278. [Google Scholar] [CrossRef] [Green Version]
  52. Cuadras, A.; Yao, J.; Quilez, M. Determination of LEDs degradation with entropy generation rate. J. Appl. Phys. 2017, 122, 145702-1–145702-7. [Google Scholar] [CrossRef]
  53. Imanian, A.; Modarres, M. A thermodynamic entropy approach to reliability assessment with applications to corrosion fatigue. Entropy 2015, 17, 6995–7020. [Google Scholar] [CrossRef] [Green Version]
  54. Imanian, A.; Modarres, M. A thermodynamic entropy-based damage assessment with applications to prognostics and health management. Struct. Heal. Monit. 2018, 17, 1–15. [Google Scholar] [CrossRef]
  55. Jang, J.Y.; Khonsari, M.M. On the evaluation of fracture fatigue entropy. Theor. Appl. Fract. Mech. 2018, 96, 351–361. [Google Scholar] [CrossRef]
  56. Liakat, M.; Khonsari, M.M. Entropic characterization of metal fatigue with stress concentration. Int. J. Fatigue 2015, 70, 223–234. [Google Scholar] [CrossRef]
  57. Osara, J.A.; Bryant, M.D. A Thermodynamic Model for Lithium-Ion Battery Degradation: Application of the Degradation-Entropy Generation Theorem. Inventions 2019, 4, 23. [Google Scholar] [CrossRef] [Green Version]
  58. Osara, J.A.; Bryant, M.D. Thermodynamics of Fatigue: Degradation-Entropy Generation Methodology for System and Process Characterization and Failure Analysis. Entropy 2019, 21, 685. [Google Scholar] [CrossRef] [Green Version]
  59. Gomez, J.; Lin, M.; Basaran, C. Damage Mechanics Modeling of Concurrent Thermal and Vibration Loading on Electronics Packaging. Multidiscip. Model. Mater. Struct. 2006, 2, 309–326. [Google Scholar] [CrossRef]
  60. Amiri, M.; Khonsari, M.M. On the role of entropy generation in processes involving fatigue. Entropy 2012, 14, 24–31. [Google Scholar] [CrossRef] [Green Version]
  61. Wang, J.; Yao, Y. An entropy-based failure prediction model for the creep and fatigue of metallic materials. Entropy 2019, 21, 1104. [Google Scholar] [CrossRef] [Green Version]
  62. Sun, F.; Zhang, W.; Wang, N.; Zhang, W. A copula entropy approach to dependence measurement for multiple degradation processes. Entropy 2019, 21, 724. [Google Scholar] [CrossRef] [Green Version]
  63. Yun, H.; Modarres, M. Measures of Entropy to Characterize Fatigue Damage in Metallic Materials. Entropy 2019, 21, 804. [Google Scholar] [CrossRef] [Green Version]
  64. Li, E.H.; Li, Y.Z.; Li, T.T.; Li, J.X.; Zhai, Z.Z.; Li, T. Intelligent analysis algorithm for satellite health under time-varying and extremely high thermal loads. Entropy 2019, 21, 983. [Google Scholar] [CrossRef] [Green Version]
  65. Sosnovskiy, L.A.; Sherbakov, S.S. On the Development of Mechanothermodynamics as a New Branch of Physics. Entropy 2019, 21, 1188. [Google Scholar] [CrossRef] [Green Version]
  66. Gunel, E.M.; Basaran, C. Damage characterization in non-isothermal stretching of acrylics. Part I: Theory. Mech. Mater. 2011, 43, 979–991. [Google Scholar] [CrossRef]
  67. Naderi, M.; Amiri, M.; Khonsari, M.M. On the thermodynamic entropy of fatigue fracture. Proc. R. Soc. A Math. Phys. Eng. Sci. 2010, 466, 423–438. [Google Scholar] [CrossRef] [Green Version]
  68. Naderi, M.; Khonsari, M.M. An experimental approach to low-cycle fatigue damage based on thermodynamic entropy. Int. J. Solids Struct. 2010, 47, 875–880. [Google Scholar] [CrossRef] [Green Version]
  69. Abdulhamid, M.F.; Basaran, C. Influence of thermomigration on lead-free solder joint mechanical properties. J. Electron. Packag. Trans. ASME 2009, 131, 011002. [Google Scholar] [CrossRef]
  70. Boltzmann, L. Ableitung des Stefan’schen Gesetzes, betreffend die Abhängigkeit der Wärmestrahlung von der Temperatur aus der electromagnetischen Lichttheorie. Ann. Phys. 1884, 258, 291–294. [Google Scholar] [CrossRef]
  71. Sharp, K.; Matschinsky, F. Translation of Ludwig Boltzmann’s paper “on the relationship between the second fundamental theorem of the mechanical theory of heat and probability calculations regarding the conditions for thermal equilibrium” Sitzungberichte der kaiserlichen akademie d. Entropy 2015, 17, 1971–2009. [Google Scholar] [CrossRef] [Green Version]
  72. Planck, M. On the Law of Distribution of Energy in the Normal Spectrum. Ann. Phys. 1901, 4, 553. [Google Scholar] [CrossRef]
  73. Lemaitre, J.; Chaboche, J.-L. Mechanics of Solid Materials; Cambridge University Press: Cambridge, UK, 1990; ISBN 0-521-32853-5. [Google Scholar]
  74. Voyiadjis, G.Z.; Faghihi, D. Thermo-mechanical strain gradient plasticity with energetic and dissipative length scales. Int. J. Plast. 2012, 30–31, 218–247. [Google Scholar] [CrossRef]
  75. Murakami, S. Continuum Damage Mechanics; Springer: Berlin/Heidelberg, Germany, 2012; ISBN 9789400726659. [Google Scholar]
  76. Carrion, P.E.; Shamsaei, N.; Daniewicz, S.R.; Moser, R.D. Fatigue behavior of Ti-6Al-4V ELI including mean stress effects. Int. J. Fatigue 2017, 99, 87–100. [Google Scholar] [CrossRef] [Green Version]
  77. Murakami, S. Mechanical modeling of material damage. J. Appl. Mech. Trans. ASME 1988, 55, 280–286. [Google Scholar] [CrossRef]
  78. Spitzig, W.A.; Sober, R.J.; Richmond, O. Pressure dependence of yielding and associated volume expansion in tempered martensite. Acta Metall. 1975, 23, 885–893. [Google Scholar] [CrossRef]
  79. Mahnken, R. Strength difference in compression and tension and pressure dependence of yielding in elasto-plasticity. Comput. Methods Appl. Mech. Eng. 2001, 190, 5057–5080. [Google Scholar] [CrossRef]
Figure 1. Coordinate system in unified mechanics theory.
Figure 1. Coordinate system in unified mechanics theory.
Entropy 22 00024 g001
Figure 2. Schematic representation of computing plastic dissipation from the engineering stress-plastic strain graph.
Figure 2. Schematic representation of computing plastic dissipation from the engineering stress-plastic strain graph.
Entropy 22 00024 g002
Figure 3. Schematics of numerical model for displacement controlled monotonic tensile loading in ABAQUS.
Figure 3. Schematics of numerical model for displacement controlled monotonic tensile loading in ABAQUS.
Entropy 22 00024 g003
Figure 4. Comparison between monotonic tensile stress-strain graphs obtained from the test data [76] and numerical model.
Figure 4. Comparison between monotonic tensile stress-strain graphs obtained from the test data [76] and numerical model.
Entropy 22 00024 g004
Figure 5. Comparison between the damage (TSI) prediction for monotonic tensile loading.
Figure 5. Comparison between the damage (TSI) prediction for monotonic tensile loading.
Entropy 22 00024 g005
Figure 6. Schematics of numerical model for displacement controlled monotonic compressive. loading in ABAQUS.
Figure 6. Schematics of numerical model for displacement controlled monotonic compressive. loading in ABAQUS.
Entropy 22 00024 g006
Figure 7. Comparison between monotonic compressive stress-strain graphs obtained from the test data [2] and numerical model.
Figure 7. Comparison between monotonic compressive stress-strain graphs obtained from the test data [2] and numerical model.
Entropy 22 00024 g007
Figure 8. Comparison between the damage prediction for monotonic compressive loading.
Figure 8. Comparison between the damage prediction for monotonic compressive loading.
Entropy 22 00024 g008
Figure 9. Analytical prediction of damage for different strain amplitudes of cyclic loading.
Figure 9. Analytical prediction of damage for different strain amplitudes of cyclic loading.
Entropy 22 00024 g009
Figure 10. Low cycle fatigue life (Nf) prediction at different strain amplitudes of cyclic loading in comparison with the test data [76].
Figure 10. Low cycle fatigue life (Nf) prediction at different strain amplitudes of cyclic loading in comparison with the test data [76].
Entropy 22 00024 g010
Figure 11. Numerical results on engineering stress-strain hysteresis loops for 1.2% strain amplitude of cyclic loading. (a) hysteresis loops at 1.2% strain amplitude for 50 cycles of loading; (b) comparative hysteresis plot for first cycle and 50th cycle of loading.
Figure 11. Numerical results on engineering stress-strain hysteresis loops for 1.2% strain amplitude of cyclic loading. (a) hysteresis loops at 1.2% strain amplitude for 50 cycles of loading; (b) comparative hysteresis plot for first cycle and 50th cycle of loading.
Entropy 22 00024 g011
Table 1. Material parameters used in the numerical model for tensile loading in Ti-6Al-4V alloy.
Table 1. Material parameters used in the numerical model for tensile loading in Ti-6Al-4V alloy.
Material ParameterValueUnit
Young’s modulus, E106GPa
Poisson’s ratio, ν0.31
Density, ρ 4540kg/m3
Critical TSI, Φ c 1
Hardening parameter, K 968.00MPa
Hardening exponent, r 0.64
Yield strength, σ y 0 992.00MPa
Molar mass, m s 0.047867kg/mol
Reference temperature, T298K
Table 2. Material parameters used in the numerical model for compressive loading in Ti-6Al-4V alloy.
Table 2. Material parameters used in the numerical model for compressive loading in Ti-6Al-4V alloy.
Material ParameterValueUnit
Young’s modulus, E118GPa
Poisson’s ratio, ν0.31
Density, ρ 4540kg/m3
Critical TSI, Φ c 1
Hardening parameter, K 550.00MPa
Hardening exponent, r 0.65
Yield strength, σ y 0 1047.00MPa
Molar mass, m s 0.047867kg/mol
Reference temperature, T298K

Share and Cite

MDPI and ACS Style

Bin Jamal M, N.; Kumar, A.; Lakshmana Rao, C.; Basaran, C. Low Cycle Fatigue Life Prediction Using Unified Mechanics Theory in Ti-6Al-4V Alloys. Entropy 2020, 22, 24. https://doi.org/10.3390/e22010024

AMA Style

Bin Jamal M N, Kumar A, Lakshmana Rao C, Basaran C. Low Cycle Fatigue Life Prediction Using Unified Mechanics Theory in Ti-6Al-4V Alloys. Entropy. 2020; 22(1):24. https://doi.org/10.3390/e22010024

Chicago/Turabian Style

Bin Jamal M, Noushad, Aman Kumar, Chebolu Lakshmana Rao, and Cemal Basaran. 2020. "Low Cycle Fatigue Life Prediction Using Unified Mechanics Theory in Ti-6Al-4V Alloys" Entropy 22, no. 1: 24. https://doi.org/10.3390/e22010024

APA Style

Bin Jamal M, N., Kumar, A., Lakshmana Rao, C., & Basaran, C. (2020). Low Cycle Fatigue Life Prediction Using Unified Mechanics Theory in Ti-6Al-4V Alloys. Entropy, 22(1), 24. https://doi.org/10.3390/e22010024

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop