Since radio telescope is the main tool for human being to search the universe secret, the astronomer reached unanimity at the 24th URSI Conference in Kyota, Japan, 1993, and proposed to construct the next generation of the large radio telescope (LT) (Nan & Peng, 2000). From then on, the astronomer of China began the project of Five-hundred meter Aperture Spherical radio Telescope (FAST) (Qiu, 1998, Li, 1998).
It is well known that Arecibo is the breakthrough of radio telescope. Its main mirror, 305m in diameter, is fixed on the karst base, and an elaborately designed feed system illuminates a part of the mirror which forms an effective aperture of the telescope with about 200m. The feed system is movable at a height of about 100 m for tracking the object to be observed. The enormous receiving area of the telescope will enable it to make many important astronomical discoveries inaccessible to lesser instruments, despite its small sky coverage (20 zenith scan angle), due to geometrical configuration, and narrow frequency bandwidth, originated from spherical aberration. An upgrade project has recently been carried out for the Arecibo telescope, in which a heavy and complex hence expensive Gregorian dual-reflector feed system is introduced for correcting the spherical aberration and a broad bandwidth is affected (Duan, 1999).
For the sake of satisfying the requirements of low cost and broad bandwidth, the project group of FAST decided to substitute the fixed spherical reflector with active reflector units. As shown in Fig. 1(a), the reflector consists of almost 2000 elementary reflector units. Fig. 1(b) shows some active reflector units and supporting mechanisms. The reflector unit is small part of spherical surface of regular hexagon and is driven by a supporting mechanism. The part of spherical reflector illuminated by the feed is continuously adjusted to fit a paraboloid of revolution in real-time, synchronous with the motion of the feed while tracking the object to be observed. As it is now free from spherical aberration, a simple, light, hence cheap feed system may be adopted to achieve broad bandwidth and full polarization.
In order to fit a paraboloid of revolution, it is necessary that every reflector units should be driven by a supporting mechanism with two rotational degrees of freedom and one translational degree of freedom (Luo et al., 2000). That means almost 6000 control nodes on the whole active reflector should be managed and controled at the same time. It is supposed to be very difficult, so a sharing strategy is derived to decrease the number of nodes, which requires three adjacent nodes combined together to share one driver. Basically, there are two types of mechanism which can fulfill the required movement for each reflector unit and fit for the sharing strategy, 3-PSS mechanism with constraint leg (Wang et al., 2006), shown in Fig. 2(a), 3-PSS+C for abbravation, and 3-PRS mechanism (Tang et al., 2007), shown in Fig. 2(b).
These mechanisms will bring errors because of the control or dimensional factor. Moreover, the fitting surface of reflector will not match exactly with the nominal paraboloid, and the sharing strategy also brings accuracy problem. In order to guarantee the highest working frequency of large spherical radio telescope, 5GHz, the fitting accuracy of active reflector should be studied systematically. Based on the kinematics of 3-PSS+C mechanism, in this chapter, one-dimensional and two-dimensional fitting accuracy on the whole active reflector is analyzed considering control errors. However, about 2000 constraint legs increase almost one quarter of the cost. Thus 3-PRS mechanism is proposed and used as supporting manipulator for reflector unit. Since 3-PRS mechanism has many problems such as parasitic motion, advanced research on kinematics with errors is necessary. Then three-dimensional fitting accuracy is analyzed based on error kinematics of 3-PRS mechanism.
2. The analysis of 3-PSS+C supporting mechanism
2.1. Supporting mechanism description
As shown in Fig. 2(a), the parallel supporting mechanism consists of a base plate, a movable platform, and four connecting legs, three of which have identical kinematic chains, PSS. Each of the three legs is composed of one fixed length link (3), and one union driven plate (5). The fixed length link (3) is connected to the movable platform (1) and the union driven plate (5) by two spherical joints (2) and (4), respectively. The union driven plate (5) is connected to the base plate (7) by a prismatic joint (6). The base plate and the movable platform are two regular triangles. The passive leg (8) connects the center points of the two regular triangles. One end of the passive leg has a 2-DOF universal joint (9), another end is fixed to the base plate (7) by a prismatic joint (10). The passive leg (8) can be extensible with the prismatic joint (10) along its axis line. Furthermore, when the supporting mechanism is assembled, the axis line of the prismatic joint (10) should pass the center of the spherical reflector. Since a supporting mechanism should be driven by three actuator legs, as shown in Fig. 2, the union driven plate (5) connects three fixed length links in order to reduce the actuator number. As a result, the number of actuators of the active reflector is equal to that of the reflector units.
From above description, one can see that the proposed mechanism is such a mechanism with n DOFs, which usually consists of n identical actuated legs with 6 DOFs and one passive leg with n DOFs connecting the movable platform and the base plate, i.e., the DOF of the mechanism is dependent on the passive leg’s DOF. For the mechanism considered in this paper, the passive leg is with three DOFs, which means that n equals to 3. The three DOFs are one translation along axis and two rotations about and axes.
2.2. Kinematics analysis
The mechanism kinematics deals with the study of the mechanism motion as constrained by the geometry of the links. Typically, the study of mechanism kinematics is divided into two parts, inverse kinematics and forward (or direct) kinematics (Wang & Tang, 2003). The inverse kinematics problem involves mapping a known pose (position and orientation) of the output platform of the mechanism to a set of input joint variables that will achieve that pose. The forward kinematics problem involves the mapping from a known set of input joint variables to a pose of the movable platform that results from those given inputs (Wang et al., 2001). Generally, as the number of closed kinematics loops in the parallel mechanism increases, the difficulty of solving the forward kinematics relationships increases while the difficulty of solving the inverse kinematics relationships decreases (Liu et al., 2001).
2.2.1. Inverse kinematics
A kinematics model of the mechanism is developed as shown in Fig. 3. The vertices of the movable platform are denoted as platform joints ), and the vertices of the base plate are denoted as ( A fixed global reference system : is located at the center of the regular triangles with the axis normal to the base plate and the axis parallel to the side . The circumcircle radius of triangles is denoted as . Another reference frame, called the top frame : , is located at the center of regular triangles . The axis is perpendicular to the movable platform and axis parallel to the side . The circumcircle radius of triangles is denoted as . Vector of fixed length links are denoted as ( ), and the link length for each legs is denoted as , where
The objective of the inverse kinematics solution is to define a mapping from the pose of the output platform in the Cartesian space to the set of actuated inputs that achieve that pose. For this analysis, the pose of the movable platform is considered known, and the position is given by the position vector and the orientation is given by a matrix . Then there are
where stands for cosine function, stands for sine function, and and are the orientational DOFs of the movable platform with respect to and axes, respectively. The coordinate of point in the frame can be described by the vector
Vectors will be defined as the position vectors of base joints in frame and
Vectors in frame can be, therefore, written as
Then the inverse kinematics of the parallel mechanism can be solved by writing following constraint equation
Hence, for a given mechanism and prescribed position and orientation of the movable platform, the required actuator inputs can be directly computed from Eq. (6), that is
2.2.2. Forward kinematics
The objective of the forward kinematics solution is to define a mapping from the known set of the actuated inputs to the unknown pose of the output platform. For the architecture with prismatic actuators, the inputs that are considered known are the lengths of the three actuator legs and . The unknown pose of the output platform is described by the position vector and angles and . Because it is very difficult to describe the direct kinematics in closed form for this type of parallel mechanism, the forward kinematics solution should be obtained by numerical methodology as following:
Decide the non-singularity workspace of the mechanism;
Give the initial value of direct kinematics solution;
Calculate the position coordinates of spherical joints, construct the nonlinear equations set by the geometry constraint relationship of fixed length links;
Solve the nonlinear equations set by Quasi-Newton method (Press et al., 1995).
From the Eq. (6), the nonlinear equations are
where the direct kinematics solutions are and
2.2.3. Velocity equation
Eq. (6) can be differentiated with respect to time to obtain the velocity equation. This leads to an equation of the form.
where is the vector of Cartesian velocities defined as
Matrices and are the forward and inverse Jacobian matrices of the mechanism and can be expressed as
where is the unit vector of , and is the element of vector with respect to axis coordinate, and are the elements of vector with respect to and axis coordinates.
2.3. Mechanism accuracy analysis
When the large spherical radio telescope works, the feed system will illuminate a working area, which is the paraboloid reflector with a three-hundred-meter aperture. The part of spherical reflector illuminated by the feed is continuously adjusted to fit a paraboloid of revolution in real-time, synchronous with the motion of the feed while tracking the object to be observed. For the fitting, the spherical surface reflector is divided into some small elementary units. When the mechanisms drive the reflector units to fit the paraboloid, the fitting surface of reflector will not match exactly with the nominal paraboloid. Moreover, the mechanism has error because of the control or dimensional factor. In this section, the mechanism accuracy is analyzed firstly.
The mechanism accuracy involves the error caused by the actuator input error and the joint error of the mechanism. The actuator input error is denoted as and the joint error is denoted as , where includes the joint error on the base platform and the input error . The output error is denoted as
From Eq. (5) and (6), the inverse kinematics equation can be written as
Differentiating Eq. (14) leads to
3. Fitting accuracy analysis of Active reflector
3.1. One dimensional fitting accuracy analysis
As shown in Fig. 4, the base active reflector of the radio telescope is a spherical surface with five-hundred-meter aperture, and the working reflector is a paraboloid with a three-hundred-meter aperture. When it works, the reflector units are driven by the parallel mechanism from the initial position to the fitting position to fit the paraboloid. Because the paraboloid is formed by the revolution of parabola, we can analyze the deviation about spherical surface and paraboloid in the reflector frame : which is built as shown in Fig. 4, where the spherical surface and the paraboloid in the frame are circular arc and parabola, respectively.
Fig. 5 shows one reflector unit which is in the initial position and fitting position, respectively. The initial position is located at the base spherical reflector surface. The deviation from the circular arc to the parabola is denoted as and symbol represents the i-th reflector unit which corresponds to the i-th mechanism. The symbol represents the supporting point of the movable platform. The explanations of other symbols used in accuracy analysis are:The supporting point while the reflector unit is in the initial position. The supporting point while the reflector unit is in the fitting position. The intersecting points of line and the parabola. The refence center in the movable platform while the reflector unit is in the initial position. The refence center in the movable platform while the reflector unit is in the fitting position. The center of spherical reflector. The radius of spherical reflector. The focal point of the paraboloid.
The absolute actuator input of the i-th mechanism is specified as , while the i-th active reflector unit is driven to fit the paraboloid. Obviously, the driven reflector unit will not match exactly with the nominal paraboloid. In order to evaluate the fitting error, as shown in Fig. 5, is defined as the center points deviation of the i-th reflector unit to the corresponding paraboloid and is equal to where the center points deviation is called as one-dimensional fitting error.
3.1.1. Parabola equation and circle equation
According to the reference (Qiu 1998), the focal length of the parabola is specified as then the parabola equation can be written as
The base spherical surface in reflector coordinate system is a circle. And the circle equation can be written as
3.1.2. Actuator input range
The coordinate of the point in the frame can be described by the vector
The equation of straight line can be written as
According to Eqs. (19) and (21), the intersecting point between line and the circle can be expressed by vector which is
Actuator input value of the i-th reflector unit can be written as
3.1.3. One-dimensional fitting error
When actuator input is specified, the fitting error can be reached. The first step is to calculate the position coordinate in the frame by the forward kinematics solution. The position vector of center point in the frame is written as
where is the rotation matrix about frame : to the frame :
where Then the fitting error is expressed as
Since the three-hundred-meter aperture paraboloid is composed of a lot of reflector units, we should analysis all the error of reflector units. When the error is studied in the reflector frame : and the side length of reflector unit is specified, the one-dimensional root-mean-square (RMS) fitting error of the paraboloid reflector with three-hundred-meter aperture is defined as
3.1.4. One-dimensional accuracy synthesis analysis
The accuracy synthesis analysis is defined as the composition RMS error that caused by the mechanism actuator input error and the fitting error. For the mechanism actuator input error has linear relationship with the value of Eq. (23) can be rewritten as
where is the actuator input error. Then Eqs. (24)- (27) can be used to calculate the one-dimensional composition RMS error
3.1.5. Simulation example
Since the position of supporting point should be limited in the range of reflector unit, as shown in Fig. 5(b), the base plate parameter of the parallel mechanism is written as
where is the side length of reflector unit and is the distance from the movable platform edge to reflector unit edge. In this work, m and
Since the paraboloid reflector with the 300m aperture is symmetry, the error can be analyzed in the range of 150m. Fig. 6(a) shows the one-dimensional fitting error when the side length of reflector unit is specified.
According to Eqs. (26) and (28), the one-dimensional RMS fitting error and composition RMS error can be drawn as shown in Fig. 6(b) for mm. In this work, we assume that the maximal input error of the mechanism is 1 mm. When the side length of reflector unit is equal to 7.5m, the one-dimensional RMS fitting error is 3.75 mm.
3.2. Two-dimensional fitting accuracy analysis
As shown in Fig. 5(a), the reflector unit fitting error is a closed region. And section 3.1 only considered the one-dimensional error. In this section, the area of the closed region will be used to analyze and evaluate the fitting error, which is called as two-dimensional fitting error. Obviously, the two-dimensional fitting error will provide more reliable index for us to analyzing the fitting accuracy of the large spherical radio telescope. Fig. 7 shows the two-dimensional fitting error of the i-th reflector unit, which is the sectional region.
3.2.1. Fitting position circle equation
When the reflector units are driven by mechanism, the circle arc equation of the i-th reflector unit will be changed in the frame . As shown in Fig. 7, the centre of circle arc is changed from to . The coordinate of in the frame is written as
The coordinate of in the frame is written as
The circle arc equation is changed to
3.2.2. Two-dimensional fitting error
According to the circle arc equation in the frame , the two-dimensional fitting error can be calculated. Firstly, as shown in Fig. 7, one point in parabola is denoted as , The equation of straight line can be written as
According to Eqs. (33) and (34), the intersecting point between line and the circle arc can be expressed by . The fitting error of given point is expressed as
The area of the closed region can be written as
which is the two-dimensional fitting error of the i-th reflector unit. Then the average error of the two-dimensional fitting error is defined as
In the end, the two-dimensional root-mean-square (RMS) fitting error of the paraboloid reflector with three-hundred-meter aperture is defined as
where is the number of reflector units that consist of three-hundred meter aperture parabola.
3.2.3. Two-dimensional accuracy synthesis analysis
The two-dimensional accuracy synthesis analysis is defined as the composition RMS error that caused by the mechanism actuator input error and the two-dimensional fitting error, which is denoted as . The Eqs. (28), (37) and (38) can be used to calculate the two-dimensional composition RMS fitting error
3.2.4. Simulation example
The two-dimensional RMS error and composition RMS fitting error are shown in Fig. 8, where all the dimensional design parameters are the same as the specified parameters in section 3.1.5. Comparing the Fig. 8 and Fig. 7, we can know although both the one-dimensional and two-dimensional RMS fitting error increase while the side length of reflector unit increases, the two-dimensional RMS error is larger than one-dimensional RMS error.
According to the fitting error requirement given by reference (Qiu, 1998), when the highest working frequency of the radio telescope is 5GHz, the reflector RMS fitting error should be less than 3.75 mm. Now, we can decide the dimensional parameters and guarantee the implementation of the working frequency by the one-dimensional or two-dimensional RMS fitting error curves. For example, according to the Fig. 7, if the side length of reflector unit is specified as 7.5m, the specified dimension of reflector units can satisfy the requirement of 5GHz work frequency. However, as shown in Fig. 8, if the two-dimensional RMS fitting error is used to evaluate the fitting accuracy, the side length of reflector unit should less than 7.0m for satisfying the requirement of 5GHz work frequency.
4. The error kinematics of 3-PRS mechanism
3-PRS mechanism has less chians which reduces its cost. The kinematics of 3-PRS mechanism has been fully analyzed (Carretero et al., 1997, Tsai & Shiau, 2003). Yet, when the mechanical manufacturing and assembling errors are brought into the model, kinematic analysis will become complicated. Therefore, analysis on parasitic motion and accuracy should be made to guarantee the application of 3-PRS mechanism as reflector unit supporting mechanism.
4.1. Kinematic modeling with errors
The magnitude of the reflector driving machine is always at meter, so input error, length error of the legs and location error of the spherical joint have little influence on motion error of the moving platform. On the other hand, the location and angle error of the rotational joint, which will be extended by the legs, will mix with parasitic motion so as to greatly affect the motion. Therefore, we introduces angle error of the rotational axis and location error of the joint point in the rotational joint as the main error resources in order to analyze kinematics of 3-PRS error model.
In the error model representation of 3-PRS mechanism, as shown in Fig. 9(a), is ideal axis vector of the rotational joint, and is actual axis vector with angle error. Similarly, is ideal vector of the rotational joint point, whereas is actual joint point vector with location error. Both and include three direction errors separately along and axis, which means that there are six errors in each leg, as shown in Fig. 9(b), in which two-dot chain line represents ideal rotational joint and real line represents actual one.
The location error vector of the rotational joint is defined as
The angle error vector of the rotational joint is defined as
Then we can find
The three components of the vector are independent while those of the vector are not since the error on the direction can be given through two parameters only. So the relationship between the components of the vector can be determined by , where are unit direction vectors. Thus, error resources are appropriately introduced and error modeling of 3-PRS mechanism is completed.
4.2. Inverse kinematics
The coordinate axes of the inertial frame fixed on the base platform are denoted by : while those of the moving frame fixed on the moving platform are denoted by : (see Fig. 9). In order to simplify the kinematic model, the origin of the inertial frame is located on the center of the base platform and axis of the inertial frame points to one of the spherical joint on the base platform. The axis is also on the plane of the base platform while axis points upward forming a right-handed orthogonal frame. The coordinate axes of the moving frame are also located on the moving platform in the same way. The rotation matrix from the coordinate axes of the moving platform to those of the base platform is denoted by which is expressed as
where , and are variables which orderly specify the rotations around the , and axis, and s represents sin, while c represents cos.is the vector from the origin of the inertial frame to the origin of the moving frame expressed in the inertial frame. is the vector expressed in the inertial frame. is the vector expressed in the inertial frame. is the vector expressed in the inertial frame. is the vector expressed in the inertial frame. is the vector from the origin of the moving frame to the i-th upper attachment point expressed in the moving frame. is the vector from the rotational point to the upper attachment point of the i-th leg. It should be noted that is constant for each leg. is the set of actuated joint variable of the 3-PRS mechanism which is the height of the rotational joint point. We can get
The inverse kinematic problem is supposed to determine the value of the actuated variables for a known position and orientation of the end-effector, that is:
In those six variables, the known numbers are three desired motions which include , and , while the unknown numbers are three parasitic motions which include , and . The parasitic motions are determined by the target motions, that is:
The structure of mechanical joint leads to two geometrical constraints which are rotation constraint and length limitation of the leg.
(a) The rotation constraint
Each attachment point of the moving platform should be restricted in the rotation plane formed by the wheeling leg. The constraint equations are
where is a constant of the rotation plane and determined by the following equation
Substituting Eq. (39) into Eq. (42)
Expressing with the elements of those vectors, we get
Substituting Eq. (40) into Eq. (41)
Expressing with the elements of those vectors, we get
(b) The leg limitation
The distance between the attachment point of the moving platform and the rotational point of the rotational joint should be constant. The constraint equations are
Substitute Eq. (39) into Eq. (47)
Expressing with the elements of those vectors, we get
In Eq. (49), there are two possible solutions for each leg, thereby yielding a total of 8 possible combinations of actuated height for a given position and orientation. In the present work, the negative square root is always selected to yield a solution where the legs are always beneath the moving platform.
Combining Eqs. (44) and (46), we get
Combining Eqs. (49) and (50), we get equations set involved with variables and three parasitic motions. We can find that if , Eqs. (49) and (50) can be solved separately. But in ordinary condition, , direct solution of the equations set will become impossible.
4.3. Arithmetic of inverse solution
In order to figure out the nonlinear equations set with the normal condition of , a numerical iterative arithmetic is proposed as follows:
Decide the initial value of the actuated height and set the loop variable
Calculate the parasitic motions of the moving platform according to Eq. (50) with
Then calculate according to Eq. (49);
If where is the acceptable error limitation, end up iterative calculation with solution of . Otherwise, then turn back to step 2.
4.4. Forward kinematics
Forward kinematic solution is supposed to determine the position and orientation of the end-effector for known actuated variables. Since we already have the inverse kinematics, Similarly with the method used in setion 2.2.2, we can also get forward solution.
5. Three-dimensional fitting accuracy analysis
5.1. Coordinate description for calculation
As shown in Fig. 10(a), the paraboloid of revolution covers a spherical surface with a 300-meter aperture, which has a coning angle of 60 degree. The paraboloid moves on the sphere to track the object in real-time. Since the axial line of the desired paraboloid always orients to the center of the base spherical surface when it works, the fitting accuracy is constant at any time when the paraboloid of revolution is at any location on the base spherical surface. Thus analysis of Fig. 10(a), where the peak of the paraboloid of revolution is located on the bottom of the base spherical surface, will be enough.
In Fig. 10(a), the global coordinate axes of the inertial frame fixed on the whole active reflector system are denoted by : . The origin of the coordinate system is located on the bottom of the base sphere while the plane is the tangent plane of the point
Fig. 10(b) is a top view of base sphere along radial direction. The axis is perpendicular with the side of the regular hexagon reflector unit. Fig. 10(b) shows the location of the three attachment points and twenty four sampling points for calculation. It is noted that the upper surface of the reflector unit is spherical which means that the sampling points are all on a spherical surface whose radius is
As shown in Fig. 11, the coordinate axes of the inertial frame fixed on the platform in initial position are denoted by : , in which the location of the axes is similar with the location in the kinematic analyses introduced in setion 4.1. The inertial frame is not moving with the reflector unit. Since the reflector unit in initial position is the tangent plane of the base sphere, axis always points to the center of the sphere. It should be noticed that the directions of axis and axis in the coordinate frame is different from those in . Similarly, the coordinate axes of the moving frame fixed on the platform in fitting position and orientation are denoted by : . The coordinate frame and are similar with the kinematic inertial and moving frame analyzed in section 4 so that the inverse and forward solution can be used.
We can divide the surface by several circular arcs. These arcs are intersections of the base sphere and the planes which are parallel with plane, and are all through the centers of the reflector units. So the fitting accuracy of all reflector units can be calculated along these arcs.is the i-th attachment on the reflector unit. is the center of the base sphere. is the radius of the base sphere. is the focal point of the ideal paraboloid of revolution. is the side length of the reflector unit. is the distance from each attachment to the border of the reflector unit on radial direction. is the vector from the origin of to on the reflector unit in the initial position expressed in the inertial frame (the subscript mean the j-th reflector unit on the t-th circular arc, the same as below). is the vector from the origin of to in the initial position expressed in the global inertial frame is the vector from to the k-th sampling point in the fitting position expressed in the moving frame , so is constant and known. is the vector from to the k-th sampling point in the fitting position expressed in the inertial frame is the vector from to the k-th sampling point in the fitting position expressed in the inertial frame is the vector from to expressed in the inertial frame which is determine by the position of current analyzing reflector unit. is the vector from to expressed in the inertial frame is the rotation matrix from coordinate frame to is the rotation matrix from coordinate frame to
5.2. Paraboloid equation and circle equation
The focal length of the paraboloid is specified as 0.467K, then the paraboloid equation can be written as
The distance from one point to the paraboloid on radial direction is
By combining the equations above, we define the function to calculate distance between one point and the paraboloid on radial direction as
5.3. Driving strategy
Driving strategy determines the method to drive the reflector unit to fit for paraboloid. In order to simplify calculation, the actuated variable is the value of actuated variable in initial position plus
5.4. Fitting accuracy calculation
Calculating the number of the arcs in the half base sphere of positive axis surface, we get:
where means the least integer which is no less than
The radius of the t-th arc is written as
The meeting point of the t-th arc to the border of the base sphere is written as
Thus the length of the t-th positive half circular arc is written as
So the number of reflector units on the positive half t-th arc can be obtained as
First, we focus on the coordinate frame and . The equation of the base sphere is
So that along the arc of even number, the position of the j-th reflector unit on the t-th arc can be written as
While along the arc of odd number, the is
Thus position of the coordinate frame which is fixed on the reflector unit in the initial position is determined by
According to the orientation of the coordinate frame , the rotation matrix to can be written as
where is rotation matrix of rotating angle of degree around vector , and is rotating vector in plane which is expressed as
The vector of sampling point can be written as
while the vector of the upper attachment in the initial position can be written as
Then, we focus on the coordinate frame and . The value of actuated variable should be determined for the unit. According to the driving strategy, the actuated variable of the reflector unit in fitting position can be obtained as
where is the initial value of the actuated variable and
According to forward solution, the Cartesian variables which specify the position and orientation from the coordinate frame to the inertial frame can be obtained as
The vector of the k-th sampling point expressed in inertial frame can be obtained by
Finally, we calculate the fitting accuracy based on the results obtained above. The fitting accuracy of one sampling point on the j-th reflector unit in the t-th arc can be written as
Synthetically, substituting Eqs. (56)-(64) into (65), we can calculate with specified t, j and k as
The RMS fitting accuracy on the whole fitting surface can be written as
5.5. Driving strategy optimization
In order to reduce the RMS fitting accuracy of the whole reflector, the current driving strategy should be optimized. Considering the real-time compensation, the optimization algorithm should be no more difficult than the current strategy.
According to the analyses above, we put forward one modified driving strategy as
The strategy will compensate the actuated variable of each reflector unit with the algebraic average among three supporting points and the center of the reflector so as to get less RMS fitting accuracy on each unit spherical surface.
5.6. Simulation example
Fig. 12(a) shows the fitting accuracy on the whole reflector range with and without optimized driving strategy when the side length is changeable, which are respectively expressed as and . All the dimensional design parameters are the same as the specified parameters in section 3.1.5.
As shown in Fig. 12(a), the fitting accuracy is reduced approximately by 40% when the driving strategy is optimized. Fig. 12(b) shows the fitting accuracy on the whole reflector range with and without errors when the optimized driving strategy is used, which are respectively expressed as and . In this work, we assume that the axis angle tolerance of the rotational joint is and the position tolerance is , and is the worst situation with these tolerances.In order to guarantee the working frequency of the large radio telescope, the side length of the reflector unit should be less than 7.3mm.
6. Conclusion and future works
In order to guarantee the usage of active reflector and achieve the highest working frequency requirement, 5GHz, in FAST, fitting accuracy of the active reflector is supposed to be analyzed. In this chapter, a novel 3 DOFs parallel mechanism, 3-PSS with constraint leg, is proposed. This mechanism can fulfill the required movement to fit a paraboloid of revolution for the active reflector. The kinematics of 3-PSS+C mechanism is studied. Based on that, the one and two-dimensional fitting accuracy are calculated and the side length limit of reflector units is evaluated as 7.0m. However, due to more expensive cost of the extra constraint chain, 3-PSS+C is not very appropriate as reflector supporting mechanism. So 3-PRS mechanism becomes more attractive and deserves to pay more attention. Then error kinematics with rotational joint tolerance is analyzed for actual application. Based on that, three-dimensional fitting accuracy is calculated with optimized driving strategy, and the side length limit turns out to be 7.3m.
The future work will still focus on the fitting accuracy not only on kinematics, but also on synthetical design, stiffness and control, as well as sharing strategy study. Further more, experiment research will be taken into account as certification for the theoretical analysis.