# Application of the Value Optimization Model of Key Factors Based on DSEM

Authors: Chao Su; Zhiming Ren
Journal: Mathematical Problems in Engineering
 (2016)

---

## Abstract

The key factors of the damping solvent extraction method (DSEM) for the analysis of the unbounded medium are the size of bounded domain, the artificial damping ratio, and the finite element mesh density. To control the simulation accuracy and computational efficiency of the soil-structure interaction, this study establishes a value optimization model of key factors that is composed of the design variables, the objective function, and the constraint function system. Then the optimum solutions of key factors are obtained by the optimization model. According to some comparisons of the results provided by the different initial conditions, the value optimization model of key factors is feasible to govern the simulation accuracy and computational efficiency and to analyze the practical unbounded medium-structure interaction.


---

## Body

## 1. Introduction

The rationality of numerical simulation in seismic loads plays an important role in the dynamic response analysis of the soil-structure interaction. With the development of science and technology, the research and application into the soil-structure interaction are of interest to those in the engineering profession and academia. Moreover, many researches on the unbounded media-structure interaction improve the rationality and accuracy of the foundations in frequency and time domains.The methods for analyzing the soil-structure interaction are currently the artificial boundary method, the boundary element method, and the damping solvent extraction method. The artificial boundary method [1, 2] intercepts the bounded domain and increases the energy transfer system at the outer boundary (such as by a damper or spring-damper). But this method is not appropriate for analyzing an irregular foundation. The boundary element method [3] proposes the principle of wave amplitude attenuation of various frequencies far afield and makes a transformation into the time domain. Likewise, this method has a low computational efficiency and does not analyze a complicated unbounded media. However, through introducing and extracting the artificial damping ratio of the bounded medium, the damping solvent extraction method [4, 5] is appropriate for the simulation of the various foundations. Meanwhile, the DSEM implements the mutual conversion in the frequency-time domain and eliminates the convolution operation in the time domain.Since the DSEM was first presented by Wolf and Song in 1994, some researches have presented gradually the perfect basis theory, the improved computational procedure, and the evaluation of key factors of the soil-structure interaction. In the computational algorithm, a series of calculation procedures had been suggested to improve the accuracy, such as the acceleration input stagger method [6], the precise step-by-step time integration scheme [7], and the subregional explicit implicit recursive method [8, 9]. In its practical application, the DSEM was applied to analyze the soil-structure interaction with fractional order [10], the nonlinear overload of the unbounded medium-arch dam [11], the impact of pile driving vibration on seawall [12], and fluid-structure dynamic interaction [13]. In the evaluation of key factors, the causes of the error of the DSEM had been analyzed and described in detail in [14–16]. Simultaneously, the key factors estimated in [17–19] have a great influence on the accuracy of the soil-structure interaction.The size of bounded domain and the finite element mesh density were randomly assigned and the artificial damping ratio was limited in a specified scope, so the key influence factors of the DSEM were roughly analyzed both in the analysis and in the application. To meet the requirements of practical engineering and get the optimum solution of key factors, this study analyzes the dimensionless dynamic stiffness coefficients of various mediums and proposes the value optimization model of key factors. The initial conditions of the optimization model are the maximum of the difference function of the dynamic stiffness coefficient, the wave amplitude attenuation, and the maximum finite element number. Numerical verification demonstrates that the value optimization model of key factors can regulate the simulation accuracy and computational efficiency and analyze the practical infinite foundation-structure interaction.


## 2. Damping Solvent Extraction Method

The damping solvent extraction method is simply reviewed in this chapter. A bounded domain of the infinite media adjacent to the soil-structure interface is discretized with a finite element meshes shown in Figure1.Figure 1
Finite element mesh of bounded domain adjacent to structure.


### 2.1. Implementation in the Frequency Domain

The procedure of the DSEM consists of the following steps for each frequency, which calculate the dynamic stiffness coefficients of the unbounded medium.Step 1.
In the selected bounded domain idealized as an assemblage of the finite elements, the artificial damping ratio in the domain and the viscous dashpots at the outer boundary serve to attenuate the vibration amplitudes of the outgoing wave and reflected wave. The coefficients of the viscous dashpots per unit surface area at the outer boundary areρcs∗ in the perpendicular direction and ρcp∗ in the tangential direction. The shear-wave velocity and the dilatation-wave velocity are affected by the factor 1+2iζ. Consider(1a)cs∗=cs1+2iζ,cp∗=cp1+2iζ,(1b)G∗=G1+2iζ,a0∗=a01+2iζ,where G and a0 are the shear modulus and dimensionless frequency, respectively.Step 2.
Eliminating all degrees of freedom, which are located in the medium except the interface, the dynamic stiffness condensed matrix of the bounded medium,S~f(ω), is transformed from the dynamic stiffness matrix Sf(ω). The dynamic stiffness matrix of the damped bounded foundation is expresses as follows:(2)Sfω=G∗Sfa0∗=G∗-a0∗2M-f+ia0∗C-f+K-f,where M-f, C-f, and K-f are the dimensionless stiffness, damping, and mass matrices.Step 3.
Since the dynamic stiffness coefficient is derived withω and the dynamic stiffness condensed matrix of the bounded medium, S~f(ω), convert to the dynamic stiffness condensed matrix of the unbounded media, S~∞(ω), the dimensionless dynamic stiffness condensed matrix, S~∞(a0∗), is written as follows:(3)S~∞a0∗=S~fa0∗,S~∞a0∗,a0∗=S~fa0∗,a0∗.
Taking the first-order Taylor expansion ofS~∞(a0∗), the dynamic stiffness matrix of the undamped unbounded medium is equal to(4)S~∞a0=S~∞a0∗+S~∞a0∗,a0∗a0-a0∗.


