finish
/clear,start

/title,Silikonzylinder
/prep7

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

et,1,solid5,0      !solid98

mp,ex,1,10e6       !E-Modul
mp,nuxy,1,0.499    !Poissonzahl
mp,dens,1,1120     !Dichte
!mp,rsvx,1,0.33e-3 !el.Widerstand
mp,c,1,1400
mp,kxx,1,0.15
mp,alpx,1,2e-4


mptemp,1,25,30,35,40,45,50
mptemp,7,55,60,65,70,75,80

mpdata,rsvx,1,1,0.33e-3,0.3301e-3,0.3303e-3,0.3308e-3,0.3317e-3,0.3325e-3    !el.Widerstand  
mpdata,rsvx,1,7,0.333e-3,0.3325e-3,0.333e-3,0.335e-3,0.3328e-3,0.324e-3
mpplot,rsvx,1

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

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

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

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

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

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

/pnum,line,1
lplot

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

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

TUNIF,Tu

/pnum,area,1
aplot

asel,s,area,,1
asel,a,area,,2
asel,a,area,,3
asel,a,area,,4
asel,a,area,,5
asel,a,area,,6
nsla,s,1
d,all,temp,Tu
alls

nsel,s,loc,z,0
d,all,ux
d,all,uy
d,all,uz
alls

nsel,s,loc,z,h
d,all,uz,-h*Stauchung  !Verschiebung
alls

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

/solu
antype,static
alls
solve

!elektr.Lösung

/prep7
upgeom,1,,,sr,rst

vplot
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  
alls

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

/solu
ANTYPE,STATIC
alls
solve

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