Optimization of Two Soil–Structure Interaction Parameters Using Dynamic Centrifuge Tests and an Analytical Approach

The response of the structure subjected to an earthquake load is greatly affected by the properties of the structure and soil so it is very important to accurately determine the characteristics of the structure and soil for analysis. However, studies on the effective profile depth where soil properties are determined, have been conducted in the presence of restricted conditions (i.e., surface foundation, linear soil properties), and without any considerations on damping. In case of the effective height of structure that affects its rocking behavior, it was only theoretically or empirically determined. In addition, most previously published studies on soil–structure interaction (SSI) focused on limited effects and parameters (e.g., rocking behavior, embedment effect, effective profile depth, spring constant, and damping coefficient) and not on comprehensive SSI parameters. Furthermore, no detailed validation procedure has been set in place which made it difficult to validate the SSI parameters. Since the effective height of structure and effective profile depth are the basis of all the input parameters of SSI analysis, it is important to validate and determine them. Therefore, in this study, the procedure used to optimize the two SSI parameters was established based on an analytical approach that considered all the possible SSI parameters that were investigated from conventional codes and studies and physical model tests. As a result of this study, the optimum values of the effective height of the structure and effective profile depth were respectively determined according to (a) the height from the bottom part of the foundation to the center of the mass of the superstructure, and according to (b) the depth at values equal to four times the radius of the foundation.


Introduction
The dynamic characteristics of the structure depend on the surrounding soil. To estimate an accurate structural response for seismic design, the soil-structure interaction (SSI) effect has been considered as the crucial effect for the seismic evaluation of the structure. There have been various studies on the soil-structure interaction (SSI) effect but the subject of each study was not comprehensive but limited to specific parameters or phenomena (e.g., spring constant and damping coefficient [1][2][3][4][5][6], effective profile depth [7], and rocking behavior [8][9][10][11]). With regard to the SSI analysis, four procedures were introduced in FEMA356 [12]: Linear static procedure (LSP), linear dynamic procedure (LDP), nonlinear static procedure (NSP), and nonlinear dynamic procedure (NDP). Among them, the static procedures (i.e., LSP and NSP) have been accepted in various standards owing to their simplicity and (NDP). Among them, the static procedures (i.e., LSP and NSP) have been accepted in various standards owing to their simplicity and practicality [12][13][14][15], and the dynamic procedures (i.e., LDP and NDP) have been used to either verify the static procedures or to obtain more detailed structural responses. Even though all the standards mentioned above (i.e., FEMA 356, ATC-40, FEMA 440, and ASCE 41-13) define the same three types of modes (i.e., structural swaying, foundation swaying, and rocking), the formula for the soil-spring constant and damping coefficient, and the two SSI parameters (i.e., the effective height of the structure and effective profile depth) were defined differently at various standards and studies, including the aforementioned standards [1][2][3][4]7,12,14,[16][17][18][19]. The two SSI parameters are very important because they determine the dynamic soil properties and the rocking potential of the soil-structure system (i.e., rocking damping coefficient and moment), but they have been derived in restricted conditions (i.e., the effective profile depth was determined based only on considerations of static soil stiffness in surface foundation condition [5,7]) and theoretically determined (i.e., the effective height was determined based on structural dynamic theory [3,10,12,13]). In addition, existing standards and prior studies have been limited in view of the following: (1) Approximate consideration of the nonlinear deformation characteristics of the soil based on peak ground acceleration [12,14,15,18], (2) determination of damping ratio of soil based only on radiation damping considerations [14,15,18], and (3) lack of detailed SSI analysis procedure and SSI parameter validation procedure with physical model tests. Therefore, in this study, the optimum SSI parameter selection procedure was established using an LDP-based analytical approach and relevant tests, whereby the SSI effects were comprehensively and appropriately considered (i.e., consideration of all the possible SSI parameters for analysis, accurate consideration of soil nonlinearity using site response analysis, and determination of total soil damping coefficient, including soil material damping). Finally, two optimum SSI parameters were determined based on a number of analyses according to the established procedure.

Soil-Structure Interaction (SSI) Parameters and Optimum SSI Parameter Selection Procedure
The SSI model is shown in Figure 1. It consists of two lumped masses and , whereby the distance between the two masses is defined as the effective height ℎ, a set of spring , and dashpot to represent the dynamic characteristics of the structure, and two sets of springs and dashpots to simulate the flexibility and damping of the soil. In the case of horizontal excitation, the total absolute displacement of the can be defined with three deformation modes, and is expressed as = + + h∅ + . In this case, is the absolute displacement of the soil, is the relative displacement between the foundation and the soil, ∅ is the rotation angle of the foundation, ℎ is the effective height of the structure, and is the relative displacement between the upper structure and the soil, or the net structural displacement. To obtain detailed structural responses, the equation of motion (EOM) of the soil-structure system that considers the aforementioned three degree-of-freedom (3DOF) is expressed by Equation (1) [8,9,20].

