Next Article in Journal
The Effects of Sn Doping MnNiFeO4 NTC Ceramic: Preparation, Microstructure and Electrical Properties
Previous Article in Journal
Research, Modelling and Prediction of the Influence of Technological Parameters on the Selected 3D Roughness Parameters, as Well as Temperature, Shape and Geometry of Chips in Milling AZ91D Alloy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Helicopter Main Rotor Blade Parametric Design for a Preliminary Aerodynamic Analysis Supported by CFD or Panel Method

1
Doctoral School, Military University of Technology, ul. gen. Sylwestra Kaliskiego 2, 00-908 Warsaw, Poland
2
Faculty of Mechatronics, Armament and Aerospace, Military University of Technology, ul. gen. Sylwestra Kaliskiego 2, 00-908 Warsaw, Poland
*
Author to whom correspondence should be addressed.
Materials 2022, 15(12), 4275; https://doi.org/10.3390/ma15124275
Submission received: 31 March 2022 / Revised: 2 June 2022 / Accepted: 13 June 2022 / Published: 16 June 2022

Abstract

:
This work is the preliminary part of a research program which is aimed at finding some new methods and design solutions for helicopter main rotor multidisciplinary optimization. The task was to develop a parametric geometric model of a single-blade main rotor applicable for varied methods of numerical aerodynamic modeling. The general analytical assumptions for the parametric main rotor design were described. The description of the main rotor blade parametric design method based on Open GRIP graphical programming was presented. Then, the parametric model of a blade was used for aerodynamic models independently developed for panel method and advanced CFD solver. The results obtained from the CFD simulations and panel analysis for main rotor aerodynamics were compared and assessed using analytical calculations. The calculations and simulations for a single-blade and completed rotor were performed for different helicopter weights and rotor pitch angles. The results of different computer aerodynamic analysis environments were compared for the possibility of their application in an optimization loop. This is preliminary work that describes only a partial problem that could be used in the future as part of a comprehensive methodology for aerodynamic and structural optimization of a helicopter rotor. As an output of the research, new options for main rotor optimization are developed. The combined parametric modeling with aerodynamic analysis, as described in this paper, provide the preliminary design for a main rotor spiral, as an element of the optimization loop.

1. Introduction

Military rotorcrafts that are currently operating are mostly constructions that were designed in the second half of the twentieth century. The constructions (in the preliminary and detail design phases) were prepared with classic mathematical calculations, eventually supported with CAD. With improved computerized options, design processes have been fully transferred to virtual environments. New aircraft structures are fully designed using computer methods; however, the structures of the existing constructions have also been modeled using engineering software, which is described in [1].
The most popular engineering design solution using CAD is a point and click method. A more comprehensive design method is to prepare a model that is generated as a result of entered parameters, which can be obtained by using scripting in the environment’s language; an example of parametric blade optimization is shown for a wind turbine blade in [2], and for a main rotor blade in [3]. The parametric design includes options to quickly change the model’s features without time-consuming 3D modeling. It is important to note that the influence of parametric change can be rapidly checked in numerous variants.
The helicopter rotor blade is an element which can be successfully adapted to parametric design, because of mathematical functions which precisely describe blade features. The parametric model is the first step to prepare the main rotor optimization process [4]. Some examples have been shown in different works: parametric airfoil design [5,6,7], structural design [8,9], and aerodynamic design [10]. A good parametric analysis of the rotor blade framework is presented in [11]. In this research, the parametric model was prepared using GRIP modeling language, which is an integral part of the Siemens NX environment. The usage of GRIP in parametric modeling was conducted in [12,13,14].
The proposed graphical programming language offers the possibility to prepare a geometric and structural model of the aircraft element. The designed part is fully described by the program script. In addition, the elements’ features can be calculated from mathematical functions, for example, polynomials for shape parametrization. After programming the structure, Open GRIP offers the possibility to analyze the inertia physiognomy of the body. The obtained data are the basis for the strength calculation. The program also grants a function for preparing forms which can be filled by the user. Therefore, some parameters can be changed during the design loop, without changing the program code. In the future, it is crucial to prepare a full parametric optimization. An effective approach to main rotor blade optimization was proposed in [8,10,15,16,17,18], and effective optimization of rotor construction was proposed in [19].
This work is part of a research program that is aimed at finding the best rotorcraft construction optimization solutions. The first step of the research program was an analysis of modern helicopter constructions and an outline of design parameters that are crucial in helicopter main rotor design. The results of the analysis showed that there is a field were rotor construction can be improved in accordance with modern combat field requirements which are dynamically changing due to fast technological improvements. The results, with charts and a comparison of parameters were published in [20].
On the basis of construction comparisons, we conducted a study of new design methods or combinations of some separate solutions. The aim of the study was to provide a versatile, accurate main rotor model that could be analyzed, validated, as well as quickly and easily changed in the optimization process. As a new approach to rotorcraft design, we combined the parametric programming CAD model generator with the CFD and panel method, with the goal of providing new possibilities for improving the construction optimization procedure. The effects and possibilities of CFD main rotor studies that were implemented in the main rotor development procedure were well presented in [21,22,23]. In addition, a CFD analysis, as a base for a strength analysis because of load prediction, was clearly shown in [24].
The development of the method started with the main rotor blade, as the main element that generates thrust. An existing rotor blade model was taken as a basis to conduct and authenticate the method. As previously described in the Introduction, the blade model was prepared using the parametric method in GRIP language. GRIP language was chosen because it is user-friendly graphical programming code for easily preparing the model with extensive built-in possibilities for quickly calculating the inertia parameters of the generated model, which could be used in further optimization steps, which is crucial when the strength of the construction is taken into consideration in the optimization loop. Further, the model was validated using the CFD and panel methods, and the thrust results were compared with the analytical calculations. The methods were also compared with each other to choose the appropriate procedure to perform the most efficient optimization. To the best of our knowledge, similar studies that combined the proposed methods for rotorcraft main rotor design and main rotor optimization procedure have not been conducted.
This paper is organized as follows: The research methods are described in Section 4, where the mathematical model, which is the basis for the parametric programming, is defined, and in the second part of this section, the CFD and panel method settings are given; the validation results are described in Section 3; and the results are discussed and concluded in Section 4.

