Next Article in Journal
A Hybrid CNN-Based Review Helpfulness Filtering Model for Improving E-Commerce Recommendation Service
Previous Article in Journal
Reliability-Based Design of Driven Piles Considering Setup Effects
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Semi-Automated Procedure to Estimate Nonlinear Kinematic Hardening Model to Simulate the Nonlinear Dynamic Properties of Soil and Rock

1
Faculty of Civil Engineering, Thuyloi University, 175 Tay Son Street, Dong Da District, Hanoi 100000, Vietnam
2
Department of Civil Engineering, Chonnam National University, Yongbong-ro 77, Buk-gu, Gwangju 61186, Korea
*
Authors to whom correspondence should be addressed.
Appl. Sci. 2021, 11(18), 8611; https://doi.org/10.3390/app11188611
Submission received: 13 August 2021 / Revised: 10 September 2021 / Accepted: 13 September 2021 / Published: 16 September 2021
(This article belongs to the Section Civil Engineering)

Abstract

:
The strain-dependent nonlinear properties of ground materials, such as shear modulus degradation (G/Gmax) and damping, are of significant importance in seismic-related analyses. However, the ABAQUS program lacks a comprehensive procedure to estimate parameters for a built-in model. In this study, a nonlinear kinematic hardening (NKH) model with three back-stress values was used, which allows better fitting to the backbone curves compared to the simplified nonlinear kinematic hardening (SNKH) model previously proposed. Instead of modeling in ABAQUS, a semi-automated procedure was implemented in MATLAB, which can predict shear stress–shear strain hysteretic loops, to find the fitting parameters to the target G/Gmax and/or damping curves. The procedure was applied for three soil and two rock samples, and the results indicate a good match between model and target backbone curves, which proves the application of the procedure and the NKH model in simulating the nonlinear properties of ground materials.

1. Introduction

Shear modulus degradation (G/Gmax) and damping (D) curves are the most important properties that need to be considered when simulating dynamic problems, such as seismic wave propagation in the ground, or dynamic soil–structure interactions (SSIs) [1,2,3,4]. Previous studies have indicated that good modeling of the nonlinear hysteretic properties of soil, including shear modulus degradation and damping, could enable wave propagation and soil–structure interaction to be accurately simulated [2,5].
Among famous commercial programs, such as ABAQUS, ANSYS, LS-DYNA, FLAC, and PLAXIS, ABAQUS has commonly been used for dynamic problems and SSI investigation, because of its capacity and convenience with numerous materials and elements for soil, structure, composite, and porous media. However, constitutive models with the implementation of G/Gmax and damping curves for soils and rocks in ABAQUS have still been quite limited. The simplified nonlinear kinematic hardening (SNKH) model was implemented in ABAQUS to simulate the strain-dependent properties of soil [5,6,7]. Although the procedure for estimating SNKH model parameters was proposed by Anastasopoulos and Gelagoti [5], with a single calibrated parameter based on the ultimate shear strength of material, the kinematic hardening parameters were varied to fit with G/Gmax obtained from experiments, without changing the shear strength at extremely large shear strains. Furthermore, the SNKH model did not accurately simulate the nonlinear dynamic properties of soil, due to the use of a simple nonlinear kinematic hardening rule. Several user-defined materials (e.g., UMAT and VUMAT) imply that more complicated hardening rules had been developed [4,8,9,10]. However, the implementation of these UMAT materials by engineers has been limited due to their complexity and difficulties, as well as the application of the nonlinear kinematic hardening model (NKH) in geotechnical engineering being limited due to the difficulty in estimating its parameters.
In this study, a semi-automated procedure implemented in MATLAB was proposed to predict NKH model parameters, with seven input values, for fitting to target G/Gmax and D curves obtained from experimental results. Three soils were used with the nonlinear G/Gmax and damping curves obtained from the resonant column test (RCT). The NKH model parameters were estimated and compared with the SNKH model results. The procedure was also used for the rock materials, proving the application of the NKH model and the semi-automated procedure for both soils and rocks.

2. Soil Sample and ABAQUS Simulation of Torsional Shear Test

2.1. Soil Sample Properties

A total of three soils were simulated in the study. Two soils, namely, a silica and a mixture of weathered soil and pebbles, were tested with the resonant column test; another popular soil, Toyoura sand, was also used in this study with the fundamental properties reported by Wu, Yamamoto [11] and Kokusho [12], as shown in Table 1. The samples prepared in laboratory RCT tests were made with a diameter of 50 mm and height of 105 mm. The silica sand was carefully poured into the cylindrical sampler at a relative density of 80%; the mixture was carefully tamped at an optimum water content, a wopt of 10.5%, defined from compaction test to obtain the highest density, as well as stiffness. The mixed soil included weathered soil and small pebbles with a mixed ratio by weight of 55% and 45%, respectively. Table 1 shows the strength properties of the two soils at similar conditions to the RCT test. The small-strain shear wave velocity, shear modulus degradation (G/Gmax), and damping (D) of the two soils were defined by performing resonant column tests (RCTs) at a confining pressure of 200 kPa, an average value of confining strain when performing RCTs. Figure 1a,b show G/Gmax and the damping curves of the two soils. Figure 1c shows the G/Gmax and D curves of the Toyoura sand, which were reported by Wichtmann and Triantafyllidis [13].

