Next Article in Journal
A Study on the Optimal Machining Parameters of the Induction Assisted Milling with Inconel 718
Next Article in Special Issue
Parallel Finite Element Model for Multispecies Transport in Nonsaturated Concrete Structures
Previous Article in Journal
Nanocellulose-Block Copolymer Films for the Removal of Emerging Organic Contaminants from Aqueous Solutions
Previous Article in Special Issue
Wear Performance Analysis of Ni–Al2O3 Nanocomposite Coatings under Nonconventional Lubrication
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Asymptotic Homogenization Applied to Flexoelectric Rods

by
David Guinovart-Sanjuán
1,
Jose Merodio
2,*,
Juan Carlos López-Realpozo
3,
Kuppalapalle Vajravelu
1,
Reinaldo Rodríguez-Ramos
3,
Raúl Guinovart-Díaz
3,
Julián Bravo-Castillero
4 and
Federico J. Sabina
4
1
Department of Mathematics, University of Central Florida, 4393 Andromeda Loop N, Orlando, FL 32816, USA
2
Departamento de Mecánica de los Medios Continuos y T. Estructuras, E.T.S. de Caminos, Canales y Puertos, Universidad Politécnica de Madrid, C.P. 28040 Madrid, Spain
3
Departamento de Matemática, Universidad de La Habana, San Lazaro y L, La Habana 10400, Cuba
4
Instituto de Investigaciones en Matemáticas Aplicadas y en Sistemas, Universidad Nacional Autónoma de México, Alcaldía Álvaro Obregón, Apartado Postal 20-126, CDMX Mexico 01000, Mexico
*
Author to whom correspondence should be addressed.
Materials 2019, 12(2), 232; https://doi.org/10.3390/ma12020232
Submission received: 6 December 2018 / Revised: 26 December 2018 / Accepted: 7 January 2019 / Published: 11 January 2019

Abstract

:
In this manuscript, the equilibrium problem for a flexoelectric one-dimensional composite material is studied. The two-scales asymptotic homogenization method is used to derive the homogenized formulation of this problem. The manuscript offers a step-by-step methodology to derive effective coefficients and to solve local problems. As an illustrative example, results reported in the literature for piezoelectric composites are obtained as a particular case of the formulation derived here. Finally, three flexoelectric/piezoelectric composites are studied to illustrate the influence of the flexoelectric property on the effective coefficients and the global behavior of the structure.

1. Introduction

In the last 35 years, the study of piezoelectric composites has increased because of their use in several areas of engineering [1], contributing to the development of new mechanical structures that are used in high-technology devices [2,3]. The wide range of applications of piezoelectric composites has advanced the development of mathematical, experimental, and computational models related to the study of the properties of these materials, see for instance References [4,5,6]. Piezoelectric composites have been extensively studied, but only in very recent years, flexoelectric materials have taken an important place in the new scientific works, see Reference [7].
Piezoelectricity is usually expressed as an interaction between mechanical strain and one of the electrical variables: the electric field, the electric displacement, or the electric polarization. In Reference [8], the author examines the consequences of considering an additional, linear, electromechanical effect: an interaction between the strain and the polarization gradient. Experiments have shown that when large amplitude mechanical disturbances propagate through a dielectric medium a voltage is developed across the ends of the sample, see Reference [9]. Some of these experiments with nonpiezoelectric elements have been shown to produce a polarization charge upon shock loading. Strain-gradient-induced polarization is known as the flexoelectric effect. According to Reference [10], flexoelectricity is an electromechanical effect of dielectric materials whereby they exhibit a spontaneous electrical polarization induced by a strain gradient (inhomogeneous deformation). The flexoelectric effect can be considered as a high-order electromechanical phenomenon with respect to the piezoelectric effect [7].
Flexoelectricity was first theoretically predicted in the 1950’s and described from a phenomenological standpoint in the 1960’s [11]. Recently, flexoelectric materials have gained prominence due to the development of new methods for manufactured materials with coupled mechanical and electrical behavior [11] and their applications in membrane structures [12,13,14].
In References [15,16], important contributions to the study of flexoelectric materials have been presented highlighting the differences between piezoelectric and flexoelectric structures. The authors emphasize the significance of the flexoelectric effect discussing several applications which include materials with large flexoelectric coefficients as well as cases in which electromechanical coupling of piezoelectricity is not present and in applications related to soft materials. Moreover, some important applications of the flexoelectric effect are focused on the study of biological membranes and the development of piezoelectric structures without using piezoelectric materials, sensing, actuating, or energy harvesting.
Many authors have studied the flexoelectric mechanical equilibrium equations (for instance References [11,17]), bringing together different approaches to solve the problem at hand. One of the most common methods in solid mechanics to approximate the solution of the equilibrium problem of a solid (elastic, piezoelectric, flexoelectric, etc.) is the two-scales asymptotic homogenization method (AHM) [18,19]. The multiscale asymptotic homogenization method has been widely used to derive the effective properties of composite materials with different mechanical properties. In References [18,20,21,22], the AHM is used to derive the effective properties of elastic composites (laminates, fibrous, wavy, etc.). These ideas are extended to the study of piezoelectric structures [23]; viscoelastic composites [24]; and thermo-piezoelectric materials [25]. To the best of the authors knowledge, a methodology to find the effective properties of flexoelectric materials has not been presented before. For that reason, this work acts to extend those methodologies (that find effective properties) to include the case of flexoelectric materials using the two-scales asymptotic homogenization method, which provides the first step to subsequently study three-dimensional curvilinear flexoelectric structures.
In general, the two-scale asymptotic homogenization method is used to find the effective properties of periodic composites, i.e., for any mechanical property p ( x ) , x R n , there is a T R n s.t. p ( x + T ) = p ( x ) for all x R n . For non-periodic media, statistical and numerical methods have been developed to find representative volume elements and the associated effective properties [26,27,28]. In Reference [29], the self-consistent method in explicit form (effective field) and implicit form (effective medium) as well as the AHM are used to derive the effective piezoelectric properties of fibrous composites with randomly positioned fibers and periodically distributed fibers, showing good agreement among the different methods.
The work is organized as follows. In Section 2, a theoretical framework for the equilibrium of a one-dimensional periodic flexoelectric composite material is presented. The asymptotic expansions of the displacement and the electric potential are introduced in Section 3, using the two-scales asymptotic homogenization method, and a methodology to derive the expression of the effective coefficients depending on the so-called local function expressions is illustrated. The process to obtain the analytic solutions of the local problems is shown. A variational formulation of the local problem is proposed and the system is solved using the classical finite element method (FEM). As an important benchmark, the model reproduces the solution of the local problems reported in Reference [30] for a one-dimensional piezoelectric structure. Finally, in Section 4, the theoretical results derived in the previous sections are used to study some numerical cases. The effective coefficients reported in References [19,30] for a two-material piezoelectric structure made of lead zirconate titanate (PZT) C91 and P-82 are derived using the model developed here, which makes use of the solution of the local problem variational formulation (FEM). In order to extend the study for the case of randomly distribute one-dimensional flexoelectric composites, a two-constituents structure Ω N is considered for which the materials properties are randomly distributed along the rod with a binomial probability. Using the concept of correlation length for the particular case of one-dimensional binomial distributed structure, as reported in Reference [31], the methodology is then extended to the case of non-periodic structures. As it is shown in Reference [29], the effective properties of random composites converge to the effective properties of periodic structures when the length of Ω N increases. In addition, the combination of flexoelectric/piezoelectric/non-piezoelectric composites are studied in order to illustrate the influence of the flexoelectricity on the effective properties. One of the composites is considered as the combination of two flexoelectric materials, barium titanate (BaTiO 3 ) and strontium titanate (SrTiO 3 ) [32,33,34]. Other composites are combinations of these two flexoelectric materials with a exclusively piezoelectric material, polyvinylidene fluoride (PVDF) [35,36], a non-piezoelectric polymer described in Reference [37] and PZT-7A reported in Reference [38]. From these material combinations important conclusions are derived with respect to the influence of the flexoelectric property in the global behavior of the composites. Finally, a comparison between the solutions of a heterogeneous and homogeneous problem for a one-dimensional flexoelectric structure is presented. The finite element method is used to solve the problems of two different examples considering the same prescribed boundary conditions but different external forces.

