/batch


voben = 2				!Voltage applied to the top of the PZT layer              	
vunten = 0				!Voltage applied to the bottom of the PZT layer

seltol,1e-6				!Selection tolerance

/prep7

					!Element types


et,1,solid98   				!Hexaedral Coupled-Field Solid

et,2,solid95  				!3-D 20-Node Structural 






/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,1.2.1E5				!Elastizitätsmodul
mp,prxy,2,0.3				!Poisson Koeffizient
mp,dens,2,7.85e-9			!Densität




/com
/com  -- Modeling
/com
BLOCK,0,440,0,40,0,2, 			!Balkenmodell

bLOCK,30,80,5,35,2,2.5, 		!Piezoelementmodell
btol,3.5e-10
vglue,all




					!Teile selektieren


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


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

VSEL,all


  					!Meshen
ESIZE,2,0,  
MSHKEY,0
MSHAPE,1,3d 

VSEL, , , ,       3 
VMESH,3

VSEL, , , ,       2   
VMESH,2









					!EInspannung
NSEL,s,loc,x,0
d,all,uz,0,,,,ux,uy


					!Kraft


Nsel,s,loc,x,440
Nsel,r,loc,y,20
Nsel,r,loc,z,2
*get,Kraftknoten,Node,,Num,min
F,Kraftknoten,Fz,4



					!Piezospannung


NSEL,S,LOC,Z,2.5
NSEL,R,LOC,X,30,80
NSEL,R,LOC,Y,5,35
CP,1,volt,all
*get,piezoknoten1,NODE,,NUM,MIN
d,piezoknoten1,volt,voben	
d,all,UX,0


NSEL,S,LOC,Z,2
NSEL,R,LOC,X,30,80
CP,2,volt,all
*get,piezoknoten2,NODE,,NUM,MIN
d,piezoknoten2,volt,vunten
d,all,UX,0


FINISH  

					!Modalanalyse

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

Solve