Research Article  Open Access
Khoa Viet Nguyen, Mai Van Cao, "Exact Receptance Function of Cracked Beams and Its Application for Crack Detection", Shock and Vibration, vol. 2019, Article ID 9485351, 19 pages, 2019. https://doi.org/10.1155/2019/9485351
Exact Receptance Function of Cracked Beams and Its Application for Crack Detection
Abstract
The receptance function is very important which interrelates the harmonic excitation and the response of a structure in the frequency domain. This paper presents the exact receptance function of cracked beams. In this work, the “receptance curvature” is defined as the second derivative of the receptance. The influence of the crack on the receptance curvature is investigated. The results show that when there are cracks, the receptance curvature is influenced significantly at crack positions. This might be useful for the detection of cracks. In this paper, the derivation of exact receptance of the beam with general boundary conditions is presented, and the numerical simulations are provided.
1. Introduction
As is known the receptance method has been applied widely in mechanical system and structural dynamics. Yang [1] presented the exact receptances of nonproportionally damped dynamic systems. Based on a decomposition of the damping matrix, an iteration procedure is developed which does not require matrix inversion. Lin and Lim [2] developed a new and effective method to derive the structural design sensitivities which include the receptance sensitivity with respect to mass modification and stiffness modification from the limited vibration test data. Mottershead [3] investigated the measured zeros from frequency response functions and their application to model assessment and updating. Gurgoze [4] was concerned with receptance matrices of viscously damped systems subject to several constraint equations. The frequency response matrix of the constrained system was established in terms of the frequency response matrix of the unconstrained system and the coefficient vectors of the constraint equations. Karakas and Gurgoze [5] presented a formulation of the receptance matrix of nonproportionally damped dynamic systems. The receptance matrix was obtained directly without using the iterations. Mottershead et al. [6] and Tehrani et al. [7] addressed the problem of polezero assignment using the receptance method in active vibration control. Bernal [8] presented a receptancebased formulation for modal scaling using mass perturbations. Bonello and Minh Hai [9] used the receptance harmonic balance technique for the computation of the vibration of a whole aeroengine model with nonlinear bearings. Albertelli et al. [10] proposed a method using the receptance coupling substructure method to improve chatter free cutting conditions prediction. Muscolino et al. [11] presented the explicit frequency response functions of discretized structures with uncertain parameters. Recently, Muscolino and Santoro [12] presented the explicit frequency response function of beams with cracks of uncertain depths in order to evaluate the main statistics as well as the upper and lower bounds of the response.
Early crack detection is extremely important in mechanical systems and engineering structures, and this issue has attracted many researchers in the last three decades. Since the cracks may influence significantly on the dynamic characteristics of structures such as natural frequencies and mode shapes, these dynamic characteristics have been investigated and applied widely for crack detection of structures. Chondros and Dimarogonas [13] and Lee and Chung [14] studied the change in natural frequencies of beams caused by a crack. Hu and Liang [15] presented an integrated method using the knowledge of changes in natural frequencies for crack detection. Zheng and Kessissoglou [16] investigated the relationship between natural frequency of a cracked beam to the depth and location of the crack. The results in these papers showed that the natural frequency of the cracked beam decreases as the crack depth increases. Gudmundson [17] studied the influence of cracks on the natural frequencies of slender structures using a flexibility matrix approach. Ruotolo and Surace [18], VakilBaghmisheh et al. [19] applied genetic algorithms for crack detection of structures using natural frequencies. In other work, Thalapil and Maiti [20] investigated the change in natural frequencies caused by longitudinal cracks for crack detection. Gillich et al. [21] proposed new mathematical expressions applying the multimodal analysis for damage detection. The proposed method makes the frequencybased damage detection method robust and not influenced by temperature variations. Along with the change in natural frequencies, the change in mode shapes caused by cracks has been investigated in many previous works. Some authors presented methods to calculate and apply the mode shape of cracked structures for crack detection purpose. Very early works by Pandey et al. [22] and Abdel Wahab [23] used the change in curvature mode shapes for determining the location of cracks in beam structures. The author of this paper [24] applied 3D finite elements to investigate the change in mode shapes at the crack positions. However, the mode shapebased methods for damage detection have the disadvantage that the mode shapes cannot be measured directly but they have to be calculated from the frequency response functions (FRF). To overcome this, Sampaio et al. [25] presented a method using the frequency response function curvature for crack detection which is the extension of the method in [22]. However, the curvatures of mode shape and FRF in [22, 25] were calculated by approximation methods, and no mathematical formulas for the curvature mode shape and FRF curvature were presented. Recently, Caddemi and Calio [26] presented the exact closedform solution for the mode shapes of the Euler–Bernoulli beam with multiple open cracks. This result can be applied to establish the exact formula of the receptance function to improve the accuracy of the methods using the receptance matrix.
In previous works, the receptance of cracked beams was derived approximately while the exact form of the receptance of cracked beams has not been presented. Mode shape and its curvature have been applied to investigate the influence of the crack on the dynamic response of structures. However, these works were based on approximate methods. The aim of this work is to establish the exact formulas of the receptance function and receptance curvature of a cracked beam using the exact mode shape. The exact receptance function provides convenient tool to calculate and study the distribution of response amplitude of the beam corresponding to harmonic force acting along the beam at any forcing frequency. The influence of cracks on the receptance curvature is investigated for the crack detection purpose. Simulation results show that the receptance curvature has significant changes at crack positions. It is concluded that only the low forcing frequency is recommended for crack detection while the high forcing frequency is not recommended. The experiment has been carried out to justify the formulas obtained in this work.
2. Derivation of the Receptance Function of a Cracked Beam
2.1. Intact Beam
In this work, the undamped simply supported Euler–Bernoulli beam is considered. The forced vibration equation of the undamped beam can be written as follows:where is the nondimensional coordinate.
The solution of equation (1) can be found in the formwhere is the n^{th} mode shape of the beam and is the timedependent amplitude which is referred to as generalized coordinate. Substituting equation (2) into equation (1) yields
Multiplying both sides of equation (3) with and integrating give
The orthogonality conditions of the beam can be expressed as
Applying this orthogonality conditions yields
If the force is harmonic exciting at the point , then and from (6) we have
Therefore, the receptance at due to the force at can be obtained as
The receptance curvature matrix is defined as the second derivative of the receptance function with respect to ξ variable as follows:
As can be seen from equation (19), the receptance curvature is dependent on the curvature of the mode shape. It is well known that when there is a crack, the mode shape is changed at the crack position but it is small when the crack size is small. However, the receptance curvature may have significant change at the crack position in comparison with that of the receptance.
2.2. Cracked Beam
The orientation of cracks changes with location and loading conditions and is not always vertical or normal to the beam central line. The orientation of crack influences the change in stiffness at the crack location and the position of the crack as presented in [27]. In this study, for simplicity reason, we consider only cracks which are open and perpendicular to the beam surface. In order to derive the exact receptance function of the cracked beam, the exact closed form of the mode shape of the simply supported Euler–Bernoulli beam with n cracks is adopted from [26] as follows:whereα_{k} is the dimensionless frequency parameter ; is the position of the ith crack, .
For simplicity reason, some notations are presented as follows:where is Heaviside function.
The parameters in (10) are presented as follows:where λ_{i} is the severity parameter of the crack and is defined as the ratio between the crack depth d_{i} and the beam height h.
From equations (10) and (11), the square of the mode shape can be calculated aswhere
To integrate the terms including Heaviside function in the square of the mode shape, following properties of Heaviside function should be noted:where is the antiderivative function of .
From equation (16), we have
Equation (17) cannot be used directly to derive the explicit formula of . However, we can express equation (17) in the formwhere δ_{ij} is the Kronecker delta.
Applying equations (14)–(18), the integral of the square of the mode shape can be derived. It is interesting that the term will vanish, and the following formula is obtained:
Substituting equations (10) and (19) into equation (8), the exact receptance of the cracked beam will be obtained.
In order to derive the exact formulas of curvature receptance, the second derivative of the mode shape needs to be calculated. It is noted that
Applying the following properties of Heaviside function and Dirac delta function [27]yields
From equations (9) and (22), the second derivative of the mode shape can be derived as follows:
Substituting equations (19) and (23) into (9), the exact receptance curvature will be determined.
The exact formulas of receptance and its curvature of cracked beams with other general boundary conditions can be obtained by the same procedure and can be seen in Appendix.
3. Numerical Simulation and Discussions
Numerical simulations of a simply supported beam with two cracks are presented in this section. Parameters of the beam are mass density ρ = 7800 kg/m^{3}; modulus of elasticity E = 2.0 × 10^{11} N/m^{2}; L = 1 m; and b = 0.02 m; h = 0.01 m. The first three natural frequencies of the intact beam are calculated and listed in Table 1.