2. Flexoelectric One-Dimensional Problem

A one-dimensional flexoelectric composite rod Ω is considered. The constitutive relations between stress σ and electrical displacement D with strain ϵ and electric field E are written in the following form (see Reference [17]),
σ ( x ) = C ( x ) ϵ ( x ) e ( x ) E ( x ) + μ ( x ) d E ( x ) d x ,
D ( x ) = e ( x ) ϵ ( x ) + μ ( x ) d ϵ ( x ) d x + κ ( x ) E ( x ) ,
where C ( x ) , e ( x ) , κ ( x ) , and μ ( x ) denote the stiffness, piezoelectric, permittivity, and flexoelectric tensors, respectively. Linear deformations of the composite, the strain, and the electric field are studied in terms of the displacement u and the electric potential ϕ as
ϵ ( x ) = d u ( x ) d x ,
E ( x ) = d ϕ ( x ) d x .
Assume the material parameters C, e, κ , μ are rapidly oscillating periodic functions along the rod with respect to the variable y, i.e.,
C ( y ) C x ε , e ( y ) e x ε ,
κ ( y ) κ x ε , μ ( y ) μ x ε ,
where y = x / ε Y is the local variable, since 0 < ε < < 1 is a very small parameter that characterizes the periodicity of the structure, and Y is the periodic cell. The composite has the structure presented in Figure 1.
To simplify the equations, material parameters and functions are considered to depend on the position variable x or local variable y, i.e., u u ( x ) , ϕ ϕ ( x ) , C C ( y ) , e e ( y ) , κ κ ( y ) , μ μ ( y ) . The equilibrium equations presented in Reference [17] for one-dimensional solids take the following expression:
d σ d x + f = d d x C d u d x + e d ϕ d x μ d 2 ϕ d x 2 + f = 0 in Ω ,
d D d x = d d x e d u d x + μ d 2 u d x 2 κ d ϕ d x = ρ in Ω ,
with boundary conditions
u ( 0 ) = u 0 , d u d x ( 0 ) = w 0 , ϕ ( 0 ) = φ 0 ,
σ ( 1 ) = S 1 , μ ( 1 ) d E d x ( 1 ) = r 1 , D ( 1 ) = τ 1 ,
where the functions f and ρ are the external forces and body charge density, respectively. On the other hand, the values u 0 , S 1 , w 0 , r 1 , φ 0 , and τ 1 are the prescribed displacement, traction, displacement derivative, high-order traction, electric potential and charge, respectively. Perfect contact conditions at the interface are considered, i.e., function displacement u is differentiable and classical stress C ϵ + e E , higher order stress μ d E d x (see Reference [17]), electric potential ϕ and electric displacement D are considered continuous functions at the interfaces of the materials. The boundary value problem (5)–(8) is cast as a linear system of third order ordinary differential equations for which, under the perfect contact conditions as mentioned above, existence and uniqueness of the solution on the interval [ 0 , 1 ] is guaranteed [39,40].

3. Asymptotic Homogenization Method

The boundary value problem (5)–(8) has rapidly oscillating coefficients. In order to approach the solution of the problem, an equivalent homogenized problem must be obtained. To derive the expression of the homogenized system, the two-scales asymptotic homogenization method is used. In Reference [19], the asymptotic expansion to the solution of the piezoelectric problem is reported as
u ( x , y ) = v 0 + k = 1 ε k N k ( y ) d k v 0 d x k + Π k ( y ) d k ϕ 0 d x k ,
ϕ ( x , y ) = ϕ 0 + k = 1 ε k Ψ k ( y ) d k v 0 d x k + Θ k ( y ) d k ϕ 0 d x k ,
where v 0 ( x ) and ϕ 0 ( x ) only depend on the global variable x and the functions N k , Π k , Ψ k , and Θ k are ε -periodic continuous functions that only depend on the local variable y, i.e., N k ( y + ε ) = N k ( y ) , Π k ( y + ε ) = Π k ( y ) , Ψ k ( y + ε ) = Ψ k ( y ) , and Θ k ( y + ε ) = Θ k ( y ) .
Now, the expansions (9)–(10) are introduced into the equilibrium Equations (5)–(6), for which one needs to consider the derivative operator
d u d x = u x + 1 ε u y , d ϕ d x = ϕ x + 1 ε ϕ y ,
where ( · ) x and ( · ) y denote the partial derivative with respect to x and y, respectively. After some simple manipulations, (5) and (6) take the following expressions:
k = 1 m = 1 ε k L k m d m v 0 d x m + P k m d m ϕ 0 d x m + f = 0 ,
k = 1 m = 1 ε k Q k m d m v 0 d x m + R k m d m ϕ 0 d x m = ρ ,
respectively, where L k m , P k m , Q k m , and R k m denote the coefficient of the m derivative of v 0 and ϕ 0 that multiplies ε k . For instance, L 11 is the coefficient of the first derivative of v 0 that multiplies ε 1 . As a result of the continuity requirements, when ε 0 , the coefficients of k = 1 are equated to zero. The non-identical zero coefficients are the so-called local problems. In what follows, we characterize these problems.

3.1. Local Problems

