• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1992.04a
• /
• pp.58-64
• /
• 1992

It is usual that underground structures are constructed within multi-layered medium. In this paper, an efficient numerical model ling of multi-layered structural systems is studied using coupled analysis of finite elements and boundary elements. The finite elements are applied to the area in which the material nonlinearity is dominated, and the boundary elements are applied to the far field area where the nonlinearity is relatively weak. In the boundary element model 1 ins of the multi-layered medium, fundamental solutions are restricted. Thus, methods which can utilize existing Kelvin and Melan solution are sought for the interior multi-layered domain problem and semi infinite domain problem. Interior domain problem which has piecewise homogeneous layers is analyzed using boundary elements with Kelvin solution; by discretizing each homogeneous subregion and applying compatibility and equilibrium conditions between interfaces. Semi-infinite domain problem is analyzed using boundary elements with Melan solution, by superposing unit stiffness matrices which are obtained for each layer by enemy method. Each methodology is verified by comparing its results which the results from the finite element analysis and it is concluded that coupled analysis using boundary elements and finite elements can be reasonable and efficient if the superposition technique is applied for the multi-layered semi-infinite domain problems.

• The Journal of the Korea institute of electronic communication sciences
• /
• v.19 no.1
• /
• pp.31-38
• /
• 2024

The digital phase-locked loop(DPLL) is one of the circuits composed of a digital detector, digital loop filter, voltage-controlled oscillator, and divider as a fundamental circuit, widely used in many fields such as electrical and circuit fields. A state estimator using various mathematical algorithms is used to improve the performance of a digital phase-locked loop. Traditional state estimators have utilized Kalman filters of infinite impulse response state estimators, and digital phase-locked loops based on infinite impulse response state estimators can cause rapid performance degradation in unexpected situations such as inaccuracies in initial values, model errors, and various disturbances. In this paper, we propose a two-layer Frobenius norm-based finite impulse state estimator to design a new digital phase-locked loop. The proposed state estimator uses the estimated state of the first layer to estimate the state of the first layer with the accumulated measurement value. To verify the robust performance of the new finite impulse response state estimator-based digital phase locked-loop, simulations were performed by comparing it with the infinite impulse response state estimator in situations where noise covariance information was inaccurate.

• The Journal of the Acoustical Society of Korea
• /
• v.17 no.3E
• /
• pp.35-42
• /
• 1998

Transient acoustic and elastic wave propagation in inhomogeneous media are studied by using the Fourier method. It is known that the fourier method has advantages in memory requirements and computing speed over conventional methods such as FDM and FEM, because the Fourier method needs less grid points for achieving the same accuracy. To verify the proposed numerical scheme, several examples having analytic solutions are considered, where two different semi-infinite media are in contact along a plane boundary. The comparisons of numerical results by the Fourier method and analytic solutions show good agreements. In addition, the fourier method is applied to a layered half-plane, in which an elastic semi-infinite medium is covered by an elastic layer of finite thickness. It is showed how to derive the analytic solutions by using the Cagniard-de Hoop method. The numerical solutions are in excellent agreements with analytic results.

• Structural Engineering and Mechanics
• /
• v.26 no.5
• /
• pp.591-615
• /
• 2007

An improvement is introduced to solve the plane problems of linear elasticity by reciprocal theorem for orthotropic materials. This method gives an integral equation with complex kernels which will be solved numerically. An artificial boundary is defined to eliminate the singularities and also an algorithm is introduced to calculate multi-valued complex functions which belonged to the kernels of the integral equation. The chosen sample problem is a plate, having a circular or elliptical hole, stretched by the forces parallel to one of the principal directions of the material. Results are compatible with the solutions given by Lekhnitskii for an infinite plane. Five different orthotropic materials are considered. Stress distributions have been calculated inside and on the boundary. There is no boundary layer effect. For comparison, some sample problems are also solved by finite element method and to check the accuracy of the presented method, two sample problems are also solved for infinite plate.

• Proceedings of the Earthquake Engineering Society of Korea Conference
• /
• 1998.10a
• /
• pp.399-406
• /
• 1998

Transient acoustic and elastic wave propagation in inhomogeneous media are studied by using the Fourier method. To verify the proposed numerical scheme, several examples having analytic solutions are considered, where two different semi-infinite media are in contact along a plane boundary. The comparisons of numerical results by the Fourier method and analytic solutions show good agreements. In addition, the Fourier method is applied to a layered half-plane, in which an elastic semi-infinite medium is covered by an elastic layer of finite thickness. It is showed how to derive the analytic solutions by using the Cagniard-de Hoop method. The numerical solutions are in excellent agreements with analytic results.

• Structural Engineering and Mechanics
• /
• v.78 no.2
• /
• pp.135-143
• /
• 2021

In this study, continuous contact problem in the functionally graded (FG) layer loaded with a FG flat punch resting on the elastic semi-infinite plane was analyzed by the finite element method (FEM). It was assumed that the shear modulus and density of the layer and punch varied according to exponentially throughout their depth. FG layer's weight was included to the problem and additionally all surfaces were considered as frictionless. Analysis of FG materials was performed with a special macro which was added to the ANSYS program. Firstly, the shear modulus of the punch was considered to be very rigid and the results of initial separation load (λcr) and distance (xcr) were compared with the analytical solution. Afterwards, results obtained from the contact analysis made according to the inhomogeneity parameters (β, γ) between FG punch-FG layer which had been unprecedented in the literature were discussed. As a result, FG punch's stress values at the punch edges where stress accumulations occurred were found to be smaller than the rigid punch. The security of the structure, longer life of the material and ease of production are directly related to the reduction of the stress values. The results obtained in this study are important in this respect. Also this work is the first study that investigates the effect of FG punch on the FG layer.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.1 no.4
• /
• pp.59-68
• /
• 1997

