Groupe : 38 crapet franck tune pierre-Damien zucchetta vincent Date : 28 Juin 2013

Bibliographie

• 
  Centrale Paris : Rapport final TMO-08.pdf
  
• 
  Centrale Paris : Fichier MATLAB
  
• 
  ESTACA : Structure/Estaca 4a Beeplane rapport.pdf (Novembre 2012)
  
• 
  ESTACA : Rapport Bee-Plane final
  
• 
  Project description for TRL2 revision 2.2
  
• 
  Technoplane.com
  
• 
  Wiki Technoplane
  
• 
  Technique de l’Ingénieur
  

Annexe 1 : Script MATLAB de la fonction ‘structure_bee_plane.m’

function [ Structure , Moteurs , numero_aile_gauche , numero_aile_droite , numero_moteur_gauche , numero_moteur_droit , numero_fuselage] = structure_bee_plane( phase_decollage , vol_avec_basket )

%Dans cette fonction sont définis les différents éléments qui composent la

%structure de la Bee ainsi que de la Basket

% L'origine est prise comme étant celle du nez de la Bee puisque le

% Bee-Plane doit être capable de voler sans basket

%% Liste de paramètres
d_moteur_fuselage=8.5; %distance du moteur par rapport à l'axe du fuselage suivant Z

d_train_atterrissage_avant=5; %distance du point d'accroche du train atterrissage avant par raport à l'axe du fuselage suivant Z

d_train_atterrissage_arriere=5; %distance du point d'accroche du train atterrissage arrière par raport à l'axe du fuselage suivant Z

masse_train_avant=150; %Masse d'un train avant complet

masse_train_arriere=150; %Masse d'un train arrière complet

hauteur_turbofan=3.5; %Un moteur de la taille 'hauteur_turbofan' représenté par un cylindre un inséré en bout de fuselage. La taille du fuselage,turbofan compris, vaut 'hauteur_fuselage_turbofan'. La taille du fuselage revient donc à la soustraction du dernier moins le premier

hauteur_fuselage_turbofan=21;

decalage_empennage_verti=1.5;

longueur_basket=25;

%% Définition des élements composant la structure du Bee-Plane par incrémentation du paramètre 'ii'

ii=1; %Numero de l'élément de la structure
%% Définition de la partie gauche de la Bee puis de la partie droite comme symétrique
%% Elements de l'Aile Gauche

Structure.Element(ii).Nom='Aile Gauche';

Structure.Element(ii).Type='Aile';

Structure.Element(ii).Masse=4000*21/20*5.2/5;

Structure.Element(ii).Incidence=2.3; %en degrés

Structure.Element(ii).Longueur=21;

Structure.Element(ii).Naca='4415';

Structure.Element(ii).GouD='G';

Structure.Element(ii).Raccord_Avion=[4,0,0];% en général : la côte est positive pour une aile gauche, négative pour une aile droite

Structure.Element(ii).Raccord_Element=[0,0,0];

Structure.Element(ii).Rotation=[0,Structure.Element(ii).Incidence*pi/180,0];%en radians : [azimut, assiette, gite]

Structure.Element(ii).C_Root=5.2;

Structure.Element(ii).Rapport_Bases=2/5;% rapport de la corde en bout d'aile sur la corde à la base de l'aile

Structure.Element(ii).Lambda_0=17*pi/180;%en radians

numero_aile_gauche=ii;%Variable muette pour définir les éléments rattachés

ii=ii+1;
Structure.Element(ii).Nom='Moteur Gauche';%Les coordonnées des ailes sont nécessaires pour placer les moteurs

Structure.Element(ii).Type='Moteur';

Structure.Element(ii).Masse=1890;

Structure.Element(ii).Raccord_Avion=[5,0,d_moteur_fuselage];

Structure.Element(ii).Raccord_Avion(1)=Structure.Element(numero_aile_gauche).Raccord_Avion(1)+sin(Structure.Element(numero_aile_gauche).Lambda_0)*Structure.Element(ii).Raccord_Avion(3);

%la ligne précédente permet d'ajuster automatiquement le point d'accroche des moteurs en fonction de l'inclinaison des ailes et de leur point d'accroche le long du fuselage, la valeur '10' renseignée précédemment ne joue donc aucun rôle

