軸の危険速度の計算

clear; close all; clc;
E = 206e9;
d = 0.1;
I = pi/64*d^4;
m = 30;
l = 2;
a = l/2;
b = l/2;
k = 6*E*I*l/(a*b*(l^2-a^2-b^2));
nc = 60/(2*pi)*sqrt(k/m);
fprintf('危険速度: %.0f [rpm]\n',nc);

実行結果

危険速度: 4294 [rpm]