1. Introduction
After the 1994 Northridge earthquake in California and the 1995 Kobe earthquake in Japan, it was observed that beam-to-column connections and baseplate connections in steel bridge piers undergo ultra-low cycle fatigue (ULCF) in such events [
1,
2,
3]. This ULCF damage has been shown to cause the progressive collapse of entire structures [
4]. The ULCF is characterized by a ductile crack that initiates at the strain concentration position, then this crack expands stably under cyclic loading, and finally, the catastrophic failure occurs in the brittle mode. ULCF has been shown to occur in the areas of strain concentration in steel bridge piers and beam-column connections under cyclic loading [
5,
6,
7]. Unlike traditional high cycle and low cycle fatigue, ULCF experiences large plastic strain amplitude and is usually characterized by few reverse loading cycles (in general less than 100). Therefore, ULCF is of great significance in the seismic design of steel structures.
The Coffin-Manson formula [
8,
9] has been widely used to predict the low cycle fatigue life of steel structures. In order to predict ULCF life, Ge et al. [
10,
11,
12] introduced a damage index to evaluate the ULCF life in steel bridge piers based on the Coffin-Manson formula and Miner’s rule [
13]. Tateishi et al. [
14] developed a new fatigue prediction model that can accurately predict the fatigue life of plain material in an extremely large strain range. Xue [
15] proposed a uniform expression to predict low cycle fatigue and ULCF by introducing an exponential function and additional material parameters. However, the above empirical models are derived under uniaxial strain conditions and cannot be applied to a multiaxial stress condition. Micro-mechanism-based models have been proposed to solve these problems in recent years and will be discussed below.
The micro-mechanical models in the literature can be classified into coupled and uncoupled models [
16]. The coupled models consider the intercoupling between material constitutive properties and damage. Mear et al. [
17] and Leblond et al. [
18] modified the Gurson–Tvergaard–Needleman (GTN) model for cyclic loading. Tong et al. [
19] proposed a model based on continuous damage mechanics (CDM) to investigate the ULCF behaviour of beam-column connections. However, in coupled models, the model parameter calibration is a complex task due to the interdependency among the parameters as well as the high calculation cost. These shortcomings impede the application of the coupled models. The uncoupled models can be efficient for crack initiation modelling. Since the uncoupled models assume independence between the material constitutive properties and damage, the parameter calibration is simpler compared to that of the coupled models, and the most accurate state-of-the-art constitutive models can be used in the uncoupled models.
Several uncoupled models have been presented in the literature. Kanvinde and Deierlein [
20] proposed the cyclic void growth model (CVGM) to predict ULCF life of structural steel based on the Rice-Tracey void growth theory [
21]. Owing to some advantages of the CVGM model, such as predicting fracture initiation at the continuous level and being suitable for multiaxial stress condition compared to empirical models, it has received extensive attention in attempting to predict the ULCF damage of steel structural members. Myers et al. [
22] and Fell et al. [
23] investigated the ULCF fracture initiation of column baseplate connections and steel frame brace components, respectively. Zhou et al. [
24] investigated the ULCF behaviour of beam-column connections, and Liao [
25] conducted ULCF fracture predictions for welded connection combined square steel pipe column and H-shaped steel beam. In general, the above research results are encouraging, and it is promising to predict the ULCF fracture initiation of steel structural members using the CVGM model. However, the results predicted by the CVGM model largely depend on the calibration values of the model parameters, especially the calibration value of the damage degradation parameter. In some literature [
25,
26], the calibration results of damage degradation parameters of the CVGM model is rather discrete. Therefore, it is significant to investigate the effect of model parameter calibration methods on the prediction accuracy of CVGM model. Additionally, it has been suggested that the Lode angle parameter should also be accounted for ULCF modelling except for stress triaxiality (
T = σ
h/σ
e), defined as the ratio of the hydrostatic stress (σ
h) and the Mises stress (σ
e), especially in the case of low stress triaxiality (T < 0.33) [
27,
28,
29]. However, high stress triaxiality (T > 0.70) and moderate stress triaxiality (0.33 < T < 0.70) are employed in the present study, and the effect of Load angle parameter can be neglected.
In order to improve the prediction accuracy of the CVGM model, tests were conducted on circular notched specimens made of Q345qC steel commonly used in the construction of steel bridges in China. A model parameter calibration method was proposed, and model parameters of the CVGM for Q345qC steel were calibrated at both high and moderate stress triaxiality based on the experimental results and finite element analysis (FEA). Comparisons were made between the experimental results and predicted results to verified the effectiveness of proposed parameter calibration method. Finally, the effect of damage degradation parameter on ULCF life prediction in steel bridge piers was discussed.
3. ULCF Fracture Initiation Prediction
In order to examine the accuracy of the CVGM model at high and moderate stress triaxiality, cyclic tests were conducted on circular notched specimens providing both high and moderate stress triaxiality. The contour and dimensions of specimens are presented in
Figure 7 and
Table 5, respectively. Loading was controlled by strain with CTF loading, the gauge length of the extensometer was 50 mm, and the loading strain is provided in
Table 5.
The two-dimensional axisymmetric finite element model was established using ABAQUS 6.14 with the element division similar to that of the aforementioned tensile loading. The Lemaitre-Chaboche hybrid hardening model [
36] was used to simulate the cyclic plastic flow of material.
Figure 11 presents an example of the force-displacement plot for a specimen, and it can be seen that the test curve is in good agreement with the FEA curve. The development of the damage index,
D, of partial specimens at the fracture initiation location is presented in
Figure 12. ULCF fracture is predicted to occur in the notched specimens when
D reaches one.
Figure 13 presents the comparison of the experimental results and predicted results, including the number of cycles,
Nf, and equivalent plastic strain,
, to fracture initiation. To evaluate the prediction accuracy of the CVGM model, relative error,
, is calculated according to Equation (13), and the results are presented in
Table 6. The results indicate that the predicted life and equivalent plastic strain by a segmentally calibrated CVGM model is nearer to the experimental results compared to the originally calibrated CVGM model. For a small part of specimens, the errors are larger than or close to 25% between experimental results and predicted results by segmentally calibrated CVGM model, which may be related to some assumptions introduced in the model, but the error is acceptable. In general, the mean value of
calculated by the segmentally calibrated CVGM model is approximately 50% smaller than that of the initially calibrated CVGM model. Therefore, separately calibrating the CVGM model at high and moderate stress triaxiality greatly improved the prediction accuracy of ULCF.
4. Effect of Damage Degradation Parameter
The different values of the damage degradation parameter for Q345qC steel were obtained at different stress triaxiality ranges. To investigate the effect of the damage degradation parameter on predicting ULCF fracture initiation of steel structural members, the single-column steel bridge pier was applied as the research object. In this study, two fillet welds were used with groove angles of 45° and 56°, as presented in
Figure 14. In
Figure 14,
h represents the height of the pier,
B represents the width of the flange,
W represents the width of the web,
t represents the thickness of the flange and the web,
bs represents the width of the stiffener,
ts represents the thickness of the stiffener,
td represents the thickness of the diaphragm and baseplate, and
a represents the spacing of the diaphragms. The geometric dimensions of the steel bridge pier are listed in
Table 7. The axial pressure,
P, and horizontal forced displacement,
δ, were applied at the top of the pier, and the axial compression ratio is 0.15. The loading pattern is shown in
Figure 15, where
δy is the horizontal yield displacement of the pier.
The numerical analyses were conducted using ABAQUS 6.14. Three-dimensional finite element models of the steel bridge pier were established, as shown in
Figure 16. The beam element (B31) was employed to simulate the upper part of the steel bridge pier. The lower part of the pier,
Ld, was simulated by the shell element (S4R) where the length of
Ld was determined according to the empirical formula for calculating the damaged domain [
40]. To obtain the stress–strain history at the bottom of the boundary between the flange and web, the solid element (C3D8R) was used. The MPC-beam connection was adopted between the beam and shell element. The shell and solid element were coupled by Shell-Solid. The most dangerous element in the point of connection between the flange and web was considered to be the calculating element, and the size of the calculating element was 0.20 mm, which is consistent with the characteristic length of Q345qC [
34]. The Lemaitre–Chaboche hybrid hardening model was used to simulate the cyclic plastic flow of material, and the model parameters used are provided in
Table 8 [
34].
The ULCF life of steel bridge piers with different groove angles was calculated using the above calibrated CVGM model, and the calculation results are shown in
Table 9. The predicted ULCF life by segmentally calibrated CVGM model varied by a maximum of 37.5% as compared to the predicted life by initially calibrated CVGM model. Therefore, based on the calculation results, the calibrated damage degradation parameter has a significant influence on the ULCF fracture prediction, and it is important to accurately calibrate the damage degradation parameters under different stress triaxiality ranges. Additionally, due to greater strain concentration, the steel bridge pier with a groove angle of 56° is prone to ULCF damage compared with the steel pier with a groove angle of 45°.
5. Conclusions
In the research presented in this paper, Q345qC steel was selected for study as it is commonly used in the construction of steel bridges in China. A model parameter calibration method that separately calibrates the damage degradation parameter at high and moderate stress triaxiality was proposed. The validity of the CVGM model calibrated by the proposed method was verified based on the tests and FEA. The effect of the damage degradation parameter on predicting the ULCF fracture initiation of steel bridge piers was investigated. Based on the conducted studies, the following conclusions can be drawn:
- (1)
The goodness of fit is improved when the damage degradation parameter is calibrated separately at high and moderate stress triaxiality. It is shown that the predicted number of cycles and equivalent plastic strain to fracture by the segmentally calibrated CVGM model agree well with the experimental results. The proposed parameter calibration method greatly improved the predictive accuracy of the CVGM model compared to the previous method.
- (2)
From the numerical simulation of steel bridge piers, the value of the damage degradation parameter has a relatively significant influence on the predicted cycles to ULCF fracture initiation. Therefore, it is significant to separately calibrate the damage degradation parameters under different stress triaxiality ranges.
- (3)
The different values of the damage degradation parameter were obtained at high and moderate stress triaxiality, and it was determined that the value of the damage degradation parameter depends on stress triaxiality. Therefore, it is clear that the relationship between damage degradation parameter and stress triaxiality deserves further study in the future.
In future work, additional experimental data are needed to determine the relationship between damage degradation parameter and stress triaxiality. Additionally, the different values of damage degradation parameter at different levels of stress triaxiality may be related to the fracture mechanism, and as such it requires further research. In general, the proposed parameter calibration method that separately calibrates parameters under different stress triaxiality ranges is promising.