Next Article in Journal
Optimization of Flow Field in Slab Continuous Casting Mold with Medium Width Using High Temperature Measurement and Numerical Simulation for Automobile Exposed Panel Production
Previous Article in Journal
Leaching Kinetics of Arsenic Sulfide-Containing Materials by Copper Sulfate Solution
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Reply

Reply to Comment on the Paper Entitled “A New Cumulative Fatigue Damage Rule Based on Dynamic Residual S-N Curve and Material Memory Concept” by Peng Z., Huang H., Zhou J. and Li Y. Published in Metals (2018; 8(6): 456)

1
School of Mechanical and Electrical Engineering, University of Electronic Science and Technology of China, Chengdu 611731, China
2
Center for System Reliability and Safety, University of Electronic Science and Technology of China, Chengdu 611731, China
*
Author to whom correspondence should be addressed.
Metals 2020, 10(1), 8; https://doi.org/10.3390/met10010008
Submission received: 9 December 2019 / Accepted: 16 December 2019 / Published: 20 December 2019
Dear Editors,
For the published paper entitled “A new cumulative fatigue damage rule based on dynamic residual S-N curve and material memory concept” (Metals (2018; 8(6): 456)), Kris Hectors and Prof. Wim De Waele give a comment about the incorrect use of the formula of fatigue life prediction. Thanks for the suggestion. The following is our reply for the comment.
According to the comment “Equation (34) implicitly assumes that the number of cycles of each load block can be scaled with the same scalar up to failure, corresponding to D i =   1 and n i =   N pre . This is however only true if the damage accumulation is described as a linear function of cycle ratios.”
This description is not true in some special cases. Considering a repeated two-level load block (with two blocks to failure), as shown in Figure 1.
It can also be viewed as a four-level fatigue loading, the specimen is fatigued by σ1 for n1 cycles, σ2 for n2 cycles, σ3 for n3 cycles, and σ4 for n4 cycles to failure, where σ 1 = σ 1 , σ 2 = σ 2 , σ 3 = σ 1 , σ 4 = σ 2 , n 1 = n 1 , n 2 = n 2 , n 3 = n 1 , n 4 = n 2 .
As is known, for a linear damage accumulation model, the damage variable Di should be a linear function of cycle ratio n i / N f i , that is:
D i     n i N f i
In Equation (1), for Miner rule, D i = n i / N fi . For a nonlinear damage model, Di is defined as a nonlinear function of the cycle ratio, such as Marco-Starkey’s model ( D i = ( n i / N fi ) a i ), Chaboche’s continuum damage mechanics model (the formula can be found in Ref [47]).
(1)
Kwofie’s model
For Kwofie’s model, Di is described as:
D i   =   n i N f i × ln ( N f i ) ln ( N f 1 )
where ln ( N fi ) / ln ( N f 1 ) is a constant without n i / N fi . Di is a linear function of cycle ratio n i / N fi and that it is a linear damage model.
Fatigue damage is accumulated by taking a summation of the segmental damage caused by each loading stress level.
(a)
The first-level loading
D 1   =   n 1 N f 1 = n 1 N f 1
(b)
The second-level loading
D 2   =   n 2 N f 2 × ln ( N f 2 ) ln ( N f 1 ) = n 2 N f 2 × ln ( N f 2 ) ln ( N f 1 )
(c)
The third-level loading
D 3   =   n 3 N f 3 × ln ( N f 3 ) ln ( N f 1 ) = n 1 N f 1 × ln ( N f 1 ) ln ( N f 1 ) = n 1 N f 1
(d)
The fourth-level loading
D 4   =   n 4 N f 4 × ln ( N f 4 ) ln ( N f 1 ) = n 2 N f 2 × ln ( N f 2 ) ln ( N f 1 )
The cumulative damage D i is expressed as:
D i   = D 1 + D 2 + D 3 + D 4 = n 1 N f 1 + n 2 N f 2 × ln ( N f 2 ) ln ( N f 1 ) + n 1 N f 1 + n 2 N f 2 × ln ( N f 2 ) ln ( N f 1 ) = 2 × n 1 N f 1 + 2 × n 2 N f 2 × ln ( N f 2 ) ln ( N f 1 )
Using the Equation (34), the predicted fatigue life is
N pre   = n i D i = n 1 + n 2 + n 3 + n 4 D i = 2 × n 1 + 2 × n 2 2 × n 1 N f 1 + 2 × n 2 N f 2 × ln ( N f 2 ) ln ( N f 1 ) = n 1 + n 2 n 1 N f 1 + n 2 N f 2 × ln ( N f 2 ) ln ( N f 1 )
Equation (8) can be scaled with a scalar C (C = 2).
(2)
A modified Kwofie’s model
Similar to Equation (2), we change the damage variable as:
D i   =   n i N f i × ln ( N f i ) ln ( N f ( i 1 ) )
Compared to Equation (2), ln ( N f 1 ) becomes ln ( N f ( i 1 ) ) , and ln ( N fi ) / ln ( N f ( i 1 ) ) is a constant without of n i / N fi . Di is still a linear function of cycle ratio n i / N f i .
(a)
The first-level loading
D 1   =   n 1 N f 1 = n 1 N f 1
(b)
The second-level loading
D 2   =   n 2 N f 2 × ln ( N f 2 ) ln ( N f 1 ) = n 2 N f 2 × ln ( N f 2 ) ln ( N f 1 )
(c)
The third-level loading
D 3   =   n 3 N f 3 × ln ( N f 3 ) ln ( N f 2 ) = n 1 N f 1 × ln ( N f 1 ) ln ( N f 2 )
(d)
The fourth-level loading
D 4   =   n 4 N f 4 × ln ( N f 4 ) ln ( N f 3 ) = n 2 N f 2 × ln ( N f 2 ) ln ( N f 1 )
The cumulative damage D i is expressed as:
D i   = D 1 + D 2 + D 3 + D 4 = n 1 N f 1 + n 2 N f 2 × ln ( N f 2 ) ln ( N f 1 ) + n 1 N f 1 × ln ( N f 1 ) ln ( N f 2 ) + n 2 N f 2 × ln ( N f 2 ) ln ( N f 1 ) = n 1 N f 1 + n 1 N f 1 × ln ( N f 1 ) ln ( N f 2 ) + 2 × n 2 N f 2 × ln ( N f 2 ) ln ( N f 1 )
Using the Equation (34), the predicted fatigue life is:
N pre   =   n i D i = n 1 + n 2 + n 3 + n 4 D i = 2 × n 1 + 2 × n 2 n 1 N f 1 + n 1 N f 1 × ln ( N f 1 ) ln ( N f 2 ) + 2 × n 2 N f 2 × ln ( N f 2 ) ln ( N f 1 ) 2 × n 1 + 2 × n 2 2 × n 1 N f 1 + 2 × n 2 N f 2 × ln ( N f 2 ) ln ( N f 1 ) = n 1 + n 2 n 1 N f 1 + n 2 N f 2 × ln ( N f 2 ) ln ( N f 1 )
Equation (15) cannot be scaled with a scalar because ln ( N f 1 ) / ln ( N f 2 ) 1 . As a result, for the linear damage model, it is uncertain whether the number of cycles of each load block can be scaled. This is attributed to ln ( N f ( i 1 ) ) in Equation (9). Although ln ( N f ( i 1 ) ) is a constant irrelevant with cycle ratio, it is changed with previous fatigue life.
(3)
Proposed model
For the proposed model, Di is described as:
D i   =   n i N f i × j = 1 i 1 ( N f j N f ( j + 1 ) ) ( k = 1 j α k ) 1
where α k = ( e n k / N fk e 1 ) / ( 1 e 1 ) .
In Equation (16), the loading effect coefficient j = 1 i 1 ( N fj N f ( j + 1 ) ) ( k = 1 j α k ) 1 (using the symbol λ i instead) is associated with fatigue lives N f 1 , N f 2 , …, N fi and previous cycle ratios n 1 / N f 1 , n 2 / N f 2 , …, n i 1 / N f ( i 1 ) but without n i / N f i . If these previous cycle ratios are determined, then λ i becomes a determined value as a constant. The damage variable Di will be a linear function of cycle ratio n i / N f i , in reality it is a linear damage model.
(a)
The first-level loading
D 1   =   n 1 N f 1 = n 1 N f 1
(b)
The second-level loading
D 2   =   n 2 N f 2 × λ 2 = n 2 N f 2 × λ 2
(c)
The third-level loading
D 3   =   n 3 N f 3 × λ 3 = n 1 N f 1 × λ 3
(d)
The fourth-level loading
D 4   =   n 4 N f 4 × λ 4 = n 2 N f 2 × λ 4
The cumulative damage D i is expressed as:
D i   = D 1 + D 2 + D 3 + D 4 = n 1 N f 1 + n 2 N f 2 × λ 2 + n 1 N f 1 × λ 3 + n 2 N f 2 × λ 4
Using the Equation (34), the predicted fatigue life is
N pre   =   n i D i = n 1 + n 2 + n 3 + n 4 D i = 2 × n 1 + 2 × n 2 n 1 N f 1 + n 2 N f 2 × λ 2 + n 1 N f 1 × λ 3 + n 2 N f 2 × λ 4 = 2 × n 1 + 2 × n 2 ( 1 + λ 3 ) × n 1 N f 1 + ( λ 2 + λ 4 ) × n 2 N f 2 n 1 + n 2 n 1 N f 1 + n 2 N f 2
Equation (22) cannot be scaled with a scalar because 1 + λ 3 2 and λ 2 + λ 4 2 . This is attributed to λ i in Equation (16). λ i is a constant without of the current cycle radio n i / N f i , it is changed and determined by previous cycle ratios from n 1 / N f 1 to n i 1 / N f ( i 1 ) . It should be pointed out that the loading effect coefficient λ i reveals a significant correlation with previous fatigue loading histories. This finding is an essential distinction with the Miner rule ( λ i = 1 ).
(4)
The reason for using Equation (34) to predict fatigue life
In the paper, the prediction of fatigue life contain two steps:
(a)
The first step:
The cumulative fatigue damage i = 1 k D i is first calculated using different fatigue damage models. It should be emphasized that i = 1 k D i denotes the predicted cumulative damage and its calculated value may not be equal to unity for different models. For the experimental value of cumulative damage, it is always defined as unity without using damage models, i.e., i = 1 k D i = 1 .
(b)
The second step:
According to the calculated values of cumulative damage obtained by different models, the fatigue life is thus calculated by Equation (34).
For the experimental result, the cumulative fatigue damage is unity. For different models, if the predicted cumulative damage is unity (i.e., D i 1 ), then the predicted fatigue life should be close to the experimental result (i.e., NpreNexp) and the corresponding damage model becomes more effective.
For the comment in Table 1, D i denotes the predicted cumulative damage by different models but not the real value of experimental result ( D i = 1 ). In theory, the predicted results cannot be equal to the experimental ones. In row (a), the specimen is fatigued by σ1 for 13,749 cycles and σ2 for 51,304 cycles to failure. This is correct because it denotes the real experimental results. However, in row (b) and row (c), the number of loading cycles of n1 and n2 are defined as the hypothetical values. The specimen is supposed first fatigued by σ1 for 18,454 (or 17,497) cycles and then σ2 for 68,860 (or 65,290) cycles. This situation is not correct (the specimen may be failure when n2 = 50,000 or other value, and n2 can only be determined by the experiments). In addition, if n1 and n2 are changed, the predicted damage (D1 and D2) should also be changed by the cumulative damage models as well as the corresponding predicted fatigue lives. It is meaningless for supposing a particular value of the loading cycles.
Moreover, for Corten’s model as stated by Equation (29), the prediction formula of fatigue life can be found in the paper entitled “A Practical Method for Determining the Corten-Dolan Exponent and Its Application to Fatigue Life Prediction” by Shun-Peng Zhu, Hong-Zhong Huang, Yu Liu, Li-Ping He and Qiang Liao published in Int. J. Turbo Jet-Engines (2012;29: 79–87), on page 81, in Equation (5), that is:
n r N g   =   i = 1 k n i N 1 ( σ i σ 1 ) d
where n r = n 1 + n 2 + + n k = i = 1 k n i , N g is the predicted fatigue life.
The predicted fatigue life can be written as:
N p r e = N g = n r i = 1 k n i N 1 ( σ i σ 1 ) d   = i = 1 k n i i = 1 k D i
This formula is the same as the Equation (34) in our paper.
As stated previously, the use of the Equation (34) for fatigue life prediction should be reasonable.

