# Error Assessment Model for the Inverse Kinematics Problem for Stewart Parallel Mechanisms for Accurate Aerospace Optical Linkage

Authors: J. J. Hernández-Gómez; I. Medina; C. R. Torres-San Miguel; A. Solís-Santomé; C. Couder-Castañeda; J. C. Ortiz-Alemán; J. I. Grageda-Arellano
Journal: Mathematical Problems in Engineering
 (2018)

---

## Abstract

In this work we develop a mathematical model to estimate the error for inverse kinematics problem for Gough-Stewart parallel mechanisms. We propose the estimation error method to include manufacture, assembly, backlash, and sensing errors. We provide the error transmission matrices for the length of each leg of the hexapod, which permits evaluation of the accuracy error in the position of each one, given a desired position and orientation of the mobile platform. We also present numerical modelling in order to estimate the accuracy of the methodology herein proposed, for specific attitude operations corresponding to performing a successful ground-LEO nanosatellite optical link. In such a case, we were able to provide the required tolerances for the actuators in order to guarantee an orientation precision requirement of the order of milliradians.


---

## Body

## 1. Introduction

Space-based laser communications (Lasercom) have the potential to transform scientific, defense, and commercial spacecraft communications platforms. Compared with traditional radio frequency (RF) communications, Lasercom offers higher bandwidth, reduced size and mass of transceivers, lower power consumption, immunity from interference, and higher gain and also avoids the significant regulatory hurdles of radio frequency allocation [1]. Nanosatellites, in particular CubeSats, have great potential to benefit from Lasercom as they are heavily constrained in size, weight, and power [2].However, the benefits of Lasercom also come with a unique set of challenges. Two of the main problems are the atmospheric loss, since Lasercom frequencies are very susceptible to atmospheric effects, and the pointing, acquisition, and tracking (PAT) subsystem. While narrow transmit beamwidths increase Lasercom link efficiency, the tradeoff is the need for very precise PAT.Lasercom systems ought to align the optical line-of-sight (LOS) very precisely, with typical systems requiring error as small as submicroradians [3]. Typical pointing accuracies for laser satellite communications (LSC) systems vary from 1 to 200 μrad; however most typical systems have a pointing accuracy within the range of 10 to 50 μrad. In the tracking phase, the current values of accuracies typically reach 1 to 3 μrad [4, 5].Lasercom links from low Earth orbit (LEO) to ground present several unique challenges for PAT. The satellite will have to slew rapidly to track the ground station, placing stress on the attitude determination and control subsystem. The fine tracking mechanism to locate a satellite-like spacecraft requires a point-ahead/look-behind calculation, as depicted in Figure1, considering the distance between Terminal 1 and Terminal 2 at t0 as the nominal orbital distance. However, slight variations of the location of Terminal 2 imply the need to consider a point-ahead or look-behind angle (θ). In order to point at Terminal 2 at t0, for instance, the point-ahead ought to be performed at time t0-Δt, where Δt is the time that the link signal takes (travelling at speed of light c) to travel from Terminal 1 to Terminal 2. Following this principle and the trigonometry from Figure 1, the point-ahead/look-behind angle would be given by(1)θ=2arctan⁡vrelΔtcΔt≈2VrelcFigure 1
Point-ahead/look-behind angle as a function of relative velocity [6].In (1), the last term arises from the fact that the link range is much greater (orbit altitude) than the change in position (Figure 1 does not represent real distances and angles). In this approximation, the point-ahead angle depends only on the relative velocity between the terminals, Vrel, and on c. Thus, the point-ahead angle given by (1) will be relatively large due to high orbital velocities. Furthermore, ground station sight of view endures less than 10 minutes for LEO satellites, so PAT subsystem ought to quickly acquire and track the ground station to maximise data transmission [6], while keeping great accuracy during the optical link.Many satellite-to-ground optical communications experiments have been tested by using the optical ground station (OGS) developed by the European Space Agency (ESA) since 2000 [7]. This includes the SILEX project that connected the OGS with the satellite ARTEMIS [8] as well as the SILEX compatible laser communication terminal named LUCE, developed by Japanese Space Exploration Agency (JAXA). The later was successfully tested in 2003 from ESA’s optical ground station in a bidirectional space to ground link with ARTEMIS [9]. Other experimental links with TerraSAR-X and NFIRE satellites were tested for only PAT performances [10].More recently, the laser communication experiment between the satellite SOCRATES and the OGS MeO—belonging to the Observatoire de la Côte d’Azur, located in France—has obtained important data from satellite-to-ground links [11]. JAXA has also been testing ground-to-satellite optical links using the NICT optical ground station in Koganei, Tokyo, since 2006 [12, 13]. Nevertheless, all of these optical aerospace communication experiments have been tested on conventional stationary observatories with two degrees of freedom and limited speed, mostly limiting tests to satellites located in higher orbits.In order to overcome the problem of atmospheric effects in optical link, we propose the development of a transportable optical ground station. Nevertheless, the mobility of such ground station would increase the requirements of accuracy in its PAT, particularly in the case of ground to LEO links to nanosatellites, due to their high orbital velocity and reduced size. So in order to achieve greater accuracy than conventional PAT systems, we propose to use a Gough-Stewart platform as PAT at ground station.Gough-Stewart platform is a parallel mechanism with six degrees of freedom, which is very popular for positioning and high accuracy operations [14–19] because of their great load capacity, as well as their particular high speeds and acceleration, while reducing vibrations. Those characteristics make Gough-Stewart manipulator ideal for optical aerospace linkage. So for a hexapod to provide such accuracy as that required for aerospace optical link from mobile ground station to LEO nanosatellites, the control of the error in the positioning and orientation of the mobile platform in the Gough-Stewart manipulator is fundamental. Nevertheless, due to the high nonlinearity of its kinematics and dynamics, it is commonly a difficult task to cast theoretical models that accurately describe errors in positioning and orientation of the platform arising, for instance, from manufacturing, assembly, backlash, and sensing. Furthermore, studies focused on determining mathematical models for error in this kind of parallel manipulators are usually casted in an empiric way [20].Many authors have explored the error and accuracy of parallel mechanisms using different techniques. Altuzarra et al. [21] provide a method to locate regions of the workspace with discontinuities for a planar parallel robot, addressing the error between the pin and the hub of a joint with clearance. Also, stiffness modelling method for overconstrained systems has been used to formulate workspaces as close as possible to a predefined velocity in flexible links [22]. Ding et al. [23] quantify kinematics stiffness error through a simplified method, which combines both analytic and numerical modelling (FEA), providing ultrahigh accurate posing applications. Chebbi et al. [24] focus their study on the effects of the clearance in the joints on the parallel robot accuracy, which yields to an analytic predictive model for the pose error. Moreover, Zhang and Gao [25] design kinematic calibration of parallel manipulators assisted though cooperative coevolutionary neural networks (CCNN) to establish a nonlinear relationship between joint variables and the measured pose of the end-effector.In the case of Gough-Stewart parallel mechanism, some error analysis methods based on differential operators [26] consider errors in the structural dimensions due to the manufacturing and assembling tolerances. Lin et al. [26] demonstrated that such structural errors greatly affect the obtained position of the mobile platform at certain angles. Ropponen and Arai [27] considered the kinematics model taking into account errors in the position of the articulations, as well as action and reaction ones, providing the accuracy in the final positioning of the mobile platform. On the other side, Wang and Ehmann [28] present an error model and accuracy analysis of Stewart platform by focusing on the sources of errors of ball joints, ball joints supports, and motion of ball joints and considering an error model for a one-DOF actuator, while Zhou et al. [29] propose a chain error propagation model based on the differential of D-H matrix, where the chains are treated as series single-opened-chains.Nowadays, each error source can be calculated separately, but up to now there is no theoretical mathematical model able to predict the final position of the legs of the actuator considering all the error sources at the same time. In this work we provide a methodology to calculate the error of position of each leg of the hexapod obtained by considering manufacture, assembly, backlash, and sensing errors, for a desired final position and orientation of the mobile platform with the required accuracy to perform attitude manoeuvres for optical aerospace linkage. This work is organised as follows: In Section2 we present both the kinematics inverse problem and the development of the error model, where we provide the global error transmission matrix for each leg of the Gough-Stewart mechanism. In Section 3 we provide a parametrisation of the orientation space, and we show a graphical solution of the model herein proposed in order to assess the finest error in actuator’s lengths, for given tolerances in manufacturing, assembling, backlash, and sensing, required to satisfy typical accuracy for aerospace optical linkage. Finally, in Section 4 we provide some interesting final remarks of the present work.