Structure.Element(ii).Raccord_Element=[0,0.9,0];

Structure.Element(ii).Rotation=[0,0,0];

Structure.Element(ii).Hauteur=3.5;

Structure.Element(ii).Rayon=0.9;

Structure.Element(ii).Barycentre=[Structure.Element(ii).Hauteur/2,0,0];

numero_moteur_gauche=ii;%Variable muette pour définir l'élément moteur droit comme symétrique par rapport plan xy

ii=ii+1;

%% Empennage Gauche
Structure.Element(ii).Nom='Jonction Empennage Gauche Horizontal';

Structure.Element(ii).Type='Aile';

Structure.Element(ii).Masse=800;

Structure.Element(ii).Incidence=0; %en degrés

Structure.Element(ii).Longueur=6;

Structure.Element(ii).Naca='4415';

Structure.Element(ii).GouD='G';

Structure.Element(ii).C_Root=3 ;

Structure.Element(ii).Rapport_Bases=.5;% rapport de la corde en bout d'aile sur la corde à la base de l'aile

Structure.Element(ii).Lambda_0=.2;%en radians

Structure.Element(ii).Raccord_Avion=[hauteur_fuselage_turbofan-Structure.Element(ii).C_Root,0.1,0];% en général : la côte est positive pour une aile gauche, négative pour une aile droite

Structure.Element(ii).Raccord_Element=[0,0,0];

Structure.Element(ii).Rotation=[0,0,0];%en radians : [azimut, assiette, gite]

numero_empennage_horiz_gauche=ii;%Variable muette pour définir l'élément Jonction Empennage Droit Horizontal comme symétrique par rapport plan xy

ii=ii+1;
Structure.Element(ii).Nom='Jonction Empennage Gauche Vertical';

Structure.Element(ii).Type='Aile';

Structure.Element(ii).Masse=400;

Structure.Element(ii).Incidence=0; %en degrés

Structure.Element(ii).Longueur=3;

Structure.Element(ii).Naca='0006';

Structure.Element(ii).GouD='G';

Structure.Element(ii).Raccord_Avion=[Structure.Element(numero_empennage_horiz_gauche).Raccord_Avion(1)+Structure.Element(numero_empennage_horiz_gauche).C_Root/2,Structure.Element(numero_empennage_horiz_gauche).Raccord_Avion(2),decalage_empennage_verti];% en général : la côte est positive pour une aile gauche, négative pour une aile droite

Structure.Element(ii).Raccord_Element=[0,0,0];

Structure.Element(ii).Rotation=[0,0,pi/2];%en radians : [azimut, assiette, gite]

Structure.Element(ii).C_Root=1.5 ;

Structure.Element(ii).Rapport_Bases=1;% rapport de la corde en bout d'aile sur la corde à la base de l'aile

