1

Satyajit
Date Sep 26, 2014 4:13:14 PM

2

Contents
1. Component 1 (comp1) ..................................................................................................................... 3
1.1. Definitions ................................................................................................................................ 3
1.2. Geometry 1 .............................................................................................................................. 4
1.3. Materials .................................................................................................................................. 6
1.4. Piezoelectric Devices (pzd) ..................................................................................................... 10
1.5. Mesh 1.................................................................................................................................... 86
2. Study 1 ........................................................................................................................................... 88
2.1. Stationary ............................................................................................................................... 88
2.2. Solver Configurations ............................................................................................................. 88
3. Results ............................................................................................................................................ 90
3.1. Data Sets ................................................................................................................................ 90
3.2. Plot Groups ............................................................................................................................ 91

3

1 Component 1 (comp1)
1.1 Definitions
1.1.1 Coordinate Systems
Boundary System 1
Coordinate system type Boundary system
Identifier sys1

Settings
Name Value
Coordinate names {t1, t2, n}
Create first tangent direction from Global Cartesian (spatial)

Base Vector System 2
Coordinate system type Base vector system
Identifier sys2

Settings
Name Value
Coordinate names {x1, x2, x3}
transposedCoord {x1, x2, x3}
Base vectors {{0, 0, 1}, {0, 1, 0}, {1, 0, 0}}

4

1.2 Geometry 1

Geometry 1
Units
Length unit m
Angular unit deg

Geometry statistics
Property Value
Space dimension 3
Number of domains 5
Number of boundaries 26
Number of edges 44
Number of vertices 24

1.2.1 Block 1 (blk1)
Position
Name Value
Position {0, 0, 0}

Axis
Name Value
Axis type z - axis

Size and shape
5

Name Value
Size {0.1, 0.03, 0.018}

1.2.2 Block 2 (blk2)
Position
Name Value
Position {0, 0, 0.008}

Axis
Name Value
Axis type z - axis

Size and shape
Name Value
Size {0.1, 0.03, 0.002}

1.2.3 Block 3 (blk3)
Position
Name Value
Position {0.055, 0, 0.008}

Axis
Name Value
Axis type z - axis

Size and shape
Name Value
Size {0.01, 0.03, 0.002}

6

1.3 Materials
1.3.1 Al

Al
Selection
Geometric entity level Domain
Selection Domains 1, 3

Material parameters
Name Value Unit
Density 2700[kg/m^3] kg/m^3
Young's modulus 70.0e9[Pa] Pa
Poisson's ratio 0.35 1

Basic Settings
Description Value
Electrical conductivity {{35.5e6[S/m], 0, 0}, {0, 35.5e6[S/m], 0}, {0, 0,
35.5e6[S/m]}}
Coefficient of thermal expansion {{23.1e-6[1/K], 0, 0}, {0, 23.1e-6[1/K], 0}, {0, 0, 23.1e-
6[1/K]}}
Heat capacity at constant pressure 904[J/(kg*K)]
Density 2700[kg/m^3]
Thermal conductivity {{237[W/(m*K)], 0, 0}, {0, 237[W/(m*K)], 0}, {0, 0,
237[W/(m*K)]}}

7

Young's modulus and Poisson's ratio Settings
Description Value
Young's modulus 70.0e9[Pa]
Poisson's ratio 0.35

1.3.2 Lead Zirconate Titanate (PZT-5H)

Lead Zirconate Titanate (PZT-5H)
Selection
Geometric entity level Domain
Selection Domain 4

Material parameters
Name Value Unit
Density 7500[kg/m^3] kg/m^3
Elasticity matrix (Ordering: xx, yy, zz, yz, xz,
xy)
{1.27205e+011[Pa], 8.02122e+010[Pa],
1.27205e+011[Pa], 8.46702e+010[Pa],
8.46702e+010[Pa], 1.17436e+011[Pa], 0[Pa],
0[Pa], 0[Pa], 2.29885e+010[Pa], 0[Pa], 0[Pa],
0[Pa], 0[Pa], 2.29885e+010[Pa], 0[Pa], 0[Pa],
0[Pa], 0[Pa], 0[Pa], 2.34742e+010[Pa]}
Pa
Coupling matrix {0[C/m^2], 0[C/m^2], -6.62281[C/m^2],
0[C/m^2], 0[C/m^2], -6.62281[C/m^2],
0[C/m^2], 0[C/m^2], 23.2403[C/m^2],
0[C/m^2], 17.0345[C/m^2], 0[C/m^2],
17.0345[C/m^2], 0[C/m^2], 0[C/m^2],
0[C/m^2], 0[C/m^2], 0[C/m^2]}
C/m^2
8

Name Value Unit
Relative permittivity {1704.4, 1704.4, 1433.6} 1

Basic Settings
Description Value
Density 7500[kg/m^3]

Strain-charge form Settings
Description Value
Compliance matrix (ordering: xx, yy, zz, yz, xz, xy) {{1.65e-011[1/Pa], -4.78e-012[1/Pa], -8.45e-
012[1/Pa], 0[1/Pa], 0[1/Pa], 0[1/Pa]}, {-4.78e-
012[1/Pa], 1.65e-011[1/Pa], -8.45e-012[1/Pa],
0[1/Pa], 0[1/Pa], 0[1/Pa]}, {-8.45e-012[1/Pa], -
8.45e-012[1/Pa], 2.07e-011[1/Pa], 0[1/Pa],
0[1/Pa], 0[1/Pa]}, {0[1/Pa], 0[1/Pa], 0[1/Pa],
4.35e-011[1/Pa], 0[1/Pa], 0[1/Pa]}, {0[1/Pa],
0[1/Pa], 0[1/Pa], 0[1/Pa], 4.35e-011[1/Pa],
0[1/Pa]}, {0[1/Pa], 0[1/Pa], 0[1/Pa], 0[1/Pa],
0[1/Pa], 4.26e-011[1/Pa]}}
Coupling matrix (ordering: xx, yy, zz, yz, xz, xy) {{0[C/N], 0[C/N], -2.74e-010[C/N], 0[C/N], 0[C/N],
-2.74e-010[C/N]}, {0[C/N], 0[C/N], 5.93e-
010[C/N], 0[C/N], 7.41e-010[C/N], 0[C/N]},
{7.41e-010[C/N], 0[C/N], 0[C/N], 0[C/N], 0[C/N],
0[C/N]}}
Relative permittivity {{3130, 0, 0}, {0, 3130, 0}, {0, 0, 3400}}
Loss factor for compliance matrix sE {{0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0}}
Loss factor for coupling matrix d {{0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0}}
Loss factor for electrical permittivity T {{0, 0, 0}, {0, 0, 0}, {0, 0, 0}}

Stress-charge form Settings
Description Value
Elasticity matrix (Ordering: xx, yy, zz, yz, xz, xy) {{1.27205e+011[Pa], 8.02122e+010[Pa],
8.46702e+010[Pa], 0[Pa], 0[Pa], 0[Pa]},
{8.02122e+010[Pa], 1.27205e+011[Pa],
8.46702e+010[Pa], 0[Pa], 0[Pa], 0[Pa]},
{8.46702e+010[Pa], 8.46702e+010[Pa],
1.17436e+011[Pa], 0[Pa], 0[Pa], 0[Pa]}, {0[Pa],
0[Pa], 0[Pa], 2.29885e+010[Pa], 0[Pa], 0[Pa]},
{0[Pa], 0[Pa], 0[Pa], 0[Pa], 2.29885e+010[Pa],
0[Pa]}, {0[Pa], 0[Pa], 0[Pa], 0[Pa], 0[Pa],
2.34742e+010[Pa]}}
9

Description Value
Coupling matrix {{0[C/m^2], 0[C/m^2], -6.62281[C/m^2], 0[C/m^2],
0[C/m^2], -6.62281[C/m^2]}, {0[C/m^2], 0[C/m^2],
23.2403[C/m^2], 0[C/m^2], 17.0345[C/m^2],
0[C/m^2]}, {17.0345[C/m^2], 0[C/m^2], 0[C/m^2],
0[C/m^2], 0[C/m^2], 0[C/m^2]}}
Relative permittivity {{1704.4, 0, 0}, {0, 1704.4, 0}, {0, 0, 1433.6}}
Loss factor for elasticity matrix cE {{0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0}, {0,
0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0}}
Loss factor for coupling matrix e {{0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0}}
Loss factor for electrical permittivity S {{0, 0, 0}, {0, 0, 0}, {0, 0, 0}}

1.3.3 Material 3

Material 3
Selection
Geometric entity level Domain
Selection Domains 2, 5

Material parameters
Name Value Unit
Young's modulus 35.3e6 Pa
Poisson's ratio 0.383 1
Density 32 kg/m^3

Basic Settings
10

Description Value
Young's modulus 35.3e6
Poisson's ratio 0.383
Density 32

1.4 Piezoelectric Devices (pzd)

Piezoelectric Devices
Selection
Geometric entity level Domain
Selection Domains 15

Equations


Settings
Description Value
Displacement field Quadratic
Compute boundary fluxes Off
Electric potential Quadratic
Value type when using splitting of complex variables {Complex, Complex}
Frame Material
Displacements control spatial frame On
11

Description Value
Reference point for moment computation {0, 0, 0}
Add mixed form pressure everywhere Off
Structural transient behavior Include inertial terms
SmeApplModeID SolidMechanics
EquationType smetffs
Activate terminal sweep Off
Reference impedance 50[ohm]
Parameter to export Z
Typical wave speed for perfectly matched layers pzd.cp

Used products
COMSOL Multiphysics
MEMS Module

Variables
Name Expression Unit Description Selection
pzd.nX nX Normal vector, X
component
Boundaries 126
pzd.nY nY Normal vector, Y
component
Boundaries 126
pzd.nZ nZ Normal vector, Z
component
Boundaries 126
pzd.nXmesh root.nXmesh Normal vector (mesh), X
component
Boundaries 126
pzd.nYmesh root.nYmesh Normal vector (mesh), Y
component
Boundaries 126
pzd.nZmesh root.nZmesh Normal vector (mesh), Z
component
Boundaries 126
pzd.nx nx Normal vector, x
component
Boundaries 126
pzd.ny ny Normal vector, y
component
Boundaries 126
pzd.nz nz Normal vector, z
component
Boundaries 126
pzd.nxmesh root.nxmesh Normal vector (mesh), x
component
Boundaries 126
12

Name Expression Unit Description Selection
pzd.nymesh root.nymesh Normal vector (mesh), y
component
Boundaries 126
pzd.nzmesh root.nzmesh Normal vector (mesh), z
component
Boundaries 126
xt d(x,TIME) m/s Mesh velocity, x component Global
yt d(y,TIME) m/s Mesh velocity, y
component
Global
zt d(z,TIME) m/s Mesh velocity, z component Global
pzd.refpntx 0 m Reference point for
moment computation, x
coordinate
Global
pzd.refpnty 0 m Reference point for
moment computation, y
coordinate
Global
pzd.refpntz 0 m Reference point for
moment computation, z
coordinate
Global
pzd.zref 50[ohm] Reference impedance Global
pzd.cref pzd.cp m/s Typical wave speed for
perfectly matched layers
Global

1.4.1 Piezoelectric Material 1

Piezoelectric Material 1
Selection
13

Geometric entity level Domain
Selection Domain 4

Equations







Settings
Settings
Description Value
Elasticity matrix (Ordering: xx, yy, zz, yz, xz, xy) From material
Elasticity matrix (Ordering: xx, yy, zz, yz, xz, xy) {{0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0}}
Coupling matrix From material
Relative permittivity From material
Relative permittivity {{0, 0, 0}, {0, 0, 0}, {0, 0, 0}}
Compliance matrix (ordering: xx, yy, zz, yz, xz, xy) {{0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0},
{0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0}}
Relative permittivity {{0, 0, 0}, {0, 0, 0}, {0, 0, 0}}
Density From material
Constitutive relation Stress - charge form

Properties from material
Property Material Property group
Elasticity matrix (Ordering: xx, yy, zz,
yz, xz, xy)
Lead Zirconate Titanate (PZT-5H) Stress-charge form
Coupling matrix Lead Zirconate Titanate (PZT-5H) Stress-charge form
Relative permittivity Lead Zirconate Titanate (PZT-5H) Stress-charge form
Density Lead Zirconate Titanate (PZT-5H) Basic

Variables
Name Expression Unit Description Selection
14

Name Expression Unit Description Selection
uTXt root.comp1.uTXTIME 1/s Tangential
gradient of u,
first time
derivative, X
component
Boundaries
1419
uTYt root.comp1.uTYTIME 1/s Tangential
gradient of u,
first time
derivative, Y
component
Boundaries
1419
uTZt root.comp1.uTZTIME 1/s Tangential
gradient of u,
first time
derivative, Z
component
Boundaries
1419
uTXtt root.comp1.uTXTIMETI
ME
1/s^2 Tangential
gradient of u,
second time
derivative, X
component
Boundaries
1419
uTYtt root.comp1.uTYTIMETI
ME
1/s^2 Tangential
gradient of u,
second time
derivative, Y
component
Boundaries
1419
uTZtt root.comp1.uTZTIMETI
ME
1/s^2 Tangential
gradient of u,
second time
derivative, Z
component
Boundaries
1419
uXt root.comp1.uXTIME 1/s Gradient of u,
first time
derivative, X
component
Domain 4
uYt root.comp1.uYTIME 1/s Gradient of u,
first time
derivative, Y
component
Domain 4
uZt root.comp1.uZTIME 1/s Gradient of u,
first time
derivative, Z
component
Domain 4
15

Name Expression Unit Description Selection
uXtt root.comp1.uXTIMETIM
E
1/s^2 Gradient of u,
second time
derivative, X
component
Domain 4
uYtt root.comp1.uYTIMETIM
E
1/s^2 Gradient of u,
second time
derivative, Y
component
Domain 4
uZtt root.comp1.uZTIMETIM
E
1/s^2 Gradient of u,
second time
derivative, Z
component
Domain 4
vTXt root.comp1.vTXTIME 1/s Tangential
gradient of v,
first time
derivative, X
component
Boundaries
1419
vTYt root.comp1.vTYTIME 1/s Tangential
gradient of v,
first time
derivative, Y
component
Boundaries
1419
vTZt root.comp1.vTZTIME 1/s Tangential
gradient of v,
first time
derivative, Z
component
Boundaries
1419
vTXtt root.comp1.vTXTIMETI
ME
1/s^2 Tangential
gradient of v,
second time
derivative, X
component
Boundaries
1419
vTYtt root.comp1.vTYTIMETI
ME
1/s^2 Tangential
gradient of v,
second time
derivative, Y
component
Boundaries
1419
vTZtt root.comp1.vTZTIMETI
ME
1/s^2 Tangential
gradient of v,
second time
derivative, Z
component
Boundaries
1419
16

Name Expression Unit Description Selection
vXt root.comp1.vXTIME 1/s Gradient of v,
first time
derivative, X
component
Domain 4
vYt root.comp1.vYTIME 1/s Gradient of v,
first time
derivative, Y
component
Domain 4
vZt root.comp1.vZTIME 1/s Gradient of v,
first time
derivative, Z
component
Domain 4
vXtt root.comp1.vXTIMETIM
E
1/s^2 Gradient of v,
second time
derivative, X
component
Domain 4
vYtt root.comp1.vYTIMETIM
E
1/s^2 Gradient of v,
second time
derivative, Y
component
Domain 4
vZtt root.comp1.vZTIMETIM
E
1/s^2 Gradient of v,
second time
derivative, Z
component
Domain 4
wTXt root.comp1.wTXTIME 1/s Tangential
gradient of w,
first time
derivative, X
component
Boundaries
1419
wTYt root.comp1.wTYTIME 1/s Tangential
gradient of w,
first time
derivative, Y
component
Boundaries
1419
wTZt root.comp1.wTZTIME 1/s Tangential
gradient of w,
first time
derivative, Z
component
Boundaries
1419
wTXtt root.comp1.wTXTIMETI
ME
1/s^2 Tangential
gradient of w,
second time
derivative, X
Boundaries
1419
17

Name Expression Unit Description Selection
component
wTYtt root.comp1.wTYTIMETI
ME
1/s^2 Tangential
gradient of w,
second time
derivative, Y
component
Boundaries
1419
wTZtt root.comp1.wTZTIMETI
ME
1/s^2 Tangential
gradient of w,
second time
derivative, Z
component
Boundaries
1419
wXt root.comp1.wXTIME 1/s Gradient of w,
first time
derivative, X
component
Domain 4
wYt root.comp1.wYTIME 1/s Gradient of w,
first time
derivative, Y
component
Domain 4
wZt root.comp1.wZTIME 1/s Gradient of w,
first time
derivative, Z
component
Domain 4
wXtt root.comp1.wXTIMETI
ME
1/s^2 Gradient of w,
second time
derivative, X
component
Domain 4
wYtt root.comp1.wYTIMETI
ME
1/s^2 Gradient of w,
second time
derivative, Y
component
Domain 4
wZtt root.comp1.wZTIMETI
ME
1/s^2 Gradient of w,
second time
derivative, Z
component
Domain 4
ut root.comp1.uTIME m/s Structural
velocity field, X
component
Domain 4
vt root.comp1.vTIME m/s Structural
velocity field, Y
component
Domain 4
18