## 2. Mathematical Modelling of Gough-Stewart Platform

Parallel manipulators, as Gough-Stewart platform, have many advantages; for example, the payload is shared between the actuators, actuators are located in the base so inertia is reduced, and parallel manipulators exhibit stiff behaviour. On the other hand, they possess some disadvantages as their small workspace, limited dexterity [30], in some cases more energy consumption [31], and a greater sensitivity to input errors [32] as compared to the serial manipulators.A problem inherent in Gough-Stewart platform design is that the performance is highly dependent on the geometry; furthermore, different performance measures (workspace, dexterity, etc.) are mutually dependent [33], which complicates the kinematics of the mechanism, exhibiting a highly nonlinear behaviour.The Gough-Stewart platform consists of 8 main elements that form a hexagonal octahedron, a fixed platform, a mobile platform, and six linear actuators. As shown in Figure2, the union of the actuators with the platforms is through articulated joints with certain number of degrees of freedom to the fixed platform. Likewise, the other end of each of actuator is connected by another type of gasket to the mobile platform.Figure 2
6-6 Gough-Stewart platform.The kinematic analysis of the Gough-Stewart platform can be done depending on the input of known variables. Forward modelling consists in the determination of orientation and position of mobile platform, for given lengths of the actuators. The inverse kinematics consists in determining the required lengths on each leg, in order to achieve a certain orientation and position of mobile platform [34, 35]. As in this work we are interested in providing a given accuracy in attitude manoeuvres for aerospace optical link, in what follows we refer only to inverse kinematics of Gough-Stewart platform.


