Reducing the Frame Vibration of Delta Robot in Pick and Place Application : An Acceleration Profile Optimization Approach

Delta robot is typically mounted on a frame and performs high speed pick and place tasks from top to bottom. Because of its outstanding accelerating capability and higher center of mass, the Delta robot can generate significant frame vibration. Existing trajectory smoothing methods mainly focus on vibration reduction for the robot instead of the frame, and modifying the frame structure increases the manufacturing cost. In this paper, an acceleration profile optimization approach is proposed to reduce the Delta robot-frame vibration.The profile is determined by the maximum jerk, acceleration, and velocity.The pick and place motion (PPM) and resulting frame vibration are analyzed in frequency domain. Quantitative analysis shows that frame vibration can be reduced by altering those dynamic motion parameters. Because the analytic model is derived based on several simplifications, it cannot be directly applied. A surrogate model-based optimization method is proposed to solve the practical issues. By directly executing the PPM with different parameters and measuring the vibration, a model is derived using Gaussian Process Regression (GPR). In order to reduce the frame vibration without sacrificing robot efficiency, those two goals are fused together according to their priorities. Based on the surrogate model, a single objective optimization problem is formulated and solved by Genetic Algorithm (GA). Experimental results show effectiveness of the proposed method. Behavior of the optimal parameters also verifies the robot-frame vibration mechanism.


