Научная статья на тему 'EFFECT OF PARAGLIDING WING DOME SHAPE ON ITS AERODYNAMIC CHARACTERISTICS'

EFFECT OF PARAGLIDING WING DOME SHAPE ON ITS AERODYNAMIC CHARACTERISTICS Текст научной статьи по специальности «Строительство и архитектура»

CC BY
158
75
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
AERODYNAMICS / PARACHUTE / DOUBLE-MEMBRANE GLIDING PARACHUTE / DISCRETE VORTEX METHOD / DOME-SHAPED PARAGLIDING PARACHUTE / AERODYNAMICS COEFFICIENTS

Аннотация научной статьи по строительству и архитектуре, автор научной работы — Ovchinnikov V. V., Petrov Yu. V., Ganiev Sh. F.

Double-membrane gliding parachutes (DGP) obtain their wide variety of application, including the solution of cargo transportation problems. This parachute is a flexible canopy, which shape is maintained by ram air. In terms of the aerodynamic performance calculation and analysis when operating, DGP is the most complex aero elastic system. The computation of DPG aerodynamic performance is only possible, utilizing the methods of nonlinear aerodynamics and the nonlinear theory of elasticity methods.This paper investigates the aerodynamic characteristics of stable geometric shapes for various gliding parachutes, taking into account their dome shape both chord-wise and span-wise. Notably, the volumetric parachute profile is modeled by its median surface. The research, conducted by the authors, showed that such an aero elastic model of DGP allows you to obtain results that reflect correctly the qualitative effects of detached and free streamline flow. To solve the problem about the airflow over a gliding parachute, considering its canopy curvature, the method of discrete vortices with closed frames is employed, which allows you to calculate the paragliding wing aerodynamic performance within a wide range of angles of attack. There is also a possibility of flow separation simulation. The ideal incompressible liquid flow over the median surface of a stable shape for a double-membrane gliding parachute is regarded. The parachute fabric porosity is not analyzed, since the upper and lower DGP panels are made of either the low permeable or non-porous fabric. In the separated flow past, the aerodynamic coefficients are identified by time averaging to its large values after computing. The DGP aerodynamic performance computation results are given at a different value of its dome shape, as in the free streamline flow as in the flow separation. The computed coefficients, that allow us to consider the influence of canopy dome shape on its aerodynamic characteristics, are obtained. The proposed technique can be used for operational estimates of aerodynamic forces while designing and planning a pipe experiment.

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Текст научной работы на тему «EFFECT OF PARAGLIDING WING DOME SHAPE ON ITS AERODYNAMIC CHARACTERISTICS»

Civil Aviation High Technologies

Vol. 24, No. 06, 2021

УДК 629.735.07

DOI: 10.26467/2079-0619-2021-24-6-54-65

EFFECT OF PARAGLIDING WING DOME SHAPE ON ITS AERODYNAMIC CHARACTERISTICS

V.V. OVCHINNIKOV1, YU.V. PETROV2, SH.F. GANIEV2

1Bauman Moscow State Technical University, Moscow, Russia 2Moscow State Technical University of Civil Aviation, Moscow, Russia

Double-membrane gliding parachutes (DGP) obtain their wide variety of application, including the solution of cargo transportation problems. This parachute is a flexible canopy, which shape is maintained by ram air. In terms of the aerodynamic performance calculation and analysis when operating, DGP is the most complex aero elastic system. The computation of DPG aerodynamic performance is only possible, utilizing the methods of nonlinear aerodynamics and the nonlinear theory of elasticity methods. This paper investigates the aerodynamic characteristics of stable geometric shapes for various gliding parachutes, taking into account their dome shape both chord-wise and span-wise. Notably, the volumetric parachute profile is modeled by its median surface. The research, conducted by the authors, showed that such an aero elastic model of DGP allows you to obtain results that reflect correctly the qualitative effects of detached and free streamline flow. To solve the problem about the airflow over a gliding parachute, considering its canopy curvature, the method of discrete vortices with closed frames is employed, which allows you to calculate the paragliding wing aerodynamic performance within a wide range of angles of attack. There is also a possibility of flow separation simulation. The ideal incompressible liquid flow over the median surface of a stable shape for a double-membrane gliding parachute is regarded. The parachute fabric porosity is not analyzed, since the upper and lower DGP panels are made of either the low permeable or non-porous fabric. In the separated flow past, the aerodynamic coefficients are identified by time averaging to its large values after computing. The DGP aerodynamic performance computation results are given at a different value of its dome shape, as in the free streamline flow as in the flow separation. The computed coefficients, that allow us to consider the influence of canopy dome shape on its aerodynamic characteristics, are obtained. The proposed technique can be used for operational estimates of aerodynamic forces while designing and planning a pipe experiment.

