Skip to main content

Aerodynamics and dynamic stability of micro-air-vehicle with four flapping wings in hovering flight

A Correction to this article was published on 28 June 2020

This article has been updated

Abstract

Recently, a novel concept of flapping Micro-Air-Vehicles (FMAVs) with four wings has been proposed, which potentially utilizes the clap-and-fling effect for lift enhancement and agile maneuvers through an adjustment of wing kinematics. However, the application of the clap-and-fling effect in the four-winged FMAVs is underexplored and the dynamic stability is still unclear. In this paper, aerodynamics and flight dynamic stability of the four-winged FMAVs are studied experimentally and numerically. Results show that the clap-and-fling effect is observed when the flapping frequency is above 18 Hz. Due to the clap-and-fling effect, the lift generation and aerodynamic efficiency are both improved, which is mainly attributed to the fling phase. Further studies show that the clap-and-fling effect becomes weaker as the wing root spacing increases and is almost absent at a wing root spacing of 1.73 chord length. In addition, a wing with an aspect ratio of 3 can increase both lift generation and efficiency due to the clap-and-fling effect. Finally, according to the dynamic stability analysis of the four-winged FMAV, the divergence speed of the lateral oscillation mode is about 4 times faster than that of the longitudinal oscillation mode. Our results can provide guidance on the design and control of four-winged FMAVs.

1 Introduction

With the rapid development of microelectronics and microfabrication technology, the concept of Micro-Air-Vehicle (MAV) [1] has been proposed by the Defense Advanced Research Projects Agency (DARPA) in the 1990s. Compared to conventional large-scale aircrafts, the MAV shows great superiority in the concealment of its movement and the agility within tiny space due to limited wing span and weight. Therefore, once the concept of MAV has been put forward, it rapidly received continuous attention of engineers due to its mass military and civil applications like intelligence reconnaissance, individual combat, and disaster relief.

Motivated by insects’ flight, flapping MAV (FMAV) is thought the most potential MAV configuration in terms of miniaturization [2]. In the last two decades, the unsteady high-lift mechanisms and flight mechanics of natural insects have been widely studied and their feasibility applied into FMAV design has also been proved [3,4,5]. Now, a series of high-lift mechanisms including delayed stall [6], wake capture [7], and clap-and-flip effect [8,9,10,11,12,13] have been clearly addressed. Flying stability and the controlling strategies used by insects have also been systematically revealed. With the theories how insects fly in mind, scientists turned their focus into how to fully apply those flying principles to a man-made FMAV and conducted systematical researches on aerodynamic, structure, power and control problems. For example, extensive experimental and numerical studies have been conducted on the effects of wing geometrical [14,15,16] (including wing aspect ratio, area, span et al.) and on both lift and power consumption. Wing passive deformations beneficial for lift augment and power efficiency enhancement have also been discovered [17,18,19,20]. Then, various MAV prototypes using gear transmission, piezoelectric as well as electromagnetic drives were built, such as Robobee [21] of Harvard University, hummingbird [22] of AeroVironment Inc., KU Beetle [23] of Konkuk University, and Colibri [24] of ULB. After flying test of both vertical take-off and hovering flight, the dynamic stability of those bionic FMAVs has been further studied to guide flight control system design so as to finally achieve maneuvering flight.

Although various FMAV prototypes are capable of free flight, rigorous requirements in size and weight not only impose great difficulty on MAV conceptual design and sub-system design but also limit MAV’s performance. Lift enhancement is always desired to carry more payloads like camera to conduct specific tasks. In addition, traditional FMAV consists of two wings which are not only used for lift and thrust generation but also taken as controlling surfaces to produce torques for maneuvering flight, which brings great difficulty in accurately controlling the flight of FMAV.

In terms to enhance lift production and reduce difficulty in the design of MAV control system, four-winged FMAV combining the concept of FMAV and rotary wing has been proposed, such as Delfly Nimble [25] and RoboBee X-wing [26] (Fig.1). A four-winged FMAV consists of two pairs of flapping wings. Lift augment was designed from the increase in the number of wings, and also benefited from the potential clap-and-fling effect of adjacent wings at the end of both upstroke and downstroke. The clap-and-fling effect was firstly proposed by Weis-Fogh [8] to explain how small wasp generated enough lift. He pointed that wing clap-and-fling effect will induce a circulation around the wing which accounts for the high lift production. This assumption was later confirmed by many scientists [9,10,11,12,13]. Although this principle was generally found used by natural tiny insects whose Reynolds (Re) number is lower than 100, researchers still explored to introduce it into MAV design. Hoang Vu Phan designed a FMAV model with a flapping amplitude of 180° and wings clapped and flinged at the end of upstroke and downstroke. They measured the averaged lift of this model and found that the clap-and-fling effect leads to a lift augment by about 16% [23]. In summary, wing that could benefit from clap-and-fling effect in lift augment has been widely accepted, however, its effect on wing aerodynamic power consumption as well as power efficiency remains unclear and detailed investigation is still needed to design a micro four-wing FMAV for larger lift production.

Fig. 1
figure 1

Typical four-winged FMAV. a Delfly Nimble, b RoboBee X-wing

With the increase in the number of flapping wings, it is possible for the four-winged FMAV to achieve maneuver flight by controlling the overall motion of the flapping wings like that used by multi-rotor drone. Even so, the dynamic stability characteristics of this type of FMAV should still be clarified firstly. In recent years, dynamic flight stability of insects has been revealed and it was found that insect hovering flight is dynamically unstable [27, 28] to enhance maneuverability. Same with insects, bionic FMAV is also of the same dynamic stability property. The four-winged FMAV combines the concepts of FMAV and rotary wing. Unlike FMAVs, multi-rotor drone is capable to stabilize its flight. Thus, it raises the question that whether the flight of the micro four-wing FMAV is inherently stable, which must be answered before the design of FMAV controlling system.

In this paper, aerodynamics and flight dynamic stability of the four-winged FMAV was studied experimentally and numerically. The paper is organized as follows: firstly, both experimental as well as numerical methods (including computational fluid dynamic (CFD) simulation and dynamic stability analysis) are introduced. Then, based on a typical model, to what extent the clap-and-fling effect could affect lift, power consumption as well as power efficiency are studied. Effect of wing geometrical parameters, i.e. wing root distance and aspect ratio, on aerodynamic performance has been further studied using CFD methods. Finally, the longitudinal and lateral dynamic stability of the four-winged FMAV model are analyzed.

