Next Article in Journal
Polydopamine-Modified Al2O3/Polyurethane Composites with Largely Improved Thermal and Mechanical Properties
Previous Article in Journal
Glass/Au Composite Membranes with Gold Nanoparticles Synthesized inside Pores for Selective Ion Transport
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Extremely-Low-Cycle Fatigue Damage for Beam-to-Column Welded Joints Using Structural Details

1
School of Civil Engineering, Hubei Engineering University, Xiaogan 432000, China
2
Hubei Key Laboratory of Roadway Bridge & Structure Engineering, Wuhan University of Technology, Wuhan 430070, China
3
School of Civil Engineering, Shandong Jianzhu University, Jinan 250101, China
*
Author to whom correspondence should be addressed.
Materials 2020, 13(7), 1768; https://doi.org/10.3390/ma13071768
Submission received: 21 February 2020 / Revised: 2 April 2020 / Accepted: 3 April 2020 / Published: 9 April 2020
(This article belongs to the Section Advanced Materials Characterization)

Abstract

:
The multiaxial fatigue critical plane method can be used to evaluate the extremely-low-cycle fatigue (ELCF) damage of beam-to-column welded joints in steel frameworks subjected to strong seismic activity. In this paper, fatigue damage models using structural detail parameters are studied. Firstly, the fatigue properties obtained from experiments are adopted to assess ELCF life for steel frameworks. In these experiments, two types of welded specimens, namely, plate butt weld (PB) and cruciform load-carrying groove weld (CLG), are designed according to the structural details of steel beam and box column joints, in which both structural details and welded factors are taken into account. Secondly, experiments are performed on three full-scale steel welded beam-to-column joints to determine the contribution of stress and/or strain to damage parameters. Finally, we introduce a modification of the most popular fatigue damage model of Fatemi and Socie (FS), modified by us in a previous study, for damage evaluation, and compare this with Shang and Wang (SW) in order to examine the applicability of the fatigue properties of PB and CLG. This study shows that the modified FS model using the fatigue properties of CLG can predict the crack initiation life and evaluate the damage of beam-to-column welded joints, and can be subsequently used for further investigation of the damage evolution law.

1. Introduction

Steel beam-to-column welded joints may undergo extremely-low-cycle fatigue (ELCF) during strong seismic activity [1,2,3], which is usually characterized by a few reverse loading cycles (in general, less than 20) with large strain amplitudes before failure. Welded joints are crucial components in antiseismic steel structures; for instance, the 1994 Northridge and 1995 Kobe earthquakes which caused severe damage to steel structures, where cracks were initiated in welded joints and eventually caused the collapse of whole structures [4,5]. The structural details of welding defects present possible new critical positions for ELCF crack initiation, which plays an important role in assessing fatigue crack initiation life and damage. The numerous available methods for seismic damage evaluation of steel frame structures are classified into two categories. The first category mainly focuses on establishing the damage evaluation model using macroparameters, such as plastic deformation [6], displacement amplitude [7,8], and rotation amplitude [9] combined with the low-cycle S-N curve, which describe the relationship between strain and the number of cycles to failure. The second category is based on cumulative damage of low-cycle fatigue, which takes plastic deformation of structures or components [10,11,12], energy dissipation [13,14,15,16,17], or the combination of these two variables [18,19,20,21,22] as a seismic damage parameter. In addition, numerous experiments have also been carried out to evaluate the structural damage of seismic excitations. However, existing studies mostly focus on the ductility and energy dissipation capability of the whole structure or structural component, and ignore the local details of welded joints and the fatigue crack initiation mechanics during the ELCF fracture process, which may lead to inaccurate damage evaluation results.
With the application of the micromechanical approach in fracture predictions of welded joints, researchers have paid attention to investigating the fatigue failure behavior [23,24,25,26,27] and severe seismic damage assessments [28,29,30] from the aspect of structural details. From these aspects, the results can be more suitable than traditional methods [26,28]. However, welded joints are always in multiaxial and nonproportional stress states due to their geometrically discontinuous properties, which result in more fatigue damage [31]. Therefore, a suitable model is necessary to correctly describe the effect of multiaxial and nonproportional stress states in fracture predictions for welded joints. The widely accepted critical plane approach is used to determine the damage parameter in a multiaxial nonproportional fatigue damage model. The critical plane approach usually uses the maximum shear plane as the critical damage plane, and takes the maximum shear strain γ max and normal strain ε n on the critical plane as the two basic damage parameters. The most popular multiaxial fatigue damage model combining strain with stress components is that proposed by Fatemi and Socie (FS) [32], which considers the nonproportional influence. Based on this study, numerous researchers have improved the nonproportionality parameter in the FS model based on their own research backgrounds [31,33,34,35], resulting in specific engineering applications. Another popular damage parameter using a path-independent parameter, Δ ε e q c r / 2 , composed of normal strain excursion ε n * and maximum shear strain γ max , was proposed by Shang and Wang (SW) [36], but has not yet been used in steel welded joints. Furthermore, the application of multiaxial fatigue damage modeling in ELCF failure prediction for steel welded joints is also limited. In conclusion, the research and application of multiaxial fatigue damage evaluation of steel beam-to-column welded joints under ELCF loading conditions still require further study.
This paper aims to propose a fatigue damage parameter for ELCF crack initiation life prediction and damage assessment for steel beam-to-column welded joints, in which the I-beam and box-column profile is considered in the joint. Firstly, ELCF tests are carried out by designing two types of welded specimens according to welded structural details for obtaining material parameters in multiaxial fatigue damage modeling. Secondly, reverse cyclic loading tests are performed on three full-scale welded joints to obtain the local stress and strain response of welded joints, and the experimental data are used to form the fatigue damage parameter. In addition, the FS model is modified from our previous study using the nonproportionality factor k * , based on both the nonproportionality of strain path Φ and material nonproportional hardening coefficient L . Finally, the modified FS model and the SW model are checked for their accuracy and reliability in assessing the fatigue damage of steel beam-to-column welded joints.

2. Experiments

In this section, ELCF tests of welded structural details and reverse cyclic loading tests of I-beam to box-column welded joints are carried out to obtain the parameters in the multiaxial fatigue damage model.

2.1. ELCF Test

