Three dimensional curved shell finite elements for heat conduction

This paper presents a finite element formulation for three dimensional curved shell heat conduction where nodal temperatures and nodal temperature gradients through the shell thickness are retained as primary variables. The three dimensional curved shell geometry is constructed using the coordinates of the nodes lying on the middle surface of the shell and the nodal point normals. The element temperature field is defined in terms of the element approximation functions, nodal temperatures and nodal temperature gradients. The weak formulation of the three dimensional Fourier heat conduction equation is constructed in the Cartesian coordinate system. The properties of the curved shell elements are then derived using the weak formulation and the element temperature approximation. The element formulation permits linear temperature distribution through the element thickness.Distributed heat flux as well as convective boundaries are permitted on all six faces of the element. The element also has internal heat generation as well as orthotropic material capability. The superiority of the formulation in terms of applications, efficiency and accuracy is demonstrated. Numerical examples are presented and comparisons are made with theoretical solutions.     

Transition finite elements with temperature and temperature gradients as primary variables for axisymmetric heat conduction

This paper presents finite element formulation for a special class of elements referred to as “transition finite elements” for axisymmetric heat conduction. The transition elements are necessary in applications requiring the use of both axisymmetric solid elements and axisymmetric shell elements. The elements permit transition from the solid portion of the structure to the shell portion of the structure. A novel feature of the formulation presented here is that nodal temperatures as well as nodal temperature gradients are retained as primary variables. The weak formulation of the Fourier heat conduction equation is constructed in the cylindrical co-ordinate system (r, z). The element geometry is defined in terms of the co-ordinates of the nodes as well as the nodal point normals for the nodes lying on the middle surface of the element. The element temperature field is approximated in terms of element approximation functions, nodal temperatures and the nodal temperature gradients. The properties of the transition elements are then derived using the weak formulation and the element temperature approximation. The formulation presented here permits linear temperature distribution through the element thickness. Convective boundaries as well as distributed heat flux is permitted on all four faces of the element. Furthermore, the element formulation also permits distributed heat flux and orthotropic material behaviour. Numerical examples are presented, first to illustrate the accuracy of the formulation and second to demonstrate its usefulness in practical applications. Numerical results are also compared with the theoretical solutions.     

Three dimensional solid-shell transition finite elements for heat conduction

This paper presents a finite element formulation for a special class of finite elements referred to as ‘Solid-Shell Transition Finite Elements’ for three dimensional heat conduction. The solid-shell transition elements are necessary in applications requiring the use of both three dimensional solid elements and the curved shell elements. These elements permit transition from the solid portion of the structure to the shell portion of the structure. A novel feature of the formulation presented here is that nodel temperatures as well as nodal temperature gradients are retained as primary variables. The element geometry is defined in terms of coordinates of the nodes as well as the nodal point normals for the nodes lying on the middle surface of the element. The temperature field with the element is approximated in terms of element approximation functions, nodal temperatures and nodal temperature gradients. The properties of the transition element are then derived using the weak formulation (or the quadratic functional) of the Fourier heat conduction equation in the Cartesian coordinate system and the element temperature approximation. The formulation presented here permits linear temperature distribution in the element thickness direction.

Convective boundaries as well as distributed heat flux is permitted on all six faces of the elements. Furthermore, the element formulation also permits internal heat generation and orthotropic material behavior. Numerical examples are presented firstly to illustrate the accuracy of the formulation and secondly to demonstrate its usefulness in practical application. Numerical results are also compared with the theoretical solutions.     

Curved shell elements for heat conduction with p-approximation in the shell thickness direction

A finite element formulation is presented for the curved shell elements for heat conduction where the element temperature approximation in the shell thickness direction can be of an arbitrary polynomial order p. This is accomplished by introducing additional nodal variables in the element approximation corresponding to the complete Lagrange interpolating polynomials in the shell thickness direction. This family of elements has the important hierarchical property, i.e. the element properties corresponding to an approximation order p are a subset of the element properties corresponding to an approximation order p + 1. The formulation also enforces continuity or smoothness of temperature across the inter-element boundaries, i.e. C⁰ continuity is guaranteed.

The curved shell geometry is constructed using the co-ordinates of the nodes lying on the middle surface of the shell and the nodal point normals to the middle surface. The element temperature field is defined in terms of hierarchical element approximation functions, nodal temperatures and the derivatives of the nodal temperatures in the element thickness direction corresponding to the complete Lagrange interpolating polynomials. The weak formulation (or the quadratic functional) of the three-dimensional Fourier heat conduction equation is constructed in the Cartesian co-ordinate space. The element properties of the curved shell elements are then derived using the weak formulation (or the quadratic functional) and the hierarchical element approximation. The element matrices and the equivalent heat vectors (resulting from distributed heat flux, convective boundaries and internal heat generation) are all of hierarchical nature. The element formulation permits any desired order of temperature distribution through the shell thickness.