Name Expression Unit Description Selection
wt root.comp1.wTIME m/s Structural
velocity field, Z
component
Domain 4
utt root.comp1.uTIMETIME m/s^2 Acceleration
field, X
component
Domain 4
vtt root.comp1.vTIMETIME m/s^2 Acceleration
field, Y
component
Domain 4
wtt root.comp1.wTIMETIM
E
m/s^2 Acceleration
field, Z
component
Domain 4
VTXt root.comp1.VTXTIME m*kg/(s^4*A) Tangential
gradient of V,
first time
derivative, X
component
Boundaries
1419
VTYt root.comp1.VTYTIME m*kg/(s^4*A) Tangential
gradient of V,
first time
derivative, Y
component
Boundaries
1419
VTZt root.comp1.VTZTIME m*kg/(s^4*A) Tangential
gradient of V,
first time
derivative, Z
component
Boundaries
1419
VTXtt root.comp1.VTXTIMETI
ME
m*kg/(s^5*A) Tangential
gradient of V,
second time
derivative, X
component
Boundaries
1419
VTYtt root.comp1.VTYTIMETI
ME
m*kg/(s^5*A) Tangential
gradient of V,
second time
derivative, Y
component
Boundaries
1419
VTZtt root.comp1.VTZTIMETI
ME
m*kg/(s^5*A) Tangential
gradient of V,
second time
derivative, Z
component
Boundaries
1419
19

Name Expression Unit Description Selection
Vt root.comp1.VTIME-
Vx*d(x,TIME)-
Vy*d(y,TIME)-
Vz*d(z,TIME)
V/s Electric
potential, first
time derivative
Domain 4
Vtt d(root.comp1.VTIME-
Vx*d(x,TIME)-
Vy*d(y,TIME)-
Vz*d(z,TIME),TIME)-
d(root.comp1.VTIME-
Vx*d(x,TIME)-
Vy*d(y,TIME)-
Vz*d(z,TIME),x)*d(x,TI
ME)-
d(root.comp1.VTIME-
Vx*d(x,TIME)-
Vy*d(y,TIME)-
Vz*d(z,TIME),y)*d(y,TI
ME)-
d(root.comp1.VTIME-
Vx*d(x,TIME)-
Vy*d(y,TIME)-
Vz*d(z,TIME),z)*d(z,TI
ME)
V/s^2 Electric
potential, second
time derivative
Domain 4
pzd.Eequ (pzd.cE44+pzd.cE55+pz
d.cE66)*(pzd.cE11+2*p
zd.cE12+2*pzd.cE13+pz
d.cE22+2*pzd.cE23+pzd
.cE33)/(pzd.cE11+2*pzd
.cE12+2*pzd.cE13+pzd.
cE22+2*pzd.cE23+pzd.c
E33+pzd.cE44+pzd.cE55
+pzd.cE66)
Pa Equivalent
Young's modulus
Domain 4
pzd.nuequ 1.5*((pzd.cE11+2*pzd.c
E12+2*pzd.cE13+pzd.cE
22+2*pzd.cE23+pzd.cE3
3)/3-
2*(pzd.cE44+pzd.cE55+
pzd.cE66)/3)/(pzd.cE11
+2*pzd.cE12+2*pzd.cE1
3+pzd.cE22+2*pzd.cE23
+pzd.cE33+pzd.cE44+pz
d.cE55+pzd.cE66)
1 Equivalent
Poisson's ratio
Domain 4
pzd.cp sqrt(((pzd.cE11+2*pzd.c
E12+2*pzd.cE13+pzd.cE
22+2*pzd.cE23+pzd.cE3
m/s Pressure-wave
speed
Domain 4
20

Name Expression Unit Description Selection
3)/9+0.4444444444444
444*(pzd.cE44+pzd.cE5
5+pzd.cE66))/pzd.rho)
pzd.rho material.rho kg/m^3 Density Domain 4
pzd.FdxX 1+uX 1 Deformation
gradient, xX
component
Domain 4
pzd.FdyX vX 1 Deformation
gradient, yX
component
Domain 4
pzd.FdzX wX 1 Deformation
gradient, zX
component
Domain 4
pzd.FdxY uY 1 Deformation
gradient, xY
component
Domain 4
pzd.FdyY 1+vY 1 Deformation
gradient, yY
component
Domain 4
pzd.FdzY wY 1 Deformation
gradient, zY
component
Domain 4
pzd.FdxZ uZ 1 Deformation
gradient, xZ
component
Domain 4
pzd.FdyZ vZ 1 Deformation
gradient, yZ
component
Domain 4
pzd.FdzZ 1+wZ 1 Deformation
gradient, zZ
component
Domain 4
pzd.Fdlx1 (1+uX)*sys2.T11+uY*sy
s2.T12+uZ*sys2.T13
1 Deformation
gradient, local,
x1 component
Domain 4
pzd.Fdly1 vX*sys2.T11+(1+vY)*sys
2.T12+vZ*sys2.T13
1 Deformation
gradient, local,
y1 component
Domain 4
pzd.Fdlz1 wX*sys2.T11+wY*sys2.
T12+(1+wZ)*sys2.T13
1 Deformation
gradient, local,
z1 component
Domain 4
21

Name Expression Unit Description Selection
pzd.Fdlx2 (1+uX)*sys2.T21+uY*sy
s2.T22+uZ*sys2.T23
1 Deformation
gradient, local,
x2 component
Domain 4
pzd.Fdly2 vX*sys2.T21+(1+vY)*sys
2.T22+vZ*sys2.T23
1 Deformation
gradient, local,
y2 component
Domain 4
pzd.Fdlz2 wX*sys2.T21+wY*sys2.
T22+(1+wZ)*sys2.T23
1 Deformation
gradient, local,
z2 component
Domain 4
pzd.Fdlx3 (1+uX)*sys2.T31+uY*sy
s2.T32+uZ*sys2.T33
1 Deformation
gradient, local,
x3 component
Domain 4
pzd.Fdly3 vX*sys2.T31+(1+vY)*sys
2.T32+vZ*sys2.T33
1 Deformation
gradient, local,
y3 component
Domain 4
pzd.Fdlz3 wX*sys2.T31+wY*sys2.
T32+(1+wZ)*sys2.T33
1 Deformation
gradient, local,
z3 component
Domain 4
pzd.gradUxX uX 1 Displacement
gradient, xX
component
Domain 4
pzd.gradUyX vX 1 Displacement
gradient, yX
component
Domain 4
pzd.gradUzX wX 1 Displacement
gradient, zX
component
Domain 4
pzd.gradUxY uY 1 Displacement
gradient, xY
component
Domain 4
pzd.gradUyY vY 1 Displacement
gradient, yY
component
Domain 4
pzd.gradUzY wY 1 Displacement
gradient, zY
component
Domain 4
pzd.gradUxZ uZ 1 Displacement
gradient, xZ
component
Domain 4
pzd.gradUyZ vZ 1 Displacement
gradient, yZ
Domain 4
22

Name Expression Unit Description Selection
component
pzd.gradUzZ wZ 1 Displacement
gradient, zZ
component
Domain 4
pzd.eXX uX 1 Strain tensor, XX
component
Domain 4
pzd.eXY 0.5*(uY+vX) 1 Strain tensor, XY
component
Domain 4
pzd.eXZ 0.5*(uZ+wX) 1 Strain tensor, XZ
component
Domain 4
pzd.eYY vY 1 Strain tensor, YY
component
Domain 4
pzd.eYZ 0.5*(vZ+wY) 1 Strain tensor, YZ
component
Domain 4
pzd.eZZ wZ 1 Strain tensor, ZZ
component
Domain 4
pzd.eX pzd.eXX 1 Strain tensor, X
component
Domain 4
pzd.eY pzd.eYY 1 Strain tensor, Y
component
Domain 4
pzd.eZ pzd.eZZ 1 Strain tensor, Z
component
Domain 4
pzd.Ldx d(pzd.gradUxX,TIME) 1/s Rate of strain
tensor, x
component
Domain 4
pzd.Ldxy 0.5*(d(pzd.gradUxY,TI
ME)+d(pzd.gradUyX,TI
ME))
1/s Rate of strain
tensor, xy
component
Domain 4
pzd.Ldxz 0.5*(d(pzd.gradUxZ,TI
ME)+d(pzd.gradUzX,TI
ME))
1/s Rate of strain
tensor, xz
component
Domain 4
pzd.Ldy d(pzd.gradUyY,TIME) 1/s Rate of strain
tensor, y
component
Domain 4
pzd.Ldyz 0.5*(d(pzd.gradUyZ,TI
ME)+d(pzd.gradUzY,TI
ME))
1/s Rate of strain
tensor, yz
component
Domain 4
pzd.Ldz d(pzd.gradUzZ,TIME) 1/s Rate of strain
tensor, z
Domain 4
23

Name Expression Unit Description Selection
component
pzd.Lwx 0 1/s Spin tensor, x
component
Domain 4
pzd.Lwxy 0.5*(d(pzd.gradUxY,TI
ME)-
d(pzd.gradUyX,TIME))
1/s Spin tensor, xy
component
Domain 4
pzd.Lwxz 0.5*(d(pzd.gradUxZ,TI
ME)-
d(pzd.gradUzX,TIME))
1/s Spin tensor, xz
component
Domain 4
pzd.Lwy 0 1/s Spin tensor, y
component
Domain 4
pzd.Lwyz 0.5*(d(pzd.gradUyZ,TI
ME)-
d(pzd.gradUzY,TIME))
1/s Spin tensor, yz
component
Domain 4
pzd.Lwz 0 1/s Spin tensor, z
component
Domain 4
pzd.el11 (sys2.T11*uX+0.5*sys2.
T12*(uY+vX)+0.5*sys2.
T13*(uZ+wX))*sys2.T11
+(0.5*sys2.T11*(uY+vX)
+sys2.T12*vY+0.5*sys2.
T13*(vZ+wY))*sys2.T12
+(0.5*sys2.T11*(uZ+wX
)+0.5*sys2.T12*(vZ+wY)
+sys2.T13*wZ)*sys2.T1
3
1 Strain tensor,
local coordinate
system, 11
component
Domain 4
pzd.el12 (sys2.T11*uX+0.5*sys2.
T12*(uY+vX)+0.5*sys2.
T13*(uZ+wX))*sys2.T21
+(0.5*sys2.T11*(uY+vX)
+sys2.T12*vY+0.5*sys2.
T13*(vZ+wY))*sys2.T22
+(0.5*sys2.T11*(uZ+wX
)+0.5*sys2.T12*(vZ+wY)
+sys2.T13*wZ)*sys2.T2
3
1 Strain tensor,
local coordinate
system, 12
component
Domain 4
pzd.el13 (sys2.T11*uX+0.5*sys2.
T12*(uY+vX)+0.5*sys2.
T13*(uZ+wX))*sys2.T31
+(0.5*sys2.T11*(uY+vX)
+sys2.T12*vY+0.5*sys2.
T13*(vZ+wY))*sys2.T32
+(0.5*sys2.T11*(uZ+wX
1 Strain tensor,
local coordinate
system, 13
component
Domain 4
24

Name Expression Unit Description Selection
)+0.5*sys2.T12*(vZ+wY)
+sys2.T13*wZ)*sys2.T3
3
pzd.el22 (sys2.T21*uX+0.5*sys2.
T22*(uY+vX)+0.5*sys2.
T23*(uZ+wX))*sys2.T21
+(0.5*sys2.T21*(uY+vX)
+sys2.T22*vY+0.5*sys2.
T23*(vZ+wY))*sys2.T22
+(0.5*sys2.T21*(uZ+wX
)+0.5*sys2.T22*(vZ+wY)
+sys2.T23*wZ)*sys2.T2
3
1 Strain tensor,
local coordinate
system, 22
component
Domain 4
pzd.el23 (sys2.T21*uX+0.5*sys2.
T22*(uY+vX)+0.5*sys2.
T23*(uZ+wX))*sys2.T31
+(0.5*sys2.T21*(uY+vX)
+sys2.T22*vY+0.5*sys2.
T23*(vZ+wY))*sys2.T32
+(0.5*sys2.T21*(uZ+wX
)+0.5*sys2.T22*(vZ+wY)
+sys2.T23*wZ)*sys2.T3
3
1 Strain tensor,
local coordinate
system, 23
component
Domain 4
pzd.el33 (sys2.T31*uX+0.5*sys2.
T32*(uY+vX)+0.5*sys2.
T33*(uZ+wX))*sys2.T31
+(0.5*sys2.T31*(uY+vX)
+sys2.T32*vY+0.5*sys2.
T33*(vZ+wY))*sys2.T32
+(0.5*sys2.T31*(uZ+wX
)+0.5*sys2.T32*(vZ+wY)
+sys2.T33*wZ)*sys2.T3
3
1 Strain tensor,
local coordinate
system, 33
component
Domain 4
pzd.eil11 0 1 Initial strain
tensor, local
coordinate
system, 11
component
Domain 4
pzd.eil12 0 1 Initial strain
tensor, local
coordinate
system, 12
component
Domain 4
pzd.eil13 0 1 Initial strain Domain 4
25

Name Expression Unit Description Selection
tensor, local
coordinate
system, 13
component
pzd.eil22 0 1 Initial strain
tensor, local
coordinate
system, 22
component
Domain 4
pzd.eil23 0 1 Initial strain
tensor, local
coordinate
system, 23
component
Domain 4
pzd.eil33 0 1 Initial strain
tensor, local
coordinate
system, 33
component
Domain 4
pzd.eiel11 0 1 Inelastic strain
tensor, local
coordinate
system, 11
component
Domain 4
pzd.eiel12 0 1 Inelastic strain
tensor, local
coordinate
system, 12
component
Domain 4
pzd.eiel13 0 1 Inelastic strain
tensor, local
coordinate
system, 13
component
Domain 4
pzd.eiel22 0 1 Inelastic strain
tensor, local
coordinate
system, 22
component
Domain 4
pzd.eiel23 0 1 Inelastic strain
tensor, local
coordinate
system, 23
Domain 4
26

