Dual-mixed variational formulation and hp finite element method for axisymmetric shell problems in elastodynamics

A new dimensionally reduced axisymmetric shell model is presented briefly for modeling time-dependent problems. This is based on the extended version of the three-field dual-mixed variational formulation of elastostatics [1, 2] to linear elastodynamics, the independent fields of which are the non-symmetric stress tensor, the displacement- and the rotation vector. An important property of the related shell model is that the classical kinematical hypotheses regarding the deformation of the normal to the shell middle surface are not used, i.e., unmodified three-dimensional constitutive equations are applied. The computational performance of the new h- and p-version axisymmetric shell finite elements is tested through a representative cylindrical shell problems. The development presented in this paper has been motivated by the fact that efficient dual-mixed hp plate and shell finite elements were managed previously to be developed for elastostatics by [1-5]. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

The stressed state of a cylindrical shell with a circular neck loaded by internal pressure

We study the stressed state of a thin cylindrical shell with a circular neck loaded by a constant internal pressure. We start from the equations obtained using Reissner's variational principle applying expansions of the components of the stress tensor and the displacement vector with respect to the coordinates along the normal to the base surface of the shell. Translated fromMatematichni Metody i Fiziko-Mekhanichni Polya, Vol. 38, 1995.     

Stabilized dual-mixed method for the problem of linear elasticity with mixed boundary conditions

We extend the applicability of the augmented dual-mixed method introduced recently in Gatica (2007), Gatica et al. (2009) to the problem of linear elasticity with mixed boundary conditions. The method is based on the Hellinger–Reissner principle and the symmetry of the stress tensor is imposed in a weak sense. The Neumann boundary condition is prescribed in the finite element space. Then, suitable Galerkin least-squares type terms are added in order to obtain an augmented variational formulation which is coercive in the whole space. This allows to use any finite element subspaces to approximate the displacement, the Cauchy stress tensor and the rotation.     

厚壳理论及其在圆柱壳中的应用

本文从Hellinger-Reissner广义变分原理出发,以位移和应力的假设为基础,建立了厚壳理论.文中把壳体的位移展开为其厚度方向的幂级数,对平行和垂直于中面的位移分别保留其级数的前四项和前三项.并假定壳体的法向挤压和横向剪切应力沿壳厚为三次曲线,使其满足上下壳面上的应力条件,利用变分原理推导出分析厚壳所需的物理方程,平衡方程和边界条件.文中对圆柱壳的情况作了实例计算,并作了光弹性实验,结果表明理论和实验符合良好.     

A Priori and a Posteriori Error Analysis of a Wavelet-Based Stabilization for the Mixed Finite Element Method

We use Galerkin least-squares terms and biorthogonal wavelet bases to develop a new stabilized dual-mixed finite element method for second-order elliptic equations in divergence form with Neumann boundary conditions. The approach introduces the trace of the solution on the boundary as a new unknown that acts also as a Lagrange multiplier. We show that the resulting stabilized dual-mixed variational formulation and the associated discrete scheme defined with Raviart–Thomas spaces are well-posed and derive the usual a priori error estimates and the corresponding rate of convergence. Furthermore, a reliable and efficient residual-based a posteriori error estimator and a reliable and quasi-efficient one are provided.     

A solid shell finite element formulation to simulate the behaviour of thin dielectric elastomer structures

To analyse the behaviour of thin structures of dielectric elastomer (DE) material a solid shell finite element is presented. The main characteristics of DEs are a non-linear hyper elastic behaviour, the quasi-incompressibility, and the ability to transform electric energy into mechanical work. Applying a voltage to thin DE structures may produce large elongation strains of 120-380%. These large strains, the efficient electro-mechanical coupling, and the light weight make DEs very attractive for the usage in actuators. Thus, there is a need for detailed research. With respect to the electro-mechanical coupling a constitutive model is presented. An electric stress tensor and a total stress tensor are introduced by considering the electrical body force and couple in the balance of linear momentum and angular momentum, respectively. The governing equations are derived and embedded in the solid shell formulation. The element formulation is based on a Hu-Washizu mixed variational principle using six independent fields: displacements, electric potential, strains, electric field, mechanical stresses, and dielectric displacements. It allows large deformations and accounts for physical nonlinearities to capture two of the main characteristics of DEs. The shell element could be applied for the modelling of arbitrary curved thin structures. The ability of the present element formulation is demonstrated in several examples. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