The material and fatigue parameters in the heat affected zone (HAZ) are usually much lower and worse than the corresponding properties of the base materials. Indeed, fatigue failure of welded steel structures may happen at the HAZ. However, according to the damage investigation of the Northridge earthquake in 1994, it was found that about 100 steel frameworks suffered from different degrees of damage at the beam-to-column joints. It was also shown that the cracks originating in the weld zone occupied comprised as much as 90% of the total number, while only about 10% originated from the base metal. There were obvious cracks initiating at the weld seam of the beam-to-column welded joints, and the cracks propagated to the members connected to the joints, causing the failure of the steel frame. It is difficult to determine the location of the HAZ and to take a sample for analysis. Therefore, we focus on the weld zone in this paper.
As shown in Figure 1, two types of welded structural detail specimens were designed to carry out ELCF tests in order to investigate the fatigue properties of the welded structural details of box-column joints in a high-rise steel frame structure. There were twelve specimens in each type. The first type, namely, plate butt (PB) weld specimens, were used to simulate the welded splicing between the beam plate or between the column plate. The second type, namely, cruciform load-carrying groove (CLG) weld specimens, were used to simulate the welded connection between the beam flange plate and column flange plate. Both of the weld forms were single-sided groove welds. The base material was Q235 low-carbon steel. The mechanical properties of tested material are listed in Table 1.
These fatigue tests were conducted on an INSTRON1342 fatigue test machine, and controlled by uniaxial constant strain amplitude using sinusoidal wave forms, as shown in Figure 2. The strain ratio, i.e., the ratio of maximum strain to minimum strain, was −1 in these tests.
The experimental results were listed in [27] as shown in Table 2. The identifier for the specimen number corresponds to the loading condition; for instance, PB04 and PB10 mean that the tests for PB specimens were controlled by total strain amplitudes of 0.004 and 0.010, respectively. The meanings of the symbols in Table 2 are the same as those in Equations (1) and (2). For either loading conditions, Δ ε 2 , Δ ε e 2 , Δ ε p 2 and Δ σ 2 are obtained based on the stable hysteresis loop curve, as shown in Figure 3.
The fatigue properties of the welded structural details, which are correlated by the Manson–Coffin equation expressed by Equation (1), are listed in Table 3. Also listed in this table are the cyclic stress–strain properties when the stable hysteresis loop curve is fitted into the Ramberg–Osgood relation expressed by Equation (2). The fatigue and cyclic stress–strain properties are fitted using General model Power1 in MATLAB, that is, f ( x ) = a x b .
Δ ε 2 = Δ ε e 2 + Δ ε p 2 = σ f E ( 2 N f ) b + ε f ( 2 N f ) c
Δ ε 2 = Δ ε e 2 + Δ ε p 2 = Δ σ 2 E + ( Δ σ 2 K ) 1 n
where E is Young’s modulus; N f is the number of cycles to failure, σ f , ε f are fatigue strength coefficient and fatigue ductility coefficient, respectively, b and c are fatigue strength exponent and fatigue ductility exponent, respectively, and Δ ε e 2 = σ f E ( 2 N f ) b , Δ ε p 2 = ε f ( 2 N f ) c , and Δ ε 2 = Δ ε e 2 + Δ ε p 2 are, respectively, elastic strain amplitude, plastic strain amplitude, and total strain amplitude. In Equation (2), Δ σ 2 is stress amplitude, K is the cyclic strength coefficient, and n is the cyclic strain hardening exponent.

2.2. Reverse Cyclic Loading Test

Three cruciform beam-to-column welded joints were designed in this experiment, in which the inflection point of the beam-to-column connection was selected. The column had box sections with a size of 300 × 300 × 16 × 16 (mm) and had a height of 2 m. The beam had an I-shaped cross-section with a size of 250 × 150 × 8 × 10 (mm) and a beam length of 1.35 m. The weld form between the beam flange and the column flange was a full penetration groove using ER50-6 gas shielded welding wire. Welding process quality inspection was performed using a digital ultrasonic flaw detector; the grade of the weld seam was I. The beam-to-column welded joint specimens and welded structural details are shown in Figure 4.
The test device is shown in Figure 5. The DH5956 dynamic signal analysis acquisition system and multiple JM3812 multifunctional static strain data acquisition systems were applied in this test, and a 200-fold ultra-high-definition electron microscope was used to observe the cracks in the beam-to-column welded joints. To facilitate the statistics of crack growth and strain response, the left and right directions of Figure 5 are marked as west and east, respectively, and the other two orthogonal directions are marked as north and south, respectively.
According to the numerical simulation results of ABAQUS finite element analysis, the cruciform I-beam to box-column welded joint is in yield state when the displacement is greater than 12 mm, which is the noted yield displacement. It can be seen from Figure 6 that at time moment t = 2.141s, when the simulation is controlled by yield displacement, the model has maximum von Mises stress at the weld roots of the upper and bottom flanges with a stress value of 255.2 MPa, which is greater than the yield strength value shown in Table 1. This indicates that the weld zone is in the stage of plastic deformation and the weld root on the beam flange plate of the welded connection is in a position of potential fatigue failure. Hence, strain rosettes are arranged in the vicinity of the dangerous position to obtain the strain response in this test, as shown in Figure 7. Moreover, a reciprocating displacement greater than the yield displacement is selected for the cyclic loading test in order to simulate the stress state of welded joints of the steel framework when undergoing strong earthquake actions. In the test, actuators are controlled synchronously and in reverse on the beam ends. The loading histories are shown in Table 4.
The bottom flange at the northwest of the second specimen, which is marked with NWB_2, was selected to analyze the strain response. This was validated by a test showing that the weld root on the beam flange plate of the welded connection is in a position where fatigue cracks are likely to initiate. From the experiment, it was observed that under cyclic loading with a constant displacement amplitude, the tiny visible crack in NWB_2 was initiated just after the 7th cycle, i.e., at the time period between the 7th and 8th cycle. Strain data are obtained by placing three strain rosettes at the three coordinate planes in the vicinity of the weld root. Three normal strains and three shear strains of the weld root are then calculated using the elastic mechanics theory. The extracted strain response time history curves of NWB_2, ranging from the beginning of the trial to crack initiation, are drawn in Figure 8.
It can be seen from Figure 8 that the sizes of the six strain components do not change in a fixed proportion, nor do they reach the peak and/or valley at the same time. There is a certain phase difference between strains, that is, the strain response has nonproportional characteristics. Furthermore, the force characteristic is usually multiaxial in the welded zone due to complex geometric characteristics. In fact, a complex component, which has a sudden change of geometric shape, such as a notch or a weld, is always in a multiaxial stress–strain state, even though it is subjected to uniaxial load. That is, the failure of a beam-to-column welded joint which undergoes severe seismic excitation is multiaxial and nonproportional ELCF. Therefore, the widely used multiaxial fatigue critical plane method was adopted to evaluate the fatigue damage of the beam-to-column weld joint in this paper.

