/PREP7 

		!CS
local,11							
local,12, , , , , ,180				!new CS			
csys,11



BLOCK,0,440,-40/2,40/2,0,2,				!beam
BLOCK,30,80,-15,15,2,2.5,				!Piezoelement
BLOCK,110,160,-15,15,2,2.5,
BLOCK,30,80,-15,15,2,-0.5,
BLOCK,110,160,-15,15,2,-0.5,
VGLUE,ALL               				!glue   all V 


!ET,1,SOLID45   			
!ET,2,SOLID5,3 			
				


				
ET,1,SOLID186   			! beam 
ET,2,SOLID226,1001
 			!piezo 

MPTEMP,,,,,,,,  			!m a t e r i j al  1
MPTEMP,1,0  
MPDATA,EX,1,,2.1E5  					! E
MPDATA,PRXY,1,,0.3				! Poisson koeficient
MPDATA,DENS,1,,0.785E-9					

				!p i e z o m a t e r i j a l   PIC 255 

MPDATA,DENS,2,,0.78E-9 					
TB,PIEZ,2 			! Piezoelectric stress matrix  [e]					 
TBDATA,3,-7.07E-6
TBDATA,6,-7.07E-6
TBDATA,9,13.82E-6
TBDATA,14,11.91E-6
TBDATA,16,11.91E-6 

MP,PERX,2,0.929		!dielectric matrix (permittivity) in x and y direction
MP,PERY,2,0.929
MP,PERZ,2,0.857 		!permittivity in z direction

TB,ANEL,2		!stiffness matrix [c]
TBDATA,1,12.29E4,7.66E4,7.02E4
TBDATA,7,12.29E4,7.02E4
TBDATA,12,9.71E4
TBDATA,16,2.32E4
TBDATA,19,2.23E4
TBDATA,21,2.23E4






				
!ET,1,SOLID186   			!beam 
!ET,2,SOLID226,1001			!piezo 






/com  -- Material matrices for PZT4 (polar axis along Y-axis): ANSYS input
/com
/com   [c11 c13 c12  0  0  0 ]        [ 0    0  e31 ]     
/com   [c13 c11 c12  0  0  0 ]        [ 0    0  e31 ]      
/com   [c12 c13 c33  0  0  0 ]        [ 0    0  e33 ]      
/com   [ 0   0   0  c44 0  0 ]        [ 0    0   0  ]
/com   [ 0   0   0   0 c55 0 ]        [ 0   e52  0  ]
/com   [ 0   0   0   0  0 c66]        [ e52  0   0  ]

!c11=12.29e4
!c13=7.66e4
!c12=7.02e4
!c33=9.71e4
!c44=2.32e4
!c55=2.23e4
!c66=2.23e4

!tb,anel,1				!Define structural table
!tbdata,1,c11,c13,c12			!Steifigkeitsmatrix [c]
!tbdata,7,c33,c13			
!tbdata,12,c11				
!tbdata,16,c44				
!tbdata,19,c44			
!tbdata,21,c66			

					!Piezoelektrische Spannungsmatrix [e]
!e31=-7.07E-6
!e33=13.82E-6
!e52=11.91E-6

!tb,piez,1				!Define Piez. table
!tbdata,3,e13				
!tbdata,6,e33
!tbdata,9,e13
!tbdata,14,e52
!tbdata,16,e52




!mp,perx,1,0.929				!Permittivity (x direction)
!mp,pery,1,0.929				!Permittivity (y direction)
!mp,perz,1,0.857				!Permittivity (z direction)
!mp,dens,1,7.8e-9			!Density





/com
/com  -- Beammaterial
/com
!mp,ex,2,2.1E5				!Elastizitätsmodul
!mp,prxy,2,0.3				!Poisson Koeffizient
!mp,dens,2,7.85e-9			!Densität




				!atribut
VSEL,S, , ,1 
VATT,1, ,1, ,

VSEL,S, , ,2,3,1 
VATT,2, ,2, ,

VSEL,S, , ,4,5,1
VATT,2, ,2,12,


VSEL,ALL
			
				!MESHING


SMRT,OFF
ESIZE,5

MSHKEY,1
MSHAPE,0,3D

VSEL,S, , ,1
Vmesh,1
VSEL,S, , ,2,3,1 
VSWEEP,ALL
VSEL,S, , ,4,5,1
CSYS,12
VSWEEP,ALL

CSYS,11


	
				!!Einspannung

FLST,2,1,5,ORDE,1   
FITEM,2,5   
!*  
/GO 
DA,P51X,ALL,0   
   



				
					!!1.Piezoelement_Oben	(x volt)

NSEL,S,LOC,Z,2.5
NSEL,R,LOC,X,30,80   
CP,1,VOLT,ALL   
*GET,knoten_1_oben,NODE, ,NUM,MIN   
D,knoten_1_oben,VOLT,4		



					!!1.Piezoelement_unten (0 Volt)
NSEL,S,LOC,Z,2
NSEL,R,LOC,X,30,80
NSEL,R,LOC,Y,-15,15
CP,2,VOLT,ALL
*GET,knoten_1_unten,NODE, ,NUM,MIN  	 
D,knoten_1_unten,VOLT,0 					


					!!2.Piezoelement_Oben	(x volt)

NSEL,S,LOC,Z,2.5 
NSEL,R,LOC,X,110,160   
CP,3,VOLT,ALL   
*GET,knoten_2_oben,NODE, ,NUM,MIN  	
D,knoten_2_oben,VOLT,4			
		
					!!2.Piezoelement_unten	(0 volt)

NSEL,S,LOC,Z,2
NSEL,R,LOC,X,110,160
NSEL,R,LOC,Y,-15,15
CP,4,VOLT,ALL
*GET,knoten_2_unten,NODE, ,NUM,MIN  	 
D,knoten_2_unten,VOLT,0 				

	
					!!3.Piezoelement_Oben	(x volt)
NSEL,S,LOC,Z,-0.5   
NSEL,R,LOC,X,30,80 
CP,5,VOLT,ALL  
*GET,knoten_3_oben,NODE,0,NUM,MIN  
D,knoten_3_oben,VOLT,4
	

					!!3.Piezoelement_unten	(0 volt)
NSEL,S,LOC,Z,0
NSEL,R,LOC,X,30,80
NSEL,R,LOC,Y,-15,15
CP,6,VOLT,ALL
*GET,knoten_3_unten,NODE, ,NUM,MIN  	  
D,knoten_3_unten,VOLT,0 
	
					!!4.Piezoelement_Oben	(x volt)
NSEL,S,LOC,Z,-0.5
NSEL,R,LOC,X,110,160 
CP,7,VOLT,ALL    
*GET,knoten_4_oben,NODE,0,NUM,MIN  
D,knoten_4_oben,VOLT,4 

					!!4.Piezoelement_unten	(0 volt)
NSEL,S,LOC,Z,0
NSEL,R,LOC,X,110,160
NSEL,R,LOC,Y,-15,15
CP,8,VOLT,ALL
*GET,knoten_4_unten,NODE, ,NUM,MIN  	
D,knoten_4_unten,VOLT,0 

allsel,all
	



					!!Kraftknoten

FLST,2,1,1,ORDE,1   
FITEM,2,32
!*  
/GO 
F,P51X,FZ,4 
!*  


finish
/SOL
ANTYPE,2
!*  
!*  
MODOPT,LANB,10  
EQSLV,SPAR  
MXPAND,10, , ,0 
LUMPM,0 
PSTRES,0
!*  
MODOPT,LANB,10,0,1600, ,OFF 

Solve

/POST1  
SET,LIST