• Transactions of the Korean Society of Mechanical Engineers
• /
• v.14 no.5
• /
• pp.1104-1111
• /
• 1990

A method to identify nonlinear elements position of a nonlinear system is presented. Nonlinear elements position can be identified by an equivalent error damping and stiffness matrices which are based on the equivalent linearization technique. The procedures of this technique are: (1) Obtain input force and system response. (2) Define error between the actual and linearized restoring forces. (3) Calculate linearized damping and stiffness coefficients to minimize the square error sum. Several examples are tested and found that these methods are very effective not only to locate the nonlinear elements position but also to identify the degree of nonlinearity qualitatively. Nonlinear type can be qualitatively identified by examining the plots of restoring force vs equivalent state values.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.33 no.11
• /
• pp.7-14
• /
• 2005

A subdomain-interface finite element method is suggested to solve a class of fully- coupled thermomechanical problems with contact boundaries. The penalty method is used for connecting subdomains that satisfy interface compatibility conditions. As a result, effective stiffness matrices are always positive definite, and computational efficiency can be improved to a considerable degree. Moreover, any complex-shaped domain can be divided into independently modeled subdomains without considering the conformity of meshes on interfaces. Using a computer code based on the present method, these advantageous features are shown through a set of numerical studies.

• Journal of the Korean Society of Hazard Mitigation
• /
• v.7 no.5
• /
• pp.27-35
• /
• 2007

This study presents a non -prismatic beam element for modeling the elastic and inelastic behavior of the steel beam, which has the post-Northridge connections in steel moment frames that are subjected to earthquake ground motions. The elastic stiffness matrix for non-prismatic members with reduced beam section (RES) connection is in the closed-form. The plasticity model is of a discrete type and is composed of a series of nonlinear hinges connected by rigid links. The hardening rules can model the inelastic behavior for monotonic and random cyclic loading, and the effects of local buckling. Verification and calibration of the model are presented in a companion paper.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.9 no.6 s.46
• /
• pp.67-73
• /
• 2005

Precise analysis of soil-structure interaction requires a proper description of soil profile. However, such approach becomes generally nearly unpractical for soil exhibiting material discontinuity and complex geometry since meshes should match that material discontinuity line. To overcome these difficulties, a different numerical integration method is adopted in this paper, which enables to integrate easily over an element with material discontinuity without regenerating mesh fellowing the discontinuity line. As a result the mesh is highly structured, loading to very regular silliness matrix. The influence of the shape of soil profile on the response is examined and it is seen that the proposed soil-structure analysis method can be easily used on soil-structure interaction problems with complicated soil profile and produce reliable results regardless of material discontinuities.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.16 no.1
• /
• pp.1-8
• /
• 2003

Topology optimization has been evolved into a very efficient conceptual design tool and has been utilized into design engineering processes. Traditional topology optimization has been using homogenization method and optimality criteria method. homogenization method provides relationship equation between structure which includes many holes and stiffness matrix in FEM. Optimality criteria method is used to update design variables while maintaining that volume fraction is uniform. Traditional topology optimization has advantage of good convergence but has disadvantage of too much convergency time. In one way to solve this problem, element removal method using the criterion of an average stress is presented. As the result of examples, it is certified that convergency time is very reduced.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.19 no.3 s.73
• /
• pp.283-293
• /
• 2006

Recently, various types of structural systems for skyscrapers are studied as the height and size of the building structures rapidly increase due to social and economical needs. The mega frame system among them, which is the structural system developed recently, is known as a suitable structural system for skyscrapers because this structural system has sufficient stiffness against the lateral forces by combination of mega members which consist of many columns and girders. Since the mega frame structure has significant numbers of elements and nodes, it takes tremendous times and computer memories to analyze and design the structures. Therefore, the exclusive structural analysis program for mega frame system is developed to reduce the efforts and time required for the analysis and design of mega frame structure. To this end, an efficient modelling technique using the characteristics of mega frame structures and an efficient analytical model, which uses a few DOFs selected by the user using the matrix condensation method, are developed in this study. Static and dynamic analyses are conducted using an example structure. The effectiveness and accuracy of the developed program we verified by the comparison between the results of the proposed method and the conventional method.