Key words: aerodynamics, parachute, double-membrane gliding parachute, discrete vortex method, dome-shaped paragliding parachute, aerodynamics coefficients.

INTRODUCTION

Double-membrane gliding parachutes (DGP) obtain their wide practical application, including the cargo transportation issues solution. This parachute is a flex-wing, which shape is maintained by ram air. From the point of view of the aerodynamic characteristics calculation and analysis during operation, DGP is the most complex aero elastic system. The calculation of aerodynamic performance for a similar system is merely possible, using the methods of nonlinear aerodynamics and the nonlinear theory of elasticity.

A large number of experimental and theoretical researches, associated with the determination of the DGP dynamic and aerodynamic characteristics1, for example, [1-9] are known. At the same time, obtaining reliable analytical results of the aerodynamic characteristics calculation is significantly complicated by the fact that DGP is actually a flex-wing, which shape is supported by the dynamic air pressure, so it represents an aero elastic system of immense complexity.

In order to identify characteristics of the given system, it is necessary to employ the methods of nonlinear aerodynamics and the nonlinear theory of elasticity, which specifies a significant computational problem. In this sphere the major researches on the basis of vortex methods were carried out by the Central Aerohydrodynamic Institute (TsAGI), the Scientific Research Institute of Parachute Construction, the Moscow State University, Zhukovsky Air Force Engineering Academy (see, for example, [10-13]).

1 The Study of a Parachute Wing in a Wind Tunnel T-101 TsAGI. (1976). Scientific and Technical Department TsAGI

№ 3415. Moscow: TsAGI, 76 p.

Vol. 24, No. 06, 2021

Civil Aviation High Technologies

In this article, the aerodynamic characteristics of various paragliding parachute stable shapes, taking into account their curvature, are explored. In this respect, the volumetric parachute profile is modeled by means of its median surface [14, 15]. The calculations, carried out by authors, showed that this DGP model allows you to obtain results reflecting correctly the qualitative effects of separation flow and free streamline flow. The DGP aerodynamic performance computational results with a different degree of a canopy curvature, as in the free streamline flow as in the separation flow, are given. The design coefficients, that enable you to consider the influence of a canopy dome shape on its aerodynamic characteristics, are obtained. The proposed technique can be used for operational assessments of aerodynamic forces during the design phase and while planning a pipe experiment. The obtained results can be useful while designing the DGP, setting up and conducting the pipe experiments.

METHODS AND METHODOLOGY OF STUDY

In this paper, in order to solve the problem of moving air current over a parachute, the method of discrete vortices with closed frames [12] is used, which permits us to calculate the aerodynamic performance of parachutes with a rather high practical accuracy. A principal assumption, demonstrated in [15], is used, according to which the DGP is modeled by a thin carrying surface. The given assumption allows you to conduct the research efficiently, both when solving the free streamline flow problems and considering the flow separation. A typical vortex diagram of the cutting out shape DGP for the closed-frame discrete vortex method is shown in Figure 1.

Fig. 1. Vortex diagram for the gliding parachute median surface

