• Structural Engineering and Mechanics
• /
• v.15 no.1
• /
• pp.83-114
• /
• 2003

A new simple 3-node triangular flat-shell element with standard nodal DOF (6 DOF per node) is proposed for the linear and geometrically nonlinear analysis of very thin to thick plate and shell structures. The formulation of element GT9 (Long and Xu 1994), a generalized conforming membrane element with rigid rotational freedoms, is employed as the membrane component of the new shell element. Both one-point reduced integration scheme and a corresponding stabilization matrix are adopted for avoiding membrane locking and hourglass phenomenon. The bending component of the new element comes from a new generalized conforming Kirchhoff-Mindlin plate element TSL-T9, which is derived in this paper based on semiLoof constrains and rational shear interpolation. Thus the convergence can be guaranteed and no shear locking will happen. Furthermore, a simple hybrid procedure is suggested to improve the stress solutions, and the Updated Lagrangian formulae are also established for the geometrically nonlinear problems. Numerical results with solutions, which are solved by some other recent element models and the models in the commercial finite element software ABAQUS, are presented. They show that the proposed element, denoted as GMST18, exhibits excellent and better performance for the analysis of thin-think plates and shells in both linear and geometrically nonlinear problems.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.14 no.4
• /
• pp.481-494
• /
• 2001

An enhanced degenerated shell finite element (FE), which has been developed for inelastic analysis of reinforced concrete structures is described in this paper. Generally, Reissner-Mindlin (RM) assumptions are adopted to develop the degenerated shell FE so that transverse shear deformation effects is considered. However, it is found that there are serious defects such as locking phenomena in RM degenerated shell FE since the stiffness matrix has been overestimated in some situations. As remedies of locking phenomena, reduced integration, incompatible mode and assumed strain method have been used. Especially, the assumed strain method has been successfully used in many FEs. But contrarily, there is a few investigation on the performance of the assumed strains in the inelastic analysis of concrete structures. Therefore, shell formulation is provided in this paper with emphasis on the terms related to the stiffness matrix based on assumed strain method and microscopic concrete material model. Finally, the performance of the present shell element is tested and demonstrated with several numerical examples. From the numerical tests, the present result shows a good agreement with experimental data or other numerical results.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.8 no.2
• /
• pp.25-32
• /
• 1988

A formulation of an isoparametric quadrilateral higher-order plate bending finite element is presented. The 8-noded 28-d.o.f. plate element has been degenerated from the three-dimensional continuum by introducing the plate assumptions and considering higher-order in-plane displacement profile. The element characteristics have been derived by the Galerkin's weighted residual method and computed by using the selective reduced integration technique to avoid shear-locking phenomenon. Several numerical examples are given to demonstrate the accuracy and versatility of the proposed quadrilateral higher-order plate bending element over the other existing plate finite elements in both static and dynamic analyses.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2000.10a
• /
• pp.11-17
• /
• 2000

,An residual stress and out-of plane deformation produced by butt welding was analyzed by four kinds of 3D-FEM programs(Thermal El-P1 Analysis) developed by authors. The magnitude of deformation of perpendicular to the welding line generated by butt welding was large when the reduced integration method was used. This was because of removal of the locking phenomenon, which it was generally known that the stiffness of the shear component of out-of-plane was largely evaluated. And the magnitude of residual stress was analyzed by using the FEM program based on a large and small deformation theory was similar to that was analyzed by the redeced integration method.

• Structural Engineering and Mechanics
• /
• v.57 no.3
• /
• pp.473-506
• /
• 2016

The aim of this work is the development of a 2D quadrilateral isoparametric finite element model, based on a layerwise approach, for the bending analysis of sandwich plates. The face sheets and the core are modeled individually using, respectively, the first order shear deformation theory and the third-order plate theory. The displacement continuity condition at the interfaces 'face sheets-core' is satisfied. The assumed natural strains method is introduced to avoid an eventual shear locking phenomenon. The developed element is a four-nodded isoparametric element with fifty two degrees-of-freedom (52 DOF). Each face sheet has only two rotational DOF per node and the core has nine DOF per node: six rotational degrees and three translation components which are common for the all sandwich layers. The performance of the proposed element model is assessed by six examples, considering symmetric/unsymmetric composite sandwich plates with different aspect ratios, loadings and boundary conditions. The numerical results obtained are compared with the analytical solutions and the numerical results obtained by other authors. The results indicate that the proposed element model is promising in terms of the accuracy and the convergence speed for both thin and thick plates.

• Structural Engineering and Mechanics
• /
• v.86 no.3
• /
• pp.291-309
• /
• 2023

