Real-Time Monitoring of the Position and Orientation of a Radio Telescope Sub-Reflector with Fiber Bragg Grating Sensors

Environmental loads linked with pointing errors, such as gravity, thermal gradients, and wind disturbances, are a serious concern for large-aperture high-frequency radio telescopes. For the purpose of maintaining the pointing performance of a telescope, a contact measurement scheme is proposed on basis of fiber Bragg grating (FBG) strain sensors that can monitor the sub-reflector shift in real time as the input data of the adjustment system. In this scheme, the relationship between the in situ strain measurement and the deformation of the supporting structure, which is the main cause of sub-reflector shift, is deduced using the inverse Finite Element Method (iFEM). Finally, experimental studies are carried out on a simple physical structure model to validate the effectiveness and accuracy of the contact measurement scheme.


Introduction
With the aim of satisfying the pointing performance of a large-aperture high-frequency radio telescope, it is important to maintain the ideal location relationship between the sub-reflector and the main reflector. Nevertheless, environmental loads such as gravity, thermal gradients, and wind disturbances may deform the shape of the supporting structure of the sub-reflector, which will seriously affect the location relationship.
In a series of three papers, the active control system has been widely used for adjusting the position and orientation of the sub-reflector [1][2][3]. In order to provide the exact input for the active control system, the position and orientation of the sub-reflector must be firstly and accurately obtained. For this purpose, a non-contact measurement scheme based on a Position Sensing Device (PSD) consisting of a laser diode and a Complementary Metal Oxide Semiconductor (CMOS) camera was developed by the Sardinia Radio Telescope (SRT) Metrology team [4,5]. In this scheme, the accuracy of the PSD was able to reach around 0.1 mm within the measurement range of ±40 mm. Because the laser-catching capability of the CMOS camera is seriously affected by ambient light, it is only during the night that this non-contact measurement scheme can be successfully performed.
To overcome the aforementioned limitation, a contact measurement technology based on sensing structural deformations from in situ strain measurements, namely shape sensing, has offered a method for structural health monitoring (SHM) that can be applied to active control systems of large structures [6][7][8]. Importantly, compared to traditional strain sensors, fiber Bragg grating (FBG) strain Sensors 2019, 19, 619 2 of 14 sensors have been extensively studied and applied for shape sensing due to their lightness, accuracy, and resistance to electromagnetic interference, radio frequency interference, and radiation [9][10][11].
For a typical radio telescope, the supporting structure comprises a fixed platform and four supporting legs; the sub-reflector is fixed on the fixed platform ( Figure 1a). Neglecting the deformation of the fixed platform (the stiffness of the fixed platform is big enough to resist environmental loads), the sub-reflector shift mainly arises from the fixed platform shift associated with the deformations of the four supporting legs caused by the weight of the sub-reflector. Thus, when the end-node shifts of the inside beams of the four supporting legs are known, the position of the fixed platform can be determined. Also, the orientation of the fixed platform can be computed from the end-node shifts with a direction cosine algorithm. Furthermore, the position and orientation of the sub-reflector can be adjusted using the Stewart platform attached to the fixed platform (Appendix A).
Sensors 2019, 19, x FOR PEER REVIEW 2 of 14 accuracy, and resistance to electromagnetic interference, radio frequency interference, and radiation [9][10][11]. For a typical radio telescope, the supporting structure comprises a fixed platform and four supporting legs; the sub-reflector is fixed on the fixed platform ( Figure 1a). Neglecting the deformation of the fixed platform (the stiffness of the fixed platform is big enough to resist environmental loads), the sub-reflector shift mainly arises from the fixed platform shift associated with the deformations of the four supporting legs caused by the weight of the sub-reflector. Thus, when the end-node shifts of the inside beams of the four supporting legs are known, the position of the fixed platform can be determined. Also, the orientation of the fixed platform can be computed from the end-node shifts with a direction cosine algorithm. Furthermore, the position and orientation of the sub-reflector can be adjusted using the Stewart platform attached to the fixed platform (Appendix A). Therefore, the present paper aims to develop a contact measurement scheme according to shape sensing technology and bare FBG strain sensors (Figure 1b) in order to measure the deformations of the inside beams of the four supporting legs in real time. The key of this scheme is to construct a suitable and robust measurement model to describe the relationship between the beam deformations and the in situ strain measurements. Also, in order to monitor the influence of the temperature variation on the strain measurement system, several FBG temperature sensors (Figure 1c) are set on the surface of the beam used to capture the current temperature ( Figure 1a). When the current temperature is obtained, the wavelength shift of the FBG strain sensor caused by temperature variation can be estimated, and then this wavelength shift can be removed from the strain measurement.
In order to sense the deformation of the beam element, a load-independent method was developed by Ko et al. according to the piece-wise continuous polynomials and the classical beam theory [12,13]. This one-dimensional scheme has been demonstrated to be sufficiently accurate for predicting deformation and somewhat less accurate for evaluating the cross-sectional twist of a beam/frame structure; however, a large number of strain sensors are required to capture in situ strain data. Bogert et al. analyzed and verified the effectiveness of the modal transformation method for reconstructing the structural deformation [14]. Despite the advantages of this method, a great many natural vibration modes are required. Another problem is that this method requires an analysis of Therefore, the present paper aims to develop a contact measurement scheme according to shape sensing technology and bare FBG strain sensors (Figure 1b) in order to measure the deformations of the inside beams of the four supporting legs in real time. The key of this scheme is to construct a suitable and robust measurement model to describe the relationship between the beam deformations and the in situ strain measurements. Also, in order to monitor the influence of the temperature variation on the strain measurement system, several FBG temperature sensors ( Figure 1c) are set on the surface of the beam used to capture the current temperature ( Figure 1a). When the current temperature is obtained, the wavelength shift of the FBG strain sensor caused by temperature variation can be estimated, and then this wavelength shift can be removed from the strain measurement.
In order to sense the deformation of the beam element, a load-independent method was developed by Ko et al. according to the piece-wise continuous polynomials and the classical beam theory [12,13]. This one-dimensional scheme has been demonstrated to be sufficiently accurate for predicting deformation and somewhat less accurate for evaluating the cross-sectional twist of a beam/frame structure; however, a large number of strain sensors are required to capture in situ strain data. Bogert et al. analyzed and verified the effectiveness of the modal transformation method for reconstructing the structural deformation [14]. Despite the advantages of this method, a great many natural vibration modes are required. Another problem is that this method requires an analysis of numerous eigenvalues and a detailed description of the elastic and inertial material properties. Therefore, a suitable algorithm for shape sensing should be [15]: (1) general enough to accommodate complex structural topologies and boundary conditions; (2) robust, stable, and accurate under a wide range of loads, material systems, inertial/damping characteristics, and inherent errors during the strain measurements; (3) sufficiently fast for applications in real time.
An algorithm framework which can fulfill the aforementioned requirements was proposed by Tessler and Spangler [16]. This algorithm, named the inverse Finite Element Method (iFEM), employs the weighted least squares variational principle to construct a moderate relationship between the strain measurements and the deformations of a structure such as a frame (truss and beam), plate, and shell, etc. On the basis of the iFEM framework and the kinematic assumptions of Timoshenko beam theory, Gherlone et al. formulated a robust inverse finite beam element for the purpose of sensing three-dimensional (3D) deformation of a beam/frame structure [15,17]. In References [18,19], the optimal placement of sensors was researched by Zhao et al. to maintain the stability and accuracy of deformation reconstruction using iFEM with an inverse beam element. The experimental validations showed that iFEM is able to accurately and effectively estimate the deformations of a three-dimensional frame/shell/plate, and even the composite structure undergoing static and/or dynamic damped harmonic excitations without any knowledge of material, inertial, loading, or damping structural properties [20].
According to the advantages of the method described above, iFEM is employed to derive the measurement model for the description of the relationship between the beam deformations (the end-node shifts of the four inside beams) and the in situ strain measurements in the present paper. To verify the effectiveness and feasibility of the measurement scheme, a simple physical structure capable of modeling the whole telescope was manufactured ( Figure 2a). The performances of the simple physical structure model are presented as follows: the whole frame structure is made of an aluminum alloy with a Young's modulus of E = 7300 MPa, a Poisson's ration of v = 0.3, and a density of ρ = 2557 kg/m 3 . The weight of the whole frame structure is 43 kg (Figure 2b). The whole frame structure is composed of four supporting legs (the length of each leg is 1000 mm, and the diameter is 20 mm), a fixed platform (Figure 2c  An algorithm framework which can fulfill the aforementioned requirements was proposed by Tessler and Spangler [16]. This algorithm, named the inverse Finite Element Method (iFEM), employs the weighted least squares variational principle to construct a moderate relationship between the strain measurements and the deformations of a structure such as a frame (truss and beam), plate, and shell, etc. On the basis of the iFEM framework and the kinematic assumptions of Timoshenko beam theory, Gherlone et al. formulated a robust inverse finite beam element for the purpose of sensing three-dimensional (3D) deformation of a beam/frame structure [15,17]. In References [18,19], the optimal placement of sensors was researched by Zhao et al. to maintain the stability and accuracy of deformation reconstruction using iFEM with an inverse beam element. The experimental validations showed that iFEM is able to accurately and effectively estimate the deformations of a threedimensional frame/shell/plate, and even the composite structure undergoing static and/or dynamic damped harmonic excitations without any knowledge of material, inertial, loading, or damping structural properties [20].
According to the advantages of the method described above, iFEM is employed to derive the measurement model for the description of the relationship between the beam deformations (the endnode shifts of the four inside beams) and the in situ strain measurements in the present paper. To verify the effectiveness and feasibility of the measurement scheme, a simple physical structure capable of modeling the whole telescope was manufactured ( Figure 2a). The performances of the simple physical structure model are presented as follows: the whole frame structure is made of an aluminum alloy with a Young's modulus of E = 7300 MPa, a Poisson's ration of v = 0.3, and a density of ρ = 2557 kg/m 3 . The weight of the whole frame structure is 43 kg (Figure 2b). The whole frame structure is composed of four supporting legs (the length of each leg is 1000 mm, and the diameter is 20 mm), a fixed platform (Figure 2c

Construction of the Measurement Model
As shown in Figure 2, each supporting leg is regarded as a straight beam member of a constant circle cross-section. The inverse finite beam element quoted from References [15,17] is used to construct the measurement model. The process of constructing the measurement model is divided into four key steps: (1) the construction of the inverse finite beam element; (2) the computation of the section strains from the strain measurements; (3) the optimization of the strain sensor placement; (4) the construction of the measurement model in the global coordinate system.

Construction of the Inverse Finite Beam Element
This step is mainly quoted from References [15,17]. Based on the assumption of Timoshenko theory, the displacement vector, which incorporates the three-dimensional Cartesian coordinates, is shown in Figure 3. The shift of a certain node, B, is presented as follows: where denotes the nodal degrees-of-freedom. In light of Equation (2), the section strains vector ( ) is obtained from the following equation:

Construction of the Measurement Model
As shown in Figure 2, each supporting leg is regarded as a straight beam member of a constant circle cross-section. The inverse finite beam element quoted from References [15,17] is used to construct the measurement model. The process of constructing the measurement model is divided into four key steps: (1) the construction of the inverse finite beam element; (2) the computation of the section strains from the strain measurements; (3) the optimization of the strain sensor placement; (4) the construction of the measurement model in the global coordinate system.

Construction of the Inverse Finite Beam Element
This step is mainly quoted from References [15,17]. Based on the assumption of Timoshenko theory, the displacement vector, which incorporates the three-dimensional Cartesian coordinates, is shown in Figure 3. The shift of a certain node, B, is presented as follows: where u x (x, y, z), v y (x, y, z), u z (x, y, z) indicates the shift of node B and (x, y, z) are the Cartesian coordinates of node B; u(x, y, z) = u(x), v(x), w(x), θ x (x), θ y (x), θ z (x) represents the kinematic variables vector denoting the change of the node along the center axis (x ∈ [0, L], y = z = 0), i.e. node A. In the finite element framework, the arbitrary-node kinematic variables vector u(x, y, z) is interpolated by the shape function N(x): where u e denotes the nodal degrees-of-freedom. In light of Equation (2), the section strains vector e(u e ) is obtained from the following equation: where the matrix B(x) contains the derivatives of the shape function N(x). However, e(u e ) cannot be directly obtained from the strain sensors. In iFEM, e(u e ) is replaced with the in situ section strains vector e ε by minimizing the weighted least squares function ϕ: where e ε is computed from the surface strain measurements. Expanding the above least squares function results in the quadratic form: where C e is a constant vector; k e and f e are defined in the following: where L is the length of the beam element; n and x i (0 ≤ x i ≤ L) are, respectively, the number and the axial coordinate of the locations where the section strains are calculated. In the following, the minimization of the function ϕ(u) in Equation (5) in terms of u e leads to the inverse finite beam element: (e) (f)

Construction of the Measurement Model
As shown in Figure 2, each supporting leg is regarded as a straight beam member of a constant circle cross-section. The inverse finite beam element quoted from References [15,17] is used to construct the measurement model. The process of constructing the measurement model is divided into four key steps: (1) the construction of the inverse finite beam element; (2) the computation of the section strains from the strain measurements; (3) the optimization of the strain sensor placement; (4) the construction of the measurement model in the global coordinate system.

Construction of the Inverse Finite Beam Element
This step is mainly quoted from References [15,17]. Based on the assumption of Timoshenko theory, the displacement vector, which incorporates the three-dimensional Cartesian coordinates, is shown in Figure 3. The shift of a certain node, B, is presented as follows: where denotes the nodal degrees-of-freedom. In light of Equation (2), the section strains vector ( ) is obtained from the following equation:

Computation of the Section Strains from the Strain Measurements
The section strains are usually computed from the strain measurements through the following equation [15]: is the in situ section strains vector at location x i along the x-axis; m is the number of the sensors used to capture the surface strains of the beam; and R is the external radius of the beam element. ε(x i , θ i , β i ) denotes the measured strain at the location (x i , θ i , β i ), which is expressed in the cylindrical coordinate system (Figure 4). T(x i , θ i , β i ) is used to define the transformational relationship between the surface strain measurements and the section strains. where ( ) = { 1 ( ), 2 ( ), … , 6 ( )|( = 1, 2, … , )} is the in situ section strains vector at location along the x-axis; is the number of the sensors used to capture the surface strains of the beam; and R is the external radius of the beam element. ( , , ) denotes the measured strain at the location ( , , ) , which is expressed in the cylindrical coordinate system ( Figure 4).
( , , ) is used to define the transformational relationship between the surface strain measurements and the section strains. When only Equation (8) is used, six strain sensors need to be placed in one section for the calculation of the section strains vector ; thus, the total number of the strain sensors used to capture the surface strain measurements is 6×n in one inverse finite beam element. Herein, n is the minimum number of the sections where the section strains are evaluated, which is different under different When only Equation (8) is used, six strain sensors need to be placed in one section for the calculation of the section strains vector e ε ; thus, the total number of the strain sensors used to capture the surface strain measurements is 6 × n in one inverse finite beam element. Herein, n is the minimum number of the sections where the section strains are evaluated, which is different under different loading cases [21]. For the end-node load, n equals 2, whereas for the uniformly distributed load, n equals 3. With regard to the environmental loads, the deformations of four supporting legs mainly arise from the variations in gravity of the sub-reflector and the fixed platform during the pitching motion of the telescope. Also, the deformations of the four supporting legs are caused by the wind disturbance exerted on the surface of the sub-reflector. These loads are regarded as the end-node forces that are exerted on the four supporting legs. Therefore, the minimum number of sections is 2; the mini-number of the strain sensors is 12 in one inverse finite beam element. It is lucky that the number of the sensors can be reduced with the use of constitutive equations (Equation (9), refer to Figure 5) and equilibrium equations (Equation (10)).
where the section forces (N, Q y and Q z ) and the moments (M x , M y and M z ) are related to the section strains e i (x) (i = 1, . . . , 6). A x = EA is the axial rigidity, where A is the area of the cross-section of the beam element. G y = k 2 y GA and G z = k 2 z GA are the shear rigidities with k 2 y and k 2 z denoting the shear correction factors, where G is the shear modulus. For the solid section, k 2 y = k 2 z = 0.887. J x = GI p is the torsional rigidity. D y = EI y . and D z = EI z . are the bending rigidities [15].
Sensors 2019, 19, x FOR PEER REVIEW 6 of 14 loading cases [21]. For the end-node load, n equals 2, whereas for the uniformly distributed load, n equals 3. With regard to the environmental loads, the deformations of four supporting legs mainly arise from the variations in gravity of the sub-reflector and the fixed platform during the pitching motion of the telescope. Also, the deformations of the four supporting legs are caused by the wind disturbance exerted on the surface of the sub-reflector. These loads are regarded as the end-node forces that are exerted on the four supporting legs. Therefore, the minimum number of sections is 2; the mini-number of the strain sensors is 12 in one inverse finite beam element. It is lucky that the number of the sensors can be reduced with the use of constitutive equations (Equation (9), refer to Figure 5) and equilibrium equations (Equation (10)).
= is the axial rigidity, where is the area of the cross-section of the beam element. = 2 and = 2 are the shear rigidities with 2 and 2 denoting the shear correction factors, where G is the shear modulus. For the solid section, 2 = 2 = 0.887. = is the torsional rigidity.
= and = are the bending rigidities [15]. When the forms of the distribution loads q x , q y , and q z are known, the forms of the section strains can be estimated [15]. For the end-node loads, the section stains e ε 1 (x i ), e ε 4 (x i ), e ε 5 (x i ), and e ε 6 (x i ) are constant, while e ε 2 (x i ) and e ε 3 (x i ) are linear. In References [18,19], the section strains are expressed as: or in a matrix equation: where p1 = [a 1 , a 2 , a 3 , a 4 , a 5 , a 6 ] T is a constant parameters vector; T1 x i defines the transfer matrix between p1 and the section strains vector e ε (x i ).
When the parameters vector p1 is solved, the arbitrary section strains vector along the x-axis is calculated through Equation (13).
Finally, the kinematic variables vector u e is calculated through Equation (7) when two section strains vectors with respect to the two different cross-sections are determined. Subsequently, the deformation of an arbitrary node on the beam element surface is evaluated through Equation (1).

Optimization of the Strain Sensor Placement
The estimation of the structural deformation from the strain measurements is ill-posed when the placement of strain sensors or the boundary condition is not appropriate [16]. For instance, when all the strain sensors are set parallel to the generatix of the beam (all the β i (i = 1, . . . ,m) are set to the same value 0 o in Equation (13)), the transformation matrix T(x i , θ i , β i )T1(x i ) in Equation (13) is singular, which causes the estimation of the beam structural deformation from the strain measurements by Equation (7) to be ill-posed. Zhao et. al, discussed the influence of the above ill-posed problem on the estimation of the beam structural deformation by Equation (7), and constructed an optimal placement model of strain sensors to maintain the stability and accuracy of the estimation of the beam structural deformation. Using the sensor placement scheme obtained from the aforementioned optimal placement model, the deformation of the beam/frame structure can be accurately and steadily reconstructed even if differences exist for strain sensor placement and strain measurement [18,19]. In this paper, the placement of strain sensors is quoted from References [18,19] (refer to Table 1). The coordinate of every sensor in Table 1 is in the coordinate system of Figure 4, which is a local coordinate system (refer to Figure 6). Table 1. Description of the optimal placement for six strain sensors; angles are expressed in degrees.

Description
Orientation (θ, β) at x i = L/5 Orientation (θ, β) at x i = 4L/5 The relationship between the local coordinate system of the j th supporting leg and the global coordinate system can be directly expressed as a direction cosine matrix:

Six Strain Sensors
where , , and ( = , , ) are the direction cosines of the local k-axis. Accordingly, the calculation model of the j th supporting leg in the global coordinate system can be obtained as follows: where ′ is used to denote the deformation of the j th supporting leg in the global coordinate system.

Verification of the Model
In this section, a static loading test is conducted on a simple physical structure to verify the effectiveness and accuracy of the contact measurement scheme proposed in this paper. The real shifts of the fixed platform under different states are obtained using a laser tracking system, which comprises a Spherically Mounted Retro reflector (SMR, Figure 7a), seven different target holders (three target holders are installed on the main reflector to establish the principal coordinate system, while the other four are installed on the fixed platform to confirm the position and orientation of the fixed platform), and a laser tracking system (LTS, API Tracker 3, Automated Precision Inc., Rockwell, MD, USA). The resolution of the LTS is 1 µ m and the accuracy is related to the distance between the LTS and the measured object (the ratio is 5 µ m/m). In our test, the distance is 2 m, and the accuracy is 10 µ m. The in situ strain is captured by the bare FBG strain sensor attached to the beam surface with alpha-cyanoacrylate glue (Figure 7b). To monitor the temperature variation, a thermometer (Figure 7c, 51 Series II, Fluke, Avery, WA, USA) is used to measure the temperature variation on the surface of the four supporting legs each half hour. In Reference [22], when the initial and current temperatures are known, the influence on the strain measurement can be estimated with Equation (16):

Construction of the Measurement Model in a Global Coordinate System
As shown in Figures 3-5, the coordinate system for the supporting leg reconstruction is a local coordinate system, where the x-axis is positioned along the centroidal axis of the beam element; y and z are the principal inertial axes of the cross-section. For the four supporting legs, there are four different local coordinate systems depicted in Figure 6, so it is necessary that these four local coordinate systems be unified into a global coordinate system. In this paper, the setup of the x-axis and the y-axis of the global coordinate system is conducted on the main reflector; the z-axis is perpendicular to the main reflector ( Figure 6).
The relationship between the local coordinate system of the j th supporting leg and the global coordinate system can be directly expressed as a direction cosine matrix: where l k , m k , and n k (k = x, y, z) are the direction cosines of the local k-axis. Accordingly, the calculation model of the j th supporting leg in the global coordinate system can be obtained as follows: where u e j is used to denote the deformation of the j th supporting leg in the global coordinate system.

Verification of the Model
In this section, a static loading test is conducted on a simple physical structure to verify the effectiveness and accuracy of the contact measurement scheme proposed in this paper. The real shifts of the fixed platform under different states are obtained using a laser tracking system, which comprises a Spherically Mounted Retro reflector (SMR, Figure 7a), seven different target holders (three target holders are installed on the main reflector to establish the principal coordinate system, while the other four are installed on the fixed platform to confirm the position and orientation of the fixed platform), and a laser tracking system (LTS, API Tracker 3, Automated Precision Inc., Rockwell, MD, USA). The resolution of the LTS is 1 µm and the accuracy is related to the distance between the LTS and the measured object (the ratio is 5 µm/m). In our test, the distance is 2 m, and the accuracy is 10 µm. The in situ strain is captured by the bare FBG strain sensor attached to the beam surface with alpha-cyanoacrylate glue (Figure 7b). To monitor the temperature variation, a thermometer (Figure 7c, 51 Series II, Fluke, Avery, WA, USA) is used to measure the temperature variation on the surface of the four supporting legs each half hour. In Reference [22], when the initial and current temperatures are known, the influence on the strain measurement can be estimated with Equation (16): where ∆λ i and λ ini(i) are, respectively, the wavelength variation affected by the temperature ∆T j and the initial wavelength of the i th FBG strain sensor attached to the surface of the j th supporting leg. α is the thermal expansion of silica; η is the thermos-optic coefficient of the FBG sensor.
where ∆λ and ( ) are, respectively, the wavelength variation affected by the temperature Δ and the initial wavelength of the i th FBG strain sensor attached to the surface of the j th supporting leg.
is the thermal expansion of silica; is the thermos-optic coefficient of the FBG sensor. During the test, the in situ strain data are obtained from the strain measurement system, which comprises FBG strain sensors (Fiber Bragg Grating| os1100, Micron Optics, Atlanta, GA, USA) and an FBG interrogator (Optical Sensing Instrument| Si 155, Micron Optics, Atlanta, GA, USA). For each supporting leg, six FBG strain sensors (the range of the initial wavelength is (1527 nm, 1560 nm)) are placed on the surface of the beam according to the optimal placement scheme (see Table 1). Because four supporting legs possess the same performances, the placements of the strain sensors on four supporting legs are the same. Thus, the total number of the strain sensors used in the test is 24.
As mentioned in Section 1, the deformations of the four supporting legs give rise to the fixed platform shift. The real end-node shifts of the four supporting legs are computed using Equation (17): where ( ) = ( , , ) is the end-node (the intersection between the center axis of the supporting leg and the fixed platform) shift of the j th supporting leg in the three-dimensional coordinate system. The superscript 'end' denotes the end-node coordinate in the current global coordinate system, whereas ' ' denotes the initial end-node coordinate in the initial global coordinate system. During the test, the in situ strain data are obtained from the strain measurement system, which comprises FBG strain sensors (Fiber Bragg Grating| os1100, Micron Optics, Atlanta, GA, USA) and an FBG interrogator (Optical Sensing Instrument| Si 155, Micron Optics, Atlanta, GA, USA). For each supporting leg, six FBG strain sensors (the range of the initial wavelength is (1527 nm, 1560 nm)) are placed on the surface of the beam according to the optimal placement scheme (see Table 1). Because four supporting legs possess the same performances, the placements of the strain sensors on four supporting legs are the same. Thus, the total number of the strain sensors used in the test is 24.
As mentioned in Section 1, the deformations of the four supporting legs give rise to the fixed platform shift. The real end-node shifts of the four supporting legs are computed using Equation (17): where (δ) j = u j , v j , w j is the end-node (the intersection between the center axis of the supporting leg and the fixed platform) shift of the j th supporting leg in the three-dimensional coordinate system. The superscript 'end' denotes the end-node coordinate in the current global coordinate system, whereas 'ini' denotes the initial end-node coordinate in the initial global coordinate system. As shown in Figure 8, the fixed platform shift can be computed by using the shifts of the end-nodes of the four supporting legs through the following equations: where iFEM(δ) j (j = 1, . . . , 4) is the end-node shift of the j th supporting leg computed from the in situ strain measurements with iFEM, while LTS(δ) j is the end-node shift of the j th supporting leg captured from the LTS. iFEM(δ) o is the shift of the fixed platform computed from iFEM(δ) j , while LTS(δ) o is the real shift of the fixed platform computed from LTS(δ) j . As shown in Figure 8, the fixed platform shift can be computed by using the shifts of the endnodes of the four supporting legs through the following equations: where iFEM( ) ( = 1, … , 4) is the end-node shift of the j th supporting leg computed from the in situ strain measurements with iFEM, while ( ) is the end-node shift of the j th supporting leg captured from the LTS. iFEM( ) is the shift of the fixed platform computed from iFEM( ) , while ( ) is the real shift of the fixed platform computed from ( ) . The in situ strain measurements are captured from the bare FBG strain sensors. Different from traditional strain gauges such as resistance strain gauges, FBG sensors capture the structure strains based on the light wavelength shifts which are caused by FBG deformations due to the tension/compressive force or the temperature variation. In Reference [22], the strain is calculated through Equation (19) without considering the effect of temperature variation.
where ( ) and ( ) are the current wavelength and the initial wavelength of the i th FBG sensor, respectively. Δ is the difference of the i th FBG sensor caused by the temperature variation, estimated using Equation (16). P e is the effective photo-elastic coefficient of the fiber; the strain measurement is expressed as a micro strain. In the test, the wavelength of the FBG sensor is defined as the initial wavelength at the initial state of the supporting structure (Figure 9a), whereas the wavelength of the FBG sensor is defined as the current wavelength when the supporting structure is leaning at a certain angle (Figures 9b-d). The typical strain measurements at the pitch angle 30° are shown in Table 2. The in situ strain measurements are captured from the bare FBG strain sensors. Different from traditional strain gauges such as resistance strain gauges, FBG sensors capture the structure strains based on the light wavelength shifts which are caused by FBG deformations due to the tension/compressive force or the temperature variation. In Reference [22], the strain is calculated through Equation (19) without considering the effect of temperature variation.
where λ end(i) and λ ini(i) are the current wavelength and the initial wavelength of the i th FBG sensor, respectively. ∆λ i is the difference of the i th FBG sensor caused by the temperature variation, estimated using Equation (16). P e is the effective photo-elastic coefficient of the fiber; the strain measurement ε i is expressed as a micro strain. In the test, the wavelength of the FBG sensor is defined as the initial wavelength at the initial state of the supporting structure (Figure 9a), whereas the wavelength of the FBG sensor is defined as the current wavelength when the supporting structure is leaning at a certain angle (Figure 9b-d). The typical strain measurements at the pitch angle 30 • are shown in Table 2.
respectively. Δ is the difference of the i th FBG sensor caused by the temperature variation, estimated using Equation (16). P e is the effective photo-elastic coefficient of the fiber; the strain measurement is expressed as a micro strain. In the test, the wavelength of the FBG sensor is defined as the initial wavelength at the initial state of the supporting structure (Figure 9a), whereas the wavelength of the FBG sensor is defined as the current wavelength when the supporting structure is leaning at a certain angle (Figures 9b-d). The typical strain measurements at the pitch angle 30° are shown in Table 2.   Table 3 presents the corresponding real shifts of the fixed platform ( ( ) ) measured from the LTS and the evaluated shifts of the fixed platform ( iFEM( ) ) computed from the in situ strain measurements with iFEM. Also, the accuracies of the evaluated shifts are assessed by means of the absolute error and the percent difference:    Table 3 presents the corresponding real shifts of the fixed platform (LTS(δ) o ) measured from the LTS and the evaluated shifts of the fixed platform (iFEM(δ) o ) computed from the in situ strain measurements with iFEM. Also, the accuracies of the evaluated shifts are assessed by means of the absolute error and the percent difference: Furthermore, the orientation of the fixed platform can be calculated from the shift of the fixed platform and the shifts of the intersections (see Figure 8) of the four supporting legs using a discrete cosine algorithm (see Table 4). Table 4. The comparison of the fixed platform orientations at different pitch angles.

Pitch Angle
x-Rotation y-Rotation z-Rotation  Table 3, it was found that the shifts of the fixed platform can be accurately evaluated from the strain measurements with iFEM at different pitch angles. Specifically, the percent differences of the shifts between iFEM calculations and LTS captures are below 6% along the maximum-shift direction (y-direction). Although the percent differences of the shifts along the other two directions (x-and zdirections) are bigger than those along the y-direction, the absolute shifts along the xand zdirections are below 0.1 mm, smaller than those along the y-direction. Moreover, the absolute differences between iFEM calculations and LTS captures are less than or equal to 0.08 mm. From Table 4, it was found that the percent differences of the orientations are below 1% for x-rotation and z-rotation. Although the percent differences of the orientations for y-rotation are bigger than those for x-rotation and z-rotation, the absolute differences of the orientations for y-rotation are below 0.02 • .

Conclusions
With aim of providing the exact adjustment values for the sub-reflector adjustment structure of a large radio telescope, a contact measurement scheme was proposed to monitor the position and orientation of the sub-reflector in real time. In the scheme, several FBG strain sensors are attached to the surfaces of the four supporting legs and the deformations of these supporting legs are computed from in situ strain measurements to evaluate the position and orientation of the fixed platform. In order to verify the accuracy and effectiveness of the contact measurement scheme, a static test was conducted on a simple physical structure. The experimental results demonstrate that the contact measurement scheme proposed in this paper can effectively and accurately evaluate the shift of the fixed platform, offering a new way to monitor the shift of a radio telescope sub-reflector. Acknowledgments: Special appreciation goes to the editor and the anonymous reviewers for their careful reviews and valuable comments.

Conflicts of Interest:
The authors declare no conflict of interest.

Appendix A
U 0 , V 0 are, respectively, defined as the initial position and initial orientation of the sub-reflector; the corresponding changes are ∆U and ∆V in the case of pitching. Thus, the new position and orientation of the sub-reflector are: where R 1 is used to define the transform matrix between the initial state and the current state of the fixed platform. In order to maintain the original position and orientation of the sub-reflector, the length of the Stewart platform leg will change through Equation (A2): where R 2 is the transform matrix between the initial state of the sub-reflector and the current state of the fixed platform. R 3 is the transform matrix between the fixed platform and the sub-reflector under the initial state. P i and B i are used to define the coordinates of the i th spherical joint and the coordinate of the i th Hooke joint, respectively ( Figure A1). Acknowledgments: Special appreciation goes to the editor and the anonymous reviewers for their careful reviews and valuable comments.

Conflicts of Interest:
The authors declare no conflict of interest.
Appendix A 0 , 0 are, respectively, defined as the initial position and initial orientation of the sub-reflector; the corresponding changes are ∆ and ∆ in the case of pitching. Thus, the new position and orientation of the sub-reflector are: where 1 is used to define the transform matrix between the initial state and the current state of the fixed platform. In order to maintain the original position and orientation of the sub-reflector, the length of the Stewart platform leg will change through Equation (A2): ∆ = ‖ 2 + 1 − ‖ − ‖ 3 + 0 − ‖ = 1,2, … ,6 where 2 is the transform matrix between the initial state of the sub-reflector and the current state of the fixed platform. 3 is the transform matrix between the fixed platform and the sub-reflector under the initial state. and are used to define the coordinates of the i th spherical joint and the coordinate of the i th Hooke joint, respectively ( Figure A1). Figure A1. The coordinate systems of the fixed platform and the sub-reflector.