3. Methods of Evaluation

In this section, two fatigue critical plane methods are used. We outline damage analysis steps and fatigue damage parameters ranging from strain response on the critical plane to nonproportional factors, and cumulative damage criteria are introduced.

3.1. Damage Analysis Steps

The multiaxial fatigue analysis method is used to analyze the extremely-low-cycle fatigue damage of beam-to-column welded joints; the steps are shown in Figure 9 and summarized as follows:
(1)
The obtained strain responses are used to determine the candidate material plane of the dangerous position based on the coordinate transformation. Then, the critical damage plane is determined by means of the average weighted method.
(2)
The cyclic rain flow counting method is applied to count the normal and shear strain response time histories on the critical plane. Then, the stress or strain components required in the SW model and the modified FS model are determined, and the load spectra of variable amplitude fatigue are compiled.
(3)
The damage parameters and fatigue properties obtained from ELCF tests are substituted into the SW model and modified FS model, respectively. The fatigue life corresponding to each strain amplitude is calculated according to the load spectrum of variable amplitude fatigue.
(4)
The fatigue damage is evaluated based on the Miner linear damage accumulation criterion, and the applicability of the two damage evaluation models is analyzed by comparing the fatigue damage assessment with the test results.

3.2. Determination of Critical Plane and Cyclic Rain Flow Counting