2.2. Torsional Shear Test Modeling

To obtain the G/Gmax and damping curves by simulation, a torsional shear (TS) test was modelled in ABAQUS. A cylindrical column was made with the same dimension as the TS test column in the laboratory, as shown in Figure 2a. The end of the column was fixed, and a torsional shear deformation was applied at the top-free end of the column as a torsional rotation with a sinusoidal form (Figure 2b). The torsional strain for a cylindrical column specimen could be defined based on the applied rotational angle (θ), radius (r), and height (l) of the cylindrical column, as follows [14,15]:
ε = r e q × θ l
In this study, the equivalent radius req was taken to be 0.82 times greater than the radius of the cylindrical column [15]. The confining pressure was applied by surface pressure load and maintained during modelling to simulate a confining stress. The value of the confining pressure was chosen to be 200 kPa to produce a similar stress condition to the experimental RCT. The peak amplitude of the torsional shear deformation applied at the top-free end was changed in ten levels to obtain the shear strain of the column, which was defined from the maximum rocking angle, in the range of 0.0001% to ~0.1%. The torsional shear stresses of elements located outside the column were captured, and the hysteretic loop of shear stress–shear strain could be obtained for each strain level. The G/Gmax and damping value could be easily defined from the measured hysteretic loops.

3. Nonlinear Kinematic Hardening Model

3.1. Simplified Nonlinear Kinematic Hardening (SNKH) Model

The nonlinear kinematic hardening model has frequently been used for soil in dynamic analyses using ABAQUS [5,6,7,16,17]. In previous studies, a simplified nonlinear kinematic hardening (SNKH) model was applied with the implied von-Mises criterion and a back stress (α), which control the translation of yield stress. The yield criterion is defined as:
F = f (σα) − σo = 0
where σo is the initial yield stress at zero plastic strain. The increment of back stress is described as follows [5,18]:
α ˙ = C 1 1 σ o ( σ α ) ε ¯ ˙ p l γ 1 α ε ¯ ˙ p l
where ε ¯ ˙ p l is the increment of equivalent plastic strain. The three parameters required for the SNKH model are the initial yield stress (σo), initial kinematic hardening modulus (C1 is equal to Young’s modulus E), and rate of decrement of kinematic hardening modulus (γ1). The model parameters could be estimated based on the assumption of equaling the stress at large plastic strain to the maximum yield stress (σy) of soil, which can be estimated from the soil strength parameters (i.e., friction angle ϕ, and cohesion, C) [5,7]. The σo value can be estimated as a fraction (λ) of the yield stress (σo = λ·σy), and the γ1 parameter could be calculated as:
γ 1 = C 1 / ( σ y σ o )
By calibrating the λ parameter, which typically ranges from 0.1 to 0.3 for soil, the model parameters could be estimated to fit to the G/Gmax curve. The mean-average-error (MAE) was used to estimate the difference between the model and target curves, as:
MAE = | A m A t |
where Am and At are the G/Gmax values obtained from the SNKH model and target curve, respectively.
Figure 1a–c show the G/Gmax and D curves with SNKH model compared to the target curves with the model parameters indicated in Table 1. The λ parameter for the three soils varied (0.15 to 0.3). Note that the main point of the calibrating method was modifying the back-stress curve by changing the σo and γ1 parameters, without changing the ultimate shear strength, which is represented by summation of the σo and C1/γ1 values [5]. In this study, the λ parameter was chosen as the value that produced the lowest MAE, as shown in Figure 1d.
The degradation of shear modulus and damping properties simulated in ABAQUS by the SNKH model partially matched the target curves. However, there were several error points: (1) the threshold shear strain ( ε t h S N K H ), the strain level that separates material behavior from elastic to nonlinear plastic that can be estimated as shear strain at G/Gmax of 0.99 [19,20] in the model was much higher than those in experiments ( ε t h R C T ); (2) the G/Gmax curve at small strain range was not well simulated; and (3) the hysteretic damping in the SNKH model at low high strain level was much smaller than in experiment, as shown in Figure 1a–c. Similar errors could also be seen in the results reported by Tsinidis and Pitilakis [6] and Tsinidis and Pitilakis [7].

3.2. Nonlinear Kinematic Hardening (NKH) Model

