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 parameter
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.
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.
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).
(1) |
(2) |
The corresponding data pairs are σ i and ϵ i - σ i/ E where i = 1, 2, 3, …….. k. Given that
(3) |
(4) |
(5) |
(6a) |
It can be deduced
(6b) |
As K and n are both nonzeros, Eq. (6)b can be further simplified to
(7) |
For a given set of σ i, ϵ i ( i = 1,2,…., k), the corresponding K and n parameters can then be determined.
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.
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:
The predicted K values show good correlation with the experimental results as summarized in Table 10.
Based on the discussion from Section 2, the proposed constitutive model of cyclic stress-strain behavior is summarized as follows:
(9) |
(10) |
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.
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.
This section presents a numerical study based upon the proposed constitutive model in Section 3.
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].
![]() | Fig.10 Finite element model |
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.
![]() | Fig.11 Hysteresis loops. |
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%.
1 |
|
2 |
|
3 |
|
4 |
|
5 |
|
6 |
|
7 |
|
8 |
|
9 |
|
10 |
|
11 |
|
12 |
|
13 |
|
14 |
|
15 |
|
16 |
|