The critical plane approach has been widely accepted in multiaxial fatigue analysis. It is the primary method for performing strain state decomposition according to six strain components in which the strain tensors on the candidate material plane are calculated, and the location of the fatigue damage critical plane by means of the weighted average method is then determined.
The z-axis is chosen parallel to the loading axis at the beam end, and the x-axis coincides with the outside normal of the surface. Then, the strain tensor at the weld root can be expressed as
ε = [ ε x x ε x y ε x z ε x y ε y y ε y z ε x z ε y z ε z z ]
A candidate material plane is characterized by the critical location with the orientation of the unit normal n to the plane defined by angles θ and ϕ [33], as shown in Figure 10. The strain tensors on the candidate material plane are calculated by computing the transformation matrices for a given set of θ and ϕ angles.
ε = M T ε M
where ε is the strain tensor on the candidate material plane and M T is the transpose of the transformation matrices M ; M is given by
Μ = [ a 11 a 12 a 13 a 21 a 22 a 23 a 31 a 32 a 33 ] = [ cos θ sin ϕ sin θ sin ϕ cos ϕ sin θ cos θ 0 cos θ cos ϕ sin θ cos ϕ sin ϕ ]
The shear strain component and normal strain component on the candidate material plane can be derived from Equations (4) and (5):
{ ε x x = 0.25 ( 1 cos 2 ϕ ) [ ε x x ( 1 + cos 2 θ ) + ε y y ( 1 cos 2 θ ) + γ x y sin 2 θ ] + 0.5 [ ε z z ( 1 + cos 2 ϕ ) + ( γ x z cos θ + γ y z sin θ ) sin 2 ϕ ] γ x y = 0.5 sin ϕ ( ε x x sin 2 θ + ε y y sin 2 θ + 2 γ x y cos 2 θ ) + cos ϕ ( ε z z sin ϕ + γ y z cos θ γ z x sin θ ) γ x z = 0.25 sin 2 ϕ [ ε x x ( 1 + cos 2 θ ) + ε y y ( cos 2 θ 1 ) + 2 ε z z 2 γ x y sin 2 θ ] cos 2 ϕ ( γ y z sin θ + γ z x cos θ )
The shear and normal strains on the candidate material plane ( θ , ϕ ) can be expressed as a function of time t as Equation (6). Then, the shear strain range acting on the candidate material plane in the current cycle can be calculated to determine the maximum shear strain range Δ γ max and the location of the Δ γ max plane. For the weld root of NWB_2, the distribution curve surfaces of the critical plane at different times are shown in Figure 11.
It can be seen from Figure 11 that the maximum shear strain range Δ γ max on the same candidate material plane changes significantly with time at 400 and 800 s. That is to say, the location of the Δ γ max plane of the fatigue failure dangerous position at the weld zone changes continuously when the beam ends are controlled by reverse cyclic loading. It is well known that normal strain and shear strain are controlled in the same phase under multiaxial proportional loading, and that the strain principal axis only changes in size but not in direction. Thus, the location of the Δ γ max plane remains unchanged at any time. However, there is a phase difference between the normal strain and the shear strain under multiaxial nonproportional loading, i.e., the strain principal axis changes not only in size but also in direction. Therefore, a new Δ γ max plane may be generated, corresponding to a new different strain state under multiaxial nonproportional loading. That is why the location of Δ γ max plane is different at any time.
In terms of fatigue damage assessment, it is necessary that the critical damage plane is representative. It is considered that the normal and shear strain on the material plane often have several different peaks during the response time history; the angle of the critical damage plane θ ¯ is determined by the weighted average method in the present study. θ ¯ is given by
θ ¯ = 1 W i = 1 n θ ( t i ) w ( t i )
where W is the sum of weights w ( t i ) and θ ( t i ) is the maximum angle of the candidate material plane in the current time. w ( t i ) is the weight function of θ ( t i ) , which indicates that the maximum shear strain is related to the damage of the material; w ( t i ) is given by
w ( t i ) = γ t , max γ min γ max γ min
where γ max , γ min is the maximum and minimum shear strain in the whole strain time history, respectively. γ t , max is the maximum shear strain at time t .
In light of the above steps, the location of the critical damage plane of the dangerous position at the weld zone of NWB_2 is calculated as θ = 86 and ϕ = 73 using MATLAB. The obtained critical plane angles are substituted into Equation (6). The maximum shear strain response on the critical plane γ c and the normal strain response on the maximum shear plane ε n , c are then determined, as shown in Figure 12.
Based on the obtained strain responses on the critical damage plane, the double rain flow counting method can be used to compile the load spectrums of variable amplitude fatigue. In practice, the number of cycles, shear strain range Δ γ , and the shear mean strain γ m are counted first, and the turning point of each cycle is also determined. The peaks of the dotted line are the turning points as shown in Figure 13a. Then, the maximum normal strain range Δ ε and the normal mean strain ε m between adjacent turning points are calculated from the normal strain response on the critical plane, as shown in Figure 13b.
The double cyclic rain flow counting results for the critical damage plane at the weld zone of NWB_2 are shown in Table 5.
As shown in Table 5, the value of 1 presented in the number of cycles means one cycle, and 0.5 means half a cycle. In terms of the extracted loading process of NWB_2, the number of rain flow counts is 14, and the total number of cycles is 10.5 weeks. The load spectra of variable amplitude fatigue obtained from the statistical cyclic counting are plotted in Figure 14.

3.3. Multiaxial Fatigue Damage Parameters

In light of the above studies, the fatigue damage parameters for fatigue life prediction are determined by the results of the cyclic rain flow counting. The damage parameters of the SW model and modified FS model are as follows:

3.3.1. Parameters of the SW model

It is indicated by the tension–torsion multiaxial fatigue test that the amplitude of ε n on the maximum shear plane is very small under proportional loading, whereas the amplitude of ε n obviously increases with the phase angle under nonproportional loading [36]. Furthermore, Shang and Wang point out that fatigue crack growth is a decohesion process along the crack tip shear band from the point of view of the micrometer, and the normal strain on the crack plane accelerates this behavior. This means that γ max and ε n on the critical plane are two important fatigue damage-controlled parameters, and the size of the normal strain excursion between adjacent turning points is one of the parameters affecting fatigue crack growth. Based on the above ideas, Shang and Wang suggest that an equivalent strain to the von Mises criterion is used as a multiaxial fatigue damage parameter:
Δ ε e q c r 2 = [ ε n * 2 + 1 3 ( Δ γ max 2 ) 2 ] 1 2
where Δ ε e q c r / 2 is the equivalent strain, ε n * is the normal strain excursion between adjacent turning points, and Δ γ max is the maximum shear strain range on the critical plane.
The multiaxial fatigue damage formula relating to the Manson–Coffin equation is expressed as Equation (10), that is, the SW model:
Δ ε e q c r 2 = σ f E ( 2 N f ) b + ε f ( 2 N f ) c
In the case of proportional loading, Equation (10) reduces the equivalent strain approach form. In the case of uniaxial loading, Equation (10) reduces the Manson–Coffin equation to a uniaxial form. Thus, Shang and Wang proposed that Equation (10) could be used as a unified fatigue damage criterion under either proportional loading (including uniaxial and multiaxial loading) or nonproportional loading.

3.3.2. Parameters of the Modified FS model

In terms of multiaxial fatigue study, researchers have made great progress in the area of damage accumulation and life prediction based on the strain components or the combined strain with stress components on the critical plane. The other popular approach is that proposed by Fatemi and Socie. They introduced a maximum normal stress σ n , max on the maximum shear strain plane to reflect the effect of additional cyclic hardening under nonproportional loading [32]. The FS damage parameter is expressed as
Δ γ max 2 ( 1 + k σ n , max σ y ) = [ ( 1 + υ e ) σ f E ( 2 N f ) b + ( 1 + υ p ) ε f ( 2 N f ) c ] ( 1 + k σ f 2 σ y ( 2 N f ) b )
where k is an empirical constant, σ y is the yield strength, and υ e , υ p are elastic Poisson ratio and plastic Poisson ratio, respectively. The other symbols are as previously stated. It is worth noting here that the stress-correlated factor, 1 2 ( 1 + k σ n , max σ y ) , found in Equation (11), is not constant. The maximum normal stress σ n , max is determined by the modified Ramberg–Osgood relation [31]:
ε n , max = σ n , max E + [ σ n , max K ( 1 + L Φ ) ] 1 n
where ε n , max is the maximum normal strain, K and n are as previously stated, and L is the nonproportional additional hardening coefficient, which depends on the material. L is determined by [37]
L = σ ¯ OP σ ¯ IP 1
where σ ¯ OP , σ ¯ IP are the equivalent stress amplitude under 90° out-of-phase loading and the equivalent stress amplitude under proportional loading, respectively. Φ is the nonproportionality of the strain path determined by [31]
Φ = 2 A θ , max A max 1
where A θ , max is the swept area of the Δ γ θ , max θ polar coordinate space and A max is the circle area with a radius of the maximum shear strain during one cycle.
It is believed that Equation (12) can account for the nonproportionality of both material additional hardening and the strain path. However, it is worth noting that not all metallic materials exhibit nonproportional additional hardening characteristics under nonproportional loading. The value of L is zero for materials without a nonproportional additional hardening effect, and the value of the modified coefficient ( 1 + L Φ ) in Equation (12) is then zero, even under nonproportional loading. This means that Equation (12) is similar to the Ramberg–Osgood relation, i.e., Equation (2). That is to say, Equation (12) cannot accurately reflect the effect of additional hardening when applied to fatigue life prediction for a material without a nonproportional additional hardening effect. To solve this problem, a new modified coefficient k * is introduced in Equation (11) to reflect the effect of additional hardening by our research team. In addition, it is suggested that σ n , max can be replaced by σ ¯ n , max . Therefore, our previously proposed modified FS model is expressed as
Δ γ max 2 ( 1 + k * σ ¯ n , max σ y ) = [ ( 1 + υ e ) σ f E ( 2 N f ) b + ( 1 + υ p ) ε f ( 2 N f ) c ] ( 1 + k * σ f 2 σ y ( 2 N f ) b )
where σ ¯ n , max is the maximum normal stress on the critical plane determined by the Ramberg–Osgood relation. The other symbols are as previously stated except k * , which is defined as the nonproportional factor, and is given by [35]
k * = 1 + Φ 2 + ( L Φ ) 2 = 1 + Φ 1 + L 2
It can be observed from Equation (16) that the value of L is zero for materials without a nonproportional additional hardening effect, whereas the effect of nonproportionality of the strain path is not ignored. That is to say, k * can account for the effect of the nonproportionality of both the strain path and material additional hardening. In other words, the proposed modified FS model for life prediction appears to be relatively accurate in adopting the modified coefficient k * , and may be used for different materials.
The nonproportionality of the strain path cannot be adequately accounted for if determined by Equation (14). This is because the value of Φ may be negative or zero when A θ , max is less than or equal to A max / 2 , resulting in limitations of application. To solve this problem, we carry out a second amendment for the nonproportional factor k * . Φ is determined by the introduction of the moment of inertia [34]. As shown in Figure 15, L ¯ is the random strain path simplified by its enveloping line, i.e., the equivalent convex path L , and the maximum chord length of L ¯ is the diameter of the maximum circle path L o . It is worth noting that the midpoint of the maximum chord length is used as the origin of the ε γ / 3 coordinate system, the x-axis represents the in-phase proportional loading path, and the y-axis represents the antiphase proportional loading path.
Φ is calculated using the following formula:
Φ = ( I / I o ) h
where I , I o , and h are respectively given by
I = S y 2 d S
I o = S o y 2 d S o
h = ( 1 S ¯ S o ) L c y c 4 Δ r max
where I is the moment of inertia of L to the x-axis in Equation (18), which indicates the deviation of the point from the in-phase proportional loading path. In other words, the moment of inertia can account for the nonproportional contribution of the strain path. Any point in area S surrounded by L contributes to nonproportional additional hardening, and y represents the distance from any point in area S to the x-axis. Similar to I , I o is the moment of inertia of L o to the x-axis. S o is the circle area with diameter of the maximum chord length. y is used to express the distance from any point in area S o to the x-axis in Equation (19). Because an equivalent convex path may appear even under different random strain path loading, the contribution of the random strain path is taken into account by the parameter h . S ¯ is the area surrounded by the random strain path L ¯ . L c y c and Δ r max are the perimeter and the maximum strain range of L ¯ , respectively.
The strain path of the dangerous position at the weld zone of NWB_2 is shown in Figure 16. The dotted line represents the original coordinate position and the new coordinate position is determined by the midpoint of the maximum chord length of the strain path marked with the blue line. Its equivalent convex path and maximum circular path are marked with the red and black lines, respectively.
It can be seen from Figure 16 that the strain path changes constantly and has obvious randomness, while the determination of Φ is convenient for engineering applications when the convex hull in the MATLAB program is adopted, and its physical significance is clear.
In light of the above study, the nonproportional parameters in the modified FS model are shown in Table 6.

3.4. Damage Accumulation Approach

Fatigue damage is an irreversible process in which material properties deteriorate continuously under reverse cyclic loading. It is also a process in which damage accumulates gradually. Generally speaking, most of the components are subjected to variable amplitude cyclic loading. Therefore, it is very important for fatigue damage evaluations to study the damage accumulation rule and fatigue failure criteria under variable amplitude cyclic loading.
It is known that the Miner linear damage accumulation rule can predict the mean life of engineering structures under random loading. The fatigue damage per cycle, D = 1 / N f , can be determined from the fatigue criteria for each parameter. Therefore, the Miner linear damage accumulation rule is used to calculate the cumulative fatigue damage in this paper, and the damage value D = 1 expresses fatigue crack initiation.

4. Discussion

The obtained fatigue damage parameters are respectively substituted into Equations (10) and (15) to calculate the fatigue lives N i for each strain amplitude. Then, the fatigue damage per cycle is determined. The calculated results of cumulative damage ranging from the 3rd to the 8th cycle are shown in Table 7.
As shown in Table 7, it can be seen that:
(1)
In terms of the same fatigue properties of either CLG or PB, ranging from the 3rd to the 8th cycle, the damage values of the SW model are much smaller than those of the modified FS model. From the comparison between the experiment and damage evaluation, it is indicated that for box-column joints, the SW model is inclined to underestimate the risk. This may be because the nonproportionality of material additional hardening and the strain path both aggravate fatigue damage under nonproportional loading [31], whereas the nonproportional factor is not taken into account in the SW model.
(2)
In terms of the same fatigue damage model, either SW or the modified FS model, the fatigue values calculated by the fatigue properties of CLG are much larger than those calculated by the fatigue properties of PB, which are more coincident with the reverse cyclic loading test results. This may be because the geometric shape and force characteristics of the CLG specimens are more similar to those of box-column joints than those of PB specimens. It is shown that the structural details, i.e., geometric shape and force characteristics, play an important role in damage evaluation; the more similar the structural details, the better the damage estimation.
(3)
In terms of the modified FS model using the fatigue properties of CLG, the calculated damage value of the dangerous position at the weld zone of NWB_2 is 1.0321 at the 7th cycle, which is slightly larger than D = 1. The Miner linear damage accumulation rule indicates that a microscopic crack will be initiated at the time when the 7th cycle is coming to an end, while it was previously stated in Section 2.2 that tiny visible crack initiation could be observed just after the 7th cycle, that is, at the time period between the 7th and 8th cycle. However, there is no contradiction between the theoretical calculation results and the experimental results. There are two reasons for this. The first is that from the point of view of a micrometer, fatigue crack is a process ranging from microcosmic to visible, which takes some time. Furthermore, the observation of a crack is limited by the test instruments. The second reason is that the error between the theoretical calculation results and the experimental results is very small, which is acceptable. Therefore, it is shown that the damage assessment from adopting this model is relatively accurate.
In summary, for ELCF damage evaluations of box-column welded joints, the method can provide good estimations in cases where the fatigue properties of CLG are adopted as the proposed modified FS model parameters.

5. Conclusions

Based on numerous existing studies, it is difficult to accurately assess the damage of a steel structure during the ELCF fracture process using the macroparameters of the whole structure or a structure component, since local structural details and fatigue crack initiation mechanics are ignored. A damage evaluation methodology based on the structural details of I-beam to box-column welded joints, extended from multiaxial fatigue critical plane theory, is proposed to assess ELCF crack initiation life and damage. Based on experimental studies associated with analytical studies, the key conclusions are as follows:
(1)
The failure of beam-to-column welded joints is multiaxial and nonproportional ELCF when subjected to severe seismic excitation.
(2)
From a comparison between the FS model containing a nonproportional factor and the SW model without this factor, it was shown that the effect of nonproportional additional hardening plays a key role in damage evaluation.
(3)
The nonproportional factor in the FS model previously proposed by us was again modified and introduced to reflect the effect of nonproportionality. This accounts for both the strain path and additional material hardening. It was found that the fatigue damage model including this proposed factor can successfully evaluate the fatigue damage of beam-to-column welded joints with reasonable accuracy.
(4)
The fatigue properties of PB and CLG when considering welded structural details were respectively adopted as the damage parameters in the modified FS model and the SW model. From the comparison between experiments and damage estimation, it was shown that the fatigue properties of CLG result in the damage evaluation being more consistent with the test results than that of PB because of the highly similar geometric shape and force characteristics.
In conclusion, the modified FS model with the fatigue properties of CLG is most in line with the test results, which verifies the applicability of the model. The presented work makes an important contribution to ELCF damage evaluations of steel frameworks.

Author Contributions

Formal analysis, E.Z.; Funding acquisition, L.H. and W.Q.; Investigation, L.H.; Methodology, L.H., W.Q. and E.Z.; Software, L.H. and E.Z.; Writing—original draft, L.H.; Writing—review & editing, L.H. and W.Q. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by [Science and Technology Research Project of Hubei Provincial Department of Education] grant number [Q20192704], [National Natural Science Foundation of China] grant number [51438002, 51378409], [Hubei Key Laboratory of Roadway Bridge and Structure Engineering (Wuhan University of Technology)] grant number [DQJJ201907]. And the APC was funded by [Q20192704].

Acknowledgments

The corresponding author also sincerely thanks Ernian Zhao for his guidance and cooperation.

Conflicts of Interest

No conflict of interest exits in the submission of this manuscript, and manuscript is approved by all authors for publication.

References

  1. Kanvinde, A.M.; Deierlein, G.G. Cyclic void growth model to assess ductile fracture initiation in structural steels due to ultra low cycle fatigue. J. Eng. Mech. 2007, 133, 701–712. [Google Scholar] [CrossRef]
  2. Mendes, L.A.; Castro, L.M. A simplified reinforcing steel model suitable for cyclic loading including ultra-low-cycle fatigue effects. Eng. Struct. 2014, 68, 155–164. [Google Scholar] [CrossRef]
  3. Nip, K.H.; Gardner, L.; Davies, C.M.; Elghazouli, A.Y. Extremely low cycle fatigue tests on structural carbon steel and stainless steel. J. Constr. Steel Res. 2010, 66, 96–110. [Google Scholar] [CrossRef]
  4. Miller, D.K. Lessons learned from the Northridge earthquake. Eng. Struct. 1998, 20, 249–260. [Google Scholar] [CrossRef]
  5. Nakashima, M.; Inoue, K.; Tada, M. Classification of damage to steel buildings observed in the 1995 Hyogoken-Nanbu earthquake. Eng. Struct. 1998, 20, 271–281. [Google Scholar] [CrossRef]
  6. Krawinkler, H.; Zohrei, M. Cumulative damage in steel structures subjected to earthquake ground motion. Comput. Struct. 1983, 16, 531–541. [Google Scholar] [CrossRef]
  7. Bernuzzi, C.; Calado, L.; Castiglioni, C.A. Ductility and load carrying capacity prediction of steel beam-to-column connections under cyclic reversal loading. J. Earthq. Eng. 1997, 1, 401–432. [Google Scholar] [CrossRef]
  8. Ballio, G.; Castiglioni, C.A. A unified approach for the design of steel structures under low or high cycle fatigue. J. Constr. Steel Res. 1995, 34, 75–101. [Google Scholar] [CrossRef]
  9. Kishiki, S.; Lee, D.S.; Yamada, S.; Ishida, T.; Jiao, Y. Low-cycle fatigue performance assessment of current Japanese steel beam-to-column connections determined by ductile fracture. Eng. Struct. 2019, 182, 241–250. [Google Scholar] [CrossRef]
  10. Kanno, R. Advances in steel materials for innovative and elegant steel structures in Japan—A. review. Struct. Eng. Int. 2016, 26, 242–253. [Google Scholar] [CrossRef]
  11. Powell, G.H.; Allahabadi, R. Seismic damage prediction by deterministic methods: Concepts and procedures. Earthq. Eng. Struct. Dyn. 1988, 16, 719–734. [Google Scholar] [CrossRef]
  12. Shi, Y.; Xiong, J.; Wang, Y.Q. Study on damage behavior of weld of beam-to-column connection in steel frameⅠ: Experiment. J. Build. Struct. 2012, 33, 48–55. [Google Scholar] [CrossRef]
  13. Darwin, D.; Nmai, C.K. Energy dissipation in RC beams under cyclic load. J. Struct. Eng. 1986, 112, 1829–1846. [Google Scholar] [CrossRef]
  14. Castiglioni, C.A.; Pucinotti, R. Failure criteria and cumulative damage models for steel components under cyclic loading. J. Constr. Steel Res. 2009, 65, 751–765. [Google Scholar] [CrossRef]
  15. Xu, L.H.; Shan, X.; Li, Z.X. Vulnerability analysis for steel frame structure due to strong earthquakes. Appl. Mech. Mater. 2011, 90–93, 1486–1489. [Google Scholar] [CrossRef]
  16. Castillo, E.; Fernández-Canteli, A. A Unified Statistical Methodology for Modeling Fatigue Damage; Springer: Dordrecht, The Netherlands, 2009. [Google Scholar] [CrossRef]
  17. Correia, J.A.F.O.; Raposo, P.; Muniz-Calvente, M.; Blasón, S.; Lesiuk, G.; De Jesus, A.M.P.; Moreira, P.M.G.P.; Calçada, R.A.B.; Canteli, A.F. A generalization of the fatigue kohout-věchet model for several fatigue damage parameters. Eng. Fract. Mech. 2017, 185, 284–300. [Google Scholar] [CrossRef]
  18. Park, Y.J.; Ang, A.H.S.; Wen, Y.K. Seismic damage analysis of reinforced concrete buildings. J. Struct. Eng. 1985, 111, 740–757. [Google Scholar] [CrossRef]
  19. Park, Y.J.; Ang, A.H.S. Mechanistic seismic damage model for reinforced concrete. J. Struct. Eng. 1985, 111, 722–739. [Google Scholar] [CrossRef]
  20. Wang, J.W.; Qu, W.L. Failure model of steel members based on deformation and energy dissipation. J. Vib. Shock 2013, 32, 71–75. [Google Scholar] [CrossRef]
  21. Dong, J.; Ma, H.; Zhang, N.; Liu, Y.; Mao, Z. Seismic damage assessment of steel reinforced recycled concrete column-steel beam composite frame joints. Earthq. Struct. 2018, 14, 73–84. [Google Scholar] [CrossRef]
  22. Zhang, Y.T.; Du, X.J.; Yang, L. Research on seismic performance assessment based on component damage for RC frame structure. J. Hunan Univ. Nat. Sci. 2016, 43, 9–21. [Google Scholar]
  23. Atzori, B.; Lazzarin, P.; Meneghetti, G. Fatigue strength of welded joints based on local, semi-local and nominal approaches. Theor. Appl. Fract. Mech. 2009, 52, 55–61. [Google Scholar] [CrossRef]
  24. Gholami, M.; Deylami, A.; Tehranizadeh, M. Seismic performance of flange plate connections between steel beams and box columns. J. Constr. Steel Res. 2013, 84, 36–48. [Google Scholar] [CrossRef]
  25. Kermajani, M.; Malek Ghaini, F.; Miresmaeili, R.; Aghakouchak, A.A.; Shadmand, M. Effect of weld metal toughness on fracture behavior under ultra-low cycle fatigue loading (earthquake). Mater. Sci. Eng. A 2016, 668, 30–37. [Google Scholar] [CrossRef]
  26. Tong, L.; Huang, X.; Zhou, F.; Chen, Y. Experimental and numerical investigations on extremely-low-cycle fatigue fracture behavior of steel welded joints. J. Constr. Steel Res. 2016, 119, 98–112. [Google Scholar] [CrossRef]
  27. Huang, L.Z.; Qu, W.L. Experimental study on low cycle fatigue for structural details of welded joints. J. Huazhong Univ. Sci. Technol. Nat. Sci. Ed. 2016, 44, 36–40. [Google Scholar] [CrossRef]
  28. Crupi, G.; Crupi, V.; Guglielmino, E.; Taylor, D. Fatigue assessment of welded joints using critical distance and other methods. Eng. Fail. Anal. 2005, 12, 129–142. [Google Scholar] [CrossRef]
  29. Čermelj, B.; Može, P.; Sinur, F. On the prediction of low-cycle fatigue in steel welded beam-to-column joints. J. Constr. Steel Res. 2016, 117, 49–63. [Google Scholar] [CrossRef]
  30. Huang, X.W.; Zhao, J. A cumulative damage model for extremely low cycle fatigue cracking in steel structure. Struct. Eng. Mech. 2017, 62, 225–236. [Google Scholar] [CrossRef]
  31. Chen, X.; Gao, Q.; Sun, X.F. Low-cycle fatigue under non-proportional loading. Fatigue Fract. Eng. Mater. Struct. 1996, 19, 839–854. [Google Scholar] [CrossRef]
  32. Fatemi, A.; Socie, D.F. A critical plane approach to multiaxial fatigue damage including out-of-phase loading. Fatigue Fract. Eng. Mater. Struct. 1988, 11, 149–165. [Google Scholar] [CrossRef]
  33. Li, J.; Zhang, Z.P.; Sun, Q.; Li, C.W. Multiaxial fatigue life prediction for various metallic materials based on the critical plane approach. Int. J. Fatigue 2011, 33, 90–101. [Google Scholar] [CrossRef]
  34. Zhong, B.; Wang, Y.R.; Wei, D.S.; Yang, L. Multiaxial fatigue life prediction based on non-proportionality of strain path. J. Aerosp. Power 2016, 31, 317–322. [Google Scholar] [CrossRef]
  35. Zhao, E.N.; Qu, W.L. A new proposal for multiaxial low-cycle fatigue life prediction under non-proportional loading. Chin. J. Theor. Appl. Mech. 2016, 37, 153–160. [Google Scholar] [CrossRef]
  36. Shang, D.G.; Wang, D.J. A new multiaxial fatigue damage model based on the critical plane approach. Int. J. Fatigue 1998, 20, 241–245. [Google Scholar] [CrossRef]
  37. Shamsaei, N.; Fatemi, A.; Socie, D.F. Multiaxial fatigue evaluation using discriminating strain paths. Int. J. Fatigue 2011, 33, 597–609. [Google Scholar] [CrossRef]
Figure 1. Design of welded specimens (mm): (a) plate butt (PB); (b) cruciform load-carrying groove (CLG).
Figure 1. Design of welded specimens (mm): (a) plate butt (PB); (b) cruciform load-carrying groove (CLG).
Materials 13 01768 g001
Figure 2. (a)Test equipment; (b) loading history.
Figure 2. (a)Test equipment; (b) loading history.
Materials 13 01768 g002
Figure 3. Cyclic stress–strain hysteresis curve.
Figure 3. Cyclic stress–strain hysteresis curve.
Materials 13 01768 g003
Figure 4. Design of specimen (mm): (a) beam-to-column welded joint; (b) welded structural details.
Figure 4. Design of specimen (mm): (a) beam-to-column welded joint; (b) welded structural details.
Materials 13 01768 g004
Figure 5. Test setup configuration.
Figure 5. Test setup configuration.
Materials 13 01768 g005
Figure 6. von Mises stress distribution nephogram.
Figure 6. von Mises stress distribution nephogram.
Materials 13 01768 g006
Figure 7. The location of strain rosettes in the specimen.
Figure 7. The location of strain rosettes in the specimen.
Materials 13 01768 g007
Figure 8. Six strain responses of NWB_2 weld. (* in the figure means the peak or valley of the strain).
Figure 8. Six strain responses of NWB_2 weld. (* in the figure means the peak or valley of the strain).
Materials 13 01768 g008
Figure 9. Extremely-low-cycle fatigue (ELCF) damage analysis steps.
Figure 9. Extremely-low-cycle fatigue (ELCF) damage analysis steps.
Materials 13 01768 g009
Figure 10. The strains on the candidate critical plane in a three-dimensional coordinate system.
Figure 10. The strains on the candidate critical plane in a three-dimensional coordinate system.
Materials 13 01768 g010
Figure 11. The location of the critical plane: (a) time = 400 s; (b) time = 800 s.
Figure 11. The location of the critical plane: (a) time = 400 s; (b) time = 800 s.
Materials 13 01768 g011
Figure 12. Strain response on the critical damage plane: (a) shear strain; (b) normal strain.
Figure 12. Strain response on the critical damage plane: (a) shear strain; (b) normal strain.
Materials 13 01768 g012
Figure 13. Double cyclic rain flow counting: (a) shear strain range; (b) normal strain range. (Red line represents strain responses; blue line represents strain range correlated to double rain flow counting method).
Figure 13. Double cyclic rain flow counting: (a) shear strain range; (b) normal strain range. (Red line represents strain responses; blue line represents strain range correlated to double rain flow counting method).
Materials 13 01768 g013
Figure 14. The load spectra of variable amplitude fatigue: (a) shear strain range; (b) normal strain range.
Figure 14. The load spectra of variable amplitude fatigue: (a) shear strain range; (b) normal strain range.
Materials 13 01768 g014
Figure 15. Random strain path and its envelope path.
Figure 15. Random strain path and its envelope path.
Materials 13 01768 g015
Figure 16. Strain path and equivalent convex path on the critical plane of NWB_2.
Figure 16. Strain path and equivalent convex path on the critical plane of NWB_2.
Materials 13 01768 g016
Table 1. Relevant mechanical parameters.
Table 1. Relevant mechanical parameters.
Young’s Modulus E ( GPa ) Yield Strength σ y ( MPa ) Ultimate Strength σ u ( MPa ) Ultimate Strain ε u ( % ) Elastic Poisson Ratio υ e Plastic Poisson Ratio υ p
200254.522415.44920.9760.30.5
Table 2. Strain amplitude used in the tests and test results.
Table 2. Strain amplitude used in the tests and test results.
Specimen No. Δ ε 2 Δ ε p 2 Δ ε e 2 Δ σ / 2 ( MPa ) 2 N f Specimen No. Δ ε 2 Δ ε p 2 Δ ε e 2 Δ σ / 2 ( MPa ) 2 N f
PB040.0040.00240.0016319464CLG040.0040.00240.0016316810
PB050.0050.00320.0018355210CLG050.0050.00330.0017329460
PB060.0060.00420.0018356140CLG060.0060.00430.0017353164
PB070.0070.00500.002036498CLG070.0070.00520.001840384
PB080.0080.00600.002039580CLG080.0080.00600.002038966
PB100.0100.00780.002241562CLG100.0100.00790.002141930
Table 3. Fatigue and cyclic stress–strain properties of welded structural details.
Table 3. Fatigue and cyclic stress–strain properties of welded structural details.
Specimen K ( MPa ) n σ f ( MPa ) b ε f c
PB11340.2085783.3−0.13970.1193−0.6739
CLG13690.2430576.24−0.08300.0273−0.3611
Table 4. Loading condition.
Table 4. Loading condition.
Specimen No.Displacement of Beam Ends (ΔU/mm)
1± 15 → ± 17 → ± 20 → ± 21 → ± 22 → ± 23 → ± 25 → ± 27
2± 17
3± 14 → ± 16 → ± 18 → ± 20
Table 5. Double cyclic rain flow counting results of the maximum damage critical plane of NWB_2.
Table 5. Double cyclic rain flow counting results of the maximum damage critical plane of NWB_2.
No.Number of Cycles γ m ( % ) Δ γ ( % ) ε m ( % ) Δ ε ( % )
10.50.3110260.6224080.1023180.217569
20.50.01498551.21449−0.074480.465991
30.50.329911.844340.2209150.636
40.5−0.55573.615560.015090.997263
51−0.609093.235080.035121.144744
61−0.6623283.102540.1416590.855743
711.092272.24082−0.035510.771647
811.205262.64228−0.007061.072557
911.070133.19876−0.067970.775517
100.50.3905355.508030.0320781.150827
110.51.374653.5398−0.105870.693131
120.50.6173752.025250.0762090.400422
130.51.056331.14734−0.039330.594464
140.51.007711.05010.2216890.072419
Table 6. Nonproportional parameters in the modified FS model.
Table 6. Nonproportional parameters in the modified FS model.
Nonproportionality of Strain Path Φ Nonproportional Additional Hardening Coefficient L Nonproportional Factor k *
0.44250.311.4425
Table 7. Comparison of fatigue damage models using fatigue properties of structural details for NWB_2.
Table 7. Comparison of fatigue damage models using fatigue properties of structural details for NWB_2.
Fatigue Damage ModelFatigue PropertiesCumulative Damage
3rd4th5th6th7th8th
SW modelCLG0.14640.20890.23270.28270.34520.5535
PB0.07510.10960.13000.16030.19480.2564
Modified FS modelCLG0.40590.60590.68920.83211.03211.3821
PB0.16580.24270.29030.35700.43390.5794

Share and Cite

MDPI and ACS Style

Huang, L.; Qu, W.; Zhao, E. Extremely-Low-Cycle Fatigue Damage for Beam-to-Column Welded Joints Using Structural Details. Materials 2020, 13, 1768. https://doi.org/10.3390/ma13071768

AMA Style

Huang L, Qu W, Zhao E. Extremely-Low-Cycle Fatigue Damage for Beam-to-Column Welded Joints Using Structural Details. Materials. 2020; 13(7):1768. https://doi.org/10.3390/ma13071768

Chicago/Turabian Style

Huang, Lizhen, Weilian Qu, and Ernian Zhao. 2020. "Extremely-Low-Cycle Fatigue Damage for Beam-to-Column Welded Joints Using Structural Details" Materials 13, no. 7: 1768. https://doi.org/10.3390/ma13071768

APA Style

Huang, L., Qu, W., & Zhao, E. (2020). Extremely-Low-Cycle Fatigue Damage for Beam-to-Column Welded Joints Using Structural Details. Materials, 13(7), 1768. https://doi.org/10.3390/ma13071768

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