2. Research Methods

2.1. Parametric Description of a Single-Blade Rotor with Aspect of Its Geometrics and Aerodynamics

The construction process of helicopter rotor blades requires consideration of fatigue, strength, stiffness, cost, and vibrations. The rotor blade works in a changeable environment, and therefore, the designer has to take into consideration all of the parameters that transform during flight. Periodic changes in the flow velocity on the blade section can modify the angle of attack and the Mach number, and as a consequence, can change the aerodynamic coefficients. Blade turns (relative to hinge axis) and blade deformation cause a periodic change in inertial forces and provoke a cyclic modification in the angles of attack and flow velocity. The flexibility of the main rotor blades and their work in a strong centrifugal force field induce an inseparable connection between loads and deformations. An analysis of the problems mentioned and also the noise level or flow disturbance requires a versatile evaluation of the main rotor blade’s work conditions. Effective blade aerodynamic modeling was proposed in [25].
The precise blade geometric design is crucial for obtaining the required main rotor power. As described in [26], the geometric parameters of the blade are the following:
-
Radius The length of the blade measured from the axis of rotation to the tip;
-
Chord The length of the blade measured from the airfoil leading edge to the trailing edge and, for a tapered blade, it is the function of local radius;
-
Airfoil The cross-sectional shape that determines the aerodynamic parameters of the blade, i.e., lift, drag and momentum coefficients;
-
Contour shape This is the final shape of the blade, which depends on the chord function and tip shape;
-
Geometric twist Variation of the airfoil angle between the chord and a plane of rotation along the radius;
-
Aerodynamic twist Variation of the airfoil shape along the radius;
-
Position and shape of trim tabs Defines the possibility to adjust the rotating blade to the plane of rotation.
The area of a blade root airfoil transfers the whole load from the rotating aerodynamic surface to the blade grip sticking out of a rotor head block. The shape of the root is dictated by the build conception and the adopted aerodynamic solutions. The geometric twist is usually settled on several degrees. The shape of the rotor tip depends on the aerodynamic problems, for example, achieving the speed of sound or the noise level. An analysis of rotor blade tips was described in [27]. According to the mentioned work, there are three main types of helicopter tip designs: BERP tip, the parabolic tip, and the swept (tapered) tip. An example of a main rotor blade is shown in Figure 1.
To evaluate the prepared parametric model, a theoretical lift calculation must be conducted. The lift can be calculated in vertical flight using the blade element theory combined with the momentum theory. The first assumptions are the mass conservation equation for rotor disc:
m ˙ = ρ A v i
and the momentum conservation equation for hoovering rotor:
T = m ˙ w
The energy conservation equation is the work of the rotor to change the rate of energy in fluid:
T v i = 1 2 m ˙ w 2
The calculation is based on the thrust coefficient, as a function of pitch angle and inflow angle, and is given [28,29]:
C T = 0 1 a σ 2 ( ϑ δ ) r 2 d r
where a is the lift curve slope, σ is the rotor solidity, ϑ is the pitch angle, δ is the inflow angle, and r is the current radius.
According to Figure 2, the inflow angle δ can be defined by:
δ = a r c t g ( W + v i Ω r )
where W is the vertical velocity, v i   is the induced velocity, and Ω is the rotor angular velocity.
To calculate the induced velocity, for required thrust in hover, the equation solved from the momentum theory is given:
v i = T 2 ρ A
where T is the rotor thrust, ρ is the air density, and A is the rotor planform area. For hover, the W velocity equals 0.
The inflow angle can be defined with the inflow ratio as λ i = δ / r , therefore, the thrust coefficient will transform to:
C T = 0 1 a σ 2 ( ϑ r 2 λ r ) d r
Rotor performance can be calculated with the combination of the blade element theory with the momentum theory. Therefore, the induced velocity can be determined for nonhomogeneous inflow distribution, by using the differential form of the momentum theory:
d C T = 4 λ i λ r d r
with the blade element theory C T equation form.
As a result, for hover where λ c = 0 , the inflow ratio for the induced velocity is:
λ i = a σ 16 [ 1 + 32 a σ ϑ r 1 ]
The σ rotor solidity equation is necessary, which, in this work, is nonuniform along the blade span as a consequence of non constant chord. The rotor solidity is calculated as:
σ = 0 1 b c ( r ) π R d r
where b is the number of blades, c(r) is the chord in function of radius, and R is the rotor radius.
The thrust coefficient is the function of radial location of the blade measured from the centre of rotation to the blade tip. It is calculated as an integral with the r limits from 0 to 1. However, to obtain better quality results, a tip loss factor can be proposed [30]:
B = 1 a c ( r = 1 ) 2 R
The tip loss factor is taken as an upper limit of the coefficient integral. The thrust losses are also a result of the root cutout. It is usually 10% to 30% of the blade radius. Including tip loss factor and root cutout, the thrust integration is:
T = 1 2 ρ a b Ω 2 R 2 r R B c ( ϑ δ ) r 2 d r
According to [22], the mean lift coefficient for the main rotor blade can be determined:
C ¯ L = 6 C T σ
The calculations described above were all prepared in MATLAB. The estimated mathematical thrust model is the basis for the CFD computation evaluation.
The program build is defined using the algorithm presented in Figure 3. The code is constructed as a loop to look for the collective angle required for a given aircraft weight.
It starts by defining the input parameters that can be obtained from first cut calculations and the mission flight conditions. In the beginning, the collective is positioned at 0° angle. Next, the position radius coordinates are defined for the twist and chord values in the given sections. The thrust is a radius position integral, and therefore, all program functions are defined as dependent on x. Using the chord and twist values, the chord and twist polynomials are interpolated. These polynomials are also used in the GRIP program to create the blade geometrics. Next, the inflow angle, the rotor solidity, and tip loss factor are given. All the specified equations are inputs to the thrust integral which is bounded by the root cutout and the tip loss factor. The blade attack angle is calculated using the pitch angle with addition of collective. The thrust equation is computed using loop. The collective angle is increased in each step until the result is similar to the hover required force for the assumed mass. The determined collective is used for the CFD analysis to check the parametric model aerodynamic parameters. The procedure described above can be applied for the hover case calculations and for vertical or horizontal flight conditions as well.
To present the possibilities of Matlab calculations for future main rotor optimization, a multi-dimensional analysis was conducted. The main rotor thrusts were calculated for different numbers of blades, twist angles, and mean chord values. The computations were made for the collective angle 7.47°. The distribution of thrust values as a function of twist angle and chord value are shown in Figure 4, Figure 5 and Figure 6. As predicted, the highest thrust value is obtained for a five-blade rotor, although the reference system is a four-blade rotor, corresponding to the rotor design of the Polish W-3 Sokol helicopter.

