In this tutorial we will learn how to simulate a 2D Film Bulk Acoustic Resonator (FBAR). We will analyse the impedance and mode shapes of the device.
You will learn:
- The basic simulation workflow in OnScale Designer
- How to set up a 2D model
- How to create a simple geometry
- How to edit material properties
- How to add circuits into a design
- How to display and post-process results
What is an FBAR?: FBAR is an acoustically isolated device consisting of a piezoelectric material sandwich between two electrodes. FBAR devices using piezoelectric films with thicknesses ranging from several micrometres down to tenth of micrometers resonate in the frequency range of roughly 100 MHz to 10 GHz. Aluminium nitride and zinc oxide are two common piezoelectric materials used in FBARs. Common applications of FBARs is radio frequency (RF) filters used in phones and other wireless applications. Filters are made from a network of resonators and are designed to remove unwanted frequencies from being transmitted whilst allowing specific frequencies to be received. These networks are usually in half-ladder, full-ladder, lattice or stacked topologies. FBARs can also be used in sensor applications in this case when a device is subject to some form of mechanical pressure its resonance frequency will shift.
Model Definition
Characteristics of the model
Model: |
FBAR |
Mesh Size: |
20 elements per wavelength |
Analysis Time: |
10 us |
Output Results: |
- Y Displacement Mode Shape at 1.967 MHz |
Material Data
Name |
Code Name |
Density |
Bulk Velocity |
Shear Velocity |
Poling |
Aluminium Nitride |
aln |
3230 kg.m-3 |
- |
- |
Y+ |
Molybdenum |
moly |
10220 kg.m-3 |
6649.79 ms-1 |
3509.29 ms-1 |
- |
Silicon, Generic |
si |
2330 kg.m-3 |
7526 ms-1 |
4346 ms-1 |
- |
Note: Two of these materials were taken from the Project Materials Database. Molybdedum is created by the user as it is not in the Database.
Why This Simulation?
It is useful to use simulation when designing FBARs because the thickness of the electrodes and piezoelectric layer determine the frequency range in which the device can operate. It allows designers to optimise their design for frequency range, quality factor and many other key performance characteristics.
This 2D FBAR model provides a simple starting structure consisting of a piezoelectric active layer (Aluminum Nitride), Molybdenum electrodes and silicon substrate material at the sides of the air cavity. In this simulation, we apply a voltage across the AlN layer and analyse the resulting impedance and mode shapes.
The Simulation Process
Let's go through the step by step tutorial and see how to simulate a 2D FBAR in OnScale!
Step 1 - Create a New Project
Open up OnScale in Designer Mode. The first step is to create a new project.
- In the Home tab of the ribbon, click New Project. The New Project window shows.
- Type a name for the project.
- 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 Model.
- Select the Advanced checkbox.
- For Distance, select μm.
- Click OK.
Step 2 - Set the Frequency of Interest
First we will specify the frequency of interest.
- Select Project Settings in the Model Tree
- Expand the Frequency of Interest property
- Tick the box
- Enter a value of 1.9e9
- Expand the Frequency of Damping property
- Tick the box
- Enter a value of 1.9e9
Step 3 - Add the Materials from the Material DB
The second step is to add the required materials for the FBAR transducer to the Project Materials database.
Molybdenum
Add electrode material
- Click + icon, next to Materials to open the Material Database
- Select Add New Material > Add Project Material
- Give the Material Description - Molybdenum
- Set Material Name to moly
- Select Material Category to be Metal
- Toggle the Damping property
- Set Density to 10220
- Set Bulk Velocity to 6649.79
- Set Shear Velocity to 3509.29
- Set Damping Units to dB/MHz/cm
- Set Bulk Attenuation to 0.1
- Set Shear Attenuation to 0.3
- Set Bulk Power Law to 1
- Set Shear Power Law to 1
- Select Save
Aln
Add piezoelectric material.
- Expand the 'Piezoelectric' materials dropdown
- Double click on 'Aluminium Nitride - aln' to add it to Project Materials
Silicon
Now add the substrate material.
- Expand 'Misc' materials dropdown
- Double click on 'Silicon, generic- si' to add it to Project Materials
Molybdenum 2
Copy electrode material
- Right click 'moly' in Project Materials and select Copy
- Rename material to moly2
- Select OK
- Select Done
Create Basic Geometric Shapes
We will make use of the geometric primitive shapes available in OnScale to build the FBAR.
Substrate
We will start by creating the substrate of the device.
- Select Rectangle button
- Select primitive_1
- Set the Material property of primitive_1 to si
- Set the X End property to 230
- Set the Y End property to 2
- Right click in the workspace and select Reset View to snap the view to the geometry
Cavity
Next we will add another rectangle to create the cavity.
- Right click primitive_1 and select Duplicate Selection
- Set the Material property of primitive_2 to void
- Set X Begin property of to 10
- Set X End property of to 220
Bottom Electrode
Next we will add another rectangle to represent the bottom electrode.
- Right click primitive_1 and select Duplicate Selection
- Set the Material property of primitive_3 to moly
- Set Y Begin property of primitive_3 to 2
- Set Y End property of primitive_3 to 2.4
Piezoelectric Layer
The piezoelectric layer is made from AlN in this example.
- Right click primitive_3 and select Duplicate Selection
- Set the Material property of primitive_4 to aln
- Set Y Begin property of primitive_3 to 2.4
- Set Y End property of primitive_3 to 3.2
Top Electrode
Lastly, we will add the top electrode.
- Right click primitive_4 and select Duplicate Selection
- Set the Material property of primitive_5 to moly2
- Set X Begin property to 15
- Set Y Begin property to to 3.2
- Set X End property to to 215
- Set Y End property to to 3.6
- Right click primitive_5 and select Duplicate Selection
- Set Y Begin property of primitive_6 to 3.6
- Set X End property to to 17
- Set Y End property to to 3.8
- Right click primitive_6 and select Duplicate Selection
- Set the Material property of primitive_7 to moly
- Set X Begin property of primitive_7 to 213
- Set X End property to to 215
Step 5 - Define a Time Function
The next step is to add a drive function. We will use a Ricker Wavelet time function with a frequency of 1.9 GHz.
- Click '+' to open the Drive Function dialogue
- Change function type from Sinusoidal to Ricker Wavelet
- Set Drive Frequency to 1.9e9
- Click Insert to close the window. A record called timefunc_1 will be added to the Model Tree
Step 6 - Define a Circuit
It is now time to define a circuit.
- Click '+' to open the Circuit dialogue
- Select the line between points 2 and 3
- Set Element to Resistor
- Set Resistance to 0.001
- Select Insert to close window. A record called circuit_1 will be added to the Model Tree
Step 7 - Choose The Right Mesh Size
It is time to set up the meshing of the model.
- Expand Mesh in the Model Tree
- Select Configuration
- Set Definitions to Advanced
- Set Elements per Wavelength to 20
- Set Mesh Velocity to Defined
- Set Value to 6000
Step 8 - Create Voltage Loads Between Electrodes and Piezo
We will now create loads on each side of the piezoelectric layer.
Create Load 1
- In the Model Tree, expand Boundary Conditions and then, beside Loads, click +.
- For Creation Mode, select Geometry Interface.
- For Geometry, select primitive_3 (moly) or click it in the model.
- For Interfacing Item, select primitive_4 (aln).
- For Load Type, select Voltage.
- For Area Scaling, type 0.0002.
- For Termination, select Ground.
- Click Create Load.
Create Load 2
The Load Definition window should still be open.
- For Geometry, select primitive_4 (aln) or click it in the model.
- For Interfacing Item, select primitive_5 (moly2).
- For Circuit, select circuit_1.
- For Termination, select timefunc_1.
- Click Create Load.
Step 9 - Define the Boundary Conditions
We need to change the boundary conditions so that the tail of the transducer is fixed and the water load has absorbing boundaries.
- Click Domain Boundaries in the Model Tree
- Set X Minimum to Fixed
- Set X Maximum to Fixed
- Set Y Minimum to Fixed
- Set Y Maximum to Free
Step 10 - Define the Type of Analysis
We will now set the model simulation time to be 10 us seconds
- Click Analysis
- Change Simulation Run Time to 1e-5
Step 11 - Define the Simulation Output Results
We will now define an output so we are able to see the displacement mode shape of the device at a 1.967 GHz.
- Click '+' to create a new output
- Set Output Type to Shape Data
- Set Array Type to Displacement
- Set Array Component to Y
- Set Frequency to 1.967e9
Step 12 - Run the Simulation 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
- Select Estimate
- Select Run
How to Get the Simulation Results?
Once the simulation has finished, the results are available in your storage to download.
- Click the Storage button to open the cloud storage
- Select your job from the dropdown menu
- Expand the simulation folder
- CTRL + select the -shape.flxdato & flxhst file, right click and select Download Selection
Choose an appropriate save location and close the cloud storage.
Step 13 - Check the Simulation Results
Switch to the Post Processor
- Click this icon to access the Post Processor to analyse simulation results
Open Results
- Click File Explorer
- Expand the job simulation folder
- Double click the data out file to open the mode shape results
- Double click the history file to open the time history results
Plot Impedance
- From Results Manager, select 'pize load2: Voltage'
- Select Impedance button
- Double click on 'Impd:load2.amp'
- Toggle Log yAxis
- Set xAxis Minimum to 2.1e9
- Set xAxis Maximum to 2.25e9
Plot Mode Shape at 1.967 GHz
- Expand Mode Shapes
- Expand Mode-1
- Select ydsp
- Right click and select Plot Shape Movie
- Select Continue
- Set Scale Factor to 0.001
- Select Play button to play mode shape animation
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 of your own.