Research Article  Open Access
Shouyan Chen, Tie Zhang, Yanbiao Zou, Meng Xiao, "Model Predictive Control of Robotic Grinding Based on Deep Belief Network", Complexity, vol. 2019, Article ID 1891365, 12 pages, 2019. https://doi.org/10.1155/2019/1891365
Model Predictive Control of Robotic Grinding Based on Deep Belief Network
Abstract
Considering the influence of rigidflexible dynamics on robotic grinding process, a model predictive control approach based on deep belief network (DBN) is proposed to control robotic grinding deformation. The rigidflexible coupling dynamics of robotic grinding is first established, on the basis of which a robotic grinding prediction model is constructed to predict the change of robotic grinding status and perform feedforward control. A rolling optimization formula derived from the energy function is also established to optimize control output in real time and perform feedback control. As the accurately model parameters are hard to obtain, a deep belief network is constructed to obtain the parameters of robotic grinding predictive model. Simulation and experimental results indicate that the proposed model predictive control approach can predict abrupt change of robotic grinding status caused by deformation and perform a feedforward and feedback based combination control, reducing control overflow and system oscillation caused by inaccurate feedback control.
1. Introduction
The deformation occurs during robotic grinding process has significant impact on robotic grinding dynamic and robotic grinding performance [1, 2]. Two ways are mainly presented in current studies to solve this problem. One way is to optimize mechanical structure of robotic machining system or increase stiffness and stability of robot machining system. The other is to adjust machining trajectory by offline planning or realtime force control according to the robot dynamic model [3].
The realtime force control approaches presented in current studies include adaptive control, fuzzy control, and control based on neural network [4]. Mendes et al. [5] proposed an adaptive fuzzy control approach, which is based on Hybrid force/motion control system, to cope with contact issues between robot and a given surface. Fu et al. [6] proposed an adaptive fuzzy force control model, which includes a speed control loop and a position control loop to control both feed rate and position of robot, to achieve stable robotic deburring control. Yen [7] proposed an adaptive control method based on recursive fuzzy wavelet neural network to optimize motion control parameters of threeaxis robot in real time.
The above feedback control approaches implement only when trajectory deviations appear, which may result in overshoot, control overspill, and system oscillations. In view of this, some scholars attempt to implement model predictive control approach to achieve feedforward compensation control [8, 9]. Many nonlinear model predictive control approaches are then proposed. Wilson [10] discussed the performances of three model predictive control approaches applied to robot system control. The three approaches are nonlinear model predictive control (nMPC) approach, PIDbased nonlinear model predictive control (PID nMPC) approach, and simplified nonlinear model predictive control (SnMPC) approach. The results of discussion indicated that the performance of nonlinear model predictive control approach is susceptible to system model errors. Some scholars try to improve nonlinear model predictive control approach by using intelligent algorithms such as neural network [11–13]. Li [14] proposed a nonlinear model control method based on neural dynamic network, where the neural dynamic network is used to obtain optimal values of the formulated constrained quadratic programming (QP) problem derived from the cost function of nonlinear model predictive control model. Zeng [15] used Gaussian radial basis function (RBF) neural networks to improve the nonlinear model predictive control approach applied in the control of nonlinear multivariable systems. Dalamagkidis [16] proposed a nonlinear model predictive control approach based on recurrent neural network to achieve the predictive control of propeller selfrotation process while unmanned aerial vehicle engine is damaged.
In this paper, a model predictive control approach based on a deep belief network (DBN) is proposed to control robotic deformation and reduce rigidflexible effect on robotic grinding dynamics. The following parts are arranged as follows: Firstly, the dynamic model of robotic grinding is established with the consideration of rigidflexible coupling effect. Based on this, the model predictive controller of robotic grinding is designed. Since the accurate parameters of robotic grinding dynamics model and model predictive controller are hard to acquire, a deep belief network is designed to access nonlinear predictive model of robotic grinding. Simulation and experiments are finally carried out to verify performance of the proposed approach.
2. RigidFlexible Coupling Dynamics of Grinding Robot
Traditional grinding dynamics model can be expressed as [17]where is the cutting force; is the system mass matrix; is the system damping; is the system dynamic stiffness; is the grinding tool position; and , are its first and second derivatives. Since the stiffness of CNC is large and the deformation is small, the grinding tool position is approximate to the planned position. However, the stiffness of robot is not sufficient which may lead to large deformation and large deviation between grinding position and planned position. Therefore, the relationship of the grinding position , the planned position , and the deformation can be expressed asSimilarly, the relationship of the force acts to robot endeffector , the cutting force , and the force caused by robotic grinding deformation is
2.1. Robotic Grinding Deformation
Robotic grinding deformation consists of extrusion deformation and periodic deformation. The extrusion deformation is caused by relative motion between grinding tool and workpiece, while the periodic deformation is caused by relative motion between blades and workpiece. Therefore the grinding deformation can be expressed as where is the extrusion deformation; is the periodic deformation. The generation of extrusion deformation is shown in Figure 1. The grinding tool is driven toward the workpiece at a feed rate of to perform cutting. According to traditional grinding theory, there are sliding and extrusion state before the actual grinding is conducted. Based on this, an assumption is made that the deformation generated when the grinding tool cutin workpiece is mostly derived from extrusion deformation, written as where is the dynamic stiffness matrix of robot grinding system; is the extrusion deformation at time . The value of extrusion deformation can be regarded as the accumulation of the difference between the feed rate and the removal speed from time to time , written as Similarly, the relationship between periodic deformation force and the periodic deformation can be expressed aswhere is the circular frequency of grinding tool; is the number of blades,; is the tool rotating speed; is the corresponding phase; is the corresponding maximum cutting deformation. Therefore, the actual grinding position can be obtained by (6) and(8): Substituting (9) and (1) into (3), the dynamic model of robot grinding can be expressed asFor the convenience of discussion, is named as robot grinding force in the following text.
3. Model Predictive Control Based on Deep Belief Network
3.1. Problem Description
During robotic grinding, the actual grinding position is expected to be approximate to the planned position However, according to (10), as the grinding robot rigidity is insufficient and the robotic grinding deformation is considerable, (11) is hard to realize. Therefore, an ideal deformation is defined as well as the corresponding ideal trajectory of robot grinding . Equation (11) can be written asAssuming that there is a corresponding relation between the grinding position and the robotic grinding force , when the system parameters (such as feed rate, grinding tool rotation speed, and material of workpiece) are constant. The control target in this paper can be converted intowhere is the ideal grinding force corresponding to the ideal grinding position .
3.2. Model Predictive Control for Robotic Grinding
Define robot grinding force deviation as Assume that robot grinding force deviation at time is affected by robot grinding force deviation and system transfer function at time . where is the control output which takes feed rate as variables; is the weighting factor of the robot grinding force deviation. Assume . Rewrite (15) according to (14)where Define , , . According to (6) there is , and (17) can be rewritten as
3.2.1. Robotic Grinding Trajectory Prediction
Substitute (19) and (15) into (14) to obtain the relationship between robotic grinding force deviation and feed ratewhere . Therefore, the robotic grinding force deviation matrix from time to time can be acquired by represents the predictive of robotic grinding force deviation at time according to the system information at time . Assume that there is a relationship between change of feed rate and input voltage: where and are weight factors. The prediction of feed rate increment from time to time can be acquired by where , respectively, represent the changes of feed rate and control voltage at time predicted according to the system information at time .
3.2.2. Control Output Optimizing
In order to optimize control output voltage, an energy function is defined aswhereSubstitute (23) into (26):DefinewhereAccording to (26), the extreme value of energy function can be obtained by . Therefore, there is extreme value when: Substitute (23) into (30):where
3.2.3. Stability Analysis
According to (22), (23), and (30), there iswhereAccording to the model predictive control principle [18], when the eigenvalues of the matrix are all within the unit circle, the control system is asymptotically stable.
3.2.4. Robot Grinding State Prediction Model
Merging (33) and (31), there iswhere and represent predictive model parameters. Then, the deep belief network is introduced to obtain predictive model parameters.
3.3. Predictive Model Based on Deep Belief Network
According to the deep belief network principle,Comparing (37) and (35), it can be found that predictive model parameter can be obtained by weights product :Therefore, an assumption is made that can be obtained by the product of . And subassumptions are made that ; ; . Based on these assumptions, the deep belief network is constructed, which is comprised of three layers, including two layers of Restricted Boltzmann Machine (RBM) and one BackPropagation (BP) layer. Two layers of RBM are used to obtain weights and , while the BP layer is used to obtain . The inputs of DBN are grinding force deviation , feed rate , increment of feed rate , and increment of control output voltage , while the outputs are and . The energy function of general RBM is [19]where are explicit layer and hidden layer; , are, respectively, network weight, explicit layer bias, and hidden layer bias. Consider as the first visible layer of first RBM and as the hidden layer first RBM, and (39) can be rewritten asThe probability that the RBM assigns to a visible vector is The contrastive divergence method is then used to obtain the update formula of DBN parameters [19]:
For the construction of BP network, a classical model is used. The relevant functions, including activation function, energy function, and update formulas, are then set based on the classical model of BP network. The structure diagram of DBN is presented in Figure 2.
The overall control flowchart is shown in Figure 3. Firstly, the function of control output voltage is obtained by minimizing energy function, which is the sufficient condition of optimal control, and the MPC model is constructed based on it. The DBN is designed according to MPC model to obtain the parameters of MPC. The training of DBN is achieved with the use of robotic grinding historical data. After identification of MPC parameters, the MPC is used to perform online closeloop control with the use of robotic grinding force deviation and system information , , .”
4. DBN Training and Simulation
Simulations are conducted to evaluate feasibility and performance of the DBN used to obtain the model predictive control parameters of robotic grinding. A BP network is constructed, and a comparison between the BP network and the proposed DBN is conducted.
According to (35) and (37), the MPC is used to obtain with the use of , which means that the MPC represents the relationship between and . Therefore, the training data is required to contain the relationship between and .
The training data of DBN is determined according to the following aspects.
(1) The parameters of training data are determined according to influence factors of robotic grinding process, including grinding depth, grinding tool feed rate, tool rotation speed, and control output voltage. As there is a functional relationship between feed rate and control output voltage when a controller is applied, this paper uses feed rate change to represent change of control voltage. The ranges of the parameters are grinding depth is from 0.6 mm to 1.2 mm; initial feed rate is from 0.5 mm/s to 1.2mm/s; feed rate change is from 0.1 /s to 1 /s; tool rotation speed is 4000 r/min. (2) Since grinding is a material removal process, it is difficult and costly to perform dozens of experiments and obtain a large amount of training data; the DBN model used in this paper is only designed for specific case (robotic grinding) as Oh and Jung did in their research [19], which can reduce the quantity demand of data for DBN training but also reduce the applicability of DBN applying in other robotic operations. Based on this, the scale and the size of learning sample are determined. The scale of learning sample batch is 40 packets and each packet has about 50004 data. 34 packets are randomly selected for training and the remaining 6 are selected for model validation. The maximum epoch is 300.
The training performances of both networks are shown in Figure 4. For the model training of robotic grinding status at time , the DBN realizes fitting at 31st epoch and BP network realizes fitting at 52nd epoch. The model training of robotic grinding status from time to is also conducted, and the training performance of both BP network and DBN is presented in Table 1. In overall, the training performances of the proposed DBN are approximate to the training performances of the classical BP network. The training results are then used to construct.

(a)
(b)
(c)
(d)
A simulation is conducted to evaluate the predictive performance of models obtained, respectively, by the proposed DBN and the BP network. The steps of the simulation are as follows. Firstly, a sample is selected randomly, and the first 10 data of the sample are taken as inputs. Then, predictions of robotic grinding state are made according to these inputs. As shown in Figure 5, the prediction results of two models are basically consistent with the actual grinding force. The average prediction deviation of grinding force is about 0.5N. The maximum prediction deviation, about 2 N, appears at the peak of grinding force during cutinto state. At 2nd second, the grinding force predicted by the BP based model changes dramatically, while the one obtained by the DBN based model is smoother. Based on the above simulations, a preconclusion can be made that the performance of predictive model obtained by DBN network is better than the one obtained by BP network.
5. Robotic Grinding Control Experiments
5.1. Robotic Grinding System
The robotic machining system consists of YASKAWA industrial robot MH24, sixaxis force sensor MEFK6D40, and grinding tool. As shown in Figure 6, the force sensor and grinding tool are installed on the endeffector of robot, while the workpiece is installed on a fixed platform. The spindle of grinding tool is perpendicular to the Z direction of industrial robot sixth axis. The world coordinate system is set according to the robot Cartesian coordinate. The relevant parameters of equipment are the sixaxis force sensor is MEFK6d40, Germany; the workpiece material is Q235 steel; the grinding tool is composed of motor, handle, and grinding head.
The force signals are collected by the force sensor and delivered to the embedded realtime control system in PC with the use of Ethercat protocol. The analogue filter frequency of force sensor is 2500 Hz, and the sampling frequency of embedded realtime control system is 1 ms. The control output is then calculated by embedded realtime control system and sent to YASKAWA robot control cabinet to modify output pulse and adjust the feed rate of endeffector. The frequency of output voltage of control system is about 100 ms.
5.2. Robotic Grinding Control Experiments
To evaluate the performance of proposed control approach, robotic grinding control experiments are conducted. The open loop robotic grinding experiment, robotic grinding experiment with fuzzyPD control, and robotic grinding experiment with model predictive control are conducted, and comparisons between these experiments are made. The robotic grinding control experiments are carried out on the steel plate plane (Q235). The rotation speed of grinding tool is 4000 r/min, and the grinding depth is 1.2 mm; the initial feed rate is 1 mm/s; the material of cutter is D1614M06.
The grinding path of the experiment is shown in Figure 7. The grinding tool moves from point A to D via point B and C. The grinding path AB is perpendicular to workpiece surface, while the path BC is tangential to workpiece surface. For the convenience of discussion, in subsequent articles, AB is regarded as cutinto state, while BC is regarded as stable grinding state.
5.2.1. Experiment with Open Loop Control and Experiment with FuzzyPD Control
(A) Experiment with Open Loop Control. The desired grinding force for the experiment is 14 N. The experiment result is shown in Figure 8. The grinding force increases dramatically to 29N during cutinto state and then decreases around 15N. The dramatically increase of grinding force, which is labeled as abrupt change, during cutinto state can be illustrated by (10) and (16). When grinding tool cuts into the workpiece, the gap between feed rate and removal rate results in a deformation and leads to an increase grinding force. After perpendicular move, the grinding tool starts to move in tangential direction and the value of removal rate in perpendicular direction is larger than the one of feed rate resulting in a decrease of grinding force.
(B) Experiment with FuzzyPD Control. To verify the effectiveness of proposed control approach, an experiment based on fuzzyPID control is conducted. The inputs of fuzzyPD controller are force deviation and its firstorder derivative and control parameters and . Their domains are , , , and . The fuzzy reasoning is “IF A1 and B1 Then C1 and D1”, while the Gaussian function is used as fuzzy membership function.
The result of grinding experiment with PD control is shown in Figure 9. The grinding force increases considerably and reaches a peak value of 25 N, which is 4 N lower than the one of open loop control. The corresponding feed rate increases to 1.4mm/s at the first time and then decrease to 0.7mm/s with a valley of 0.4mm/s. During this process, a fluctuation of feed rate appears at second and the fluctuation amplitude of feed rate is around 0.1mm/s. This is because when the grinding force is close to the target force, the control parameters swift according to the fuzzy rules which results in a fluctuation of control output. This is a common problem of fuzzy control as fuzzy rules are often defined according to the artificial experience which is easy to result in control oscillations caused by switching of control rate. In view of this, a model predictive control approach is applied to robotic grinding control.
5.2.2. Robotic Grinding Experiment with Model Predictive Control
The robotic grinding experiment with model predictive control is then conducted to reduce the grinding force deviation and feed rate fluctuation that occurs in cutinto state. The prediction time of model predictive control is set as , and the experiment result is shown in Figure 10. The peak value of robotic grinding force is 18N, which is 7N lower and the corresponding feed rate curve changes from 1.4mm/s to 0.7mm/s with a valley of 0.3mm/s. The changes of robotic grinding force and feed rate shown in Figure 10 are smoother compared with the ones shown in Figure 9. This is because that the model predictive control approach can predict the future grinding deviation based on acquired information and result in a control compensation to reduce the coming up force deviation and feed rate fluctuation.
To further verify the performance of proposed control approach, another robotic grinding experiment based on model predictive control is conducted of which the prediction time is set as . As shown in Figure 11, the grinding force increases smoothly from 0N to 14N and the fluctuation of grinding force is around 2N. The abrupt change of grinding force, which occurs during cutinto state as shown in Figures 10 and 9, is completely eliminated. The corresponding feed rate changes from 1 mm/s to 0.75 mm/s without an abrupt increase as shown in Figures 9 and 10. This is because when the prediction time is adjusted to , the control system can have sufficient time to adjust feed rate smoothly, which can lead to sufficient control compensation.
To explore the flexibility of the proposed control approach, the grinding depth is changed to 1.4mm, while the grinding depth range of study sample is from 0.6mm to 1.2mm. The experiment result with open loop control is shown in Figure 12(a), while the experiment results with model predictive control are shown in Figures 12(b) and 12(c). As shown in Figure 12(a), the grinding force increases dramatically to about 33 N during cutinto state and then decreases to around 19N. When model predictive control approach is implemented, the grinding force increases smoothly from 0N to 19N without an abrupt change during cutinto state. The corresponding feed rate curve shown in Figure 12(c) is also smooth, which decreases from 1 mm/s to 0.7mm/s with a small valley of 0.65 mm/s at second.
(a)
(b)
(c)
Compared with fuzzyPD control approach, the proposed MPC approach can avoid mutation of grinding force and reduce the grinding force deviation when prediction interval is large enough, as shown in Figure 12. Also the feed rate change of robotic grinding based on MPC is smoother compared with the one based on fuzzyPD control. The conclusion is the same while the comparison is made with the experiments based on adaptive PID approach which presented in a previous research of writers [20]. The comparison shows that adaptive PID control can reduce the mutation of grinding force but not able to avoid it, as the adaptive PID control reacts only when the mutation comes out.
6. Conclusion
In this paper, a model predictive control approach based on deep belief network is proposed to control the deviation of grinding force caused by robotic deformation. The robotic grinding dynamic impacted by rigidflexible coupling effect is analyzed and the deep belief network is used to study the robotic grinding dynamic to acquire predict model. Based on this, the model predictive control approach can predict the change of grinding force deviation and achieve a combination of feedforward and feedback control.
The experimental results show that the prediction model of robot grinding based on the deep belief network can accurately predict the change of robot grinding state in time and in time . Compared with open loop control and fuzzyPD control, the model predictive control approach can perform a compensate control to beforehand reduce the force deviation. The elimination of the abrupt change of grinding force caused by deformation as shown in Figures 11 and 12 verifies the performance of proposed control approach.
Since current researches are mostly focusing on the feedback control of robot machining, while feedforward compensation control is rarely applied to robotic machining, the contributions of this paper are as follows: (1) the model predictive control approach is adjusted according to robotic grinding model within the consideration of coupling effect between rigidflexible system and control system. The MPC is used to predict the mutation of grinding force and perform forward control to avoid grinding force mutation. The proposed approach has both the advantages of feedforward control and feedback control, which can reduce the control oscillations and overflow caused by system delay and feedback control. (2) A deep belief network is designed according to the MPC model of robotic grinding control system. Since the MPC model is nonlinear and the parameters are difficult to obtain, the use of deep belief network can benefit the obtainment of MPC model parameters. Therefore, this study provides model reference and data support for the research of nonlinear control method of robot machining process based on intelligent model.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
References
 S. Mousavi, V. Gagnol, B. C. Bouzgarrou, and P. Ray, “Stability optimization in robotic milling through the control of functional redundancies,” Robotics and ComputerIntegrated Manufacturing, vol. 50, pp. 181–192, 2018. View at: Publisher Site  Google Scholar
 M. Cordes and W. Hintze, “Offline simulation of path deviation due to joint compliance and hysteresis for robot machining,” The International Journal of Advanced Manufacturing Technology, vol. 90, no. 14, pp. 1075–1083, 2017. View at: Publisher Site  Google Scholar
 A. Brunete, E. Gambao, J. Koskinen et al., “Hard material smallbatch industrial machining robot,” Robotics and ComputerIntegrated Manufacturing, vol. 54, pp. 185–199, 2018. View at: Publisher Site  Google Scholar
 R. Mei and C. Yu, “Adaptive neural output feedback control for uncertain robot manipulators with input saturation,” Complexity, vol. 2017, Article ID 7413642, 12 pages, 2017. View at: Publisher Site  Google Scholar
 N. Mendes and P. Neto, “Indirect adaptive fuzzy control for industrial robots: a solution for contact applications,” Expert Systems with Applications, vol. 42, no. 22, pp. 8929–8935, 2015. View at: Publisher Site  Google Scholar
 F.Y. Hsu and L.C. Fu, “Intelligent robot deburring using adaptive fuzzy hybrid position/force control,” IEEE Transactions on Robotics and Automation, vol. 16, no. 4, pp. 325–335, 2000. View at: Publisher Site  Google Scholar
 V. T. Yen, W. Y. Nan, and P. van Cuong, “Recurrent fuzzy wavelet neural networks based on robust adaptive sliding mode control for industrial robot manipulators,” Neural Computing and Applications, vol. 18, pp. 1–14, 2018. View at: Google Scholar
 T. Faulwasser, T. Weber, P. Zometa, and R. Findeisen, “Implementation of nonlinear model predictive pathfollowing control for an industrial robot,” IEEE Transactions on Control Systems Technology, vol. 25, no. 4, pp. 1505–1511, 2017. View at: Publisher Site  Google Scholar
 Y. G. Xi, D. L. Li, and S. Lin, “Model predictive control — status and challenges,” Acta Automatica Sinica, vol. 39, no. 3, pp. 222–236, 2013. View at: Google Scholar
 J. Wilson, M. Charest, and R. Dubay, “Nonlinear model predictive control schemes with application on a 2 link vertical robot manipulator,” Robotics and ComputerIntegrated Manufacturing, vol. 41, no. C, pp. 23–30, 2016. View at: Publisher Site  Google Scholar
 H. Peng, J. Wu, G. Inoussa, Q. Deng, and K. Nakano, “Nonlinear system modeling and predictive control using the RBF netsbased quasilinear ARX model,” Control Engineering Practice, vol. 17, no. 1, pp. 59–66, 2009. View at: Publisher Site  Google Scholar
 H. Peng, K. Nakano, and H. Shioya, “Nonlinear predictive control using neural netsbased local linearization arx model—stability and industrial application,” IEEE Transactions on Control Systems Technology, vol. 15, no. 1, pp. 130–143, 2007. View at: Publisher Site  Google Scholar
 N. Nikdel, P. Nikdel, M. A. Badamchizadeh, and I. Hassanzadeh, “Using neural network model predictive control for controlling shape memory alloybased manipulator,” IEEE Transactions on Industrial Electronics, vol. 61, no. 3, pp. 1394–1401, 2014. View at: Publisher Site  Google Scholar
 Z. Li, W. Yuan, Y. Chen, F. Ke, X. Chu, and C. L. P. Chen, “Neuraldynamic optimizationbased model predictive control for tracking and formation of nonholonomic multirobot systems,” IEEE Transactions on Neural Networks and Learning Systems, pp. 1–10, 2018. View at: Google Scholar
 X. Zeng, H. Peng, and J. Wu, “An improved multivariable RBFARX modelbased nonlinear model predictive control approach and application,” Journal of Central South University (Science and Technology), vol. 46, no. 10, pp. 3710–3717, 2015. View at: Google Scholar
 K. Dalamagkidis, K. P. Valavanis, and L. A. Piegl, “Nonlinear model predictive control with neural network optimization for autonomous autorotation of small unmanned helicopters,” IEEE Transactions on Control Systems Technology, vol. 19, no. 4, pp. 818–831, 2011. View at: Publisher Site  Google Scholar
 C. Zhang, Machinery Dynamics, Higher Education Press, Beijing, China, 2nd edition, 2007.
 H. Chen, Model_Predictive_Control, Science Press, Beijing, China, 1st edition, 2013.
 H. Oh, J. H. Jung, B. C. Jeon, and B. D. Youn, “Scalable and unsupervised feature engineering using vibrationimaging and deep learning for rotor system diagnosis,” IEEE Transactions on Industrial Electronics, vol. 65, no. 4, pp. 3539–3549, 2018. View at: Publisher Site  Google Scholar
 S. Chen and T. Zhang, “Force control approaches research for robotic machining based on particle swarm optimization and adaptive iteration algorithms,” Industrial Robot: An International Journal, vol. 45, no. 1, pp. 141–151, 2018. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2019 Shouyan Chen et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.