Soil-Spring Constants and Damping Coefficients
In the SSI system, the soil-spring constant and damping coefficient represent the characteristics of the foundation and lower ground, and impose a considerable influence on the system's response. Regarding the soil-spring constant and damping coefficient, Wolf [1] and Richart and Lysmer [2] suggested the use of the frequency-independent spring constant and radiation damping coefficient formulas that comprised simple parameters. These formulas have been used previously [14,19,21]. Gazetas [3] and Stewart et al. [4] proposed a precise frequency-dependent spring constant and damping coefficient formula that considered the soil impedance. However, to utilize the formulas, it is necessary to consider numerous parameters and numerous corrections for each frequency subject to earthquake loading. In FEMA356 [12], the frequency-independent soil stiffness and correction factor formulas-that considered the geometry and embedment depth of the rectangular footing in detail-were presented, but the damping coefficient formula was not presented. Therefore, two types of simple frequency-independent soil-spring constant and radiation damping coefficient formulas have been presented in Table 1 (i.e., the formulas suggested by Wolf [1] and Richart and Lysmer [2]) were adopted in this study.
Where B is the width of the foundation, L is the length of the foundation, G is the shear modulus of the soil, υ is the Poisson's ratio of the soil, and r x and r φ are the equivalent radii of the foundation expressed as √ BL/π and 4 BL 3 /3π, respectively. The parameters β x and β φ are a function of L/B, and the values in case of square footing (i.e., L/B = 1) can be obtained from the suggested chart, where β x and β φ are equal to 1 and 0.5, respectively.

Strain-Dependent Soil Properties
Linear SSI analyses (i.e., LSP and LDP) require constant soil properties regardless of depth and strain, but the actual shear wave velocity of soil V s (V s ) and damping ratio of soil ξ g (ξ g ) vary with depth and strain level. In conventional standards, the reduction factor for V s is determined and used to convert the strain dependent V s to a strain-independent factor, but a reduction factor is expressed as a function of the peak ground acceleration at the ground surface without considering the frequency contents of motion and depth. Moreover, there is no correction factor for ξ g in the standards [12,14,15,18]. To overcome the above limitation and to obtain strain-independent V s and ξ g values (i.e., equivalent V s and ξ g values) accurately during the excitation, one-dimensional equivalent linear-site-response analyses (SRAs) were performed in this study. Although equivalent V s and ξ g values are converted to their strain-independent forms, they are still non-uniformly distributed as a function of depth. Therefore, it is necessary to define an effective profile depth (Z p ). In this way, depth-independent, equivalent V s and ξ g values have to be obtained as the average values within a depth Z p . The averaged values of V s and ξ g that consider Z p can be obtained by Equation (2), as follows, where (V s ) i is the shear wave velocity of the ith soil layer, ξ g i is the material damping ratio of the ith soil layer, and ∆Z i is the thickness of the ith soil layer, respectively. Stewart et al. [7] regarded the static soil-spring constants obtained from the impedance solutions by Wong and Luco [5] as reference values, and repeatedly calculated the soil-spring constants at various profile depths. Note that Z p is 0.75 times the radius of the foundation (r), wherein the residual between the reference and the calculated value is minimized. In addition to prior research publications, in the recommended provision of national earthquake hazard reduction program (NEHRP) [16], 4 r and 1.5 r were proposed as the respective values of Z p for swaying and rocking behaviors, respectively. However, in previous studies, the embedment effect and damping ratio of soil were not considered. Therefore, in this study, three types of scenarios of 0.75 r, 2 r, and 4 r, were considered to evaluate the optimum effective profile depth and necessary considerations (i.e., embedment effect and ξ g ), and were included in the analytical approach for the evaluation of the optimum effective profile depth.

Embedment Correction for Spring Constant and Radiation Damping Coefficient
When the foundation was embedded, the natural frequency of the soil-structure system increased owing to an increase in the soil stiffness, and the overall response of the system decreased owing to an increase in radiation damping [22]. In view of this phenomenon, the formulas of the correction factor that considered the geometry and embedment depth of the rectangular footing in detail was presented in FEMA356 [12]. In contrast, Whitman [6] proposed a simple embedment correction factor (i.e., correction factor for both stiffness and damping) as a function of the embedment depth (D), radius of foundation (r), and Poisson's ratio (ν) of the soil. Table 2 shows the correction factor for embedment suggested by Whitman. According to the Whitman's suggestion, the damping correction factor was suggested for the radiation damping ratio of the soil but the damping matrix in Equation (1) needed a damping coefficient (C x and C φ ). Therefore, the correction factor formula was converted to Equation (3) and was used in this study, wherein C x(emb) and C ∅(emb) are expressed as α x √ η x c x and α φ √ η φ c φ , respectively.

Determination of Soil Damping Based on Radiation and Material Damping Considerations
The total damping coefficient formula consists of the material damping and the radiation damping coefficients [1,3,19]. Assuming that the structure is rigid (k s = ∞) and that the foundation cannot rock (k r = ∞), or is only allowed to rock (k h = ∞), the natural frequency of each case follows 2 h = k h /m t and 2 r = k r /I 0 [1], where m t and I 0 are the total mass and total mass moment of inertia of the structure, respectively. Equation (4) expresses a form of the total damping coefficient based on the consideration of the embedment effect. Substituting the natural frequency equation in each mode in Equation (4) estimates the swaying and rocking total damping coefficients, as listed in Table 3, wherein C gx and C gφ are the respective swaying and rocking material damping coefficients. The total damping coefficient was expressed as function of C h and C r to distinguish them from the subscripts x and φ of the radiation damping coefficient. Accordingly, the soil-spring constant was also denoted by k h and k r for the swaying and rocking modes.

