!****************************************************************
!**  Design Parameter                                          **
!****************************************************************

/com,...........................Geometrie........................

breite=0.0054
hoehe=0.0002
laenge=0.03
l_einsp=0.004
h_inv=0.0001

/com,........................material properties

/com,........................

s11 =  5.78e-12		      ! [m/N]
s12 = -1.01e-12		      ! 
s13 = -1.47e-12
s14 = -1.02e-12
s33 =  5.02e-12
s44 = 17.00e-12
s56 =    2*s14
s66 = 13.60e-12
d22 = 20.80e-12		      ! [m/V]] 
d31 = -0.85e-12			! 
d33 =  6.00e-12
d15 = 69.20e-12
d16 =   -2*d22
eps1=    44			      ! AULD
eps3=    29
density=4700	    		! AULD 
d=0.01
spannung_unten=0
nmode=1
startfreq_modal=200
endfreq_modal=300
harm_nmode=1
*dim,frequency,,nmode

/com,============================================================
/com,==                         PRE-PROCESSING                 ==
/com,============================================================
/com,
/prep7

/com,............................element types

et,1,solid226,1001
et,2,circu94,0
r,2,100000

/com,............................geometry

block,0,laenge,0,breite,0,hoehe
lsel,s,tan1,z,1
lsel,a,tan1,z,-1
lesize,all,hoehe/4

lsel,s,tan1,x,1
lsel,a,tan1,x,-1
lesize,all,laenge/30

lsel,s,tan1,y,1
lsel,a,tan1,y,-1
lesize,all,breite/10

/com,............................material properties

tb,anel,1,,,1
tbdata, 1,s11,s12,s13, s14,  0,  0
tbdata, 7,    s11,s13,-s14,  0,  0  
tbdata,12,        s33,   0,  0,  0
tbdata,16,             s66,  0,s56
tbdata,19,                 s44,  0
tbdata,21,                     s44

tb,piez,1,,,1
tbdata, 1,  0,-d22,d31
tbdata, 4,  0, d22,d31
tbdata, 7,  0,   0,d33
tbdata,10,d16,   0,  0
tbdata,13,  0, d15,  0
tbdata,16,d15,   0,  0

mp,perx,1,eps1
mp,pery,1,eps1
mp,perz,1,eps3
mp,dens,1,density

/com,.............................meshing

mat,1
type,1
numstr,node,10
numstr,elem,10
vmesh,all
eplo

/com,.............................components

nsel,s,loc,x,0,l_einsp
nsel,r,loc,z,0
nplo
cm,n_einsp,node

/com,******************************************
/com,*********   SOLUTION 1 - MODAL open   ***
/com,******************************************

/solu
antype,modal
modopt,lanb,nmode,startfreq_modal,endfreq_modal
outres,all,all
mxpand,nmode

/com,........................boundaries

cmsel,s,n_einsp
d,all,ux,0
d,all,uy,0
d,all,uz,0

nsel,s,loc,z,0
d,all,volt,spannung_unten

nsel,s,loc,z,hoehe	
cp,1,volt,all

alls
/pbc,all,,1
eplo

solve

/com,*********************************************
/com,*******   POSTPROCESSING 1 - MODAL OPEN   ***
/com,*********************************************

/post1

*do,i,1,nmode
 set,1,i
 pldisp
 *get,f,mode,i,freq
 frequency(i) = f			! [Hz]
*enddo

/com,*********************************************
/com,*******   PREPROCESSING 2 - HARMONIC      ***
/com,*********************************************

/prep7

n,1,laenge/3,breite/2,hoehe*6
n,2,laenge/3,breite/2,-hoehe*7
n_bot=2
n_top=1
real,2
type,2
e,n_bot,n_top

/com,********************************************
/com,*********   SOLUTION 2 - harmonic  *********
/com,********************************************

/solu
antype,harmic 
hropt,full
hrout,off
outpr,nsol
outres,all,all
harfrq,frequency(harm_nmode)*0.9,frequency(harm_nmode)*1.1
nsubst,15

/com,...........................damping..........................
/com,                           
   
mdamp,1,d

/com,...........................Randbedingungen

nsel,s,loc,x,laenge
f,all,fz,-0.001

alls
/pbc,u,,1
eplo
solve

/com,********************************************
/com,*******   POSTPROCESSING 2 -harmonic  ******
/com,********************************************

/post26
numvar,200
/com,

/com,...........................misc.............................
/com,                           2 Kontroll-Knoten finden

/erase
*dim,nt,array,2
nt(1)=node(0.025,0.003,-0.0021 )
nt(2)=node(0.025,0.003,0.0027)
nsol,11,nt(1),u,z,U_1
nsol,12,nt(2),u,z,U_2
nsel,all
/grid,1
/axlab,x,Frequenz [Hz]
/xrange,280,300
/axlab,y,Verschiebung in z-Richtung [m]
plvar,11,12

/erase
/win,1,top
esol,13,e_resistor,,smisc,1,VOLT_Resistor
esol,14,e_resistor,,smisc,2,CURRENT_Resistor
esol,15,e_resistor,,nmisc,1,POWER_Resistor
store
/grid,1
/axlab,x,Frequenz [Hz]
/xrange,280,300
/axlab,y,allerlei [x]

plvar,13,14,15

/noerase
/win,1,bot
nsol,16,n_top,volt,VOLT_Resistor_top
nsol,17,n_bot,volt,VOLT_Resistor_bot
/grid,1
/axlab,x,Frequenz [Hz]
/xrange,280,300
/axlab,y,Spannung [V]

plvar,16,17

/erase
/win,1,full

/eof