Name Expression Unit Description Selection
component
pzd.eiel33 0 1 Inelastic strain
tensor, local
coordinate
system, 33
component
Domain 4
pzd.eel11 (sys2.T11*uX+0.5*sys2.
T12*(uY+vX)+0.5*sys2.
T13*(uZ+wX))*sys2.T11
+(0.5*sys2.T11*(uY+vX)
+sys2.T12*vY+0.5*sys2.
T13*(vZ+wY))*sys2.T12
+(0.5*sys2.T11*(uZ+wX
)+0.5*sys2.T12*(vZ+wY)
+sys2.T13*wZ)*sys2.T1
3-pzd.eil11-pzd.eiel11
1 Elastic strain
tensor, local
coordinate
system, 11
component
Domain 4
pzd.eel12 (sys2.T11*uX+0.5*sys2.
T12*(uY+vX)+0.5*sys2.
T13*(uZ+wX))*sys2.T21
+(0.5*sys2.T11*(uY+vX)
+sys2.T12*vY+0.5*sys2.
T13*(vZ+wY))*sys2.T22
+(0.5*sys2.T11*(uZ+wX
)+0.5*sys2.T12*(vZ+wY)
+sys2.T13*wZ)*sys2.T2
3-pzd.eil12-pzd.eiel12
1 Elastic strain
tensor, local
coordinate
system, 12
component
Domain 4
pzd.eel13 (sys2.T11*uX+0.5*sys2.
T12*(uY+vX)+0.5*sys2.
T13*(uZ+wX))*sys2.T31
+(0.5*sys2.T11*(uY+vX)
+sys2.T12*vY+0.5*sys2.
T13*(vZ+wY))*sys2.T32
+(0.5*sys2.T11*(uZ+wX
)+0.5*sys2.T12*(vZ+wY)
+sys2.T13*wZ)*sys2.T3
3-pzd.eil13-pzd.eiel13
1 Elastic strain
tensor, local
coordinate
system, 13
component
Domain 4
pzd.eel22 (sys2.T21*uX+0.5*sys2.
T22*(uY+vX)+0.5*sys2.
T23*(uZ+wX))*sys2.T21
+(0.5*sys2.T21*(uY+vX)
+sys2.T22*vY+0.5*sys2.
T23*(vZ+wY))*sys2.T22
+(0.5*sys2.T21*(uZ+wX
)+0.5*sys2.T22*(vZ+wY)
1 Elastic strain
tensor, local
coordinate
system, 22
component
Domain 4
27

Name Expression Unit Description Selection
+sys2.T23*wZ)*sys2.T2
3-pzd.eil22-pzd.eiel22
pzd.eel23 (sys2.T21*uX+0.5*sys2.
T22*(uY+vX)+0.5*sys2.
T23*(uZ+wX))*sys2.T31
+(0.5*sys2.T21*(uY+vX)
+sys2.T22*vY+0.5*sys2.
T23*(vZ+wY))*sys2.T32
+(0.5*sys2.T21*(uZ+wX
)+0.5*sys2.T22*(vZ+wY)
+sys2.T23*wZ)*sys2.T3
3-pzd.eil23-pzd.eiel23
1 Elastic strain
tensor, local
coordinate
system, 23
component
Domain 4
pzd.eel33 (sys2.T31*uX+0.5*sys2.
T32*(uY+vX)+0.5*sys2.
T33*(uZ+wX))*sys2.T31
+(0.5*sys2.T31*(uY+vX)
+sys2.T32*vY+0.5*sys2.
T33*(vZ+wY))*sys2.T32
+(0.5*sys2.T31*(uZ+wX
)+0.5*sys2.T32*(vZ+wY)
+sys2.T33*wZ)*sys2.T3
3-pzd.eil33-pzd.eiel33
1 Elastic strain
tensor, local
coordinate
system, 33
component
Domain 4
pzd.evol uX+vY+wZ 1 Volumetric strain Domain 4
pzd.J pzd.FdxX*pzd.FdyY*pzd
.FdzZ+pzd.FdxY*pzd.Fd
yZ*pzd.FdzX+pzd.FdxZ*
pzd.FdyX*pzd.FdzY-
pzd.FdxX*pzd.FdyZ*pzd
.FdzY-
pzd.FdxY*pzd.FdyX*pzd
.FdzZ-
pzd.FdxZ*pzd.FdyY*pzd
.FdzX
1 Volume ratio Domain 4
pzd.Sl11 pzd.cE11*pzd.eel11+2*
pzd.cE16*pzd.eel12+2*
pzd.cE15*pzd.eel13+pz
d.cE12*pzd.eel22+2*pz
d.cE14*pzd.eel23+pzd.c
E13*pzd.eel33-
pzd.eES11*(sys2.T11*p
zd.EX+sys2.T12*pzd.EY
+sys2.T13*pzd.EZ)-
pzd.eES21*(sys2.T21*p
zd.EX+sys2.T22*pzd.EY
N/m^2 Second Piola-
Kirchhoff stress,
local coordinate
system, 11
component
Domain 4
28

Name Expression Unit Description Selection
+sys2.T23*pzd.EZ)-
pzd.eES31*(sys2.T31*p
zd.EX+sys2.T32*pzd.EY
+sys2.T33*pzd.EZ)+pzd.
Sil11
pzd.Sl12 pzd.cE16*pzd.eel11+2*
pzd.cE66*pzd.eel12+2*
pzd.cE56*pzd.eel13+pz
d.cE26*pzd.eel22+2*pz
d.cE46*pzd.eel23+pzd.c
E36*pzd.eel33-
pzd.eES16*(sys2.T11*p
zd.EX+sys2.T12*pzd.EY
+sys2.T13*pzd.EZ)-
pzd.eES26*(sys2.T21*p
zd.EX+sys2.T22*pzd.EY
+sys2.T23*pzd.EZ)-
pzd.eES36*(sys2.T31*p
zd.EX+sys2.T32*pzd.EY
+sys2.T33*pzd.EZ)+pzd.
Sil12
N/m^2 Second Piola-
Kirchhoff stress,
local coordinate
system, 12
component
Domain 4
pzd.Sl13 pzd.cE15*pzd.eel11+2*
pzd.cE56*pzd.eel12+2*
pzd.cE55*pzd.eel13+pz
d.cE25*pzd.eel22+2*pz
d.cE45*pzd.eel23+pzd.c
E35*pzd.eel33-
pzd.eES15*(sys2.T11*p
zd.EX+sys2.T12*pzd.EY
+sys2.T13*pzd.EZ)-
pzd.eES25*(sys2.T21*p
zd.EX+sys2.T22*pzd.EY
+sys2.T23*pzd.EZ)-
pzd.eES35*(sys2.T31*p
zd.EX+sys2.T32*pzd.EY
+sys2.T33*pzd.EZ)+pzd.
Sil13
N/m^2 Second Piola-
Kirchhoff stress,
local coordinate
system, 13
component
Domain 4
pzd.Sl22 pzd.cE12*pzd.eel11+2*
pzd.cE26*pzd.eel12+2*
pzd.cE25*pzd.eel13+pz
d.cE22*pzd.eel22+2*pz
d.cE24*pzd.eel23+pzd.c
E23*pzd.eel33-
pzd.eES12*(sys2.T11*p
zd.EX+sys2.T12*pzd.EY
N/m^2 Second Piola-
Kirchhoff stress,
local coordinate
system, 22
component
Domain 4
29

Name Expression Unit Description Selection
+sys2.T13*pzd.EZ)-
pzd.eES22*(sys2.T21*p
zd.EX+sys2.T22*pzd.EY
+sys2.T23*pzd.EZ)-
pzd.eES32*(sys2.T31*p
zd.EX+sys2.T32*pzd.EY
+sys2.T33*pzd.EZ)+pzd.
Sil22
pzd.Sl23 pzd.cE14*pzd.eel11+2*
pzd.cE46*pzd.eel12+2*
pzd.cE45*pzd.eel13+pz
d.cE24*pzd.eel22+2*pz
d.cE44*pzd.eel23+pzd.c
E34*pzd.eel33-
pzd.eES14*(sys2.T11*p
zd.EX+sys2.T12*pzd.EY
+sys2.T13*pzd.EZ)-
pzd.eES24*(sys2.T21*p
zd.EX+sys2.T22*pzd.EY
+sys2.T23*pzd.EZ)-
pzd.eES34*(sys2.T31*p
zd.EX+sys2.T32*pzd.EY
+sys2.T33*pzd.EZ)+pzd.
Sil23
N/m^2 Second Piola-
Kirchhoff stress,
local coordinate
system, 23
component
Domain 4
pzd.Sl33 pzd.cE13*pzd.eel11+2*
pzd.cE36*pzd.eel12+2*
pzd.cE35*pzd.eel13+pz
d.cE23*pzd.eel22+2*pz
d.cE34*pzd.eel23+pzd.c
E33*pzd.eel33-
pzd.eES13*(sys2.T11*p
zd.EX+sys2.T12*pzd.EY
+sys2.T13*pzd.EZ)-
pzd.eES23*(sys2.T21*p
zd.EX+sys2.T22*pzd.EY
+sys2.T23*pzd.EZ)-
pzd.eES33*(sys2.T31*p
zd.EX+sys2.T32*pzd.EY
+sys2.T33*pzd.EZ)+pzd.
Sil33
N/m^2 Second Piola-
Kirchhoff stress,
local coordinate
system, 33
component
Domain 4
pzd.Sil11 0 N/m^2 Initial stress
tensor, local
coordinate
system, 11
component
Domain 4
30

Name Expression Unit Description Selection
pzd.Sil12 0 N/m^2 Initial stress
tensor, local
coordinate
system, 12
component
Domain 4
pzd.Sil13 0 N/m^2 Initial stress
tensor, local
coordinate
system, 13
component
Domain 4
pzd.Sil22 0 N/m^2 Initial stress
tensor, local
coordinate
system, 22
component
Domain 4
pzd.Sil23 0 N/m^2 Initial stress
tensor, local
coordinate
system, 23
component
Domain 4
pzd.Sil33 0 N/m^2 Initial stress
tensor, local
coordinate
system, 33
component
Domain 4
pzd.u_ttX d(root.comp1.uTIME,TI
ME)
m/s^2 Acceleration, X
component
Domain 4
pzd.u_ttY d(root.comp1.vTIME,TI
ME)
m/s^2 Acceleration, Y
component
Domain 4
pzd.u_ttZ d(root.comp1.wTIME,TI
ME)
m/s^2 Acceleration, Z
component
Domain 4
pzd.accX d(root.comp1.uTIME,TI
ME)
m/s^2 Total
acceleration, X
component
Domain 4
pzd.accY d(root.comp1.vTIME,TI
ME)
m/s^2 Total
acceleration, Y
component
Domain 4
pzd.accZ d(root.comp1.wTIME,TI
ME)
m/s^2 Total
acceleration, Z
component
Domain 4
pzd.u_tX root.comp1.uTIME m/s Velocity, X Domain 4
31

Name Expression Unit Description Selection
component
pzd.u_tY root.comp1.vTIME m/s Velocity, Y
component
Domain 4
pzd.u_tZ root.comp1.wTIME m/s Velocity, Z
component
Domain 4
pzd.disp sqrt(real(u)^2+real(v)^2
+real(w)^2)
m Total
displacement
Domain 4
pzd.SX (sys2.invT11*pzd.Sl11+
sys2.invT12*pzd.Sl12+s
ys2.invT13*pzd.Sl13)*(s
ys2.gSub11*sys2.invT1
1+sys2.gSub12*sys2.inv
T12+sys2.gSub13*sys2.i
nvT13)+(sys2.invT11*pz
d.Sl12+sys2.invT12*pzd
.Sl22+sys2.invT13*pzd.
Sl23)*(sys2.gSub21*sys
2.invT11+sys2.gSub22*
sys2.invT12+sys2.gSub2
3*sys2.invT13)+(sys2.in
vT11*pzd.Sl13+sys2.inv
T12*pzd.Sl23+sys2.invT
13*pzd.Sl33)*(sys2.gSu
b31*sys2.invT11+sys2.g
Sub32*sys2.invT12+sys
2.gSub33*sys2.invT13)
N/m^2 Second Piola-
Kirchhoff stress,
X component
Domain 4
pzd.SXY (sys2.invT11*pzd.Sl11+
sys2.invT12*pzd.Sl12+s
ys2.invT13*pzd.Sl13)*(s
ys2.gSub11*sys2.invT2
1+sys2.gSub12*sys2.inv
T22+sys2.gSub13*sys2.i
nvT23)+(sys2.invT11*pz
d.Sl12+sys2.invT12*pzd
.Sl22+sys2.invT13*pzd.
Sl23)*(sys2.gSub21*sys
2.invT21+sys2.gSub22*
sys2.invT22+sys2.gSub2
3*sys2.invT23)+(sys2.in
vT11*pzd.Sl13+sys2.inv
T12*pzd.Sl23+sys2.invT
13*pzd.Sl33)*(sys2.gSu
b31*sys2.invT21+sys2.g
Sub32*sys2.invT22+sys
N/m^2 Second Piola-
Kirchhoff stress,
XY component
Domain 4
32

Name Expression Unit Description Selection
2.gSub33*sys2.invT23)
pzd.SXZ (sys2.invT11*pzd.Sl11+
sys2.invT12*pzd.Sl12+s
ys2.invT13*pzd.Sl13)*(s
ys2.gSub11*sys2.invT3
1+sys2.gSub12*sys2.inv
T32+sys2.gSub13*sys2.i
nvT33)+(sys2.invT11*pz
d.Sl12+sys2.invT12*pzd
.Sl22+sys2.invT13*pzd.
Sl23)*(sys2.gSub21*sys
2.invT31+sys2.gSub22*
sys2.invT32+sys2.gSub2
3*sys2.invT33)+(sys2.in
vT11*pzd.Sl13+sys2.inv
T12*pzd.Sl23+sys2.invT
13*pzd.Sl33)*(sys2.gSu
b31*sys2.invT31+sys2.g
Sub32*sys2.invT32+sys
2.gSub33*sys2.invT33)
N/m^2 Second Piola-
Kirchhoff stress,
XZ component
Domain 4
pzd.SY (sys2.invT21*pzd.Sl11+
sys2.invT22*pzd.Sl12+s
ys2.invT23*pzd.Sl13)*(s
ys2.gSub11*sys2.invT2
1+sys2.gSub12*sys2.inv
T22+sys2.gSub13*sys2.i
nvT23)+(sys2.invT21*pz
d.Sl12+sys2.invT22*pzd
.Sl22+sys2.invT23*pzd.
Sl23)*(sys2.gSub21*sys
2.invT21+sys2.gSub22*
sys2.invT22+sys2.gSub2
3*sys2.invT23)+(sys2.in
vT21*pzd.Sl13+sys2.inv
T22*pzd.Sl23+sys2.invT
23*pzd.Sl33)*(sys2.gSu
b31*sys2.invT21+sys2.g
Sub32*sys2.invT22+sys
2.gSub33*sys2.invT23)
N/m^2 Second Piola-
Kirchhoff stress,
Y component
Domain 4
pzd.SYZ (sys2.invT21*pzd.Sl11+
sys2.invT22*pzd.Sl12+s
ys2.invT23*pzd.Sl13)*(s
ys2.gSub11*sys2.invT3
1+sys2.gSub12*sys2.inv
T32+sys2.gSub13*sys2.i
N/m^2 Second Piola-
Kirchhoff stress,
YZ component
Domain 4
33

Name Expression Unit Description Selection
nvT33)+(sys2.invT21*pz
d.Sl12+sys2.invT22*pzd
.Sl22+sys2.invT23*pzd.
Sl23)*(sys2.gSub21*sys
2.invT31+sys2.gSub22*
sys2.invT32+sys2.gSub2
3*sys2.invT33)+(sys2.in
vT21*pzd.Sl13+sys2.inv
T22*pzd.Sl23+sys2.invT
23*pzd.Sl33)*(sys2.gSu
b31*sys2.invT31+sys2.g
Sub32*sys2.invT32+sys
2.gSub33*sys2.invT33)
pzd.SZ (sys2.invT31*pzd.Sl11+
sys2.invT32*pzd.Sl12+s
ys2.invT33*pzd.Sl13)*(s
ys2.gSub11*sys2.invT3
1+sys2.gSub12*sys2.inv
T32+sys2.gSub13*sys2.i
nvT33)+(sys2.invT31*pz
d.Sl12+sys2.invT32*pzd
.Sl22+sys2.invT33*pzd.
Sl23)*(sys2.gSub21*sys
2.invT31+sys2.gSub22*
sys2.invT32+sys2.gSub2
3*sys2.invT33)+(sys2.in
vT31*pzd.Sl13+sys2.inv
T32*pzd.Sl23+sys2.invT
33*pzd.Sl33)*(sys2.gSu
b31*sys2.invT31+sys2.g
Sub32*sys2.invT32+sys
2.gSub33*sys2.invT33)
N/m^2 Second Piola-
Kirchhoff stress,
Z component
Domain 4
pzd.PxX pzd.FdxX*pzd.SX+pzd.F
dxY*pzd.SXY+pzd.FdxZ*
pzd.SXZ
N/m^2 First Piola-
Kirchhoff stress,
xX component
Domain 4
pzd.PyX pzd.FdyX*pzd.SX+pzd.F
dyY*pzd.SXY+pzd.FdyZ*
pzd.SXZ
N/m^2 First Piola-
Kirchhoff stress,
yX component
Domain 4
pzd.PzX pzd.FdzX*pzd.SX+pzd.F
dzY*pzd.SXY+pzd.FdzZ*
pzd.SXZ
N/m^2 First Piola-
Kirchhoff stress,
zX component
Domain 4
pzd.PxY pzd.FdxX*pzd.SXY+pzd.
FdxY*pzd.SY+pzd.FdxZ*
pzd.SYZ
N/m^2 First Piola-
Kirchhoff stress,
xY component
Domain 4
34

Name Expression Unit Description Selection
pzd.PyY pzd.FdyX*pzd.SXY+pzd.
FdyY*pzd.SY+pzd.FdyZ*
pzd.SYZ
N/m^2 First Piola-
Kirchhoff stress,
yY component
Domain 4
pzd.PzY pzd.FdzX*pzd.SXY+pzd.
FdzY*pzd.SY+pzd.FdzZ*
pzd.SYZ
N/m^2 First Piola-
Kirchhoff stress,
zY component
Domain 4
pzd.PxZ pzd.FdxX*pzd.SXZ+pzd.
FdxY*pzd.SYZ+pzd.FdxZ
*pzd.SZ
N/m^2 First Piola-
Kirchhoff stress,
xZ component
Domain 4
pzd.PyZ pzd.FdyX*pzd.SXZ+pzd.
FdyY*pzd.SYZ+pzd.FdyZ
*pzd.SZ
N/m^2 First Piola-
Kirchhoff stress,
yZ component
Domain 4
pzd.PzZ pzd.FdzX*pzd.SXZ+pzd.
FdzY*pzd.SYZ+pzd.FdzZ
*pzd.SZ
N/m^2 First Piola-
Kirchhoff stress,
zZ component
Domain 4
pzd.sx pzd.SX N/m^2 Stress tensor, x
component
Domain 4
pzd.sxy pzd.SXY N/m^2 Stress tensor, xy
component
Domain 4
pzd.sxz pzd.SXZ N/m^2 Stress tensor, xz
component
Domain 4
pzd.sy pzd.SY N/m^2 Stress tensor, y
component
Domain 4
pzd.syz pzd.SYZ N/m^2 Stress tensor, yz
component
Domain 4
pzd.sz pzd.SZ N/m^2 Stress tensor, z
component
Domain 4
pzd.sl11 pzd.Sl11 N/m^2 Stress tensor,
local coordinate
system, 11
component
Domain 4
pzd.sl12 pzd.Sl12 N/m^2 Stress tensor,
local coordinate
system, 12
component
Domain 4
pzd.sl13 pzd.Sl13 N/m^2 Stress tensor,
local coordinate
system, 13
component
Domain 4
pzd.sl22 pzd.Sl22 N/m^2 Stress tensor, Domain 4
35