In analysis of underground structures, the effects of artificial boundary conditions are considered as one of the major reasons for differences from experimental results. These phenomena can be overcome by using the boundary elements which satisfy the multi-layered half space conditions. The fundamental solutions of multi-layered half-space for boundary element method is formulated satisfying the transmission and reflection of waves at each layer interface and radiation conditions at bottom layer. The governing equations can be obtained from the displacements at each layer which are expressed in terms of harmonic functions. All types of waves can be included using the complete response from semi-infinite integrals with respect to horizontal wavenumbers using expansion of Fourier series and Hankel transformation. Two dimensional Green's functions are derived from cylindrical Navier equations and potentials performing infinite integration in y-direction. In this case, it is effective to transform into two dimensional problem using semi-analytical integration and sinusoidal Bessel function. Some verifications are given to show the accuracy and efficiency of the developed method, and numerical examples to demonstrate the dynamic behavior of underground with various properties.

• Interaction and multiscale mechanics
• /
• v.3 no.1
• /
• pp.95-110
• /
• 2010

In recent years, study of the static response of pavements to moving vehicle and aircraft loads has received significant attention because of its relevance to the design of pavements and airport runways. The static response of beams resting on an elastic foundation and subjected to moving loads was studied by several researchers in the past. However, most of these studies were limited to steady-state analytical solutions for infinitely long beams resting on Winkler-type elastic foundations. Although the modelling of subgrade as a continuum is more accurate, such an approach can hardly be incorporated in analysis due to its complexity. In contrast, the two-parameter foundation model provides a better way for simulating the underlying soil medium and is conceptually more appealing than the one-parameter (Winkler) foundation model. The finite element method is one of the most suitable mathematical tools for analysing rigid pavements under moving loads. This paper presents an improved solution algorithm based on the finite element method for the static analysis of rigid pavements under moving vehicular or aircraft loads. The concrete pavement is discretized by finite and infinite beam elements, with the latter for modelling the infinity boundary conditions. The underlying soil medium is modelled by the Pasternak model allowing the shear interaction to exist between the spring elements. This can be accomplished by connecting the spring elements to a layer of incompressible vertical elements that can deform in transverse shear only. The deformations and forces maintaining equilibrium in the shear layer are considered by assuming the shear layer to be isotropic. A parametric study is conducted to investigate the effect of the position of moving loads on the response of pavement.

• Journal of the Mineralogical Society of Korea
• /
• v.10 no.2
• /
• pp.97-104
• /
• 1997

Based on the method of determination for relative stability of each phase from the difference among the interaction parameters of the phases consisting the mixed layer, the types of interactions between layers were specified and interaction parameter between layers in ordered domain was analytically derived as a function parameter between layers in ordered domain was analytically derived as a function of not only temperature and mole fraction of layers but also ordering parameter. Interaction parameter between the different layers in ordered phase, L is as follows:{{{{ {L }_{1 } (X,Q,T)= { C} over { Q} -4(1-2Q) { L}^{2 } - { RT} over {2} ln { 1} over {2 } - { 2RT} over { { X}_{ s} } ln { { 4QX}`_{s } ^{2 } } over {(1- { X}_{s }- { QX}_{s })( { X}_{s }- {QX }_{s } ) } }}}}L2 is the interaction parameter between ordered and disordered phase in domain and is the mole fraction of the domain which represent the infinite length of mixed layer mineral and Q and C are the reaction progress parameter and arbitrary constant, respectively. This equation was used for the I/S mixed layer clay minerals to infer the relative stability of R1 type I/S mixed layer in the temperature range from 373K to 450K. The result of calculation suggest that, owing to the decrease in interaction parameter with increasing temperature. The interaction parameter decreases more rapidly with decreasing mole fraction of smectite in domain, which is consistent with the fact that the probability of finding the series smectite layer is lo in the domain with small mole fraction of smectite layers in natural system.

• Journal of KSNVE
• /
• v.6 no.5
• /
• pp.671-677
• /
• 1996

This paper presents the analytical method for predicting turbulence- induced noise in the multilayered cylinder composed of an outer hose, an inner fluid and an internal core. It is assumed that an infinite axisymmetric cylinder is located horizontally in water with free stream velocity and the turbulent boundary layer (TBL) surrounding the outer hose is fully developed and homogeneous. The transfer function at the core surface due to the propagation of the pressure fluctuation within the TBL is formulated using the linearized Navier-Stockes equation for solid and fluid. In the estimation of the energy spectrum of wall pressure fluctuation, the empirical formula proposed by Strawderman based on the Corcos model is used. A general algorithm for the calculation of the pressure level at the surface of a core, that is, turbulence- induced noise, is presented. Through the detailed numerical simulation, it is found that the major noise mechanism is the propagation of the bulge wave along hose.