Helicopter aeroelastic analysis is highly complex and multidisciplinary in nature; the flexibility of main rotor blades is coupled with aerodynamics, dynamics and control systems. A key component of an aeroelastic analysis is the vehicle trim procedure. Trim requires calculation of the main rotor and tail rotor controls and the vehicle attitude which cause the six steady forces and moments about the helicopter center of gravity to be zero. Trim simulates steady level flight of the helicopter. The trim equations are six nonlinear equations which depend on blade response and aerodynamic forcing through finite element analysis. Simulating the behavior of the helicopter in flight requires the solution of this system of nonlinear algebraic equations with unknowns being pilot controls and vehicle attitude angles. The nonlinear solution procedure is prone to slow convergence and occasional divergence causing problems in optimization and stochastic simulation studies. In this thesis, an attempt is made to efficiently solve the nonlinear equations involved in helicopter trim. Typically, nonlinear equations in mathematical physics and engineering are solved by linearizing the equations and forming various iterative procedures, then executing the numerical simulation. Helicopter aeroelasticity involves the solution of systems of nonlinear equations in a computationally expensive environment. The Newton method is typically used for the solution of these equations. Due to the expensive nature of each aeroelastic analysis iteration, Jacobian calculation at each iteration for the Newton method is not feasible for the trim problems. Thus, the Jacobian is calculated only once about the initial trim estimate and held constant thereafter. However, Jacobian modifications and updates can improve the performance of the Newton method. A comparative study is done in this thesis by incorporating different Jacobian update methods and selecting appropriate damping schemes for solving the nonlinear equations in helicopter trim. A modified Newton method with varying damping factor, Broyden rank-1 update and BFGS rank-2 update are explored using the Jacobian calculated at the initial guess. An efficient and robust approach for solving the strongly coupled nonlinear equations in helicopter trim based on the modified Newton method is developed. An appropriate initial estimate of the trim state is needed for successful helicopter trim. Typically, a guess from a simpler physical model such as a rigid blade analysis is used. However, it is interesting to study the impact of other starting points on the helicopter trim problem. In this work, an attempt is made to determine the control inputs that can have considerable effect on the convergence of trim solution in the aeroelastic analysis of helicopter rotors by investigating the basin of attraction of the nonlinear equations (set of initial guess points from which the nonlinear equations converge). It is illustrated that the three main rotor pitch controls of collective pitch, longitudinal cyclic pitch and lateral cyclic pitch have significant contribution to the convergence of the trim solution. Trajectories of the Newton iterates are shown and some ideas for accelerating the convergence of trim solution in the aeroelastic analysis of helicopter are proposed.

New helicopter rotor designs are desired that offer increased efficiency, reduced vibration, and reduced noise. This problem is multidisciplinary, requiring knowledge of structural dynamics, aerodynamics, and aeroacoustics. Rotor optimization requires achieving multiple, often conflicting objectives. There is no longer a single optimum but rather an optimal trade-off space, the Pareto Frontier. Rotor Designers in industry need methods that allow the most accurate simulation tools available to search for Pareto designs. Computer simulation and optimization of rotors have been advanced by the development of "comprehensive" rotorcraft analysis tools. These tools perform aeroelastic analysis using Computational Structural Dynamics (CSD). Though useful in optimization, these tools lack built-in high fidelity aerodynamic models. The most accurate rotor simulations utilize Computational Fluid Dynamics (CFD) coupled to the CSD of a comprehensive code, but are generally considered too time consuming where numerous simulations are required like rotor optimization. An approach is needed where high fidelity CFD/CSD simulation can be routinely used in design optimization. This thesis documents the development of physics based rotor simulation frameworks. A low fidelity model uses a comprehensive code with simplified aerodynamics. A high fidelity model uses a parallel processor capable CFD/CSD methodology. Both frameworks include an aeroacoustic simulation for prediction of noise. A synergistic process is developed that uses both frameworks together to build approximate models of important high fidelity metrics as functions of certain design variables. To test this process, a 4-bladed hingeless rotor model is used as a baseline. The design variables investigated include tip geometry and spanwise twist. Approximation models are built for high fidelity metrics related to rotor efficiency and vibration. Optimization using the approximation models found the designs having maximum rotor efficiency and minimum vibration. Various Pareto generation methods are used to find frontier designs between these two anchor designs. The Pareto anchors are tested in the high fidelity simulation and shown to be good designs, providing evidence that the process has merit. Ultimately, this process can be utilized by industry rotor designers with their existing tools to bring high fidelity analysis into the preliminary design stage of rotors.

