• Structural Engineering and Mechanics
• /
• v.40 no.2
• /
• pp.191-213
• /
• 2011

This paper presents an improved 8-node shell element for the analysis of plates and shells. The finite element, based on a refined first-order shear deformation theory, is further improved by the combined use of assumed natural strain method. We analyze the influence of the shell element with the different patterns of sampling points for interpolating different components of strains. Using the assumed natural strain method with proper interpolation functions, the present shell element generates neither membrane nor shear locking behavior even when full integration is used in the formulation. Further, a refined first-order shear deformation theory, which results in parabolic through-thickness distribution of the transverse shear strains from the formulation based on the third-order shear deformation theory, is proposed. This formulation eliminates the need for shear correction factors in the first-order theory. Numerical examples demonstrate that the present element perform better in comparison with other shell elements.

• East Asian mathematical journal
• /
• v.34 no.5
• /
• pp.601-614
• /
• 2018

In this paper, we introduce an extrapolated higher order characteristic finite element method to approximate solutions of nonlinear Sobolev equations with a convection term and we establish the higher order of convergence in the temporal and the spatial directions with respect to $L^2$ norm.

• Structural Engineering and Mechanics
• /
• v.3 no.4
• /
• pp.383-390
• /
• 1995

A 2-D four-noded finite element which contains a ${\lambda}$ singularity is developed. The new element is compatible with quadratic standard isoparametric elements. The element is tested on two different examples. In the first example, an edge crack problem is analyzed using two different meshes and different integration orders. The second example is a crack perpendicular to the interface problem which is solved for different material properties and in turn different singularity order ${\lambda}$. The results of those examples illustrate the efficiency of the proposed element.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.4 no.1
• /
• pp.67-77
• /
• 2000

In this paper we consider finite element mathods for nonlinear parabolic problems defined in ${\Omega}{\subset}R^d$ ($d{\leq}4$). A new initial approximation is taken. Optimal order error estimates in $L_p$ for $2{\leq}p{\leq}{\infty}$ are established for arbitrary order finite element. One order superconvergence in $W^{1,p}$ for $2{\leq}q{\leq}{\infty}$ are demonstrated as well.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2003.04a
• /
• pp.261-268
• /
• 2003

A first order shear deformable Loop-subdivision triangular element which can handle transverse shear deformation of moderately thick shell is developed. The developed element is general since it includes the effect of transverse shear deformation and has standard six degrees of freedom per node.(three translations and three rotations) The quartic box-spline function is employed as interpolation basis function. Numerical examples for the benchmark problems are analyzed in order to assess the performance of the newly developed subdivision shell element. Both in the uniform and in the distorted mesh configurations.

• Structural Engineering and Mechanics
• /
• v.65 no.3
• /
• pp.213-222
• /
• 2018

The purpose of this paper is to study free vibration analysis of thick plates resting on Winkler foundation using Mindlin's theory with shear locking free fourth order finite element, to determine the effects of the thickness/span ratio, the aspect ratio, subgrade reaction modulus and the boundary conditions on the frequency paramerets of thick plates subjected to free vibration. In the analysis, finite element method is used for spatial integration. Finite element formulation of the equations of the thick plate theory is derived by using higher order displacement shape functions. A computer program using finite element method is coded in C++ to analyze the plates free, clamped or simply supported along all four edges. In the analysis, 17-noded finite element is used. Graphs are presented that should help engineers in the design of thick plates subjected to earthquake excitations. It is concluded that 17-noded finite element can be effectively used in the free vibration analysis of thick plates. It is also concluded that, in general, the changes in the thickness/span ratio are more effective on the maximum responses considered in this study than the changes in the aspect ratio.

• Proceedings of the Korean Society For Composite Materials Conference
• /
• 2004.10a
• /
• pp.69-72
• /
• 2004

A new three-node triangular shell element based on higher order zig-zag theory is developed for laminated composite shells with multiple delaminations. The present higher order zig-zag shell theory is described in a general curvilinear coordinate system and in general tensor notation. All the complicated curvatures of surface including twisting curvatures can be described in an exact manner in the present shell element because this element is based on geometrically exact surface representation. The displacement field of the proposed finite element includes slope of deflection, which requires continuity between element interfaces. Thus the nonconforming shape function of Specht's three-node triangular plate bending element is employed to interpolate out-of-plane displacement. The present element passes the bending and twisting patch tests in flat surface configurations. The developed element is evaluated through the eigenvalue problems of composite cylindrical shells with multiple delaminations. Through the numerical examples it is demonstrated that the proposed shell element is efficient because it has minimal degrees of freedom per node. The present shell element should serve as a powerful tool in the prediction of natural frequency and modes of multi-layered thick laminated shell structures with arbitrary-shaped multiple delaminations.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2008.04a
• /
• pp.20-25
• /
• 2008

The finite element analysis for porous media is severe job because constituents have different physical peoperties, and element's continuity and stability should be considered. Thus, we propose the new mixed finite element method in order to overcome the problems. In this method, multi time step, remeshing step, and sub iteration step are introduced. The multi time step and remeshing step make it possible to satisfy a stability and an accuracy during sub iteration in which global time is determined. Finally, the proposed method is compared with the ABAQUS(2007) software and exact solution(Schiffman 1967) through two dimensional consolidation model.

• Structural Engineering and Mechanics
• /
• v.17 no.6
• /
• pp.735-749
• /
• 2004

This work presents a direct integration scheme, based on a fourth order finite difference approach, for elastodynamics. The proposed scheme was chosen as an alternative for attenuating the errors due to the use of the central difference method, mainly when the time-step length approaches the critical time-step. In addition to eliminating the spurious numerical oscillations, the fourth order finite difference scheme keeps the advantages of the central difference method: reduced computer storage and no requirement of factorisation of the effective stiffness matrix in the step-by-step solution. A study concerning the stability of the fourth order finite difference scheme is presented. The Finite Element Method and the Boundary Element Method are employed to solve elastodynamic problems. In order to verify the accuracy of the proposed scheme, two examples are presented and discussed at the end of this work.

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