Numerical modelling of piezoelectric actuators exposed to hydrogen

Modern fuel injectors have been developed based on piezoelectric stack actuators. Performance and durability of actuators in a hydrogen environment are important considerations in the development of hydrogen injectors. 2D plane stress and 3D models for analysis of coupled diffusion and thermo-electromechanical response of actuators are presented. Chemical potential, electric field and temperature gradients are taken as driving forces for hydrogen transport. The explicit Euler finite difference method is used to solve the nonlinear diffusion governing equation. The finite element method is used for time-dependent analysis of fully coupled mechanical, electric and thermal fields. The diffusion process and thermo-electromechanical deformations are coupled through the dependence of piezoelectric properties on hydrogen concentration. Experimental results for the piezoelectric coefficient d₃₃ of PZT ceramics exposed to different hydrogen concentrations are used. A comparison of a fully coupled 2D model with 2D and 3D models with reduced coupling is made to examine the significance of coupling and computational efficiency. Selected numerical results are presented for time histories of hydrogen concentration, temperature and stroke of an idealized actuator unit cell to obtain a preliminary understanding of the performance of actuators exposed to hydrogen.