屈服条件由应力张量的二次齐次函数与一次齐次函数之和来表达的极限分析定理

本文建立了由应力张量σ_ij的二次齐次函数与一次齐次函数的和来表达其屈服条件的刚理想塑性体的极限分析变分原理,它可用于岩土力学的极限分析问题,并把屈服条件为应力张量σ_ij 的二次齐次函数或一次齐次函数来表达的情况作为其特例.     

A mixed finite element formulation for piezoelectric shells

A finite element formulation to analyze piezoelectric shell problems is presented. A reference surface of the shell is modelled with a four node element. Each node possesses six mechanical degrees of freedom, three displacements and three rotations, and one electric degree of freedom, which is the difference of the electric potential in thickness direction. The formulation is based on the mixed field variational principle of Hu-Washizu. The independent fields are displacements u , electric potential φ, strains E , electric field E , stresses S and dielectric displacements D . The mixed formulation allows an interpolation of the strains and the electric field in thickness direction. Accordingly a three-dimensional material law is incorporated in the variational formulation. It is remarked that no simplification regarding the constitutive law is assumed. The formulation allows the consideration of arbitrary constitutive relations. The normal zero stress condition and the normal zero dielectric displacement condition are enforced by the independent stress and dielectric displacement fields. They are defined as zero in thickness direction. The present shell element fulfills the important patch tests: the in-plane, bending and shear test. Some numerical examples demonstrate the applicability of the present piezoelectric shell element. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A dual-dual mixed formulation for nonlinear exterior transmission problems

We combine a dual-mixed finite element method with a Dirichlet-to-Neumann mapping (derived by the boundary integral equation method) to study the solvability and Galerkin approximations of a class of exterior nonlinear transmission problems in the plane. As a model problem, we consider a nonlinear elliptic equation in divergence form coupled with the Laplace equation in an unbounded region of the plane. Our combined approach leads to what we call a dual-dual mixed variational formulation since the main operator involved has itself a dual-type structure. We establish existence and uniqueness of solution for the continuous and discrete formulations, and provide the corresponding error analysis by using Raviart-Thomas elements. The main tool of our analysis is given by a generalization of the usual Babuska-Brezzi theory to a class of nonlinear variational problems with constraints.

     

只含两个独立变量的扁壳大挠度问题修正的海林格—赖斯内变分泛函

本文首先用海林格-赖斯内变分原理建立任意形状扁壳大挠度问题的泛函，然后用修正的变分原理导出适合于有限单元法的变分泛函表达式．泛函中只包含应力函数F和挠度W两个独立交量．其中也导出了在边界上用上述两个变量表示的中面位移的表达式．推导中考虑了边界的曲率，所以适用于任意形状的边界．     

Dynamic stability of a porous cylindrical shell

This paper is devoted to a closed cylindrical shell made of a porous-cellular material. The mechanical properties vary continuously on the thickness of a shell. The mechanical model of porosity is as described as presented by Magnucki, Stasiewicz. A shell is simply supported on edges. On the ground of assumed displacement functions the deformation of shell is defined. The displacement field of any cross section and linear geometrical and physical relationships are assumed in cylindrical coordinate system. The components of deformation and stress state were found. Using the Hamilton's principle the system of differential equations of dynamic stability is obtained. The forms of unknown functions are assumed and the system of a differential equations is reduced to a simple ordinary equation of dynamic stability of shell (Mathieu's equation). The derived equation are used for solving a problem of dynamic stability of porous-cellular shell with intensity of load directed in generators of shell. The critical loads are derived for a family of porous shells. The unstable space of family porous shells is found. The influence a coefficient of porosity on the stability regions in Figures is presented. The results obtained for porous shell are compared to a homogeneous isotropic cylindrical shell. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Coupled multiscale finite element analysis of shell structures

