Research Article  Open Access
Alessandro Peloni, Daniele Barbera, Susanna Laurenzi, Christian Circi, "Dynamic and Structural Performances of a New Sailcraft Concept for Interplanetary Missions", The Scientific World Journal, vol. 2015, Article ID 714371, 14 pages, 2015. https://doi.org/10.1155/2015/714371
Dynamic and Structural Performances of a New Sailcraft Concept for Interplanetary Missions
Abstract
Typical square solarsail design is characterised by a central hub with fourquadrant sails, conferring to the spacecraft the classical Xconfiguration. One of the critical aspects related to this architecture is due to the large deformations of both membrane and booms, which leads to a reduction of the performance of the sailcraft in terms of thrust efficiency. As a consequence, stiffer sail architecture would be desirable, taking into account that the rigidity of the system strongly affects the orbital dynamics. In this paper, we propose a new solarsail architecture, which is more rigid than the classical Xconfiguration. Among the main pros and cons that the proposed configuration presents, this paper aims to show the general concept, investigating the performances from the perspectives of both structural response and attitude control. Membrane deformations, structural offset, and sail vibration frequencies are determined through finite element method, adopting a variable pretensioning scheme. In order to evaluate the manoeuvring performances of this new solarsail concept, a 35degree manoeuvre is studied using a feedforward and feedback controller.
1. Introduction
Solar sailing is a promising technology, which allows planning missions otherwise impracticable using traditional propulsion systems. Therefore, the solarsailing concept opens up new avenues for scientific discoveries in many fields of astronautic science, from materials engineering to flight dynamics [1–5]. Due to the continuous and propellantfree thrust, solar sails are mainly studied for orbits with high requirements, such as Earth Polesitter orbits [6–8], orbits at the EarthMoon libration points [9, 10], or new kinds of orbits around the Earth, as the Taranis orbits [11, 12]. Attitude control and optimal steering laws to improve sailcraft performances have recently been studied in several works as well [13–17].
Sailcrafts have a very large and complex structure, typically formed by four petal membranes, which are tensioned to form the square shape using deployable ultrathin composite booms [18–22]. Alternative methods for deploying and tensioning the membrane were also investigated and tested, such as in the case of IKAROS, the firstlaunched solarsail demonstrator [23]. In that case, the four trapezoidal membranes are linked together using spaced strips, which facilitate the folding of the membrane. The solar sail was then deployed and kept extended in a flat shape by the centrifugal force due to the spin of the sailcraft itself. In all mentioned cases, even in the IKAROS demonstrator that did not use booms, the solar sail is visually and physically divided into four membranes and a central hub, which gives the typical Xconfiguration to the spacecraft.
This work proposes a new approach to the solarsail design with a different kind of configuration, in which the classical central bus is divided into four hubs displaced at the corners of the square sail. With this configuration, the sail tensioning can be controlled more easily and the tensioning motors, if any, can be directly placed on the hubs. The membrane tensioning is an important topic, due to the formation of wrinkles or bubbles born from the tensioning or thermal load. This feature is crucial for the masking, shadowing, and thermal issue that may afflict sail performances. From an attitude control point of view, this configuration allows the attitude control thrusters, if any, to be directly mounted on the hubs rather than on the top of a flexible boom. Therefore, the thruster is more likely to be in the nominal position. This feature and the stiffer nature of the architecture itself entail the sail being flatter than in the Xconfiguration. Moreover, this type of architecture can be considered as the unit part of a bigger modular solar sail. On the other hand, the dislocated nature of this configuration increases the moments of inertia of the sailcraft, with a possible decrease of the attitude control efficiency.
This paper investigates the structural and the dynamics performances of the proposed novel configuration of the solar sail. To compare both the structural and the attitude control performances of this new architecture with those available in literature, a 40 m side sail has been considered for the study.
The paper is organised as follows: the new sail configuration is presented in Section 2; a finite element model is reported in Section 3 for the evaluation of maximum outofplane displacement, vibration frequencies, and calculation of the offset between the centre of mass and the centre of pressure. The offset value is used in Section 4 to investigate the solarsail attitude control performances for a 35degree deep space manoeuvre.
2. SolarSail Concept Configuration
The solarsail geometry proposed in this study is a classical square, with the booms on the perimeter of the membrane and the mass of the satellite divided into four parts, collocated on the square’s corners and joined at the booms’ end. The opening sequence from the closedshape launch configuration to the deployed one (Figure 1) is helped by the strain energy stored in the booms. The deployment velocity is a function of the booms’ shapes and the parameters of the deployment mechanism [24, 25].
Table 1 shows the main sailcraft’s characteristics, according to [26], while Figure 2 shows the solarsail reference frame taken into account. The sail mass in Table 1 is computed by considering the same film as in [26], in which the sail has a mass of 6 kg. The boom mass is calculated in the same way. The origin of the body reference frame is set on the geometric centre, roll axis is perpendicular to the sail plane, and pitch and yaw axes are the transverse axes parallel to the booms. is the cone angle between the Sunline direction and the roll axis, is the position vector of centre of mass with respect to the Sun, and is the position vector of the origin of body reference frame with respect to the Sun. The primary attitude control system is based on centreofmass (CM)/centreofpressure (CP) offset due to the shift of 4 ballast masses along the sailcraft’s perimeter. The steadystate offset between the centre of mass and the centre of pressure of the proposed configuration, due to the membrane tensioning, is smaller than the one used in literature for the classical Xshaped solar sail [27–29], as will be discussed in Section 3.3. Because sliding masses do not affect rotations along roll axis, a PulsedPlasmaThrusters (PPTs) system is introduced as well. In the scheme presented there are two pairs of PPTs mounted on two opposite satellites, but four pairs of thrusters can be utilized for redundancy or if greater torque on roll axis is required. It is important to underline that thrusters are mounted on a satellite at the corners of the square sail, instead of at the end of the booms as in Xconfiguration. This means a greater ease of assembly and a major stability of the structure. The total satellite mass considered is about 150 kg [26] and the four dislocated satellite buses have equal masses.