A number of numerical examples are presented to demonstrate the superiority, efficiency and accuracy of the present formulation and the results are also compared with the analytical solutions. For the first three examples, the h-approximation results are also presented for comparison purposes.     

Isoparametric line and transition finite elements with temperature and temperature gradients as primary variables for two dimensional heat conduction

This paper presents isoparametric line and transition finite element formulation for two dimensional heat conduction. The element properties are derived using weak formulation of the Fourier heat conduction equation and the element approximation where nodal temperatures and the nodal temperature gradients are retained as primary variables. The formulation permits linear temperature distribution through the element thickness. Distributed heat flux as well as convective boundaries are permitted on all four faces of the elements. Furthermore, the elements can have internal heat generation as well as orthotropic material properties. The superiority of the formulation in terms of efficiency and accuracy is demonstrated. Numerical examples are presented to illustrate their applications, and a comparison is made with theoretical results.     

p-Version axisymmetric shell element for geometrically nonlinear analysis

This paper presents a p-version geometrically nonlinear (GNL) formulation based on total Lagrangian approach for a three-node axisymmetric curved shell element. The approximation functions and the nodal variables for the element are derived directly from the Lagrange family of interpolation functions of order p_ξ and p_η. This is accomplished by first establishing one-dimensional hierarchical approximation functions and the corresponding nodal variable operators in the ξ and η directions for the three- and one-node equivalent configurations that correspond to p_ξ + 1 and p_η+ 1 equally spaced nodes in the ξ and η directions and then taking their products. The resulting element approximation functions and the nodal variables are hierarchical and the element approximation ensures C⁰ continuity. The element geometry is described by the coordinates of the nodes located on the middle surface of the element and the nodal vectors describing top and bottom surfaces of the element.

The element properties are established using the principle of virtual work and the hierarchical element approximation. In formulating the properties of the element complete axisymmetric state of stresses and strains are considered, hence the element is equally effective for very thin as well as extremely thick shells. The formulation presented here removes virtually all of the drawbacks present in the existing GNL axisymmetric shell finite element formulations and has many additional benefits. First, the currently available GNL axisymmetric shell finite element formulations are based on fixed interpolation order and thus are not hierarchical and have no mechanism for p-level change. Secondly, the element displacement approximations in the existing formulations are either based on linearized (with respect to nodal rotation) displacement fields in which case a true Lagrangian formulation is not possible and the load step size is severely limited or are based on nonlinear nodal rotation functions approach in which case though the kinematics of deformation is exact but additional complications arise due to the noncummutative nature of nonlinear nodal rotation functions. Such limitations and difficulties do not exist in the present formulation. The hierarchical displacement approximation used here does not involve traditional nodal rotations that have been used in the existing shell element formulations, thus the difficulties associated with their use are not present in this formulation.

Incremental equations of equilibrium are derived and solved using the standard Newton method. The total load is divided into increments, and for each increment of load, equilibrium iterations are performed until each component of the residuals is within a preset tolerance. Numerical examples are presented to show the accuracy, efficiency and advantages of the preset formulation. The results obtained from the present formulation are compared with those available in the literature.     

A robust non-linear solid shell element based on a mixed variational formulation

The paper is concerned with a geometrically non-linear solid shell finite element formulation, which is based on the Hu-Washizu variational principle. For the approximation of the independent displacement, stress and strain fields, the strain field is additively decomposed into two parts. Due to the fact that one part of the strain field is interpolated in the same manner as proposed by the enhanced assumed strain (EAS) method, it is denoted as EAS field. The other strain field is approximated with the same interpolation functions as the stress field. In contrast to the EAS concept the approximation spaces of the stresses and the enhanced assumed strains are not orthogonal. Consequently the stress field is not eliminated from the finite element equations. For the displacements tri-linear shape functions are considered. Shear locking and curvature thickness locking are treated using assumed natural strain interpolations. A static condensation leads to a simple low order hexahedral solid shell element. Numerical tests show that the present model is very robust and allows larger load steps than an EAS solid shell element.     

A general isoparametric finite element program SDRC SUPERB