### 2.1. Inverse Modelling

For Gough-Stewart parallel mechanism, the inverse kinematics consists in determining the length of each leg, given the position and orientation of the mobile platform. It is of great importance when one needs the mobile platform to control its position and orientation, particularly when one needs it to follow a determined trajectory with certain accuracy. In what follows, we shall follow the notation of Gosselin [36].If the fixed base radium isRB and the angular separation between connection articulations in the base platform is ϕB, then b→i is the position vector of the articulation of the ith leg from a coordinate system whose origin is attached to the centre of the fixed base (see Figure 3) and is given by(2)b→i=RBcos⁡θiRBsin⁡θi0,i=1,…,6,where θi is the angular position of the ith articulation leg in the fixed base and is given by(3)θ1θ2θ3θ4θ5θ6=ϕB2π3-ϕB2π3+ϕB4π3-ϕB4π3+ϕB-ϕB.Figure 3
Considered 6-6 Gough-Stewart platform. Base platform in blue, mobile platform in magenta, and actuators in green. Base attached system of reference in red and mobile attached system of reference in cyan.In the same fashion, if the mobile base radium isRp and the angular separation between connection articulations in the mobile platform is ϕp, then p→i′ is the position vector of the articulation of the ith leg from a primed coordinate system whose origin is attached to the centre of the mobile base (see Figure 3) and is given by(4)p→i′=Rpcos⁡ηiRpsin⁡ηi0,i=1,…,6,where ηi is the angular position of the ith articulation leg in the mobile base (primed coordinate system) and is given by(5)η1η2η3η4η5η6=ϕP2π3-ϕP2π3+ϕP4π3-ϕP4π3+ϕP-ϕP.If the position vector of the mobile platform’s centre is given by(6)r→=xyz,measured from the nonprimed coordinate system whose origin lies in the centre of the fixed platform, and if the mobile platform orientation is determined by the rotation matrix(7)Q=q11q12q13q21q22q23q31q32q33,then the position vector of the articulation in the mobile platform of the ith leg measured from the nonprimed coordinate system is given by(8)p→i=r→+Qp→i′,i=1,…,6.Subtractingb→i from both sides of the later equation yields (9)p→i-b→i=r→+Qp→i′-b→i.The vectorp→i-b→i points from the fixed to the mobile articulation of the ith leg, so its length ci is given by its Euclid’s norm:(10)ci=p→i-b→i=r→+Qp→i′-b→i,i=1,…,6.The general solution to the inverse kinematics problem can be rewritten as(11)ci=Ui2+Vi2+Wi2,i=1,…,6,where(12)Ui=x+q11Rpcos⁡ηi+q12Rpsin⁡ηi-RBcos⁡θi,(13)Vi=y+q21Rpcos⁡ηi+q22Rpsin⁡ηi-RBsin⁡θiand(14)Wi=z+q31Rpcos⁡ηi+q32Rpsin⁡ηi.Therefore, for a given (or desired) position of the centre of the mobile platformr→ and orientation Q, the required length of each leg is governed by (11).


### 2.2. Leg Length Error Modelling