The following local problems have to be considered. The system of local problems (LQ):
L 11 d d y C + C d N 1 d y + e d Ψ 1 d y = 0 ,
L 12 d d y μ d 2 Ψ 2 d y 2 + 2 μ d Ψ 1 d y μ d 2 Ψ 1 d y 2 = 0 ,
Q 11 d d y e + e d N 1 d y κ d Ψ 1 d y = 0 ,
Q 12 d d y μ + μ d 2 N 2 d y 2 + 2 μ d N 1 d y + μ d 2 N 1 d y 2 = 0 ,
with interface conditions
C + C d N 1 d y + e d Ψ 1 d y = μ d 2 Ψ 2 d y 2 + 2 μ d Ψ 1 d y = e + e d N 1 d y κ d Ψ 1 d y = μ + μ d 2 N 2 d y 2 + 2 μ d N 1 d y = 0 ,
where · = ( · ) ( + ) ( · ) ( ) denotes the jump at the interface.
Furthermore, the system of local problems (PR):
P 11 d d y e + C d Π 1 d y + e d Θ 1 d y = 0 ,
P 12 d d y μ + μ d 2 Θ 2 d y 2 + 2 μ d Θ 1 d y μ d 2 Θ 1 d y 2 = 0 ,
R 11 d d y κ + κ d Θ 1 d y e d Π 1 d y = 0 ,
R 12 d d y μ d 2 Π 2 d y 2 + 2 μ d Π 1 d y + μ d 2 Π 1 d y 2 = 0 ,
with interface conditions
e + C d Π 1 d y + e d Θ 1 d y = μ + μ d 2 Θ 2 d y 2 + 2 μ d Θ 1 d y = κ + κ d Θ 1 d y e d Π 1 d y = μ d 2 Π 2 d y 2 + 2 μ d Π 1 d y = 0 .
The LQ system is derived from the coefficients multiplying the derivatives of v 0 in each equation and it relates the local functions N 1 , Ψ 1 , N 2 , and Ψ 2 . Similarly, the PR system is obtained from the coefficients of the derivatives of ϕ 0 and it relates the local functions Θ 1 , Π 1 , Θ 2 , and Π 2 .

3.2. Methodology to Solve Local Problems Using a System of Linear Equations

To obey conditions (13)–(20), the LQ and PR systems need to be solved. More specifically, the functions d N 1 d y , d Ψ 1 d y , d 2 N 2 d y 2 , d 2 Ψ 2 d y 2 , d Π 1 d y , d Θ 1 d y , d 2 Π 2 d y 2 , and d 2 Θ 2 d y 2 need to be obtained. Later, we show that the local functions N 1 , N 2 , Ψ 1 , Ψ 2 , Π 1 , Π 2 , Θ 1 , and Θ 1 are not necessary to find the homogenized problem.
The system of local problems is solved using ideas described in Reference [20]. Here, we summarize the methodology. Let us focus on the LQ system. Integrating Equations (13) and (15) the following system is derived:
C e e κ d N 1 d y d Ψ 1 d y = λ 1 C λ 3 e ,
where λ α λ α ( x ) , α = 1 , 3 . Thus, the solutions of the system are
d N 1 d y = κ e 2 + C κ λ 1 + e e 2 + C κ λ 3 1 ,
d Ψ 1 d y = e e 2 + C κ λ 1 C e 2 + C κ λ 3 .
To find expressions for λ α in terms of the materials parameters C, e, κ and μ , one needs to use the average operator with respect to the local variable. This average operator is given by
· = 1 V Y Y ( · ) d y ,
where V Y = | Y | . From the periodicity of N 1 and Ψ 1 and taking the average operator in both sides of the Equations (21) and (22), the following system for λ α is obtained:
κ e 2 + C κ λ 1 + e e 2 + C κ λ 3 = 1 ,
e e 2 + C κ λ 1 C e 2 + C κ λ 3 = 0 .
It follows that λ 1 and λ 3 are solutions of the system (24) and (25). Considering C, e, κ , and μ to be constant functions along each material, one gets that d 2 N 1 d y 2 = d 2 Ψ 1 d y 2 = 0 . Then, Equations (14) and (16) are transformed into
d d y μ d 2 Ψ 2 d y 2 + 2 μ d Ψ 1 d y = 0 ,
d d y μ + μ d 2 N 2 d y 2 + 2 μ d N 1 d y = 0 ,
respectively. Finally, solutions of the LQ system are given by (21), (22) and
d 2 N 2 d y 2 = 2 k Δ λ 1 2 e Δ λ 3 + 1 μ λ 4 + 1 ,
d 2 Ψ 2 d y 2 = 2 e Δ λ 1 + 2 C Δ λ 3 + 1 μ λ 2 ,
where Δ = e 2 + C κ and λ i , i = 1 , 2 , 3 , 4 , satisfy the system of equations
κ / Δ 0 e / Δ 0 2 κ / Δ 0 2 e / Δ 1 / μ e / Δ 0 C / Δ 0 2 e / Δ 1 / μ 2 C / Δ 0 λ 1 λ 2 λ 3 λ 4 = 1 1 0 0 .
Similarly, solutions of the PR system (17)–(20) are obtained. Expressions for these solutions are given in Appendix A.

3.3. Variational Formulation of the LQ System. Finite Element Method (FEM)

The finite element method is a strong and convenient numerical tool to solve systems of ordinary or partial differential equations. It is an alternative method to approximate the solution for the LQ systems (13)–(16). This is the procedure. Let us take Y = [ 0 , γ ) [ γ , 1 ] , where γ ( 0 , 1 ) . Consider the test functions v i , i = 1, 2, 3, 4, multiplying the LQ problem. Integrating along Y, using integration by parts and taking into account the continuity of the function at the interface, the following variational formulation for the LQ system is derived:
Y C + C d N 1 d y + e d Ψ 1 d y v 1 y d y = 0 ,
Y μ d 2 Ψ 2 d y 2 + 2 μ d Ψ 1 d y v 2 y d y = 0 ,
Y e + e d N 1 d y κ d Ψ 1 d y v 3 y d y = 0 ,
Y μ + μ d 2 N 2 d y 2 + 2 μ d N 1 d y v 4 y d y = 0 .
Periodic conditions for the local functions and perfect contact at the interface are considered. Now, using Equations (31)–(34) and the standard finite element method [41] on the domain Y, the local functions are approximated.

3.4. Effective Coefficients and Homogenized Problem

