• Journal of the Society of Naval Architects of Korea
• /
• v.30 no.1
• /
• pp.145-156
• /
• 1993

An improved analysis model for material nonlinearity induced by elasto-plastic deformation and damage including large strain response was proposed. The elasto-plastic-damage constitutive model based on the continuum damage mechanics approach was adopted to overcome limitations of the conventional plastic theory, which can manage the anisotropic tonsorial damages evolved during time-independent plastic deformation process of materials. Updated Lagrangian finite element formulation for elasto-plastic damage coupling problem including large deformation, large rotation and large strain problems was completed to develop a numerical model which can predict all kinds of structural nonlinearities and damage rationally. Finally, a finite element analysis code for the 2-dimensional plane problem was developed and the applicability and validity of the numerical model was investigated through some numerial examples. Calculations showed reasonable results in both geometrical nonlinear problem due to large deformation and material nonlinearity including the damage effect.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2011.04a
• /
• pp.760-763
• /
• 2011

본 논문의 목적은 효율적인 3절점 판 유한요소를 개발하는 것이다. Hellinger-Reissner 범함수에 근거한 혼합정식화(mixed formulation)를 사용한다. 잠김현상을 일으키는 전단변형률장을 독립적으로 분리한 후, MITC(Mixed Interpolation of Tensorial Components)방법을 이용하여 대체전단변형률장(assumed transverse shear strain field)을 구성한다. 추가적으로 회전된 반변기저벡터(contravariant base vector)로 정의된 근사전단변형률장(approximated transverse shear strain field)에 미지수(unknowns)를 도입하여 혼합정식화를 완성시키고 정적응축(static condensation)을 통해 최종 강성행렬을 구성한다. 거짓영에너지모드시험(spurious zero energy mode test), 조각시험(patch test), 등방성시험(isotropic test) 등을 실시하였으며, 4변 완전구속 정사각형 판 구조물과 60도 기울어진 단순지지 판 구조물 등 예제들을 해석하여 MITC3판 유한요소와 수렴성능을 비교하였다.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.20 no.3
• /
• pp.303-309
• /
• 2007

For the propagation of elastic waves in unbounded domains, absorbing boundary conditions at the fictitious numerical boundaries have been proposed. In this paper we focus on both first and second order paraxial boundary conditions(PBCs) in the framework of variational approximations which are based on paraxial approximations of the scalar and elastic wave equations. We propose a penalty function method for the treatment of PBCs and apply these into finite element analysis. The numerical verification of the efficiency is carried out through comparing PBCs with Lysmer-Kuhlemeyer's boundary conditions.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.16 no.2
• /
• pp.165-172
• /
• 2003

In this study, a new and efficient harmonic axisymmetric shell element for static and dynamic analysis Is proposed. The present element considering shear strain is based on a modified mixed variational principle in which the independent unknowns are only the Quantities prescribable at the shell edges. Unlike existing hybrid-mixed axisymmetric shell elements, the present element introduces additional nodeless degrees for displacement field Interpolation In order to enhance the numerical performance. The stress parameters are eliminated by the stationary condition and the nodeless degrees are condensed out by the Guyan reduction. Through several numerical examples, the hybrid-miked shell element with the additional nodeless degrees and the consistent stress parameters is shown to be efficient and yield very accurate results for static and vibration analysis.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.28 no.5
• /
• pp.485-490
• /
• 2015

This paper presents development of a geometrically nonlinear triangular planar element including rotational degrees of freedom, based on the co-rotational(CR) formulation. The CR formulation is one of the efficient geometrically nonlinear formulations and it is based on the assumptions on small strain and large rotation. In this paper, modified CR formulation is suggested for the developemnt of a triangular planar element. The present development is validated regarding the static and time transient problems. The present results are compared with the results predicted by the previous researchers and those obtained by the existing commercial software.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.31 no.5
• /
• pp.55-62
• /
• 2003

Subdomain-interface variational formulation is presented to solve a class of dynamic contact problems of composite structures. The penalty method is used for imposing inequality constraints on contact surfaces and for connecting finite element subdomains that satisfy interface compatibility conditions. As a result, any complex-shaped domain can be easily divided into independently modeled subdomains without considering the conformity of meshes on interfaces. Some advantageous features of the present method are shown through a set a numerical studies with a developed computer code.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.12 no.2
• /
• pp.119-128
• /
• 1999

본 연구는 등방탄성체에 고속충격하중이 작용하는 경우에 대한 유한요소해석에 관한 것으로 대상구조는 여러 가지 모양의 2차원 평면구조를 택하였다. Galerkin 방법을 이용하여 유한요소 정식화하였으며 직접시간적분법에 의해 수치해를 구하였다. 본 해석에서는 균열이 없는 평판으로 수치해와 이론해를 비교하여 수치해의 신뢰성을 확인하였으며, 0°, 30°, 45°경사 균열이 없는 평판에 적용한 3가지 예를 분석하였다. 수치해석 결과는 이론해의 결과와 상호 잘 일치하였다.

• Journal of Korean Association for Spatial Structures
• /
• v.3 no.4 s.10
• /
• pp.11-17
• /
• 2003

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.13 no.1
• /
• pp.37-49
• /
• 2000

This paper is for the first essential study on the development of unified finite element formulations for solving problems related to the dynamics/thermoelastics behavior of solids. In the first part of formulations, the finite element method is based on the introduction of a new quantity defined as heat displacement, which allows the heat conduction equations to be written in a form equivalent to the equation of motion, and the equations of coupled thermoelasticity to be written in a unified form. The equations obtained are used to express a variational formulation which, together with the concept of generalized coordinates, yields a set of differential equations with the time as an independent variable. Using the Laplace transform, the resulting finite element equations are described in the transform domain. In the second, the Laplace transform is applied to both the equation of heat conduction derived in the first part and the equations of motions and their corresponding boundary conditions, which is referred to the transformed equation. Selections of interpolation functions dependent on only the space variable and an application of the weighted residual method to the coupled equation result in the necessary finite element matrices in the transformed domain. Finally, to prove the validity of two approaches, a comparison with one finite element equation and the other is made term by term.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.26 no.2A
• /
• pp.383-390
• /
• 2006

In this paper, a stochastic field that is compatible with Monte Carlo simulation is suggested for an expansion-based stochastic analysis scheme of weighted integral method. Through investigation on the way of affection of stochastic field function on the displacement vector in the series expansion scheme, it is noticed that the stochastic field adopted in the weighted integral method is not compatible with that appears in the Monte Carlo simulation. As generally recognized in the field of stochastic mechanics, the response variability is not a linear function of the coefficient of variation of stochastic field but a nonlinear function with increasing variability as the intensity of uncertainty is increased. Employing the stochastic field suggested in this study, the response variability evaluated by means of the weighted integral scheme is reproduced with high precision even for uncertain fields with moderately large coefficient of variation. Besides, despite the fact that only the first-order expansion is employed, an outstanding agreement between the results of expansion-based weighted integral method and Monte Carlo simulation is achieved.