• Structural Engineering and Mechanics
• /
• v.17 no.6
• /
• pp.751-766
• /
• 2004

This paper extends the use of the hierarchic degenerated shell element to geometric non-linear analysis of composite laminated skew plates by the p-version of the finite element method. For the geometric non-linear analysis, the total Lagrangian formulation is adopted with moderately large displacement and small strain being accounted for in the sense of von Karman hypothesis. The present model is based on equivalent-single layer laminate theory with the first order shear deformation including a shear correction factor of 5/6. The integrals of Legendre polynomials are used for shape functions with p-level varying from 1 to 10. A wide variety of linear and non-linear results obtained by the p-version finite element model are presented for the laminated skew plates as well as laminated square plates. A numerical analysis is made to illustrate the influence of the geometric non-linear effect on the transverse deflections and the stresses with respect to width/depth ratio (a/h), skew angle (${\beta}$), and stacking sequence of layers. The present results are in good agreement with the results in literatures.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.15 no.3
• /
• pp.491-499
• /
• 2002

A p-version finite element model based on degenerate shell element is proposed tot the analysis of orthotropic laminated plates. In the nonlinear formulation of the model, the total Lagrangian formulation is adopted with large deflection and moderate rotation being accounted tot in the sense of yon Karman hypothesis. The material model is based on the Huber-Mises yield criterion and Prandtl-Reuss flow rule in accordance with the theory of strain hardening yield function, which is generalized lot anisotropic materials by introducing the parameters of anisotropy. The model is also based on extension of equivalent-single layer laminate theory(ESL theory) with shear deformation, leading to continuous shear strain at the interface of two layers. The integrals of Legendre polynomials are used for shape functions with p-level varying from 1 to 10. Gauss-Lobatto numerical quadrature is used to calculate the stresses at the nodal points instead of Gauss points. The validity of the proposed P-version finite element model is demonstrated through several comparative points of iew in terms of ultimate load, convergence characteristics, nonlinear effect, and shape of plastic tone.

• Selected Papers of The Society of Naval Architects of Korea
• /
• v.1 no.1
• /
• pp.91-100
• /
• 1993

A displacement-based finite element method Is presented for the geometrically nonlinear analysis of eccentrically stiffened plates. A nonlinear degenerated shell element and a nonlinear degenerated eccentric isoparametric beam (isobeam) element are formulated on the basis of Total Agrangian and Updated Lagrangian descriptions. In the formulation of the isobeam element, some additional local decrees of freedom are implementd to describe the stiffener's local plate buckling modes. Therefore this element can be effectively employed to model the eccentric stiffener with fewer D.O.F's than the case of a degenerated shell element. Some detailed buckling and nonlinear analyses of an eccentrically stiffened plate are performed to estimate the critical buckling loads and the post buckling behaviors including the local plate buckling of the stiffeners discretized with the degenerated shell elements and the isobeam elements. The critical buckling loads are found to be higher than the analytical plate buckling load but lower than Euler buckling load of the corresponding column, i.e, buckling strength requirements of the Classification Societies for the stiffened plates.

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

• Journal of the Earthquake Engineering Society of Korea
• /
• v.8 no.5 s.39
• /
• pp.15-24
• /
• 2004

The purpose of this study is to investigate the inelastic behavior and ductility capacity of reinforced concrete bridge piers including P-delta effects. A computer program, named RCAHEST (Reinforced Concrete Analysis in Higher Evaluation System Technology), for the analysis of reinforced concrete structures was used. Material nonlinearity is taken into account by comprising tensile, compressive and shear models of cracked concrete and a model of reinforcing steel. The smeared crack approach is incorporated. In addition to the material nonlinear properties, the algorithm for large displacement problem that may give an additional deformation has been formulated using total Lagrangian formulation. The proposed numerical method for the inelastic behavior and ductility capacity of reinforced concrete bridge piers is verified by comparison with reliable experimental results.

• Journal of the Korean Society of Fisheries and Ocean Technology
• /
• v.28 no.1
• /
• pp.79-92
• /
• 1992

This is analyzed using the finite element method which is appling excellent isoparametric curve element in the aspect of large usages of dynamic responses in which is regarding geometric and material nonlinear of a large scale shell structure of an airplane, a submarine, a ship, and an ocean structure. The solution of dynamic equations is got by direct integration method using time-stepping procedure and regarding Central Difference Method of the both solutions. But because formal matrix factorization is not necessary in each time step and it does not take less time to compute relatively, this method must be regarded very few time steps on the condition. Axisymmatric shell problems are inspected using 8 node Isoparametric element in this paper. Partial axisymmatric spherical shell is used as a model to analyze axisymmatric nonlinear dynamic behavior regarding. Total Lagrangian formulation in geometric nonlinear behavior and elastio-viscoplastic in material nonlinear behavior.

