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

• Geosystem Engineering
• /
• v.21 no.6
• /
• pp.309-317
• /
• 2018

We employ a three-dimensional indirect boundary element method (BEM) to simulate temperature change around an underground liquefied natural gas storage cavern. The indirect BEM (IBEM) uses fictitious heat source strength on boundary elements as basic variables which are solved from equations of boundary conditions and then used to compute the temperature change at other points in the considered problem domain. The IBEM requires evaluation of singular integration for temperature change due to heat conduction from a constant heat source on a planar (triangular) region. The singularity can be eliminated by a semi-analytical integration scheme. However, it is found that the semi-analytical integration scheme yields sharp temperature gradient for points close to vertices of triangle. This affects the accuracy of heat flux, if they are evaluated by finite difference method at these points. This difficulty can be overcome by a combination of using a direct numerical integration for these points and the semi-analytical scheme for other points distance away from the vertices. The IBEM and the hybrid integration scheme have been verified with an analytic solution and then used to the application of the underground storage.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.20 no.12
• /
• pp.3792-3803
• /
• 1996

Since the skew beam element has two curvatures which are a curvature and a torsion, spatial behavior of curved beam which cannot be included in one plane can be anlayzed by emploting the skew beam element. The $C^{0}$-continuous skew beam element shows the stiffness locking phenomenon when full integration is employed. The locking phenomenpn is characterized by two typical phenomena ; one is the much smaller displacement thant the exact one and theother is the undelation phenomenon is stress distribution. In this paper, we examine how unmatched coefficient in the constrained energy brings about the locking by Reduced Minimization theory. We perform the numerical ones. These comparisons show that uniformly full integration(UFI), which employs full integration for the constrained energy, entails the locking phenomenon. But the use of uniformly reduced integration(URI) of selectively reduced integration(SRI), which employs reduced integration for constrained energy, does not produce the significant errors of displacements of the undulation phenomenon in stress distribution since they do not entails the locking, Additionally, the error due to the approximated parameters for describing the geometry of skew beam is examined.d.

• Structural Engineering and Mechanics
• /
• v.21 no.5
• /
• pp.539-551
• /
• 2005

Finite elements based on isoparametric formulation are known to suffer spurious stiffness properties and corresponding stress oscillations, even when care is taken to ensure that completeness and continuity requirements are enforced. This occurs frequently when the physics of the problem requires multiple strain components to be defined. This kind of error, commonly known as locking, can be circumvented by using reduced integration techniques to evaluate the element stiffness matrices instead of the full integration that is mathematically prescribed. However, the reduced integration technique itself can have a further drawback - rank deficiency, which physically implies that spurious energy modes (e.g., hourglass modes) are introduced because of reduced integration. Such instability in an existing stiffness matrix is generally detected by means of an eigenvalue test. In this paper we show that a knowledge of the dimension of the solution space spanned by the column vectors of the strain-displacement matrix can be used to identify the instabilities arising in an element due to reduced/selective integration techniques a priori, without having to complete the element stiffness matrix formulation and then test for zero eigenvalues.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2008.04a
• /
• pp.172-176
• /
• 2008

In this paper, three dimensional linear conforming variable-finite elements are presented with the aid of a smoothed integration (a class of stabilized conforming nodal integration), for mnltiscale mechanics problems. These elements meet the desirable properties of an interpolation such as the Kronecker delta condition, the partition of unity condition and the positiveness of interpolation function. The necessary condition of linear exactness is fully relaxed by employing the smoothed integration, which renders us to meet the linear exactness in a straightforward manner. This novel element description extend the category of the conventional finite elements space to ration type function space and give the flexibility on the number of nodes of element which are fixed in the conventional finite elements. Several examples are provided to show the convergence and the accuracy of the proposed elements, and to demonstrate their potential with emphasis on the multiscale mechanics problems such as global/local analysis, nonmatching contact problems, and modeling of composite material with defects.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1997.10a
• /
• pp.87-94
• /
• 1997