Effective Height of Structure h
In the soil-structure system, the effective height (h) of the structure affects its rocking response of the system. Stewart et al. [17] defined h as the distance from the foundation to the centroid of the inertial force in relation to the fundamental mode. In FEMA440 [14] and FEMA450 [18], the full height was considered as the value of h of the one-story structure, and the distance from the foundation to the center point of the first modal shape was set to h in multistory structures. In the Electric Power Research Institute (EPRI) training module [19], h was calculated based on the moment equilibrium to the center of mass of the rigid foundation. According to the recent research by Gavras et al. [10], the value of h of the footing-flexible column-bridge deck system was defined as the distance from the base of the footing to the center of the deck. However, a limited number of studies have verified or validated the optimum effective height of the structure. In this study, three effective heights (h) were considered to choose the optimum h, whereby the two heights were suggested by conventional standards, and the other height satisfied the total mass moment of inertia (I 0 ) of the structure. The three effective heights used in this study adhered to the following order: Height from the bottom of the foundation to the center of mass of the superstructure (h base to m s ) < height compatible to the total mass moment of inertia of the structure (h MMI ) < height needed to satisfy moment equilibrium (h moment ). Additionally, h MMI and h moment can be obtained as indicated below, Sustainability 2020, 12, 7113 6 of 18 where I 0 is the total mass moment of inertia of the upper, middle, and lower structures, and I f is the mass moment of inertia based on the consideration of the effective mass m f and the geometry of foundation, M 0 is the summation of the moment of the upper, middle, and lower structures, and M f is the moment that considers the effective mass m f and the geometry of the foundation.
2.3. Procedure Used to Identify the Optimum Z p and h of the Structure Figure 2 summarizes the procedure used to identify the optimum parameter using an analytical approach and test. The objective of the procedure was the identification of the optimum Z p and h. As an analytical approach, the time domain SSI analysis by the state space equation (SSE) was adopted and dynamic centrifuge tests were performed at the 20 g level.
Sustainability 2020, 12, x FOR PEER REVIEW 6 of 18 the foundation to the center of mass of the superstructure (ℎ ) < height compatible to the total mass moment of inertia of the structure (ℎ ) < height needed to satisfy moment equilibrium (ℎ ). Additionally, ℎ and ℎ can be obtained as indicated below, where is the total mass moment of inertia of the upper, middle, and lower structures, and is the mass moment of inertia based on the consideration of the effective mass and the geometry of foundation, is the summation of the moment of the upper, middle, and lower structures, and is the moment that considers the effective mass and the geometry of the foundation. Figure 2 summarizes the procedure used to identify the optimum parameter using an analytical approach and test. The objective of the procedure was the identification of the optimum and ℎ. As an analytical approach, the time domain SSI analysis by the state space equation (SSE) was adopted and dynamic centrifuge tests were performed at the 20 g level.

Dynamic Centrifuge Tests
The dynamic centrifuge model tests were performed to obtain the seismic responses of the soilstructure system and to evaluate the optimum SSI parameters. As depicted in Figure 3b, the experimental model consists of a soft sandy soil, a shallow foundation, and a Single Degree of Freedom (SDOF) structure. The beam-type centrifuge facility in the Korea Advanced Institute of Science and Technology (KAIST) and KAIST Analysis center for Research Advancement (KARA) was used in this study, and the radius and maximum payload capacity of the facility were 5 m and 240 g·tons, respectively [8]. This facility was equipped with a shaking table that operated in flight conditions that could generate maximum horizontal acceleration of 40 g in model scale. The equivalent shear beam (ESB) soil container that consists of 10 aluminum alloy rectangular frames was mounted on the shaking table. The dynamic performance of the ESB box corresponded to that of ground motion and it was validated by Lee et al. [23]. The internal dimensions of the ESB box were 0.49 m × 0.49 m × 0.6 m at the model scale, and correspond to the dimensions of 9.8 m × 9.8 m × 12.0 m at the 20 g prototype scale. Figure 3b shows the ESB box soil-structure system mounted on the shaking table. All the dimensions were expressed in the form of a prototype scale according to the scaling law of the centrifuge model test [24]. In this study, eight cases (= four structures × two input motions/structure) of dynamic centrifuge tests were conducted.

Dynamic Centrifuge Tests
The dynamic centrifuge model tests were performed to obtain the seismic responses of the soil-structure system and to evaluate the optimum SSI parameters. As depicted in Figure 3b, the experimental model consists of a soft sandy soil, a shallow foundation, and a Single Degree of Freedom (SDOF) structure. The beam-type centrifuge facility in the Korea Advanced Institute of Science and Technology (KAIST) and KAIST Analysis center for Research Advancement (KARA) was used in this study, and the radius and maximum payload capacity of the facility were 5 m and 240 g·tons, respectively [8]. This facility was equipped with a shaking table that operated in flight conditions that could generate maximum horizontal acceleration of 40 g in model scale. The equivalent shear beam (ESB) soil container that consists of 10 aluminum alloy rectangular frames was mounted on the shaking table. The dynamic performance of the ESB box corresponded to that of ground motion and it was validated by Lee et al. [23]. The internal dimensions of the ESB box were 0.49 m × 0.49 m × 0.6 m at the model scale, and correspond to the dimensions of 9.8 m × 9.8 m × 12.0 m at the 20 g prototype scale. Figure 3b shows the ESB box soil-structure system mounted on the shaking table. All the dimensions were expressed in the form of a prototype scale according to the scaling law of the centrifuge model test [24]. In this study, eight cases (=four structures × two input motions/structure) of dynamic centrifuge tests were conducted.