The error in the length of theith leg, in virtue of (11), is given by(15)dci=1ciUidUi+VidVi+WidWi,where the error in functions Ui, Vi, and Wi is given by the following vector relation: (16)dF→iT=MidE→iT,where dF→i=dUi,dVi,dWi, Mi is the error transmission matrix for the ith leg, and dE→i is a vector that takes into account the global error on each element that determines the length of ith leg. There are thirteen elements per leg that may affect the accuracy of Gough-Stewart platform, so(17)dE→i=dx,dy,dz,dRB,dRP,dθi,dηi,dq11,dq12,dq21,dq22,dq31,dq32.d x, dy, and dz measure the error in position while dqij∀i,j=1,…,3, measure the error in attitude determination. The errors in the radii of the base and the platform (dRB and dRP, respectively), as well as in their angular positions (dθi and dηi∀i=1,…,6), are the addition of those errors due to manufacture, assembly, backlash, and sensing.The error transmission matrix for the length of theith leg is therefore given by(18)MiT=100010001-cos⁡θi-sin⁡θi0q11cos⁡ηi+q12sin⁡ηiq21cos⁡ηi+q22sin⁡ηiq31cos⁡ηi+q32sin⁡ηiRBsin⁡θiRBcos⁡θi0RPq12cos⁡ηi-q11sin⁡ηiRPq22cos⁡ηi-q21sin⁡ηiRPq32cos⁡ηi-q31sin⁡ηiRPcos⁡ηi00RPsin⁡ηi000RPcos⁡ηi00RPsin⁡ηi000RPcos⁡ηi00RPsin⁡ηi.


## 2.1. Inverse Modelling

For Gough-Stewart parallel mechanism, the inverse kinematics consists in determining the length of each leg, given the position and orientation of the mobile platform. It is of great importance when one needs the mobile platform to control its position and orientation, particularly when one needs it to follow a determined trajectory with certain accuracy. In what follows, we shall follow the notation of Gosselin [36].If the fixed base radium isRB and the angular separation between connection articulations in the base platform is ϕB, then b→i is the position vector of the articulation of the ith leg from a coordinate system whose origin is attached to the centre of the fixed base (see Figure 3) and is given by(2)b→i=RBcos⁡θiRBsin⁡θi0,i=1,…,6,where θi is the angular position of the ith articulation leg in the fixed base and is given by(3)θ1θ2θ3θ4θ5θ6=ϕB2π3-ϕB2π3+ϕB4π3-ϕB4π3+ϕB-ϕB.Figure 3
Considered 6-6 Gough-Stewart platform. Base platform in blue, mobile platform in magenta, and actuators in green. Base attached system of reference in red and mobile attached system of reference in cyan.In the same fashion, if the mobile base radium isRp and the angular separation between connection articulations in the mobile platform is ϕp, then p→i′ is the position vector of the articulation of the ith leg from a primed coordinate system whose origin is attached to the centre of the mobile base (see Figure 3) and is given by(4)p→i′=Rpcos⁡ηiRpsin⁡ηi0,i=1,…,6,where ηi is the angular position of the ith articulation leg in the mobile base (primed coordinate system) and is given by(5)η1η2η3η4η5η6=ϕP2π3-ϕP2π3+ϕP4π3-ϕP4π3+ϕP-ϕP.If the position vector of the mobile platform’s centre is given by(6)r→=xyz,measured from the nonprimed coordinate system whose origin lies in the centre of the fixed platform, and if the mobile platform orientation is determined by the rotation matrix(7)Q=q11q12q13q21q22q23q31q32q33,then the position vector of the articulation in the mobile platform of the ith leg measured from the nonprimed coordinate system is given by(8)p→i=r→+Qp→i′,i=1,…,6.Subtractingb→i from both sides of the later equation yields (9)p→i-b→i=r→+Qp→i′-b→i.The vectorp→i-b→i points from the fixed to the mobile articulation of the ith leg, so its length ci is given by its Euclid’s norm:(10)ci=p→i-b→i=r→+Qp→i′-b→i,i=1,…,6.The general solution to the inverse kinematics problem can be rewritten as(11)ci=Ui2+Vi2+Wi2,i=1,…,6,where(12)Ui=x+q11Rpcos⁡ηi+q12Rpsin⁡ηi-RBcos⁡θi,(13)Vi=y+q21Rpcos⁡ηi+q22Rpsin⁡ηi-RBsin⁡θiand(14)Wi=z+q31Rpcos⁡ηi+q32Rpsin⁡ηi.Therefore, for a given (or desired) position of the centre of the mobile platformr→ and orientation Q, the required length of each leg is governed by (11).


