3D mixed finite elements for curved, flat piezoelectric structures

Martin Meindlhumer,Astrid Pechstein
DOI: https://doi.org/10.1080/19475411.2018.1556186
2018-08-02
Abstract:The Tangential-Displacement Normal-Normal-Stress (TDNNS) method is a finite element method that was originally introduced for elastic solids and later extended to piezoelectric materials. It uses tangential components of the displacement and normal components of the normal stress vector as degrees of freedom for elasticity. For the electric field, the electric potential is used. The TDNNS method has been shown to provide elements which do not suffer from shear locking. Therefore thin structures (e.g. piezoelectric patch actuators) can be modeled efficiently. Hexahedral and prismatic elements of arbitrary polynomial order are provided in the current contribution. We show that these elements can be used to discretize curved, shell-like geometries by curved elements of high aspect ratio. The order of geometry approximation can be chosen independently from the polynomial order of the shape functions. We present two examples of curved geometries, a circular patch actor and a radially polarized piezoelectric semi-cylinder. Simulation results of the TDNNS method are compared to results gained in ABAQUS. We obtain good results for displacements and electric potential as well as for stresses, strains and electric field when using only one element in thickness direction.
Numerical Analysis
What problem does this paper attempt to address?
The main problem that this paper attempts to solve is to develop an effective finite - element method for simulating piezoelectric structures with complex geometries (especially thin, curved or shell - like structures) and avoid common numerical locking problems (such as shear locking). Specifically, the authors introduce and extend the Tangential - Displacement Normal - Normal - Stress (TDNNS) method to deal with the mechanical and electrical coupling problems in piezoelectric materials. The following are the key points of the paper: 1. **Avoiding Shear Locking**: - Traditional finite - element methods are prone to the shear - locking phenomenon when simulating thin structures, resulting in inaccurate calculation results. The TDNNS method effectively avoids this problem by using tangential displacement and normal stress as degrees of freedom. 2. **High - Order Geometric Approximation**: - The paper proposes a method that can accurately approximate complex geometries (such as curved piezoelectric structures) while maintaining a low number of elements. This makes it possible to model curved structures with high aspect ratios. 3. **Multi - Physical Field Coupling**: - Piezoelectric materials involve the coupling between mechanical deformation and electric fields. The TDNNS method can handle these two physical fields simultaneously and ensure the accuracy of the solution. 4. **Efficient Computation**: - By using the TDNNS method, good results of displacement, electric potential, stress, strain and electric field can be obtained even when only one element is used in the thickness direction. This greatly reduces the computational cost. 5. **Eigenvalue Problem Solving**: - The paper also discusses how to use the inverse iteration method to solve the eigenvalue problem of the coupled electromechanical system, which is very important for studying vibration characteristics. ### Formula Summary - **Stress - Strain Relationship**: \[ \sigma = C_E : \varepsilon - e^T \cdot E \] \[ D = e : \varepsilon + \epsilon_\varepsilon \cdot E \] - **Variational Formula**: \[ \int_{\Omega} (-S_E : \sigma + d^T \cdot \nabla \phi) : \delta \sigma \, d\Omega + \langle \varepsilon(u), \delta \sigma \rangle + \langle \varepsilon(\delta u), \sigma \rangle + \int_{\Omega} (d : \sigma - \epsilon_\sigma \cdot \nabla \phi) \cdot \nabla \delta \phi \, d\Omega = \int_{\Omega} f \cdot \delta u \, d\Omega \] - **Inverse Iteration Eigenvalue Problem**: \[ A \cdot q_k = \lambda_k M \cdot q_k \] \[ A \cdot q^{n + 1}_k = \lambda^n_k M \cdot q^n_k \] \[ \lambda(q) = \frac{q^T \cdot A \cdot q}{q^T \cdot M \cdot q} \] Through these improvements, the TDNNS method provides a powerful tool for the efficient and accurate simulation of piezoelectric structures.