 Open Access
 Total Downloads : 83
 Authors : R. Jagadeesh , S. Madhu Mathi , Dr. K. Vijayaraja
 Paper ID : IJERTV8IS040062
 Volume & Issue : Volume 08, Issue 04 (April – 2019)
 Published (First Online): 02042019
 ISSN (Online) : 22780181
 Publisher Name : IJERT
 License: This work is licensed under a Creative Commons Attribution 4.0 International License
Design and Analysis of Flapping Wing Unmanned Aerial Vehicle

Jagadeesh
Department of Aeronautical Engineering KCG College of Technology
Chennai, India

Madhu Mathi
Department of Aeronautical Engineering KCG College of Technology
Chennai, India
Dr. K. Vijayaraja
Department of Aeronautical Engineering KCG College of Technology
Chennai, India
Abstract This project presents the design of a Flapping wing UAV which is inspired by various bird mechanism and its action during flight. In this project, the real actions have been tried to convert into a perfect mechanism in order to get a stable flight manuvering. Biomemic plays a major role in this design. The design has been made CATIA V5 with all the parameters calculated according to the bird selected. The moving mesh analysis is completed with the help of ANSYS 19.0. This UAV can be widely used for surveillance for civilian and military applications.
KeywordsCATIA; ANSYS; flapping; UAV; lift; dynamic; meshing;
Figure 1 Flight Action of Flapping Wing UAV

INTRODUCTION
The flapping wing aircraft which will produce lift and thrust by the flapping mechanism. Using Biomimic various Ornithopter, designs have been suggested for civil and military applications especially for the purpose of surveillance. In this paper, the highly aerodynamic design for flapping wing UAV with advanced specifications have been created. The flapping wing UAV made which can be used for surveillance or reconnaissance of a particular target and also for a specific environment without its own consciousness. The Ornithopter uses battery power, gear mechanism, which enables to increase the number of flaps. We are bringing down the specifications of various birds and trying it to convert it into perfect real time mechanism. As of the first initiative, the basic and operating principles of flight have been studied in order to understand the flapping mechanism.

For an airplane/ bird to stay at a constant height, Lift force upwards = Weight force downwards

For an airplane/bird to stay at a constant speed, Thrust force forwards = opposing force of Drag
Ease of Use
Figure 2 Up stroke And Down stroke action
Figure 3 Flight Operating Principle
The parameter has been calculated for the design. After the analysis of 2D design, the 3D design has been created with the respect.


PARAMETERS CALCULATED FOR DESIGNING THE FLAPPING WING UAV

Design Parameters
Table I
S. NO
DESIGN PARAMETERS
DENOTION
UNITS
VALUE
1
LENGTH
L
M
0.23
2
WING SPAN
B
M
0.58
3
MASS
M
Kg
0.476
4
WING AREA
S
m2
0.0745
5
ASPECT RATIO
AR
———
4.515
6
TORQUE
T
N
0.358
7
WING LOADING
W/S
Kg/ m2
10.836

Flight Parameters

FREQUENCY
Figure 4. Flapping frequency
f 1.08 (m1/ 3 g1/ 2 b1 S1/ 4 1/ 3)
f (Small Birds) 116.3 m1/6
f (Large Birds) 28.7 m1/3
The frequencies for the parameters designed are
f = 2.6 Hz
f = 2.6 Hz

FLIGHT SPEED (U)
The relationship between flight speed & the mass of the bird can be given,
U 4.77 m1/6
Where m is the mass of the bird and U is the Flight speed.
U = 4.212 m/s
U = 4.212 m/s

FLAPPING ANGLE ()
Flapping angle varies as sinusoidal function. and its rate is given by following equations.
(t) = 22.538o
(t) = 22.538o
(t) = max cos2 f t

FLAPPING RATE (r)
The equation for flapping rate is given as
r(t) = 2 f t sin 2 f t
r(t) = 374.943 Hz.s
r(t) = 374.943 Hz.s

PITCHING RATE ()
The angle of climb (slope) is the angle between the horizontal axis and the aircraft flight path.
Figure 6. Pitching angle
Figure 7. Flapping action
The pitch is the angle between the horizontal axis and the longitudinal axis of the aircraft.
(t) = (r(i)/B) o cos (2 f t + ) (t) = 12.5o


Velocity Parameters

VERTICAL COMPONENET OF RELATIVE WIND VELOCITY (Vy)
Vysin( ) (r(i). .cos( ) ) (0.75.c. .cos( )))
Vy = 2.427m/s

HORIZONTAL COMPONENET OF RELATIVE WIND VELOCITY (Vx)
The horizontal component of relative wind velocity is given by
VxU cos( ) (0.75.c. .sin( )) Vx m/s
Figure 5 Flapping Angle

