Next Article in Journal
Smartphone-Based Activity Recognition for Indoor Localization Using a Convolutional Neural Network
Previous Article in Journal
Tell the Device Password: Smart Device Wi-Fi Connection Based on Audio Waves
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

1
Key Laboratory of Electronic Equipment Structure Design of Ministry of Education, Xidian University, Xi’an 710071, China
2
Xinjiang Observatory, National Astronomical Observatories, Chinese Academy of Sciences, Urumqi 830011, China
*
Author to whom correspondence should be addressed.
Sensors 2019, 19(3), 619; https://doi.org/10.3390/s19030619
Submission received: 25 December 2018 / Revised: 19 January 2019 / Accepted: 23 January 2019 / Published: 1 February 2019
(This article belongs to the Section Physical Sensors)

Abstract

:
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.

1. 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 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).
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), a main reflector (Figure 2d), and a base (Figure 2e). The sub-reflector structure is composed of five iron plates. The weight of each iron plate (Figure 2f) is 7.5 kg and the total weight of the sub-reflector is 37.5 kg.

2. 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.

2.1. 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:
{ u x ( x ,   y , z ) = u ( x ) + z θ y ( x ) y θ z ( x ) v y ( x , y , z ) = v ( x ) z θ x ( x ) u z ( x , y , z ) = w ( x ) + y θ x ( x )
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 ) :
u ( x ,   y , z ) = N ( x ) u e
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:
e ( u e ) = B ( x ) u e
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 φ :
φ ( u ) = e ( u ) e ε 2
where e ε is computed from the surface strain measurements. Expanding the above least squares function results in the quadratic form:
φ ( u ) = 1 2 ( u e ) T k e u e ( u e ) T f e + C e
where C e is a constant vector; k e and f e are defined in the following:
k k e = L n i = 1 n [ B k T ( x i ) B k ( x i ) ] f k e = L n i = 1 n [ B k T ( x i ) e k ε ( x i ) ]
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:
k e u e = f e

2.2. Computation of the Section Strains from the Strain Measurements

The section strains are usually computed from the strain measurements through the following equation [15]:
ε ( x i , θ i , β i ) = e 1 ε ( x i ) ( c β 2 v s β 2 ) + e 2 ε ( x i ) ( c β 2 v s β 2 ) s θ R + e 3 ε ( x i ) ( c β 2 v s β 2 ) c θ R + e 4 ε ( x i ) c β s β c θ e 5 ε ( x i ) c β s β s θ + e 6 ε ( x i ) c β s β R = [ c β 2 v s β 2 , ( c β 2 v s β 2 ) s θ R , ( c β 2 v s β 2 ) c θ R , c β s β c θ , c β s β s θ , c β s β R ] × e ε ( x i ) = T ( x i , θ i , β i ) × e ε ( x i ) with   c β cos β i , s β s i n β i , c θ c o s θ i , s θ s i n θ i
where e ε ( x i ) = { e 1 ε ( x i ) ,   e 2 ε ( x i ) , ,   e 6 ε ( x i ) | ( i = 1 ,   2 ,   ,   m ) } 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.
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)).
N = A x e 1 Q y = G y e 5 Q z = G z e 4 M x = J x e 6 M y = D y e 2 M z = D z e 3
where the section forces ( N ,   Q y   a n d   Q z ) and the moments ( M x ,   M y   a n d   M z ) are related to the section strains e i ( x ) (i = 1, …,6). A x = E A is the axial rigidity, where A is the area of the cross-section of the beam element. G y = k y 2 G A and G z = k z 2 G A are the shear rigidities with k y 2 and k z 2 denoting the shear correction factors, where G is the shear modulus. For the solid section, k y 2 = k z 2 = 0.887. J x = G I p is the torsional rigidity. D y = E I y . and D z = E I z . are the bending rigidities [15].
N x + q x = 0 Q y x + q y = 0 Q z x + q z = 0 M x x = 0 M y x Q z = 0 M z x Q y = 0
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:
e 1 ε ( x i ) = a 1 e 2 ε ( x i ) = a 2 x i + a 4 e 4 ε ( x i ) = m 1 a 2 e 6 ε ( x i ) = a 6 e 3 ε ( x i ) = a 3 x i + a 5 e 3 ε ( x i ) = a 3 x i + a 5
or in a matrix equation:
e ε ( x i ) = { e 1 ε ( x i ) , e 2 ε ( x i ) , , e 6 ε ( x i ) } T = [ 1 0 0 0 0 0 0 x i 0 0 0 0 0 0 x i 0 1 0 0 m 1 0 0 0 0 0 0 m 1 0 0 0 0 0 0 0 0 1 ] × [ a 1 , a 2 , a 3 , a 4 , a 5 , a 6 ] T = T 1 ( x i ) × p 1 , with   m 1 = E I z G k y 2 A ,   m 2 = E I y G k z 2 A
where p 1 = [ a 1 ,   a 2 ,   a 3 ,   a 4 ,   a 5 ,   a 6 ] T is a constant parameters vector; T 1 x i defines the transfer matrix between p 1 and the section strains vector e ε ( x i ) .
When the parameters vector p 1 is solved, the arbitrary section strains vector along the x-axis is calculated through Equation (13).
e ε ( x k ) = T 1 ( x j ) × ( T ( x i , θ i , β i ) T 1 ( x i ) ) 1 × ε ( x i ,   θ i ,   β i )      ( i = 1 ,   ,   6 ;   k = 1 ,   2 )
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).