2 Method and material

2.1 Lift-measurement platform and experimental model

Averaged lift (L) of the four-winged FMAV model was measured experimentally using a fast lift measurement method to reveal the clap-and-fling effect on averaged lift production, as shown in Fig. 2. The flapping mechanism was fixedly connected to the optical plate through the support column. The regulated power supply provided steady electricity for the flapping mechanism. The laser frequency measurement (measurement accuracy is 0.01HZ) was fixedly connected above the wing and the laser shed onto the top surface of the flapping mechanism at wing root. When the wing flapped and went through the laser, the laser was reflected back by the white resin, and the digital counter was increased. By counting how many times the wing flapped under the measurement during one second, we computed the flapping frequency. A precision electronic scale with measurement accuracy of 0.01 g was used for averaged lift measurement. It was not directly connected with the FMAV model and was just placed below the wing, and the distance between the wing and the electronic scale was about one span. When the wing flapped, the downwash air flowed to the plate. Therefore, the force measured by the precision electronic scale was equal to the lift generated by the flapping wing. It has been tested that this method measures lift with an error of about 10% [22].

Fig. 2
figure 2

Experimental apparatus and model (a) the fast lift measurement platform, (b) experimental test model, (c) Wing structure, (d) the flapping wings with clap-and-fling motion

Two pairs of flapping wings of the four-winged FMAV were completely symmetrical in structure and function, so only one pair of flapping wings was investigated in the experiment. The 3D printed FMAV model (Fig. 2b) was driven by a 7.4 V hollow cup motor, which was decelerated by a gear set with the reduction ratio of 20. Wing flapping amplitude (ϕf) of the model was about 90° and the maximum flapping frequency (n) was about 30 Hz. The wing (Fig. 2c) was made of a leading-edge carbon rod (1 mm diameter), a reinforcing carbon rod (0.5 mm diameter carbon rod) and a vertical carbon rod (1 mm diameter carbon rod) covered by a polyimide film with the thickness of 0.025 mm. The span of wing was 9 cm and the chord length of the wing was 4 cm. During wing flapping, the wing membrane was deformed with the acting of rod inertia force and aerodynamic pressure to achieve angle of attack. Using a high-speed camera, the experiment captured the motion of the wing at the end of the upstroke (Fig. 2d) when the flapping frequency was 25 Hz, at which time a significant clap-and-fling motion can be observed. In addition to measuring the lift of the paired wings, we also measured the lift of a single wing (\( \overline{L} \)S), which was not contributed to clap-and-fling effect. Comparing 2 \( \overline{L} \)S and \( \overline{L} \) at the same flapping frequency, the clap-and-fling effect on averaged lift can be revealed.

2.2 CFD method

The experimental measurement can only obtain the averaged lift and it is difficult to analyze the clap-and-fling effect on lift, energy consumption and aerodynamic efficiency. Therefore, a simplified model was built according to the experimental test model and simulated using CFD method to comprehensively explore the influence of clap-and-fling motion on the aerodynamic performance (lift, energy consumption and lift efficiency) of the four-winged FMAV.

The wing was simplified to a rectangular wing that flaps horizontally. The geometry of the wing was determined by aspect ration AR (defined as the ratio of the length R to the chord length c), the distance between the body center and the wing root CR (CR = 0.7c), and the root spacing (d) between the adjacent wing roots at the end of upstroke. In order to define wing kinematics, the inertial coordinate system og-xgygzg, the body-fixed coordinate system ob-xbybzb, and the wing-fixed coordinate system o-xyz were introduced (Seen in Fig. 3). The og-xgyg plane of the inertial coordinate system was horizontal. The origin of the body-fixed coordinate systems ob was located at the center of the body. The ob-xbyb plane of the body-fixed coordinate system coincided with wing flapping plane, and obxb was symmetric axis of the body in vertical plane and pointed backward. The wing-fixed coordinate system o-xyz moved with wing flaps with oz axis normal to the flapping stroke plane and always pointed upward. ox was tangent to the flapping trajectory of the wing, pointing in the direction of downstroke. The flapping angle (ϕ) and angle of attack (α) were defined using a sinusoidal function as Eqs. 1 and 2, respectively.

$$ \dot{\phi}=\frac{\pi }{T}{\phi}_f\sin \left(\frac{2\pi }{T}t\right) $$
(1)
$$ \dot{\alpha}=\frac{\pi -2\overline{\alpha}}{\varDelta {\tau}_r}\left(\cos \left(\frac{4\pi }{\varDelta {\tau}_r}t\right)-1\right) $$
(2)

where T was the non-dimensional flapping period, \( \overline{\alpha} \) was the angle of attack in the middle of upstroke, and Δτr (Δτr = 0.25 T) was time for wing performing a flipping motion. In the CFD simulation, ϕf = 90° and \( \overline{\alpha} \) =40°. Refer to wing kinematics of those insects using clap-and-fling effect for high lift production, the pitch axis varied during the clap-and-fling motion. When the FMAV clapped its wings, the pitch axis was the leading edge and it turned into the trailing edge after wing entered fling phase. To eliminate the changes in wing vertical displacement caused by the changes in the pitch axis, an additional vertical velocity of wing was introduced as presented in Eq. 3.

$$ {V}_z=\frac{2c\left(1-\sin \alpha \right)}{T-2\varDelta {\tau}_r} $$
(3)
Fig. 3
figure 3

The four-winged FMAV model sketch (a-c) Definition of coordinate systems and kinematics of four-winged FMAV, (d) as well as clap-and-fling motion schematic diagram

Based on wing kinematics defined, Re was defined based on the velocity (U = 2ϕfnr2) at the second moment of wing area (r2) of the flapping wing area as Re = cU/ν, where ν was viscosity coefficient. In this paper, Re was set at 3725 based on the experimental measurement of the four-winged FMAV model.

