Numerical modeling of single-layered piezoelectric elements

Abstract

The focus of this report is on numerical simulation of piezoelectric materials. One example of such a material is quartz. By applying a pressure load to a piezoelectric material, electricity is generated – the so-called direct piezoelectric effect. In the opposite case, by applying an electric field, the material deforms (mechanically) – the so-called inverse piezoelectric effect. A mathematical model for linear piezoelectricity is described. This includes a model for linear elasticity, as well as a model for linear electrostatics. Piezoelectricity is expressed mathematically by combining the two above mentioned models. The explicit coupling between the elasticity problem and the electrostatics problem is through (modified) constitutive laws. Two simulators are implemented; one in Diffpack and one in MSC.Marc. In this way we are able to compare numerical results, and then verify the implementation in both software packages. In addition, this increases the general understanding of the problem. Simulation in Diffpack requires some low-level programming, but gives large possibilities for adjustments and assessments. Moreover, it provides access to all programming details needed. MSC.Marc is on the other hand to a larger extent an application software package. An advantage of MSC.Marc is that one is able to use robust software for advanced problems relatively fast. A main disadvantage is that most of the implementation details are hidden, which gives, among other things, uncertainty regarding the choice of solution method, what presumptions are made for the calculations, and the interpretation of quantities that are returned (results). For such a complex and comprehensive software package as MSC.Marc, detailed documentation, that easily can be searched, is a major challenge to work out. Two test examples are simulated in both software tools for different piezoelectric materials. In the first case we simulate the piezoelectric effect, whereas in the second case we consider the inverse piezoelectric effect. Comparison of the numerical results shows very good agreement in the two software packages for the material model cases included in this report.

Denne rapporten tar for seg numerisk modellering av piezoelektriske materialer. Et eksempel på et slikt materiale er kvarts. Piezoelektriske materialer har den egenskapen at ved å utsette disse for en trykklast, genereres elektrisitet – såkalt direkte piezoelektrisk effekt. I motsatt tilfelle, ved å påføre materialet et elektrisk felt, deformeres det (mekanisk) – såkalt invers piezoelektrisk effekt. En matematisk modell for lineær piezoelektrisitet er beskrevet. Dette inkluderer en modell for lineær elastisitet, så vel som en modell for lineær elektrostatikk. Piezoelektrisitet uttrykkes matematisk ved en kombinasjon av de to ovennevnte modellene. Den eksplisitte koblingen mellom elastisitetsproblemet og elektrostatikkproblemet er gjennom (modifiserte) konstitutive lover. To simulatorer er implementert; en i Diffpack og en i MSC.Marc. Dette gjør oss i stand til å sammenlikne numeriske resultater, og på den måten verifisere implementasjonen i begge programvarepakkene. I tillegg gir dette økt generell forståelse av problemet. Simulering i Diffpack krever en del lavnivå programmering, men gir store muligheter for egne tilpasninger og vurderinger. I tillegg har vi fullt innsyn i alle programmeringsdetaljer vi har behov for. MSC.Marc er på den andre siden i større grad en applikasjonsprogramvarepakke. En fordel med MSC.Marc er at man kan benytte robust programvare til avanserte problemer relativt raskt. En betydelig ulempe er at mange av implementeringsdetaljene er skjult, noe som blant annet gir usikkerhet knyttet til løsningsmetode, hvilke forutsetninger som er gjort i beregningene og hvilke størrelser som returneres (resultater). For en så komplisert og omfattende programvarepakke som MSC.Marc, er for øvrig detaljert nok dokumentasjon, som det samtidig er mulig å finne frem i, en stor utfordring å utarbeide. To testeksempler er simulert i begge programvarepakker for forskjellige piezoelektriske materialer. I det ene testproblemet simulerer vi piezoelektrisk effekt, mens vi i det andre testproblemet tar for oss invers piezoelektrisk effekt. Det er meget godt samsvar mellom de numeriske resultatene i de to programvarepakkene for de materialmodellene som er vist i denne rapporten.