This paper presents a design procedure to achieve a flapping wing mechanism for a micro-air vehicle that coordinates both the wing swing and wing pitch with one actuator. The mechanism combines a planar four-bar linkage with a spatial four-bar linkage attached to the input and output links forming a six-bar linkage. The planar four-bar linkage was designed to control the wing swing trajectory profile and the spatial four-bar linkage was designed to coordinate the pitch of the wing to the swing movement. Tolerance zones were specified around the accuracy points, which were then sampled to generate a number of design candidates. The result was 29 designs that achieve the desired coordination of wing swing and pitch, and a prototype was constructed.
Introduction
Flapping wing micro-air vehicles are winged aircraft generally smaller than 15 cm that mimic insect flight. The wings of these aircraft generally swing in a planar path without control of the wing pitch angle. This pitch angle is set by air drag forces that shape a membrane that forms the wing. See, for example, the aerovironment nanohummingbird [1].
Recent research by Taha et al. [2,3] shows that coordinating the wing pitch angle with the wing swing movement improves flight aerodynamics. This paper presents the synthesis procedure for a one degree-of-freedom flapping wing mechanism that coordinates the wing swing and pitch angles to achieve a specified performance.
This new flapping wing mechanism uses a spatial four-bar linkage to control the wing pitch that is built on a planar four-bar linage that controls the swing movement. The result is a six-bar mechanism, Fig. 1, that uses a single input to coordinate both the wing swing and pitch angles.