To find the effective coefficients, the average operator is used on the coefficients of ε 0 in Equations (11) and (12), i.e., one deals with L 0 m , P 0 m , Q 0 m , and R 0 m . Taking into account the periodicity of the local functions as well as C, μ , e, and κ , the non-identical zeros are
L 02 = C + C d N 1 d y + e d Ψ 1 d y ,
L 03 = 2 μ d Ψ 1 d y + μ d 2 Ψ 2 d y 2 ,
P 02 = e + C d Π 1 d y + e d Θ 1 d y ,
P 03 = μ + 2 μ d Θ 1 d y + μ d 2 Θ 2 d y 2 .
Q 02 = e + e d N 1 d y κ d Ψ 1 d y ,
Q 03 = μ + 2 μ d N 1 d y + μ d 2 N 2 d y 2 ,
R 02 = κ + κ d Θ 1 d y e d Π 1 d y ,
R 03 = 2 μ d Π 1 d y + μ d 2 Π 2 d y 2 .
Substituting the expressions of the solutions of the LQ and PR systems, which are given in (21), (22), (28), (29), and the Appendix Equations (A1)–(A4), into the Equations (35)–(42), the following relations are obtained
L 03 R 03 0 ; P 03 Q 03 ; P 02 Q 02 .
Therefore, the effective coefficients of Equations (11) and (12) are
C ^ = C + C d N 1 d y + e d Ψ 1 d y ,
e ^ = e + C d Π 1 d y + e d Θ 1 d y e + e d N 1 d y κ d Ψ 1 d y ,
κ ^ = κ + κ d Θ 1 d y e d Π 1 d y ,
μ ^ = μ + 2 μ d Θ 1 d y + μ d 2 Θ 2 d y 2 μ + 2 μ d N 1 d y + μ d 2 N 2 d y 2 .
The original problem (5) and (6), by means of the homogenization technique, is transformed into the following homogenized problem (in terms of the effective coefficients):
d d x C ^ d v 0 d x + e ^ d ϕ 0 d x μ ^ d 2 ϕ 0 d x 2 + f = 0 in Ω ,
d d x e ^ d v 0 d x + μ ^ d 2 v 0 d x 2 κ ^ d ϕ 0 d x = ρ in Ω ,
v 0 ( 0 ) = u 0 , d v 0 d x ( 0 ) = w 0 , ϕ 0 ( 0 ) = φ 0 ,
σ 0 ( 1 ) = S 1 , μ ^ ( 1 ) E 0 d x ( 1 ) = r 1 , D 0 ( 1 ) = τ 1 ,
where σ 0 , E 0 , and D 0 are the effective stress, effective electric field, and effective electrical displacement, respectively.

4. Analysis of Numerical Results

4.1. Model Validation with the Particular Case of Piezoelectric Materials

To validate the model, our goal is to reproduce the results of the one-dimensional piezoelectric composite analyzed in Reference [30]. The constituents of this one-dimensional piezoelectric structure are lead zirconate titanate (PZT) C-91 and P-82. The materials properties are given in Table 1 of Reference [30]. For this particular case, the flexoelectric coefficient μ associated with the components of the structure is considered zero. Therefore, the equilibrium problems (5)–(8) take the expression (1) and (2) of Reference [30]. Using the asymptotic expansion (9) and (10) and following the methodology described in the previous section, the non-vanishing local problems are (13), (15), (18), and (20), which corresponds with Equations (23) and (24) in Reference [30]. Finally, the expressions of the effective coefficients for this piezoelectric composite are given by Equations (43), (44), and (45).
In Table 1, the values reported for the effective coefficients C ^ , e ^ , κ ^ of this piezoelectric composite considered in References [19,30] are compared with the results obtained using the present model (Section 3) and the results obtained using FEM (solving the variational formulation of the LQ and PR systems). The volume fraction of C-91, for the numerical calculation, is considered to be V 1 (0, 0.4 , 0.8 , 1). To avoid singularities of solutions of LQ and PR systems, the value of flexoelectric property for the piezoelectric composites is taken as μ = 10 20 . Suffice to say that the results reported in Reference [19,30] are reproduced using both the AHM (present model) and FEM (variational formulation) methods. Therefore, the piezoelectric composite results are obtained as a particular case of the flexoelectric composite formulation developed here.

4.2. Flexoelectric Composite Rod without Ideal Periodicity

It is interesting to study the the above-described methodology for the case of random composites, i.e., composites without periodicity. A two-elements flexoelectric rod where the constituents are L 1 (BaTiO 3 ) and L 2 (SrTiO 3 ) with a fixed length of 1 unit is considered. Assume Ω N = [ 0 , N ] , for every 0 n < N , n N , the region [ n , n + 1 ] is occupied by the layer L 1 with probability p = 0.7146 . The numbers N 1 and N 2 represent the amount of components L 1 and L 2 in Ω N , respectively, and N = N 1 + N 2 . The average operator is considered
· N = 1 N Ω N ( · ) d x = N 1 N ( · ) ( 1 ) + N 2 N ( · ) ( 2 ) .
This average operator is a particular case of the Equation (2.2.1) of Reference [31] for one-dimensional structures. From the results described in Section 3 and considering the average operator (51), the effective properties (43)–(46) are easily computed for a structure without periodicity. To compare the behavior of the effective properties for random composites with the ones for periodic structures, a periodic composite Ω is considered. The constituents of Ω are BaTiO 3 with volume fraction p = 0.7146 and SrTiO 3 with volume fraction 1 p . In Table 2, the material properties of BaTiO 3 and SrTiO 3 are given. Figure 2 illustrates the behavior of three experiments with random composites considering a probability p = 0.7146 for L 1 and a composite with periodicity and a volume fraction of p = 0.7146 for BaTiO 3 for N [ 1 , 2 25 ] . The three random composites were randomly generated with MATLAB using the so-called rand function.
Figure 2 shows that the values of effective properties for non-periodic composites coincide with the values of the effective coefficients for a periodic structure when the number of layers increase. It follows that the methodology is valid for structures with no periodicity, making it possible to approximate the effective properties for random composites considering the periodic structures. Similar results are obtained in Reference [42], for the particular case of an elastic Fibonacci laminate composite.

4.3. Influence of the Flexoelectric Property on the Effective Coefficients