2.2. Modeling the Blade Geometrics and Aerodynamics Using Proprietary Software Applications

The calculations of the blade were performed for an existing main rotor blade of a Polish military W-3 SOKOL helicopter. The model could be easily modified with different parameters. The technical specifications are shown in Table 1, however, some details such as the blade and tip shape (which are provided from design documentation) cannot be published. The simulation was performed for the aerodynamic calculations, therefore, it was modeled as an empty shell only to provide the blade shape.
Figure 7 shows the original algorithm dedicated to the rotor blade parametric modeling. The rotor blade parametric models were prepared in Siemens NX Open Grip language. It is a programming language that can draw a CAD model with commands and is able to execute advanced customized operations in a more effective way than interactive NX operations. Numerous interactive operations can be executed using GRIP language. It is possible to draw geometric objects using control of parameters or load data from a file to modify the geometrics. Exemplary program code for parametric modeling of the blade geometrics is presented in Listing 1.
Listing 1. Part of GRIP code with blade surface generation.
….
SPL3(K+1)=SPLINE/PT0(1..N-1) $$ upper surface
SPL4(K+1)=SPLINE/PT1(1..N-1) $$ lower surface
dfi=fi(4)*cr+fi(3)*cr+fi(2)*cr+fi(1)

MAT=MATRIX/XYROT,-dfi
SPL1(K+1)=TRANSF/MAT,SPL3(K+1)
SPL2(K+1)=TRANSF/MAT,SPL4(K+1)
$$ LN(K+1)=TRANSF/MAT,LN1(K+1)

PT2(K+1)=POINT/ENDOF,XSMALL,SPL1(K+1)
PT3(K+1)=POINT/ENDOF,XLARGE,SPL1(K+1)
PT4(K+1)=POINT/ENDOF,XLARGE,SPL2(K+1)
LN(K+1)=SPLINE/PT3(K+1),PT4(K+1)

DELETE/PT0,PT1,SPL3,SPL4,LN1

IFTHEN/K<M
K=K+1
JUMP/L10:
ELSE
JUMP/L40:
ENDIF

L40:
A=&POINT(PT3(1))
B=&POINT(PT4(1))

SPLC(1)=SPLINE/PT2(1..K+1)
SPLC(2)=SPLINE/PT3(1..K+1)
SPLC(3)=SPLINE/PT4(1..K+1)
SSRF(1)=BSURF/MESH,SPLC(1..2),WITH,SPL1(1..K+1),TYPE,3,TOLER,.01,.01
SSRF(2)=BSURF/MESH,SPLC(1),SPLC(3),WITH,SPL2(1..K+1),TYPE,3,TOLER,.01,.01

