Published July 31, 2024 | Version v2
Dataset Open

Randomised material parameter impedance dataset of piezoelectric rings (RaPIDring)

  • 1. ROR icon Paderborn University

Description

Version notes

This updated version uses larger ranges for the material parameters to account for a wider range of samples. Further, an erroneous boundary condition in the simulation model used to generate the data set in the previous version has been fixed.

Description of the dataset

This dataset contains the results of 301008 finite-element simulations of the complex, frequency-dependent electrical impedance of a piezoelectric ring with randomised material parameters. Each impedance consists of 2000 samples in the frequency domain up to 8 MHz. We assume the sample to be dielectric, thus the impedance a frequency 0 Hz is infinite. The piezoelectric ring has an outer radius of 6.35 mm an inner radius of 2.6 mm and a thickness of 1 mm. The transverse isotropic material parameters are sampled from independent uniform distributions with ranges that are intended to represent the behaviour of different piezoceramic materials. The parameters of the Rayleigh damping model (alpha_M and alpha_K) are sampled from a logarithmic distribution to account for the larger parameter range.

Parameter   Min Max Unit     Description
c11 110 165 GPa Elastic stiffness
c12 60   150 GPa Elastic stiffness
c13 60   100   GPa Elastic stiffness
c33 110   140 GPa Elastic stiffness
c44 18   30   GPa Elastic stiffness
eps11   3 12 nF/m Dielectric permittivity
eps33   4  8 nF/m Dielectric permittivity
e15 6 20 C/m^2 Piezoelectric coupling
e31 1 7.5 C/m^2 Piezoelectric coupling
e33 8 20 C/m^2 Piezoelectric coupling
alpha_M 0.2 150     1/ms Mass-proportional damping
alpha_K 10 1000 ps   Stiffness-proportional damping
density 7600 7850 kg/m^3 Density

Files

The dataset contains the following files:

  • dataset.hdf5: The main dataset file in HDF5 format. Refer to the next section on how to load the dataset.
  • simulation.xml: The file describing the parameters of the finite element simulation. This file can be used with openCFS along with the mesh and material file to simulate the behaviour of the piezoelectric ceramic in the frequency domain. This simulation file was used with randomised material parameters to generate the dataset.
  • material.xml: The material file used for the simulation with exemplary material parameters values.
  • ring.geo: The gmsh geometry file with the axisymmetric representation of the piezoelectric ring.
  • ring.msh: The mesh file generated with gmsh using the ring.geo file.

Loading the dataset

The dataset is stored as a HDF5 file, which can be opened with all libraries that support that format, e.g. in Python using the h5py library:

import h5py

# Open the dataset in read mode.
file = h5py.File("dataset.hdf5", "r")

# Impedances as a 301008 x 2000 array of complex numbers.
impedances = file["impedances"]
# Material parameter values as a 301008 x 13 array of real numbers.
parameters = file["parameters"]
# Frequency vector of the impedance with length 2000.
frequencies = file["meta"]["frequencies"]
# 13 strings with the identifiers of the material parameters.
parameter_labels = file["meta"]["parameter_labels"]

Simulating impedances

To generate a result for the electrical impedance using the supplied simulation files, download and install openCFS and call the executable with the simulation.xml, but omit the file extension, e.g.:

cfsbin.exe simulation

The path to the executable of openCFS will depend on your operating system and installation. Running the simulation will result in the creation of several files and folders. Among those files will be the result for the electric charge on one of the electrodes of the sample, which will be placed in the 'history' subfolder. We can determine the current by taking the time derivative of the charge and already know the voltage because we excited the piezoceramic with an electric potential of 1 V in the simulation. Because the simulation is conducted in the frequency regime, all we have to do is to divide voltage by current to get the frequency dependent electrical impedance. The loading of the result and calculation of the impedance is implemented in the following Python script as an example:

import numpy as np

# Load result file for electric charge
result_path = 'history/simulation-elecCharge-surfRegion-ground.hist'
data = np.loadtxt(result_path)

frequency = data[:, 0]
# Convert polar representation from file to complex numbers.
charge = data[:, 1] * np.exp(1j * 2 * np. pi / 360 * data[:, 2])

# Excitation potential is 1 V in simulation.
potential = 1
# Determine impedance by applying Z = V / I = V / (j omega Q).
impedance = potential / (1j * 2 * np.pi * frequency * charge)

Files

simulation.xml

Files (9.7 GB)

Name Size Download all
md5:0366527557b030789fd5458c8f20a891
9.7 GB Download
md5:6e195cc169e74b2b90a9c88da9f1cfab
1.5 kB Preview Download
md5:42289f4c8a2860b0d17f065e9689c92c
483 Bytes Download
md5:b19edddfe7fa9f961ea1d79b26863900
11.7 kB Download
md5:6180a6583fcb002bbc997a17e4f4475d
2.5 kB Preview Download