### 2.2. Implementation in the Time Domain

The medium-structure interaction force is obtained by the DSEM in the time domain, so the motion equation of the damped bounded medium takes the form(5)Mmm00Mmmu¨mu¨b+CmmCmbCbmCbbfu˙mu˙b+MmmMmbMbmMbbumub=0Rbtwith(6)K=K-+ζ2M-,C=C-+2ζM-,M=M-,where M-, C-, and K- are the mass, damping, and stiffness matrices, respectively, and Rb is the interaction force at the soil-structure interface. The subscript b denotes the nodes of the bounded medium at the soil-structure interface and m denotes the remaining nodes. Furthermore, the interaction force of the soil-structure was written in [19]:(7)Rbt=Mbbu¨b+Kbb-ζCbbub+Kbmum+ζKbbvm,where u¨b, u˙b, and ub are the displacement, velocity, and acceleration vectors and the unknown um and vm vectors are provided by the following equations:(8)Mmmu¨m+Cmmu˙m+Kmmum=-Kmbub-Cmbu˙b,Mmmv¨m+Cmmv˙m+Kmmvm=2Mmmu˙m+Cmmum+Cmbub.


## 2.1. Implementation in the Frequency Domain

The procedure of the DSEM consists of the following steps for each frequency, which calculate the dynamic stiffness coefficients of the unbounded medium.Step 1.
In the selected bounded domain idealized as an assemblage of the finite elements, the artificial damping ratio in the domain and the viscous dashpots at the outer boundary serve to attenuate the vibration amplitudes of the outgoing wave and reflected wave. The coefficients of the viscous dashpots per unit surface area at the outer boundary areρcs∗ in the perpendicular direction and ρcp∗ in the tangential direction. The shear-wave velocity and the dilatation-wave velocity are affected by the factor 1+2iζ. Consider(1a)cs∗=cs1+2iζ,cp∗=cp1+2iζ,(1b)G∗=G1+2iζ,a0∗=a01+2iζ,where G and a0 are the shear modulus and dimensionless frequency, respectively.Step 2.
Eliminating all degrees of freedom, which are located in the medium except the interface, the dynamic stiffness condensed matrix of the bounded medium,S~f(ω), is transformed from the dynamic stiffness matrix Sf(ω). The dynamic stiffness matrix of the damped bounded foundation is expresses as follows:(2)Sfω=G∗Sfa0∗=G∗-a0∗2M-f+ia0∗C-f+K-f,where M-f, C-f, and K-f are the dimensionless stiffness, damping, and mass matrices.Step 3.
Since the dynamic stiffness coefficient is derived withω and the dynamic stiffness condensed matrix of the bounded medium, S~f(ω), convert to the dynamic stiffness condensed matrix of the unbounded media, S~∞(ω), the dimensionless dynamic stiffness condensed matrix, S~∞(a0∗), is written as follows:(3)S~∞a0∗=S~fa0∗,S~∞a0∗,a0∗=S~fa0∗,a0∗.
Taking the first-order Taylor expansion ofS~∞(a0∗), the dynamic stiffness matrix of the undamped unbounded medium is equal to(4)S~∞a0=S~∞a0∗+S~∞a0∗,a0∗a0-a0∗.


## 2.2. Implementation in the Time Domain

The medium-structure interaction force is obtained by the DSEM in the time domain, so the motion equation of the damped bounded medium takes the form(5)Mmm00Mmmu¨mu¨b+CmmCmbCbmCbbfu˙mu˙b+MmmMmbMbmMbbumub=0Rbtwith(6)K=K-+ζ2M-,C=C-+2ζM-,M=M-,where M-, C-, and K- are the mass, damping, and stiffness matrices, respectively, and Rb is the interaction force at the soil-structure interface. The subscript b denotes the nodes of the bounded medium at the soil-structure interface and m denotes the remaining nodes. Furthermore, the interaction force of the soil-structure was written in [19]:(7)Rbt=Mbbu¨b+Kbb-ζCbbub+Kbmum+ζKbbvm,where u¨b, u˙b, and ub are the displacement, velocity, and acceleration vectors and the unknown um and vm vectors are provided by the following equations:(8)Mmmu¨m+Cmmu˙m+Kmmum=-Kmbub-Cmbu˙b,Mmmv¨m+Cmmv˙m+Kmmvm=2Mmmu˙m+Cmmum+Cmbub.


## 3. Constraint Function

The constraint functions regarded as the essential components are crucial if one wishes to solve the optimization model. Therefore, the references and analysis of the constraint functions are described in detail in this chapter.


### 3.1. Constraint Function of Artificial Damping Ratio