The commercial software FLUENT was used to calculate the aerodynamic forces and power consumption. Considering that the four flapping wings were symmetrically arranged along the longitudinal and lateral directions of the FMAV respectively, only one flapping wing was calculated and the influence of the adjacent two flapping wings was simulated by setting symmetric boundary conditions. The simulation used dynamic mesh technology to realize the unsteady motion of the wing. The wing surface was constructed by 80,000 structural grids with unstructured grid wrapped. To improve the grid quality, the unstructured grid zone was further divided into two parts: the inner part was dense unstructured grid with a radius of 20c around the wing surface and the outer was a sparse unstructured gird with a radius of 40c. (seen in Fig. 4).

Fig. 4
figure 4

Mesh used for CFD simulation (a) Computational domain mesh, (b) Side view of computational domain

The 3D pressure-based coupled solver with laminar model was used for the transient simulation of the flapping wing. We defined the lift normal to the wing flapping plane and drag tangent to the flapping trajectory of the wing. Lift (L) was always positive in a flapping cycle, while drag (D) was positive in upstroke and negative in downstroke. Aerodynamic power (P = M·ω) can be further obtained based on the aerodynamic moment (M) of the wing around the origin and the angular velocity (ω) of the wing. The lift coefficient, the drag coefficient and the power coefficient CL, CD and CP were further calculated as \( {C}_L=\frac{L}{0.5\rho {U}^2S},{C}_D=\frac{D}{0.5\rho {U}^2S},{C}_P=\frac{P}{0.5\rho {U}^3S} \), where S was wing area. We used \( {\overline{C}}_L \), \( {\overline{C}}_P \) to denote the average lift coefficient and aerodynamic power coefficient during one stroke cycle, and \( {\overline{C}}_D \) was defined as average drag coefficient during half stroke cycle. Based on \( {\overline{C}}_L \), \( {\overline{C}}_P \), aerodynamic efficiency η was computed as \( \eta ={{\overline{C}}_L}^{\frac{3}{2}}/{\overline{C}}_P \) [29].

In order to test the independence of the grid used in the numerical simulation, we have verified the grid density, the boundary domain size and the time step separately. During grid-independence test, the parameters of the flapping wing were selected as Re = 3725, d = 0.28c and AR = 3. Under the conditions of 20 times far field and 800 steps/cycle, mesh 1 (220,000 grids), mesh 2 (360,000 grids), and mesh 3 (650,000 grids) were selected for grid density verification. The results are shown in Table 1. It can be seen from the comparison of the data that the calculation results of the medium density grid (mesh 2) and the dense grid (mesh 3) were relatively close, and the results of the sparse grid (mesh 1) calculation were quite different. Therefore, in order to ensure the calculation accuracy and improve the calculation efficiency at the same time, a 360,000 grid was selected for subsequent calculation. The same method was used to verify the boundary domain size and calculation time step. Finally, the numerical simulation of this study was carried out using a grid with a far field size of 20c and a time step of 0.005 T.

Table 1 Grid-independence test parameters

In addition, in order to verify the difference between the single-wing simulation and the full simulation, the study also compared the results of symmetric processing with that of full model as shown in Table 2. It can be seen from the table that the results obtained by the two methods were basically the same. Therefore, it was feasible to use the symmetric processing to perform the aerodynamic simulation of the multi-wing wing model.

Table 2 Aerodynamic coefficient comparison between symmetric and full model

To validate the CFD method used in this paper, CFD simulation toward a flapping wing with Re = 16,100, AR = 3, ϕ = 120°and α = 45°was conducted and compared with the results calculated by Hao et al. [30]. Instantaneous CL and CD in a flapping cycle from two methods were shown in Fig. 5. As can be seen from Fig. 5, instantaneous CL and CD from these two methods were approximately the same despite some slight differences in the reversal phase. Besides, the stroke-averaged aerodynamic coefficient calculated using the method in this paper was also nearly the same as that calculated by Hao et al. with an error less than 2%. In summary, the CFD simulation method used in this work was reasonable.

Fig. 5
figure 5

Comparison of (a) CL and (b) CD coefficient from reference and ANSYS Fluent

2.3 Dynamic stability model

Based on the idea of insect dynamic stability analysis [27, 31], a six-degree-of-freedom rigid body dynamics equation was established based on the average aerodynamic force and moment (Fig. 3a). The longitudinal and lateral motion equations of the aircraft were as Eq. 4.

$$ {\displaystyle \begin{array}{l}m\left(\dot{u}+ qw\right)=X- mg\sin \theta \\ {}m\dot{v}=Y+ mg\sin \varphi \\ {}m\left(\dot{w}- qu\right)=Z- mg\cos \theta \\ {}{I}_{xx}\dot{p}=L\\ {}{I}_{yy}\dot{q}=M\\ {}{I}_{zz}\dot{r}=N\\ {}\dot{\theta}=q\\ {}\dot{\varphi}=p\end{array}} $$
(4)

where, u, v, w, p, q, r were the corresponding projections of the centroid velocity and angular velocities of the FMAV on the xb, yb, zb axis, respectively. θ was the angle between the obxb axis and the ox axis, and φ was the angle between the obyb axis and the oy axis. Ixx, Iyy and Izz were the moment of inertia around the xb, yb, zb axis, respectively, whose parameters were from the experimental model (Ixx = 5.634 × 10− 5kg·m2, Iyy = 84.296 × 10− 6kg·m2, Izz = 4.810 × 10− 5kg·m2). X, Y, Z, L, M, N were the corresponding projections of aerodynamic forces and moments on the xb, yb, zb axis. The small perturbation equation (shown as Eq.5) can be obtained by linearizing Eq. 4 in the equilibrium state and ignoring the small amount above the second order.

