Ziya Peng·Ping Yang,2·Kang Hu ·Yueling Song
Abstract Welding residual stress in the engineering structure has a non-negligible influence on crack propagation,and crack closure is a significant factor affecting the crack propagation.Based on the elastoplastic finite element method and crack closure theory,we studied crack closure and residual compressive stress field of butt-welded plates under constant amplitude loading and overloading regarding the stress ratio,maximum load,overload ratio,and number of overloads.The results show that the welding residual tensile stress can decrease the crack closure because of a decrease in the residual compressive stress in the wake zone,but the effect is gradually reduced with increased stress ratio or maximum load.And the combined effect of welding residual tensile stress and overload can produce a stronger retardation effect on crack propagation.
Keywords Welding residual stress .Crack closure .Overloading .Crack propagation .Butt-welded plate
With the development of large ships,low-cycle fatigue fractures in ship hull structures have become more serious and are a key problem to address in ship development(Han et al.2007). When the crack extends to a certain length under cyclic loading, the hull structure will be failure. Therefore,accurately predicting the life of crack propagation is significant. Crack closure occurs in the process of crack propagation and is an essential factor affecting the crack growth rate.The closing of crack faces leads to reduction in the stress intensity factor range that reduces the crack growth rate.The hull is made up of a large number of welded structures,resulting in welding residual stress in the structure, which has a non-negligible influence on crack propagation. The combined effect of welding residual stress and residual stress caused by the external load will lead to a change in the crack closure,which affects the crack propagation rate.There have been few investigations into the effect of the welding residual stress on the butt-welded plate under constant amplitude loading and overloading.
Elber (1971) first proposed the concept of crack closure,suggesting that the crack faces might be in contact even under cyclic tensile loading.The contact leads to a decrease of the stress intensity factor range, which provides a new understanding of fatigue crack propagation. Antunes et al. (2015 2008 2010)performed a large amount of numerical simulation calculations investigating the influence of stress ratio, crack length,cycle number,and mesh size on crack closure,and also verified the validity of crack closure. Many results have shown that the plasticity (Matos and Nowell 2008; Cochran et al.2011)is still the primary factor to induce crack closure,and the residual compressive stress produced by the residual plastic deformation in the wake zone has been widely studied.Ding et al.(2013)focused on the effect of loading history on crack propagation.The study shows that the residual compressive stress in the wake zone is an important factor to promote crack closure during loading.Li et al. (2014) carried out the tensile and compression overloading analysis of MT specimens by numerical simulation. The results conclude that the increase of residual compressive stress in the wake zone caused by the tensile overloading promotes crack closure.Compression overloading will have a reverse effect on crack closure,however,the effect is small.
To study the effect of the welding residual stress on crack propagation in engineering structures,relevant researchers introduced welding residual stress as the initial residual stress into a model for analysis. Lv and Wang et al. (2016) introduced the initial residual stress of shot peening into the model and investigated its influence on crack propagation and crack closure.The study shows that the crack closure is caused by the combined effect of external loading and initial residual compressive stress in shot peening, and the distribution of initial residual stress in shot peening also plays an important role in the crack growth rate. Wang et al. (2014) considered the effect of the initial residual stress and established an improved crack propagation model of stiffened plates. The results are in good agreement with the experimental results and can be used in engineering practice.In some studies(Terada 2011;Tan et al.2007;Bao and Zhang 2010;Bao et al.2010),the effect of the initial residual stress is mainly expressed in the form of the stress intensity factorKre,but these works do not investigate the generation process of the welding residual stress. Therefore, some researchers (Lee and Chang 2011)have simulated the process of welding by numerical simulation,and then they treated the residual stress obtained by the simulation as the initial stress in the model, which is a more reasonable method to investigate the effect of welding residual stress on crack propagation. Of course, the welding residual stress will be released to a certain extent and appear as redistribution during the crack propagation process, which will affect the crack growth rate.This effect is more obvious under different overloads (Liljedahl et al. 2008; Huang and Moan 2007; Li et al. 2007). There are two primary methods for predicting the crack growth rate considering the welding residual stress.One is to determine the effective stress ratio Reffby the superposition principle to explain the effect of the welding residual stress. The other method is calculating the open stress intensity factor Kopand effective stress intensity factor range ΔKeff under the combined effect of the external load and welding residual stress field based on the principle of crack closure. Liljedahl et al. (2009) used the residual stress approach and crack closure approach to assess the crack growth rate,and compared the finite element results with the experimental results.It was found that the results obtained by the residual stress approach were in good agreement with those from the experiment under low loading levels,and the results obtained by the crack closure approach were in good agreement with those from the experiment under high loading level.
The discussed studies show that the welding residual stress has a non-negligible influence on crack propagation.The concept analysis of fatigue crack closure can reflect the influence of residual stress on the crack propagation law more profoundly(Wang et al.1999;Kim et al.1991;Larue and Daniewicz 2007).The existing research on crack closure and crack propagation law of a butt-welded plate considering the welding residual stress under constant amplitude loading and overloading is limited. In this paper, the butt-welded plate with a center penetrating crack is considered,and the welding residual stress is taken as the initial stress in the model investigating the influence of the stress ratio,maximum load,overload ratio,and number of overloads based on the elastoplastic finite element method and crack closure theory. The results show that the residual compressive stress field formed by the combined effect of welding residual stress and constant amplitude loading and overloading have different effects on the crack closure under various loading cases, which provides a new way to study the complex mechanism of crack closure and crack propagation law for the butt-welded plate.
In the weight function method,the superposition principle is usually used to superimpose the stress intensity factors caused by the external load and welding residual stress to predict the crack propagation life.However,it is assumed that the superposition principle is only valid if the residual stresses are not affected by the presence and growth of a fatigue crack.Due to the closing of crack faces and release of welding residual stress during the crack propagation,the crack closure mechanism has been widely used to study the combined effect of the residual stress field caused by the external load and welding residual stress field (Lv et al. 2016; Wang et al. 2016); the effective stress intensity factor range is expressed as:
where Kopis the stress intensity factor corresponding to the combined effect of the external load and welding residual stress.
In a previous study,the crack closure parameter U is often used to describe the crack closure level:
To facilitate the study of the distribution of welding residual stress, the analysis model adopted in this paper is shown in Figure 1.The sizes of the butt-welded plate are 2L=160 mm,2W=80 mm, t=3 mm, a0=4 mm. The yield stress is 345 MPa and the applied maximum cyclic loading is 320 MPa. The load direction is applied along the weld line(longitudinal direction) (Bao et al. 2010; 2010). During the analysis, the geometry of the weld is not considered under the assumption of a constant plate thickness t and a crack length exceeding the width of the weld.
Figure 1 Geometric model
Because the model is welded by two flat plates with symmetry, 1/4 of the model was built in Figure 2 to simplify the analysis.To better simulate crack closure and crack propagation,an analytical rigid line is established at the bottom of the model as a master surface,and the bottom surface of the model is used as a slave surface.Because the plate thickness is much smaller than the length and width of the plate, the model is assumed to be analyzed under plane stress. The boundary condition is taken as follows:a fixed constraint in the x-direction is applied along the y-axis direction. To prevent mutual embedding of materials, an analytical rigid line is defined along the x-axis direction, and a fixed constraint in the y-direction is applied in the region where the crack stops expanding.
This paper uses node release technology to simulate the process of crack propagation.Due to the influence of calculation time and releasing node time,crack growth is simulated by successive node release at minimum load once every two cycles of loading (Matos and Nowell 2009). Crack opening loads are determined by monitoring the contact state of the first node behind the crack tip(Borrego et al.2012;Dong et al.2019).To improve the accuracy of the opening values when the crack actually opens,linear interpolation is used between the remote loads corresponding to the last contact in the first node behind the crack tip and the following load increment(Allison 1988). The model is meshed using the plane stress four-node reduction integral element(CPS4R).The expanded area adopts a refined mesh with a mesh size of 0.05×0.05 mm to ensure that the cyclic plastic zone contains at least 3 to 4 elements(Lei 2008;Li and Li 2006).
Figure 2 Finite element model
The material properties parameters are σy=345 MPa, E=206GPa,and v=0.3.Because plasticity induced crack closure(PICC)is a plastic deformation-based phenomenon,it is very important to accurately simulate the yield and hardening behavior of the materials.The Chaboche combination model can better capture the basic cyclic plastic response of materials such as the Bauschinger effect,plastic shakedown,ratcheting resulting from asymmetric cycles of stress, and mean stress relaxation resulting from an asymmetric strain cycle.Therefore,the Chaboche combinatorial constitutive model is adopted to characterize the nonlinear behavior of the materials,and the model parameters shown in Table 1 are selected from literature(Dong et al.2016).
In Table 1,v is Poisson’s ratio,Q is the maximum change of the yield surface,b is the change rate of the yield surface with increasing plastic strain, c1, c2, and c3are the initial kinematic hardening moduli, and γ1, γ2, and γ3are the rate of strain increase.
The welding residual stress should be self-balanced.However,in the experimental measurement,the measured tensile residual stress is often greater than the compressive residual stress leading to the unbalanced stress field. This stress field will influence the calculation results(Bao et al.2010).The purpose of this paper is not to simulate the welding residual stress distribution but to better study the redistribution and release of welding residual stress during crack propagation and its effect on crack propagation under constant and variable amplitude loading (Liljedahl et al. 2008). Therefore, this paper adopts the empirical simplification formula of the welding residual stress proposed by Tada and Paris (1983). Due to the mode-I crack propagation,the longitudinal residual stress has a more important role than the transverse residual stress.To simplify the calculation,only the distribution of the longitudinal residual stress is considered:
Table 1 Chaboche composite model material parameters
where σ0yis a parameter defining the maximum value of the tensile residual stress and cois the distance from the y-axis at which the residual stress field changes from positive to negative,i.e.,from tensile to compression.Because the stress at the weld is close to or exceeds the material yield stress,σ0ycan be taken as 345 MPa, and cohas a typical value of 10 mm(George and Ioannis 2009).Considering the calculation efficiency,it is not advisable to choose a larger crack propagation length.Through a large number of finite element calculations,this paper finally chooses to expand only Δa=5mm (The reasons for selection of coinclude the following. This paper mainly studies the effect of the ideal residual stress distribution on crack propagation but does not study the specific welding process. In addition, because only the effect of the residual tensile stress is investigated, it is desirable that cois larger than the sum of the crack extension length and initial crack length (c0>Δa+a0:10mm >5mm+2mm)). The residual stress distribution is shown in Figure 3.
The residual stress file is written in the FORTRAN program and introduced into the calculation model as initial stress by ABAQUS definition user subroutine SIGINI. Since the introduced initial stress cannot fully satisfy the selfbalancing condition of the residual stress and boundary conditions,a static analysis is needed to make the residual stress in the test material become the equilibrium stress state.Figure 4 shows the introduced residual stress distribution before and after equilibrium and Eq.4.The difference between the residual stress distribution of pre-equilibrium, post-equilibrium,and the formula is very small.Therefore,the residual stress distribution after equilibrium is used as the initial residual stress in the calculation model.
Selected literature (Lv et al. 2016) was used to verify the validity of the finite element method used in this paper.The test model is a four-point bending specimen. First, the finite element model shown in Figure 5 is established by ABAQUS software.This work also uses the node release technique simulating crack growth.The introduced welding residual stress is realized by the ABAQUS definition user subroutine. The crack opening load is determined by the displacement state of the first node behind the crack tip.As shown in Figure 6,the results calculated by finite element analysis in this paper show very good agreement with literature.Therefore,the finite element method used has reliability.
Figure 3 The distribution of welding residual stress
Figure 4 Stress distribution before and after equilibrium
To study the effect of the welding residual stress on crack propagation under constant amplitude loading, the analysis was performed under a stress ratio of R=0.Figure 7 a and b show the curves of crack opening load and crack closure parameter with and without welding residual stress, indicating the welding residual stress has a non-negligible influence on crack propagation.
The crack closure weakened in the welding tensile stress region primarily because the residual tensile stress is able to offset some part of the residual compressive stress formed by the external load during unloading. Additionally, the difference between crack closure parameters in the two different cases continuously decreases as the crack length increases,which is probably related to the release of residual stress.This can be considered that the effect on crack propagation is also gradually reducing.
Figure 5 Finite element model and boundary condition(Lv et al.,2016)
Figure 6 Comparison of the result
To investigate the evolution of the residual compressive stress field in the plastic wake zone, the paper uses σrre(2),σre(2)to represent the size of the residual compressive stress field with or without welding residual stress for the crack extension length Δa=2mm (σrre,σre are the sum of all the compressive stress values behind the crack tip at the minimum load with or without welding residual stress, respectively).Figure 7 c and d show the stress distribution near the crack tip for Δa=2mm.It can be seen that the residual compressive stress field in the plastic wake zone is reduced at the minimum load due to the welding residual tensile stress (|σrre(2)|=999MPa <|σre(2)|=2784MPa). Because the residual compressive stress field is the primary factor to affect crack closure, the decreased residual compressive stress will weaken the crack closure to a certain extent,resulting in the increase in crack growth rate.The work also extracts accumulative plastic strain at the crack tip for Δa=2mm, and the accumulative plastic strain value with welding residual stress is larger than the situation without residual stress (εpre=0.0373 >εp=0.021),which indicates that the welding residual tensile stress can increase the plastic deformation at the crack tip. This conclusion agrees with the literature(Li et al.2011).
The effect of the stress ratio on crack propagation is also significant.The crack propagation in four loading cases with different stress ratios(R=0,0.1,0.2,and 0.3)was studied to investigate the crack closure parameter U and residual stress field distribution.
Figure 7 The effect of welding residual stress on Pop,U,and stress distribution
Based on these results, the crack closure will weaken with an increase in the stress ratio, which is mainly related to the decrease in the residual compressive stress field in the wake zone behind the crack tip. Figure 8c shows that the residual compressive stress field will be reduced with the increase in the stress ratio (|σre(0.1)|=1974MPa <|σre(0)|=2993MPa,|σrre(0.1)|=1023MPa<|σrre(0)|=1223MPa). The residual compressive stress field formed by the external load and welding residual tensile stress is smaller than with only the external load(|σrre(0.1)|<|σre(0.1)|, |σrre(0)|<|σre(0)|).
Quantitative analysis was used to investigate the effect of the stress ratio on the welding residual stress.S0,S0.1,S0.2,and S0.3are used to represent the area between crack closure parameter curves with or without welding residual stress under different stress ratios(R=0,0.1,0.2,0.3,respectively),which are the effect of the welding residual stress on crack closure and crack propagation.As shown in Table 2,the effect of the welding residual tensile stress on crack closure gradually reduced with an increase in stress ratio(Uav,Uav(re)),resulting in a decrease in the effect of the welding residual stress on the crack growth rate (S0>S0.1>S0.2>S0.3), which is closely related to the release and redistribution of welding residual stress during crack propagation.
The maximum load also has a non-negligible effect on crack propagation. Cases with maximum load (Pmax=320, 370,420 MPa) under R=0.1 were analyzed. As shown in Figure 9a, the increased maximum load decreases the crack closure without welding residual stress primarily because the load can offset part of the residual stress on the crack faces(|σre(320)|=1973MPa >|σre(370)|=1672MPa), as shown in Figure 9b. Alternatively, the law of crack closure is different with welding residual stress.When Δa >3 mm,the crack closure decreases with an increase in the maximum load.Figure 9c shows the stress field distribution at the minimum load for Δa=5mm, illustrating the compressive stress field has some reduction (|σre(320)|=1589MPa >|σre(370)|=1386MPa). However, the variation law is not obvious for Δa= <3mm, which is mainly related to the
Fig.8 The effect of the stress ratio on crack closure
Table 2 Effect of the stress ratio
Where Uav(re) and Uav represent the average value of the crack parameter with or without the welding residual stress, respectively, i.e., the level of crack closure. Sidescribes the effect of stress ratio on crack closure and crack propagation with welding residual stress (i=0, 0.1 0.2, 0.3, …, …, n)distribution of the welding residual stress. The same quantitative method was performed to investigate the influence of the maximum load for the welding residual stress. As shown in Table 3, overall, the increased maximum load will reduce the influence of the welding residual stress on crack closure (S320>S370>S420), indicating that the increased maximum load may promote the release and redistribution of welding residual stress.
The stress intensity factor range is a significant factor affecting crack propagation.The stress intensity factor can be calculated by Feng et al.(2012):
The effective stress intensity factor ranges with or without welding residual stress are obtained by Eq. (7). As shown in Figure 10 a and b, the effective stress intensity factor range with welding residual stress is larger than the case without welding residual stress, indicating the welding residual tensile stress can affect the crack growth rate to some degree. When the stress ratio increases, the effective stress intensity factor range will decrease, indicating a decrease in the crack growth rate.Additionally,the effective stress intensity factor range still increases with an increase in the maximum load.However,the effect caused by the increased maximum load on welding residual stress appears to be a reduction.
Figure 9 The effect of the maximum load on crack closure
Table 3 The effect of the maximum load
Figure 10 The variety of effective stress intensity factor range under different constant amplitude loads
The influence of overload ratio, number of overloads, and other factors were used to study the effect of the welding residual stress on crack propagation under overloading. The main purpose is to consider the effect of overload on the welding residual stress redistribution that result in complex stress changes at the crack tip,thus affecting the crack growth rate.The crack closure parameter and change of residual compressive stress field with or without welding residual stress are investigated for ROL=1.2(where ROL=POL/Pmax,POLis the maximum load during the overload) and Δa=1mm in this section.
As shown in Figure 11a,crack closure appears to decrease first after the overload,then gradually increase,finally recovering to constant amplitude loading. Figure 11b shows that residual compressive stress field in the plastic wake zone will be enhanced after the overload (Δa=1.3m, |σre(1.2)|=2 3 4 5 M P a >|σ r e(1)| = 1 6 6 7 M P a, |σ r r e(1.2)| =2024MPa >|σrre(1)|=931MPa), and the residual compressive stress field after the overload will relatively reduce under the effect of the welding residual tensile stress (|σrre(1.2)|=2024MPa <|σre(1.2)|=2345MPa), which is the main reason for the decrease of crack closure and increase of the crack growth rate. In addition, the accumulative plastic strain(Δa=1mm) will increase at the time of overload(εpre(1.2) = 0.0644 >εpre(1) = 0.0349, εp(1.2) =0.0394 >εp(1)=0.0171),and the accumulative plastic strain with welding residual stress is larger(εpre(1.2)>εp(1.2)),indicating that the overload and welding residual tensile stress will increase the plastic deformation at the crack tip.
Figure 11 U and stress distribution near the crack tip under overloading
A similar literature method was used to quantify the effect of the welding residual stress on crack propagation under different overloads(Chen et al. 2018).S(re)and S0 are used to represent the influence of the retardation effect during the overload with or without the welding residual stress in Fig.11a, respectively. The values can be obtained by numerical integration.The increase in S0and S(re)implies that the retardation effect is enhanced,and it can be expressed as:
where S1and S2represent the acceleration and deceleration effects of crack propagation without welding residual stress,respectively.S1(re)and S2(re)represent the acceleration and deceleration effects of crack propagation with welding residual stress, respectively. Δa1, Δa2, and Δa3are described in Figure 11a.
As shown in Figure 11a (S(re)=0.0864 >S0=0.0139), the retardation effect with the welding residual tensile stress caused by the overload is larger, which is similar to previous reports(Daneshpour et al.2012).The residual compressive stress field shown in Table 4 has a large proportion of change from the position of pre-overload to the strongest retardation effect after overload due to the combined effect of residual tensile stress and overload.This leads to the fact that the decreased amplitude of ΔU is larger(the amplitude of ΔU refers to the difference between the minimum value of Umin(overload)after the overload and Uconstant,and the definition of Uconstant is the intersection of the overload curve and constant amplitude loading curve when the overload curve decreased,as shown in Figure 11a;when Urecovers to the level of constant amplitude loading,the retardation effect will be also larger.
Table 4 The change in the residual compressive stress field
In Table 4, σwrs is the welding residual stress, “×” indicates the case without welding residual stress, “√” indicates the case with welding residual stress,a1is the position before the overload, a2is the position of the strongest retardation effect after overload,σre(1),σrre(1)indicate the residual compressive stress field without or with welding residual for the Δa=a1,and σre(2),σrre(2)indicate the residual compressive stress field without or with welding residual for Δa=a2.
Figure 12 The effect of the overload ratio on crack closure
To investigate the effect of the overload ratio on welding residual stress, the analysis is carried out for four cases with different overload ratios (ROL=1.2, 1.3, 1.5, and 1.7).The relationship between the crack closure parameter and crack propagation length with or without welding residual stress under different overload ratios is shown in Figure 12 a and b. Based on the quantitative result of the retardation effect shown in Table 5,the retardation effect is gradually enhanced with the increase in overload ratio in two different cases,and the retardation effect with welding residual stress is larger(S(re)>S0,ΔUre >ΔU0).In addition,with the increase in the overload ratio,the value of U is more difficult to recover to the level of constant amplitude loading after overloading.The difference between the overloading and constant amplitude loading gradually increases indicating the increased overload ratio may somewhat reduce the effect of welding residual stress.
To investigate the effect of number of overloads on welding residual stress, the analysis is performed for four cases with different number of overloads (N=1, 3, 7,and 11). Based on the quantitative result of the retardation effect shown in Table 6, the retardation effect is gradually enhanced with an increase in the number of overloads, and the retardation effect with the welding residual stress is larger (S(re)>S0, ΔUre >ΔU0) than the case without welding residual stress (Figure 13).Although the retardation effect increases with the increased number of overloads, the trend of the crack closure parameters is relatively close in the late stage, which indirectly indicates that the increased number of overloads do not have an obvious effect on the release of welding residual stress.
To express the law of crack propagation during overload,the effective stress intensity factor range is obtained by Eq.(7).As shown in Figure 14,the crack parameter U suddenly increases and then gradually decreases,finally returning to the level of constant amplitude loading.When the overload ratio and number of overloads increase, the retardation effect is relatively enhanced resulting in the reductionin the crack growth rate.The retardation effect is relatively stronger with the welding residual stress. The increased overload ratio may also reduce the effect of welding residual stress to a certain extent. In contrast, the effect of the number of overloads is not very obvious.
Table 5 The effect of the retardation effect
Table 6 The effect of the retardation effect
Figure 13 The effect of the number of overloads on crack closure
Figure 14 The variety of effective stress intensity factor range under different overloads
Welding residual stress has a non-negligible influence on crack propagation. Based on the crack closure method, the effect of the welding residual stress on crack propagation was investigated under constant amplitude loading and overloading,and the influencing factors of stress ratio,maximum load,overload ratio,and number of overloads discussed.The factors were quantitatively analyzed to understand the effect of the welding residual stress on the crack propagation law.The following conclusions were drawn:
1) Welding residual tensile stress can offset part of the residual compressive stress caused by the external loading,reducing the total residual compressive stress field in the plastic wake zone, which can reduce the crack closure thereby increasing the crack growth rate.
2) The effect of the welding residual tensile stress on the crack closure continuously decreased with the increase in the stress ratio under constant amplitude loading.The ability of the welding residual tensile stress to offset the residual compressive stress in the wake zone tends to decrease,which is related to the release and redistribution of the welding residual stress during crack propagation.
3) The effect of the maximum load is somewhat different with or without welding residual stress. Under no welding residual stress,the crack closure decreased with an increase in the maximum load.Alternatively,the variable law of the crack closure is not very obvious. Only when the crack length extends to a certain extent, the crack closure will decrease with the increase in the maximum load.
4) Under the combined effect of the overload and welding residual tensile stress,the retardation effect will increase to some degree primarily because the residual compressive stress field has a significant increase resulting in a large decrease of the crack closure parameter. When U recovers to the level of constant amplitude loading, the retardation effect also gradually increased.
5) The retardation effect caused by overload is enhanced with the increase in the overload ratio and number of overloads,and the effect with the welding residual tensile is particularly obvious.When the overload ratio increases,it will promote the release of the welding residual stress and reduce its influence on the crack propagation to a certain extent. In contrast, the effect of the number of overloads is not very obvious.
6) The effective stress intensity factor range can reflect the crack growth rate well. In the case of welding residual tensile stress,the crack growth rate is relatively large.A decreased maximum load or increased stress ratio will reduce the effective stress intensity factor range,resulting in a reduction in the crack growth rate.For the case of the overload,the retardation effect is relatively enhanced with the increase in the overload ratio and number of overloads.The retardation effect caused by the overload under welding residual stress is greater.
Funding This study was supported by the National Natural Science Foundation of China(Grant No.51779198).
Journal of Marine Science and Application2020年2期