Contents     Functions     Index          PDF

m_piezo, p_piezo

• 1 : Simplified 3D piezoelectric properties
• 2 : General 3D piezo
• 3 : General 3D piezo, e matrix
• Shell element properties

Purpose

Material function for piezoelectric solids

Syntax

 mat= m_piezo('database name')
 pl = m_piezo('dbval MatId name');

Sample calls are given at the end of this page and in section 6.6.

Material

The m_piezo Const and BuildConstit commands support integration constant building for piezo electric volumes integrated in the standard volume elements. Element properties are given by p_solid entries, while materials formats are detailed here. Theoretical details are given in  section 6.1.5.

The elastic constitutive law is declared using ElasMatId referencing an elastic material declared as a row of model.pl with this identifier, see m_elastic for input formats. The piezoelectric constants can be declared using the following sub-types

1 : Simplified 3D piezoelectric properties

[ProId Type ElasMatId d31 d32 d33 eps1T eps2T eps3T EDType]

These simplified piezoelectric properties can be used for PVDF, but also for PZT if shear mode actuation/sensing is not considered (d₂₄=d₁₅=0). For EDType==0 on assumes d is given. For EDType==1, e is given. Note that the values of є^T (permitivity at zero stress) should be given (and not є^S).

2 : General 3D piezo

[ProId Type ElasMatId d_1:18 epsT_1:9]

d_1:18 are the 18 constants of the [d] matrix, and epsT_1:9 are the 9 constants of the [є^T] matrix. One reminds that strains are stored in order xx,yy,zz,yz,zx,yx.

3 : General 3D piezo, e matrix

[ProId Type ElasMatId e_1:18 epsT_1:9]

e_1:18 are the 18 constants of the [d] matrix, and epsT_1:9 are the 9 constants of the [є^T] matrix.

Examples are given in section 6.6.

Shell element properties

Piezo shell elements with electrodes are associated with element property rows of the form

[ProId Type UnderlyingProId ElNodeId1 LayerId1 UNU1 ElNodeId2...]

Type=fe_mat('p_piezo','SI',1)

UnderlyingProId : ProId of underlying element p_shell composite entry piezo layers must be associated with piezo electric material properties as detailed above. ElNodId1 : NodeId for electrode 1 LayerId : layer number as declared in the composite entry UNU1 : unused property for now (angle for polarization)

Examples are given in section 6.6.

Shell strain is defined by the membrane, curvature and transverse shear as well as the electric field components. It is assumed that in each piezoelectric layer i=1...n, the electric field takes the form E^→= (0    0    E_zi). E_zi is assumed to be constant over the thickness h_i of the layer and is therefore given by E_zi=−Δ φ_i/h_i where Δ φ_i is the difference of potential between the electrodes at the top and bottom of the piezoelectric layer i. It is also assumed that the piezoelectric principal axes are parallel to the structural orthotropy axes.

The strain state of a piezoelectric shell takes the form

There are thus n additional degrees of freedom Δ φ_i, n being the number of piezoelectric layers in the laminate shell

The constitutive law associated to this strain state is given by :

where D_zi is the electric displacement in piezoelectric layer (assumed constant and in the z-direction), z_mi is the distance between the midplane of the shell and the midplane of piezoelectric layer i, and G_i, H_i are given by :

. denotes the direction of polarization. If the piezoelectric is used in extension mode, the polarization is in the z-direction, therefore H_i =0 and G_i ={

}_i . If the piezoelectric is used in shear mode, the polarization is in the x or y-direction, therefore G_i=0, and H_i = {

}_i or H_i = {

}_i . It turns out however that the hypothesis of a uniform transverse shear strain distribution through the thickness is not satisfactory, a more elaborate shell element would be necessary. Shear actuation should therefore be used with caution.

[R_s]_i and [R]_i are rotation matrices associated to the angle θ of the piezoelectric layer.

 model=femesh('testquad4'); 
 model.pl=m_elastic('dbval 1 steel','dbval 2 steel', ...
 'dbval 12 ',[12 fe_mat('m_elastic','SI',1) 65e9 .3 7800]) % PZT elastic prop
 model.pl=m_piezo(model.pl,'dbval 3 -elas 12 Sample_ULB');
 model.il=p_shell('dbval 100 kirchoff .05 -f5','dbval 1 kirchoff .05 -f5');
 % 3 layer composite, .01 mm piezo, .2 mm steel, .01mm piezo
 model.il=p_shell(model.il,'dbval 111 laminate 3 1e-3 0 1 .05 0 3 1e-3 0');
 % piezo shell : 3 layer
 model.il=p_piezo(model.il,'dbval 110 shell 111 1682 1 0 1683 3 0');

 sdtdef('diag',2) % show constitutive properties
 [Case,model.DOF]=fe_mknl('init',model);fe_c(model.DOF,.21,'dofs')
 sdtdef('diag',0);