• Journal of Theoretical and Applied Mechanics
• /
• v.2 no.1
• /
• pp.31-50
• /
• 1996

In this paper, a finite element analysis based on the local approach concept to fracture in the continuum damage mechanics is performed to analyze ductile fracture in two dimensional quasi-static state. First an isotropic damage model based on the generalized concept of effective stress is proposed for structural materials in the context of large deformation. In this model, the stiffness degradation is taken as a measure of damage and so, the fracture phenomenon can be explained as the critical deterioration of stiffness at a material point. The modified Riks' continuation technique is used to solve incremental iterative equations. Crack propagation is achieved by removing critically damaged elements. The mesh size sensitivity analysis and the simulation of the well known shearing mode failure in plane strain state are carried out to verify the present formulation. As numerical examples, an edge cracked plate and the specimen with a circular hole under plane stress are taken. Load-displacement curves and successively fractured shapes are shown. From the results, it can be concluded that the proposed model based on the local approach concept in the continuum damage mechanics may be stated as a reasonable tool to explain ductile fracture initiation and crack propagation.

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

An efficient numerical procedure for material and geometric nonlinear analysis of reinforced concrete shells under monotonically increasing loads through their elastic, inelastic and ultimate load ranges is developed by using the finite element method. The 8-node Serendipity isoparametric element developed by the degeneration approach including the transverse shear deformation is used. A layered approach is used to represent the steel reinforcement and to discretize the concrete behavior through the thickness. The total Lagrangian formulation based upon the simplified Von Karman strain expressions is used to take into account the geometric nonlinearity of the structure. The material nonlinearities are taken into account by comprising the tension, compression, and shear models of cracked concrete and a model for reinforcement in the concrete; and also a so-called smeared crack model is incorporated. The steel reinforcement is assumed to be in a uniaxial stress state and is modelled as a smeared layer of equivalent thickness. This method will be verified a useful tool to account for geometric and material nonlinearities in detailed analysis of reinforced concrete concrete shells of general form through numerical examples of the sequential paper( ).

• Magazine of the Korea Concrete Institute
• /
• v.8 no.6
• /
• pp.223-234
• /
• 1996

The purpose of this paper is to present an analysis method by using the finite element method which can exactly analyze load-deflection relationships, crack propagations. and stresses and strains of reinforcements, tendons, and concrete in behaviors of elastic. inelastic and ultimate ranges of reinforced and prestressed concrete slabs under monotonically increasing loads. For t h i s purpose, the m a t e r i a l and geometric nonlinearities are taken into account in this study. The total Lagrangian formulation based upon the simplified Von Karman strain expressions is used to take into account the geometric nonlinearities of the structure. The material nonlinearities are taken into account by comprising the tension, compression. and shear models of cracked concrete and models for reinforcements and tendons in the concrete : and also a so-called smeared crack model is incorporated. The reinforcements and t,endons are assumed to be in a uniaxial stress state and are modelled as smeared layers of equivalent thickness. For the verification of application and validity of the method proposed in this paper, several numerical examples are analyzcd and compared with experimental results. As a result, this method can successfully predict the nonlinear and inelastic behaviors throughout the fracture of reinforced and prestressed concrete slabs.

• Coupled systems mechanics
• /
• v.1 no.4
• /
• pp.361-379
• /
• 2012

This paper presents a fully coupled three-dimensional solver for the analysis of interaction between pulsatile flow and large deformation structure. A partitioned time marching algorithm is employed for the solution of the time dependent coupled discretised problem, enabling the use of highly developed, robust and well-tested solvers for each field. Conservative transfer of information at the fluid-structure interface is combined with an effective multi-predict-correct iterative scheme to enable implicit coupling of the interacting fields at each time increment. The three-dimensional unsteady incompressible fluid is solved using a powerful implicit time stepping technique and an ALE formulation for moving boundaries with second-order time accurate is used. A full spectrum of total variational diminishing (TVD) schemes in unstructured grids is allowed implementation for the advection terms and finite element shape functions are used to evaluate the solution and its variation within mesh elements. A finite element dynamic analysis of the highly deformable structure is carried out with a numerical strategy combining the implicit Newmark time integration algorithm with a Newton-Raphson second-order optimisation method. The proposed model is used to predict the wave flow fields of a particular flow-induced vibrational phenomenon, and comparison of the numerical results with available experimental data validates the methodology and assesses its accuracy. Another test case about three-dimensional biomedical model with pulsatile inflow is presented to benchmark the algorithm and to demonstrate the potential applications of this method.