Name Expression Unit Description Selection
local coordinate
system, 22
component
pzd.sl23 pzd.Sl23 N/m^2 Stress tensor,
local coordinate
system, 23
component
Domain 4
pzd.sl33 pzd.Sl33 N/m^2 Stress tensor,
local coordinate
system, 33
component
Domain 4
pzd.sdevx pzd.sx-
(pzd.sx+pzd.sy+pzd.sz)/
3
N/m^2 Deviatoric stress
tensor, x
component
Domain 4
pzd.sdevxy pzd.sxy N/m^2 Deviatoric stress
tensor, xy
component
Domain 4
pzd.sdevxz pzd.sxz N/m^2 Deviatoric stress
tensor, xz
component
Domain 4
pzd.sdevy pzd.sy-
(pzd.sx+pzd.sy+pzd.sz)/
3
N/m^2 Deviatoric stress
tensor, y
component
Domain 4
pzd.sdevyz pzd.syz N/m^2 Deviatoric stress
tensor, yz
component
Domain 4
pzd.sdevz pzd.sz-
(pzd.sx+pzd.sy+pzd.sz)/
3
N/m^2 Deviatoric stress
tensor, z
component
Domain 4
pzd.SdevX pzd.SX-
(pzd.SX+pzd.SY+pzd.SZ)
/3
N/m^2 Deviatoric
second Piola-
Kirchhoff stress,
X component
Domain 4
pzd.SdevXY pzd.SXY N/m^2 Deviatoric
second Piola-
Kirchhoff stress,
XY component
Domain 4
pzd.SdevXZ pzd.SXZ N/m^2 Deviatoric
second Piola-
Kirchhoff stress,
XZ component
Domain 4
36

Name Expression Unit Description Selection
pzd.SdevY pzd.SY-
(pzd.SX+pzd.SY+pzd.SZ)
/3
N/m^2 Deviatoric
second Piola-
Kirchhoff stress,
Y component
Domain 4
pzd.SdevYZ pzd.SYZ N/m^2 Deviatoric
second Piola-
Kirchhoff stress,
YZ component
Domain 4
pzd.SdevZ pzd.SZ-
(pzd.SX+pzd.SY+pzd.SZ)
/3
N/m^2 Deviatoric
second Piola-
Kirchhoff stress,
Z component
Domain 4
pzd.I1s pzd.sl11+pzd.sl22+pzd.s
l33
N/m^2 First principal
invariant of
stress
Domain 4
pzd.I2s 0.5*(pzd.I1s^2-
pzd.sl11^2-
2*pzd.sl12^2-
2*pzd.sl13^2-
pzd.sl22^2-
2*pzd.sl23^2-
pzd.sl33^2)
kg^2/(m^2*s^4) Second principal
invariant of
stress
Domain 4
pzd.I3s pzd.sl11*pzd.sl22*pzd.s
l33+2*pzd.sl12*pzd.sl2
3*pzd.sl13-
pzd.sl11*pzd.sl23^2-
pzd.sl12^2*pzd.sl33-
pzd.sl22*pzd.sl13^2
kg^3/(m^3*s^6) Third principal
invariant of
stress
Domain 4
pzd.II2s 0.5*((pzd.sl11-
(pzd.sl11+pzd.sl22+pzd.
sl33)/3)^2+2*pzd.sl12^
2+2*pzd.sl13^2+(pzd.sl
22-
(pzd.sl11+pzd.sl22+pzd.
sl33)/3)^2+2*pzd.sl23^
2+(pzd.sl33-
(pzd.sl11+pzd.sl22+pzd.
sl33)/3)^2)
kg^2/(m^2*s^4) Second invariant
of stress deviator
Domain 4
pzd.II3s (pzd.sl11-
(pzd.sl11+pzd.sl22+pzd.
sl33)/3)*(pzd.sl22-
(pzd.sl11+pzd.sl22+pzd.
sl33)/3)*(pzd.sl33-
kg^3/(m^3*s^6) Third invariant of
stress deviator
Domain 4
37

Name Expression Unit Description Selection
(pzd.sl11+pzd.sl22+pzd.
sl33)/3)+2*pzd.sl12*pz
d.sl23*pzd.sl13-
(pzd.sl11-
(pzd.sl11+pzd.sl22+pzd.
sl33)/3)*pzd.sl23^2-
pzd.sl12^2*(pzd.sl33-
(pzd.sl11+pzd.sl22+pzd.
sl33)/3)-(pzd.sl22-
(pzd.sl11+pzd.sl22+pzd.
sl33)/3)*pzd.sl13^2
pzd.thetaL atan2(sqrt(max(0.1481
4814814814814*pzd.II
2s^3-
pzd.II3s^2,eps)),pzd.II3s
)/3
rad Lode angle Domain 4
pzd.eelvol pzd.eel11+pzd.eel22+p
zd.eel33
1 First invariant of
elastic strain
tensor
Domain 4
pzd.eeldev11 pzd.eel11-
(pzd.eel11+pzd.eel22+p
zd.eel33)/3
1 Deviatoric elastic
strain tensor,
local coordinate
system, 11
component
Domain 4
pzd.eeldev12 pzd.eel12 1 Deviatoric elastic
strain tensor,
local coordinate
system, 12
component
Domain 4
pzd.eeldev13 pzd.eel13 1 Deviatoric elastic
strain tensor,
local coordinate
system, 13
component
Domain 4
pzd.eeldev22 pzd.eel22-
(pzd.eel11+pzd.eel22+p
zd.eel33)/3
1 Deviatoric elastic
strain tensor,
local coordinate
system, 22
component
Domain 4
pzd.eeldev23 pzd.eel23 1 Deviatoric elastic
strain tensor,
local coordinate
system, 23
Domain 4
38

Name Expression Unit Description Selection
component
pzd.eeldev33 pzd.eel33-
(pzd.eel11+pzd.eel22+p
zd.eel33)/3
1 Deviatoric elastic
strain tensor,
local coordinate
system, 33
component
Domain 4
pzd.II2eel 0.5*((pzd.eel11-
(pzd.eel11+pzd.eel22+p
zd.eel33)/3)^2+2*pzd.e
el12^2+2*pzd.eel13^2+
(pzd.eel22-
(pzd.eel11+pzd.eel22+p
zd.eel33)/3)^2+2*pzd.e
el23^2+(pzd.eel33-
(pzd.eel11+pzd.eel22+p
zd.eel33)/3)^2)
1 Second invariant
of deviatoric
elastic strain
tensor
Domain 4
pzd.pm -
(pzd.sl11+pzd.sl22+pzd.
sl33)/3
N/m^2 Pressure Domain 4
pzd.curlUX wY-vZ 1 Curl of
displacement, X
component
Domain 4
pzd.curlUY -wX+uZ 1 Curl of
displacement, Y
component
Domain 4
pzd.curlUZ vX-uY 1 Curl of
displacement, Z
component
Domain 4
pzd.mises (3*pzd.II2s)^0.5 N/m^2 von Mises stress Domain 4
pzd.tresca pzd.sp1-pzd.sp3 N/m^2 Tresca stress Domain 4
pzd.RFx reacf(u) N Reaction force, x
component
Domain 4
pzd.RFy reacf(v) N Reaction force, y
component
Domain 4
pzd.RFz reacf(w) N Reaction force, z
component
Domain 4
pzd.RMx pzd.RFz*(y-
pzd.refpnty)-
pzd.RFy*(z-pzd.refpntz)
N*m Reaction
moment, x
component
Domain 4
pzd.RMy -pzd.RFz*(x-
pzd.refpntx)+pzd.RFx*(
N*m Reaction
moment, y
Domain 4
39

Name Expression Unit Description Selection
z-pzd.refpntz) component
pzd.RMz pzd.RFy*(x-
pzd.refpntx)-
pzd.RFx*(y-pzd.refpnty)
N*m Reaction
moment, z
component
Domain 4
pzd.Tax pzd.sx*pzd.nx+pzd.sxy*
pzd.ny+pzd.sxz*pzd.nz
N/m^2 Traction
(force/area), x
component
Boundaries
1419
pzd.Tay pzd.sxy*pzd.nx+pzd.sy*
pzd.ny+pzd.syz*pzd.nz
N/m^2 Traction
(force/area), y
component
Boundaries
1419
pzd.Taz pzd.sxz*pzd.nx+pzd.syz
*pzd.ny+pzd.sz*pzd.nz
N/m^2 Traction
(force/area), z
component
Boundaries
1419
pzd.Ws 0.5*((pzd.Sl11+pzd.Sil1
1)*pzd.eel11+2*(pzd.Sl
12+pzd.Sil12)*pzd.eel1
2+2*(pzd.Sl13+pzd.Sil1
3)*pzd.eel13+(pzd.Sl22
+pzd.Sil22)*pzd.eel22+
2*(pzd.Sl23+pzd.Sil23)*
pzd.eel23+(pzd.Sl33+pz
d.Sil33)*pzd.eel33)
J/m^3 Elastic strain
energy density
Domain 4
pzd.Ws_tot root.comp1.pzd.pzm1.i
nt8(pzd.Ws)
J Total elastic
strain energy
Global
VXt root.comp1.VXTIME m*kg/(s^4*A) Gradient of V,
first time
derivative, X
component
Domain 4
VYt root.comp1.VYTIME m*kg/(s^4*A) Gradient of V,
first time
derivative, Y
component
Domain 4
VZt root.comp1.VZTIME m*kg/(s^4*A) Gradient of V,
first time
derivative, Z
component
Domain 4
VXtt root.comp1.VXTIMETI
ME
m*kg/(s^5*A) Gradient of V,
second time
derivative, X
component
Domain 4
VYtt root.comp1.VYTIMETIM m*kg/(s^5*A) Gradient of V,
second time
Domain 4
40

Name Expression Unit Description Selection
E derivative, Y
component
VZtt root.comp1.VZTIMETIM
E
m*kg/(s^5*A) Gradient of V,
second time
derivative, Z
component
Domain 4
Vxt root.comp1.VxTIME-
Vxx*d(x,TIME)-
Vxy*d(y,TIME)-
Vxz*d(z,TIME)
m*kg/(s^4*A) Gradient of V,
first time
derivative, x
component
Domain 4
Vyt root.comp1.VyTIME-
Vyx*d(x,TIME)-
Vyy*d(y,TIME)-
Vyz*d(z,TIME)
m*kg/(s^4*A) Gradient of V,
first time
derivative, y
component
Domain 4
Vzt root.comp1.VzTIME-
Vzx*d(x,TIME)-
Vzy*d(y,TIME)-
Vzz*d(z,TIME)
m*kg/(s^4*A) Gradient of V,
first time
derivative, z
component
Domain 4
Vxtt d(root.comp1.VxTIME-
Vxx*d(x,TIME)-
Vxy*d(y,TIME)-
Vxz*d(z,TIME),TIME)-
d(root.comp1.VxTIME-
Vxx*d(x,TIME)-
Vxy*d(y,TIME)-
Vxz*d(z,TIME),x)*d(x,TI
ME)-
d(root.comp1.VxTIME-
Vxx*d(x,TIME)-
Vxy*d(y,TIME)-
Vxz*d(z,TIME),y)*d(y,TI
ME)-
d(root.comp1.VxTIME-
Vxx*d(x,TIME)-
Vxy*d(y,TIME)-
Vxz*d(z,TIME),z)*d(z,TI
ME)
m*kg/(s^5*A) Gradient of V,
second time
derivative, x
component
Domain 4
Vytt d(root.comp1.VyTIME-
Vyx*d(x,TIME)-
Vyy*d(y,TIME)-
Vyz*d(z,TIME),TIME)-
d(root.comp1.VyTIME-
Vyx*d(x,TIME)-
Vyy*d(y,TIME)-
m*kg/(s^5*A) Gradient of V,
second time
derivative, y
component
Domain 4
41

Name Expression Unit Description Selection
Vyz*d(z,TIME),x)*d(x,TI
ME)-
d(root.comp1.VyTIME-
Vyx*d(x,TIME)-
Vyy*d(y,TIME)-
Vyz*d(z,TIME),y)*d(y,TI
ME)-
d(root.comp1.VyTIME-
Vyx*d(x,TIME)-
Vyy*d(y,TIME)-
Vyz*d(z,TIME),z)*d(z,TI
ME)
Vztt d(root.comp1.VzTIME-
Vzx*d(x,TIME)-
Vzy*d(y,TIME)-
Vzz*d(z,TIME),TIME)-
d(root.comp1.VzTIME-
Vzx*d(x,TIME)-
Vzy*d(y,TIME)-
Vzz*d(z,TIME),x)*d(x,TI
ME)-
d(root.comp1.VzTIME-
Vzx*d(x,TIME)-
Vzy*d(y,TIME)-
Vzz*d(z,TIME),y)*d(y,TI
ME)-
d(root.comp1.VzTIME-
Vzx*d(x,TIME)-
Vzy*d(y,TIME)-
Vzz*d(z,TIME),z)*d(z,TI
ME)
m*kg/(s^5*A) Gradient of V,
second time
derivative, z
component
Domain 4
pzd.cE11 material.cE11 Pa Elasticity matrix,
11 component
Domain 4
pzd.cE21 material.cE21 Pa Elasticity matrix,
21 component
Domain 4
pzd.cE31 material.cE31 Pa Elasticity matrix,
31 component
Domain 4
pzd.cE41 material.cE41 Pa Elasticity matrix,
41 component
Domain 4
pzd.cE51 material.cE51 Pa Elasticity matrix,
51 component
Domain 4
pzd.cE61 material.cE61 Pa Elasticity matrix,
61 component
Domain 4
42

Name Expression Unit Description Selection
pzd.cE12 material.cE12 Pa Elasticity matrix,
12 component
Domain 4
pzd.cE22 material.cE22 Pa Elasticity matrix,
22 component
Domain 4
pzd.cE32 material.cE32 Pa Elasticity matrix,
32 component
Domain 4
pzd.cE42 material.cE42 Pa Elasticity matrix,
42 component
Domain 4
pzd.cE52 material.cE52 Pa Elasticity matrix,
52 component
Domain 4
pzd.cE62 material.cE62 Pa Elasticity matrix,
62 component
Domain 4
pzd.cE13 material.cE13 Pa Elasticity matrix,
13 component
Domain 4
pzd.cE23 material.cE23 Pa Elasticity matrix,
23 component
Domain 4
pzd.cE33 material.cE33 Pa Elasticity matrix,
33 component
Domain 4
pzd.cE43 material.cE43 Pa Elasticity matrix,
43 component
Domain 4
pzd.cE53 material.cE53 Pa Elasticity matrix,
53 component
Domain 4
pzd.cE63 material.cE63 Pa Elasticity matrix,
63 component
Domain 4
pzd.cE14 material.cE14 Pa Elasticity matrix,
14 component
Domain 4
pzd.cE24 material.cE24 Pa Elasticity matrix,
24 component
Domain 4
pzd.cE34 material.cE34 Pa Elasticity matrix,
34 component
Domain 4
pzd.cE44 material.cE44 Pa Elasticity matrix,
44 component
Domain 4
pzd.cE54 material.cE54 Pa Elasticity matrix,
54 component
Domain 4
pzd.cE64 material.cE64 Pa Elasticity matrix,
64 component
Domain 4
pzd.cE15 material.cE15 Pa Elasticity matrix,
15 component
Domain 4
43

Name Expression Unit Description Selection
pzd.cE25 material.cE25 Pa Elasticity matrix,
25 component
Domain 4
pzd.cE35 material.cE35 Pa Elasticity matrix,
35 component
Domain 4
pzd.cE45 material.cE45 Pa Elasticity matrix,
45 component
Domain 4
pzd.cE55 material.cE55 Pa Elasticity matrix,
55 component
Domain 4
pzd.cE65 material.cE65 Pa Elasticity matrix,
65 component
Domain 4
pzd.cE16 material.cE16 Pa Elasticity matrix,
16 component
Domain 4
pzd.cE26 material.cE26 Pa Elasticity matrix,
26 component
Domain 4
pzd.cE36 material.cE36 Pa Elasticity matrix,
36 component
Domain 4
pzd.cE46 material.cE46 Pa Elasticity matrix,
46 component
Domain 4
pzd.cE56 material.cE56 Pa Elasticity matrix,
56 component
Domain 4
pzd.cE66 material.cE66 Pa Elasticity matrix,
66 component
Domain 4
pzd.eES11 material.eES11 C/m^2 Coupling matrix,
11 component
Domain 4
pzd.eES21 material.eES21 C/m^2 Coupling matrix,
21 component
Domain 4
pzd.eES31 material.eES31 C/m^2 Coupling matrix,
31 component
Domain 4
pzd.eES12 material.eES12 C/m^2 Coupling matrix,
12 component
Domain 4
pzd.eES22 material.eES22 C/m^2 Coupling matrix,
22 component
Domain 4
pzd.eES32 material.eES32 C/m^2 Coupling matrix,
32 component
Domain 4
pzd.eES13 material.eES13 C/m^2 Coupling matrix,
13 component
Domain 4
pzd.eES23 material.eES23 C/m^2 Coupling matrix,
23 component
Domain 4
44