In this section, the influence of the flexoelectric property on the effective coefficients is studied. A bi-material one-dimensional composite is considered. In Table 2, the properties of the flexoelectric materials barium titanate (BaTiO 3 ) and strontium titanate (SrTiO 3 ), the active piezoelectric materials polyvinylidene fluoride (PVDF) and PZT-7A, and soft non-piezoelectric polymer are shown. To illustrate the effective coefficients of flexoelectric/piezoelectric structures, several combinations of flexoelectric composite with other well-known piezoelectric materials are considered.
In Figure 3, the influence of the flexoelectricity in combination with an active polymer PVDF which has piezoelectric properties on the effective characteristics is studied using three composites. The first composite is made of the flexoelectric material BaTiO 3 and the piezoelectric active polymer PVDF; the second one is the composition of the flexoelectric composite SrTiO 3 and PVDF; the materials of the third composite are the piezoelectric constituents PZT-7A and PVDF (see Reference [38]). With the values given in Table 2 and the methodology described in Section 3, the effective coefficients for each composite are obtained. Figure 3 illustrates the behavior of the effective coefficients for different volume fractions of the composites. The effective coefficient C ^ continues increasing for the three combinations of materials. In addition, as is shown in the figure, the influence of the flexoelectric materials on the piezoelectric and the dielectric effective coefficients is remarkable. For both BaTiO 3 -PVDF and SrTiO 3 -PVDF, the piezoelectric effective coefficient e ^ increases until V 1 = 0.9 and V 1 = 0.8 , respectively, from which the property decays. On the other hand, for PZT-PVDF, the piezoelectric effective tensor decreases for all values of V 1 . Similarly to C ^ , the dielectric effective coefficient κ ^ increases with the volume fraction for the three combinations of materials, although for the composites formed by a flexoelectric material this property is greater than that of the PZT-PVDF composite. The effective flexoelectric coefficient μ ^ varies only for the composites with a flexoelectric component (either BaTiO 3 or SrTiO 3 ).
To study the influence of non-flexoelectric and non-piezoelectric constituents in the effective properties for composites with a flexoelectric component, three structures are studied. One of the structures is the combination of the flexoelectric materials BaTiO 3 and SrTiO 3 ; the constituents for the second structure are SrTiO 3 and a non-piezoelectric polymer; the last structure is made of strontium titanate (SrTiO 3 ) and a piezoelectric material PZT-7A. In Figure 4, the effective properties of these three composites for different volume fractions are shown. An important detail to highlight is that for all the effective properties in Figure 4, the values of the effective coefficients of the SrTiO 3 -polymer are lower than the values of the effective coefficients of the other two combinations. The non-piezoelectric polymer properties diminish the values of the composite effective coefficients, compared with the other composites. On the other hand, the combination of the two flexoelectric materials reinforce the values of C ^ , e ^ , κ ^ , and μ ^ , compared with the composite made of the non-piezoelectric polymer. It is shown in Figure 4 that the effective piezoelectric property e ^ is almost zero for the combination of SrTiO 3 with a non-piezoelectric polymer. A similar situation is shown for the flexoelectric parameter, i.e., e ^ remains close to zero for the composites that involve a non-flexoelectric component (either SrTiO 3 -polymer or SrTiO 3 -PZT). Therefore, it can be concluded that for composites made of flexoelectric material with non-flexoelectric material, the global behavior of the effective coefficients of the structure is similar to a non-flexoelectric material. A parallel situation occurs when considering the composition of a piezoelectric with a non-piezoelectric structure as it is illustrated in Figure 4. Hence, for composites made of a piezoelectric material with a non-piezoelectric material the global behavior of the effective coefficients of the structure is similar to a non-piezoelectric material.
The dependence of the effective coefficients C ^ , e ^ , and κ ^ on the local functions d N 1 d y , d Ψ 1 d y , d Θ 1 d y , d Π 1 d y is shown in Equations (43)–(45). From Equations (21), (22), the Appendix Equations (A1), (A3), and the solutions of the systems (30), (A5), the non-influence of the μ parameter in the local functions is derived. Therefore, the effective coefficients C ^ , e ^ , and κ ^ are not affected by the flexoelectric component in the case of non-flexoelectric materials. From the homogenized problem (47)–(50), it can be followed that the influence of the effective flexoelectric component in the effective stress and electric displacement cannot be underestimated.
The investigation of the flexoelectric effect goes beyond the theoretical analysis of the materials properties to include real measurements [43]. Recent papers have established detailed descriptions of their experiments to determine the flexoelectric properties of PVDF. More precisely, References [44] and [45] describe experiments to find the effective flexoelectric components μ 1123 and μ 2312 of the fourth order tensor μ i j k l for polyvinylidene fluoride. On account of the great difficulty in obtaining some experimental measurements, the relationship between the strain gradient and torque is deduced theoretically and further verified with finite element analysis. The approach is applied to test responses in bars machined from bulk polyvinylidene fluoride [45]. On the other hand, in Reference [46], the flexoelectricity of prototypical semicrystalline polymer, α -phase PVDF, films are investigated. The paper presents a step-by-step description of the experiment highlighting the direct flexoelectric effect in the α -phase PVDF films.

4.4. Solution of the Heterogeneous and Homogenized Problems

In this section, the homogenized problem (47)–(50) is solved for a flexoelectric one-dimensional structure. In order to compare the approximation between the solution of the heterogeneous problem and the homogenized problem, the finite element method is used. Consider Ω a two-element structure made of barium titanate (BaTiO 3 ) and strontium titanate (SrTiO 3 ) periodically distributed along Ω . In order to illustrate the example, the prescribed conditions for the heterogeneous (5)–(8) and the homogeneous (47)–(50) problems are given as follows: u 0 = 1 , w 0 = 1 , φ 0 = 0 , S 1 = 0 , r 1 = 0 , and τ 1 = 0 . Two different volume fractions for BaTiO 3 are considered V 1 = 0.5 , case 1, and V 1 = 0.8 , case 2. The external forces are taken as f ( x ) = e x for case 1 and f ( x ) = sin ( x ) for case 2. The body charge density ρ is considered identically zero along Ω . The values of the effective coefficients of the homogenized problem are computed using the methodology described in Section 3 and are shown in Table 3.
In Figure 5, a comparison of the stress obtained from the solution of the heterogeneous problem (5)–(8) and the homogeneous equivalent problem (47)–(50) for ε = 1 / 5 is shown. As a result of the continuity of the coefficients of the homogeneous problem, the stress function has a smooth behavior along Ω . On the other hand, the discontinuity of the coefficients in the heterogeneous problem leads to the obstacles in the solution, see Reference [47].

5. Conclusions

In this manuscript, the equilibrium equation for a flexoelectric composite has been studied. The two-scales asymptotic homogenization method is used to find the homogenized problem formulation. The work presents the expressions of the effective coefficients for the case of a flexoelectric structure. The local problems are derived along with their corresponding variational formulations. The manuscript offers a detailed description to solve the local problems for one-dimensional structures. The effective coefficients reported in References [19,30] for piezoelectric composite are computed here as particular cases of both the present model (Section 3) and the finite element method. The methodology is used to find the effective properties of random structures. The results illustrate the convergence of the effective coefficients for random structures to the case of a composite with ideal periodicity. In addition, the influence of the flexoelectric parameter on the global behavior of composites has been studied. The effective properties of three bi-materials (composites) made of combinations of BaTiO 3 , SrTiO 3 , PVDF, a soft non-piezoelectric polymer, and PZT-7A were obtained. The influence of the flexoelectric parameter ( μ ) on the effective stress and electric displacement was shown. As an example, a comparison between the solutions of the heterogeneous and homogeneous problem for a flexoelectric structure is presented. A good approximation between the solutions can be appreciated when ε is zero. The results presented in this manuscript can be extended to more complex structures with wavy behavior and three-dimensional composites.

Author Contributions

The individual contributions of each author are the following: conceptualization, D.G.-S., J.C.L.-R. and R.R.-R.; methodology, D.G.-S., K.V. and J.M.; software, R.G.-D. and F.J.S.; validation, D.G.-S., J.C.L.-R. and R.G.-D.; formal analysis, D.G.-S. and J.B.-C.; writing—original draft preparation, D.G.-S., R.R.-R. and K.V.; writing—review and editing, J.M., J.B.-C. and F.J.S.

Funding

This research received no external funding.

Acknowledgments