$$ \left[\begin{array}{c}\begin{array}{c}\begin{array}{c}\delta {u}^{+}\\ {}\delta {w}^{+}\end{array}\\ {}\delta {q}^{+}\\ {}\delta \theta \\ {}\delta {v}^{+}\end{array}\\ {}\delta {p}^{+}\\ {}\delta {r}^{+}\\ {}\delta {\phi}^{+}\end{array}\right]=\left[\begin{array}{cc}\begin{array}{cccc}\frac{{X_u}^{+}}{m^{+}}& \frac{{X_w}^{+}}{m^{+}}& \frac{{X_q}^{+}}{m^{+}}& -{g}^{+}\\ {}\frac{{Z_u}^{+}}{m^{+}}& \frac{{Z_w}^{+}}{m^{+}}& \frac{{Z_q}^{+}}{m^{+}}& 0\\ {}\frac{{M_u}^{+}}{{I_{yy}}^{+}}& \frac{{M_w}^{+}}{{I_{yy}}^{+}}& \frac{{M_q}^{+}}{{I_{yy}}^{+}}& 0\\ {}0& 0& 1& 0\end{array}& 0\\ {}0& \begin{array}{cccc}\frac{{Y_v}^{+}}{m^{+}}& \frac{{Y_p}^{+}}{m^{+}}& \frac{{Y_r}^{+}}{m^{+}}& {g}^{+}\\ {}\frac{{L_v}^{+}}{{I_{xx}}^{+}}& \frac{{L_p}^{+}}{{I_{xx}}^{+}}& \frac{{L_r}^{+}}{{I_{xx}}^{+}}& 0\\ {}\frac{N_v}{{I_{zz}}^{+}}& \frac{N_p}{{I_{\mathrm{z}z}}^{+}}& \frac{N_r}{{I_{zz}}^{+}}& 0\\ {}0& 1& 0& 0\end{array}\end{array}\right]\left[\begin{array}{c}\begin{array}{c}\begin{array}{c}\delta {u}^{+}\\ {}\delta {w}^{+}\end{array}\\ {}\delta {q}^{+}\\ {}\delta \theta \\ {}\delta {v}^{+}\end{array}\\ {}\delta {p}^{+}\\ {}\delta {r}^{+}\\ {}\delta {\phi}^{+}\end{array}\right] $$
(5)

where:

$$ {\displaystyle \begin{array}{l}{X}^{+}=\frac{X}{0.5\rho {U}^2{S}_t},{Y}^{+}=\frac{Y}{0.5\rho {U}^2{S}_t},{Z}^{+}=\frac{Z}{0.5\rho {U}^2{S}_t}\\ {}{L}^{+}=\frac{L}{\rho {U}^2{S}_tc},{M}^{+}=\frac{M}{\rho {U}^2{S}_tc},{N}^{+}=\frac{N}{\rho {U}^2{S}_tc}\\ {}\delta {u}^{+}=\frac{\delta u}{U},\delta {w}^{+}=\frac{\delta w}{U},\delta {v}^{+}=\frac{\delta v}{U},\delta {p}^{+}=\delta p\cdotp T,\delta {q}^{+}=\delta q\cdotp T,\delta {r}^{+}=\delta r\cdotp T,{t}^{+}=\frac{t}{T}\\ {}{m}^{+}=\frac{m}{\rho {U}^2{S}_tT},{g}^{+}=\frac{g{t}_w}{U},{I_{xx}}^{+}=\frac{I_{xx}}{0.5{U}^2{S}_t\mathrm{c}{T}^2},{I_{yy}}^{+}=\frac{I_{yy}}{0.5{U}^2{S}_tc{T}^2},{I_{zz}}^{+}=\frac{I_{zz}}{0.5{U}^2{S}_tc{T}^2}\end{array}} $$

Where m was the mass of the two pairs of FMAV and m = 62 g. Because the longitudinal motion equations and lateral motion equations can be decoupled, we dealed with these equations separately.

For the aerodynamic derivatives, we can get the accurate results of this parameters but at a high cost of calculation sources and calculation time, which did not satisfy requirements in the analysis of dynamic stability especially in conceptual design phase. Therefore, the quasi-steady aerodynamic estimation method was adopted to calculate the aerodynamic derivatives. A detailed description and verification of this method can be found in the literature of these researches [27, 28, 31].

3 Results

3.1 Aerodynamic performance of a typical micro four-wing FMAV

Averaged lift of a typical micro four-winged FMAV model was firstly measured to verify the effect of the clap-and-fling effect on lift production. To get rid of the error coming from wing structure and experimental measurement, five pairs of wing are measured and the experiments are repeated three times for each trial. Averaged lift (\( \overline{L} \)) produced by two wings with and without the clap-and-fling effect at different flapping frequencies is presented in Fig. 6. It can be seen from Fig. 6 that whether the clap-and-fling effect could bring lift augment depends on wing flapping frequency. When the flapping frequency is lower than a critical value, which is 18 Hz in this experiment, there is no significant change in lift production. Once the flapping frequency is higher than 18 Hz, lift produced by wing with clap-and-fling effect motion becomes larger than that without the motion and the lift difference caused by clap-and-fling effect almost increases linearly with wing flapping frequency. In the studied range of flapping frequency, the lift is increased by 11.2% due to clap-and-fling effect when f = 25 Hz. As shown in Fig.2d, the clap-and-fling effect of two adjacent wings is mainly caused by wing’s flexible deformation. At a flapping frequency smaller than 18 Hz, the flexible deformation of the flapping wing is not significant and as a result, clap-and-fling effect is also weak. As wing flapping frequency increases, wing passive deformation increases and there is obvious clap-and-fling motion of wing trailing edge at the end of upstroke, which consequently leads to a higher lift production at this phase.

Fig. 6
figure 6

Comparison between wing with or without clap-and-fling effect

To further test the clap-and-fling effect on lift, power consumption and aerodynamic efficiency at high flapping frequency, a CFD simulation is conducted. Note that for the experiment model, wing always pitches around its leading edge. Differently, for the model used in CFD simulation, wing pitching axis during the reverse phase is changed between leading edge and trailing edge, which resembles that of tiny insect based on the observation results, to fully utilize the clap-and-fling effect. Parameters used in this simulation are as follows: Re = 3725, d = 0.28c, AR = 3, ϕ = 90°. \( {\overline{C}}_L \), \( {\overline{C}}_D \), \( {\overline{C}}_P \) and η computed are presented in Table 3. As shown in Table 3, clap-and-fling effect leads to wing \( {\overline{C}}_L \) increasing by 19% from 2.044 to 2.436 at a cost of a larger power consumption, which increases by 23.9%. As the increase in lift is more pronounced, η also increases by about 4.8%. It suggests that wing clap-and-fling effect is not only beneficial for lift augment but for power efficiency enchantment.

Table 3 Average aerodynamic coefficient and efficiency of wing with or without clap-and-fling effect

