The efficiency of photovoltaic power generation is critically dependent on the surface cleanliness of the solar panel. Accumulated dust and debris can significantly reduce light transmittance, leading to power output losses ranging from 10% to 25% or even higher in arid environments. Maintaining the pristine condition of solar panel arrays, especially in large-scale utility installations, is therefore not merely an operational concern but an economic imperative. While manual cleaning and water spraying are common, they are labor-intensive, inefficient, and water-consumptive. Automated cleaning systems have emerged as a vital solution. A key component of many such robotic cleaners is a mechanical arm or boom that positions the cleaning head across the surface of the solar panel. The design of this structural member directly influences the system’s overall weight, energy consumption, dynamic stability, and ultimately, its practicality and cost-effectiveness.
This article details the structural optimization process for the cross-arm of an automated solar panel cleaning device. The primary objective was to achieve a substantial mass reduction while preserving the structural integrity and dynamic performance necessary for reliable operation. The methodology integrates topology optimization—a mathematical approach for determining the optimal material layout within a given design space—with detailed finite element analysis (FEA) for static, torsional, and modal performance validation.
1. Design Overview of the Solar Panel Cleaning System
The developed cleaning system is designed for deployment on a tracked mobile platform, enabling it to navigate along rows of tilted solar panel arrays. The system comprises two main subsystems: the cleaning head assembly and the supporting boom mechanism. The cleaning head, constructed from lightweight industrial aluminum profiles, houses a motor-driven rotating brush and a water spray nozzle system. It is attached to the distal end of a cantilevered cross-arm.
The boom mechanism is essentially a planar four-bar linkage, consisting of the main cross-arm, a support arm, and a driven arm. This linkage is actuated to extend and retract the cleaning head across the solar panel surface. A critical feature of this kinematic design is that it maintains a constant orientation of the cleaning head, parallel to the plane of the solar panel, throughout its sweeping motion. This ensures consistent contact and cleaning performance.
The initial prototype, while functionally sound, revealed a significant drawback: excessive weight. The cross-arm, being a long cantilevered structure, was identified as the heaviest single component, accounting for approximately half the mass of the entire boom assembly. This heavy mass increases the load on the mobile platform, demands more powerful actuators, and can lead to higher inertial forces and vibrations during operation. Therefore, a targeted lightweight design of the cross-arm was undertaken to improve the system’s overall efficiency and performance.
2. Topology Optimization and Conceptual Redesign
Topology optimization is a powerful computational design tool that solves the problem of distributing a limited amount of material within a predefined design domain to achieve optimal structural performance for given loads and constraints. The fundamental formulation for a stiffness-based topology optimization, aiming to minimize compliance (maximize stiffness), can be expressed as:
$$
\begin{aligned}
& \underset{\rho}{\text{minimize}} & & c(\boldsymbol{\rho}) = \mathbf{U}^T \mathbf{K} \mathbf{U} = \sum_{e=1}^{N} E_e(\rho_e) \mathbf{u}_e^T \mathbf{k}_0 \mathbf{u}_e \\
& \text{subject to} & & \frac{V(\boldsymbol{\rho})}{V_0} = f \\
& & & \mathbf{K} \mathbf{U} = \mathbf{F} \\
& & & 0 < \rho_{\min} \leq \rho_e \leq 1, \quad e = 1, \ldots, N
\end{aligned}
$$
where \(c\) is the structural compliance, \(\boldsymbol{\rho}\) is the vector of design variables (element densities), \(\mathbf{U}\) and \(\mathbf{F}\) are the global displacement and force vectors, \(\mathbf{K}\) is the global stiffness matrix, \(E_e(\rho_e)\) is the Young’s modulus as a function of element density (often via a Solid Isotropic Material with Penalization, SIMP, model), \(\mathbf{u}_e\) and \(\mathbf{k}_0\) are the element displacement vector and stiffness matrix, \(V\) and \(V_0\) are the material volume and design domain volume, \(f\) is the prescribed volume fraction, and \(\rho_{\min}\) is a small lower bound to avoid numerical singularity.
For the cross-arm, the design domain was its original maximum envelope volume. The optimization constraints were defined as follows:
- Boundary Conditions: The four mounting hole surfaces at the rear (proximal end) of the cross-arm were fixed (all degrees of freedom constrained).
- Loading Conditions: Two primary load cases were considered to simulate the forces from the solar panel cleaning head.
- Vertical Bending: A total force of 1200 N (simulating the weight of the cleaning head and ancillary components) was applied downward on the two side faces at the distal end of the cross-arm (600 N per side).
- Horizontal Friction: A force of 300 N was applied horizontally, opposite the direction of travel, to simulate the friction between the rotating brush and the solar panel surface.
- Objective: Minimize structural compliance (maximize stiffness) under the combined loading.
- Constraint: Limit the material volume to 65% of the original design domain (a 35% mass reduction target).
The topology optimization solution, represented by a material density distribution, revealed areas of low stress contribution where material could be removed. The result highlighted six primary regions suitable for material reduction. Interpreting this density map requires engineering judgment to create manufacturable and structurally sound geometries, avoiding sharp corners that induce stress concentrations.
Based on this interpretation, three distinct conceptual redesign schemes were proposed. Schemes 1, 2, and 3 applied similar modifications to five of the six regions but differed significantly in their treatment of the sixth, largest region near the center of the arm’s length. Scheme 1 proposed a conservative removal with a single vertical rib for support. Scheme 2 proposed a more aggressive removal with two vertical ribs. Scheme 3 proposed a removal similar in extent to Scheme 1 but omitted the connecting rib entirely, creating a larger open cavity. The key parameters of the initial concepts are summarized below.
| Scheme | Description of Key Modification (Region 6) | Theoretical Mass Removal (vs. Original) | Primary Design Intent |
|---|---|---|---|
| Original | Solid, unmodified structure. | 0% | Baseline for comparison. |
| Scheme 1 | Moderate opening with one central supporting rib. | 26.5% | Balance weight reduction with stiffness retention. |
| Scheme 2 | Large opening with two supporting ribs. | 29.5% | Maximize weight reduction. |
| Scheme 3 | Moderate opening with no supporting rib. | 26.8% | Maximize weight reduction with simpler geometry. |
3. Comparative Finite Element Analysis of Design Schemes
3.1 Static Bending Analysis
Each proposed scheme was modeled in detail, and a linear static FEA was performed under the same vertical bending load (1200 N downward) and horizontal friction load (300 N) used in the topology optimization. The material was structural steel. The results for maximum deformation, maximum von Mises stress, and final mass were compared against the original design.
The deformation \(\mathbf{u}\) under static load is governed by the linear equation:
$$
\mathbf{K} \mathbf{u} = \mathbf{f}
$$
where \(\mathbf{K}\) is the stiffness matrix of the discretized model and \(\mathbf{f}\) is the applied force vector. The von Mises stress \(\sigma_v\) is calculated from the stress tensor \(\boldsymbol{\sigma}\) to assess yield criteria for ductile materials:
$$
\sigma_v = \sqrt{\frac{(\sigma_{11}-\sigma_{22})^2 + (\sigma_{22}-\sigma_{33})^2 + (\sigma_{33}-\sigma_{11})^2 + 6(\sigma_{12}^2+\sigma_{23}^2+\sigma_{31}^2)}{2}}
$$
| Performance Metric | Design Scheme | |||
|---|---|---|---|---|
| Original | Scheme 1 | Scheme 2 | Scheme 3 | |
| Mass (kg) | 65.02 | 47.78 | 45.82 | 47.57 |
| Mass Reduction | — | 17.24 kg (26.5%) | 19.20 kg (29.5%) | 17.45 kg (26.8%) |
| Max. Displacement X (mm) | 2.409 | 2.9745 | 3.6833 | 4.0337 |
| Max. Displacement Y (mm) | 0.152 | 0.2579 | 0.4583 | 0.382 |
| Max. Displacement Z (mm) | 1.513 | 1.655 | 2.2665 | 2.177 |
| Total Max. Displacement (mm) | 2.838 | 3.1674 | 4.3123 | 4.564 |
| Increase in Total Displacement | — | +0.3294 mm (+11.6%) | +1.4743 mm (+51.9%) | +1.726 mm (+60.8%) |
| Max. von Mises Stress (MPa) | 45.579 | 52.153 | 79.889 | 55.611 |
| Stress Increase | — | +6.574 MPa (+14.4%) | +34.31 MPa (+75.3%) | +10.032 MPa (+22.0%) |
The analysis clearly shows the trade-off between mass reduction and structural stiffness. While all schemes achieved significant weight savings, Scheme 1 offered the best compromise, with a modest 11.6% increase in deflection and a 14.4% increase in stress. Scheme 2, despite achieving the greatest mass reduction, suffered a disproportionate loss of stiffness, with deflection increasing by over 50% and stress by 75%. Scheme 3 performed worse than Scheme 1 in stiffness, demonstrating the importance of the central supporting rib in resisting bending loads applied by the solar panel cleaning head.
3.2 Torsional Analysis
The horizontal friction force applied at the distal end of the cantilevered arm creates a torsional moment about the arm’s longitudinal axis. This torsion is a critical load case as it can lead to twisting and affect the cleaning head’s alignment on the solar panel. The torsional stiffness was evaluated by applying horizontal forces of 600 N and 900 N (simulating higher friction scenarios) at the cleaning head mounting point.
The shear stress \(\tau\) due to pure torsion in a thin-walled section can be approximated by Bredt’s formula:
$$
\tau = \frac{T}{2 A_m t}
$$
where \(T\) is the applied torque, \(A_m\) is the area enclosed by the median line of the cross-section, and \(t\) is the wall thickness. While the cross-arm is not a simple thin-walled tube, this principle highlights that removing material from the central region (reducing \(A_m\)) directly reduces torsional stiffness.
| Performance Metric | 600 N Horizontal Load | 900 N Horizontal Load | ||||||
|---|---|---|---|---|---|---|---|---|
| Orig. | Sch.1 | Sch.2 | Sch.3 | Orig. | Sch.1 | Sch.2 | Sch.3 | |
| Max. Displacement (mm) | 11.234 | 13.2653 | 13.933 | 14.7656 | 16.851 | 19.898 | 20.899 | 22.148 |
| Displacement Increase | — | +18.1% | +24.0% | +31.4% | — | +18.1% | +24.0% | +31.4% |
| Max. Stress (MPa) | 38.452 | 56.461 | 88.905 | 66.024 | 57.678 | 84.692 | 133.36 | 99.036 |
The torsional analysis reinforced the findings from the bending analysis. Scheme 1 consistently showed the smallest increase in deformation and stress among the lightweight options. The percentage increase in displacement remained constant across load levels for each scheme, confirming linear elastic behavior. The maximum stress in all schemes remained below the yield strength of structural steel, but Scheme 2 again showed the largest relative increase, indicating a more pronounced vulnerability to torsional loading during solar panel cleaning operations.
4. Dynamic Modal Analysis of the Optimal Design
Modal analysis determines the inherent vibrational characteristics of a structure—its natural frequencies and mode shapes—which are crucial for assessing dynamic stability and avoiding resonance. The undamped free vibration is governed by the eigenvalue equation:
$$
\left( \mathbf{K} – \omega_i^2 \mathbf{M} \right) \boldsymbol{\phi}_i = \mathbf{0}
$$
where \(\omega_i\) is the \(i\)-th natural circular frequency (\(\omega_i = 2\pi f_i\)), \(\boldsymbol{\phi}_i\) is the corresponding mode shape (eigenvector), and \(\mathbf{M}\) is the mass matrix. The effective modal mass for a given direction indicates how much each mode participates in a dynamic response in that direction.
A comparative modal analysis between the Original design and the selected optimal design (Scheme 1) was conducted, with the proximal mounting holes fixed. The first six natural frequencies and mode shapes were extracted.
| Mode Number | Original Design Frequency (Hz) | Scheme 1 Frequency (Hz) | Frequency Change |
|---|---|---|---|
| 1 | 31.848 | 28.685 | -9.9% |
| 2 | 33.197 | 33.894 | +2.1% |
| 3 | 146.97 | 68.354 | -53.5% |
| 4 | 168.59 | 118.54 | -29.7% |
| 5 | 348.73 | 203.19 | -41.7% |
| 6 | 481.9 | 252.59 | -47.6% |
The reduction in mass and changes in stiffness distribution caused shifts in the natural frequencies. The first two fundamental bending modes (lateral and vertical) saw relatively minor changes. However, the higher-order modes, which involve more complex bending and torsional deformations, experienced significant frequency drops, particularly Modes 3, 5, and 6. This is a direct consequence of removing material from the central region, which reduces both torsional and local bending stiffnesses that resist these higher-order deformations.
The primary mode shapes for both designs were identified from the effective modal mass participation. The first mode is dominated by lateral bending (X-direction translation, Rz rotation). The second mode is dominated by vertical bending (Z-direction translation, Rx rotation). While the general character of the first four mode shapes remained similar between the two designs, the specific deformation patterns in Modes 5 and 6 showed some alterations due to the new geometry. Crucially, the first natural frequency of Scheme 1 at 28.7 Hz is sufficiently high and distinct from common low-frequency excitation sources in a solar panel cleaning environment (e.g., engine vibration, drive motor frequencies), minimizing the risk of resonant vibration during operation.
5. Discussion and Implications for Solar Panel Cleaning Systems
The systematic design optimization process yielded Scheme 1 as the optimal lightweight cross-arm. It achieves a 26.5% mass reduction (17.24 kg) while limiting the increase in static deflection to 11.6% and maintaining all stresses within a safe margin. Although its higher-order natural frequencies decreased, the fundamental frequencies critical for dynamic interaction remain acceptably high.
The success of this optimization has several positive implications for the solar panel cleaning robot:
- Reduced Energy Consumption: A lighter arm requires less energy to accelerate and decelerate during its reciprocating motion across the solar panel, and places a lower static load on the mobile platform’s actuators.
- Potential for Smaller Actuators: The reduced inertial forces may allow for the use of smaller, less expensive motors and drives in the boom mechanism.
- Improved Portability and Deployment: The overall system weight reduction facilitates transportation and setup in large-scale solar panel farms.
- Material Cost Savings: The design uses less raw material, directly reducing manufacturing cost.
The dynamic analysis, particularly the drop in higher-mode frequencies, points to an area for further investigation. Future work should involve a harmonic or transient dynamic analysis to simulate the response of the optimized arm to realistic time-varying loads, such as those induced by the rotating brush engaging with the solar panel surface or by the mobile platform traversing uneven terrain. This would validate that the vibrational amplitudes under operational conditions remain within acceptable limits to ensure cleaning effectiveness and structural longevity.
6. Conclusion
This article presented a comprehensive methodology for the lightweight design of a critical component in an automated solar panel cleaning system. By employing topology optimization to guide the conceptual redesign and rigorous finite element analysis for performance validation, a significant mass reduction of over 26% was achieved for the cross-arm without compromising its core structural integrity. The selected optimal design balances the imperative for weight reduction with the necessary stiffness to withstand the static and dynamic loads encountered during the cleaning of solar panel arrays. The comparative analysis of static displacement, stress, and modal parameters provides a clear engineering rationale for the design choice and establishes a foundation for further dynamic performance evaluation. This work demonstrates that strategic structural optimization is a powerful tool for enhancing the efficiency, economy, and performance of robotic systems dedicated to maintaining the energy output of solar panel installations.