• International Journal of Railway
• /
• v.3 no.4
• /
• pp.117-122
• /
• 2010

In transportation geotechnical engineering, stress-strain behavior of earth structures has been analyzed by numerical simulations with the implemented plasticity constitutive model. It is a fact that many advanced plasticity constitutive models on predicting the mechanical behavior of soils have been developed as well as experimental research works for geotechnical applications in the past decades. In this study, recently developed, a unified constitutive model for both clay and sand, which is referred to as CASM (clay and sand model), was compared with a classical constitutive model, Cam-Clay model. Moreover, integration methods of stress-strain rate equations using CASM were presented for simulation of undrained and drained triaxial compression tests. As a conclusion, it was observed that semi-implicit integration method has more improved accuracy of capturing strain rate response to applied stress than explicit integration by the multiple correction and iteration.

• Proceedings of the KIEE Conference
• /
• 1987.07a
• /
• pp.105-109
• /
• 1987

A solution which combines ray and aperture integration(AI) techniques is presented for the problem of electromagnetic plane wave scattering by an open-ended, perfectly-conducting, semi-infinite parallel plate waveguide with a thin, uniform layer of lossy or absorbing material on its inner walls, and with a simple planar termination inside. Numerical results are given for the fields outside the waveguide.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.25 no.9
• /
• pp.1376-1383
• /
• 2001

Meshless methods, developed in various ways over the past decade, have been attractive as new computational methods in that they do not need mesh generation in analyzing procedure. But most of these methods were not truly meshless methods because background meshes were required for the spatial integration of a weak form. Accordingly, in this paper, nodal integration for truly meshless methods has been studied, and an improvement scheme is proposed. To improve stabilization and accuracy, which are the weak points in previous nodal integration methods, the integration area is transformed to circle and then numerically integrated. This method does not need any adding term for stabilization in the variational formulation and then simplifies the integration procedure. Numerical test results show that the proposed method is more accurate, stable, and reasonable than the existed nodal integration methods.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.10
• /
• pp.2640-2652
• /
• 1994

The stiffness of 6-node isotropic element is stiffer than that of 8-node isotropic element of same configuration. This phenomenon was called 'Relative Stiffness Stiffening Phenomenon'. In this paper, an equation of sampling point modification which correct this phenomenon was derived for the composite plate, as well as an equation for an isotropic plate. The relative stiffness stiffening phenomena of an isotropic plate element could be corrected by modifying Gauss sampling points in the numerical integration of stiffness matrix. This technique could also be successfully applied to the static analyses of composite plate modeled by the 3-dimensional 16-node elements. We predicted theoretical errors of stiffness versus the number of layers that result from the reduction of numerical integration order. These errors coincide very well with the actual errors of stiffness. Therefore, we can choose full integration of reduced integration based upon the permissible error criterion and the number of layers by using the thoretically predicted error.

• International Journal of Aeronautical and Space Sciences
• /
• v.7 no.1
• /
• pp.13-26
• /
• 2006

When the stagnation temperature of a perfect gas increases, the specific heat for constant pressure and ratio of the specefic heats do not remain constant any more and start to vary with this temperature. The gas remains perfect: its state equation remains always valid, with exception that it will be named by calorically imperfect gas. The aim of this research is to develop the relations of the necessary thermodynamics and geometrical ratios. and to study the supersonic flow at high temperature. lower than the threshold of dissociation. The results are found by the resolution of nonlinear algebraic equations and integration of complex analytical functions where the exact calculation is impossible. The dichotomy method is used to solve the nonlinear equation. and the Simpson algorithm for the numerical integration of the found integrals. A condensation of the nodes is used. Since. the functions to be integrated have a high gradient at the extremity of the interval of integration. The comparison is made with the calorifcally perfect gas to determine the error made by this last. The application is made for the air in a supersonic nozzle.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2004.04a
• /
• pp.333-340
• /
• 2004

This paper deals with geometric nonlinear analyses using a new meshfree technique which improves the numerical integration accuracy. The new method called the Petrov-Galerkin natural element method (PGNEM) is based on the Voronoi diagram and the Delaunay triangulation which is based on the same concept used for conventional natural element method called the Bubnov-Galerkin natural element method (BGNEM). But, unlike BGNEM, the test shape function is differently chosen from the trial shape function. In the linear static analysis, it is ensured that the numerical integration error of the PGNEM is remarkably reduced. In this paper, the PGNEM is applied to large deformation problems, and the accuracy of the proposed numerical technique is verified through the several examples.

• Structural Engineering and Mechanics
• /
• v.75 no.3
• /
• pp.357-367
• /
• 2020

The objective of this article is investigation of dynamic response of thick multilayer functionally graded (FG) beam under generalized dynamic forces. The plane stress problem is exploited to describe the constitutive equation of thick FG beam to get realistic and accurate response. Applied dynamic forces are assumed to be sinusoidal harmonic, sinusoidal pulse or triangle in time domain and point load. Equations of motion of deep FG beam are derived based on the Hamilton principle from kinematic relations and constitutive equations of plane stress problem. The numerical finite element procedure is adopted to discretize the space domain of structure and transform partial differential equations of motion to ordinary differential equations in time domain. Numerical time integration method is used to solve the system of equations in time domain and find the time responses. Numerical parametric studies are performed to illustrate effects of force type, graduation parameter, geometrical and stacking sequence of layers on the time response of deep multilayer FG beams.

• The Korean Journal of Applied Statistics
• /
• v.13 no.2
• /
• pp.477-488
• /
• 2000

A distance measure between two Latin-hypercube designs is defined and its expected value is computed. It was computed by using mathematical statistics, numerical analysis (multidimensional numerical integration), Monte-carlo method, and the theory of asymptotic normal distribution. For the comparison of two Latin-hypercube designs with same structure but different randomness, the difference of expected values of response function and information mass of experimental designs are considered. These methods may be useful in comparison between two general experimental designs.

• Journal of Mechanical Science and Technology
• /
• v.17 no.1
• /
• pp.64-76
• /
• 2003

In this work, a new penalty formulation is proposed for the analysis of Mindlin-Reissner plates by using the element-free Galerkin method. A penalized weak form for the Mindlin-Reissner Plates is constructed through the exterior penalty method to enforce the essential boundary conditions of rotations as well as transverse displacements. In the numerical examples, some typical problems of Mindlin-Reissner plates are analyzed, and parametric studies on the order of integration and the size of influence domain are also carried out. The effect of the types of background cells on the accuracy of numerical solutions is observed and a proper type of background cell for obtaining optimal accuracy is suggested. Further, optimal order of integration and basis order of Moving Least Squares approximation are suggested to efficiently handle the irregularly distributed nodes through the triangular type of background cells. From the numerical tests, it is identified that unlike the finite element method, the proposed element-free Galerkin method with penalty technique gives highly accurate solution without shear locking in dealing with Mindlin-Reissner plates.

• Journal of Mechanical Science and Technology
• /
• v.19 no.spc1
• /
• pp.481-486
• /
• 2005

This paper introduces a numerical method for integration of the linear and nonlinear differential dynamic equation of motion. The variation of acceleration in two time steps is approximated as a combination of both trigonometric cosine and hyperbolic cosine functions with weighted coefficient. From which all necessary formulae are elaborated for the direct integration of the governing equation. A number of linear and nonlinear dynamic problems with various degrees of freedom are analysed using both the suggested method and Newmark method for the comparison. The numerical results show high advantages and effectiveness of the new method.