2.3. 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 ) T 1 ( 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).

2.4. Construction of the Measurement Model in a Global Coordinate System

As shown in Figure 3, Figure 4 and Figure 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 jth supporting leg and the global coordinate system can be directly expressed as a direction cosine matrix:
Λ j = [ l x m x n x l y m y n y l z m z n z ] i ( j = 1 , 2 , 3 , 4 )
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 jth supporting leg in the global coordinate system can be obtained as follows:
k j e u j e = f j e ( j = 1 ,   2 ,   3 ,   4 ) with   k j e = Φ j T k j e Φ j ,   f j e = Φ j T f j e , Φ j = [ Λ j 0 0 Λ j 0 0 0 0 0 0 0 0 Λ j 0 0 Λ j ]
where u j e is used to denote the deformation of the jth supporting leg in the global coordinate system.

3. 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):
Δ λ i = ( α + η ) Δ T j × λ i n i ( i )
where Δ λ i and λ i n i ( i ) are, respectively, the wavelength variation affected by the temperature Δ T j and the initial wavelength of the ith FBG strain sensor attached to the surface of the jth 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):
( δ ) j = δ j e n d δ j i n i ( j = 1 ,   2 , 3 , 4 )
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 jth supporting leg in the three-dimensional coordinate system. The superscript ‘end’ denotes the end-node coordinate in the current global coordinate system, whereas ‘ i n i ’ 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:
L T S ( δ ) o = j = 1 4 L T S ( δ ) j / 4 ,     iFEM ( δ ) o = j = 1 4 iFEM ( δ ) j / 4
where iFEM ( δ ) j   ( j = 1 ,   ,   4 ) is the end-node shift of the jth supporting leg computed from the in situ strain measurements with iFEM, while L T S ( δ ) j is the end-node shift of the jth supporting leg captured from the LTS. iFEM ( δ ) o is the shift of the fixed platform computed from iFEM ( δ ) j , while L T S ( δ ) o is the real shift of the fixed platform computed from L T S ( δ ) j .
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.
ε i = 1 K * ( λ e n d ( i ) λ i n i ( i ) Δ λ i λ i n i )   w i t h   K = 1 P e
where λ e n d ( i ) and λ i n i ( i ) are the current wavelength and the initial wavelength of the ith FBG sensor, respectively. Δ λ i is the difference of the ith 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.
Table 3 presents the corresponding real shifts of the fixed platform ( L T S ( δ ) 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:
Diff ( δ ) o = | iFEM ( δ ) o A P I ( δ ) o | % Diff ( δ ) o = 100 × [ | iFEM ( δ ) o A P I ( δ ) o | | A P I ( δ ) o | ]
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).
From 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 z- directions) are bigger than those along the y-direction, the absolute shifts along the x- and z- directions 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°.

4. 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.

Author Contributions

Y.Z. carried out the majority of work on the paper, constructing the measurement model for monitoring the sub-reflector shift in real time, verifying measurement scheme, and writing the article. J.D. constructed the finite element model of the sub-reflector. H.B. proposed this new measurement method for monitoring the shift of the sub-reflector in real time. Q.X. gave advice on the writing of the article.

Funding

This work is financially supported by the National Natural Science Foundation of China through Grants 51775401 and 51675398, the National Key Basic Research Program of China through Grant 2015CB857100, the CAS “Light of West China” Program (No. 2016-QNXZ-A-7), and the joint Research Fund in Astronomy (No. U1731135) through cooperation between the National Science Foundation of China (NSFC) and the Chinese Academy of Sciences (CAS).

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:
U 1 = R 1 U 0 + Δ U ,   V 1 = V 0 Δ V
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):
Δ L i = R 2 P i + U 1 B i R 3 P i + U 0 B i   i = 1 , 2 , , 6
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 ith spherical joint and the coordinate of the ith Hooke joint, respectively (Figure A1).
Figure A1. The coordinate systems of the fixed platform and the sub-reflector.
Figure A1. The coordinate systems of the fixed platform and the sub-reflector.
Sensors 19 00619 g0a1

