Abstract

The interaction between urban river-canyon topography and the river-side building is investigated by using a whole analytic model of a semicircle river-canyon and a shear wall supported by a semicircle rigid foundation embedded in a homogenous half-space. The closed-form analytical solution for system response is presented based on the wave function expansion method. The analysis focuses on the effects of the canyon-building interaction on system response. The strength of the interaction between the river-canyon topography and the building changes periodically as the distance between the canyon and the structure increases, leading to the interaction having beneficial or harmful effects on the building’s seismic response. The foundation peak response of the building can be amplified by about 10%, and the peak of the building relative response can be amplified by about 40%. The distribution of canyon-structure spacing with strong or weak interaction is closely related to the dynamic characteristics of the building and the incident angle of the wave. When designing buildings along the river, the building and canyon should be analyzed as a whole model to determine whether the location of the building is in a position with strong interaction with the river-canyon. The model in this paper may be useful for obtaining insight into the effects of canyon-structure interaction and interpreting the observed response in buildings and seismic response estimation in general.

1. Introduction

This paper investigates the interaction between urban river-canyon and river-side building by using an analytical model. The work in this article is an extension of Trifunac’s researches on the seismic response analysis of a single building [1] and a single canyon [2]. The problem studied in this paper involves the seismic interaction between urban local irregular topography and nearby buildings, which is currently termed as site-city interaction [35].

In the past, the amplification effect of canyon topography on ground motions received widespread attention [624]. However, scholars considered this amplification effect of canyon topography on building responses only through studying the ground motion characteristics of canyon sites rather than directly seeing the canyon and the nearby buildings as a whole model to study the interaction. Actually, when a strong earthquake occurs, seismic waves scatter on the canyon boundary and the surface of the building foundation. Their scattering wave fields will couple with each other, forming the complex canyon-structure interaction system. Therefore, the ground motion characteristics of the canyon site and the building seismic response are no longer determined only by the seismic response characteristics of the canyon site. The correct results may not be obtained if only the canyon topography and the surrounding soil are taken out for analysis. It is reasonable to analyze the river-canyon, the soil, and the buildings along the river as a whole model.

This paper proposes a two-dimensional overall analytical model of the interaction between the canyon topography and the buildings along the river. This model consists of a 2D infinite long shear wall supported by a rigid foundation along a river-canyon embedded in a homogenous half-space, subjected by incident SH wave. The closed-form analytical solution of the canyon-structure interaction system response is obtained. The closed-form analytical solution is further used to reveal the interaction mechanism of the canyon-structure system.

In the next section, the methodology is presented, followed by the verification through comparison with the published results. Then, the results in the frequency domain are presented and the anti-plane building response and canyon motion are discussed. Finally, the main findings and the conclusions are summarized.

2. Methodology

2.1. Analytical Model

The dynamic interaction model of canyon-structure, as shown in Figure 1, consists of a river-canyon and a shear wall, infinitely long in the y-direction, standing on a rigid foundation embedded in a half-space. With the points of and as the origin, there are two sets of coordinate systems. The rectangular coordinates and the corresponding cylindrical coordinates are drawn. The y-axis of the rectangular coordinate system is perpendicular to the paper (points to the reader). The river-canyon and the rigid foundations are semicircular with the centers at points of and , respectively, and the radii are a and b, respectively. The half-space and the shear wall are assumed to be linearly elastic, uniform, and isotropic, and the boundary between the rigid foundation and the half-space is assumed to be tightly connected. The half-space is marked with shear modulus , shear wave velocity , and mass density . Within the unit length along the y-axis, the mass of the rigid foundation is . The shear wall has shear wave velocity , height , undamped natural frequency , and mass per unit length in the y-direction. The distance between the canyon center and the foundation center is . The excitation is an incident plane SH wave with circular frequency , incidence angle measured from the vertical z-axis, and unit amplitude.

To illustrate the canyon deformation, Figure 2 shows 5 observation points on the canyon which are chosen to compute their motions. These are point 1 , point 2 , point 3 , point 4 , and point 5 , with point 3 being at the bottom of the canyon.

2.2. Governing Equations and Solution

In the local coordinate system (--) of the canyon, the motion in the half-space along the y-axis direction caused by the incident SH wave can be written aswhere cx and cz represent the phase velocities of waves along the x- and z-axes, respectively

The free-field motion consists of the incident wave and the reflected wave on the surface of half-space.

