Frontiers of Structural and Civil Engineering

ISSN 2095-2430

ISSN 2095-2449(Online)

CN 10-1023/X

邮发代号 80-968

2019 Impact Factor: 1.68

Cyclic stress-strain behavior of structural steel with yield-strength up to 460 N/mm<sup>2</sup>
Cyclic stress-strain behavior of structural steel with yield-strength up to 460 N/mm2
Yiyi CHEN1, Wei SUN2, Tak-Ming CHAN3
1. State Key Laboratory of Disaster Reduction in Civil Engineering, Tongji University, Shanghai 200092, China
2. Architecture Design & Research Institute of Tongji University (Group) Co. Ltd., Shanghai 200092, China
3. School of Engineering, University of Warwick, Coventry CV4 7AL, UK
t.m.chan@warwick.ac.uk
Abstract

This paper presents a constitutive model based on Ramberg-Osgood equation to describe the hysteresis material behavior of structural carbon steel with nominal yield strength between 235 to 420 N/mm2. The proposed model was calibrated against a series of cyclic material tests with strain amplitude varying from 0.5% to 2.0%. A simple relationship between the modular parameterK and the yield strengthfy was proposed. The calibrated Ramberg-Osgood model revealed excellent agreement with the experimental results and captured further the experimental behavior of test specimens with nominal yield strength of 460 N/mm2. The proposed constitutive model was also adopted in conjunction with the combined kinematic/isotropic materials description in ABAQUS to mimic a full scale experimental test under cyclic loading. The numerical results revealed close agreement with the experimental observations.

Keyword: constitutive model; finite element analysis; hysteresis curve; Ramberg-Osgood equation
Introduction

Ramberg and Osgood [ 1] first proposed an analytical formulation to describe the nonlinear behavior of material under compression and tension. Based on Ramberg and Osgood relations, numerous studies have been conducted to modify and improve the relationship to cover a wider range of strain for nonlinear materials under monotonic loads [ 2 6] and to describe the stress-strain relationship at elevated temperature [ 7, 8]. Concerning hysteresis behavior, Skelton et al. [ 9] described the relationships between the various Masing-type models and the Ramberg-Osgood stress-strain response. A new link was established between the Ramberg-Osgood strength and cyclic hardening parameters. Broggiato et al. [ 10] also described the use of the Chaboche nonlinear kinematic hardening model for cyclic plasticity assessment. The tuned Chaboche model for elastoplastic material description was employed in the finite element simulation which was capable of following the cyclic elastoplastic response very closely up to large accumulated plastic strain levels. However, limited research [ 11] is conducted in the development of hysteresis stress-strain models based on Ramberg-Osgood relationships for structural carbon steel. This paper presents the calibration of the Ramberg-Osgood equations against a series of cyclic materials tests [ 12] and proposes a constitutive model to describe the material hysteresis performance.

Calibration

A series of experimental tests on steel coupons under monotonic and cyclic loadings were conducted by Chen et al. [12]. Mean measured dimensions and key results from direct tensile tests are summarized in Table 1. The specimens are categorized according to the material grade as shown in columns 1 and 2, measured yield strength, measured ultimate strength, ratio of yield to ultimate strengths and measured elongation are shown in columns 3 to 6 respectively. Further experimental details can be found in Chen et al. [ 12].

Tab.1 Mean measured dimensions and key results from direct tensile tests

For the cyclic material tests, three different loading protocols 1) cyclic ascend 2) cyclic alternate and 3) cyclic tensile were employed as shown in Fig. 1. As shown in Fig. 1(a), cyclic ascend is a loading protocol with an equal strain increment of+/‒0.005 for each five-cycle loading phase up to+/- 0.02 strain amplitude. cyclic alternate refers to the loading scheme as shown in Fig. 1(b) which is similar to cyclic ascend, except the second and third loading phases were interchanged. The loading scheme of cyclic tensile is indicated in Fig. 1(c). All loadings are in tension and it oscillates incrementally up to 0.08 strain amplitude. Previous research [ 12] revealed that the loading protocols in particular between cyclic ascend and cyclic alternate affected the hysteresis behavior but the effect is minimal. This paper primarily focuses on the hysteresis response from cyclic ascend and cyclic alternate loading protocols.

Fig.1 Cyclic loading protocols. (a) Cyclic ascend; (b) cyclic alternate; (c) cyclic tensile

Ramberg-Osgood model

The most commonly used expression to describe nonlinear material stress-strain behavior is that proposed by Ramberg and Osgood [ 1] as shown in Eq. (1).

