外圧を受ける円管 ラーメの式 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'); 回転…