Author Contributions

Formal analysis, J.Z.; Funding acquisition, Y.-F.L.; Methodology, Z.P.; Supervision, H.-Z.H.; Validation, Y.-F.L.; Writing-original draft, Z.P.; Writing-review & editing, H.-Z.H. All authors have read and agreed to the published version of the manuscript.

Acknowledgments

This study was sponsored by the National Natural Science Foundation of China under Grant No. 51775090.

Conflicts of Interest

The authors declare no conflict of interest.
Figure 1. Repeated two-level load block (two blocks).
Figure 1. Repeated two-level load block (two blocks).
Metals 10 00008 g001
Table 1. Predicted lifetime calculated by Peng et al. for their non-linear damage model and two block loading sequences (a), with the corresponding damage (b) and the corrected lifetime prediction (c).
Table 1. Predicted lifetime calculated by Peng et al. for their non-linear damage model and two block loading sequences (a), with the corresponding damage (b) and the corrected lifetime prediction (c).
High–low loading sequence σ 1 = 485 M P a σ 2 = 400 M P a
n 1 / n 2 n 1 n 2 n i α 1 D i N p r e
(a) N p r e Peng et al.0.26813,74951,30465,0530.65010.745087,314
(b) Corresponding damage0.26818,45468,86087,3140.54911.0687
(c) Correct N p r e 0.26817,49765,29082,7870.56891.000082,786
Low–high loading sequence σ 1 = 400 M P a σ 2 = 485 M P a
(a) N p r e Peng et al.2.341109,31046,693156,0030.16531.1264138,502
(b) Corresponding damage2.34197,04741,455138,5020.23091.0221
(c) Correct N p r e 2.34194,48240,359134,8410.24531.0000134,847