## 2.2. Leg Length Error Modelling

The error in the length of theith leg, in virtue of (11), is given by(15)dci=1ciUidUi+VidVi+WidWi,where the error in functions Ui, Vi, and Wi is given by the following vector relation: (16)dF→iT=MidE→iT,where dF→i=dUi,dVi,dWi, Mi is the error transmission matrix for the ith leg, and dE→i is a vector that takes into account the global error on each element that determines the length of ith leg. There are thirteen elements per leg that may affect the accuracy of Gough-Stewart platform, so(17)dE→i=dx,dy,dz,dRB,dRP,dθi,dηi,dq11,dq12,dq21,dq22,dq31,dq32.d x, dy, and dz measure the error in position while dqij∀i,j=1,…,3, measure the error in attitude determination. The errors in the radii of the base and the platform (dRB and dRP, respectively), as well as in their angular positions (dθi and dηi∀i=1,…,6), are the addition of those errors due to manufacture, assembly, backlash, and sensing.The error transmission matrix for the length of theith leg is therefore given by(18)MiT=100010001-cos⁡θi-sin⁡θi0q11cos⁡ηi+q12sin⁡ηiq21cos⁡ηi+q22sin⁡ηiq31cos⁡ηi+q32sin⁡ηiRBsin⁡θiRBcos⁡θi0RPq12cos⁡ηi-q11sin⁡ηiRPq22cos⁡ηi-q21sin⁡ηiRPq32cos⁡ηi-q31sin⁡ηiRPcos⁡ηi00RPsin⁡ηi000RPcos⁡ηi00RPsin⁡ηi000RPcos⁡ηi00RPsin⁡ηi.


## 3. Tolerance Calculations for Given Expected Position and Orientation Errors for Aerospace Optical Link

For an aerospace optical link, it is necessary to consider the communications laser properties, such as power, wavelength, and divergence angle. In order to achieve a correct position and orientation to establish an accurate link between an optical ground station and a satellite orbiting at LEO orbit (for example, 350 km of altitude), the laser divergence angle is a critical parameter. The latter actually intervenes in the computation of the required tolerances for the linear actuators of a Gough-Stewart platform. Figure4 depicts the divergence angle of a laser.Figure 4
Laser beam divergence angle.In order to assess the required tolerances of the platform for the specific requirement of orientation precision for PAT at a LEO orbit ofh=350 km of altitude and considering that the movement of any pointing mechanism is best described by using spherical coordinates, we ought to parametrise rotation matrix (7) in terms of two spherical angles, named α and β, as follows: (19)Q=Rotz;αRotx;β=cos⁡α-cos⁡βsin⁡αsin⁡αsin⁡βsin⁡αcos⁡αcos⁡β-cos⁡αsin⁡β0sin⁡βcos⁡β,where α∈[0,2π) and β∈[0,π/2], which correspond to the z and x related Euler angles. Thus, the subtended area dA by a solid angle dΩ from PAT to the CubeSat, since the altitude remains unchanged along all the pointing process, is given by (see Figure 5)(20)dA=h2dΩ=h2sin⁡βdβdα.Figure 5
Solid angle subtended from PAT to CubeSat.Considering that a satellite has a constant velocity, the accuracy in bothα and β is the same; i.e., dβ=dα; it is given by(21)dββ=dAh2sin⁡β.The areadA ought to be estimated as follows: considering that the high power communications laser ray has a divergence angle of γ=4 mrad, then(22)dA=πh2tan2⁡γ2is the required area to be matched by the PAT system. Thus, by inserting (22) in (20), the equation for accuracy in α is finally given by(23)dββ=πtan2⁡γ/2sin⁡β=3.545mrad.


### 3.1. Assessment of Required Leg’s Tolerances in Positioning