References

  1. Pisanu, T.; Buffa, F.; Poppi, S.; Marongiu, P.; Serra, G.; Vargiu, G.P.; Concu, R. The SRT inclinometer for monitoring the rail and the thermal gradient effects on the alidade structure. In Proceedings of the SPIE Astronomical Telescopes + Instrumentation, Montréal, QC, Canada, 22–27 June 2014; Volume 91454R. [Google Scholar]
  2. Ries, P.; Hunter, T.R.; Constantikes, K.T.; Brandt, J.J.; Ghigo, F.D.; Mason, B.S.; Prestage, R.M.; Ray, J.; Schwab, F.R. Measuring and Correcting Wind-Induced Pointing Errors of the Green Bank Telescope Using an Optical Quadrant Detector. Publ. Astron. Soc. Pac. 2011, 123, 682–691. [Google Scholar] [CrossRef]
  3. Jiang, H.; Kesteven, M.; Wilson, W.; Li, C. A distributed control system for a radio telescope with six-meter hexapod mount. In Proceedings of the IEEE International Conference on Control and Automation, Edinburgh, UK, 9–11 December 2009; pp. 2003–2010. [Google Scholar]
  4. Pisanu, T.; Poppi, S. Architecture of the metrology for the SRT. Proc. SPIE 2012, 8444, 84442E1-9. [Google Scholar]
  5. Pisanu, T.; Buffa, F.; Concu, R.; Marongiu, P.; Pili, M.; Poppi, S.; Serra, G.; Urru, E.; Vargiu, G. A PSD (position sensing device) to map the shift and tilt of the SRT secondary mirror. Proc. SPIE 2014, 9145, 91454U1-9. [Google Scholar]
  6. Derkevorkian, A.; Masri, S.F.; Alvarenga, J.; Boussalis, H.; Bakalyar, J.; Richards, W.L. Strain-based deformation shape-estimation algorithm for control and monitoring applications. AIAA J. 2013, 51, 2231–2240. [Google Scholar] [CrossRef]
  7. Pak, C. Wing shape sensing from measured strain. AIAA J. 2016, 54, 1064–1073. [Google Scholar] [CrossRef]
  8. Akl, W.; Poh, S.; Baz, A. Wireless and distributed sensing of the shape of morphing structures. Sens. Actuators A Phys. 2006, 140, 94–102. [Google Scholar] [CrossRef]
  9. Sante, R.D. Fiber optic sensors for structural health monitoring of aircraft composite structures: Recent advances and applications. Sensors 2015, 15, 18666–18713. [Google Scholar] [CrossRef] [PubMed]
  10. Amanzadeh, M.; Aminossadati, S.M.; Kizil, M.S.; Rakić, A.D. Recent Developments in Fiber Optic Shape Sensing. Measurement 2018, 128, 119–137. [Google Scholar] [CrossRef]
  11. Vorathin, E.; Hafizi, Z.M.; Ghani, S.A.C.; Lim, K.S. Real-time monitoring system of composite aircraft wings utilizing Fibre Bragg Grating sensor. IOP Conf. Ser. Mater. Sci. Eng. 2016, 152, 012024. [Google Scholar] [CrossRef]
  12. Ko, W.L.; Richards, W.L.; Fleischer, V.T. Applications of the Ko displacement theory to the deformed shape predictions of the doubly-tapered Ikhana wing. In NASA/TP-2009-214652; NASN Dryden Flight Research Center: Edwards, CA, USA, 2009. [Google Scholar]
  13. Jutte, C.V.; Ko, W.L.; Stephens, C.A.; Al, E. Deformed Shape Calculation of a Full-Scale Wing Using Fiber Optic Strain Data from a Ground Loads Test. In NASA Langley Research Center, Rept, TP-215975; NASA Langley Research Center: Hampton, VA, USA, 2011. [Google Scholar]
  14. Bogert, P.; Haugse, E.; Gehrki, R. Structural shape identification from experimental strains using a modal transformation technique. In Proceedings of the 44th AIAA/ASME/ASCE/AHS Structures, Structural Dynamics, and Materials Conference, Norfolk, VA, USA, 7–10 April 2003; p. 1626. [Google Scholar]
  15. Gherlone, M.; Cerracchio, P.; Mattone, M.; Di Sciuva, M.; Tessler, A. Shape sensing of 3D frame structures using an inverse finite element method. Int. J. Solids Struct. 2012, 49, 3100–3112. [Google Scholar] [CrossRef]
  16. Tessler, A.; Spangler, J.L. A variational principle for reconstruction of elastic deformations in shear deformable plates and shells. In NASA Langley Research Center TM-212445; NASA Langley Research Center: Hampton, VA, USA, 2003. [Google Scholar]
  17. Gherlone, M.; Cerracchio, P.; Mattone, M.; Di Sciuva, M.; Tessler, A. An inverse finite element method for beam shape sensing: Theoretical framework and experimental validation. Smart Mater. Struct. 2014, 23, 1–13. [Google Scholar] [CrossRef]
  18. Zhao, Y.; Du, J.; Bao, H.; Xu, Q. Optimal Sensor Placement Based on Eigenvalues Analysis for Sensing Deformation of Wing Frame Using iFEM. Sensors 2018, 18, 2424. [Google Scholar] [CrossRef] [PubMed]
  19. Zhao, Y.; Du, J.; Bao, H.; Xu, Q. Optimal Sensor Placement for Inverse Finite Element Reconstruction of Three-Dimensional Frame Deformation. Int. J. Aerospace Eng. 2018, 2018, 1–10. [Google Scholar] [CrossRef]
  20. Gherlone, M.; Cerracchio, P.; Mattone, M. Shape sensing methods: Review and experimental comparison on a wing-shaped plate. Prog. Aerosp. Sci. 2018, 99, 14–26. [Google Scholar] [CrossRef]
  21. Zhao, Y.; Bao, H.; Duan, X.; Fang, H. The Application Research of Inverse Finite Element Method for Frame Deformation Estimating. Int. J. Aerosp. Eng. 2017, 2017, 1–8. [Google Scholar]
  22. Werneck, M.M.; Allil, R.C.S.B.; Ribeiro, B.A.; de Nazaré, F.V.B. A Guide to Fiber Bragg Grating Sensors. In Current Trends in Short-and Long-Period Fiber Gratings; INTECH: Vienna, Austria, 2013; pp. 1–24. [Google Scholar]
