This paper aims to study the hysteretic characteristics of the steel plates artificially corroded by neutral salt spray. Salt spray was applied to accelerate the corrosion on the steel plates; specimens of varying degrees of corrosion were obtained in this manner. And each specimen was subject to cyclic loading test to get the hysteretic curve. Then the experimental results were extensively discussed, focusing on strength and ductility, hysteretic energy, the skeleton curve, and unloading and loading curve. After that, the hysteretic constitutive model of corroded steel was established based on the first time loading criterion, unloading criterion, cycle skeleton criterion, and reloading curve criterion. The result of the experiment showed that, with the increase of the degree of corrosion, the mechanical properties and seismic energy dissipation performance of seismic energy of the steel decreased; the deterioration of ductility got aggravated. On the other hand, the skeleton curve and the Ramberg-Osgood model were well matched, and the coefficient of circular enhancement showed a decreasing trend; the variation of cyclic hardening exponent did not have an obvious pattern. Meanwhile, the hysteretic constitutive model of corroded steel and the results of the experiment were well matched.
1. Introduction
Steel structure is widely used in large bridges, offshore oil platforms, industrial constructions, and residential buildings [1–3], as it has excellent mechanical properties and it is not expensive with good performance. Particularly its good ductility makes it easier to be the first choice for seismic design. However, if it is exposed to corrosive environment for a long time such as soil, air, acid rain, and seawater, it is easy to be corroded [4–7].
In order to study the phenomenon of the metal corrosion while being exposed in the real circumstance, some test methods are adopted to accelerate corrosion, such as exposing metal in hot and humid environment, using salt spray, using salt composite spray, and alternately exposing metal in wet and dry environments. Research had shown that [8–10] the accelerated corrosion test could represent the long-term natural exposure test. Moreover, it has shorter test cycle, and the test result can be reproduced. Nonetheless, it is difficult to simulate the natural environment completely due to the relatively big difference between the natural environment and the test conditions.
Corrosion will cause the mechanical properties of materials to degrade. Many scholars, who studied the deterioration of mechanical properties of corroded steel [11–15], had found that not only would the cross section of corroded steel reduce and the surface features change, but also the yield strength, the ultimate strength, and the ductility of corroded steel would reduce as well; and among them, the decrease of ductility was the most obvious one. In [16], author cut one corroded H shape steel and one noncorroded H shape steel which were actual engineering components and studied the hysteretic behaviors of bending under low cyclic loading. The results showed that the corroded specimen cracked faster than the noncorroded one, and the corroded specimen had fewer hysteretic cycles than the noncorroded one under the same loading system. In [17], the formulation and verification of a cyclic stress-strain relationship of reinforcing bars were presented. The cyclic loops followed Giuffré-Menegotto-Pinto [18] equations with some modifications to account for the effect of buckling. A complete path-dependent cyclic constitutive model was then obtained by combining the equations of the two monotonic envelopes and the cyclic loops. Mendes and Castro [19] presented a new constitutive model for the simulation of reinforcing steel bars used in common reinforced concrete structures and it was designed to be used for general loading cases. The model included the well-known Guiffré-Menegotto-Pinto [18] softened branch, although new expressions were proposed for the evolutions of the curvature-related parameter and of the yield surface. Shi et al. [20] analyzed cyclic performance of high-strength structural steel and established an appropriate constitutive relationship. The study showed that the Q460D was similar to the ordinary strength steel [21] in some aspects, such as plasticity, cyclic hardening or softening, and Bauschinger effect [22]. Both experimental and modeling results showed that the responses of high-strength structural steel under cyclic loading and monotonic loading were different. The necking and fracture behavior would occur in advance for the former because the accumulated damage reduced the ductility of steel.
However, studying the damage accumulation of corroded steel and the cyclic plastic constitutive model of corroded steel under low cyclic loading are not covered in other researches. The hysteretic loop of steel under low cyclic loading belongs to low-cycle fatigue, and corrosion has a great and obvious effect on plastic deformation, crack growth, and fracture toughness of steel under low cyclic loading. Therefore, in order to study the hysteretic behaviors of the steel with different degrees of corrosion in the neutral salt spray environment and in order to study the cyclic plastic constitutive model, this paper will carry out the following experiment method. Firstly, the Q235 steel plate would undergo the process of accelerated corrosion. Secondly, cyclic loading tests would be carried out to analyze the hysteretic behavior of steel. Thirdly, the Ramberg-Osgood model [22] would be applied to compare the cyclic skeleton curves of steel in different corrosion conditions, and the corresponding parameters would be analyzed. Finally, a hysteretic constitutive model would be established, which could provide basic data for the study of the seismic behavior of steel structures.
Neutral salt spray accelerated corrosion test was conducted to obtain the corroded specimens, based on the standard of GB/T 10125-2012 [23] and GB/T 24517-2009 [24], where chemical reagents (NaCl) and distilled water were adopted, the NaCl solution concentration was 50 mg/L, the pH value was 6.2–7.2, and the diameter of the nozzle and the rate of mist spray of the spraying equipment were 0.5 mm–1.5 mm and 0.5–1.5 L/(min·m2), respectively. There were eight sets of specimens named A01–A71 (280 mm × 50 mm × 8 mm [25]). All specimens for corrosion were placed individually on an exposure test fixture which are parallel at an angle of 45° with regard to the vertical direction. The temperature and humidity were the same as the natural environment. In order to ensure the uniform corrosion of both sides of the test pieces, salt spray would be used on both sides of the test pieces every 2-3 days. All sample times (corrosion time): A11: 30 days, A21: 110 days, A31: 150 days, A41: 250 days, A51: 310 days, A61: 440 days. The specimens were retrieved and cleared away from the corrosion products with dilute hydrochloric acid and distilled water and then kept in CaO desiccator until the tensile test. In order to achieve the mass-loss rates of the corroded specimens, they were weighed with analytical balance before and after the accelerated corrosion test.
2.2. Cyclic Load Experiment
The specific sizes are shown in Figure 1(a). Loading device is the INSTRON Model 1343 which is a universal tension, compression, and torsion fatigue testing machine (Figure 1(c)). The displacement gauge distance was 50 mm, and its measuring range was 50% in the direction of tension and 10% in the pressure direction. The antibuckling device was designed to prevent the lateral buckling deformation of the plate under the action of larger pressure strain. The strain control was used to carry out large strain cyclic loading. Loading cycle was two times per level. Loading was applied with a 0.3% increase in amplitude. All of the cyclic loading waveforms were triangular wave and the strain loading rate was 0.00025/s. When the pressure strain reached 2.1%, the cyclic load was no longer applied, and then the specimen was subject to the uniaxial tensile load until the specimen cracked. The cyclic loading system is shown in Figure 1(b).
Size of experimental components, device, and system for loading.
3. Experiment Result3.1. Hysteretic Curve
It is difficult to measure the real strength of steel due to the uneven surface caused by corrosion. Therefore, the nominal stress (namely, the ratio of actual load to design area) was adopted to objectively reflect the effects of corrosion on the hysteretic behavior of steel. The hysteretic loading curves of all specimens are shown in Figure 2.
Hysteretic curve of specimen.
Figure 2 shows that steel hysteresis curves are divided into three parts including skeleton curve, unloading curve, and softening curve (reloading curve). The irregular changes in individual curves may be related to the antibuckling device and the displacement meter. The first loading curve of the steel extends along the monotonic curve, which is uniaxial tension effect. The unloading of the steel is in accordance with the elastic linear unloading, and the unloading stiffness is the same as that of the initial elastic modulus Es, but with the increase of the strain, the elastic modulus of the unloading curve at all levels slightly decreases. The loading curve of the steel is the same as that of the loading curve of the same direction, which is in accordance with the model of the peak value. The peak stress of cyclic loading at the next stage of the steel has a greater dependence on the upper level. Steel was pulled into the plastic deformation stage, and stress is unloaded to zero. And then steel will undergo reverse loading to yield. The new yield stress decreases, which shows Bauschinger effect [22]. With the increase of strain, the stress of the steel increases, the effect of cyclic hardening occurs, and the strengthening effect of the steel is more obvious in the pulling direction than that in the pressing direction. Under the same strain amplitude, the hysteretic curves are not exactly the same as those in the target strain, and the stress in the last lap is greater than that in the previous lap. All hysteresis loops are plump, which shows that the corroded steel can also achieve full hysteresis loops, indicating seismic performance and energy dissipation capacity.
The difference between monotonic curve and hysteresis curve is shown in Figure 3. In cyclic loading scenario, the hardening or softening effect caused the change of stress-strain relationship until its stabilization. It shows that steel hysteresis and monotonic constitutive curves were quite different; for example, hardening phenomenon is quite obvious and the damage accumulation due to cyclic effect always leads the steel ductility to deteriorate, and the rate of stress decrease after peak point was significantly increased [21].
Comparison chart between monotonous and hysteresis curves.
The main parameters of mechanical performance are shown in Table 1. The yielding effect of the steel under cyclic loading is not obvious, and so the corresponding strain of the steel is not listed. The ultimate strength is named fu; the ultimate strain is named εu; and the breaking elongation is named δ. The hysteretic energy is obtained by calculating the area of the hysteresis loop. It reflects the seismic energy dissipation capacity of steel in the material level.
Summary on mechanical properties of the specimen.
Serial number
Corrosion rate/%
fu/MPa
εu
δ
Hysteresis laps
Hysteretic energy/MPa
A01
0
382.5
0.173
0.298
14
253.737
A11
5.08
374.75
0.171
0.278
14
236.679
A21
7.49
360
0.169
0.267
14
226.987
A31
7.89
364.75
0.167
0.265
14
226.391
A41
9.40
354.25
0.164
0.254
14
215.232
A51
10.47
356.95
0.166
0.256
14
225.762
A61
14.80
328.5
0.163
0.233
14
208.254
A71
18.82
306.85
0.149
0.208
14
207.808
Table 1 shows that, with the increase of the degree of corrosion, the ultimate strength of steel under cyclic loading decreases. The strain corresponding to the ultimate strength and the elongation also shows a decreasing trend. During the process of the test, some of the corroded specimens cracked when they reached the maximum strength, which showed that the ductility of corroded steel was severely degraded under cyclic loading. The changes of hysteresis energy with different rates of weight loss are shown in Figure 4.
Hysteretic energy change with weight loss rate.
Compare the rate of corrosion with the hysteretic energy:(1)J=64.499×exp-ρw13.541+189.905,where J is the hysteretic energy and ρw is the weight loss rate.
It can be seen that the dissipation performance of seismic energy decreases exponentially with the increase of the degree of corrosion under the same amount of cycles.
3.2. Skeleton Curve
The cyclic stress-strain skeleton curves of steel with different degrees of weight loss are shown in Figure 5. It can be seen that the steel has an obvious cyclic hardening behavior under cyclic loading but does not have a yield platform. Compared with the noncorroded steel, the strength of corroded steel appears to be lower. And the greater the degree of corrosion is, the more obvious the decrease of the strength becomes.
Hysteresis curve on varying degrees of corrosion steel.
The maximum normal stress and normal strain obtained from the results of the experiment as well as the Ramberg-Osgood model [22] were used to be compared with the hysteretic skeleton curve of steel with different degrees of corrosion. Mathematical expressions of Ramberg-Osgood model [22] are as follows:(2)Δε2=Δεe2+Δεp2=Δσ2Es+Δσ2K′1/n′,where Δε is the total strain amplitudes; Δεe is the elastic strain amplitude; Δεp is the shaping strain amplitude; E is the elastic modulus of steel; K′ is the cyclic hardening coefficient; n′ is the cyclic hardening index; Δσ is the stable stress amplitude.
In order to make it easier to calculate and analyze, the mathematical expressions were simplified as follows:(3)ε=σEs+σK′1/n′.
Formula (3) was used to fit the cyclic loading skeleton curve. See Table 2 for the values of the cyclic hardening coefficient K′ and of the cyclic hardening index n′. The fitting curve is shown in Figure 6. Figure 6 shows that Ramberg-Osgood model [22] fitted well with test data.
Skeleton curve model.
Number
A01
A11
A21
A31
A41
A51
A61
A71
Es (105)/MPa
1.970
1.856
1.801
1.792
1.758
1.734
1.636
1.545
K′/MPa
816.813
626.500
592.927
584.254
551.420
567.530
545.944
534.919
n′
0.211
0.164
0.155
0.153
0.153
0.157
0.170
0.201
Fitting curve of skeleton model.
3.3. Parametric Degradation of Cyclic Skeleton Curve
Formula (3) shows that the right side of the model is composed of two parts. The first one shows the elastic characteristics of steel, which is linear expression; the second one shows the hardening effect of steel under cyclic loading, which is exponential function. The cyclic hardening coefficient K′ and the cyclic hardening exponent n′ reflect the degree of the hardening of steel under cyclic loading [22].
The cyclic enhancement coefficient K′ under different degrees of corrosion is shown in Figure 6. It can be seen that, with the increase of the rate of weight loss, the steel circulation strengthening coefficient decreases, whereas the variation of cyclic hardening exponent is small and does not have an obvious pattern. According to the definition of the exponential function, the decrease of the cyclic hardening coefficient indicates that, with the increasing degree of weight loss of the steel, the strengthening effect of the steel under cyclic loading is weakened. According to Figure 7, the relationship between the coefficient of cyclic hardening and the weight loss rate of steel plate is as follows:(4)K′=85.926×exp-ρw0.046+531.090.
Relationship between recycling intensified factor and weight loss rate.
The change of pattern of the cyclic hardening index n′ is not obvious. The noncorroded steel plate is 0.211. The steel plate with corrosion rate of 19% is 0.201. When the weight loss rate is within 18%, the circular enhancement index n′ is 0.16 n, in which the value of n can be 0.95–1.05.
4. Constitutive Model of Steel
Figure 2 shows that the hysteretic curve is mainly composed of three parts, the skeleton curve, the unloading curve, and the curve of the loading curve. The first loading curve of the steel extends along the monotonic curve, which is uniaxial tension effect. The unloading of the steel is in accordance with the elastic linear unloading, and the unloading stiffness is the same as that of the initial elastic modulus Es. See Sections 2.2 and 3.1 for the first loading and the changes of the peak value of the steel. The unloading stiffness (initial elastic modulus Es) is shown in Table 2.
4.1. Cyclic Skeleton Criterion
The Ramberg-Osgood model can relatively accurately simulate the shape of the circular skeleton curve; however, in this model, stress is the independent variable and strain is given as a function of stress, and it is difficult to obtain an explicit analytical expression of stress as a function of strain. This fact significantly increases calculation work when a strain value is given (e.g., measured from experiments) to initiate the program. According to [20, 21], and combining with the experimental data of this paper, a more efficient two-stage cyclic skeleton curve is proposed in the following and its performance is evaluated by the comparison with the experimental results in Figure 8. In the proposed curve, the first stage is elastic stage before yield and the second stage is the cyclic hardening stage, as given in (5)σ=σy;ε≤εya+b×ec×ε;ε>εy.
Model for skeleton curve of A11.
In the formula, a, b, and c are used to control the shape of the circular skeleton curve. As the yield effect is not obvious under cyclic loading, the yield strength of the steel is not easy to determine. However, in the past researches on the behavior of the steel plate under monotonic tension, the corresponding strain values of the yield strength of the noncorroded steel plate and the corroded steel plate were all between 0.0013 and 0.0012. This paper took the stress which corresponded to the strain 0.0012 as the yield strength. Take A11 as an example. The comparisons of the tests and models are shown in Figure 8. It can be seen that the proposed form of function simulates the cyclic stress characteristics well under different rates of weight loss of Q235, which can reflect the curve of the circular skeleton.
Cyclic hardening parameters are shown in Table 3. According to the experimental data of corrosion rate and hysteresis curve, and the relationship between the parameters of the cyclic skeleton curve and the weight loss rate is shown in (6)σy=-23.898×exp-ρw-12.612+260.006;a=234.421×exp-ρw4.961+305.723;b=-148.302-169.314×0.651ρw.
Cyclic hardening parameters.
Number
A01
A11
A21
A31
A41
A51
A61
A71
σy
236.451
222.67
216.132
215.046
210.950
208.048
179.942
154.496
a
541.571
381.042
358.927
349.929
348.703
345.251
311.33
306.544
b
-317.442
-171.496
-153.727
-148.409
-149.861
-149.231
-142.583
-159.821
c
-26.699
-56.401
-67.515
-74.386
-56.205
-64.317
-71.826
-44.761
The pattern of the variation of coefficient C is not obvious. The value of coefficient C of the noncorrosive test pieces is −26.7. Under different rates of weight loss, the value of C is −62.2 K, and the K is within the range of 0.70–1.20.
4.2. Reload Criterion
Reloading curve mentioned in the calculation model [21] was applied and modified according to the results of the experiment, as shown in Figure 9.
Reloading curve model.
As shown in Figure 9, the CA section is the unloading curve, and the AB section is the reloading curve for tension. The point A is the starting point of reloading curve (εAi,0), and εAi is defined as the corresponding values of displacement system; the point B is the target point of reloading curve for tension (εBi,σBi), defined as the maximum value for each cycle of the loop, referring to frame curve criterion; the point C is the target point of reloading curve for compression (εBi,σBi), defined as the maximum value for each cycle of the loop, referring to frame curve criterion; Es is unloading slope (see Table 2); Ek is the angle between line a and ε axis. Assume that X point is required; then line b is Ek(ε-εAi), line e is Es(ε-εAi), and line c+d is (Es-Ek)(ε-εAi). Assume that the ratio of line d to line c+d is η; then line d is η(Es-Ek)(ε-εAi).
The formula of stress and strain of X point is shown in (7)σ=Esε-εA+σA-Es-Ekε-εAη,Ek=σB-σAεB-εA,where (8)η=a-bε-εA/εB-εA+b.
According to the results of the experiment, the coefficients a, b have no obvious patterns. The value of a is from 1.05 to 1.4, and the value of b is from 0.05 to 0.6.
4.3. Experiment and Model Comparison
According to the several criteria above and the results of the hysteretic curve test, the criterion of the development of hysteretic curve is shown in Figure 10. The target reaches the yield strength A1 when the curve is loaded on the 0-A1 curve for the first time (see Sections 2.2 and 3.1). According to the criterion of cyclic skeleton curves (see Section 3.1), it reaches A2. After reaching the unloading point of the steel, target will be unloaded from the elastic line to 01, assuming the unloading stiffness and the initial elastic modulus Es are the same (see Table 2). Then the target reaches the peak value B2 in the first lap of the compression cycle (see Section 3.2). After being unloaded from the straight line, the target will reach -01. The unloading stiffness and the initial elastic modulus Es are the same (see Table 2). After the target being pulled up to the maximum of A2, the second lap will start, and it will repeat the steps above (i.e., unloaded to 01, pressed to B2, and unloaded to -01). Finally the target will move along the reloading curve and pass through A2 to the unloading point A3 (see Sections 3.1 and 3.2); after elastic unloading, the target will reach the starting point of reloading curve for compression. Then the whole process above will be repeated.
The specific description.
By comparing the curve of cyclic constitutive model and hysteretic curve of steel plate under cyclic loading, A11, for example (Figure 11), had shown that the model matched the results of the experiment well. The steel loop constitutive model presented in this paper not only could predict the hysteretic behavior of the steel under different rates of weight loss but also could predict the residual strength of the specimen.
Comparison chart between model of the hysteresis criterion and experimental data of A11.
5. Conclusion
Both noncorroded steel and corroded steel had shown behaviors like cyclic hardening, Bauschinger effect, and peak point effect under cyclic loading. The hysteresis loop was relatively plump. With the increase of the degree of weight loss, the mechanical properties of the steel were decreased, and the deterioration of ductility was aggravated.
The Ramberg-Osgood model matched the skeleton curve of corroded steel under cyclic loading well. With the increase of the degree of weight loss, the coefficient of steel circular enhancement showed a decreasing trend, whereas the variation of cyclic hardening exponent was small and did not have an obvious pattern.
In this paper, the hysteretic constitutional relationship of corroded steel was established based on the first time loading criterion, unloading criterion, cycle skeleton criterion, and reloading curve criterion. The model matched the results of the experiment well. Therefore, the steel loop constitutive model presented in this paper not only could predict the hysteretic behavior of steel under different rates of weight loss but also could predict the residual strength of the specimen.
Competing Interests
The authors declare that they have no competing interests.
Acknowledgments
The financial support is provided by the National Natural Science Foundation of China under Grant no. 51378417. This paper is also supported by Program for Changjiang Scholars and Innovative Research Team in University (IRT13089).
ANSI/AISC1999Chicago, Ill, USAAmerican Institute of Steel ConstructionKřivýaV.KonečnýbP.Real material properties of weathering steels used in bridge structures201357624633KravanjaS.ŽulaT.Cost optimization of industrial steel building structures201041344245010.1016/j.advengsoft.2009.03.005ZBL1183.901472-s2.0-72549089747YatesT. J. S.CooteA. T.ButlinR. N.The effect of acid deposition on buildings and building materials198821202610.1016/0950-0618(88)90005-02-s2.0-45549118262LaqueF. L.1975New York, NY, USAJohn Wiley & SonsWhiteE. E.WatkinsK. O.Design for prevention of corrosion in buildings and structures1994914161MelchersR. E.Effect of small compositional changes on marine immersion corrosion of low alloy steels20044671669169110.1016/j.corsci.2003.10.004MelchersR. E.Modeling of marine immersion corrosion for mild and low-alloy steels—part 1: phenomenological model200359431933410.5006/1.32775642-s2.0-0037648521MelchersR. E.WellsT.Models for the anaerobic phases of marine immersion corrosion20064871791181110.1016/j.corsci.2005.05.0392-s2.0-33744790872WangZ.-X.YuG.-A.HanW.Investigation on interrelation of indoor accelerated corrosion and atmospheric exposure corrosion of steels20041627073LiH. B.JiangZ. H.FengH.ZhangS. C.LiL.HanP. D.MisraR. D. K.LiJ. Z.Microstructure, mechanical and corrosion properties of friction stir welded high nitrogen nickel-free austenitic stainless steel20158429129910.1016/j.matdes.2015.06.1032-s2.0-84941276996SultanaS.WangY.SobeyA. J.WhartonJ. A.ShenoiR. A.Influence of corrosion on the ultimate compressive strength of steel plates and stiffened panels2015969510410.1016/j.tws.2015.08.0062-s2.0-84940383774WangY.WhartonJ. A.ShenoiR. A.Ultimate strength analysis of aged steel-plated structures exposed to marine corrosion damage: a review201486426010.1016/j.corsci.2014.04.0432-s2.0-84903490051KhedmatiM. R.RoshanaliM. M.NouriZ. H. M. E.Strength of steel plates with both-sides randomly distributed with corrosion wastage under uniaxial compression201149232534210.1016/j.tws.2010.10.0022-s2.0-78650706111XuS.-H.WangY.-D.Estimating the effects of corrosion pits on the fatigue life of steel plate based on the 3D profile201572274110.1016/j.ijfatigue.2014.11.0032-s2.0-84912138947BruneauM.ZahraiS. M.Effect of severe corrosion on cyclic ductility of steel1997123111478148610.1061/(ASCE)0733-9445(1997)123:11(1478)2-s2.0-0031279767DhakalR. P.MaekawaK.Path-dependent cyclic stress–strain relationship of reinforcing bar including buckling200224111383139610.1016/s0141-0296(02)00080-92-s2.0-0036837707CEB1996London, UKThomas Telford ServicesMendesL. A. M.CastroL. M. S. S.A simplified reinforcing steel model suitable for cyclic loading including ultra-low-cycle fatigue effects20146815516410.1016/j.engstruct.2014.02.0312-s2.0-84897472105ShiG.WangM.BaiY.WangF.ShiY.WangY.Experimental and modeling study of high-strength structural steel under cyclic loading20123711310.1016/j.engstruct.2011.12.0182-s2.0-84862816109ShiY. J.WangM.WangY. Q.Experimental and constitutive model study of structural steel under cyclic loading20116781185119710.1016/j.jcsr.2011.02.0112-s2.0-79955483678RambergW.OsgoodW. R.Description of stress-strain curves by three parameters1943902National Advisory Committee for AeronauticsGB/T 10125Beijing, ChinaStandards Press of China2012 (Chinese)GB/TCorrosion of metals and alloys. Outdoors exposure test methods for periodic water spray24517Beijing, ChinaStandards Press of China2009 (Chinese)GB/T2975Beijing, ChinaStandards Press of China1998 (Chinese)