Analyzing the effect of residual stresses on the fatigue life of colddrawn steel wire specimens

Danielle Lopes Zaidan, Paulo Pedro Kenedi, Luís Felipe Guimarães de SouzaPostgraduate Program in Mechanical Engineering and Materials Technology - PPEMM, CEFET-RJ, Braz̧il danielle.zaidan@aluno.cefet-rj.br, http://orcid.org/0009-0008-7141-8275paulo.kenedi@cefet-rj.br, https:// orcid.org/ 0000-0001-5563-843Xluis.souza@cefet-rj.br, https://orcid.org/0000-0002-6669-8451

Introduction

Fatigue and residual stress are responsible for many mechanical failures. Although there are many areas where fatigue can be applied, such as civil engineering, Laterza et al. [1] or bioengineering, Teoh [2], the focus of this manuscript will be set on mechanical engineering. Most of the residual stress literature references are experimental, as in Tabatabaeian et al. [3]. In this manuscript, however, the residual stress approach will be developed in an analytical way. These two important failure causes, one widely known as fatigue and the other lesser known as residual stress, will be combined to form a compound analytical model. A step-by-step explanation will be given for the residual stress analytic approach, leaving the fatigue basics available in the Appendix. Residual stress can develop in numerous ways in manufacturing processes, such as welding or machining. Any procedure that partially yields a mechanical part cross-section can generate residual stress. Textbooks about the mechanics of solids, such as Crandall et al. [4], present, as a basic subject, the residual stress distribution in static beams submitted to one semi-cycle bending moment, for an elastic perfectly plastic material. Although analytical approaches can be used to gain insight into cross-sectional residual stress distribution, they usually handle only simple geometries, such as rectangular or circular cross-sections. Stok and Halilovic [5] show an elastoplastic analytical solution for bending beams with rectangular cross-sections. Rimovskis and Sabaliauskas[6] analyzed, for rectangular and circular beam cross-sections, the drift of the neutral axis relative to the centroidal axis of the beam as a function of dissimilar yielding parameters in tension and compression. Castro and Meggiolaro [7] and Jirasek and Bazant[8] presented improvements in residual stress analytical models. De Castro et al. [9] show different cases of analytical models' applications to estimate the residual stress distribution of partially yielded structures of rectangular cross sections. Experimental techniques were used to estimate residual stress on surfaces, as in Schajer [10]. Atienza and Elices [11] analyze the residual stress in cold-drawn wires experimentally and numerically, through the influence of residual stresses on the tensile test, where the residual stress was relieved based on mechanical and/or thermomechanical treatments. Yet, the availability of the combination of fatigue and residual stress, in technical literature, is quite uncommon. One exception is the interesting manuscript of McClung [12], that provides a broad literature survey of the current understanding of significant residual stress issues for fatigue lifetime for shot peening, cold expansion of holes, welding, and machining. Also, see Webster et al. [13], James et al. [14] and Vaara et al. [15] for additional references.
This article proposes an analytical model to describe the effect of residual stress on the fatigue life of partially yielded colddrawn steel wires, modeled as a bi-linear material, of bulged rectangular cross-section. This research is a development of former references, such as Zaidan et al. [16]. The results obtained in this manuscript can contribute to improving the knowledge about the effect of residual stress on fatigue performance of partially yielded cold-drawn steel wires, extensively used, for instance, in the offshore industry.

Analytical model

The compound analytical model, which intends to reproduce the experimental part of the manuscript, is divided into residual stress and fatigue approaches, which refer, respectively, to the experimental part of phases 1 and 2 . In phase 1, a monotonic displacement is imposed on the specimen, as shown in Fig.1.b. A circular forming tool acts transversally on the specimen, supported by a 3 -point bending gripper, yielding partially its cross-section. After spring-back, as shown in Fig.1.c, a residual stress distribution pattern is generated in its cross-section. Phase 1 is described by the residual stress approach of this analytical model, which is considered a relevant contribution, represented by Eqns. (1-11). In phase 2 a variable transverse load P ( t ) P ( t ) P(t)\mathrm{P}(\mathrm{t})P(t) is applied on the specimen that has already passed by phase 1 , turned upside down, as shown in Fig. 1.d. The variable transverse load is applied until the specimen completely fails. The phase 2 is described by the fatigue approach of this analytical model. The applied fatigue theory is the classical one, so it was placed in the Appendix section, as shown in Eqns. (1A-7A). The compound analytical model, as shown in Eqns. (12-22), merges two different phenomena, residual stress and fatigue, to produce a comprehensive analysis of the effect of residual stresses on the fatigue life of colddrawn steel wire specimens. Mathcad software is used to perform calculations.
Figure 1: (a) The forming tool, with radius 1 1 _(1)_{1}1, touches the specimen, (b) a transversal displacement δ δ delta\deltaδ is applied towards the specimen, (c) the specimen, after spring back, presents radius 2 2 _(2){ }_{2}2, (d) the specimen is turned upside down and the fatigue loading is applied, (e) the specimen bulged rectangle cross-section.
Where radius 1 1 _(1)_{1}1 is the forming tool radius. Radius 2 2 _(2){ }_{2}2 is the radius of curvature of the specimen after spring back. Radius 3 3 _(3){ }_{3}3 is the radius of the specimen's lateral cross-section. L L LLL is the length of the specimen between the 3-point bending gripper external rollers. δ δ delta\deltaδ is transversal displacement applied towards the specimen. P ( t ) P ( t ) P(t)\mathrm{P}(\mathrm{t})P(t) is the fatigue loading. b 1 b 1 b_(1)\mathrm{b}_{1}b1 and b 2 b 2 b_(2)\mathrm{b}_{2}b2 are, respectively, the minor and the major width of the specimen's cross-section, y y yyy is a vertical distance and d y d y dyd ydy is a differential of vertical distance. And c c ccc is the semi-height and b b bbb is the height of the specimen's cross-section. y y y y y_(y)\mathrm{y}_{\mathrm{y}}yy is the elastic-plastic border.

Residual stress basics