Structure.Element(ii).Lambda_0=0.7;%en radians Compris entre ]-pi/2;pi/2[

numero_empennage_verti_gauche=ii;%Variable muette pour définir l'élément Jonction Empennage Droit Vertical comme symétrique par rapport plan xy

ii=ii+1;
Structure.Element(ii).Nom='Empennage Gauche';

Structure.Element(ii).Type='Aile';

Structure.Element(ii).Masse=700*3.5/1.5*5/6;

Structure.Element(ii).Incidence=0; %en degrés

Structure.Element(ii).Longueur=5;

Structure.Element(ii).Naca='4415';

Structure.Element(ii).GouD='G';

Structure.Element(ii).Rotation=[0,0,0];%en radians : [azimut, assiette, gite]

Structure.Element(ii).C_Root=3.5;

Structure.Element(ii).Rapport_Bases=0.7/2.3;% rapport de la corde en bout d'aile sur la corde à la base de l'aile

Structure.Element(ii).Lambda_0=atan(Structure.Element(ii).C_Root*(1-Structure.Element(ii).Rapport_Bases)/Structure.Element(ii).Longueur);%en radians, est modifié directement dans l'analyse paramétrique pour que le bord de fuite reste droit
if mod(Structure.Element(numero_empennage_verti_gauche).Lambda_0,2*pi)==0

ffr=1;

else

ffr=Structure.Element(numero_empennage_verti_gauche).Longueur/tan(pi/2-Structure.Element(numero_empennage_verti_gauche).Lambda_0);

end

Structure.Element(ii).Raccord_Avion=[Structure.Element(numero_empennage_verti_gauche).Raccord_Avion(1)+ffr-Structure.Element(ii).C_Root/2,Structure.Element(numero_empennage_verti_gauche).Raccord_Avion(2)+Structure.Element(numero_empennage_verti_gauche).Longueur,0];% en général : la côte est positive pour une aile gauche, négative pour une aile droite

Structure.Element(ii).Raccord_Element=[0,0,0];

numero_empennage_gauche=ii;%Variable muette pour définir l'élément Empennage droit comme symétrique par rapport plan xy

ii=ii+1;
%% Partie Droite de la Bee comme symétrique de la gauche
%% Elements de l'Aile Droite

% L'aile Droite et ses éléments associés, sont définis comme le symétrique

% des éléments de l'Aile Gauche
%Rentrer les paramètres des ailes dans la partie Aile Gauche
Structure.Element(ii).Nom='Aile Droite';

Structure.Element(ii).Type='Aile';

Structure.Element(ii).Masse=Structure.Element(numero_aile_gauche).Masse;

Structure.Element(ii).Incidence=Structure.Element(numero_aile_gauche).Incidence;

Structure.Element(ii).Longueur=Structure.Element(numero_aile_gauche).Longueur;

Structure.Element(ii).Naca=Structure.Element(numero_aile_gauche).Naca;

Structure.Element(ii).GouD='D';

Structure.Element(ii).Raccord_Avion=Structure.Element(numero_aile_gauche).Raccord_Avion; Structure.Element(ii).Raccord_Avion(3)=-Structure.Element(ii).Raccord_Avion(3);% en général : la côte est positive pour une aile gauche, négative pour une aile droite

Structure.Element(ii).Raccord_Element=[0,0,0];

Structure.Element(ii).Rotation=Structure.Element(numero_aile_gauche).Rotation;%en radians : [azimut, assiette, gite]

Structure.Element(ii).C_Root=Structure.Element(numero_aile_gauche).C_Root;

Structure.Element(ii).Rapport_Bases=Structure.Element(numero_aile_gauche).Rapport_Bases;% rapport de la corde en bout d'aile sur la corde à la base de l'aile

Structure.Element(ii).Lambda_0=Structure.Element(numero_aile_gauche).Lambda_0; %en radians

numero_aile_droite=ii;%Variable muette pour définir les éléments rattachés

ii=ii+1;
%Les coordonnées des ailes sont nécessaires pour placer les moteurs. Le

%moteur droit est donc défini après l'aile droite

Structure.Element(ii).Nom='Moteur Droit';

Structure.Element(ii).Type='Moteur';

Structure.Element(ii).Masse=Structure.Element(numero_moteur_gauche).Masse;

Structure.Element(ii).Raccord_Avion=Structure.Element(numero_moteur_gauche).Raccord_Avion;

Structure.Element(ii).Raccord_Avion(3)=-Structure.Element(ii).Raccord_Avion(3);

Structure.Element(ii).Raccord_Element=[0,0.9,0];

Structure.Element(ii).Rotation=Structure.Element(numero_moteur_gauche).Rotation;

Structure.Element(ii).Hauteur=Structure.Element(numero_moteur_gauche).Hauteur;

Structure.Element(ii).Rayon=Structure.Element(numero_moteur_gauche).Rayon;

Structure.Element(ii).Barycentre=Structure.Element(numero_moteur_gauche).Barycentre;

numero_moteur_droit=ii;

ii=ii+1;
%% Empennage Droit

Structure.Element(ii).Nom='Empennage Droit';

Structure.Element(ii).Type='Aile';

Structure.Element(ii).Masse=Structure.Element(numero_empennage_gauche).Masse;

Structure.Element(ii).Incidence=Structure.Element(numero_empennage_gauche).Incidence; %en degrés

Structure.Element(ii).Longueur=Structure.Element(numero_empennage_gauche).Longueur;

Structure.Element(ii).Naca=Structure.Element(numero_empennage_gauche).Naca;

Structure.Element(ii).GouD='D';

Structure.Element(ii).Raccord_Avion=Structure.Element(numero_empennage_gauche).Raccord_Avion;

Structure.Element(ii).Raccord_Avion(3)=-Structure.Element(ii).Raccord_Avion(3);% en général : la côte est positive pour une aile gauche, négative pour une aile droite

Structure.Element(ii).Raccord_Element=[0,0,0];

Structure.Element(ii).Rotation=Structure.Element(numero_empennage_gauche).Rotation;%en radians : [azimut, assiette, gite]

Structure.Element(ii).C_Root=Structure.Element(numero_empennage_gauche).C_Root;

Structure.Element(ii).Rapport_Bases=Structure.Element(numero_empennage_gauche).Rapport_Bases;% rapport de la corde en bout d'aile sur la corde à la base de l'aile

Structure.Element(ii).Lambda_0=Structure.Element(numero_empennage_gauche).Lambda_0;%en radians, est modifié directement dans l'analyse paramétrique pour que le bord de fuite reste droit

ii=ii+1;
Structure.Element(ii).Nom='Jonction Empennage Droit Horizontal';

Structure.Element(ii).Type='Aile';

Structure.Element(ii).Masse=Structure.Element(numero_empennage_horiz_gauche).Masse;

Structure.Element(ii).Incidence=Structure.Element(numero_empennage_horiz_gauche).Incidence; %en degrés

Structure.Element(ii).Longueur=Structure.Element(numero_empennage_horiz_gauche).Longueur;

Structure.Element(ii).Naca=Structure.Element(numero_empennage_horiz_gauche).Naca;

Structure.Element(ii).GouD='D';

Structure.Element(ii).Raccord_Avion=Structure.Element(numero_empennage_horiz_gauche).Raccord_Avion;% en général : la côte est positive pour une aile gauche, négative pour une aile droite

Structure.Element(ii).Raccord_Avion(3)=-Structure.Element(ii).Raccord_Avion(3); %Z=-Z

Structure.Element(ii).Raccord_Element=[0,0,0];

Structure.Element(ii).Rotation=Structure.Element(numero_empennage_horiz_gauche).Rotation;%en radians : [azimut, assiette, gite]

Structure.Element(ii).C_Root=Structure.Element(numero_empennage_horiz_gauche).C_Root ;

Structure.Element(ii).Rapport_Bases=Structure.Element(numero_empennage_horiz_gauche).Rapport_Bases;% rapport de la corde en bout d'aile sur la corde à la base de l'aile

Structure.Element(ii).Lambda_0=Structure.Element(numero_empennage_horiz_gauche).Lambda_0;%en radians

numero_empennage_horiz_droit=ii;

ii=ii+1;
Structure.Element(ii).Nom='Jonction Empennage Droit Vertical';

Structure.Element(ii).Type='Aile';

Structure.Element(ii).Masse=Structure.Element(numero_empennage_verti_gauche).Masse;

Structure.Element(ii).Incidence=Structure.Element(numero_empennage_verti_gauche).Incidence; %en degrés

Structure.Element(ii).Longueur=Structure.Element(numero_empennage_verti_gauche).Longueur;

Structure.Element(ii).Naca=Structure.Element(numero_empennage_verti_gauche).Naca;

Structure.Element(ii).GouD='D';

Structure.Element(ii).Raccord_Avion=Structure.Element(numero_empennage_verti_gauche).Raccord_Avion;

Structure.Element(ii).Raccord_Avion(3)=-Structure.Element(ii).Raccord_Avion(3);% en général : la côte est positive pour une aile gauche, négative pour une aile droite

Structure.Element(ii).Raccord_Element=[0,0,0];

Structure.Element(ii).Rotation=Structure.Element(numero_empennage_verti_gauche).Rotation;%en radians : [azimut, assiette, gite]

Structure.Element(ii).Rotation(3)=-Structure.Element(ii).Rotation(3);

Structure.Element(ii).C_Root=Structure.Element(numero_empennage_verti_gauche).C_Root;

Structure.Element(ii).Rapport_Bases=Structure.Element(numero_empennage_verti_gauche).Rapport_Bases;% rapport de la corde en bout d'aile sur la corde à la base de l'aile

Structure.Element(ii).Lambda_0=Structure.Element(numero_empennage_verti_gauche).Lambda_0;%en radians

ii=ii+1;
%%

if phase_decollage==true

Structure.Element(ii).Nom='Train Atterrissage Avant Gauche';

Structure.Element(ii).Type='Train Atterrissage';

Structure.Element(ii).Masse=masse_train_avant;

Structure.Element(ii).GouD='G';%Train gauche ou droit

Structure.Element(ii).Rayon_Support=0.15;%rayon du cylindre qui lie la roue à l'aile

Structure.Element(ii).Longueur_Support=3;%longueur du cylindre

Structure.Element(ii).Masse_Support=100;

Structure.Element(ii).Masse_Roue=50;

Structure.Element(ii).Rayon_Roue=0.5;%Rayon de la roue

Structure.Element(ii).Epaisseur_Roue=0.2;%Epaisseur de la roue

Structure.Element(ii).Raccord_Avion=[5,0,d_train_atterrissage_avant];

%Ligne suivante à définir avec le groupe 'Train d'atterrissage'

Structure.Element(ii).Raccord_Avion(1)=Structure.Element(numero_aile_gauche).C_Root/2+Structure.Element(numero_aile_gauche).Raccord_Avion(1)+sin(Structure.Element(numero_aile_gauche).Lambda_0)*Structure.Element(ii).Raccord_Avion(3);

Structure.Element(ii).Raccord_Element=[0,0,0];

Structure.Element(ii).Rotation=[0,0,0];

Structure.Element(ii).Coeff_Lift=0;

numero_train_avant_gauche=ii;%Variable muette pour définir l'élément Train Aterrissage Avant droit droit comme symétrique par rapport plan xy

ii=ii+1;

Structure.Element(ii).Nom='Train Atterrissage Arriere Gauche';

Structure.Element(ii).Type='Train Atterrissage';

Structure.Element(ii).Masse=masse_train_arriere;

Structure.Element(ii).GouD='G';%Train gauche ou droit

Structure.Element(ii).Rayon_Support=0.15;%rayon du cylindre qui lie la roue à l'aile

Structure.Element(ii).Longueur_Support=3;%longueur du cylindre

Structure.Element(ii).Masse_Support=100;

Structure.Element(ii).Masse_Roue=50;

Structure.Element(ii).Rayon_Roue=0.5;%Rayon de la roue

Structure.Element(ii).Epaisseur_Roue=0.2;%Epaisseur de la roue

Structure.Element(ii).Raccord_Avion=[24.1,0,d_train_atterrissage_arriere];

Structure.Element(ii).Raccord_Avion(1)=Structure.Element(numero_empennage_horiz_gauche).C_Root/2+Structure.Element(numero_empennage_horiz_gauche).Raccord_Avion(1)+sin(Structure.Element(numero_empennage_horiz_gauche).Lambda_0)*Structure.Element(ii).Raccord_Avion(3);

Structure.Element(ii).Raccord_Element=[0,0,0];

Structure.Element(ii).Rotation=[0,0,0];

Structure.Element(ii).Coeff_Lift=0;

numero_train_arriere_gauche=ii;%Variable muette pour définir l'élément Train Atterrissage Arriere Droit comme symétrique par rapport plan xy

ii=ii+1;

Structure.Element(ii).Nom='Train Atterrissage Avant Droit';

Structure.Element(ii).Type='Train Atterrissage';

Structure.Element(ii).Masse=Structure.Element(numero_train_avant_gauche).Masse;

Structure.Element(ii).GouD='D';%Train gauche ou droit

Structure.Element(ii).Rayon_Support=Structure.Element(numero_train_avant_gauche).Rayon_Support;%rayon du cylindre qui lie la roue à l'aile

Structure.Element(ii).Longueur_Support=Structure.Element(numero_train_avant_gauche).Longueur_Support;%longueur du cylindre

Structure.Element(ii).Masse_Support=Structure.Element(numero_train_avant_gauche).Masse_Support;

Structure.Element(ii).Masse_Roue=Structure.Element(numero_train_avant_gauche).Masse_Roue;

Structure.Element(ii).Rayon_Roue=Structure.Element(numero_train_avant_gauche).Rayon_Roue;%Rayon de la roue

Structure.Element(ii).Epaisseur_Roue=Structure.Element(numero_train_avant_gauche).Epaisseur_Roue;%Epaisseur de la roue

Structure.Element(ii).Raccord_Avion=Structure.Element(numero_train_avant_gauche).Raccord_Avion;

Structure.Element(ii).Raccord_Avion(3)=-Structure.Element(ii).Raccord_Avion(3); %La coordonée Z_droit = - Z_gauche

Structure.Element(ii).Raccord_Element=[0,0,0];

Structure.Element(ii).Rotation=Structure.Element(numero_train_avant_gauche).Rotation;

Structure.Element(ii).Coeff_Lift=Structure.Element(numero_train_avant_gauche).Coeff_Lift;

numero_train_avant_droit=ii;

ii=ii+1;

Structure.Element(ii).Nom='Train Atterrissage Arriere Droit';

Structure.Element(ii).Type='Train Atterrissage';

Structure.Element(ii).Masse=Structure.Element(numero_train_arriere_gauche).Masse;

Structure.Element(ii).GouD='D';%Train gauche ou droit

Structure.Element(ii).Rayon_Support=Structure.Element(numero_train_arriere_gauche).Rayon_Support;%rayon du cylindre qui lie la roue à l'aile

Structure.Element(ii).Longueur_Support=Structure.Element(numero_train_arriere_gauche).Longueur_Support;%longueur du cylindre

Structure.Element(ii).Masse_Support=Structure.Element(numero_train_arriere_gauche).Masse_Support;

Structure.Element(ii).Masse_Roue=Structure.Element(numero_train_arriere_gauche).Masse_Roue;

Structure.Element(ii).Rayon_Roue=Structure.Element(numero_train_arriere_gauche).Rayon_Roue;%Rayon de la roue

Structure.Element(ii).Epaisseur_Roue=Structure.Element(numero_train_arriere_gauche).Epaisseur_Roue;%Epaisseur de la roue

Structure.Element(ii).Raccord_Avion=Structure.Element(numero_train_arriere_gauche).Raccord_Avion;

Structure.Element(ii).Raccord_Avion(3)=-Structure.Element(ii).Raccord_Avion(3);

Structure.Element(ii).Raccord_Element=[0,0,0];

Structure.Element(ii).Rotation=Structure.Element(numero_train_arriere_gauche).Rotation;

Structure.Element(ii).Coeff_Lift=Structure.Element(numero_train_arriere_gauche).Coeff_Lift;

ii=ii+1;

else

% Ajout de la masse des trains d'atterrissage dans les ailes lors de la

% phase de vol

Structure.Element(numero_aile_gauche).Masse=Structure.Element(numero_aile_gauche).Masse+masse_train_avant;%Aile Gauche

Structure.Element(numero_aile_droite).Masse=Structure.Element(numero_aile_droite).Masse+masse_train_avant;%Aile Droite

Structure.Element(numero_empennage_horiz_gauche).Masse=Structure.Element(numero_empennage_horiz_gauche).Masse+masse_train_arriere;%Jonction Empennage Gauche Horizontal

Structure.Element(numero_empennage_horiz_droit).Masse=Structure.Element(numero_empennage_horiz_droit).Masse+masse_train_arriere;%Jonction Empennage Droit Horizontal

end
%% Assemblage Fuselage + Turbofan

Structure.Element(ii).Nom='Fuselage Bee';

Structure.Element(ii).Type='Fuselage';

Structure.Element(ii).Masse=10000;

Structure.Element(ii).Raccord_Avion=[0,0,0];

Structure.Element(ii).Raccord_Element=[0,-1,0];

Structure.Element(ii).Rotation=[0,0,0];

Structure.Element(ii).Hauteur_Cylindre_Principal=hauteur_fuselage_turbofan-hauteur_turbofan;

Structure.Element(ii).Forme_Base_Cylindre='Elliptique';

Structure.Element(ii).Longueur_Caracteristique_Base_1=1;

Structure.Element(ii).Longueur_Caracteristique_Base_2=.9;

Structure.Element(ii).Barycentre=[Structure.Element(ii).Hauteur_Cylindre_Principal/2,0,0];

numero_fuselage=ii;

ii=ii+1;
Structure.Element(ii).Nom='Turbofan';

Structure.Element(ii).Type='Moteur';

Structure.Element(ii).Masse=1000; %Valeur prise pur le moment arbitrairement

Structure.Element(ii).Raccord_Avion=[Structure.Element(numero_fuselage).Hauteur_Cylindre_Principal,0,0];

Structure.Element(ii).Rotation=[0,0,0];

Structure.Element(ii).Hauteur=2.6;

Structure.Element(ii).Rayon=0.87;

Structure.Element(ii).Raccord_Element=[0,-Structure.Element(ii).Rayon,0];

Structure.Element(ii).Barycentre=[Structure.Element(ii).Hauteur/2,0,0];

ii=ii+1;

%% Structure

%Déterminer la corde moyenne des ailes

% Utile pour calculer barycentre Bee

Structure=corde_moyenne(Structure);

%% Barycentre des masses

%Appel à la fonction bar_train(Structure) pour déterminer le barycentre

%du train d'atterrissage en prenant pour repère le bord supérieur du train

%d'atterrissage.

Structure=bar_train(Structure);
% Appel à la fonction bary_ailes(Structure) pour déterminer le barycentre

% des ailes en prennant pour origine du repère le bord supérieur gauche de

% l'aile.

Structure=bary_ailes(Structure);

% Ensuite il faut déterminer les coordonnées de ce barycentre dans le

% repère associé à l'avion. Pour cela, on a besoin de savoir où se raccorde

% l'élément à la structure et de connaître les différents angles de

% rotation.

Structure=coord_bary_repere_avion(Structure);
% Maintenant que l'on dispose de tous les barycentres des éléments, on peut

% déterminer le barycentre total de l'avion dans le repère de l'avion. On

% fait donc appel à la fonction barycentre_global(Structure).

Structure=barycentre_global(Structure);

%% Définition de la Basket en fonction du barycentre précédent
if vol_avec_basket==true

Structure.Element(ii).Nom='Basket';

Structure.Element(ii).Type='Fuselage';

Structure.Element(ii).Masse=45000;

Structure.Element(ii).Raccord_Element=[0,0,0];

Structure.Element(ii).Rotation=[0,0,0];

Structure.Element(ii).Hauteur_Cylindre_Principal=longueur_basket;

Structure.Element(ii).Forme_Base_Cylindre='Elliptique';

Structure.Element(ii).Longueur_Caracteristique_Base_1=3.5;

Structure.Element(ii).Longueur_Caracteristique_Base_2=1.5;

Structure.Element(ii).Barycentre=[Structure.Element(ii).Hauteur_Cylindre_Principal/2,0,0];

Structure.Element(ii).Raccord_Avion=[Structure.Barycentre(1)-Structure.Element(ii).Barycentre(1),-Structure.Element(ii).Longueur_Caracteristique_Base_2,0];

ii=ii+1;

end

% Moteurs

Moteurs.Element(1).nom='Moteur Droit';

Moteurs.Element(1).Type='TP400';

Moteurs.Element(1).Puissance=8*10.^6;% en W

Moteurs.Element(2).nom='Moteur Gauche';

Moteurs.Element(2).Type='TP400';

Moteurs.Element(2).Puissance=8*10.^6;% en W;
if phase_decollage==true

Moteurs.Element(3).nom='Turbofan';

Moteurs.Element(3).Type='TP400';

Moteurs.Element(3).Puissance=8*10.^6;% en W

end
end

Annexe 2 : Schéma de calcul de la fonction ‘structure_bee_plane.m’

Paramètres du Bee-Plane

Définition de la Bee

Barycentre de la Bee

Définition de la Basket

Définition des moteurs

Sorties :

• 
  Structure
  
• 
  Moteurs
  
• 
  Numéro des éléments
  

Entrées:

• 
  Avec / sans Basket
  
• 
  Phase de décollage / atterrissage ou non
  
•