• 한국전산유체공학회:학술대회논문집
• /
• 2001.05a
• /
• pp.38-44
• /
• 2001

An efficient Gauss-Seidel time integration scheme is developed for solving the Euler and Navier-Stokes equations on unstructured meshes. Roe's FDS is used for the explicit residual computations and van Leer's FVS for evaluating implicit flux Jacobian. To reduce the memory requirement to a minimum level, off-diagonal flux Jacobian contributions are repeatedly calculated during the Gauss-Seidel sub-iteration process. Computational results based on the present scheme show that approximately $15\%$ of CPU time reduction is achieved while maintaining the memory requirement level to $50-60\%$ of the original Gauss-Seidel scheme.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.11
• /
• pp.2922-2931
• /
• 1994

We propose a modified 6-node element, where the sampling point of Gauss quadrature moved in the thickness direction. The modified 6-node element has been applied to static problems and forced motion analyses. In this study, this method is extended to the finite element analysis of the natural frequencies of two dimensional problems. We also propose a modified 16-node element for three dimensional problems, which behaves much like a 20-node element with smaller degree of freedom. The modified 6-node and 16-node elements have been applied to the modal analyses of beams and plates, respectively. The results agree well with the results of the 8-node or 20-node element models.

• Journal of IKEEE
• /
• v.19 no.2
• /
• pp.219-224
• /
• 2015

Electrical impedance tomography is an imaging technique to reconstruct the internal conductivity distribution based on applied currents and measured voltages in a domain of interest. In this paper, a modified Gauss-Newton method is proposed for conductivity image reconstruction. In the proposed method, the dimension of the inverse term is reduced by replacing the number of elements with the number of measurement data in the conductivity updating equation of the conventional Gauss-Newton method. Therefore, the computation time is greatly reduced as compared to the conventional Gauss-Newton method. Moreover, the regularization parameter is selected by computing the minimum-maximum from the diagonal components of the Jacobian matrix at every iteration. The numerical experiments with several scenarios were carried out to evaluate the reconstruction performance of the proposed method.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.12 no.4
• /
• pp.642-651
• /
• 1988

The present paper develops finite element procedures to calculate displacements, strains and stresses in initially stressed elastic solids subjected to static or time-dependent loading conditions. As a point of departure, we employ Hamilton's principle to obtain nonlinear equations of motion characterizing the displacement in a solid. The equations of motion reduce to linear equations of motion if incremental stresses are assumed to be infinitesimal. In the case of linear problem, finite element solutions are obtained by Newmark's direct integration method and by modal analysis. An analytic solution is referred to compare with the linear finite element solution. In the case of nonlinear problem, finite element solutions are obtained by Newton-Raphson iteration method and compared with the linear solution. Finally, the effect of the order of Gauss-Legendre numerical integration on the nonlinear finite element solution, has been investigated.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.11 no.3
• /
• pp.47-58
• /
• 1991

This paper presents dynamic infinite elements for soil-structure interaction analysis. In order to discretize the far field of the unbounded soil media, axisymmetric infinite elements which are capable of propagating multi-waves are proposed. An efficient numerical integration scheme for constructing the element characteristic matrices of the infinite elements in developed based on Gauss-Laguerre quadrature. The efficiency of the infinite elements is demonstrated by comparing the computed impedances of rigid circular footings on an elastic half space and on a layered half spaces with those obtained analytically.

• Transactions of the Korean Society of Machine Tool Engineers
• /
• v.13 no.2
• /
• pp.81-87
• /
• 2004

This paper deals with a new approach to tolerance optimization problems. Optimal tolerance allotment problems can be formulated as stochastic optimization problems. Most schemes to solve the stochastic optimization problems have been found to exhibit difficulties in multivariate integration of the probability density function. As a typical example of stochastic optimization the optimal tolerance allotment problem has the same difficulties. In this stochastic model, manufacturing system is represented by Gauss-Markov stochastic process and the manufacturing unit availability is characterized for realistic optimization modeling. The new algorithm performed robustly for a large deviation approximation. A significant reduction in computation time was observed compared to the results obtained in previous studies.

• Journal of Mechanical Science and Technology
• /
• v.17 no.4
• /
• pp.538-544
• /
• 2003

Although many methodologies exist for determining the constrained equations of motion, most of these methods depend on numerical approaches such as the Lagrange multiplier's method expressed in differential/algebraic systems. In 1992, Udwadia and Kalaba proposed explicit equations of motion for constrained systems based on Gauss's principle and elementary linear algebra without any multipliers or complicated intermediate processes. The generalized inverse method was the first work to present explicit equations of motion for constrained systems. However, numerical integration results of the equation of motion gradually veer away from the constraint equations with time. Thus, an objective of this study is to provide a numerical integration scheme, which modifies the generalized inverse method to reduce the errors. The modified equations of motion for constrained systems include the position constraints of index 3 systems and their first derivatives with respect to time in addition to their second derivatives with respect to time. The effectiveness of the proposed method is illustrated by numerical examples.

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

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.4
• /
• pp.990-1004
• /
• 1995

A modified 16-node element for composite plate has been proposed and compared with the 20-node element to check the validity of it. The fields of numerical inspection include mode analysis and specific damping analysis. By symetrizing the conventional unsymmetric damping matrix in the analysis of specific damping capacity, we could compute the specific damping capacity and make a program, effectively. In addition, we could predict the errors caused by reduction of integration order in thickness direction depending upon the number of layers.

• Computational Structural Engineering
• /
• v.10 no.2
• /
• pp.139-148
• /
• 1997

A finite element method is programmed to analyse the nonlinear behavior of axisymmetric structures. The lst order Mindlin shell theory which takes into account the transversal shear deformation is used to formulate a conical two node element with six degrees of freedom. To evade the shear locking phenomenon which arises in Mindlin type element when the effect of shear deformation tends to zero, the reduced integration of one point Gauss Quadrature at the center of element is employed. This method is the Updated Lagrangian formulation which refers the variables to the state of the most recent iteration. The solution is searched by Newton-Raphson iteration method. The tangent matrix of this method is obtained by a finite difference method by perturbating the degrees of freedom with small values. For the moment this program is limited to the analyses of non-linear elastic problems. For structures which could have elastic stability problem, the calculation is controled by displacement.