MPE Mathematical Problems in Engineering 1563-5147 1024-123X Hindawi 10.1155/2017/9125734 9125734 Research Article Identification of Random Dynamic Force Using an Improved Maximum Entropy Regularization Combined with a Novel Conjugate Gradient http://orcid.org/0000-0002-2026-9598 Ren ChunPing 1 Wang NengJian 1 Liu ChunSheng 2 Chen Xinkai 1 College of Mechanical and Electrical Engineering Harbin Engineering University Harbin 150001 China hrbeu.edu.cn 2 School of Mechanical Engineering Heilongjiang University of Science and Technology Harbin 150022 China usth.edu.cn 2017 30102017 2017 02 05 2017 07 08 2017 13 08 2017 30102017 2017 Copyright © 2017 ChunPing Ren 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.

We propose a novel mathematical algorithm to offer a solution for the inverse random dynamic force identification in practical engineering. Dealing with the random dynamic force identification problem using the proposed algorithm, an improved maximum entropy (IME) regularization technique is transformed into an unconstrained optimization problem, and a novel conjugate gradient (NCG) method was applied to solve the objective function, which was abbreviated as IME-NCG algorithm. The result of IME-NCG algorithm is compared with that of ME, ME-CG, ME-NCG, and IME-CG algorithm; it is found that IME-NCG algorithm is available for identifying the random dynamic force due to smaller root mean-square-error (RMSE), lower restoration time, and fewer iterative steps. Example of engineering application shows that L-curve method is introduced which is better than Generalized Cross Validation (GCV) method and is applied to select regularization parameter; thus the proposed algorithm can be helpful to alleviate the ill-conditioned problem in identification of dynamic force and to acquire an optimal solution of inverse problem in practical engineering.

National Natural Science Foundation of China 51674106 51274091
1. Introduction

The inverse problem, which is a typical mathematical problem, has been widely applied in many fields of practical engineering, such as vibration isolation, health monitoring, and fault diagnosis . The primary concern of mathematical workers and engineering technicians in the inverse problem derives from the development of the theory of solving the problems of universal existence and ill-conditioned problems in practical engineering applications in the past few years . Identification of dynamic force between pick and coal-seam has been a major concern in failure analysis of rotary machinery; once accurate force sources are obtained, it is possible to apply some advanced algorithms to analyze the characteristics of the dynamic force signals. However, we can not obtain the accurate dynamic force directly by the traditional methods due to the ill-posedness of systems . Aiming to reduce the above problems, it is necessary to obtain a method for solving ill-conditioned problem.

Accordingly, regularization techniques, which can treat the ill-conditioned problem, have been utilized in many fields . In the past few years, maximum entropy (ME) regularization techniques have been put forward one after another, which has been applied successfully in the wide areas of image reconstruction, signal processing, force identification problem, and so on , and the advantages of using ME regularization technique are as follows: firstly, some important information from the incomplete data can be extracted; then the probability distribution of the constraint is hidden in it; and, finally, entropy term can be changed according to different objects . Hence, ME regularization technique mentioned above can lay a firm foundation for the steady and accurate solutions of ill-conditioned problem.

However, ME regularization technique is not perfect, which has been made better in many forms, and the means of improvement are generally concentrated in the entropy function combined with some optimization algorithms to solve regularization solution . Mohammad-Djafari et al.  applied the improved entropy function combining the Lagrange method to process mass spectrometry data problems. Haffner et al.  discussed an improved maximum entropy modeling to handle with multiclass categorization problems using Newton iteration algorithm. In , by combining maximum entropy and Bayesian method, a new maximum entropy regularization is constructed to distribute force acting on boundary structure. Myrheim and Rue  presented a novel algorithm combining the ME method with the preconditioned conjugate gradient method to reconstruct image. In , the maximum entropy function and the conjugate gradient method were utilized, which can be used to solve absolute value equation problems. However, as for the improved entropy function mentioned, we consider that it is a numerical optimization algorithm, and its disadvantage is that the rate of convergence is not high, which is only suitable for simple objects, but the identified results are not ideal for complex objects. In fact, the entropy function can be redesigned according to the properties of different objects, and identification results of different objects can be achieved when we use different entropy functions combined with some different optimization algorithms .

