In this tutorial, you will learn how to set up a simple 2D transducer model with a backing and matching layer and extract outputs such as impedance, directivity plots and transmit sensitivities.
You will learn:
- The Basic Simulation Workflow in OnScale Designer
- How to set up a 2D axi-symmetric model
- How to create a simple geometry
- How to simulate piezoelectric materials
- How to display and post-process your results
- How to use Extrapolation tool
Model Definition
Characteristics of the model
Model: |
PZT Disc dimensions: 10 mm x 2 mm Backing dimensions: 10 mm x 5 mm Matching Layer dimensions: 10 mm x 1/4 Wavelength mm Water of dimensions: 15mm x 20mm |
Mesh Size: |
0.1 mm |
Output Results: |
- Impedance - Directivity at 1 MHz - Transmit Sensitivity |
Material Data
Name | Water at 25C | CTS 3203HD | Vantico HY1300/CY1301 | Backing (20% VF) |
Code Name | watr | pmt3 | hard | back20 |
Density | 1000 kg.m-3 | 7820 kg.m-3 | 1149 kg.m-3 | 4800 kg.m-3 |
Bulk Velocity | 1496 ms-1 | 4708.36 ms-1 | 2536 ms-1 | 1800 ms-1 |
Shear Velocity | 0 ms-1 | 1687.891 ms-1 | 1179 ms-1 | 962 ms-1 |
Note: Material Data in OnScale are generally defined using the bulk velocity and the shear velocity parameters instead of the more traditional Elastic Modulus and Poisson's Ratio. You can check this page if you want to understand the relation between those parameters.
Why This Simulation?
Transducer sensors often are constructed using multiple layers to try maximize key performance indicators (KPIs) such as bandwidth. Common ways are to introduce a heavy backing layer and matching layers.
The full design is often tested under specific environmental conditions to measure real-world performance to validate the design.
Important: Make Sure that you select the correct Axis of Symmetry from the beginning in the New project window because it is impossible to change afterwards in the designer mode.
The Simulation Process
Let's go through the step by step tutorial and see how to simulate this in OnScale!
Step 1 - Create a New Project
- Click New Project to open the New Project window.
- Name the project Transducer_2D.
- If desired, change the save location and/or project file name by clicking … beside Project File.
- For Analysis, select Mechanical Dynamic.
-
For Model Type, select 2D Axi-Symmetric Model (Y).
Warning: You can't change the axis of symmetry later, so be sure to select Y here and not X!
- Select the Advanced checkbox.
- For Distance, select mm.
- Click OK.
Step 2 - Add the Materials from the Material DB
First we will add the materials needed from the material database. We will use pmt3 and water in this tutorial.
- Click Project Materials to open the database
- Expand Fluid and add watr to the Project Materials (double click)
- Expand Epoxy and add hard and back20
- Expand Piezoelectric and add pmt3
- Expand pmt3 and change the poling direction to Y+
- Click Done
Step 3 - Create Basic Geometry Shapes
Note: After making changes to X and Y right click the workspace and select Reset View
- Click Rectangle
- Change Material to watr
- Set End (mm): X (mm) = 15
- Set End (mm): Y (mm) = 15
- Right click primitive_1 and select Duplicate Selection
- Click primitive_2
- Change material from watr to back20
- Set End (mm): X(mm) = 10
- Set End (mm): Y(mm) = 5
- Right click primitive_2 and select Duplicate Selection
- Click primitive_3
- Change material from back20 to pmt3
- Set Begin (mm): Y(mm) = 5
- Set End (mm): Y(mm) = 7
- Right click primitive_3 and select Duplicate Selection
- Click primitive_4
- Change material from pmt3 to hard
- Set Begin (mm): Y(mm) = 7
- Set End (mm): Y(mm) = 7.634
Step 4 - Define a Time Function
We will now add a Ricker Wavelet drive function for later use as out loads require a time function be set.
- Click '+' to open the Define Input Time function window
- Change to Ricker Wavelet
- Click Insert to close the window. A record called timefunc_1 will be added to the window
Step 5 - Choose the right Mesh Size
We will now change the mesh settings to use 15 elements per wavelength
- Select Configuration
- Set Definitions to Wavelength Based
- Set Elements Per Wavelength to 15
Step 6 - Create the two electrodes
We'll create two loads to represent the two electrodes. We'll drive the top electrode and ground the bottom electrode.
Create Electrode 1
- In the Model Tree, expand Boundary Conditions and then, beside Loads, click +.
- For Creation Mode, select Geometry Interface.
- For Geometry, select primitive_3 (pmt3) or click it in the model.
- For Interfacing Item, select primitive_4 (hard).
- For Load Type, select Voltage.
- For Termination, select timefunc_1.
- Click Create Load.
Create Electrode 2
The Load Definition window should still be open.
- For Geometry, select primitive_3 (pmt3) or click it in the model.
- For Interfacing Item, select primitive_2 (back20).
- For Termination, select Ground.
- Click Create Load.
Step 7 - Define the Boundary Conditions
We will need to change the X minimum boundary condition to Symmetry as this model is symmetrical across that axis.
- Click Domain Boundaries
- Change the X Minimum boundary condition to Symmetry. All others will be absorbing
- Change remaining boundary conditions to Absorbing
Step 8 - Define the Analysis Simulation Time
We will now set the model simulation time to be 5e-5 seconds
- Click Analysis
- Change Simulation Run Time (s) to 2e-05
Step 9 - Define the Output Results
We will now define 1 output, a extrapolation boundary to allow us to generate outputs using the Extrapolation tool
Output Result 1: Extrapolation Data
- Click '+' this will create a new output
- Change Output Type to Extrapolation Data
- Change Cutting Plane to Y
- Change Cutting Point (mm) to 8.5
Step 10 - Run on the Cloud
At this point the model is completely set up and it can now be run on the cloud.
- Click Run on Cloud
- The option to rename your job. This is how it will appear in the storage
- Click Estimate
- Click Run
How to Get the Simulation Results?
The simulation results will need to be downloaded from the cloud storage in order to analyse the results in the post processor. More experience users may also be able to process Time Histories in Review.
- Click Storage this opens the window shown above
- Locate the job
- Click Download
- Click Download all
Choose an appropriate save location when the file explorer pops up and click Select Folder to close the window.
Step 11 - Check the Simulation Results
Switch to the Post Processor
- Click this icon to access the Post Processor
Open Results, Calculate and Plot Conductance
- In the File Explorer, Locate the Flex History File (.flxhst) & double Left Click
- Select the Charge record on load1
- Select Real/Imag from the drop down menu
- Click Admittance
- Two new records will appear, double click on Admt:load1.re to plot
- The real part of Admittance (Conductance) will be plotted allowing you to analyze the transducer performance in greater detail.
Using the Extrapolation Tool
- In the Tools tab, Select Extrapolation Tool
- Click Open
- Select the downloaded extrapolation file (.flxext)
- Click Open
- You will be prompted to open the flxmdl file. Click Yes
Directivity (Radial Plot)
- Select Radial Plot
- Change Y1 origin of directivity to 7.634e-3 m
- Change Radius (m) to 1
- Change Angular Spread to 180
- Change No. of points on line to 181
- Change Propagation Vector to Y
- Change Drive Frequency (Hz) to 1e6
- Click Start to calculate
Radial Plot Result
Transmit Sensitivity (TVR)
- Select TVR Calculation
- Select User Defined
- Alter Y reference value to 7.634e-3 m
- Change Propagation direction to Y
- Click Start to calculate
TVR Response
Try For Yourself
Now that we have introduced you to the tutorial try have a play around with some of the settings, add some other outputs, or use this model as a starting point for your own.