Dynamic tracking performance of indoor global positioning system: An experimental and theoretical study

Transcription

1 Research Article Dynamic tracking performance of indoor global positioning system: An experimental and theoretical study Advances in Mechanical Engineering 015, Vol. 7(10) 1 8 Ó The Author(s) 015 DOI: / aime.sagepub.com Gang Zhao, Xishuang Jing, Liangjie Chen, Zhanlei Sun, Pengfei Han and Pengfei Zhang Abstract The automation level has been improved rapidly with the introduction of large-scale measurement technologies, such as indoor global positioning system, into the production process among the fields of car, ship, and aerospace due to their excellent measurement characteristics. In fact, the objects are usually in motion during the real measurement process; however, the dynamic measurement characteristics of indoor global positioning system are much limited and still in exploration. In this research, we focused on the dynamic tracking performance of indoor global positioning system and then successfully built a mathematical model based on its measurement principles. We first built single and double station system models with the consideration of measurement objects movement. Using MATLAB simulation, we realized the dynamic measurement characteristics of indoor global positioning system. In the real measurement process, the experimental results also support the mathematical model that we built, which proves a great success in dynamic measurement characteristics. We envision that this dynamic tracking performance of indoor global positioning system would shed light on the dynamic measurement of a motion object and therefore make contribution to the automation production. Keywords igps, dynamic measurement, mathematical model, experimental study Date received: 18 December 014; accepted: 1 September 015 Academic Editor: Elsa de Sa Caetano Introduction With the rapid development of modern science and technology, manufacturing process has been greatly improved and become with the features of informatization, automatization, intelligentialization, and integration. 1, In assembly process of manufacturing industry, people often use some new technologies to improve the automatization level. It is widely considered that the quality, precision, and efficiency of the measurement technology are significantly important in the automation assembly process since they determine the success of the assembly process. 3 As a traditional and popular large-scale measurement technology, laser tracker is widely used in the assembly process of manufacturing industry, 4,5 yet the small measuring range and single State Key Laboratory of Virtual Reality Technology and Systems, School of Mechanical Engineering and Automation, Beihang University, Beijing, P.R. China Corresponding authors: Liangjie Chen, State Key Laboratory of Virtual Reality Technology and Systems, School of Mechanical Engineering and Automation, Beihang University, Beijing , P.R. China. ljchen390711@163.com Zhanlei Sun, State Key Laboratory of Virtual Reality Technology and Systems, School of Mechanical Engineering and Automation, Beihang University, Beijing , P.R. China. sunzhanlei@buaa.edu.cn Creative Commons CC-BY: This article is distributed under the terms of the Creative Commons Attribution 3.0 License ( which permits any use, reproduction and distribution of the work without further permission provided the original work is attributed as specified on the SAGE and Open Access pages ( open-access-at-sage).

2 Advances in Mechanical Engineering point measuring of laser tracker still limit its industrial application. As a new large-scale measurement technology, indoor global positioning system (igps) provides a novel way to realize automation assembly due to its high precision, high efficiency, high reliability, and multi-point tracking ability at the same time. 6,7 It has been demonstrated that igps is used in various automation assembly processes among aerospace, automotive, and ship-building structures (e.g. the jigless assembly of aircraft components and the positioning of automatic robots). 8 It is also reported that the igps in some cases presented better results than a laser tracker. 9 In the past decades, lots of researches on static measurement characteristics of igps have been conducted. For instance, DA Maisano et al. 10 clearly stated that measurement accuracy of igps could improve if the number of transmitters increases, while the change will not be too much if the number of transmitters is four or more. In 009, JE Muelaner et al. 11 discussed the uncertainty of angle measurement for a single transmitter and gave out the relationship diagram. It was proved by R Schmitt et al. 1 that the layout of transmitters could significantly influence its measurement accuracy and the C Type Layout presented better precision. In fact, the objects are normally in motion during the assembly process which causes measurement deviation of igps based on its operational principles. For these applications, the objects need to be dynamically tracked in order to monitor the whole process. Most recently, some scholars have started to investigate the accuracy in igps dynamic measurement. For example, Claudia Depenthal 13 claimed that the dynamic measurement accuracy of igps would decrease if the measurement object is in rotational movement. Zheng Wang et al. 9 considered that the velocity direction of the measurement object could influence the increment direction of deviation. Unfortunately, there is little work on the qualitative relationship between the measurement deviation and the object s velocity, and the dynamic tracking performance of igps system still remains to be a big challenge. In this article, we focused on the dynamic tracking performance of igps in automation production process through theoretical simulation and experimental study. The emphasis of this article lied in the estimation of dynamic measurement accuracy by exploring the quantitative relationship between the object s velocity and the measurement angles relative to transmitters. First, we proposed a single station system model to compensate the angle deviation when the object is in motion. Then, we used the angles information relative to two transmitters to obtain the object s position by building a double station system model. Most importantly, the experimental study demonstrated that the deviation has a linear relationship with the object s velocity, which was in accordance with the MATLAB simulation, and both of the theoretical simulation and experimental study results indicated the great success of igps applied in dynamic measurements. We consider that the results of our reported work will improve the measurement accuracy of igps for tracking dynamic objects and guide the direction of practical industrial production when high precision of igps is needed for assembly process. Mathematics model An igps system typically consists of two or more transmitters, a control center, a number of wired/wireless sensors, and a scale bar. During the measurement process, transmitters, which operate as reference points, continually generate three signals: an infrared lightemitting diode (LED) strobe and two infrared laser fanned beams that rotate in the head of the transmitter. Each transmitter presents two measurement values to each sensor: the horizontal (azimuth, u) and the vertical angles (elevation, u), as shown in Figure 1(a). The position calculation engine (PCE) can calculate these measurement values related to transmitters trough pulsing signals, which are gained from the sensor being placed on the surface of the object. If the sensor is located in the line sight of two or more transmitters, whose coordinate locations have been solved before starting measuring, and the position of sensor can be calculated based on the measurement principle of triangulation, as shown in Figure 1(b). Considering that two fanned lasers sweeping through the measuring volume with a predetermined speed at the beginning of the measurement process and the strobe illuminating the volume, a sensor that is placed in the volume can receive two kinds of optical signals (flat light and synchronous light) and convert these optical signals into pulsing signals. During a measurement cycle, a sensor creates a pulsing signal to record the beginning moment (t 0 ) on receiving a synchronous signal and create two timing pulsing signals to record the time (t 1 and t ) when it receives two flat signals. Based on these pulsing signals and the postures of two flat planes, the azimuth angle and elevation angle of a sensor related to a transmitter can be calculated as presented in Figure. Since these two infrared laser fanned beams do not reach the sensor simultaneously, time delay will cause measurement deviation if the sensor which is set as a measurement point is in motion. In order to obtain the accurate position, a single station system model is proposed to compensate the deviation of two measurement values related to each transmitter. When a transmitter launches a synchronous signal, the intersecting line of the light plane and transmitter s horizontal plane is defined as x-axis, the rotation axis is defined as z axis, and the y axis can be determined by the right-hand rule.

3 Zhao et al. 3 Figure 1. Measurement principle of igps. (a) The transmitter presents two measurement values: the azimuth angle (u) and elevation angle (u). When the coordinate system of the transmitter is fixed, the position of sensor can be located in a line related to the transmitter based on these two angles. (b) If two or more transmitters are placed in the measuring volume, the accurate position of sensor is at the line intersection. Figure. Two infrared laser fanned beams and an infrared LED strobe of a transmitter. (a) A transmitter continually operates an infrared LED strobe and two infrared laser fanned beams during its working process. The angle between fan 1 (L 1 ) and fan (L )is about u off in the horizontal plane of the transmitter. (b) During a measurement cycle, a sensor receives three signals from each transmitter. When it receives a signal of the infrared LED strobe, it creates a pulse to record the time of t 0. The time of t 1 and t is to record the moment when the sensor receives the signals of fan 1 and fan. This measurement process can be seen in Figure 3 for a measurement cycle. If we bypass its motion, the direction vector of the measurement point can be given by equation (1) as follows 6 r = 4 r x r y r z = n 1 ðt 1 Þ3 n ðt Þ= i j k cos c 1 sin f 1 cos c 1 cos f 1 sin c 1 cos c sin f cos c cos f sin c ð1þ where f 1 and f are the angles between the fan and the x-axis when the point receives the signals of fan 1 and fan, and c 1 or c is the intersection angle which is between the fan and the z-axis. However, the point is often in motion during the measurement process. In this situation, when the point receives the signal of two infrared laser fanned beams, it is not at the original position of t 0 which is needed to be measured. Thus, the values of f 1 and f should be compensated considering the motion of the point. If the point that is far away from the original point in coordinate system (R) is in motion of the velocity (v), the accurate rotation angle of fan 1 aligned with the x-axis can be calculated by

4 4 Advances in Mechanical Engineering Figure 3. Measurement process for a single station. (a) At the beginning moment of t 0, the normals to the planes of fan 1 and fan are n 1 (t 0 ) and n (t 0 ). (b) Through t 1, fan 1 has rotated u 1 about the z axis and its normal turns n 1 (t 1 ). (c) Then the normal to fan becomes n (t ) at t. (d) The point lies on the first fan at t 1 and the second fan at t, and its direction vector related to the origin coordinate can be calculated through cross product. f 1 = v ðt 1 t 0 Þ Df ðþ where Df is the compensation angle when the velocity of the point is considered and can be expressed as Df = v m ðt 1 t 0 Þ ð3þ where v m is the angular velocity of the measurement point aligned with the rotation axis of the transmitter, which is positive value when its rotation direction is the same as that of the transmitter. Since the product of the velocity (v) and time (t) is much less than R, equation (3) can be transformed into Df = v m ðt 1 t 0 Þ= a v R ð t 1 t 0 Þ ð4þ where a is the coefficient, which is determined by the angle between the normal to fan 1 and the direction of velocity. And the values of f 1 and f can be given by f 1 = v ðt 1 t 0 Þ a v R ð t 1 t 0 Þ f = v ðt t 0 Þ b v R t ð5þ ð t 0 Þ f off where b is the coefficient, which is determined by the angle between the normal to fan and the direction of velocity, and f off is the angle between fan 1 and fan in the horizontal plane of the transmitter. Finally, if we consider the influence of the velocity, the azimuth angle (u) and elevation angle (u) of the measurement point related to the transmitter can be calculated by u = v R ðaðt 1 t 0 Þ+ bðt t 0 ÞÞ + v ðð t 1 t 0 Þ+ ðt t 0 ÞÞ f off u = arctan 0 sin v R ð a t 1 t ð Þ+ bðt t 0 ÞÞ 4 + v ðð t 1 t 0 Þ+ ðt t 0 ÞÞ f off 4 tan c ð6þ 1 A ð7þ

5 Zhao et al. 5 calculating the azimuth angle (u B ) and elevation angle (u B ) in the single station system of transmitter B, whose infrared LED strobe is supposed to sweep the point later. In this situation, the values of f 1 and f can be described as 8 < f 1 = v ðt 1 t 0 Þ a v R ðð t 1 t 0 Þ+ DtÞ : f = v ðt t 0 Þ b v R ðð t t 0 Þ+ Dt Þ f off where Dt is the time interval between the start time of these two transmitters. From equation (8), the position of M can be calculated by Figure 4. General double station system model. In the double station system model, there are two transmitters (A and B) in the measuring volume, and the position of M can be calculated accurately based on this model. From discussion above, we can see that the single station system model can obtain the point s relatively position to one transmitter; however, the accurate position must be determined by at least two transmitters. Therefore, we proposed a double station system model which has compensated angle deviation to calculate the accurate position of the point in the following section. For the double station system model, there are two transmitters (A and B) in the measuring volume, and the position of M can be calculated accurately because each coordinate value contains two components. For example, the x-axis value consists of x 1 and x. As presented in Figure 4, the distance between A and B is d in the horizontal direction, l in the y-axis direction, and h in the vertical direction. Thus, the components of the coordinate value of M can be given by 8 x 1 = ð l d tan u BÞ d tan u A tan u B x = ð l d tan u AÞ + d tan u A tan u B y 1 = ð l d tan u BÞ >< l 1 cot u A tan u B y = ð l d tan u ð8þ AÞ cot u B tan u A 1 + l ðl d tan u z 1 = B Þ cos u A ðsin u A cos u A tan u B Þ h ðl d tan u z = A Þ >: cos u B ðcos u B tan u A sin u B Þ + h where u A and u B are the azimuth angles related to transmitters A and B, and u A and u B are the elevation angles, which can be gained from equations (6) and (7). As the start time of these two transmitters are different, special processing should be carried out when x = x 1 + x, y = y 1 + y, z = z 1 + z ð9þ During the measurement process, the position of M is considered to be (x 0, y 0, z 0 ) T when it is static. However, if the velocity of M is v, the position of M is (x, y, z) T under the same measurement condition at this moment. Thus, the deviation of the coordinate values can be given by e = ðx, y, zþ T ðx 0, y 0, z 0 Þ T ð10þ To be noted, a sensor, which represents the measurement point, can receive signals from several transmitters in the measuring volume, and the system will calculate its position based on the signals from two random transmitters. In this consideration, the position of M may have several values. The system will combine those values together and obtain its final coordinate value. Results and discussion To investigate the dynamic tracking performance of igps, we first give the relationship between the dynamic measurement accuracy of igps and velocity of the measurement point using MATLAB simulation. For the working conditions of a single station system model, the point which is in motion with the velocity of v is R far away from the original point of one transmitter coordinated system, and the results are presented in Figure 5. Worth noting is that in the process of simulation, R is equal to 3 m, while v varies from 00 to 00 mm/s, t 0 is equal to 0 s, t 1 is equal to s, t is equal to s, and v is equal to 89.6p. The static values of the azimuth angle (u) and elevation angle (u) are the standard data which need to be obtained when the velocity of the point is 0 mm/s. However, these measurement values will change if the point is in motion, as illustrated in Figure 5(a). The azimuth angle is in a linear relationship with the velocity of the point in the range of velocity varying from 00 to 00 mm/s. If the velocity is 00 mm/s, the azimuth

6 6 Advances in Mechanical Engineering Figure 5. Simulation results of a single station system model. (a) The azimuth angle (u) in the single station system model is in a linear relationship with the velocity of the point, while its increment direction is in the opposite direction of the velocity. (b) The elevation angle (u) is in a linearly negative correlation relationship with velocity. angle is 67.93, and its value will decrease to when the velocity accelerates to 00 mm/s. As illustrated in Figure 5(b), the relationship between the elevation angle and the velocity of the point presents to be similar with the relationship of azimuth angle (u), which is shown in Figure 5(a). When the velocity of the point varies from 00 to 00 mm/s, the elevation angle decreases proportionally with the minimum value of 0.9 and the maximum value of From equations (6) and (7), we consider that if the linear relationship in the single station system is obtained, the deviation of azimuth angle and elevation angle can be compensated. It is noted that if there are several transmitters in the measuring volume and their measurement values are compensated, the accurate position of the point can be calculated by two arbitrary transmitters through the triangulation principle. Thus, in the following section, we built a double station system model to solve this problem with the help of MATLAB simulation. It is assumed that the distance between transmitters A and B is mm in horizontal direction, mm in the y-axis direction, and mm in vertical direction. The position of A is ( , , and ) and the position of B is ( , , ), and the result is summarized in Figure 6. For the transmitter A, fan 1 sweeps the point after (t 1 t 0 ), and fan sweeps the same point after (t t 0 ), where t 0 is equal to 0s,t 1 is equal to s, t is equal to s, and f 1 is equal to 44.8 Hz. For the transmitter B, fan 1 sweeps the point after (t 4 t 3 ), and fan sweeps the same point after (t 5 t 3 ), where t 3 is equal to s, t 4 is equal to s, t 5 is equal to s, and f is equal to 43.0 Hz. In the measuring volume, the velocity of the measurement point varies from 00 to 00 mm/s. Figure 6. Relationship between the velocity and deviation in the double station system model. When the velocity of the point is 0 mm/s, the deviation is 0 mm without considering the influence of other factors. The deviation is in a linear relationship with the velocity, and it increases if the absolute value of velocity rises in the MATLAB simulation. As illustrated in Figure 6, the deviation in the double station system model is in a linear relationship with the point velocity in the range from 00 to 0 mm/s or from 0 to 00 mm/s. The deviation is considered to be 0 mm when the velocity is 0 mm/s without considering influence of other factors, and the deviation will enhance with the increase in velocity absolute value. Moreover, while the velocity varies from 00 to 00 mm/s, the maximum value of deviation reaches 0.08 mm. We also conducted an experiment to verify the double station system model, as presented in Figure 7. During the experimental process, the min-vector bar which was set on the slider of the linear membrane

7 Zhao et al. 7 Figure 7. Experimental results of dynamic tracking performance of igps. (a) The igps sensor and transmitter used in our experiments. (b) Experimental results of the double station system model. The black dots are the experimental results done by two transmitters and the linear membrane group, and the blue line is the theoretical calculation. It is shown that the relationship between the deviation and the velocity is in a positive correlation, which verifies the established double station system model. group moved with a set velocity varying from 16.5 to 16.5 mm/s, while the two transmitters (Figure 7(a)) were used to measure the position of the min-vector bar. The distance of the slider moving in a second was measured by igps, and the standard distance was given by the control system of linear membrane group. Then we compared the distance with the standard value and obtained the deviation. The experimental results proved a positive correlation between the deviation and the min-vector bar s velocity as shown in Figure 7(b), which was in accordance with the double station system model that we established in section Mathematics model. Therefore, if the parameters, including the object s velocity, time signals, and rotating frequencies, are known, the deviation of each position in the igps measuring volume can be theoretically compensated. Conclusion As an alternative solution for describing dynamic tracking performance of igps, we propose a method to compensate the measurement results of igps. We first built a set of mathematical models including single and double station system models. It has been verified that the relationship between measurement angles and velocity of the object is in linear correlation when the velocity varies from 00 to 00 mm/s in the single station system model with the help of MATLAB simulation. Based on the single station system model, the measurement values of transmitters in the measuring volume can be compensated. Then we chose two random transmitters to calculate the coordinate deviation of the object in double station system model. And we proved that the deviation is in a linear relationship with the velocity of the object in the range from 00 to 0 mm/s or from 0 to 00 mm/s. Therefore, the coordinate of each point in the measuring volume can be theoretically compensated based on the double station system model during the measurement process of igps. We also performed an experiment to evaluate the reasonability of our theoretical study, and the results indicated a positive correlation between the deviation and the velocity, which provided additional support to the mathematical model that we built. We consider that this work we reported would improve the measurement accuracy of igps for tracking dynamic objects and shed light on the industrial application of igps in automation assembly. Declaration of conflicting interests The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article. Funding This work was supported by the Fund of National Engineering and Research Center for Commercial Aircraft Manufacturing (SAMC13-JS-15-07). References 1. Chen SB, Ye Z and Fang G. Intelligentized technologies for welding manufacturing. In: Kiet Tieu A, Zhu H and Zhu Q (eds) Materials science forum (Vols ). Switzerland: Trans Tech Publications Ltd, 013, pp Zhou J. Digitalization and intelligentization of manufacturing industry. Adv Manuf 013; 1: Mei Z and Maropoulos PG. Review of the application of flexible, measurement-assisted assembly technology in aircraft manufacturing. Proc IMechE, Part B: J Engineering Manufacture. Epub ahead of print 7 February 014. DOI: /

8 8 Advances in Mechanical Engineering 4. Liu Z, Xie Y, Xu J, et al. Laser tracker based robotic assembly system for large scale peg-hole parts. In: 014 IEEE 4th annual international conference on cyber technology in automation, control, and intelligent systems (CYBER), Hong Kong, 4 7 June 014, pp New York: IEEE. 5. Zhang F, Qu X and Shi G. A dimensional measurement method based on a laser tracker for a crane assembly. Laser Eng 013; 6: Norman AR, Scho nberg A, Gorlach IA, et al. Validation of igps as an external measurement system for cooperative robot positioning. Int J Adv Manuf Tech 013; 64: Muelaner JE, Wang Z, Martin O, et al. Verification of the indoor GPS system, by comparison with calibrated coordinates and by angular reference. J Intell Manuf 01; 3: Sharke P. Measuring across space and time. Mech Eng 003; 15: Wang Z, Mastrogiacomo L, Franceschini F, et al. Experimental comparison of dynamic tracking performance of igps and laser tracker. Int J Adv Manuf Tech 011; 56: Maisano DA, Jamshidi J, Franceschini F, et al. A comparison of two distributed large-volume measurement systems: the mobile spatial co-ordinate measuring system and the indoor global positioning system. Proc IMechE, Part B: J Engineering Manufacture 009; 3: Muelaner JE, Wang Z, Jamshidi J, et al. Study of the uncertainty of angle measurement for a rotary-laser automatic theodolite (R-LAT). Proc IMechE, Part B: J Engineering Manufacture 009; 3: Schmitt R, Nisch S, Demeester F, et al. Performance evaluation of igps for industrial applications. In: 010 international conference on indoor positioning and indoor navigation (IPIN), Zu rich, September 010, pp.1 8. New York: IEEE. 13. Depenthal C. Path tracking with IGPS. In: 010 international conference on indoor positioning and indoor navigation (IPIN), Zürich, September 010, pp.1 6. New York: IEEE.