In SNKH model, the G/Gmax and damping curves were not well modelled; therefore, in this study, the nonlinear kinematic hardening (NKH) model, which was basically developed from the SNKH model, was recommended. The NKH model has been a very common model used in mechanical engineering, and several researchers have proposed the procedure to estimate model parameters for metal materials based on experimental results [21,22,23,24]. Using a similar yield surface function and the von-Mises failure criterion to the SNKH model, the NKH model was upgraded to three kinematic hardening moduli (i.e., C1, C2, and C3), and three rates of decrement of kinematic hardening modulus (i.e., γ1, γ2, and γ3, respectively). With the initial yielding stress, σo, the model had seven input parameters, resulting in difficulties in finding appropriate parameters for ground materials. The back stress that included the three components, α1, α2, and α3 with their evolution could be estimated from each pair of kinematic parameters (i.e., C1-γ1, C2-γ2, and C3-γ3, respectively) by Equation (3).
Figure 3a shows the back-stress curves for silica sand in the SNKH and NKH models with the plastic shear strain using the equation proposed by Mahmoudi and Pezeshki-Najafabadi [21]. In the NKH model, the total back stress is the summation of three back-stress components. It can be seen that with the increasing γpeq or shear strain, each back-stress component including α1, α1, and α3, respectively, sequentially prevail in the total back-stress (α in NKH model), and the total back-stress curve can be divided into three segments by the limitation of each back stress component. Therefore, each pair of Ci and γi parameters controls the G/Gmax curve at each level of shear strain, as indicated in Figure 3b. A manual process was performed to estimate the seven input parameters for the NKH model. Firstly, the initial yield stress, σo, was estimated from G and the threshold shear strain ( ε t h R C T ), which is defined from the G/Gmax curve. Compared to those in the SNKH model, the σo value in the NKH model was slightly reduced, because the SNKH model produced a larger threshold strain. Due to the reduction in the σo, three pairs of kinematic hardening parameters were firstly chosen to have slightly higher back stress at all strain range (Figure 3a). For each trial, the G/Gmax curve was obtained, and each pair of kinematic hardening parameters was modified following the instruction indicated in Figure 3b. By variation in its slope, the target G/Gmax curves are divided into three parts, which are controlled by each pair of kinematic parameters, C1-γ1, C2-γ2, and C3-γ3, respectively. For each part of the G/Gmax curve with its pair of kinematic hardening parameters, the reduction in C and increment in γ resulted in the increase in slope of the G/Gmax curve, and vice versa, as shown in Figure 3b with blue and green arrows. This is because the reduction in C and/or increase in γ in each segment results in reduction in limitation of back stress; therefore, G reduces in that segment.
The estimation process for kinematic hardening parameters should be performed in order from C1-γ1 to C2-γ2, and finally C3-γ3; and as a result, the G/Gmax curve could gradually be fitted to the experimental curve. Figure 4 illustrates the results of the nonlinear curve by using the NKH model with parameters for the three soils obtained after approximately 20 trials with a considerable time consumption. The NKH model perfectly simulated the G/Gmax curve at all ranges of shear strain level. In particular, the results for damping obtained from the NKH model were significantly improved compared to the SNKH model in Figure 1a,b, due to the reasonable threshold shear strain, and the back-stress curves with more kinematic hardening parameters. Anastasopoulos and Gelagoti [5] reported the requirements for the yield stress of material by controlling back-stress at large strain, which were implemented in the SNKH model. Figure 3a typically shows the back-stress curve of silica sand in the NKH model, and at extremely large shear strain, the back-stress value was slightly higher than that in the SNKH model, because of the small reduction in the initial yield stress, σo. Therefore, the requirement of the ultimate shear strength of soil could be allowed.

4. Semi-Automated Procedure to Estimate the NKH Model Parameters for Soils and Rocks

4.1. MATLAB Algorithm Estimating Shear Stress–Shear Strain Loop