RELATIVE VELOCITY (Vrel)
The relative velocity for the given parameters is given by
Vrel (Vx2Vy2)1/2
Vrel 5.34 m/s

RELATIVE ANGLE ()
Relative angle between the two velocity components and the effective angle of attack
= tan1 (Vx / Vy) = 62.97o


ANGLE OF ATTACK

EFFECTIVE ANGLE OF ATTACK (eff)
Effective Angle of Attack is that part of a given angle of attack that lies between the chord of an airfoil and the effective airflow.
Figure 8. Effective angle of attack
Effective airflow is a line representing the resultant velocity of the disturbed airflow.
eff
eff

RELATIVE ANGLE OF ATTACK ()
Figure 9. Relative angle of attack
= 6O

FLOW RELATIVE ANGLE OF ATTACK (I)
Figure 10. Angle of attack
I = 2.4O

CYCLE ANGLE () The cycle angle is given by
= 58o

SECTION MEAN PITCH ANGLE (I) The section mean pitch angle is given by
i = 60o


FORCE PARAMETERS

HORIZONTAL COMPONENT OF FORCE( Fhor) The horizontal component of force is given by
FhordLcsin.cosdNnc sin().cosdDdcos.cos Fhor = 1.232 N

VERTICAL COMPONENT OF FORCE ( Fhor) FverdLccos.cosdNnc cos().cos.cos
Fver = 0.2777N

NORMAL FORCE( Nnc)
The vertical component of force is given by
Nnc = 0.306N


LIFT PARAMETERS

THEODORSON LIFT DEFICIENCY (C (K))
Theodorsen Lift Deficiency factor which is a function of reduced frequency k
C

LIFT COEFFICIENT DUE TO CIRCULATION (Clc)
The section lift coefficient due to circulation (KuttaJoukowski condition) for flat plate is given by
Clc = 0.17

SECTIONAL LIFT (LC)
The section lift can thus be calculated by
LC = 0.157N

INSTANTANEOUS LIFT (L) The instantaneous lift is given by
L = 0.565 N

LIFT ALONG SPAN (LX)
The lift along the span is given by
LX = 1.105N


THRUST PARAMETERS

INSTANTANEOUS THRUST (T) The instantaneous thrust is given by
T = 0.204 N

THRUST ALONG SPAN (TX) The thrust along span is given by
Tx = 1.16
Pin = 0.9 W
b) POWER OUTPUT (Pout)
Pout = 0.72 W
b) EFFICENCY ()


DRAG PARAMETERS

INDUCED DRAG COEFFICIENT (Cdi)
Cdi = 0.002468

INDUCED DRAG ( Di)
Di = 0.001708 N

PARASITE DRAG COEFFICIENT (Cdp)
Cdp = 11.1257

PARASITE DRAG (Dp)
Dp = 2.887 N

TOTAL DRAG (Dd)
Dd = Dp + Di Dd = 2.888 N

SKIN FRICTIONAL DRAG COEFFICIENT (Cf)
Cf = 2.529

PRESSURE DIFFERENCE BETWEEN UPPER AND LOWER SURFACE
Pressure (symbol: p or P) is the force applied perpendicular to the surface of an object per unit area over which that force is distributed.
= 82.22

DESIGN

SELECTION OF AEROFOIL
The aerofoil to be selected must be similar to the shape of the bird wing.
Figure 12. NACA 4309
So on analyzing different aerofoil and finding similarities the NACA 4309 aerofoil is been selected

DESIGN OF AEROFOIL
The design of aerofoil in CATIA V5 by importing the coordinates of aerofoil from excel sheet.
Figure 13. Aerofoil

CONNECTING RODS
The connecting rod which
Figure 11. Pressure difference
Gauge pressure is the pressure relative to the ambient pressure.
connects the two gears and two sides of the assembled wing structure which creates the flapping motion.
J. POWER
p = 0.507 BAR
Figure 14. Connecting rod
a) POWER INPUT (Pin)

WING SKELETON STRUCTURE
The aerofoil and different crank that are used to connect the aerofoil to construct the skeleton structure of the wing which is driven by the gear mechanism.
Figure 15. Wing skeleton structure

RUNNING GEAR
The running gear is the main gear which is attached to the wing skeleton structure which is connected to the connecting rod to make the flapping motion.
Figure 19. Top body
H. ASSEMBLED WING STRUCTURE ALONG WITH GEAR The wing structure after the assembling of aerofoil and connecting rod is then connected with running gear and pinion.
Figure 20. Assembled gear with wing

PINION
Figure 16. Running gear

ASSEMBLING OF OTHER PARTS
After assembling the gear and pinion to the wing skeleton structure the body structure is merged with the wing structure.
The pinion gear is the gear which drives the main gear. The
gear is attached to the key which is driven by rotating the key.
Figure 21. Assembled structure