Replacing the DGP profile with its median surface makes it possible to obtain results, correctly reflecting the qualitative effects of separation flow and free streamline flow and can be applied when studying the major consistent patterns of airflow over parachute finite span wings. We will explore the flow by ideal incompressible liquid over the median surface of the double-membrane paragliding parachute of a stable shape. The parachute fabric permeability is not analyzed, since the DGP upper and lower panels are made of either low permeable or non-porous fabric, so that at the first stage of studies their permeability can be neglected (if necessary, it is taken into account by well-known modifications of the task boundary conditions [13]). In the separation flow, the aerodynamic coefficients are identified by time averaging to its large values after computing.

RESEARCH RESULTS AND THEIR DISCUSSION

The cutting out shape of the double-membrane gliding parachute changes significantly in the air flow (fig. 2), its wing takes a stable dome shape of span ZH with the camber height f which makes 10-20% of the initial span l. Obviously, the stated deformation influences the integral wing aerodynamic performance. This paper proposes a technique that allows you to assess the influence of the canopy dome shape on its aerodynamic characteristics and give formulae for a curved wing recalculation.

Civil Aviation High Technologies

Vol. 24, No. 06, 2021

Fig. 2. The dome shape of the wing of the gliding parachute

As it has been shown in [15], the aerodynamic performance of the parachute cutting out shape is calculated with acceptable accuracy, applying the method of discrete vortices with closed frames as in the free streamline flow (small angles of attack) as in the separation flow.

Here we are going to illustrate how it will be possible to obtain the wing aerodynamic performance with the known camber height by a computational scheme for a cutting out shape.

First, let us explore the case where the parachute wing has a small extension A. Generally speaking, for the most of practicable DGPs, the elongation does not exceed 2. (So, for a gliding parachute system of PO-9 series A = 1.24, for PO-9-2 and PO-9-7 A = 1.7). The sail-wings, which elongation value can reach 3^4, are the exception. The following symbols are used for analytical studies in the given article. The dimensionless aerodynamic performance, assigned to the area in the design of the curved wing 5H, will be designated as cyand cy as the symbol for the aerodynamic characteristics, assigned to the area in the design of the original wing S (cutting out shape). Then we have

Y

Y

Cy qSИ ' qS

(1)

where q is the air pressure, Y is the DGP normal force. When describing the parameters the subscript

Y

V indicates that one or another characteristic refers to a curved wing, for example, cVH = — c is the

* qS

coefficient of normal force for a curved wing, attributed to the non-curved wing area.

It is known that for small extension wings, the slope of lift curve is expressed as [16]

ya

жЯ

T

к Я

(2)

In the area where the aerodynamic performance is linear, it is possible to write (hereinafter, to reduce the entry, the subscript "a" in the designations for the lift coefficient is omitted):

_ ^и ~ _ ^и

суи — ~суи — ~суа•

(3)

Substituting c" into (3) according to (2) and taking into account that AH = —, we get

V

5и ¿и %

Слги — k Œ — k — Œ.

Уи S Би S

(4)

а

Том 24, № 06, 2021_Научный Вестник МГТУ ГА

Vol. 24, No. 06, 2021 Civil Aviation High Technologies

For non-deformed wing

-у ~ ^y" ~ ,v s

Hence follows the following ratio

i2

cy = Cya = k — a. (5)

тИт)2- <6>

This formula is true for the wings of small elongation at small angles of attack. In a similar way, for wings of very large elongation which curve slope of lift can be considered weakly dependent on A value, we obtain:

cyu = ^H (7)

Cy I

Let us analyze the curvature effect on parachute wing drag. We will first make an analysis for small extension wings in the linear region. Taking into account the dependence for the drag coefficient (polar equation):

Cx = cXo + Ac£ (8)

where cXo is the wing profile drag, ACy is the drag dependent on lift. So, we have the following for a curved wing:

2

r = (p + A c2 A = r + X ?2 AA £ = r + a £ r2 (9)

V ^o^ H^yu/ 5 lX0h ^-^Hcyu\sJ 5 h s '-yu-

Since at small angles of attack you can take cy = c£a, so, A = —. Then, considering (2), we get:

Cy

cXH = cX0H + "¡¿X~CyH = cX0H + cyu, (10)

^X = CX0 + . (11)

Since the value of drag at zero lift depends only on the profile shape, so, cxqh = cxq . Then

. S 2 r I S r2 S r2 I S r2 2

+тсуи Сх°+ш2СУи и2 cy+ki2 ъ у ф2-

= 1 + ■

S 2 I S 2 Сх0 , 1

0 kl2 Uy 0 kl2 Uy s c2

kl

2 ^y

To convert the denominator, we use (11), and consider (6) in the numerator, after that we finally have:

tt^+^^+TMI)2). (12)

