Abstract
The main purpose of this paper is to identify the dynamic forces between the conical pick and the coal-seam. According to the theory of time domain method, the dynamic force identification problem of the system is established. The direct problem is described by Green kernel function method. The dynamic force is expressed by a series of functions superposed by impulses, and the dynamic response of the structure is expressed as a convolution integral form between the input dynamic force and the response of Green kernel function. Because of the ill-conditioned characteristics of the structure matrix and the influence of measurement noise in the process of dynamic force identification, it is difficult to deal with this problem by the usual numerical method. In present content, a novel improved Tikhonov regularization method is proposed to solve ill-posed problems. An engineering example shows that the proposed method is effective and can obtain stable approximate solutions to meet the engineering requirements.
1. Introduction
Various mechanical and electrical products are often subjected to a variety of dynamic loads, including vibration, impact, noise, and thermal environment [1]. It is important for the safety and accurate design of structures to accurately determine these loads. In many cases, it is difficult to directly measure the external loads of large structures such as missiles, aircraft, and offshore platforms under the action of wind waves or alternating excitations. Dynamic load identification is based on the dynamic response of the measured system and the known dynamic characteristics of the system to solve the dynamic loads on the structure. In practical engineering, it is very difficult to directly measure dynamic loads. Therefore, load identification is an indirect method to obtain the required load [2–4].
At present, load identification methods mainly include frequency domain method and time domain method [5–7]. Frequency domain method was put forward earlier. Its basic idea is to identify excitation spectrum by response spectrum. It is mainly realized by inversion of frequency response function between excitation and response. However, the ill-conditioned problem of coefficient matrix and singular value decomposition problem are often encountered in matrix inversion [8–10]. Based on the back analysis of the complex convolution relationship between load and response, the time domain method directly determines the time history of dynamic force, which has a good application prospect in engineering. For instance, [11] adopted the singular value decomposition (SVD) method to predict hourly load and peak load for the next selected time span. Reference [12] presented a new approach for solving specific classes of inverse source identification problems. Reference [13] used a modified inverse patch transfer function (IPTF) method to reconstruct the normal velocities of the target source in a noisy environment, and [14] used neural network in combination with genetic programming to implement the load identification of electric equipment. Reference [15] proposed the finite element and wavelet-based method for reconstructing the moving force. In [16], a new and effective algorithm for optimization problems, GPSA, is utilized in civil engineering. Reference [17] presented a single-point method for the identification of prevailing disturbing loads in power systems. But we may encounter some difficulties. It is very difficult to use the above methods for large-scale and complex load identification problems. Meanwhile, in practical engineering problems, it is necessary to identify the random dynamic force on pick and coal-seam structure. Unfortunately, they are complex inverse problems with inherent ill-posedness.
Fortunately, there are many popular regularization methods to deal with the ill-posed problem in practical engineering applications [18–22]. It is well known that the Tikhonov regularization method [23], as a traditional technique, has been widely used in various mathematical and engineering problems to deal with ill-posedness [24], for instance, parameter reconstruction, force source reconstruction problem, and so on [25–29]. Reference [30] applied the Tikhonov regularization method to reconstruct the dynamic force. Reference [31] demonstrated that the Tikhonov regularization method can solve the inverse problem in a stable manner despite the presence of noisy data. Reference [32] identified the shock load on an electroelastic bimorph disk using the Tikhonov regularization method. However, the Tikhonov regularization method is not completely perfect; there are numerous unavoidable limitations and disadvantages, which is described as follows: (a) the approximate solution provided by Tikhonov regularization method is too smooth; (b) the approximate solution provided by Tikhonov regularization method may lack some details that the desired real solution might possess. Therefore, it is very necessary to discuss other regularization methods to improve the limitations and disadvantages of the regularization method mentioned above.
Considering the limitations of the above methods in a certain case, we propose a novel improved Tikhonov regularization method to offer a stable solution for the ill-posed problems in practical engineering applications. Compared with the conventional methods, the uniqueness and improvement of the proposed method are obtained in this paper.
First and foremost, the paper is different from the conventional methods, and the uniqueness and improvement of this paper can be described as follows: (a) a novel improved Tikhonov regularization method can make the error of the regularization solution reach the asymptotic optimal order; (b) a novel improved Tikhonov regularization method has better convergence, and the detailed information of the solution is clearly enhanced; (c) the method we propose can effectively overcome the smoothness of the solution of the ill-posed problem; (d) the proposed method not only eliminates the influence of large singular values but also ensures the filtering of small singular values; (e) the small perturbation of right-end data will not affect the accuracy of inversion results, which fully reflects the good stability of the improved method.
In this study, our main purpose is to effectively identify the dynamic force between conical pick and coal-seam by using a novel improved Tikhonov regularization method. The rest of this study is structured as follows. In Section 2, the inverse problem model is established. In Section 3, a novel improved Tikhonov regularization method is proposed. In Section 4, the stability of the proposed method is proved. In Section 5, an engineering application example is discussed. Some valuable conclusions are described in Section 6.
2. Establishment of Inverse Problem Model
The inverse problems encountered in engineering can be expressed by the following [33, 34]:where y(t) is the response which can be displacement, velocity, acceleration, strain, etc. k(t) is the corresponding Green’s function that is the kernel of impulse response. x(t) is the identified dynamic force.
Using the rectangular formula, we can get the discrete form of (1). where denotes the time interval and it can be rewritten as follows:So Hence, (2) can be simply written asor equivalently,Where , and denote structural dynamic response, Green’s function matrix, and input force at time , respectively.
Because the condition number of K is very large and the response signal Y measured in practice always contains noise, a very small disturbance can result in a huge solution deviation. So dynamic force identification is a typical ill-conditioned problem, and it is inappropriate to directly invert K in (3).
The above characteristics of ill-conditioned problems do not mean that ill-conditioned problems are not solvable but that traditional linear algebraic methods cannot be directly applied to solve such problems.
Therefore, in order to make the solution meaningful, the approximate solution is usually obtained by regularization technique.
3. The Regularization Method for Solving Inverse Problem
3.1. Preliminaries
As we all know, we have the following lemma for regularized filter functions.
Lemma 1 (see [35]). Let and be Hilbert spaces, be a compact operator, be a singular system for the linear operator , and the function be defined: ; then we have the following properties:, there is a constant , such that and ;Operator is defined by the following formula:which is a regularization method, and .
The function with the above properties is called a regularized filter function of .
For the famous Tikhonov regularization method, the Tikhonov regularization solution is the minimum value of Tikhonov functional; a singular system with compact operators can be expressed as is the Tikhonov regularized filter function, and , , . Next, we will get a new regularization method by giving a new regularization filter function.
3.2. A Novel Improved Tikhonov Regularization Method
According to the regularization theory, the regularization operator can be constructed by regularization filter function , which provides a theoretical basis for establishing a new regularization method.
In this paper, a novel improved regularized filter function is constructed as follows:where , , and
The method proposed in this paper can not only ensure that the large singular values are not corrected but also filter the small singular values without affecting the accuracy of the solutions.
Theorem 2. The function defined by formula (8) is a kind of regularized filter function.
Proof. ① When , ;
when , we all know Therefore, we get② When , , such thatwhen , .
When , we haveWhen , .
Letting , we haveThen such thatHence, , we have .
And ; we can obtain③ When ,
It is known from Theorem 2which is a kind of regularized filter functions and its corresponding regularization operator is In order to express the advanced filtering operator of the proposed method more clearly, we give another representation of its filtering algorithm.
Substituting formula (13) into formula (14), we can get Introduce the singular value decomposition (SVD) [36]:in which and are orthogonal matrices and Substituting formula (27) into formula (26), we can get Then, we have in which Hence, we can haveComparing the Tikhonov regularization method with the improved Tikhonov regularization method, according to formula (31), we can see that, in part , change from to 1, , this eliminates the effect of large singular values.
When change from to , ; this not only guarantees the filtering of small singular values but also improves the approximation order of solutions.
4. Stability Analysis of the Proposed Method
Theorem 3. Suppose that the solution of satisfies , and ; if choosing the regularization parameter ( is a positive constant), then there are the following error estimates:
Proof. As we all know, the error estimates of true solution and regular solution can be expressed as follows:It is known from Theorem 2 that we have By using the singular system of operator K, then we havewhen and , such that , hence, we can obtain When and , then we haveAnd So we can obtainIf , thenHence, we can getFrom the above proofs, we can see that the method proposed in this paper can make the error of the regularization solution reach the asymptotically optimal order.
5. An Example of Engineering Application
In this section, the proposed method is used to an engineering example of the identification problem of dynamic force acting on conical pick and coal-seam structure.
5.1. The Establishment of Experimental System
The structure diagram of conical pick and coal-seam is made in Figure 1. The parameters are as follows: the coal-seam was 2000 mm × 800 mm × 1750 mm. The density of conical pick and coal-seam is 7800 kg/m3 and 1607.98 kg/m3, respectively. The elastic modulus of conical pick and coal-seam is 600 MPa and 322.09 MPa, respectively. The Poisson ratio of conical pick and coal-seam is set as 0.25 and 0.25, respectively [37].