No description available.

L’une des principales sources d’inconfort dans un hélicoptère sont les vibrations transmises par le rotor à la structure de l’appareil. En vol d’avancement, des efforts aérodynamiques cycliques sont subis par l’ensemble des pales en tête rotor et génèrent de très fortes vibrations basse fréquence (aux alentours des 17Hz) transmises aux passagers via la boîte de transmission principale puis le fuselage lui-même. Afin de garantir le confort des membres d’équipage et des passagers, de nombreux systèmes antivibratoires ont été conçus. Ces systèmes sont généralement passifs car la majorité de l’énergie vibratoire transmise à la structure se situe à une fréquence unique ωc correspondant à bΩ avec b le nombre de pales et Ω la fréquence de rotation du rotor. Cependant, les appareils modernes évoluent et le régime rotor jusqu’alors fixe durant toutes les phases de vol varie à présent pour des préoccupations de performances et de consommation (variation de l’ordre de +/-10% autour de bΩ). Cette nouvelle contrainte dans la conception des hélicoptères rend pertinente la technologie des systèmes antivibratoires actifs, pouvant s’adapter à la sollicitation en termes d’amplitude et fréquence. Lors de ces travaux de thèse, la suspension passive SARIB de Airbus Helicopters basée sur le principe du DAVI (Dynamic Antiresonant Vibration Isolator) est modifiée afin d’être rendue active par ajout d’une partie actuation/commande. La théorie des lois et algorithmes de contrôle utilisés dans ces travaux, est présentée en détail afin de poser solidement les bases du contrôle actif du prototype de suspension conceptualisé ici à savoir le contrôle FXLMS (adaptatif) et le contrôle optimal LQG. Afin de simuler le fonctionnement du système, un modèle tridimensionnel de la suspension active est construit, couplé à la structure souple de l’hélicoptère (NH90). Sur ce modèle sont alors appliquées les différentes lois de commande introduites auparavant et leurs performances comparées dans différents cas de chargement en tête rotor et surtout pour différentes fréquences de sollicitation. De même, pour chaque algorithme, différentes localisations des capteurs d’erreur sont étudiées afin de converger vers une configuration optimale. Les simulations démontrent que l’algorithme FXLMS feedforward est très bien adapté au contrôle des perturbations harmoniques et permet de réduire très significativement le niveau vibratoire du plancher cabine, sans réinjection parasite dans le reste de la structure. Une comparaison de l’efficacité du SARIB actif avec les systèmes d’absorbeurs en cabine est ensuite effectuée pour démontrer la pertinence d’utiliser le principe du DAVI comme base d’un système actif. Les travaux de cette thèse traitent également des essais réalisés en laboratoire sur le prototype échelle 1 de la suspension SARIB active avec contrôle FXLMS. / One of the main causes of discomfort in helicopters are the vibrations transmitted from the rotor to the structure. In forward flight, the blades are submitted to cyclic aerodynamic loads which generate low frequency (around 17Hz) but high energy mechanical vibrations. These vibrations are transmitted from the rotor to the main gearbox, then to the structure and finally to the crew and passengers. In order to maintain acceptable comfort for crew members and passengers, a lot of antivibration devices have been developed since the last 30 years. These systems are generally passive because most of the mechanical energy transmitted to the structure is at only one frequency ωc which is equal to the product bΩ with b the number of blades and Ω the rotor rotational speed. However, modern helicopters evolve and the rotor rpm, which has always been considered as fixed during flight is now a function of time, depending on the flight phases in order to increase performances and reduce energy consumption (variation bandwidth of Ω +/- 10%). This new constraint on the design of helicopters makes the active antivibration technology completely relevant with its capacity to adapt in terms of amplitude and frequency to the perturbation. During this thesis, the passive suspension called SARIB from Airbus Helicopters, based on the DAVI principle (Dynamic Antiresonant Vibration Isolator) is modified in order to implement active components and command (actuation). The theory of the control algorithms used in this thesis is presented in detail in order to define the theoretical tools of the active DAVI control which are : FXLMS control (adaptive control) and LQG (optimal control). To simulate the complete system, a 3D multibody model of the active suspension has been set up, coupled to a the flexible structure of a NH90 (Airbus Helicopters). On this model are applied the different control algorithms presented before and their performances are compared for different loads with variable frequency on the rotor hub. In the same way, different locations for the error sensors in the structure are studied to find the optimal control configuration. The simulations show that the FXLMS algorithm is well suited for the control of harmonic perturbations and reduce significantly the dynamic acceleration level on the cabin floor, without parasite reinjection on other parts of the structure. A comparison of the active SARIB with classical cabin vibration absorbers is also made in terms of efficiency in order to show the advantages of using the DAVI system as a base for an active antivibration device. Finally, this thesis also presents the experiments realized in the dynamics laboratory of Airbus Helicopters on a 1:1 scale prototype of the active SARIB suspension with FXLMS control. The results demonstrate the efficiency of the active suspension architecture and control algorithms.

