Mechanics & Industry
Volume 22, 2021
|Number of page(s)||8|
|Published online||08 March 2021|
On the way to fault detection method in moving load dynamics problem by modified recurrent neural networks approach
Department of Mechanical Engineering, Vardhaman College of Engineering, Hyderabad, India
2 Department of Mechanical Engineering, National Institute of Technology, Rourkela, India
Accepted: 30 January 2021
Parameters identification on structure subjected to moving load can be predicted by using the accurate and reliable data. The concepts of recurrent neural networks (RNNs) approach have been used in parameters (crack locations and severities) identifications in structure subjected to moving load in the present methodology. This methodology has incorporated the knowledge based Elman's recurrent neural networks (ERNNs) and Jordan's recurrent neural networks (JRNNs) jointly for the identification of parameters. This approach has been addressed as the inverse problem for predicting the locations and quantification of cracks in the structure in a supervised manner. The Levenberg-Marquardt's back propagation algorithm is implemented to train the proposed networks. To check the robustness of the present method, Numerical studies followed by Finite Element Analysis (FEA) and experimental verifications (Forward problems) are presented as a case study by considering a multi-cracked simply supported structure under a moving mass. The estimated crack locations and severities obtained from the proposed RNNs model converge well with those of FEA and experiments. From the demonstration of the case study, it concludes that the proposed analogy can identify and quantify the crack locations and severities effectively.
Key words: Crack locations / crack severities / L-M algorithm / ERNNs / JRNNs
© S.P. Jena and D.R. Parhi, Hosted by EDP Sciences 2021
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Vibration based fault detection approaches are global methods which have been considered to spot the positions and severities of cracks in structures. The fault detection methods have achieved a wide range of attentions from engineering and scientific communities due to the unexpected structural failure. A consistent and efficient fault detection approach is essential to preserve the protection and reliability of structures. The fault or crack detection approaches have been carried out by using various conventional and intelligent methods for the past few years. In contrast, attention in the literature survey was raised significantly over the recent years in neural networks (NNs) based fault detection approaches. The NNs are quietly used to improve the abilities of conventional methods. Numerous studies have been carried out by various researchers and engineers to spot the positions and severities of cracks in structure by using the vibration based dynamical characteristics of structural responses. A brief literature study on vibration-based condition monitoring is carried out here with more inclusive reviews.
Bagherahmadi and Seyedpoor  proposed a proficient damage index approach for multiple damages in structure using the strain energy and frequency response function matrix of the structure. They have considered the change in strain energy of the structure as damage index parameter. The concepts of natural frequencies and mode shapes analysis mechanism are applied by Chasalevris and Papadopoulos  for detection of multiple cracks on structure. Chang et al.  developed a structural health monitoring approach based on ANNs modeling to interpret the damage parameters using the modal properties of the structures. The proposed approach has been demonstrated successfully to a 7DOF building structure for damage identification. Duimitriu  proposed a fault detection approach on railway damper vehicle suspension system by considering the vertical bogie acceleration.
Eroglu and Tufekc  have developed a FEM model for an edge cracked Timoshenko straight beam. They have solved this problem by using genetic algorithm in an inverse manner. Hou et al.  have developed Genetic Algorithm (GA) based Optimal Sensor Placement (OSP) approach for damage recognition in structure with l1-regularization damage detection. They have conducted experiments for a cantilever beam and three-storey frame structures to check the convergence of the developed method. Hakim et al.  have developed a damage identification procedure in the domain ANNs by taking into account an I-beam structure. Jena et al. [8–11] have carried out numerous studies like theoretical, FEA and experimental to determine the responses of different types of structures subjected to transit mass. They have determined the structure's response at various damage and moving load parameters. In the later part, Jena and Parhi [12,13] have also developed damage detection procedures for structure under traverse load as inverse methods using the concepts of different types of rule-based recurrent neural networks. (ERNNs and JRNNs).
Mao and Lu  have explored an effective method to determine the critical speed and resonance condition of bridge structure subjected to moving trains with different speeds. Their method is very much effective to study the effects of resonance in accordance with resonance speed. Pala and Reis  have investigated the effects of inertial, centripetal and Coriolis's forces on the response of a cracked beam subjected to a moving load. Pathirage et al.  have developed an auto encoder based approach for structural damage identification approach in a steel frame structure. The concepts of NNs and deep learning approaches are applied here for pattern recognination problem. Yang et al.  have explored a method based on the combination of ANNs and dynamics strain mechanism to identify the damage parameters in moving load problem. Yang et al.  have developed a static algorithm by implementing the moving load induced response and bending theory mechanism to locate damage in structure. Wu et al.  have developed a 3-dimensional finite element model to determine the dynamic response of a concrete structure due to the excitation of a moving mass. They have calculated the stress and deflection of the structure at the critical positions of the moving load. Zhang et al.  presented a noble non-model based approach for crack detection in bridge structure under vehicular load based on the concept of phase trajectory alteration of multi-type induced vibration measurement. The proposed mechanism is very useful to recognize the shear failure in composite bridge structure.
So far from the literature studies, the authors have got the attentions that numerous works have been carried out for structural health monitoring problem using different techniques as forward and inverse methodologies. Even if many works are also carried using NNs based approach, but using the concepts of RNNs methodologies are few in the applications of mechanical engineering era like structural dynamics. Again as per the concern of authors' knowledge, the literatures are very scanty on damaged detection in structure under moving mass using modified RNNs approach. The novelty of the proposed work is to develop a structural damage detection procedure using the combined model of the modified ERNNs and JRNNs approaches. The results obtained from the modified RNNs (combined model of the modified ERNNs and JRNNs) model are corroborated with those of FEA and experimental works. The entire damage detection procedure has been carried out in a supervised manner. The numerical, FEA and experimental works which are already validated with each other earlier by the authors [8–10] are considered as forward problems, while the modified RNNs analogy is considered as inverse method in the present formulation. The present method has been also compared with author's previous works [12,13].
In the present investigation, a simply supported beam (Fig. 1) with multiple cracks subjected to a traversing load has been considered for the analysis. The Euler Bernoulli's beam analogy is applied in the present problem without inclusion of the damping effect of the beam. The problem formulation has already described by the authors [8–10] in their previous work.
δ = Dirac delta function.
= Applied force.
γ = vt = Position of the traversing mass at any moment of time ‘t’.
L1, L2, L3 = L1,2,3 = Locations of the first, second and third cracks from the left supported end of the cracked simply supported structure respectively.
d1, d2, d3 = d1,2,3 = Deepness of the first, second and third cracks respectively.
M = Mass of the traversing mass, ν = Traversing speed.
EI = Flexural rigidity of the structure.
m = Beam mass per unit length. g = gravitational constant, L = Length of the beam.
The solution of the above equation (1) has been written in series from i.e.(2)where Yn(x) = Shape function of the structure without the transit mass.
Tn(t) = Amplitude function, n = Number modes of vibrating structure.
For the evaluation of Y(x), the equation (2) has been expressed as-(3)
A numerical example has been formulated to determine the responses of the simply supported beam subjected to the traversing load. The forward problem (Numerical, FEA and experimental works) has been considered in the same way of Jena and Parhi [8–10] as of their earlier approach and responses of the structure are also verified using FEA and experimentation methods by the authors in their previous works. The transient dynamic analysis method (Newmark-β integration method) has already been carried out to determine the responses of the simply supported structure under transit mass in ANSYS WORKBENCH 15 domain. The laboratory test method was also conducted for the proposed work. The test procedures (Fig. 2) remain same as the procedures of the author's  previous work. The frictional force between the moving mass and the beam is neglected here. The responses (amplitudes) of the beam are recorded at various desired locations. The responses of the structure are also observed from the test method. The FEA model of the present problem is shown in Figure 3 (Tab. 1).
Schematic presentation of cracked simply supported beam under traversing mass.
Laboratory test model for simply supported structure.
FEA model for the moving load problem for η1,2,3 = 0.331, 0.545, 0.703. α1,2,3 = 0.13, 0.28, 0.41. M = 1.13 kg
Frequency ratio for the cracked beam at different modes.
To corroborate all the adopted method, a numerical example has been formulated.
Beam type-Rectangular (Mild steel specimen)
Beam size-145 cm × 4.5 cm × 0.5 cm.
Moving speed- 5.23 and 6.41 m/s.
Moving mass- 1.13 and 2.6 kg.
Relative crack locations-.
Relative crack depth-
The natures of the cracks on the structure are open and transverse and the detailed analysis of cracks are already analyzed by Jena and Parhi . The cracks are open and transverse in nature, the load is assumed to be acting on the transverse direction only. From the analysis of Figure 4, the sudden rise in dynamic amplification on the beam response is occurring due to only the presence of cracks. From the observed dynamic amplification, the positions and severities of cracks are determined by FEA and experimental analyses. The severities of cracks are quantified using the concepts of natural frequencies and mode shapes analysis mechanism using FEA approach by Chasalevris and Papadopoulos  and Bagherahmadi and Seyedpoor . Based upon the analysis of the response, the crack positions and depth are observed.
Graph for crack detection at M = 1.13 kg, v = 5.12 m/s, η1,2,3 = 0.358, 0.503, 0.634, α1,2,3 = 0.358, 0.503, 0.634.
RNNs are those kinds of networks where the output data from the preceding step can be fed as input data to the present step. In case of conventional NNs, the data from input and output are independent with each other. But in RNNs, the data are dependent with each other and can be used whenever it's required. RNNs have additional memory (context layer) which remembers the information that has been required later. In the present analogy, a modified RNNs based approach has been developed to identify the positions of cracks and determine the severities of crack in a supervised manner. The modified RNNs structure is the combination of the modified ERNNs and JRNNs structures. The architecture of the proposed RNNs structure has been represented in Figure 5. Based upon the architectural structure, the equations are developed for the modified RNNs. The different terminologies used in the network structure are described here. Primarily, the proposed RNNs model is trained with the training data parameters like relative deflections of the beam, number of cracks, speed of the transit mass, length of the structure, mass of the transit mass as input data and crack locations & severities as output data.
RD = Relative deflection at the crack position = Deflection of the beam at the position of crack to that of the beam without crack with the same moving speed and mass.
RD-1, RD-2 and RD-3 are the relative deflections of the structure at the first, second and third crack positions respectively.
n and n ′ are the number of cracks in the input and output layers of the network model respectively.
W1 = RD-1. W2 = RD-2. W3 = RD-3.
W4 = Length of the beam (L).
W5 = Speed of the moving mass (v).
W6 = Magnitude of the moving mass (M).
W7 = Number of cracks.
φ1 = First relative crack location (η1).
φ2 = First relative crack depth (α1).
φ3 = Second relative crack location (η2).
φ4 = Second relative crack depth (α2).
φ5 = Third relative crack location (η3).
φ6 = Third relative crack depth (α3).
φ7 = No of cracks present in the structure (n ′).
where, i = 1, 2 ... N, ‘N’ is the total input nodes.
r = 1, 2...r1, ‘r1’ is the total number of neurons in the context layers, the context layers 1 & 2 having same number of neurons.
r = 1, 2 ... r2, ‘r2’ is the total number of neurons in the context layers, the context layers 3 & 4 having same number of neurons.
k = 1, 2 ... k, ‘k’ is the number of nodes in the output layer.
‘j1= j2 = j3=1, 2 ... j’, ‘j’ is the number of neurons in each of the hidden layers, that is, first, second and third hidden layers which is constant throughout the layer.
‘β’ is the self-recurrent value for each node in the output and context layers that varies from 0 to 1.
Z−1 is the unit delay.
w is weight of the connection.
U1−7 and U ′ 1−7 are the context unit values in the context layers-3 and 4 respectively.
V1−13 and V ′ 1−13 are the context unit values in the context layers-1 and 2 respectively.
φc 1 t and φc 1 t −1 are the output values of the nodes at the context layer-1 at time index ‘t’ and ‘t−1’ respectively.
φc 2 t and φc 2 t −1 are the output values of the nodes at the context layer-2 at time index ‘t’ and ‘t−1’ respectively.
φc 3 t and φc 3 t −1 are the output values of the nodes at the context layer-3 at time index ‘t’ and ‘t−1’ respectively.
φc 4 t and φc 4 t −1 are the output values of the nodes at the context layer-4 at time index ‘t’ and ‘t−1’ respectively.
φj t are φj t −1 are the output values of the nodes at the hidden layer at time index ‘t’ and ‘t−1’ respectively.
φk t and φk t − 1 are the output values of the nodes at the output layer at time index ‘t’ and ‘t−1’ respectively.
‘t−1’ is the time index that is deferred by one-time step due to the feedback connections which is represented by Z−1 in the network architecture.
‘w’ is the weight of connection.
By analyzing the modified RNNs structure (Fig. 5), the mathematical equations for the structure are extended as follows:(5) (6) (7) (8)
The applied activation functions in the hidden and output layers are denoted by f (.) and g (.) respectively. The ‘tan-sigmoid’ activation function has been applied for the hidden layers, while that of output layer is ‘purelin’. The error functions (ε) of the output neurons are estimated approximately by employing the efficient weight factors relation, that is, w new = w old + ϕΔw, where ‘ϕ’, the learning rate coefficient that is varies from 0 to 1. The training procedures of the proposed network model are carried out using sum square error function mechanism. The L-M (Levenberg-Marquardt) back propagation algorithm by Yu and Wilamski  is implemented to the proposed network to estimate the crack locations and depth of the structures. The training of the network model has been conducted with the same procedures of authors [12,13] earlier works and Yu and Wilamski . Due to the fast and steady convergence properties, the L.M back propagation algorithm is applied for the proposed network model. This algorithm is the combination of Gauss-Newton and steepest descent methods. The proposed algorithm adopts the mechanism of the speed improvement of the Gauss-Newton and steadiness of steepest descent methods. The algorithm accomplishes a computational solution to the problem to minimize the non-linear function. The equation of the L-M back propagation mechanism  is given-(13)
Where ‘J’ is the Jacobian matrix that is calculated from the Gauss-Newton approach. ‘I’ is the identity matrix. ‘ξ’ is the combination coefficient. The value of ‘ξ’ approximate to zero, then equation (13) will follow the mechanism of Gauss-Newton and, if that of ‘ξ’ is quiet greater value, then that will follow the mechanism of steepest descent methods.
ξ = (1/ν), ‘ν’ is the step size or training constant.
e = Error vector = ϕdesired − ϕactual.
The proposed network training procedures are conducted using the rules of sum square error function. As per the rules of the sum square error function, if the present error value is less than the preceding one, it seems that the quadratic approximation on the sum square error is executing properly, then the value of ‘ξ’ is reduced to minimize the implication of gradient descent section. But, if the present error value is greater than the preceding one, then the value of ‘ξ’ is enhanced.
The present RNNs model is trained with 750 numbers of patterns with different parameters of the moving load dynamic problem. From the 750 numbers of patterns, 650 patterns are used for training procedures of the network while 100 patterns are for testing procedures of the networks.
While training the network model, the number of neurons and hidden layers are selected on the basis of cross-validation analogy. The model data (neurons, layers) has been determined from the available training data set or patterns (650 numbers). The remaining 100 patterns are reserved to test the model. The number of neurons in each hidden layer is found to be 13, while the number of hidden layer to be 3. The numbers of neurons in the context layers (3 & 4) are 7 while those in context layers (1 & 2) are 13 respectively. It's due to the reason that the data or information can be copied and transformed exactly from the respective preceding layers and forwarded to the proceeding layer. The training method has been conducted using the L-M algorithm. The patterns generated for the training process is exemplified in Table 2. The patterns are generated for both for the damaged and healthy structures. In Table 2, in the input parameters, the values of RD characterize to 1 indicates that there is no crack at that particular location which leads the values of are η and α zero in the output parameters. The present methodology can identify up to three numbers of cracks.
Modified architecture of RNNs model.
Training Patterns for modified RNNs model.
In the present investigation, a parameter identification (crack locations and severities) problem has been developed using modified recurrent neural networks in a supervised manner. The modified RNNs blend the combination of modified ERNNs and JRNNs. The L-M back propagation algorithm has been implemented to train and test the networks. The present formulation is based to explore an indirect approach for fault detection in moving load dynamic problem in the domain RNNs. To check the convergence of the developed method (Modified RNNs based approach), a numerical study has been carried out followed by FEA and laboratory corroboration. To substantiate all the adopted method, the same numerical example has been considered which has been analyzed in the numerical formulation part. The detailed analysis of the problem in the forward manner has been earlier discussed by the authors earlier. The severities of cracks are found out by analyzing FEA model. 750 numbers of patterns are generated for the proposed analysis, where 650 patterns are used for training process while rests of these are used for testing procedures.
The L-M back propagation methodology has been applied to train the modified RNNs structure. The present methodology is carried out in a supervised way as an inverse problem. 900 numbers of iterations are carried out. The results (relative crack locations and severities) found out from each of the analogies (FEA, experiment and RNNs) are corroborated with the numerical values and represented in Tables 3 and 4. The percentages of error or deviation values are calculated with the exemplified numerical values. The percentage of error values of FEA and experiments with numerical values are about 1.9% & 2.8% for relative crack locations while those for relative crack depth are about 2.1% & 3.0% respectively.
A graphical representation between sum square error values and number of iteration for both relative crack locations and severities are shown in Figure 6. The results obtained from the modified RNNs analogy vary with approximately 3.5% with the numerical values for relative crack locations while the results for crack severities vary about 3.9% approximately with the numerical values which are considered as good convergent. In the earlier works, the authors have also applied the analogy of the modified ERNNs  and JRNNs  based approaches individually with L-M back propagation algorithm.
The authors have also compared the present results with their earlier works [12,13] and found that the present approach yields better results. It has been observed that the modified recurrent neural networks approach yield better accuracy as comparison to the individual ERNNs and JRNNs analogies. So the proposed method can be used efficiently to identify, locate and quantify cracks on structure subjected to travelling mass.
Comparison of results for relative crack locations (η) with various methods.
Comparison of results for relative crack depth (α) with various methods.
Number of iterations vs. sum square error values (%) for relative crack depth and locations.
The present study has quietly elaborated a structural health monitoring problem in the domain of moving load dynamics problem using modified RNNs based approach. The developed RNNs based analogy has been carried out in a supervised manner. The numerical, FEA and experimental studies are considered as the forward problem while the RNNs based approach as inverse or indirect approach in the present investigation. The proposed RNNs structure has been trained by implementing the mechanism of L-M back propagation algorithm. The existence, relative crack locations and depth of the simply supported structure under travelling mass are found out by both the forward and inverse approaches. The obtained results (α and η) from RNNs methodology are compared with the numerical problem values and found to be about 3.9% and 3.5% for relative crack depth and positions respectively which seems to be convergent. The proposed method can be used effectively for crack detection in structure under moving load.
- S.A. Bagherahmadi, S.M. Seyedpoor, Structural damage detection using a damage probability index based on frequency response function and strain energy concept, Structural Engineering and Mechanics 67 , 327–336 (2018) [Google Scholar]
- A.C. Chasalevris, C.A. Papadopoulos, Identification of multiple cracks in beams under bending, Mechanical Systems and Signal Processing 20, 1631–1673 (2016) [Google Scholar]
- C.-M. Chang, T.-K. Lin, C.-W. Chang, Applications of neural network models for structural health monitoring based on derived modal properties, Measurement 129 , 457–470 (2018) [Google Scholar]
- M. Duimitriu, Fault detection of damper in railway vehicle suspension based on cross-correlation analysis of bogie accelerations, Mechanics & Industry 20 , 102 (2019) [EDP Sciences] [Google Scholar]
- U. Eroglu, E. Tufekci, Exact solution based finite element formulation of cracked beams for crack detection, International Journal of Solids and Structures 96 , 240–253 (2016) [Google Scholar]
- R. Hou, Y. Xia, Q. Xia, X. Zhou, Genetic algorithm based optimal sensor placement for L1‐regularized damage detection, Structural Control and Health Monitoring 26 , 1–14 (2018) [Google Scholar]
- S.J.S. Hakim, H.A. Razak, S.A. Ravanfa, Fault diagnosis on beam-like structures from modal parameters using artificial neural networks, Measurement 76 , 45–61 (2015) [CrossRef] [Google Scholar]
- S.P. Jena, D.R. Parhi, D. Mishra, Comparative study on cracked beam with different types of cracks carrying moving mass, Structural Engineering and Mechanics 56 , 797–811 (2015) [Google Scholar]
- S.P. Jena, D.R. Parhi, Parametric study on the response of cracked structure subjected to moving mass, Journal of Vibration Engineering & Technologies 5 , 11–19 (2017) [Google Scholar]
- D.R. Parhi, S.P. Jena, Dynamic and experimental analysis on response of multi-cracked structures carrying transit mass, Proceedings of the Institution of Mechanical Engineers, Part 0: Journal of Risk and Reliability 231 , 25–35 (2017) [Google Scholar]
- S.P. Jena, D.R. Parhi, Dynamic response and analysis of cracked beam subjected to transit mass, International Journal of Dynamics and Control 6 , 961–972 (2018) [Google Scholar]
- S.P. Jena, D.R. Parhi, Fault detection in cracked structure under moving load using RNNs based approach, Scientia Iranica: Transaction of Mechanical Engineering-B (2019) [Google Scholar]
- S.P. Jena, D.R. Parhi, Structural damage detection in moving load problem using JRNNs based approach, Journal of Theoretical and Applied Mechanics 57 , 665–676 (2019) [Google Scholar]
- L. Mao, Y. Lu, Critical speed and resonance criteria of railway bridge response to moving trains, Journal of Bridge Engineering 18 , 131–141 (2013) [Google Scholar]
- Y. Pala, M. Reis, Dynamic response of a cracked beam under a moving mass load, Journal of Engineering Mechanics 139 , 1229–1238 (2013) [Google Scholar]
- C.S. Nadith Pathirage, J. Li, L. Li, H. Hao, C. Liu, P. Ni, Structural damage identification based on autoencoder neural networks and deep learning, Engineering Structures 172 , 13–28 (2018) [Google Scholar]
- H. Yang, W. Yan, H. He, Parameters Identification of moving load using ANN and dynamic strain, Shock and Vibration 16 , 1–13 (2016) [Google Scholar]
- Q. Yang, J.K. Liu, B.X. Sun, C.F. Liang, Damage localization for beam structure by moving load, Advances in Mechanical Engineering 9 , 1–6 (2017) [CrossRef] [Google Scholar]
- J. Wu, J. Liaug, S. Adhikari, Dynamic response of concrete pavement structure with asphalt isolating layer under moving loads, Journal of Traffic and Transportation Engineering 1 , 439–447 (2014) [Google Scholar]
- W. Zhang, J. Lib, H. Haob, H. Mac, Damage detection in bridge structures under moving loads with phase trajectory change of multi-type vibration measurements, Mechanical Systems and Signal Processing 87 , 410–425 (2017) [Google Scholar]
- H. Yu, B.M. Wilamski, Levenberg-Marquardt Training, The Industrial Electronics Handbook, Intelligent Systems, 2nd Edition, 12, CRC Press, New York, 2011, pp. 1–16 [Google Scholar]
Cite this article as: S.P. Jena, D.R. Parhi, On the way to fault detection method in moving load dynamics problem by modified recurrent neural networks approach, Mechanics & Industry 22, 10 (2021)
FEA model for the moving load problem for η1,2,3 = 0.331, 0.545, 0.703. α1,2,3 = 0.13, 0.28, 0.41. M = 1.13 kg
|In the text|
Graph for crack detection at M = 1.13 kg, v = 5.12 m/s, η1,2,3 = 0.358, 0.503, 0.634, α1,2,3 = 0.358, 0.503, 0.634.
|In the text|
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.