Within our framework, the experimental system was shown in Figure 2. The experimental system includes test-bed, motor, reducer, coupling, cutting mechanism, platform, coal-seam structure, the force-measured device, force sensor, the signal amplifier, and the DaspV10 intelligent data acquisition and signal processing system [38, 39].

5.2. The Principle of Force Measurement
The cutting resistance of the pick when cutting the coal-seam is transmitted through the gear sleeve, and its size is measured by the force sensor at the back end. The force direction of the sensor is defined as the axial load Fz in accordance with the axis of the pick, and the force direction measured is defined as the radial load Fy perpendicular to the axis of the pick.
The force state of the pick in the test force measuring device is shown in Figure 3. Z is the cutting resistance, Y is the propulsion resistance, f is the friction resistance between the support structure and the pick sleeve, beta is the tangential installation angle of the pick, O is the support point of the pick sleeve, l1 is the distance between the tip of the tooth and the support point, and l2 is the distance between the sensor and the support point. According to Figure 3, the force balance and moment balance equations of the pick are obtained. Setting , so we have and are called friction coefficient and structural dimension coefficient, respectively.

Experimental conditions: the installation angle of the pick is , the cutting impedance of coal - seam is set to 180 kN/m, the speed of the cutting arm is set to 40.8 r/min, and the traction speed is set to 0.8 m/min. The measured axial and radial loads of the pick are shown in Figures 4(a) and 4(b). According to the load curves of Figures 4(a) and 4(b), the cutting force is converted from (19) at the sampling discrete points, as shown in Figure 4(c), and the measured displacement response is shown in Figure 5.

