2次元軸対称問題の解

外圧を受ける円管

ラーメの式

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');

集中力を受ける半無限板

対向集中力を受ける半無限板