Figure 7a plots the instantaneous aerodynamic forces and the power consumption in a flapping period. As can be seen from Fig. 7, CL has a periodic change over time with a period at 0.5 T, while CD seems changes without period which is caused by the definition of L and D. In fact, CD in upstroke and downstroke has the same absolute value while the sign of CD is changed. For both wings with or without clap-and-fling effect, there is a larger contribution of aerodynamics observed due to the rapid pitching of the wing. The increase in lift, drag and power consumption caused by the clap-and-flip effect mainly occurs when wings fling (e.g. t = 0~0.125 T and 0.5~0.625 T). In the four-winged FMAV, the clap-and-fling motion repeated twice for a specific wing at the end of both downstroke and upstroke. Consequently, there are two significant increment in instantaneous lift in a flapping period.

Fig. 7
figure 7

Aerodynamics, pressure distribution and flow field of the flapping wing with or without clap-and-fling effect in typical time. a-c lift, drag, aerodynamic power coefficient changes during a stroke cycle, (d) pressure distribution of typical time

In order to further analyze the clap-and-fling effect on the aerodynamic forces production at Re~ο (103), vorticity field around wing with and without clap-and-fling motion at three typical times is compared and analyzed (shown in Fig. 7d). In the fling phase (e.g. t = 0.06 T), a significantly stronger attachment leading edge vortex (LEV) and tip vortex (TV) is formed in the upper surface of the wing due to a rapid pitching. The zone occupied by the LEV and TV is of low pressure, which results in a large normal suction acting on the wing leading to large lift and drag. As seen in Fig. 7d, at this time, LEV and TV formed around the wing with clap-and-fling motion is much intensive with larger strength, so there is much larger lift and drag produced and instantaneous CL, CD and CP increases by 48.3%, 53.6% and 14.5%, respectively. When the wing translates (e.g. t = 0.32 T), the fling motion of wings has already finished and the clap-and-fling effect becomes so weak and the instantaneous structure of the vorticity field around the wing with or without the clap-and-fling effect is basically the same, so the instantaneous aerodynamic force is basically the same at this time. When the wings move to the end of upstroke and the clap motion begins (e.g. t = 0.42 T), the LEV remains attached to the lower surface of the wing. Although a downwash flow is induced, the intensity of the LEV is not influenced significantly. As a result, the wing aerodynamics and power consumption of wing with clap-and-fling effect motion slightly changed and CL, CD as well as CP increase by 12%, 15% and 6%, respectively, compared to its counterpart.

3.2 Influence of geometric parameters on aerodynamics

This section explored the effect of wing geometrical parameters on the aerodynamic performance enhancement led by clap-and-fling effect. Influence of wing root spacing and aspect ratio of the wing was studied using CFD simulation.

3.2.1 Wing root spacing

Three wing roots spacing, that is d = 0.28c, 0.43c, 1.73c, are considered with Re, AR, ϕ unchanged. \( {\overline{C}}_L \), \( {\overline{C}}_D \), \( {\overline{C}}_P \) and η computed at different d are shown in Table 4. CL, CD and CP are plotted in Fig.8. Results of wing without clap-and-cling effect are also presented for comparison. It can be seen from Table 4 that within the studied range of d, wing lift and power efficiency could always be improved by the clap-and-fling effect. When the d is 0.28c, increase in \( {\overline{C}}_L \) can be maximized and its value increases from 2.044 to 2.503 (by 23%). Despite a significant increase in \( {\overline{C}}_P \) (increased by 19%), aerodynamic efficiency in lift production is still slightly increased by 7.4%. As wing root spacing increases, clap-and-fling motion between two adjacent wings is weakened, and lift augment and power efficiency caused by the clap-and-fling effect are all gradually reduced. When the wing root spacing increases larger than 1.73c, the aerodynamic performance of these two types of flapping wing is basically the same, which means the clap-and-fling effect is negligible. As can be seen from Fig. 8, CL, CD and CP produced during two reversal phases reduced with the increment of d (seen from t = 0.06 T, 0.42 T, 0.56 T and 0.92 T), especially in the fling phase (e.g. t = 0.06 T). At t = 0.06 T, wing lift coefficient of d = 0.28c, 0.43c and 0.73c is increased by 55.5%, 32.7% and 3.3% compared to wing without clap-and-fling effect, while it only increased by 13.2%, 8.3% and 0.5% at t = 0.42 T (i.e. clap phase). In summary, a smaller wing root spacing is suggested to fully utilize the “clap-and-fling” effect for a better aerodynamic performance.

Table 4 Averaged aerodynamic coefficient and power coefficient under different wing root spacing
Fig. 8
figure 8

Instantaneous aerodynamic coefficient and power coefficient under different wing root spacing. (a) CL , (b) CD , (c) CP

3.2.2 Aspect ratio

As the aspect ratio of insect wings is in the range of 1~5, therefore three typical aspect ratio, that is AR = 1, 3 and 5, are considered. Please note that the chord length is unchanged with aspect ratio. To remain Re = 3725, the non-dimensional flapping frequency is changed with aspect ratio. Wing with small aspect ratio flaps fast than that with large one. \( {\overline{C}}_L \), \( {\overline{C}}_D \), \( {\overline{C}}_P \) and η of wing with different aspect ratio are shown in Table 5 while time courses of CL, CD and CP in a flapping period are presented in Fig. 9.

Table 5 Average aerodynamic coefficient and power coefficient under different aspect ratio
Fig. 9
figure 9

Instantaneous aerodynamic coefficient and power coefficient under different AR (a) AR = 1, (b) AR = 3, (c) AR = 5

As can be seen from Table 5 that at a small AR (about 1), the clap-and-fling effect plays a reverse role in lift production while power consumption is slightly enlarged. \( {\overline{C}}_L \) of wing with clap-and-fling effect reduced from 6.145 to 5.952 by 3.4% compared to single wing. At a constant Re, wing with small aspect ratio flaps at high frequency and larger aerodynamic contribution is observed in fast reversal phase at the end of both upstroke and downstroke. Compared to single wing, increment in instantaneous lift of wing with clap-and-fling effect is greatly reduced compared to cases with AR > 1. After wings fling and undergo translation, loss in lift further offsets the lift augment from clap-and-fling effect. Taking a complete flapping period into account, wing’s averaged lift and power efficiency decreases.