The dynamic stiffness coefficientS∞a0 and the dimensionless dynamic stiffness coefficient S¯∞a0 of the undamped unbounded medium are, respectively, given by(9)S∞a0=K∞ka0+ia0ca0,(10)S¯∞a0=S∞a0K∞=1-a02=ka0+ia0ca0.Ifa0≤1, ka0=1-a02 and ca0=0.Ifa0>1, ka0=0 and ca0=1-1/a02.On the basis of the plural damping coefficient correspondence principle, the dimensionless dynamic stiffness coefficient of the damped unbounded medium is deduced as(11)S¯∞a0∗=1-a0∗2with (12)a0∗=ω-iζbc1=a0-iζ¯,where a0∗ and a0 are the dimensionless frequencies of the damped and undamped bounded medium.The dimensionless dynamic stiffness of the damped bounded medium is determined by the size of the bounded domain,l/b, and the artificial damping ratio ζ-:(13)S¯a0∗=1-a0∗21-c1/c21+c1/c2,where c1 and c2 are the vibration amplitudes of the incoming wave and outgoing wave.The dynamic stiffness coefficient of the undamped unbounded medium is calculated by taking the first-order Taylor expansion to eliminate the artificial damping ratio:(14)S¯∞a0=S¯a0∗+S¯a0∗,a0∗a0-a0∗.The various values of the size of the bounded medium and the artificial damping ratio listed in Table1 are used to analyze the dynamic stiffness coefficients of the undamped unbounded, damped unbounded, damped bounded, and simulating undamped unbounded medium.Table 1
Values of bounded domain size and artificial damping ratio.
Case C1 C2 C3 C4 C5 C6 C7 C8 l / b 3 3 3 3 1 3 5 7 ζ ¯ 0.01 0.1 0.5 1 0.3 0.3 0.3 0.3When the dynamic stiffness coefficient of the undamped unbounded medium is regarded as the accurate solution, the dynamic stiffness coefficients determined by(l/b=3,ζ-=0.2) shown in Figure 2 indicate that (1) the larger artificial damping ratio increases the difference between the dynamic stiffness coefficient of the undamped and damped unbounded medium; (2) the larger size of the bounded medium reduces the amplitude of the dynamic stiffness coefficient of the damped unbounded medium; (3) the smaller artificial damping ratio and the larger size of the bounded medium lead to the more accurate dynamic stiffness coefficient of the undamped unbounded medium. Therefore, the constraint function of artificial damping ratio is the difference function between the dimensionless dynamic stiffness coefficient of the undamped and damped unbounded medium.Figure 2
Dimensionless dynamic stiffness coefficients ofl/b=3, ζ-=0.2.


### 3.2. Constraint Function of Size of the Bounded Medium

Suppose that the artificial damping ratioζ-=ζ(b/cs) is zero and the size of bounded medium is l/b; e-iwl/cs is the phase angle caused by radiating the shear wave from the interface to the outer boundary of the undamped bounded domain. If the artificial damping ratio is nonzero, the frequency ω is replaced by w-iζ and the phase angle is replaced by e-ζl/c1·e-iwl/c1. Thus, e-ζl/c1 is the attenuation of the wave amplitude in the damped bounded domain.Due to the addition of the artificial damping ratio into the bounded medium, the wave is reflected by the outer boundary and radiated to the interface. The wave amplitude attenuation function of the initial displacement wave excited at the interface is written as(15)ϕζ,lb=exp⁡-2ζlbbc1=exp⁡-2ζ¯lc1.If the wave amplitude attenuation at the interface is known, the constraint function of the size of bounded medium is transformed from function (15). Consider (16)lb=-ln⁡ϕ02ζ¯.


### 3.3. Constraint Function of Finite Element Mesh Density

In the finite element discretization of the bounded medium, the finite element mesh of only one row, denoted as the “base row,” is merely required.Sr denotes the dynamic stiffness matrix of one row of the finite element mesh with length r (Figure 3). The same is true for the condensed matrix obtained by eliminating all degrees of freedom except those related to interior nodes:(17)Sαbω=Sbαω,S~αbω=S~bαω.Figure 3
First row of finite element mesh.The formula relatesS~r1 for exterior nodes to the corresponding matrix of interior nodes S~r0:(18)S~r1ω=S~r0κω,in which κ=r1/r0=γ2.The dashpots are assigned at the exterior nodes of the base row. BothS~r0γ2n-2ω and S~r1γ2n-4ω are calculated by (18). After removing the dashpots, either S~r0γ2n-2ω or S~r1γ2n-4ω is assembled at the exterior nodes of the base row for frequency γ2n-4ω. S~r0γ2n-4ω is calculated by condensing again into the dynamic stiffness matrix for interior nodes. S~fω is obtained by repeating the last process n-2 and the proportionality factor γ must satisfy(19)r0+1/2Δ+γΔr0=γ2.


## 3.1. Constraint Function of Artificial Damping Ratio

The dynamic stiffness coefficientS∞a0 and the dimensionless dynamic stiffness coefficient S¯∞a0 of the undamped unbounded medium are, respectively, given by(9)S∞a0=K∞ka0+ia0ca0,(10)S¯∞a0=S∞a0K∞=1-a02=ka0+ia0ca0.Ifa0≤1, ka0=1-a02 and ca0=0.Ifa0>1, ka0=0 and ca0=1-1/a02.On the basis of the plural damping coefficient correspondence principle, the dimensionless dynamic stiffness coefficient of the damped unbounded medium is deduced as(11)S¯∞a0∗=1-a0∗2with (12)a0∗=ω-iζbc1=a0-iζ¯,where a0∗ and a0 are the dimensionless frequencies of the damped and undamped bounded medium.The dimensionless dynamic stiffness of the damped bounded medium is determined by the size of the bounded domain,l/b, and the artificial damping ratio ζ-:(13)S¯a0∗=1-a0∗21-c1/c21+c1/c2,where c1 and c2 are the vibration amplitudes of the incoming wave and outgoing wave.The dynamic stiffness coefficient of the undamped unbounded medium is calculated by taking the first-order Taylor expansion to eliminate the artificial damping ratio:(14)S¯∞a0=S¯a0∗+S¯a0∗,a0∗a0-a0∗.The various values of the size of the bounded medium and the artificial damping ratio listed in Table1 are used to analyze the dynamic stiffness coefficients of the undamped unbounded, damped unbounded, damped bounded, and simulating undamped unbounded medium.Table 1
Values of bounded domain size and artificial damping ratio.
Case C1 C2 C3 C4 C5 C6 C7 C8 l / b 3 3 3 3 1 3 5 7 ζ ¯ 0.01 0.1 0.5 1 0.3 0.3 0.3 0.3When the dynamic stiffness coefficient of the undamped unbounded medium is regarded as the accurate solution, the dynamic stiffness coefficients determined by(l/b=3,ζ-=0.2) shown in Figure 2 indicate that (1) the larger artificial damping ratio increases the difference between the dynamic stiffness coefficient of the undamped and damped unbounded medium; (2) the larger size of the bounded medium reduces the amplitude of the dynamic stiffness coefficient of the damped unbounded medium; (3) the smaller artificial damping ratio and the larger size of the bounded medium lead to the more accurate dynamic stiffness coefficient of the undamped unbounded medium. Therefore, the constraint function of artificial damping ratio is the difference function between the dimensionless dynamic stiffness coefficient of the undamped and damped unbounded medium.Figure 2
Dimensionless dynamic stiffness coefficients ofl/b=3, ζ-=0.2.


