外圧を受ける円管
ラーメの式
clear; close all; a=1; b=2; p=1; r = linspace(a,b,20); sr = -b^2*p/(b^2-a^2)*(1-a^2./r.^2); ss = -b^2*p/(b^2-a^2)*(1+a^2./r.^2); plot(r,sr); title('\sigma_r'); figure(); plot(r,ss); title('\sigma_\theta');
回転円板
clear; close all; v = 0.3; b = 1; rho = 1; omega = 1; r = linspace(0,b,20); sr = (3+v)*rho*omega^2*b^2/8*(1-r.^2/b^2); ss = (3+v)*rho*omega^2*b^2/8*(1-(1+3*v)/(3+v)*r.^2/b^2); plot(r,sr); title('\sigma_r'); figure(); plot(r,ss); title('\sigma_\theta');
集中力を受ける半無限板