!Piezo-Biegebalken
!michel23, 15.04.2010


!/BATCH

/PREP7

/TITLE, Monomorpher Piezo-Biegebalken




!Elementtypen, Keyoption 1001 für Piezo
ET, 1,  SOLID226, 1001
ET, 2,  SOLID186
KEYOPT,1,1,1001 				
KEYOPT,1,2,0					
KEYOPT,1,4,0
	

!Anisotropic elastic Matrix				
TB,ANEL,1,1,0,1
TBTEMP,0
TBDATA,,14.2e-12,-3.7e-12,-6.5e-12,0,0,0  
TBDATA,,14.2e-12,-6.5e-12,0,0,0,20.6e-12   
TBDATA,,0,0,0,43.0e-12,0,0  
TBDATA,,43.0e-12,0,35.8e-12,,,  	

						
!Piezo matrix					
TB,PIEZ,1,,,1
TBDATA,,0,0,-315e-12,0,0,-315e-12
TBDATA,,0,0,640e-12,0,895e-12,0
TBDATA,,895e-12,0,0,0,0,0,0 


!Dielektrizitätszahlen (DPER  —  Anisotropic electric permittivity) 
TB,DPER,1,,,1   
TBDATA,,4750,4750,4500,0,0,0
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,DENS,1,,8100 


!zwischenmaterial -> hier: CFK
MPTEMP,,,,,,,,  			
MPTEMP,0  				
MPDATA,EX,2,,2.3e11 			
MPDATA,PRXY,2,,0.29      
MPTEMP,,,,,,,,  
MPTEMP,0  
MPDATA,DENS,2,,1800
MPTEMP,,,,,,,,
MPTEMP,0
MPDATA,DMPR,2,,0.02



!Geometrie erstellen
BLOCK,  0,    45E-3,    0, 7.2E-3,          0,   0.25E-3,
BLOCK,  0,    45E-3,    0, 7.2E-3,    0.25E-3,   0.53E-3,

VSEL,ALL
VGLUE,ALL

VSEL,,,,1
VATT,1,,1,0

VSEL,,,,3
VATT,2,,2,0




!Meshen

LESIZE,17, , ,20, , , , ,1 
lsel,s,line,,9
lsel,a,line,,21
lesize,all,,,1
allsel
vmesh,all
lsel,all




!Elektrode oben:
asel,s,,,1
nsla,s,1
cp,1,volt,all
*GET,oben,NODE,0,NUM,MIN
ALLSEL



!Elektrode unten:
asel,s,,,2
nsla,s,1
cp,2,volt,all
*GET,unten,NODE,0,NUM,MIN
ALLSEL



!Feste einspannung, Auslenkung A
A=5e-4

da,5,ux,0
da,5,uy,0
da,5,uz,A

da,16,ux,0
da,16,uy,0
da,16,uz,A



!Widerstand anschliessen
D=1000

!R,1,D,
!ET,4,CIRCU94,0,0
!TYPE,4
!E,oben,unten


/SOL

!spannung anlegen und erden

!D,oben,volt,200
D,unten,volt,0


ANTYPE,3
HARFRQ,288,292,
!HROUT,ON
NSUBST,1,   
KBC,1
SOLVE   


/POST26
NUMVAR,200


FINISH