A Piezoelectric Finite Shell Element for the Geometrically Nonlinear Analysis of Sensors

A geometrically nonlinear finite element formulation to analyze piezoelectric shell structures is presented. The formulation is based on the mixed field variational functional of Hu–Washizu. Within this variational principle the independent fields are displacements, electric potential, strains, electric field, stresses and dielectric displacements. The mixed formulation allows an interpolation of the strains and the electric field through the shell thickness, which is an essential advantage when using a three dimensional material law. It is remarked that no simplification regarding the constitutive relation is assumed. The normal zero stress condition and the normal zero dielectric displacement condition are enforced by the independent resultant stress and resultant dielectric displacement fields. The shell structure is modeled by a reference surface with a four node element. Each node possesses six mechanical degrees of freedom, three displacements and three rotations, and one electrical degree of freedom, which is the difference of the electric potential through the shell thickness. The developed mixed hybrid shell element fulfills the in–plane, bending and shear patch tests, which have been adopted for coupled field problems. A numerical investigation of a smart antenna demonstrates the applicability of the piezoelectric shell element under the consideration of geometrical nonlinearity. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A Finite Element Formulation for the Nonlinear Analysis of Piezoelectric Three-Dimensional Beam Structures

A finite element formulation for a three-dimensional piezoelectric beam which includes geometrical and material nonlinearities is presented. To account for the piezoelectric effect, the coupling between the mechanical stress and the electrical displacement is considered. Based on the Timoshenko theory, an eccentric beam formulation is introduced which provides an efficient model to analyze piezoelectric structures. The geometrically nonlinear assumption allows the calculation of large deformations including buckling analysis. A quadratic approximation of the electric potential through the cross section of the beam ensures the fulfilment of the charge conservation law exactly. This assumption leads to a finite element formulation with six mechanical and five electrical degrees of freedom per node. To take into account the typical ferroelectric hysteresis phenomena, a nonlinear material model is essential. For this purpose, the phenomenological Preisach model is implemented into the beam formulation which provides an efficient determination of the remanent part of the polarization. The applicability of the introduced beam formulation is discussed with respect to available data from literature. (© 2005 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A Mixed Finite Element Formulation for Piezoelectric Materials

A geometrically nonlinear finite element formulation to analyze multi-field problems as they arise e.g. in piezoelectric or magnetostrictive materials is presented. Here we focus on piezoelectric problems. The formulation is based on a Hu-Washizu functional considering six independent fields. These are displacements u , electric potential ϕ, strains E , stresses S , electric field , and the electric flux density . The finite element approximation leads to an 8-node hexahedral element with u and ϕ as nodal degrees of freedom. The fields E , S , , and are interpolated on element level by employing some internal degrees of freedom. These fields do not require continuity across the element boundaries, thus the internal degree of freedoms are eliminated on element level by a static condensation. The geometrically non-linear theory allows large deformations and accounts for stability problems. To fulfill the charge conservation law in bending dominated situations exactly a quadratic approximation of the electric potential is necessary. This leads in general to additional nodal degrees of freedom, which is circumvented by the presented formulation by employing appropriate interpolations of and . Numerical examples show that the locking effect which arise in low order elements are significantly reduced and that the element provides good accuracy with respect to experimental data. (© 2005 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Simulation of Piezoelectric Smart Lightweight Structures by the Finite Element Method for Acoustic Control

Piezoelectric patches are widely used as distributed actuators and sensors for active vibration control. Active control of sound in this way is mainly interesting in the low frequency range, where passive methods are less effective. To simulate such electromechanical‐acoustical systems a virtual overall model based on the FEM is presented in the paper. Different approaches for model reduction with respect to controller design purposes are discussed. As a simple test example the numerical model of an acoustic box is studied. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

壳体的唯象理论及其有限元分析方法

本文从具有不同抗拉与抗压性能的各种材料的宏观特性出发,建立了以描述不同弹性模量材料性质的壳体唯象理论,并由此给出了一种有效的分析方法以求解由抗拉与拉压性能不同的材料所制造的壳体结构强度和变形问题.     

Implementation of a Free-locking Composite Piezoelectric Shell Element

In this contribution, a free–locking shell element with piezoelectric ceramic layers is presented. The material law and the constitutive equations of the coupled electro–mechanical problem are derived in convective coordinates, so that the presented element accounts for geometrical nonlinearities [2]. To avoid zero–energy mode and locking effects of the thin shell structure, assumed natural strain (ANS) method [1] is implemented. Finally, some numerical examples demonstrate the ability of this element to solve linear and nonlinear problems. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

周期结构区域中压电问题的双尺度有限元方法

近年来纤维压电复合材料的力电性能预测已发展为一个重要的研究领域．对力电耦合周期结构的复合材料问题,通过引入匹配的边界层得到了电势与位移解的新型双尺度有限元计算方法,建立了电势与位移的双尺度耦合关系,分析了双尺度有限元解的误差.数值算例验证了方法的有效性．     

弹性结构有限元控制系统

本文讨论了经有限元方法处理后的弹性结构系统的可控、可观测、镇定等问题.所得的结论与用分布参量系统模型所得的结论一致,但却便于用计算机计算且方法简单.在一、中研究了系统的可控与可观测的问题,给出了易于用计算机判别的条件.在二、中对于采用线性反馈镇定弹性体的问题进行了仔细的讨论,指出对弹性结构系统而言,若系统完全可控仅用位移反馈可以任意配置振动频率但却无法镇定系统,而仅用速度反馈虽可以进行镇定但镇定能力是有限的,对于在系统运动方程中包含刚体运动成分的情形也作了研究.在三、中对梁的控制问题用有限元进行了处理,指出直梁作为一个系统可以分解为拉压、扭转和两个方向弯曲这四个互不关联的子系统,它们的可控与可观测问题可以分别进行讨论.最后对折线型刚架的可控与可观测的问题也作了探讨.     

A New Finite Element Space for Expanded Mixed Finite Element Method

In this article, we propose a new finite element space Λ$_h$ for the expanded mixed finite element method (EMFEM) for second-order elliptic problems to guarantee its computing capability and reduce the computation cost. The new finite element space Λ$_h$ is designed in such a way that the strong requirement V$_h\subset$Λ$_h$ in [9] is weakened to {v$_h\in$V$_h$; divv$_h$=0}$\subset$Λ$_h$ so that it needs fewer degrees of freedom than its classical counterpart. Furthermore, the new Λ$_h$ coupled with the Raviart-Thomas space satisfies the inf-sup condition, which is crucial to the computation of mixed methods for its close relation to the behavior of the smallest nonzero eigenvalue of the stiff matrix, and thus the existence, uniqueness and optimal approximate capability of the EMFEM solution are proved for rectangular partitions in $\mathbb{R}^d, d=2,3$ and for triangular partitions in $\mathbb{R}^2$. Also, the solvability of the EMFEM for triangular partition in $\mathbb{R}^3$ can be directly proved without the inf-sup condition. Numerical experiments are conducted to confirm these theoretical findings.     

摄动有限元法解一般轴对称壳几何非线性问题

本文在处理几何非线性问题时,利用在变分方程中引入振动过程,得到各级变分摄动方程,并通过有限元法求解.由于有限元法能成功地处理各种复杂边界条件、几何形状的力学问题,摄动法又可将非线性问题转化为线性问题求解.若结合这两种方法的优点,将能够解决大量复杂的非线性力学问题.并能够消除单独使用有限元法或摄动法求解复杂非线性问题所出现的困难. 本文应用摄动有限元法求解了一般轴对称壳的几何非线性问题.     

半球壳振动的有限元法分析

本文建立了用有限元法分析球壳的动力学方程.计算分析了不同边界条件下壳体的振动频率.     

有限元体系的M—P逆结构拓扑变化法

本文应用Ｍｏｏｒｅ＿Ｐｅｎｒｏｓｅ逆理论，给出了一套当结构拓扑发生变化时结构静力响应的计算公·式其特点是不需建立求解线性方程组     

Efficient Finite Element Analysis of Inelastic Structures with Iterative Solvers

This article treats the efficient solution of linear systems of equations which arise during the iterative process within the finite element analysis of inelastic structures. For the finite element analysis high order time integration methods and diagonally implicit Runge–Kutta methods (DIRK), in combination with an inexact Multilevel–Newton algorithm (MLNA) are applied. Up to 80% of the total computation time is spent by the solver for the linear systems, which suggests investigating this process. Two simple strategies to speed up the solution of the linear systems are described. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Computation of the Three-Dimensional Stress State in Composite Shell Structures with Mixed Finite Elements

Being able to compute the complete three-dimensional stress state in layered composite shell structures is essential in order to examine complicated interlaminar failure modes such as delamination. We lay out a mixed finite element formulation with independent displacements, rotations, stress resultants and shell strains. A mixed hybrid shell element with 4 nodes and 5 or 6 nodal degrees of freedom is developed, so that the element formulation can also be used for problems with shell intersections. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

一种高精度的裂纹奇异单元

基于广义伽辽金法的多变量有限元算法,增加了连续体力学有限元模型建立的灵活性.本文利用它,通过数值试验的对比建立了一种高精度的含奇异性的裂纹单元,并对多变量奇异元的构成进行了探讨.     

A Convergence Proof for Constrained Finite Element Solutions

A convergence proof is given for the finite-element solutionof the infinite dimensional quadratic programming problem ofminimizing a quadratic functional subject to linear constraints.The proof for the unconstrained problem is briefly reviewed,and then extended to the constrained case. Only the first partof the proof is given, in which necessary conditions for convergenceare derived for the specific problem and its finite-elementapproximation. The final step of proving that any problem doesobey these conditions will depend on the specific problem, butit is shown that if the finite element formulation is pointwiseconvergent and the unconstrained problem is convergent, thenso too will be the constrained problem.     

A Least-Squares Mixed Finite Element for Quasi-Incompressible Elastodynamics

In the present work, a mixed finite element based on a modified least-squares formulation is proposed. Here, we consider the time-dependent equations for quasi-incompressible elastodynamics under small strain assumptions. The main goal is to obtain an accurate approximation of both displacements and stresses in particular for the lowest-order element. Basis for the element formulation is a weak form resulting from a least-squares method. The L₂-norm minimization of the time-discretized residuals of the given first-order system leads to a functional depending on approximations for displacements and stresses. By introducing a time-independent displacement test function, a weak form is derived. A numerical example concerning quasi-incompressible elasticity shows the performance of the approach for the lowest-order element RT₀P₁as. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A Simple Finite Element Method for Meteorological Problems

A Galerkin method is applied to simple two dimensional equationsimportant in meteorological problems. The construction of thespace of trial functions for the Galerkin method is done usingthe "finite element" method, where the functions are definedas polynomials on individual elements and values are matchedon element boundaries. This method is applied to passive advectionproblems and to a non-linear gravity wave problem. The resultsare compared with those obtained by finite difference methodsand the computation time for given accuracy is shown to be atleast as short using the finite element method as with finitedifferences. Sharp local gradients are especially well handled.Extension of this approach to irregular grids and the possibleuse of higher order polynomials are proposed.     

A Finite Element Algorithm for Electromechanical Contact Problems

In the last hundred years a lot of work is done in describing and measuring the influence of the pressure on the resistance and wear in electrical contacts. But up to now there exists a lack of knowledge in predicting and optimizing the behavior of electrical contacts with numerical simulation tools considering the pressure dependency. The present work provides a new constitutive model for the contact interface in the case of current flow and a new friction law including electrical wear phenomena. Additionally numerical investigations are made to compare the numerical results with experimental data. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

线弹性结构非平稳随机振动分析的有限元方法

当前结构分析的有效方法是有限单元法,对于结构动力学问题,将变位、应力等物理量通过Fou-rier变换进行谱分解,在谱分解的形式下推求动力刚度矩阵,这样所得的矩阵和有关方程不能用结构的随机振动问题常用的振型分解法求解.本文提出了一个普遍化的求解方法.文中考虑如地震、风震等外载是如下非平稳随机过程:P(t)={P_i(t)},P_i(t)=α_i(t)P_i0(t),α_i(t)是巳知的时间函数,P_i0(t)是平稳随机过程.本文将有限单元法所得的离散化方程进行Fourier变换,利用随机过程谱分解的正交增量性质推导了激励谱和反应谱之间关系的公式.用这些公式可以寻求反应的互功率谱密度矩阵,再根据反应的统计量进行结构的安全度分析.在本文提出的计算方法中,当α_i(t)=1(i=1.,2,…,n)时方法可以简化为求解平稳过程的特殊情况.在实际应用中可以根据地震、风震记录所得的功率谱密度矩阵,按本文方法用计算机对高层、高耸、大跨度等结构问题进行分析,为了说明计算方法的特点,文中首先考虑单自由度情况,其次考虑多自由度情况,列出几个重要统计量的计算公式,并对数值计算方法和安全度分析作了讨论.