SDRC SUPERB is a general purpose finite element program that performs linear static, dynamic and steady state heat conduction analyses of structures made of isotropic and/or orthotropic elastic materials having temperature dependent properties. The finite element library of SUPERB contains isoparametric plane stress, plane strain, flat plate, curved shell, solid type curved shell and solid elements in addition to conventional beam and spring elements. Linear, quadratic and cubic interpolation functions are available for all isoparametric elements. Independent parameters such as displacements and temperatures are obtained from SUPERB using the stiffness method of analysis. The remaining dependent parameters, such as stresses and strains, are evaluated at element gauss points and extrapolated to nodal locations. Averaged values are given as final output. The graphic capabilities of SUPERB consists of geometry and distorted geometry plotting, and stress, strain and temperature contouring. Contours are plotted at user defined cutting planes for solids and at top, middle or bottom surfaces for plate and shell types of structures.In the first part of this paper, the program capabilities of SUPERB are summarized. Extrapolation techniques used for determining dependent nodal parameters and for contour plotting are explained in the second part of the paper. Behavior of standard, wedge and transition type isoparametric elements and the effect of interpolation function orders on accuracy are discussed in the third part. The results of several illustrative problems are included.     

An asymptotic variational formulation for dynamic analysis of multilayered anisotropic plates

A multilevel variational formulation for dynamic analysis of multilayered anisotropic plates is developed within the framework of three-dimensional elasticity. By means of asymptotic expansions the Hellinger-Reissner functional for the elastodynamic problem is decomposed into a series of functionals with which a computational model can be constructed. In the formulation multiple time scales are introduced so that the secular terms can be eliminated systematically in obtaining a uniform expansion leading to valid asymptotic solution. Modifications to the approximation of various orders are determined by considering the solvability conditions of the higher-order equations. The model is adaptive, when combined with the finite element method, it has many appealing features, including that the displacements and transverse stresses may be interpolated independently, that the nodal degree-of-freedom (DOF) at each level is less than that of Kirchhoff plates, and that the mass and stiffness matrices generated at the leading-order level are always used at subsequent levels. Above all, the solution is three-dimensional in effect yet requires only two-dimensional interpolation. The through-thickness variations of the field variables are determined analytically with no need of interpolating in the thickness direction.     

Element free method for static and free vibration analysis of spatial thin shell structures

The implementation of the element free Galerkin method (EFG) for spatial thin shell structures is presented in this paper. Both static deformation and free vibration analyses are considered. The formulation of the discrete system equations starts from the governing equations of stress resultant geometrically exact theory of shear flexible shells. Moving least squares approximation is used in both the construction of shape functions based on arbitrarily distributed nodes as well as in the surface approximation of general spatial shell geometry. Discrete system equations are obtained by incorporating these interpolations into the Galerkin weak form. The formulation is verified through numerical examples of static stress analysis and frequency analysis of spatial thin shell structures. For static load analysis, essential boundary conditions are enforced through penalty method and Lagrange multipliers while boundary conditions for frequency analysis are imposed through a weak form using orthogonal transformation techniques. The EFG results compare favorably with closed-form solutions and that of finite element analyses.     

A piezoelectric solid shell element based on a mixed variational formulation for geometrically linear and nonlinear applications

The paper is focused on a piezoelectric solid shell finite element formulation. A geometrically nonlinear theory allows large deformations and includes stability problems. The formulation is based on a variational principle of the Hu–Washizu type including six independent fields: displacements, electric potential, strains, electric field, mechanical stresses and dielectric displacements. The element has eight nodes with displacements and the electric potential as nodal degrees of freedom. A bilinear distribution through the thickness of the electric field is assumed to obtain correct results in bending dominated situations. The presented element is able to model arbitrary curved shells and incorporates a 3D-material law. Numerical examples demonstrate the ability of the proposed model to analyze piezoelectric devices.     

Thermal stress analysis of laminated composite plates and shells

In this paper, Semiloof shell finite element formulation has been extended to thermal stress analysis of laminated plates and shells. The accuracy of the formulation has been verified using sample problems available in the literature. Thermal stresses in cross-ply and angle-ply laminated plates and shells subjected to thermal gradients across the thickness are presented for different boundary conditions, taking into account the temperature dependence of the material properties. The behaviour of laminates under thermal load is found to be different from that under mechanical loads in certain respects.     

A posteriori error analysis of an augmented fully-mixed formulation for the stationary Boussinesq model

In this paper we undertake an a posteriori error analysis along with its adaptive computation of a new augmented fully-mixed finite element method that we have recently proposed to numerically simulate heat driven flows in the Boussinesq approximation setting. Our approach incorporates as additional unknowns a modified pseudostress tensor field and an auxiliary vector field in the fluid and heat equations, respectively, which possibilitates the elimination of the pressure. This unknown, however, can be easily recovered by a postprocessing formula. In turn, redundant Galerkin terms are included into the weak formulation to ensure well-posedness. In this way, the resulting variational formulation is a four-field augmented scheme, whose Galerkin discretization allows a Raviart–Thomas approximation for the auxiliary unknowns and a Lagrange approximation for the velocity and the temperature. In the present work, we propose a reliable and efficient, fully-local and computable, residual-based a posteriori error estimator in two and three dimensions for the aforementioned method. Standard arguments based on duality techniques, stable Helmholtz decompositions, and well-known results from previous works, are the main underlying tools used in our methodology. Several numerical experiments illustrate the properties of the estimator and further validate the expected behavior of the associated adaptive algorithm.     