In our study, a fast mathematical algorithm that can identify dynamic force between pick and coal-seam is determined. The technique we propose is new to our research field. Significant improvements are made in comparison with our previous work.

First, the force identification procedure is transformed into an unconstrained optimization problem due to being suboptimal for conventional regularization methods. In this work, a novel optimization mathematical algorithm, combining a novel improved maximum entropy (IME) regularization technique and a novel conjugate gradient (NCG) method (for short IME-NCG), is proposed to simplify the optimization process and improve the applicability of the algorithm.

Second, as a fast technique of IME-NCG, the optimized objective function is redesigned using the 1-norm estimation and the elongated entropy function in order to overcome the numerical instability.

Third, the optimal solution (stable solution) of inverse problem in practical engineering is studied. The technique proposed does not require any a priori information on the model for the random dynamic force identification. Only the displacement response measured by experiment is adequate for the identification algorithm.

Fourth, global convergence and stability of the solution for identified force model are proved to improve the precision and detailed information of random dynamic force identification.

Fifth, the L-curve criteria, combined with the IME regularization technique, are applied to select the suitable regularization parameter for dynamic force identification in complex vibration system.

Finally, the applicability of the proposed technique is effectively demonstrated through an engineering application example. It is found that the performance of random dynamic force identification using the combination of IME-NCG algorithm and L-curve is better than that of other previous works, which can effectively reduce the ill-posedness of random dynamic force identification.

In the present work, we aim at identifying random dynamic force between pick and coal-seam combining a novel improved maximum entropy (IME) regularization technique and a novel conjugate gradient (NCG) method. The structure of the paper is arranged as follows. In Section 2, the identification model of random dynamic force is established. In Section 3, the IME-NCG algorithm is proposed. In Section 4, we present the engineering application and results discussion. Some important conclusions are obtained in Section 5.

2. Model of Identified Force

An abridged general view of random dynamic force between pick and coal-seam structure is plotted in Figure 1. For deterministic MDOF structures, the kinetic equation can be expressed by the following form :(1)Myt+Cyt+Kyt=Sx,y,z,t,where M, C, and K represent the mass matrix, damping matrix, and stiffness matrix, respectively; S(x,y,z,t) denotes the dynamic force column vector in different directions; y(t),y(t), and y(t) are displacement, velocity, and acceleration response vector, respectively.

A schematic of dynamic forces between pick and coal-seam.

The displacement of the structure can be described as (2)yt=i=1nφiqit=ΦQ,where Φ is matrix of mode and Φ=φ1,φ2,,φn; Q is the displacement vector and time function in the generalized coordinates; Q=q1,q2,,qnT.

By means of (1) and (2), the dynamic equations represented by the physical coordinates are transformed into the modal coordinate system, and the n decoupled modal equations are obtained in the following form: (3)qi+2ξiωiqi+ωi2qi=FitMi,in which ξi is the i order modal damping ratio; Mi,Fi(t) are the generalized mass and the generalized force, respectively, and Mi=φiTMφi and Fi(t)=φiTS(t).

The initial velocity and initial displacement of the system are assumed to be x0 and x0, respectively. So qi in (2) can be expressed as follows: (4)qit=q1it+q2it.

It is well known that q1i(t) can be obtained by the homogeneous equation of (1), which can be described as follows :(5)q1it=e-ξiωitqi0tcosωidt+qi0+ξiωiqi0tωidsinωidt,where qi0 and qi0 denote the corresponding values in modal coordinates, which can be described as follows: (6)qi0=φiTMx0Mi,

q 2 i can be expressed in the following form due to the fact it has nothing to do with the initial conditions of the system:(7)q2it=1Miωid0tFiτe-ξiωit-τsinωidt-τdτ,

The displacement of the system can be expressed as follows:(8)y=y0+y,in which y0=i=1nφiq1i and y=i=1nφiq2i.

The system force is identified according to the known displacement; the displacement is represented as follows: (9)y=y-y0.

The following expressions can be obtained through (5) and (7) which are introduced into (9):(10)y=0ti=1nφiφiTMiωidSiτe-ξiωit-τsinωidt-τdτ.

In general, y(t) can be described as a convolution integral of the modal force and the corresponding Green kernel function: (11)yt=0tGt-τSτdτ.

Hence, Green’s kernel function can be described by comparing (10) and (11): (12)Gt=i=1nφiφiTMiωide-ξiωitsinωidt.Let y(t)=y1(t)ym(t)T, S(t)=S1(t)Sm(t)T, and G(t)=g1(t)gm(t)T then t (11) was discretized, and it is transformed into a matrix form that can be expressed as follows:(13)y1y2y3ym=g1000g2g100g3g2g100gmgm-1gm-20s1s2s3sm,or simply noted as(14)Y=GS.

Owing to the responses containing large noise, so (14) is expressed by the following form:(15)Yc=GSid+e,where Sid denotes the identified dynamic force, e is the unknown noise, and e=l·std(y(t))·rand, l denotes a parameter level, std() denotes the standard deviation, and rand denotes a random number which ranges from −1 to 1.

To identify Sid, we need to get Yc and G; moreover, it was found that (15) is a typical ill-conditioned problem. In some cases, small varieties in the right-hand side e of (15) may lead to dramatic varieties in the solution . In other words, they are very susceptible to error, in the inverse problem, this may directly lead to greater deviations of the identified results , and thus it is not realistic to solve (15) using traditional mathematical methods .

3. Algorithm

From the above analysis it can be seen that identification problem of dynamic force is ill-conditioned, whose solution is usually astaticism . Therefore, it is necessary to find a way to achieve a stable solution. The regularization techniques are good method to alleviate the ill-conditioned problem. In the past few years, maximum entropy (ME) regularization techniques have been put forward one after another, which have been applied successfully in the wide area of practical engineering problem . However, ME regularization technique is not very perfect; it has some disadvantages. Hence, in this work, an improved maximum entropy (IME) regularization technique is studied.

3.1. Maximum Entropy (ME) Regularization Method

It is generally known that maximum entropy (ME) regularization technique has been widely adopted in the past few years. The idea of the ME method is shown that the solving of (15) is transformed into that of an optimization problem :(16)minJSid=GSid-Yc2+λi=1nSilnSi,where λ represents the regularization parameter.

Let Ω(Si)=i=1nSiln(Si), and it is called the entropy item. Thus, (16) can be replaced by the following:(17)minJSid=GSid-Yc2+λΩSi.

3.2. An Improved Maximum Entropy (IME) Regularization Method

In some circumstances, such as complex identified objects, ME regularization technique is not very perfect, and it has some disadvantages of poor quality of identification. Hence, in this paper, the entropy item of ME regularization technique is redesigned using the 1-norm estimation and the elongated entropy function according to the ill-posedness characteristics of dynamic force identification, which was named as an improved maximum entropy (IME) regularization technique, was proposed to process the dynamic force.

In this section, a function is used to design the entropy item which is expressed as follows according to 1-norm estimation method:(18)ΩSi=i=1n1+ln1+Sipiln1+ln1+Sipi,where represents 1 norm and 0<pi<2.

And S can be expressed by  (19)SS2+ξ1/2.Hence, (18) is described by(20)ΩSii=1n1+ln1+Si2+ξpi/2ln1+ln1+Si2+ξpi/2,where ξ denotes the predetermined parameter.

Hence, a novel improved entropy item for dynamic force identification can be obtained, and it can be described as follows:(21)minJSid=GSid-Yc2+λi=1n1+ln1+Si2+ξpi/2ln1+ln1+Si2+ξpi/2.

It is obvious that the proposed entropy item is different with that of ME regularization technique; it mainly considers the 1-norm estimation and uses a better entropy item to make regular solutions more stable, thus which may be more suitable for the force identification of complex objects.

3.3. Selection of the Regularization Parameter

In this paper, we discuss L-curve method and Generalized Cross Validation (GCV) due to the validity and the classical character in choosing the regularization parameter [35, 42].

According to the regularization theory, the regularization parameter λ for an improved maximum entropy (IME) regularization plays an important role in applying the numerical experiment process. First, the L-curve criterion is introduced to select an optimal regularization parameter.

By reading the relevant literature, we know that the L-curve criterion is a very powerful way for selecting a suitable value of the regularization parameter for many meaningful engineering and mathematical problems.

It is well known that the L-curve is represented as a log-log plot (logGS-Yc2,logS2) of the norm of the residual versus the corresponding norm of a regularized solution.

The curve shows a typical “L” shape, and the optimal value for the regularization parameter is considered to be the one that corresponds to the corner of the curve .

Typically, the GCV function can be defined in the following form :(22)minGCVλ=GSλ-Yc22trI-GGλI2,where tr() denotes trace of matrix.

3.4. Solving of Objective Function

Equation (21) is considered as a class of unconstrained optimization problems and a series of optimization algorithms, which can treat such unconstrained optimization problem. In this work, a novel conjugate gradient (NCG) method was proposed to minimize (21). It has always been known that iterative form of conjugate gradient (CG) method are described as follows :(23)Sk+1=Sk+αkhk,where αk denotes the step length, which can be computed using standard Armijo search conditions(24)αk=maxρj,j=0,1,2,,fSk+αkhkfSk+δαkgkThk,where ρ0,1 and δ(0,0.5).

h k denotes the search direction(25)hk=-gk,ifk=0-gk+βkhk-1,ifk1,where gk denotes the gradient of J(Sk) at the point Sk, gk=J(Sk), and βk is conjugate parameter.

A novel search direction hknew based on [46, 47] PRP conjugate gradient method is determined as follows: (26)hknew=-gk,ifk=0-gk+βkPRPzkhk-1new,ifk1,where (27)zk=minθk,1,θk=c1gk-12Lkhk-1new2,Lk=maxgk-gk-1vk-1,Lk-1,c00,1,vk-1=αk-1hk-1new,L0>0,βkPRP=gkTgk-gk-1gk-12.

Algorithm 1.

Step 1. Given an initial point S0, c1(0,1), and L0>0, set h0=-g0 and k=0.

Step 2. If gk=0, then terminate; else go to the next step.

Step 3. αk is computed based on (24).

Step 4. Set Sk+1=Sk+αkhknew.

Step 5. hknew is computed based on (26).

Step 6. Set k=k+1 and go to Step 2.

Theorem 2.

If the search direction hknew is calculated by (26), and the step length αk meets the standard Armijo search conditions, for all k0, and then one has(28)gkThknew-1-c1gk2.

Proof.

When k=0, we can obtain (29)gkThknew=-gk2.Then Theorem 2 is proved.

When k1, because the step length αk is calculated by formula (24), αk1, and then we can obtain(30)gkThknew=-gk2+βkPRPminθk-1,1gkThknew-gk2+gkgk-gk-1gk-12θk-1gkThknew-gk2+Lkvk-1hknewgk-12θk-1gk2-gk2+Lkhknew2gk-12θk-1gk2=-1-c1gk2.Hence, the conclusion is established.

3.5. Global Convergence

Under the following assumptions, the global convergence of the algorithm is established.

Assumption A.

(i) The objective function is bounded in the following level:(31)Ω=sRnfsfs0.(ii) The gradient g is Lipschitz continuous; that is, there exists a constant L>0 such that, for any x,yΩ, we obtain(32)gx-gyLx-y.

Lemma 3.

Suppose that Assumption A holds, the step length αk is calculated by (24) and the search direction hknew is calculated by (26), and there is a constant c2>0, then(33)αk-gkThknewhknew2,(34)k=0gk4hknew2<.

Proof.

Set c2=min1-c1,ρ1-δ/L, because the step length αk is calculated by (24), then αk1.

When αk=1, according to Theorem 2, we can obtain the following:(35)gkThknew-1-c1gk2.Because -gkThknewgkhknew, then we obtain(36)hknew1-c1gk.Then one has(37)hknew21-c1gkhknew-1-c1gkThknew.We have(38)αk=1-1-c1gkThknewhknew2.Hence, (33) is proved.

When αk<1, we know that ρ-1αk can not meet (24) according to the definition of the step length, and then(39)fsk+ρ-1αkhknew>fsk+δρ-1αkgkThknew.According to the mean value theorem, we have(40)fsk+ρ-1αkhknew-fsk=ρ-1αkgsk+tkρ-1αkhknew2hknewρ-1αkgkThknew+Lρ-2αk2hknew2,where tk0,1.

According to the formula, we obtain(41)δρ-1αkgkThknewρ-1αkgkThknew+Lρ-2αk2hknew2.So we have(42)αk-ρ1-δgkThknewLhknew2c2gkThknewhknew2.Hence, (33) is proved.

We can see that fk is monotone bounded according to (24) and Assumption A, so we can obtain the following conclusion based on (24) and (33) and Lemma 3:(43)k=0kδc21-c12gk4hknew2k=0kδc2-gkThknew2hknew2k=0k-δαkgkThknewk=0kfk-fk+1<.Hence, (34) is proved.

Theorem 4.

Suppose that the step length αk and the search direction hknew are generated by IME-NCG algorithm, and then one has that(44)limkgk=0.

Proof.

Owing to the fact that the step length αk and the search direction hknew are generated by (24) and (26), respectively, then(45)hknewgk+βkminθk-1,1hk-1newgk+gkgk-gk-1gk-12c1gk-12hknew2Lkhknewgk+gkgk-gk-1gk-12c1gk-12sk-1hknew2gk-gk-1hknew=gk+c1αkgk1+c1gk.So we can obtain(46)gk4dk2gk4dk21+c12=gk21+c12.Hence, (34) is proved.

4. The Engineering Application 4.1. Experimental Set-Up

In this paper, a diagram of experimental set-up is organized in Figure 2. The size of the coal-seam structure specimen manufactured artificially was 2000mm×800mm×1750 mm, the pick is mounted on the rotary cutting arm, it can be driven by the reducer and torque, and its power driven by the motor is rated at 55 kW . The cutting bench achieves free forward and backward movement by the hydraulic pressure drive control system. System output responses can be measured by the corresponding sensors, then converted by a signal amplifier, and finally recorded by V10Dasp data vibration signal acquisition system [49, 50].

A diagram of experimental set-up.

4.2. The Engineering Application

All codes were written in MATLAB 7.0 and run on a HP with 2.0 GB RAM and Windows 7 operating system. Stop the iteration if criterion is defined as gk10-5 being satisfied or run time more than 500 seconds. And the step length αk is calculated by (24), and the parameters are obtained as follows: δ=0.45, ρ=0.12, L0=10-4, c1=1-10-6,pi=1.3, i=1,2,,n, and ξ=2×10-10. For performance analysis of dynamic force identification techniques used, the performance measurement metrics for dynamic force are described as follows: restoration time, iterative steps, and root mean-square-error (RMSE), where RMSE is defined by(47)RMSE=1Ni=1NSid-Sactual2,where Sactual is the actual measured random dynamic force and Sid is the identified random dynamic force.

In this section, the material parameters of experimental conditions are shown in Tables 1 and 2.

The material parameters.

Materials Compressive strength (MPa) Elastic modulus (MPa) Poisson ratio Density (kg/m3)
Coal-seam 2.68 322.09 0.25 1607.98
Pick 600 0.25 7800

The experimental conditions.

Installation angle (°) Maximum cutting thickness (mm) Hauling speed (r/min) Rotary arm speed (m/min)
35 15 0.6 25

It is generally known that direct numerical treatment fails to solve the present inverse problem in the engineering application. Therefore, the proposed method is to be applied to identify dynamic force. However, the key problem is how to select the suitable regularization parameter to obtain the optimal solution. Moreover, the suitable regularization parameter obtained by L-curve is really important for numerical experiment process.

First, the graphical illustration of L-curve was shown in Figure 3(a) according to L-curve criterion. Let it be noted that it is very possible to obtain the optimal solution by combining the L-curve criterion with the proposed regularization method.

The selection of the suitable regularization parameter from the L-curve.

The graphical illustration of L-curve

λ = 10−1

λ = 10−2

λ = 10−5

λ = 10−7

It is clearly seen from Figure 3(a) that when the regularization parameter is small, the graph shows very clearly that the norm of regularized solution S2 is large, but the norm of the corresponding residual GS-Yc2 is small; we can speculate that it is likely to be caused by measurement errors. However, when the regularization parameter is large, it is found that the norm of regularized solution S2 is a poor approximation; on the contrary, the norm of the corresponding residual GS-Yc2 is large.

The corner of the L-curve indicates this transition, owing to the fact it shows a compromise between the minimization of the norm of the residual GS-Yc2 and norm of the solution S2. The selection of the regularization parameter of IME regularization can be not completely guaranteed to be suitable for all systems with the ill-posedness. However, many examples show that the L-curve criterion is a powerful method for determining a suitable value of the regularization parameter for many significant engineering and mathematical problems.

In order to further discuss the influence of L-curve on numerical experiments, then make sure that the optimal regularization parameter is selected from the L-curve. According to the L-curve criterion, the value of different regularization parameters are selected from Figure 3(a) to find the optimal regularization parameter. The measured displacement response was shown in Figure 4. And the identified results were formed in Figures 3(b), 3(c), 3(d), and 3(e).

Measured displacement response.

Figures 3(d) and 3(e) show that the identified result is extremely unsatisfactory due to small regularization parameter. However, we can see from Figure 3(b) that even if when the regularization parameter is large, that is, beyond a certain threshold, the identified result is far away from what is normal. At the same time, the identified result of the value at the corner of the L-curve is shown in Figure 3(c). Figure 3(c) clearly shows that the most accurate results are obtained under the optimal regularization parameter (λ = 10−2) compared with other cases of regularization parameters. So we conclude that the values at the corner of the L-curve are selected as the most suitable regularization parameter.

As we can see from the above discussion, the regularization parameter was obtained by L-curve method, and λL-curve=0.01. Then identified dynamic force can be given based on measured displacement response using the identified modal model and the proposed algorithm, and the identified results were formed in Figure 5 by assembling the analysis results of ME, ME-CG, ME-NCG, IME-CG, and IME-NCG algorithm.

Identified results of different algorithms at λL-curve=0.01.

ME algorithm

ME-CG algorithm

ME-NCG algorithm

IME-CG algorithm

IME-NCG algorithm

It was obtained that the random dynamic force is easily identified with application of the above identification algorithms, and the identified random dynamic force is obviously observed using the regularization parameter selected by L-curve method as seen in Figures 5(a), 5(b), 5(c), 5(d), and 5(e), respectively.

Figure 5(a) shows that the identified result is obtained using ME algorithm, Figure 5(b) shows that the identified result is obtained using ME-CG algorithm, identified result using ME-NCG algorithm is obtained from Figure 5(c), Figure 5(d) shows that the identified result is obtained using IME-CG algorithm, and the identified result using IME-NCG algorithm is shown in Figure 5(e).

Moreover, it can be shown that different identification algorithms can commendably identify the random dynamic force. However, to some extent, the performance of random dynamic force identification has a certain difference.

We analyzed the different identification algorithms at λL-curve=0.01. And restoration time, iterative steps, and RMSE were used to evaluate the quality of dynamic force identification.

In Table 3, it was found that the proposed IME-NCG algorithm is associated with minimum RMSE, minimum restoration time, and minimum iterative steps in comparison to ME, ME-CG, ME-NCG, and IME-CG algorithm, which proves that the IME-NCG algorithm is superior to the other algorithms in the identification of random dynamic force at λL-curve=0.01, which is closer to the actual dynamic force.

The measurement metrics comparisons using L-curve method.

Method ME algorithm ME-CG algorithm ME-NCG algorithm IME-CG algorithm ME-NCG algorithm
Restoration time (s) 36.0152 25.4632 18.6788 10.1442 3.0421
Iterative steps 30 25 19 15 8
RMSE 0.4702 0.3775 0.1867 0.0828 0.0783

However, as is well known, different regularization parameters are obtained by different regularization parameter selection methods, which will lead to different identified results. In this paper, the GCV is applied to compute the regularization parameter compared with L-curve method.

Figure 6 shows that the regularization parameter was obtained by GCV method, when λGCV=0.006, GCV(λ) reaches the minimum, and the identified results were formed in Figure 7 by assembling the analysis results of different algorithms at λGCV=0.006.

Selection of regularization parameters using GCV method.

Identified results of different algorithms at λGCV=0.006.

ME algorithm

ME-CG algorithm

ME-NCG algorithm

IME-CG algorithm

IME-NCG algorithm

In Figures 7(a), 7(b), 7(c), 7(d), and 7(e), it can be shown that five different identification algorithms can identify the random dynamic force, and it was found that the proposed IME-NCG algorithm is associated with minimum RMSE, minimum restoration time, and minimum iterative steps in comparison to the other algorithm from Table 4; hence, we can draw a conclusion that the IME-NCG algorithm is closer to the actual force at λGCV=0.006, which shows that the IME-NCG algorithm is superior to the other algorithms whether the regularization parameter is selected by using L-curve or GCV method.

The measurement metrics comparisons using GCV method.

Method ME algorithm ME-CG algorithm ME-NCG algorithm IME-CG algorithm ME-NCG algorithm
Restoration time (s) 40.1442 28.3664 20.4112 17.1100 5.0021
Iterative steps 35 28 22 18 10
RMSE 0.6047 0.4241 0.2827 0.0919 0.0806

From Tables 3 and 4, we can see that RMSE, restoration time, and iteration steps using L-curve method are lower than those of using GCV method whether using IME-NCG algorithm or the other algorithms.

Hence, we can conclude that the IME-NCG algorithm combined with L-curve method is available to identify the dynamic force.

5. Conclusions

In this paper, the identification model for random dynamic force between pick and coal-seam is established through the relationship between the identified force vector and the displacement responses measured by machinery vibration system. Then the entropy item of maximum entropy (ME) regularization is redesigned using the 1-norm estimation method and the elongated entropy function according to the ill-conditioned characteristics of force identification, which was named as an improved maximum entropy (IME) regularization technique, was proposed to process the random dynamic force. Then the force identification problem is transformed into an unconstrained optimization problem and a novel conjugate gradient (NCG) method was applied to solve the objective function, which was called IME-NCG algorithm.

The results show that the proposed IME-NCG algorithm is associated with minimum RMSE, minimum restoration time, and minimum iterative steps in comparison to ME, ME-CG, ME-NCG, and IME-CG algorithm, and L-curve method is introduced which is better than GCV method and is applied to select the regularization parameter. Thus we can safely conclude that the combination of IME-NCG algorithm and L-curve method can reduce the ill-posedness of random dynamic force identification and achieve the optimal solution of inverse problem in practical engineering.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported by the Chinese National Natural Science Foundation (Contract nos. 51674106 and 51274091).