Analytical Approach Based on State Space Equation (SSE)
A time-domain analysis by the State Space Equation (SSE) has been used for the estimation of the dynamic structural response or system identification [11,25,26], and the analytical approach in this study was implemented using MATLAB R2014b. Equation (1) can be expressed as Equation (7) and the equation was converted to continuous-time SSE (Equation (8)

Analytical Approach Based on State Space Equation (SSE)
A time-domain analysis by the State Space Equation (SSE) has been used for the estimation of the dynamic structural response or system identification [11,25,26], and the analytical approach in this study was implemented using MATLAB R2014b. Equation (1) can be expressed as Equation (7) and the equation was converted to continuous-time SSE (Equation (8)), wherein Y in {X} consists of u s , u x , and φ. The discrete-time SSE of Equation (9) can be constructed by considering the general solution of Equation (8) [27][28][29]. In the discrete SSE, a sufficiently small-time interval (∆t) setting was required because the constant U in ∆t must be satisfied, and the stability of the system also needs to be secured [29]. Therefore, in this study, we used ∆t = 5 ms. [ Based on Equation (9) T was obtained in a step-by-step manner, and the relative acceleration response was then calculated according to Equation (10). Finally, the absolute acceleration response of the superstructure and the foundation ..

Properties of Structure in Tests and Analytical Approach
Four SDOF structural models made of steel were used. Each structural model consisted of (1) a lumped mass at the top position (i.e., upper structure, m u ), (2) two thin plates (i.e., middle plates, m m ) below the lumped mass considered to show shear deformation, and (3) a direct mat foundation (i.e., lower structure, m f ) at the bottom. To obtain the characteristics of the structure for the time domain SSI analyses, impact hammer tests were performed on four structures with the exception of all the lower structures subject to the fixed-base condition. Accordingly, acceleration time history (ATH) data were acquired. The natural frequency ( f ni ) of each system was then obtained based on the fast Fourier transformation (FFT), and the relationship between f ni and added mass (m add(i−1) ) of each model structure was formulated. Finally, the effective lateral stiffness (k s ) and effective mass of the superstructure (m s ) were estimated based on the nonlinear curve fitting of the data with the use of the natural frequency formula of SDOF shown in Figure 3a. As the massless beam was considered in the 3DOF time domain SSI analysis, m f was calculated by subtracting the m s from the total mass of the structure (m t ) that was defined as m t = m u + m m + m l = m s + m f . Given that the damping ratio of the structure cannot be obtained analytically, it is determined through free vibration tests [30]. The scaling factor N (prototype/model) is 20 because the centrifugal acceleration in each test is 20 g. The obtained  Table 4. The width of the square footing is 2 m, and the effective radius r (= r x r φ ) was calculated as 1.135 m. use of the natural frequency formula of SDOF shown in Figure 3a. As the massless beam was considered in the 3DOF time domain SSI analysis, was calculated by subtracting the from the total mass of the structure ( ) that was defined as = + + = + . Given that the damping ratio of the structure cannot be obtained analytically, it is determined through free vibration tests [30]. The scaling factor N (prototype/model) is 20 because the centrifugal acceleration in each test is 20 g. The obtained results are listed in Table 4. The width of the square footing is 2 m, and the effective radius (= ≒ ) was calculated as 1.135 m. Note: As the mass of the superstructure increases, the natural frequency clearly decreases, but the tendency of damping ratio increase is unclear.

Soil Properties and Ground Condition in Tests and Analytical Approaches
Dry silica sand layers were prepared with a relative density ( ) of 55% by using sand rainer. The properties of silica sand are listed in Table 5, and the profiles are shown in Figure 4a. In this case, USCS is the unified soil classification, is the median particle size, is the specific gravity, is the maximum void ratio, is the minimum void ratio, is the dry unit weight of the soil, is the internal friction angle of the soil, and is the calculated coefficient of the earth pressure at rest based on Jacky's formula (i.e., = 1 − ). Kim et al. [31] obtained the − ′ profile by performing in-flight bender element tests on the silica sand with the relative densities of 40% and 80%. In this study, the − profile was obtained by interpolating the results of Kim et al. [31]. Finally, the − profile was acquired based on the = / relationship. Resonant column tests were performed with mean effective stress values ′ = 25, 50, and 100 kPa, and the results were used in one-dimensional equivalent SRAs, as shown in Figure 4b, wherein / − is the modulus reduction curve of the soil, − is the damping ratio curve of the soil and ′ is the mean effective stress that can be calculated as (1 + 2 ) /3.
use of the natural frequency formula of SDOF shown in Figure 3a. As the massless beam was considered in the 3DOF time domain SSI analysis, was calculated by subtracting the from the total mass of the structure ( ) that was defined as = + + = + . Given that the damping ratio of the structure cannot be obtained analytically, it is determined through free vibration tests [30]. The scaling factor N (prototype/model) is 20 because the centrifugal acceleration in each test is 20 g. The obtained results are listed in Table 4. The width of the square footing is 2 m, and the effective radius (= ≒ ) was calculated as 1.135 m. Note: As the mass of the superstructure increases, the natural frequency clearly decreases, but the tendency of damping ratio increase is unclear.

Soil Properties and Ground Condition in Tests and Analytical Approaches
Dry silica sand layers were prepared with a relative density ( ) of 55% by using sand rainer. The properties of silica sand are listed in Table 5, and the profiles are shown in Figure 4a. In this case, USCS is the unified soil classification, is the median particle size, is the specific gravity, is the maximum void ratio, is the minimum void ratio, is the dry unit weight of the soil, is the internal friction angle of the soil, and is the calculated coefficient of the earth pressure at rest based on Jacky's formula (i.e., = 1 − ). Kim et al. [31] obtained the − ′ profile by performing in-flight bender element tests on the silica sand with the relative densities of 40% and 80%. In this study, the − profile was obtained by interpolating the results of Kim et al. [31]. Finally, the − profile was acquired based on the = / relationship. Resonant column tests were performed with mean effective stress values ′ = 25, 50, and 100 kPa, and the results were used in one-dimensional equivalent SRAs, as shown in Figure 4b, wherein / − is the modulus reduction curve of the soil, − is the damping ratio curve of the soil and ′ is the mean effective stress that can be calculated as (1 + 2 ) /3.
use of the natural frequency formula of SDOF shown in Figure 3a. As the massless beam was considered in the 3DOF time domain SSI analysis, was calculated by subtracting the from the total mass of the structure ( ) that was defined as = + + = + . Given that the damping ratio of the structure cannot be obtained analytically, it is determined through free vibration tests [30]. The scaling factor N (prototype/model) is 20 because the centrifugal acceleration in each test is 20 g. The obtained results are listed in Table 4. The width of the square footing is 2 m, and the effective radius (= ≒ ) was calculated as 1.135 m. Note: As the mass of the superstructure increases, the natural frequency clearly decreases, but the tendency of damping ratio increase is unclear.

Soil Properties and Ground Condition in Tests and Analytical Approaches
Dry silica sand layers were prepared with a relative density ( ) of 55% by using sand rainer. The properties of silica sand are listed in Table 5, and the profiles are shown in Figure 4a. In this case, USCS is the unified soil classification, is the median particle size, is the specific gravity, is the maximum void ratio, is the minimum void ratio, is the dry unit weight of the soil, is the internal friction angle of the soil, and is the calculated coefficient of the earth pressure at rest based on Jacky's formula (i.e., = 1 − ). Kim et al. [31] obtained the − ′ profile by performing in-flight bender element tests on the silica sand with the relative densities of 40% and 80%. In this study, the − profile was obtained by interpolating the results of Kim et al. [31]. Finally, the − profile was acquired based on the = / relationship. Resonant column tests were performed with mean effective stress values ′ = 25, 50, and 100 kPa, and the results were used in one-dimensional equivalent SRAs, as shown in Figure 4b, wherein / − is the modulus reduction curve of the soil, − is the damping ratio curve of the soil and ′ is the mean effective stress that can be calculated as (1 + 2 ) /3. use of the natural frequency formula of SDOF shown in Figure 3a. As the massless beam was considered in the 3DOF time domain SSI analysis, was calculated by subtracting the from the total mass of the structure ( ) that was defined as = + + = + . Given that the damping ratio of the structure cannot be obtained analytically, it is determined through free vibration tests [30]. The scaling factor N (prototype/model) is 20 because the centrifugal acceleration in each test is 20 g. The obtained results are listed in Table 4. The width of the square footing is 2 m, and the effective radius (= ≒ ) was calculated as 1.135 m. Note: As the mass of the superstructure increases, the natural frequency clearly decreases, but the tendency of damping ratio increase is unclear.

Soil Properties and Ground Condition in Tests and Analytical Approaches
Dry silica sand layers were prepared with a relative density ( ) of 55% by using sand rainer. The properties of silica sand are listed in Table 5, and the profiles are shown in Figure 4a. In this case, USCS is the unified soil classification, is the median particle size, is the specific gravity, is the maximum void ratio, is the minimum void ratio, is the dry unit weight of the soil, is the internal friction angle of the soil, and is the calculated coefficient of the earth pressure at rest based on Jacky's formula (i.e., = 1 − ). Kim et al. [31] obtained the − ′ profile by performing in-flight bender element tests on the silica sand with the relative densities of 40% and 80%. In this study, the − profile was obtained by interpolating the results of Kim et al. [31]. Finally, the − profile was acquired based on the = / relationship. Resonant column tests were performed with mean effective stress values ′ = 25, 50, and 100 kPa, and the results were used in one-dimensional equivalent SRAs, as shown in Figure 4b, wherein / − is the modulus reduction curve of the soil, − is the damping ratio curve of the soil and ′ is the mean effective stress that can be calculated as (1 + 2 ) /3. Note: As the mass of the superstructure increases, the natural frequency clearly decreases, but the tendency of damping ratio increase is unclear.

Soil Properties and Ground Condition in Tests and Analytical Approaches
Dry silica sand layers were prepared with a relative density (D r ) of 55% by using sand rainer. The properties of silica sand are listed in Table 5, and the V s profiles are shown in Figure 4a. In this case, USCS is the unified soil classification, D 50 is the median particle size, G s is the specific gravity, e max is the maximum void ratio, e min is the minimum void ratio, γ d is the dry unit weight of the soil, φ is the internal friction angle of the soil, and K 0 is the calculated coefficient of the earth pressure at rest based on Jacky's formula (i.e., K 0 = 1 − sinφ). Kim et al. [31] obtained the V s − σ v profile by performing in-flight bender element tests on the silica sand with the relative densities of 40% and 80%. In this study, the V s − σ v profile was obtained by interpolating the results of Kim et al. [31]. Finally, the V s − Z profile was acquired based on the Z = σ v /γ d relationship. Resonant column tests were performed with mean effective stress values σ m = 25, 50, and 100 kPa, and the results were used in one-dimensional equivalent SRAs, as shown in Figure 4b, wherein G/G max − logγ is the modulus reduction curve of the soil, ξ g − logγ is the damping ratio curve of the soil and σ m is the mean effective stress that can be calculated as (1 + 2K 0 )σ v /3.

Input Motions in Tests and analyses
In this study, three input motions were used for tests and analyses: (1) Hachinohe earthquake motion with abundant low-frequency contents, (2) Northridge earthquake motion with abundant high-frequency contents, and (3) synthesized motion that consisted of 2, 6, and 10 Hz sinusoidal waves. Figure 5 is normalized Fourier amplitude spectra showing the frequency characteristics of each representative input motion. The characteristics of each motion type, including peak ground acceleration (PGA), are summarized in Table 6, wherein contains every frequency of the Fourier amplitude spectrum that ranges from 0.25 to 20 Hz (Δ ≤ 0.05 Hz), is the Fourier amplitude, and the parameter is the reciprocal of the mean period ( ) suggested by Rathje et al. [33] and used as a representative frequency motion index.
Base motions were the measured outcrop motions at the outer bottom of the centrifuge test box and were used to SRAs. Free-field motions (FFMs) were measured as within motions at half the embedment depth of the foundation (i.e., Z = 0.3 m) and were used for analytical approach (i.e., 3DOF SSI analyses).

Input Motions in Tests and Analyses
In this study, three input motions were used for tests and analyses: (1) Hachinohe earthquake motion with abundant low-frequency contents, (2) Northridge earthquake motion with abundant high-frequency contents, and (3) synthesized motion that consisted of 2, 6, and 10 Hz sinusoidal waves. Figure 5 is normalized Fourier amplitude spectra showing the frequency characteristics of each representative input motion. The characteristics of each motion type, including peak ground acceleration (PGA), are summarized in Table 6, wherein f i contains every frequency of the Fourier amplitude spectrum that ranges from 0.25 to 20 Hz (∆ f ≤ 0.05 Hz), A i is the Fourier amplitude, and the parameter f m is the reciprocal of the mean period (T m ) suggested by Rathje et al. [33] and used as a representative frequency motion index.
Base motions were the measured outcrop motions at the outer bottom of the centrifuge test box and were used to SRAs. Free-field motions (FFMs) were measured as within motions at half the embedment depth of the foundation (i.e., Z = 0.3 m) and were used for analytical approach (i.e., 3DOF SSI analyses).

Site-Response Analyses (SRAs) and Strain-Independent Soil Properties
One dimensional equivalent linear SRAs, One of the methods that can be used to obtain the

Site-Response Analyses (SRAs) and Strain-Independent Soil Properties
One dimensional equivalent linear SRAs, One of the methods that can be used to obtain the acceleration response of the soil layer and to obtain the equivalent linear soil properties, were conducted to obtain strain-independent V s and ξ g values with the program STRATA that was developed by Kottke and Rathje [34]. Figure 6a shows that the ATHs from SRAs match well those from the test results at Z = 0.3 m, and Figure 6b,c shows the equivalent V s and ξ g -depth results (i.e., the strain-independent V s depth and ξ g -depth relationships) at their final effective shear strain values.  Figure 7 shows the variations of and that are converted to their depth-independent forms using Equation (2) Figure 7 shows the variations of V s and ξ g that are converted to their depth-independent forms using Equation (2), and are based on three types of Z p scenarios. Regarding the effective radius r in square footing, r x = 0.56 B in the swaying mode is similar to r φ = 0.57 B in the rocking mode, and the effective profile depth Z p is normalized with one effective radius r (i.e., r = 0.565 B). In addition, Zp/r = 0.75, 2, and 4, correspond to Z p = 0.85, 2.3, and 4.5 m, respectively.  Figure 7 shows the variations of and that are converted to their depth-independent forms using Equation (2), and are based on three types of scenarios. Regarding the effective radius in square footing, = 0.56 B in the swaying mode is similar to = 0.57 B in the rocking mode, and the effective profile depth is normalized with one effective radius (i.e., = 0.565 B). In addition, / = 0.75, 2, and 4, correspond to = 0.85, 2.3, and 4.5 m, respectively.

Determination of Analysis CASES Considering the SSI Parameters
The soil-spring constant and damping coefficient depend on V s and ξ g as well as on the foundation conditions (e.g., geometry of foundation, embedded depth, and other subsoil properties) as follows • k h(emb) and k r(emb) = f (ρ, V s , geometry o f f oundation, ν, D) • C h(emb) and C r(emb) = f ρ, V s , geometry o f f oundation, ν, D, ξ g , m t or I 0 V s and ξ g are determined by the effective profile depth and I 0 determined by the effective height affects C g∅ . Given that h and Z p affect the soil-spring constants and damping coefficients, the number of the analysis cases were set as follows: • In the case of the optimum h value: 24 soil-spring constants and damping coefficients = one formula × eight input motions × three effective heights × one effective profile depth • In the case of the optimum Z p value: 48 soil-spring constants and damping coefficients = two formulas × eight input motions × one effective height × three effective profile depths

Quantification of Differences between Test and Analysis Results
To evaluate the differences in phase, amplitude, frequency contents, and amplification characteristics between the test and analysis results, the mean squared errors (MSEs) for ATH and response spectrum (RS) of each superstructure were obtained in the forms of MSE TH and MSE RS , respectively, and their sums (i.e., MSE i ) were used to comprehensively evaluate the conformity of the analysis to the test results. In Equations (13) to (15), a test(k) and a anal(k) were the acceleration Sustainability 2020, 12, 7113 12 of 18 responses of the test and analysis at a certain time t = k, and SA test( j) and SA anal( j) were the spectral acceleration values of the test and analyses at a certain frequency f = j, respectively. When the value of MSE i is small, the overall agreement between the test and analysis results is better.

Evaluation of the Optimum h
To evaluate the optimum effective height (h), 24 analyses were conducted in total (i.e., 24 cases = four structures × two input motions/structure × three effective heights), where Z p /r = 4 (i.e., Z p = 4.5 m). The Richart and Lysmer formula was considered for the analyses. In each case, the minimum h value corresponds to the base to m s case, and the maximum h value to the moment equilibrium case. As h changes, I 0 and m s h also changes in all the analysis cases. The analysis cases are summarized in Table 7, wherein no case matched the I 0 and m s h condition simultaneously. The shaded area denoted the case matching of each property of the prototype structure (i.e., I 0 or m s h).  Figure 8 shows the ATHs and ∑MSE values of the test and analysis results. According to the results of the analysis, as h increases, ∑MSE increases. Based on the tendency and calculation result of ∑MSE , ∑MSE at ℎ yielded a minimum error of 0.106 on average, ∑MSE at h yielded an error of 0.127 on average, and ∑MSE at ℎ yielded the greatest error of 0.229 on average. Given that ∑MSE at ℎ was estimated to be lower than that at full height, ℎ was considered as the optimum effective height.

Evaluation of the Optimum Effective Profile Depth
In the evaluation of the optimum effective profile depth ( ), 48 analyses were conducted in total (i.e., 48 cases = four structures × two input motions/structure × three effective profile depths × two

Evaluation of the Optimum Effective Profile Depth Z p
In the evaluation of the optimum effective profile depth (Z p ), 48 analyses were conducted in total (i.e., 48 cases = four structures × two input motions/structure × three effective profile depths × two types of formulas), wherein the embedment correction factor suggested by Whitman [6] and h base to m s was used in common. The analysis cases are summarized in Table 8. Figure 9 shows the MSE i value of each analysis case, and Table 9 lists the maximum acceleration responses of the superstructure, wherein the shaded areas represent the analysis cases that best match the test results (i.e., the case with the minimum MSE i ). Although two analysis cases (i.e., structure 3-Northridge II and structure 4-Synthesized II) show that the maximum acceleration responses were smaller than those of the corresponding test results, the differences were not considerable.
Specifically, the maximum acceleration responses of the conducted analyses were 93.2% and 99.0% compared with those of the corresponding test results. The average values of MSE i in the case at which Wolf's formula was used was 0.098, and the average value of MSE i in the case at which Richart and Lysmer's formula was used was 0.106. Thus, the difference of MSE i between the results of the two formulas was not significant. In terms of the maximum acceleration response, the results generated by Richart and Lysmer's formula were found to be 0.8-12.4% more conservative than those obtained based on Wolf's formula, as listed in Table 9. In terms of the input motion, the average value of MSE i was 0.137 for the Hachinohe motion, 0.124 for Northridge motion, and 0.045 for synthesized motion. These findings show that the analysis results based on the use of the seismic wave were associated with a larger error compared with that obtained based on the analysis with the use of synthesized motion. In all the tested cases, the acceleration response of the superstructure at Z p /r = 0.75 yielded a higher discrepancy between tests and analysis outcomes (i.e., the maximum MSE i value in all the cases of Z p /r = 0.75 was 2.403) but the responses at Z p /r = 4 showed good agreement (i.e., maximum MSE i in all the cases of Z p /r = 4 was 0.425). types of formulas), wherein the embedment correction factor suggested by Whitman [6] and ℎ was used in common. The analysis cases are summarized in Table 8.  Figure 9 shows the ∑MSE value of each analysis case, and Table 9 lists the maximum acceleration responses of the superstructure, wherein the shaded areas represent the analysis cases that best match the test results (i.e., the case with the minimum ∑MSE ). Although two analysis cases (i.e., structure 3-Northridge II and structure 4-Synthesized II) show that the maximum acceleration responses were smaller than those of the corresponding test results, the differences were not considerable. Specifically, the maximum acceleration responses of the conducted analyses were 93.2% and 99.0% compared with those of the corresponding test results. The average values of ∑MSE in the case at which Wolf's formula was used was 0.098, and the average value of ∑MSE in the case at which Richart and Lysmer's formula was used was 0.106. Thus, the difference of ∑MSE between the results of the two formulas was not significant. In terms of the maximum acceleration response, the results generated by Richart and Lysmer's formula were found to be 0.8-12.4% more conservative than those obtained based on Wolf's formula, as listed in Table 9. In terms of the input motion, the average value of ∑MSE was 0.137 for the Hachinohe motion, 0.124 for Northridge motion, and 0.045 for synthesized motion. These findings show that the analysis results based on the use of the seismic wave were associated with a larger error compared with that obtained based on the analysis with the use of synthesized motion. In all the tested cases, the acceleration response of the superstructure at /r = 0.75 yielded a higher discrepancy between tests and analysis outcomes (i.e., the maximum ∑MSE value in all the cases of /r = 0.75 was 2.403) but the responses at /r = 4 showed good agreement (i.e., maximum ∑MSE in all the cases of /r = 4 was 0.425).    Figure 10 shows the variations of the ATHs and RS of the superstructure during Hachinohe motion. The tendency to increase or decrease the maximum acceleration of the superstructure as the effective profile depth increases is unclear, but as the effective profile depth increases, MSE i clearly decreases, thus yielding the lowest value at Z p /r = 4. Therefore, based on all the results in this study, it can be concluded that 4 r is the optimum effective profile depth.   Figure 10 shows the variations of the ATHs and RS of the superstructure during Hachinohe motion. The tendency to increase or decrease the maximum acceleration of the superstructure as the effective profile depth increases is unclear, but as the effective profile depth increases, ∑MSE clearly decreases, thus yielding the lowest value at /r = 4. Therefore, based on all the results in this study, it can be concluded that 4 r is the optimum effective profile depth.

Conclusions
In this study, dynamic centrifuge tests and LDP-based analyses (i.e., time domain SSI analyses by SSE) were performed to evaluate the optimum effective height and effective profile depth conditions proposed in the conventional standards and prior research publications. Four structures and three ground motions were used in the centrifuge tests, and three effective heights and three effective profile depth conditions were considered as the SSI analysis cases in addition to the aforementioned test conditions. The main results of this research are summarized as follows.

1.
In this study, the applicability of the SSI parameters suggested by various standards and studies was discussed, and the optimum SSI parameter selection procedure that (a) comprehensively considered the SSI parameters, (b) adopted an analytical approach and a physical model test, was suggested. Based on the established procedure, the optimum values of two controversial SSI parameters (i.e., the effective height and effective profile depth) were determined • Unlike the conventional standards that apply a simplified reduction factor for the initial shear wave velocity profile and do not apply any corrections in initial damping ratio profile, one dimensional equivalent linear site response analyses were performed to accurately obtain the equivalent strain-independent shear wave velocity (V s ) and damping ratio (ξ g ) of the soil. The equivalent V s and ξ g values that varied with depth obtained herein were converted to depth independent V s and ξ g values based on considerations of the effective profile depth (Z p ).

•
Unlike the conventional research efforts that ignored soil material damping and indirectly determined soil damping based on the effective period lengthening ratio, the total soil damping was obtained directly by the addition of soil material damping to soil radiation damping. In addition, the stiffness and total damping of soil were determined based on embedded foundation conditions. • Unlike the conventional SDOF SSI analysis that was based on the RS, this study adopted a 3DOF time domain SSI analysis based on structural translation, foundation translation, and rocking behavior considerations to accurately obtain structural responses.

2.
The effective height of the structure affected the rocking behavior of the soil-structure system (i.e., m s h and C g∅ ). In this study, applicability of the following three effective height scenarios were evaluated based on the following test and analysis results: (1) Height from the bottom of the foundation to the center of the mass of the superstructure, and (2) height compatible to the total mass moment of inertia, and (3) height to satisfy moment equilibrium. As the effective height increased within the effective height range used in the analysis, the differences between test and analysis results increased. Consequently, in all the cases, the height from the bottom of the foundation to the center of the mass of the superstructure was the optimum effective height with the lowest MSE i value.

3.
The optimum effective profile depth (Z p ) was determined. This was used to average the dynamic soil property that varied as a function of depth. Three scenarios were considered in the analyses: (1) Z p /r = 0.75, (2) Z p /r = 2, and (3) Z p /r = 4. As a result, 4 r was found to be the optimum effective profile depth with the lowest difference between the test and analysis results. In addition, the results associated with earthquake motions yielded higher errors than those of synthesized motions that consisted of sinusoidal waves at the frequencies of 2, 6, and 10 Hz. The maximum acceleration responses by Richart and Lysmer's formula was more conservative than those obtained with Wolf's formula.