The funding of Proyecto Nacional de Ciencias Básicas 2013–2015 (Project No. 7515) are gratefully acknowledged. Thanks to Departamento de Matemáticas y Mecánica IIMAS-UNAM and FENOMEC for their support and Ramiro Chávez Tovar and Ana Pérez Arteaga for computational assistance. David Guinovart-Sanjuan thanks Ariel Ramírez-Torres and Serge Dumont for computational assistance. The authors want to thank the reviewers for the time spent in reading the manuscript and for the helpful comments

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
BaTiO 3 Barium titane
SrTiO 3 Strontium titanate
PVDFPolyvinylidene fluoride
PZTLead zirconate titanate
AHMtwo-scales Asymptotic Homogenization method
FEMFinite Element method

Appendix A. Solutions of the (PR) System

The solutions of the (PR) system (17)–(20) are
d Π 1 d y = k Δ λ 5 e Δ λ 7 ,
d 2 Π 2 d y 2 = 2 k Δ λ 5 + 2 e Δ λ 7 + 1 μ λ 8 ,
d Θ 1 d y = e Δ λ 5 + C Δ λ 7 1 ,
d 2 Θ 1 d y 2 = 2 e Δ λ 5 2 C Δ λ 7 + 1 μ λ 6 + 1 ,
where Δ = e 2 + C κ and λ i for i = 5 , 6 , 7 , 8 , satisfy the following system of equations
κ / Δ 0 e / Δ 0 2 κ / Δ 0 2 e / Δ 1 / μ e / Δ 0 C / Δ 0 2 e / Δ 1 / μ 2 C / Δ 0 λ 5 λ 6 λ 7 λ 8 = 0 0 1 1 .

