Coupling Modal Analysis with the BEM for the Transient Response of Bar Structures Interacting with Three-Dimensional Soil Profiles

Abstract This work investigates the transient response of bar structures interacting with three-dimensional soil profiles. The structures are modeled by the Finite Element Method (FEM) and the soil models are described by a three-dimensional Boundary Element Formulation (BEM) in the frequency domain. A classic modal analysis is performed on the structure in terms of the relative displacements with respect to the soil. The dynamic response of the structure is coupled to the soil response, aiming to obtain frequency response functions (FRFs) of the soil-structure system. A new set of modal parameters are extracted from the FRFs of the coupled system. These new parameters allow for the synthesis of a set of orthogonal differential equations in the time domain. These equations are integrated by a classical numerical scheme resulting in the transient response of the structure interacting with the supporting soil. It is shown that for soil profiles that present eigenfrequencies, the system modal basis must be expanded to properly include the soil dynamics. The cases of a structure interacting with a homogeneous half-space and with a horizontal layer over a rigid stratum are considered. The results presented for both soil models are consistent. Graphical Abstract


INTRODUCTION
The dynamic response of structures interacting with the supporting soil profiles, called Dynamic Soil-Structure Interaction (DSSI), has received the attention of many researchers throughout the last four decades and still is a topic of continuous research.The most significant characteristic of DSSI problems is related to the unbounded dimensions that the soil profile presents.The existence of, at least, one unbounded dimension introduces an effect known as geometric or radiation damping, which is related to the energy that is withdrawn from the structure-foundation system in the form of outgoing and non-reflected waves.So, any attempt to model the dynamics of soils or unbounded domains by numerical methods must be able to take into account the radiation damping, or the Sommerfeld radiation condition.Typical domain-type numerical methods, such as the Finite Element Method (FEM) of the Finite Difference Method (FDM) present finite meshes and are not able to, automatically, consider the geometric damping.The FEM has established itself as the most versatile and used numerical simulation method in solid mechanics.So many attempts have been made to include radiation damping into FEM schemes presenting finite meshes.Examples of these strategies are in the inclusion of "Infinite Elements", the development of Dirichlet to Newman (DtN) and of Perfectly Matched Layers (PML) schemes to model the unbounded, wave propagating domain (Mesquita and Pavanello, 2005;Zhang et al., 2019).
In the last decades, the Boundary Element Method (BEM) has established itself as the most efficient and accurate method to model problems of wave propagation in unbounded domains and as such the modeling of the dynamic response of unbounded soil profiles.When formulated with the proper "fundamental solution" the BEM only requires the discretization of the surface of the domain being analyzed, reducing the meshing efforts by one dimension.More important for the present study is that the BEM can automatically take into account the radiation condition when dealing with the dynamics of unbounded domains (Dominguez, 1993).
There is a significant amount of research based on the BEM devoted to describing the dynamic response of structures interacting with soil profiles.In these dynamic analyses, the approaches of frequency-domain (stationary) or time-domain (transient) methodologies also apply.Frequency-domain analysis for DSSI problems is well established.The most efficient scheme to analyze DSSI problems is to model the soil by the BEM and the structure by the FEM (Spyrakosa and Xu, 2003;Mehdizadeh et al., 2021).Time domain analysis using the BEM had a start with the works of Mansur andBrebbia (1982a, 1982b) and has seen a constant evolution ever since.Some important contributions were made by Schanz and Antes (1997), in which the "Operational Quadrature Method" was used to obtain transient BEM solutions.Yerli et al. (1998) used a coupling of finite and infinite elements within the realm of the FEM to obtain the transient solution of unbounded domains.A direct coupling of the FEM with the BEM in the time domain, with the prescription of kinematic compatibility and equilibrium conditions at the soil-structure interface was presented by von Estorff and Prabuki (1990).The coupling of the FEM-BEM domains by a staggered solution was proposed by Rizos and Wang (2002).A time-domain iterative coupled scheme for BEM-FEM was presented by Soares et al. (2004).More recently, the Perfectly Matched Layer (PML) has been used in the FEM methodology to model the transient analysis of structures interacting with horizontally layered soil profiles (Zhang et al., 2019).All the previously mentioned efforts that describe the transient response of structures interacting with unbounded soil profiles are computationally very expensive, which makes the simulation of complex realistic problems almost unattainable.Moreover, the obtained transient responses tend to grow inaccurate for large time periods of analysis.
The aim of the present article is to report the development of an alternative strategy to synthesize the transient response of structures interacting with soils.It builds upon an idea presented by Wu and Smith (1995) in which a modal analysis of the structure was performed with respect to structure displacements relative to the foundation-soil degrees of freedom.The analysis was performed in the frequency domain.The system excitation was in the form of waves impinging upon the foundation.After the coupling of the structure with the soil-foundation subsystem, a set of frequency response functions (FRFs) could be synthesized for the coupled soil-structure system.In that analysis, it was possible to choose arbitrarily the number of structural eigenfrequencies and eigenmodes that would be used to synthesize the FRFs of the coupled system.This is important because in the dynamic response of buildings, the higher structural modes do not play a significant role in the dynamic response.To obtain a transient response to an earthquake excitation Wu and Smith (1995) did multiply the spectrum of the earthquake excitation with the FRFs of the coupled system and performed an inverse FFT (IFFT) operation.So, in this work the FFT algorithm was used to obtain the transient response and the system excitation was given by prescribed incoming waves.Louzada et al. (2019aLouzada et al. ( , 2019b) ) extended the work of Wu and Smith (1995) to consider external excitations acting directly on the structure.The work also considered the influence of distinct supporting mechanisms on the structural response.Examples of the supporting models were the fixed base, a linear spring, the soil modeled as a half-space, a layer over a rigid bedrock and also a structure supported by a pile embedded in the soil.The transient responses were obtained by the use of the IFFT on the stationary, frequency, responses.
Latin American Journal of Solids and Structures, 2023, 20(6), e502 3/16 The work presented by Ferraz (2021a) was also based on the FRFs of the coupled soil-structure system and was also able to consider any arbitrary number of modes to describe the structural response but did not use the FFT algorithm to obtain the transient response.From the FRFs of the coupled system, a new set of modal parameters were extracted to build a novel modal basis.This new modal basis was used to synthesize a new set of orthogonal differential equations of motion in the time domain.The direct integration of this uncoupled set of differential equations rendered the transient response of the structure considering the soil influence for the case of external excitations applied directly to the structure.
The present article enlarges the previously described works by allowing the soil model to be more complex and present, itself, eigenfrequencies.The typical model of an unbounded domain presenting eigenfrequencies and eigenmodes is a horizontal layer supported by a rigid bedrock.Now both subsystems present eigenfrequencies.And it is no longer possible to reproduce the dynamics of the coupled system by only considering the number of eigenfrequencies of the original structural system.An expanded modal basis is required to describe the coupled dynamics of both systems.This article uses this expanded modal basis concept to describe the transient dynamic behavior of structures interacting with layer over bedrock.The parameters of the extended modal basis are extracted from the FRFs of the coupled soil-structure system.The results are compared with those stemming from structure interaction with the homogeneous half-space.
Section 2 describes the statement of the problem, the BEM models used to describe the soil response as well as the structural equations of motion in terms of displacements relative to the foundation degrees of freedom.The soil-structure coupling procedure, the expansion of the modal basis and the methodology used to extract modal parameters from the FRFs are described next.Section 3 is dedicated to present numerical examples and conclusions are presented in section 4.

STATEMENT OF THE PROBLEM
The typical problem addressed in the present article is presented in Figure 1a, and consists of a linear structure with N degrees of freedom (DOFs) interacting with a soil profile.The structure is modelled using the Finite Element Method (FEM) considering bar elements.Only the displacement in the vertical direction, z u , is considered.The soil model depicted in Figure 1a is the classical homogeneous three-dimensional half-space.Nevertheless, in the present work, a horizontal soil layer with thickness, s h , supported by a rigid bedrock is also analyzed.The three-dimensional soil response in the frequency domain, ( ) ω , is obtained by the Direct Boundary Element Method (DBEM) following the work developed by Carrion et al. (2007).The three-dimensional version of the DBEM allows, in principle, the modeling of soil profiles with arbitrary geometry (Carrion et al., 2007).It is also assumed that the structure interacts with a rigid and massless foundation at the soil-foundation interface.The considered excitation are external forces, ext F , applied to the foundation degrees of freedom.
Latin American Journal of Solids and Structures, 2023, 20(6), e502 4/16 The system presented in Figure 1a is subdivided into two subsystems, respectively, the structure (Subsystem I) and the soil with the rigid and massless foundation (Subsystem II), as shown in Figure 1b.The interface forces between the structure and the soil-foundation systems are, respectively 1 I F and 2 I F .The responses of each subsystem are obtained according to the formulations described in the sequence.

Soil-foundation formulation
The three-dimensional Direct Boundary Element Method (DBEM) with constant quadrilateral elements is used to derive the soil response in the frequency domain, ( ) ω .The soil model is characterized by its density, S ρ , Young modulus, S E , Poisson ratio, S ν , and the internal damping factor, S η .A frequency independent, constant, damping factor, S η = cte, is considered for the soil in the present analysis (Beskos, 1987).The three-dimensional (3D) Boundary Element formulation in the frequency domain is based on the "frequency domain full space fundamental solution" (Dominguez, 1993).The BEM implementation described in (Carrion, 2002) allows for the dynamic response of 3D soils under applied surface tractions and also for the interaction of the soil with rigid surface or embedded foundations.Different soil profiles can be modeled, such as the homogeneous half-space, layer over a horizontal bedrock or layer on a non-horizontal bedrock, among others (Carrion et al., 2007).In the next paragraphs, the formulation of a rigid and massless foundation interacting with a 3D soil layer supported by a rigid bedrock is described.The formulation for rigid foundations interacting with the half-space can be obtained by a simplification of the case shown in this article (Carrion, 2002).
Figure 2 shows a rigid and massless foundation with dimensions 2 2 a a × bonded at the surface of a horizontal soil layer with thickness s h and supported by a rigid bedrock.The boundary at the soil-foundation interface is f Γ , the soil-free surface is 1 s Γ and the interface soil-rigid bedrock is 2 s Γ .At this first step is it convenient to consider that Γ =Γ ∪ Γ , which represents the soil boundaries that are not interacting with the rigid foundation interface.The Boundary Element equations can be formulated in terms of the surface displacements, U , and surface tractions, T .Assuming that the surface displacements, U , and tractions can also be subdivided according to ( , )   S S S U T Γ and ( , ) , the discretized BEM equations in matrix form can be written as (Carrion, 2002) (1) The displacements at the soil-foundation interface, { } f U , can be related to the rigid body degrees of freedom of the massless and rigid foundation, The tractions of the soil-foundation interface,{ } f T , are related to the vector of the resulting external forces, , applied to the rigid and massless foundation through a matrix of equilibrium equations,[ ]  and Structures, 2023, 20(6), e502 5/16 The problem can be rewritten in terms of the rigid body displacements of the rigid foundation, 0 { } U , and the applied external forces, { } F (Carrion, 2002) In equation ( 4), if the external forces are known, it is possible to determine the foundation rigid body displacements, 0 U , and the remaining soil-surface displacements, S U , the tractions at the soil-foundation interface, f T , and the tractions at the remaining soil surface, S T .Assuming the tractions, 1 { } s T , on the free surface, 1 s Γ , and the displacements, , on the soil bedrock foundation, 2 s Γ , are known, equation ( 4) can be rearranged as (5) The boundary conditions for the tractions on the free surface of the soil, 1 s Γ , are assumed to be zero, 1 { } 0 s T = .The soil is considered completely bonded at the soil-bedrock surface, 2 s Γ , and the corresponding boundary condition is Considering these boundary conditions, equation ( 5) can be written as This system of equations ( 6) can be solved to find the components of the rigid body displacements, 0 { } U , for a given frequency, ( ) ω .To obtain the rigid body displacements, the vector of the external forces,{ } F , assumes a sequence of unit value loads, i.e., { } { 1 0 0 0 0 0} and so on.The columns of the displacements resulting from the 6-unit value external force vectors may be organized in matrix form, to generate a dynamic flexibility (compliance) matrix, [ ( )] N ω , of the rigid and massless foundation interacting with the soil (Carrion, 2002).This frequency dependent flexibility matrix relates the vector of the external forces applied to the rigid foundation,{ } F , to the rigid body displacements of the foundation, 0 In equation ( 7) a is half the length of the foundation's side and S G is the shear modulus of the soil.Figures 3a and  3b show the real and imaginary parts of the flexibility functions, ( ) ω , relating the external vertical excitation, z F , to the vertical rigid foundation degree of freedom, z u , for the case of a homogeneous half-space (Figure 3a) and a horizontal layer over bedrock (Figure 3b).The layer depth is  (Carrion, 2002).
These figures show a large difference between the dynamic flexibility of a rigid and massless foundation interacting with the homogeneous half-space (Figure 3a) and with the horizontal layer over a bedrock (Figure 3b).The fact that the horizontal layer has a finite height, s h , introduces natural frequencies in the vertical direction.This point will be addressed again in the formulations and results that will follow.