(a) Axial load

(b) Radial load

(c) Cutting force

Under the above experimental conditions, it can be seen that the cutting force obtained by conversion is similar to that of axial load, and the radial load has little effect on its variation. That is to say, the cutting force Z is proportional to the axial load Fz. Therefore, the measured axial load can reflect the size and variation of cutting resistance and can be approximately characterized by the measured axial load when analyzing the characteristics of cutting resistance.
5.3. Result Analysis
Adding a random distributed perturbation to displacement response data, we obtain where indicates the noise level of the measurement displacement response and the generates arrays of random numbers whose elements are normally distributed with mean 0, variance , and standard deviation . returns an array of random entries that is of the same size as [40]. The noise level can be measured in the sense of root mean square error (RMSE) according to This process is known as Morozov’s discrepancy principle [41, 42].
Relative error (RE)
In order to consider the effect of noise level, we set . All algorithmic programs are set in MATLAB 7.0 and Windows 7 operating system. The identified results of dynamic force are achieved by using different regularization methods in this paper, which are shown in Figures 6, 7, 8, and 9, respectively.

(a) Tikhonov regularization method

(b) Iterative Tikhonov regularization method
![(c) Ref. [33] method](https://static.hindawi.com/articles/mpe/volume-2019/6095184/figures/6095184.fig.006c.jpg)
(c) Ref. [33] method

(d) Present method

(a) Tikhonov regularization method

(b) Iterative Tikhonov regularization method
![(c) Ref. [33] method](https://static.hindawi.com/articles/mpe/volume-2019/6095184/figures/6095184.fig.007c.jpg)
(c) Ref. [33] method

(d) Present method

(a) Tikhonov regularization method

(b) Iterative Tikhonov regularization
![(c) Reference [33] method](https://static.hindawi.com/articles/mpe/volume-2019/6095184/figures/6095184.fig.008c.jpg)
(c) Reference [33] method

(d) Present method

(a) Tikhonov regularization method

(b) Iterative Tikhonov regularization method
![(c) Reference [33] method](https://static.hindawi.com/articles/mpe/volume-2019/6095184/figures/6095184.fig.009c.jpg)
(c) Reference [33] method

(d) Present method
Figures 6(a), 6(b), 6(c), and 6(d) illustrate that the identified force is obtained by adopting various different regularization methods under , respectively. Figures 7(a), 7(b), 7(c), and 7(d) show the identified force by adopting various different regularization methods under , respectively. The identified force can be seen from Figures 8(a), 8(b), 8(c), and 8(d), which gives that the identified force under , and Figures 9(a), 8(b), 8(c), and 9(d) offer the identified force by adopting various different methods under , respectively.
It can be seen from Figures 6, 7, 8, and 9 that the various different methods can identify the dynamic force. On the contrary, the effect of dynamic force identification is different.
To explore the quality of dynamic force identification by adopting various different methods, we give the value of RE and CC and iterative steps, which are shown in Tables 1, 2, 3, and 4. It is easy to see from Table 1 that RE and iterative steps of Figure 6(d) are smaller than that of Figures 6(a), 6(b), and 6(c), On the contrary, CC value presents opposite result. Table 2 shows that RE and iterative steps of Figures 7(a), 7(b), and 7(c) are larger than that of Figure 7(d), and the CC value presents opposite result. Table 3 shows that RE and iterative steps of Figure 8(d) are smaller than that of Figures 8(a), 8(b), and 8(c), but the CC value is the opposite. Table 4 offers that RE and iterative steps of Figure 9(d) are smaller than that of Figures 9(a), 9(b), and 9(c); nevertheless, the CC value presents opposite result.
Hence, compared with the other methods, we find from Tables 1–4 that the proposed method has the smallest RE, the largest CC, and the least iterative steps, so the optimal value is .
In addition, with the increase of , although RE and iterative steps increase, they do not affect the identification accuracy; that is to say, the small perturbation of the right-hand data will not affect the accuracy of the inversion results, which fully reflects that the present method has a good stability.
Through the analysis, we draw the following conclusions that the present method is successful to identify dynamic force acting on conical pick and coal-seam.
6. Conclusions
In this paper, a novel improved Tikhonov regularization method is proposed to solve the ill-posed problem. Under the current framework, a new regularized filter function is constructed by using the singular system of compact operator, and we prove the stability of the proposed method effectively.
An engineering example indicates that the proposed method has the smallest RE, the largest CC, and the least iterative steps and can obtain an efficient approximation of the actual force compared with the Tikhonov regularization method, iterative Tikhonov regularization method, and Reference [33] method. Meanwhile, the small perturbation of the right-hand data will not affect the accuracy of the inversion results, which fully reflects the present method has a good stability.
Consequently, we can give a new conclusion that the proposed method can offer a unified stable solution for solving the ill-posed problem. Furthermore, we are confident that the proposed method will provide a reference for the next research work.
Data Availability
No data were used to support this study.
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) and Research and Development of Scientific and Technological Achievements in Provincial Universities of Heilongjiang Education Department (Contract no. TSTAU-C2018005).