Published: 19 October 2018

Contact parameters identification of a cannon cradle and its bushing based on FDA and ELM

Hao Wen1
Baolin Hou2
1, 2School of Mechanical Engineering, Nanjing University of Science and Technology, Nanjing, China
Corresponding Author:
Baolin Hou
Views 22
Reads 10
Downloads 1194

Abstract

An identification method based on functional data analysis (FDA) and extreme learning machine (ELM) is presented to identify the contact parameters of a cannon cradle and its bushing. A virtual prototype of the cannon is built in ADAMS. The response curves of muzzle vertical acceleration with different contact parameters of the cradle and its bushing are obtained by simulation experiments and used for FDA as sample data. Features of the sample data are extracted by FDA and functional principle component analysis (FPCA), and the features and contact parameters are used to train the ELM. Simulation data and test data are used to verify the proposed method. The presented method is also proved to be feasible and effective by comparing actual muzzle vertical acceleration curve and the muzzle vertical acceleration curve from the virtual prototype with respect to the test data identification results.

1. Introduction

A cradle is the carrier of a cannon recoiling part, which rotates around the trunnion and transfers firing load to other carriages. Changes in contact parameters of the cradle and its bushing will affect mechanism kinematic accuracy during firing and firing accuracy of the cannon. To establish a virtual prototype of the cannon, it is significant to infer the contact parameters of the cradle and the front and back bushing. Due to a cannon is a complex and multi-parameter system, determination of key parameters is a core problem in the modeling process, but these parameters are difficult to measure and can only be obtained by identification.

Since it is difficult to find the analytic expressions of the measurable response to the parameters to be identified, optimization method based on analytic expression cannot be used to identify the contact parameters of the cradle and its bushing [1]. Machine learning (ML) can estimate the dependence of data based on the known sample, so ML can predict and judge the unknown data. Actual experiments can obtain high-quality samples, but for a cannon system, a large number of experiments require a lot of cost and resources and it is easier to obtain samples by virtual simulation.

In this paper, a virtual prototype of the cannon is established, and the response curves of muzzle vertical acceleration are obtained by sampling and simulation experiments of the parameters to be identified and used for FDA as sample data. Using FDA and FPCA to extract the sample data, the extracted features and parameters to be identified are used for ELM training as training samples. Finally, the simulation data and test data are used to identify the parameters and the presented method is proved to be feasible and effective by comparing real muzzle vertical acceleration curve and the muzzle vertical acceleration curve from the virtual prototype with respect to the test data identification results.

2. Modeling and simulation of virtual prototype

The recoiling part of a cannon on zero firing angle is mainly affected by recoil brake force, recuperator force, sealing device friction and cradle rail friction. According to the topological structure of the cannon, the dynamics model is established in ADAMS as shown in Fig. 1.

According to the research purpose, the contact parameters to be identified of the cradle and its bushing are the following five parameters: contact stiffness Kf and Kr, damping C, and clearance Df and Dr. Consulting engineering experience, the range of Kf is [100000, 200000] N⋅mm, the range of Kr is [150000, 300000] N⋅mm, the range of C is [10, 50], and the range of Df and Dr is [0.5, 1.5] mm. The five parameters are sampled within the distribution range, and virtual simulation experiments are carried out according to the sampling results. A total of 100 groups of muzzle vertical acceleration curves in counter-recoiling movement are obtained as shown in Fig. 2.

Fig. 1Virtual prototype of the cannon

Virtual prototype of the cannon

Fig. 2Sample curves of muzzle vertical acceleration

Sample curves of muzzle vertical acceleration

3. Features extraction based on FDA and FPCA

3.1. The FDA and FPCA algorithm

The basic idea of FDA is to treat the observed data as a whole, expressed as a smooth curve or a continuous function, and analyze it from a functional perspective. Assume that the ith observation sample contains a series of observation values yi1, yi2, ..., yin, it can be transformed to be a function xit by FDA, where t is the argument. Generally, the non-periodic data are expanded by the B-spline basis function system. Its functional form is represented by a linear combination of basis functions by letting c indicate the vector of length K and ϕ as the functional vector whose elements are the basis functions ϕk:

1
x=cΤϕ=ϕTc.

For making the estimated curves more stable, a roughness penalty PENm(x) is introduced:

2
PENmx=Dmxs2ds=DmcTϕs2ds=cTDmϕsDmϕTscds=cTRc,

where Dm is the m-order derivative and R=Dmϕ(s)DmϕT(s)ds.

Define H as a weighted matrix and λ as a smoothing parameter, we obtain:

3
PENSSEm(yc)=(y-Φc)TH(y-Φc)+λcTRc.

The expression for the estimated coefficient vector is:

4
c^=(ΦTHΦ+λR)-1ΦTHy.

FPCA is an expansion of Principle Component Analysis to Hilbert space. According to the derivation of Ramsay [4], the eigenfunction must satisfy the following equation:

5
v(s,t)ξ(t)dt=ρξs,

where ρ is the eigenvalue and v(s,t)=1Ni=1Nxi(s)xi(t) is the covariance function.

Let ξ indicate the eigenfunction and define an integral transform V, Vξ=v(,t)ξ(t)dt, then Eq. (5) can be expressed as:

6
Vξ=ρξ.

The principle component function also needs to be smoothed. Maximizing the variance of the samples with a roughness penalty:

7
PCAPSVξ=varξxidtξ2+λ×PEN2ξ.

For sample function, the form of basis function combination is xi(t)=k=1Kcikϕk(t). Let A be the covariance matrix of vectors ci and define W=ϕϕT, then Eq. (7) can be written as:

8
PCAPSV=bTWAWbbTWb+λbTRb.

The eigenequation corresponding to Eq. (8) is given by:

9
WAWb=ρW+λRb.

Now, preforming a Choleski factorization W+λR=LLT and defining S=L-1, the Eq. (9) can be written as:

10
SWAWSTLTb=ρLTb.

Defining u=LTb, the Eq. (10) can be written as:

11
SWAWSTu=ρu.

This is a eigenvalue problem, and u, b and the eigenfunction can be carried out in turns.

3.2. Features extraction of muzzle vertical acceleration curves

In this paper, a 4-order B-spline basis function and a 2-order roughness penalty function with smoothing coefficient λ= 50000 are used. Fig. 3(a) shows the sample data for the muzzle vertical acceleration after this functional procedure. Fig. 3(b) shows the first ten principal component functions.

To ensure that the sum of ratio of the principle component functions is greater than 90 %, we select the first 10 principal component functions. For every sample, after it is made to be functional data, the principal component scores can be acquired by calculating the inner product of the functional data and principal component functions, and the principal component scores are the eigenvalues we need.

4. Parameter identification

4.1. ELM algorithm

ELM is a typical single-hidden layer feedforward neural network [9]. Assume any independent sample xi,ti, xi=[xi1,xi2,...,xin]TRn, ti=[ti1,ti2,...,tim]TRm, mathematical model of ELM with n input nodes and m output nodes is:

12
i=1N~βigi(xj)=i=1N~βig(wixj+bi)=oj, j=1,2,...,N,

where N~ is the number of hidden layer nodes, g(x) is the activation function, wi is the weight connectiing the ith input node and hidden node, βi is the weight connecting the ith hidden node and output node, and bi is the threshold of the ith hidden node.

Fig. 3The results of FDA and FPCA

The results of FDA and FPCA

a) Muzzle vertical acceleration after function procedure

The results of FDA and FPCA

b) Proportion of first ten principal components

When the number of hidden nodes is equal to training samples, ELM can approximate the training sample with zero error. There are βi, wi, bi and output layer matrix H which make:

13
Hβ=T,H(w1,...,wN~,b1,...,bN~,x1,...,xN)=g(w1x1+b1)g(wN~x1+bN~)g(w1xN+b1)g(wN~xN+bN~)N×N~β=β1 β2 ... βN~N~×mT, T=t1 t2 ... tNN×mT.,

β can be acquired by solving the least squares solution of Eq. (13):

14
β^=H+T,

where H+ is the Moore-Penrose generalized inverse of the output matrix.

The identification accuracy and generalization performance of ELM are evaluated by average relative error E and determination coefficient R2, respectively:

15
E=y^i-yiyi, i=1,2,...,n, R2=li=1ly^iyi-i=1ly^ii=1lyi2[li=1ly^i2-(i=1ly^i)2][li=1lyi2-(i=1lyi)2].

4.2. Testing of muzzle vertical acceleration

The artificial recoil method is used to test the muzzle vertical acceleration by using the wireless acceleration node A104EX in the wireless test system. Fig. 4 shows the components of the wireless test system. The sensors are installed on the muzzle brake as shown in Fig. 5. The test is repeated 5 times at a firing angle of 0°, and the test results are shown in the Fig. 6.

Fig. 4Hardware components

Hardware components

Fig. 5The installation diagram of the sensors

The installation diagram of the sensors

Fig. 6Test curves of muzzle vertical acceleration

Test curves of muzzle vertical acceleration

4.3. Identification results and analysis

In this paper, 100 groups of sample curves are obtained, of which 95 groups are used for ELM training, and the other 5 groups are used to verify the identification effect of ELM. Table 1 shows the identification results of five contact parameters, and the first five times are the identification results of five groups of simulation data and the last is the identification results of test data. The average relative error E and determination coefficient R2 of the identification results are shown in Table 2. The identification results of Ef and Eb are better and the identification results of Kf and Kb can also meet the requirements. The identification results of C are relatively poor, and the probable cause is that C has less effect on the muzzle vertical acceleration, or the number of the samples used to train ELM is insufficient.

Table 1Identification results

No. of identification
C
Df (mm)
Dr (mm)
Kf (×105 N⋅mm)
Kr (×105 N⋅mm)
1
Real value
21.2857
0.8438
1.2438
1.7786
2.9587
Identification value
22.3214
0.8306
1.2275
1.7571
2.9852
2
Real value
17.6429
1.1844
1.1281
1.3293
2.9913
Identification value
17.6071
1.1975
1.1206
1.3607
2.8201
3
Real value
24.5357
0.5594
0.9234
1.3929
1.6229
Identification value
27.1286
0.5625
0.9313
1.2857
1.9583
4
Real value
21.6929
0.8825
0.6094
1.3536
2.9737
Identification value
24.2143
0.9281
0.5813
1.3821
2.8256
5
Real value
27.7143
0.9625
1.0063
1.6393
1.5703
Identification value
25.8714
96563
1.0156
1.7071
1.7862
6
Real value
\
\
\
\
\
Identification value
14.5357
1.3924
0.8263
2.2214
1.7903

The parameters identified by the test data are substituted into the virtual prototype model for simulation, and the test curve is compared with the simulated curve. As shown in Fig. 7, the two curves are more consistent, which proves the reliability of the identification results.

Table 2Average relative error and determination coefficient

Parameter
Average relative error
Determination coefficient
C
8.01 %
0.7681
Df
1.40 %
0.9902
Dr
1.47 %
0.9962
Kf
3.13 %
0.9035
Kr
9.12 %
0.9749

Fig. 7Comparison of identification curve and test curve

Comparison of identification curve and test curve

5. Conclusions

In this paper, a virtual prototype of a cannon is established in ADAMS. The muzzle vertical acceleration curves are obtained by the sampling of the parameters to be identified and simulation experiments and used for FDA as sample data. Using FDA and FPCA to extract the sample data, the extracted features and parameters to be identified are used for ELM training. The simulation data and test data are used to identify the contact parameters of the cradle and bushing. The presented method is proved to be feasible and effective by comparing real muzzle vertical acceleration curve and the muzzle vertical acceleration curve from the virtual prototype with respect to the test data identification results.

References

  • Zhao Q. Q., Hou B. L. Parameter identification of a shell transfer arm using FDA and optimized ELM. Chinese Journal of Engineering, Vol. 39, 4, p. 611-618.
  • Zong Weiwei, Huang Guang Bin Face recognition based in extreme learning machine. Neurocomputing, Vol. 74, Issue 16, 2011, p. 2541-2551.
  • Ramsay J. O. When the data are functions. Psychometrika, Vol. 47, Issue 4, 1982, p. 379.
  • Ramsay J. O., Hooker G., Graves S. Functional Data Analysis with R and MATLAB. Second Edition, Springer Science + Business Media Inc., New York, 2005.
  • Li M. Analysis of water data function based on principle component analysis method. Natural Sciences Edition, Journal of Hefei University, Vol. 24, Issue 4, 2014, p. 21-25.
  • Zipunnikov V., Caffo B., Yousem D. M., et al. Functional principal component model for high-dimensional brain imaging. NeuroImage, Vol. 58, Issue 4, 2011, p. 772-784.
  • Gao X. X., Sun H. G., Hou B. L., et al. Optimization design of positioning precision for howitzer shell transfer arm with parameter uncertainty. Acta Armamentarii, Vol. 35, Issue 6, 2014, p. 776-781.
  • Huang G. B., Zhao Q. Y., Siew C. K. Extreme learning machine: theory and applications. Neurocomputing, Vol. 70, Issue 1, 2006, p. 489-501.

About this article

Received
07 September 2018
Accepted
21 September 2018
Published
19 October 2018
SUBJECTS
Measurements in engineering
Keywords
parameter identification
functional data analysis
extreme learning machine
cradle