In this paper, a coupled two-scale shell model is presented. A variational formulation and associated linearisation for the coupled global-local boundary value problem is derived. The discretisation of the shell is performed with quadrilaterals, whereas the local boundary value problems at the integration points of the shell are discretised using 8-noded or 27-noded brick elements, or solid shell elements. The coupled boundary value problem is simultaneously solved within a Newton iteration scheme. Solutions for small strain problems are computed within the so-called FE² method. In an important test, the correct material matrix for the stress resultants assuming linear elasticity and a homogeneous continuum is verified. Examples show that the developed two-scale model is able to analyse the global and local mechanical behaviour of heterogeneous shell structures. (© 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Dual-mixed finite element approximation of Stokes and nonlinear Stokes problems using trace-free velocity gradients

In this work a finite element method for a dual-mixed approximation of Stokes and nonlinear Stokes problems is studied. The dual-mixed structure, which yields a twofold saddle point problem, arises in a formulation of this problem through the introduction of unknown variables with relevant physical meaning. The method approximates the velocity, its gradient, and the total stress tensor, but avoids the explicit computation of the pressure, which can be recovered through a simple postprocessing technique. This method improves an existing approach for these problems and uses Raviart-Thomas elements and discontinuous piecewise polynomials for approximating the unknowns. Existence, uniqueness, and error results for the method are given, and numerical experiments that exhibit the reduced computational cost of this approach are presented.     

Estimates of the deviations from the exact solutions for variational inequalities describing the stationary flow of certain viscous incompressible fluids

This paper is concerned with computable and guaranteed upper bounds of the difference between exact solutions of variational inequalities arising in the theory of viscous fluids and arbitrary approximations in the corresponding energy space. Such estimates (also called error majorants of functional type) have been derived for the considered class of nonlinear boundary‐value problems in (Math. Meth. Appl. Sci. 2006; 29:2225–2244) with the help of variational methods based on duality theory from convex analysis. In the present paper, it is shown that error majorants can be derived in a different way by certain transformations of the variational inequalities that define generalized solutions. The error bounds derived by this techniques for the velocity function differ from those obtained by the variational method. These estimates involve only global constants coming from Korn‐ and Friedrichs‐type inequalities, which are not difficult to evaluate in case of Dirichlet boundary conditions. For the case of mixed boundary conditions, we also derive another form of the estimate that contains only one constant coming from the following assertion: the L² norm of a vector‐valued function from H¹(Ω) in the factor space generated by the equivalence with respect to rigid motions is bounded by the L² norm of the symmetric part of the gradient tensor. As for some ‘simple’ domains such as squares or cubes, the constants in this inequality can be found analytically (or numerically), we obtain a unified form of an error majorant for any domain that admits a decomposition into such subdomains. Copyright © 2009 John Wiley & Sons, Ltd.     

A simplified approach to estimating the stress state of flexible cylindrical shells under nonaxisymmetric load

An approach is proposed to approximate numerical analysis of the stress state of a geometrically nonlinear cylindrical shell under an unsymmetric load. An approximate technique is also proposed to reduce the solution of a two-dimensional nonlinear problem of a cylindrical shell with unsymmetric parameters to the solution of axisymmetric problems. A numerical example is considered. Translated from Vychislitel'naya i Prikladnaya Matematika, No. 58, pp. 61–66, 1986.     

A method for the approximate solution of quasi-static problems for hardening bodies

A method for the approximate solution of quasi-static problems for hardening elastoplastic bodies is proposed. The constitutive relation of the model is taken in the form of a variational inequality. An approximate solution of the initial problem is constructed in time steps and, by means of the finite element method, is reduced to the solution of a system of two variational inequalities in corresponding finite-dimensional space. It is shown that the solution of this system is equivalent to finding the saddle point of the corresponding quadratic functional. To find the saddle point, Udzawa's algorithm is used, by means of which the process of finding the velocity vector and stress tensor reduces to the successive calculation of these quantities: the velocity vector is determined from the variational inequality corresponding to the equilibrium equation, and the stress tensor is determined from the variational inequality corresponding to the constitutive relation. The latter inequality is reduced to a certain non-linear equation containing the operation of projection onto a closed convex set corresponding to the elastic strains of the medium. In turn, the solution of the non-linear equation is constructed using the method of successive approximations. To illustrate the use of the proposed method, the one-dimensional problem of the quasi-static deformation of a cylindrical tube under a load applied to its internal surface is considered.     

A variational method of deriving the equations of the non-linear mechanics of liquid crystals

The non-linear equations of the dynamics of liquid crystals [1], derived previously by the Poisson brackets method, are derived from the Hamilton-Ostrogradskii variational principle. The variational problem of an unconditional extremum of the action functional in Lagrange variables is investigated. The difference between the volume densities of the kinetic and free energy of the liquid crystal is used as the Lagrangian. It is shown that the variational equations obtained are equivalent to the differential laws of conservation of momentum and the kinetic moment of the liquid crystal in Euler variables, while the Ericksen stress tensor and the molecular field are defined in terms of the derivatives of the free energy.     

叠层圆柱厚壳混合状态方程的弱形式及其边值问题

从Helinger-Reisner变分原理出发,在柱坐标系中,导出圆柱壳的弱形式混合状态方程和边界条件,联用状态空间法给出叠层柱壳的解析解,此法使得求解该类问题的形式得以扩大和统一。     

A novel numerical solution strategy for solving nonlinear free and forced vibration problems of cylindrical shells

Nonlinear vibration analysis of circular cylindrical shells has received considerable attention from researchers for many decades. Analytical approaches developed to solve such problem, even not involved simplifying assumptions, are still far from sufficiency, and an efficient numerical scheme capable of solving the problem is worthy of development. The present article aims at devising a novel numerical solution strategy to describe the nonlinear free and forced vibrations of cylindrical shells. For this purpose, the energy functional of the structure is derived based on the first-order shear deformation theory and the von–Kármán geometric nonlinearity. The governing equations are discretized employing the generalized differential quadrature (GDQ) method and periodic differential operators along axial and circumferential directions, respectively. Then, based on Hamilton's principle and by the use of variational differential quadrature (VDQ) method, the discretized nonlinear governing equations are obtained. Finally, a time periodic discretization is performed and the frequency response of the cylindrical shell with different boundary conditions is determined by applying the pseudo-arc length continuation method. After revealing the efficiency and accuracy of the proposed numerical approach, comprehensive results are presented to study the influences of the model parameters such as thickness-to-radius, length-to-radius ratios and boundary conditions on the nonlinear vibration behavior of the cylindrical shells. The results indicate that variation of fundamental vibrational mode shape significantly affects frequency response curves of cylindrical shells.     

A nonlinear surface-stress-dependent model for vibration analysis of cylindrical nanoscale shells conveying fluid

A nonlinear surface-stress-dependent nanoscale shell model is developed on the base of the classical shell theory incorporating the surface stress elasticity. Nonlinear free vibrations of circular cylindrical nanoshells conveying fluid are studied in the framework of the proposed model. In order to describe the large-amplitude motion, the von Kármán nonlinear geometrical relations are taken into account. The governing equations are derived by using Hamilton's principle. Then, the method of multiple scales is adopted to perform an approximately analytical analysis on the present problem. Results show that the surface stress can influence the vibration characteristics of fluid-conveying thin-walled nanoshells. This influence becomes more and more considerable with the decrease of the wall thickness of the nanoshells. Furthermore, the fluid speed, the fluid mass density, the initial surface tension and the nanoshell geometry play important roles on the nonlinear vibration characteristics of fluid-conveying nanoshells.