SSRF(3) = RLDSRF/SPLC(2),,SPLC(3)

SSRF(4) =RLDSRF/SPL1(1),,SPL2(1)
SSRF(5) =RLDSRF/SPL1(K+1),,SPL2(K+1)
BLD(1)=SEW/SSRF(1..5)
….
The program prepared for generating the main rotor blade was easy to modify, and the blade parameters could be changed without needing to repeat an interactive method. The blade radius, chord, chord distribution, twist, and airfoil shape could be easily changed to generate a model with new features. The process of parametric model generation is presented in Figure 8 and Figure 9.

2.3. Aerodynamic Modeling of a Single-Blade and Complete Main Rotor

The generated model was implemented in two computer environments to evaluate the properties of the obtained blade shape. The main rotor blade was implemented in a CFD software, i.e., ANSYS Fluent and an aerodynamic modeling software, i.e., DARcorporation Flightstream. The domain for CFD was prepared in ANSYS Spaceclaim. The air density in the simulations was 1.25 kg/m3 and it was constant.
The implementation was conducted to assess the model and to prepare its applicability in further research; however, it was also conducted to compare two different tools for aerodynamic rotorcraft analysis with mathematical calculation. With the results of the CFD analysis, a further strength analysis is possible. In addition, the CFD analysis is a crucial step in the optimization process of a modern main rotor blade design.

2.3.1. One Blade CFD Modeling in ANSYS Fluent

In ANSYS Fluent, a one blade CFD model was prepared. The fluid domain was prepared using Ansys Spaceclaim. It was generated with a quarter of a sphere and half of a cylinder. The air enclosure were 20 m long and high, and the width of the enclosure was 10 m. A hybrid mesh was built for the computational area using the ANSYS Mesh module. Body and face sizing were used with inflation to generate the correct mesh. The boundary layer was modeled using the full thickness option. It was set at 20 mm for 25 layers with a growth ratio of 1.2. The obtained y+ value was 8. The maximum element size was set at 1000 mm, with a size reduction closer to the blade of 20 mm. The 3D bodies were transformed into a tetrahydra mesh with prismatic components within the boundary layer. The mesh consisted of almost 2,837,051 elements. The mesh is shown in Figure 10.
The boundary conditions are presented in Figure 11. The blade is fixed in the enclosure which is modeled by combining half of a hemisphere with a cylinder. In addition, to obtain correct body sizing, a smaller cuboid enclosure was made around the blade, and by naming the object faces, an air flow direction was projected.
To imitate the main rotor blade working environment, velocity change, pitch, and inflow angle were implemented with mathematical expressions. A velocity change was realized by an expression that raised the inflow speed with the span. To model a change in the inflow angle, the inflow directions on the X and Y axes were expressed using the cosine of the X axis and sine of the Y axis. The inflow angle was calculated from the inflow ratio with the Matlab program. Table 2 presents the inflow input function to model the air flow over the blade. The polynomial was a 4th grade polynomial, the higher grades gave similar results. The functions values changed with the span.
The viscosity model which was used for simulations was k-omega SST. It was applied because it is recommended for rotating machinery.

2.3.2. One Blade Panel Modeling in FligtStream

In the second aerodynamic environment, a main rotor blade model was also prepared. The parametric blade model was easily implemented. Software was used to import the CAD files from commercial software, and therefore, the tool could be used in the optimization loop.
In the Flightstream, the mesh is generated only on the studied element. It does not required from user preparing an enclosure and environment. The mesh is generated automatically by the software. Because of the simplicity of the software, the inflow angle and the change of the pitch angle were modeled with the blade position in the reference frame. Velocity change was simulated by the build “shear” freestream option and prepared input .txt file with the velocity magnitude across the span. The advantage of the solver is the fact that the model for conducting a simulation is ready to use three–four times faster than in the Fluent environment. However, rotorcrafts are complex constructions that work in versatile conditions, and therefore, some working states may not be imitated. The solver was setup with steady viscous parameters. The solver setup is shown in Figure 12.

2.3.3. Main Rotor CFD Modeling in ANSYS Fluent

In the second stage of the research, the blades were multiplied and positioned in a rotational plane. The fluid domain were also modeled in Spaceclaim. A mesh was built for the computational area using the Mesh Ansys module. The boundary conditions are presented in Figure 12. Body sizing and inflation were used to generate the correct mesh. The boundary layer was modeled using the full thickness option. The air enclosure was 75 m long and wide, and the height of the enclosure was 25 m. The maximum element size was set at 2000 mm, with a size reduction closer to the blade of 20 mm. The boundary layer was modeled using the full thickness option. It was set at 20 mm for 25 layers with a growth ratio of 1.2. The obtained y+ value was 8. The 3D bodies were transformed into a triangular mesh. The mesh consisted of nearly 5,300,000 elements.
The blades were placed in the enclosure which was modeled with cuboid. In addition, to obtain correct body sizing, a smaller cylinder enclosure was made around the blades. By naming the object faces, an air flow direction was projected. The analysis was performed with mesh motion, and therefore, there was no inlet, because the air flow was provoke by the rotating blades. The boundary conditions are presented on Figure 13.
The blades were set at the pitch angle calculated for the required thrust. The models were set to rotate with the main rotor angular speed. The inflow angle was assumed to be generated in the solver with the rotational movement. The viscosity model for the simulation that produced the most accurate results with the lowest time consumed was realizable k-epsilon with scalable wall function.