As AR becomes larger than 1, lift augment and power efficiency can be got by clap-and-fling motion. Furthermore, compared to case at AR = 5, wing with an aspect ratio of 3 can benefit more in lift augment. As seen from Fig. 9, the instantaneous aerodynamic forces of wings with or without clap-and-fling effect at reversal phase both decreased significantly as wing flapping velocity decreased with large AR. From instantaneous aerodynamic coefficient and power coefficient shown in Fig. 8, the peak value of CL is found about 1.5 times larger than that produced by the single wing when AR ≥ 3, causing lift augment in those cases at AR > 1. Further comparison of the instantaneous aerodynamic forces produced by wing with AR = 3 and AR = 5 shows that in case of a large aspect ratio, there is a small amount of lift loss (about 5.8%) occurred during the translation phase compared to single wing while there is a lift increment (about 6%) occurred during the translation when AR = 3, so the overall lift augment of large aspect ratio wing is slightly smaller than that of wing with medium aspect ratio. To summarize, a medium aspect ratio wing is more recommended to design a micro four-wing FMAV when Re is kept constant.

3.3 Dynamic stability of a four-winged FMAV under typical design conditions

The dynamic stability analysis is based on the small perturbations from equilibrium state, so it is necessary to estimate the kinematic parameters of the FMAV in equilibrium state. Here, the quasi-steady aerodynamic method is applied to achieve the equilibrium state of aerodynamic forces and moment. The equilibrium state parameters are roughly as follows: \( \overline{\alpha} \) =45°, ϕf = 90°. In order to solve the small perturbation equation, the stability derivatives need to be calculated firstly in both longitudinal and lateral equations. The calculation method is the same as that used by Sun and Xiong [27] and Sun et al. [28]. Here is an example of the longitudinal U-series aerodynamic derivatives (Xu, Zu and Mu): it is necessary to take a number of effective small disturbances (usually within 10% of the reference velocity), and set w and q to zero, then calculate the aerodynamic forces and moments curve with u. The slope of the curve at the equilibrium state is thus the corresponding aerodynamic derivative. The other aerodynamic derivatives in longitudinal and lateral are computed similarly. Finally, the non-dimensional aerodynamic derivatives of the four-winged FMAV can be obtained as shown in Table 6.

Table 6 Non-dimensional aerodynamic derivatives

Bring the non-dimensional aerodynamic derivative into the system matrix, we can get longitudinal and lateral equations as Eq. 6 and 7, respectively.

$$ \left[\begin{array}{c}\delta {u}^{+}\\ {}\delta {w}^{+}\\ {}\delta {q}^{+}\\ {}\delta \theta \end{array}\right]={A}_0\left[\begin{array}{c}\delta {u}^{+}\\ {}\delta {w}^{+}\\ {}\delta {q}^{+}\\ {}\delta \theta \end{array}\right]=\left[\begin{array}{cccc}-0.049& 0& 0& -0.071\\ {}-0.002& -0.043& 0& 0\\ {}0.014& 0& -0.089& 0\\ {}0& 0& 1& 0\end{array}\right]\left[\begin{array}{c}\delta {u}^{+}\\ {}\delta {w}^{+}\\ {}\delta {q}^{+}\\ {}\delta \theta \end{array}\right] $$
(6)
$$ \left[\begin{array}{c}\delta {v}^{+}\\ {}\delta {p}^{+}\\ {}\delta {r}^{+}\\ {}\delta {\phi}^{+}\end{array}\right]={A}_1\left[\begin{array}{c}\delta {v}^{+}\\ {}\delta {p}^{+}\\ {}\delta {r}^{+}\\ {}\delta {\phi}^{+}\end{array}\right]=\left[\begin{array}{cccc}-0.034& -0.004& 0& 0.050\\ {}-0.041& -0.063& 0& 0\\ {}0& 0& -0.654& 0\\ {}0& 1& 0& 0\end{array}\right]\left[\begin{array}{c}\delta {v}^{+}\\ {}\delta {p}^{+}\\ {}\delta {r}^{+}\\ {}\delta {\phi}^{+}\end{array}\right] $$
(7)

Then the techniques of eigenvalue and eigenvector analysis are applied to solve the longitudinal equations and the lateral equations of the four-winged FMAV, and the corresponding eigenvalues of longitudinal stability matrix A0 and lateral stability matrix A1 are shown in Table 7, respectively. It can be seen that the longitudinal mode and the lateral mode have three modes respectively: mode 1 (unstable oscillation mode), mode 2 (slow subsidence mode) and mode 3 (fast subsidence mode).

Table 7 Eigenvalues of system matrix

The period of longitudinal mode 1(Tosc) and the time to double the initial perturbation (tdouble) (or half the initial perturbations, thalf) can be calculated by Eq. 8.

$$ {\displaystyle \begin{array}{l}{T}_{osc}=2\pi /\hat{\omega}\\ {}{t}_{double}\;\mathrm{or}\;{t}_{half}=0.693/\mid \lambda \mid \left(\lambda\ is\ real\right)\\ {}{t}_{double}\;\mathrm{or}\;{t}_{half}=0.693/\mid \hat{n}\mid \left(\lambda\ is\ complex\right)\end{array}} $$
(8)

where, \( \hat{n} \) and \( \hat{\omega} \) are the real and imaginary parts of the eigenvalue, respectively. Thus, the time constants (dimensionless by the stroke period) of the above modes are shown in Table 8. For the oscillating divergent modes of the FMAV, the longitudinal perturbation needs 93.65 strokes to diverge twice, while the lateral perturbation needs only 21.3 strokes to diverge twice. So, the lateral control response of the four-winged FMAV needs to be faster.

Table 8 Time constants of longitudinal and lateral mode

The eigenvector of each mode determines the magnitude and phase relationship of state variables. The eigenvectors of the longitudinal stability matrix A0 and the lateral stability matrix A1 are shown in Table 9. The longitudinal and lateral eigenvectors are normalized by δθ and δφ, respectively. As can be seen from Table 9, the longitudinal mode 1 is dominated by horizontal forward motion (u) and pitch motion (θ and q), and the vertical motion is about one order smaller in magnitude and can be ignored in this mode. Since the phase difference of u and q is less than 90°, the motion is horizontal forward motion with nose-up and horizontal backward motion accompanied by nose-down. Longitudinal mode 2 is dominated by horizontal forward motion and pitch motion, too. But the phase difference between u and q is 180°, indicating that there is a nose-down angular velocity in nose-up state and an nose-up angular velocity in the nose-down state, which means there will be a tendency to overcome the disturbance and the movement tends to be stable. Longitudinal mode 3 is dominated by ascending and sinking motions. The lateral mode 1 is dominated by the side slip motion (v) and the roll motion (p), and the yaw motion is of a smaller magnitude. Since the phase difference of v and p is greater than 90 degrees, the motion is a right horizontal movement with rolling left or left horizontal movement with rolling right. The lateral mode 2 is dominated by yaw motion. The lateral flight mode 3 is dominated by the side motion (v) and the roll motion (p), too. But the phases of δv+, δp + and δϕ+ indicate that the motion is a translation motion to left with a right roll angle accompanied by rolling left. Rolling angular velocity is opposite to the roll angle.

