Shell and hexahedral curvilinear finite elements for the analysis of piezoceramics ferroelectric and ferroelastic behaviors
Développement d'éléments finis ferroélectriques et ferroélastiques de type solide et coque curvilignes
Résumé
Piezoceramics like lead zirconate titanate or PZT can produce an electric potential when they are subjected to a mechanical stress and deform in the presence of an electric field. This electromechanical coupling can be described by linear constitutive equations for moderate loadings. Nevertheless, this coupling becomes highly non linear when piezoceramics are subjected to high electromechanical loadings due to the electric polarization switching. In this thesis work, a phenomenological material constitutive model that describe the electric polarization ferroelectric switching (by an electric field) and ferroelastic switching (by a mechanical stress) is proposed. To describe the loading history, two internal variables are considered and two electric and mechanical loading surfaces are defined to indicate the onset of domain switchings. A bi-dimensional version of this model is developed to study thin piezoelectric structures. The phenomenological model 2D and 3D versions are implicitly integrated by adopting the return-mapping algorithm. Two shell and hexahedral first-order finite elements are then formulated and implemented into the commercial finite element code Abaqus via the user subroutine UEL (User ELement).
Les céramiques piézoélectriques, comme le Titatano-Zirconate de Plomb (PZT), peuvent produire une tension électrique quand elles sont soumises à une contrainte mécanique et, inversement, se déforment sous l'effet d'un champ électrique. Ce couplage électromécanique peut être décrit par des équations de comportement linéaires pour des chargements modérés. Cependant, au-delà de certaines valeurs de champ électrique ou de contrainte mécanique, ce couplage devient fortement non linéaire à cause des phénomènes de réorientation de la polarisation électrique. Dans ce travail de thèse, un modèle phénoménologique, qui tient compte des réorientations ferroélectrique (par un champ électrique) et ferroélastique (par une contrainte mécanique) de la polarisation électrique, est proposé. Deux variables internes sont considérées pour décrire l'histoire du chargement et deux surfaces de charges électrique et mécanique sont définies pour déterminer les débuts des écoulements ferroélectrique et ferroélastique. Une version bi-dimensionnelle de ce modèle est développée également pour faire l'étude des structures piézoélectriques minces. Les deux versions 2D et 3D du modèle phénoménologique sont intégrées implicitement en adoptant la méthode de retour radial (prédiction/correction). Deux éléments finis coque et hexaédrique de premier ordre, qui intègrent ce modèle phénoménologique non linéaire, sont ensuite développés et implémentés dans le code de calcul par éléments finis Abaqus via la routine utilisateur UEL (User ELement).
| Origine | Fichiers produits par l'(les) auteur(s) |
|---|