This research investigates on-blade partial span active plain trailing edge flaps (TEFs)with an aim to alleviate the helicopter vibrations. Among all the available smart materials, piezoelectric stack actuator(PEA)has shown its strong candidature for full scale rotor systems. Although, PEAs are quite robust in operation, however, they exhibit rate dependent hysteresis phenomenon and can generate only very small displacements. Dynamic hysteresis is a complex phenomenon which, if not modeled, can lead to drift in the vibration predictions. In this research, a comprehensive experimental analysis is performed on a commercially available piezostack actuator, APA-500L, which is well suited for full scale applications. Rate dependent hysteresis loops are obtained for helicopter operational frequencies. Nonlinear rate-dependent hysteresis loops are modeled using conic section approach and the results are validated with experimental data. Dynamic hysteresis exhibited by the PEA is further cascaded with the helicopter aeroelastic analysis and its effect on helicopter vibration predictions is investigated. PEAs generate high force but are limited by small translational motions. A linear to rotary motion amplification mechanism is required to actuate the TEF for vibration alleviation. A smart flap is designed and developed using computer-aided-design models. A rotor blade test section is fabricated and a lever-fulcrum mechanism (AM-1) is developed for a feasibility study. Smart flap actuation is demonstrated on the rotor blade test section. The conventional motion amplification devices contain several linkages, which are potential sites for structural failure. A novel pinned-pinned post-buckled beam linear-to-rotary motion amplifier (AM-2) is designed and developed to actuate the flaps. A new design of linear-to-linear amplification mechanism (LX-4) is developed and is employed in conjunction with AM-2 to increase the flap angles by an order of magnitude. An analytical model is developed using Mathieu-Hill type differential equations. Static and dynamic tests are conducted on a scaled flap model. Helicopter aeroelastic simulations show substantial reduction in hub loads using AM-2 mechanism. To further enhance the flap angles, an optimization study is performed and optimal beam dimensions are obtained. A new technique is also proposed to actively bias the flaps for both upward and downward motion. Critical flap design parameters, such as flap span, flap chord and flap location influences the flap power requirement and vibration objective function significantly. A comprehensive parametric investigation is performed to obtain the best design of TEFs at various advance ratios. Although, parametric study equips the designer with vital information about various critical system parameters, however, it is a computationally expensive exercise especially when used with large comprehensive helicopter aero elastic codes. A formal optimization procedure is employed to obtain the optimal flap design and location. Surrogate models are developed using design of experiments based on response surface methodology. Two new orthogonal arrays are proposed to construct the second order polynomial response surfaces. Pareto analysis is employed in conjunction with a newly developed computationally efficient evolutionary multi-objective bat algorithm. Optimal flap design and flap locations for dual trailing edge flaps are obtained for mutually conflicting objectives of minimum vibration levels and minimum power requirement to actuate the flaps.