Introduction
Delta robot is a 3/4 Degree-of-Freedoms(DOFs) parallel robot widely used in manufacturing industries affording the pick and place tasks [1][2][3].A typical sorting work cell is shown in Figure 1.Delta robot is usually mounted on a big heavy frame to operate workpieces from top to bottom.The frame is usually customized according to the environmental conditions, such as workpiece size or conveyor arrangement.Typically, a frame is 2 meters high and 1.5 meters wide, which covers work range of the Delta robot.Thick square steels are welded together to provide a relatively static basis.Other types of base/frames also exist which are summarized in [4].
Although lightweight in mass and small in inertia, the Delta robot can still generate huge reaction force because of its excellent acceleration capability.Moreover, Center of Mass (CoM) of the frame is higher than normal robot base.Therefore, stiffness of the frame is an important factor affecting system noise level, reliability, and precision.Pick and Place Motion (PPM) is the primary behavior of Delta robot.As shown in Figure 1, it moves between two locations along a door shaped trajectory, picking and placing the workpieces.When the robot moves fast, there will be huge reaction force and the frame is forced to vibrate.When the robot moves cyclically, the vibrations may superimpose together and cause significant side effect.Severe frame vibration also affects stability of the frame mounted industrial cameras.Therefore, reducing the frame vibration is a nonnegligible problem when installing new Delta robots.
Smoothing the trajectory and strengthening the frame are two feasible solutions.A lot of trajectory smoothing methods have been proposed to reduce vibration of the Delta robot, such as using special functions to smooth the piecewise linear motions and optimizing the dynamic parameters such as jerk or acceleration.Dai and Sheng propose to use fifth Bézier curve to blend the adjacent positions and linear motions, which can generate a smooth trajectory with continuous second-order derivative and curvature in real-time [5].Huang adopts the elliptical trajectory with modified sine motion profile to plan the pick and place operation, which can guarantee smoothness of the torque.Because the motion planning algorithm only involves two parameters, it can be efficiently implemented.However, no quantitative analysis is performed for the robot vibrations [6].
Kuo studies the mathematical modeling problem of the Delta robot with flexible links.Kineto-elastodynamics and the finite element method are used to derive a mathematical model.Natural frequency analysis is performed to demonstrate the importance of flexible link compensation for precise robot motion [7].However, the frequency analysis is not applied to motion optimization problem and the method is too complex for practical application.In order to improve the tracking accuracy, Liu studies the trajectory planning problem for Delta robot with joint friction and jerk constraints [8,9].Smooth and nearly time optimal trajectory can be generated using the proposed solution.The trajectory is modeled using quadratic and quantic polynomial splines in Cartesian space and septuple polynomial splines in joint space.Zhang and Wang study the trajectory planning and optimization method to obtain smoother and more efficient trajectories for Delta robot [10,11].Béarée introduces input shaping technique into the jerk-limited profile, derives the damped-jerk profile to cancel the residual vibration of an undamped system, and proposes several practical tuning rules for the equivalent input filter parameters [12].The dominating flexible mode analysis closely relates to system physical parameters; therefore, it is difficult to apply in the vibration suspension problem of Delta robot and frame.
The causes of frame and robot vibration are different.Robot vibration originates from discontinuity of the drive torque [13,14], nonlinearity of the robot structure [15][16][17], or flexibility of the links [7].Frame vibration is excited by the force coupled through the static base.It originates from motion of the robot.No matter how smooth the robot motion is, there is always excitation force imposing on the frame.The only ways are modifying either the frame structure or the robot motion.
A commonly used solution is to modify the frame structure.Increasing the frame mass and stiffness or integrating reinforcement beam can all alter frame's natural frequency characteristics.The designing process involves 3D modeling, Finite Element Analysis (FEA) on the stiffness, deformation, and frequency characteristics of the new structure [4].Apparently, it is not an economic solution because the frame should be specifically designed for each application.More materials and costs are required.Even though, for high speed pick and place motions, the vibration is still significant.
The other method is to modify the robot motion.However, it is restricted by task and physical constraints.A feasible way is to customize the dynamic motion parameters for each application.This can also change frequency characteristics of the excitation signal imposed on the robot frame.Consequently, the frame vibration can be reduced without modifying the hardware design, software architecture, and sacrificing the efficiency.Furthermore, because the proposed method does not change the underlying motion smoothing algorithms, there is no side effect on the robot vibration.
Vibration signal contains rich information about the state of the robot/machine and hence can be used to assist the highlevel diagnosis and decision-making process.A lot of model and signal based methods have been proposed to solve the fault detection problem for machines [18][19][20] and robots [21,22].In this paper, the vibration signals are used to find the optimal dynamic motion parameters.
There are two objectives for the optimization problem: reducing the frame vibration and keeping efficiency.It is noticed that they have different priorities.Efficiency is the primary goal.Therefore, in this paper, a specific fitness function is designed to fuse them together safely.Although a lot of multiobjective optimization algorithms have been proposed [23,24], single objective optimization problems are easier to handle because there is no trade-off between each objective.
The paper has the following contributions.(1) A frequency domain analysis is performed on the PPM acceleration profile.The internal relation between dynamic motion parameters and distribution/ amplitude of the harmonic components is established, which proves the possibility of reducing the frame vibration by optimizing the acceleration profile.(2) A surrogate model-based parameter optimization method is proposed to solve the practical vibration reduction problem.The method uses a nonparametric regression tool: Gaussian Process Regression (GPR) to generate the surrogate model which makes it easier to be applied in other acceleration profiles.The method fuses two objectives into one single fitness function based on their priority differences.It can be easily integrated in existing optimization algorithms.   1 shows a commonly used configuration.The workpieces are supplied through conveyors and the robot picks them from conveyor and places them on target locations (boxes, plates, or other conveyors).Typically, in order to facilitate tracking of the conveyor, moving direction of the conveyor is aligned with X-axis of the Delta robot.Therefore, pick and place motion (PPM) is mainly within YZ plane as shown in Figure 1.In the paper, the conveyor is assumed to align with X-axis for convenience.Results of aligning with Y-axis are similar.

Problem Formulation
Pick and Place Motion (PPM) is a door shaped trajectory consisted of one horizontal, two vertical, and two smooth transition segments as shown in Figure 2 [25].The trajectory formed by points P0→P3 in Figure 2 is a typical PPM, where P0 is the pick-up location and P3 is the put-down location.P1 and P2 define the intermediate transition positions.The robot picks workpiece up at P0 and moves it to P3 along the predefined trajectory.During this process, the robot accelerates and decelerates cyclically and swiftly.Huge reaction forces are imposed on the robot and the frame.Therefore, smoothness of the PPM greatly affects the vibrational amplitude.
When the robot moves along PPM cyclically and swiftly, the reaction force will make the frame vibrate.This will generate noise and affect the robot accuracy.The problem is to reduce the frame vibration while keeping system efficiency based on existing hardware and software configuration.

Acceleration Profile and Frame Vibration.
Efficiency and frequency of the robot motion are controlled by the acceleration profile, which refers to how the velocity and acceleration change in the trajectory tracking control process.The jerklimited profile is a smooth command pattern used by modern motion systems such as machine-tools and industrial robots.Jerk is time derivative of acceleration.The introduced jerk can eliminate the discontinuities of acceleration and reduce the residual vibration.Figure 3 shows an example plot of jerk-limited profile.By using the bang-bang control law, a time optimal accelerating trajectory is derived.A jerk-limited profile equals to a trapezoidal or a triangular acceleration profile.
A point to point motion profile can be expressed using the following equation:

𝑎 (𝑡 | 𝑆, 𝐽
where  is the stroke,  max ,  max ,  max are the dynamic motion parameters for acceleration profile, Δ is the time when the robot is in uniform velocity state, which is proportional to , and Δ 1 , Δ 2 are the accelerating time.
A detailed discussion on the jerk-limited acceleration profile can be found in [12,26].Besides the symmetric profile, there are also many variations such as the asymmetric acceleration profile in [27].No matter which type the robot uses, the profile all can be parameterized.For the standard jerk-limited acceleration profile, the dynamic parameters  max ,  max ,  max greatly affect the robot motion efficiency and smoothness.
Typically they are chosen empirically to achieve small cycle time or to reduce the robot vibration.Significant frame vibration originates from mismatch of the robot and frame frequency characteristics.Before installation, the frame can be reinforced to match the application.However, this increases the material and labor costs.The acceleration profile is the key to change the robot motion frequency characteristics.Therefore, an economic  and general applicable solution is to customize the dynamic parameters according to specific application based on the existing frame configuration.Instead of empirically setting those parameters, it is necessary to find an optimal solution for efficient and lower vibration motion.The following sections further discuss this problem from mechanism to solutions.

Frequency Domain Analysis of PPM and Frame Vibration
When Delta robot performs the pick and place tasks, acceleration and deceleration occur cyclically.Consequently, a periodic reaction force drives the frame to vibrate.The robot motion can be seen as the input while the frame vibration is the output.The robot-frame union is a black box with certain frequency response characteristics.To understand the vibration excitation mechanism and find methods to reduce the vibration, it is necessary to analyze the robot motion and robot-frame union in frequency domain.

Spectrum of the Pick and Place Motion. As shown in
Figure 1, in pick and place applications, Delta robots usually move between two locations cyclically.In the door shaped trajectory, horizontal segment is the dominant one which covers most of the stroke.Therefore, it can be treated as a cyclic point to point motion for simplicity (P1→P2, P2→P1).When applying trapezoidal acceleration profile, the corresponding motion is illustrated in Figure 4(a), where the thin solid line is the acceleration.
Considering the point to point motion (1) and symmetricity, acceleration profile of the cyclic motion can be modeled using which is a piecewise linear function as shown in Figure 4(a), where  is an elementary point to point motion whose period is ;  is the stroke between initial and target positions;  max ,  max ,  max are the maximum velocity, acceleration, and jerk.
Because the movement between P0 and P3 are symmetric, they both cover half of the period.Therefore, one can get where   is the base frequency and   is the corresponding angular frequency.From Figure 4(a) and ( 5), it is able to find that the acceleration is actually an odd harmonic function.
Considering the Fourier transformation with  =   , there is If  is odd, − − = 1, which will lead to an identity.If  is even, − − = −1, which will lead to There is (  ) = 0. Therefore, for the cyclic PPM motion, only odd harmonics such as   , 3  , 5  , ⋅ ⋅ ⋅ exist.Figure 4(b) shows Fast Fourier Transform (FFT) results for the acceleration profile shown in Figure 4(a).It is clearly seen that the base frequency is 0.44Hz and the harmonics are 1.32Hz, 2.2Hz, which is consistent with the theoretical analysis.4(a), (4), and (1), it is able to calculate the amplitude-frequency characteristics of the acceleration profile.Noticing that the jerk profile is simper (piecewise constant) than the acceleration profile (piecewise linear), the Fourier transformation can be performed on jerk based on the facts that

Amplitude-Frequency Characteristics Analysis. According to Figures
where () refers to the jerk profile.The derivation process of the differential property (8) for Fourier transformation can be found in [28].The Fourier expansion coefficient for acceleration profile is where Assuming  is odd, using Euler formula to expand above equation and merge similar items, one gets Substituting ( 12) into (10), one gets According to ( 13), (2), and (3), one can find that the dynamic motion parameters  max ,  max ,  max are the only factors affecting the acceleration profile, cycle time, and amplitude-frequency characteristics.Larger  max and smaller cycle time  tend to increase the amplitude.However, items in the brackets can also be reduced by properly choosing those dynamic motion parameters.This indicates that it is possible to reduce the frame vibration by tuning the acceleration profile instead of modifying the hardware, and a feasible way is to optimize dynamic motion parameters according to given constraints.
An illustrative example is shown in Figure 5.The cycle time  is set to be 0.5s, i.e., locations of the harmonic components are fixed.Different acceleration profiles are generated by changing  max = [10000, 50000]mm/s 2 and  max = [1000000, 2000000]mm/s 3 as shown in Figure 5(a).
Amplitudes of base frequency, 3rd and 5th harmonics are plotted in Figure 5(b).It is seen that amplitudes vary according to the dynamic motion parameters.Globally, the amplitude is proportional to  max .However, local minimums also exist which can generate lower vibration without lowering production efficiency.This is consistent with (13).

Measures to Reduce Frame Vibration.
Different from the robots, the robot frames usually have no fixed models.They are customized for each application to meet the specific conveyor width and height, workspace arrangement, and workpiece size requirements.Although weighs hundreds kilograms, it can still vibrate with some natural frequency.This frequency is decided by the frame material, structure, dimension, and mounting type.When external excitation signal is close to the natural frequency, the frame generates significant vibration.
Frequency response of the frame can be described using where () is the amplitude-frequency characteristics of robot-frame union whose natural frequency is   .() and () are the robot motion and frame vibration signals in frequency domain.For a given robot-frame union, () is determined.Without changing the hardware, the only way is to tuning the robot motion ().
The robot motion can be tuned in two aspects.
(i) Harmonics distribution For the PPM motion, only odd harmonics exist in the acceleration profile.The proximity between (2−1)  and   is defined by It is easy to verify that 0 ≤   ≤   .  = 0 will lead to resonance.Bigger   can decrease the vibration.The harmonics distribution is closely related to cycle time  because   = 1/.One can alter  to tune the distribution to avoid the natural frequency.
(ii) Natural Frequency Amplitude The other method is to decrease the amplitude of (  ).According to (13), this can be realized by increasing , decreasing  max , and optimizing the dynamic motion parameters.However, it is hardly possible to directly find such parameters because (  ) and  are internally related in a complex reverse way.
Remark 1.As shown in Figure 1, typically X-axis of Delta robot is aligned with the conveyor and the pick and place motion mainly lies in the Y-Z plane.Because the vertical motion is usually very small compared to the horizontal motion, the PPM motion is mainly in Y direction.Therefore, the acceleration profile and the vibration model are both considered as one dimensional in this paper.
The above theoretical analysis reveals the mechanism for frame vibration and shows two possible methods to reducing the frame vibration by tuning the acceleration profile of PPM.However, it is not suitable for practical application because of the following four reasons.(1) Real PPM motion has vertical segments and smooth transition segments.(2) Formula of acceleration profile may change for certain parameters.(3) Real PPM may derivate from Y-axis.(4) There is a lack of robot-frame union model.

Modeling and Parameter Optimization
Because of the complex robot-frame structure, it is hardly possible to derive an accurate analytic model for the excitation and vibration process.In previous section, robot motion and frame vibration are analyzed in frequency domain.Guidelines for reducing the frame vibration are identified.However, there still lacks quantitative methods to optimize such dynamic motion parameters.
This paper proposes a surrogate model-based parameter optimization solution.Instead of using the analytic model, a surrogate model is obtained by performing regression on real datasets.The modeling and optimization process is shown in Figure 6.
The whole process consisted of three parts: data acquisition, signal processing and modeling, and dynamic motion parameter optimization.

Data Acquisition.
To obtain a surrogate model for (14), all related variables should be able to be extracted from the recorded datasets.According to (14), robot motion and frame vibration are both relevant.Therefore, acceleration profile of robot motion and frame vibration is recorded simultaneously.The recorded signals are time series () = {  (),   ()}, where   () is the robot motion while   () is the frame vibration.
The arguments are p = { max ,  max ,  max }, which are known in advance when performing the evaluation experiments.Therefore, the recorded dataset is In order to minimize the effect of random factors, the robot repeats the pick and place tasks several times for each parameter candidate p.

Signal
Processing and Feature Extraction.Performances of the dynamic motion parameters are evaluated from two aspects: efficiency of the robot motion and amplitude of the frame vibration.In this paper, efficiency is measured by cycle time .Vibration is measured using maximum acceleration â.Both indices can be calculated from the datasets.
The surrogate numerical model is obtained by performing regression over the datasets.The dynamic motion parameters are  max ,  max , and  max and the performance indices are â, .However, the unit and value range of those variables cause serious problems.Typically, their values have significant differences, ranging from 0-2000000.Without normalization, some of the variables may be submerged by others and the modeling process may fail.To treat each variable equally, a normalization operation is required.Take  as example; the normalizing process follows where T refers to the normalized ; () and () stand for the mean and standard division of .This operation can make all the variables zero mean and with identical variance.The original variable can be recovered easily through Same operations can be performed on â and  max ,  max , and  max .After normalization, inputs and outputs of the model are converted to Jmax , Ãmax , Ṽmax , and T, ã.The latent models can be represented using where   ,   are the underlying functions.

Gaussian Process Regression Based Modeling.
Considering the uncertainties in datasets, numbers of samples, and nonlinearity of the underlying function, Gaussian Process Regression (GPR) is used to modeling the datasets.GPR is a nonparametric tool that can handle the modeling problem with noisy observations and system uncertainties [29].It has been applied in robot process modeling and parameter optimization [30].
For the modeling problem (21), there is The goal is to derive a model to approximate the underlying function   (x) and   (x).Because there are many random factors in the robot system and the observations are also noisy, they can be seen as two Gaussian Processes and can be modeled using GPR.
For a Gaussian Process (x), a set of multivariate Gaussian random variables F = {(x 1 ), (x 2 ),, ⋅ ⋅ ⋅ , (x N )} can be defined over X, where {x  ∈ X}  =1 are the inputs and (x i ) are Shock and Vibration the corresponding outputs.X is the parameter space which is defined over R  .Instead of deriving an explicit form for (x), it is approximated by its covariance function (x, x  ): (x) ∼ GP(0, (x, x  , ) where x and x  are two arbitrary variables in X.
The goal of model construction is to find a covariance function (x, x  , ) which fits the data set (X, y) best.Suppose (x) is a candidate latent function ( in short), the marginal likelihood of  and  given the data set (X, y) can be derived by performing marginalization over the function   (y | X, , ) = ∫  (y | X, , , )  ( | X, , )  (23) where  is the hypothesis on the structure of the covariance function,  are the hyperparameters, X, y are the sample data sets, and is the posterior probability.
In order to simplify the modeling process, a log operation is performed to convert multiplication into addition.The Gaussian assumption makes it possible to derive an analytic solution for the marginal log likelihood log Then the hyperparameters  can be found by maximizing the marginal log likelihood function.
This optimization problem is solved by Conjugate Gradient Algorithm.After training, the model is constructed and can be used for prediction.
The joint distribution of the modeling dataset (X, y) and test dataset (X * , y * ) is where (X  , X  ) (X  and X  refer to X and X * ) is the covariance matrix whose element   at  ℎ row and  ℎ column equals to (x  i , x  j ).By deriving the conditional distribution, the predictive function is obtained where (y * ) is the predicted mean for y * and V(y * ) is the variance.
In summery, a covariance function should be specified for the GPR model.In the proposed data acquisition procedure, the pick and place motion are repeated for several times for each parameter combination.The random noise in the extracted feature is averaged.Therefore the following covariance function is used: where two Squared Exponential covariance functions with Automatic Relevance Determination (SEARD) items are used to provide feasibility in different length scales   .Considering the amplitude   , there are 8 hyperparameters.Based on (29) and the obtained dataset (22), the optimal hyperparameters are found by optimizing (26).Finally two models are created.
T =   ( J , Ã , Ṽ ) ã =   ( J , Ã , Ṽ ) . ( These models exactly recover the original datasets and provide smoothly weighted interpolations on the unevaluated dynamic motion parameters.It is able to predict the performance of the dynamic motion parameters using (28).
Based on the GPR surrogate model, one can perform global optimization to reduce frame vibration.For efficiency consideration,  must be smaller than the threshold value   .Therefore, the first constraint is The other goal is to minimize the frame vibration, which can be expressed using In general multiobjective optimization problems, usually different objectives may conflict with each other.However, for this problem, those two goals, keeping efficient (31) and reducing vibration (32), have different priorities and requirements.It follows the following two rules: (1) If a solution cannot meet the cycle time constraint (31), it is bad no matter how small the vibration is.(2) If two solutions both satisfy  1 <  2 <   , the one with smaller vibration is better.
In other words, the optimal solution is neither the fastest nor the vibration-free one.

Unified Fitness Function.
These two objectives have different priorities.Keeping efficiency is the primary goal while minimizing the vibration is the secondary one.Simply summing them together ( =  1  +  2  ) is not applicable because it confuses each index and leads to undesired extreme results.
To avoid handling multiobjective optimization problem, the following fitness function is proposed to cooperate the above two rules into one goal, where step() is defined as Clearly, the proposed rule (33) mixes the two objectives together in a right way.Because  is normalized within range [-1, 1], the first item 2 × step( T − T ) dominates (33) if step() = 1.The second item ã cannot surpass the first item, which satisfies rule (1).Meanwhile, parameter candidates with  ≤   is treated equally and their differences are evaluated by the second item ã, which satisfies rule (2).
Equation (34) divided the parameters into two classes.For parameters violates rule (1), (33) maps them in the range of [1,3].For parameters satisfies rule (1), (33) maps them in the range of [-1, 1].Therefore, the "good" parameters are distinguishable.Existence of the solution to this problem can be guaranteed.
The only hyperparameter in the fitness function is   .It is assigned based on practical requirements.The Delta robot can operate at very high speed in pick and place applications.Cycle time of state-of-the-art Delta robot is 0.2s-0.6s(varies according to particular testing configuration), i.e., 2Hz-5Hz.However, not all applications need such high speed.It is necessary to optimize the robot behavior according to real requirements, such as the pick and place path configuration and the workpiece supply rate.The first one is identified as  0 ,  1 ,  2 ,  3 which are modeled implicitly in (30).The second one is identified as required cycle time   and introduces a constraint  <   for the optimization process.

Parameter Optimization.
The fitness function (33) formulates the following single objective optimization problem: With the surrogate model and the fitness function, it is able to embed the model into an optimization tool such as Genetic Algorithm (GA) [31].The grid search method also works but there is a trade-off between accuracy and efficiency.Figure 7 shows diagram of the optimization process.
The GA has been well implemented in MATLAB [32].The real encoding GA is chosen and the population size is set as 50."Stochastic uniform", "scatter", and "mutationadaptfeasible" are used to realize the selection, crossover, and mutation operation.The crossover rate is 0.8 and error tolerance is 1 −6 .Starting from a group of randomly initialized parameter candidates, the evaluation and regeneration process are iteratively executed until the terminate criterion condition is satisfied.
Remark 2. Because GPR is a nonparametric regression tool that heavily relies on datasets and its computational complexity is O( 3 ), it cannot be applied online when  is large.On the other hand, accuracy of the surrogate model cannot be guaranteed if size of the datasets is too small.The signal acquisition, modeling, and optimization process is performed offline (before actual production) to find better dynamic motion parameters which are suitable for practical production requirements (PPM trajectory, cycle time).
In this paper, GA is used to solve (35).Other optimization methods such as Conjugate Gradient Algorithm can also be integrated.Vibrations can be measured by several kinds of devices, such as laser tracker, IMU, strain sensor, and even radio signals [33].In this paper, a nine axis IMU module is used to measure the frame vibration.The used module is based on the InvenSense MPU9250 chip which consisted of a 3-axis gyroscope, a 3-axis accelerometer, and a 3-axis magnetometer.It is an easy-to-use and low cost IMU device with acceptable speed and precision.Two IMU modules are integrated.The #1 IMU sensor is mounted on the robot frame as shown in Figure 8  Figure 10 gives the overall flowchart of the signal acquisition process.

Experiments and Discussion
Each parameter candidate is repeated for four times to reduce effect of noises and random factors.Figure 11 shows a typical vibration datasets.
Each sample lasts for 6 seconds and contains 4 pick and place cycles.12.The horizontal axis is cycle time while the vertical axis is the vibration.The blue curve is a fitted polynomial for comparison.From this figure, one can find that the vibration is inversely proportional to the cycle time in general.This is consistent with (13) because  is in the numerator.Meanwhile, significant variances exist in the vibration amplitude even for same cycle time.Take  = 0.65s, for example, the maximum vibration is nearly double of the minimum one.This indicates the possibility of decreasing the frame vibration while maintaining the efficiency.
Figure 13 shows a slice of the generated GPR model with   = 3000mm/s.From Figure 13, the saturation phenomenon (cycle time cannot be further reduced if velocity is limited) and vibration variability ("peaks" and "valleys") are clearly demonstrated.The goal of this paper is to find those "valleys" and refine the system performance.

Optimization.
For the optimization problem defined using (33) and (35), Genetic Algorithm (GA) is used to perform the optimization process.Given desired cycle time   , the corresponding optimal dynamic parameters can be found.GA is configured with the settings presented in Section 4.4.3.The optimization process is shown in Figure 14.
It can be seen that the algorithm converges after about 10 generations.The optimal fitness values are all within range [-1, 1], which indicates the optimal parameters all satisfy the cycle time constraint (31).After optimization and recovery, the optimal parameters and the corresponding performance indices are listed in Table 2.
Applying these optimized dynamic motion parameters on the experiment platform and recording the vibration, one can find that the predicted values are close to the real ones.The results verify effectiveness of the proposed modeling and optimization method.Waveform of the optimal parameter for   =0.8s is shown in Figure 18.   15.
From this figure, one can find that the resonance frequency lies around 7.6Hz, from 6.76Hz to 9.2Hz.

Spectrum of the Robot Motion Acceleration Profile.
According to Figure 12, cycle time ranges from 0.56s to 1.2s, i.e., the base frequency is less than 1.8Hz and the 5th harmonics component lies in the range of [6.76, 9.2].
For cases with  >   , the 5th harmonics are outside the resonance range.Although the 7th harmonics are within the resonance range, typically, their amplitudes are smaller than the 5th one.Therefore, the frame vibration is also small.For cases with  <   , because of resonance, the vibration increases significantly.The frame vibration is dominated by the amplitudes of the 5th harmonics.The above conclusion is consistent with the results given in Figure 12.
To verify the influence of 5th harmonics, a comparative example is given in Figure 16.
Two groups of dynamics parameters [  ,   ,   ] =[2000,50000,2500000] (Case B:Blue solid line) and [3000, 50000, 1750000](Case R:Red dash line) are examined.Both cases have same cycle time 0.65s.The harmonic distributions are identical (because some of the accelerations are beyond the sensor range, spurious harmonics also exist in the spectrum).
x 10 4 ) x 10 4 (b) Vibration Comparing the amplitude of each harmonics, Case B is larger than Case R except the 5th harmonics.However, frame vibration of Case B is smaller than Case R, which proves importance of the amplitude of the close-to-naturalfrequency.
Figure 17 shows the spectrum of the optimal parameters listed in Table 2.One can find that both of the above two conditions are satisfied.The 5th and 7th harmonics are away from the resonance frequency and also the amplitude is low.Therefore, the frame vibration is minimized.
The waveform of those three groups of parameters is shown in Figure 18.It is clearly seen that the dynamic motion parameter optimization process indeed reduces the frame vibration.

General Applicability of the Proposed Methods.
The central topic of this paper is to find a way to optimize dynamic   motion parameters of PPM motion so as to reduce the induced robot base/frame vibration.As shown in Figure 6, the proposed method uses two vibration datasets obtained from robot body and frame to build a model characterizing the relation between dynamic motion parameters and system performance.Then based on the model, an optimization problem can be formulated and solved to find the best parameter.
From Figure 6, one can find that only the "evaluating" process is related to the actual robot system.Although the derivation and experiments are based on Delta robot, it can be easily applied to other robots with similar PPM motions, such as SCARA robot and palletizing robot.However, because cycle time of non-Delta robots is typically longer and they are usually mounted on the ground/table, the frame vibration is not so obvious.
It is easy to apply the proposed method to a different robot.The only modification is to write a simple robot program that uses different parameters to perform PPM motion or just add several instructions to the existing program to modify the dynamic parameters.Then the dataset can be recorded automatically and processed offline.

Conclusion
This paper focuses on the frame vibration suppression problem for Delta robot in pick and place application.By analyzing acceleration profile of the pick and place motion in frequency domain, distribution and amplitude of harmonics in the robot motion are identified.They are proved to be internally related to dynamic motion parameters such as jerk, acceleration, and velocity.Guidelines for reducing the vibration are found.In order to solve the practical vibration reduction problem, a surrogate model-based optimization solution is proposed.The model is based on Gaussian Process Regression (GPR) and can be constructed from the recorded acceleration signals on real robot platform.In order to achieve the goal of keeping efficiency and reducing vibration, their priority differences are utilized.They are fused together in a single fitness function.This single objective optimization problem is solved by Genetic Algorithm (GA).The proposed dynamic motion parameter optimization method is flexible and generally applicable and can work with other PPM generation algorithms or acceleration profiles.
The future works focus on two main aspects: developing an embed measurement unit with general applicable interface and integrating them on different robot platforms; exploring the possibility of online dynamic optimization problem.

Figure 1 :
Figure 1: The Delta robot based pick and place work cell.

Figure 2 :
Figure 2: The door shape pick and place motion.

Figure 3 :
Figure 3: An example plot of jerk-limited acceleration profile.

Figure 4 :
Figure 4: An illustrative example of the trapezoidal acceleration profile and the corresponding spectrum for cyclic motion.

Figure 5 :
Figure 5: An example illustrating the relation between the dynamic motion parameters and amplitude-frequency characteristics.

Figure 6 :
Figure 6: Diagram of the signal processing, modeling, and optimization process.

4. 4 .
Surrogate Model-Based Dynamic Motion Parameter Optimization 4.4.1.Objectives and Priorities.The goal of the optimization problem is to reduce the frame vibration without sacrificing robot efficiency.There are two indices evaluating the performance: the cycle time  and the frame vibration amplitude .

Figure 7 :
Figure 7: The surrogate model-based parameter optimization problem.
to detect the frame vibration signals.It is located between the frame center and the frame boundary.The location is rough because what really matters is the relative amplitude of the vibration.The #2 IMU sensor is attached on the moving platform to measure the reference driving signals.Y-axis of both sensors is aligned to the Y-axis of the robot, which is the primary direction of the pick and place motion.Both sensors transfer the data at 200Hz back to computer through USB-UART interfaces.
Figure 11(a) is the acceleration profile of the robot's moving platform.Figure 11(b) is the frame vibration

Figure 11 :Figure 12 :
Figure 11: Typical vibration signals recorded by the IMU sensors.(a) Acceleration of moving platform along Y-axis; (b) frame vibration along Y-axis.

Figure 13 :
Figure 13: Slice of the generated model.
Spectrum of the frame vibration

Figure 16 :
Figure 16: Comparative results of two cases with same cycle time  = 0.65s.However, the frame vibration amplitudes have significant difference: 0.7353 m/s 2 and 1.2286 m/s 2 .

Figure 17 :Figure 18 :
Figure 17: Spectrum of the acceleration profile generated by the optimal dynamic motion parameters.
2.1.Delta Robot in Pick and Place Application.Delta robots are typically used in pick and place applications.Figure

Table 1 :
Motion related parameters.(mm/s)(mm/s 2 )   (mm/s 3 ) Both signals are in Y-axis, i.e., the main direction of the pick and place motion.From this figure, it is able to identify that the robot performs the pick and place motion for four times.During the motion process, the frame is forced to vibrate.Even after the robot stops at 3.5s, the frame still vibrates for nearly 3 seconds.Dynamic motion parameters that may affect the smoothness and vibration are listed in Table1.The maximum, minimum, and interval of each variable are presented.Totally 9 × 5 × 7 = 315 experiments were conducted.5.2.Experimental Results and Analysis5.2.1.Data Processing and Modeling.After signal processing, performance indices of all the 315 data sets are extracted. Ccle time  and vibration amplitude â are plotted in Figure

Table 2 :
The optimized dynamic parameters.
5.3.1.Frame Frequency Response Characteristics.From Figure11, it is able to find that when the external excitation disappears, the frame still vibrate for several seconds.Applying the Fourier transformation, spectrum of the robot motion and frame vibration can be obtained.Dividing the corresponding frequency component, one gets the frequency response chart as shown in Figure