# Error Modeling and Compensation of Circular Motion on a New Circumferential Drilling System

Authors: Qiang Fang; Yu-Chao Li; Shao-Hua Fei; Yong-Wei Han; Ying-Lin Ke
Journal: Mathematical Problems in Engineering
 (2015)

---

## Abstract

A new flexible circumferential drilling system is proposed to drill on the fuselage docking area. To analyze the influence of the circular motion error to the drilling accuracy, the nominal forward kinematic model is derived using Denavit-Hartenberg (D-H) method and this model is further developed to model the kinematic errors caused by circular positioning error and synchronization error using homogeneous transformation matrices (HTM). A laser tracker is utilized to measure the circular motion error of the two measurement points at both sides. A circular motion compensation experiment is implemented according to the calculated positioning error and synchronization error. Experimental results show that the positioning error and synchronization error were reduced by 65.0% and 58.8%, respectively, due to the adopted compensation, and therefore the circular motion accuracy is substantially improved. Finally, position errors of the two measurement points are analyzed to have little influence on the measurement result and the validity of the proposed compensation method is proved.


---

## Body

## 1. Introduction

In large aircraft fuselage assembly process, fuselage sections are supported by three-coordination-system numerical positioners [1, 2] to align with each other and there are lots of fastener holes to be drilled on the docking area with wide range and high difficulty to join the sections together. Thus automated precision drilling system becomes an essential integral part of aircraft numerical assembly technology. However, traditional approaches for automated drilling of the fuselage sections which employ large truss equipment require high investment and large space during implementation [3]. While this kind of technology remains effective for most applications, the aerospace industry is in demand of developing a low-cost, highly flexible, and lightweight special hole drilling equipment for fuselage sections. The aerospace manufacturing industry has successfully developed such equipment. Boeing Company developed the first Flex Track drilling system; then Electroimpact Company developed a variety of different Flex Track drilling systems to satisfy all kinds of applications, achieving hole accuracy of 0.025 mm. The newly developed Flex Track drilling system in recent years has been successfully applied to the assembly of Boeing’s B777 and B787 fuselage sections and Airbus’s A380 wing box drilling system [4–6]. The domestic institutes have caught tightly the international steps on the theory research and equipment development of the Flex Track drilling system; among them AVIC Beijing Aeronautical Manufacturing Technology Research Institute developed BAA300 system; this system can ensure the fine accuracy and quality of the holes [7] but is limited to theoretical research and prototype verification; the system has not been put into the practical application. The main advantage of the Flex Track drilling system over traditional drilling system is that Flex Track system uses vacuum cups to install the tracks of the drilling unit directly on the fuselage sections instead of requiring large-scale machine tools like infrastructure. Since the tracks are flexible, they are able to adapt to different machined work pieces, but the precision of the system becomes vulnerable due to the geometry properties of the surface, the structural rigidity, and other factors. Besides, tracks of the Flex Track drilling system need to be installed repeatedly and the geometric error of the splice region increases the difficulty of the implementation.Considering the demand of some fuselage sections splice region hole drilling, based on the features of the existing Flex Track drilling system, Zhejiang University developed a new type of circumferential automated drilling system which has been successfully applied to the production. The circumferential drilling system includes double tracks as the basic system support. The inner circular rail consisting of ten segments is installed on the fuselage surface by adjustable contact pads, and the outer arc guide is used as the support of the drilling actuator. The drilling actuator can move along the outer arc guide and drilling holes within a certain drilling station while the arc guide can move along the circular rails to switch a total of twelve drilling stations. Thus, the circumferential drilling system can complete the whole splice region drilling task through one installation. Closed structure of the circular rails and the double rigid tracks of the system effectively guarantee the system stiffness and constrain fuselage deformation when machined. The adjustable contact pads make the circular rails adapt to large surface curvature, nonround feature, and nonequivalence segment of the drilling surface [7]. The new circumferential automated drilling system introduced in this paper was fully upgraded over the one in literature [8, 9] with improved functions and reliability, achieving automatic drilling stations switching and dual-motor synchronous drive for the circular movement axis.Drilling actuator’s position of the circumferential drilling system with respect to the fuselage section is calibrated by the reference holes’ position measured by vision-based measurement unit on the arc-based drilling unit. The positions of the holes to be drilled can be planned under the drilling system coordinate system like a 2D flat problem. While the normality sensor guarantees the normal accuracy of the drilled holes, the position accuracy mainly depends on the accuracy of the circular movement axis and longitudinal movement axis. The circular movement axis differs from traditional revolute axis that it is driven by dual-motor through two gear and arc rack transmission pairs. Due to large motion radius and low rail stiffness, the circular movement axis only engages semiclosed loop control. Therefore, in order to guarantee high ultimate drilling accuracy, it is necessary to improve the circular motion accuracy of the drilling actuator on the arc guide. Electroimpact Company adds a nondrive axis equipped with an encoder to measure the angular difference between the movement along the vacuum track and the longitudinal movement axis in the Flex Track drilling system for B777 fuselage sections assembly in order to eliminate the synchronization error and realize the docking drilling on the round table facet [10–12]. In this solution, the accuracy of effector is determined by the rail on one side and effector’s torsional rigidity in the direction of drilling axis is comprised. Thus, it is not suitable for the circumferential drilling system.Error compensation technology [13–15] is applied to CNC machine tools first and can be utilized to improve the circular motion accuracy of the circumferential drilling system. Three main error sources [16, 17] of CNC machine tools are geometric and kinematic errors, thermal errors, and cutting force induced errors, among which geometric and kinematic errors typically account for more than 40% with good repeatability. The circumferential drilling system has a lower movement velocity compared to CNC machine tools, and other movement axes remain steady while drilling holes so that they produce lower thermal errors and cutting forces errors. Therefore, compensation for geometric errors can effectively improve the accuracy of movement of the drilling actuator. Error compensation technique includes error modeling, error measurement, and error compensation. As the theoretical foundation of error measurement and compensation, establishing accurate comprehensive error model can effectively evaluate the influence of each error component, the normal ways of modeling are geometric modeling, Homogeneous coordinate transformation matrix, and modeling based on multibody system theory. The measurement system with laser tracker is of the advantages of fast and high accuracy, and common equipment for error detection fits to this system to measure the points in a large space. Considering the error after measuring and identifying, the opposite offset has been added into the software in advance, in order to do the compensation of the system.In Section2 the structure of the proposed circumferential drilling system is analyzed. The nominal forward kinematic model is derived and this model is further developed to model the kinematic errors caused by circular positioning error and synchronization error using the homogeneous coordinate transformation method in Section 3. The compensation method is introduced according to Section 4. To verify the compensation method, the compensation experiment has been performed with the utilization of a laser tracker and the experiment result is discussed in Section 5.