cx~cxO

Научный Вестник МГТУ ГА_Том 24, № 06, 2021

Civil Aviation High Technologies Vol. 24, No. 06, 2021

Formula (12) shows how the drag coefficient of the curved wing changes compared to the non-curved wing. The second addend appears when the wing has the drag equal to cx0 (for ideal liquid, it is observed when the wing curves along the chord). It is clear that the wing camber along the chord slightly attenuates the effect of curvature. If cx0 _ 0 (there is no chord curvature), then the dome-shaped wing drag coefficient changes in the same way as the coefficient of lift.

Conducting the similar argumentation for the wing of big elongation, we get:

схи _ — + Cx0 (1 _1и.) (13)

Су l Су I

The given formulae, obtained for small angles of attack, were later checked by numerical calculations under various flow conditions. A straight wing was taken as the initial carrying surface, coinciding in the shape and size with the lower surface of the parachute PO-9 (A = 1.24). The aerodynamic characteristics cxa, cya of this wing were calculated under its deformation over the cylindrical surface (fig. 2).

Figure 3 shows the dependences of cya coefficient of this wing from the angle of attack a in the free streamline flow for various deformations. Figure 3 shows clearly the slope of the lift curve of the

deformed wing decreases when its span varies from 100 to 70% from the initial value and with the in- f

crease of camber height f = - from 0 to 20% of span l respectively. The wing drag coefficient without curvature along the chord depends on the canopy dome shape as well as the coefficient of lift.

Fig. 3. Dependence of wing coefficient cya on the angle of attack a in the free streamline flow

for different values of span ln

Figure 4 shows the same dependencies, as Figure 2 does, after dividing them by (y) - the recalculation coefficient from formula (6). Now it is obvious that all the dependencies practically coincide, which indicates the possibility of recalculating the aerodynamic characteristics of the dome-shape canopy of small elongation according to formula (6) in the free streamline flow within the range of a < 10° angles of attack.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Vol. 24, No. 06, 2021

Civil Aviation High Technologies

Fig. 4. Dependence of wing coefficient cya on an angle of attack a in the free streamline flow for different values of span lH after recalculation, using formula (6)

In case the wing surface is additionally curved along the chord (Figure 5, where PO-9 parachute median surface, which is the surface of two curvatures, is shown), in the free streamline flow, it is possible to recalculate the aerodynamic characteristics in conformity with formulae (6), (13), as it is illustrated in Figures 6 and 7.

/„ -1 0,95.

0,90 ,

0**

5 10 15 20 25 a> °

Fig. 6. Lift coefficient of PO-9 parachute median surface considering its dome shape (wing with double curvature)

Civil Aviation High Technologies

Vol. 24, No. 06, 2021

-ya

1,2

0,8

0,4

lK = 0,95 ^ \

0,90 ___

Суаи ~ С^и/ 1 0 cya

10

15

20

25

a,

Fig. 7. To the possibility of recalculation of the lift coefficient for a dome-shaped wing with double curvature

In some cases, it is necessary to take into account the possibility of the flow separation from the surface of the parachute wings. Figure 8 shows the dependencies of cya coefficients in case of the separation flow of deformed wings when the span decreases from 100% to 85% from the original value due to their deformation. During the calculations, it was found out that the recalculation of coefficients, using formula (6), does not lead to positive results in the separation flow.

¿и = 1 0,95 ^^^^^. —^ ■

0,90

■ 0,85

10 15 20 25 30 35 a,0

Fig. 8. Lift coefficient for the dome-shaped wing in the separation flow

Figure 9 represents the dependencies of aerodynamic performance of the deformed wings after dividing by -U coefficient. You can see that the curves are almost congruent.