Table 9 Eigenvectors of the system matrix

4 Discussion

In order to further explore the dynamic stability influence of the four-winged FMAV, we compared both the longitudinal mode and lateral mode of four-winged FMAV with those of some insects and a pair of flapping wings MAV (as shown in Table 10 and Table 11, respectively) [28, 29, 31]. It can be found that they have the same longitudinal modes. The longitudinal modes of the four-winged FMAV, insects and a pair of flapping wings are divided into three types: the unstable oscillatory divergence mode, the stable fast subsidence mode and the stable slow subsidence mode. For the control design of aircraft, we are more concerned with the unstable oscillatory divergence mode. The time to double the initial perturbation is determined by the real part of eigenvalues. It can be inferred that the tdouble of four-winged FMAV is about one order larger in magnitude than that of the KU Beetle, thus, leaving enough margin for the response time of the control system. From the longitudinal eigenvalues, we can see the mode of four-winged FMAV is more like the mode of insects. Thus, the four-winged FMAV can have more possibility to maneuver as the natural insects.

Table 10 Longitudinal modal comparison of FMAVs and insects
Table 11 Lateral modal comparison of four-winged FMAV and insects

While the four-winged FMAV’s lateral mode is not the same as some insects (i.e. BB, DF, HM). Insects like BB, DF, HM have one stable oscillatory mode, one stable subsidence mode and one unstable divergence mode. But the four-winged FMAV has one unstable oscillatory divergence mode and two stable subsidence mode, the same as its longitudinal mode. The modal structure of the four-winged FMAV is more appropriate to the modal structure of HB. Although the oscillatory subsidence mode of four-winged FMAV is unstable, its time to double the initial perturbation is smaller than BB, DF and HM. Thus, we can call the unstable oscillatory divergence mode as a nearly neutrally stable oscillatory mode [29]. It means that we can control this unstable mode in time through the control system.

Considering the longitudinal and lateral mode comprehensively, the four-winged FMAV is more advantageous for the control system, because it provided enough response time for the control system.

5 Conclusion

A four-winged FMAV, which utilizes the clap-and-fling effect to improve its aerodynamic performance, is proposed in this research. Experimental measurements of wing kinematics are conducted and numerical simulations are employed to unfold the effect of clap-and-fling effect on the four-winged FMAV. A dynamic stability analysis is then carried out using the modal analysis method. Our findings are concluded as follows:

  1. 1.

    For a flapping wing with an aspect ratio of 3, the clap-and-fling effect is observed when the flapping frequency is above 18 Hz. Within the parameter space of our study, the clap-and-fling effect results a linear lift enhancement over the flapping frequency. Due to the clap-and-fling effect, the lift enhancement is around 11.2% at a flapping frequency of 25 Hz.

  2. 2.

    The effect of clap-and-fling effect on the lift enhancement of the four-winged FMAV is intensified as the wing root spacing decreases. Moreover, the lift enhancement is mainly attributed to the fling phase. The clap-and-fling effect is almost absent when the wing root spacing is above 1.73 chord length.

  3. 3.

    A further study on the aspect ratio (AR) shows that the lift enhancement due to the clap-and-fling effect is observed at AR = 3 or 5. In addition, the improvements of cycle-averaged lift generation and efficiency for a flapping wing with AR = 3 are 19.2% and 4.8%, respectively. However, as a result of the clap-and-fling effect, the cycle-averaged lift generation and efficiency of a flapping wing with AR = 1 is reduced by 3.1% and 5.6%, respectively.

  4. 4.

    The longitudinal and lateral dynamic stability modes of the four-winged FMAV are akin to insects and can be separated into three components: the unstable oscillatory divergence mode, the stable fast subsidence mode and the stable slow subsidence mode. The divergence speed of the lateral oscillatory divergence mode is about 4 times faster than that of the longitudinal unstable mode. To be specific, the longitudinal oscillatory divergence mode of the four-winged FMAV is doubled over 93.65 strokes, while the lateral oscillatory divergence modes doubles their magnitudes over 21.3 strokes.

Our study can thus contribute to the design of four-winged FMAVs to better exploit the clap-and-fling effect, thereby extending its flight endurance under a heavy payload. Our dynamic stability analysis can shed light on the control of four-winged FMAVs.

Availability of data and materials

All data obtained or analyzed in the course of this study are included in this published paper.

Change history

  • 28 June 2020

    An amendment to this paper has been published and can be accessed via the original article.

Abbreviations

CFD:

Computational Fluid Dynamic

DARPA:

Defense Advanced Research Projects Agency

FMAVs:

Flapping Micro-air-vehicles

LEV:

Leading Edge Vortex

MAV:

Micro-air-vehicles

TV:

Tip Vortex

