• Proceedings of the Earthquake Engineering Society of Korea Conference
• /
• 2000.10a
• /
• pp.81-88
• /
• 2000

The elastic wave equation is solved using the finite-difference method in 3D space to simulate the seismic wave propagation. It is based on the velocity-stress formulation of the equation of motion on a staggered grid. The nonreflecting boundary conditions are used to attenuate the wave field close to the numerical boundary. To satisfy the stress-free conditions at the free-surface boundary, a new formulation combining the zero-stress formalism with the vacuum one is applied. The effective media parameters are employed to satisfy the traction continuity condition across the media interface. With use of the moment-tensor components, the wide range of source mechanism parameters can be specified. The numerical experiments are carried out in order to test the applicability and accuracy of this scheme and to understand the fundamental features of the wave propagation under the generalized elastic media structure. Computational results show that the scheme is sufficiently accurate for modeling wave propagation in 3D elastic media and generates all the possible phases appropriately in under the given heterogeneous velocity structure. Also the characteristics of the ground motion in an sedimentary basin such as the amplification, trapping, and focusing of the elastic wave energy are well represented. These results demonstrate the use of this simulation method will be helpful for modeling the ground motion of seismological and engineering purpose like earthquake hazard assessment, seismic design, city planning, and etc..

• Journal of Ocean Engineering and Technology
• /
• v.21 no.6
• /
• pp.26-33
• /
• 2007

Finite element simulation of elastic wave propagation in a concrete plate was carried out to investigate its modeling and damage detection procedures. For the numerical stability three criteria were introduced and tested. With a proper element size and time increment, two different kinds of damage scenarios (crack and deterioration) were applied to verify the feasibility of the finite element simulation. It is shown that the severities of those damages are sensitive to the received displacement signals.

• Smart Structures and Systems
• /
• v.18 no.2
• /
• pp.335-354
• /
• 2016

This contribution presents an extended one-dimensional theory for piezoelectric beam-type structures with non-ideal electrodes. For these types of electrodes the equipotential area condition is not satisfied. The main motivation of our research is originated from passive vibration control: when an elastic structure is covered by several piezoelectric patches that are linked via resistances and inductances, vibrational energy is efficiently dissipated if the electric network is properly designed. Assuming infinitely small piezoelectric patches that are connected by an infinite number of electrical, in particular resistive and inductive elements, one obtains the Telegrapher's equation for the voltage across the piezoelectric transducer. Embedding this outcome into the framework of Bernoulli-Euler, the final equations are coupled to the wave equations for the longitudinal motion of a bar and to the partial differential equations for the lateral motion of the beam. We present results for the wave propagation of a longitudinal bar for several types of electrode properties. The frequency spectra are computed (phase angle, wave number, wave speed), which point out the effect of resistive and inductive electrodes on wave characteristics. Our results show that electrical damping due to the resistivity of the electrodes is different from internal (=strain velocity dependent) or external (=velocity dependent) mechanical damping. Finally, results are presented, when the structure is excited by a harmonic single force, yielding that resistive-inductive electrodes are suitable candidates for passive vibration control that might be of great interest for practical applications in the future.

• Earthquakes and Structures
• /
• v.10 no.1
• /
• pp.91-104
• /
• 2016

In this paper we investigate the existence of SH-waves in fiber-reinforced layer placed over a heterogeneous elastic half-space. The heterogeneity of the elastic half-space is caused by the exponential variations of density and rigidity. As a special case when both the layers are homogeneous, our derived equation is in agreement with the general equation of Love wave. Numerically, it is observed that the velocity of SH-waves decreases with the increase of heterogeneity and reinforced parameters. The dimensionless phase velocity of SH-waves increases with the decreases of dimensionless wave number and shown through figures.

• Structural Engineering and Mechanics
• /
• v.87 no.2
• /
• pp.165-172
• /
• 2023

This paper presents a novel method for modeling elastic wave propagation in long beams. The proposed method derives a solution for the transient transverse displacement of the beam's neutral axis without assuming the separation of variables (SV). By mapping the governing equation from the space domain to the frequency domain using Fourier transformation (FT), the transverse displacement function is determined as a convolution integral of external loading functions and a combination of trigonometric and Fresnel functions. This method determines the beam's response to general loading conditions as a linear combination of the analytical response of a beam subjected to an abrupt localized loading. The proposed solution method is verified through finite element analysis (FEA) and wave propagation patterns are derived for tone burst loading with specific frequency contents. The results demonstrate that the proposed solution method accurately models wave dispersion, reduces computational cost, and yields accurate results even for high-frequency loading.

• Geophysics and Geophysical Exploration
• /
• v.11 no.2
• /
• pp.93-98
• /
• 2008

Abstract: Surface topography has a significant influence on seismic wave propagation in a reflection seismic exploration. Effects of surface topography on two-dimensional elastic wave propagation are investigated through modeling using a weighted-averaging (WA) finite-element method (FEM), which is computationally more efficient than conventional FEM. Effects of air layer on wave propagation are also investigated using flat surface models with and without air. To validate our scheme in modeling including topography, we compare WA FEM results for irregular topographic models against those derived from conventional FEM using one set of rectangular elements. For the irregular surface topography models, elastic wave propagation is simulated to show that breaks in slope act as a new source for diffracted waves, and that Rayleigh waves are more seriously distorted by surface topography than P-waves.

• 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.

• Geophysics and Geophysical Exploration
• /
• v.6 no.2
• /
• pp.77-86
• /
• 2003

We designed a new time-domain, finite-difference, elastic wave modeling technique, based on a displacement formulation. which yields nearly correct solutions to Lamb's problem. Unlike the conventional, displacement-based, finite-difference method using a node-based grid set (where both displacements and material properties such as density and Lame constants are assigned to nodal points), in our new finite-difference method, we use a cell-based grid set (where displacements are still defined at nodal points but material properties within cells). In the case of using the cell-based grid set, stress-free conditions at the free surface are naturally described by the changes in the material properties without any additional free-surface boundary condition. Through numerical tests, we confirmed that the new second-order finite differences formulated in the cell-based grid let generate numerical solutions compatible with analytic solutions unlike the old second-order finite-differences formulated in the node-based grid set.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.35 no.8
• /
• pp.957-964
• /
• 2011

This investigation presents a finite element method to obtain the transmission properties of bulk elastic waves in piezoelectric band gap structures(phonon crystals) for varying frequencies and modes. To this end, periodic boundary conditions are imposed on a three-dimensional model while both in-plane and out-of-plane modes are included. In particular, the mode decoupling characteristics between in-plane and out-of-plane modes are identified for each electric poling direction and the results are incorporated in the finite element modeling. Through numerical simulations, the proposed modeling method was found to be a useful, effective one for analyzing the wave characteristics of various types of piezoelectric phononic band gap structures.

• Proceedings of the Korean Geotechical Society Conference
• /
• 2005.03a
• /
• pp.302-307
• /
• 2005

Love wave and Rayleigh wave are the major elastic waves belonging to the category of the surface wave. The fact that Love wave is not contaminated by P-wave which makes Love wave superior to Rayleigh wave and other body waves. Therefore, the information that Love wave carries is more distinct and clearer than the information of Rayleigh wave. Based on theoretical research, the joint inversion analysis which is used both Love wave dispersion information and Rayleigh wave dispersion information was proposed. Purpose of the joint inversion analysis is to improve accuracy and convergency of inversion results utilizing that frequency contribution of each wave is different. This analysis technique is consisted of the forward modeling using transfer matrix, the sensitivity matrix determined to the ground system and DLSS(Damped Least Square Solution) as a inversion technique. The application of this analysis was examined through the field test.