References

  1. Ting, R.Y. Evaluation of new piezoelectric composite materials for hydrophone applications. Ferroelectrics 1986, 67, 143–157. [Google Scholar] [CrossRef]
  2. Gururaja, T.R.; Schulze, W.A.; Cross, L.E.; Newnham, R.E.; Auld, B.A.; Wang, Y.J. Piezoelectric Composite Materials for Ultrasonic Transducer Applications. Part I: Resonant Modes of Vibration of PZT Rod-Polymer Composites. IEEE Trans. Sonics Ultrason. 1985, 32, 481–498. [Google Scholar] [CrossRef] [Green Version]
  3. Akdogan, E.K.; Allahverdi, M.; Safari, A. Piezoelectric composites for sensor and actuator applications. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 2005, 52, 746–775. [Google Scholar] [CrossRef] [PubMed]
  4. Lin, X.; Zhou, K.; Zhang, X.; Zhang, D. Development, modeling and application of piezoelectric fiber composites. Trans. Nonferrous Met. Soc. China 2013, 23, 98–107. [Google Scholar] [CrossRef]
  5. Zhu, J.; Wang, Z.; Zhu, X.; Yang, B.; Fu, C. Theoretical and experimental study on the effective piezoelectric properties of 1-3 type cement-based piezoelectric composites. Materials 2018, 11, 1698. [Google Scholar] [CrossRef] [PubMed]
  6. Fang, X.Q.; Huang, M.J.; Liu, J.X.; Feng, W.J. Dynamic effective property of piezoelectric composites with coated piezoelectric nano-fibers. Compos. Sci. Technol. 2014, 98, 79–85. [Google Scholar] [CrossRef]
  7. Yudin, P.V.; Tagantsev, A.K. Fundamentals of flexoelectricity in solids. Nanotechnology 2013, 24, 432001. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Mindlin, R. Polarization gradient in elastic dielectrics. Int. J. Solids Struct. 1968, 4, 637–642. [Google Scholar] [CrossRef]
  9. Harris, P. Mechanism for the Shock Polarization of Dielectrics. J. Appl. Phys. 1965, 36, 739–741. [Google Scholar] [CrossRef]
  10. Zubko, P.; Catalan, G.; Tagantsev, A.K. Flexoelectric Effect in Solids. Annu. Rev. Mater. Res. 2013, 43, 387–421. [Google Scholar] [CrossRef] [Green Version]
  11. Enakoutsa, K.; Corte, A.D.; Giorgio, I. A model for elastic flexoelectric materials including strain gradient effects. Math. Mech. Solids 2016, 21, 242–254. [Google Scholar] [CrossRef]
  12. Mohammadi, P.; Liu, L.; Sharma, P. A theory of flexoelectric membranes and effective properties of heterogeneous membranes. J. Appl. Mech. 2013, 81, 011007. [Google Scholar] [CrossRef]
  13. Deng, Q.; Liu, L.; Sharma, P. Flexoelectricity in soft materials and biological membranes. J. Mech. Phys. Solids 2014, 62, 209–227. [Google Scholar] [CrossRef]
  14. Narita, F.; Shindo, Y.; Mikami, M. Analytical and experimental study of nonlinear bending response and domain wall motion in piezoelectric laminated actuators under ac electric fields. Acta Mater. 2005, 53, 4523–4529. [Google Scholar] [CrossRef]
  15. Krichen, S.; Sharma, P. Flexoelectricity: A Perspective on an Unusual Electromechanical Coupling. J. Appl. Mech. 2016. [Google Scholar] [CrossRef]
  16. Majdoub, M.S.; Sharma, P.; Cagin, T. Enhanced size-dependent piezoelectricity and elasticity in nanostructures due to the flexoelectric effect. Phys. Rev. B 2008, 77, 125424. [Google Scholar] [CrossRef]
  17. Nguyen, B.; Zhuang, X.; Rabczuk, T. Numerical model for the characterization of Maxwell-Wagner relaxation in piezoelectric and flexoelectric composite material. Comput. Struct. 2018, 208, 75–91. [Google Scholar] [CrossRef]
  18. Meguid, S.; Kalamkarov, A. Asymptotic homogenization of elastic composite materials with a regular structure. Int. J. Solids Struct. 1994, 31, 303–316. [Google Scholar] [CrossRef]
  19. Otero, J.; Castillero, J.; Ramos, R. Homogenization of heterogeneous piezoelectric medium. Mech. Res. Commun. 1997, 24, 75–84. [Google Scholar] [CrossRef]
  20. Tsalis, D.; Chatzigeorgiou, G.; Charalambakis, N. Homogenization of structures with generalized periodicity. Compos. Part B Eng. 2012, 43, 2495–2512. [Google Scholar] [CrossRef]
  21. Chen, C.M.; Kikuchi, N.; Rostam-Abadi, F. An enhanced asymptotic homogenization method of the static and dynamics of elastic composite laminates. Comput. Struct. 2004, 82, 373–382. [Google Scholar] [CrossRef]
  22. Guinovart-Díaz, R.; Rodríguez-Ramos, R.; Bravo-Castillero, J.; Sabina, F.; Otero, J.; Maugin, G. A recursive asymptotic homogenization scheme for multi-phase fibrous elastic composites. Mech. Mater. 2005, 37, 1119–1131. [Google Scholar] [CrossRef]
  23. Castillero, J.; Otero, J.; Ramos, R.; Bourgeat, A. Asymptotic homogenization of laminated piezocomposite materials. Int. J. Solids Struct. 1998, 35, 527–541. [Google Scholar] [CrossRef]
  24. Yi, Y.M.; Park, S.H.; Youn, S.K. Asymptotic homogenization of viscoelastic composites with periodic microstructures. Int. J. Solids Struct. 1998, 35, 2039–2055. [Google Scholar] [CrossRef]
  25. Fantoni, F.; Bacigalupo, A.; Paggi, M. Multi-field asymptotic homogenization of thermo-piezoelectric materials with periodic microstructure. Int. J. Solids Struct. 2017, 120, 31–56. [Google Scholar] [CrossRef] [Green Version]
  26. Kanit, T.; Forest, S.; Galliet, I.; Mounoury, V.; Jeulin, D. Determination of the size of the representative volume element for random composites: Statistical and numerical approach. Int. J. Solids Struct. 2003, 40, 3647–3679. [Google Scholar] [CrossRef]
  27. Drygaś, P.; Mityushev, V. Effective elastic properties of random two-dimensional composites. Int. J. Solids Struct. 2016, 97-98, 543–553. [Google Scholar] [CrossRef]
  28. Gotlib, V.A.; Sato, T.; Beltzer, A.I. Neural computing of effective properties of random composite materials. Comput. Struct. 2001. [Google Scholar] [CrossRef]
  29. Levin, V.; Sabina, F.; Bravo-Castillero, J.; Guinovart-Díaz, R.; Rodríguez-Ramos, R.; Valdiviezo-Mijangos, O. Analysis of effective properties of electroelastic composites using the self-consistent and asymptotic homogenization methods. Int. J. Eng. Sci. 2008, 46, 818–834. [Google Scholar] [CrossRef]
  30. López-Realpozo, J.; Caballero-Pérez, R.; Rodríguez-Ramos, R.; Guinovart-Díaz, R.; Bravo-Castillero, J.; Camacho-Montes, H.; Espinosa-Almeyda, Y.; Sabina, F.; Chi-Vinh, P. Analysis of mechanical and electrical imperfect contacts in piezoelectric composites. Mech. Res. Commun. 2018, 93, 96–102. [Google Scholar] [CrossRef]
  31. Snarskii, A.; Bezsudnov, I.; Sevrukov, V.; Morozovskiy, A.; Malinsky, J. Transport Processes in Macroscopically Disordered Media: From Mean Field Theory to Percolation; Springer: New York, NY, USA, 2016. [Google Scholar]
  32. Yang, T.; Hu, T.; Liang, X.; Shen, S. On band structures of layered phononic crystals with flexoelectricity. Arch. Appl. Mech. 2018. [Google Scholar] [CrossRef]
  33. Della, C.N.; Shu, D. On the performance of 1–3 piezoelectric composites with a passive and active matrix. Sens. Actuators A Phys. 2007, 140, 200–206. [Google Scholar] [CrossRef]
  34. Erba, A.; El-Kelany, K.E.; Ferrero, M.; Baraille, I.; Rérat, M. Piezoelectricity of SrTiO3: An ab initio description. Phys. Rev. B 2013, 88, 035102. [Google Scholar] [CrossRef]
  35. Holmes-Siedle, A.; Wilson, P.; Verrall, A. PVdF: An electronically-active polymer for industry. Mater. Des. 1983, 4, 910–918. [Google Scholar] [CrossRef]
  36. Jain, A.; Sharma, A.K.; Jain, A. Dielectric and piezoelectric properties of PVDF/PZT composites: A review. Polym. Eng. Sci. 2015, 55, 1589–1616. [Google Scholar] [CrossRef] [Green Version]
  37. Berger, H.; Kari, S.; Gabbert, U.; Rodríguez-Ramos, R.; Bravo-Castillero, J.; Guinovart-Díaz, R. A comprehensive numerical homogenisation technique for calculating effective coefficients of uniaxial piezoelectric fibre composites. Mater. Sci. Eng. A 2005, 412, 53–60. [Google Scholar] [CrossRef]
  38. Berger, H.; Kari, S.; Gabbert, U.; Rodriguez-Ramos, R.; Bravo-Castillero, J.; Guinovart-Diaz, R. Calculation of effective coefficients for piezoelectric fiber composites based on a general numerical homogenization technique. Compos. Struct. 2005, 71, 397–400. [Google Scholar] [CrossRef]
  39. Klaasen, G.A. Existence Theorems for boundary value problems of n-th order ordinary differential equations. Rocky Mt. J. Math. 1973, 3, 457–472. [Google Scholar] [CrossRef]
  40. Maranganti, R.; Sharma, N.D.; Sharma, P. Electromechanical coupling in nonpiezoelectric materials due to nanoscale nonlocal size effects: Green’s function solutions and embedded inclusions. Phys. Rev. B 2006, 74, 014110. [Google Scholar] [CrossRef]
  41. Guinovart-Sanjuán, D.; Rizzoni, R.; Rodríguez-Ramos, R.; Guinovart-Díaz, R.; Bravo-Castillero, J.; Alfonso-Rodríguez, R.; Lebon, F.; Dumont, S.; Sevostianov, I.; Sabina, F.J. Behavior of laminated shell composite with imperfect contact between the layers. Compos. Struct. 2017, 176, 539–546. [Google Scholar] [CrossRef]
  42. Guinovart-Sanjuán, D.; Rodríguez-Ramos, R.; Guinovart-Díaz, R.; Bravo-Castillero, J.; Sabina, F.; Merodio, J.; Lebon, F.; Dumont, S.; Conci, A. Effective properties of regular elastic laminated shell composite. Compos. Part B Eng. 2016, 87, 12–20. [Google Scholar] [CrossRef] [Green Version]
  43. Lu, J.; Lv, J.; Liang, X.; Xu, M.; Shen, S. Improved approach to measure the direct flexoelectric coefficient of bulk polyvinylidene fluoride. J. Appl. Phys. 2016, 119, 094104. [Google Scholar] [CrossRef]
  44. Liu, K.; Zhang, S.; Xu, M.; Wu, T.; Shen, S. The research of effective flexoelectric coefficient along 1123 direction in polyvinylidene fluoride. J. Appl. Phys. 2017, 121, 174104. [Google Scholar] [CrossRef]
  45. Zhang, S.; Liu, K.; Xu, M.; Shen, H.; Chen, K.; Feng, B.; Shen, S. Investigation of the 2312 flexoelectric coefficient component of polyvinylidene fluoride: Deduction, simulation, and mensuration. Sci. Rep. 2017, 7, 3134. [Google Scholar] [CrossRef] [PubMed]
  46. Baskaran, S.; He, X.; Chen, Q.; Fu, J.Y. Experimental studies on the direct flexoelectric effect in α-phase polyvinylidene fluoride films. Appl. Phys. Lett. 2011, 98, 242901. [Google Scholar] [CrossRef]
  47. Bakhvalov, N.; Panasenko, G. Homogenisation: Averaging Processes in Periodic Media: Mathematical Problems in the Mechanics of Composite Materials; Mathematics and its Applications; Springer: Dordrecht, The Netherlands, 2012. [Google Scholar]
