• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.18 no.2
• /
• pp.123-131
• /
• 2005

According to ow previous study, we confirmed That the Petrov-Galerkin natural element method(PG-NEM) completely resolves the numerical integration inaccuracy in the conventional Bubnov-Galerkin natural element method(BG-NEM). This paper is an extension of PG-NEM to two-dimensional geometrically nonlinear problem. For the analysis, a linearized total Lagrangian formulation is approximated with the PS-NEM. At every load step, the grid points ate updated and the shape functions are reproduced from the relocated nodal distribution. This process enables the PG-NEM to provide more accurate and robust approximations. The representative numerical experiments performed by the test Fortran program, and the numerical results confirmed that the PG-NEM effectively and accurately approximates The large deformation problem.

• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 2004.11a
• /
• pp.179-182
• /
• 2004

This paper deals with a new mathematical formulation of nonlinear wave profile based on Banach fixed point theorem. As application of the formulation and its solution procedure, some numerical solutions was presented in this paper and nonlinear equation was derived. Also we introduce a new operator for iteration and getting solution. A numerical study was accomplished with Stokes' first-order solution and iteration scheme, and then we can know the nonlinear characteristic of Stokes' high-order solution. That is, using only Stokes' first-oder(linear) velocity potential and an initial guess of wave profile, it is possible to realize the corresponding high-oder Stokian wave profile with tile new numerical scheme which is the method of iteration. We proved the mathematical convergence of tile proposed scheme. The nonlinear strategy of iterations has very fast convergence rate, that is, only about 6-10 iterations arc required to obtain a numerically converged solution.

• Steel and Composite Structures
• /
• v.4 no.6
• /
• pp.489-507
• /
• 2004

The flexibility of the connection between steel and concrete largely influences the global behaviour of the composite beam. Therefore the way the connection is modelled is the key issue in its structural analysis. Here we present a new strain-based finite element formulation in which we consider non-linear material and contact models. The computational efficiency and accuracy of the formulation is proved with the comparison of our numerical results with the experimental results of Abdel Aziz (1986) obtained in a full-scale laboratory test. The shear connectors are assumed to follow a non-linear load-slip relationship proposed by Ollgaard et al. (1971). We introduce the notion of the generalized slip, which offers a better physical interpretation of the behaviour of the contact and gives an additional material slip parameter. An excellent agreement of experimental and numerical results is obtained, using only a few finite elements. This demonstrates that the present numerical approach is appropriate for the evaluation of behaviour of planar composite beams and perfect for practical calculations.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.27 no.1
• /
• pp.37-43
• /
• 2014

The extended framework of Hamilton's principle provides a new rigorous weak variational formalism for a broad range of initial boundary value problems in mathematical physics and mechanics in terms of mixed formulation. Based upon such framework, a new variational numerical method of linear elasticity is provided for the classical single-degree-of-freedom dynamical systems. For the undamped system, the algorithm is symplectic with respect to the time step. For the damped system, it is shown to be accurate with good convergence characteristics.

• Structural Engineering and Mechanics
• /
• v.82 no.3
• /
• pp.401-416
• /
• 2022

The aim of this study is to examine the frictionless double receding contact problem for two functionally graded (FG) layers pressed with a uniformly distributed load and resting on a homogeneous half plane (HP) using analytical and numerical methods. The FG layers are made of a non-homogeneous material with an isotropic stress-strain law with exponentially varying properties. It is assumed that the contact at the FG layers and FG layer-HP interface is frictionless. The body force of the FG layers and homogeneous HP are ignored in the study. Firstly, an analytical solution for the contact problem has been realized using the theory of elasticity and the Fourier integral transform techniques. Then, the problem modeled and two-dimensional analysis was carried out by using the ANSYS package program based on FEM. Numerical results for contact lengths and contact pressures between FG layers and FG layer-HP were provided for various dimensionless quantities including material inhomogeneity, distributed load width, the shear module ratio, and the heights of the FG layers for both methods. The results obtained using FEM were compared with the results found using the analytical formulation. It was found that the results obtained from analytical formulation were in perfect agreement with the FEM study.