• Journal of the Korean Society for Precision Engineering
• /
• v.27 no.2
• /
• pp.109-116
• /
• 2010

It is important to make a mechanical structure precisely and reasonably in predicting the dynamic characteristics, controlling the vibration, and designing the structure dynamics. In finite element analysis model updating is appropriate as the design parameter is used to analyze the dynamic system. The errors can be contained from the physical parameters and the element modeling. From the dynamic test, more precise dynamic characteristics can be obtained. In this paper, model updating algorithm is developed using frequency difference between experiment and calculation. Modal frequencies are obtained by experiment and finite element analysis for beams with various cross section and shapes which have added masses and holes in the middle. For plates with and without groove, experiment and analyses are carried out by applying free boundary conditions as well. Mass and stiffness matrices are updated by comparing test and analytical modal frequencies. The result shows that the updated frequencies become closer to the test frequencies in case that both matrices are updated. An improved analytical model is obtained by changing model parameters such that the discrepancy between test and finite element frequencies is minimized. For beam and plate models updating of mass and stiffness matrices can improve the dynamical behavior of the model by acting on the physical parameters such as masses and stiffness.

• Journal of Korean Society of Steel Construction
• /
• v.10 no.1 s.34
• /
• pp.47-61
• /
• 1998

The behavior of horizontally curved I-girder bridges is complex because the flexural and torsional behavior of curved girders are coupled due to their initial curvature. Also, the behavior is affected by cross beams. To investigate the behavior of horizontally curved I-girder bridges, it is necessary to consider curved girders with cross beams. In order to perform free vibration analyses of horizontally curved I-girder bridges, a finite element formulation is presented here and a finite element analysis program is developed. The formulation that is presented here consists of curved and straight beam elements, including the warping degree of freedom. Based on the theory of thin-walled curved beams, the shape functions of the curved beam elements are derived from homogeneous solutions of the static equilibrium equations. Third-order hermits polynomials are used to form the shape functions of the straight beam elements. In the finite element analysis program, global stiffness and mass matrix are composed, based on the Cartesian coordinate system. The Gupta method is used to efficiently solve the eigenvalue problem. Comparing the results of several examples here with those of previous studies, the formulation presented is verified. The validity of the program developed is shown by comparing results with those analyzed by the shell element.

• Journal of the Korean Geotechnical Society
• /
• v.28 no.4
• /
• pp.67-78
• /
• 2012

A three-dimensional approximate computer-based method, YSPR (Yonsei Piled Raft), was developed for analysis of behavior of piled raft foundations. The raft was modeled as a flat shell element having 6 degrees of freedom at each node and the pile was modeled as a beam-column element. The behaviors of pile head and soil were controlled by using $6{\times}6$ stiffness matrix. To model the non-linear behavior, the soil-structure interaction between soil and pile was modeled by using nonlinear load-transfer curves (t-z, q-z and p-y curves). Comparison with previous model and FEM analysis showed that YSPR gave similar load-displacement behaviors. Comparison with field measurement also indicated that YSPR gave a reasonable result. It was concluded that YSPR could be effectively used in analysis and design of piled raft foundations.

• Journal of the Korean Society for Railway
• /
• v.15 no.6
• /
• pp.588-596
• /
• 2012

This paper presents a new semi-analytical annular Mindlin plate element with which out-of-plane natural vibration of thick disks can be analyzed simply, efficiently, and accurately through FEM by including effects of rotary inertia and transverse shear deformation. Using static deformation modes which are exact solutions of equilibrium equations of annular Mindlin plate, the element interpolation functions, stiffness and mass matrices corresponding to each number of nodal diameters are derived. The element is capable of representing out-of-plane rigid-body motions exactly and free from shear locking. Natural frequencies of uniform and multi-step disks with or without concentric ring support are analyzed by applying the presented element. Such results are compared with theoretical predictions of previous works or FEA results obtained by using two-dimensional shell element to investigate the convergence and accuracy of the presented element.