## 3.2. Constraint Function of Size of the Bounded Medium

Suppose that the artificial damping ratioζ-=ζ(b/cs) is zero and the size of bounded medium is l/b; e-iwl/cs is the phase angle caused by radiating the shear wave from the interface to the outer boundary of the undamped bounded domain. If the artificial damping ratio is nonzero, the frequency ω is replaced by w-iζ and the phase angle is replaced by e-ζl/c1·e-iwl/c1. Thus, e-ζl/c1 is the attenuation of the wave amplitude in the damped bounded domain.Due to the addition of the artificial damping ratio into the bounded medium, the wave is reflected by the outer boundary and radiated to the interface. The wave amplitude attenuation function of the initial displacement wave excited at the interface is written as(15)ϕζ,lb=exp⁡-2ζlbbc1=exp⁡-2ζ¯lc1.If the wave amplitude attenuation at the interface is known, the constraint function of the size of bounded medium is transformed from function (15). Consider (16)lb=-ln⁡ϕ02ζ¯.


## 3.3. Constraint Function of Finite Element Mesh Density

In the finite element discretization of the bounded medium, the finite element mesh of only one row, denoted as the “base row,” is merely required.Sr denotes the dynamic stiffness matrix of one row of the finite element mesh with length r (Figure 3). The same is true for the condensed matrix obtained by eliminating all degrees of freedom except those related to interior nodes:(17)Sαbω=Sbαω,S~αbω=S~bαω.Figure 3
First row of finite element mesh.The formula relatesS~r1 for exterior nodes to the corresponding matrix of interior nodes S~r0:(18)S~r1ω=S~r0κω,in which κ=r1/r0=γ2.The dashpots are assigned at the exterior nodes of the base row. BothS~r0γ2n-2ω and S~r1γ2n-4ω are calculated by (18). After removing the dashpots, either S~r0γ2n-2ω or S~r1γ2n-4ω is assembled at the exterior nodes of the base row for frequency γ2n-4ω. S~r0γ2n-4ω is calculated by condensing again into the dynamic stiffness matrix for interior nodes. S~fω is obtained by repeating the last process n-2 and the proportionality factor γ must satisfy(19)r0+1/2Δ+γΔr0=γ2.


## 4. Value Optimization Model of Key Factors

A value optimization model of key factors is proposed to control the simulation accuracy and computational efficiency of the infinite soil-structure interaction. Then the optimum solutions of key factors of the DSEM are obtained by the value optimization.


### 4.1. Design Variable

Because the optimum values of key factors are provided by solving the value optimization model, the design variables of optimization model include the size of the bounded medium, the artificial damping ratio, the mesh size of finite element, and the proportionality factor.


### 4.2. Objective Function

The relationship between the size of bounded medium and the finite element mesh density is expressed as(20)lb=22Δb1-γ2n1-γ.Substituting (19) into (20), the objective function is transformed from (20).


### 4.3. Constraint Function

The constraint function system is composed of the constraint functions of the size of the bounded medium, the artificial damping ratio, and the finite element mesh density as described in Section2.The difference function between the dimensionless dynamic stiffness coefficient of the damped and undamped unbounded foundation is expressed as follows:(21)ΔSa0=S¯∞a0∗-S¯∞a0=1-a0∗2-1-a02.The results shown in Figure 4 indicate that the smaller artificial damping ratio results in the smaller difference of the dynamic stiffness coefficient and the maximum of the difference function is calculated at the point given by a0=1 (MaxS=ΔS(1)).Figure 4
Difference function of dimensionless dynamic stiffness coefficients.Taking the accuracy and computational cost of the soil-structure interaction into account, the initial conditions of the optimization model consist of the maximum of the difference function of the dynamic stiffness coefficient, the wave amplitude attenuation, and the maximum finite element number in the bounded medium. Therefore, the value optimization model of key factors based on the DSEM is expressed as follows:(1)
Initial conditions:MaxS, ϕ0, MaxN.(2)
Design variables:ζ-, l/b, Δ/b, γ.(3)
Objective function:(22)γ=exp⁡12nlnlb-1.(4)
Constraint functions:(23)absΔS1,ζ¯≤MaxS=ΔS1,lb=-lnϕ02ζ¯,1+12bΔ+γΔ=γ2,n≤MaxN.


## 4.1. Design Variable

Because the optimum values of key factors are provided by solving the value optimization model, the design variables of optimization model include the size of the bounded medium, the artificial damping ratio, the mesh size of finite element, and the proportionality factor.


## 4.2. Objective Function

The relationship between the size of bounded medium and the finite element mesh density is expressed as(20)lb=22Δb1-γ2n1-γ.Substituting (19) into (20), the objective function is transformed from (20).


## 4.3. Constraint Function

