• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2009.10a
• /
• pp.321-321
• /
• 2009

This paper presents quasi-zero-stiffness (QZS) characteristic of a passive isolator using flexures under compression force. The passive isolator consists of a positive stiffness element (a vertical coil spring) and a negative stiffness element (flexures under compression force), and their proper combination of the positive and negative stiffness elements can produce both substantial static and zero dynamic stiffness, so called QZS. Firstly, a nonlinear dimensionless expression of a flexure under compression force is derived. A dynamic model of the passive isolator is developed and numerical simulations of its time and frequency response are performed. Then, undesirable nonlinear vibration is quantified using a period doubling bifurcation diagram and a Poincare's map of the isolator under forced excitation. Finally, experiments are performed to validate the QZS characteristic of the passive isolator.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2007.04a
• /
• pp.487-492
• /
• 2007

The enhancement of the service life of damaged or cracked structures is a major issue for researchers and engineers. The hierarchical void element with the integrals of Legend polynomials is used to characterize the fracture behavior of unpatched crack as well as repaired crack with bonded composite patches by computing the stress intensity factors and stress contours at the crack tip. The numerical approach is based on the v-version degenerate shell element including the theory of anisotropic laminated composites. Since the equivalent single layer approach is adopted in this study, the proposed element is necessary to represent a discontinuous crack part as a continuum body with zero stiffness of materials. Thus the aspect ratio of this element to represent the crack should be extremely slender. The sensitivity of numerical solution with respect to energy release rate, displacement and stress has been tested to show the robustness of hierarchical void element as the aspect ratio is increased up to 2000. The stiffness derivative method and displacement extrapolation method have been applied to calculate the stress intensity factors of Mode I problem.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2006.04a
• /
• pp.589-596
• /
• 2006

An improved numerical method to obtain the exact element stiffness matrix is newly proposed to perform the spatially coupled elastic and stability analyses of non-symmetric thin-walled beam-columns with two-types of elastic foundation. This method overcomes drawbacks of the previous method to evaluate the exact stiffness matrix for the spatially coupled stability analysis of thin-walled beam-column. This numerical technique is firstly accomplished via a generalized eigenproblem associated with 14 displacement parameters by transforming equilibrium equations to a set of first order simultaneous ordinary differential equations. Then exact displacement functions are constructed by combining eigensolutions and polynomial solutions corresponding to non-zero and zero eigenvalues, respectively. Consequently an exact stiffness matrix is evaluated by applying the member force-deformation relationships to these displacement functions.

• Journal of Korean Society of Steel Construction
• /
• v.14 no.4
• /
• pp.509-518
• /
• 2002

In order to perform the spatial buckling and static analysis of the nonsymmetric thin-walled beam-column element, improved exact static stiffness matrices were evaluated using equilibrium equation and force-deformation relationships. This numerical technique was obtained using a generalized linear eigenvalue problem, by introducing 14 displacement parameters and system of linear algebraic equations with complex matrices. Unlike the evaluation of dynamic stiffness matrices, some zero eigenvalues were included. Thus, displacement parameters related to these zero eigenvalues were assumed as polynomials, with their exact distributions determined using the identity condition. The exact displacement functions corresponding to three loadingcases for initial stress-resultants were then derived, by consistently combining zero and nonzero eigenvalues and corresponding eigenvectors. Finally, exact static stiffness matrices were determined by applying member force-displacement relationships to these displacement functions. The buckling loads and displacement of thin-walled beam were evaluated and compared with analytic solutions and results using ABAQUS' shell element or straight beam element.

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

• Magazine of the Korean Society of Agricultural Engineers
• /
• v.27 no.2
• /
• pp.38-46
• /
• 1985

The reinforced concrete farm silos on the elastic foundatin are widely used in agricultural engineering because of their superior structural performance, economy and attractive appearance. Various methods for the analysis and design of farm silo, such as the analytical method, the finite difference method, and the finite element methods, can be used. But the analytical procedure can not be applied for the intricate conditions in practice. Therefore lately the finite element method has been become in the structural mechanics. In this paper, a method of finite element analysis for the cylindrical farm silo on ffness matrix for the elastic foundation governed by winkler's assumption. A complete computer programs have been developed in this paper can be applicable not only to the shell structures on elastic foundation but also to the arbitrary three dimensional structures. Assuming the small deflection theory, the membrane and plate bending behaviours of flat plate element can be assumed mutually uncoupled. In this case, the element has 5 degrees of freedom per node when defined in the local coordinate system. However, when the element properties are transformed to the global coordinates for assembly, the 6th degree of freedom should be considered. A problem arises in this procedure the resultant stiffness in the 6th degree of freedom at this node will be zero. But this singularity of the stiffness matrix can be eliminated easily by merely replacing the zero diagonal by dummy stiffness.

• 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
• /
• 1996.04a
• /
• pp.182-189
• /
• 1996

The focus of this study is on the problem of the design of structure of undetermined topology. This problem has been regarded as being the most challenging of structural optimization problems, because of the difficulty of allowing topology to change. Conventional approaches break down when element sizes approach to zero, due to stiffness matrix singularity. In this study, a novel nonlinear Programming formulation of the topology Problem is developed and examined. Its main feature is the ability to account for topology variation through zero element sizes. Stiffness matrix singularity is avoided by embedding the equilibrium equations as equality constraints in the optimization problem. Although the formulation is general, two dimensional plane elasticity examples are presented. The design problem is to find minimum weight of a plane structure of fixed geometry but variable topology, subject to constraints on stress and displacement. Variables are thicknesses of finite elements, and are permitted to assume zero sizes. The examples demonstrate that the formulation is effective for finding at least a locally minimal weight.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.13 no.3
• /
• pp.317-322
• /
• 1989

Effective computational strategies have been proposed in the evaluation of stiffness matrices of rigid-plastic finite element method widely used in simulation of metal forming processes. The stiffness matrices are expressed as the sum of stiffness matrices evaluated by reduced integration and Liu's stabilization matrices which control the occurrence os zero-energy mode due to excessive reduced integration. The proposed method has been applied to the solution of fundamental 3-dimensional problems. The results clarified that the deformed mesh configuration was remarkably stabilized and computation speed attained about 3 times as fast as that of conventional 3-dimensional analyses. Furthermore, computation speed increases by a factor 60 when parallel computation is introduced. This speed has a tendency to increase as the total degree of freedom increases. As a result, this rigid-plastic finite element method enables us to analyze real 3-dimensional forming processes with practically acceptable computation time.

• Structural Engineering and Mechanics
• /
• v.39 no.4
• /
• pp.541-558
• /
• 2011

This paper presents an analysis of stochastic eigenvalue problem of plane bar structures. Particular attention is paid to the effect of spatial variations of the flexural properties of the structure on the first four eigenvalues. The problem of spatial variations of the structure properties and their effect on the first four eigenvalues is analyzed in detail. The stochastic eigenvalue problem was solved independently by stochastic finite element method (stochastic FEM) and Monte Carlo techniques. It was revealed that the spatial variations of the structural parameters along the structure may substantially affect the eigenvalues with quite wide gap between the two extreme cases of zero- and full-correlation. This is particularly evident for the multi-segment structures for which technology may dictate natural bounds of zero- and full-correlation cases.