The assessment of the required leg’s tolerancesdci for a specified accuracy in position and orientation for a typical PAT operation, for a given Gough-Stewart platform, can be done as shown in Figure 6.Figure 6
Conceptual design of Gough-Stewart platform’s payload for optical ground station for LEO satellite link, which consists essentially of a communications laser for the uplink and a telescope for downlink.As the optical ground station (OGS) for LEO satellite link (see Figure6) consists essentially of a telescope to establish downlink and a laser device for uplink, these elements constitute the payload of Gough-Stewart platform in this work. Considering a 33 kg, 10-inch diameter telescope with a base of 25 inches built into the Gough-Stewart platform with a long exposure adjustment mirror and a fibber to free-space communications laser of 4.23×1.37 inches embedded to the telescope, the required mobile platform is of RP=1 m and ϕP=π/12 rad. Thus, a base platform of RB=1.5 m and ϕB=π/6 rad is optimal to perform the required PAT operations.This design was achieved by deciding that the (maximal) radius of moving platform (RP) in order to mount the almost 65 cm base of the telescope should provide a margin of at least 10 cm from the base of the telescope to the (minimal radius) edge of the moving platform. This margin is enough for fixing elements of the base of the telescope (without compromising the mechanical integrity of the platform), as well as for placing electronic elements as cables and circuitry for the operation of payload. Then, with considerations of distribution of the weight of payload between the actuators and the constrained workspace (there is a maximum azimuth angle established by ITU for optic links [37], which constrains the workspace of the parallel manipulator), the fixed platform was sized (RB). The angular design of joints placement (ϕP=π/12 rad and ϕB=π/6) responds to the distribution of the weight load between actuators. This joint distribution also guarantees the mobility of passive universal joints in both fixed and moving platforms in the constrained workspace.Once the platform is specified, errors in the radii of the base and the moving platform, as well as in the angular positions of passive universal joints, due to manufacture, assembly, backlash, and sensing, can be observed in Table1. For this case, we shall consider a position of the mobile platform of (0,0,45) cm, according to (6).Table 1
Total errors of the platform due to manufacture, assembly, backlash, and sensing for optical LEO link application.
Error Value d θ 1 (mrad) 1.745 d θ 2 (mrad) 1.747 d θ 3 (mrad) 1.744 d θ 4 (mrad) 1.743 d θ 5 (mrad) 1.746 d θ 6 (mrad) 1.745 d η 1 (mrad) 1.746 d η 2 (mrad) 1.748 d η 3 (mrad) 1.743 d η 4 (mrad) 1.744 d η 5 (mrad) 1.747 d η 6 (mrad) 1.745 d R B (mm) 0.15 d R P (mm) 0.1 d x (cm) 0.1 d y (cm) 0.1 d z (cm) 0.1Considering accuracy in orientation given by (23) and errors in Table 1, (15) for i=1,…,6 turns into scalar functions of the orientation angles; i.e., dci=dci(α,β). Thus, it is possible to observe (see Figure 7) in the region of the parameter space (0≤α≤2π and 0≤β≤2π/9  (this maximum azimuth angle is the one established by ITU for optic links [37], which constrains the workspace of the parallel manipulator)) the points that require the finest precision in the length of the legs in order to provide the required precision in positioning and orientation. As for the lack of (probably expected) symmetry in the graphs of Figure 7, the scale of the vertical axes of such graphics ought to be taken into account. The depicted changes in dci∀i=1,…,6 are very small across the parallel robot workspace, so although important and not negligible they are due to the error configuration of current platform (see Table 1). Furthermore, the distribution of weight along the moving platform is not completely equal; for instance, as mentioned previously, besides the telescope there are other elements as a high power laser device.Figure 7
Graphical assessment of the required length error of each leg of the Gough-Stewart platform as function of the orientation anglesα and β. Functions plotted taking the modules of each term in (15).
(a)
d c 1 ( α , β ) (b)
d c 2 ( α , β ) (c)
d c 3 ( α , β ) (d)
d c 4 ( α , β ) (e)
d c 5 ( α , β ) (f)
d c 6 ( α , β )Analytically, the precision of each leg can be determined as follows. First of all, critical points ought to be identified for each leg. Asdci(α,β):R2→R is a scalar function of two variables, the critical points of the ith leg are those in which the first derivative of the function (its gradient) is zero [38]. In other words, these are the set of points Xij=(αij,βij) such that(24)∇dciXij=0.Then, in order to identify which of those points{Xi,j} are the sought global minima, the evaluation of the second derivatives of dci is required. For a scalar multivariable function, this requires that(25)HdciXij>0and(26)∂dci∂α2>0or∂dci∂β2>0where H is the determinant of the Hessian matrix of the system [38], formed by all the possible second partial derivatives of the function. For the case of the platform considered in this work, the results for the global minima and their values for each leg are shown in Table 2. From these results, we can observe that the most accurate actuator ought to be leg 5, which ought to possess an error in the length of dc5=±0.186 mm.Table 2
Minimum value for the tolerances of each leg,dci, in the orientation domain α∈[0,2π) and β∈[0,2π/9].
Leg Minimum tolerance value Orientation d c i (cm) α (rad) β (rad) 1 0.1416 0.0000 0.0000 2 0.1977 0.3829 0.0000 3 0.1217 0.3185 0.6979 4 0.1794 0.0000 0.6898 5 0.0186 0.2940 0.4846 6 0.1982 0.0531 0.6981We ought to remark that the obtained results for each leg (see Table2) correspond to the orientation which requires most accuracy for it. Also that the reported accuracy for each actuator is different due to the differences in specific errors (Table 1) as well as difference in the distribution of load across the moving platform. On the other side, as we set the accuracy of all actuators as the minimum (of those reported in Table 2), this guarantees that when the ith actuator requires its maximum accuracy, the others are working far from their tolerance. Thus, for instance, if one is to consider actuators with an error of dci=±0.15 mm or dci=±0.10 mm, the Gough-Stewart platform would be able to accomplish its pointing within the required orientation error of dα=dβ=3.5449124284mrad (see (23)), with security factors of 19.35% and 46.24%, respectively.


