Computer homogenization of porous piezoceramics of different ferrohardness with random porous structure and inhomogeneous polarization field

The article is concerned with the homogenization problems, in which the effective moduli of porous piezoceramic composites are determined taking into account the inhomogeneity of the polarization field. The homogenization problems are solved by the finite element method in the framework of the theory of effective moduli and the Hill energy principle using the ANSYS package. To this end, in static problems of electroelasticity, the displacements and electric potential, which are linear in spatial variables, are specified on the boundary of a representative volume to provide constant stress and electric induction fields for a homogeneous reference medium. After solving a set of boundary value problems under different boundary conditions and determining the volume-averaged stress components and the electric induction vector, a complete set of effective moduli for the piezoelectric composite is calculated. A representative volume of the piezocomposite is created in the form of a regular finite element mesh consisting of cubic elements. Pores in the representative volume are assumed to be filled with a piezoelectric material with extremely small moduli. Finite elements with pore properties are selected according to a simple random algorithm. The inhomogeneous polarization field is found by solving an electrostatic problem, in which the polarization process in the representative volume is modeled based on a simplified linear formulation. The local coordinate systems for individual finite elements of the composite matrix are specified by the directions of the polarization vectors. In the following, when solving the problems of electroelasticity, these local coordinate systems associated with the elements of the piezoelectric matrix allow recalculating the material properties according to the formulas of transformation of the tensor components as the coordinate systems rotate. In addition, consideration is given to different models describing the change in the moduli of the material from an unpolarized state to a polarized one as a function of the polarization vector. Computational experiments were carried out for three types of piezoceramics: soft ferroelectric piezoceramics PZT-5H, piezoceramics PZT-4 of medium ferrohardness, and piezoceramics PZT-8 with higher degree of ferrohardness. The dependences of the effective moduli on porosity are compared for different laws of polarization inhomogeneity and different kinds of piezoceramic material of the composite matrix.