// Determinação da ordem do sistema de massa de fluido no tanque
//Parâmetros
hf = 2 //altura do líquido no tanque em [m]
l = 17.121 //Comprimento na direção do movimento em [m]
mf = 71000 //massa do líquido no tanque em [kg]
g = 9.8 //aceleração da gravidade em [m/s^2]
n = 5 //numero de modos

visc_cin = 0.64 * 10 ^ ( - 6 ) //viscosidade cinemática da gasolina a 20ºC em [m^2/s]

//Funções
function m_i = massa_i ( mf , i , hf , l )
m_i = mf * 8 * tanh (( 2 * i + 1 ) * %pi * hf / l ) / ( %pi ^ 3 * ( 2 * i + 1 ) ^ 3 * hf /l)
endfunction

function m_0 = massa_0 ( m , mf )
//m é lista com as n massa_i
n = length ( m )
m_0 = mf
for i = 1 : n
m_0 = m_0 - m ( i )
end
endfunction

function ki = rigidez_i ( hf , mf , g , l , i )
ki = ( 8 * mf * g * ( tanh (( 2 * i + 1 ) * %pi * hf / l )) ^ 2 ) / ( hf * ( 2 * i + 1 ) ^ 2 )
endfunction

function ci = amortecimento_i ( mi , fn , df )
ci = 2 * mi * fn * df
endfunction

//Resolução
DF = sqrt ( visc_cin / ( l ^ ( 3 / 2 ) * g ))
//definição dos parametros
//listas para cada um dos i modos
m = zeros ( n )
k = zeros ( n )
fn = zeros ( n )
c = zeros ( n )
for i = 1 : n
m ( i ) = massa_i ( mf,i,hf,l )
end
m0 = massa_0 ( m,mf )
for i = 1 : n
k ( i ) = rigidez_i ( hf,mf,g,l,i )
fn ( i ) = sqrt ( k ( i ) / m ( i ))
c ( i ) = amortecimento_i ( m ( i ) ,fn ( i ) ,DF )
end

// Resultados
disp ( m, 'massas' )
disp ( m0, 'm0' )
disp ( hf, 'alturas' )
disp ( hf, 'h0' )
disp ( k, 'coeficientes de mola' )
disp ( fn, 'frequencias naturais' )
disp ( c, 'coeficientes de amortecimento' )