This paper addresses the finite element modeling of functionally graded porous (FGP) beams for free vibration and buckling behaviour cases. The formulated finite element is based on simple and efficient higher order shear deformation theory. The key feature of this formulation is that it deals with Euler-Bernoulli beam theory with only three unknowns without requiring any shear correction factor. In fact, the presented two-noded beam element has three degrees of freedom per node, and the discrete model guarantees the interelement continuity by using both C⁰ and C¹ continuities for the displacement field and its first derivative shape functions, respectively. The weak form of the governing equations is obtained from the Hamilton principle of FGP beams to generate the elementary stiffness, geometric, and mass matrices. By deploying the isoparametric coordinate system, the derived elementary matrices are computed using the Gauss quadrature rule. To overcome the shear-locking phenomenon, the reduced integration technique is used for the shear strain energy. Furthermore, the effect of porosity distribution patterns on the free vibration and buckling behaviours of porous functionally graded beams in various parameters is investigated. The obtained results extend and improve those predicted previously by alternative existing theories, in which significant parameters such as material distribution, geometrical configuration, boundary conditions, and porosity distributions are considered and discussed in detailed numerical comparisons. Determining the impacts of these parameters on natural frequencies and critical buckling loads play an essential role in the manufacturing process of such materials and their related mechanical modeling in aerospace, nuclear, civil, and other structures.

• Journal of the Korea institute for structural maintenance and inspection
• /
• v.8 no.2
• /
• pp.247-258
• /
• 2004

The composite shell element is developed for the solution of undamped forced vibration problem of composite curved panels. The finite element used in the current study is an 4-node enhanced assumed shell element with six degrees of freedom per node. The composite shell element is free of both shear and membrane locking phenomenon by using the enhanced assumed strain(EAS) method. A modification to the first-order shear deformation shell theory is proposed, which results in parabolic thorough-thickness distribution of the transverse shear strains and stresses. It eliminates the need for shear correction factors in the first order theory. Newmark's direct integration technique is used for carrying out the integration of the equation motion, to obtain the repones history. Parametric studies of curved composite panels are carried out for forced vibration analysis by geometrical shapes and by laminated composite; such as fiber orientation, stacking sequence.

• Steel and Composite Structures
• /
• v.42 no.1
• /
• pp.23-32
• /
• 2022

In this paper, a nonlinear numerical method to solve the large deflection problem is introduced. And the non-dimensional load-deflection behavior of functionally graded carbon nanotube-reinforced composite (FG-CNTRC) plates is parametrically investigated. The large deflection problem is formulated according to the von Kármán nonlinear theory and the (1,1,0)^* hierarchical model, and it is approximated by 2-D natural element method (NEM). The shear locking phenomenon is suppressed by the selectively reduced integration method. The nonlinear matrix equations are solved by combining the incremental loading scheme and the Newton-Raphson iteration method. The proposed method is validated from the benchmark experiments, where the propose method shows an excellent agreement with the reference methods. The nonlinear behavior of FG-CNTRC plates is evaluated in terms of the non-dimensional load-deflection curve, and it is parametrically investigated with respect to the existence/non-existence and gradient pattern of CNTs, the width-to-thickness and aspect ratios of plates and the type of boundary conditions. The non-dimensional central deflection is significantly reduced when CNTs and added, and it decreases with the volume fraction of CNTs. But, it shows a uniform increase in proportion to the width-to-thickness and aspect ratios. Both the gradient pattern of CNTs and the type of boundary conditions do also show the remarkable effects.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.19 no.1 s.71
• /
• pp.93-103
• /
• 2006

It is absolutely essential that safety assessment of the containment buildings during service life because containment buildings are last barrier to protect radioactive substance due to the accidents. Therefore, this study describes an enhanced degenerated shell finite element(FE) which has been developed for nonlinear FE analysis of reinforced concrete(RC) containment buildings with elasto-plastic material model. For the purpose of the material nonlinear analysis, Drucker-Prager failure criteria is adapted in compression region and material parameters which determine the shape of the failure envelop are derived from biaxial stress tests. Reissner-Mindlin(RM) assumptions are adopted to develop the degenerated shell FE so that transverse shear deformation effects is considered. However, it is found that there are serious defects such as locking phenomena in RM degenerated shell FE since the stiffness matrix has been overestimated in some situations. Therefore, shell formulation is provided in this paper with emphasis on the terms related to the stiffness matrix based on assumed strain method. Finally, the performance of the present shell element to analysis RC containment buildings is tested and demonstrated with several numerical examples. From the numerical tests, the present results show a good agreement with experimental data or other numerical results.