finish
/clear,start

/title,Silikonzylinder
/prep7

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

et,1,solid62       !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

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!elastische Verformungsoption

TB,Melas,1,1,3
TBTEMP,25

TBPT,Defi,0.25,2.5e+6   !Dehnung %, Spannung MPa
TBPT,Defi,1.25,3e+6
TBPT,Defi,1.5,4.5e+6

TBPLOT,Melas,1

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

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

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

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

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

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

/pnum,line,1
lplot

CM,_Y,VOLU  
VSEL, , , ,       1 
CM,_Y1,VOLU 
CHKMSH,'VOLU'   
CMSEL,S,_Y  
!*  
VMESH,_Y1   
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  


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

TUNIF,Tu

/pnum,area,1
aplot

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,trans
kbc,1
time,1
nsub,5
nlgeom,on
alls
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  
alls

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

/solu
antype,trans
kbc,1
time,1
nsub,5
neqit,10
nlgeom,on
alls
solve

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