Name Expression Unit Description Selection
pzd.eES33 material.eES33 C/m^2 Coupling matrix,
33 component
Domain 4
pzd.eES14 material.eES14 C/m^2 Coupling matrix,
14 component
Domain 4
pzd.eES24 material.eES24 C/m^2 Coupling matrix,
24 component
Domain 4
pzd.eES34 material.eES34 C/m^2 Coupling matrix,
34 component
Domain 4
pzd.eES15 material.eES15 C/m^2 Coupling matrix,
15 component
Domain 4
pzd.eES25 material.eES25 C/m^2 Coupling matrix,
25 component
Domain 4
pzd.eES35 material.eES35 C/m^2 Coupling matrix,
35 component
Domain 4
pzd.eES16 material.eES16 C/m^2 Coupling matrix,
16 component
Domain 4
pzd.eES26 material.eES26 C/m^2 Coupling matrix,
26 component
Domain 4
pzd.eES36 material.eES36 C/m^2 Coupling matrix,
36 component
Domain 4
pzd.epsilonrS11 material.epsilonrS11 1 Relative
permittivity, 11
component
Domain 4
pzd.epsilonrS21 material.epsilonrS21 1 Relative
permittivity, 21
component
Domain 4
pzd.epsilonrS31 material.epsilonrS31 1 Relative
permittivity, 31
component
Domain 4
pzd.epsilonrS12 material.epsilonrS12 1 Relative
permittivity, 12
component
Domain 4
pzd.epsilonrS22 material.epsilonrS22 1 Relative
permittivity, 22
component
Domain 4
pzd.epsilonrS32 material.epsilonrS32 1 Relative
permittivity, 32
component
Domain 4
pzd.epsilonrS13 material.epsilonrS13 1 Relative Domain 4
45

Name Expression Unit Description Selection
permittivity, 13
component
pzd.epsilonrS23 material.epsilonrS23 1 Relative
permittivity, 23
component
Domain 4
pzd.epsilonrS33 material.epsilonrS33 1 Relative
permittivity, 33
component
Domain 4
pzd.epsilonS11 epsilon0_const*pzd.eps
ilonrS11
F/m Electrical
permittivity, 11
component
Domain 4
pzd.epsilonS12 epsilon0_const*pzd.eps
ilonrS12
F/m Electrical
permittivity, 12
component
Domain 4
pzd.epsilonS13 epsilon0_const*pzd.eps
ilonrS13
F/m Electrical
permittivity, 13
component
Domain 4
pzd.epsilonS22 epsilon0_const*pzd.eps
ilonrS22
F/m Electrical
permittivity, 22
component
Domain 4
pzd.epsilonS23 epsilon0_const*pzd.eps
ilonrS23
F/m Electrical
permittivity, 23
component
Domain 4
pzd.epsilonS33 epsilon0_const*pzd.eps
ilonrS33
F/m Electrical
permittivity, 33
component
Domain 4
pzd.TXX sys2.T11 1 Coordinate
system in
piezoelectric
domain, XX
component
Domain 4
pzd.TYX sys2.T21 1 Coordinate
system in
piezoelectric
domain, YX
component
Domain 4
pzd.TZX sys2.T31 1 Coordinate
system in
piezoelectric
domain, ZX
component
Domain 4
46

Name Expression Unit Description Selection
pzd.TXY sys2.T12 1 Coordinate
system in
piezoelectric
domain, XY
component
Domain 4
pzd.TYY sys2.T22 1 Coordinate
system in
piezoelectric
domain, YY
component
Domain 4
pzd.TZY sys2.T32 1 Coordinate
system in
piezoelectric
domain, ZY
component
Domain 4
pzd.TXZ sys2.T13 1 Coordinate
system in
piezoelectric
domain, XZ
component
Domain 4
pzd.TYZ sys2.T23 1 Coordinate
system in
piezoelectric
domain, YZ
component
Domain 4
pzd.TZZ sys2.T33 1 Coordinate
system in
piezoelectric
domain, ZZ
component
Domain 4
pzd.DrX 0 C/m^2 Remanent
electric
displacement
field, X
component
Domain 4
pzd.DrY 0 C/m^2 Remanent
electric
displacement
field, Y
component
Domain 4
pzd.DrZ 0 C/m^2 Remanent
electric
displacement
Domain 4
47

Name Expression Unit Description Selection
field, Z
component
pzd.PX sys2.invT11*(pzd.eES11
*pzd.eel11+2*pzd.eES1
6*pzd.eel12+2*pzd.eES
15*pzd.eel13+pzd.eES1
2*pzd.eel22+2*pzd.eES
14*pzd.eel23+pzd.eES1
3*pzd.eel33+(pzd.epsil
onS11-
epsilon0_const)*(sys2.T
11*pzd.EX+sys2.T12*pz
d.EY+sys2.T13*pzd.EZ)+
pzd.epsilonS12*(sys2.T
21*pzd.EX+sys2.T22*pz
d.EY+sys2.T23*pzd.EZ)+
pzd.epsilonS13*(sys2.T
31*pzd.EX+sys2.T32*pz
d.EY+sys2.T33*pzd.EZ))
+sys2.invT12*(pzd.eES2
1*pzd.eel11+2*pzd.eES
26*pzd.eel12+2*pzd.eE
S25*pzd.eel13+pzd.eES
22*pzd.eel22+2*pzd.eE
S24*pzd.eel23+pzd.eES
23*pzd.eel33+pzd.epsil
onS12*(sys2.T11*pzd.E
X+sys2.T12*pzd.EY+sys
2.T13*pzd.EZ)+(pzd.eps
ilonS22-
epsilon0_const)*(sys2.T
21*pzd.EX+sys2.T22*pz
d.EY+sys2.T23*pzd.EZ)+
pzd.epsilonS23*(sys2.T
31*pzd.EX+sys2.T32*pz
d.EY+sys2.T33*pzd.EZ))
+sys2.invT13*(pzd.eES3
1*pzd.eel11+2*pzd.eES
36*pzd.eel12+2*pzd.eE
S35*pzd.eel13+pzd.eES
32*pzd.eel22+2*pzd.eE
S34*pzd.eel23+pzd.eES
33*pzd.eel33+pzd.epsil
onS13*(sys2.T11*pzd.E
X+sys2.T12*pzd.EY+sys
2.T13*pzd.EZ)+pzd.epsil
C/m^2 Polarization, X
component
Domain 4
48

Name Expression Unit Description Selection
onS23*(sys2.T21*pzd.E
X+sys2.T22*pzd.EY+sys
2.T23*pzd.EZ)+(pzd.eps
ilonS33-
epsilon0_const)*(sys2.T
31*pzd.EX+sys2.T32*pz
d.EY+sys2.T33*pzd.EZ))
+pzd.DrX
pzd.PY sys2.invT21*(pzd.eES11
*pzd.eel11+2*pzd.eES1
6*pzd.eel12+2*pzd.eES
15*pzd.eel13+pzd.eES1
2*pzd.eel22+2*pzd.eES
14*pzd.eel23+pzd.eES1
3*pzd.eel33+(pzd.epsil
onS11-
epsilon0_const)*(sys2.T
11*pzd.EX+sys2.T12*pz
d.EY+sys2.T13*pzd.EZ)+
pzd.epsilonS12*(sys2.T
21*pzd.EX+sys2.T22*pz
d.EY+sys2.T23*pzd.EZ)+
pzd.epsilonS13*(sys2.T
31*pzd.EX+sys2.T32*pz
d.EY+sys2.T33*pzd.EZ))
+sys2.invT22*(pzd.eES2
1*pzd.eel11+2*pzd.eES
26*pzd.eel12+2*pzd.eE
S25*pzd.eel13+pzd.eES
22*pzd.eel22+2*pzd.eE
S24*pzd.eel23+pzd.eES
23*pzd.eel33+pzd.epsil
onS12*(sys2.T11*pzd.E
X+sys2.T12*pzd.EY+sys
2.T13*pzd.EZ)+(pzd.eps
ilonS22-
epsilon0_const)*(sys2.T
21*pzd.EX+sys2.T22*pz
d.EY+sys2.T23*pzd.EZ)+
pzd.epsilonS23*(sys2.T
31*pzd.EX+sys2.T32*pz
d.EY+sys2.T33*pzd.EZ))
+sys2.invT23*(pzd.eES3
1*pzd.eel11+2*pzd.eES
36*pzd.eel12+2*pzd.eE
S35*pzd.eel13+pzd.eES
C/m^2 Polarization, Y
component
Domain 4
49

Name Expression Unit Description Selection
32*pzd.eel22+2*pzd.eE
S34*pzd.eel23+pzd.eES
33*pzd.eel33+pzd.epsil
onS13*(sys2.T11*pzd.E
X+sys2.T12*pzd.EY+sys
2.T13*pzd.EZ)+pzd.epsil
onS23*(sys2.T21*pzd.E
X+sys2.T22*pzd.EY+sys
2.T23*pzd.EZ)+(pzd.eps
ilonS33-
epsilon0_const)*(sys2.T
31*pzd.EX+sys2.T32*pz
d.EY+sys2.T33*pzd.EZ))
+pzd.DrY
pzd.PZ sys2.invT31*(pzd.eES11
*pzd.eel11+2*pzd.eES1
6*pzd.eel12+2*pzd.eES
15*pzd.eel13+pzd.eES1
2*pzd.eel22+2*pzd.eES
14*pzd.eel23+pzd.eES1
3*pzd.eel33+(pzd.epsil
onS11-
epsilon0_const)*(sys2.T
11*pzd.EX+sys2.T12*pz
d.EY+sys2.T13*pzd.EZ)+
pzd.epsilonS12*(sys2.T
21*pzd.EX+sys2.T22*pz
d.EY+sys2.T23*pzd.EZ)+
pzd.epsilonS13*(sys2.T
31*pzd.EX+sys2.T32*pz
d.EY+sys2.T33*pzd.EZ))
+sys2.invT32*(pzd.eES2
1*pzd.eel11+2*pzd.eES
26*pzd.eel12+2*pzd.eE
S25*pzd.eel13+pzd.eES
22*pzd.eel22+2*pzd.eE
S24*pzd.eel23+pzd.eES
23*pzd.eel33+pzd.epsil
onS12*(sys2.T11*pzd.E
X+sys2.T12*pzd.EY+sys
2.T13*pzd.EZ)+(pzd.eps
ilonS22-
epsilon0_const)*(sys2.T
21*pzd.EX+sys2.T22*pz
d.EY+sys2.T23*pzd.EZ)+
pzd.epsilonS23*(sys2.T
C/m^2 Polarization, Z
component
Domain 4
50

Name Expression Unit Description Selection
31*pzd.EX+sys2.T32*pz
d.EY+sys2.T33*pzd.EZ))
+sys2.invT33*(pzd.eES3
1*pzd.eel11+2*pzd.eES
36*pzd.eel12+2*pzd.eE
S35*pzd.eel13+pzd.eES
32*pzd.eel22+2*pzd.eE
S34*pzd.eel23+pzd.eES
33*pzd.eel33+pzd.epsil
onS13*(sys2.T11*pzd.E
X+sys2.T12*pzd.EY+sys
2.T13*pzd.EZ)+pzd.epsil
onS23*(sys2.T21*pzd.E
X+sys2.T22*pzd.EY+sys
2.T23*pzd.EZ)+(pzd.eps
ilonS33-
epsilon0_const)*(sys2.T
31*pzd.EX+sys2.T32*pz
d.EY+sys2.T33*pzd.EZ))
+pzd.DrZ
pzd.Dx epsilon0_const*pzd.Ex+
(spatial.invF11*pzd.PX+
spatial.invF12*pzd.PY+s
patial.invF13*pzd.PZ)*s
patial.detInvF
C/m^2 Electric
displacement
field, x
component
Domain 4
pzd.Dy epsilon0_const*pzd.Ey+
(spatial.invF21*pzd.PX+
spatial.invF22*pzd.PY+s
patial.invF23*pzd.PZ)*s
patial.detInvF
C/m^2 Electric
displacement
field, y
component
Domain 4
pzd.Dz epsilon0_const*pzd.Ez+
(spatial.invF31*pzd.PX+
spatial.invF32*pzd.PY+s
patial.invF33*pzd.PZ)*s
patial.detInvF
C/m^2 Electric
displacement
field, z
component
Domain 4
pzd.Ex -Vx V/m Electric field, x
component
Domain 4
pzd.Ey -Vy V/m Electric field, y
component
Domain 4
pzd.Ez -Vz V/m Electric field, z
component
Domain 4
pzd.normE sqrt(realdot(pzd.Ex,pzd.
Ex)+realdot(pzd.Ey,pzd.
V/m Electric field
norm
Domain 4
51

Name Expression Unit Description Selection
Ey)+realdot(pzd.Ez,pzd.
Ez))
pzd.normD sqrt(realdot(pzd.Dx,pzd
.Dx)+realdot(pzd.Dy,pz
d.Dy)+realdot(pzd.Dz,pz
d.Dz))
C/m^2 Electric
displacement
field norm
Domain 4
pzd.normJ sqrt(sys2.gSub11*reald
ot(root.comp1.pzd.Jx1,
sys2.gSub11*root.comp
1.pzd.Jx1+sys2.gSub12*
root.comp1.pzd.Jx2+sys
2.gSub13*root.comp1.
pzd.Jx3)+sys2.gSub12*r
ealdot(root.comp1.pzd.
Jx2,sys2.gSub21*root.c
omp1.pzd.Jx1+sys2.gSu
b22*root.comp1.pzd.Jx
2+sys2.gSub23*root.co
mp1.pzd.Jx3)+sys2.gSu
b13*realdot(root.comp
1.pzd.Jx3,sys2.gSub31*
root.comp1.pzd.Jx1+sys
2.gSub32*root.comp1.
pzd.Jx2+sys2.gSub33*r
oot.comp1.pzd.Jx3)+sys
2.gSub21*realdot(root.
comp1.pzd.Jx1,sys2.gSu
b11*root.comp1.pzd.Jx
1+sys2.gSub12*root.co
mp1.pzd.Jx2+sys2.gSub
13*root.comp1.pzd.Jx3
)+sys2.gSub22*realdot(
root.comp1.pzd.Jx2,sys
2.gSub21*root.comp1.
pzd.Jx1+sys2.gSub22*r
oot.comp1.pzd.Jx2+sys
2.gSub23*root.comp1.
pzd.Jx3)+sys2.gSub23*r
ealdot(root.comp1.pzd.
Jx3,sys2.gSub31*root.c
omp1.pzd.Jx1+sys2.gSu
b32*root.comp1.pzd.Jx
2+sys2.gSub33*root.co
mp1.pzd.Jx3)+sys2.gSu
b31*realdot(root.comp
1.pzd.Jx1,sys2.gSub11*
A/m^2 Current density
norm
Domain 4
52

Name Expression Unit Description Selection
root.comp1.pzd.Jx1+sys
2.gSub12*root.comp1.
pzd.Jx2+sys2.gSub13*r
oot.comp1.pzd.Jx3)+sys
2.gSub32*realdot(root.
comp1.pzd.Jx2,sys2.gSu
b21*root.comp1.pzd.Jx
1+sys2.gSub22*root.co
mp1.pzd.Jx2+sys2.gSub
23*root.comp1.pzd.Jx3
)+sys2.gSub33*realdot(
root.comp1.pzd.Jx3,sys
2.gSub31*root.comp1.
pzd.Jx1+sys2.gSub32*r
oot.comp1.pzd.Jx2+sys
2.gSub33*root.comp1.
pzd.Jx3))
pzd.EX -VX V/m Electric field, X
component
Domain 4
pzd.EY -VY V/m Electric field, Y
component
Domain 4
pzd.EZ -VZ V/m Electric field, Z
component
Domain 4
pzd.tEX -VTX V/m Tangential
electric field, X
component
Boundaries
1419
pzd.tEY -VTY V/m Tangential
electric field, Y
component
Boundaries
1419
pzd.tEZ -VTZ V/m Tangential
electric field, Z
component
Boundaries
1419
pzd.JX 0 A/m^2 Current density,
X component
Domain 4
pzd.JY 0 A/m^2 Current density,
Y component
Domain 4
pzd.JZ 0 A/m^2 Current density,
Z component
Domain 4
pzd.unTeX -
0.5*dnX*(up(pzd.Dx)*u
p(pzd.Ex)+up(pzd.Dy)*u
p(pzd.Ey)+up(pzd.Dz)*u
Pa Maxwell upward
electric surface
stress tensor, X
component
Boundaries
1419
53

Name Expression Unit Description Selection
p(pzd.Ez))+up(pzd.Dx)*(
up(pzd.Ex)*dnX+up(pzd
.Ey)*dnY+up(pzd.Ez)*d
nZ)
pzd.unTeY -
0.5*dnY*(up(pzd.Dx)*u
p(pzd.Ex)+up(pzd.Dy)*u
p(pzd.Ey)+up(pzd.Dz)*u
p(pzd.Ez))+up(pzd.Dy)*(
up(pzd.Ex)*dnX+up(pzd
.Ey)*dnY+up(pzd.Ez)*d
nZ)
Pa Maxwell upward
electric surface
stress tensor, Y
component
Boundaries
1419
pzd.unTeZ -
0.5*dnZ*(up(pzd.Dx)*u
p(pzd.Ex)+up(pzd.Dy)*u
p(pzd.Ey)+up(pzd.Dz)*u
p(pzd.Ez))+up(pzd.Dz)*(
up(pzd.Ex)*dnX+up(pzd
.Ey)*dnY+up(pzd.Ez)*d
nZ)
Pa Maxwell upward
electric surface
stress tensor, Z
component
Boundaries
1419
pzd.dnTeX -
0.5*unX*(down(pzd.Dx)
*down(pzd.Ex)+down(p
zd.Dy)*down(pzd.Ey)+d
own(pzd.Dz)*down(pzd
.Ez))+down(pzd.Dx)*(do
wn(pzd.Ex)*unX+down(
pzd.Ey)*unY+down(pzd.
Ez)*unZ)
Pa Maxwell
downward
electric surface
stress tensor, X
component
Boundaries
1419
pzd.dnTeY -
0.5*unY*(down(pzd.Dx)
*down(pzd.Ex)+down(p
zd.Dy)*down(pzd.Ey)+d
own(pzd.Dz)*down(pzd
.Ez))+down(pzd.Dy)*(do
wn(pzd.Ex)*unX+down(
pzd.Ey)*unY+down(pzd.
Ez)*unZ)
Pa Maxwell
downward
electric surface
stress tensor, Y
component
Boundaries
1419
pzd.dnTeZ -
0.5*unZ*(down(pzd.Dx)
*down(pzd.Ex)+down(p
zd.Dy)*down(pzd.Ey)+d
own(pzd.Dz)*down(pzd
.Ez))+down(pzd.Dz)*(do
Pa Maxwell
downward
electric surface
stress tensor, Z
component
Boundaries
1419
54

Name Expression Unit Description Selection
wn(pzd.Ex)*unX+down(
pzd.Ey)*unY+down(pzd.
Ez)*unZ)
pzd.unTX 0 Pa Maxwell upward
surface stress
tensor, X
component
Boundaries
1419
pzd.unTY 0 Pa Maxwell upward
surface stress
tensor, Y
component
Boundaries
1419
pzd.unTZ 0 Pa Maxwell upward
surface stress
tensor, Z
component
Boundaries
1419
pzd.dnTX pzd.dnTeX Pa Maxwell
downward
surface stress
tensor, X
component
Boundaries
1419
pzd.dnTY pzd.dnTeY Pa Maxwell
downward
surface stress
tensor, Y
component
Boundaries
1419
pzd.dnTZ pzd.dnTeZ Pa Maxwell
downward
surface stress
tensor, Z
component
Boundaries
1419
pzd.We 0.5*(pzd.Dx*pzd.Ex+pz
d.Dy*pzd.Ey+pzd.Dz*pz
d.Ez)
J/m^3 Electric energy
density
Domain 4
pzd.W 0.5*(pzd.Dx*pzd.Ex+pz
d.Dy*pzd.Ey+pzd.Dz*pz
d.Ez)
J/m^3 Energy density Domain 4
pzd.intWe pzd.int7(0.5*(pzd.Dx*p
zd.Ex+pzd.Dy*pzd.Ey+p
zd.Dz*pzd.Ez))
J Total electric
energy
Global

