• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.11 no.3
• /
• pp.1-10
• /
• 1991

The paper is concerned with the analysis of laminated composite 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. An incremental total Lagrangian formulation is presented which allows the calculation of arbitrarily large displacements. The resulting non-linear equilibrium equations are solved by the Newton-Raphson method. The versatility and accuracy of this improved degenerated shell element are demonstrated by solving several numerical examples.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.4 no.3
• /
• pp.53-61
• /
• 1984

A nonlinear formulation of a beam finite element(NB6) on the total Lagrangian mode for the geometrically nonlinear analysis of two-dimensional elastic framed structures is presented. The NB6 beam element has been degenerated from the three-dimensional continuum by introducing the deep beam assumptions and consists of three reference nodes and three relative nodes. The element characteristics are derived by discretizing the beam equations of motion using the Galerkin weighted residual method and are reduced-integrated repeatedly for each loading step by the Newton-Raphson iteration techpique. Several numerical examples are given to demonstrate the accuracy and versatility of the proposed nonlinear NB6 beam element.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.8 no.3
• /
• pp.1-10
• /
• 1988

A higher order plate bending finite element using cubic in-plane displacement profiles is proposed for geometrically nonlinear analysis of thin and thick plates. The higher order plate bending element has been derived from the three dimensional plate-like continuum by discretization of the equations of motion by Galerkin weighted residual method, together with enforcing higher order plate assumptions. Total Lagrangian formulation has been used for geometrically nonlinear analysis of plates and consistent linearization by Newton-Raphson method has been performed to solve the nonlinear equations. The element characteristics have been computed by, selective reduced integration technique using Gauss quadrature to avoid shear locking phenomenon in case of extremely thin plates. Several numerical examples were solved with FEAP macro program to demonstrate versatility and accuracy of the present higher order plate bending element.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.13 no.2
• /
• pp.209-220
• /
• 2000

For non-linear analysis of stiffened shell structures, the total Lagrangian formulation is presented based upon the degenerated shell element. Geometrically correct formulation is developed by updating the direction of normal vectors and taking into account second order rotational terms in the incremental displacement field. Assumed strain concept is adopted in order to overcome shear locking phenomena and to eliminate spurious zero energy mode. The post-buckling behaviors of stiffened shell structures are traced by modeling the stiffener as a shell element and considering general transformation between the main structure and the stiffener at the connection node. Numerical examples to demonstrate the accuracy and the effectiveness of the proposed shell element are presented and compared with references' results.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.7 no.2
• /
• pp.21-28
• /
• 1987

The paper describes the analysis of large displacement problems including instability phenomena. The element used in this is a degenerated isoparametric shell element with eight nodes. Total Lagrangian formulation has been adopted in this study using Newton-Raphson iteration method with incremental load. The linear stability analyses performed usually for the initial position can be repeated at several advanced fundamental states on the non-linear buckling path. Thus a current estimate of the failure load is given. The numerical examples of a cylindrical panel under uniform load, simply supported plate under axial load, and clamped plate under uniform load are carried out. The examples applying degenerated isoparametric elements to bifurcation buckling and nonlinear collapse problems are also performed.

• Structural Engineering and Mechanics
• /
• v.54 no.3
• /
• pp.433-451
• /
• 2015

This paper focuses on large deflection static behavior of edge cracked simple supported beams subjected to a non-follower transversal point load at the midpoint of the beam by using the total Lagrangian Timoshenko beam element approximation. The cross section of the beam is circular. The cracked beam is modeled as an assembly of two sub-beams connected through a massless elastic rotational spring. It is known that large deflection problems are geometrically nonlinear problems. The considered highly nonlinear problem is solved considering full geometric non-linearity by using incremental displacement-based finite element method in conjunction with Newton-Raphson iteration method. There is no restriction on the magnitudes of deflections and rotations in contradistinction to von-Karman strain displacement relations of the beam. The beams considered in numerical examples are made of Aluminum. In the study, the effects of the location of crack and the depth of the crack on the non-linear static response of the beam are investigated in detail. The relationships between deflections, end rotational angles, end constraint forces, deflection configuration, Cauchy stresses of the edge-cracked beams and load rising are illustrated in detail in nonlinear case. Also, the difference between the geometrically linear and nonlinear analysis of edge-cracked beam is investigated in detail.