## 2. Structure of the Circumferential Drilling System

The circumferential drilling system is newly proposed, highly flexible, and lightweight numerical drilling equipment for the demand of drilling on the fuselage sections splice region based on a comprehensive analysis of the traditional drilling machine and contemporary Flex Track system. The main structure of the circumferential drilling system consists of bottom auxiliary carriage, circular rails, holders, and arc-base drilling unit, equipped with vacuum cleaner, cooling system, and controller cart, as shown in Figure1.Figure 1
Structure of circumferential drilling machine.The bottom auxiliary carriage is the main support for circular rails. The circular rails are installed on the fuselage sections by contact pads, and the arc-base drilling unit is mounted on the circular rails.Arc-base drilling unit is used for drilling in the splice region within an arc’s span, containing six axes of motion, as shown in Figure2. The X-axis along the arc guide is driven by the dual-motor with the constraint of the arc guide. The arc guide occupies range of 53° to ensure finishing the drilling task of the entire circumferential region in twelve drilling stations and keeping space for station switching. The Y-axis moves along the longitudinal direction. The X-axis and Y-axis are used to control the position of the holes. The A-axis offers rotation about the tangent line of the X-axis and the B-axis offers rotation about the Y-axis. The A-axis and B-axis both control the normality of the drilled holes. The two Z-axes move along the drilling direction. The Z1-axis with larger motion range is used to avoid obstacles and make space for tool exchange, while the feeding Z2-axis can guarantee precise countersink depth by controlling the relative position with the pressure foot. The drilling unit is mounted with the normality sensor and the vision-based measurement unit.Figure 2
Arc-base drilling unit.Since the drilling actuator can only drill holes within an arc’s span of 30°, the arc-base drilling unit needs to be fixed on the circular rails at different angle to finish the drilling task. The automated drilling station switching process is shown in Figure3. When the drilling unit drills holes within current station, the arc guide is locked on the circular rails by one group of electromagnetic latches. When the arc-base drilling unit needs to switch to the next station, the drilling actuator moves to specified position and the other group of electromagnetic latches locks the drilling actuator on the circular rails; with the first group of electromagnetic latches released, the arc guide moves along the circular rails driven by the dual-motor to switch drilling station. Therefore, the circumferential drilling system can finish drilling task of the entire splice region automatically.Figure 3
Drilling station switch process.