55

Shape functions
Name Shape
function
Unit Description Shape
frame
Selection
u Lagrange
(Quadratic)
m Displacement field, X component Material Domain 4
v Lagrange
(Quadratic)
m Displacement field, Y component Material Domain 4
w Lagrange
(Quadratic)
m Displacement field, Z component Material Domain 4
V Lagrange
(Quadratic)
V Electric potential Material Domain 4
V Lagrange
(Quadratic)
V Electric potential Spatial Domain 4

Weak expressions
Weak expression Integration frame Selection
-pzd.Sl11*test(pzd.el11)-2*pzd.Sl12*test(pzd.el12)-
2*pzd.Sl13*test(pzd.el13)-pzd.Sl22*test(pzd.el22)-
2*pzd.Sl23*test(pzd.el23)-pzd.Sl33*test(pzd.el33)
Material Domain 4
-pzd.PX*test(VX)-pzd.PY*test(VY)-pzd.PZ*test(VZ) Material Domain 4
epsilon0_const*(-VX*test(-VX)-VY*test(-VY)-VZ*test(-VZ)) Material Domain 4

1.4.2 Free 1

Free 1
Selection
56

Geometric entity level Boundary
Selection Boundaries 23, 5, 8, 1013, 15, 18, 20, 2326

1.4.3 Zero Charge 1

Zero Charge 1
Selection
Geometric entity level Boundary
Selection Boundaries 1415, 1819

Equations

Variables
Name Expression Unit Description Selection
VTXt root.comp1.VTXTIME m*kg/(s^4*A) Tangential
gradient of V,
first time
derivative, X
component
Boundaries
14, 19
VTYt root.comp1.VTYTIME m*kg/(s^4*A) Tangential
gradient of V,
first time
derivative, Y
component
Boundaries
14, 19
VTZt root.comp1.VTZTIME m*kg/(s^4*A) Tangential
gradient of V,
first time
Boundaries
14, 19
57

Name Expression Unit Description Selection
derivative, Z
component
VTXtt root.comp1.VTXTIMETIME m*kg/(s^5*A) Tangential
gradient of V,
second time
derivative, X
component
Boundaries
14, 19
VTYtt root.comp1.VTYTIMETIME m*kg/(s^5*A) Tangential
gradient of V,
second time
derivative, Y
component
Boundaries
14, 19
VTZtt root.comp1.VTZTIMETIME m*kg/(s^5*A) Tangential
gradient of V,
second time
derivative, Z
component
Boundaries
14, 19
Vt root.comp1.VTIME V/s Electric
potential, first
time derivative
Boundaries
14, 19
Vtt root.comp1.VTIMETIME V/s^2 Electric
potential,
second time
derivative
Boundaries
14, 19
pzd.nD 0 C/m^2 Surface charge
density
Boundaries
1415, 1819

Shape functions
Name Shape function Unit Description Shape frame Selection
V Lagrange
(Quadratic)
V Electric potential Material Boundaries 14, 19

58

1.4.4 Initial Values 1

Initial Values 1
Selection
Geometric entity level Domain
Selection Domains 15

Settings
Settings
Description Value
Displacement field {0, 0, 0}
Structural velocity field {0, 0, 0}
Electric potential 0

Variables
Name Expression Unit Description Selection
pzd.uInitx 0 m Initial value of displacement, x
component
Domains 15
pzd.uInity 0 m Initial value of displacement, y
component
Domains 15
pzd.uInitz 0 m Initial value of displacement, z
component
Domains 15
pzd.utInitx 0 m/s Initial value of structural velocity, x
component
Domains 15
pzd.utInity 0 m/s Initial value of structural velocity, y Domains 15
59

Name Expression Unit Description Selection
component
pzd.utInitz 0 m/s Initial value of structural velocity, z
component
Domains 15

1.4.5 Linear Elastic Material 1

Linear Elastic Material 1
Selection
Geometric entity level Domain
Selection Domains 13, 5

Equations



Settings
Settings
Description Value
Young's modulus From material
Poisson's ratio From material
Elasticity matrix {{0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0}, {0, 0, 0, 0,
0, 0}, {0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0}}
Elasticity matrix, Voigt notation {{0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0}, {0, 0, 0, 0,
60

Description Value
0, 0}, {0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0}}
Density From material
Solid model Isotropic
Specify Young's modulus and Poisson's ratio
Standard (XX, YY, ZZ, XY, YZ, XZ)
Nearly incompressible material Off
Energy dissipation 0
Calculate dissipated energy 0

Properties from material
Property Material Property group
Young's modulus Al Young's modulus and Poisson's ratio
Poisson's ratio Al Young's modulus and Poisson's ratio
Density Al Basic
Young's modulus Material 3 Basic
Poisson's ratio Material 3 Basic
Density Material 3 Basic

Variables
Name Expression Unit Description Selection
uTXt root.comp1.uTXTIME 1/s Tangential
gradient of u,
first time
derivative, X
component
Boundaries
114, 1617,
1926
uTYt root.comp1.uTYTIME 1/s Tangential
gradient of u,
first time
derivative, Y
component
Boundaries
114, 1617,
1926
uTZt root.comp1.uTZTIME 1/s Tangential
gradient of u,
first time
derivative, Z
component
Boundaries
114, 1617,
1926
uTXtt root.comp1.uTXTIMETIME 1/s^2 Tangential
gradient of u,
second time
Boundaries
114, 1617,
61

Name Expression Unit Description Selection
derivative, X
component
1926
uTYtt root.comp1.uTYTIMETIME 1/s^2 Tangential
gradient of u,
second time
derivative, Y
component
Boundaries
114, 1617,
1926
uTZtt root.comp1.uTZTIMETIME 1/s^2 Tangential
gradient of u,
second time
derivative, Z
component
Boundaries
114, 1617,
1926
uXt root.comp1.uXTIME 1/s Gradient of u,
first time
derivative, X
component
Domains 13,
5
uYt root.comp1.uYTIME 1/s Gradient of u,
first time
derivative, Y
component
Domains 13,
5
uZt root.comp1.uZTIME 1/s Gradient of u,
first time
derivative, Z
component
Domains 13,
5
uXtt root.comp1.uXTIMETIME 1/s^2 Gradient of u,
second time
derivative, X
component
Domains 13,
5
uYtt root.comp1.uYTIMETIME 1/s^2 Gradient of u,
second time
derivative, Y
component
Domains 13,
5
uZtt root.comp1.uZTIMETIME 1/s^2 Gradient of u,
second time
derivative, Z
component
Domains 13,
5
vTXt root.comp1.vTXTIME 1/s Tangential
gradient of v,
first time
derivative, X
component
Boundaries
114, 1617,
1926
vTYt root.comp1.vTYTIME 1/s Tangential Boundaries
62

Name Expression Unit Description Selection
gradient of v,
first time
derivative, Y
component
114, 1617,
1926
vTZt root.comp1.vTZTIME 1/s Tangential
gradient of v,
first time
derivative, Z
component
Boundaries
114, 1617,
1926
vTXtt root.comp1.vTXTIMETIME 1/s^2 Tangential
gradient of v,
second time
derivative, X
component
Boundaries
114, 1617,
1926
vTYtt root.comp1.vTYTIMETIME 1/s^2 Tangential
gradient of v,
second time
derivative, Y
component
Boundaries
114, 1617,
1926
vTZtt root.comp1.vTZTIMETIME 1/s^2 Tangential
gradient of v,
second time
derivative, Z
component
Boundaries
114, 1617,
1926
vXt root.comp1.vXTIME 1/s Gradient of v,
first time
derivative, X
component
Domains 13,
5
vYt root.comp1.vYTIME 1/s Gradient of v,
first time
derivative, Y
component
Domains 13,
5
vZt root.comp1.vZTIME 1/s Gradient of v,
first time
derivative, Z
component
Domains 13,
5
vXtt root.comp1.vXTIMETIME 1/s^2 Gradient of v,
second time
derivative, X
component
Domains 13,
5
vYtt root.comp1.vYTIMETIME 1/s^2 Gradient of v,
second time
derivative, Y
Domains 13,
5
63

Name Expression Unit Description Selection
component
vZtt root.comp1.vZTIMETIME 1/s^2 Gradient of v,
second time
derivative, Z
component
Domains 13,
5
wTXt root.comp1.wTXTIME 1/s Tangential
gradient of w,
first time
derivative, X
component
Boundaries
114, 1617,
1926
wTYt root.comp1.wTYTIME 1/s Tangential
gradient of w,
first time
derivative, Y
component
Boundaries
114, 1617,
1926
wTZt root.comp1.wTZTIME 1/s Tangential
gradient of w,
first time
derivative, Z
component
Boundaries
114, 1617,
1926
wTXtt root.comp1.wTXTIMETIME 1/s^2 Tangential
gradient of w,
second time
derivative, X
component
Boundaries
114, 1617,
1926
wTYtt root.comp1.wTYTIMETIME 1/s^2 Tangential
gradient of w,
second time
derivative, Y
component
Boundaries
114, 1617,
1926
wTZtt root.comp1.wTZTIMETIME 1/s^2 Tangential
gradient of w,
second time
derivative, Z
component
Boundaries
114, 1617,
1926
wXt root.comp1.wXTIME 1/s Gradient of w,
first time
derivative, X
component
Domains 13,
5
wYt root.comp1.wYTIME 1/s Gradient of w,
first time
derivative, Y
Domains 13,
5
64

Name Expression Unit Description Selection
component
wZt root.comp1.wZTIME 1/s Gradient of w,
first time
derivative, Z
component
Domains 13,
5
wXtt root.comp1.wXTIMETIME 1/s^2 Gradient of w,
second time
derivative, X
component
Domains 13,
5
wYtt root.comp1.wYTIMETIME 1/s^2 Gradient of w,
second time
derivative, Y
component
Domains 13,
5
wZtt root.comp1.wZTIMETIME 1/s^2 Gradient of w,
second time
derivative, Z
component
Domains 13,
5
ut root.comp1.uTIME m/s Structural
velocity field, X
component
Domains 13,
5
vt root.comp1.vTIME m/s Structural
velocity field, Y
component
Domains 13,
5
wt root.comp1.wTIME m/s Structural
velocity field, Z
component
Domains 13,
5
utt root.comp1.uTIMETIME m/s^2 Acceleration
field, X
component
Domains 13,
5
vtt root.comp1.vTIMETIME m/s^2 Acceleration
field, Y
component
Domains 13,
5
wtt root.comp1.wTIMETIME m/s^2 Acceleration
field, Z
component
Domains 13,
5
pzd.Eequ pzd.E Pa Equivalent
Young's
modulus
Domains 13,
5
pzd.nuequ pzd.nu 1 Equivalent
Poisson's ratio
Domains 13,
5
65

Name Expression Unit Description Selection
pzd.cp sqrt(pzd.E*(1+pzd.nu/(1-
2*pzd.nu))/((1+pzd.nu)*pzd.r
ho))
m/s Pressure-wave
speed
Domains 13,
5
pzd.rho material.rho kg/m^3 Density Domains 13,
5
pzd.FdxX 1+uX 1 Deformation
gradient, xX
component
Domains 13,
5
pzd.FdyX vX 1 Deformation
gradient, yX
component
Domains 13,
5
pzd.FdzX wX 1 Deformation
gradient, zX
component
Domains 13,
5
pzd.FdxY uY 1 Deformation
gradient, xY
component
Domains 13,
5
pzd.FdyY 1+vY 1 Deformation
gradient, yY
component
Domains 13,
5
pzd.FdzY wY 1 Deformation
gradient, zY
component
Domains 13,
5
pzd.FdxZ uZ 1 Deformation
gradient, xZ
component
Domains 13,
5
pzd.FdyZ vZ 1 Deformation
gradient, yZ
component
Domains 13,
5
pzd.FdzZ 1+wZ 1 Deformation
gradient, zZ
component
Domains 13,
5
pzd.Fdlx1 1+uX 1 Deformation
gradient, local,
x1 component
Domains 13,
5
pzd.Fdly1 vX 1 Deformation
gradient, local,
y1 component
Domains 13,
5
pzd.Fdlz1 wX 1 Deformation
gradient, local,
Domains 13,
5
66

Name Expression Unit Description Selection
z1 component
pzd.Fdlx2 uY 1 Deformation
gradient, local,
x2 component
Domains 13,
5
pzd.Fdly2 1+vY 1 Deformation
gradient, local,
y2 component
Domains 13,
5
pzd.Fdlz2 wY 1 Deformation
gradient, local,
z2 component
Domains 13,
5
pzd.Fdlx3 uZ 1 Deformation
gradient, local,
x3 component
Domains 13,
5
pzd.Fdly3 vZ 1 Deformation
gradient, local,
y3 component
Domains 13,
5
pzd.Fdlz3 1+wZ 1 Deformation
gradient, local,
z3 component
Domains 13,
5
pzd.gradUxX uX 1 Displacement
gradient, xX
component
Domains 13,
5
pzd.gradUyX vX 1 Displacement
gradient, yX
component
Domains 13,
5
pzd.gradUzX wX 1 Displacement
gradient, zX
component
Domains 13,
5
pzd.gradUxY uY 1 Displacement
gradient, xY
component
Domains 13,
5
pzd.gradUyY vY 1 Displacement
gradient, yY
component
Domains 13,
5
pzd.gradUzY wY 1 Displacement
gradient, zY
component
Domains 13,
5
pzd.gradUxZ uZ 1 Displacement
gradient, xZ
component
Domains 13,
5
67

