% Parâmetros de Sistema Físico

clc
clearvars
close all

% dados dimensionais
l=1; % comprimento da torre
h=1; % espessura da torre
b=1; % largura da torre

% dados inerciais
mm=1; % massa do motor
mh=1; % massa da hélice
ms=1; % massa do suporte do motor

% dados materiais
ro=1; % densidade do material da torre
E=1;  % módulo de elasticidade do material da torre

% parâmetros intermediários
I=b*h^3/12;  % momento de inércia de área da torre
k=3*E*I/l^3; % constante de rigidez da torre
mt=ro*l*h*b; % massa da torre
m=mm+mh+ms+0.23*mt; % massa total

% parâmetros de interesse
fn=sqrt(k/m)/(2*pi); % frequência natural do sistema

disp(['A frequência natural do sistema é ',num2str(fn),' Hz.']);
disp(' ')