• Geophysics and Geophysical Exploration
• /
• v.15 no.2
• /
• pp.57-65
• /
• 2012

Various numerical methods in simulation of seismic wave propagation have been developed. Recently an innovative numerical method called as the Spectral Element Method (SEM) has been developed and used in wave propagation in 3-D elastic media. The SEM that easily implements the free surface of topography combines the flexibility of a finite element method with the accuracy of a spectral method. It is generally used a weak formulation of the equation of motion which are solved on a mesh of hexahedral elements based on the Gauss-Lobatto-Legendre integration rule. Variational formulations of velocity-stress motion are newly modified in order to implement ADE-PML (Auxiliary Differential Equation of Perfectly Matched Layer) in wave propagation in 3-D elastic media, because a general weak formulation has a difficulty in adapting CFS (Complex Frequency Shifted) PML (Perfectly Matched Layer). SEM of Velocity-Stress motion having ADE-PML that is very efficient in absorbing waves reflected from finite boundary is verified with simulation of 1-D and 3-D wave propagation.

• Proceedings of the KIEE Conference
• /
• 1988.07a
• /
• pp.950-953
• /
• 1988

The equations of motion for linearly elastic bodies undergoing large displacement motion are derived. This produces a set of equations which are efficient to numerically integrate. The equations for the elastic bodies are formulated and simplified to provide as much efficiency as possible in their numerical solution. A futher efficiency is obtained through the use of floating reference frame. The equation are presented in two forms for numerical integration. 1) Explicit numerical integration 2) Implicit numerical integration. In this paper, there was used the numerical integration. The implicit numerical integration is extended to solved second order equation, futher reducing the numerical effort required. The formulation given is seen to be occulate and is expected to be efficient for many types of problems.

• Structural Engineering and Mechanics
• /
• v.52 no.2
• /
• pp.239-260
• /
• 2014

Non-stationary random vibration of linear structures with uncertain parameters is investigated in this paper. A time-domain explicit formulation method is first presented for dynamic response analysis of deterministic structures subjected to non-stationary random excitations. The method is then employed to predict the random responses of a structure with given values of structural parameters, which are used to fit the conditional expectations of responses with relation to the structural random parameters by the response surface technique. Based on the total expectation theorem, the known conditional expectations are averaged to yield the random responses of stochastic structures as the total expectations. A numerical example involving a frame structure is investigated to illustrate the effectiveness of the present approach by comparison with the power spectrum method and the Monte Carlo simulation method. The proposed method is also applied to non-stationary random seismic analysis of a practical arch bridge with structural uncertainties, indicating the feasibility of the present approach for analysis of complex structures.

• Structural Engineering and Mechanics
• /
• v.11 no.5
• /
• pp.547-564
• /
• 2001

The numerical formulation of a two-phase interface element appropriate for porous lining system is presented. The formulation is isoparametric and can be applied both for 2-D and 3-D analysis. Biot's theory is utilized as the basis for the development of the element constitutive theory. In order to be capable of simulating the reinforcing characteristics of some geotextiles utilized as lining system, a reinforcement component has also been implemented into the formulation. By employing this specially developed interface finite element, the influence of soil consolidation on the stress distribution along the lining system of a reservoir and a landfill has been investigated.

• Structural Engineering and Mechanics
• /
• v.50 no.3
• /
• pp.403-417
• /
• 2014

The present paper addresses a general formulation for the thermo elastic analysis of a functionally graded cylindrical shell subjected to external loads. The shear deformation theory and energy method is employed for this purpose. This method presents the final relations by using a set of second order differential equations in terms of integral of material properties along the thickness direction. The proposed formulation can be considered for every distribution of material properties, whether functional or non functional. The obtained formulation can be used for manufactured materials or structures with numerical distribution of material properties which are obtained by using the experiments. The governing differential equation is applied for two well-known functionalities and some previous results are corrected with present true results.