## 3.1. Assessment of Required Leg’s Tolerances in Positioning

The assessment of the required leg’s tolerancesdci for a specified accuracy in position and orientation for a typical PAT operation, for a given Gough-Stewart platform, can be done as shown in Figure 6.Figure 6
Conceptual design of Gough-Stewart platform’s payload for optical ground station for LEO satellite link, which consists essentially of a communications laser for the uplink and a telescope for downlink.As the optical ground station (OGS) for LEO satellite link (see Figure6) consists essentially of a telescope to establish downlink and a laser device for uplink, these elements constitute the payload of Gough-Stewart platform in this work. Considering a 33 kg, 10-inch diameter telescope with a base of 25 inches built into the Gough-Stewart platform with a long exposure adjustment mirror and a fibber to free-space communications laser of 4.23×1.37 inches embedded to the telescope, the required mobile platform is of RP=1 m and ϕP=π/12 rad. Thus, a base platform of RB=1.5 m and ϕB=π/6 rad is optimal to perform the required PAT operations.This design was achieved by deciding that the (maximal) radius of moving platform (RP) in order to mount the almost 65 cm base of the telescope should provide a margin of at least 10 cm from the base of the telescope to the (minimal radius) edge of the moving platform. This margin is enough for fixing elements of the base of the telescope (without compromising the mechanical integrity of the platform), as well as for placing electronic elements as cables and circuitry for the operation of payload. Then, with considerations of distribution of the weight of payload between the actuators and the constrained workspace (there is a maximum azimuth angle established by ITU for optic links [37], which constrains the workspace of the parallel manipulator), the fixed platform was sized (RB). The angular design of joints placement (ϕP=π/12 rad and ϕB=π/6) responds to the distribution of the weight load between actuators. This joint distribution also guarantees the mobility of passive universal joints in both fixed and moving platforms in the constrained workspace.Once the platform is specified, errors in the radii of the base and the moving platform, as well as in the angular positions of passive universal joints, due to manufacture, assembly, backlash, and sensing, can be observed in Table1. For this case, we shall consider a position of the mobile platform of (0,0,45) cm, according to (6).Table 1
Total errors of the platform due to manufacture, assembly, backlash, and sensing for optical LEO link application.
Error Value d θ 1 (mrad) 1.745 d θ 2 (mrad) 1.747 d θ 3 (mrad) 1.744 d θ 4 (mrad) 1.743 d θ 5 (mrad) 1.746 d θ 6 (mrad) 1.745 d η 1 (mrad) 1.746 d η 2 (mrad) 1.748 d η 3 (mrad) 1.743 d η 4 (mrad) 1.744 d η 5 (mrad) 1.747 d η 6 (mrad) 1.745 d R B (mm) 0.15 d R P (mm) 0.1 d x (cm) 0.1 d y (cm) 0.1 d z (cm) 0.1Considering accuracy in orientation given by (23) and errors in Table 1, (15) for i=1,…,6 turns into scalar functions of the orientation angles; i.e., dci=dci(α,β). Thus, it is possible to observe (see Figure 7) in the region of the parameter space (0≤α≤2π and 0≤β≤2π/9  (this maximum azimuth angle is the one established by ITU for optic links [37], which constrains the workspace of the parallel manipulator)) the points that require the finest precision in the length of the legs in order to provide the required precision in positioning and orientation. As for the lack of (probably expected) symmetry in the graphs of Figure 7, the scale of the vertical axes of such graphics ought to be taken into account. The depicted changes in dci∀i=1,…,6 are very small across the parallel robot workspace, so although important and not negligible they are due to the error configuration of current platform (see Table 1). Furthermore, the distribution of weight along the moving platform is not completely equal; for instance, as mentioned previously, besides the telescope there are other elements as a high power laser device.Figure 7
Graphical assessment of the required length error of each leg of the Gough-Stewart platform as function of the orientation anglesα and β. Functions plotted taking the modules of each term in (15).
(a)
d c 1 ( α , β ) (b)
d c 2 ( α , β ) (c)
d c 3 ( α , β ) (d)
d c 4 ( α , β ) (e)
d c 5 ( α , β ) (f)
d c 6 ( α , β )Analytically, the precision of each leg can be determined as follows. First of all, critical points ought to be identified for each leg. Asdci(α,β):R2→R is a scalar function of two variables, the critical points of the ith leg are those in which the first derivative of the function (its gradient) is zero [38]. In other words, these are the set of points Xij=(αij,βij) such that(24)∇dciXij=0.Then, in order to identify which of those points{Xi,j} are the sought global minima, the evaluation of the second derivatives of dci is required. For a scalar multivariable function, this requires that(25)HdciXij>0and(26)∂dci∂α2>0or∂dci∂β2>0where H is the determinant of the Hessian matrix of the system [38], formed by all the possible second partial derivatives of the function. For the case of the platform considered in this work, the results for the global minima and their values for each leg are shown in Table 2. From these results, we can observe that the most accurate actuator ought to be leg 5, which ought to possess an error in the length of dc5=±0.186 mm.Table 2
Minimum value for the tolerances of each leg,dci, in the orientation domain α∈[0,2π) and β∈[0,2π/9].
Leg Minimum tolerance value Orientation d c i (cm) α (rad) β (rad) 1 0.1416 0.0000 0.0000 2 0.1977 0.3829 0.0000 3 0.1217 0.3185 0.6979 4 0.1794 0.0000 0.6898 5 0.0186 0.2940 0.4846 6 0.1982 0.0531 0.6981We ought to remark that the obtained results for each leg (see Table2) correspond to the orientation which requires most accuracy for it. Also that the reported accuracy for each actuator is different due to the differences in specific errors (Table 1) as well as difference in the distribution of load across the moving platform. On the other side, as we set the accuracy of all actuators as the minimum (of those reported in Table 2), this guarantees that when the ith actuator requires its maximum accuracy, the others are working far from their tolerance. Thus, for instance, if one is to consider actuators with an error of dci=±0.15 mm or dci=±0.10 mm, the Gough-Stewart platform would be able to accomplish its pointing within the required orientation error of dα=dβ=3.5449124284mrad (see (23)), with security factors of 19.35% and 46.24%, respectively.


