Modélisation du comportement dynamique des rotors d'hélicoptèresL'objectif de la thèse est l'étude et la modélisation du comportement dynamique et aérodynamique du rotor de l'hélicoptère en considérant à la fois les forces d'inertie et les forces aérodynamiques et en tenant compte des déformations élastiques des pales. L'algorithme de couplage proposé permet d'effectuer le calcul transitoire avec échange de données entre les solveurs fluide et structure à chaque pas de temps.La particularité de cette étude est l'utilisation du modèle aérodynamique de la ligne active, qui représente les forces de pale appliquées au fluide par des termes sources. Ces termes sources sont répartis dans les cellules de maillage à l'emplacement de la pale. Ainsi, la rotation, la torsion et le battement de la pale peuvent être représentés sans aucune déformation du maillage. Un avantage de la ligne active est que la simulation utilise un nombre réduit de nœuds, car des conditions aux limites «lois des parois» ne doivent pas être modélisées.Le cas d'un petit rotor expérimental d'hélicoptère est étudié en vol d'avancement. Les solveurs de fluide et de structure sont couplés pour calculer le comportement aérodynamique et dynamique du rotor. Pour ce faire, un algorithme de couplage faible en série décalé est appliqué. Les calculs itératifs sont contrôlés par un code spécialement conçu. Au début de chaque itération, le code calcule et répartit les termes sources dans le domaine fluide. A la fin du pas de temps, le code exécute le solveur de calcul de structure pour calculer un seul pas de temps. Ce solveur calcule le déplacement de la pale sous l'effet des forces aérodynamiques, élastiques et d'inertie et renvoi les résultats au solveur fluide. Les déplacements de la pale calculés servent de référence pour le solveur fluide au pas de temps suivant, pour distribuer les termes sources. Le calcul s'arrête lorsque le critère de convergence est vérifié.Afin de valider le cas simulé, des expérimentations sont réalisées en soufflerie. La puissance et la poussée aérodynamique du rotor sont mesurées. La Vélocimétrie par images de particules (PIV) est utilisée pour obtenir le champ de vitesse autour du rotor. Les mesures PIV à phase bloqué dans des plans azimutaux ont permis de reconstituer le champ d'écoulement 3D. La comparaison entre les résultats numériques et les expériences montre un bon accord et permet de valider la méthode de couplage proposée. / MODELING THE DYNAMIC BEHAVIOR OF HELICOPTER ROTORThe aim of the thesis is the investigation and modeling of dynamic and aerodynamic behavior of helicopter rotor considering both inertial and aerodynamic forces and taking into account the elastic deformation of the blades. The proposed coupling algorithm allows the transient calculations with data exchange between the fluid and structure solvers at each time step.The particularity of this research is the use of an actuator line aerodynamic model, which represents the blade forces applied to the fluid as source terms. These source terms are distributed in the grid cells where the blade is located. Thus the rotation, flapping and torsion of the blade can be represented without any grid deformation. An advantage of the actuator line is that the simulation uses a reduced number of nodes, because the “wall” boundary conditions do not need to be modeled.The case of small experimental helicopter rotor is studied in forward flight. The fluid and structure solvers are coupled to calculate aerodynamic and dynamic behavior of the rotor. For this purpose, a loosely coupling serial staggered algorithm is applied. The iterative calculations are controlled by specially developed code. At the beginning of each iteration, this code calculates and distributes the source terms in the fluid domain. At the end of the time step, the code runs the structural solver to execute a single time step. This solver calculates the blade displacement under aerodynamic, elastic and inertial forces, and the results are returned to the fluid solver. The calculated blade displacements serve as reference in the next fluid step to distribute the source terms. The calculation stops when the convergence criteria are met.In order to validate the simulated case, measurements are carried on in the wind tunnel. The power and aerodynamic thrust of the rotor are measured. Particle Image Velocimetry (PIV) is used to obtain the velocity field around the rotor. Phase locked measurement in azimuth planes enabled to reconstruct 3D flow field. The comparison between numerical results and experiments shows good agreement and permits to validate the proposed coupling method.
| Identifer | oai:union.ndltd.org:theses.fr/2013ENAM0055 |
| Date | 17 October 2013 |
| Creators | Velkova, Cvetelina Vladimirova |
| Contributors | Paris, ENSAM, Massouh, Fawaz, Todorov, Michael, Dobrev, Ivan |
| Source Sets | Dépôt national des thèses électroniques françaises |
| Language | French |
| Detected Language | French |
| Type | Electronic Thesis or Dissertation, Text |
Page generated in 0.0022 seconds