# Development of a Prototype Robotic System for Radiosurgery with Upper Hemispherical Workspace

Authors: Sun Young Noh; Kyungmin Jeong; Yong-chil Seo; Chang-hoi Kim; Jongwon Park; Yoo Rark Choi; Sung Uk Lee; Yeong-Geol Bae; Seungho Kim
Journal: Journal of Healthcare Engineering
 (2017)

---

## Abstract

This paper introduces a specialized robotic system under development for radiosurgery using a small-sized linear accelerator. The robotic system is a 5-DOF manipulator that can be installed above a patient to make an upper hemispherical workspace centered in a target point. In order to determine the optimal lengths of the link, we consider the requirements for the workspace of a linear accelerator for radiosurgery. A more suitable kinematic structure than conventional industrial manipulators is proposed, and the kinematic analysis is also provided. A graphic simulator is implemented and used for dynamic analysis. Based on those results, a prototype manipulator and its control system are under development.


---

## Body

## 1. Introduction

The aim of stereotactic radiosurgery systems such as Linac, Gamma Knife, and CyberKnife is to destroy tumor tissue while preserving adjacent normal tissue using ionizing radiation rather than excision with a blade [1, 2]. Thus, they should be able to accurately focus radiation beams from as many different angles as possible to converge on one target tumor [3, 4].A Gamma Knife is used for treating not only brain tumors but also vascular and functional pathologies. And it typically contains 201 Co-60 radiation sources, each placed in a circular array, and emits gammy ray through a target point in the patient’s brain. Thus, the Gamma Knife has over 200 beam delivery angles [5, 6].Compared with the Gamma Knife, a Linac system uses X-rays generated from a linear accelerator, and the accelerator is mechanically rotated around the patient, in a full or partial circle to change the delivery angle aiming a target point. And the couch where the patient is lying can also be moved in small linear or angular steps. Thus, the combination of the movements of the gantry and of the couch can make more beam delivery angles. But they are confined within an approximately 2-dimensional space [7, 8]. A CyberKnife system also uses a linear accelerator as a radiation source, but relatively small and light-weighed. It uses a 6-DOF industrial robot manipulator to move the X-ray source to aim arbitrary points from different angles [9, 10]. So, CyberKnife systems are recognized as one of the most versatile radiosurgery systems having a broad range of the beam delivery angle in a 3-dimensional space [11, 12].But even though the degree of freedom of the manipulator used in a CyberKnife system is 6, 5-DOF is sufficient to aim an arbitrary 3-dimensional point from an arbitrary delivery angle because the change in the rolling angle of the beam about its axis has no effects [13, 14].And besides, whereas the required radius of workspace for the radiosurgery is about 1 meter and the maximum moving speed of the end-effector is less than 1 m/s, the maximum radius of the workspace in the CyberKnife system is approximately 3 meters and its moving speed is also higher than required. That means that the performance of the joint-driving motors of industrial robot manipulators in the market is usually higher than required, and the driving motors are oversized. It may cause the system to be heavy-weighted and make its cost high. Furthermore, in cases when the robot manipulator is installed on the floor beside the couch, there are limitations for increasing the range of the beam delivery angle.So, this paper introduces a specialized robotic system under development for radiosurgery using a small-sized linear accelerator. The robotic system is a 5-DOF manipulator installed above a patient to make an upper hemispherical workspace centered in a target point.In Section2, the basic structure and the design process for determining the lengths of its links are described. Section 3 provides the kinematic analysis of the manipulator, and Section 4 explains a graphic simulation for dynamic analysis and its results. The manufactured manipulator under development and the structure of control system are shown in Section 5.


## 2. Design of Lengths of Robot Arm Links

As mentioned above section, the manipulator for radiosurgery needs 5-DOF joints in minimum to deliver a linear accelerator, so the lengths of the links should be optimized for its operation. And also, the manipulator workspace takes into consideration the position of the patient and is designed to avoid contact with the patient. This is achieved by creation of a safety zone around the patient and the treatment couch.In Figure1, the center of a tumor tissue to be treated is assumed to be located at the origin of xyz coordinates. The x- and y-axes define a horizontal plane, and the z-axis is vertical.Figure 1
Workspace for radiosurgery.Letxe,ye,ze refer to the position of X-ray source inside the linear accelerator.If SAD (source-to-axis distance) is betweendmin and dmax, and the workspace is limited to the upper hemisphere, then the required workspace where the end point should reach at is
(1)dmin≤xe2+ye2+ze2≤dmax,ze≥0.At any point in the workspace, the accelerator should direct radiation beams to a small region around the origin.In our study, a serial manipulator with revolute joints is considered as shown in Figure2(a).Figure 2
Proposed configuration of robot.
(a)
(b)
(c)The joint 1 (J1) is aligned with the z-axis and defines a central sagittal plane, and by using joints 2 (J2), 3 (J3), and 4 (J4) parallel to each other, the accelerator is directed to the origin at any point in that plane. Both joints 4 and 5 (J5) are used to direct the beam to the neighboring region.The center of joint 2 is shifted above the origin withh to avoid collision with a patient.One extreme constraint with the lengths of the links is the case whenϕ=−π/2,d=dmax where the angle ϕ is defined as shown in Figure 2(a), and the reach of the arm is maximum as shown in Figure 2(b). To access such a posture in its full reach, the lengths of the links satisfy the following equation:
(2)L2+L3=h2+dmax+L42.Another case is whenϕ=π/2,d=dmin→ϕ=−π,d=dmin. In that case, the following equation satisfies the posture:
(3)h+L2−L3−L4=dmin.The link’s lengthsL2,L3,L4 satisfy (2) and (3), and the required workspace of the arm is obtained.Whendmin,dmax,L4 are given, L2 and L3 can be determined as follows:
(4)L2=12h2+dmax+L42−h+L4+dmin,L3=12h2+dmax+L42+h−L4−dmin.Table1 shows L2 and L3 determined from (4) when dmin=0.6m,dmax=1.0m,L4=0.1m optimal length according to H.Table 1
Optimal length.
H
L
2
L
3
0.5
0.8
0.5
0.6
0.78
0.58
0.7
0.75
0.65
0.8
0.73
0.73
0.9
0.7
0.8
1.0
0.69
0.89
1.1
.67
0.97