3. Structural Analysis
In this study, we performed nonlinear static analysis, based on Finite Element Method (FEM), adopting the commercial code ABAQUS. The aims of this investigation were the determination of the membrane outofplane deflections caused by the solar pressure, the natural modes of the structures, and the disturbing momentum due to offset between the CP and the origin of the body reference frame, varying the tension applied to the corner of the solar sail. Particular attention was paid when calculating the offset value, which will be used in the dynamics analysis to evaluate the manoeuvring performances of the novel square configuration.
The solar sail is a large thin membrane structure with the bending stiffness negligible compared to the inplane stiffness; thus the membrane cannot carry compressive stress. The flat square shape is controlled by the use of tensioning loads applied at the membrane corners. When the tensioned membrane is exposed to the solar pressure, outofplane large deformations occur and wrinkles can be formed with the origin sets in the corner membrane. The numerical simulation of the wrinkles’ formation is still an open issue, since the wrinkles alter the membrane shape and thus the final performances of the solar sail. Authors investigated numerically the wrinkles’ amplitude tension load dependency, considering the membrane truncated at the corners in order to avoid stress concentrations in those locations [30]. However, the analysis of the wrinkles’ behaviour is not the focus of this current study and will be analysed in a separate work.
In this study, we investigate the structural response of a flat, square, thinfilm membrane tensioned at the corners. Considering the FEM model, the easiest way of pretensioning the membrane is to apply a tension load on the nodes positioned at the four corners. However, this approach can produce difficulties on the convergence of the numerical solution as consequence of the singularities which arise when a single force is applied to a single node of the FEM model. To avoid these singularities, Sleight and Muheim [31] proposed the use of a virtual tension obtained by applying a fictitious thermal load on the sail tensioning cable. This solution, which is a pure mathematical expedient, creates a defined displacement of the corners that correspond to a tensioning load. Furthermore, the extreme dimensions of the solar sails can add convergence problems of the nonlinear analysis. To overcome these problems, the authors in [31] used the cables, modelled as truss elements, when applying the fictitious negative temperatures and producing a membrane tensioning.
Similarly, in our FEM model, we applied a negative temperature to the cable nodes to induce the shrinkage of the wires which connect the membrane and the booms. This shrinkage strains the membrane with wellknown tensioning load. In addition, the prestresses produced in the membrane help the stabilisation of the analysis and the convergence of the solution [32].
The structural simulation was composed of three steps. The first one was a linear displacementtemperature coupled step, during which a negative temperature difference was imposed on the Kevlar cables in order to prestress the membrane. The second step consisted of a nonlinear quasistatic analysis, where the load was applied linearly. In this step, the volumeproportional damping factor was considered to help the convergence of the quasistatic problem. The third step regarded the determination of the natural modes using the Subspace algorithm [33].
At the end of each simulation, the nodal translations and rotations were processed by a Matlab script to calculate the centre of pressure, which was then used to determine the offset between the centre of mass (CM) and centre of pressure (CP).
The considerations concerning the determination of the CP position started from the following equation: where is the projected area of the element, is the solar radiation unit vector, is the normal to the element, and is the position vector of the area element . The normal to the surface element, , varies as a consequence of the nodal rotation; hence an opportune rotation matrix was required to be calculated.
In our system, (1) can be rewritten in discretised form aswhere is the nodal inplane translation. The use of triangular membrane elements avoids outofplane rotations and, as a consequence, the normal vector of the element is constant because of rigid deformation. The offset calculation was performed applying a solar radiation pressure (SRP) force with an angle of 35 degrees between the normal vector to the solarsail plane and the unit vector of the solar radiation.
3.1. Finite Element Model
The geometry of the solar sail was schematised in three main parts: the solarsail membrane given by a square plate with 40 m of edge; the booms represented by four large wires; the tensioning cables given by four small wires. Booms were localised at 0.25 m from the membrane edges, and the tensioning cables were positioned at the corners. The membrane flatness of the classical Xconfiguration is usually increased by dividing the entire sail into several strips. This helps the deployment and also reduces the stresses along the booms [34]. On the other hand, both packaging and jointing of the membrane are challenging tasks and require particular attention. The solarsail architecture presented in this work does not require the strip strategy described above.
The material properties adopted in this work were taken from literature [31, 35] and are summarised in Table 2.

