• Structural Engineering and Mechanics
• /
• v.50 no.3
• /
• pp.349-364
• /
• 2014

Perfectly matched layers are employed in time harmonic analysis of dam-foundation systems. The Lysmer boundary condition at the truncation boundary of the PML region has been incorporated in the formulation of the dam-foundation FE model (including PML). The PML medium is defined in a way that the formulation of the system can be transformed into time domain. Numerical experiments show that applying Lysmer boundary conditions at the truncation boundary of the PML area reduces the computational cost and make the PML approach a more efficient technique for the analysis of dam-foundation systems.

• Structural Engineering and Mechanics
• /
• v.51 no.4
• /
• pp.685-705
• /
• 2014

This paper presents a new formulation for forward scalar wave simulations in semi-infinite media. Perfectly-Matched-Layers (PMLs) are used as a wave absorbing boundary layer to surround a finite computational domain truncated from the semi-infinite domain. In this work, a hybrid formulation was developed for the simulation of scalar wave motion in two-dimensional PML-truncated domains. In this formulation, displacements and stresses are considered as unknowns in the PML domain, while only displacements are considered to be unknowns in the interior domain. This formulation reduces computational cost compared to fully-mixed formulations. To obtain governing wave equations in the PML region, complex coordinate stretching transformation was introduced to equilibrium, constitutive, and compatibility equations in the frequency domain. Then, equations were converted back to the time-domain using the inverse Fourier transform. The resulting equations are mixed (contain both displacements and stresses), and are coupled with the displacement-only equation in the regular domain. The Newmark method was used for the time integration of the semi-discrete equations.

• Earthquakes and Structures
• /
• v.3 no.3_4
• /
• pp.563-580
• /
• 2012

We modify the formulation of a recently developed absorbing boundary condition (ABC), the perfectly matched discrete layers (PMDL), to incorporate the excitation coming from the exterior such as earthquake waves. The modified formulation indicates that the effect of the exterior excitation can be incorporated into PMDL ABCs (traditionally designed to treat only interior excitation) simply by applying appropriate forces on the nodes connected to the first PMDL layer. Numerical results are presented to clearly illustrate the effectiveness of the proposed method.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.28 no.2
• /
• pp.123-130
• /
• 2015

This paper presents a new formulation for transient simulations of microwave propagation in heterogeneous unbounded domains. In particular, perfectly-matched-layers(PMLs) are introduced to allow for wave absorption at artificial boundaries used to truncate the infinite extent of the physical domains. The development of the electromagnetic PML targets the application to engineering mechanics problems such as structural health monitoring and inverse medium problems. To formulate the PML for plane electromagnetic waves, a complex coordinate transformation is introduced to Maxwell's equations in the frequency-domain. Then the PML-endowed partial differential equations(PDEs) for transient electromagnetic waves are recovered by the application of the inverse Fourier transform to the frequency-domain equations. A mixed finite element method is employed to solve the time-domain PDEs for electric and magnetic fields in the PML-truncated domain. Numerical results are presented for plane microwaves propagating through concrete structures, and the accuracy of solutions is investigated by a series of error analyses.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.26 no.4
• /
• pp.223-231
• /
• 2013

This paper presents a time-domain Gauss-Newton full-waveform inversion method for the material profile reconstruction in heterogeneous semi-infinite solid media. To implement the inverse problem in a finite computational domain, perfectly-matchedlayers( PMLs) are introduced as wave-absorbing boundaries within which the domain's wave velocity profile is to be reconstructed. The inverse problem is formulated in a partial-differential-equations(PDE)-constrained optimization framework, where a least-squares misfit between measured and calculated surface responses is minimized under the constraint of PML-endowed wave equations. A Gauss-Newton-Krylov optimization algorithm is utilized to iteratively update the unknown wave velocity profile with the aid of a specialized regularization scheme. Through a series of one-dimensional examples, the solution of the Gauss-Newton inversion was close enough to the target profile, and showed superior convergence behavior with reduced wall-clock time of implementation compared to a conventional inversion using Fletcher-Reeves optimization algorithm.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.26 no.4
• /
• pp.213-221
• /
• 2013