Vol. 24, No. 06, 2021

Civil Aviation High Technologies

cya 0,9

0,6

0,3

0

10 15 20 25 30 35 a,°

Fig. 9. To the possibility of recalculation of the lift coefficient of dome-shaped wings in the separation flow

It turns out that in the separation flow, the conversion calculation of the aerodynamic coefficients can be conducted, using formula (7) obtained for wings of large elongation. It can be explained by an insignificant influence of edges effects compared to the interaction effect of the vortex sheets in the separation flow over the parachute wing. Figure 10 and Figure 11 illustrate the formula application for the conversion calculation of the drag force coefficient (13), derived for a large elongation wing to recalculating the similar coefficient for a small elongation wing in the separation flow. It is clear that the use of this formula is quite justified.

Thus, the numerical simulation justifies the use of formulae (7) and (13) to recalculate the aerodynamic characteristics obtained for the parachute cutting shape on the dome-shaped wing in the flow separation.

0,8 0,6

0,4 0,2 0

10 15 20 25 30 35 a, °

Fig. 10. Drag coefficient of a dome-shaped wing in the separation flow

Суаи = QJl) Cya

0,95

----

0,90 0,85 = =

Civil Aviation High Technologies

Vol. 24, No. 06, 2021

С 0,8

0,6 0,4

0,2 О

Fig. 11. To the possibility of recalculating the drag force coefficient of the dome-shaped wings in the separation flow

CONCLUSION

This paper provides the investigation results of influence of the median surface curvature (dome shape) of a double-membrane gliding parachute on its basic aerodynamic performance (normal force coefficient cya and drag coefficient cxa) at various values of the parachute dome, its elongation and angles of attack. The results were obtained, using the DGP mathematical model that was developed by the authors, as in the free streamline flow as in the flow separations from the parachute surface.

It is shown that in the free streamline flow at the high angles of attack and the parachute significant curvature (ZH up to 0,7), the normal force coefficient decrease comes to up to 50%. When taking into account the flow separation, the flow pattern changes significantly and cya coefficient decrease does not exceed 10%.

As a result of the calculated data and the analytical studies processing, allowing us to estimate the aerodynamic performance value, taking into account parachute curvature, were obtained. Such an approach enables you to utilize the simplified conversion formulae while designing the DGP and avoid solving the non-linear problems of the elasticity theory, taking into account quite accurately (the error does not exceed 3-5%) the curvature of wing span-wise.

Thus, the authors proposed the solution algorithm for aerodynamic problems when conducting the investigation of parachute optimal shapes and sizes: at the first stage, aerodynamic performance for the flat median surface of the parachute, as in the free streamline flow as in the separation flow, are calculated. Then the parachute dome shape is taken into account by means of utilizing the proposed conversion calculation coefficients. It is obvious that the DGP design the final problems should be solved in a general setting, using complete mathematical models.

REFERENCES

1. Nicolaides, J.D. (1971). Parafoil wind tunnel tests. Technical Rept. AFFDL TR-70-146. University of Notre Dame, 206 p. DTIC_AD0731564

2. Boffadossi, M. and Savorgnan, F. (2016). Analysis on aerodynamic characteristics of a paraglider airfoil. Aerotecnica Missili & Spazio, vol. 95, issue 4, p. 211-218. DOI: 10.1007/BF03404729

Vol. 24, No. 06, 2021

Civil Aviation High Technologies

3. Kulhanek, R. (2019). Identification of a degradation of aerodynamic characteristics of a paraglider due to its flexibility from flight test. Aircraft Engineering and Aerospace Technology, vol. 91, no. 6, p. 873-879. DOI: 10.1108/AEAT-06-2018-0162

4. Wannan, W., Qinglin, S., Shuzhen, L. and others. (2018). Accurate calculation of aerodynamic coefficients of parafoil airdrop system based on computational fluid dynamic. International Journal of Advanced Robotic Systems, vol. 15, issue 2, p. 1-16. DOI: 10.1177/1729881418766190