Name Expression Unit Description Selection
pzd.gradUyZ vZ 1 Displacement
gradient, yZ
component
Domains 13,
5
pzd.gradUzZ wZ 1 Displacement
gradient, zZ
component
Domains 13,
5
pzd.eXX uX 1 Strain tensor,
XX component
Domains 13,
5
pzd.eXY 0.5*(uY+vX) 1 Strain tensor,
XY component
Domains 13,
5
pzd.eXZ 0.5*(uZ+wX) 1 Strain tensor,
XZ component
Domains 13,
5
pzd.eYY vY 1 Strain tensor,
YY component
Domains 13,
5
pzd.eYZ 0.5*(vZ+wY) 1 Strain tensor,
YZ component
Domains 13,
5
pzd.eZZ wZ 1 Strain tensor, ZZ
component
Domains 13,
5
pzd.eX pzd.eXX 1 Strain tensor, X
component
Domains 13,
5
pzd.eY pzd.eYY 1 Strain tensor, Y
component
Domains 13,
5
pzd.eZ pzd.eZZ 1 Strain tensor, Z
component
Domains 13,
5
pzd.Ldx d(pzd.gradUxX,TIME) 1/s Rate of strain
tensor, x
component
Domains 13,
5
pzd.Ldxy 0.5*(d(pzd.gradUxY,TIME)+d(
pzd.gradUyX,TIME))
1/s Rate of strain
tensor, xy
component
Domains 13,
5
pzd.Ldxz 0.5*(d(pzd.gradUxZ,TIME)+d(
pzd.gradUzX,TIME))
1/s Rate of strain
tensor, xz
component
Domains 13,
5
pzd.Ldy d(pzd.gradUyY,TIME) 1/s Rate of strain
tensor, y
component
Domains 13,
5
pzd.Ldyz 0.5*(d(pzd.gradUyZ,TIME)+d(
pzd.gradUzY,TIME))
1/s Rate of strain
tensor, yz
component
Domains 13,
5
68

Name Expression Unit Description Selection
pzd.Ldz d(pzd.gradUzZ,TIME) 1/s Rate of strain
tensor, z
component
Domains 13,
5
pzd.Lwx 0 1/s Spin tensor, x
component
Domains 13,
5
pzd.Lwxy 0.5*(d(pzd.gradUxY,TIME)-
d(pzd.gradUyX,TIME))
1/s Spin tensor, xy
component
Domains 13,
5
pzd.Lwxz 0.5*(d(pzd.gradUxZ,TIME)-
d(pzd.gradUzX,TIME))
1/s Spin tensor, xz
component
Domains 13,
5
pzd.Lwy 0 1/s Spin tensor, y
component
Domains 13,
5
pzd.Lwyz 0.5*(d(pzd.gradUyZ,TIME)-
d(pzd.gradUzY,TIME))
1/s Spin tensor, yz
component
Domains 13,
5
pzd.Lwz 0 1/s Spin tensor, z
component
Domains 13,
5
pzd.el11 uX 1 Strain tensor,
local coordinate
system, 11
component
Domains 13,
5
pzd.el12 0.5*(uY+vX) 1 Strain tensor,
local coordinate
system, 12
component
Domains 13,
5
pzd.el13 0.5*(uZ+wX) 1 Strain tensor,
local coordinate
system, 13
component
Domains 13,
5
pzd.el22 vY 1 Strain tensor,
local coordinate
system, 22
component
Domains 13,
5
pzd.el23 0.5*(vZ+wY) 1 Strain tensor,
local coordinate
system, 23
component
Domains 13,
5
pzd.el33 wZ 1 Strain tensor,
local coordinate
system, 33
component
Domains 13,
5
pzd.eil11 0 1 Initial strain
tensor, local
Domains 13,
69

Name Expression Unit Description Selection
coordinate
system, 11
component
5
pzd.eil12 0 1 Initial strain
tensor, local
coordinate
system, 12
component
Domains 13,
5
pzd.eil13 0 1 Initial strain
tensor, local
coordinate
system, 13
component
Domains 13,
5
pzd.eil22 0 1 Initial strain
tensor, local
coordinate
system, 22
component
Domains 13,
5
pzd.eil23 0 1 Initial strain
tensor, local
coordinate
system, 23
component
Domains 13,
5
pzd.eil33 0 1 Initial strain
tensor, local
coordinate
system, 33
component
Domains 13,
5
pzd.eiel11 0 1 Inelastic strain
tensor, local
coordinate
system, 11
component
Domains 13,
5
pzd.eiel12 0 1 Inelastic strain
tensor, local
coordinate
system, 12
component
Domains 13,
5
pzd.eiel13 0 1 Inelastic strain
tensor, local
coordinate
system, 13
component
Domains 13,
5
70

Name Expression Unit Description Selection
pzd.eiel22 0 1 Inelastic strain
tensor, local
coordinate
system, 22
component
Domains 13,
5
pzd.eiel23 0 1 Inelastic strain
tensor, local
coordinate
system, 23
component
Domains 13,
5
pzd.eiel33 0 1 Inelastic strain
tensor, local
coordinate
system, 33
component
Domains 13,
5
pzd.eel11 uX-pzd.eil11-pzd.eiel11 1 Elastic strain
tensor, local
coordinate
system, 11
component
Domains 13,
5
pzd.eel12 0.5*(uY+vX)-pzd.eil12-
pzd.eiel12
1 Elastic strain
tensor, local
coordinate
system, 12
component
Domains 13,
5
pzd.eel13 0.5*(uZ+wX)-pzd.eil13-
pzd.eiel13
1 Elastic strain
tensor, local
coordinate
system, 13
component
Domains 13,
5
pzd.eel22 vY-pzd.eil22-pzd.eiel22 1 Elastic strain
tensor, local
coordinate
system, 22
component
Domains 13,
5
pzd.eel23 0.5*(vZ+wY)-pzd.eil23-
pzd.eiel23
1 Elastic strain
tensor, local
coordinate
system, 23
component
Domains 13,
5
pzd.eel33 wZ-pzd.eil33-pzd.eiel33 1 Elastic strain
tensor, local
coordinate
Domains 13,
5
71

Name Expression Unit Description Selection
system, 33
component
pzd.evol uX+vY+wZ 1 Volumetric
strain
Domains 13,
5
pzd.J pzd.FdxX*pzd.FdyY*pzd.FdzZ
+pzd.FdxY*pzd.FdyZ*pzd.Fdz
X+pzd.FdxZ*pzd.FdyX*pzd.Fd
zY-
pzd.FdxX*pzd.FdyZ*pzd.FdzY
-
pzd.FdxY*pzd.FdyX*pzd.FdzZ
-
pzd.FdxZ*pzd.FdyY*pzd.FdzX
1 Volume ratio Domains 13,
5
pzd.Sl11 pzd.D11*pzd.eel11+2*pzd.D
14*pzd.eel12+2*pzd.D16*pz
d.eel13+pzd.D12*pzd.eel22+
2*pzd.D15*pzd.eel23+pzd.D
13*pzd.eel33+pzd.Sil11
N/m^2 Second Piola-
Kirchhoff stress,
local coordinate
system, 11
component
Domains 13,
5
pzd.Sl12 pzd.D14*pzd.eel11+2*pzd.D
44*pzd.eel12+2*pzd.D46*pz
d.eel13+pzd.D24*pzd.eel22+
2*pzd.D45*pzd.eel23+pzd.D
34*pzd.eel33+pzd.Sil12
N/m^2 Second Piola-
Kirchhoff stress,
local coordinate
system, 12
component
Domains 13,
5
pzd.Sl13 pzd.D16*pzd.eel11+2*pzd.D
46*pzd.eel12+2*pzd.D66*pz
d.eel13+pzd.D26*pzd.eel22+
2*pzd.D56*pzd.eel23+pzd.D
36*pzd.eel33+pzd.Sil13
N/m^2 Second Piola-
Kirchhoff stress,
local coordinate
system, 13
component
Domains 13,
5
pzd.Sl22 pzd.D12*pzd.eel11+2*pzd.D
24*pzd.eel12+2*pzd.D26*pz
d.eel13+pzd.D22*pzd.eel22+
2*pzd.D25*pzd.eel23+pzd.D
23*pzd.eel33+pzd.Sil22
N/m^2 Second Piola-
Kirchhoff stress,
local coordinate
system, 22
component
Domains 13,
5
pzd.Sl23 pzd.D15*pzd.eel11+2*pzd.D
45*pzd.eel12+2*pzd.D56*pz
d.eel13+pzd.D25*pzd.eel22+
2*pzd.D55*pzd.eel23+pzd.D
35*pzd.eel33+pzd.Sil23
N/m^2 Second Piola-
Kirchhoff stress,
local coordinate
system, 23
component
Domains 13,
5
pzd.Sl33 pzd.D13*pzd.eel11+2*pzd.D
34*pzd.eel12+2*pzd.D36*pz
d.eel13+pzd.D23*pzd.eel22+
2*pzd.D35*pzd.eel23+pzd.D
33*pzd.eel33+pzd.Sil33
N/m^2 Second Piola-
Kirchhoff stress,
local coordinate
system, 33
component
Domains 13,
5
72

Name Expression Unit Description Selection
pzd.Sil11 0 N/m^2 Initial stress
tensor, local
coordinate
system, 11
component
Domains 13,
5
pzd.Sil12 0 N/m^2 Initial stress
tensor, local
coordinate
system, 12
component
Domains 13,
5
pzd.Sil13 0 N/m^2 Initial stress
tensor, local
coordinate
system, 13
component
Domains 13,
5
pzd.Sil22 0 N/m^2 Initial stress
tensor, local
coordinate
system, 22
component
Domains 13,
5
pzd.Sil23 0 N/m^2 Initial stress
tensor, local
coordinate
system, 23
component
Domains 13,
5
pzd.Sil33 0 N/m^2 Initial stress
tensor, local
coordinate
system, 33
component
Domains 13,
5
pzd.u_ttX d(root.comp1.uTIME,TIME) m/s^2 Acceleration, X
component
Domains 13,
5
pzd.u_ttY d(root.comp1.vTIME,TIME) m/s^2 Acceleration, Y
component
Domains 13,
5
pzd.u_ttZ d(root.comp1.wTIME,TIME) m/s^2 Acceleration, Z
component
Domains 13,
5
pzd.accX d(root.comp1.uTIME,TIME) m/s^2 Total
acceleration, X
component
Domains 13,
5
pzd.accY d(root.comp1.vTIME,TIME) m/s^2 Total
acceleration, Y
component
Domains 13,
5
73

Name Expression Unit Description Selection
pzd.accZ d(root.comp1.wTIME,TIME) m/s^2 Total
acceleration, Z
component
Domains 13,
5
pzd.u_tX root.comp1.uTIME m/s Velocity, X
component
Domains 13,
5
pzd.u_tY root.comp1.vTIME m/s Velocity, Y
component
Domains 13,
5
pzd.u_tZ root.comp1.wTIME m/s Velocity, Z
component
Domains 13,
5
pzd.disp sqrt(real(u)^2+real(v)^2+real
(w)^2)
m Total
displacement
Domains 13,
5
pzd.SX pzd.Sl11 N/m^2 Second Piola-
Kirchhoff stress,
X component
Domains 13,
5
pzd.SXY pzd.Sl12 N/m^2 Second Piola-
Kirchhoff stress,
XY component
Domains 13,
5
pzd.SXZ pzd.Sl13 N/m^2 Second Piola-
Kirchhoff stress,
XZ component
Domains 13,
5
pzd.SY pzd.Sl22 N/m^2 Second Piola-
Kirchhoff stress,
Y component
Domains 13,
5
pzd.SYZ pzd.Sl23 N/m^2 Second Piola-
Kirchhoff stress,
YZ component
Domains 13,
5
pzd.SZ pzd.Sl33 N/m^2 Second Piola-
Kirchhoff stress,
Z component
Domains 13,
5
pzd.PxX pzd.FdxX*pzd.SX+pzd.FdxY*p
zd.SXY+pzd.FdxZ*pzd.SXZ
N/m^2 First Piola-
Kirchhoff stress,
xX component
Domains 13,
5
pzd.PyX pzd.FdyX*pzd.SX+pzd.FdyY*p
zd.SXY+pzd.FdyZ*pzd.SXZ
N/m^2 First Piola-
Kirchhoff stress,
yX component
Domains 13,
5
pzd.PzX pzd.FdzX*pzd.SX+pzd.FdzY*p
zd.SXY+pzd.FdzZ*pzd.SXZ
N/m^2 First Piola-
Kirchhoff stress,
zX component
Domains 13,
5
pzd.PxY pzd.FdxX*pzd.SXY+pzd.FdxY*
pzd.SY+pzd.FdxZ*pzd.SYZ
N/m^2 First Piola-
Kirchhoff stress,
Domains 13,
5
74

Name Expression Unit Description Selection
xY component
pzd.PyY pzd.FdyX*pzd.SXY+pzd.FdyY*
pzd.SY+pzd.FdyZ*pzd.SYZ
N/m^2 First Piola-
Kirchhoff stress,
yY component
Domains 13,
5
pzd.PzY pzd.FdzX*pzd.SXY+pzd.FdzY*
pzd.SY+pzd.FdzZ*pzd.SYZ
N/m^2 First Piola-
Kirchhoff stress,
zY component
Domains 13,
5
pzd.PxZ pzd.FdxX*pzd.SXZ+pzd.FdxY*
pzd.SYZ+pzd.FdxZ*pzd.SZ
N/m^2 First Piola-
Kirchhoff stress,
xZ component
Domains 13,
5
pzd.PyZ pzd.FdyX*pzd.SXZ+pzd.FdyY*
pzd.SYZ+pzd.FdyZ*pzd.SZ
N/m^2 First Piola-
Kirchhoff stress,
yZ component
Domains 13,
5
pzd.PzZ pzd.FdzX*pzd.SXZ+pzd.FdzY*
pzd.SYZ+pzd.FdzZ*pzd.SZ
N/m^2 First Piola-
Kirchhoff stress,
zZ component
Domains 13,
5
pzd.sx pzd.SX N/m^2 Stress tensor, x
component
Domains 13,
5
pzd.sxy pzd.SXY N/m^2 Stress tensor, xy
component
Domains 13,
5
pzd.sxz pzd.SXZ N/m^2 Stress tensor, xz
component
Domains 13,
5
pzd.sy pzd.SY N/m^2 Stress tensor, y
component
Domains 13,
5
pzd.syz pzd.SYZ N/m^2 Stress tensor, yz
component
Domains 13,
5
pzd.sz pzd.SZ N/m^2 Stress tensor, z
component
Domains 13,
5
pzd.sl11 pzd.Sl11 N/m^2 Stress tensor,
local coordinate
system, 11
component
Domains 13,
5
pzd.sl12 pzd.Sl12 N/m^2 Stress tensor,
local coordinate
system, 12
component
Domains 13,
5
pzd.sl13 pzd.Sl13 N/m^2 Stress tensor,
local coordinate
system, 13
component
Domains 13,
5
75

Name Expression Unit Description Selection
pzd.sl22 pzd.Sl22 N/m^2 Stress tensor,
local coordinate
system, 22
component
Domains 13,
5
pzd.sl23 pzd.Sl23 N/m^2 Stress tensor,
local coordinate
system, 23
component
Domains 13,
5
pzd.sl33 pzd.Sl33 N/m^2 Stress tensor,
local coordinate
system, 33
component
Domains 13,
5
pzd.sdevx pzd.sx-
(pzd.sx+pzd.sy+pzd.sz)/3
N/m^2 Deviatoric
stress tensor, x
component
Domains 13,
5
pzd.sdevxy pzd.sxy N/m^2 Deviatoric
stress tensor, xy
component
Domains 13,
5
pzd.sdevxz pzd.sxz N/m^2 Deviatoric
stress tensor, xz
component
Domains 13,
5
pzd.sdevy pzd.sy-
(pzd.sx+pzd.sy+pzd.sz)/3
N/m^2 Deviatoric
stress tensor, y
component
Domains 13,
5
pzd.sdevyz pzd.syz N/m^2 Deviatoric
stress tensor, yz
component
Domains 13,
5
pzd.sdevz pzd.sz-
(pzd.sx+pzd.sy+pzd.sz)/3
N/m^2 Deviatoric
stress tensor, z
component
Domains 13,
5
pzd.SdevX pzd.SX-
(pzd.SX+pzd.SY+pzd.SZ)/3
N/m^2 Deviatoric
second Piola-
Kirchhoff stress,
X component
Domains 13,
5
pzd.SdevXY pzd.SXY N/m^2 Deviatoric
second Piola-
Kirchhoff stress,
XY component
Domains 13,
5
pzd.SdevXZ pzd.SXZ N/m^2 Deviatoric
second Piola-
Kirchhoff stress,
Domains 13,
5
76

