finish
/clear,start

/title,Silikonzylinder
/prep7

!Materialparameter!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

et,1,solid226,101

mp,ex,1,1e9        !E-Modul
mp,nuxy,1,0.499    !Poissonzahl
mp,rsvx,1,1/1e8    !el.Widerstand  
mp,dens,1,1.12E-9  !Dichte

!Piezoresistive Konstanten
p11=6.5E-5     
p12=-1.1E-5   
p44=138.1E-5 

TB,PZRS,1                    ! piezoresistive stress matrix
TBDATA,1,p11,p12,p12
TBDATA,7,p12,p11,p12
TBDATA,13,p12,p12,p11
TBDATA,22,p44

!el. und mech Parameter!!!!!!!!!!!!!!!!!

Strom=1e-3     !eingeprägter Strom
Druck=1e-4     !Einheit [Pa]
Stauchung=0.05  !Stauchung in % 

!Geometrie!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

h=0.01
r=0.005
mw=10
k,1,,,,

pcirc,0,r,0,360
voffset,1,h,20

!Vernetzung!!!!!!!!!!!!!!!!!!!!!!

/pnum,line,1
lplot

mshape,1,3,D
mshkey,0
smrtsize,8
vmesh,1

!Einspannung!!!!!!!!!!!!!!!!!!!!!!!

/pnum,area,1
aplot

asel,s,area,,1
nsla,s,1
!d,all,ux
!d,all,uy
d,all,uz
alls

asel,s,area,,2
nsla,s,1
d,all,ux
d,all,uy
d,all,uz,-h*Stauchung   !Verschiebung
alls

!mech. Lösung!!!!!!!!!!!!!!!!!!!!!!!!

/solu
antype,static
nlgeom,on
solve

!elektr.Lösung

/prep7
upgeom,1,,,sr,rst

nsel,s,loc,z,h-h*Stauchung  !Fläche 2 auswählen
cp,1,volt,all               ! alle el. Freiheitsgrade  
*get,node1,node,,num,min
f,node1,amps,-strom  
!sf,all,CHRGS,10
alls

nsel,s,loc,z,0
CP,2,VOLT,ALL       
d,all,volt,0
alls
vplot

/solu
ANTYPE,STATIC
!CNVTOL,AMPS,1,1.E-3   !typ.Werte 1A Toleranz 1mA
solve

/post1
plnsol,volt
*get,spannung,node,node1,volt
widerst=spannung/strom