References

  1. McMichael, J. M. (1997). Micro air vehicles-toward a new dimension in flight. http://www.arpa.2.gov/tto/MAV/mav_auvsi.html

  2. Karasek M (2014) Robotic hummingbird: design of a control mechanism for a hovering flapping wing micro air vehicle. Universite libre de Bruxelles

    Google Scholar 

  3. Lentink D, Dickinson MH (2009) Rotational accelerations stabilize leading edge vortices on revolving fly wings. J Exp Biol 212(16):2705–2719

    Article  Google Scholar 

  4. Van Truong T, Le TQ, Park HC, Byun D (2017) Experimental and numerical studies of beetle-inspired flapping wing in hovering flight. Bioinspir Biomim 12(3):036012

    Article  Google Scholar 

  5. Zheng Y, Wu Y, Tang H (2016) An experimental study on the forewing–hindwing interactions in hovering and forward flights. Int J Heat Fluid Flow 59:62–73

    Article  Google Scholar 

  6. Ellington CP, Van Den Berg C, Willmott AP, Thomas AL (1996) Leading-edge vortices in insect flight. Nature 384(6610):626

    Article  Google Scholar 

  7. Harbig RR, Sheridan J, Thompson MC (2014) The role of advance ratio and aspect ratio in determining leading-edge vortex stability for flapping flight. J Fluid Mech 751:71–105

    Article  MathSciNet  Google Scholar 

  8. Weis-Fogh T (1975) Unusual mechanisms for the generation of lift in flying animals. Sci Am 233(5):80–87

    Article  Google Scholar 

  9. Cheng X, Sun M (2019) Revisiting the clap-and-fling mechanism in small wasp Encarsia Formosa using quantitative measurements of the wing motion. Phys Fluids 31(10):101903

    Article  Google Scholar 

  10. Weis-Fogh T (1973) Quick estimates of flight fitness in hovering animals, including novel mechanisms for lift production. J Exp Biol 59(1):169–230

    Google Scholar 

  11. Lighthill MJ (1973) On the Weis-Fogh mechanism of lift generation. J Fluid Mech 60(1):1–17

    Article  MathSciNet  Google Scholar 

  12. Maxworthy T (1979) Experiments on the Weis-Fogh mechanism of lift generation by insects in hovering flight. Part 1. Dynamics of the ‘fling’. J Fluid Mech 93(1):47–63

    Article  Google Scholar 

  13. Wakeling JM, Ellington CP (1997) Dragonfly flight. II. Velocities, accelerations and kinematics of flapping flight. J Exp Biol 200(3):557–582

    Google Scholar 

  14. Carr ZR, Chen C, Ringuette MJ (2013) Finite-span rotating wings: three-dimensional vortex formation and variations with aspect ratio. Exp Fluids 54(2):1444

    Article  Google Scholar 

  15. Phillips N, Knowles K, Bomphrey RJ (2015) The effect of aspect ratio on the leading-edge vortex over an insect-like flapping wing. Bioinspir Biomim 10(5):056020

    Article  Google Scholar 

  16. Nan Y, Karásek M, Lalami ME, Preumont A (2017) Experimental optimization of wing shape for a hummingbird-like flapping wing micro air vehicle. Bioinspir Biomim 12(2):026010

    Article  Google Scholar 

  17. Young J, Walker SM, Bomphrey RJ, Taylor GK, Thomas AL (2009) Details of insect wing design and deformation enhance aerodynamic function and flight efficiency. Science 325(5947):1549–1552

    Article  Google Scholar 

  18. Wu P, Stanford BK, Sällström E, Ukeiley L, Ifju PG (2011) Structural dynamics and aerodynamics measurements of biologically inspired flexible flapping wings. Bioinspir Biomim 6(1):016009

    Article  Google Scholar 

  19. Du G, Sun M (2010) Effects of wing deformation on aerodynamic forces in hovering hoverflies. J Exp Biol 213(13):2273–2283

    Article  Google Scholar 

  20. Zhou C, Zhang Y, Wu J (2020) Effect of flexibility on unsteady aerodynamics forces of a purely plunging airfoil. Chin J Aeronaut 33(3). https://doi.org/10.1016/j.cja.2019.08.002

  21. Pérez-Arancibia NO, Ma KY, Galloway KC, Greenberg JD, Wood RJ (2011) First controlled vertical flight of a biologically inspired microrobot. Bioinspir Biomim 6(3):036009

    Article  Google Scholar 

  22. Keennon M, Klingebiel K, Won H (2012) Development of the nano hummingbird: a tailless flapping wing micro air vehicle. In: 50th AIAA aerospace sciences meeting including the new horizons forum and aerospace exposition, p 588

    Google Scholar 

  23. Phan HV, Kang T, Park HC (2017) Design and stable flight of a 21 g insect-like tailless flapping wing micro air vehicle with angular rates feedback control. Bioinspir Biomim 12(3):036006

    Article  Google Scholar 

  24. Roshanbin A, Altartouri H, Karasek M, Preumont A (2017) COLIBRI: a hovering flapping twin-wing robot. Int J Micro Air Vehicles 9(4):270–282

    Article  Google Scholar 

  25. Karásek M, Muijres FT, De Wagter C, Remes BD, de Croon GC (2018) A tailless aerial robotic flapper reveals that flies use torque coupling in rapid banked turns. Science 361(6407):1089–1094

    Article  MathSciNet  Google Scholar 

  26. Fuller SB (2019) Four wings: an insect-sized aerial robot with steering ability and payload capacity for autonomy. IEEE Robotics Automation Lett 4(2):570–577

    Article  Google Scholar 

  27. Sun M, Xiong Y (2005) Dynamic flight stability of a hovering bumblebee. J Exp Biol 208(3):447–459

    Article  Google Scholar 

  28. Sun M, Wang J, Xiong Y (2007) Dynamic flight stability of hovering insects. Acta Mech Sinica 23(3):231–246

    Article  MathSciNet  Google Scholar 

  29. Jianghao W, Chao Z, Yanlai Z (2016) Aerodynamic power efficiency comparison of various micro-air-vehicle layouts in hovering flight. AIAA J 55(4):1265–1278

  30. Hao J., Wu J., and Zhang Y. (2019). Aerodynamic performance of a passive pitching model on bionic flapping wing micro air vehicles. Applied bionics and biomechanics. 2019:1504310

  31. Mou X, Sun M (2012) Dynamic flight stability of a model hoverfly in inclined-stroke-plane hovering. J Bionic Eng 9(3):294–303

    Article  Google Scholar 

Download references

Acknowledgments

The authors would like to thank Dr. Chen Long and Mr. Gu Mancang for revising this manuscript.

Funding

This research was supported by National Natural Science Foundation of China (NSFC, No. 11672022 and No. 11902017).

Author information

Authors and Affiliations

Authors

Contributions

Our group worked together on the topic and the research output is coming from our joint effort. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Chao Zhou.

Ethics declarations

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Cheng, C., Wu, J., Zhang, Y. et al. Aerodynamics and dynamic stability of micro-air-vehicle with four flapping wings in hovering flight. Adv. Aerodyn. 2, 5 (2020). https://doi.org/10.1186/s42774-020-0029-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s42774-020-0029-0

Keywords