The constraint function system is composed of the constraint functions of the size of the bounded medium, the artificial damping ratio, and the finite element mesh density as described in Section2.The difference function between the dimensionless dynamic stiffness coefficient of the damped and undamped unbounded foundation is expressed as follows:(21)ΔSa0=S¯∞a0∗-S¯∞a0=1-a0∗2-1-a02.The results shown in Figure 4 indicate that the smaller artificial damping ratio results in the smaller difference of the dynamic stiffness coefficient and the maximum of the difference function is calculated at the point given by a0=1 (MaxS=ΔS(1)).Figure 4
Difference function of dimensionless dynamic stiffness coefficients.Taking the accuracy and computational cost of the soil-structure interaction into account, the initial conditions of the optimization model consist of the maximum of the difference function of the dynamic stiffness coefficient, the wave amplitude attenuation, and the maximum finite element number in the bounded medium. Therefore, the value optimization model of key factors based on the DSEM is expressed as follows:(1)
Initial conditions:MaxS, ϕ0, MaxN.(2)
Design variables:ζ-, l/b, Δ/b, γ.(3)
Objective function:(22)γ=exp⁡12nlnlb-1.(4)
Constraint functions:(23)absΔS1,ζ¯≤MaxS=ΔS1,lb=-lnϕ02ζ¯,1+12bΔ+γΔ=γ2,n≤MaxN.


## 5. Numerical Verification

### 5.1. Verification in the Frequency Domain

#### 5.1.1. Problem Chosen

A cross section shown in Figure5 is the intercepted plane of a rectangular rigid strip foundation of width 2b and depth b embedded in a homogeneous half plane with shear modulus G, Poisson ratio v, and artificial damping ratio ζ. The DOFs, vertical (V), horizontal (H), and rocking (R), describe the plane motion of the foundation. The procedure analyzes the dynamic stiffness matrix, S∞(a0), of the rigid strip on the damped unbounded medium in the frequency domain.Figure 5
Cross section of a strip foundation embedded in a homogeneous half plane.LetPn and Δn, n∈(V,H,R), respectively, denote the amplitudes of the harmonic force and displacement; the vectors are related by the dynamic stiffness matrix:(24)PVPHPRb=S~∞a0ΔVΔHbΔR=SVVa0000SHHa0SHRa00SRHa0SRRa0ΔVΔHbΔR.The dynamic stiffness coefficients of the bounded medium are decomposed into the spring coefficientk(a0) and the damping coefficient c(a0):(25)Sjja0=1Gkjja0+ia0cjja0j∈V,H,R.


#### 5.1.2. Numerical Evaluation

The discretization of the finite element mesh is shown in Figure6, in which the eight-node isoparametric finite element is employed. The finer mesh is utilized near the strip foundation and it becomes coarser away from that with a proportionality factor γ. The material constants of the foundation are Poisson ratio v=0.25 and natural damping ratio β=0.05. The various initial conditions of the value optimization model and the corresponding optimum solutions of key factors are listed in Table 2. The dynamic stiffness coefficient solved by case C1 is seen as the accurate soil-structure interaction.Table 2
Characteristics of the finite element discretization of bounded medium.
Case Initial condition Optimal value of key factor Max ⁡ S ϕ 0 Max ⁡ N ζ ¯ l / b n Δ / b γ C1 0.32 0 2000 0.05 3000 2000 0.0028 1.002  C2 1.12 1.00E − 04 25 0.6 7.68 25 0.062 1.044 C3 0.9 1.00E − 04 25 0.4 11.51 25 0.074 1.052 C4 0.64 1.00E − 04 25 0.2 23.03 25 0.094 1.066 C5 0.32 1.00E − 04 25 0.1 46.06 25 0.114 1.080  C6 1.12 1.00E − 04 25 0.6 7.68 25 0.062 1.044 C7 1.12 1.00E − 06 25 0.6 11.51 25 0.074 1.052 C8 1.12 1.00E − 13 25 0.6 24.95 25 0.096 1.068 C9 1.12 1.00E − 48 25 0.6 46.06 25 0.114 1.095  C10 1.12 1.00E − 13 6 0.6 23.03 6 0.429 1.304 C11 1.12 1.00E − 13 12 0.6 23.03 12 0.201 1.142 C12 1.12 1.00E − 13 18 0.6 23.03 18 0.131 1.093 C13 1.12 1.00E − 13 25 0.6 23.03 25 0.094 1.066Figure 6
Finite element mesh of bounded medium.As the wave amplitude attenuation and the maximum finite element number in the bounded medium are held constant, the maximum of the difference function of the dynamic stiffness coefficient is a variable value. According to the comparison of the dynamic stiffness coefficient of cases C1, C2, C3, C4, and C5 shown in Figure7, the smaller maximum of the difference function of the dynamic stiffness coefficient promotes the accuracy of the dimensionless dynamic stiffness coefficient of the unbounded medium.Figure 7
Effect of maximum of difference function of dynamic stiffness coefficient on dynamic stiffness coefficient.As the maximum of the difference function of the dynamic stiffness coefficient and the maximum finite element number in the bounded medium remain unchangeable, the wave amplitude attenuation is a variable value. According to the comparison of the dynamic stiffness coefficient of cases C1, C6, C7, C8, and C9 shown in Figure8, the larger wave amplitude attenuation improves the accuracy of the dimensionless dynamic stiffness coefficient of the unbounded medium.Figure 8
Effect of wave amplitude attenuation on dynamic stiffness coefficient.As the maximum of the difference function of the dynamic stiffness coefficient and the wave amplitude attenuation are changeless, the maximum finite element number in the bounded medium is a variable value. According to the comparison of the dynamic stiffness coefficient of cases C1, C10, C11, C12, and C13 shown in Figure9, the larger maximum finite element number reduces the computational efficiency but enhances the simulation precision of the dimensionless dynamic stiffness coefficient of the unbounded medium.Figure 9
Effect of maximum finite element number on dynamic stiffness coefficient.


### 5.2. Verification in the Time Domain

