!Piezo-Biegebalken
!michel23, 09.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,
BLOCK,  0,    45E-3,    0, 7.2E-3,    0.53E-3,   0.78E-3,




!Zusammenkleben
VSEL,ALL
VGLUE,ALL

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

VSEL,,,,4
VATT,2,,2,0

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

VSEL,ALL



!Meshen
!Anzahl Elemente in der Länge
LESIZE,29, , ,30, , , , ,1 

!Anzahl Elemente in der Dicke
lsel,s,line,,9
lsel,a,line,,37
lsel,a,line,,41
lesize,all,,,1
allsel

vmesh,all





!Elektroden aussen:
asel,s,,,1	
asel,a,,,14				
nsla,s,1					
cp,1,volt,all					
*GET,aussen,NODE,0,NUM,MIN			
ALLSEL


!Elektroden innen:
asel,s,,,2
asel,a,,,8
nsla,s,1
cp,2,volt,all	
*GET,innen,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,22,ux,0
da,22,uy,0
da,22,uz,A

da,26,ux,0
da,26,uy,0
da,26,uz,A


!Feste Einspannung Alternative
!ASEL,S,,,22
!ASEL,A,,,26
!ASEL,A,,,5
!NSLA,S,0
!CP,3,UX,ALL
!CP,3,UY,ALL
!CP,3,UZ,ALL
!*GET,EINSPANNUNG,NODE,0,NUM,MIN
!D,EINSPANNUNG,UX,0
!D,EINSPANNUNG,UY,0
!D,EINSPANNUNG,UZ,A




!Wiederstand der Größe D zwischen Elektroden aussen und innen
/PREP7

D=1000

!get the number of solid elements
*get,elementcount,elem,,count 

R,1,D,
ET,4,CIRCU94,0,0
TYPE,4
E,aussen,innen



!Solution
/SOL

!spannung anlegen			

D,aussen,volt,0
!D,innen,volt,0


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


/POST26
NUMVAR,200


FINISH