## 3. Circular Motion Error Modeling

Error modeling is not only an effective way to analyze the causes of error but also the theoretical basis for improving the drilling accuracy through error compensation. The arc-base drilling unit contains six movement axes:X-axis, A-axis, Y-axis, B-axis, Z1-axis, and Z2-axis according to the order of the kinematic chain. Denavit-Hartenberg (D-H) method is adopted to derive the nominal forward kinematic model of the drilling machine [8], as shown in Figure 4.Figure 4
Kinematic structure of drilling unit.There are six joint axes in the kinematic model of the circumferential drilling system, denoted as axisi, i=1,2,…,6, where axis i describes motion between link i-1 and link i and frame {i} is attached rigidly to link i. Under D-H notation, the Z^-axis of frame {i} called Z^i is coincident with the motion direction of axis i and X^i points along ai from axis i to axis i+1, where ai is the length of the common perpendicular between axis i and axis i+1.For convenience, base frame{0} attached to the base of the circumferential drilling system is chosen coincident to frame {1} when axis 1 is in zero position. The origin of the frame {0} is located right between the centers of the two circular rails. Axis 1 is the revolute X-axis and frame {1} is rotated from frame {0} by the joint variable θ1. The pose of frame {1} with respect to frame {0} is(1)T10=Rzθ1=cθ1-sθ100sθ1cθ10000100001,where s(·) and c(·) denote sin(·) and cos(·), respectively.Axis 2 is the rotationalA-axis about the tangent line of the X-axis. Frame {2} can be acquired in three steps:(1)
Translate frame{1} along X^1 by a1 (radius of the X-axis).(2)
Rotate aroundX^1 by −90°.(3)
Rotate aroundZ^1 by joint variable θ2.The pose of frame{2} with respect to frame {1} is(2)T21=Txa1Rx-90°Rzθ2=cθ2-sθ20a10010-sθ2-cθ2000001.Axis 3 is the translationalY-axis along the longitudinal direction. Frame {3} can be acquired in two steps:(1)
Rotate frame{2} around X^2 by 90°.(2)
Translate alongZ^2 by joint variable d3.The pose of frame{3} with respect to frame {2} is(3)T32=Rx90°Tzd3=100000-1-d301000001.Axis 4 is the rotationalB-axis about the Y-axis. Frame {4} can be acquired in three steps:(1)
Translate frame{3} along X^3 by a3 (distance between A-axis and B-axis).(2)
Rotate aroundZ^3 by −90°.(3)
Rotate aroundZ^3 by joint variable θ4.The pose of frame{4} with respect to frame {3} is(4)T43=Txa3Rz-90°Rzθ4sθ4cθ40a3-cθ4-sθ40000100001.Axis 5 and axis 6 are theZ1-axis and the Z2-axis moving along the drilling direction, respectively. Frame {5} can be acquired in three steps:(1)
Rotate frame{4} around Z^4 by 180°.(2)
Rotate aroundX^4 by 90°.(3)
Translate alongZ^4 by joint variable d5.The pose of frame{5} with respect to frame {4} is(5)T54=Rz180°Rx90°Tzd5=-1000001d501000001.And frame{6} is translated from frame {5} along Z^5 by joint variable d6. The pose of frame {6} with respect to frame {5} is(6)T65=Tzd6=10000100001d60001.In order to determine the coordinate of the tool center point with respect to the base frame{0}, frame {7} is attached to the tool with fixed position offset (x7,y7,z7) with respect to frame {6}. The pose of frame {7} with respect to frame {6} is(7)T76=100x7010y7001z70001.Link parameters of the arc-base drilling unit are shown in Table1.Table 1
Link parameters of the arc-base drilling unit.
Linki
a
i
-
1
α
i
-
1
d
i
θ
i
Joint variable
1
0
0°
0
θ1
θ1: ±16°
2
a1 (2985 mm)
−90°
0
θ2
θ2: ±10°
3
0
90°
d3
0°
d3: 450 mm
4
a3 (170 mm)
180°
0
θ3
θ3: ±15°
5
0
90°
d5
180°
d5: 500 mm
6
0
180°
d6
0°
d6: 300 mmThe nominal forward kinematic model of the circumferential drilling system can be derived by multiplying coordinate transformation matrix (1) to (7) as follows:(8)T70=∏i=17Tii-1.When the drilling actuator moves along the arc guide, there are six geometric error components in all. Since circular positioning error and synchronization error can only be compensated by the dual-motor, the other errors are small and can only be reduced by other axes. Hence, the circular dual-drive positioning error and synchronization error are considered to establish comprehensive error model to do compensation on the dual-motor. The parameters of the circular movement axis are shown in Table2, where R is the radius of the arc rack, L is the distance between the two racks, and D is the travel distance of the arc guide.Table 2
Parameters of circular movement axis.
Name
R
L
D
Reduction ratio
Motor type
Gear parameters
Value
2900 mm
942 mm
1600 mm
25
AKM32E-ANC2C-00
d
=
40 mm, m=2Notex1 and x2 as the actual displacements of the drilling actuator on the two racks, respectively. The circular positioning error causes angle error Δγ around Z^1 of frame {1}:(9)Δγ=x1+x2/2-xcmdR.The circular synchronization error causes angle errorΔα around X^1 of frame {1}:(10)Δα=arctan⁡x1-x2L.Therefore, after taking circular positioning error and synchronization error into consideration, transformT10 is replaced by T10′:(11)T10′=T10·RZ^1Δγ·RX^1Δα,where(12)RZ^1Δγ=cΔγ-sΔγ00sΔγcΔγ0000100001(13)RX^1Δα=10000cΔα-sΔα00sΔαcΔα00001.Thus, the forward kinematic model of the circumferential drilling system considering circular motion error is(14)T70′=T10′∏i=27Tii-1.And the tool center positioning error is(15)E=T70′-T70.Therefore, according to the circumferential drilling system forward kinematic model and parameters, given circular positioning error and synchronization error, tool center positioning error can be calculated to evaluate the influence to the ultimate drilling accuracy.