Figure 1. (a) Radio telescope and sensor placement, (b) bare fiber Bragg grating (FBG) strain sensor, and (c) FBG temperature sensor.
Figure 1. (a) Radio telescope and sensor placement, (b) bare fiber Bragg grating (FBG) strain sensor, and (c) FBG temperature sensor.
Sensors 19 00619 g001
Figure 2. The simple physical model of the radio telescope and its size: (a) Telescope model; (b) The size of the whole structure; (c) The size of the fixed platform; (d) The size of the main reflector; (e) The size of the pedestal; (f) The size of the iron plate.
Figure 2. The simple physical model of the radio telescope and its size: (a) Telescope model; (b) The size of the whole structure; (c) The size of the fixed platform; (d) The size of the main reflector; (e) The size of the pedestal; (f) The size of the iron plate.
Sensors 19 00619 g002aSensors 19 00619 g002b
Figure 3. Beam geometry and kinematic variables.
Figure 3. Beam geometry and kinematic variables.
Sensors 19 00619 g003
Figure 4. Location and coordinates of a strain sensor placed on an external beam surface [18].
Figure 4. Location and coordinates of a strain sensor placed on an external beam surface [18].
Sensors 19 00619 g004
Figure 5. Beam section forces and moments [18].
Figure 5. Beam section forces and moments [18].
Sensors 19 00619 g005
Figure 6. The coordinate system.
Figure 6. The coordinate system.
Sensors 19 00619 g006
Figure 7. The test system of the simple physical structure: (a) Spherically Mounted Retro reflector (SMR) and target holder; (b) Bare FBG strain sensor attached to the beam surface; (c) Thermometer; (d) Test system.
Figure 7. The test system of the simple physical structure: (a) Spherically Mounted Retro reflector (SMR) and target holder; (b) Bare FBG strain sensor attached to the beam surface; (c) Thermometer; (d) Test system.
Sensors 19 00619 g007
Figure 8. The top view of the fixed platform posture (the observation point is perpendicular to the main reflector): (a) The ideal posture; (b) The shift posture.
Figure 8. The top view of the fixed platform posture (the observation point is perpendicular to the main reflector): (a) The ideal posture; (b) The shift posture.
Sensors 19 00619 g008
Figure 9. The side views of the supporting structure under different states: (a) Initial state; (b) Pitch angle: 30°; (c) Pitch angle: 45°; (d) Pitch angle: 60°.
Figure 9. The side views of the supporting structure under different states: (a) Initial state; (b) Pitch angle: 30°; (c) Pitch angle: 45°; (d) Pitch angle: 60°.
Sensors 19 00619 g009aSensors 19 00619 g009b
Table 1. Description of the optimal placement for six strain sensors; angles are expressed in degrees.
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 = 4 L / 5
Six Strain Sensors(−120,0), (0,0), (120,0)(−120,0), (0,45), (120,0)
Table 2. The typical strain measurements at a pitch angle of 30°.
Table 2. The typical strain measurements at a pitch angle of 30°.
Sensor LocationSupporting Leg 1Supporting Leg 2Supporting Leg 3Supporting Leg 4
L / 5 , −120°,0°0.0000790.000112−0.000172−0.000166
L / 5 , 0°,45°0.0000190.0000280.0000490.000053
L / 5 , 120°,0°−0.000272−0.000190.0001320.000107
4 L / 5 , −120°,0°−0.000335−0.0002950.0002630.000261
4 L / 5 , 0°,45°−0.000024−0.0000150.000001−0.000046
4 L / 5 , 120°,0°0.0004290.000327−0.000395−0.000289
Table 3. The shifts measured using API Tracker 3 and the shifts computed from the strain data using the inverse Finite Element Method (iFEM).
Table 3. The shifts measured using API Tracker 3 and the shifts computed from the strain data using the inverse Finite Element Method (iFEM).
Pitch Angle Shift in x-DirectionShift in y-DirectionShift in z-Direction
30° L T S ( δ ) o −0.01 mm−0.62 mm0.08 mm
iFEM ( δ ) o 0.03 mm−0.59 mm0.0 mm
Diff ( δ ) o 0.04 mm0.03 mm0.08 mm
% Diff ( δ ) o 400%4.8%100%
45° L T S ( δ ) o −0.02 mm−1.82 mm0.03 mm
iFEM ( δ ) o 0.05 mm−1.73 mm0.01 mm
Diff ( δ ) o 0.07 mm0.09 mm0.02 mm
% Diff ( δ ) o 350%4.9%66.7%
60° L T S ( δ ) o 0.03 mm−2.06 mm−0.04 mm
iFEM ( δ ) o 0.08 mm−1.95 mm−0.04 mm
Diff ( δ ) o 0.05 mm0.11 mm0.0 mm
% Diff ( δ ) o 166.7%5.3%0.0%
Table 4. The comparison of the fixed platform orientations at different pitch angles.
Table 4. The comparison of the fixed platform orientations at different pitch angles.
Pitch Angle x-Rotationy-Rotationz-Rotation
30° L T S ( θ ) 1.626°−0.02°0.0°
iFEM ( θ ) 1.619°−0.013°0.0°
Diff ( θ ) 0.007°0.007°0.0°
% Diff ( θ ) 0.4%35%0.0%
45° L T S ( θ ) 1.722°−0.04°0.21°
iFEM ( θ ) 1.709°−0.027°0.226°
Diff ( θ ) 0.013°0.013°0.016°
% Diff ( θ ) 0.8%32.5%1.9%
60° L T S ( θ ) 1.721°−0.066°0.5°
iFEM ( θ ) 1.704°−0.048°0.504°
Diff ( θ ) 0.017°0.018°0.004°
% Diff ( θ ) 0.98%27.3%0.8%

Share and Cite

MDPI and ACS Style

Zhao, Y.; Du, J.; Xu, Q.; Bao, H. Real-Time Monitoring of the Position and Orientation of a Radio Telescope Sub-Reflector with Fiber Bragg Grating Sensors. Sensors 2019, 19, 619. https://doi.org/10.3390/s19030619

AMA Style

Zhao Y, Du J, Xu Q, Bao H. Real-Time Monitoring of the Position and Orientation of a Radio Telescope Sub-Reflector with Fiber Bragg Grating Sensors. Sensors. 2019; 19(3):619. https://doi.org/10.3390/s19030619

Chicago/Turabian Style

Zhao, Yong, Jingli Du, Qian Xu, and Hong Bao. 2019. "Real-Time Monitoring of the Position and Orientation of a Radio Telescope Sub-Reflector with Fiber Bragg Grating Sensors" Sensors 19, no. 3: 619. https://doi.org/10.3390/s19030619

APA Style

Zhao, Y., Du, J., Xu, Q., & Bao, H. (2019). Real-Time Monitoring of the Position and Orientation of a Radio Telescope Sub-Reflector with Fiber Bragg Grating Sensors. Sensors, 19(3), 619. https://doi.org/10.3390/s19030619

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop