• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1995.04a
• /
• pp.28-34
• /
• 1995

Behavior of underground structural systems is usually complicated because of various unknown parameters. In order to construct those structural systems safely and economically, exact identification of the system parameters and accurate analysis of the system behaviors are essentially required. In this study, a forward analysis program, which is able to eliminate numerical errors due to far field boundary effect, is developed by coupling finite and boundary elements. In this coupled analysis, boundary elements are used in the semi-infinite domain where stress variation is small, and finite elements in the stress concentration region where material nonlinearity should be considered. Then, a back analysis program which can identify the system parameters is developed using the direct method to be combined with the forward analysis program. The elastic modulus and initial stress, which are most important in the description of the behavior of underground structures, are taken as the system parameters. A simple example is examined 0 show that the method can be used effectively.

• Journal of Mechanical Science and Technology
• /
• v.14 no.3
• /
• pp.331-337
• /
• 2000

A numerical procedure for contact analysis and calculating subsurface stress was developed. The procedure takes the advantage of signal processing technique in frequency domain to achieve shorter computing time. Boussinesq's equation was adopted as a response function in contact analysis. The validity of this procedure was proved by comparing the numerical results with the exact solutions. The fastness of this procedure was also compared with other algorithm.

• Proceedings of the Earthquake Engineering Society of Korea Conference
• /
• 2002.09a
• /
• pp.123-130
• /
• 2002

The ground-borne vibration from pile driving is causing many public discontents. However, because of the fact that the characteristics of wave propagation and attenuation are not well understood, systematic and effective vibration reduction measures can not be taken. This paper attempts to study the propagation of stress waves induced by the pile driving. To simulate the wave propagation in a semi-infinite domain, the so-called absorbing boundaries are incorporated in the finite element method and a series of numerical simulations is performed. Numerical results show that the surface displacement and velocity increase first and then decrease as the pile penetration depth becomes larges.

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

• Tribology and Lubricants
• /
• v.38 no.3
• /
• pp.73-83
• /
• 2022

This paper is concerned with an analysis of a surface edge crack emanated from a sharp contact edge. For a geometrical model, a square wedge is in contact with a half plane whose materials are identical, and a surface perpendicular crack initiated from the contact edge exists in the half plane. To analyze this crack problem, it is necessary to evaluate the stress field on the crack line which are induced by the contact tractions and pseudo-dislocations that simulate the crack, using the Bueckner principle. In this Part I, the stress filed in the half plane due to the contact is re-summarized using an asymptotic analysis method, which has been published before by the author. Further focus is given to the stress field in the half plane due to a pseudo-edge dislocation, which will provide a stress solution due to a crack (i.e. a continuous distribution of edge dislocations) later, using the Burgers vector. Essential result of the present work is the corrective functions which modify the stress field of an infinite domain to apply for the present one which has free surfaces, and thus the infiniteness is no longer preserved. Numerical methods and coordinate normalization are used, which was developed for an edge crack problem, using the Gauss-Jacobi integration formula. The convergence of the corrective functions are investigated here. Features of the corrective functions and their application to a crack problem will be given in Part II.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.15 no.5
• /
• pp.499-508
• /
• 2004

An efficient method of moments(MoM), which can lead to the analytical evaluation of the matrix elements, is proposed to analyze microstrip structures. The present method is formulated in conjunction with use of a new closed-form spatial-domain Green's functions which are derived by use of the integral formula for semi-infinite integrals of Bessel functions. It is observed that the computational efficiency such as the amount of calculation and computation speed has been improved due to the present MoM scheme by a factor of about 4 in comparison with the previous method. To validate the proposed method, several numerical examples are presented.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.9 no.1
• /
• pp.43-54
• /
• 1997

