• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.7 s.94
• /
• pp.1719-1730
• /
• 1993

This paper presents the theoretical natural frequencies of out-of-plane deformable ring based on the variables such as out-of-plane deflection, torsional rotation and shear rotation. Based on the same variables, a finite element eigen analysis is carried out by using the $C^0$-continuous, isoparametric element which has three nodes per element and three degrees-of-freedom at each node. Numerical experiments are peformed to find the integration scheme which produces accurate natural frequencies, natural modes and correct rigid body motion. The uniformly reduced integration and the selective reduced integration give more accurate numerical frequencies than the uniformly full integration, but the uniformly reduced integration produces incorrect rigid body motion while selective reduced integration does correct one. Therefore, the ring element based on the three variables which employes selective reduced integration is recommended to avoid spurious modes, to alleviate the error due to shear locking and to produce correct rigid body motion, simultaneously.

• Journal of the Korean Society for Precision Engineering
• /
• v.12 no.3
• /
• pp.32-41
• /
• 1995

In the analysis of metal forming processes by the finite element method, there are many numerical instabilities such as element locking, hourglass mode and shear locking. These instabilities may have a bad effect upon accuracy and convergence. The present work is concerned with improvement of stability and efficiency in two-dimensional rigid-plastic finite element method using various type of elemenmts and numerical intergration schemes. As metal forming examples, upsetting and backward extrusion are taken for comparison among the methods: various element types and numerical integration schemes. Comparison is made in terms of stability and efficiency in element behavior and computational efficiency and a new scheme of adaptive directional reduced integration is introduced. As a result, the finite element computation has been stabilized from the viewpoint of computational time, convergency, and numerical instability.

• Proceedings of the KSME Conference
• /
• 2004.11a
• /
• pp.453-458
• /
• 2004

This study proposes a new two-level condensation scheme for the construction of a reduced system. In the first step, the candidate area is selected for the construction of the reduced system by energy estimation in element-level. In the second step, primary degrees of freedom are selected by sequential elimination from the candidate degrees of freedom linked to the selected elements. Numerical examples demonstrate that the proposed method saves the computational cost effectively and provides a reduced system which predicts the eigenvalues accurately. Moreover, the well-constructed reduced system can present the reliable behavior of the structure under arbitrary dynamic loads comparing to that of global system. Time integration in a reduced system can save the computing time remarkably. Through a few numerical examples, the efficiency and reliability of the proposed scheme are verified.

• Structural Engineering and Mechanics
• /
• v.85 no.1
• /
• pp.89-104
• /
• 2023

This work presents a proposal for employing reduced numerical integration in the formulation of the 4-node quadrilateral solid finite element. The use of these low-order integration rules leads to numerical instabilities such as those producing the hourglass effect. The proposed procedure allows evaluating a given constitutive model only in one integration point, achieving an attractive computational cost reduction and, also, successfully controls the hourglass effect. A validation of the proposal is included and discussed throughout the paper. To show the efficiency of the proposal, several application examples of masonry structures are studied and discussed. To represent the non-linear mechanical behaviour of masonry a plastic-damage model is implemented within the application of this sub-integration scheme. Also, in order to have a full and computationally efficient strategy to determine the behaviour of masonry structures, involving its evolution to collapse, a homogenization technique with a macro-modeling approach is used. The methodology discussed throughout this paper demonstrates a substantial computational cost reduction and an improved approximation of the non-linear problem evidenced by a reduction of up to 85% of the computational time for some cases.

• Transactions of Materials Processing
• /
• v.4 no.4
• /
• pp.365-374
• /
• 1995

Numerical simulation of blade forming is carried out as stretch forming by an elasto-plastic finite element method. The method adopts a Lagrangian formulation, which incorporates large deformation and rotation, with a penalty method to treat the contact boundary condition. Numerical integration is done with a directional reduced integration scheme to avoid shear locking. The numerical results demonstrates various final shapes of blades which depend on the variation of the stretching force. The strain distributions in deformed blades are also obtained with the variation of the stretching force.

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

Based on numerical reduced minimization theory, new anisoparametric 3-node elements for out-of-plane curved beam are developed. The elements are designed to be free from spurious constraints. In this paper, the effect of the Jacobian upon numerical solution is analyzed and predicted through reduced minimization analysis of anisoparametric 3-node elements with different Jacobian assumption. The prediction is verified by numerical tests for circular and spiral out-of-plane deformable curved beam models. This paper proposes two kinds of 3-node elements with 7-DOF; one element employs 2-point integration for all strains, and the other element uses 3-point integration with a constant Jacobian within element for calculation of shear strain.

• Journal of the Korean Institute of Illuminating and Electrical Installation Engineers
• /
• v.29 no.10
• /
• pp.89-96
• /
• 2015

Explicit numerical integration methods for power system transient stability simulation require very small time steps to avoid numerical instability. The EXST1 exciter model is a primary source of fast dynamics in power system transients. In case of the EXST1, the required small integration time step for entire system simulation increases the computational demands in terms of running time and storage. This paper presents a practical exciter model reduction approach which allows the increase of the required step size and thus the method can decrease the computational demands. The fast dynamics in the original EXST1 are eliminated in the reduced exciter model. The use of a larger time step improves the computational efficiency. This paper describes the way to eliminate the fast dynamics from the original exciter model based on linear system theory. In order to validate the performance of the proposed method, case studies with the GSO-37 bus system are provided. Comparisons between the original and reduced models are made in simulation accuracy and critical clearing time.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2004.10a
• /
• pp.19-26
• /
• 2004

This study constructs the reduced system by two-level condensation scheme. This scheme consists of two steps. First step selects the candidate area for the primary degrees of freedom by energy estimation in element level. In the second step, the primary degrees of freedom are selected by the sequential elimination scheme. The efficiency and reliability of this scheme is shown through the prediction of eigenvalues of a few numerical examples. Time integration in the reduced system can save the computing time effectively. The well-constructed reduced system can present the accurate behavior of the structure under arbitrary dynamic loads so much as the global system. Through the numerical example, the efficiency and reliability of the proposed scheme will be demonstrated.

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

• Structural Engineering and Mechanics
• /
• v.8 no.6
• /
• pp.591-605
• /
• 1999

An accurate and efficient shell element is presented. The stiffness of the shell element is decomposed into two parts with one part corresponding to stretching and bending deformation and the other part corresponding to shear deformation of the shell. Both parts of the stiffness are calculated with reduced integration rules, thereby improving computational efficiency. Shear strains are averaged on the reference surface such that neither locking phenomena nor any zero energy mode can occur. The satisfactory behaviour of the element is demonstrated in several numerical examples.