! == Material parameters for piezo
! -- Elastic coefficients
! Allgemeine Definition eines isotroen Werkstoffes
! über Anpassung der Koeffizienten entsteht ein transversal isotroper
! Werkstoff (keine Zug-Scherungskopplung!)
! Alle Werte in N/m^2
c11      = 14.7e10    
c12      = 8.1e10   
c13      = 8.1e10
c16	 = 0
c14	 = 0
c15	 = 0
c22	 = c11
c23	 = c13
c26	 = 0
c24	 = 0
c25	 = 0 
c33      = 13.2e10
c36	 = 0
c34	 = 0
c35	 = 0
c44      = 3.1e10
c45	 = 0
c55	 = c44
c66	 = c44
c64	 = 0
c65	 = 0

! Nur Eingabe des oberen Dreiecks, da auch für anisotrope Materialien gillt:
! cij=cji. deswegen nur unteres oder oberes Dreieck notwendig
tb,ANEL,0	! ,1 bedeutet Eingabe als Nachgiebigkeitsmatrix, 0 ist Steifigkeitsmatrix
tbdata, 1,  c11, c12, c13,    c16, 	c14,    c15
tbdata, 7,       c22, c23,    c26,	c24,	c25
tbdata,12,            c33,    c36,      c34,    c35
tbdata,16,                    c66,      c64,    c65   
tbdata,19,                              c44,    c45   
tbdata,21,                                      c55

! -- Piezoelectric coefficients, Coulombs/m^2
e11 = 0
e21 = 0
e31 = -3.88*1e10
e12 = 0
e22 = 0
e32 = e31 
e13 = 0
e23 = 0
e33 = 0
e33 = 13.91*1e10
e16 = 0
e26 = 0
e36 = 0
e14 = 0
e24 = 10.34*1e10
e34 = 0
e15 = e24
e25 = 0
e35 = 0  
!	0 -- 
!	Piezoelectric stress matrix [e] (used as supplied)
!	1 -- 
!	Piezoelectric strain matrix [d] (converted to [e] form before use)

tb,PIEZ,0
tbdata, 1,  e11,    e21,   e31
tbdata, 4,  e12,    e22,   e32
tbdata, 7,  e13,    e23,   e33
tbdata,10,  e16,    e26,   e36     
tbdata,13,  e14,    e24,   e34     
tbdata,16,  e15,    e25,   e35


! alt für solid 226: dielectricity matrix in F/m
! für solid 5 eingabe als absolutwerte nötig
mp,PERx,1,812                  
mp,PERy,1,812 
mp,PERz,1,894


! -- Density, kg/m3 ist ungefähr so schwer wie Stahl
mp,DENS,1,7500