This paper presents an approximate analytical solution to a two-region one-dimensional model for the charging process of stratified thermal storage tanks with variable inlet temperature in the presence of momentum-induced mixing. Based on the superposition principle, an arbitrary-varying inlet temperature is decomposed into inherent discontinuous steps and continuous intervals approximated as a finite number of piecewise linear functions. This approximation allows the temperature of the upper perfectly-mixed layer to be expressed in terms of constant, linear and exponential functions with respect to time. Applying the Laplace transform technique to the model equation for the lower thermocline layer subject to each of three representative interfacial conditions yields compact-form solutions, a linear combination of which constitutes the final temperature profile. A systematic method for deriving solutions to the plug-flow problem having polynomial-type boundary conditions is also established. The effect of adiabatic exit boundary on solution behaviors proves to be negligible under the actual working conditions, which justifies the assumption of semi-infinite domain introduced in the solution procedure. Finally, the approximate solution is validated by comparing it with an exact solution obtained for a specific variation of inlet temperature. Excellent agreements between them suffice to show the necessity and utility of this work.

• Geophysics and Geophysical Exploration
• /
• v.11 no.2
• /
• pp.153-160
• /
• 2008

Seismic modeling is used to simulate wave propagation in the earth. Although the earth's subsurface is usually semi-infinite, we cannot handle the semi-infinite model in seismic modeling because of limited computational resources. For this reason, we usually assume a finite-sized model in seismic modeling. In that case, we need to eliminate the edge reflections arising from the artificial boundaries introducing a proper boundary condition. In this study, we changed three kinds of boundary conditions (sponge boundary condition, Clayton and Engquist's absorbing boundary condition, and Higdon's transparent boundary condition) so that they can be applied in elastic wave modeling for anisotropic media. We then apply them to several models whose Poisson's ratios are different. Clayton and Engquist's absorbing boundary condition is unstable in both isotropic and anisotropic media, when Poisson's ratio is large. This indicates that the absorbing boundary condition can be applied in anisotropic media restrictively. Although the sponge boundary condition yields good results for both isotropic and anisotropic media, it requires too much computational memory and time. On the other hand, Higdon's transparent boundary condition is not only inexpensive, but also reduce reflections over a wide range of incident angles. We think that Higdon's transparent boundary condition can be a method of choice for anisotropic media, where Poisson's ratio is large.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.14 no.3
• /
• pp.11-22
• /
• 2010

An analysis method is proposed for the transient linear or nonlinear analysis of dynamic interactions between a flexible dam body and reservoir impounding compressible water under earthquake loadings. The coupled dam-reservoir system consists of three substructures: (1) a dam body with linear or nonlinear behavior; (2) a semi-infinite fluid region with constant depth; and (3) an irregular fluid region between the dam body and far field. The dam body is modeled with linear and/or nonlinear finite elements. The far field is formulated as a displacement-based transmitting boundary in the frequency domain that can radiate energy into infinity. Then the transmitting boundary is transformed for the direct coupling in the time domain. The near field region is modeled as a compressible fluid contained between two substructures. The developed method is verified and applied to various earthquake response analyses of dam-reservoir systems. Also, the method is applied to a nonlinear analysis of a concrete gravity dam. The results show the location and severity of damage demonstrating the applicability to the seismic evaluation of existing and new dams.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.25 no.6
• /
• pp.549-558
• /
• 2012

A topology optimization method for phononic crystals is developed for the design of sound barriers, using the level set approach. Given a frequency and an incident wave to the phononic crystals, an optimal shape of periodic inclusions is found by minimizing the norm of transmittance. In a sound field including scattering bodies, an acoustic wave can be refracted on the obstacle boundaries, which enables to control acoustic performance by taking the shape of inclusions as the design variables. In this research, we consider a layered structure which is composed of inclusions arranged periodically in horizontal direction while finite inclusions are distributed in vertical direction. Due to the periodicity of inclusions, a unit cell can be considered to analyze the wave propagation together with proper boundary conditions which are imposed on the left and right edges of the unit cell using the Bloch theorem. The boundary conditions for the lower and the upper boundaries of unit cell are described by impedance matrices, which represent the transmission of waves between the layered structure and the semi-infinite external media. A level set method is employed to describe the topology and the shape of inclusions. In the level set method, the initial domain is kept fixed and its boundary is represented by an implicit moving boundary embedded in the level set function, which facilitates to handle complicated topological shape changes. Through several numerical examples, the applicability of the proposed method is demonstrated.