Figure 1 presents the 3D graphs of normalized receptance matrices at the first two natural frequencies of the intact beam. As can be seen from Figure 1(a), when the forcing frequency is equal to the first natural frequency, the receptance is zero at the ends of the beam and then increases and reaches the maximum when the response point moves from the ends to the middle of beam. When the forcing frequency is equal to the second natural frequency, the receptance matrix is minimum at the middle of the beam and becomes maximum at positions L/4 and 3L/4 as depicted in Figure 1(b). When the forcing frequency is equal to the third natural frequency, there are three peaks at L/6, 3L/6, and 5L/6 as illustrated in Figure 1(c). The distribution of response amplitude of the beam can be investigated clearer when the force position is fixed. Figure 2 illustrates the 2D graphs of receptance when the force position is at 0.58L. As can be seen from this figure, the receptance curves are smooth except the minimum position. The positions of maxima and minima of these receptance graphs coincide with the positions of nodes of the corresponding mode shapes. These results imply that when the beam is excited at natural frequencies, the distribution of response amplitude along the beam can be predicted easily by using the receptance matrices as well as using the mode shapes.
(a)
(b)
(c)
(a)
(b)
(c)
However, when the forcing frequency is in between any two natural frequencies, the distribution of response amplitude is more complicated. For example, as presented in Figure 3(a), when the forcing frequency is 400 Rad/s, between the first and the second frequencies, the receptance matrix has two main peaks in the middle area of the beam and two small peaks close to the ends. The shape of this receptance matrix suggests that it is a hybrid shape between the shapes of receptance of mode 1 and mode 2 where the two main peaks correspond to mode 1 and two other peaks correspond to mode 2. This means that in this case, mode 1 and mode 2 are both excited, but mode 1 is dominant. When the forcing frequency is 950 Rad/s, between the second and the third natural frequencies, the receptance matrix has six main peaks and two other small peaks as can be observed from Figure 3(b). In this case, modes 2 and 3 are both excited, but mode 3 is dominant. These results show that the distribution of response amplitude is complicated at arbitrary frequency and cannot be predicted using mode shapes, but it can be predicted easily by using the 3D graph of the receptance matrix.
(a)
(b)
3.1. Influence of the Crack on the Receptance
When there are cracks, the mode shapes will have changes at the crack positions leading to the changes in the receptance at the crack positions. These changes are dependent on the crack depths and the lengthtoheight ratio of the beam. Simulation results show that the changes in the receptance are small when the crack depth is small, and they only become significant when the crack depth is large. In our simulation, when the lengthtoheight ratio is fixed, the changes in the receptance are very difficult to be detected visually when the crack depth is smaller than 40% of the beam height. Figure 4 depicts the distortions of the receptance matrix at the first natural frequency when the crack depth is up to 50% of the beam height. As can be seen from this figure, the receptance matrix is distorted at the crack positions. The distorted positions coincide with the crack positions are at 0.4L and 0.76L as can be observed from Figure 4(b).
(a)
(b)
3.2. Influence of the Crack on the Receptance Curvature
In this section, in order to investigate the influence of the crack on the receptance curvature matrices, two cracks with the same depths are made at arbitrary positions of 0.4L and 0.76L from the left end of the beam. Three crack depth levels ranging from 10% to 30% have been applied. The receptance and receptance curvature matrices are calculated at 100 points spaced equally on the beam while the force moves along these points.
In order to amplify influence of the crack, power 4 of the receptance curvature is presented.
Figures 5–7 present the graphs of normalized receptance curvature matrices of the cracked beam with different levels of the crack depth when the forcing frequencies are at the first three natural frequencies, respectively. As can be seen from these figures, there are sharp peaks in the receptance curvature matrices at crack positions. In order to determine exactly the sharp peak positions, the receptance curvatures along the beam with a fixed force at the position of 0.58L are presented as shown in Figures 8–10. As can be observed from these figures, the positions of the sharp peaks are clearly detected at 0.4L and 0.76L which are the same with the crack positions. It is noted from Figures 5–10 that when the crack position is closer to the maxima of the receptance curvature matrices, the heights of sharp peaks are greater. When the crack position is far from the maxima of the receptance curvature, the heights of sharp peaks in receptance curvature are smaller. As can be observed from Figures 8–10, when the crack depth increases from 10% to 30%, the height of sharp peaks in receptance curvature increases and becomes clearer.
(a)
(b)
(a)
(b)
(a)
(b)
(a)
(b)
(a)
(b)
(a)
(b)
3.3. Influence of Forcing Frequency
When forcing frequency is equal to 146 Rad/s which is close to the first natural frequency, there are clear peaks at the crack position of 0.4L and 0.76L as presented in Figure 11(b). However, in higherorder modes, if the crack position coincides with a minimum of the receptance, the vibration at this position may be zero, so no vibration at this point can be measured; thus, the influence of the crack cannot be measured as well. If the crack position is close to the minimum of the receptance curvature, the influence of the crack will be small. Figure 11(a) presents the overview of the receptance curvature matrix obtained when the forcing frequency is 950 Rad/s. As can be observed from this figure, the peaks caused by the crack are disappeared at some minima of the receptance curvature. The receptance extracted when the force is fixed at 0.8L with the frequency of force 950 Rad/s is presented by the dotted line in Figure 11(b). The crack at 0.4L coincides with the minimum of the receptance curvature, so no peak at this position is presented. Meanwhile, the crack at 0.76L causes a small peak at this position in the receptance curvature since this crack is close to another minimum of receptance curvature. Thus, the receptance curvature corresponding to the low forcing frequency is recommended for crack detection while the receptance curvatures corresponding to the high frequency are not recommended. This conclusion has not been reported in previous works yet.
(a)
(b)
3.4. Influence of the Boundary Conditions
The proposed method can be applied for different boundary conditions since this method is just used to detect the disruption in the receptance curvature. However, the influence of the cracks on the receptance curvature will be different with different boundary conditions. Figure 12 presents the receptance curvatures obtained at the first natural frequency of the beam with clampedfree and clampedclamped conditions. As can be seen from this figure, two peaks at crack positions are presented in both cases. For the cantilever beam, the peak at 0.76L which is close to the free end of the beam is smaller than the peak at 0.4L which is close to the fixed end. For the clampedclamped beam, the peak at 0.76L which is close to the minimum of the receptance curvature is smaller than the peak far from the minimum as mentioned above.
(a)
(b)
4. Experiment Validations
The experimental setup is depicted in Figure 13. The simply supported steel beam with the same parameters presented in Section 3 has been tested. Two sawcut cracks perpendicular to the beam surface with the same depths ranging from 10% to 30% were made at 0.4L and 0.76L. The widths of the cracks are 1.2 mm, and these cracks work as open cracks. The forcing frequencies are 147.65 Rad/s and 146 Rad/s which are the first natural frequencies measured from the intact beam and the cracked beam, respectively. The position of the exciter is selected at 0.58L which is close to the middle of the beam to excite the first vibration mode. This position is also the midpoint of two cracks to limit the influence of exciter on these cracks. The beam is excited by the vibration exciter Bruel & Kjaer 4808, and the response is measured by the polytec laser vibrometer PVD100. The receptance is measured at 99 points spaced equally along the beam when the forcing frequency is set to the first natural frequency of the cracked beam. The receptance curvature is then calculated as the second derivative of the receptance.
It is noted that the excitation amplitude is important to obtain highquality data and should be chosen suitably. If the excitation amplitude is too small, the vibration signal will be polluted by the measurement noise. If the amplitude is too large, the vibration signal will be nonlinear, so the proposed method cannot be applied since the theory in this paper is only applied for linear problems.
In order to reduce the presence of measurement noise and amplify the peaks at crack positions, the fourth power of the receptance curvature is presented instead of the receptance curvature. The reason is that by presenting the fourth power of the receptance curvature, the low values including noise in the receptance curvature will approach zero earlier than the large values at crack positions. Figure 14 depicts the measured receptance and the fourth power of receptance curvature graphs measured at the first natural frequency with the crack depths ranging from 10% to 30%. As can be seen from this figure, the simulations and the experimental results are in very good agreement. The influence of the cracks on the receptance is small and cannot be detected visually. However, when the crack depth is as small as 10%, two peaks in the receptance curvature are presented at crack positions although it is polluted due to the measurement noise. When the crack depth is greater than 10%, the significant peaks are presented clearly at the positions of the cracks. These experimental results prove the correctness of the mathematical formulas and the numerical simulations presented in this work.
(a)
(b)
(c)
(d)
(e)
(f)
5. Conclusions
In this paper, exact formula of the simply supported cracked beam is presented. The graphs of receptance matrices are used to present the overall picture of the vibration amplitude distribution of the beam. When the forcing frequency is equal to the natural frequency, the vibration amplitude distribution is similar to the vibration mode shape. However, the vibration amplitude distribution is complicated when the forcing frequency is arbitrary.
Receptance curvaturebased techniques have been applied recently for crack detection. However, no exact mathematical expression has been established so far. In this study, the exact formula of the receptance curvature of the cracked beam is established and applied to investigate the crack influence on the dynamic characteristics of the beam. When there are cracks, the receptance curvature is changed significantly at the crack positions, and it can be used for crack detection. Only receptance curvatures corresponding to low forcing frequencies can be applied efficiently for crack detection while the receptance curvatures corresponding to the high frequencies are not recommended. This result has not been reported yet since no detailed investigation of receptance curvature graphs at different forcing frequencies has been carried out.
The experimental results and numerical simulations of receptance and its curvature are in very good agreement. This justifies that the formulas of the receptance and receptance curvature derived in this paper are reliable. Again, the experimental results confirm that the receptance curvature can be applied for crack detection. In this experiment, the crack can be detected when the depth is as small as 10% of the beam height with measurement noise.
Appendix
Cantilever Beam
ClampedClamped Beam
FreeFree Beam
Data Availability
All data included in this study are available upon request by contacting the corresponding author.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
This research was funded by the Vietnam National Foundation for Science and Technology Development (NAFOSTED) under grant no. 107.022017.300.
References
 B. Yang, “Exact receptances of nonproportionally damped dynamic systems,” Journal of Vibration and Acoustics, vol. 115, no. 1, pp. 47–52, 1993. View at: Publisher Site  Google Scholar
 R. M. Lin and M. K. Lim, “Derivation of structural design sensitivities from vibration test data,” Journal of Sound and Vibration, vol. 201, no. 5, pp. 613–631, 1997. View at: Publisher Site  Google Scholar
 J. E. Mottershead, “On the zeros of structural frequency response functions and their sensitivities,” Mechanical Systems and Signal Processing, vol. 12, no. 5, pp. 591–597, 1998. View at: Publisher Site  Google Scholar
 M. Gürgöze, “Receptance matrices of viscously damped systems subject to several constraint equations,” Journal of Sound and Vibration, vol. 230, no. 5, pp. 1185–1190, 2000. View at: Publisher Site  Google Scholar
 A. Karakas and M. Gurgoze, “A novel formulation of the receptance matrix of nonproportionally damped dynamic systems,” Journal of Sound and Vibration, vol. 264, pp. 733–740, 2003. View at: Publisher Site  Google Scholar
 J. E. Mottershead, M. G. Tehrani, S. James, and Y. M. Ram, “Active vibration suppression by polezero placement using measured receptances,” Journal of Sound and Vibration, vol. 311, no. 3–5, pp. 1391–1408, 2008. View at: Publisher Site  Google Scholar
 M. G. Tehrani, L. Wilmshurst, and S. J. Elliott, “Receptance method for active vibration control of a nonlinear system,” Journal of Sound and Vibration, vol. 332, no. 19, pp. 4440–4449, 2013. View at: Publisher Site  Google Scholar
 D. Bernal, “A receptance based formulation for modal scaling using mass perturbations,” Mechanical Systems and Signal Processing, vol. 25, no. 2, pp. 621–629, 2011. View at: Publisher Site  Google Scholar
 P. Bonello and P. Minh Hai, “A receptance harmonic balance technique for the computation of the vibration of a whole aeroengine model with nonlinear bearings,” Journal of Sound and Vibration, vol. 324, no. 12, pp. 221–242, 2009. View at: Publisher Site  Google Scholar
 P. Albertelli, M. Goletti, and M. Monno, “A new receptance coupling substructure analysis methodology to improve chatter free cutting conditions prediction,” International Journal of Machine Tools and Manufacture, vol. 72, pp. 16–24, 2013. View at: Publisher Site  Google Scholar
 G. Muscolino, R. Santoro, and A. Sofi, “Explicit frequency response functions of discretized structures with uncertain parameters,” Computers & Structures, vol. 133, pp. 64–78, 2014. View at: Publisher Site  Google Scholar
 G. Muscolino and R. Santoro, “Explicit frequency response function of beams with crack of uncertain depth,” Procedia Engineering, vol. 199, pp. 1128–1133, 2017. View at: Publisher Site  Google Scholar
 T. G. Chondros and A. D. Dimarogonas, “Identification of cracks in welded joints of complex structures,” Journal of Sound and Vibration, vol. 69, no. 4, pp. 531–538, 1980. View at: Publisher Site  Google Scholar
 Y.S. Lee and M.J. Chung, “A study on crack detection using eigenfrequency test data,” Computers & Structures, vol. 77, no. 3, pp. 327–342, 2000. View at: Publisher Site  Google Scholar
 J. Hu and R. Y. Liang, “An integrated approach to detection of cracks using vibration characteristics,” Journal of the Franklin Institute, vol. 330, no. 5, pp. 841–853, 1993. View at: Publisher Site  Google Scholar
 D. Y. Zheng and N. J. Kessissoglou, “Free vibration analysis of a cracked beam by finite element method,” Journal of Sound and Vibration, vol. 273, no. 3, pp. 457–475, 2004. View at: Publisher Site  Google Scholar
 P. Gudmundson, “The dynamic behaviour of slender structures with crosssectional cracks,” Journal of the Mechanics and Physics of Solids, vol. 31, no. 4, pp. 329–345, 1983. View at: Publisher Site  Google Scholar
 R. Ruotolo and C. Surace, “Damage assessment of multiple cracked beams: numerical results and experimental validation,” Journal of Sound and Vibration, vol. 206, no. 4, pp. 567–588, 1997. View at: Publisher Site  Google Scholar
 M.T. VakilBaghmisheh, M. Peimani, M. H. Sadeghi, and M. M. Ettefagh, “Crack detection in beamlike structures using genetic algorithms,” Applied Soft Computing, vol. 8, no. 2, pp. 1150–1160, 2008. View at: Publisher Site  Google Scholar
 J. Thalapil and S. K. Maiti, “Detection of longitudinal cracks in long and short beams using changes in natural frequencies,” International Journal of Mechanical Sciences, vol. 83, pp. 38–47, 2014. View at: Publisher Site  Google Scholar
 G.R. Gillich, H. Furdui, M. Abdel Wahab, and Z.I. Korka, “A robust damage detection method based on multimodal analysis in variable temperature conditions,” Mechanical Systems and Signal Processing, vol. 115, pp. 361–379, 2019. View at: Publisher Site  Google Scholar
 A. K. Pandey, M. Biswas, and M. M. Samman, “Damage detection from changes in curvature mode shapes,” Journal of Sound and Vibration, vol. 145, no. 2, pp. 321–332, 1991. View at: Publisher Site  Google Scholar
 M. M. Abdel Wahab and G. De Roeck, “Damage detection in bridges using modal curvatures: application to a real damage scenario,” Journal of Sound and Vibration, vol. 226, no. 2, pp. 217–235, 1999. View at: Publisher Site  Google Scholar
 K. V. Nguyen, “Mode shapes analysis of a cracked beam and its application for crack detection,” Journal of Sound and Vibration, vol. 333, no. 3, pp. 848–872, 2014. View at: Publisher Site  Google Scholar
 R. P. C. Sampaio, N. M. M. Maia, and J. M. M. Silva, “Damage detection using the frequencyresponsefunction curvature method,” Journal of Sound and Vibration, vol. 226, no. 5, pp. 1029–1042, 1999. View at: Publisher Site  Google Scholar
 S. Caddemi and I. Caliò, “Exact closedform solution for the vibration modes of the EulerBernoulli beam with multiple open cracks,” Journal of Sound and Vibration, vol. 327, no. 3–5, pp. 473–489, 2009. View at: Publisher Site  Google Scholar
 R. F. Hoskins, Delta Function, Horwood Publishing Limited, Sawston, Cambridge, England, 2nd edition, 2009.
Copyright
Copyright © 2019 Khoa Viet Nguyen and Mai Van Cao. 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.