ϵ=σE+(σK)n(1)
where, σ and ϵ are engineering stress and strain, respectively, E is the material Young’s modulus which was taken as 2 × 105 N/mm2 in this study, n is a strain hardening exponent and K is a modular parameter. Grouping terms, Eq. (1) can be simplified to
y^=(σK)n(2)
where y^=ϵ-σE.

The corresponding data pairs are σ i and ϵ i - σ i/ E where i = 1, 2, 3, …….. k. Given that

yi=ϵi-σiE(3)
y^i=(σiK)n(4)
and adopting the least-square method, the target is to minimize the following function to obtain an optimum pair of modular parameter K and strain-hardening exponent n.
f(K,n)=i=1k(yi-y^i)2(5)
From
{f(K,n)K=0f(K,n)n=0(6a)

It can be deduced

{nKn+1i=1kσin(yi-y^i)=01Kni=1kσin(lnσi-lnK)(yi-y^i)=0(6b)

As K and n are both nonzeros, Eq. (6)b can be further simplified to

{i=1kσin(yi-y^i)=0i=1kσinlnσi(yi-y^i)=0(7)

For a given set of σ i, ϵ i ( i = 1,2,…., k), the corresponding K and n parameters can then be determined.

Cyclic stress-strain curves

From Chen et al. [ 12], typical hysteresis loops are shown in Fig. 2. The hysteresis loops for Group 1 specimens under cyclic ascend (1-1) and cyclic alternate (1-2) loading protocol are shown in Figs. 2(a) and 2(b) respectively.

In general, the hysteresis loops are bi-symmetric, therefore, to simplify the calibration process against Ramberg-Osgood equations, the hysteresis loops in the first quadrant (I) are calibrated and the calibration results are shown in Tables 2 to 8.

Tab.2 Group 1 specimens Ramberg-Osgood parameters
Tab.3 Group 2 specimens Ramberg-Osgood parameters
Tab.4 Group 3 specimens Ramberg-Osgood parameters
Tab.5 Group 4 specimens Ramberg-Osgood parameters
Tab.6 Group 5 specimens Ramberg-Osgood parameters
Tab.7 Group 6 specimens Ramberg-Osgood parameters
Tab.8 Group 7 specimens Ramberg-Osgood parameters

Tables 2 to 8 show the calibrated parameters for Group 1 to 7 specimens under cyclic ascend and cyclic alternate loading protocols at different amplitudes. The coefficient of determination r2 between the experimental data and proposed hysteresis models are also tabulated.

As shown in the above tables, all coefficient of determination r2 are close to 1 showing close relation between experimental data and Ramberg-Osgood equations. Under different strain amplitudes, for example in Group 7 specimens (Table 8), the modular parameter K remains fairly constant while the strain-hardening exponent, n increases with the strain amplitude.

Fig.2 Hysteresis loops for Group 1 specimens. (a) Cyclic Ascend (1-1); (b) cyclic alternate (1-2)

Proposed K fy relationship

A summary of modular parameters for hysteresis loops under cyclic ascend loading protocol is shown in Table 9.

Tab.9 Modular parameters for hysteresis loops
Tab.10 Modular parameter and yield strength relation

As shown in the table, under different imposed strain amplitude, the variation of modular parameter K is minimal, for example for specimen 2-1, the maximum modular parameter is 1071 MPa while the minimum modular parameter is 1003 MPa, so an average value was used. It is also observed that as the yield stress increases, the modular parameter K displays an increasing trend. It is therefore proposed a linear relationship between modular parameter K and yield strength as follows: K=1.79fy+532.

The predicted K values show good correlation with the experimental results as summarized in Table 10.

Constitutive model
Primary skeleton stress-strain curve

Based on the discussion from Section 2, the proposed constitutive model of cyclic stress-strain behavior is summarized as follows:

ϵ={σE+(σK)n,σ0σE-(-σK)n,σ<0(9)
where K and n can be evaluated based on the methodology as shown in Section 2.1 for cyclic ascend loading protocol.

Cyclic stress-strain curve

ϵ={σ-σuE+(σ-σuK)n+ϵu,loadingieσ-σu>0σ-σuE-(σu-σK)n+ϵu,unloadingieσ-σu<0(10)
where σu, ϵu are the unloading reference point at each cycle, E is taken as 2 × 105 MPa, K can be obtained from the relationship proposed in Eq. (8) and n can be defined by the approach in Halford and Morrow [ 13]. Hence, to define full hysteresis stress-strain behavior of structural carbon steel, two key parameters which are K and n are needed. The parameters K and n also depend on the strength of steel, and the latter changes significant with the increase of maximum strain.

Based on Eqs. (9) and (10), the derived hysteresis stress-strain curves were compared with the experimental results. For Group 1 specimens, as shown before, the corresponding parameters are K = 1067 MPa (primary skeleton), n = 4.34 (primary skeleton), E = 2 × 105 MPa and K = 1033 MPa (hysteresis behavior as shown in Table 10). For specimens under cyclic ascend loading protocol, the straining history is 0→0.005→-0.005→0.01→-0.01→0.015→-0.015→0.02→-0.02→0.02 and for specimens under cyclic alternate loading protocol, the straining history is 0→0.005→-0.005→0.015→-0.015→0.015→-0.01→0.01→-0.01→0.02→-0.02→0.02. The comparison between the proposed constitutive model and the experimental results are shown in Figs. 3 and 4.

The above figures revealed the proposed hysteresis models correlate closely with the experimental data. The proposed material stress-strain model can also capture the effect of different loading protocols. Figure 5 shows the hysteresis loops for Group 1 specimens under cyclic ascend (1-1) and cyclic alternate (1-2). The difference is indicated at strain amplitude of 0.01, the cyclic peak values for specimen 1-2 is greater than 1-1, which illustrates the effect of loading protocol but as discussed in Chen et al. [ 12], the effect is insignificant toward seismic applications.

In general, the calibrated Ramberg-Osgood equation shows excellent agreement with the experimental results on the primary skeleton and the cyclic stress-strain curves.

Fig.3 Hysteresis loops for Group 1 specimen under cyclic ascend loading protocol

Fig.4 Hysteresis loops for Group 1 specimen under cyclic alternate loading protocol

Fig.5 The effect of loading protocols on hysteresis loops

Application to Q460 specimens

The proposed constitutive model in Section 3 was calibrated by experimental coupons with nominal yield strengths between 235 to 420 MPa. To investigate their applicability to a higher grade of steel, this section presents the comparison between the proposed constitutive model and the experimental result on coupons with nominal yield strength of 460 MPa (Q460 specimens). The direct tensile tests from the additional Group 8 specimens were conducted by the Authors on Grade Q460E materials with the measured yield strength of 491 MPa, ultimate strength of 611 MPa, yield-ultimate strength ratio of 0.8 and the elongation of 22.7%. Other Group 9 specimens on Grade Q460C materials were tested by Shi et al. [ 14] with the reported yield strength of 466 MPa.

Based on the proposed constitutive models (Eqs. (8)‒(10)), Figs. 6 and 7 illustrate the comparison between the proposed models and the experimental observations which reveal close agreement between the two.

Fig.6 Hysteresis loops for Group 8 specimen under cyclic ascend loading protocol

Fig.7 Hysteresis loops for Group 9 specimen under cyclic ascend loading protocol

Numerical validation

This section presents a numerical study based upon the proposed constitutive model in Section 3.

Validation experiment

The validation experiment refers to a cyclic test of tubular column of steel grade Q235B with nominal yield strength of 235 N/mm2. The geometry of the test arrangement is shown in Fig. 8. The bottom part of the specimen was anchored to the strong floor facilities at the Structures laboratory of Tongji University. The loading was applied through the top part of the specimen under displacement-controlled loading scheme. The loading scheme (Displacement ratio Δ/Δ1 vs. number of cycle n) is shown in Fig. 9 which is similar to the cyclic ascend loading protocol as shown in Fig. 1(a). As shown in Fig. 9, Δ1 (0.7 mm) is the yield displacement. Further experimental details are presented in Zhao et al. [ 15].

Numerical simulation

A numerical modeling investigation, employing the finite element (FE) package ABAQUS [ 16], was carried out to replicate the experimental results based on the proposed constitutive model in Section 3. The element chosen for the FE model was a general purpose linear brick element, with reduced integration, designated as C3D8R in the ABAQUS element library. The mesh density was carefully chosen by carrying out a mesh convergence study with the aim of achieving accurate results while minimizing computational effort. The meshed finite element model is shown in Fig. 10.

The true material stress-strain relationships were generated from the engineering stress-strain curves obtained from the tensile coupon tests and cyclic coupons [ 12]. A kinematic/isotropic combined model [ 16] was employed in which the kinematic hardening parameters employed were determined through the full stress-strain response of specimen 1-1 from Eqs. (8) and (9) under four different strain amplitudes. Boundary conditions were applied to the model fixed end support conditions at both ends but allowing free vertical displacement at the top. The vertical displacement at top end was monitored throughout the analysis. Geometric and material nonlinearity was incorporated in the models which also enabled the cyclic behavior to be traced. The numerical derived hysteresis curves are shown in Fig. 11 and are compared with the corresponding experimental hysteresis results. Replication of test results was found to be satisfactory with the numerical models able to successfully capture the initial primary stress-strain behavior and general hysteresis response.

Conclusions

This paper proposes a constitutive model based on Ramberg-Osgood equations to describe the primary skeleton stress-strain behavior and the cyclic stress-strain response for structural carbon steel with yield strength up to 420 N/mm2. The proposed model was calibrated against a series of cyclic material tests with nominal yield strength between 235 to 420 N/mm2 which revealed close agreement between the proposed model and the experimental results. A linear relation between modular parameter K and yield strength fy was also proposed. Two key parameters K and n were defined which were used to describe the full stress-strain response under cyclic load. The proposed model was also able to capable the material hysteresis performance of coupons with nominal yield strength of 460 N/mm2. The proposed constitutive model was then employed in a numerical validation which has showed excellent agreement between the numerical and the experimental results. This research forms part of an on-going research program at Tongji University to investigate the material hysteresis response of structural carbon steel with strain amplitude up to 10%.

References
1 Ramberg W, Osgood W R. Description of stress-strain curves by three parameters. National Advisory Committee for Aeronautics, Technical Note No. 902, 1943 [Cited within: 2]
2 Mirambell E, Real E. On the calculation of deflections in structural stainless steel beams: An experimental and numerical investigation. Journal of Constructional Steel Research, 2000, 54(1): 109133 [Cited within: 1] [JCR: 1.327]
3 Rasmussen K J R. Full-range stress-strain curves for stainless steel alloys. Journal of Constructional Steel Research, 2003, 59(1): 4761 [Cited within: 1] [JCR: 1.327]
4 Gardner L, Nethercot D A. Experiments on stainless steel hollow sections. Part 1: Material and cross-sectional behaviour. Journal of Constructional Steel Research, 2004, 60(9): 12911318 [Cited within: 1] [JCR: 1.327]
5 Gardner L, Ashraf M. Structural design for nonlinear metallic materials. Engineering Structures, 2006, 28(6): 926934 [Cited within: 1] [JCR: 1.713]
6 Quach W M, Teng J G, Chung K F. Three-stage full-range stress-strain model for stainless steels. Journal of Structural Engineering, 2008, 134(9): 15181527 [Cited within: 1]
7 Chen J, Young B. Effects of elevated temperatures on the mechanical properties of cold-formed steel. Proceedings in the International symposium on cold formed steel metal structures, 2004, 139158 [Cited within: 1]
8 Chen J, Young B. Stress-strain curves for stainless steel at elevated temperatures. Engineering Structures, 2006, 28(2): 229239. DOI:10.1016/j.engstruct.2005.07.005 [Cited within: 1] [JCR: 1.713]
9 Skelton R P, Maier H J, Christ H J. The Bauschinger effect, Masing model and the Ramberg-Osgood relation for cyclic deformation in metals. Materials Science and Engineering, 1997, 238(2): 377390 [Cited within: 1]
10 Broggiato G B, Campana F, Cortese L. The Chaboche nonlinear kinematic hardening model: Calibration methodology and validation. Meccanica, 2008, 43(2): 115124 [Cited within: 1] [JCR: 1.747]
11 Shi G, Wang M, Bai Y, Wang F, Shi Y J, Wang Y Q. Experimental and modeling study of high-strength structural steel under cyclic loading. Engineering Structures, 2012, 37: 113 [Cited within: 1] [JCR: 1.713]
12 Chen Y Y, Sun W, Chan T M. Effect of loading protocols on the hysteresis behaviour of structural carbon steel with yield strength up to 460 N/mm2. Advances in Structural Engineering, 2013, 16(4): 707719 [Cited within: 6] [JCR: 0.489]
13 Halford G R, Morrow J. Low-cycle fatigue in torsion. Proceedings of the American Society for Testing and Materials, 1962, 62: 695709 [Cited within: 1]
14 Shi G, Ban H Y, Shi Y J, Wang Y Q. Engineering application and recent research progress on high strength steel structures. Industrial Construction, 2012, 42(1): 17 (in Chinese) [Cited within: 1] [CJCR: 0.566]
15 Zhao X Z, Ge L J, Chen Y Y, He M X. Research on the effect of assembling weld form on mechanical behavior of box brace. Journal of Building Structures, 2008, 29(3): 8895 (in Chinese) [Cited within: 1] [CJCR: 1.074]
16 ABAQUS. Analysis User’s Manual I–V. Version 6. 7. USA: ABAQUS, Inc, 2007 [Cited within: 2]