• Computational Structural Engineering
• /
• v.3 no.3
• /
• pp.113-123
• /
• 1990

The paper is concerned with the elasto-plastic and geometrically nonlinear analysis of shell structures using an improved degenerated shell element. In the formulation of the element stiffness, the combined use of three different techniques was made. They are; 1) an enhanced interpolation of transverse shear strains in the natural coordinate system to overcome the shear locking problem ; 2) the reduced integration technique in in-plane strains to avoid the membrane locking behavior ; and 3) selective addition of the nonconforming displacement modes to improve the element performances. This element is free of serious shear/membrane locking problems and undesirable compatible/commutable spurious kinematic deformation modes. In the formulation for plastic deformation, the concept of a layered element model is used and the material is assumed von Mises yield criterion. An incremental total Lagrangian formulation is presented which allows the calculation of arbitrarily large displacements and rotations. The resulting non-linear equilibrium equations are solved by the Netwon-Raphson method combined with load or displacement increment. The versatility and accuracy of this improved degenerated shell element are demonstrated by solving several numerical examples.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.27 no.1
• /
• pp.45-56
• /
• 1990

This study deals with the nonlinear strength analysis of laminated composite cylindrical shells to find the optimum structure of pressure vessel. By applying the F.E.M. using the 8-node degenerated Isoparametric shell element and Total Lagrangian formulation and being adopted Newton-Raphson method with incremental load as a solution scheme. the optimum structure is found from the viewpoint of minimum displacement. As a results of linear analysis on the 9 cases of laminated structure, $[50^{\circ}/-50^{\circ}]$ composition of the shell laminate give the minimum deflection. In case of the nonlinear analysis by applying Quadratic Failure Criteria on laminated combination $[{\theta}^{\circ}/-{\theta}^{\circ}]$, shell laminate structure of ${\theta}=50^{\circ}$ under external uniform pressure was founded as a optimum structure and ${\theta}=50^{\circ}$ for the case of external and axial loading combined.

• Journal of Korean Society for Quality Management
• /
• v.9 no.2
• /
• pp.31-36
• /
• 1981

This paper deals with the reliability optimization of parallel - in - series system subject to several linear constraints. The model of nonlinear constrained optimization is transformed to a saddle point problem by using Lagrange multipliers. Then Newton - Raphson method is used to solve the resulting problem and these step - by - step solution procedures are programmed in Basic Level II of micro - computer TRS-80. An example which has two linear constraints is solved and the results are analyzed.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.12
• /
• pp.3007-3019
• /
• 1993

For the same configuration, the stiffness of 6-node two dimensional isoparametric is stiffer than that of 8-node two dimensional isoparametric element. This phenomenon may be called 'Relative Stiffness Stiffening Phenomenon.' In this paper, the relative stiffness stiffening phenomenon was studied, and could be corrected by modifying the position of Gauss integral points used in the numerical integration of the stiffness matrix. For the same deformation (bending) energy of 6-node and 8-node two dimensional isoparametric elements, Gauss integral points of 6-node element have to move closer, in comparison with those of 8-node element, in the case of numerical integration along the thickness direction.

• Journal of the Society of Naval Architects of Korea
• /
• v.39 no.1
• /
• pp.28-37
• /
• 2002

In this study nonlinear dynamic behaviors of towed tow-tension cables are numerically analysed. In the case of a taut cable analysis, a bending stiffness term is usually neglected due to its minor effect but it plays an important role in a low-tension cable analysis. A low-tension cable may experience large displacements due to relatively small restoring forces and thus the effects of fluid and geometric non-linearities become predominant. The bending stiffness and non-linearity effects are considered in this work. In order to obtain dynamic behaviors of a towed low-tension cable, three-dimensional nonlinear dynamic equation is described and discretized by employing a finite difference method. An implicit method and Newton-Raphson iteration are adopted for the time integration and nonlinear solutions. For the calculation of huge size of matrices. block tri-diagonal matrix method is applied, which is much faster than the well-known Gauss-Jordan method in two point boundary value problems. Some case studies are carried out and the results of numerical simulations are compared with those of a in-house program of WHOI Cable with good agreements.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.15 no.2
• /
• pp.523-536
• /
• 1991

This dissertation analyzes the running mechanism of flexible and thin tape above rotating head through the numerical simulation and the experiment. The scope of analysis is confined to the phenomena of two dimensional elasto hydrodynamic lubrication between the protruded bump on a rotating cylinder and the running tape. This model is based on the elastic deformation equation of plate and shell and Reynolds equation. Finite difference method is employed as a numerical technique to calculate (1) the distribution of pressure between the running tape and rotating bump and (2) the vertical deformation of elastic thin tape over he rotating bump under hydrodynamic pressure. In numerical analyses, the effects of bump size on flying characteristics of the tape were evaluated and examined considering the influence of tension and stiffness of tape.