## 4. Error Compensation Methodology

Error compensation refers to the method of artificially creating a new error to neutralize the systematic error of the current movement. The circular movement axis of the circumferential drilling system generates circular paths. Due to small central angle of the path, arc fitting will introduce large errors. Hence, chord lengths of the adjacent points are used to calculate positioning error and synchronization error. Also, the average chord lengthcavg is used as the ideal chord length:(16)cavg=∑i=1Nc1i+∑i=1Nc2i2N,where c1i and c2i are the chord lengths of the adjacent points on the two rails, respectively, and N is the quantity of chords of each rail.Circular positioning errorΔn is calculated by(17)Δn=∑i=1nc1i+c2i2-cavg·n,where n=1,2,…,N is the number of measurement points.Circular synchronization errorξn is calculated by(18)ξn=c1i-c2ii=n.Compensation value is calculated as follows:(19)xcmd1n′=xcmdn·cavg·n∑i=1nc1ixcmd2n′=xcmdn·cavg·n∑i=1nc2i,where xcmd1n′ and xcmd2n′ are the compensation value of the two rails, respectively, and xcmdn is the demand value before compensation. The compensation value between measurement points is calculated by linear interpolation.


## 5. Experiment and Discussion

In order to verify the compensation method, a compensation experiment is perforated on the circumferential drilling system as shown in Figure5. A Leica laser tracker is adopted to measure the circular positioning error and synchronization error when the drilling actuator is moving along the arc guide, as shown in Figure 6.Figure 5
The circumferential drilling machine.Figure 6
Compensation experiment.Here is the experiment procedure:(1)
Select a proper drilling station and pick two measurement points symmetrically on both sides of the drilling actuator near the circular rails.(2)
Within the whole 1600 mm range of theX-axis, measure both sides of the measurement points every 50 mm movement.(3)
Repeat(2) three times and take an average to gain 32 pairs of point coordinates approximate on two parallel arc segments.(4)
Use 31 pairs of chord lengths of the adjacent points to calculate positioning error and synchronization error for compensation value calculation.(5)
Repeat(3) and (4) using compensation value and compare the positioning error and synchronization error with the ones before compensation.Error compensation implementation intends to set the tool or work piece position to generate relative motion in the reverse error direction to eliminate the original error. Software compensation technology will set offset on the original command value in advance which is economic and convenient to implement. The circumferential drilling system movement axes use Kollmorgen motion solutions and the drilling control system is developed with the provided MPI (motion programming interface) in C++/MFC environment. Hence, compensation program can be easily developed.Do compensation experiment with compensation value calculated according to formula (19); the comparison of the positioning error and synchronization error before and after compensation is shown in Figures 7 and 8, respectively. The circular positioning error and synchronization error were reduced by 65.0% and 58.8%, respectively, after compensation. Thus, the circular motion accuracy is largely improved due to the compensation.Figure 7
Positioning error comparison.Figure 8
Synchronous error comparison.In general, only one measurement point is utilized when using a laser tracker to measure position. In order to measure synchronization error, two measurement points are adopted here. It is essential to analyze the influence of the position errors of the two measurement points, as shown in Figure9.Figure 9
Measurement error analysis.Without loss of generality, assume that one of the measurement points is in ideal position, the other measurement point has spatial position errorδ=δx,δy,δz in frame {1}, where δx and δy cause deviation along the radius direction and generate additional synchronization error. δz along the longitude direction changes distance to the plane X1O1Y1 and magnify synchronization error by (L1-δz)/L1, where L1 is the nominal distance between two measurement points. In summary, the measured synchronization error ξn′ and the actual synchronization error ξn have the following relationship:(20)ξn′=∑i=1nc1i+c2iδxcosα-δysinα2R1+L1-δz·ξnL1,where R1 is the nominal motion radius of the measurement points and α is the angle of the measurement points in plane X1O1Y1. Since α is very small and L1≫δz, formula (20) can be simplified to(21)ξn′-ξn=∑i=1nc1i+c2iδx2R1.Supposeδx=0.2mm; according to formula (21), measurement error induced synchronization error is less than 0.1 mm. Therefore, δ has little influence on the measurement result so that the measurement and compensation can effectively improve the circular motion accuracy within a certain precision.


## 6. Conclusion

This paper describes a newly developed circumferential drilling system to satisfy fuselage section splice region drilling demand and derives the nominal forward kinematic model. Comprehensive error model is established by using the homogeneous coordinate transformation considering circular positioning error and synchronization error. The laser tracker is utilized to measure the circular motion error and the compensation experiment is implemented. The experimental result shows that the positioning error and synchronization error were reduced by 65.0% and 58.8%, respectively. At last, position errors of the two measurement points are analyzed to have little influence on the measurement result and the proposed compensation method is effective.



---

*Source: 286796-2015-09-27.xml*