Share and Cite

MDPI and ACS Style

Peng, Z.; Huang, H.-Z.; Zhou, J.; Li, Y.-F. Reply to Comment on the Paper Entitled “A New Cumulative Fatigue Damage Rule Based on Dynamic Residual S-N Curve and Material Memory Concept” by Peng Z., Huang H., Zhou J. and Li Y. Published in Metals (2018; 8(6): 456). Metals 2020, 10, 8. https://doi.org/10.3390/met10010008

AMA Style

Peng Z, Huang H-Z, Zhou J, Li Y-F. Reply to Comment on the Paper Entitled “A New Cumulative Fatigue Damage Rule Based on Dynamic Residual S-N Curve and Material Memory Concept” by Peng Z., Huang H., Zhou J. and Li Y. Published in Metals (2018; 8(6): 456). Metals. 2020; 10(1):8. https://doi.org/10.3390/met10010008

Chicago/Turabian Style

Peng, Zhaochun, Hong-Zhong Huang, Jie Zhou, and Yan-Feng Li. 2020. "Reply to Comment on the Paper Entitled “A New Cumulative Fatigue Damage Rule Based on Dynamic Residual S-N Curve and Material Memory Concept” by Peng Z., Huang H., Zhou J. and Li Y. Published in Metals (2018; 8(6): 456)" Metals 10, no. 1: 8. https://doi.org/10.3390/met10010008

APA Style

Peng, Z., Huang, H. -Z., Zhou, J., & Li, Y. -F. (2020). Reply to Comment on the Paper Entitled “A New Cumulative Fatigue Damage Rule Based on Dynamic Residual S-N Curve and Material Memory Concept” by Peng Z., Huang H., Zhou J. and Li Y. Published in Metals (2018; 8(6): 456). Metals, 10(1), 8. https://doi.org/10.3390/met10010008

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