Finite element of the heat conduction equation with temperature dependent coefficients

A finite element formulation is derived for the equation of heat conduction with temperature dependent conductivity and heat capacity. The derivation of the finite element model is based on a variational formulation of the heat conduction equation which, together with generalized coordinates, yields an equation of similar form to Lagrange's equation in mechanics. The obtained equation is of special interest for applying the finite element method to solve problems with temperature-dependent properties.     

考虑刚-柔-热耦合的板结构多体系统的动力学建模

对热载荷作用下中心刚体与大变形薄板多体系统的动力学建模问题进行研究.基于Kirchhoff假设,从格林应变和曲率与绝对位移的非线性关系式出发,推导了非线性广义弹性力阵,用绝对节点坐标法建立了大变形矩形薄板的有限元离散的动力学变分方程.为了考虑刚体姿态运动、弹性变形和温度变化的相互耦合作用,推导了热流密度与绝对节点坐标之间的关系式.引入系统的运动学约束方程,建立了中心刚体-矩形板多体系统的考虑刚-柔-热耦合的热传导方程和带拉格朗日乘子的第一类拉格朗日动力学方程.为了有效地提高计算效率,将改进的中心差分法和广义-α法相结合,求解热传导方程和动力学方程,差分后的方程通过牛顿迭代法耦合求解.对刚-柔耦合和刚-柔-热三者耦合两种模型的仿真结果进行比较表明,刚体运动对温度梯度和热变形的影响显著.此外,本文建模方法考虑了几何非线性项,因此也考虑了热膨胀引起的轴向变形对横向变形的影响.     

Stability of cylindrical shell panels subjected to follower forces based on a mixed finite element formulation

A mixed finite element formulation is developed from a weak variational priniciple. This formulation is applied to stability analysis of cylindrical shell structures subjected to follower loading. Bilinear trial functions are used for all field variables. The rectangular curved elements presented here satisfy the continuity requirements for the field variables at the element interface. Two examples of a cantilevered cylindrical shell panel under different kinds of loading are solved.     

A generalized geometrically nonlinear formulation with large rotations for finite elements with rotational degrees of freedoms

A generalized geometrically nonlinear formulation using total Lagrangian approach is presented for the finite elements with translational as well as rotational degrees of freedoms. An important aspect of the formulation presented here is that the restriction on the magnitude of the nodal rotations is eliminated by retaining true nonlinear nodal rotation terms in the definition of the element displacement field and the consistent derivation of the element properties based on this displacement field. The general derivation and the formulation steps are applicable to any element with translational and rotational nodal degrees of freedoms. The specific forms of the formulation for axisymmetric shells, two-dimensional isoparametric beams, curved shells, two-dimensional transition elements and solid-shell transition elements can be easily derived by considering the explicit forms of the nonlinear nodal rotations for the element at hand. The specific forms of this formulation have already been well tested and applied to various two- and three-dimensional elements, the results for some of which are presented here. Currently it is being applied to the three-dimensional isoparametric beam elements.     

A family of special-purpose elements for analysis of ribbed and reinforced shells

A family of super-parametric special-purpose finite elements for analysis of ribbed and reinforced concrete shells is introduced. Any shell element may comprise an arbitrary number of curved ribs and/or reinforcing bars. The finite element formulation is conceived as an extension of Ahmad's thick shell element. The displacements and deformations of the ribs and/or reinforcing bars are consequently derived from the customary displacement definition of the thick shell elements. The formulation properly takes into account the excentricity of the ribs and/or reinforcing bars with respect to the middle surface of the plate or shell. Examples shown at the end of the paper illustrate the great efficiency of the concept in practical applications.     

Nonlinear shell models with seven kinematic parameters

We discuss design of nonlinear finite rotation shell model with seven kinematic displacement-like parameters, which are: three displacements of the middle surface, two rotations of the shell director, and two through-the-thickness stretching parameters. From the theoretical side we examine several possibilities for constructing the enriched kinematic field, which leads to different higher-order 7-parameter shell formulations. From the finite element implementation side a shell director interpolation is identified which eliminates the “curvature thickness locking”. Numerical examples are presented in order to compare different formulations and to illustrate the performance of the developed finite elements.