Deformation-compensated modeling of flexible material processing based on T-S fuzzy neural network and fuzzy clustering
Yaohua Deng^{1} , Sicheng Chen^{2} , Jiayuan Chen^{3} , Jiayi Wu^{4} , Qiaofen Zhang^{5} , Liming Wu^{6} , Hui Chen^{7}
^{1, 2, 3, 4, 5, 6, 7}Guangdong University of Technology, Guangzhou, Guangdong, China
^{1}Corresponding author
Journal of Vibroengineering, Vol. 16, Issue 3, 2014, p. 1455-1463.
Received 19 April 2014; received in revised form 23 April 2014; accepted 25 April 2014; published 15 May 2014
JVE Conferences
According to the factors that influence flexible material processing (FMP), the deformation compensation modeling method based on T-S fuzzy neural network is proposed. This method combines T-S fuzzy reasoning with a fuzzy neural network. Firstly, fuzzy clustering is introduced to extract fuzzy membership functions and the fitness of fuzzy rules of T-S fuzzy neural network antecedent from the past processing data. Secondly, with the steepest descent method, back-propagation iteration is used to calculate the connection weights of the network. The processing of experiments shows that T-S fuzzy neural network modeling is superior to typical T-S model. The angle error and the straightness error processed by NTS-FNN is 40.4 %, 28.8 % lower than those of STS-FNN. The minimum processing time processed by NTS-FNN is lower by 46.1 % than that of STS-FNN.
Keywords: flexible material path processing, deformation-compensated modeling, T-S fuzzy neural network, fuzzy clustering.
1. Introduction
The flexible material processing (FMP) refers to the process that takes various complex graphical processing on a flexible material (combined by multi-layer soft material), and an uneven three-dimensional pattern strayed on the surface of a flexible material [1, 2].
FMP deformation can be caused by multiple factors. Generally, several relevant variables will be regulated in order to reduce processing errors. Essentially, the compensation control of FMP deformation is a multi-input and a multi-output process. At present, the intelligent modeling of FMP is mainly based on the method of regression analysis, fuzzy reasoning and neural network [3, 4]. The method of Neural Network PLSR (NNPLSR) [5]; Kernel PLSR (KPLSR), which is suitable for real-time process modeling, is studied in [6-7]. Regarding the methods of modeling artificial intelligence, in 1993, a Japanese scholar Yasukawa proposed a simplified T-S fuzzy model. Its parameter of conclusion is expressed by a single value, which greatly simplifies the identification process [8]. In 1996, another Japanese scholar Sugeno proposed a T-S fuzzy model of a continuous recognition algorithm, the model parameter adjustment is achieved by the least squares method of weighted round robin [9]. In 2002, S. G. Cao from Australia proposed a fuzzy dynamic model, which consisted of multiple linear equations that smoothly connected global models through fuzzy membership functions. The accuracy of approximation can be arbitrary to any continuous nonlinear functions, which are defined on the compact set [10]. In 2011, Uros Zuperl from Slovenia studied the neural network model cutting force prediction of 3D spherical milling process [11].
Through the analysis of the literature above, it can be concluded that fuzzy, neural network modeling method has a strong non-linear expression. It is suitable for dynamic and complex uncertainty process modeling. A primitive angle, the interpolation rate, the processing direction angle and the flexible parts clamping way influence FMP deformation. Furthermore, some problems exist in the actual machining process, which requires deformation compensation prediction in real-time.
This paper proposes the FMP deformation compensation modeling, which effectively combines T-S fuzzy reasoning and a fuzzy neural network. It also introduces fuzzy clustering in order to achieve antecedent network model identification [12]. Based on fuzzy clustering, it adaptively adjusts the input space clustering center, the radius and the number of clusters, completes the input space fuzzy classification, determines the membership function and the rule fitness of data points, and improves the model training speed and approximation accuracy.
2. FMP deformation compensation modeling based on T-S fuzzy neural network
As the consequent network structure of a typical T-S model is relatively simple, if a new hidden layer is added to the consequent network, it will contribute to improve the universal approximation capacity of the model. Fig. 1 shows a new T-S fuzzy neural network extended model diagram of the standard T-S model. An antecedent network and a consequent network construct the new T-S fuzzy neural network model. The antecedent network is used to match the premise of T-S fuzzy rules. Meanwhile the consequent network generates the fuzzy rules. Then, the non-linear relationship modeling between the compensation output $S=\{{s}_{1},{s}_{2},\dots ,{s}_{r}\}$ and the deformation factors $X=\{{x}_{1},{x}_{2},\dots ,{x}_{n}\}$ is established through weighted calculations of antecedent and consequent networks.
Fig. 1. Construction of T-S fuzzy neural network model after extension
2.1. FMP deformation compensation model constructed by T-S fuzzy neural network
In Fig. 1, the T-S fuzzy neural network structure of FMP deformation compensation includes the function of a structuring membership, the rules of fitness of an antecedent network, and the calculations of the weight parameter of a consequent network [4].
2.1.1. Constructed antecedent network through clustering division of input data
The purpose of the division of FMP input data (FMP deformation influence factors) is to divide the input data fuzzy level, and to identify the coordinates of the center of each clustering and the width of each division district. Introducing the front part of the antecedent of T-S fuzzy neural network extracts input membership functions and the fitness of the fuzzy rules. Fuzzy clustering of the input data is mainly for working out the fuzzy partition matrix and clustering centers. Fuzzy clustering distance measurements are shown by Mahalanobis distance by calculating the minimum objective function (Eq. (1)) of $c$ groups of fuzzy clustering division. This is done in order to obtain the data point ${x}_{k}$’s clustering centers and the fuzzy partition matrix:
In Eq. (1): $X=\{{x}_{1},{x}_{2},...,{x}_{n}\}$ is the $n$-dimensional input data set. $U=\left[{u}_{ik}\right]$ is the fuzzy partition of $X\left({u}_{k}\right)$ is a $k$-th sample that about $i$-th class membership degree), and $\sum _{i=1}^{c}{u}_{ik}=\text{1}$, $k=$1, 2,…, $N$. ${D}_{ik{M}_{i}}^{2}$ is the square inner product from the mid-point of $n$-dimensional data space ${x}_{k}$ to clustering center ${v}_{i}$ (Eq. (2)):
In Eq. (2), ${M}_{i}=\mathrm{d}\mathrm{e}\mathrm{t}{\left({F}_{i}\right)}^{\frac{1}{n}}{F}_{i}^{-1}$, where ${F}_{i}$ is a positive definite symmetric matrix:
In Eq. (1), the minimum objective function is the basis for input data fuzzy clustering. Meanwhile, $\sum _{i=1}^{c}{u}_{ik}=\text{1}$, $k=$1, 2,…, $N\text{.}$ Eq. (1) use Lagrange multiplier method to get the objective function and constrains. Then, introduced Lagrange multiplier $\lambda $:
The requirement for getting the extreme of $\stackrel{-}{L}\left(X;U,V,\lambda \right)$ is:
Therefore, if the values of data set $X$, the categories number of clustering $c$, and the fuzzy weighting exponent $\beta $ are known, the data partitioning, the optimal fuzzy classification matrix and the cluster center all can be calculated by Eq. (5) [13].
Suppose the fuzzy partition matrix of input data is $U=\left[{u}_{ik}\right]$ after the fuzzy clustering division, the fuzzy category of $c$ set is ${G}_{i}$ (1$\le i\le c$), then categories’ center of ${G}_{i}\text{,}$${v}_{iq}\text{,}$ and the corresponding variance ${\sigma}_{iq}^{2}$ can be:
If the input data space division requires a high fuzzy clustering, in each ${G}_{i}$ category data ${q}_{k}$ is very close to the category component of ${v}_{i,n+1}$ ($i=\text{1,}\text{}\text{2,...,}\text{}c$) of the clustering center ${V}_{i}={[{v}_{i,1},{v}_{i,2},...,{v}_{i,n+1}]}^{T}$ (variance is ${\sigma}_{iq}^{2}\approx \text{0}$), the clustering center corresponds to ${G}_{i}$, then the category value ${q}_{k}$, which is the shortest distance from the categories’ center ${v}_{iq}$ can be used as the decision function $D\left({G}_{i}\right)$:
Corresponding, relative to ${G}_{i}$ membership function of ${x}_{k}$ is presented in Eq. (8), where $i=\text{1,}\text{}\text{2,...,}\text{}c$, $j=\text{1,}\text{}\text{2,...,}\text{}n$, $|{v}_{iq}-{v}_{iq}^{\text{'}}|$ is the width of the corresponding input data division region. ${v}_{iq}^{\text{'}}$ is the value of the cluster center, which is closest to the center of the $i$-th cluster. $\gamma $ reflects the coefficient of the membership reduction speed when the input sample is away from the cluster center:
The physical meaning of the membership function ${Z}_{ji}\left({x}_{kj}\right)$ can be described as a measuring of a similar relationship between the input sample ${x}_{k}$ and the fuzzy category ${G}_{i}$’s prototype. If ${x}_{k}$ is away from the prototype, ${Z}_{ji}\left({x}_{kj}\right)$ is closer to 0. If ${x}_{k}$ is closer to prototype, ${Z}_{ji}\left({x}_{kj}\right)$ moves near 1 [4].
For each fuzzy category ${G}_{i}$ and its decision category, ${R}_{i}$ reflects one of the decision rules (${R}_{i}:$ if ${x}_{k}\in {G}_{i}$, then $d\left({x}_{k}\right)=D\left({G}_{i}\right)$, $i=\text{1,}\text{}\text{2,...,}\text{}c$). If the multidimensional fuzzy set ${G}_{i}$ is projected to the whole input data space, then the decision rule ${R}_{i}^{\text{'}}$ is as follows:
Respectively, the rule fitness of ${x}_{k}$, relative to the rules of ${R}_{i}^{\text{'}}$, is expressed by Eq. (10):
Eqs. (8) and (10) are the calculation formula of the membership function and rule fitness when constructing the T-S fuzzy neural network model. Therefore the constructing of T-S fuzzy neural network model follows the below steps: 1) The execution of the fuzzy clustering algorithm in input data space to get category ${G}_{i}$ ($\text{1}\le i\le c$); 2) The category center ${v}_{iq}$, corresponding to the variance ${\sigma}_{iq}^{2}$ of all fuzzy categories is calculated (by using Eq. (6)); 3) By using the quality indicators parameter $\xi $, an excellent divided fuzzy category is selected; 4) The fuzzy category decision function $D\left({G}_{i}\right)$ after a filter is calculated (by using Eq. (7)); 5) The membership ${Z}_{ji}\left({x}_{kj}\right)$ of sample data ${x}_{k}$ with respect to the fuzzy category ${G}_{i}$ is calculated; 6) Fuzzy rule sets are generated according to Eq. (9); rule fitness ${\upsilon}_{i}$ are computed according to Eq. (10).
$\xi $ is introduced as the parameter of the quality standard. Firstly, it assigns one small value to $c$, calculating all variance of fuzzy division categories. If ${\sigma}_{iq}^{2}>\xi $, $c=c+1$, then the fuzzy partition is redone and looped until the condition ${\sigma}_{iq}^{2}\le \xi $ is met. Only then the corresponding fuzzy categories are accepted to construct a membership of a fuzzy neural network model.
2.1.2. Calculation of consequent network parameters
The value of the consequent network parameters ${w}_{ij}^{r}$ and ${w}_{jk}^{r}$ are calculated based on the steepest descent learning algorithm. Supposed, ${S}_{r}={\left[{s}_{1},{s}_{2},\dots ,{s}_{r}\right]}^{T}\text{,}$${S}_{r}^{\text{'}}={\left[{s}_{1}^{\text{'}},{s}_{2}^{\text{'}},\dots ,{s}_{r}^{\text{'}}\right]}^{T}$ represent the actual and the expected outputs of the consequent neural network. Then, the network output minimum variance of error cost function is ${E}_{mse}=\frac{1}{2}\sum _{r=1}^{p}{\left({s}_{r}^{\text{'}}-{s}_{r}\right)}^{2}$.
Moreover, supposed $\mathrm{\Delta}{E}_{mse}=\frac{\partial {E}_{mse}}{\partial w}$, $W$ represents the connection weights between nodes, $u\in \text{(0,1)}$ is learning rate, neurons coefficient iterative function each layer of consequent neural network is $W(t+1)=W\left(t\right)-u\mathrm{\Delta}{E}_{mse}$, and calculated by the error back-propagation algorithm, the calculation formulas of ${w}_{jk}^{r}$ and ${w}_{ij}^{r}$ are:
3. Processing test experiments
Real measure data of FMP is taken to test the performance of the new deformation compensation model T-S fuzzy neural network (below called NTS-FNN), which was discussed above. The analysis and the calculation is carried out with Intel i5 processor, 2.4 GHz frequency, 2 GB RAM, Microsoft Windows 7 Professional operating system, and MATLAB R2008 software platform.
3.1. Modeling of FMP deformation compensation
The primitive angle ${x}_{1}$ (°); the feed depth ${x}_{2}$ (mm); primitive types (circle or line) ${x}_{3}$ (circle: 1, line: 2); the processing step ${L}_{step}$ (mm) ${x}_{4}$ as the deformation influence factors (which are the input parameters of NTS-FNN) are taken, and ${s}_{1}$, ${s}_{2}$ are taken ass the feed compensation parameters of the $X$, $Y$ direction (which are the output of NTS-FNN) in order to construct NTS-FNN model. This mainly includes: 1) Construction of the antecedent network membership function ${Z}_{ji}\left({x}_{kj}\right)$, the fuzzy rule fitness ${\upsilon}_{i}$; 2) Calculation of ${w}_{jk}^{1}$,$\mathrm{}{w}_{ij}^{1}$,$\mathrm{}{w}_{jk}^{2}$, ${w}_{ij}^{2}$ of consequent network parameter.
Fig. 2. Experiment platform of quilting CNC
Choose five pieces of flexible parts, which length and width is 200 mm×120 mm and the material is polyurethane foam (the modulus of elasticity $E=$0.3652 MPa, Poisson’s ratio $\mu =$0.25). Flexible parts have different thickness. The track machining experiment was conducted on the quilting CNC platform (Fig. 2), and the track angle machining error was controlled within ±2.0 %. 265 samples were gained (flexible parts’ thickness (mm) is 5, 15, 20, 25, 35, the machining track is parallelogram, cure angle is 15°, 30°, 45°, 60°) as test data of NST-FNN.
Assigned $\beta =\text{2}$, $\epsilon =\text{0.00001,}$ maximum iteration number ${t}_{max}=$400, $c=$2, $\xi =\text{1.3}\text{.}$ According to Eq. (8), assigned $\gamma =\text{2}$, membership function ${Z}_{ji}\left({x}_{kj}\right)$ and the fuzzy rule fitness ${\upsilon}_{i}$ of ${x}_{1}$-$\mathrm{}{x}_{4}$ are gained. Furthermore, the number of nodes of each layer of the antecedent network is determined: from 1 to 4 is 4, 40, 10, 10. Finally, according to the Hecht-Nielsen method, it is confirmed that each layer neurons nodes of antecedent network is 4-40-10-10, consequent network is 5-11-10-2.
233 samples are chosen from the original 265 samples to train consequent network. 32 samples are left for the accuracy test, at the learning efficiency initial value $u=$0.4, the minimum expected error ${e}_{min}=$0.001. The training step is 600, the gain consequent network parameter is ${w}_{ij}^{2}$:
Furthermore,$\mathrm{}{w}_{ij}^{1}$, $\mathrm{}{w}_{jk}^{1}$, $\mathrm{}{w}_{jk}^{2}$ is calculated.
3.2. Processing experiments
The processing track angle error ${f}_{\alpha}$, the straightness error ${f}_{l}$, and the primitive minimum processing time ${t}_{p}$ are chosen as the machining performance test indicators. Moreover, performance indicators are compared to the standard-based fuzzy neural network prediction model STS-FNN [4].
Fig. 3. Processing track geometry
The pattern in Fig. 3 is used as the processing track, polyurethane sponge is used as a flexible material (the modulus of elasticity $E=$0.3652 MPa, Poisson’s ratio is 0.25), length, width and height of the flexible material is 150 mm×100 mm×15 mm. The processing experiment is without precompensation (NC), the feedback compensation control based on STS-FNN or NTS-FNN is conducted separately.
1) Process of errors ${f}_{\alpha}$, ${f}_{l}$
Under the control of STS-FNN and NTS-FNN respectively, 14 groups processing samples are received. After data treatment, 3 groups of samples are removed (those with the starting points and end points not closed). Finally, 11 groups samples are selected to respectively measure the processing track angle errors and the straightness errors.
Fig. 4 shows a comparative graph of processing errors of the processing track angle error ${f}_{\alpha}$, the straightness error ${f}_{l}$ (that without precompensation (NC), and with compensation by STS-FNN, NTS-FNN). Processing track angle error: ${f}_{a-NC}^{}$, ${f}_{a-STSFNN}^{}$, ${f}_{a-NTSFNN}^{}$; straightness error: ${f}_{l-NC}^{}$, ${f}_{l-STSFNN}^{}$, ${f}_{l-NTSFNN}^{}$.
Fig. 4. Comparative graph of processing errors
a)${f}_{a-NC}^{}$, ${f}_{a-STSFNN}^{}$, ${f}_{a-NTSFNN}^{}$
b)${f}_{l-NC}^{}$, ${f}_{l-STSFNN}^{}$, ${f}_{l-NTSFNN}^{}$
Furthermore, more calculations are made: the processing track angle error of the average without precompensation (NC), compensation by STS-FNN, NTS-FNN are: ${f}_{a-NC}^{\text{'}}=$5.09°, ${f}_{a-STSFNN}^{\text{'}}=\text{4.06\xb0,}$${f}_{a-NTSFNN}^{\text{'}}=\text{2.42\xb0;}$ average straightness errors are: ${f}_{l-NC}^{\text{'}}=$0.44 mm, ${f}_{l-STSFNN}^{\text{'}}=$0.34 mm, ${f}_{l-ATSFNN}^{\text{'}}=$0.25 mm.
2) Primitive minimum processing time ${t}_{p}$
From the testing results shown in Fig. 5, the time-consuming of completing track processing is as below.
The average primitive minimum processing time of process without pre-compensation (NC) ${t}_{p-NC}=$2.86 s; with feedback compensation control based on STS-FNN ${t}_{p-STSFNN}=$5.66 s; based on NTS-FNN ${t}_{p-ATS-FNN}=$3.05 s.
From the processing experiment results, we can conclude that the processing track on the flexible material of thickness is 15 mm, ${f}_{a-NTSFNN}^{\text{'}}$ is 40.4 %, 52.55 % smaller than ${f}_{a-STSFNN}^{\text{'}}$, ${f}_{a-NC}^{\text{'}}$; ${f}_{l-NTSFNN}^{\text{'}}$ decreased 28.8 %, 44.45 % relative to $\mathrm{}{f}_{l-STSFNN}^{\text{'}}$, ${f}_{l-NC}^{\text{'}}$; ${t}_{p-NTSFNN}$ decreased 46.1 % relative to ${t}_{p-STSFNN}$, and increased by 6.65 % more than ${t}_{p-NC}$.
Fig. 5. Comparative graph of primitive minimum processing time
4. Conclusion
1) In the T-S fuzzy neural network modeling process, pre-optimized input spaces (by using the fuzzy clustering method) ensure an antecedent network membership function, simplify the fuzzy rule number, help save fuzzy reasoning time, and improve model training speed. The testing result shows that the primitive minimum processing time processed by NTS-FNN is lower by 46.1 % than that processed by STS-FNN.
2) FMP deformation compensation modeling method includes the advantage of fuzzy clustering and T-S fuzzy neural network method. Compensated by NTS-FNN, an angle error, the straightness error is 40.4 %, which is 28.8 % lower than these of STS-FNN.
3) The consequent network and the antecedent network are relatively independent. Therefore they can easily realize embedded parallel computing. This is conducive to the application of the model in the prediction of the actual FMP compensation.
Acknowledgements
This work is partially supported by National Natural Science Foundation of China (No. 51205069), Natural Science Foundation of Guangdong Province China (No. S2013010013288) and Combination project of Industry-university research of Foshan Shunde of Guangdong Province China (No. 2013CXY02). Thank you for all your help.
References
- Yaohua Deng, Sicheng Chen, Bingjing Li, Jiayuan Chen, Liming Wu Study and testing of processing trajectory measurement method of flexible workpiece. Mathematical Problems in Engineering, Vol. 2013, 2013. [Search CrossRef]
- Yaohua Deng, Bingjing Li, Sicheng Chen, Jiyuan Chen Deformation forecast of flexible material process by spline finite element method and application. International Journal on Smart Sensing and Intelligent Systems, Vol. 6, Issue 1, 2013, p. 333-351. [Search CrossRef]
- Sun Zengqi, Xu Hongbin Fuzzy-neural network based on T-S model. Journal of Tsinghua University, Vol. 37, Issue 3, 1997, p. 76-80. [Search CrossRef]
- Deng Yaohua Liu Guixiong ATS-FNN-based modeling and simulation for compensation prediction of fwp machining deformation. Journal of South China University of Technology, Natural Science Edition, Vol. 40, Issue 3, 2012, p. 146-151. [Search CrossRef]
- J. M. Karthikeya Udayagiri V. R.,Taleb Moazzeni, et al. Detection algorithms for the Nano Nose. 19th International Conference on Systems Engineering, 2008, p. 399-404. [Search CrossRef]
- Kyungpil Kim, Jong-Min Lee, In-Beum Lee A novel multivariate regression approach based on kernel partial least squares with orthogonal signal correction. Chemometrics and Intelligent Laboratory System, Vol. 79, 2005, p. 22-30. [Search CrossRef]
- Bai Yifeng, Xiao Jian, Yu Long Kernel partial least-squares regression. International Joint Conference on Neural Networks, 2006, p. 1231-1238. [Search CrossRef]
- Sugeno M., Yasukawa T. A fuzzy logic based approach to qualitative modeling. IEEE Trans on Fuzzy Systems, Vol. 1, Issue 1, 1993, p. 7-31. [Search CrossRef]
- Dan A. Ralescu, Michio Sugeno Fuzzy integral representation. Fuzzy Sets and Systems, Vol. 84, 1996, p. 127-133. [Search CrossRef]
- Feng G., Cao S. G., Rees N. W. Stable adaptive control of fuzzy dynamic systems. Fuzzy Sets and Systems, Vol. 131, 2002, p. 217-224. [Search CrossRef]
- Zuperl U. F. Cus neural control strategy of constant cutting force system in end milling. Robotics and Computer Integrated Manufacturing, Vol. 27, Issue 3, 2011, p. 485-493. [Search CrossRef]
- Francisco de A. T. de Carvalho, Camilo P. Tenório Fuzzy K-means clustering algorithms for inter valued data based on adaptive quadratic distances. Fuzzy Sets and Systems, Vol. 161, Issue 23, 2010, p. 2978-2999. [Search CrossRef]
- Zhang Dong-B., Wang Yao-Nan, Huang Hui-Xian Fuzzy rough model based rough neural network modeling. Acta Automatica Sinica, Vol. 34, Issue 8, 2008, p. 1016-1023. [Search CrossRef]