Name Expression Unit Description Selection
XZ component
pzd.SdevY pzd.SY-
(pzd.SX+pzd.SY+pzd.SZ)/3
N/m^2 Deviatoric
second Piola-
Kirchhoff stress,
Y component
Domains 13,
5
pzd.SdevYZ pzd.SYZ N/m^2 Deviatoric
second Piola-
Kirchhoff stress,
YZ component
Domains 13,
5
pzd.SdevZ pzd.SZ-
(pzd.SX+pzd.SY+pzd.SZ)/3
N/m^2 Deviatoric
second Piola-
Kirchhoff stress,
Z component
Domains 13,
5
pzd.I1s pzd.sl11+pzd.sl22+pzd.sl33 N/m^2 First principal
invariant of
stress
Domains 13,
5
pzd.I2s 0.5*(pzd.I1s^2-pzd.sl11^2-
2*pzd.sl12^2-2*pzd.sl13^2-
pzd.sl22^2-2*pzd.sl23^2-
pzd.sl33^2)
kg^2/(m^2*s^4) Second
principal
invariant of
stress
Domains 13,
5
pzd.I3s pzd.sl11*pzd.sl22*pzd.sl33+2
*pzd.sl12*pzd.sl23*pzd.sl13-
pzd.sl11*pzd.sl23^2-
pzd.sl12^2*pzd.sl33-
pzd.sl22*pzd.sl13^2
kg^3/(m^3*s^6) Third principal
invariant of
stress
Domains 13,
5
pzd.II2s 0.5*((pzd.sl11-
(pzd.sl11+pzd.sl22+pzd.sl33)/
3)^2+2*pzd.sl12^2+2*pzd.sl1
3^2+(pzd.sl22-
(pzd.sl11+pzd.sl22+pzd.sl33)/
3)^2+2*pzd.sl23^2+(pzd.sl33
-
(pzd.sl11+pzd.sl22+pzd.sl33)/
3)^2)
kg^2/(m^2*s^4) Second
invariant of
stress deviator
Domains 13,
5
pzd.II3s (pzd.sl11-
(pzd.sl11+pzd.sl22+pzd.sl33)/
3)*(pzd.sl22-
(pzd.sl11+pzd.sl22+pzd.sl33)/
3)*(pzd.sl33-
(pzd.sl11+pzd.sl22+pzd.sl33)/
3)+2*pzd.sl12*pzd.sl23*pzd.
sl13-(pzd.sl11-
(pzd.sl11+pzd.sl22+pzd.sl33)/
kg^3/(m^3*s^6) Third invariant
of stress
deviator
Domains 13,
5
77

Name Expression Unit Description Selection
3)*pzd.sl23^2-
pzd.sl12^2*(pzd.sl33-
(pzd.sl11+pzd.sl22+pzd.sl33)/
3)-(pzd.sl22-
(pzd.sl11+pzd.sl22+pzd.sl33)/
3)*pzd.sl13^2
pzd.thetaL atan2(sqrt(max(0.148148148
14814814*pzd.II2s^3-
pzd.II3s^2,eps)),pzd.II3s)/3
rad Lode angle Domains 13,
5
pzd.eelvol pzd.eel11+pzd.eel22+pzd.eel
33
1 First invariant
of elastic strain
tensor
Domains 13,
5
pzd.eeldev11 pzd.eel11-
(pzd.eel11+pzd.eel22+pzd.ee
l33)/3
1 Deviatoric
elastic strain
tensor, local
coordinate
system, 11
component
Domains 13,
5
pzd.eeldev12 pzd.eel12 1 Deviatoric
elastic strain
tensor, local
coordinate
system, 12
component
Domains 13,
5
pzd.eeldev13 pzd.eel13 1 Deviatoric
elastic strain
tensor, local
coordinate
system, 13
component
Domains 13,
5
pzd.eeldev22 pzd.eel22-
(pzd.eel11+pzd.eel22+pzd.ee
l33)/3
1 Deviatoric
elastic strain
tensor, local
coordinate
system, 22
component
Domains 13,
5
pzd.eeldev23 pzd.eel23 1 Deviatoric
elastic strain
tensor, local
coordinate
system, 23
component
Domains 13,
5
pzd.eeldev33 pzd.eel33- 1 Deviatoric Domains 13,
78

Name Expression Unit Description Selection
(pzd.eel11+pzd.eel22+pzd.ee
l33)/3
elastic strain
tensor, local
coordinate
system, 33
component
5
pzd.II2eel 0.5*((pzd.eel11-
(pzd.eel11+pzd.eel22+pzd.ee
l33)/3)^2+2*pzd.eel12^2+2*
pzd.eel13^2+(pzd.eel22-
(pzd.eel11+pzd.eel22+pzd.ee
l33)/3)^2+2*pzd.eel23^2+(pz
d.eel33-
(pzd.eel11+pzd.eel22+pzd.ee
l33)/3)^2)
1 Second
invariant of
deviatoric
elastic strain
tensor
Domains 13,
5
pzd.pm -
(pzd.sl11+pzd.sl22+pzd.sl33)/
3
N/m^2 Pressure Domains 13,
5
pzd.curlUX wY-vZ 1 Curl of
displacement, X
component
Domains 13,
5
pzd.curlUY -wX+uZ 1 Curl of
displacement, Y
component
Domains 13,
5
pzd.curlUZ vX-uY 1 Curl of
displacement, Z
component
Domains 13,
5
pzd.mises (3*pzd.II2s)^0.5 N/m^2 von Mises
stress
Domains 13,
5
pzd.tresca pzd.sp1-pzd.sp3 N/m^2 Tresca stress Domains 13,
5
pzd.RFx reacf(u) N Reaction force,
x component
Domains 13,
5
pzd.RFy reacf(v) N Reaction force,
y component
Domains 13,
5
pzd.RFz reacf(w) N Reaction force,
z component
Domains 13,
5
pzd.RMx pzd.RFz*(y-pzd.refpnty)-
pzd.RFy*(z-pzd.refpntz)
N*m Reaction
moment, x
component
Domains 13,
5
pzd.RMy -pzd.RFz*(x-
pzd.refpntx)+pzd.RFx*(z-
N*m Reaction
moment, y
Domains 13,
5
79

Name Expression Unit Description Selection
pzd.refpntz) component
pzd.RMz pzd.RFy*(x-pzd.refpntx)-
pzd.RFx*(y-pzd.refpnty)
N*m Reaction
moment, z
component
Domains 13,
5
pzd.Tax pzd.sx*pzd.nx+pzd.sxy*pzd.n
y+pzd.sxz*pzd.nz
N/m^2 Traction
(force/area), x
component
Boundaries
114, 1617,
1926
pzd.Tay pzd.sxy*pzd.nx+pzd.sy*pzd.n
y+pzd.syz*pzd.nz
N/m^2 Traction
(force/area), y
component
Boundaries
114, 1617,
1926
pzd.Taz pzd.sxz*pzd.nx+pzd.syz*pzd.
ny+pzd.sz*pzd.nz
N/m^2 Traction
(force/area), z
component
Boundaries
114, 1617,
1926
pzd.Ws 0.5*((pzd.Sl11+pzd.Sil11)*pz
d.eel11+2*(pzd.Sl12+pzd.Sil1
2)*pzd.eel12+2*(pzd.Sl13+pz
d.Sil13)*pzd.eel13+(pzd.Sl22
+pzd.Sil22)*pzd.eel22+2*(pz
d.Sl23+pzd.Sil23)*pzd.eel23+
(pzd.Sl33+pzd.Sil33)*pzd.eel
33)
J/m^3 Elastic strain
energy density
Domains 13,
5
pzd.Ws_tot root.comp1.pzd.lemm1.int10
(pzd.Ws)
J Total elastic
strain energy
Global
pzd.E material.E Pa Young's
modulus
Domains 13,
5
pzd.nu material.nu 1 Poisson's ratio Domains 13,
5
pzd.K pzd.E/(3-6*pzd.nu) N/m^2 Bulk modulus Domains 13,
5
pzd.G 0.5*pzd.E/(1+pzd.nu) N/m^2 Shear modulus Domains 13,
5
pzd.lambLame pzd.E*pzd.nu/((1+pzd.nu)*(1
-2*pzd.nu))
N/m^2 Lam
parameter
Domains 13,
5
pzd.muLame 0.5*pzd.E/(1+pzd.nu) N/m^2 Lam
parameter
Domains 13,
5
pzd.cs sqrt(0.5*pzd.E/((1+pzd.nu)*p
zd.rho))
m/s Shear-wave
speed
Domains 13,
5
pzd.D11 pzd.E*(1+pzd.nu/(1-
2*pzd.nu))/(1+pzd.nu)
Pa Elasticity
matrix, 11
component
Domains 13,
5
80

Name Expression Unit Description Selection
pzd.D21 pzd.E*pzd.nu/((1+pzd.nu)*(1
-2*pzd.nu))
Pa Elasticity
matrix, 21
component
Domains 13,
5
pzd.D31 pzd.E*pzd.nu/((1+pzd.nu)*(1
-2*pzd.nu))
Pa Elasticity
matrix, 31
component
Domains 13,
5
pzd.D41 0 Pa Elasticity
matrix, 41
component
Domains 13,
5
pzd.D51 0 Pa Elasticity
matrix, 51
component
Domains 13,
5
pzd.D61 0 Pa Elasticity
matrix, 61
component
Domains 13,
5
pzd.D12 pzd.E*pzd.nu/((1+pzd.nu)*(1
-2*pzd.nu))
Pa Elasticity
matrix, 12
component
Domains 13,
5
pzd.D22 pzd.E*(1+pzd.nu/(1-
2*pzd.nu))/(1+pzd.nu)
Pa Elasticity
matrix, 22
component
Domains 13,
5
pzd.D32 pzd.E*pzd.nu/((1+pzd.nu)*(1
-2*pzd.nu))
Pa Elasticity
matrix, 32
component
Domains 13,
5
pzd.D42 0 Pa Elasticity
matrix, 42
component
Domains 13,
5
pzd.D52 0 Pa Elasticity
matrix, 52
component
Domains 13,
5
pzd.D62 0 Pa Elasticity
matrix, 62
component
Domains 13,
5
pzd.D13 pzd.E*pzd.nu/((1+pzd.nu)*(1
-2*pzd.nu))
Pa Elasticity
matrix, 13
component
Domains 13,
5
pzd.D23 pzd.E*pzd.nu/((1+pzd.nu)*(1
-2*pzd.nu))
Pa Elasticity
matrix, 23
component
Domains 13,
5
pzd.D33 pzd.E*(1+pzd.nu/(1- Pa Elasticity
matrix, 33
Domains 13,
81

Name Expression Unit Description Selection
2*pzd.nu))/(1+pzd.nu) component 5
pzd.D43 0 Pa Elasticity
matrix, 43
component
Domains 13,
5
pzd.D53 0 Pa Elasticity
matrix, 53
component
Domains 13,
5
pzd.D63 0 Pa Elasticity
matrix, 63
component
Domains 13,
5
pzd.D14 0 Pa Elasticity
matrix, 14
component
Domains 13,
5
pzd.D24 0 Pa Elasticity
matrix, 24
component
Domains 13,
5
pzd.D34 0 Pa Elasticity
matrix, 34
component
Domains 13,
5
pzd.D44 0.5*pzd.E/(1+pzd.nu) Pa Elasticity
matrix, 44
component
Domains 13,
5
pzd.D54 0 Pa Elasticity
matrix, 54
component
Domains 13,
5
pzd.D64 0 Pa Elasticity
matrix, 64
component
Domains 13,
5
pzd.D15 0 Pa Elasticity
matrix, 15
component
Domains 13,
5
pzd.D25 0 Pa Elasticity
matrix, 25
component
Domains 13,
5
pzd.D35 0 Pa Elasticity
matrix, 35
component
Domains 13,
5
pzd.D45 0 Pa Elasticity
matrix, 45
component
Domains 13,
5
82

Name Expression Unit Description Selection
pzd.D55 0.5*pzd.E/(1+pzd.nu) Pa Elasticity
matrix, 55
component
Domains 13,
5
pzd.D65 0 Pa Elasticity
matrix, 65
component
Domains 13,
5
pzd.D16 0 Pa Elasticity
matrix, 16
component
Domains 13,
5
pzd.D26 0 Pa Elasticity
matrix, 26
component
Domains 13,
5
pzd.D36 0 Pa Elasticity
matrix, 36
component
Domains 13,
5
pzd.D46 0 Pa Elasticity
matrix, 46
component
Domains 13,
5
pzd.D56 0 Pa Elasticity
matrix, 56
component
Domains 13,
5
pzd.D66 0.5*pzd.E/(1+pzd.nu) Pa Elasticity
matrix, 66
component
Domains 13,
5

Shape functions
Name Shape
function
Unit Description Shape
frame
Selection
u Lagrange
(Quadratic)
m Displacement field, X component Material Domains 1
3, 5
v Lagrange
(Quadratic)
m Displacement field, Y component Material Domains 1
3, 5
w Lagrange
(Quadratic)
m Displacement field, Z component Material Domains 1
3, 5

Weak expressions
Weak expression Integration frame Selection
-pzd.Sl11*test(pzd.el11)-2*pzd.Sl12*test(pzd.el12)-
2*pzd.Sl13*test(pzd.el13)-pzd.Sl22*test(pzd.el22)-
2*pzd.Sl23*test(pzd.el23)-pzd.Sl33*test(pzd.el33)
Material Domains 13, 5

83

1.4.6 Fixed Constraint 1

Fixed Constraint 1
Selection
Geometric entity level Boundary
Selection Boundaries 1, 4, 7

Equations

Settings
Settings
Description Value
Apply reaction terms on All physics (symmetric)
Use weak constraints Off

Shape functions
Constraint Constraint force Shape function Selection
-u test(-u) Lagrange (Quadratic) Boundaries 1, 4, 7
-v test(-v) Lagrange (Quadratic) Boundaries 1, 4, 7
-w test(-w) Lagrange (Quadratic) Boundaries 1, 4, 7

84

1.4.7 Electric Potential 1

Electric Potential 1
Selection
Geometric entity level Boundary
Selection Boundary 16

Equations

Settings
Settings
Description Value
Electric potential 20
Apply reaction terms on All physics (symmetric)
Use weak constraints Off

Variables
Name Expression Unit Description Selection
pzd.nD unX*(down(pzd.Dx)-
up(pzd.Dx))+unY*(down(pzd.Dy)-
up(pzd.Dy))+unZ*(down(pzd.Dz)-
up(pzd.Dz))
C/m^2 Surface charge
density
Boundary 16
pzd.V0 20 V Electric
potential
Boundary 16

85

Shape functions
Constraint Constraint force Shape function Selection
pzd.V0-V test(pzd.V0-V) Lagrange (Quadratic) Boundary 16

1.4.8 Ground 1

Ground 1
Selection
Geometric entity level Boundary
Selection Boundary 17

Equations

Settings
Settings
Description Value
Apply reaction terms on All physics (symmetric)
Use weak constraints Off

Variables
Name Expression Unit Description Selection
pzd.nD unX*(down(pzd.Dx)-
up(pzd.Dx))+unY*(down(pzd.Dy)-
up(pzd.Dy))+unZ*(down(pzd.Dz)-
up(pzd.Dz))
C/m^2 Surface charge
density
Boundary 17
86

Name Expression Unit Description Selection
pzd.V0 0 V Electric
potential
Boundary 17

Shape functions
Constraint Constraint force Shape function Selection
pzd.V0-V test(pzd.V0-V) Lagrange (Quadratic) Boundary 17

1.5 Mesh 1
Mesh statistics
Property Value
Minimum element quality 0.1305
Average element quality 0.6686
Tetrahedral elements 5586
Pyramid elements 0
Prism elements 0
Hexahedral elements 0
Triangular elements 1938
Quadrilateral elements 0
Edge elements 264
Vertex elements 24


Mesh 1
87

1.5.1 Size (size)
Settings
Name Value
Maximum element size 0.01
Minimum element size 0.0018
Curvature factor 0.6
Resolution of narrow regions 0.5
Maximum element growth rate 1.5

1.5.2 Free Tetrahedral 1 (ftet1)
Selection
Geometric entity level Remaining

88

2 Study 1
2.1 Stationary
Study settings
Property Value
Include geometric nonlinearity Off

Mesh selection
Geometry Mesh
Geometry 1 (geom1) mesh1

Physics selection
Physics Discretization
Piezoelectric Devices (pzd) physics

2.2 Solver Configurations
2.2.1 Solver 1
Compile Equations: Stationary (st1)
Study and step
Name Value
Use study Study 1
Use study step Stationary

Dependent Variables 1 (v1)
General
Name Value
Defined by study step Stationary

Initial values of variables solved for
Name Value
Solution Zero

Values of variables not solved for
Name Value
Solution Zero

89

Displacement field (Material) (comp1.u) (comp1_u)
General
Name Value
Field components {comp1.u, comp1.v, comp1.w}

Electric potential (comp1.V) (comp1_V)
General
Name Value
Field components comp1.V

Stationary Solver 1 (s1)
General
Name Value
Defined by study step Stationary

Log
Stationary Solver 1 in Solver 1 started at 26-Sep-2014 16:48:03.
Linear solver
Number of degrees of freedom solved for: 25806.
Symmetric matrices found.
Scales for dependent variables:
Displacement field (Material) (comp1.u): 4e-011
Electric potential (comp1.V): 1
Iter Damping Stepsize #Res #Jac #Sol LinErr LinRes
1 1.0000000 0.78 1 1 1 1.3e-009 2.3e-014
Stationary Solver 1 in Solver 1: Solution time: 5 s
Physical memory: 1.16 GB
Virtual memory: 1.27 GB
Fully Coupled 1 (fc1)
General
Name Value
Linear solver Direct 1

Direct 1 (d1)
General
Name Value
Solver PARDISO

90

3 Results
3.1 Data Sets
3.1.1 Solution 1
Solution
Name Value
Solution Solver 1
Component Save Point Geometry 1


Data set: Solution 1
91

3.2 Plot Groups
3.2.1 Displacement (pzd)

Surface: Total displacement (m)
3.2.2 Potential (pzd)

Slice: Electric potential (V)