2.3.4. Main Rotor Panel Modeling in FligtStream

Since the modeling was performed using CFD software, the panel environment was also tested to provide results for parametric modeling in the optimization procedure. The model building was more time-consuming and it was easy to make a mistake, because the pitch angle and inflow angle were modeled by the blade’s position. The blades were multiplied and positioned in a rotational plane. The mesh was generated the same way as it was generated for the one blade case. The cases were calculated with an unsteady solver. The rotational movement was modeled with motion settings, where the angular velocity and rotation axis could be defined.

3. Results

The evaluation results are shown in Table 1 and Table 2. In both cases, the parametric model has no issues when implemented. Calculations of the aerodynamic thrust were conducted for selected design weight values corresponding to those from the range of operating weights of the W-3 Sokol helicopter (i.e., between 4900 and 6400 kg).
Table 3 shows the comparison of the results of Fluent and Flightstream stationary blade modeling with mathematical calculations for a given thrust and the required collective angle for that thrust. The collective angles were calculated using thrust Equation (9) using Matlab procedure, combining the blade element theory with the momentum theory. According to the obtained results, shown in Figure 14, CFD modeling resulted in a high accuracy of expected resultant forces; both solvers produced similar results. However, for lower thrust, Flightstream, which is mainly designated to steady simulations, delivered a more accurate result.
During the preparation of the stationary simulation, it is crucial to model the exact inflow conditions. The main advantage of this type of imitation is the short time for the solver calculations. As compared with a four blade rotary model, it provides the results about 9–10 times faster. It also offers the possibility of checking the different collective values and inflow angle using one model and solver settings. This type of simulation is a good choice for evaluating the first cut blade loads.
Table 4 presents the comparison of results from Fluent and Flightstream for the main rotor model with modeled rotational movement. The mathematical assumptions are the same as in the steady model. The results presented in Figure 15 show that, for the rotary model, definitely better results are obtained with the Fluent solver. The Flightstream environment gives the most accurate findings with the unsteady solver, nevertheless, they is a 20% difference btween the mathematical and CFD solver model. The advantage of the Flighstream simulation is the duration of calculation. For the Ansys software, the calculations with mesh preparation last about 3 h (for a 30 threads computer station); with the DARcorporation solver, the simulation with meshing is ready in up to 30 min (with 8 threads). Despite the results, the Flightstream simulation can be used, when a designer needs to take a first look at the rotor aerodynamics behavior. However, to obtain accurate findings, a CFD simulation needs to be conducted.
The rotary model provides the calculations and forces that can be used in further analysis. It can be applied to the strength model and can calculate the required mass parameters to sustain the blade loads. Combining the models’ mass loads from rotary motion can also be simulated.

4. Conclusions

As stated in the Introduction, this work is part of an aerospace construction optimization research program. It is the second step in preparing the optimization loop for rotorcraft constructions.
In this research, it was confirmed that graphical programming language with a user-friendly syntax is a good solution for this type of geometric modeling. The commands are intuitive and easy to use. The code is short and the blade model is generated within a few seconds. In addition, the language offers an inertia analysis for solids, which is planned to be used in future research on the optimal blade design. The generated model was versatile and easy to adjust in different environments. There were no compatibility problems with the CFD solvers.
The parametric model behaved as it was assumed. The CFD analysis confirmed that the generated shape had the aerodynamic features that were preferred during the design phase. The mathematical assumptions of the inflow angle were confirmed with the rotary CFD model.
This research shows a new approach to rotor blade inflow modeling. The stationary blade with complex air flow analytical model gives correct results and consumes less computer power and time. Therefore, it is possibile to implement this type of simulation in the preliminary design of a main rotor, which is useful in optimization studies. Combining the stationary blade with parametric modeling is the first stage of optimization to evaluate a construction and to enhance the first cut shape. This method provides an initial shape in a short period of time. Therefore, even when the calculations are combined with mechanical analysis (FSI), the analysis time is reduced. Finally, the iterative procedure for optimization studies is proposed in Figure 16. In the first phase, the CFD model can be replaced with a panel code environment. To summarize, the results of this study can serve as a basis for developing a main rotor aerostructural model (combining aerodynamic properties and finite element structure) for further strength analysis, where both the aerodynamic and weight loads can be considered. In future work, in accordance with the research program, the CFD and FE model will be combined using scripting methods to prepare an optimization loop to find a design shape that fulfils modern combat field requirements and offers the best mass properties.

Author Contributions

Conceptualization, J.K., R.R. and S.K.; methodology, J.K.; software, J.K.; validation, R.R. and J.K.; formal analysis, J.K.; writing—original draft preparation, J.K.; writing—review and editing, R.R.; visualization, J.K.; supervision, S.K.; project administration, S.K. All authors have read and agreed to the published version of the manuscript.

Funding

The methods and results presented in this paper have been obtained during research works conducted within the university research project entitled “Methods of optimal design of aircraft and helicopters in application to improve their design and aerodynamic properties”. This work was financed by the Military University of Technology (Warsaw, PL), in 2022, under research project UGB-782.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Program code is available at: https://github.com/jakubkocjan/parametricrotordesign (accessed on 30 March 2022).