This paper introduces a mesh continuation scheme for a one-dimensional inverse medium problem to reconstruct the spatial distribution of elastic wave velocities in heterogeneous semi-infinite solid domains. To formulate the inverse problem, perfectly-matched-layers(PMLs) are introduced as wave-absorbing boundaries that surround the finite computational domain truncated from the originally semi-infinite extent. To tackle the inverse problem in the PML-truncated domain, a partial-differential-equations(PDE)-constrained optimization approach is utilized, where a least-squares misfit between calculated and measured surface responses is minimized under the constraint of PML-endowed wave equations. The optimization problem iteratively solves for the unknown wave velocities with their updates calculated by Fletcher-Reeves conjugate gradient algorithms. The optimization is performed using a mesh continuation scheme through which the wave velocity profile is reconstructed in successively denser mesh conditions. Numerical results showed the robust performance of the mesh continuation scheme in reconstructing target wave velocity profile in a layered heterogeneous solid domain.

• Geophysics and Geophysical Exploration
• /
• v.15 no.3
• /
• pp.121-128
• /
• 2012

Absorbing boundary conditions are used to mitigate undesired reflections that can arise at the model's truncation boundaries. We apply a complex frequency shifted perfectly matched layer (CFS-PML) to elastic wave modeling in the frequency domain. Modeling results show that the performance of our implementation is superior to other absorbing boundaries. We consider the coefficients of CFS-PML to be optimal when the kinetic energy becomes to the minimum, and propose the modified CFS-PML that has the CFS-PML coefficient ${\alpha}_{max}$ defined as a function of frequency. Results with CFS-PML and modified CFS-PML are significantly improved compared with those of the classical PML technique suffering from large spurious reflections at grazing incidence.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.34 no.4
• /
• pp.221-230
• /
• 2021

This paper proposes full waveform inversion (FWI) for estimating the physical properties of a layered half-space. An FWI solution is obtained using a genetic algorithm (GA), which is a well-known global optimization approach. The dynamic responses of a layered half-space subjected to a harmonic vertical disk load are measured and compared with those calculated using the estimated physical properties. The responses are calculated using the thin-layer method, which is accurate and efficient for layered media. Subsequently, a numerical model is constructed for a layered half-space using mid-point integrated finite elements and perfectly matched discrete layers. An objective function of the global optimization problem is defined as the L2-norm of the difference between the observed and estimated responses. A GA is used to minimize the objective function and obtain a solution for the FWI. The accuracy of the proposed approach is applied to various problems involving layered half-spaces. The results verify that the proposed FWI based on a GA is suitable for estimating the material properties of a layered half-space, even when the measured responses include measurement noise.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.37 no.2
• /
• pp.133-141
• /
• 2024

Considering fluid-structure-soil interactions, a finite-element model for a liquid-storage tank is presented and the nonlinear earthquake response analysis is formulated. The tank structure is modeled considering shell elements with geometric and material nonlinearities. The fluid is represented by acoustic elements and combined with the structure using interface elements. To consider the soil-structure interactions, the near- and far-field regions of soil are modeled with solid elements and perfectly matched discrete layers, respectively. This approach is applied to the seismic fragility analysis of a 200,000 kL liquid-storage tank. The fragility curve is observed to be influenced by the amplification and filtering of rock outcrop motions at the site when the soil-structure interactions are considered.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.29 no.4
• /
• pp.317-325
• /
• 2016

In this study, nonlinear earthquake responses of a soil-structure interaction(SSI) system which is subjected to a three-directional ground motion are examined. The structure and the near-field region of soil, where the geometry is irregular, the material properties are heterogeneous, and nonlinear dynamic responses are expected, are modeled by nonlinear finite elements. On the other hand, the infinite far-field region of soil, which has a regular geometry and homogeneous material properties and dynamic responses is assumed linearly elastic, is represented by three-dimensional perfectly matched discrete layers which can radiate elastic waves into infinity efficiently. Nonlinear earthquake responses of the system subjected to a three-directional ground motion are calculated with the numerical model. It is observed that the dynamic responses of a SSI system to a three-directional motion have a predominant direction according to the characteristics of the ground motion. The responses must be evaluated using precise analysis methods which can consider nonlinear behaviors of the system accurately. The the method employed in this study can be applied easily to boundary nonlinear problems as well as material nonlinear problems.