!/solu
!MODALNA
!antype, modal
!MODOPT,LANB,10
!EQSLV,SPAR
!MXPAND,10, , ,0
!LUMPM,0
!STRES,0
!MODOPT,LANB,10,0,0, ,OFF
!/STATUS,SOLU
!SOLVE


fini
/clear
/prep7

!Stałe Materiałowe piezo
cE11=10.76e10
cE12=10.04e10
cE13=6.312e10
cE33=6.385e10
cE44=1.926e10
CE66=(cE11-cE12)/2

e31=-9.6
e33=15.1
e15=12
Epsilon11=1110
Epsilon33=852
ro=7760
Epsilon0=8.85e-12


!wprowadzenie materiału piezo
TB,Anel,1
tbdata,1,ce11,ce12,ce13
TBDATA,7,ce11,ce13
TBDATA,12,ce33
TBDATA,16,ce66
TBDATA,19,ce44
TBDATA,21,ce44

MP,dens,1,ro

TB,Piez,1
TBdata,3,e31
TBdata,6,e31
TBdata,9,e33
TBdata,14,e15
TBdata,16,e15

MP,PERX,1,epsilon11
MP,PERY,1,epsilon11
MP,PERZ,1,epsilon33
MP,LSST,1,0.013
DMPRAT,0.015

!materiał dociskający
mp,ex,2,16000000000
mp,prxy,2,0.45
mp,dens,2,11340


ET,1,solid226,1001
ET,2,solid95,1



xcorner=0
ycorner=0
lenght=25e-3
width=2.5e-3
height=10e-3

BLC4, xcorner, ycorner,lenght,width,height !top
BLC4, xcorner, ycorner,lenght,-width,height !bot

!masa dociskająca
BLC4,0,2.5e-3,25e-3,8e-3,10e-3

LOCAL,11,,,,,,-90
LOCAL,12,,,,,,90

mat,1 /Sets the element material attribute pointer.
type,1 / Sets the element type attribute pointer.
esize,1e-3

ESYS,12
vmesh,1
ESYS,11
vmesh,2

mat,2
type,2
esize,1e-3

vmesh,3

nummrg,node / Merges coincident or equivalently defined items.
numcmp,node / Compresses the numbering of defined items.


nsel,s,loc,z,width
*get,gora,node,0,num,min
cp,1,volt,all
nsel,all

nsel,s,loc,z,-width
*get,dol,node,0,num,min
cp,2,volt,all
nsel,all


nsel,s,loc,z,0
*get,srodek,node,nie 0,num,min
cp,3,volt,all
nsel,all
 
!d,gora,volt,0
!d,dol,volt,0
!d,srodek,volt,1

!rezystor
REZYST=100
et,3,CIRCU94,0
R,1,100
R,2,0,100,10
N,500001,20e-3,-10e-3,-10e-3
TYPE,3
REAL,1
E,srodek,500001
D,500001,VOLT,0



DA,9,UX,0
DA,9,Uy,0
DA,9,Uz,0


*DEL,_FNCNAME   
*DEL,_FNCMTID   
*DEL,_FNCCSYS   
*SET,_FNCNAME,'Przysp'  
*SET,_FNCCSYS,0
! /INPUT,AccelFunc.func,,,1
*DIM,%_FNCNAME%,TABLE,6,8,1,,,,%_FNCCSYS%   
!   
! Begin of equation: 9.8*sin(8*3.14159265358979310*{TIME})  
*SET,%_FNCNAME%(0,0,1), 0.0, -999   
*SET,%_FNCNAME%(2,0,1), 0.0
*SET,%_FNCNAME%(3,0,1), 0.0
*SET,%_FNCNAME%(4,0,1), 0.0
*SET,%_FNCNAME%(5,0,1), 0.0
*SET,%_FNCNAME%(6,0,1), 0.0
*SET,%_FNCNAME%(0,1,1), 1.0, -1, 0, 4, 0, 0, 0  
*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 3.14159265358979310, 0, 0, -1   
*SET,%_FNCNAME%(0,3,1),   0, -3, 0, 1, -1, 3, -2
*SET,%_FNCNAME%(0,4,1), 0.0, -1, 0, 1, -3, 3, 1
*SET,%_FNCNAME%(0,5,1), 0.0, -1, 9, 1, -1, 0, 0
*SET,%_FNCNAME%(0,6,1), 0.0, -2, 0, 9.8, 0, 0, -1   
*SET,%_FNCNAME%(0,7,1), 0.0, -3, 0, 1, -2, 3, -1
*SET,%_FNCNAME%(0,8,1), 0.0, 99, 0, 1, -3, 0, 0
! End of equation: 9.8*sin(4*3.14159265358979310*{TIME})
!-->


ACEL,,%Przysp%

!/solu
!Antype,static
!solve

/SOLU                      ! Enter SOLUTION
ANTYPE,TRANS               ! Transient analysis

DMPRAT,0.01


nsubst,40
time,0.5

outres,nsol,all
outres,esol,all
outres,misc,all


tintp,,0.25,0.5,0.5
solve