Equations of motion for the Structure in relative coordinates
Consider the vector of the total displacements of the structural system, , and the vertical displacement of the rigid foundation, ( ) , respectively, can be written as where, is the vector of external forces applied to the rigid and massless foundation.Equation ( 8) can be rearranged in terms of the relative structural displacement,{ In order to obtain the stationary, frequency response of Subsystem I, equation ( 9) must be transformed from the time domain to the frequency domain with the application of the Fourier Transform leading to Displacements in the frequency domain, ( ) ω , are described in capital letters.In the present study, the damping matrix is assumed to be proportional to the mass and stiffness matrix , with the proportional coefficients µ and β (Caughey, 1960;Chopra, 2012).Next, a classical modal analysis is performed on the structural system described in equation ( 10).The solution of the linear eigenvalue problem will deliver the structural eigenvalues (natural frequencies), i ω , (i=1, n), the modal damping coefficients, i ξ , (i=1, n) and the eigenmodes (vibration modes),{ } i φ , (i=1, n) for the original structure in terms of the relative displacements.The eigenvectors can be organized into the modal matrix,[ ] Φ (Chopra, 2012;Fu and He, 2001).The modal matrix allows the transformation of the vector with structural physical degrees of freedom expressed in relative coordinates, ω , into the vector of modal coordinates, also in relative coordinates,{ By substituting equation ( 11) in equation ( 10) and pre-multiplying the resulting expression by the transposed of the modal matrix, [ ] T Φ , leads to Latin American Journal of Solids and Structures, 2023, 20(6), e502 7/16 Consider the orthogonal conditions for classical modal analysis of linear systems related to the modal matrix, [ ] Φ (Fu and He, 2001) In which [ ] I is the identity matrix, n ω and n ξ represent, respectively, the system n-th natural frequency and modal damping factor.Applying the transformations (13) to equations ( 12) results in the uncoupled set of equations in terms of the modal response of the structure subsystem, in relative coordinates, and in the frequency domain where the transfer function, [ ( )] H ω , and the generalized modal load coefficient, { ( )} ω Γ , are defined, respectively, by (Ferraz, 2021;Wu and Smith, 1995) The relative physical displacement of the structure, { ( )} relz U ω , can be recovered using the equation ( 11) Equations ( 17) express the relative displacement of the structure degrees of freedom,{ ω , as a function of the external excitation, { ( )} ext F ω , and the rigid foundation displacement, ( ) bz U ω .The soil influence is not yet incorporated into the structural response.One important characteristic of this formulation is that it is possible to choose an arbitrary number of structural modes in equation ( 15) to represent the dynamics of the structure (Ferraz, 2021;Wu and Smith, 1995).For structures with a large number of degrees of freedom and modes, the formulation allows investigating how many structural modes need to be considered to obtain a structural response within a given accuracy.

Coupling soil and structures responses
The objective of this session is to connect the soil and structural response formulations in the frequency domain, which will lead to the synthesis of Frequency Response Functions (FRFs) of the coupled system.
After defining the dynamic response of the soil, equation ( 7), and of the structural subsystem, equation ( 17), the coupling between the subsystems is carried out by the balance of forces and kinematic compatibility at the soil-structure interface, f Γ .First, a balance of forces in Subsystem I (structure) is performed using the total structural degrees of freedom, { ( )} z u t , in the time domain In equation ( 18), 1 ( ) F t is the interface force between the rigid foundation and the 1-st structural degree of freedom (see Figure 1b).Rewriting equation ( 18) in matrix form will lead to The equilibrium equation ( 19) can also be transformed to the frequency domain by the Fourier Transform Now consider, in equation ( 7), only the vertical soil-foundation compliance, ( ) ω relating the vertical displacement of the rigid and massless foundation interacting with the soil, ( ) bz U ω to the vertical foundation-structure interface force, 2 ( ) I F ω (see Figure 1b) The compliance, ( ) , can be inverted to produce a vertical dynamic stiffness of the soil-foundation, ( ) Considering equations ( 21) and ( 22), the displacement of the rigid and massless foundation interacting with the soil in the frequency domain, ( ) bz U ω , can be related to the external force applied to the foundation-structure interface, 2 ( ) Equilibrium conditions at the rigid foundation-structure interface prescribe that 1 2 ( ) ( ) 0 Equations ( 20), ( 23) and ( 24) can be rearranged resulting in On the other hand, the total response of the structure degrees of freedom, { ( )} z U ω , is given by the sum of the relative displacements, { ( )} relz U ω , and the displacement from the soil, ( Reorganizing equations ( 17), ( 25) and ( 26) will result in Collecting the terms where, H ω , as can be seen by analyzing equation ( 29).The upper index 'dssi' has been added to the displacement to indicate that it already contains the soil-structure effect.The excitation remains the external forces applied directly to the n foundation degrees of freedom, { ( )} ext F ω .
To determine the total structure response with soil influence, { ( )} given by equation ( 17) where, ( ) H ω .Equation ( 30) allows the construction of modified frequency response functions (FRFs) for the structure, using an arbitrary number of structural modes, and already considering the soil-structure interaction effects.As will be explained in the next sessions, structural modal parameters may be extracted from these modified FRFs and used to reconstruct a system of orthogonal equations of motion in the time domain in terms of modal parameters.

Modal parameters extraction from modified FRFs
The Rational Fraction Polynomial Method (RFPM) was chosen to extract the modal parameters (natural frequencies, damping factors and modal forms) from the FRFs containing soil-structure interaction soil effects given in equation ( 30).The RFPM is a classic extraction method presented in Ewins (2000), and a more detailed formulation can be seen in Richardson andFormenti (1982) andFerraz (2021).This method is based on the approximation of the FRF curve by means of complex orthogonal polynomials, represented in equation (32).From this approximation, the relationships between the poles and residuals of the partial fractions and the modal quantities are determined.
In equation ( 32) N represents the number of modes in the FRF to be adjusted.The modal quantities extracted from the formulation (32) are eigenfrequencies including the soil-foundation effects, dssi i ω , as well as the modal damping coefficients, dssi i ξ , and eigenvectors, { } dssi i φ .Results of the extracted modified modal quantities using this methodology will be presented in section 3.

Transient response by modal superposition
Starting from the modal quantities, represented by the natural frequencies, dssi , a set of orthogonal, uncoupled, equations of motion in the time domain, in terms of the modal coordinates, { ( )} dssi z q t , may be constructed to represent the dynamic response of the system under analysis (Fu and He, 2001;Ferraz, 2021) Latin American Journal of Solids and Structures, 2023, 20(6), e502 10/16 Equation ( 33) represents a set of uncoupled ordinary differential equations that can be directly integrated in time domain, delivering the response in terms of the modal coordinates including the soil-structure effects, { ( )} dssi z q t .The displacement response in physical coordinates, { ( )} dssi z u t , can be obtained through equation ( 11) In the present work, the Newmark method (Chopra, 2012) was chosen to perform the numerical integration of the equations of motion (33), due to its unconditional stability for any time step when the parameters  = 0.25 and  = 0.5 are used (Chopra, 2012).

Expanded modal basis
In soil models with unbounded dimensions, like the half-space, there no natural frequencies or normal vibration modes.On the other hand, for soil profiles with a limited dimension, such as a layer of finite depth over a rigid stratum or bedrock, there are natural frequencies and normal vibration modes along the finite dimension.For soil profiles that present eigenfrequencies, it might be necessary to expand the original modal basis of the structure to include the new eigenfrequencies that may appear on the FRFs of the coupled soil-structure system due to soil response.
This influence can be seen in the example of Figure 4, which simulates a soil-structure system in which the soil is represented by a mass-spring-damper system with 2 dofs and the structure with 3 dofs.Note that the FRFs of each separate subsystem only show resonances referring to their number of dofs.However, when the coupling is performed, 5 resonances can be observed in the FRF of the coupled system.So, frequency response functions obtained by equations (30) may present more resonances than the structure's N original degrees of freedom, depending on the characteristics of the soil profile supporting the foundation and structure.For these cases, it is necessary to expand the modal basis to include a larger number of system degrees of freedom S , which is the sum of the original N structural dofs and the new degrees of freedom added to the system by the inclusion of the soil response.The equation transforming the expanded modal displacements, 1 { ( )} d s si z s q t × , with S degrees of freedom is given by The expanded modal basis is given by [ ] dssi n s ×

Φ
. The orthogonality conditions (13) must also be changed to reflect the expansion of the modal basis It should be noted that this expanded system has now S eigenfrequencies and S eigenmodes.Based on equations (36) an updated set of orthogonal equations of motion in the time domain can be synthesized The formulation resulting in equation ( 37) indicates that a new, expanded, system of equations of motion with S degrees of freedom is synthesized to cope with the modifications in the time domain structural response to include the soil effect.

Systems and input data
For the numerical analysis, two distinct models are investigated.These are depicted in Figures 5a and 5b.The structure is composed of bar elements with 9 dofs and is subjected to an excitation force applied to the 9-th dof.The soil models are the homogeneous half-space (Figure 5a) and a horizontal layer of depth, = , with a being the foundation half-width (Figure 5b).The soil properties are given in Table 1.The soil dynamic flexibility,    For comparison purposes, the classic modal analysis was carried out only for the structure, considering a fixed base, obtaining the natural frequencies and damping factors for the 9 modes, as shown in Table 3.
Table 3 Modal parameters of the structure (fixed base).

Modified FRFs for the structure interacting with the soil models
Using the methodology summarized in equation ( 30), a set of 9 FRFs for the structure interacting with both soil profiles is shown in Figure 6.The black dashed line corresponds to the case of the homogeneous half-space and the red line correspond to the layer over a rigid stratum.

Extraction of modal data from both models
The Rational Fraction Polynomial Method (RFPM) was applied to this set of 9 FRFs to extract eigenfrequencies, i ω , and modal damping coefficients, i ξ .Table 4 contains the values of the 9 natural frequencies and damping factors extracted for the case of the half-space.Table 5 shows the natural frequencies and damping factors for the case of the horizontal layer.It should be noted that in this second case, 12 eigenfrequencies were extracted from the FRFs by the RFP Method.As expected, the eigenfrequencies of the soil layer do influence the dynamic behavior of the soil-structure system that needs an expanded modal basis to describe its dynamics.
A comparison among the modal data from Table 3 (structure on a rigid base), Table 4 (structure on the half-space) and Table 5 (structure on a layer over bedrock) reveals that the presence of the soil introduces higher damping coefficients (Tables 4 and 5) than the structure on a fixed base (Table 3).This is consistent with the additional damping introduced in the system by the geometric damping mechanism.

Transient responses of systems with soil influence
To determine the transient responses of the structures, the previously extracted modal parameters are used as input data for the motion equations (37).The excitation force was applied at the last dof of the structure (i=9) and its time dependency is given by: The motion equations were numerically integrated by the Newmark method, using the parameters given in Table 6. Figure 7 shows the transient responses of both systems for each dof of the structure.The validation of the present approach for the soil as a half-space model was presented in Ferraz et al. (2021b).The extension to the model of a layer over bedrock, presented in the current article is a new result.Nevertheless, the response for the structure supported by a horizontal layer over bedrock is consistent with the expected results.The amount of geometric damping of the layer is notably smaller than the one presented by the homogeneous half-space, due to the wave reflections at the rigid base.The transient response for the structure interacting with the layer shown in Figure 7 for all 9 structural dofs present a much smaller overall damping.The layer over bedrock is also slightly more rigid than the half-space resulting in a small increase in the frequency of the structure response.

CONCLUSION
This article presented a methodology to obtain the transient response of linear structures interaction with soil profiles.The soil response in the frequency domain was synthesized by the Boundary Element Method.A rigid and massless foundation interacting with the homogenous half-space and with a horizontal layer over a rigid stratum were the considered soil models.
Latin American Journal of Solids and Structures, 2023, 20(6), e502 15/16 The dynamic response of the structure was based on a modal description of the system, which allows to consider an arbitrary number of degrees of freedom.After the coupling of the structure with the soil, in the frequency domain, a set of Frequency Response Functions (FRFs) for the structure were synthesized in which the effects of the soil-structure interaction were already accounted for.From these modified FRFs, modal parameters were extracted in order to build a set of uncoupled differential equations that describe the coupled transient behavior of the structure.An expansion of the structural modal base was presented to properly describe the dynamics of soil-structures in which the soil model presents eigenfrequencies.The layer over a bedrock represents such a system.The results obtained for the transient response of the structures were consistent with the expected behavior of the considered cases.The methodology can be applied to arbitrary soil profiles, provided the frequency response of the soil model is available.

Figure 1
Figure 1 Soil-structure system model and system subdivision.

Figure 2
Figure 2 Rigid and massless foundation at the surface of a 3D horizontal layer over bedrock

Figure 3
Figure 3 Comparison flexibilities N uzFz (ω) for the half-space (a) and for a layer over bedrock (b).

.
bz u t .A vector of the structural displacements, { } r el z u , relative to the rigid foundation vertical displacement vertical can be defined as: Using the relative displacements, the time domain equations of motion of the structure with mass matrix, [ ] M , damping and stiffness matrices [ ] C and [ ] K by expression (28) must be added to the vector of the structure relative displacements, compliance matrix relating the total displacement of the structure degrees of freedom, { ( )} dssi z U ω , under an external excitation, { ( )} ext F ω , applied at the structure degrees of freedom and already considering the soil effects, ( ) s K ω , for an arbitrary number of modes to describe the structure dynamics, [ ( )] -th mode and the matrix of modal shapes, [ ] dssi Φ

Figure 4
Figure 4 Soil influence on structure FRFs.
System I and II is  shown, respectively in Figures (3a) and (3b).Table2presents the properties of the structure, including the mass and stiffness values for each element, along with the coefficients of proportional damping µ and β .

Figure 5
Figure 5 Soil-structure systems with different soil profiles.

Figure 6
Figure 6 Structure FRFs with soil influence (System I x System II).

Figure 7
Figure 7 Transient responses (System I x System II).
Coupling Modal Analysis with the BEM for the Transient Response of Bar Structures Interacting with Three-Dimensional Soil Profiles Amauri Coelho Ferraz et al.
Latin American Journal of Solids Coupling Modal Analysis with the BEM for the Transient Response of Bar Structures Interacting with Three-Dimensional Soil Profiles Amauri Coelho Ferraz et al.Equation (28)furnishes the vertical displacement of the rigid and massless foundation interacting with the soil, , and the dynamics of the structure described in terms of modal quantities, [ ( )] Latin American Journal of Solids and Structures, 2023, 20(6), e502 9/16 s K ω

Table 1
Soil properties.

Table 2
Structure properties.

Table 4
Modal parameters of System I.

Table 5
Modal parameters of System II.
Coupling Modal Analysis with the BEM for the Transient Response of Bar Structures Interacting with Three-Dimensional Soil Profiles Amauri Coelho Ferraz et al.Latin American Journal of Solids and Structures, 2023, 20(6), e502 14/16