The finite element model was obtained by discretising the geometry with different element types. In particular, the sail membrane was modelled using membrane triangular elements (M3D3) with constant thickness, whereas the linear beam elements (B31) with circular cross section were used for the booms [25, 34, 36]. The tensioning cables were modelled using thermally coupled truss elements (T3D2T) as required by the displacementtemperature coupled step.
Tensioning force was calculated using a Matlab script, which implements (3) considering the Kevlar cable as an isotropic material:where is the crosssection area of the cable, is the Kevlar Young modulus, is the thermal expansion coefficient, and is the imposed temperature.
Assuming that the sail membrane is perfectly reflective, the total pressure load due to SRP is . The pressure was applied in the normal direction with respect to the sail membrane. The boundary conditions were applied on the nodes at the corners of the square membrane, where a beam element and a truss element interlock. In particular, multipoint constraints (MPC) were used to connect the abovementioned nodes to a master node. In this case, the MPC type is a beam providing a rigid link between the master node and the slave ones. The master nodes may have different degrees of freedom, which can vary during the simulation. In our analysis, all master nodes were pinned during the pretensioning step, whereas the translation along the perpendicular direction with respect to the membrane plane and the relative rotation were constrained during the SRP loading step.
Before starting the simulations, a mesh convergence analysis was accomplished in order to investigate the mesh influence on both analysis and final results. In particular, we studied the effects of the number of the membrane elements, considering the maximum displacement, which occurs over the sail during the solarpressure action as a triggering parameter. As discussed above, a pretension load is required to achieve enough stiffness, since the membrane bending rigidity is negligible. The minimum pretensioning load was calculated by a trialanderror technique, establishing that 0.6 N for each tensioning cable is the minimum load necessary to stabilise the membrane in flat configuration. Figure 3 shows the trend of maximum displacement of the membrane as function of the number of elements. The graph shows that over 13,000 elements the value of maximum displacement is constant; thus we adopted this number of membrane elements for the analysis.
3.2. Numerical Results
In this section, we present the results of the maximum outofplane deformation and the first three vibration frequencies due to solar radiation pressure, considering several tensioning cases. In Table 3 the results of the structural analysis for the maximum outofplane displacements achieved varying the tensioning force at the corners of the membrane are reported. As discussed in Section 3, the tensioning force on the membrane corners was obtained by inducing the shrinkage of the cables. This was achieved by imposing a negative temperature as boundary condition at the nodes of the truss elements. However, it is worth to note that the temperature boundary conditions reported in Table 3 are a fictitious thermal load, which is a numerical expedient to generate an effective tension applied to the membrane corners (see (3)).

As shown in Table 3, the first tension load at which the analysis converges is given by 0.6 N. As expected, the final outofplane displacement decreases with an increase of the pretensioning load. Figure 4 shows the outofplane displacements obtained with the minimum tensioning force (on the left) and the maximum tensioning force (on the right), in the case of a maximum thrust during a 35degree manoeuvre. The two displacement distributions are similar, but the maximum displacement value is largely reduced in case of tensioning at 18 N. The load required to tension the sail membrane properly is an important design key factor, because it affects the wrinkles’ formation on the sail membrane [30], and it is limited by structural stability of the booms.
The last part of the FEM analysis investigates the natural vibration modes of the solarsail membrane for the different pretensioning cases. The results are reported in Table 4, where it can be observed that the vibration modes shift to higher frequencies with the increasing of the pretensioning load. These results can be related to the membrane stress state, which is directly influenced by the intensity of the tensioning load. An increase of such loads corresponds to an increase of the vibration frequency. The total displacements associated with the first four modes of the 0.6 N and 18 N tensioning cases are shown in Figures 5 and 6, respectively. Comparing the images in Figures 5 and 6, it can be observed that the shapes of the modes encountered in the 0.6 N tensioning case are different from those of the 18 N tensioning one. This difference may be explained by the loose state of the membrane in case of tension at 0.6 N. Further, the 0.6 N tension load is the lowest value of the tension force to reach the numerical convergence of the solution in the static analysis, but this value may add some uncertainties to the dynamic analysis. In fact, we noted that, in the case of pretensioning at 0.6 N, the first four modes of the solarsail membrane presented shapes similar to ones reported in Figure 6, which are the typical shape modes of a square membrane.