5. Uddin, M.N. and Mashud, M. (2010). Wind tunnel test of a paraglider (flexible) wing canopy. International Journal of Mechatronics Engineering IJMME-IJENS, vol. 10, no. 3, p. 7-13.

6. Mashud, M. and Umemura, A. (2006). Experimental investigation on aerodynamic characteristics of a paraglider wing. The Japan Society for Aeronautical and Space Sciences, vol. 49, issue 163, p. 9-17. DOI: 10.2322/tjsass.49.9

7. Belloc, H. (2015). Wind tunnel investigation of a rigid paraglider reference wing. Journal of Aircraft, vol. 52, no. 2, p. 703-708. DOI: 10.2514/1.C032513

8. Liang, W., Sun, Q., Tao, J. and others. (2017). Dynamic modeling of parafoil based on CFD simulation and least square. Beijing Ligong Daxue Xuebao/Transaction Beijing Institute Technology, vol. 37, no. 2, p. 157-162. DOI: 10.15918/j.tbit1001-0645.2017.02.009

9. Ortega, E., Flores, R. and Pons-Prats, J. (2016). Ram-air parachute simulation with panel methods and staggered coupling. Journal of Aircraft, vol. 54, no. 2, p. 1-18. DOI: 10.2514/1.C033677

10. Aparinov, A.A., Aparinov, V.A. and Setukha, A.V. (2013). Primeneniye vikhrevykh metodov k zadacham dinamiki parashyutov [Vortex method application to problems of parachute dynamics]. Trudy XVI Mezhdunarodnogo simpoziuma MZDOZMF-2013 [Proceedings of the XVI International Symposium MZDOZMF-2013]. Kharkov: KhNU, p. 41-44. (in Russian)

11. Rysev, O.V., Ponomarev, A.T., Vasil'yev, M.I., Vishnyak, A.A., Dneprov, I.V. and Moseev, Yu.V. (1996). Parashyutnyye sistemy [Parachute systems]. Moscow: Nauka, 288 p. (in Russian)

12. Lyalin, V. and Morozov, V. (2005). Comparative efficiency analyses of the parachute system for space vehicle return blocks. 18th AIAA Aerodynamic Decelerator Systems Technology Conference and Seminar. AIAA-2005-1682. Germany, Munich, 23-26 May. DOI: 10.2514/6.2005-1682

13. Lyalin, V., Morozov, V. and Ponomarev, A. (2005). Mathematical modelling of parachute loading taking into account the canopy shape changes. 18th AIAA Aerodynamic Decelerator Systems Technology Conference and Seminar. AIAA-2005-1676. Germany: Munich, 23-26 May. DOI: 10.2514/6.2005-1676

14. Pisarev, I.V. and Setukha, A.V. (2014). Transferring the boundary conditions to the middle surface for the numerical solution of a boundary value problem in the linear wing theory. Numerical Methods and Programming, vol. 15, no. 1, p. 109-120. (in Russian)

15. Ovchinnikov, V.V. and Petrov, Yu.V. (2018). The aerodynamic chracteristics calculation methodology of two-shell paragliders. Civil Aviation High Technologies, vol. 21, no. 3, p. 91-100. DOI: 10.26467/2079-0619-2018-21-3-91-100 (in Russian)

16. Ovchinnikov, V.V. and Petrov, Yu.V. (2017). Chislennyye metody issledovaniya aerouprugosti letatelnykh apparatov: Monografiya [Numerical methods for the study of aircraft aeroe-lasticity: a Monograph]. Moscow: Izdatelskiy dom Academii imeni N.Ye. Zhukovskogo, 160 p. (in Russian)

INFORMATION ABOUT THE AUTHORS

Valery V. Ovchinnikov, Doctor of Technical Sciences, Professor, Head of the Fundamentals of Physics Chair (Specialized Educational and Scientific Center-2), Bauman Moscow State Technical University, vvovchinnikov@bk.ru.

Civil Aviation High Technologies

Vol. 24, No. 06, 2021