Figure 5 shows the algorithm of the program based on the ABAQUS calculation procedures [18]. For a given incremental shear strain Δεk, the von-Mises yield criterion was initially checked with the original predicted elastic stress (Δσprk) by the assumption of pure elastic strain. Then, the iteration process was applied to determine the correct incremental plastic strain (Δεp), incremental back stress (Δα), and incremental total stress (Δσ) at a step. In the iteration process, new plastic strain was estimated based on the flow rule [18], and the back-stress values either for the positive (αp) or negative (αn) phase could be calculated as follows [21]:
{ α p = i = 1 3 [ 2 3 C i γ i + ( α i n 2 3 C i γ i ) e x p ( γ i Δ ε p ) ] ( p o s i t i v e   p h a s e ) α n = i = 1 3 [ 2 3 C i γ i + ( α i p + 2 3 C i γ i ) e x p ( γ i Δ ε p ) ] ( n e g a t i v e   p h a s e )
where αip and αin are the maximum and minimum values of the back-stress component during positive and negative loading, respectively, and Δεp is the incremental plastic strain at the step. The total shear stress could be calculated from the incremental elastic strain (Δεe), which was calculated from the incremental strain and plastic strain at that step (Δεe = Δεk − Δεp), and the yield function, F, could be defined following Equation (2). The process of iteration finishes when F reaches the tolerated value [18].
The procedure was implemented in MATLAB, which can estimate shear stress–strain loops with a much faster speed compared to ABAQUS. The results obtained from this procedure were then validated by comparing the stress, hysteretic loops, G/Gmax, and damping curves obtained in the ABAQUS simulation. Figure 6a,b show a comparison of the stresses obtained in ABAQUS and MATLAB with 10 cycles of sinusoidal shear strains, and the identical results between hysteretic loops obtained in both ABAQUS and MATLAB were typical for the silica and Toyoura sands. The shear modulus (G) and damping (D) were then estimated and compared with the ABAQUS results for the silica and Toyoura sands, as shown in Figure 6c. Similar results were also obtained for the mixture soil. The good matching of the results in both the G/Gmax and D curves proves that the algorithm could predict the soil nonlinear properties in a wide range of shear strain of 0.0001% to 1%.

4.2. Semi-Automated Procedure to Estimate the NKH Model Parameters

In the previous section, the NKH model parameters, including the initial yield stress σo and three pairs of kinematic hardening parameters (i.e., C1-γ1, C2-γ2, and C3-γ3), were found for the three soils by manually using the trial-and-error method, which consumed a lot of time and effort. Therefore, in this section, an algorithm is proposed with the main purpose of rapidly finding the NKH model parameters to match the target G/Gmax and damping curves obtained from experimental results based on the basic properties of soil, such as the Young’s modulus (E) or shear modulus (G), density (ρ), and Poisson’s ratio (ν). Figure 7 shows the flowchart of the semi-automated procedure estimating the NKH model parameters. The outputs of the process are the seven parameters of the NKH model previously mentioned from the input, including E or G, ρ, and ν, and the target G/Gmax and damping curves. The procedure could be divided into two parts: (1) finding the initial parameters from the G/Gmax curve; and (2) updating the parameters with consideration of both the G/Gmax and damping curves with a weighting factor (w).
There are three steps to estimating the initial parameter from the G/Gmax curve. The first parameter, initial yield stress, σo, was estimated from G and the threshold shear strain ( ε t h R C T ), as shown in Figure 7, following the method mentioned above. At a shear strain larger than ε t h R C T , the G of the material was reduced compared to Gmax, and plastic strain and back stress are produced. In the SNKH and NKH models, the yield stress at a strain level is the summation of the initial yield stress and the back stress. Therefore, the back-stress target curve with the plastic shear strain could be estimated for a target G/Gmax curve in the second step, as shown in Figure 7. The target G/Gmax curve was discretized to approximately 200 data points in the logarithm scale for the shear strain space. At the shear strain εk with increment shear strain Δεk, the shear stress, σ, could be estimated by the shear modulus G calculated from the G/Gmax value. Plastic shear stress could be defined from the shear stress and the elastic shear stress, which is equal to σ/Gmax. Back stress could also be calculated from the σ and σo values. Finally, by applying the equation of evolution of back stress with plastic strain (i.e., the positive phase of Equation (6)), the initial parameters (i.e., C1, γ1, C2, γ2, C3, and γ3) could be predicted. Figure 8 compares the G/Gmax and damping curves with the initial parameters with the target curves for the different materials. The results show slight differences between the G/Gmax curve obtained by the model with the target in the shear strain range (0.0005% to 0.01%), especially for the Toyoura sand in Figure 8c. Therefore, these initial parameters needed to be upgraded.
In the second part of the procedure, the model parameters were modified to better fit the G/Gmax and/or damping curves. The G/Gmax and damping curves were divided into two or three segments by manually assigning the strain level at a slope-changing point, as shown in Figure 3b. Then, each pair of kinematic hardening parameters, including C1-γ1, C2-γ2, and C3-γ3, respectively, were modified with a reducing changing rate to fit each part of the G/Gmax and/or damping curves. The changing rate, which could be manually input into MATLAB, indicates the change in the percent of the parameter (i.e., 20%, 10%, and 5%), and the degree of change should depend on the MAE value. The weighting factor (w) was used to estimate the error considering both G/Gmax and damping with MAEG and MAED, respectively, as follows:
MAE = w·MAEG + (1 − w)·MAED
Note that in fitting to the damping curve, the new target damping was considered by subtracting the damping (i.e., damping obtained from the RCT) to the initial damping value, which was defined as damping at the threshold strain. New back-stress parameters for each segment of G/Gmax or damping curves could be obtained with the lowest MAE value at that segment, following the method proposed by Anastasopoulos and Gelagoti [5]. The convergence of the procedure could be reached if the MAE value was stable, with slight modification in the rate of change. Figure 9a–c show the results of the G/Gmax and damping curves for the three soils with parameters obtained by the proposed procedure with w factors of 1, which means that only G/Gmax was considered during the modification process. Acceptable results were obtained in G/Gmax compared between the model and target curves, which show considerable improvement compared to the initial parameters in Figure 8a–c. By performing the algorithm in the MATLAB compiler, the results of the parameters were obtained in a very short time, compared to the manual trial-and-error method mentioned above. The results for the damping improved, even though the underrated damping was seen at small shear strain and damping was overestimated at large shear strains. Figure 9d compares the back-stress curves of Toyoura sand for three cases, including the SNKH model, the NKH model with parameters obtained from the semi-automated procedure (NKH_G), and from modifying the SNKH curves through trial-and-error (NKH_Trial). The back-stress value in large plastic strain value obtained by the semi-automated method is reasonable compared to those of the trial-and-error method and the SNKH model. The same results were also obtained for the silica sand and mixture soil. It could be concluded that the NKH model and its parameters obtained by the algorithm match the requirement for the shear strength of the material proposed by Anastasopoulos and Gelagoti [5].
Figure 10 shows the results of the shear modulus degradation and damping curves for the three soils with parameters obtained by the semi-automated procedure with a w factor of 0.95, which means that both the G/Gmax and damping were considered. The results for damping were much improved for all materials compared to those in Figure 9a–c with acceptable shear modulus degradation curves. The value of the weighting factor could be increased to improve the accuracy of the model in simulating the shear modulus of the material.

4.3. Application of the NKH Model and the Semi-Automated Procedure for Rocks

Although the application of the nonlinear model for rock material during dynamic analyses has recently been seen frequently with the increase in the structure’s size and weight, such as in nuclear power plants (NPPs), high-rise buildings, and gravity dams, a few nonlinear models did exist for rocks in ABAQUS. Drucker–Prager, an isotropic hardening model recommended for rock materials, requires yield stress and plastic strain data, which have to calibrated with complicated techniques [25]. Moreover, in rock engineering, the Young’s modulus, shear wave velocity, Poisson’s ratio, and rock quality are the common properties provided during the geotechnical survey, whereas the cohesion and friction angle are not very common. Nonlinear dynamic properties for intact or jointed rocks have been proposed in several reports [26,27,28,29]. EPRI [30] proposed referent G/Gmax and damping curves for rock material, depending on its depth. Therefore, by using the semi-automated procedure proposed above, the parameters for the NKH model can be found for rock based on common fundamental properties and referent dynamic curves.
Table 2 shows the ground profile for Wolsung NPP located in Korea. There were 11 rock layers with shear wave velocities varying from about 549 m/s at the surface to around 2804 m/s at a depth of 100 m. The semi-automated procedure was applied to define the NKH model parameters for layers 3 and 5 at depths of 6 m and 17 m, respectively; based on these depths, the referent G/Gmax and damping curves are presented in EPRI [30]. Figure 11 shows the G/Gmax and damping curves for finding the NKH model for the two rocks, considering only the shear modulus (i.e., the weighting factor, w, is 1). Good results were obtained for G/Gmax for the two rocks, proving the application of the NKH model with the semi-automated procedure for both soil and rock materials.

5. Conclusions

The nonlinear kinematic hardening (NKH) model was used in this study to simulate the nonlinear strain-dependent dynamic properties of soils and rocks. Seven model parameters were estimated either by trial-and-error or by using a semi-automated procedure with the modeling of torsional shear test in either ABAQUS or MATLAB. The main conclusions are as follows:
  • The simplified nonlinear kinematic hardening (SNKH) model could simulate the strain-dependent nonlinear properties of soil with a simple procedure of calibrating model parameters using the yield strength of the material. However, limitations of the model were illustrated, such as the threshold shear strain, shear modulus, and damping at low strain level;
  • Parameters for the nonlinear kinematic hardening (NKH) model were estimated by using the proposed method of estimation based on the G/Gmax curve with initial values defined from the SNKH back-stress curve. The NKH model with seven parameters improved on the limitations of the SNKH model, and its simulated nonlinear shear modulus degradation and damping curves of soil were much better than the frequently used SNKH model;
  • The trial-and-error procedure was proposed to estimate NKH model parameters, with one drawback of considerable time consumed, which has resulted in the very limited implementation of the NKH model in geotechnical engineering;
  • A semi-automated procedure was proposed to estimate the NKH model parameters based on the basic properties of soil and the nonlinear dynamic curves. By implementation with the help of MATLAB, the parameters of the NKH model could easily be estimated in a very short time, and considering either the shear modulus degradation or the damping of material. The procedure was also applied for rock material, with good results being obtained;
  • In this study, the NKH model could accurately simulate the shear modulus degradation, but the damping was not well simulated, especially at the high shear strain level, which would be not appropriate for analyses of a strong earthquake. Moreover, a more automated procedure should be further developed in combination with ABAQUS for estimating NKH model parameters which would be more applicable for researchers and engineers. Additionally, a built-in model combined with yielding criteria for geotechnical material needs to be developed further for ABAQUS that has been available in other commercial programs.

Author Contributions

Conceptualization, J.-M.K., C.L. and V.-L.N.; methodology, V.-L.N. and E.-h.L.; software, V.-L.N. and E.-h.L.; validation, V.-L.N. and C.L.; formal analysis, V.-L.N.; investigation, V.-L.N. and E.-h.L.; resources, J.-M.K.; data curation, V.-L.N. and E.-h.L.; writing—original draft preparation, V.-L.N.; writing—review and editing, V.-L.N. and C.L.; visualization, V.-L.N. and C.L.; supervision, C.L. and J.-M.K.; project administration, J.-M.K.; funding acquisition, J.-M.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Ministry of Trade, Industry & Energy (MOTIE) of Korea (No.20193110100020).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data sharing is not applicable to this article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lashin, I.; Hussein, M.N.; Karray, M. Assessment of Conventional Interpretation Methods of RC Results Based on 3D Numerical Simulations. Int. J. Geomech. 2018, 18, 04018160. [Google Scholar] [CrossRef]
  2. Rayhani, M.H.; El Naggar, M.H. Numerical Modeling of Seismic Response of Rigid Foundation on Soft Soil. Int. J. Geomech. 2008, 8, 336–346. [Google Scholar] [CrossRef]
  3. Rayhani, M.T.; El Naggar, M.H. Physical and Numerical Modeling of Seismic Soil-Structure Interaction in Layered Soils. Geotech. Geol. Eng. 2012, 30, 331–342. [Google Scholar] [CrossRef]
  4. Chen, G.; Ruan, B.; Zhao, K.; Chen, W.; Zhuang, H.; Du, X.; Khoshnevisan, S.; Juang, C.H. Nonlinear Response Characteristics of Undersea Shield Tunnel Subjected to Strong Earthquake Motions. J. Earthq. Eng. 2020, 24, 351–380. [Google Scholar] [CrossRef]
  5. Anastasopoulos, I.; Gelagoti, F.; Kourkoulis, R.; Gazetas, G. Simplified Constitutive Model for Simulation of Cyclic Response of Shallow Foundations: Validation against Laboratory Tests. J. Geotech. Geoenviron. Eng. 2011, 137, 1154–1168. [Google Scholar] [CrossRef] [Green Version]
  6. Tsinidis, G.; Pitilakis, K.; Trikalioti, A.D. Numerical simulation of round robin numerical test on tunnels using a simplified kinematic hardening model. Acta Geotech. 2014, 9, 641–659. [Google Scholar] [CrossRef]
  7. Tsinidis, G.; Pitilakis, K.; Anagnostopoulos, C. Circular tunnels in sand: Dynamic response and efficiency of seismic analysis methods at extreme lining flexibilities. Bull. Earthq. Eng. 2016, 14, 2903–2929. [Google Scholar] [CrossRef]
  8. Kammerer, A.; Coleman, J.; Hashash; Johnson, J.; Kennedy, R.; Whittaker, A.; Yu, C.-C. Nonlinear Soil-Structure-Interaction Analysis in Support of Seismic Design and Probabilistic Risk Assessment of Nuclear Facilities; Idaho National Laboratory: Idaho Falls, ID, USA, 2018.
  9. Ma, Z.; Liao, H.; Dang, F.; Cheng, Y. Seismic slope stability and failure process analysis using explicit finite element method. Bull. Eng. Geol. Environ. 2020, 80, 1287–1301. [Google Scholar] [CrossRef]
  10. Li, F.; Chen, G. Nonlinear Seismic Response Characteristics of CAP1400 Nuclear Island Structure on Soft Rock Sites. Sci. Technol. Nucl. Install. 2020, 2020, 8867026. [Google Scholar] [CrossRef]
  11. Wu, Y.; Yamamoto, H.; Yao, Y. Numerical study on bearing behavior of pile considering sand particle crushing. Geomech. Eng. 2013, 5, 241–261. [Google Scholar] [CrossRef]
  12. Kokusho, T. Cyclic Triaxial Test of Dynamic Soil Properties for Wide Strain Range. Soils Found. 1980, 20, 45–60. [Google Scholar] [CrossRef]
  13. Wichtmann, T.; Triantafyllidis, T. Effect of Uniformity Coefficient on G/Gmax and Damping Ratio of Uniform to Well-Graded Quartz Sands. J. Geotech. Geoenviron. Eng. 2013, 139, 59–72. [Google Scholar] [CrossRef] [Green Version]
  14. Paterson, M.S.; Olgaard, D.L. Rock deformation tests to large shear strains in torsion. J. Struct. Geol. 2000, 22, 1341–1358. [Google Scholar] [CrossRef]
  15. Kim, D.S.; Stokoe, K.H.; Hudson, W.R. Deformational Characteristics of Soils at Small to Intermediate Strains from Cyclic Tests. Ph.D. Thesis, Center for Transportation Research, Bureau of Engineering Research, University of Texas at Austin, Austin, TX, USA, August 1991. [Google Scholar]
  16. Gazetas, G.; Adamidis, O.; Kontourupi, T. Nonlinear rocking stiffness of foundations. Soil Dyn. Earthq. Eng. 2013, 47, 83–91. [Google Scholar] [CrossRef]
  17. Kourkoulis, R.S.; Lekkakist, P.C. Suction caisson foundations for offshore wind turbines subjected to wave and earthquake loading: Effect of soil–foundation interface. Géotechnique 2014, 64, 171–185. [Google Scholar] [CrossRef] [Green Version]
  18. Simulia. Abaqus/Standard Version 6.18 Analysis User’s Manual; Pawtucket: Rhode Island, NJ, USA, 2018. [Google Scholar]
  19. Ichii, K.; Mikami, T. Cyclic threshold shear strain in pore water pressure generation in clay in situ samples. Soils Found. 2018, 58, 756–765. [Google Scholar] [CrossRef]
  20. Lee, C.; Shin, H.; Lee, J.-S. Behavior of sand–rubber particle mixtures: Experimental observations and numerical simulations. Int. J. Numer. Anal. Methods Geomech. 2014, 38, 1651–1663. [Google Scholar] [CrossRef]
  21. Mahmoudi, A.H.; Pezeshki-Najafabadi, S.M.; Badnava, H. Parameter determination of Chaboche kinematic hardening model using a multi objective Genetic Algorithm. Comput. Mater. Sci. 2011, 50, 1114–1122. [Google Scholar] [CrossRef]
  22. Moslemi, N.; Stokoe, K.H.; Hudson, W.R. Evaluation of Sensitivity and Calibration of the Chaboche Kinematic Hardening Model Parameters for Numerical Ratcheting Simulation. Appl. Sci. 2019, 9, 2578. [Google Scholar] [CrossRef] [Green Version]
  23. Hales, R.; Rayhani, M.H.; El Naggar, M.H. A Code of Practice for the determination of cyclic stress-strain data. Mater. High Temp. 2002, 19, 165–185. [Google Scholar] [CrossRef]
  24. Rezaiee-Pajand, M.; Sinaie, S. On the calibration of the Chaboche hardening model and a modified hardening rule for uniaxial ratcheting prediction. Int. J. Solids Struct. 2009, 46, 3009–3017. [Google Scholar] [CrossRef] [Green Version]
  25. Atrian, A.; Hussein, M.N.; Karray, M. A novel approach to calibrate the Drucker–Prager Cap model for Al7075 powder. Arch. Appl. Mech. 2018, 88, 1859–1876. [Google Scholar] [CrossRef]
  26. Perino, A.; Barla, G. Resonant Column Apparatus Tests on Intact and Jointed Rock Specimens with Numerical Modelling Validation. Rock Mech. Rock Eng. 2015, 48, 197–211. [Google Scholar] [CrossRef]
  27. Sebastian, R.; Sitharam, T.G. Column Tests and Nonlinear Elasticity in Simulated Rocks. Rock Mech. Rock Eng. 2018, 51, 155–172. [Google Scholar] [CrossRef]
  28. Hosseini, R.; Paterson, M.S.; Olgaard, D.L. Free-Free Resonant Column Testing of Rock Cores from Two Spillways and Adjacent Exposed Rock Areas of a Dam. In Geo-Congress 2020: Geotechnical Earthquake Engineering and Special Topics; American Society of Civil Engineers: Reston, VA, USA, 2020; pp. 409–417. [Google Scholar]
  29. Chong, S.; Kim, J.; Cho, G. Rock Mass Dynamic Test Apparatus for Estimating the Strain-Dependent Dynamic Properties of Jointed Rock Masses. Geotech. Test. J. 2014, 37, 311–318. [Google Scholar] [CrossRef]
  30. EPRI. Seismic Evaluation Guidance: Screening, Prioritization and Implementation Details (SPID) for the Resolution of Fukushima Near-Term Task Force Recommendation 2.1: Seismic; Rep. No. 1025287; EPRI: Palo Alto, CA, USA, 2012. [Google Scholar]
Figure 1. G/Gmax and damping curve obtained from RCTs and the SNKH model for (a) silica sand; (b) mixture soil; and (c) Toyoura sand; and (d) MAE value for three soils. The arrows show threshold strain in (ac), and the lowest lambda value in (d).
Figure 1. G/Gmax and damping curve obtained from RCTs and the SNKH model for (a) silica sand; (b) mixture soil; and (c) Toyoura sand; and (d) MAE value for three soils. The arrows show threshold strain in (ac), and the lowest lambda value in (d).
Applsci 11 08611 g001aApplsci 11 08611 g001b
Figure 2. (a) Cylindrical column, and (b) torsional shear strain.
Figure 2. (a) Cylindrical column, and (b) torsional shear strain.
Applsci 11 08611 g002
Figure 3. (a) Back-stress curves for silica sand, and (b) the estimation of NKH model parameters for each part of the G/Gmax curve.
Figure 3. (a) Back-stress curves for silica sand, and (b) the estimation of NKH model parameters for each part of the G/Gmax curve.
Applsci 11 08611 g003
Figure 4. G/Gmax and damping curves in the NKH model for (a) silica sand, (b) mixture soil, and (c) Toyoura sand.
Figure 4. G/Gmax and damping curves in the NKH model for (a) silica sand, (b) mixture soil, and (c) Toyoura sand.
Applsci 11 08611 g004
Figure 5. Algorithm for estimating back stress and plastic strain in MATLAB.
Figure 5. Algorithm for estimating back stress and plastic strain in MATLAB.
Applsci 11 08611 g005
Figure 6. Comparison of (a) shear stress estimation, (b) shear stress–shear strain hysteretic loop, and (c) G/Gmax and damping curves between ABAQUS and MATLAB typically for silica sand (left) and Toyoura sand (right).
Figure 6. Comparison of (a) shear stress estimation, (b) shear stress–shear strain hysteretic loop, and (c) G/Gmax and damping curves between ABAQUS and MATLAB typically for silica sand (left) and Toyoura sand (right).
Applsci 11 08611 g006
Figure 7. Flowchart for estimating NKH model parameters.
Figure 7. Flowchart for estimating NKH model parameters.
Applsci 11 08611 g007
Figure 8. G/Gmax and damping curves with initial parameters of the NKH model for (a) silica sand, (b) mixture soil, and (c) Toyoura sand.
Figure 8. G/Gmax and damping curves with initial parameters of the NKH model for (a) silica sand, (b) mixture soil, and (c) Toyoura sand.
Applsci 11 08611 g008
Figure 9. G/Gmax and damping curves with parameters of the NKH model (weighting factor, w = 1) for (a) silica sand, (b) mixture soil, and (c) Toyoura sand, and (d) back-stress curve for Toyoura sand.
Figure 9. G/Gmax and damping curves with parameters of the NKH model (weighting factor, w = 1) for (a) silica sand, (b) mixture soil, and (c) Toyoura sand, and (d) back-stress curve for Toyoura sand.
Applsci 11 08611 g009
Figure 10. G/Gmax and damping curves with parameters of the NKH model (weighting factor, w = 0.95) for (a) silica sand, (b) mixture soil, and (c) Toyoura sand.
Figure 10. G/Gmax and damping curves with parameters of the NKH model (weighting factor, w = 0.95) for (a) silica sand, (b) mixture soil, and (c) Toyoura sand.
Applsci 11 08611 g010
Figure 11. G/Gmax and damping curves with the NKH model for two rocks in NPP profile: (a) layer 3 and (b) layer 5.
Figure 11. G/Gmax and damping curves with the NKH model for two rocks in NPP profile: (a) layer 3 and (b) layer 5.
Applsci 11 08611 g011
Table 1. Soil properties and parameters of SNKH and NKH models.
Table 1. Soil properties and parameters of SNKH and NKH models.
PropertiesSilica SandMixtureToyoura Sand
Specific gravity, Gs2.652.672.65
Dry density (t/m3)1.541.821.62
Young modulus, E (MPa)74.26284.86291.52
Poison ratio0.2430.2370.227
Cohesion (kN/m2)0.10.10.1
Friction angle (deg)36.638.544
Shear wave velocity (m/s)139.27250.15270.89
SNKH model
Initial yield stress, σo (kPa)15.4933.7030.08
C1 (kPa)74,257.42284,861.68291,520.28
γ1845.962113.253230.54
NKH model (weighting factor w = 1)
σo (kPa)1.566.170.36
C11,007,085.811,826,119.563,901,266.00
γ1101,045.1579,421.78132,583.39
C2107,118.11371,966.54363,633.29
γ24774.634700.15574.23
C3098,866.0127,380.88
γ3111
NKH model (weighting factor w = 0.95)
C1652,591.601,725,682.992,808,911.52
γ161,112.1154,324.5057,276.02
C2107,118.11361,551.47261,815.97
γ23819.715414.5630,10.08
C3067,624.3537,380.88
γ3111
Table 2. Fundamental property profiles for Wolsung NPP.
Table 2. Fundamental property profiles for Wolsung NPP.
Layer
No
Depth
(m)
Vs
(m/s)
υρ
(t/m3)
G
(GPa)
E
(GPa)
13549.00.4182.0810.631.78
29554.60.4162.0810.641.81
3151010.60.3612.1452.195.96
42111520.352.1612.877.74
5271173.30.352.1612.978.03
6331781.90.3562.3217.3719.99
7391798.50.3372.3217.5120.08
8451815.20.3382.3217.6520.46
9511832.10.3382.3217.7920.85
10571848.80.3392.3217.9321.25
11632804.20.332.48119.5151.90
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ngo, V.-L.; Lee, C.; Lee, E.-h.; Kim, J.-M. Semi-Automated Procedure to Estimate Nonlinear Kinematic Hardening Model to Simulate the Nonlinear Dynamic Properties of Soil and Rock. Appl. Sci. 2021, 11, 8611. https://doi.org/10.3390/app11188611

AMA Style

Ngo V-L, Lee C, Lee E-h, Kim J-M. Semi-Automated Procedure to Estimate Nonlinear Kinematic Hardening Model to Simulate the Nonlinear Dynamic Properties of Soil and Rock. Applied Sciences. 2021; 11(18):8611. https://doi.org/10.3390/app11188611

Chicago/Turabian Style

Ngo, Van-Linh, Changho Lee, Eun-haeng Lee, and Jae-Min Kim. 2021. "Semi-Automated Procedure to Estimate Nonlinear Kinematic Hardening Model to Simulate the Nonlinear Dynamic Properties of Soil and Rock" Applied Sciences 11, no. 18: 8611. https://doi.org/10.3390/app11188611

APA Style

Ngo, V. -L., Lee, C., Lee, E. -h., & Kim, J. -M. (2021). Semi-Automated Procedure to Estimate Nonlinear Kinematic Hardening Model to Simulate the Nonlinear Dynamic Properties of Soil and Rock. Applied Sciences, 11(18), 8611. https://doi.org/10.3390/app11188611

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