Every step of the residual stress part of the analytical model is explicitly presented and commented on. The value of the bending moment M ( ϱ i ) M ( ϱ i ) M(ϱi)\mathrm{M}(\varrho \mathrm{i})M(ϱi) imposed by the forming tool is inside the range M y < M ( ϱ i ) < M p M y < M ( ϱ i ) < M p M_(y) < M(ϱi) < M_(p)\mathrm{M}_{\mathrm{y}}<\mathrm{M}(\varrho \mathrm{i})<\mathrm{M}_{\mathrm{p}}My<M(ϱi)<Mp. The maximum elastic moment M y M y M_(y)\mathrm{M}_{\mathrm{y}}My and the plastic moment M p M p M_(p)\mathrm{M}_{\mathrm{p}}Mp equations are presented.
(01) M y = 2 { 0 c y y c S y 0 2 radius 3 2 y 2 d y } radius 3 = c 2 + ( b 1 2 ) 2 (02) M p = 2 { 0 c y ( S y 0 + E t ( y ρ i ε y ) ) 2 radius 3 2 y 2 d y } (03) ρ i = radius 1 + c radius 2 = c | ε a | ρ y = c ε y ρ p = c ε u E t = S u t 0 S y 0 ε u ε y (01) M y = 2 0 c y y c S y 0 2 radius 3 2 y 2 d y radius 3 = c 2 + b 1 2 2 (02) M p = 2 0 c y S y 0 + E t y ρ i ε y 2 radius 3 2 y 2 d y (03) ρ i = radius 1 + c radius 2 = c ε a ρ y = c ε y ρ p = c ε u E t = S u t 0 S y 0 ε u ε y {:[(01)M_(y)=2*{int_(0)^(c)y*(y)/(c)*S_(y0)*2*sqrt(radius_(3)^(2)-y^(2))dy}quadradius_(3)=sqrt(c^(2)+((b_(1))/(2))^(2))],[(02)M_(p)=2*{int_(0)^(c)y*(S_(y0)+E_(t)*((y)/(rho_(i))-epsi_(y)))*2*sqrt(radius_(3)^(2)-y^(2))dy}],[(03)rho_(i)=radius_(1)+c quadradius_(2)=(c)/(|epsi_(a)|)quadrho_(y)=(c)/(epsi_(y))quadrho_(p)=(c)/(epsi_(u))quadE_(t)=(S_(ut0)-S_(y0))/(epsi_(u)-epsi_(y))]:}\begin{align*} & M_{y}=2 \cdot\left\{\int_{0}^{c} y \cdot \frac{y}{c} \cdot S_{y 0} \cdot 2 \cdot \sqrt{\operatorname{radius}_{3}^{2}-y^{2}} d y\right\} \quad \operatorname{radius}_{3}=\sqrt{c^{2}+\left(\frac{b_{1}}{2}\right)^{2}} \tag{01}\\ & M_{p}=2 \cdot\left\{\int_{0}^{c} y \cdot\left(S_{y 0}+E_{t} \cdot\left(\frac{y}{\rho_{i}}-\varepsilon_{y}\right)\right) \cdot 2 \cdot \sqrt{\operatorname{radius}_{3}^{2}-y^{2}} d y\right\} \tag{02}\\ & \rho_{i}=\operatorname{radius}_{1}+c \quad \operatorname{radius}_{2}=\frac{c}{\left|\varepsilon_{a}\right|} \quad \rho_{y}=\frac{c}{\varepsilon_{y}} \quad \rho_{p}=\frac{c}{\varepsilon_{u}} \quad E_{t}=\frac{S_{u t 0}-S_{y 0}}{\varepsilon_{u}-\varepsilon_{y}} \tag{03} \end{align*}(01)My=2{0cyycSy02radius32y2dy}radius3=c2+(b12)2(02)Mp=2{0cy(Sy0+Et(yρiεy))2radius32y2dy}(03)ρi=radius1+cradius2=c|εa|ρy=cεyρp=cεuEt=Sut0Sy0εuεy
Where S y 0 S y 0 S_(y0)\mathrm{S}_{\mathrm{y} 0}Sy0 is the material yielding resistance and S ut 0 S ut 0 S_(ut0)\mathrm{S}_{\mathrm{ut} 0}Sut0 is the material's ultimate resistance, both, as received, E t E t E_(t)\mathrm{E}_{\mathrm{t}}Et is the plastic tangent modulus, ρ i , ρ y ρ i , ρ y rho_(i),rho_(y)\rho_{i}, \rho_{y}ρi,ρy, and ρ p ρ p rho_(p)\rho_{p}ρp are the radius of curvature of the specimen, respectively, at loading, at the beginning of yielding, and the full cross-section yielding. The ε a ε a epsi_(a)\varepsilon_{\mathrm{a}}εa is the strain at point a a aaa of Fig 1.a, ε y ε y epsi_(y)\varepsilon_{\mathrm{y}}εy is the yielding strain, and ε u ε u epsi_(u)\varepsilon_{\mathrm{u}}εu is the ultimate strain. The elastic-plastic border y y y y y_(y)y_{y}yy equation is presented:
(04) y y = ε y ρ i ρ p < ρ i < ρ y (04) y y = ε y ρ i ρ p < ρ i < ρ y {:(04)y_(y)=epsi_(y)*rho_(i)quadrho_(p) < rho_(i) < rho_(y):}\begin{equation*} y_{y}=\varepsilon_{y} \cdot \rho_{i} \quad \rho_{p}<\rho_{i}<\rho_{y} \tag{04} \end{equation*}(04)yy=εyρiρp<ρi<ρy
The equation of the bending moment M ( ρ i ) M ρ i M(rho_(i))M\left(\rho_{i}\right)M(ρi) in function of the imposed radius of curvature of the specimen:
M ( ρ i ) = 2 { 0 y y y ( y ε y ρ i S y 0 ) 2 radius 3 2 y 2 d y + y y c y ( S y 0 + E t ( y ρ i ε y ) ) 2 radius 3 2 y 2 d y } (05) M y < M ( ρ i ) < M p M ρ i = 2 0 y y y y ε y ρ i S y 0 2  radius  3 2 y 2 d y + y y c y S y 0 + E t y ρ i ε y 2  radius  3 2 y 2 d y (05) M y < M ρ i < M p {:[M(rho_(i))=2*{int_(0)^(y_(y))y*((y)/(epsi_(y)*rho_(i))*S_(y0))*2*sqrt(" radius "_(3)^(2)-y^(2))dy+int_(y_(y))^(c)y*(S_(y0)+E_(t)*((y)/(rho_(i))-epsi_(y)))*2*sqrt(" radius "_(3)^(2)-y^(2))dy}],[(05)M_(y) < M(rho_(i)) < M_(p)]:}\begin{align*} M\left(\rho_{i}\right)= & 2 \cdot\left\{\int_{0}^{y_{y}} y \cdot\left(\frac{y}{\varepsilon_{y} \cdot \rho_{i}} \cdot S_{y 0}\right) \cdot 2 \cdot \sqrt{\text { radius }_{3}^{2}-y^{2}} d y+\int_{y_{y}}^{c} y \cdot\left(S_{y 0}+E_{t} \cdot\left(\frac{y}{\rho_{i}}-\varepsilon_{y}\right)\right) \cdot 2 \cdot \sqrt{\text { radius }_{3}^{2}-y^{2}} d y\right\} \\ & M_{y}<M\left(\rho_{i}\right)<M_{p} \tag{05} \end{align*}M(ρi)=2{0yyy(yεyρiSy0)2 radius 32y2dy+yycy(Sy0+Et(yρiεy))2 radius 32y2dy}(05)My<M(ρi)<Mp
The stress σ σ sigma\sigmaσ and strain E E E\mathcal{E}E evolution of the specimen point a a aaa (shown in Fig. 1.a) is presented in Fig. 2.a, whose values are shown in Eqn. (06), Castro and Meggiolaro [6].
(06) σ 0 = 0 ε 0 = 0 σ 1 = S y 0 ε 1 = ε y σ 2 = S y 0 + E t ( y ρ i ε y ) ε 2 = ( y y c ) 1 ε y σ 3 = S y 0 + E t ( y ρ i ε y ) + σ s b ( c ) ε 3 = [ ( y y c ) 1 α ] ε y (06) σ 0 = 0 ε 0 = 0 σ 1 = S y 0 ε 1 = ε y σ 2 = S y 0 + E t y ρ i ε y ε 2 = y y c 1 ε y σ 3 = S y 0 + E t y ρ i ε y + σ s b ( c ) ε 3 = y y c 1 α ε y {:(06){:[sigma_(0)=0,epsi_(0)=0,],[sigma_(1)=S_(y0),epsi_(1)=epsi_(y),],[sigma_(2)=S_(y0)+E_(t)*((y)/(rho_(i))-epsi_(y)),epsi_(2)=((y_(y))/(c))^(-1)*epsi_(y)],[sigma_(3)=S_(y0)+E_(t)*((y)/(rho_(i))-epsi_(y))+sigma_(sb)(-c),epsi_(3)=[((y_(y))/(c))^(-1)-alpha]*epsi_(y)]:}:}\begin{array}{lll} \sigma_{0}=0 & \varepsilon_{0}=0 & \\ \sigma_{1}=S_{y 0} & \varepsilon_{1}=\varepsilon_{y} & \\ \sigma_{2}=S_{y 0}+E_{t} \cdot\left(\frac{y}{\rho_{i}}-\varepsilon_{y}\right) & \varepsilon_{2}=\left(\frac{y_{y}}{c}\right)^{-1} \cdot \varepsilon_{y} \\ \sigma_{3}=S_{y 0}+E_{t} \cdot\left(\frac{y}{\rho_{i}}-\varepsilon_{y}\right)+\sigma_{s b}(-c) & \varepsilon_{3}=\left[\left(\frac{y_{y}}{c}\right)^{-1}-\alpha\right] \cdot \varepsilon_{y} \tag{06} \end{array}(06)σ0=0ε0=0σ1=Sy0ε1=εyσ2=Sy0+Et(yρiεy)ε2=(yyc)1εyσ3=Sy0+Et(yρiεy)+σsb(c)ε3=[(yyc)1α]εy
where, α = M ( ϱ i ) / M y α = M ϱ i / M y alpha=M(ϱ_(i))//M_(y)\alpha=\mathrm{M}\left(\varrho_{\mathrm{i}}\right) / \mathrm{M}_{\mathrm{y}}α=M(ϱi)/My is a multiplicative factor. σ s b ( c ) σ s b ( c ) sigma_(sb)(-c)\sigma_{s b}(-c)σsb(c) is the springback stress at specimen point a a aaa, of Fig 1.a, estimated by Eqn. (10). The Eqns. (07) and (08) represent, respectively, in Fig. 2.a, the material behavior (light gray dash lines) and the loading/unloading equations (black lines for point a a aaa, and gray lines for point b b bbb ), where points a a aaa and b b bbb are shown in Fig 1.a.
(07) σ s ( ε ) = { S y 0 + E t ( ε + ε y ) , ε < ε y S y 0 ε y ε , ε y ε ε y S y 0 + E t ( ε ε y ) , ε > ε y (07) σ s ( ε ) = S y 0 + E t ε + ε y , ε < ε y S y 0 ε y ε , ε y ε ε y S y 0 + E t ε ε y , ε > ε y {:(07)sigma_(s)(epsi)={[-S_(y0)+E_(t)*(epsi+epsi_(y))",",epsi < -epsi_(y)],[(S_(y0))/(epsi_(y))*epsi",",-epsi_(y) <= epsi <= epsi_(y)],[S_(y0)+E_(t)*(epsi-epsi_(y))",",epsi > epsi_(y)]:}:}\sigma_{s}(\varepsilon)=\left\{\begin{array}{lc} -S_{y 0}+E_{t} \cdot\left(\varepsilon+\varepsilon_{y}\right), & \varepsilon<-\varepsilon_{y} \tag{07}\\ \frac{S_{y 0}}{\varepsilon_{y}} \cdot \varepsilon, & -\varepsilon_{y} \leq \varepsilon \leq \varepsilon_{y} \\ S_{y 0}+E_{t} \cdot\left(\varepsilon-\varepsilon_{y}\right), & \varepsilon>\varepsilon_{y} \end{array}\right.(07)σs(ε)={Sy0+Et(ε+εy),ε<εySy0εyε,εyεεySy0+Et(εεy),ε>εy
(08) σ ( ε ) = { S y 0 ε y ε ε 0 ε ε 1 S y 0 + E t ( ε ε y ) ε 1 ε ε 2 σ 3 σ 2 ε 3 ε 2 ε + σ 2 ε 3 σ 3 ε 2 ε 3 ε 2 ε 2 ε ε 3 (08) σ ( ε ) = S y 0 ε y ε ε 0 ε ε 1 S y 0 + E t ε ε y ε 1 ε ε 2 σ 3 σ 2 ε 3 ε 2 ε + σ 2 ε 3 σ 3 ε 2 ε 3 ε 2 ε 2 ε ε 3 {:(08)sigma(epsi)={[(S_(y0))/(epsi_(y))*epsi,epsi_(0) <= epsi <= epsi_(1)],[S_(y0)+E_(t)*(epsi-epsi_(y)),epsi_(1) <= epsi <= epsi_(2)],[(sigma_(3)-sigma_(2))/(epsi_(3)-epsi_(2))*epsi+(sigma_(2)*epsi_(3)-sigma_(3)*epsi_(2))/(epsi_(3)-epsi_(2)),epsi_(2) <= epsi <= epsi_(3)]:}:}\sigma(\varepsilon)= \begin{cases}\frac{S_{y 0}}{\varepsilon_{y}} \cdot \varepsilon & \varepsilon_{0} \leq \varepsilon \leq \varepsilon_{1} \tag{08}\\ S_{y 0}+E_{t} \cdot\left(\varepsilon-\varepsilon_{y}\right) & \varepsilon_{1} \leq \varepsilon \leq \varepsilon_{2} \\ \frac{\sigma_{3}-\sigma_{2}}{\varepsilon_{3}-\varepsilon_{2}} \cdot \varepsilon+\frac{\sigma_{2} \cdot \varepsilon_{3}-\sigma_{3} \cdot \varepsilon_{2}}{\varepsilon_{3}-\varepsilon_{2}} & \varepsilon_{2} \leq \varepsilon \leq \varepsilon_{3}\end{cases}(08)σ(ε)={Sy0εyεε0εε1Sy0+Et(εεy)ε1εε2σ3σ2ε3ε2ε+σ2ε3σ3ε2ε3ε2ε2εε3
In Fig. 2.a, in the final state of point a a aaa, labeled with the number 3, there is a change in the stress signal from its former state, labeled with the number 2, while the strain signal remains the same signal. Also, in Fig. 2.a the abscissa axis was maintained as strain (not divided by ε y 0 ε y 0 epsi_(y0)\varepsilon_{\mathrm{y} 0}εy0 ) to show that at the end of phase 1 , point a a aaa has a residual strain of 0.06 (or 6 % 6 % 6%6 \%6% ). The Eqns. (0911) were used to build the various cross-section stress distributions in Fig. 2.b.
(09) σ curv ( y ) = { [ S y 0 + E t ( y ρ i ε y ) ] y y < y c ( y y y ) S y 0 y y y y y [ S y 0 + E t ( y ρ i ε y ) ] c y < y y (10) σ s b ( y ) = M ( ρ i ) y I , c y c where I = 2 c c y 2 ( radius 3 2 y 2 b 1 2 ) d y + b 1 b 3 12 (11) σ res _ c ( y ) = σ curv ( y ) + σ s b ( y ) , c y c (09) σ curv  ( y ) = S y 0 + E t y ρ i ε y y y < y c y y y S y 0 y y y y y S y 0 + E t y ρ i ε y c y < y y (10) σ s b ( y ) = M ρ i y I , c y c  where  I = 2 c c y 2  radius  3 2 y 2 b 1 2 d y + b 1 b 3 12 (11) σ res  _ c ( y ) = σ curv  ( y ) + σ s b ( y ) , c y c {:[(09)sigma_("curv ")(y)={[-[S_(y0)+E_(t)*((y)/(rho_(i))-epsi_(y))],y_(y) < y <= c],[-((y)/(y_(y)))*S_(y0),-y_(y) <= y <= y_(y)],[[S_(y0)+E_(t)*((-y)/(rho_(i))-epsi_(y))],-c <= y < -y_(y)]:}],[(10)sigma_(sb)(y)=(M(rho_(i))*y)/(I)","-c <= y <= c quad" where "I=2*int_(-c)^(c)y^(2)*(sqrt(" radius "_(3)^(2)-y^(2))-(b_(1))/(2))*dy+(b_(1)*b^(3))/(12)],[(11)sigma_("res "_c)^(**)(y)=sigma_("curv ")(y)+sigma_(sb)(y)","quad-c <= y <= c]:}\begin{align*} & \sigma_{\text {curv }}(y)= \begin{cases}-\left[S_{y 0}+E_{t} \cdot\left(\frac{y}{\rho_{i}}-\varepsilon_{y}\right)\right] & y_{y}<y \leq c \\ -\left(\frac{y}{y_{y}}\right) \cdot S_{y 0} & -y_{y} \leq y \leq y_{y} \\ {\left[S_{y 0}+E_{t} \cdot\left(\frac{-y}{\rho_{i}}-\varepsilon_{y}\right)\right]} & -c \leq y<-y_{y}\end{cases} \tag{09}\\ & \sigma_{s b}(y)=\frac{M\left(\rho_{i}\right) \cdot y}{I},-c \leq y \leq c \quad \text { where } I=2 \cdot \int_{-c}^{c} y^{2} \cdot\left(\sqrt{\text { radius }_{3}^{2}-y^{2}}-\frac{b_{1}}{2}\right) \cdot d y+\frac{b_{1} \cdot b^{3}}{12} \tag{10}\\ & \sigma_{\text {res } \_c}^{*}(y)=\sigma_{\text {curv }}(y)+\sigma_{s b}(y), \quad-c \leq y \leq c \tag{11} \end{align*}(09)σcurv (y)={[Sy0+Et(yρiεy)]yy<yc(yyy)Sy0yyyyy[Sy0+Et(yρiεy)]cy<yy(10)σsb(y)=M(ρi)yI,cyc where I=2ccy2( radius 32y2b12)dy+b1b312(11)σres _c(y)=σcurv (y)+σsb(y),cyc
where, σ curv ( y ) , σ s b ( y ) σ curv  ( y ) , σ s b ( y ) sigma_("curv ")(y),sigma_(sb)(y)\sigma_{\text {curv }}(y), \sigma_{s b}(y)σcurv (y),σsb(y) and σ res_c ( y ) σ res_c  ( y ) sigma_("res_c ")^(**)(y)\sigma_{\text {res_c }}^{*}(y)σres_c (y) are, respectively, loading stress distribution, the spring-back stress distribution, and the residual stress distribution, all within phase 1 . I 1 . I 1.I1 . I1.I is the area moment of inertia for a straight specimen.
Fig. 2.a shows, the σ ( ε ) / S y 0 σ ( ε ) / S y 0 sigma(epsi)//S_(y0)\sigma(\varepsilon) / \mathrm{S}_{\mathrm{y} 0}σ(ε)/Sy0 vs ε ε epsi\varepsilonε evolution diagram of specimen points a a aaa and b b bbb during phase 1, and Fig. 2.b shows the specimen cross-section parameterized stress distribution in phase 1 1 1-1-1 loading ( σ curv ( y ) / S y 0 ) σ curv  ( y ) / S y 0 (sigma_("curv ")(y)//S_(y0))\left(\sigma_{\text {curv }}(\mathrm{y}) / \mathrm{S}_{\mathrm{y} 0}\right)(σcurv (y)/Sy0), spring back ( σ sb ( y ) / S y 0 ) σ sb ( y ) / S y 0 (sigma_(sb)(y)//S_(y0))\left(\sigma_{\mathrm{sb}}(\mathrm{y}) / \mathrm{S}_{\mathrm{y} 0}\right)(σsb(y)/Sy0), and residual ( σ res_c ( y ) / S y 0 ) σ res_c  ( y ) / S y 0 (sigma_("res_c ")^(**)(y)//S_(y0))\left(\sigma_{\text {res_c }}^{*}(y) / \mathrm{S}_{\mathrm{y} 0}\right)(σres_c (y)/Sy0).
Figure 2: (a) Specimen phase 1 graph ( σ ( ε ) / S y 0 σ ( ε ) / S y 0 sigma(epsi)//S_(y0)\sigma(\varepsilon) / \mathrm{S}_{\mathrm{y} 0}σ(ε)/Sy0 vs ε ε epsi\varepsilonε ) for points a a aaa and b b bbb, and (b) Specimen cross-section parameterized stress distribution in phase 1 1 1-1-1 loading ( σ curv ( y ) / S y 0 ) σ curv  ( y ) / S y 0 (sigma_("curv ")(y)//S_(y0))\left(\sigma_{\text {curv }}(\mathrm{y}) / \mathrm{S}_{\mathrm{y} 0}\right)(σcurv (y)/Sy0), spring back ( σ sb ( y ) / S y 0 ) σ sb ( y ) / S y 0 (sigma_(sb)(y)//S_(y0))\left(\sigma_{\mathrm{sb}}(\mathrm{y}) / \mathrm{S}_{\mathrm{y} 0}\right)(σsb(y)/Sy0), and residual ( σ res_c ( y ) / S y 0 ) σ res_c  ( y ) / S y 0 (sigma_("res_c ")^(**)(y)//S_(y0))\left(\sigma_{\text {res_c }}^{*}(y) / \mathrm{S}_{\mathrm{y} 0}\right)(σres_c (y)/Sy0).
Note that Fig.2.b axes are dimensionless to improve conceptual understanding. For instance, when σ curv ( y ) / S y 0 = 1 σ curv  ( y ) / S y 0 = 1 sigma_("curv ")(y)//S_(y0)=1\sigma_{\text {curv }}(\mathrm{y}) / \mathrm{S}_{\mathrm{y} 0}=1σcurv (y)/Sy0=1, the material begins to yield. Also, an interesting point for mechanical parts severely yielded is that the residual stresses at the elastic-plastic border ( y / c = ± y y / c y / c = ± y y / c y//c=+-y_(y)//c\mathrm{y} / \mathrm{c}= \pm \mathrm{y}_{\mathrm{y}} / \mathrm{c}y/c=±yy/c ) are bigger than the residual stresses at the external surfaces ( y / c = ± 1 y / c = ± 1 y//c=+-1\mathrm{y} / \mathrm{c}= \pm 1y/c=±1 ).

Combining the residual stress and the fatigue approaches

At the end of phase 1, the specimen cross-section stress distribution σ res _ c ( y ) σ res  _ c ( y ) sigma_("res "_c)^(**)(y)\sigma_{\text {res } \_c}^{*}(y)σres _c(y) is shown in Fig. 2.b. Before beginning phase 2 the specimen is turned upside down, obeying Eqn. (12) to convert σ res _ c ( y ) σ res  _ c ( y ) sigma_("res "_c)^(**)(y)\sigma_{\text {res } \_c}^{*}(y)σres _c(y) of Fig. 2.b into σ res _ ( y ) σ res  _ ( y ) sigma_("res "_)(y)\sigma_{\text {res } \_}(y)σres _(y) of Fig. 3.
(12) σ r e s c ( y ) = σ r e s c ( y ) , c y c (12) σ r e s c ( y ) = σ r e s c ( y ) , c y c {:(12)sigma_(res_(-)c)(y)=-sigma_(res_(-)c)^(**)(y)","quad-c <= y <= c:}\begin{equation*} \sigma_{r e s_{-} c}(y)=-\sigma_{r e s_{-} c}^{*}(y), \quad-c \leq y \leq c \tag{12} \end{equation*}(12)σresc(y)=σresc(y),cyc
Combining the residual stress and the fatigue stress loading distributions, using Eqns. (12) and Eqns. (4A.a) and (5A.a), results in Eqns. (13) and (14), as shown, graphically, in Fig. 3.
(14) σ f max ( y ) = σ r e s c ( y ) + σ fad_ max ( y ) , c y c σ f min ( y ) = σ r e s c ( y ) + σ fad_ min ( y ) , c y c (14) σ f max ( y ) = σ r e s c ( y ) + σ fad_  max ( y ) , c y c σ f min ( y ) = σ r e s c ( y ) + σ fad_  min ( y ) , c y c {:(14){:[sigma_(f_(-)max)(y)=sigma_(res_(-)c)(y)+sigma_("fad_ "max)(y)",",-c <= y <= c],[sigma_(f_(-)min)(y)=sigma_(res_(-)c)(y)+sigma_("fad_ "min)(y)",",-c <= y <= c]:}:}\begin{array}{ll} \sigma_{f_{-} \max }(y)=\sigma_{r e s_{-} c}(y)+\sigma_{\text {fad_ } \max }(y), & -c \leq y \leq c \\ \sigma_{f_{-} \min }(y)=\sigma_{r e s_{-} c}(y)+\sigma_{\text {fad_ } \min }(y), & -c \leq y \leq c \tag{14} \end{array}(14)σfmax(y)=σresc(y)+σfad_ max(y),cycσfmin(y)=σresc(y)+σfad_ min(y),cyc
Fig. 3 shows the specimen cross-section stress distribution resulting from Eqns. (4А.a, 5A.a, 12, 13, 14). Both, σ fad_ max ( y ) / S y 0 σ fad_ max  ( y ) / S y 0 sigma_("fad_ max ")(y)//S_(y0)\sigma_{\text {fad_ max }}(y) / S_{y 0}σfad_ max (y)/Sy0 and σ fad_ min ( y ) / S y 0 σ fad_ min  ( y ) / S y 0 sigma_("fad_ min ")(y)//S_(y0)\sigma_{\text {fad_ min }}(y) / S_{y 0}σfad_ min (y)/Sy0 of Fig. 3.a and 3.b, respectively, have elastic stress distribution, but for the σ f max ( y ) / S y 0 σ f max ( y ) / S y 0 sigma_(f_(-)max)(y)//S_(y0)\sigma_{f_{-} \max }(y) / S_{y 0}σfmax(y)/Sy0 distribution, shown in Fig.3.a, it seems to have some yielding at external specimen surfaces. The main reason for applying such severe combined loading is to impose short lives on the specimens, intending to finish the experimental part of this research in a reasonable period.
Figure 3: Phase 2 - Specimen cross-section loading stress distribution: (a) maximum and (b) minimum.

Experimental setup

The experimental part of this research was done in LAMAT Laboratory (https://dippg.cefet-rj.br/ppemm/index.php/pt/laboratorios/57-estrutura/296-laboratorio-de-materiais-lamat) of the Federal Center for Technological Education Celso Suckow da Fonseca - CEFET/RJ (https://www.cefet-rj.br/), in Rio de Janeiro, Brazil. It was used a servo-hydraulic material testing machine INSTRON 8801 of 100 kN of capacity, with a three-point bending grip setup with a phase 2 specimen positioned, as shown in Fig.4.a. The frequency of the fatigue test was limited to 10 Hz , due to the difficulties in adjusting the PID feedback control for such flexible specimens. The fatigue test was made in load control, and an excessive displacement interlock was set to stop the test immediately after the specimen had failed. Fig. 4.b shows an example of the experimental graph output displacement vs load of a specimen in phase 2. Fig. 4.c shows the two specimens' groups ( a a a\boldsymbol{a}a and b b b\boldsymbol{b}b ) after phase 1 , where the curvatures were imposed, before being used in phase 2 of the fatigue test.
Each specimen has 25 cm length, and cross-section nominal values of 14 mm × 6 mm 14 mm × 6 mm 14mmxx6mm14 \mathrm{~mm} \times 6 \mathrm{~mm}14 mm×6 mm, which were taken from the manufacturer's wire spools of cold-drawn steel wires. Figs. 4.d and 4.e show examples of specimens' final fractures, from an external and an internal point of view.
In Fig. 4.b the abscissa axis shows the displacement of the piston rod, and the ordinate axis shows the load cell measurements during the fatigue test. The test began with the cycles positioned at the right part of the graphic. As the test heads towards the end, the distance between cycles increases until the specimen final failure, when the excessive displacement activates the interlock, ending the test. In the macrography of Fig. 4.e, it is possible to see a fatigue crack propagation region (the smooth region at the cross-section superior area) and a final failure region (the rough surface at the cross-section central area).

RESULTS AND DISCUSSION

I n n nnn the experimental part of this research, an unexpected behavior of the tested specimens occurred. The specimens submitted to phases 1 and 2 had surprisingly long lives. The experimental results are summarized from the experimental approach of the original research developed in Zaidan [17], which is not entirely available because of a confidentiality agreement. Nevertheless, the results presented in the text are sufficient to support the development of the proposed analytical model.
It was initially assumed that the specimens taken from the manufacturer's wire spools of cold-drawn steel wires, did not have appreciable residual stress. As the specimens were taken from the manufacturer's wire spools, they were marked on the upper side with a yellow marker. The specimens were divided into groups a a a\boldsymbol{a}a, and b b b\boldsymbol{b}b. These groups were related to how the specimens were positioned on the 3 -point bending gripper, with the wire's original yellow marking above or below. The objective was to verify if the specimen side, with the yellow marking top side or the bottom side, affects the specimen's
fatigue performance. Fig. 5 shows the theoretical SN curve ( SN teo 0 % SN teo  0 % SN_("teo ")0%\mathrm{SN}_{\text {teo }} 0 \%SNteo 0% ), with experimental results of groups a a a\boldsymbol{a}a, and b b b\boldsymbol{b}b, each group with three specimens, and the notable point O F O F OFO FOF.
Figure 4: Experimental setup - (a) Servo-hydraulic material testing machine INSTRON of 100 kN of capacity, (b) Experimental graph output displacement vs load of a specimen in phase 2 (c) Groups a a a\boldsymbol{a}a and b b b\boldsymbol{b}b specimens after phase 1 . Specimen failed after phase 2 2 2-2-2 (d) external view and (e) internal view (macrography).
The fatigue results for both specimen groups a a a\boldsymbol{a}a and b b b\boldsymbol{b}b eventually show that they were quite similar. So, it was taken the average of the number of cycles generated by the specimens of both groups, N exp N exp  N_("exp ")N_{\text {exp }}Nexp , shown in Fig. 5 . Note that the number of cycles that the specimens failed N exp N exp  N_("exp ")N_{\text {exp }}Nexp  exceeded the number of cycles estimated by the proposed fatigue model N teo N teo  N_("teo ")N_{\text {teo }}Nteo . With the aid of Eqns. (7A.a) and (7A.b), the equivalent alternate stress σ a e q σ a e q sigma_(a_(-)eq)\sigma_{a_{-} e q}σaeq and the number of cycles until fatigue failure, Nteo are estimated, which correspond to point OF of Fig. 5. In effect, the point of coordinates ( N exp , σ a e q / S y 0 N exp  , σ a e q / S y 0 N_("exp "),sigma_(a_(-)eq)//S_(y0)N_{\text {exp }}, \sigma_{a_{-} e q} / S_{y 0}Nexp ,σaeq/Sy0 ) is not on the theoretical SN curve ( SN teo 0 % SN teo  0 % SN_("teo ")0%\mathrm{SN}_{\text {teo }} 0 \%SNteo 0% ), described by Eqn. (2A). An explanation for this behavior is proposed next.
Figure 5: Theoretical SN curve ( SN teo 0 % SN teo  0 % SN_("teo ")0%\mathrm{SN}_{\text {teo }} 0 \%SNteo 0% ), with experimental results.
It was tentatively solved by adjusting the S ut 0 S ut 0 S_(ut0)\mathrm{S}_{\mathrm{ut} 0}Sut0 for a higher value, considering | 6 % | | 6 % | |6%||6 \%||6%| residual strain present in specimen surfaces after phase 1, as shown in point a a aaa of Fig. 2.a. This S ut S ut S_(ut)\mathrm{S}_{\mathrm{ut}}Sut correction generated a new theoretical SN curve (SN SN teo 6 % SN teo  6 % SN_("teo ")6%\mathrm{SN}_{\text {teo }} 6 \%SNteo 6% ) of Eqn. (2A) approximating to the experimental results, but not sufficiently. Moving the SN curve more in the "right" direction, a new SN curve was proposed, named the hypothetical SN curve ( SN hyp 6 % SN hyp  6 % SN_("hyp ")6%\mathrm{SN}_{\text {hyp }} 6 \%SNhyp 6% ), which passes through the fatigue experimental points. To describe the SN hyp 6 % SN hyp  6 % SN_("hyp ")6%\mathrm{SN}_{\text {hyp }} 6 \%SNhyp 6%, shown in Fig. 6, Eqn. (15) is proposed
(15) S f b y p = { S u t ( N ) ( 1 3 log ( f ) ) , 10 0 N 10 3 a h y p ( N ) b h p p , 10 3 < N < 10 6 a h y p ( 10 6 ) b h p p , N 10 6 a h y p = σ a e q _ 6 ( N e x p ) b h p p b h y p = b t e o (15) S f b y p = S u t ( N ) 1 3 log ( f ) , 10 0 N 10 3 a h y p ( N ) b h p p , 10 3 < N < 10 6 a h y p 10 6 b h p p , N 10 6 a h y p = σ a e q _ 6 N e x p b h p p b h y p = b t e o {:(15)S_(f)^(byp)={[S_(ut)*(N)^(((1)/(3)*log(f)))",",10^(0) <= N <= 10^(3)],[a_(hyp)*(N)^(b_(hpp))",",10^(3) < N < 10^(6)],[a_(hyp)*(10^(6))^(b_(hpp))",",N >= 10^(6)]quada_(hyp)=(sigma_(a_(-)eq_6))/((N_(exp))^(b_(hpp)))quadb_(hyp)=b_(teo):}:}S_{f}^{b y p}=\left\{\begin{array}{lr} S_{u t} \cdot(N)^{\left(\frac{1}{3} \cdot \log (f)\right)}, & 10^{0} \leq N \leq 10^{3} \tag{15}\\ a_{h y p} \cdot(N)^{b_{h p p}}, & 10^{3}<N<10^{6} \\ a_{h y p} \cdot\left(10^{6}\right)^{b_{h p p}}, & N \geq 10^{6} \end{array} \quad a_{h y p}=\frac{\sigma_{a_{-} e q \_6}}{\left(N_{e x p}\right)^{b_{h p p}}} \quad b_{h y p}=b_{t e o}\right.(15)Sfbyp={Sut(N)(13log(f)),100N103ahyp(N)bhpp,103<N<106ahyp(106)bhpp,N106ahyp=σaeq_6(Nexp)bhppbhyp=bteo
where S f b y p S f b y p S_(f)^(byp)S_{f}^{b y p}Sfbyp is the hypothetical fatigue resistance, and a hyp a hyp  a_("hyp ")\mathrm{a}_{\text {hyp }}ahyp  and b hyp b hyp  b_("hyp ")\mathrm{b}_{\text {hyp }}bhyp  are hypothetical constants. Note that equivalent alternate stress σ a e q _ 6 σ a e q _ 6 sigma_(a_(-)eq_6)\sigma_{a_{-} e q \_6}σaeq_6 of Eqn. (15.b) uses the same Eqn. (7A.a), but with two differences: a) the material has a 6 % 6 % 6%6 \%6% residual strain, instead of 0 % 0 % 0%0 \%0%, and b ) the loading combines residual and fatigue stresses, instead of only fatigue stresses.
Considering that the average experimental number of cycles N exp N exp  N_("exp ")\mathrm{N}_{\text {exp }}Nexp  is the only undeniable experimental result, two different approaches were proposed to explain the difference between N teo N teo  N_("teo ")\mathrm{N}_{\text {teo }}Nteo  and N exp N exp  N_("exp ")\mathrm{N}_{\text {exp }}Nexp . The two notable points related to the two different proposed approaches are shown in Fig. 6.
The first approach requires the existence of the new SN curve, called SN SNy hyp 6 % SNy hyp  6 % SNy_("hyp ")6%\mathrm{SNy}_{\text {hyp }} 6 \%SNyhyp 6%, which is described by Eqn. (15). This approach receives the name of "resistance increase". One explanation for this proposition can be reached through the modification of the miscellaneous-effects factor k f k f k_(f)\mathrm{k}_{\mathrm{f}}kf of Eqn. (1A), substituting the value 1 by a value greater than one, rising the S e S e S_(e)\mathrm{S}_{\mathrm{e}}Se value, and therefore generating the a hyp a hyp  a_("hyp ")a_{\text {hyp }}ahyp  and b hyp b hyp  b_("hyp ")b_{\text {hyp }}bhyp  values of SN hyp 6 % SN hyp  6 % SN_("hyp ")6%\mathrm{SN}_{\text {hyp }} 6 \%SNhyp 6% curve. The point ( N exp , σ a eq R I / S y 0 N exp  , σ a eq  R I / S y 0 N_("exp "),sigma_(a_(-)"eq ")^(RI)//S_(y0)N_{\text {exp }}, \sigma_{a_{-} \text {eq }}^{R I} / S_{y 0}Nexp ,σaeq RI/Sy0 ) is represented by a black triangle RI on the SN hyp 6 % SN hyp  6 % SN_("hyp ")6%\mathrm{SN}_{\text {hyp }} 6 \%SNhyp 6% curve in Fig. 6:
(16) σ a e q R I = a h y p ( N e x p ) b b p p (16) σ a e q R I = a h y p N e x p b b p p {:(16)sigma_(a_(-)eq)^(RI)=a_(hyp)*(N_(exp))^(b_(bpp)):}\begin{equation*} \sigma_{a_{-} e q}^{R I}=a_{h y p} \cdot\left(N_{e x p}\right)^{b_{b p p}} \tag{16} \end{equation*}(16)σaeqRI=ahyp(Nexp)bbpp
The second approach maintains the original SN curve, called SN SN teo 0 % SN teo  0 % SN_("teo ")0%\mathrm{SN}_{\text {teo }} 0 \%SNteo 0%, but reduces σ a e q σ a e q sigma_(a_(-)eq)\sigma_{a_{-} e q}σaeq originally estimated by Eqn. (7A.a) to σ a e q S R σ a e q S R sigma_(a_(-)eq)^(SR)\sigma_{a_{-} e q}^{S R}σaeqSR of Eqn. (17). This approach receives the name "stress reduction". This effect can appear, for instance, by the existence of superficial residual compressive stress in function of the manufacturing process, such as cold rolling. This approach doesn't modify the fatigue limit as the former approach but actuates only diminishing the effective tensile stresses that is acting in the specimen's cross-section area. See Gothivarekar et al. [18] for a schematic cross-section residual stress distribution representation after a cold rolling process. The point ( N exp , σ a e q S R / S y 0 N exp  , σ a e q S R / S y 0 N_("exp "),sigma_(a_(-)eq)^(SR)//S_(y0)N_{\text {exp }}, \sigma_{a_{-} e q}^{S R} / S_{y 0}Nexp ,σaeqSR/Sy0 ) is represented by a black inverted triangle SR on the SN SN teo 0 % SN teo  0 % SN_("teo ")0%\mathrm{SN}_{\text {teo }} 0 \%SNteo 0% curve in Fig. 6:
(17) σ a e q S R = a t e o ( N e x p ) b t e o (17) σ a e q S R = a t e o N e x p b t e o {:(17)sigma_(a_(-)eq)^(SR)=a_(teo)*(N_(exp))^(b_(teo)):}\begin{equation*} \sigma_{a_{-} e q}^{S R}=a_{t e o} \cdot\left(N_{e x p}\right)^{b_{t e o}} \tag{17} \end{equation*}(17)σaeqSR=ateo(Nexp)bteo
From the mechanical engineering point of view, this second approach, named "stress reduction", is more reasonable. Fig. 6 resumes, in a graphical form, the proposed SN curves and the notable points.
Figure 6: Proposed SN curves and the notable points.
The RI and SR stresses are calculated, respectively, through Eqns. (16) and (17), positioned in different SN curves, respectively, on SN hyp 6 % SN hyp  6 % SN_("hyp ")6%\mathrm{SN}_{\text {hyp }} 6 \%SNhyp 6% and on SN teo 0 % SN teo  0 % SN_("teo ")0%\mathrm{SN}_{\text {teo }} 0 \%SNteo 0%, but have the same number of cycles, N exp N exp  N_("exp ")\mathrm{N}_{\text {exp }}Nexp . Comparing OF point (only fatigue) of Fig 5, with the SR or RI points (both for fatigue plus residual stress) of Fig. 6, there is a perceptible difference in cycles to failure between N teo N teo  N_("teo ")\mathrm{N}_{\text {teo }}Nteo  and N exp N exp  N_("exp ")\mathrm{N}_{\text {exp }}Nexp , of almost an order of magnitude. The compressive superficial residual stress generated by the wire lamination could extend the specimen's life, even though the specimen has been submitted, in phase 1, to severe tensile residual stress. In summary, it is proposed that the lowering σ a e q R I σ a e q R I sigma_(a_(-)eq)^(RI)\sigma_{a_{-} e q}^{R I}σaeqRI to a new equivalent alternate stress σ a e q S R σ a e q S R sigma_(a_(-)eq)^(SR)\sigma_{a_{-} e q}^{S R}σaeqSR is feasible by the existence of superficial compressive residual stress caused by the lamination process, which can be estimated as:
(18) σ a e q l a m = σ a e q S R σ a e q R I (18) σ a e q l a m = σ a e q S R σ a e q R I {:(18)sigma_(a_(-)eq)^(lam)=sigma_(a_(-)eq)^(SR)-sigma_(a_(-)eq)^(RI):}\begin{equation*} \sigma_{a_{-} e q}^{l a m}=\sigma_{a_{-} e q}^{S R}-\sigma_{a_{-} e q}^{R I} \tag{18} \end{equation*}(18)σaeqlam=σaeqSRσaeqRI
Besides the superficial compressive residual stress caused by the lamination process, estimated by Eqn. (18), it is possible to propose an estimation of the effect on the cross-section residual stress distribution resulting from the specimen cold rolling process:
(19) σ r e s l ( y ) = { σ a e q l a m y l y c ( c y l y l ) σ a e q l a m y l y y l where, y l = 0.8 c σ a e q l a m y l y c (19) σ r e s l ( y ) = σ a e q l a m y l y c c y l y l σ a e q l a m y l y y l  where,  y l = 0.8 c σ a e q l a m y l y c {:(19)sigma_(res_(l))(y)={[sigma_(a_(-)eq)^(lam),y_(l) <= y <= c],[-((c-y_(l))/(y_(l)))*sigma_(a_(eq))^(lam),y_(l) <= y <= -y_(l)quad" where, "y_(l)=0.8*c],[sigma_(a_(-)eq)^(lam),-y_(l) <= y <= -c]:}:}\sigma_{r e s_{l}}(y)= \begin{cases}\sigma_{a_{-} e q}^{l a m} & \mathrm{y}_{l} \leq y \leq c \tag{19}\\ -\left(\frac{c-\mathrm{y}_{l}}{\mathrm{y}_{l}}\right) \cdot \sigma_{a_{e q}}^{l a m} & \mathrm{y}_{l} \leq y \leq-\mathrm{y}_{l} \quad \text { where, } \mathrm{y}_{l}=0.8 \cdot \mathrm{c} \\ \sigma_{a_{-} e q}^{l a m} & -\mathrm{y}_{l} \leq y \leq-c\end{cases}(19)σresl(y)={σaeqlamylyc(cylyl)σaeqlamylyyl where, yl=0.8cσaeqlamylyc
The cross-section semi-height, for tensile stress is yl and for compressive stress is ( c yl c yl c-yl\mathrm{c}-\mathrm{yl}cyl ). Its value was chosen arbitrarily since this information is relatively scarce in technical literature, usually presented only qualitatively, as in Gothivarekar et al. [18].
Adding Eqn. (19) to Eqns. (13) and (14), generates the maximum and minimum specimen's cross-section stress distribution, respectively, Eqns. (20) and (21):
(21) σ f max ( y ) = σ f max ( y ) + σ r s l ( y ) c y c σ f min ( y ) = σ f min ( y ) + σ r e l ( y ) c y c (21) σ f max ( y ) = σ f max ( y ) + σ r s l ( y ) c y c σ f min ( y ) = σ f min ( y ) + σ r e l ( y ) c y c {:(21){:[sigma_(f_(-)max)^(**)(y)=sigma_(f_(-)max)(y)+sigma_(rs_(l))(y),-c <= y <= c],[sigma_(f_(-)min)^(**)(y)=sigma_(f_(-)min)(y)+sigma_(re_(l))(y),-c <= y <= c]:}:}\begin{array}{ll} \sigma_{f_{-} \max }^{*}(y)=\sigma_{f_{-} \max }(y)+\sigma_{r s_{l}}(y) & -c \leq y \leq c \\ \sigma_{f_{-} \min }^{*}(y)=\sigma_{f_{-} \min }(y)+\sigma_{r e_{l}}(y) & -c \leq y \leq c \tag{21} \end{array}(21)σfmax(y)=σfmax(y)+σrsl(y)cycσfmin(y)=σfmin(y)+σrel(y)cyc
Finally, adding all the residual stress (from the curvature in phase 1 and from the cold rolling process):
(22) σ r e s ( y ) = σ r e s c ( y ) + σ r e s l ( y ) c y c (22) σ r e s ( y ) = σ r e s c ( y ) + σ r e s l ( y ) c y c {:(22)sigma_(res)(y)=sigma_(res_(-)c)(y)+sigma_(res_(l))(y)quad-c <= y <= c:}\begin{equation*} \sigma_{r e s}(y)=\sigma_{r e s_{-} c}(y)+\sigma_{r e s_{l}}(y) \quad-c \leq y \leq c \tag{22} \end{equation*}(22)σres(y)=σresc(y)+σresl(y)cyc
Fig. 7.a and 7.b exhibit the graphical representation, respectively, for Eqns. (20) and (21), and for Eqn. (22).
Figure 7: Specimen cross-section stress distribution: (a) in Phase 2 and (b) after Phase 2.
Fig. 7.a shows the effect of the lamination in the specimen's cross-section loading stress distribution during phase 2. Fig. 7.b shows the sum of all residual stresses, that exist after phase 2. Comparing Fig. 7.a with Fig. 3.a, there is a reduction of superficial tensile stress at extrados ( y = c y = c y=-cy=-cy=c ) and increasing the compressive stress at intrados ( y = + c y = + c y=+cy=+cy=+c ). Also, note in Fig. 7.b, that the residual stresses at the elastic-plastic border ( y / c = ± y y / c y / c = ± y y / c y//c=+-y_(y)//c\mathrm{y} / \mathrm{c}= \pm \mathrm{y}_{\mathrm{y}} / \mathrm{c}y/c=±yy/c ) are bigger than at the external surfaces ( y / c = ± 1 y / c = ± 1 y//c=+-1\mathrm{y} / \mathrm{c}= \pm 1y/c=±1 ).

Conclusions

An analytical model was built in detail, combining fatigue and residual stresses, to check the effect of residual stress on the fatigue performance of partially yielded cold-drawn steel wire specimens. All equations and their graphical representations were explicitly presented. Initially, it was assumed that the specimens, taken from the manufacturer's wire spools, did not have appreciable residual stress. So, a curvature was imposed on specimens to generate
a specified distribution of residual stresses. In sequence, the specimens were subjected to time-variable experimental bending loading until fatigue failure occurred. It was realized that the number of cycles the specimens spent up to failure was nearly one order of magnitude larger than those theoretically predicted. It was recognized that the residual stresses imposed by curvature on specimens were insufficient to remove the former beneficial compressive residual stresses, probably generated by cold rolling manufacturing. Two approaches were designated, with success, to adjust the analytical model to encompass the number of cycles difference. An analytical estimation of the level of the former superficial compressive residual stresses was also achieved. These results can contribute to improving the knowledge about the effectiveness of residual stress in the fatigue performance of partially yielded cold-drawn steel wires, including the initially unexpected beneficial effects promoted by the cold rolling fabrication process.

References

[1] Laterza, M., D'amato, M. and Casamassima, V.M. (2017). Stress-Life Curves Method for Fatigue Assessment of Ancient Brick Arch Bridges. International Journal of Arch. Heritage. 11(6), pp. 843-858.
DOI: https://doi.org/10.1080/15583058.2017.1315621.
[2] Teoh, S.H. (2000). Fatigue of biomaterials: a review. International Journal of Fatigue, 22, pp. 825-837.
[3] Tabatabaeian, A., Ghasemi, A.R., Shokrieh, M.M., Marzbanrad, B., Baraheni, M., Fotouhi, M. (2022). Advanced Engineering Materials; 24(3), pp. 2100786. DOI: https://doi.org/10.1002/adem. 202100786.
[4] Crandall, S.H., Dahl, N.C. and Lardner, T.J. (1978). An Introduction to the Mechanics of Solids. Second Edition with SI units, McGraw Hill International Editions.
[5] Stok, B. and Halilovic M. (2008). Analytical solutions in elastoplastic bending of beams with rectangular cross-section. Applied Mathematical Modelling; 33, pp. 1749-1760. DOI: https://doi.org/10.1016/j.apm.2008.03.011.
[6] Rimovskis, S. and Sabaliauskas, A. (2012). Analysis of Rectangular and Circular Cross-section Power Hardening Elements under Pure Bending. International Journal of Materials Engineering, 2(6), pp.84-89.
DOI: https://doi.org/10.5923/j.ijme.20120206.03.
[7] Castro, J.T.P. and Meggiolaro, M.A. (2016). Fatigue Design Techniques, Volume 2: Low-Cycle and Multiaxial Fatigue, 1. ed., Scotts Valley, CA 95066, USA: CreateSpace.
[8] Jirásek, M. and Bazant, Z.P. (2002). Inelastic Analysis of Structures. John Wiley and Sons Ltd.
[9] De Castro, F.A., Kenedi, P.P., Vignolli, L.L. and Riagusoff, I.I.T. (2021) Residual stress distribution in built-in beams. Proceedings of the Institution of Mechanical Engineers Part L-Journal of Materials-Design and Applications. 235(1), pp. 216-231. DOI: https://doi.org/10.1177/1464420720958617.
[10] Schajer, G.S. (2013). Practical Residual Stress Measurement Methods. Wiley.
[11] Atienza, I., and Elices, M. (2003). Influence of residual stresses in the tensile test of cold-drawn wires. Materials and Structures, 36, pp. 548-552.
[12] McClung, R. C. (2007). A literature survey on the stability and significance of residual stresses during fatigue. Fatigue & Fracture of Engineering Materials and Structures, 30(3), pp.173-205.
DOI: https://doi.org/10.1111/j.1460-2695.2007.01102.x.
[13] Webster, G.A. and Ezielo, A.N. (2001). Residual stress distributions and their influence on fatigue lifetimes. International Journal of Fatigue. 23 (1), pp. 375-383. DOI: https://doi.org/10.1016/S0142-1123(01)00133-5.
[14]James, M.N., Hughes, D.J., Chen, Z., Lombard, H., Hattingh, D.G., Asquith, D., Yates, J.R., Webster, P.J. (2007) Residual stresses and fatigue performance. Engineering Failure Analysis. 14(2), pp. 384-395.
DOI: https://doi.org/10.1016/j.engfailanal.2006.02.011.
[15] Vaara, J., Kunnari, A. and Frondelius, T. (2020). Engineering Failure Analysis. 110, pp. 104379. DOI: https://doi.org/10.1016/j.engfailanal.2020.104379.
[16] Zaidan, L.D., Buarque, Y. C., Kenedi, P.P., de Souza, L. F. G. (2022). Fatigue and residual stresses - a deleterious combination. In: Proceedings of the International Symposium on Solid Mechanics, MECSOL 2022, Campinas - SP, Brazil.
[17] Zaidan, L.D. (2023). Influência da Tensão Residual na vida à Fadiga de Estruturas Parcialmente Plastificadas, Master Thesis (in Portuguese), Centro Federal de Educação Tecnológica Celso Suckow da Fonseca - CEFET/RJ, Brazil.
[18] Gothivarekar, S., Coppieters, S., Van de Velde, A., Debruyne, D. (2020). Advanced FE model validation of cold-forming process using DIC: Air bending of high-strength steel. Int J Mater Form; 13, pp. 409-421.
DOI: https://doi.org/10.1007/s12289-020-01536.
[19] Budynas, R.G. and Nisbett, J.K. (2015). Shigley's Mechanical Engineering Design. McGraw-Hill, 10th Edition.

Appendix

Fatigue basics


he fatigue part of the analytical model is explicitly presented and commented on, based on the classical approach, as in Budynas et al. [19]. The fatigue limit S e S e S_(e)S_{e}Se is estimated as:
(1A) S e = k a k b k c k d k e k f 0.5 S u t = ( 4.51 ( S u t ) 0.265 ) [ 0.808 b 2 h 7.62 ] 0.107 1 1 1.1 0.5 S u t (1A) S e = k a k b k c k d k e k f 0.5 S u t = 4.51 S u t 0.265 0.808 b 2 h 7.62 0.107 1 1 1.1 0.5 S u t {:(1A)S_(e)=k_(a)*k_(b)*k_(c)*k_(d)*k_(e)*k_(f)*0.5*S_(ut)=(4.51*(S_(ut))^(-0.265))*[(0.808*sqrt(b_(2)*h))/(7.62)]^(-0.107)*1*1*1.1*0.5*S_(ut):}\begin{equation*} S_{e}=k_{a} \cdot k_{b} \cdot k_{c} \cdot k_{d} \cdot k_{e} \cdot k_{f} \cdot 0.5 \cdot S_{u t}=\left(4.51 \cdot\left(S_{u t}\right)^{-0.265}\right) \cdot\left[\frac{0.808 \cdot \sqrt{b_{2} \cdot h}}{7.62}\right]^{-0.107} \cdot 1 \cdot 1 \cdot 1.1 \cdot 0.5 \cdot S_{u t} \tag{1A} \end{equation*}(1A)Se=kakbkckdkekf0.5Sut=(4.51(Sut)0.265)[0.808b2h7.62]0.107111.10.5Sut
where, k a , k b , k c , k d , k e k a , k b , k c , k d , k e k_(a),k_(b),k_(c),k_(d),k_(e)\mathrm{k}_{\mathrm{a}}, \mathrm{k}_{\mathrm{b}}, \mathrm{k}_{\mathrm{c}}, \mathrm{k}_{\mathrm{d}}, \mathrm{k}_{\mathrm{e}}ka,kb,kc,kd,ke and k f k f k_(f)\mathrm{k}_{\mathrm{f}}kf are the modifying factors. b 2 b 2 b_(2)b_{2}b2 and h h hhh are used in mm and S ut S ut S_(ut)\mathrm{S}_{\mathrm{ut}}Sut are used in MPa. The theoretical SN curves SN teo 0 % SN teo  0 % SN_("teo ")0%\mathrm{SN}_{\text {teo }} 0 \%SNteo 0% and SN teo 6 % SN teo  6 % SN_("teo ")6%\mathrm{SN}_{\text {teo }} 6 \%SNteo 6%, shown in Fig. 6, is generated by Eqn. (2A):
(2A) S f = { S u t N ( 1 3 log ( f ) ) 10 0 N < 10 3 a t e o N b t o o 10 3 N 10 6 S e N > 10 6 a t e o = ( f S u t ) 2 S e b t e o = 1 3 log ( f S u t S e ) (2A) S f = S u t N 1 3 log ( f ) 10 0 N < 10 3 a t e o N b t o o 10 3 N 10 6 S e N > 10 6 a t e o = f S u t 2 S e b t e o = 1 3 log f S u t S e {:(2A)S_(f)={[S_(ut)*N^(((1)/(3)*log(f))),10^(0) <= N < 10^(3)],[a_(teo)*N^(b_(too)),10^(3) <= N <= 10^(6)],[S_(e),N > 10^(6)]quada_(teo)=((f*S_(ut))^(2))/(S_(e))b_(teo)=-(1)/(3)*log((f*S_(ut))/(S_(e))):}:}S_{f}=\left\{\begin{array}{ll} S_{u t} \cdot N^{\left(\frac{1}{3} \cdot \log (f)\right)} & 10^{0} \leq N<10^{3} \tag{2A}\\ a_{t e o} \cdot N^{b_{t o o}} & 10^{3} \leq N \leq 10^{6} \\ S_{e} & N>10^{6} \end{array} \quad a_{t e o}=\frac{\left(f \cdot S_{u t}\right)^{2}}{S_{e}} b_{t e o}=-\frac{1}{3} \cdot \log \left(\frac{f \cdot S_{u t}}{S_{e}}\right)\right.(2A)Sf={SutN(13log(f))100N<103ateoNbtoo103N106SeN>106ateo=(fSut)2Sebteo=13log(fSutSe)
S f S f S_(f)\mathrm{S}_{\mathrm{f}}Sf is the theoretical fatigue resistance, f f fff is the fatigue strength fraction, N N NNN is the number of fatigue cycles, a teo a teo  a_("teo ")\mathrm{a}_{\text {teo }}ateo  and b teo b teo  b_("teo ")\mathrm{b}_{\text {teo }}bteo  are theoretical constants. Note that in Eqns. (1A) and (2A), for the as-received straight specimen, with no residual stress, the S e , S f S e , S f S_(e),S_(f)\mathrm{S}_{\mathrm{e}}, \mathrm{S}_{\mathrm{f}}Se,Sf, and S ut S ut S_(ut)\mathrm{S}_{\mathrm{ut}}Sut are substituted by S e 0 , S f 0 S e 0 , S f 0 S_(e0),S_(f0)\mathrm{S}_{\mathrm{e} 0}, \mathrm{~S}_{\mathrm{f} 0}Se0, Sf0, and S ut 0 S ut 0 S_(ut0)\mathrm{S}_{\mathrm{ut} 0}Sut0, respectively. This material is generically labeled as 0 % 0 % 0%0 \%0%. For the curved specimen at the beginning of phase 2, a 6 % 6 % 6%6 \%6% superficial residual strain is built, as shown in Fig. 2.a. Then, S e , S f S e , S f S_(e),S_(f)\mathrm{S}_{\mathrm{e}}, \mathrm{S}_{\mathrm{f}}Se,Sf, and S ut S ut  S_("ut ")\mathrm{S}_{\text {ut }}Sut  are substituted by S e 6 , S ff S e 6 , S ff S_(e6),S_(ff)\mathrm{S}_{\mathrm{e} 6}, \mathrm{~S}_{\mathrm{ff}}Se6, Sff, and S ut 6 S ut 6 S_(ut6)\mathrm{S}_{\mathrm{ut} 6}Sut6. This material is generically labeled as 6 % 6 % 6%6 \%6%. The fatigue stress distributions σ fad_max ( y ) σ fad_max  ( y ) sigma_("fad_max ")(y)\sigma_{\text {fad_max }}(\mathrm{y})σfad_max (y) and σ fad_min ( y ) σ fad_min  ( y ) sigma_("fad_min ")(y)\sigma_{\text {fad_min }}(\mathrm{y})σfad_min (y), at the middle of the specimen, can be estimated by:
(3~A) M min = P min L 4 M max = P max L 4 (4A) σ fad_ max ( y ) = M max y b 2 b exc ( radius 2 c ) σ x max = σ fad_ max ( c ) (5A) σ fad_ min ( y ) = R σ fad_ max ( y ) σ x min = σ fad_ min ( c ) (6A) exc = ρ i radius radius = b / [ ln ( ρ i + h 2 ρ i h 2 ) ] σ x m = σ x min + σ x max 2 σ x a = σ x max σ x min 2 (3~A) M min = P min L 4 M max = P max L 4 (4A) σ fad_ max ( y ) = M max y b 2 b exc radius 2 c σ x max = σ fad_ max ( c ) (5A) σ fad_ min ( y ) = R σ fad_ max ( y ) σ x min = σ fad_ min ( c ) (6A)  exc  = ρ i  radius   radius  = b / ln ρ i + h 2 ρ i h 2 σ x m = σ x min + σ x max 2 σ x a = σ x max σ x min 2 {:[(3~A)M_(min)=P_(min)*(L)/(4)quadM_(max)=P_(max)*(L)/(4)],[(4A)sigma_("fad_"_(-)max)(y)=(M_(max)*y)/(b_(2)*b*exc*(radius_(2)-c))quadsigma_(x_(-)max)=sigma_("fad_"_(-)max)(-c)],[(5A)sigma_("fad_"_(-)min)(y)=R*sigma_("fad_"_(-)max)(y)quadsigma_(x_(-)min)=sigma_("fad_"_(-)min)(-c)],[(6A)" exc "=rho_(i)-" radius "quad" radius "=b//[ln((rho_(i)+(h)/(2))/(rho_(i)-(h)/(2)))]],[sigma_(x_(m))=(sigma_(x_(-)min)+sigma_(x_(-)max))/(2)quadsigma_(x_(a))=(sigma_(x_(-)max)-sigma_(x_(-)min))/(2)]:}\begin{align*} & M_{\min }=P_{\min } \cdot \frac{L}{4} \quad M_{\max }=P_{\max } \cdot \frac{L}{4} \tag{3~A}\\ & \sigma_{\text {fad_}_{-} \max }(y)=\frac{M_{\max } \cdot y}{b_{2} \cdot b \cdot \operatorname{exc} \cdot\left(\operatorname{radius}_{2}-c\right)} \quad \sigma_{x_{-} \max }=\sigma_{\text {fad_}_{-} \max }(-c) \tag{4A}\\ & \sigma_{\text {fad_}_{-} \min }(y)=R \cdot \sigma_{\text {fad_}_{-} \max }(y) \quad \sigma_{x_{-} \min }=\sigma_{\text {fad_}_{-} \min }(-c) \tag{5A}\\ & \text { exc }=\rho_{i}-\text { radius } \quad \text { radius }=b /\left[\ln \left(\frac{\rho_{i}+\frac{h}{2}}{\rho_{i}-\frac{h}{2}}\right)\right] \tag{6A}\\ & \sigma_{x_{m}}=\frac{\sigma_{x_{-} \min }+\sigma_{x_{-} \max }}{2} \quad \sigma_{x_{a}}=\frac{\sigma_{x_{-} \max }-\sigma_{x_{-} \min }}{2} \end{align*}(3~A)Mmin=PminL4Mmax=PmaxL4(4A)σfad_max(y)=Mmaxyb2bexc(radius2c)σxmax=σfad_max(c)(5A)σfad_min(y)=Rσfad_max(y)σxmin=σfad_min(c)(6A) exc =ρi radius  radius =b/[ln(ρi+h2ρih2)]σxm=σxmin+σxmax2σxa=σxmaxσxmin2
where exc and radius are well-known variables of curved beams theory (because the specimens in phase 2 are curved). P min P min  P_("min ")\mathrm{P}_{\text {min }}Pmin  and P max P max  P_("max ")\mathrm{P}_{\text {max }}Pmax  are the minimum and maximum fatigue transversal loads, respectively. Consequently, M min M min  M_("min ")\mathrm{M}_{\text {min }}Mmin  and M max M max  M_("max ")\mathrm{M}_{\text {max }}Mmax  are, respectively, the minimum, and the maximum fatigue bending moments in the middle of the specimen. σ x _ min , σ x _ max , σ x _ m σ x _ min  , σ x _ max  , σ x _ m sigma_(x_"min "),sigma_(x_"max "),sigma_(x_m)\sigma_{\mathrm{x} \_ \text {min }}, \sigma_{\mathrm{x} \_ \text {max }}, \sigma_{\mathrm{x} \_\mathrm{m}}σx_min ,σx_max ,σx_m and σ x _ a σ x _ a sigma_(x_a)\sigma_{\mathrm{x} \_\mathrm{a}}σx_a are respectively, the minimum, the maximum, the average, and the alternate fatigue bending stress at y = c y = c y=-c\mathrm{y}=-\mathrm{c}y=c (point b b bbb of Fig.1.d), in the middle of the specimen. R R RRR is the stress ratio.
Note that the shear stresses are null, τ x y m = τ x y a = 0 τ x y m = τ x y a = 0 tau_(xy_(m))=tau_(xy_(a))=0\tau_{x y_{m}}=\tau_{x y_{a}}=0τxym=τxya=0, as the torsion stress is null, and the shear stresses caused by the transversal load are disregarded. So, using the von Mises criterion, it is possible to show that the Mises average stress is
σ m = σ x m σ m = σ x m sigma_(m)=sigma_(x_(-)m)\sigma_{m}=\sigma_{x_{-} m}σm=σxm and the Mises alternate stress σ a = σ x a σ a = σ x a sigma_(a)=sigma_(x_(-)a)\sigma_{a}=\sigma_{x_{-} a}σa=σxa. Eqn. (7A.a) converts the stresses σ m σ m sigma_(m)\sigma_{m}σm and σ a σ a sigma_(a)\sigma_{a}σa to one unique equivalent alternate stress σ a e q σ a e q sigma_(a_(-)eq)\sigma_{a_{-} e q}σaeq, based on Goodman Criterion, which is a Fatigue failure criterion which includes, in addition to the alternate stress, the average ones. The estimation of the theoretical finite life of the specimens N teo N teo  N_("teo ")\mathrm{N}_{\text {teo }}Nteo  is done:
(7A) σ a e q = ( 1 1 ( σ m S u t ) ) σ a N teo = ( σ a e q a teo ) 1 / b teo for 10 3 N 10 6 (7A) σ a e q = 1 1 σ m S u t σ a N teo  = σ a e q a teo  1 / b teo   for  10 3 N 10 6 {:(7A)sigma_(a_(-)eq)=((1)/(1-((sigma_(m))/(S_(ut)))))sigma_(a)quadN_("teo ")=((sigma_(a_(-)eq))/(a_("teo ")))^(1//b_("teo "))quad" for "10^(3) <= N <= 10^(6):}\begin{equation*} \sigma_{a_{-} e q}=\left(\frac{1}{1-\left(\frac{\sigma_{m}}{S_{u t}}\right)}\right) \sigma_{a} \quad N_{\text {teo }}=\left(\frac{\sigma_{a_{-} e q}}{a_{\text {teo }}}\right)^{1 / b_{\text {teo }}} \quad \text { for } 10^{3} \leq N \leq 10^{6} \tag{7A} \end{equation*}(7A)σaeq=(11(σmSut))σaNteo =(σaeqateo )1/bteo  for 103N106
When Eqn. (7A) uses the mechanical properties of an as-received material (labeled as 0 % 0 % 0%0 \%0% ), with no residual stress, S e , S f S e , S f S_(e),S_(f)\mathrm{S}_{\mathrm{e}}, \mathrm{S}_{\mathrm{f}}Se,Sf, and S ut S ut S_(ut)\mathrm{S}_{\mathrm{ut}}Sut are substituted by S e 0 , S f 0 S e 0 , S f 0 S_(e0),S_(f0)\mathrm{S}_{\mathrm{e} 0}, \mathrm{~S}_{\mathrm{f} 0}Se0, Sf0, and S ut 0 S ut 0 S_(ut0)\mathrm{S}_{\mathrm{ut} 0}Sut0, respectively. The point ( N teo , σ a e q / S y 0 N teo  , σ a e q / S y 0 N_("teo "),sigma_(a_(-)eq)//S_(y0)N_{\text {teo }}, \sigma_{a_{-} e q} / S_{y 0}Nteo ,σaeq/Sy0 ) is on the SN teo 0 % SN teo  0 % SN_("teo ")0%\mathrm{SN}_{\text {teo }} 0 \%SNteo 0% curve, which is graphically shown on Fig.5.