/prep7

vtop = 1000				!Voltage applied to the top of the PZT layer              	
vbot = 0

/com
/com  -- Element types
/com 
et,1,solid5,3 				!Hexaedral Coupled-Field Solid; UX, UY, UZ, VOLT
et,2,solid95   				!3-D 20-Node Structural 


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

c11 = 107.6e9
c12 = 63.1e9
c13 = 63.8e9
c33 = 100.4e9
c44 = 19.6e9
c66 = 22.2e9

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

e13 = -9.6
e33 = 15.1
e15 = 12

tb,piez,1					!Define Piez. table
tbdata,2,e13				!Input Piezoelectric stress matrix [e]
tbdata,5,e33
tbdata,8,e13
tbdata,10,e15
tbdata,15,e15

mp,perx,1,804.6				!Permittivity (x direction)
mp,pery,1,659.7				!Permittivity (y direction)
mp,perz,1,804.6				!Permittivity (z direction)
mp,dens,1,7760				!Density

/com
/com  -- Material properties for Beam
/com
mp,ex,2,190E9				!Young modulus
mp,prxy,2,0.3				!Poisson coefficient
mp,dens,2,7850				!Density


!Balken modelieren

K,1,0,-0.5,2
K,2,10,-0.5,2
K,3,10,-0.5,-2
K,4,0,-0.5,-2

K,5,0,0.5,2
K,6,10,0.5,2
K,7,10,0.5,-2
K,8,0,0.5,-2

L,1,2 
L,2,3 
L,3,4
L,4,1

L,5,6
L,6,7
L,7,8
L,8,5

L,1,5
L,2,6
L,3,7
L,4,8

V,1,2,3,4,5,6,7,8


lesize,1,,,5
lesize,2,,,3
lesize,9,,,2

!Piezoelement modelieren

K,9,2.5,0.5,1
K,10,7.5,0.5,1
K,11,7.5,0.5,-1
K,12,2.5,0.5,-1

K,13,2.5,0.75,1
K,14,7.5,0.75,1
K,15,7.5,0.75,-1
K,16,2.5,0.75,-1

L,9,10 
L,10,11 
L,11,12
L,12,9

L,13,14
L,14,15
L,15,16
L,16,13

L,9,13
L,10,14
L,11,15
L,12,16

V,9,10,11,12,13,14,15,16


lesize,13,,,5
lesize,14,,,3
lesize,21,,,2


/com
/com  -- Mesh - Beam
/com
vsel,s,,,1
type,2
mat,2
mshape,0,3d				!Volume mesh using Hex.
mshkey,1				!Mapped meshing
vmesh,all
allsel

/com
/com  -- Mesh - Piezo
/com
vsel,s,,,2				!Piezomodell auswählen
type,1
mat,1
mshape,0,3d				!Volume mesh using Hex. 
mshkey,1				!Mapped meshing
vmesh,all
vsel,all


!Boundary Conditions Beam

esel,s,type,,2
nsel,r,loc,x,0
d,all,ux,0
d,all,uy,0
d,all,uz,0
allsel

!Boundary Conditions Piezoelemts

esel,s,type,,1
allsel,below,elem
nsel,r,loc,y,0.5
d,all,volt,vbot

esel,s,type,,1
allsel,below,elem
nsel,r,loc,y,0.75
d,all,volt,vtop
allsel



!Erstellen der Kontakflächen


et,3,conta174
keypot,3,12,3
et,4,targe170


asel,s,,,6
nsla,s,1
type,4
mat,3
real,2
esurf

asel,s,,,7
nsla,s,1
type,3
esurf

allsel

asel,s,,,6
nsla,s,1
type,3
mat,3
real,3
esurf

asel,s,,,7
nsla,s,1
type,4
esurf
allsel


finish


/solu							!Specify static analysis type
antype,static
pstres,on
solve
finish

/solu							!Specify dynamic analysis type
antype,modal
upcoord,1.0,on
pstres,on
modopt,lanb,4
mxpand,4
solve
finish