3.3. Determination of the Structural Offset
The disturbing offset was determined using a Matlab script for each tensioning case. In this calculation, the model was simplified assuming that the booms could withstand the axial loads due to the sail tensioning and that the wrinkles of the membrane were negligible. This approach allowed us to investigate the maximum tensioning load required to reduce as much as possible the offset. The offset position on the sail plane is greatly influenced by the SRP modelling. The solar vector is given by two components, one along the negative direction and one along the negative direction. Because of the symmetry of the structures, the same conclusion can be obtained considering direction instead of direction. Results of these calculations are reported in Table 5, where it can be noted that an increased load grants an offset reduction due to the increasing of the inplane membrane stiffness. In particular, the solarsail configuration proposed in this work has a smaller offset than the one for the classical Xshape configuration reported in literature [26].

The percentage change of the offsets extrapolated from the FEM analysis was taken into account for the attitude dynamics analysis. In particular, the dynamics analysis considered only the worst case offset scenario, which is the minimum applicable tensioning force (0.6 N) and is the 0.0128% of the sail edge size. The new sail concept is stiffer and guarantees a smaller offset than the Xshape one, which is 0.25% of the sail edge size [26].
4. Dynamics Sailcraft Performances
Due to dislocated mass, the proposed sailcraft’s configuration is characterised by moments of inertia greater than the classical one. Therefore the study of the performances for an attitude manoeuvre is essential to understand whether this architecture can be a valid flight configuration. A 35degree manoeuvre is taken into account to compare these performances with literature data. Note that, due to slow dynamics, classical interplanetary missions require smallamplitude manoeuvres per day [37, 38], while a fast manoeuvre is required only for “nonclassical” interplanetary missions [39, 40]. A body reference frame is considered for modelling sailcraft attitude dynamics, as described in Section 2.
The primary attitude dynamics control is performed with four ballast masses (), of which masses 1 and 3 shift only along axis and masses 2 and 4 shift only along axis (Figure 7). The coordinates of the masses are reported in (4), where is the sail side length and The sailcraft attitude dynamics is as follows [41]:where is the angular momentum of the system referred to the origin of body reference , is the total sailcraft mass, is the absolute acceleration of point , is the torque referred to , and is the position vector of the centre of mass from the reference point , as shown inThe total angular momentum is given by the angular momentum of the bus, booms, and membrane and the angular momentum of masses for attitude control . In particular isand iswhere is the matrix of inertia due to th control mass. The term which appears in (7)(8) is the angular velocity vector of the spacecraft expressed asThe absolute acceleration of the reference point is expressed by the absolute acceleration of the CM aswhere is the gravitational force vector and, considering the simplified solarsail force model [18], the solar radiation pressure force vector is , where is the value of SRP at 1 AU , is the sail area, is the solarsail efficiency factor, and is the angle between Sunline direction and roll axis, as described in Section 2. The torque relative to the origin is given by the gravitational force and SRP force, as shown inwhere is the Sun gravitational constant, is the distance between CM and the th point of the sailcraft, and is the disturbance offset vector. Since ballast masses 1, 3 and 2, 4 are coupled, (5) can be rewritten through the following scalar:where are SRP force components in body reference frame and are the components of expressed in body reference as and depend on the set of rotations chosen for the attitude representation. Let us call , respectively, the roll, pitch, and yaw angles of the spacecraft relative to the orbital reference frame, obtained by a rotational sequence of from the orbital to the body reference frame. The kinematics equations arewhereis the force per unit mass in the orbital angular momentum direction . In order to achieve the desired manoeuvre, a combination of feedforward and feedback control, as described in [42, 43], is used. In Sections 4.1 and 4.2 feedforward and feedback methods are briefly presented.
4.1. The Feedforward Controller
Feedforward control is based on a parameterisation of a desired manoeuvre, expressed as a nthorder polynomial in the generic angle . The order of polynomial depends on the boundary conditions and it should not be too big to reduce wandering phenomena. A seventhorder polynomial has been considered in this study as follows:where is the desired angle of manoeuvre and is the nondimensional time. is the final time after the manoeuvre. In order to find the coefficients in (15), the boundary conditions for are listed in Table 6.