Figure 1. One-dimensional flexoelectric rod and the corresponding periodic cell.
Figure 1. One-dimensional flexoelectric rod and the corresponding periodic cell.
Materials 12 00232 g001
Figure 2. Effective properties and structures of three random flexoelectric composites with ideal periodicity.
Figure 2. Effective properties and structures of three random flexoelectric composites with ideal periodicity.
Materials 12 00232 g002
Figure 3. Values of different effective coefficients for three composites: i) flexoelectric material BaTiO 3 and the piezoelectric active polyvinylidene fluoride (PVDF) polymer; ii) composition of the flexoelectric composite SrTiO 3 and PVDF; and iii) the piezoelectric constituents PZT-7A and PVDF vs. values of volume fraction.
Figure 3. Values of different effective coefficients for three composites: i) flexoelectric material BaTiO 3 and the piezoelectric active polyvinylidene fluoride (PVDF) polymer; ii) composition of the flexoelectric composite SrTiO 3 and PVDF; and iii) the piezoelectric constituents PZT-7A and PVDF vs. values of volume fraction.
Materials 12 00232 g003
Figure 4. Values of different effective coefficient of three composites made of the combination of two flexoelectric material SrTiO 3 -BaTiO 3 , a non-piezoelectric polymer with SrTiO 3 , and piezoelectric PZT-7A with a flexoelectric material SrTiO 3 vs. values of volume fraction.
Figure 4. Values of different effective coefficient of three composites made of the combination of two flexoelectric material SrTiO 3 -BaTiO 3 , a non-piezoelectric polymer with SrTiO 3 , and piezoelectric PZT-7A with a flexoelectric material SrTiO 3 vs. values of volume fraction.
Materials 12 00232 g004
Figure 5. Values of the stress function σ , solutions of the heterogeneous (5)–(8) and homogeneous (47)–(50) problems for different external forces functions: f ( x ) = e x and f ( x ) = sin ( x ) .
Figure 5. Values of the stress function σ , solutions of the heterogeneous (5)–(8) and homogeneous (47)–(50) problems for different external forces functions: f ( x ) = e x and f ( x ) = sin ( x ) .
Materials 12 00232 g005
Table 1. Comparison of the effective coefficients C ^ , e ^ , κ ^ for a piezoelectric composite computed using the methodology described in Section 3, the finite element method (FEM), and the results reported in References [19,30].
Table 1. Comparison of the effective coefficients C ^ , e ^ , κ ^ for a piezoelectric composite computed using the methodology described in Section 3, the finite element method (FEM), and the results reported in References [19,30].
Effective Coefficient C ^ (10 9 N/m 2 )
V 1 Present ModelFEMRef. [30]Ref. [19]
0.0118.300118.300118.299118.300
0.4116.953116.953116.952116.953
0.8115.401115.401115.400115.401
1.0114.500114.500114.499114.500
Effective Coefficient e ^ (C/m 2 )
0.026.40026.40026.39926.400
0.424.80624.80624.80524.806
0.822.63422.86322.63322.634
1.021.20021.20021.19921.200
Effective Coefficient κ ^ (10 10 F/m)
0.0110.900110.900110.899110.900
0.4148.056148.056148.054148.056
0.8200.835200.836200.832200.835
1.0236.600236.600236.595236.600
Table 2. Material properties of barium titanate (BaTiO 3 ), strontium titanate (SrTiO 3 ), PVDF, non-piezoelectric polymer, and PZT-7A are given. C, e, κ , and μ denote the elastic, piezoelectric, dielectric and flexoelectric tensors, respectively. The materials with μ = 0 are called non-flexoelectric (PVDF, polymer, PZT-7A).
Table 2. Material properties of barium titanate (BaTiO 3 ), strontium titanate (SrTiO 3 ), PVDF, non-piezoelectric polymer, and PZT-7A are given. C, e, κ , and μ denote the elastic, piezoelectric, dielectric and flexoelectric tensors, respectively. The materials with μ = 0 are called non-flexoelectric (PVDF, polymer, PZT-7A).
MaterialsC (10 9 N/m 2 )e (C/m 2 ) κ (10 10 F/m) μ (10 4 C/m)
BaTiO 3 162 [32]17.36 [33]4000 [32]5 [32]
SrTiO 3 350 [32]8.82 [34]300 [32]1 [32]
PVDF2.7 [35]30 [36]10 [36]0
Polymer3.86 [37]0 [37]0.7965 [37]0
PZT-7A131.4 [38]9.522 [38]0.372 [38]0
Table 3. Coefficients of the homogenized flexoelectric composite made of barium titanate (BaTiO 3 ) and strontium titanate (SrTiO 3 ) for V 1 = 0.5 and V 1 = 0.8 .
Table 3. Coefficients of the homogenized flexoelectric composite made of barium titanate (BaTiO 3 ) and strontium titanate (SrTiO 3 ) for V 1 = 0.5 and V 1 = 0.8 .
Materials C ^ (10 9 N/m 2 ) e ^ (C/m 2 ) κ ^ (10 10 F/m) μ ^ (10 4 C/m)
V 1 = 0.5 221.86646.8629560.041.6667
V 1 = 0.8 181.83748.02881158.112.7778

Share and Cite

MDPI and ACS Style

Guinovart-Sanjuán, D.; Merodio, J.; López-Realpozo, J.C.; Vajravelu, K.; Rodríguez-Ramos, R.; Guinovart-Díaz, R.; Bravo-Castillero, J.; Sabina, F.J. Asymptotic Homogenization Applied to Flexoelectric Rods. Materials 2019, 12, 232. https://doi.org/10.3390/ma12020232

AMA Style

Guinovart-Sanjuán D, Merodio J, López-Realpozo JC, Vajravelu K, Rodríguez-Ramos R, Guinovart-Díaz R, Bravo-Castillero J, Sabina FJ. Asymptotic Homogenization Applied to Flexoelectric Rods. Materials. 2019; 12(2):232. https://doi.org/10.3390/ma12020232

Chicago/Turabian Style

Guinovart-Sanjuán, David, Jose Merodio, Juan Carlos López-Realpozo, Kuppalapalle Vajravelu, Reinaldo Rodríguez-Ramos, Raúl Guinovart-Díaz, Julián Bravo-Castillero, and Federico J. Sabina. 2019. "Asymptotic Homogenization Applied to Flexoelectric Rods" Materials 12, no. 2: 232. https://doi.org/10.3390/ma12020232

APA Style

Guinovart-Sanjuán, D., Merodio, J., López-Realpozo, J. C., Vajravelu, K., Rodríguez-Ramos, R., Guinovart-Díaz, R., Bravo-Castillero, J., & Sabina, F. J. (2019). Asymptotic Homogenization Applied to Flexoelectric Rods. Materials, 12(2), 232. https://doi.org/10.3390/ma12020232

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop