# An Alignment Method of Human-Robot Collaboration Based on the Six-Dimensional Force/Torque Dynamic Measurement for Large-Scale Components

Authors: Ke Wen; Fuzhou Du; Jiabo Zhang; JiZhi Yang
Journal: Journal of Robotics
 (2018)

---

## Abstract

A Stewart parallel robot (SPR) is a promising choice for alignment or assembly of components that are large or heavy. This paper presents a method for human-robot collaboration, for positioning and orientation of large components. Use of interactive force measurements is important for human-robot collaboration. It is based on six-dimensional force/torque (F/T) measurements. First, the six-dimensional F/T data are calculated based on the six-actuator SPR geometry and screw theory. Second, the effects of gravity forces (dynamic gravity compensation) are considered, and a method to offset their effects is explained. Third, force estimation experiments were performed using an S-type force sensor and known applied test forces. Finally, the F/T-driven feedback was tested for the alignment of a large-scale component. The experimental results show that the calculated six-dimensional F/T can accurately track the force applied to a large and/or heavy component by a human worker. It can also accurately predict the F/T required to compensate for inertial forces and components’ weight. Thus, the alignment method of human-robot collaboration based on the six-dimensional force/torque dynamic measurements for large-scale components is correct and effective.


---

## Body

## 1. Introduction

With the constant increase in the difficulty of assembly and the numbers of large-scale satellite cabins, rocket cabins, ship cabins, and other large and/or heavy assemblies, the assembly requirements for large-size and heavy-weight components present obstacles. These can include safety risks for workers. Currently, assembly methods based on manual operations have many problems such as physical difficulty for workers, high collision risks, inadequate precision, and low productivity.Modern industrial practice is moving towards intelligence, diversification, and personalization. Human-robot collaborations have become hot spots in the field of robotics. Such robots allow people to enter their work spaces and collaborate to complete operations. Robots and humans no longer need to work separately. This type of collaborative process not only takes advantage of the large load capacity, high stability, and high-precision position and orientation (P&O) capabilities of robots but also complements the capabilities of people for on-site observation and task flexibility [1, 2]. The team of Mitsantisuk introduced a new force sensor less haptic human-robot collaboration system based on 3D parallel mechanism delta robot, and the controller of delta robot was designed with considering of gravity compensation [3]. The team of Geravand introduced two examples of reaction schemes for collaboration, with the user pushing/pulling the robot at any point of its structure (e.g., for manual guidance) or with a compliant-like robot behaviour in response to forces applied by the human [4].An industrial robot’s characteristics of high load capacity, high-precision P&O control, and good stability mean it can offer stable and accurate positioning of heavy components. In recent years, parallel robots have been widely used for P&O of the large components and assemblies, due to these advantages as well as high stiffness and rapid motion [5, 6]. Stewart parallel robots (SPRs) are the most frequently used type of parallel robot structure. They are suitable for machining and manufacturing [7, 8], surgical operations [9], flight simulators [10], flexible and precise assembly of aircraft sections [11], P&O of spacecraft on the ground [12], and low-impact alignments [13]. An SPR is composed of a moving platform and a base platform, which are connected by six extensible links with spherical/universal joints. In its operational range, 6-Degrees-of-Freedom (DOF) motion of the moving platform can be achieved by the motions of the six limbs as a group [14, 15]. The SPR’s limbs’ position and their force measurements are also very important [16]. In 2016, the team of Du from Beihang University presented an analytical algorithm of six-dimensional F/T by measuring limbs’ position and their force based on a SPR [17]. It laid the foundation for the research of this paper.In the process of human-robot collaboration, methods of compensating for the weight of an unknown load are a research hot spot. In general, the mass and center of gravity of the load are not known in advance, so it is necessary to estimate them online. In 2015, the team of Liu from South China University of Technology analyzed the center of gravity of the load according to the relationship between force and torque and then used this calculation to compensate for the effect of load center of gravity on the required screw torque [18]. In 2016, Mitsantisuk of Thai Agricultural University proposed a solution for human-robot collaboration using a delta parallel robot. A disturbance observer was used at the drive joint of the delta robot, to assist in judging the load characteristics by observing additional torque at each drive joint beyond the torque required for motion of the unloaded robot. They also compensated for gravity [3]. In 2017, Boisclair of Laval University in Canada studied a method combining of gravity compensation and the six-dimensional force/torque (F/T) sensor with a cylindrical Halbach array. The cylindrical Halbach array was installed at the joint of an industrial robot and generated joint torque through a magnetic field to compensate the additional torque of the load precisely [19]. In 2017, using six-dimensional F/T sensors and a set of calibration and calculation methods, the team of Zhang from Beijing Institute of Spacecraft Environment Engineering eliminated the effect of sensor zero and load gravity on force perception and accurately obtained the external force and torque applied at the load end of the robot [20, 21].As mentioned above, human-robot collaboration is becoming more and more important in the assembly of large and heavy components, and SPR’s are widely used in this type of process. In this paper, an SPR design is described which can adjust the P&O of the component and dynamically measure interactive forces. Then, the dynamic gravity deviations are considered and accurately compensated for. Lastly, the large-scale alignment experiments were carried out, using an SPR incorporating the six-dimensional F/T dynamic measurements.


## 2. Large-Scale Alignment System Using an SPR and a Human

The large-scale alignment system with a Stewart parallel robot (SPR) for large components is shown in Figure1. It includes an SPR, fixed platform, control cabinet, assembly fixtures, movable component, datum component, and an operator. The SPR is designed based on a 6-SPS (six limbs, each with an actuator, spherical hinge-prismatic pair-spherical hinge) Stewart parallel robot and six force sensors, which can adjust the P&O of the moving platform and dynamically measure six-dimensional F/T from the interactive force. It contains a moving platform, stationary platform (base platform), force sensors, AC servo motors, and branched chains. One of the SPR’s six force sensors is placed in each limb to measure the limbs’ actuation force to calculate the six-dimensional F/T data.Figure 1
The large-scale alignment system with an SPR for large components.The principles for positioning and aligning components are described in this paper. They are based on the measured six-dimensional F/T data and the input forces from a human operator on a component. These principles are shown by the flowchart in Figure2.Figure 2
The human-SPR component alignment flowchart based on the measured six-dimensional F/T data.


## 3. Six-Dimensional Force/Torque Dynamic Measurement of the SPR

### 3.1. Method of Measuring Six-Dimensional F/T Based on the SPR

As presented in Figure3, the SPR consists of a moving platform and a stationary platform, which are connected to each other with six limbs, adjustable in length through sliding joints. In the operation range, the 6-DOF motion of the moving platform could be achieved by the motions of the six limbs. Force sensors are placed in each limb to measure the limbs’ actuation forcefi. Encoders are used to measure the limbs’ lengthli (i = 1, …, 6). The Cartesian coordinate system ofo0-x0y0z0 is located in the center of the top surface of the base platform, while the Cartesian coordinate system ofo1-x1y1z1 is located in the center of the bottom surface of the moving platform. The centers of the spherical joints are denoted asAi andBi.[FsMs]T=[Fx,Fy,Fz,Mx,My,Mz]T is expressed the external load which is decomposed ino1-x1y1z1.Ai is the coordinate ino1-x1y1z1. B’i is the coordinate ino1-x1y1z1 andBi is the coordinate ino0-x0y0z0. The two vectors are different descriptions of the same point; hence they can be related using a linear transformation as follows:(1)Bi′=R3×3-1·Bi3×1-M3×1where R3×3 represents a rotation matrix and M3×3 represents a translation matrix. Once the distance between Ai and Bi (limb lengthli) is set, the P&O parameters {xyzαβγ} betweeno1-x1y1z1 ando0-x0y0z0 could be solved by the newton iteration method. Among the P&O parameters, {xyz} is the position vector ofo1-x1y1z1, with respect too0-x0y0z0, and {αβγ} is the orientation angle ofo1-x1y1z1, with respect too0-x0y0z0. Thus, R3×3 and M3×3 can be expressed as(2)R=cos⁡αcos⁡βcos⁡αsin⁡βsin⁡γ-sin⁡αcos⁡γcos⁡αsin⁡βcos⁡γ+sin⁡αsin⁡γsin⁡αcos⁡βsin⁡αsin⁡βsin⁡γ+cos⁡αcos⁡βsin⁡αsin⁡βcos⁡γ-cos⁡αsin⁡γ-sin⁡βcos⁡βsin⁡γcos⁡βcos⁡γ(3)M=xyzTFigure 3
The theoretical model for a six-dimensional F/T based on the SPR.The force equilibrium equation could be defined in theo1-x1y1z1 using the screw theory as(4)Fs+∈Ms=∑i=16fi$iwhere $i is the unit screw along the ith leg and could be obtained by $i=Si,S0iT.Si is the unit vector offi with respect too1-x1y1z1.S0i is the torque vector ofSi with respect too1-x1y1z1. According to (1),Si andS0i could be calculated as(5)Si=Ai-BiAi-BiS0i=Ai×SiEquation (4) can be rewritten in the form of matrix equation as (6)F=Gfwhere(7)F=FsMsT=FxFyFzMxMyMzT(8)f=f1f2f3f4f5f6T(9)G=S1S2S3S4S5S6S01S02S03S04S05S06Hence, the external load[FsMs]T can be calculated by (6).


### 3.2. Dynamic Gravity Compensation

During the assembly process, the moving platform, assembly fixtures, and components are relatively heavy and bulky, so their barycenter and gravity deviations, which are caused by manufacturing errors and installation errors, will lead to wrong calculation results of the six-dimensional F/T. The dynamic gravity compensation is needed.The influential factors of the calculation results, which are the barycenter and gravity of the moving platform, assembly fixtures, and other components, cannot be ignored. This paper considers them as a rigid system. Equation (6) can be thus rewritten as(10)F+SGS0GW=GfwhereW is the dimensionless value of the gravity, so it is not a vector.SG is the gravity unit vector of the said rigid system, so it is a 3-column vector.S0G is the torque vector ofSG with respect too1-x1y1z1, soS0G is also a 3-column vector.When the external load is zero, the six-dimensional F/T is caused by the gravity of the rigid system (as in Figure4). The coordinateC=[x y z]T indicates the barycenter of the rigid system ino1-x1y1z1. According to the least square principle,C andW can be solved by the six-dimensional F/T under three different P&O sets. As an example, four measurements are performed here:(11)Mx1Mx2Mx3Mx4My1My2My3My4Mz1Mz2Mz3Mz4=0-zyz0-x-yx0Fx1Fx2Fx3Fx4Fy1Fy2Fy3Fy4Fz1Fz2Fz3Fz4Figure 4
Schematic diagram of gravity of the rigid system ino1-x1y1z1.The resolving process ofC andW from (11) is similar to that of the generalized inverse matrix of(12)Fx1Fx2Fx3Fx4Fy1Fy2Fy3Fy4Fz1Fz2Fz3Fz4.When P&O changes,C andWof the rigid system do not differ ino1-x1y1z1, and the direction S=00-1T of gravity does not vary ino0-x0y0z0, meaning(13)SG=R-1·SR-1·SS0G=C×SGwhich could serve for the solution of (10).For the preparation of a six-dimensional F/T measurement, experiments without external loads were carried out first, and the measured F/Ts could be used for the calculations ofC andW, using (11), after whichC is substituted into (13) for the vector[SGS0G]T. Then, for an arbitrary external load, the six-dimensional F/T ino1-x1y1z1 could be obtained using (10) and the gravity of the rigid system can be dynamically compensated.


## 3.1. Method of Measuring Six-Dimensional F/T Based on the SPR

As presented in Figure3, the SPR consists of a moving platform and a stationary platform, which are connected to each other with six limbs, adjustable in length through sliding joints. In the operation range, the 6-DOF motion of the moving platform could be achieved by the motions of the six limbs. Force sensors are placed in each limb to measure the limbs’ actuation forcefi. Encoders are used to measure the limbs’ lengthli (i = 1, …, 6). The Cartesian coordinate system ofo0-x0y0z0 is located in the center of the top surface of the base platform, while the Cartesian coordinate system ofo1-x1y1z1 is located in the center of the bottom surface of the moving platform. The centers of the spherical joints are denoted asAi andBi.[FsMs]T=[Fx,Fy,Fz,Mx,My,Mz]T is expressed the external load which is decomposed ino1-x1y1z1.Ai is the coordinate ino1-x1y1z1. B’i is the coordinate ino1-x1y1z1 andBi is the coordinate ino0-x0y0z0. The two vectors are different descriptions of the same point; hence they can be related using a linear transformation as follows:(1)Bi′=R3×3-1·Bi3×1-M3×1where R3×3 represents a rotation matrix and M3×3 represents a translation matrix. Once the distance between Ai and Bi (limb lengthli) is set, the P&O parameters {xyzαβγ} betweeno1-x1y1z1 ando0-x0y0z0 could be solved by the newton iteration method. Among the P&O parameters, {xyz} is the position vector ofo1-x1y1z1, with respect too0-x0y0z0, and {αβγ} is the orientation angle ofo1-x1y1z1, with respect too0-x0y0z0. Thus, R3×3 and M3×3 can be expressed as(2)R=cos⁡αcos⁡βcos⁡αsin⁡βsin⁡γ-sin⁡αcos⁡γcos⁡αsin⁡βcos⁡γ+sin⁡αsin⁡γsin⁡αcos⁡βsin⁡αsin⁡βsin⁡γ+cos⁡αcos⁡βsin⁡αsin⁡βcos⁡γ-cos⁡αsin⁡γ-sin⁡βcos⁡βsin⁡γcos⁡βcos⁡γ(3)M=xyzTFigure 3
The theoretical model for a six-dimensional F/T based on the SPR.The force equilibrium equation could be defined in theo1-x1y1z1 using the screw theory as(4)Fs+∈Ms=∑i=16fi$iwhere $i is the unit screw along the ith leg and could be obtained by $i=Si,S0iT.Si is the unit vector offi with respect too1-x1y1z1.S0i is the torque vector ofSi with respect too1-x1y1z1. According to (1),Si andS0i could be calculated as(5)Si=Ai-BiAi-BiS0i=Ai×SiEquation (4) can be rewritten in the form of matrix equation as (6)F=Gfwhere(7)F=FsMsT=FxFyFzMxMyMzT(8)f=f1f2f3f4f5f6T(9)G=S1S2S3S4S5S6S01S02S03S04S05S06Hence, the external load[FsMs]T can be calculated by (6).


## 3.2. Dynamic Gravity Compensation

During the assembly process, the moving platform, assembly fixtures, and components are relatively heavy and bulky, so their barycenter and gravity deviations, which are caused by manufacturing errors and installation errors, will lead to wrong calculation results of the six-dimensional F/T. The dynamic gravity compensation is needed.The influential factors of the calculation results, which are the barycenter and gravity of the moving platform, assembly fixtures, and other components, cannot be ignored. This paper considers them as a rigid system. Equation (6) can be thus rewritten as(10)F+SGS0GW=GfwhereW is the dimensionless value of the gravity, so it is not a vector.SG is the gravity unit vector of the said rigid system, so it is a 3-column vector.S0G is the torque vector ofSG with respect too1-x1y1z1, soS0G is also a 3-column vector.When the external load is zero, the six-dimensional F/T is caused by the gravity of the rigid system (as in Figure4). The coordinateC=[x y z]T indicates the barycenter of the rigid system ino1-x1y1z1. According to the least square principle,C andW can be solved by the six-dimensional F/T under three different P&O sets. As an example, four measurements are performed here:(11)Mx1Mx2Mx3Mx4My1My2My3My4Mz1Mz2Mz3Mz4=0-zyz0-x-yx0Fx1Fx2Fx3Fx4Fy1Fy2Fy3Fy4Fz1Fz2Fz3Fz4Figure 4
Schematic diagram of gravity of the rigid system ino1-x1y1z1.The resolving process ofC andW from (11) is similar to that of the generalized inverse matrix of(12)Fx1Fx2Fx3Fx4Fy1Fy2Fy3Fy4Fz1Fz2Fz3Fz4.When P&O changes,C andWof the rigid system do not differ ino1-x1y1z1, and the direction S=00-1T of gravity does not vary ino0-x0y0z0, meaning(13)SG=R-1·SR-1·SS0G=C×SGwhich could serve for the solution of (10).For the preparation of a six-dimensional F/T measurement, experiments without external loads were carried out first, and the measured F/Ts could be used for the calculations ofC andW, using (11), after whichC is substituted into (13) for the vector[SGS0G]T. Then, for an arbitrary external load, the six-dimensional F/T ino1-x1y1z1 could be obtained using (10) and the gravity of the rigid system can be dynamically compensated.


## 4. Experimental Results and Discussion

The component alignment system using an SPR to align large components is designed as shown in Figure5. It includes the SPR, fixed platform, control cabinet, assembly fixtures, movable component, datum component, and operator. The parameters of the SPR are listed in Table 1.Table 1
The parameters of the SPR.
Parameter name Value Moving range along x direction ±50 mm Moving range along y direction ±50 mm Moving range along z direction ±50 mm Rotation range along x direction ±5° Rotation range along y direction ±5° Rotation range along z direction ±5° Maximum load 800 kg Accuracy 0.05 mm Overall size 1.5m×1.5m×1mFigure 5
The large-scale alignment system with a SPR supporting a large component.


### 4.1. Measurement of Known Applied Force

An S-type force sensor (with an accuracy of ±0.03 N) is used to measure the magnitude of applied force in real time. The operator applies the force to the SP or component to be moved, and the force is transmitted through the SPR to the S-type force sensor. The value of the six-dimensional F/T data can then be calculated and compared with the measurement value of the S-type force sensor. Datapoints were collected every 0.05 s during the experiment. The results are presented in Figure6.Figure 6
Comparison between the measured forces and calculated forces.To measure an applied force, three steps were followed in this experiment. The first step is to compensate for the weight of the SPR and component, and the dynamic forces due to motion. When no external forces are applied to the SPR, any four groups of P&O parameters are selected as inputs to the SPR. Then, the four groups of the six-dimensional F/T can be calculated. The weight and center of gravity of the rigid body (moving platform, assembly fixture, and movable component) are inputs to the dynamic gravity compensation method. The weight is 1969.1017 N. The barycentric coordinate of the rigid body is-36.7749,-6.2296,-208.3247T ino1-x1y1z1. Second, the mean values of the six-dimensional F/T data are calculated. Due to measurement errors that caused the six-dimensional F/T data to fluctuate, the mean values of 100 groups of the six-dimensional T/F data are taken as the reference zero values. Finally, the six-dimensional F/T data are calculated. Figure 6 compares actual force data, as applied by the operator, to the input (interactive) force calculated by the system. These F/T data allows the system to estimate the interactive force between a person and the SPR or load component.The results compared in Figure6 show that the SPR system can accurately track the actual (measured) interactive force. The measured force and the calculated force are in good agreement, with an average error magnitude of 14.70 N (the difference between measured force and calculated force is taken as the absolute value, and then the average error magnitude is taken). Thus, these experimental results verify that this method of using six-dimensional SPR F/T data can track the operator’s actual force inputs.


### 4.2. F/T-Driven Alignment of Large Components

The operator applies the force to the SPR, and the SPR adjusts the P&O of large component to follow the intentions of the operator. Thus, the F/T feedback-driven alignment of large components can be implemented as in Figure2. It is described as follows. First, some parameters are initialized, and the threshold values of six-dimensional F/T are obtained by experimental optimization. The applied threshold value of the force magnitude was 60 N, and the applied threshold value of the torque magnitude was 25 N∗m. Second, the operator visually judged the direction of the desired motion, and the external force was applied to the SPR by the operator. Third, the six-dimensional F/T tensor was calculated in real time by the dynamic force and weight compensation. This interactive force is taken as an input to the SPR. Fourth, the process of the F/T-driven alignment of large components was then executed though the direction and magnitude of interactive force and torque inputs. At last, the SPR adjusts the P&O of large component to follow the intentions of the operator. The experiment and the result are presented in Figure 7.Figure 7
P&O and F/T curves over time.In Figure7(a), when the calculated force Fx is greater than the threshold value 60N, the SPR moved along the -x directions, until it stops at about -40 mm. At this location, the alignment process was successfully completed. Small changes in other degrees of freedom indicate that the components have high assembly accuracy and the adjustment range is small aside from the -x movement (as shown in Figure 7(b)). After the alignment process is completed, the maximum force value in the six-dimensional F/T is 35 N, and the maximum torque value is 12 N∗m (as shown in Figure 7(a)). These experimental results demonstrate that the calculation method of the six-dimensional F/T based on the SPR is correct, and the F/T-driven alignment of large components is effective.


## 4.1. Measurement of Known Applied Force

