Numerical simulation of fluid-saturated piezoelectric porous media¶
We consider a porous piezoelectric medium which consists of a piezoelectric matrix, embedded metallic electrodes (conductors) and fluid-filled inclusions. These components are arranged in a periodic lattice so that the medium can be generated by copies of the reference unit cell, see Fig. 1. The mechanical behavior of such a structure can be described using the two-scale asymptotic homogenization method, see [RohanLukes2018].
The mechanical properties of the piezoelectric solid are given by the following constitutive equations
where is the Cauchy stress tensor, is the electric displacement, is the strain tensor, is the displacement field, is the electric field and is the electric potential. On the right hand side of (1), we have the fourth-order elastic tensor (), the third-order tensor (), which couples mechanical and electric quantities, and the permeability tensor . The superscript denotes the quantities oscillating within the heterogeneous structure with the period equal to the size of the periodic unit.
The quasi-static problem of the piezoelectric medium is given by the following equilibrium equations
end by the mass conservation equation for the k th fluid inclusion occupying domain
where is the inclusion pressure and is the fluid compressibility.
We apply the standard homogenization techniques to the above problem. It results in the limit model for , where is the scale parameter relating the microscopic and macroscopic length scales. The homogenization process leads to local microscopic problems, defined within a reference periodic cell, and to the global problem describing the behavior of the homogenized medium at the macroscopic level. The global problem involves the homogenized material coefficients which are evaluated using the solutions of the local problems. Due to linearity of the problem, the microscopic and macroscopic problems are decoupled.
As we assume given potentials in each of the electrode networks, the dielectric properties must be appropriately rescaled in order to preserve the finite electric field in the limit:
The local microscopic responses of the piezoelectric structure are given by the following sub-problems which are solved within the periodic reference cell , see Fig. 1, that is decomposed similarly to the decomposition of domain :
Find , such that for all , and for any
Find , such that for all , satisfying
Find , such that for all , satisfying
Find , such that for all , and for any ( is the number of conductors)
The microscopic sub-problems are solved with the periodic boundary conditions and on for , is the interface between the matrix part and -th conductor . Functions , and are equal to zero on .
With the characteristic responses obtained by solving local sub-problems, the homogenized material coefficients , , , and can be evaluated using the following expressions:
The global macroscopic problem is defined in terms of the homogenized coefficients as: Find the macroscopic displacements and such that for all and
We assumed that the volume electric charge is equal to zero, otherwise we would need extra coefficients and right hand side terms.
To run the numerical simulation, download the archive, unpack it in the main SfePy directory and type:
This invoke the
simply.py script which calculates the macroscopic
problem (10) and calls the homogenization engine that solves the
local subproblems (5)–(6), evaluates the homogenized
coefficients (9) and finally performs the reconstruction of the
solutions at the microscopic level. See [CimrmanLukesRohan2019] for more details related
to the SfePy homogenization engine.
The macroscopic sample is fixed on its left side, so that no displacements are allowed, see Fig. 2 left. The defromation is induced due to piezoelectric effect, as the responce to the prescribed electric potentials , , see (10).
The deformed shape of the sample, pressure field and the magnitude of macroscopic strain are depicted in Fig. 3. The reconstructed strain and electric fields for a given in the part of the macroscopic domain are shown in Fig. 4.
Rohan E., Lukeš V. Homogenization of the fluid-saturated piezoelectric porous media. International Journal of Solids and Structures, 147:110-120, 2018, DOI:10.1016/j.ijsolstr.2018.05.017
Cimrman R., Lukeš V., Rohan E. Multiscale finite element calculations in Python using SfePy. Advances in Computational Mathematics, 45(4):1897-1921, 2019, DOI:10.1007/s10444-019-09666-0