A transient horizontal displacement is excited at the interface of the discretization of the finite element shown in Figure10 to test the feasibility of the value optimization model of key factors in the time domain. The material constants of the medium are shear modular G=0.32 GPa, Poisson’s ratio v=0.25, a density of ρ=2000 kg/m3, shear-wave velocity cs = 400 m/s, dilatational wave velocity cp=693 m/s, and time step Δt=0.001 s.Figure 10
Finite element mesh of bounded media.The transient excitation is a prescribed horizontal displacement at the center of the rigid base (with zero values for the vertical and rocking motions)(26)ubt=u021-cos⁡2πtT0≤t≤2T0t>2Twith period T=8b/cs.The initial conditions and the optimum solutions by the value optimization listed in Table3 are realized by the damping solvent extraction method. The interaction force based on the case A1 is regarded as the exact interaction force of the foundation-structure.Table 3
Various initial conditions and the optimum values of key factors.
Case (Max⁡S, ϕ0, Max⁡N) ( l / b, ζ¯, n, Δ/b, γ)opt A1 (0.1, 0, 2000) (3000, 0.05, 2000, 0.0028, 1.002) A2 (1, 0.1, 25) (2.4, 0.48, 25, 0.035, 1.025) A3 (1, 0.05, 25) (3.2, 0.48, 25, 0.044, 1.031) A4 (1, 0.01, 25) (4.8, 0.48, 25, 0.051, 1.036) A5 (0.5, 0.05, 25) (12.4, 0.12, 25, 0.075, 1.053) A6 (0.3, 0.05, 25) (37.4, 0.04, 25, 0.107, 1.076)Since the maximum finite element number in the bounded medium remains unchanged, the computational efficiency of the unbounded foundation-structure interaction is fixed. The interaction force shown in Figure11 is calculated by the various initial conditions of the value optimization model of key factors. When the maximum value of the difference function of the dynamic stiffness coefficients is a variable value (such as in cases A1, A2, A3, and A4), the larger wave amplitude attenuation improves the precision of the interaction force; nevertheless, the over-large attenuation has a negligible effect. When the wave amplitude attenuation is changeless (such as in cases A1, A3, A5, and A6), the smaller maximum of the difference function of the dynamic stiffness coefficient increases the accuracy of the soil-structure interaction force.Figure 11
Time history of horizontal interaction force varying maximum of difference function of dynamic stiffness coefficient and wave amplitude attenuation.In summary, the implementation of the unbounded medium-structure interaction in the time-frequency domain illustrates that the alterable initial conditions of the value optimization model of key factors can improve the simulation accuracy and computational efficiency of the interaction, which are made up with the smaller maximum finite element number, the larger wave amplitude attenuation, and the smaller maximum of the difference function between the dimensionless dynamic stiffness coefficients of the damped and undamped unbounded foundation.


## 5.1. Verification in the Frequency Domain

### 5.1.1. Problem Chosen

A cross section shown in Figure5 is the intercepted plane of a rectangular rigid strip foundation of width 2b and depth b embedded in a homogeneous half plane with shear modulus G, Poisson ratio v, and artificial damping ratio ζ. The DOFs, vertical (V), horizontal (H), and rocking (R), describe the plane motion of the foundation. The procedure analyzes the dynamic stiffness matrix, S∞(a0), of the rigid strip on the damped unbounded medium in the frequency domain.Figure 5
Cross section of a strip foundation embedded in a homogeneous half plane.LetPn and Δn, n∈(V,H,R), respectively, denote the amplitudes of the harmonic force and displacement; the vectors are related by the dynamic stiffness matrix:(24)PVPHPRb=S~∞a0ΔVΔHbΔR=SVVa0000SHHa0SHRa00SRHa0SRRa0ΔVΔHbΔR.The dynamic stiffness coefficients of the bounded medium are decomposed into the spring coefficientk(a0) and the damping coefficient c(a0):(25)Sjja0=1Gkjja0+ia0cjja0j∈V,H,R.


### 5.1.2. Numerical Evaluation

The discretization of the finite element mesh is shown in Figure6, in which the eight-node isoparametric finite element is employed. The finer mesh is utilized near the strip foundation and it becomes coarser away from that with a proportionality factor γ. The material constants of the foundation are Poisson ratio v=0.25 and natural damping ratio β=0.05. The various initial conditions of the value optimization model and the corresponding optimum solutions of key factors are listed in Table 2. The dynamic stiffness coefficient solved by case C1 is seen as the accurate soil-structure interaction.Table 2
Characteristics of the finite element discretization of bounded medium.
Case Initial condition Optimal value of key factor Max ⁡ S ϕ 0 Max ⁡ N ζ ¯ l / b n Δ / b γ C1 0.32 0 2000 0.05 3000 2000 0.0028 1.002  C2 1.12 1.00E − 04 25 0.6 7.68 25 0.062 1.044 C3 0.9 1.00E − 04 25 0.4 11.51 25 0.074 1.052 C4 0.64 1.00E − 04 25 0.2 23.03 25 0.094 1.066 C5 0.32 1.00E − 04 25 0.1 46.06 25 0.114 1.080  C6 1.12 1.00E − 04 25 0.6 7.68 25 0.062 1.044 C7 1.12 1.00E − 06 25 0.6 11.51 25 0.074 1.052 C8 1.12 1.00E − 13 25 0.6 24.95 25 0.096 1.068 C9 1.12 1.00E − 48 25 0.6 46.06 25 0.114 1.095  C10 1.12 1.00E − 13 6 0.6 23.03 6 0.429 1.304 C11 1.12 1.00E − 13 12 0.6 23.03 12 0.201 1.142 C12 1.12 1.00E − 13 18 0.6 23.03 18 0.131 1.093 C13 1.12 1.00E − 13 25 0.6 23.03 25 0.094 1.066Figure 6
Finite element mesh of bounded medium.As the wave amplitude attenuation and the maximum finite element number in the bounded medium are held constant, the maximum of the difference function of the dynamic stiffness coefficient is a variable value. According to the comparison of the dynamic stiffness coefficient of cases C1, C2, C3, C4, and C5 shown in Figure7, the smaller maximum of the difference function of the dynamic stiffness coefficient promotes the accuracy of the dimensionless dynamic stiffness coefficient of the unbounded medium.Figure 7
Effect of maximum of difference function of dynamic stiffness coefficient on dynamic stiffness coefficient.As the maximum of the difference function of the dynamic stiffness coefficient and the maximum finite element number in the bounded medium remain unchangeable, the wave amplitude attenuation is a variable value. According to the comparison of the dynamic stiffness coefficient of cases C1, C6, C7, C8, and C9 shown in Figure8, the larger wave amplitude attenuation improves the accuracy of the dimensionless dynamic stiffness coefficient of the unbounded medium.Figure 8
Effect of wave amplitude attenuation on dynamic stiffness coefficient.As the maximum of the difference function of the dynamic stiffness coefficient and the wave amplitude attenuation are changeless, the maximum finite element number in the bounded medium is a variable value. According to the comparison of the dynamic stiffness coefficient of cases C1, C10, C11, C12, and C13 shown in Figure9, the larger maximum finite element number reduces the computational efficiency but enhances the simulation precision of the dimensionless dynamic stiffness coefficient of the unbounded medium.Figure 9
Effect of maximum finite element number on dynamic stiffness coefficient.