An S-type force sensor (with an accuracy of ±0.03 N) is used to measure the magnitude of applied force in real time. The operator applies the force to the SP or component to be moved, and the force is transmitted through the SPR to the S-type force sensor. The value of the six-dimensional F/T data can then be calculated and compared with the measurement value of the S-type force sensor. Datapoints were collected every 0.05 s during the experiment. The results are presented in Figure6.Figure 6
Comparison between the measured forces and calculated forces.To measure an applied force, three steps were followed in this experiment. The first step is to compensate for the weight of the SPR and component, and the dynamic forces due to motion. When no external forces are applied to the SPR, any four groups of P&O parameters are selected as inputs to the SPR. Then, the four groups of the six-dimensional F/T can be calculated. The weight and center of gravity of the rigid body (moving platform, assembly fixture, and movable component) are inputs to the dynamic gravity compensation method. The weight is 1969.1017 N. The barycentric coordinate of the rigid body is-36.7749,-6.2296,-208.3247T ino1-x1y1z1. Second, the mean values of the six-dimensional F/T data are calculated. Due to measurement errors that caused the six-dimensional F/T data to fluctuate, the mean values of 100 groups of the six-dimensional T/F data are taken as the reference zero values. Finally, the six-dimensional F/T data are calculated. Figure 6 compares actual force data, as applied by the operator, to the input (interactive) force calculated by the system. These F/T data allows the system to estimate the interactive force between a person and the SPR or load component.The results compared in Figure6 show that the SPR system can accurately track the actual (measured) interactive force. The measured force and the calculated force are in good agreement, with an average error magnitude of 14.70 N (the difference between measured force and calculated force is taken as the absolute value, and then the average error magnitude is taken). Thus, these experimental results verify that this method of using six-dimensional SPR F/T data can track the operator’s actual force inputs.


## 4.2. F/T-Driven Alignment of Large Components

The operator applies the force to the SPR, and the SPR adjusts the P&O of large component to follow the intentions of the operator. Thus, the F/T feedback-driven alignment of large components can be implemented as in Figure2. It is described as follows. First, some parameters are initialized, and the threshold values of six-dimensional F/T are obtained by experimental optimization. The applied threshold value of the force magnitude was 60 N, and the applied threshold value of the torque magnitude was 25 N∗m. Second, the operator visually judged the direction of the desired motion, and the external force was applied to the SPR by the operator. Third, the six-dimensional F/T tensor was calculated in real time by the dynamic force and weight compensation. This interactive force is taken as an input to the SPR. Fourth, the process of the F/T-driven alignment of large components was then executed though the direction and magnitude of interactive force and torque inputs. At last, the SPR adjusts the P&O of large component to follow the intentions of the operator. The experiment and the result are presented in Figure 7.Figure 7
P&O and F/T curves over time.In Figure7(a), when the calculated force Fx is greater than the threshold value 60N, the SPR moved along the -x directions, until it stops at about -40 mm. At this location, the alignment process was successfully completed. Small changes in other degrees of freedom indicate that the components have high assembly accuracy and the adjustment range is small aside from the -x movement (as shown in Figure 7(b)). After the alignment process is completed, the maximum force value in the six-dimensional F/T is 35 N, and the maximum torque value is 12 N∗m (as shown in Figure 7(a)). These experimental results demonstrate that the calculation method of the six-dimensional F/T based on the SPR is correct, and the F/T-driven alignment of large components is effective.


## 5. Conclusions

The SPR can adjust the P&O of the component and dynamically measure forces from interaction between a robot and a human. The experimental results show that the forces due to weight can be reliably compensated for. The contributions of the paper are summarized as follows:( 1 ) A method of measuring six-dimensional F/T, based on the SPR characteristics and screw theory, is proposed for the accurate estimation of the external forces during alignment of large-scale components. Forces due to gravity are taken into consideration, and a precise weight compensation method is provided.( 2 ) Known applied force measurement experiments are performed on the SPR using an S-type force sensor. The six-dimensional F/T can accurately track the change of the measurement force and the average error is about 14.70 N. The experiment of F/T-driven alignment for large component was successfully completed.( 3 ) Interactive force measurements are important for human-robot collaboration. The proposed method fulfills the requirement for accurate force estimation. It provides satisfactory compensation for forces due to gravity, in order to accurately separate forces due to the interaction of a human with a robot. The measurements from the six-dimensional F/T based on the SPR could be further improved with more accurate measurements of the forces and the limb lengths’ variation with motion.



---

*Source: 2914871-2018-09-02.xml*