In recent years, unmanned aerial vehicles (UAVs) have gained widespread adoption in both military and civilian domains due to their compact size, high maneuverability, diverse payload capabilities, and low cost. Among these, consumer-grade quadrotor drones have emerged as a mainstream choice for recreational flight experiences and aerial photography, owing to their lightweight design, portability, stability in hover, and affordability. The quadrotor drone system typically comprises a flight controller, airframe structure, power system, and data communication link. While extensive research has focused on flight control algorithms, structural design and analysis for quadrotor drones remain relatively underexplored and lack depth. This study aims to address this gap by presenting a comprehensive structural design, analysis, and optimization methodology for a consumer-grade quadrotor drone, leveraging finite element methods to ensure lightweight, high-strength, and durable performance.
Our primary objective is to develop a full carbon fiber composite structure for a quadrotor drone that meets specific performance criteria: an outer diameter under 0.6 meters, a takeoff weight between 1.8 kg and 2.0 kg (with a payload capacity of at least 300 g), a flight time exceeding 25 minutes, and robust structural integrity. Specifically, under a 2 g load factor, the structure must exhibit a buckling critical load to maximum service load ratio greater than 1.5 and an ultimate static load to maximum service load ratio greater than 2.0. Additionally, the design must accommodate integration of flight control, power, and communication systems while facilitating ease of maintenance. To achieve this, we propose a structured design workflow, as illustrated in Figure 1, which involves initial concept development, finite element modeling, static strength and stability analysis, iterative optimization, and natural vibration analysis.
The initial structural design for the quadrotor drone features an X-symmetric configuration, with the motor mounts integrated into the upper panel to streamline the layout. The upper and lower panels are connected via link plates to form a beam-like structure, enhancing overall stiffness and strength. Motors are positioned at the ends of the arms on the upper panel, with motors 1 and 2 rotating counterclockwise and motors 3 and 4 rotating clockwise. For simplicity in analysis, we consider only the thrust forces generated by the motors and propellers, neglecting torque effects due to their minimal impact. The finite element model was constructed using a combination of four-node rectangular elements and three-node triangular elements at critical junctions, such as the arm-to-center plate connections, to balance computational efficiency and accuracy. Multi-point constraints (MPC) were employed to simulate the interaction between panels and link plates, ensuring realistic load transfer.
Material selection is crucial for achieving lightweight and high-strength properties. We opted for T-300 3K/934 carbon fiber composite, which offers excellent specific strength and stiffness. The material properties, determined through experimental testing, are summarized in Table 1. The density is 1.43 g/cm³, with a single ply thickness of 0.22 mm. To validate our modeling approach, we conducted uniaxial tensile tests on carbon fiber specimens and compared experimental strain data with finite element predictions. The error was less than 5%, confirming the reliability of our model. The layup sequences for the panels are as follows: upper and lower panels use a symmetric sequence of [0°/45°/90°/45°/0°] with reference to the xy-coordinate system, while the link plates employ a [0°/45°/90°/45°/0°] layup oriented along their longitudinal axis.
| Property | Symbol | Value | Unit |
|---|---|---|---|
| Longitudinal Modulus | $$E_1$$ | 62,400 | MPa |
| Transverse Modulus | $$E_2$$ | 62,400 | MPa |
| Longitudinal Tensile Strength | $$F_{1}^{tu}$$ | 628 | MPa |
| Longitudinal Compressive Strength | $$F_{1}^{cu}$$ | 655 | MPa |
| Transverse Tensile Strength | $$F_{2}^{tu}$$ | 607 | MPa |
| Transverse Compressive Strength | $$F_{2}^{cu}$$ | 621 | MPa |
| Shear Modulus | $$G_{12}$$ | 3,400 | MPa |
| Poisson’s Ratio | $$\nu_{12}$$ | 0.33 | – |
| Shear Strength | $$F_{sbs}^{12}$$ | 82.8 | MPa |
We defined two load cases to simulate real-world operating conditions for the quadrotor drone. Load Case 1 represents the takeoff phase, where the lower panel’s central region is fixed, and equal thrust loads are applied at each arm, summing to the maximum takeoff weight of 2 kg. Load Case 2 mimics hover flight, with simply supported conditions at motor mounts and a concentrated load of 2 kg at the center of gravity. These scenarios allow us to assess stress distribution and deformation under critical operational states. The finite element analysis revealed that stress concentrations primarily occur at the junctions between the arms and center plates, as well as near constraint areas. For instance, in Load Case 1, the maximum displacement was 1.05 mm, with an ultimate load factor of 5.90 times the maximum service load and a buckling critical load factor of 12.90. Similarly, Load Case 2 showed a maximum displacement of 1.52 mm, with factors of 6.50 and 11.18, respectively. Under a 2 g overload condition, these factors remain well above design requirements, indicating that the initial structure is overly conservative and heavier than necessary, prompting optimization.
Structural optimization was conducted to minimize weight while maintaining strength and stability margins. We iteratively adjusted the number of carbon fiber plies and modified the geometric configuration of the panels. The optimization process introduced support columns between the upper and lower panels to reduce deformation and address conflicts with component installation space. Table 2 summarizes the optimization data for various models under a 2 g load factor, where the maximum service load is 4 kg. The strength margin is calculated as:
$$\text{Strength Margin} = \frac{\text{Ultimate Load}}{\text{Maximum Service Load}} – 1$$
and the stability margin as:
$$\text{Stability Margin} = \frac{\text{Buckling Critical Load}}{\text{Maximum Service Load}} – 1$$
| Model | Ply Count | Weight (g) | Max Displacement (mm) | Strength Margin | Stability Margin | Flight Time (min) |
|---|---|---|---|---|---|---|
| Initial | 5 | 137.6 | 2.10 | 1.99 | 5.14 | 26.69 |
| Optimization I | 4 | 114.0 | 4.01 | 1.01 | 1.54 | 27.05 |
| Optimization II (3-ply) | 3 | 90.4 | 9.27 | 0.29 | 0.09 | 27.41 |
| Optimization II (with columns) | 3 | 90.4 | 4.65 | 1.93 | 0.26 | 27.41 |
| Final Design | 3 | 108.6 | 2.49 | 1.57 | 1.09 | 27.13 |
The final optimized design for the quadrotor drone incorporates three plies of carbon fiber, with an arm span of 510 mm. Aluminum alloy support columns are added to enhance stiffness, and the arm geometry is refined for better aerodynamic efficiency. This configuration achieves a maximum displacement of 2.49 mm, a strength margin of 1.57, and a stability margin of 1.09, all satisfying the design targets. Notably, the weight is reduced by 21% to 108.6 g compared to the initial model, and the payload capacity is approximately three times the structural weight, significantly improving flight time and performance. The stress distributions in the optimized quadrotor drone structure show reduced concentrations, with peak stresses occurring at strategic reinforcement points, ensuring durability under operational loads.
Natural vibration analysis is essential to prevent resonance, which could compromise the quadrotor drone’s stability and integrity. We performed free modal analysis to determine the inherent dynamic characteristics of the optimized structure. The first four natural frequencies and corresponding mode shapes are listed in Table 3 and described below. The fundamental frequency is 27.863 Hz, corresponding to a torsional vibration about the axis of the arms aligned with the Y-direction. The second mode at 31.867 Hz involves torsion about the X-direction arm axis, while the third mode at 41.553 Hz represents torsional vibration about the vertical Z-axis through the center. The fourth mode at 92.129 Hz exhibits bending vibrations, with X-direction arms deflecting downward and Y-direction arms upward. Given that typical motor operational frequencies range from 50 to 70 Hz, which lies between the third and fourth natural frequencies, the quadrotor drone structure is unlikely to experience resonance, ensuring safe and stable flight.
| Mode Order | Natural Frequency (Hz) | Vibration Mode Description |
|---|---|---|
| 1 | 27.863 | Torsion about Y-direction arm axis |
| 2 | 31.867 | Torsion about X-direction arm axis |
| 3 | 41.553 | Torsion about vertical Z-axis |
| 4 | 92.129 | Bending with X-arms down and Y-arms up |
To further elucidate the structural behavior, we can express key mechanical relationships using mathematical formulations. For instance, the stress components in the composite laminate can be derived from classical lamination theory. The constitutive equation for a laminate under plane stress is given by:
$$\begin{bmatrix} \sigma_1 \\ \sigma_2 \\ \tau_{12} \end{bmatrix} = \begin{bmatrix} Q_{11} & Q_{12} & 0 \\ Q_{12} & Q_{22} & 0 \\ 0 & 0 & Q_{66} \end{bmatrix} \begin{bmatrix} \epsilon_1 \\ \epsilon_2 \\ \gamma_{12} \end{bmatrix}$$
where $$Q_{ij}$$ are the reduced stiffness coefficients, which depend on material properties and ply orientation. For a ply at angle $$\theta$$, the transformed stiffness matrix $$[\bar{Q}]$$ is calculated as:
$$\bar{Q}_{11} = Q_{11} \cos^4\theta + 2(Q_{12} + 2Q_{66}) \sin^2\theta \cos^2\theta + Q_{22} \sin^4\theta$$
$$\bar{Q}_{12} = (Q_{11} + Q_{22} – 4Q_{66}) \sin^2\theta \cos^2\theta + Q_{12}(\sin^4\theta + \cos^4\theta)$$
$$\bar{Q}_{22} = Q_{11} \sin^4\theta + 2(Q_{12} + 2Q_{66}) \sin^2\theta \cos^2\theta + Q_{22} \cos^4\theta$$
$$\bar{Q}_{66} = (Q_{11} + Q_{22} – 2Q_{12} – 2Q_{66}) \sin^2\theta \cos^2\theta + Q_{66}(\sin^4\theta + \cos^4\theta)$$
These equations enable precise prediction of stress states in the quadrotor drone’s carbon fiber components, facilitating accurate strength evaluation. Additionally, the optimization objective function for weight minimization can be formulated as:
$$\min W = \rho \sum_{k=1}^{n} t_k A_k$$
subject to constraints:
$$\sigma_i \leq \sigma_{\text{allowable}}, \quad i = 1,2,\dots,m$$
$$\delta_{\text{max}} \leq \delta_{\text{limit}}$$
$$P_{\text{cr}} \geq \eta P_{\text{service}}$$
where $$W$$ is the total weight, $$\rho$$ is material density, $$t_k$$ and $$A_k$$ are the thickness and area of the k-th ply, $$\sigma_i$$ are stress components, $$\delta_{\text{max}}$$ is maximum displacement, $$P_{\text{cr}}$$ is buckling critical load, and $$\eta$$ is the safety factor. This formulation guided our iterative design adjustments, ensuring that the quadrotor drone meets all performance criteria without unnecessary mass.
The integration of finite element analysis with optimization algorithms has proven highly effective for advancing quadrotor drone design. By simulating various load cases and material configurations, we can rapidly prototype and refine structures in a virtual environment, reducing development time and cost. For the quadrotor drone in this study, the optimized carbon fiber layout not only enhances mechanical performance but also contributes to energy efficiency, as evidenced by the extended flight time. Moreover, the modal analysis provides insights into dynamic response, allowing for tuning of control parameters to mitigate vibrations during flight. As quadrotor drones continue to evolve for applications in delivery, surveillance, and environmental monitoring, such structural optimization techniques will be pivotal in achieving robust, lightweight, and reliable designs.
In conclusion, this study demonstrates a systematic approach to designing and optimizing a carbon fiber composite structure for a consumer-grade quadrotor drone using finite element methods. The final design achieves a weight of 108.6 g, a 21% reduction from the initial model, while satisfying strength and stability requirements under 2 g overload conditions. The natural vibration analysis confirms that resonance risks are minimal, ensuring operational safety. Future work could explore advanced materials like hybrid composites or additive manufacturing techniques to further push the boundaries of quadrotor drone performance. Overall, our methodology offers a scalable framework for developing efficient and durable aerial platforms, contributing to the growing field of UAV technology.