• Journal of the Society of Naval Architects of Korea
• /
• v.28 no.2
• /
• pp.307-317
• /
• 1991

A displacement-based finite element method is presented for the geometrically nonlinear analysis of eccentrically stiffened plates. The nonlinear degenerated shell and eccentric isobeam(isoparametric beam) elements are formulated on the basis of total Lagrangian and updated Lagrangian descriptions. To describe the stiffener's local plate buckling mode, some additional local degrees of freedom are used in the eccentric isobeam element. The eccentric isobeam element can be affectively employed to model the eccentric stiffener just like the case of the degenerated shell element. A detailed nonlinear analysis including the effects of stiffener's eccentricity is performed to estimate the critical load and the post buckling behaviour of an eccentrically stiffened plate. The critical buckling loads are found higher than analytic plate buckling load but lower than Euler buckling load which are the buckling strength requirements of classification society.

• Steel and Composite Structures
• /
• v.30 no.6
• /
• pp.493-516
• /
• 2019

This paper is dedicated to nonlinear static and free vibration analysis of Uniform Distributed Carbon Nanotube Reinforced Composite (UD-CNTRC) structures under in-plane loading. The authors have suggested an efficient six-node triangular element. Mixed Interpolation of Tensorial Components (MITC) approach is employed to alleviate the membrane locking phenomena. Moreover, the behavior of the well-known LST element is considerably improved by applying an additional linear interpolation on the strain fields. Based on the rule of mixture, the properties of CNTRC are obtained. In this study, only the uniform distributed CNTs are employed through the thickness direction of element. To achieve the natural frequencies and shape modes, the eigenvalue problem is also solved. Using Total Lagrangian Principles, large amplitude free vibration is considered based on the first normalized mode shape of structure. Different well-known plane problem benchmarks and some proposed ones are studied to validate the accuracy and capability of authors' formulations. In addition, the effects of length to the height ratio of beam, CNT's characteristics, support conditions and normalized amplitude parameter on the linear and nonlinear vibration parameters are investigated.

• Coupled systems mechanics
• /
• v.6 no.4
• /
• pp.399-415
• /
• 2017

This paper deals with the nonlinear static deflections of functionally graded (FG) porous under thermal effect. Material properties vary in both position-dependent and temperature-dependent. The considered nonlinear problem is solved by using Total Lagrangian finite element method within two-dimensional (2-D) continuum model in the Newton-Raphson iteration method. In numerical examples, the effects of material distribution, porosity parameters, temperature rising on the nonlinear large deflections of FG beams are presented and discussed with porosity effects. Also, the effects of the different porosity models on the FG beams are investigated in temperature rising.

• The Korean Journal of Rheology
• /
• v.11 no.1
• /
• pp.18-27
• /
• 1999

Predicting the deformation behaviors of sheets in thermoforming processes has been a daunting challenge due to the strong nonlinearities arising from very large deformations, mold-polymer contact condition and hyperelasticity constitutive equations. Nonlinear numerical analysis is always required to face this challenge especially for realistic processing conditions. In this study a 3-D algorithm and the membrane approximation are developed for thermoforming processes. The constitutive equation is expressed in terms of the 2nd Piola-Kirchhoff stress tensor and the Cauchy-Green deformation tensor. The 2-term Mooney-Rivlin model is used for the material model equation. The algorithm is established by the finite element formulation employing the total Lagrangian coordinate. The deformation behavior and the stress distribution results of 3-D algorithm with various point boundary conditions are compared to those of the membrane approximation algorithm. Also, the slip boundary condition and the no-slip boundary condition are applied for the systems that have molds. Finally, the effect of sheet temperatures on the final thickness distribution is investigated for the ABS material.