fini
/clear
*abbr,AUFBAU2D ohne Fluid,/input,2D-Scheibe,txt
/prep7


et,1,Plane223,1001,,1		!Piezo
et,2,Plane82,,,1		!Kleber
et,3,Plane82,,,1		!Si






!----Morganparameter
!---Y-Richtung Polarisierung
MP,PERX,1,1850
MP,PERY,1,2250
MP,PERZ,1,1850
MP,DENS,1,8000



!----Glue---
mp,ex,2,2.7e9
mp,nuxy,2,0.25
mp,dens,2,1054


!	Piezoelektrische Matrix
!----------------------------------  Welche Parameter wählen???


GSV=1
!---Morganparameter 
d31=-225e-12*GSV
d33=545e-12*GSV
d15=515e-12*GSV


!-3D Fall
TB,PIEZ,1,,,1			!-> d-Matrix verwenden und Material 1
tbdata,2,d31				!Input Piezoelectric stress matrix [e]
tbdata,5,d33
tbdata,8,d31
tbdata,10,d15
tbdata,15,d15

1-2D Fall
!TBDATA,2,d31			! 2D Matrix wird auf 4 Werte reduziert
!TBDATA,5,d31		
!TBDATA,8,d33
!TBDATA,10,d15		

!---Werte von Morgan

s11=15.5e-12
s12=-4.96e-12
s13=-7.22e-12
s33=20.7e-12
s44=39e-12
s55=39e-12
s66=2*(s11-s12)


TB,ANEL,1,,,1
!TBDATA,1,s11,s12,s13
!TBDATA,7,,s33,s23
!TBDATA,13,,,s11
!TBDATA,19,,,,s66
TBDATA,1,s11,s12,s13	
TBDATA,7,s11,s13		!s12,s31,s44,s55,s66 aus Schaumburg S.204, PZT-5A
TBDATA,12,s33			!die Größe von s11, s33 sind etwa gleich, aber s33 aus PZT-5H Morgan Werte
TBDATA,16,s44			!PZT503 Morgan	
TBDATA,19,s55			!PZT503 Morgan	
TBDATA,21,s66


!----Silicon-----------------------------------------

mp,ex,3,129.198e9
mp,prxy,3,0.3			!Poisson coefficient
mp,dens,3,2329	


openingx=400e-6
lip=605e-6
endlip=655e-6
actorend=3.50e-3
end=3.65e-3			! X-RICHTUNGSPARAMETER	

openingz=40e-6		
gap=0
silicon=100e-6			! Z/Y_RICHTUNGSPARAMETER
kleber=110e-6

actorlayer=15e-6		! Dicke, nicht Absolutwert!! 




local,11
local,12,,,,,180
csys,11				!Koordinatensysteme
		
k,3,0,gap
k,4,0,silicon
k,5,0,kleber
k,6,0,kleber+1*actorlayer
k,7,0,kleber+2*actorlayer
k,8,0,kleber+3*actorlayer
k,9,0,kleber+4*actorlayer
k,10,0,kleber+5*actorlayer
k,11,0,kleber+6*actorlayer
k,12,0,kleber+7*actorlayer
k,13,0,kleber+8*actorlayer
k,14,0,kleber+9*actorlayer		!X=0; y/z = variabel

k,24,actorend,gap			
k,25,actorend,silicon			!actorend
k,26,actorend,kleber		
k,27,actorend,kleber+1*actorlayer
k,28,actorend,kleber+2*actorlayer
k,29,actorend,kleber+3*actorlayer
k,30,actorend,kleber+4*actorlayer
k,31,actorend,kleber+5*actorlayer
k,32,actorend,kleber+6*actorlayer
k,33,actorend,kleber+7*actorlayer
k,34,actorend,kleber+8*actorlayer
k,35,actorend,kleber+9*actorlayer

k,38,end,gap				!end	
k,39,end,silicon



a,3,24,25,4


a,24,38,39,25
a,4,25,26,5

a,5,26,27,6
a,6,27,28,7
a,7,28,29,8
a,8,29,30,9
a,9,30,31,10
a,10,31,32,11
a,11,32,33,12
a,12,33,34,13
a,13,34,35,14

allsel
aglue,all		!alle vergluen
allsel


x=1
lesize,4,,,10*x
lesize,10,,,1*x

allsel
lsel,s,,,13
lsel,a,,,16
lsel,a,,,19
lsel,a,,,22
lsel,a,,,25
lsel,a,,,28
lsel,a,,,31
lsel,a,,,34
lsel,a,,,37
lesize,all,,,1.5*x
allsel

lesize,1,,,350*x
lesize,5,,,15*x
allsel

asel,s,,,1
asel,a,,,2
type,3
mat,3
mshape,0,2D
mshkey,1
amesh,all
allsel			!Silizium	

asel,s,,,3
type,2
mat,2
mshape,0,2D
mshkey,1
amesh,all
allsel			!Klebeschicht

asel,a,,,5
asel,a,,,7
asel,a,,,9
asel,a,,,11
type,1 $ esys,12
mat,1			!Piezo


mshape,0,2D
mshkey,1
amesh,all

allsel

asel,s,,,4
asel,a,,,6
asel,a,,,8
asel,a,,,10
asel,a,,,12
type,1 $ esys,11		!normal (GND)
mat,1
	

mshape,0,2D
mshkey,1
amesh,all
allsel



!----------------------------Randbedingungen----------------------

allsel
nsel,s,loc,x,end
D,all,Ux,0
D,all,Uy,0
allsel				!Fest eingespannt am Rand

nsel,s,loc,x,0
D,all,Ux,0
!Option des Symm.			!Achsensymmetrisch
!dl,all,,symm		
allsel



!----------------------------Spannungen---------------------------

spannung=-15


lsel,s,,,36
lsel,a,,,30
lsel,a,,,24
lsel,a,,,18
lsel,a,,,12
nsll,s,1
d,all,volt,spannung
allsel

lsel,s,,,33
lsel,a,,,27
lsel,a,,,21
lsel,a,,,15
lsel,a,,,9
nsll,s,1
d,all,volt,0
allsel

/solu
solve