BODY DESIGN
Figure 17. Pinion
The design of body parts is designed according to the shape of the birds analyzed.
Figure 18 Side Body
Figure 22. Assembled structure side view
Figure 23. Assembled structure cross view
Figure 24. Assembled structure front view
J. SHEET METAL APPLICATION

Mesh Deformation Y Velocity
Figure 27 Mesh Y Velocity

Relative X Velocity
Figure 28 Relative X Velocity

Relative Y Velocity
Figure 24 Applied Skin Over The Wing


RESULTS
A. Dynamic velocity
F. Vorticity
Figure 29 Relative Y Velocity
Figure 25 Velocity
B. Mesh Deformation X Velocity

Tangential Velocity
Figure 30 Vorticity
Figure 26 Mesh X Velocity Figure 31 Tangential Velocity

Radial velocity M. Dynamic Pressure

X Velocity
Figure 32 Radial Velocity
N. Total Pressure
Figure 37 Dynamic Pressure
J. Y Velocity
Figure 33 X Velocity

Density
Figure 38 Total Pressure
Figure 34 Y Velocity

Temperature
Figure 39 Density

Velocity Angle
Figure 35 Velocity Angle

Cell Reynolds number

Enthalpy
Figure 40 Temperature
Figure 36 Reynolds Number
Figure 41 Enthalpy

Wall Temperature
Figure 42 Wall Temperature


CONCLUSION

In this paper the flapping wing UAV has been designed and analysed to evaluate its aerodynamic performance in the airflow. Using GAMBIT, the 2D structure has been designed and analysed by FLUENT. The results show that this structure have good aerodynamic design with standard flow parameters. With the reference of 2D design, the 3D structure has been with all its interior parts and gear systems. The results obtained from the FLUENT analysis clearly proves and ensures that it can work better with the flow over it.

Entropy
REFERENCES

Total Energy

Internal Energy

Total Enthalpy
Figure 43 Entropy
Figure 44 Total Energy
Figure 45 Internal Energy

Zachary John Jackowski, Design and Construction of Autonomous Ornithopter, June 7, 2009

Akshay Bhargava, Apoorv Bhargava, Saurabh Kr. Jangid , Ornithopter Development and Operation, March 24,2015

M Afzaal Malik, Farooq Ahmad, Effect of Different Design Parameters On Lift, Thrust and Drag of an Ornithopter, July 2,2010

Jithin M P, Jude Olivero, Naveen M Nair, Sarath R ,Vinod Yeldho Baby, Design of Ornithopter Flapping Wing Mechanism , April 04,2018

Alphan Canga, Michael Delia, Alexander Hyman, Angela Nagelin,
Fluidic Muscle Ornithopter, March 27, 2015

BoHsun Chen, LiShu Chen, Yueh Lu, ZihJie Wang and PeiChun Lin,
Design of a Butterfly Ornithopter, July 7,2016

K. P. Preethi Manohari Sai, K. Bharadwaj, K. Ravi Teja, Kalyan Dagamoori,
Kasiraju Venkata Sai Tarun, Vishnu Vijayan , Design, Fabrication and Testing Of Flapping Wing Micro Air Vehicle, January 16, 2016

S. Adjerid and J. E. Flaherty, A moving FEM with error estimation and refinenment for
1D time dependent PDEs, SIAM J. Numer. Anal. 23 (1986), 778{795

A movingmesh FEM with local refinement for parabolic PDEs, Comput. Meth.
Appl. Mech. Engrg. 55 (1986), 3{26

D. A. Anderson, Adaptive mesh schemes based on grid speeds, AIAA Paper 831931 (1983), 311

Equidistribution schemes, poisson generators, and adaptive grids, Appl.
Math. Comput. 24 (1987), 211{227

D.C. Arney and J.E. Flaherty, A twodimensional mesh moving technique for time dependent partial differential equations, J. Comput. Phys. 67 (1986), 124{144

B.N. Azarenok, Realization of a secondorder Godunovs scheme, Comput. Meth. in Appl.
Mech. and Engin. 189 (2000), 1031{1052

Variational barrier method of adaptive grid generation in hyperbolic problems of
gas dynamics, SIAM J. Numer. Anal. 40 (2002), 651{682

B.N. Azarenok, S.A. Ivanenko, and T. Tang, Adaptive mesh redistribution method based on Godunovs scheme, Commun. Math. Sci. 1 (2003), 152{179

B.N. Azarenok and T. Tang, Reactive flow calculations on moving meshes, 2005, To appear in J. Comput. Phys

M.J. Baines, Moving Finite Elements, Oxford University Press, 1994

Grid adaptation via node movement, Appl. Numer. Math. 26 (1998), 77{9}
Figure 46 total enthalpy