• Proceedings of the Korean Society for Bio-Environment Control Conference
• /
• 2001.04b
• /
• pp.31-33
• /
• 2001

잠열축열 그린하우스 난방 시스템의 난방특성을 분석하기 위하여 이에 대한 열평형 이론을 정립하고 수치해석에 의하여 컴퓨터 시뮬레이션 모델을 개발하고자 잠열축열 그린하우스 난방 시스템의 열저항 회로망을 구성하였다. 그리고 그린하우스의 피복재, 내부 공기, 토양표면, 잠열 축열재와의 열평형 방정식을 구성하였으며, Newton-Raphson반복법을 이용하여 수치해석을 하였고, 실험 분석을 통하여 수치해의 타당성을 검증하였다. 시뮬레이션 모델을 위하여 C언어를 사용하였으며, 겨울철 (11월-2월)의 기후 조건이 유사한 여러 날을 선정하여 온도, 태양강도, 상대습도, 토양 수분함량 등을 자료로 하여 모델링을 하였다. 여기에 사용된 토양 조건은 사양토로 건조한 상태를 유지하였다. 이상과 같은 분석에 의하여 그린하우스내 경시적 공기온도 변화와 열전달 현상의 실험치와 이론분석 결과가 잘 일치하고 있음을 알 수 있었다.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.32 no.3B
• /
• pp.175-183
• /
• 2012

In this study, Galerkin scheme and SU/PG scheme of Petrov-Galerkin family were applied to the shallow water equations and a finite element model for shallow water flow was developed. Numerical simulations were conducted in several flumes with convection-dominated flow condition. Flow simulation of channel with slender structure in the water course revealed that Galerkin and SU/PG schemes showed similar results under very low Fr number and Re number condition. However, when the Fr number increased up to 1.58, Galerkin scheme did not converge while SU/PG scheme produced stable solutions after 5 iterations by Newton-Raphson method. For the transcritical flow simulation in diverging channel, the present model predicted the hydraulic jump accurately in terms of the jump location, the depth slope, and the flow depth after jump, and the numerical results showed good agreements with the hydraulic experiments carried out by Khalifa(1980). In the oblique hydraulic jump simulation, in which convection-dominated supercritical flow (Fr=2.74) evolves, Galerkin scheme blew up just after the first iteration of the initial time step. However, SU/PG scheme captured the boundary of oblique hydraulic jump accurately without numerical oscillation. The maximum errors quantified with exact solutions were less than 0.2% in water depth and velocity calculations, and thereby SU/PG scheme predicted the oblique hydraulic jump phenomena more accurately compared with the previous studies (Levin et al., 2006; Ricchiuto et al., 2007).

• Journal of the Korean Society of Safety
• /
• v.11 no.4
• /
• pp.143-150
• /
• 1996

This is an explicit-Implicit, finite element analysis for linear as well as nonlinear hygrothermal stress problems. Additional features, such as moisture diffusion equation, crack element and virtual crack extension(VCE ) method for evaluating J-integral are implemented in this program. The Linear Elastic Fracture Mechanics(LEFM) Theory is employed to estimate the crack driving force under the transient condition for and existing crack. Pores in materials are assumed to be saturated with moisture in the liquid form at the room temperature, which may vaporize as the temperature increases. The vaporization effects on the crack driving force are also studied. The Ideal gas equation is employed to estimate the thermodynamic pressure due to vaporization at each time step after solving basic nodal values. A set of field equations governing the time dependent response of porous media are derived from balance laws based on the mixture theory Darcy's law Is assumed for the fluid flow through the porous media. Perzyna's viscoplastic model incorporating the Von-Mises yield criterion are implemented. The Green-Naghdi stress rate is used for the invariant of stress tensor under superposed rigid body motion. Isotropic elements are used for the spatial discretization and an iterative scheme based on the full newton-Raphson method is used for solving the nonlinear governing equations.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.36 no.6
• /
• pp.947-954
• /
• 2016

A simplified analysis method is first proposed in order to determine an optimized initial shape of cable-stayed bridges including all unstrained element lengths without using complicated nonlinear FE analysis. The unstrained-length based FE method is then presented using the unstrained lengths by the simplified analysis. To demonstrate validity and accuracy of the proposed method, Incheon bridge model having the fabrication camber is constructed and initial shaping analysis is performed using the presented method and commercial finite element analysis program, MIDAS. Resultantly it is shown that the initial solutions by the proposed algorithm are well optimized and in good agreement with those by MIDAS except for axial displacements of the main member.