## 4. Conclusions

In this paper we have developed a mathematical error assessment model for the inverse kinematics for Gough-Stewart parallel mechanisms, which we apply in a case of aerospace optical linkage. Although the model herein developed was intensely derived assuming a 6-6 Gough-Stewart platform, the methodology herein depicted is easily extended to other configurations. Furthermore, the modelling developed in this work has been cast in such a way that we provide transmission error matrices for each leg length, as they are commonly required in the field of instrumentation.The model, for which we show a graphical and an analytic solutions, has been validated in a scenario with high accuracy requirements in orientation as the establishment of a ground-LEO satellite optical link. Although the presented graphical solution heavily relies on posing the error in each leg’s lengths as functions of the orientation angles, i.e.,dci=dci(α,β)∀i=1,…,6, the analytic solution is valid even in cases of spanned solution spaces. For instance, in a more general case, the position of the mobile platform should not be restricted as in the example performed in this paper, yielding a 5-dimensional space in which global mimima ought to be found; i.e., dci=dci(α,β,x,y,z)∀i=1,…,6. Moreover, more general expressions for rotation matrix (7) than those used to exemplify the usage of the method herein depicted (see (19)) can be considered, farther increasing the parametric space. Nevertheless, although in such a case no graphical solution is achievable, analytic solution can be equally performed. Furthermore, also numerical assessment of first and second derivatives in order to find global minima of dci=dci(α,β,x,y,z)∀i=1,…,6 (see (25) and (26)) can be performed to solve the system. This latter fact shows the generality and validity of the method herein proposed.



---

*Source: 9592623-2018-10-03.xml*