According to the boundary conditions in Table 6, the parameterised manoeuvre angle is expressed by the following seventhorder polynomial:In order to design simply the feedforward controller, the following assumptions are made [42]:(i)Inertia matrix is diagonal and constant, not affected by the position of control masses.(ii)SRP force is the only force acting on the sailcraft.(iii)The centre of mass lies on sail plane .With the hypothesis above, Euler equation is simply given bywhere is the control torque:Using (17)(18) and defining as the angle between Euler axis and pitch axis, the feedforward control law for masses is given byTo comply with geometrical boundaries, the dynamical scaling of the manoeuvre time is used [44]:where and are the maximum shift of control masses required along the axis and axis, respectively.
4.2. The Feedback Controller
The pitch/yaw feedback control is based on the error between the desired manoeuvre and the one carried out by the sailcraft. The feedback control logic is in ProportionalIntegralDerivative (PID) form as below [26]:where is the error between desired angle and real one and , , and are the derivative, proportional, and integral gain, respectively. This control can be decoupled in each axis and gains can be determined as SingleInputSingleOutput (SISO) problem. Of course, whenever tuning a gain, the system response changes and the best gains are iteratively set.
A rate limiter and a saturation limit are added to simple PID controller, because of mechanical and geometrical boundaries. These boundaries arewhere is the actuator time constant taken into account, according to the value in [26]. The roll feedback control is performed with onoff controllers that work when the tolerance on roll angle is exceeded. The controller chosen is a set of 4 PPTs positioned coupled on 2 opposite satellites and the average thrust of each PPT chosen is 150 μN [45]. For this study, a required roll angle of 0 degrees with a tolerance of ±0.1 degrees is set, so that the thrusters switch on when the roll angle exceeds this value.
4.3. Numerical Results
Numerical simulations are carried out in order to verify the performances of the proposed configuration. As reported in previous sections, a feedforward controller was used to generate the desired manoeuvre over time and a feedback controller with PID logic was set to control the nonmodelled trends in feedforward controller. The characteristics of the sailcraft are those reported in Table 1, while velocities and accelerations of masses in (12) are considered null [26, 28]. The attitude is represented by the rotational matrix which transforms body into orbital reference frame:In order to evaluate the performances of the proposed configuration, a 35degree manoeuvre in a circular planar Earthlike orbit is considered. In Figure 8(a) a 35degree yaw manoeuvre with disturbance offset on axis is schematically represented; in Figure 8(b) a 35degree pitch manoeuvre with disturbance offset on axis is schematically represented.
(a)
(b)
First, a 35degree yaw manoeuvre without offset () is performed. Time histories of Euler angles and control masses over time are shown in Figures 9 and 10.
The manoeuvre is completed in less than 3 hours and roll and pitch angles do not change from initial conditions. Figure 10 shows that the feedforward control ensures the success of the manoeuvre, though the real position of the sliding masses differs slightly from the predicted one, as can be seen around 2.5 hours, due to the nonmodelled forces in the feedforward controller. However, no evident differences between real angles and predicted ones are visible by Figure 9.
The analysis with offset takes into account only the maximum offset calculated in Section 3.3, in order to have a worstcase analysis. Figures 11 and 12 show a 35degree yaw manoeuvre with 0.005 m offset (red curve) and with the literature one of 0.1 m (black curve).
Similar to the case without offset, roll and pitch angles do not change from initial conditions and real and predicted angles overlap. The manoeuvre with the offset presented in this paper is performed in less than 3 hours, as well as in the case without offset. On the other hand, the case with the literature offset performs the manoeuvre in more than 3.5 hours. The reason for this gap can be found in the differences between the steadystate positions of control masses. The steadystate position of each sliding mass can be simply obtained by (19), as shown inAs shown in Figure 12 and in (25), the steadystate position of the sliding masses on axis is m with the literature offset of 0.1 m; on the other hand, for an offset value of 0.005 m, the steadystate position of the control masses on axis is only .
As shown from Figures 9 to 12, no roll control is necessary during a 35degree yaw manoeuvre with offset on axis, because the roll angle is null during the entire manoeuvre. Figures 13 and 14 show that for a 35degree pitch manoeuvre with a disturbance offset on the same axis a roll control is necessary.
Figure 13 shows that in a pitch manoeuvre the roll angle decreases slowly, so that only after about 5 hours does the roll controller act, due to the threshold set to 0.1 degrees. For missions that require different attitude accuracy, this threshold can be set to different values. Excluding the time history of roll angle, the manoeuvre is completed as well as in the previous examples.
Comparing these manoeuvres with those in literature, a 35degree yaw manoeuvre presented in [26] and a 30degree pitch manoeuvre presented in [28] are both completed in about 2 hours. The manoeuvre presented in [42] for a perfectly reflecting solar sail is much faster than the one presented here, due to the higher moments of inertia of the proposed architecture. As mentioned above, high moments of inertia are a disadvantage for an attitude manoeuvre and one of the aims of this study is to ensure the manoeuvrability of this configuration. However, in interplanetary trajectories a fast manoeuvre is not always required and the configuration presented has performances compatible with the most challenging interplanetary missions.
5. Conclusions
In this work, a new solarsailing architecture has been proposed and structural and dynamics performances have been investigated. The structural analysis has shown that the new solar sail is stiffer than the centralhub sail. This characteristic produces smaller outofplane deformation and, as a consequence, a reduction of the disturbing centreofmass/centreofpressure offset. Furthermore, the proposed configuration can easily manage this offset with an opportune sail tensioning load. The applied tensioning load depends on the booms’ stiffness and it can be foreseen that it affects greatly the membrane wrinkling and that it increases with the load. In the proposed configuration, the tensioning motors, the sliding masses motors, and PPTs can find a good location in the corner hubs, whereas in the centralhub sail they are positioned at the end of the booms. This can allow the designing of a real shape and vibration membrane control, reducing the offset with respect to the offset reported in literature for the centralhub sail. Furthermore, the cornertensioned configuration produces more important effects such as the reduction of the masking, the shadowing, and the local thermal problems that can be controlled during the flight. In addition, it is worth noting that in the classical Xconfiguration the effective area for a 40 m side sail is 1200 m^{2} instead of 1600 m^{2}, due to the membrane deformation and the fact that around the booms there is no membrane. On the contrary, in the proposed configuration the disposition of the booms along the perimeter and the deformation of the membrane allow one to consider a major effective area, closer to nominal value of 1600 m^{2}. Starting from these results, numerical simulations of the attitude manoeuvres demonstrated that the proposed architecture gives good performances, despite the large moments of inertia. It was shown that a 35degree manoeuvre can be completed in less than 3 hours, according to the usual requirements for interplanetary missions. The new solarsail concept is proposed for an interplanetary mission, but advantages of this configuration are useful also for planetary missions. In general, the principal constraint that must be taken into account is the velocity of change attitude manoeuvres. In fact, for a planetary mission, around the nodal line a fast change on the sail attitude can be necessary; this event can be critical for this configuration and the performances must be evaluated case by case. On the contrary the higher moment of inertia can be favourable for missions where the sail attitude must remain constant (but this is valid also for the interplanetary case). In all cases studied, the disturbing torque, caused by structural offset, determines the steadystate positions of the sliding masses. As a consequence, the small offset value of this sail configuration guarantees a great increase in manoeuvrability.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
References
 M. MacDonald and C. R. McInnes, “Solar sail science mission applications and advancement,” Advances in Space Research, vol. 48, no. 11, pp. 1702–1716, 2011. View at: Publisher Site  Google Scholar
 C. Circi, “Mars and mercury missions using solar sails and solar electric propulsion,” Journal of Guidance, Control, and Dynamics, vol. 27, no. 3, pp. 496–498, 2004. View at: Publisher Site  Google Scholar
 B. Dachwald and W. Seboldt, “Multiple nearEarth asteroid rendezvous and sample return using first generation solar sailcraft,” Acta Astronautica, vol. 57, no. 11, pp. 864–875, 2005. View at: Publisher Site  Google Scholar
 G. Mengali, A. A. Quarta, C. Circi, and B. Dachwald, “Refined solar sail force model with mission application,” Journal of Guidance, Control, and Dynamics, vol. 30, no. 2, pp. 512–520, 2007. View at: Publisher Site  Google Scholar
 G. Aliasi, G. Mengali, and A. A. Quarta, “Passive control feasibility of collinear equilibrium points with solar balloons,” Journal of Guidance, Control, and Dynamics, vol. 35, no. 5, pp. 1657–1661, 2012. View at: Publisher Site  Google Scholar
 M. Ceriotti and C. R. McInnes, “Systems design of a hybrid sail polesitter,” Advances in Space Research, vol. 48, no. 11, pp. 1754–1762, 2011. View at: Publisher Site  Google Scholar
 M. Ceriotti and C. R. McInnes, “Generation of optimal trajectories for Earth hybrid pole sitters,” Journal of Guidance, Control, and Dynamics, vol. 34, no. 3, pp. 847–859, 2011. View at: Publisher Site  Google Scholar
 J. Heiligers, M. Ceriotti, C. R. McInnes, and J. D. Biggs, “Displaced geostationary orbit design using hybrid sail propulsion,” Journal of Guidance, Control, and Dynamics, vol. 34, no. 6, pp. 1852–1866, 2011. View at: Publisher Site  Google Scholar
 J. Simo and C. R. McInnes, “Solar sail orbits at the EarthMoon libration points,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 12, pp. 4191–4196, 2009. View at: Publisher Site  Google Scholar  MathSciNet
 J. Simo and C. R. McInnes, “Asymptotic analysis of displaced lunar orbits,” Journal of Guidance, Control, and Dynamics, vol. 32, no. 5, pp. 1666–1670, 2009. View at: Publisher Site  Google Scholar
 P. Anderson and M. Macdonald, “Extension of highly elliptical earth orbits using continuous lowthrust propulsion,” Journal of Guidance, Control, and Dynamics, vol. 36, no. 1, pp. 282–292, 2013. View at: Publisher Site  Google Scholar
 P. Anderson and M. Macdonald, “Static, highly elliptical orbits using hybrid lowthrust propulsion,” Journal of Guidance, Control, and Dynamics, vol. 36, no. 3, pp. 870–880, 2013. View at: Publisher Site  Google Scholar
 G. Mengali and A. A. Quarta, “Solar sail trajectories with piecewiseconstant steering laws,” Aerospace Science and Technology, vol. 13, no. 8, pp. 431–441, 2009. View at: Publisher Site  Google Scholar
 B. Dachwald, “Optimization of verylowthrust trajectories using evolutionary neurocontrol,” Acta Astronautica, vol. 57, no. 2–8, pp. 175–185, 2005. View at: Publisher Site  Google Scholar
 T. Ingrassia, V. Faccin, A. Bolle, C. Circi, and S. Sgubini, “Solar sail elastic displacement effects on interplanetary trajectories,” Acta Astronautica, vol. 82, no. 2, pp. 263–272, 2013. View at: Publisher Site  Google Scholar
 C. Circi, “Threeaxis attitude control using combined gravitygradient and solar pressure,” Proceedings of the Institution of Mechanical Engineers, Part G: Journal of Aerospace Engineering, vol. 221, no. 1, pp. 85–90, 2007. View at: Publisher Site  Google Scholar
 G. Mengali and A. A. Quarta, “Nearoptimal solarsail orbitraising from low earth orbit,” Journal of Spacecraft and Rockets, vol. 42, no. 5, pp. 954–958, 2005. View at: Publisher Site  Google Scholar
 C. R. McInnes, Solar Sailing: Technology, Dynamics and Mission Applications, Springer Praxis Publishing, Chichester, UK, 2004.
 L. Johnson, R. Young, E. Montgomery, and D. Alhorn, “Status of solar sail technology within NASA,” Advances in Space Research, vol. 48, no. 11, pp. 1687–1694, 2011. View at: Publisher Site  Google Scholar
 L. Johnson, M. Whorton, A. Heaton, R. Pinson, G. Laue, and C. Adams, “NanoSailD: a solar sail demonstration mission,” Acta Astronautica, vol. 68, no. 56, pp. 571–575, 2011. View at: Publisher Site  Google Scholar
 V. Lappas, N. Adeli, L. Visagie et al., “CubeSail: a low cost CubeSat based solar sail demonstration mission,” Advances in Space Research, vol. 48, no. 11, pp. 1890–1901, 2011. View at: Publisher Site  Google Scholar
 U. Geppert, B. Biering, F. Lura, J. Block, M. Straubel, and R. Reinhard, “The 3step DLR–ESA gossamer road to solar sailing,” Advances in Space Research, vol. 48, no. 11, pp. 1695–1701, 2011. View at: Publisher Site  Google Scholar
 Y. Tsuda, O. Mori, R. Funase et al., “Achievement of IKAROS—Japanese deep space solar sail demonstration mission,” Acta Astronautica, vol. 82, no. 2, pp. 183–188, 2013. View at: Publisher Site  Google Scholar
 C. Sickinger, L. Herbeck, and E. Breitbach, “Structural engineering on deployable CFRP booms for a solar propelled sailcraft,” Acta Astronautica, vol. 58, no. 4, pp. 185–196, 2006. View at: Publisher Site  Google Scholar
 A. Stabile and S. Laurenzi, “Coiling dynamic analysis of thinwalled composite deployable boom,” Composite Structures, vol. 113, no. 1, pp. 429–436, 2014. View at: Publisher Site  Google Scholar
 B. Wie and D. Murphy, “Solarsail attitude control design for a sail flight validation mission,” Journal of Spacecraft and Rockets, vol. 44, no. 4, pp. 809–821, 2007. View at: Publisher Site  Google Scholar
 B. Wie, D. Murphy, M. Paluszek, and S. Thomas, “Robust attitude control systems design for solar sail, part 1: propellantless primary ACS,” in Proceedings of the AIAA Guidance, Navigation, and Control Conference and Exhibit, Providence, RI, USA, 2004. View at: Publisher Site  Google Scholar
 S. N. Adeli, V. J. Lappas, and B. Wie, “A scalable busbased attitude control system for Solar Sails,” Advances in Space Research, vol. 48, no. 11, pp. 1836–1847, 2011. View at: Publisher Site  Google Scholar
 B. Wie, “Solar sail attitude control and dynamics, part 1,” Journal of Guidance, Control, and Dynamics, vol. 27, no. 4, pp. 526–535, 2004. View at: Publisher Site  Google Scholar
 W. Wong and S. Pellegrino, “Wrinkled membranes III: numerical simulations,” Journal of Mechanics of Materials and Structures, vol. 1, no. 1, pp. 63–95, 2006. View at: Publisher Site  Google Scholar
 D. W. Sleight and D. M. Muheim, “Parametric studies of square solar sails using finite element analysis,” in Proceedings of the 45th AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics & Materials Conference, AIAA20041509, Palm Springs, Calif, USA, April 2004. View at: Publisher Site  Google Scholar
 S. C. Gajbhiye, S. H. Upadhayay, and S. P. Harsha, “Free vibration analysis of flat thin membrane,” International Journal of Engineering Science and Technology, vol. 4, no. 8, pp. 3942–3948, 2012. View at: Google Scholar
 I. Abaqus, Abaqus Analysis User's Manual, 2010.
 J. M. Fernandez, V. J. Lappas, and A. J. DatonLovett, “Completely stripped solar sail concept using bistable reeled composite booms,” Acta Astronautica, vol. 69, no. 12, pp. 78–85, 2011. View at: Publisher Site  Google Scholar
 D. W. Sleight, T. Mann, D. Lichodziejewski, and B. Derbes, “Structural analysis and test comparison of a 20meter inflationdeployed solar sail,” in Proceedings of the 47th AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics and Materials Conference, pp. 1408–1423, May 2006. View at: Google Scholar
 S. Laurenzi, D. Barbera, and M. Marchetti, “Buckling design of boom structures by FEM analysis,” in Proceedings of the 63rd International Astronautical Congress (IAC '12), pp. 6367–6371, Naples, Italy, October 2012. View at: Google Scholar
 A. Bolle and C. Circi, “Solar sail attitude control through inplane moving masses,” Proceedings of the Institution of Mechanical Engineers, Part G: Journal of Aerospace Engineering, vol. 222, no. 1, pp. 81–94, 2008. View at: Publisher Site  Google Scholar
 G. Colasurdo and L. Casalino, “Optimal control law for interplanetary trajectories with nonideal solar sail,” Journal of Spacecraft and Rockets, vol. 40, no. 2, pp. 260–265, 2003. View at: Publisher Site  Google Scholar
 G. Mengali, A. A. Quarta, D. Romagnoli, and C. Circi, “H2reversal trajectory: a new mission application for highperformance solar sails,” Advances in Space Research, vol. 48, no. 11, pp. 1763–1777, 2011. View at: Publisher Site  Google Scholar
 G. Vulpetti, “3D highspeed escape heliocentric trajectories by allmetallicsail lowmass sailcraft,” Acta Astronautica, vol. 39, no. 1–4, pp. 161–170, 1996. View at: Publisher Site  Google Scholar
 B. Wie, “Solar sail attitude control and dynamics, part 2,” Journal of Guidance, Control, and Dynamics, vol. 27, no. 4, pp. 536–544, 2004. View at: Publisher Site  Google Scholar
 D. Romagnoli and T. Oehlschlägel, “High performance two degrees of freedom attitude control for solar sails,” Advances in Space Research, vol. 48, no. 11, pp. 1869–1879, 2011. View at: Publisher Site  Google Scholar
 C. Scholz, D. Romagnoli, B. Dachwald, and S. Theil, “Performance analysis of an attitude control system for solar sails using sliding masses,” Advances in Space Research, vol. 48, no. 11, pp. 1822–1835, 2011. View at: Publisher Site  Google Scholar
 B. Siciliano, L. Sciavicco, L. Villani, and G. Oriolo, Robotics—Modelling, Planning and Control, Advanced Textbooks in Control and Signal Processing, Springer, London, UK, 2009.
 B. Wie, D. Murphy, M. Paluszek, and S. Thomas, “Robust attitude control systems design for solar sails, part 2: micropptbased secondary acs,” in Proceedings of the AIAA Guidance, Navigation, and Control Conference and Exhibit, Providence, RI, USA, 2004. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2015 Alessandro Peloni et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.