The circular plate resting on Boussinesq's half-space model under axisymmetric loading is studied by a finite element procedure to evaluate the distribution of contact pressure between plate and elastic half-space. The displacement of half-space due to axisymmetric surface loading can be evaluated by double integration of Boussinesq's solution. On that case the analytical integration can be executed for the radial direction but the analytical integration for the circumferential direction is impossible and the numerical integration should be considered. With the radial integration we can get non-dimensional function. Then the numerical integration for the formula is executed for the circumferential direction and the results are approximated 5th order Polynomials by using the least square method. With these 5th order approximate formula, the flexibility matrix of half-space is constructed as the coefficient matrix of nodal contact pressure by the finite element procedures. Iteration procedures are attempted by using this method to determine the separated region.

• Coupled systems mechanics
• /
• v.2 no.4
• /
• pp.389-410
• /
• 2013

Three-dimensional Lagrangian fluid finite element is applied to seismic response analysis of an oil storage tank with a floating roof. The fluid element utilized in the present analysis is formulated based on the displacement finite element method considering only volumetric elasticity and its element stiffness matrix is derived by using one-point integration method in order to avoid volumetric locking. The method usually adds a rotational penalty stiffness to satisfy the irrotational condition for fluid motion and modifies element mass matrices through the projected mass method to suppress spurious hourglass-mode appeared in compensation for one-point integration. In the fluid element utilized in the present paper, a small hourglass stiffness is employed. The fluid and structure domains for the objective oil storage tank are modeled by eight-node solid elements and four-node shell elements, respectively, and the transient response of the floating roof structure or the free surface are evaluated by implicit direct time integration method. The results of seismic response analyses are compared with those by other method and the validation of the present analysis using three-dimensional Lagrangian fluid finite elements is shown.

• Proceedings of the KSME Conference
• /
• 2001.06c
• /
• pp.402-407
• /
• 2001

Sensitivity analysis scheme is developed in the elasto-plastic finite element method with explicit time integration using direct differentiation method. The direct differentiation is concerned with the time integration, constitutive relation, shell element with reduced integration and the contact scheme. Sensitivity analysis results are mainly examined with the highly nonlinear and quasi-static problem with the complicated contact condition. The result shows stable sensitivity especially in the sheet metal forming analysis.

• Proceedings of the KSME Conference
• /
• 2001.11a
• /
• pp.324-329
• /
• 2001

A novel method for non-matching interfaces on the boundaries of the finite elements in partitioned domains is presented by introducing interface elements in this paper. The interface element method (IEM) satisfies the continuity conditions exactly through interfaces without recourse to the Lagrange multiplier technique. The moving least square (MLS) approximation in the present study is implemented to construct the shape functions of the interface elements. Alignment of the boundaries of sub-domains in the MLS approximation and integration domains provides a consistent numerical integration due to one form of rational functions in an integration domain. The compatibility of displacements on the boundaries of the finite elements and the interface elements is always preserved in this method, and the completeness of the shape functions of the interface elements guarantees the convergence of numerical solutions. The numerical examples show that the interface element method is a useful tool for the analysis of a partitioned system and for a global-local analysis.

• Korean Journal of Computational Design and Engineering
• /
• v.3 no.1
• /
• pp.48-56
• /
• 1998

Finite element analysis is very important in the design and analysis of mechanical engineering. The process of FEA encompasses shape modeling, mesh generation, matrix solving and post-processing. Some of these processes can be tightly integrated with the current software architectures and data sharing mode. However, complete integration of all the FEA process itself and the integration to the manufacturing processes is almost impossible in the current practice. The barriers to this problem are inconsistent data format and the enterprise-wise software integration technology. In this research, the information model based on STEP AP209 was chosen for handling finite element analysis data. The international standard for the FEA data can bridge the gap between design, analysis and manufacturing processes. The STEP-based FEA system can be further tightly integrated to the distributed software and database environment using CORBA technology. The prototype FEA system DICESS is implemented to verify the proposed concepts.