//Dimensions de la lame (en m)

Long=0.25;    
larg=0.03;   
ep=0.003/6; // 3 mm est assez important, 0.5 mm me parait suffisant

S = ep*larg;    

Iz=((larg*(ep)^3))/12;
Io = ((larg*ep)*((larg)^2 + (ep)^2))/12;

//Propriétés matériau

E = 69e9 // polyéthylène : 0.5e9; //Pa (N/m^2)    // Acier : 210e9   Alu : 69e9
rau=2700;  //kg/m^3        // Acier : 7500-8100   Alu : 2700
v = 0.346;    //Alu : 0.346     Acier : 0.27-0.30

m=rau*S*Long; //kg

//Bande de pulsation accéléromètre en rad/s

Rmin = 0
Rmax = (1600*(2*%pi))

BP = Rmax - Rmin;

//
// Echantillonnage : nombre de points de mesure par période

npt_ech=20;

//Nombre de pulsations propres à calculer : 

k=3;


//Calcul des pulsations propres longitudinales : 

WL = [];

for i=1:k
    
    WL(i) = (2*i -1)*(%pi/2)*(1/Long)*sqrt(E/rau)
    if(Rmin <= npt_ech*WL(i) && npt_ech*WL(i) <= Rmax)
        printf('\n La pulsation propre longitudinale est dans la bande de pulsation du capteur')
    else
        printf('\n La pulsation propre longitudinale n est pas dans la bande de pulsation du capteur')
    end 
end

//Calcul des pulsations propres transversales : 

WT = [];

w1 = (1.87)^2*sqrt((E*Iz)/(m*(Long^3)))

WT(1) = w1;

    if(Rmin <= npt_ech*WT(1) && npt_ech*WT(1) <= Rmax)
        printf('\n La pulsation propre transversale est dans la bande de pulsation du capteur')
    else
        printf('\n La pulsation propre transversale n est pas dans la bande de pulsation du capteur')
    end 
    
for i=2:k
    
    WT(i) = (((2*i-1)*(%pi/2))^2)*sqrt((E*Iz)/(m*(Long^3)))
    if(Rmin <= npt_ech*WT(i) && npt_ech*WT(i) <= Rmax)
        printf('\n La pulsation propre transversale est dans la bande de pulsation du capteur')
    else
        printf('\n La pulsation propre transversale n est pas dans la bande de pulsation du capteur')
    end 
end


//Calcul des pulsations propres angulaires : 

G = E/(2*(1+v));

J = (larg*(ep)^3)*((16/3)-3.36*(ep/larg)*(1 - (ep^4)/(12*larg^4))); //Approximation à 4%

WA = [];

for i=1:k
    
    WA(i) = (2*i -1)*(%pi/2)*(1/Long)*sqrt((G*J)/(rau*Io))
    if(Rmin <= npt_ech*WA(i) && npt_ech*WA(i) <= Rmax)
        printf('\n La pulsation propre angulaire est dans la bande de pulsation du capteur')
    else
        printf('\n La pulsation propre angulaire n est pas dans la bande de pulsation du capteur')
    end 
end


//Charge critique en flambement : 

//Position masse sur la lame

L1 = Long; 
L2 = Long/2;
L3 = Long/3;

F = ((%pi^2)*(E*Iz))/(2*L1)^2;

k_sec=0.5 // coefficient de sécurité

M = k_sec*F/9.8; //masse critique en la longueur selectionnée (en kg)

disp(M/m)