Yury V. Petrov, Doctor of Technical Sciences, Professor, Head of the Engineering Mechanics and Engineering Graphics Chair, Moscow State Technical University of Civil Aviation, yu.petrov@m stuca.aero.

Shamil F. Ganiev, Candidate of Technical Sciences, Associate Professor of the Chair of Flight and Life Safety, Moscow State Technical University of Civil Aviation, shamgan@mstuca.aero.

ВЛИЯНИЕ АРОЧНОСТИ ПАРАШЮТНОГО КРЫЛА НА ЕГО АЭРОДИНАМИЧЕСКИЕ ХАРАКТЕРИСТИКИ

12 2 В.В. Овчинников , Ю.В. Петров , Ш.Ф. Ганиев

1 Московский государственный технический университет имени Н.Э. Баумана,

г. Москва, Россия

2Московский государственный технический университет гражданской авиации,

г. Москва, Россия

Двухоболочковые планирующие парашюты (ДПП) находят широкое практическое применение, в том числе для решения задач грузоперевозки. Данный парашют является мягким крылом, форма которого поддерживается набегающим потоком воздуха. С точки зрения расчета и анализа аэродинамических характеристик в процессе эксплуатации ДПП является сложнейшей аэроупругой системой. Вычисление аэродинамических характеристик такой системы возможно только с привлечением методов нелинейной аэродинамики и нелинейной теории упругости. В данной работе исследуются аэродинамические характеристики установившихся геометрических форм различных планирующих парашютов с учетом их арочности как по хорде, так и по размаху. При этом объемный профиль парашюта моделируется его срединной поверхностью. Проведенные авторами исследования показали, что такая аэроупругая модель ДПП позволяет получать результаты, правильно отражающие качественные эффекты отрывного и безотрывного обтекания. Для решения задачи об обтекании планирующего парашюта потоком воздуха с учетом его арочности используется метод дискретных вихрей с замкнутыми рамками, который позволяет вычислять аэродинамические характеристики парашютов в широком диапазоне углов атаки, а также имеется возможность моделирования отрыва потока. Рассматривается обтекание срединной поверхности установившейся формы двухоболочкового планирующего парашюта потоком идеальной несжимаемой жидкости. Проницаемость ткани парашюта не учитывается, так как верхнее и нижнее полотнища ДПП выполняются либо из слабо проницаемой, либо из непроницаемой ткани. При отрывном обтекании аэродинамические коэффициенты определяются путем усреднения по времени после расчета до его больших значений. Приводятся результаты расчетов аэродинамических характеристик ДПП при различном значении его арочности как при безотрывном обтекании, так и при наличии отрывов. Получены расчетные коэффициенты, позволяющие учесть влияние арочности парашюта на его аэродинамические характеристики. Предлагаемая методика может использоваться для оперативных оценок аэродинамических сил на этапе проектирования и при планировании трубного эксперимента.

Ключевые слова: аэродинамика, парашют, двухоболочковый планирующий парашют, метод дискретных вихрей, арочность парашюта, аэродинамические коэффициенты.

СПИСОК ЛИТЕРАТУРЫ

1. Nicolaides J.D. Parafoil wind tunnel tests. Technical Rept. AFFDL TR-70-146. University of Notre Dame, June 1971. 206 p. DTIC_AD0731564

2. Boffadossi M., Savorgnan F. Analysis on aerodynamic characteristics of a paraglider airfoil // Aerotecnica Missili & Spazio. 2016. Vol. 95, iss. 4. P. 211-218. DOI: 10.1007/BF03404729

3. Kulhanek R. Identification of a degradation of aerodynamic characteristics of a paraglider due to its flexibility from flight test // Aircraft Engineering and Aerospace Technology. 2019. Vol. 91, no. 6. P. 873-879. DOI: 10.1108/AEAT-06-2018-0162

4. Wannan W., Qinglin S., Shuzhen L. and others. Accurate calculation of aerodynamic coefficients of parafoil airdrop system based on computational fluid dynamic // International Journal of Advanced Robotic Systems. 2018. Vol. 15, iss. 2. P. 1-16. DOI: 10.1177/1729881418766190