Acknowledgments

The present work was conducted as part of the Aerospace Constructions research program. The authors gratefully acknowledge their support.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

AMain rotor disc area
aMain rotor blade lift curve slope
BTip loss factor
bNumber of blades
C L Lift coefficient
C T Thrust coefficient
cBlade chord
m ˙ Mass flux through the rotor disc
RMain rotor blade total radius
r Main rotor blade current radius position
T Main rotor thrust
v i Main rotor induced velocity in hover
wRotor induced velocity in far wake
W Climb velocity
δ Inflow angle
ρ Air density
λ Total inflow ratio
λ c Climb inflow ratio
λ i Induced inflow ratio
σ Rotor solidity
ϑ Pitch angle
Ω Rotor angular velocity
Abbreviations
C A D Computer-aided design
CAMComputer-aided manufacturing
CAEComputer-aided engineering
CFDComputational fluid dynamics
GRIPGraphic integrated programming
FSIFluid structure interaction
BERPBritish engineering rotor program
MATLABMatrix laboratory-multi-paradigm programming language and numerical computing environment
SIEMENS NXCAD/CAM/CAE software developed from Unigraphics system

References

  1. Stalewski, W.; Zalewski, W. Performance Improvement of Helicopter Rotors through Blade Redesigning. Aircr. Eng. Aerosp. Technol. 2019, 91, 747–755. [Google Scholar] [CrossRef]
  2. Sagimbayev, S.; Kylyshbek, Y.; Batay, S.; Zhao, Y.; Fok, S.; Lee, T.S. 3D Multidisciplinary Automated Design Optimization Toolbox for Wind Turbine Blades. Processes 2021, 9, 581. [Google Scholar] [CrossRef]
  3. Ma, L.; Zhao, Q.; Zhang, K.; Zhang, X.; Zhao, M. Aeroelastic Analysis and Structural Parametric Design of Composite Rotor Blade. Chin. J. Aeronaut. 2020, 34, 336–349. [Google Scholar] [CrossRef]
  4. Fitzgibbon, T.A.; Woodgate, M.A.; Barakos, G.N. Practical Considerations in Rotor Design Optimisation. In Proceedings of the 47th European Rotorcraft Forum, Virtual, 7–9 September 2021. [Google Scholar]
  5. Lim, J.W. Application of Parametric Airfoil Design for Rotor Performance Improvement. In Proceedings of the 44th European Rotorcraft Forum 2018, ERF 2018, Delft, The Netherlands, 18–21 September 2018; Volume 1. [Google Scholar]
  6. Vu, N.A.; Lee, J.W. Aerodynamic Design Optimization of Helicopter Rotor Blades Including Airfoil Shape for Forward Flight. Aerosp. Sci. Technol. 2015, 42, 106–117. [Google Scholar] [CrossRef]
  7. Allen, L.D.; Lim, J.W.; Haehnel, R.B.; Dettwiller, I.D. Rotor Blade Design Framework for Airfoil Shape Optimization with Performance Considerations. In Proceedings of the AIAA SciTech 2021 Forum, Online, 11–15 January 2021. [Google Scholar]
  8. Tixadou, E. Optimisation of the structure of a helicopter blade in view of reducing cabin vibration. In Proceedings of the 47th European Rotorcraft Forum, Virtual, 7–9 September 2021. [Google Scholar]
  9. Гребеников, А.Г.; Дьяченко, Ю.В.; Коллеров, В.В.; Коцюба, В.Ю.; Малков, И.В.; Урбанович, В.А.; Воронько, И.А. Analysis of Design and Technological Features of the Main Rotor Blades of Heavy Transport Helicopters. Open Inf. Comput. Integr. Technol. 2021, 93, 59–103. [Google Scholar] [CrossRef]
  10. Bailly, J.; Bailly, D.; Didier, B. Multifidelity Aerodynamic Optimization of a Helicopter Rotor Blade. AIAA J. 2019, 57, 3132–3144. [Google Scholar] [CrossRef]
  11. Roy, A.M. Finite Element Framework for Efficient Design of Three Dimensional Multicomponent Composite Helicopter Rotor Blade System. Eng 2021, 2, 69–79. [Google Scholar] [CrossRef]
  12. Ryazanov, A.I. Automated 3D Modeling of Working Turbine Blades. Russ. Eng. Res. 2016, 36, 751–754. [Google Scholar] [CrossRef]
  13. Shabliy, L.S.; Dmitrieva, I.B. Conversion of the Blade Geometrical Data from Points Cloud to the Parametric Format for Optimization Problems. ARPN J. Eng. Appl. Sci. 2014, 9, 1849–1853. [Google Scholar]
  14. Grabowik, C.; Kalinowski, K.; Kempa, W.; Paprocka, I. A Method of Computer Aided Design with Self-Generative Models in NX Siemens Environment. In Proceedings of the IOP Conference Series: Materials Science and Engineering, Mamaia, Romania, 17–20 June 2015; Volume 95. [Google Scholar]
  15. Stalewski, W. Design and Aerodynamic Optimization of Rotorcraft (Polish Translated); Institute of Aviation: Warsaw, Poland, 2017. [Google Scholar]
  16. Stalewski, W.; Zalewski, W. Aerodynamic Design and Optimisation of Main Rotors for Light Rotorcrafts. In Proceedings of the 45th European Rotorcraft Forum 2019, ERF 2019, Warsaw, Poland, 17–20 September 2019; Volume 1. [Google Scholar]
  17. Tamer, A.; Yücekayali, A.; Ortakaya, Y. Design and Optimization of Helicopter Rotor for Minimum Power Required. In Proceedings of the 6th AIAC Ankara International Aeropsace Conference, Ankara, Turkey, 14–16 September 2011. [Google Scholar]
  18. Xie, J.; Xie, Z.; Zhou, M.; Qiu, J. Multidisciplinary Aerodynamic Design of a Rotor Blade for an Optimum Rotor Speed Helicopter. Appl. Sci. 2017, 7, 639. [Google Scholar] [CrossRef] [Green Version]
  19. Slavik, S.; Klesa, J.; Brabec, J. Propeller Selection by Means of Pareto-Optimal Sets Applied to Flight Performance. Aerospace 2020, 7, 21. [Google Scholar] [CrossRef] [Green Version]
  20. Kachel, S.; Rogólski, R.; Kocjan, J. Review of Modern Helicopter Constructions and an Outline of Rotorcraft Design Parameters. Problems of Mechatronics. Armament Aviat. Safety Eng. 2021, 12, 27–52. [Google Scholar]
  21. Pölzlbauer, P.; Kümmel, A.; Desvigne, D.; Breitsamter, C. Numerical Investigation of an Optimized Rotor Head Fairing for the RACER Compound Helicopter in Cruise Flight. Aerospace 2021, 8, 66. [Google Scholar] [CrossRef]
  22. Samad, A.; Tagawa, G.B.S.; Morency, F.; Volat, C. Predicting Rotor Heat Transfer Using the Viscous Blade Element Momentum Theory and Unsteady Vortex Lattice Method. Aerospace 2020, 7, 90. [Google Scholar] [CrossRef]
  23. Lei, Y.; Huang, Y.; Wang, H. Aerodynamic Performance of an Octorotor SUAV with Different Rotor Spacing in Hover. Processes 2020, 8, 1364. [Google Scholar] [CrossRef]
  24. Spyropoulos, N.; Papadakis, G.; Prospathopoulos, J.M.; Riziotis, V.A. Investigating the Level of Fidelity of an Actuator Line Model in Predicting Loads and Deflections of Rotating Blades under Uniform Free-Stream Flow. Appl. Sci. 2021, 11, 12097. [Google Scholar] [CrossRef]
  25. Truong, K.V. Modeling Aerodynamics, Including Dynamic Stall, for Comprehensive Analysis of Helicopter Rotors. Aerospace 2017, 4, 21. [Google Scholar] [CrossRef] [Green Version]
  26. Wstęp do Konstrukcji Śmigłowców: Praca Zbiorowa, 2nd ed.; Wydawnictwa Komunik. i Łączn.: Warszawa, Poland, 2002; ISBN 8320611571.
  27. Brocklehurst, A.; Barakos, G.N. A Review of Helicopter Rotor Blade Tip Shapes. Prog. Aerosp. Sci. 2013, 56, 35–74. [Google Scholar] [CrossRef]
  28. Seddon, J.; Newman, S. Basic Helicopter Aerodynamics; John Wiley & Sons, Ltd.: Hoboken, NJ, USA, 2011. [Google Scholar]
  29. Johnson, W. Helicopter Theory; Dover Publications: New York, NY, USA, 1994; Volume 12. [Google Scholar]
  30. Bramwell, A.R.S.; Done, G.; David, B. Bramwell ’s Helicopter Dynamic; Elsevier: Amsterdam, The Netherlands, 2001. [Google Scholar]