## 3. Kinematic Analysis

This section provides analysis for forward and inversed kinematics of the radiosurgery manipulator based on the proposed structure and the lengths of links.


### 3.1. Forward Kinematic Analysis

Figure3 shows the outline of the manipulator and its coordinates in D-H convention. The base coordinate system is located on the mounting plate of the manipulator, and the origin of the coordinate systems 3 and 4 is coincident with each other. The joints 2 and 3 are supported by two links in both sides to reduce deflections. And the dummy joint 6 is not really actuated but added to follow conventional kinematic analysis procedure for 6-DOF manipulators.Figure 3
D-H parameters layout for manipulator.Table2 provides the D-H parameters in mm for length and in degree for angle.Table 2
D-H parameters.
Joint number (i)
a
i-1 (m)
∝i-1 [deg]
di [m]
θ
i [deg]
1
0
−90
−1.604
90
2
0.8
0
0
−90
3
0.73
0
0
180
4
0
90
0
0
5
0
−90
0
90
6
0
180
−0.567
0Using the D-H parameters, homogeneous transform matrices are given by(5)iTi−1=cosθi−sinθi0ai−1sinθicosαi−1cosθicosαi−1−sinαi−1−sinαi−1disinθisinαi−1cosθisinαi−1cosαi−1cosαi−1di0001.By multiplying all homogeneous transform matrices sequentially from the left to the right, the final homogeneous matrix is obtained as follows:(6)T10T21T32T43T54T65=T60=nxoxaxPxnyoyayPynzozazPx.


### 3.2. Inverse Kinematic Analysis

To calculate joint angles from a given target position and orientation, the central position of the accelerator, the origin of the coordinate system 3 is derived at first by the following equation:(7)Pm=PxPyPz′−L6×WL6=θ6.Like a usually 6-DOF manipulator, there are 4 sets of solution forθ1~θ3 in maximum depending to its configurations such as the left or the right arm, elbow up or down. There are two sets of θ4 and θ5 that are determined to satisfy the given orientation for each set of θ1~θ3.Thus, 8 sets of solution can be found if the given end point is inside the workspace. Among them, the joint angle limit criteria are applied to screen the solutions and the solution that is nearest from the current joint configuration is selected as the final solution.A high continuous path can make smooth and efficient movements. To improve the continuity of the path, two successive points are interpolated with spline curves that consist of 5th-order polynomial equation and satisfy given points, velocities, and accelerations as follows:(8)St=a+bt+ct2+dt3+et4+ft5.


## 3.1. Forward Kinematic Analysis

Figure3 shows the outline of the manipulator and its coordinates in D-H convention. The base coordinate system is located on the mounting plate of the manipulator, and the origin of the coordinate systems 3 and 4 is coincident with each other. The joints 2 and 3 are supported by two links in both sides to reduce deflections. And the dummy joint 6 is not really actuated but added to follow conventional kinematic analysis procedure for 6-DOF manipulators.Figure 3
D-H parameters layout for manipulator.Table2 provides the D-H parameters in mm for length and in degree for angle.Table 2
D-H parameters.
Joint number (i)
a
i-1 (m)
∝i-1 [deg]
di [m]
θ
i [deg]
1
0
−90
−1.604
90
2
0.8
0
0
−90
3
0.73
0
0
180
4
0
90
0
0
5
0
−90
0
90
6
0
180
−0.567
0Using the D-H parameters, homogeneous transform matrices are given by(5)iTi−1=cosθi−sinθi0ai−1sinθicosαi−1cosθicosαi−1−sinαi−1−sinαi−1disinθisinαi−1cosθisinαi−1cosαi−1cosαi−1di0001.By multiplying all homogeneous transform matrices sequentially from the left to the right, the final homogeneous matrix is obtained as follows:(6)T10T21T32T43T54T65=T60=nxoxaxPxnyoyayPynzozazPx.


## 3.2. Inverse Kinematic Analysis

To calculate joint angles from a given target position and orientation, the central position of the accelerator, the origin of the coordinate system 3 is derived at first by the following equation:(7)Pm=PxPyPz′−L6×WL6=θ6.Like a usually 6-DOF manipulator, there are 4 sets of solution forθ1~θ3 in maximum depending to its configurations such as the left or the right arm, elbow up or down. There are two sets of θ4 and θ5 that are determined to satisfy the given orientation for each set of θ1~θ3.Thus, 8 sets of solution can be found if the given end point is inside the workspace. Among them, the joint angle limit criteria are applied to screen the solutions and the solution that is nearest from the current joint configuration is selected as the final solution.A high continuous path can make smooth and efficient movements. To improve the continuity of the path, two successive points are interpolated with spline curves that consist of 5th-order polynomial equation and satisfy given points, velocities, and accelerations as follows:(8)St=a+bt+ct2+dt3+et4+ft5.


## 4. Graphic Simulator and Dynamic Analysis

In this study, a MATLAB-based open source code, ARTE (A Robotics Toolbox for Education) [15], is used and modified to simulate the motion of the radiosurgery manipulator under development in 3-dimensional graphic virtual environments.To input geometric model of the manipulator to the simulator, 3D CAD models are converted to STL files in ASCII format.And such dynamic parameters as masses, mass moments of inertia, and centers of masses for each link are also extracted from the 3D CAD files to enable dynamic simulation capability implemented in the toolbox.Figure4 shows the manipulator and a patient lying on a couch in the 3D graphic simulation environment.Figure5 shows the control panel window that displays current joint angles and target inputs to move the manipulator. It is modified from the teach program packaged in ARTE.To command the manipulator to move, the azimuth angle, the elevation angle, and the distant from the origin are directed by operators. The desired target postures can be saved to text files for later dynamic simulations.As an example of trajectory, six postures are given as shown in Figure6 using a function named MoveL() and the joints 4 and 5 are rotated slightly about 50 in some posture using a function named MoveB() that independently moves the joints without moving all the other joints.Figure 4
The 3D graphic simulation environment.Figure 5
The control panel window.Figure 6
An example trajectory six postures.Figure 7
The six postures in the 3D graphic environments.Figure 8
Joint angle and static torque.
(a)
Static joint angle
(b)
Static joint torqueFigure 9
Simulation results.
(a)
The velocity
(b)
The acceleration
(c)
The torque
(d)
The powerFigure 10
S velocity-torque diagram.Figure 11
Prototype manipulator.Figure 12
Graphical use interface.Figure7 shows the six postures in the 3D graphic environments, and Figure 8 shows the static joint torques.As a dynamic simulation, the linear velocity of the end point is given 0.1 m/s, and the velocity, the acceleration, the torque, and the power of each joint are provided in Figure9.Figure10 shows the joint velocity, the joint torque, and the power in velocity-torque diagram for each joint to help motor selection for performing the simulated motion.As shown in Table3, the maximum torque and power for the joint 2 is 6271 Nm and 2317 W.Table 3
Maximum torque and power.
Joint 1
Joint 2
Joint 3
Joint 4
Joint 5
Max. torque
289.61
6271
2446
250.60
87.17
Max. power
59.41
2318
718.32
100.7
6.06


## 5. Manufactured Robot and Control System

The motors and gearboxes are selected based on the dynamic simulation results, and the manipulator has been manufactured as shown in Figure11.As for controlling motors, MR-J4-B series servo amplifiers manufactured by Mitsubishi electric company are selected. They have an interface to a battery-backup optical encoder, and so there is no need to start initialization process at every power-up stage.Those of 5 servo amplifiers are linked through SSCNET optical fibers to MR-MC 241 multiaxis position control board mounted in a PCI slot of a desktop computer operated in MS Windows 7. And also, we consider two types of safety features of the robotic system when danger exists on machines and systems they have to be immediately shut down, in order to protect people, machines, and systems. At first, stoppers are attached to the side of the joint 2 of the manipulator. It can be sure to take preventive method to the manipulator against falling to patients when danger exists on machines. Second, an emergency stop button and a signal emergency stop button in control system are designed to happen only in emergency situations. Pressing the emergency stop button, brakes, is activated in the manipulator axis, and the robot stops. The motor power turns off at this timing.Figure12 shows the test GUI programed in MS Visual C++ environment. This program has included several functions implemented in the simulator.


## 6. Conclusion

This paper introduces a radiosurgery manipulator system that has 5-DOF and a more suitable workspace for its operation. This is achieved by creation of a safety zone around the patient and the treatment couch. In order to design the manipulator, kinematic and dynamic analyses are conducted and implemented in a 3D graphic simulator based on one of MATLAB-based open source robotics toolboxes, ARTE. Resulted from some motion simulations for the manipulator, a prototype manipulator is manufactured, and its control system is implemented.



---

*Source: 4264356-2017-07-24.xml*