## 5.1.1. Problem Chosen

A cross section shown in Figure5 is the intercepted plane of a rectangular rigid strip foundation of width 2b and depth b embedded in a homogeneous half plane with shear modulus G, Poisson ratio v, and artificial damping ratio ζ. The DOFs, vertical (V), horizontal (H), and rocking (R), describe the plane motion of the foundation. The procedure analyzes the dynamic stiffness matrix, S∞(a0), of the rigid strip on the damped unbounded medium in the frequency domain.Figure 5
Cross section of a strip foundation embedded in a homogeneous half plane.LetPn and Δn, n∈(V,H,R), respectively, denote the amplitudes of the harmonic force and displacement; the vectors are related by the dynamic stiffness matrix:(24)PVPHPRb=S~∞a0ΔVΔHbΔR=SVVa0000SHHa0SHRa00SRHa0SRRa0ΔVΔHbΔR.The dynamic stiffness coefficients of the bounded medium are decomposed into the spring coefficientk(a0) and the damping coefficient c(a0):(25)Sjja0=1Gkjja0+ia0cjja0j∈V,H,R.


## 5.1.2. Numerical Evaluation

The discretization of the finite element mesh is shown in Figure6, in which the eight-node isoparametric finite element is employed. The finer mesh is utilized near the strip foundation and it becomes coarser away from that with a proportionality factor γ. The material constants of the foundation are Poisson ratio v=0.25 and natural damping ratio β=0.05. The various initial conditions of the value optimization model and the corresponding optimum solutions of key factors are listed in Table 2. The dynamic stiffness coefficient solved by case C1 is seen as the accurate soil-structure interaction.Table 2
Characteristics of the finite element discretization of bounded medium.
Case Initial condition Optimal value of key factor Max ⁡ S ϕ 0 Max ⁡ N ζ ¯ l / b n Δ / b γ C1 0.32 0 2000 0.05 3000 2000 0.0028 1.002  C2 1.12 1.00E − 04 25 0.6 7.68 25 0.062 1.044 C3 0.9 1.00E − 04 25 0.4 11.51 25 0.074 1.052 C4 0.64 1.00E − 04 25 0.2 23.03 25 0.094 1.066 C5 0.32 1.00E − 04 25 0.1 46.06 25 0.114 1.080  C6 1.12 1.00E − 04 25 0.6 7.68 25 0.062 1.044 C7 1.12 1.00E − 06 25 0.6 11.51 25 0.074 1.052 C8 1.12 1.00E − 13 25 0.6 24.95 25 0.096 1.068 C9 1.12 1.00E − 48 25 0.6 46.06 25 0.114 1.095  C10 1.12 1.00E − 13 6 0.6 23.03 6 0.429 1.304 C11 1.12 1.00E − 13 12 0.6 23.03 12 0.201 1.142 C12 1.12 1.00E − 13 18 0.6 23.03 18 0.131 1.093 C13 1.12 1.00E − 13 25 0.6 23.03 25 0.094 1.066Figure 6
Finite element mesh of bounded medium.As the wave amplitude attenuation and the maximum finite element number in the bounded medium are held constant, the maximum of the difference function of the dynamic stiffness coefficient is a variable value. According to the comparison of the dynamic stiffness coefficient of cases C1, C2, C3, C4, and C5 shown in Figure7, the smaller maximum of the difference function of the dynamic stiffness coefficient promotes the accuracy of the dimensionless dynamic stiffness coefficient of the unbounded medium.Figure 7
Effect of maximum of difference function of dynamic stiffness coefficient on dynamic stiffness coefficient.As the maximum of the difference function of the dynamic stiffness coefficient and the maximum finite element number in the bounded medium remain unchangeable, the wave amplitude attenuation is a variable value. According to the comparison of the dynamic stiffness coefficient of cases C1, C6, C7, C8, and C9 shown in Figure8, the larger wave amplitude attenuation improves the accuracy of the dimensionless dynamic stiffness coefficient of the unbounded medium.Figure 8
Effect of wave amplitude attenuation on dynamic stiffness coefficient.As the maximum of the difference function of the dynamic stiffness coefficient and the wave amplitude attenuation are changeless, the maximum finite element number in the bounded medium is a variable value. According to the comparison of the dynamic stiffness coefficient of cases C1, C10, C11, C12, and C13 shown in Figure9, the larger maximum finite element number reduces the computational efficiency but enhances the simulation precision of the dimensionless dynamic stiffness coefficient of the unbounded medium.Figure 9
Effect of maximum finite element number on dynamic stiffness coefficient.