Figure 1. The geometry of the main rotor blade: (a) Rotor shape (c, chord and c0, mounting holes spacing); (b) blade parameters (αS(r/R), geometric twist angle in dependence on nondimensional radius; A and B, airfoil change zone; d, blade thickness) (figure prepared by authors).
Figure 1. The geometry of the main rotor blade: (a) Rotor shape (c, chord and c0, mounting holes spacing); (b) blade parameters (αS(r/R), geometric twist angle in dependence on nondimensional radius; A and B, airfoil change zone; d, blade thickness) (figure prepared by authors).
Materials 15 04275 g001
Figure 2. Blade element aerodynamic forces (figure prepared by the authors).
Figure 2. Blade element aerodynamic forces (figure prepared by the authors).
Materials 15 04275 g002
Figure 3. The algorithm developed in Matlab for determining the distribution of the blade’s parameters and its aerodynamic thrust.
Figure 3. The algorithm developed in Matlab for determining the distribution of the blade’s parameters and its aerodynamic thrust.
Materials 15 04275 g003
Figure 4. Three-blade rotor thrust distribution as a function of twist and mean chord.
Figure 4. Three-blade rotor thrust distribution as a function of twist and mean chord.
Materials 15 04275 g004
Figure 5. Four-blade rotor thrust distribution as a function of twist and mean chord.
Figure 5. Four-blade rotor thrust distribution as a function of twist and mean chord.
Materials 15 04275 g005
Figure 6. Five-blade rotor thrust distribution as a function of twist and mean chord.
Figure 6. Five-blade rotor thrust distribution as a function of twist and mean chord.
Materials 15 04275 g006
Figure 7. Diagram of the algorithm for parametric programming of the main rotor blade geometrics.
Figure 7. Diagram of the algorithm for parametric programming of the main rotor blade geometrics.
Materials 15 04275 g007
Figure 8. Main rotor blade shape curves.
Figure 8. Main rotor blade shape curves.
Materials 15 04275 g008
Figure 9. Model generated with defined curves.
Figure 9. Model generated with defined curves.
Materials 15 04275 g009
Figure 10. Rotor blade mesh with boundary layer: (a) Mesh; (b) boundary layer.
Figure 10. Rotor blade mesh with boundary layer: (a) Mesh; (b) boundary layer.
Materials 15 04275 g010
Figure 11. Single-blade model boundary conditions.
Figure 11. Single-blade model boundary conditions.
Materials 15 04275 g011
Figure 12. One blade panel code model with solver setup.
Figure 12. One blade panel code model with solver setup.
Materials 15 04275 g012
Figure 13. Four blade model boundary conditions.
Figure 13. Four blade model boundary conditions.
Materials 15 04275 g013
Figure 14. Line graphs of aerodynamic thrust dependencies for the one blade model.
Figure 14. Line graphs of aerodynamic thrust dependencies for the one blade model.
Materials 15 04275 g014
Figure 15. Line graphs of aerodynamic thrust dependencies for the complete rotor model.
Figure 15. Line graphs of aerodynamic thrust dependencies for the complete rotor model.
Materials 15 04275 g015
Figure 16. Graphical proposal of design spiral with results of the studies.
Figure 16. Graphical proposal of design spiral with results of the studies.
Materials 15 04275 g016
Table 1. Main rotor technical specifications.
Table 1. Main rotor technical specifications.
Rotor radius7.85 m
Blade mean chord0.44 m
Blade airfoilNACA 23015
Geometric twist−12°
Root cutout0.21R
Root incidence angle5.34°
Table 2. Example of the inflow functions for the inflow modeling in the FLUENT rotor blade analysis.
Table 2. Example of the inflow functions for the inflow modeling in the FLUENT rotor blade analysis.
Inflow X direction function for 5.96 pitch anglecos(PI/180 × (5.96 − (−1.40420898843306 × 10−5 × (z/7.85[m])3 + 0.00154624103083247 × (z/7.85[m])2 − 0.0591574921598430 × (z/7.85[m]) + 1.07777970150355)))
Inflow Y direction function for 5.96 pitch anglesin(PI/180 × (5.96 − (−1.40420898843306 × 10−5 × (z/7.85[m])3 + 0.00154624103083247 × (z/7.85[m])2 − 0.0591574921598430 × (z/7.85[m]) + 1.07777970150355)))
Table 3. Results obtained from analyses of the one blade model.
Table 3. Results obtained from analyses of the one blade model.
Simulation Input ParameterFluentFlightstream
Weight (kg)Analytical Thrust (N)Collective (°)CFD Thrust (N)Δ (%)Panel Method Thrust (N)Δ (%)
640015,6967.4715,5910.67%155201.12%
600014,7157.0614,8701.05%148500.92%
560013,7346.6614,0412.24%141503.03%
520012,7536.2613,2423.83%126021.18%
490012,0175.9612,5924.78%121330.96%
Table 4. Results obtained from analyses of the complete rotor model.
Table 4. Results obtained from analyses of the complete rotor model.
Simulation Input ParameterFluentFlightstream
Weight (kg)Analytical Thrust (N)Collective (°)CFD Thrust (N)Δ (%)Panel Method Thrust (N)Δ (%)
640062,7847.4764,4722.69%7418618.16%
600058,8607.0658,5040.60%7111220.82%
560054,9366.6655,8621.69%6891725.45%
520051,0126.2552,0442.02%6320523.90%
490048,0695.9549,1752.30%5946223.70%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kocjan, J.; Kachel, S.; Rogólski, R. Helicopter Main Rotor Blade Parametric Design for a Preliminary Aerodynamic Analysis Supported by CFD or Panel Method. Materials 2022, 15, 4275. https://doi.org/10.3390/ma15124275

AMA Style

Kocjan J, Kachel S, Rogólski R. Helicopter Main Rotor Blade Parametric Design for a Preliminary Aerodynamic Analysis Supported by CFD or Panel Method. Materials. 2022; 15(12):4275. https://doi.org/10.3390/ma15124275

Chicago/Turabian Style

Kocjan, Jakub, Stanisław Kachel, and Robert Rogólski. 2022. "Helicopter Main Rotor Blade Parametric Design for a Preliminary Aerodynamic Analysis Supported by CFD or Panel Method" Materials 15, no. 12: 4275. https://doi.org/10.3390/ma15124275

APA Style

Kocjan, J., Kachel, S., & Rogólski, R. (2022). Helicopter Main Rotor Blade Parametric Design for a Preliminary Aerodynamic Analysis Supported by CFD or Panel Method. Materials, 15(12), 4275. https://doi.org/10.3390/ma15124275

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop