kill(all);

tensione circonferenziale al bordo interno

stt_ri(ri,re,pi,pe):=(pi*(ri^2+re^2)-2*pe*re^2)/(re^2-ri^2);

stt_ri(10,20,1,0),numer;

tensione ideale, noto srr_ri=-pi e trt_ri=0

sidVM_ri(ri,re,pi,pe):=sqrt( stt_ri(ri,re,pi,pe)^2 + (-pi)^2 - (-pi)*stt_ri(ri,re,pi,pe) );

sidVM_ri(10,20,1,0),numer;

spostamento al bordo interno

ur_ri(ri,re,pi,pe,E,nu) := ri/E*(stt_ri(ri,re,pi,pe)-nu*(-pi));

ur_ri(10,20,1,0,70000,0.3),numer;

forzamento albero-mozzo

pf(ri,rm,re,E,dr):= E*dr/rm/( (re^2+rm^2)/(re^2-rm^2)+(rm^2+ri^2)/(rm^2-ri^2) );

pf(0,10,20,70000,0.05);