In the area near the canyon and foundation, the wave field is composed of three parts, which are (A) the free-field motion expressed by equation (3), (B) the scattering field of waves caused by the canyon boundary at , and (C) the scattering field generated by the rigid foundation boundary at . Since (B) and (C) represent the wavefield radiating outward in the half-space, the form of their final solution must be the same. The time factor is omitted in the following solutions.

All the wave motions , , and must satisfy the wave differential equationwhere stands for the three wave motions , , and , respectively. The three wave motions should satisfy the boundary conditions given by

In the local coordinate system (--) of the canyon, the free-field motion in the half-space along the y-axis caused by the incident SH wave can be expanded into Fourier-Bessel series as

Similarly, in the local coordinate system (--), the free-field motion can also be expanded into

Here, Jn(x) is the Bessel function of the first kind with argument x and order n, and .

The scattering wave motions near the canyon and the two foundations can be written as

Here, and are the scattering wave motions from the canyon and the foundation, respectively, and is the Hankel function of the second kind with argument x and order n.

Equations (8)–(11) can be proved to satisfy equation (4) and the boundary condition equation (5). An, Bn, Cn, and Dn are the unknown complex coefficients that should be determined by the boundary conditions equations (6) and (7). Δ is the effective input motion of the rigid foundation to the superstructure.

Because equations (10) and (11) are expressed in different local cylindrical coordinate systems, equations (10) and (11) should be transformed between any two local cylindrical coordinate systems as shown in Figure 1. Based on Graf’s Addition Theorem [25, 26], equations (10) and (11) can be transformed between any two cylindrical coordinate systems in Figure 1. The details of the derivation of the coordinate transformation formula are omitted.

Putting equations (8), (10), and (11) into equation (6), we obtain the following equations:

Submitting equations (9)– (11) into equation (7), we obtain the following equations:

Here, are presented as follows:where, in equation (18e), replacing C with H stands for coordinate transformation for Hankel function.

The foundation effective input motion Δ can be determined by the dynamic equilibrium equation of the rigid foundation.where, , the action per unit length of the shear wall on the foundation, is given by equation (20). is the action of the soil on the foundations, given by equation (21).

By submitting equations (20) and (21) into equation (19), it can be shown that

The unknown complex coefficients An, Bn, Cn, Dn, and Δ are determined by solving the linear equations (12)–(17) and (22) together, and the closed-form analytical solution of the problem can be obtained. Particularly, the relative responses of the shear wall to the foundation can be expressed as

To compare with the published results, Δ and are normalized by the surface displacement amplitude of free-field ground motion as

2.3. Impedance Function of the Rigid Foundation

According to the definition of impedance function, the foundation has unit displacement response () at this time. The wave field in the half-space only has canyon scattering field and foundation scattering field , but does not include free-field wave motion,

The boundary conditions are

Compared with the derivation in Section 2.2, the unknown coefficients , , , and can be obtained by solving the linear equations (12)∼(17) and (22) but removing the terms that correspond to the free-field. After that, the foundation impedance function Kyy can be obtained by integrating the stress on the foundation boundary.

After the normalization by the shear modulus of the half-space, the foundation impedance function Kyy can be further written as follows:where and are the dimensionless stiffness and radiation damping coefficients.

3. Method Verification

3.1. Dimensionless Parameters

The results in this paper can be seen as functions of the governing dimensionless parameters, which are expressed in terms of the characteristic length of the problem, a, and reference wave velocity, β. The dimensionless frequency is defined as [27]where is canyon width and is the wavelength of the shear waves in the half-space. Further, based on [28], the dimensionless parameter which describes the relative flexibility of the shear wall to surrounding soil is defined as

The dimensionless parameter has values 0 ≤ ε < ∞. For a given ε, the first-order vibration frequency of the shear wall is at . For a very stiff structure and/or very soft soil, ε = 0 (corresponding to a rigid shear wall added as a rigid mass to the foundation), and for a very flexible structure and/or rigid soil, ε ⟶ ∞.

Other dimensionless parameters used in this paper are the dimensionless distance between the structure and the canyon, , and the ratios, and , where is the mass of the soil per unit length in the y-direction that has been replaced by the foundation.

3.2. Verification

To show the verification of our solution, we set the canyon radius equal to the foundation radius and the dimensionless distance between the canyon and the foundation tends to be infinite, to verify whether our results of the river-side building and the canyon are consistent with the published results that correspond to single building [1] and single canyon [2].

Figures 3 and 4 show the comparison between our results and the published results when . The figures demonstrate that our results agree well with the published results, indicating the verification of our solution.

