Piezoelectric Unit Cube¶
All the files in this tutorial can be
This example illustrates the linear piezoelectric material, by considering a simple unit cube made from piezoelectric material.
The faces of the cube are called like the compass, N (north), E (east), S (south), W (west) for the +y, +x, -y, -x normal faces, and T (top) and B (bottom) for the +z and -z normal faces, respectively.
The lines are named by the adjoining faces, e.g.
L_NW for the line created by the intersection of north-surface
The same is done for corner points, e.g.
We set statically determined displacement BCs on three nodes, and apply an electric field in negative z-direction.
Sketch of the domain +------+. |`. | `. | `+--+---+ | | | | +---+--+ | `. | `. | `+------+
Mesing and CFS computation can be done by running
- Type 'cubit' on the terminal.
- On the Cubit GUI open up the journal file
UnitCube.jouand run it.
- A geometry of an unit cube is created.
- Export the geometry as an ANSYS-cdb mesh file
- You can also save the created geometry as
UnitCube.cub5to open it directly on Cubit (optional).
UnitCube.cdb was created this way.
Simulation with CFS¶
To start the computation run the following command in the terminal
cfs -p piezo_static.xml job
job can be any name you choose for the simulation.
CFS will write some output on the terminal, and produce the two files
job.info.xml, which contains some details about the run, and
job.cfs in the
results_hdf5 directory, which you can view with ParaView.
Open ParaView and choose File -> Load State ... to load the provided visualization pipeline
Is the deformation expected? You could also compute the results analytically from the linearized piezoelectric material law.
Additionally, you can modify the example to answer the following questions:
- Switch to fully constraint BCs and compute the stress state. Is the result expected?
- What happens when you change the orientation of the polarization?