## 5.2. Verification in the Time Domain

A transient horizontal displacement is excited at the interface of the discretization of the finite element shown in Figure10 to test the feasibility of the value optimization model of key factors in the time domain. The material constants of the medium are shear modular G=0.32 GPa, Poisson’s ratio v=0.25, a density of ρ=2000 kg/m3, shear-wave velocity cs = 400 m/s, dilatational wave velocity cp=693 m/s, and time step Δt=0.001 s.Figure 10
Finite element mesh of bounded media.The transient excitation is a prescribed horizontal displacement at the center of the rigid base (with zero values for the vertical and rocking motions)(26)ubt=u021-cos⁡2πtT0≤t≤2T0t>2Twith period T=8b/cs.The initial conditions and the optimum solutions by the value optimization listed in Table3 are realized by the damping solvent extraction method. The interaction force based on the case A1 is regarded as the exact interaction force of the foundation-structure.Table 3
Various initial conditions and the optimum values of key factors.
Case (Max⁡S, ϕ0, Max⁡N) ( l / b, ζ¯, n, Δ/b, γ)opt A1 (0.1, 0, 2000) (3000, 0.05, 2000, 0.0028, 1.002) A2 (1, 0.1, 25) (2.4, 0.48, 25, 0.035, 1.025) A3 (1, 0.05, 25) (3.2, 0.48, 25, 0.044, 1.031) A4 (1, 0.01, 25) (4.8, 0.48, 25, 0.051, 1.036) A5 (0.5, 0.05, 25) (12.4, 0.12, 25, 0.075, 1.053) A6 (0.3, 0.05, 25) (37.4, 0.04, 25, 0.107, 1.076)Since the maximum finite element number in the bounded medium remains unchanged, the computational efficiency of the unbounded foundation-structure interaction is fixed. The interaction force shown in Figure11 is calculated by the various initial conditions of the value optimization model of key factors. When the maximum value of the difference function of the dynamic stiffness coefficients is a variable value (such as in cases A1, A2, A3, and A4), the larger wave amplitude attenuation improves the precision of the interaction force; nevertheless, the over-large attenuation has a negligible effect. When the wave amplitude attenuation is changeless (such as in cases A1, A3, A5, and A6), the smaller maximum of the difference function of the dynamic stiffness coefficient increases the accuracy of the soil-structure interaction force.Figure 11
Time history of horizontal interaction force varying maximum of difference function of dynamic stiffness coefficient and wave amplitude attenuation.In summary, the implementation of the unbounded medium-structure interaction in the time-frequency domain illustrates that the alterable initial conditions of the value optimization model of key factors can improve the simulation accuracy and computational efficiency of the interaction, which are made up with the smaller maximum finite element number, the larger wave amplitude attenuation, and the smaller maximum of the difference function between the dimensionless dynamic stiffness coefficients of the damped and undamped unbounded foundation.


## 6. Application

To demonstrate the adaptation of the value optimization model of key factors in practical engineering, the rock foundation-gravity dam shown in Figure12 is applied to analyze the dynamic response of the dam to harmonic ground motion. The geometric properties of the gravity dam are the height, 132 m, and the width, which is 14 m at the crest and 110 m at the bottom. Taking the cross section of the medium-dam as an example, the dynamic response analysis of the gravity dam is calculated as a plane-strain problem. The material parameters of the gravity dam are Young’s modulus of 20 GPa, Poisson’s ratio of 0.18, and a density of 2600 kg/m3. The material parameters of the rock foundation are Young’s modulus of 2.4 GPa, Poisson’s ratio of 0.333, a density of 2100 kg/m3, and the characteristic length of 40 m. The constants of the displacement input harmonic wave are the amplitude, 0.2 m, and the cycle, 0.4 s. Meanwhile, the various initial conditions of the value optimization model of key factors based on the DSEM are implemented to analyze the dynamic response of the gravity dam.Figure 12
Finite element mesh of the gravity dam-bounded medium.The displacement function of the dam at the crest shown in Figure13 indicates that the simulation precision and computational efficient of the displacement response of the gravity dam can be enhanced by reducing the maximum of the difference function of the dynamic stiffness coefficient and raising the wave amplitude attenuation at the interface and the maximum of finite element number in the bounded medium.Figure 13
Dynamic magnification of crest displacements under seismic excitationT=0.4 s.


## 7. Conclusion

Considering the simulation accuracy and computational efficiency of the soil-structure interaction, this study’s aim was to establish the value optimization model of key factors (such as the artificial damping ratio, the size of the bounded medium, and the finite element mesh density). From the impact assessment of the key factors on the dimensionless dynamic stiffness coefficient and the analysis of the interrelation among the key factors, the value optimization model of key factors based on the DSEM was proposed. In additions, the initial conditions of the value optimization model were composed of the wave amplitude attenuation, the maximum finite element number, and the maximum of the difference function between the dimensionless dynamic stiffness coefficient of the undamped and damped unbounded medium.The implementations of the dimensionless dynamic stiffness coefficient and the interaction force demonstrated that the value optimization model of key factors was capable of controlling the simulation accuracy and computational efficiency of the unbounded medium-structure interaction in the frequency and time domain. Finally, the displacement history function at the crest of the gravity dam on the rock foundation indicated that the value optimization model of key factors based on the DSEM was appropriate for analyzing the interaction of the practical unbounded medium-structure.



---

*Source: 3240720-2016-07-12.xml*