Vol. 24, No. 06, 2021

Civil Aviation High Technologies

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

5. Uddin M.N., Mashud M. Wind tunnel test of a paraglider (flexible) wing canopy // International Journal of Mechatronics Engineering IJMME-IJENS. 2010. Vol. 10, no. 3. P. 7-13.

6. Mashud M., Umemura A. Experimental investigation on aerodynamic characteristics of a paraglider wing // The Japan Society for Aeronautical and Space Sciences. 2006. Vol. 49, iss. 163. P. 9-17. DOI: 10.2322/tjsass.49.9

7. Belloc H. Wind tunnel investigation of a rigid paraglider reference wing // Journal of Aircraft. 2015. Vol. 52, no. 2. P. 703-708. DOI: 10.2514/1.C032513

8. Liang W., Sun Q., Tao J. and others. Dynamic modeling of parafoil based on CFD simulation and least square // Beijing Ligong Daxue Xuebao/Transaction Beijing Institute Technology. 2017. Vol. 37, no. 2. P. 157-162. DOI: 10.15918/j.tbit1001-0645.2017.02.009

9. Ortega E., Flores R., Pons-Prats J. Ram-air parachute simulation with panel methods and staggered coupling // Journal of Aircraft. 2016. Vol. 54, no. 2. P. 1-18. DOI: 10.2514/1.C033677

10. Апаринов А.А., Апаринов В.А., Сетуха А.В. Применение вихревых методов к задачам динамики парашютов // Труды XVI Международного симпозиума МЗД0ЗМФ-2013. Харьков: ХНУ, 2013. С. 41-44.

11. Рысев О.В. Парашютные системы / О.В. Рысев, А.Т. Пономарев, М.И. Васильев, А.А. Вишняк, И В. Днепров, Ю.В. Мосеев. М.: Наука, 1996. 288 с.

12. Lyalin V., Morozov V. Comparative efficiency analyses of the parachute system for space vehicle return blocks // 18th AIAA Aerodynamic Decelerator Systems Technology Conference and Seminar. AIAA-2005-1682. Germany: Munich, 23-26 May 2005. DOI: 10.2514/6.2005-1682

13. Lyalin V., Morozov V., Ponomarev A. Mathematical modelling of parachute loading taking into account the canopy shape changes // 18th AIAA Aerodynamic Decelerator Systems Technology Conference and Seminar. AIAA-2005-1676. Germany: Munich, 23-26 May 2005. DOI: 10.2514/6.2005-1676

14. Писарев И.В., Сетуха А.В. Снесение граничного условия на срединную поверхность при численном решении краевой задачи линейной теории крыла // Вычислительные методы и программирование. 2014. Т. 15, № 1. С. 109-120.

15. Овчинников В.В., Петров Ю.В. Методика расчета аэродинамического нагружения двухоболочковых планирующих парашютов // Научный Вестник МГТУ ГА. 2018. Т. 21, № 3. С. 91-100. DOI: 10.26467/2079-0619-2018-21-3-91-100

16. Овчинников В.В., Петров Ю.В. Численные методы исследования аэроупругости летательных аппаратов: монография. М.: ИД Академии имени Н.Е. Жуковского, 2017. 160 с.

СВЕДЕНИЯ ОБ АВТОРАХ

Овчинников Валерий Валерьевич, доктор технических наук, профессор, заведующий кафедрой «основы физики» (СУНЦ-2) МГТУ им. Н.Э. Баумана, vvovchinnikov@bk.ru.

Петров Юрий Владимирович, доктор технических наук, профессор, заведующий кафедрой технической механики и инженерной графики МГТУ ГА, yu.petrov@mstuca.aero.

Ганиев Шамиль Фангалиевич, кандидат технических наук, доцент кафедры безопасности полетов и жизнедеятельности МГТУ ГА, shamgan@mstuca.aero.

Поступила в редакцию 01.07.2021 Received 01.07.2021

Принята в печать 25.11.2021 Accepted for publication 25.11.2021

i Надоели баннеры? Вы всегда можете отключить рекламу.