4. Numerical Results and Analysis

4.1. Model Parameters

Figures 510, respectively, show foundation impedance function, foundation responses, the shear wall’s relative response, and the canyon response, at different canyon-structure spacing. For easy comparison with the published results, the responses of single structure [1] and single river-canyon [2] are also shown in each subfigure. The parameters used in this section are as follows. The radius of the canyon is equal to the foundation radius; that is, a = b. The dimensionless distance between the canyon and the river-side building (referred to as “canyon-structure spacing” for short)  = 2.5, 3, 4, 5, and 10. The ratios of the foundation mass densities to the half-space mass density are 1/2, which are equivalent to  = 1/2. The mass ratio of the shear wall to the foundation per unit length is  = 1, 2, and 4. The stiffness parameters of the shear walls are ε = 0, 2, and 4. The incident angles of the waves are γ = 0° (vertical incident) and 90° (horizontal incident).

4.2. Foundation Impedance Function

As shown in Figure 5, the stiffness coefficient and damping coefficient of the foundation impedance function are highly related to the canyon-structure spacing () and the frequency of the incident wave (). First, the stiffness coefficient of the foundation impedance function will fluctuate strongly with different canyon-structure spacing and the frequency of the incident wave. For the closer distance between the canyon and the structure ( = 2.5, 3, 4, and 5), as the frequency of the incident wave increases, the stiffness coefficient changes more widely. The change in the damping coefficient of the foundation impedance function is similar to the stiffness coefficient . When the frequency of the incident wave is small, the damping coefficient is roughly the same as that of the single rigid foundation [29]. As the frequency of the incident wave increases, the damping coefficient will also fluctuate around the damping coefficient of the single rigid foundation impedance function (). However, the amplitude of the fluctuation of the damping coefficient is smaller than that of the stiffness coefficient .

4.3. Building Responses for Different Canyon-Structure Spacing

Figure 6 shows how the foundation response varies with the canyon-structure spacing () under different incident angles (). Observing Figure 6, we can see that, under different canyon-structure spacing, the difference between the peak of foundation response and the peak of single foundation response is small. This can be explained as follows. First, when the frequency of the incident wave is small, the wavelength of the wave is longer, resulting in strong wave diffraction. Therefore, even if the foundation is close to the canyon, the incident wave cannot feel well the existence of the canyon due to wave diffraction, which leads to a weak interaction between the canyon and the foundation. The foundation vibrates like a single foundation, so the peak of the foundation response has a small difference from that of the single foundation. Secondly, when the distance between the canyon and the structure is large, the system interaction is also naturally weaker. From this, the two factors that determine the strength of the interaction between the canyon and the foundation can be known: the canyon-structure spacing and the wavelength of the incident wave (that is, the frequency of the incident wave). The proof for this is as follows: for example, when η = 0.5, the wavelength of the incident wave is . For canyon-structure spacing  = 2.5, 3, 4, and 5, the wavelength is the same order of magnitude and close to the canyon-structure spacing. Therefore, foundation response corresponding to these locations under η = 0.5 deviates significantly from that of the single rigid foundation.

It can be seen from Figure 7 that the canyon-structure interaction has a great influence on the peak of structure relative response, which is mainly reflected in the case of structural stiffness . As shown in Figure 7(a), when and , for different superstructure masses  = 1, 2, and 4, the peaks of the structure relative responses are 10.910, 5.377, and 2.648, respectively, while the corresponding peaks of the single structure relative response are 7.825, 3.857, and 1.893, respectively, which are enlarged by 39.4%, 39.4%, and 39.9%, respectively. Besides, observing the results under different incident angles (), it can be found that when the wave is incident obliquely (), the peak of the structure relative response has a small difference from the corresponding result of the single building. However, when the wave is incident perpendicularly (), the difference is large. The influence of the canyon-structure interaction on the structure relative response is mainly reflected in the case of vertical incidence. It can be concluded from this that another factor that determines the strength of the interaction between the canyon and the river-side building is the incident angle of the wave. The influence of the incident angle () on the system interaction cannot be ignored. The deviation magnitude of the foundation response in Figure 6 and the structure relative response in Figure 7 from the single foundation result under the horizontal incidence are significantly smaller than those of the vertical incidence of the wave. This phenomenon can be called the “shielding effect.” In other words, when the wave is incident obliquely, the left canyon has the effect of shielding the incident wave reaching the foundation, leading to weak foundation response and structure relative response.

4.4. Peak of Structure Relative Response under Different Canyon-Structure Spacing

Figures 8 and 9, respectively, show the variation of the foundation peak response and the shear wall’s peak relative response with the canyon-structure spacing (). As shown in Figures 8 and 9, both the foundation peak response and the shear wall’s peak relative response show an attenuation trend as the canyon-structure spacing increases, but the foundation peak response and the shear wall’s peak relative response are always greater than those results of single building model [1] ( the gray line in Figures 8 and 9). This phenomenon is still obvious even when the distance between the canyon and the structure is very large ( = 40∼50). This illustrates that even when the canyon and the structure are far away, there still is an obvious interaction between them. The peak of the foundation response can be enlarged by approximately 10% (corresponding results of ε = 2 and  = 2 in Figure 8), and the peak of the shear wall’s relative response can be amplified by approximately 39.4% (corresponding results of ε = 2 and  = 1 and 2 in Figure 9).

The degree to which the color line deviates from the gray line in Figure 9 can be regarded as the degree of interaction between the river-canyon and the building. It can be seen from Figure 9 that the strength of the interaction between the river-canyon topography and the building changes periodically as the distance between the canyon and the structure increases, leading to the interaction having beneficial or harmful effects on the building’s seismic response. In many valley-structure intervals, the results of the valley-structure interaction model are greater than the results of single building model [1]. The distribution of valley-structure spacing with strong or weak interaction is closely related to the dynamic characteristics of the building itself and the incident angle of the wave. Therefore, when designing buildings along the river, the building and canyon should be modeled and analyzed as a whole to determine whether the location of the building is in a position with strong interaction with the river-canyon.

4.5. Canyon Response

Figure 10 shows the responses of canyon observation points under different incident angles and canyon-structure spacing. When the waves are vertical incident, because point 5 of the canyon is the closest to the building, the response of point 5 is the strongest, especially when the building is close to the river (D/a = 2.5, 3, 4, and 5). When the wave is incident obliquely, since point 1 is on the side of the incoming wave, the response of point 1 becomes the strongest, while the dynamic response at point 5 is less than that at point 1 because the river-canyon shields the wave.

Another important phenomenon is that canyon response is affected by the dynamic characteristics of the buildings, especially for small canyon-structure spacing (D/a = ≤5). For example, the response of point 5 becomes zero at the fundamental frequency of the building. This phenomenon is very similar to the foundation response; that is, at the fundamental frequency of the building, the foundation response tends to be zero. As mentioned in Section 3.1, for given ε, the first-order vibration of the structure is at ηb1 = 1/(2ε). In Figure 10, since the structural stiffness parameters are ε = 4, the first-order vibration frequency of the building is ηb1 = 0.125. It can be seen from Figure 10 that, for D/a = <5 and η = 0.125, the response of point 5 tends to be zero. Even when the buildings are all far from the canyon (D/a = 10), this trend can still be observed.

5. Conclusions

This paper established a 2D model for the anti-plane dynamic interaction between the urban river-canyon and the river-side building and presented a closed-form analytical solution for the system response based on the wave function expansion method. The main features of the interaction between the river-canyon and the building along the river are investigated by the analytical solution. The conclusions and findings are as follows:(1)The canyon-structure interaction depends not only on the distances between the canyon and the river-side building but also on the frequency of the incident wave and the incident angle. When the wave is incident obliquely, the canyon on the incoming wave side has the shielding effect on the incident wave.(2)The strength of the interaction between the river-canyon topography and the building changes periodically as the distance between the canyon and the structure increases, leading to the interaction having beneficial or harmful effects on the building’s seismic response. The results in this paper show that the foundation peak response of the building can be amplified by about 10%, and the peak of the building relative response can be amplified by about 40%.(3)The distribution of valley-structure spacing with strong or weak interaction is closely related to the dynamic characteristics of the building and the incident angle of the wave. When designing buildings along the river, the building and canyon should be modeled and analyzed as a whole to determine whether the location of the building is in a position with strong interaction with the river-canyon.(4)The canyon response of the canyon is closely related to the frequency of the incident wave, the canyon-structure spacing, the incident angle of the wave, and the dynamic characteristics of the building along the river.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that the work is original research that has not been published previously, and not under consideration for publication elsewhere. No conflicts of interest exist in the submission of this paper.

Authors’ Contributions

This study was approved by all authors for publication.

Acknowledgments

This study was supported by the National Natural Science Foundation of China under Grants 51808290 and U1839202, which are gratefully acknowledged.