The six-bar mechanism is composed of the spatial RSSR and planar four-bar mechanisms. The frames of reference used to define the rotation axes are shown.
In what follows, we present the design methodology, solve for an example design, and present a prototype for this flapping wing mechanism.
Literature Review
The aerovironment nanohummingbird [1], and the Harvard RoboBee [4] have created flying flapping wing micro-air vehicle's whose wings pitch due to the aerodynamic drag induced by the wing's planar swinging motion. Our design provides a linkage, which coordinates both the wing's swing and pitch follow specific functions. Yan et al. [5], Conn et al. [6], and Balta et al. [7] use parallel four-bar linkages and a spatial linkage to coordinate the swing and pitch of the wing, while we achieve this using a single six-bar mechanism. Recently, the Robo Raven [8] and the Bat Bot [9] have utilized a different flying mode in which extreme deformation of large wings is used to achieve flight.
A four-bar linkage that coordinates input and output angles is known as a function generator, Svoboda [10] and Freudenstein [11]. In our case, this linkage transforms a constant rotation of an input crank to an oscillating swing movement of the wing. The synthesis theory for function generators of this type is given by Brodell and Soni [12].
The control of the wing pitch is achieved by an RSSR function generator, see Denavit and Hartenberg [13] and Suh and Radcliffe [14]—R denotes a revolute or hinged joint and S denotes a spherical or ball joint. Our design approach uses the specified angle data to position the input and output cranks of the spatial four-bar linkage, and then solves the constraint equations for the SS link that connects the two, as described by Innocenti [15] and McCarthy and Soh [16].
The connection of the input and output links of a planar four-bar linkage with a spatial RSSR linkage creates a spatial six-bar linkage. The synthesis of spatial six-bar linkages has focused on the control of end-effector position and orientation rather than function generation. See, for example, Sandor et al. [17], who designed an RSSR-SS spatial six-bar and Chiang et al. [18], who designed an RSCC-RRS spatial six-bar, where C denotes a cylindric joint. Recent research includes the design of a spatial six-bar linkage to guide the end effector to trace a spatial path, see Chung [19].
This paper is the first to explore the use of a spatial six-bar linkage to coordinate two outputs to a given input.
Wing Swing and Wing Pitch Requirements
A study of flapping wing designs for micro-air vehicles [7] shows that planar crank-rocker linkages with a passive position of the wing pitch provide effective wing performance. However, Yan et al. [20] show that coordinated control of the wing pitch and wing swing movement improves the aerodynamics of a micro-air vehicle.
where the driving angle is , and ω is the flapping frequency, Fig. 2.
Synthesis of the Wing Swing Mechanism
In order to design the wing swing mechanism, we follow Brodell and Soni [12] to obtain a four-bar linkage that has a time ratio of one. This ensures the speed of the flapping movement is the same in both the forward and back strokes.
where and r4 are the lengths of the ground link, the input crank, the coupler link, and the follower link, respectively.
Our goal is to achieve the four-bar linkage output that matches the wing swing function at θ = 0 and where and the minimum of γ is at . Thus, and . The result is = −0.045 or −2.6 deg and k0 = −0.045 or −2.6 deg, as seen in Fig. 3.
Synthesis of the Pitch Mechanism
In order to control the pitch of the flapping wing, we introduce an RSSR linkage connecting the input and output links of the swing mechanism that will orient the pitch of the wing during the flapping movement.
The planar wing swing mechanism NDEG, shown in Fig. 1, is positioned in the ground frame W such that the fixed pivot G of the link EGC is on the z-axis of W. The axes of the four joints of the planar four-bar NDEG are also directed along the z-axis of W. The rotation γ of EGC provides the swing of the wing. The output crank CB of the RSSR linkage NABC controls the pitch of the wing around the axis of C. The output of the wing swing mechanism NDEG is an input to the wing pitch mechanism NABC. The dimension g is selected by the designer.
The location of the fixed pivot of input crank DNA is selected by the designer so that the ground link GN has the length r1. The input rotation of DNA about the axis of N drives both the planar wing swing mechanism and spatial wing pitch mechanism.
Our goal is to determine the dimensions of the linkage by solving the SS constraint equations for seven values of the wing pitch function, , see Ref. [16]. The values, , are chosen using Chebyshev spacing along the wing pitch function [21].
The synthesis equations coordinate the rotation of the link DNA, EGC, and BC links so that they simultaneously satisfy, and .
These equations can be reduced to a degree 20 polynomial in terms of one of the following u, v, w, x, y, or z, [15,16], and thus the system of equations has a maximum of 20 unique solutions. This set of equations was solved using Mathematica, which yielded values for the points and . The solutions found are all possible solutions for the set of specified dimensions and angles.
Analysis of the Wing Pitch Mechanism
In order to evaluate the linkage obtained from the synthesis routine, we analyze RSSR Wing pitch mechanism at each input crank position. The goal is to verify that the mechanism moves smoothly through the specified wing swing and wing pitch movements. Every solution from the synthesis is analyzed.
where g and t are determined by the designer.
Solve for using Eq. (5).
The solutions that pass through all precision points on a single branch pass the branching analysis. The passing solutions are then checked for continuity. The range of the input crank angles is divided into 200 intermediate angles, and the links are checked to ensure that they do not change length at the intermediary positions. The solutions are plotted for the entire range of motion of the input crank. An example of a solution that fails in continuity but passes branching is shown in Fig. 4.
Design Methodology
Our design methodology consists of three steps: (i) the required wing swing movement is used to design a planar four-bar function generator; (ii) the coordinate wing pitch movement is used to formulate and solve the SS chain constraints for inverted task positions within specified tolerance zones; and (iii) each solution is analyzed to evaluate its continuous movement through the task positions. This process is iterated with randomly selected task positions within the tolerance zones. Successful design candidates are collected, ranked, and evaluated by the designer.
Example Flapping Wing Mechanism
The precision points of θ and were chosen by applying Chebyshev Spacing on Eq. (1). The calculated input crank angles(θ) are then substituted into the planar linkage equation, see Eq. (6), to obtain the wing swing angle . The values of the tolerance zones for each precision point and are selected by the designer. The values are seen in Table 1.
Table of initial parameter angles for seven precision points
Wing pitch angle requirements (rad) | |||||||
---|---|---|---|---|---|---|---|
i | 1 | 2 | 3 | 4 | 5 | 6 | 7 |
θi | 0 | 0.15 | 1.30 | 3.14 | 4.99 | 6.13 | 6.28 |
0 | ±0.26 | ±0.26 | ±0.26 | ±0.26 | ±0.26 | 0 | |
1.57 | 1.41 | 0.56 | 1.57 | 2.58 | 1.73 | 1.57 | |
0 | ±0.26 | ±0.26 | ±0.26 | ±0.26 | ±0.26 | 0 | |
0.52 | 0.54 | 1.34 | 2.62 | 1.41 | 0.55 | 0.52 |
Wing pitch angle requirements (rad) | |||||||
---|---|---|---|---|---|---|---|
i | 1 | 2 | 3 | 4 | 5 | 6 | 7 |
θi | 0 | 0.15 | 1.30 | 3.14 | 4.99 | 6.13 | 6.28 |
0 | ±0.26 | ±0.26 | ±0.26 | ±0.26 | ±0.26 | 0 | |
1.57 | 1.41 | 0.56 | 1.57 | 2.58 | 1.73 | 1.57 | |
0 | ±0.26 | ±0.26 | ±0.26 | ±0.26 | ±0.26 | 0 | |
0.52 | 0.54 | 1.34 | 2.62 | 1.41 | 0.55 | 0.52 |
Solution for planar crank rocker for given input parameters
Dimensions for planar linkage | |||||
---|---|---|---|---|---|
GN | ND | DE | θ0 | k0 | |
5.00 | 0.39 | 4.99 | 0.45 | −0.04 | −0.04 |
Dimensions for planar linkage | |||||
---|---|---|---|---|---|
GN | ND | DE | θ0 | k0 | |
5.00 | 0.39 | 4.99 | 0.45 | −0.04 | −0.04 |
To align the planar and spatial linkages, the input crank axis of the spatial linkage was placed such that t = −5 in and the output crank axis such that g = 0. Knowing and for and solving Eq. (13) for gives the coordinates of and in the global frame W. The θi values were selected in the range of and the values were selected in the range of . The values were recalculated using Eq. (18) and the values of θi.
where is the wing pitch function from the solution, and n is the number of samples.
Design number 1 was selected because it has the lowest RMS error of the designs with a link ratio less than 10, see Table 3. Results with link ratios less than 10 improved the ease of manufacturing. Its associated inputs are in Table 4 and the output is plotted in Fig. 5.
Table of coordinate solutions
Ranked design candidates | ||||||||
---|---|---|---|---|---|---|---|---|
A coordinates | B coordinates | |||||||
ε | κ | u | v | w | x | y | z | |
0.89 | 2.86 | 2.25 | 0.57 | −3.32 | −1.16 | −1.22 | −1.23 | |
2 | 1.39 | 2.92 | 7.64 | −3.24 | −4.31 | 2.73 | 0.12 | 1.67 |
3 | 1.29 | 3.28 | −7.06 | −4.66 | −4.14 | 1.21 | −2.64 | 0.22 |
4 | 1.37 | 3.42 | −3.98 | −1.87 | −4.12 | 0.72 | −1.72 | 0.47 |
5 | 1.35 | 12.9 | −6.42 | −0.52 | −5.46 | 0.15 | 0.05 | −0.64 |
6 | 1.03 | 16.6 | −3.50 | −0.02 | −4.68 | 0.01 | −0.19 | 0.31 |
Ranked design candidates | ||||||||
---|---|---|---|---|---|---|---|---|
A coordinates | B coordinates | |||||||
ε | κ | u | v | w | x | y | z | |
0.89 | 2.86 | 2.25 | 0.57 | −3.32 | −1.16 | −1.22 | −1.23 | |
2 | 1.39 | 2.92 | 7.64 | −3.24 | −4.31 | 2.73 | 0.12 | 1.67 |
3 | 1.29 | 3.28 | −7.06 | −4.66 | −4.14 | 1.21 | −2.64 | 0.22 |
4 | 1.37 | 3.42 | −3.98 | −1.87 | −4.12 | 0.72 | −1.72 | 0.47 |
5 | 1.35 | 12.9 | −6.42 | −0.52 | −5.46 | 0.15 | 0.05 | −0.64 |
6 | 1.03 | 16.6 | −3.50 | −0.02 | −4.68 | 0.01 | −0.19 | 0.31 |
Table of randomized input parameters
Selected task position for solution 4 | |||||||
---|---|---|---|---|---|---|---|
θ | 0. | 0.40 | 1.17 | 3.35 | 5.24 | 6.09 | 6.28 |
1.57 | 1.19 | 0.80 | 1.83 | 2.76 | 1.81 | 1.57 | |
γ | 0.52 | 0.65 | 1.24 | 2.58 | 1.18 | 0.56 | 0.52 |
Selected task position for solution 4 | |||||||
---|---|---|---|---|---|---|---|
θ | 0. | 0.40 | 1.17 | 3.35 | 5.24 | 6.09 | 6.28 |
1.57 | 1.19 | 0.80 | 1.83 | 2.76 | 1.81 | 1.57 | |
γ | 0.52 | 0.65 | 1.24 | 2.58 | 1.18 | 0.56 | 0.52 |
The mechanism was modeled in SolidWorks, see Figs. 6–8. Figure 9 compares the SolidWorks output, the generated configuration output, and the required wing pitch angle. The SolidWorks model shows some discrepancy due to minor adjustments to the structure that allow it to be physically constructed. The dimensions used are shown in Table 5.
Table of link lengths for the spatial mechanism
Solution 4: dimensions for spatial linkage | |||||
---|---|---|---|---|---|
NA | AB | BC | CG | θs | |
4.48 | 6.57 | 1.44 | 1.27 | 2.70 | 3.66 |
Solution 4: dimensions for spatial linkage | |||||
---|---|---|---|---|---|
NA | AB | BC | CG | θs | |
4.48 | 6.57 | 1.44 | 1.27 | 2.70 | 3.66 |
The physical model of the linkage was built, shown in Figs. 10 and 11. The wing span of the flapping mechanism is 34 in. The straight links and end caps were machined from aluminum or brass tubing. The remaining links and base were made from polycarbonate using additive manufacturing. The revolute joints used brass tubes as bushings while the spherical joints used a compliant joint that used a wire rope to connect the links. The right and left wings are driven by a gear train that connects a motor to their input cranks. The mechanism was driven at a rate of approximately 2–3 Hz, but was not driven any faster due to its size. The focus of this model is to confirm the movement of the proposed linkage. Future research will pursue developing methods to miniaturize and construct a model with a wing span of 5 in.
Conclusion
This paper presents a design methodology for an RSSR mechanism built on the planar four-bar mechanism to coordinate the wing pitch angle with the swing movement of a flapping wing micro-air vehicle. The resulting mechanism is a spatial six-bar linkage that transforms a rotational input to coordinated wing pitch angle and swing angles. The design procedure is demonstrated for wing swing and wing pitch profiles recommended for micro-air vehicles. It results in 29 different designs that meet the requirements, one of which is presented in detail together with a prototype.
Acknowledgment
Benjamin Liu's preparation of the geometric models is gratefully acknowledged.
Funding Data
Division of Civil, Mechanical and Manufacturing Innovation (Grant No. 1636017).