Access the full text.

Sign up today, get DeepDyve free for 14 days.

Acoustics
, Volume 1 (4) – Sep 27, 2019

/lp/multidisciplinary-digital-publishing-institute/quantification-of-bolt-tension-by-surface-acoustic-waves-an-1B3L0UQGHe

- Publisher
- Multidisciplinary Digital Publishing Institute
- Copyright
- © 1996-2021 MDPI (Basel, Switzerland) unless otherwise stated Disclaimer The statements, opinions and data contained in the journals are solely those of the individual authors and contributors and not of the publisher and the editor(s). MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations. Terms and Conditions Privacy Policy
- ISSN
- 2624-599X
- DOI
- 10.3390/acoustics1040046
- Publisher site
- See Article on Publisher Site

acoustics Article Quantiﬁcation of Bolt Tension by Surface Acoustic Waves: An Experimentally Veriﬁed Simulation Study 1 , 2 1 , Hani Alhazmi and Rasim Guldiken * Mechanical Engineering Department, University of South Florida, Tampa, FL 33620, USA; haalhazmi@mail.usf.edu Mechanical Engineering Department, Umm Al-Qura University, Makkah 21421, Saudi Arabia * Correspondence: guldiken@usf.edu; Tel.: +1-813-974-5628 Received: 30 August 2019; Accepted: 26 September 2019; Published: 27 September 2019 Abstract: Quantifying bolt tension and ensuring that bolts are appropriately tightened for large-scale civil infrastructures are crucial. This study investigated the feasibility of employing the surface acoustic wave (SAW) for quantifying the bolt tension via ﬁnite element modeling. The central hypothesis is that the real area of contact in a bolted joint increases as the tension or preload is increased, causing an acoustical signature change. The experimentally veriﬁed 3-D simulations were carried out in two steps: A preload was ﬁrst applied to the bolt body to simulate the realistic behavior of bolted joint; and the SAW propagation was then excited on the top surface of the plate to reﬂect from the bolted joint. The bolt tension value was varied between 4 and 24 kN (properly tightened bolt) in the steps of 4 kN to study the eect of the bolt tension. The results indicate an increased reﬂected wave amplitude and a gradual phase shift, up to 0.5 s, as the bolt tension increased. Furthermore, the result shows that the distance between the ﬁrst reﬂected wave and the source becomes shorter as the preload increases, as hypothesized. A 1.9 mm dierence in the distance between the maximum and minimum preload was observed. As part of this study, the simulation results were also compared with the experimental results, and a good agreement between the simulation and experiments was demonstrated. Keywords: bolted joint; surface acoustic wave; bolt tension; ﬁnite element modeling; the real area of contact 1. Introduction The health of bolted joints depends mainly on the status of the threaded fastener since the fastener is responsible for holding two structural components together. Therefore, any failure occurring in the fasteners may lead to a catastrophic failure in the entire engineering structure or machine. One common problem that bolted joints commonly experience is the loosening of a bolt, indicating a decrease in the tension force on the bolt, which is equivalent to the clamping force on the joints. A bolt loosening can take place due to various causes such as plastic deformation of the bolted joint components, the slip between the joint parts, and wear of joint parts [1]. The bolt tension should be monitored periodically to prevent failure in bolted joints due to bolt loosening. Various useful techniques were developed over the last few decades to quantify the bolt tension. These techniques can be classified into two main categories: the techniques necessitating bolt contact, and the contactless bolt measurement techniques. The torque wrench technique is the most common method to estimate the bolt tension through applying a wide range of torque values to clamp the bolt. However, the measurement accuracy of this standard methodology can be low due to disregarding the friction between the bolt thread and the nut [2,3]. Another example of a measurement technique that operates by contacting the bolt is the acoustoelastic effect based method [4]. This technique requires attaching an Acoustics 2019, 1, 794–807; doi:10.3390/acoustics1040046 www.mdpi.com/journal/acoustics Acoustics 2019, 1 795 ultrasonic transducer to the bolt head to generate a combination of longitudinal and transverse ultrasonic waves that travel to the bolt end and then reflects to the transducer. By processing the time of flight (TOF) of wave data, the effective bolt length can be estimated which then enables the measurement of the tension force on the bolt based on the acoustoelastic constant of the material. Among the contactless bolt tension quantiﬁcation techniques, one standard methodology is the piezoelectric active sensing method [5]. In this technique, two piezoelectric patches are attached to the top surfaces of both plates of bolted joints. The role of the ﬁrst piezoelectric transducer is to generate the ultrasonic wave that travels through the interface contact between the two plates, while the second piezoelectric transducer is responsible for receiving the propagated wave. The energy of the received wave signal depends on the clamping force at the interface, which is then analyzed to quantify the bolt tension and identify bolt loosening. In our previous study, the authors experimentally investigated an alternative method to quantify bolt tension without contacting the bolt using a synthetic phased-array surface acoustic wave (SAW) operating at 5 MHz [6]. The SAWs were generated by a linear transducer array that is located only on one surface of bolted joints. The central hypothesis was that the real area of contact in a bolted joint increases as the tension or preload is increased, which causes an acoustical signature that can be easily measured with the SAWs synthetic phased array [6]. The phased array method is utilized to acquire high-resolution ultrasonic images. The reﬂection of SAWs created by the real contact area between the bolt head and the plate surface is then used to quantify the tension on the bolt. It was observed that the distance between the transducer and the reﬂected wave shortens as the bolt tension increases. This study takes the investigation of this technique to the next level by exploring the use of ultrasonic waves to quantify the bolt tension via simulating the surface acoustic wave (SAW) propagation over the bolted joint surface. The 3D simulations were carried out by a commercially available Finite Element modeling software (ANSYS 18.1). The simulations were conducted in two main steps: First, applying a preload to the bolt body to simulate the real behavior of bolted joint under the real conditions; second, exciting the SAWs signal on the top surface of the bolted joint to reflect from the interface between the bolt head and the top surface. For a thorough investigation, six separate simulation studies were carried out. Each simulation has different preload values to understand the influence of preload on the incident wave amplitude. The preload value was varied between 4 and 24 kN in the steps of 4 kN. 2. Fundamentals of Ultrasonic Waves The ultrasonic waves have been widely utilized in nondestructive testing and structural health monitoring due to their abilities to detect defects as well as to inspect large engineering structures. The ultrasonic waves can exist in various modes and frequencies and can travel a long distance in both ﬂat and curved structures with low energy loss [7,8]. The ultrasonic elastic waves in solids can be classiﬁed into bulk elastic waves and guided elastic waves. Unlike bulk ultrasonic waves, guided elastic waves are dispersive, which indicates that the phase velocity of wave relies on their frequency. As a result, the group velocity of bulk waves is not equal to their phase velocity. The guided ultrasonic elastic waves can propagate in solids as elastic surface waves (Rayleigh wave) which represents a combination of shear and longitudinal waves, Lamb waves, or shear horizontal waves. Lamb waves are the most complicated type of all since they have two dierent modes, symmetric or antisymmetric. Lamb waves can be excited in thin plates (ﬁnite plates), where the plate thickness is much smaller than the wavelength. The surface waves propagate along the free surface of the solid where the plate thickness is large compared to the wavelength, and the surface particles move in an elliptical path. The amplitude of the Rayleigh wave decays in the direction of the plate thickness [7–9]. In a perfectly elastic isotropic solid, the longitudinal (C ) and shear velocity (C ) are related to L s the material properties which can be expressed in Equations (1) and (2). The approximation of the Rayleigh wave velocity (C ), which depends on the shear wave speed and the material properties, can be expressed in Equation (3) [7]: Acoustics 2019, 2 FOR PEER REVIEW 3 (1) 𝐶 = Acoustics 2019, 1 796 𝜌(1 − 𝑣) 𝐶 = (2) C = (1) (1 ) C = (2) μ 0.87 + 1.12 𝐶 = ( ) (3) ρ 1+ 0.87 + 1.12 C = ( ) (3) where E is the elastic modulus, ρ is the density, µ is the modulus o 1 + f rigidity and ν is the Poisson ratio. where E is the elastic modulus, is the density, is the modulus of rigidity and is the Poisson ratio. 3. Simulation Work 3. Simulation Work 3.1. Initial Simulation Studies 3.1. Initial Simulation Studies In order to understand and verify the behavior of SAW propagation on a flat plate, a 40 mm thick, In13 or0 m derm to long, understand and 10 and 0 m veri m wi fy the de p behavior late wi of thSA out W any ho propagation les or b ono altﬂat s was s plate,im a 40 ulat mm ed in a thick, 130 mm long, and 100 mm wide plate without any holes or bolts was simulated in a commercially commercially available finite element software, (ANSYS Inc., Release18.1, Canonsburg, PA, USA). The plate available is taken ﬁnite element as steesoftwar l with Young e, (ANSYS ′s modu Inc., luRelease18.1, s of 200 Gpa,Canonsbur Poisson′s rat g, P io o A, f USA). 0.33, sh The ear plate modis ules o taken f as steel with Young’s modulus of 200 Gpa, Poisson’s ratio of 0.33, shear modules of 76.9 Gpa, and 76.9 Gpa, and density of 7850 kg/m [10]. For more consistent meshing control, the top surface was sl density iced to of 3m 7850 m inkg depth, a /m [10 n]. d then the two For more consistent slices were meshing consocontr lidate ol, d w the ith top a bsurface onded cwas ontasliced ct to be to a3 s o mm ne in depth, and then the two slices were consolidated with a bonded contact to be as one part [11]. part [11]. The incorporation of slices not only allows the localization of a dense mesh for the top sur The face incorporation and less dense of slices mesh not alon only g wallows ith the rem the localization aining thickness of a dense of the pla mesh te, for but i thet top also al surface lows and for less dense mesh along with the remaining thickness of the plate, but it also allows for reducing the reducing the computational time to solve the finite element (FE) model since the number of nodes is sign computational ificantly decrea timesed. to solve Initia the lly, ﬁnite a sing element le nod(FE) al po model int on t since he upper le the number ft sur offnodes ace was is se signiﬁcantly lected to represent a single decreased. Initially tr,ansd a single ucer nodal to exci point te the wa on the ve. However, the upper left surface auth was ors observed tha selected to repr t the genera esent a single ted wa transducer ve is circto ula excite r whithe ch lwave. eads to u However ndesira , the ble authors and non- observed realistic r that eflthe ectigenerated ons from the wave pla is te s ciricular des. Thes which e ref leads lectito ons i undesirable nteracted wi and th non-r the in ealistic cident rwa eﬂections ve, maki fr ng om the resu the plate lts sides. more comp These lica reﬂections ted to anal interacted yze. For thi with s reason, twel the incident ve wave, evenlymaking distributhe ted p results oint sources more complicated near the uppe to r lanalyze. eft end of the pla For this ter we eason, re setwelve lected to excite evenly distributed SAW where point easour ch point so ces nearur the ce repre upper s left ents end a si ofngle transducer. Th the plate were selected e ge tone excite rated wa SAWve wher is a e st each raight point line sour when it ce repr leesents aves th ae source single transducer and becom . The es a curv generated ed liwave ne as is it t arstraight avels far line from when the it source leaves , similar to our the source and experimental observation becomes a curved line as s. However it travels, far large plates nec from the sour essit ce, a similar te addit to ional point so our experimental urces to excite a particular w observations. However a,ve lar mode an ge plates d re necessitate duce the no additional ise levelpoint [12]. In a sour d ces ditito on to excite twel a particular ve transmiwave tting point sources, even mode and reduce the ly noise spaced three level [12 rece ]. Inivers wer additioneto introduced twelve transmitting to the model. The point sour se ces, areevenly nodal spaced points tthr hatee ar re eceivers located wer at t ehintr e to oduced p surface to o the f tmodel. he plate. These Both the l are nodal eft a points nd right end of the pla that are located at tthe e ha top d a fi surface xed boundar of the plate. y con Both dition. the Fi left gure and 1 il right lustrat end esof ththe e geomet plate had rical det a ﬁxed ails and boundar boundary condition. y condit Figur ions e of 1 tillustrates he initial verificat the geometrical ion studie details s. and boundary conditions of the initial veriﬁcation studies. Figure 1. Schematic representation of the top surface of the plate, including the locations of the three receivers and twelve transmitting point sources. All dimensions are in mm. Acoustics 2019, 1 797 The main goal of the initial simulation studies was to obtain the Rayleigh surface wave velocities at each receiver and validate them against the theoretical value reported elsewhere. This was achieved by obtaining the zero-crossing time of the wave at each receiver. Another goal of the initial simulation studies was to obtain the length of time that the surface wave needed to travel from the wave source to the right end of the plate and verify it with the theoretical time. By knowing the entire time needed for the input signal, unwanted reﬂections from the far end of the plate can be avoided, and the computational time needed can be reduced [13]. Most FEA commercial software lacks a perfectly matched layer (PML), which is essential for absorbing or attenuating propagating waves and hence preventing them from being reﬂected in the study area. 3.2. Preload Analysis The bolt pretension simulation study aims to replicate the real behavior of bolted joints after applying a particular preload value on the bolt body. Typically, the preload on a bolt occurs when applying the torque to either the bolt head or the nut. There is a directly proportional correlation between the torque and preload on a bolt: The higher the applied torque to bolt, the higher the preload is on the bolt and vice versa. However, the preload magnitude not only depends on the applied torque, but also on other important factors such as the bolt diameter, the friction coecient between the bolt body and plate surfaces, and bolt thread geometry [14]. The preload and torque relationship is described by the following equation: T = KFd (4) where T is the applied torque, F is the preload, d is the bolt diameter, and K is the nut factor. The estimated range of K value is 0.10~0.2 depending on the type of lubricant, size of bolt, and bolt material [14]. In this study, K = 0.18 was utilized [6]. To avoid the plastic deformation on the bolt, the applied preload on a bolt should not exceed the maximum allowable limit. Thus, the maximum allowable preload can be estimated using Equation [14]: F = 0.75A 0.85S (5) i t y where F is the maximum allowable preload (proof load), A is tensile stress area which is a function i t of thread pitch, S is yield strength of the bolt, and d is the nominal bolt diameter. In this study, a standard medium carbon steel bolt and nut with 10 mm nominal diameter were imported from the commercially available CAD software library, Inventor Autodesk 2018. Based on the bolt material 2 2 and nominal diameter, the value of S and A are 640 N/mm [15] and 57.1 mm respectively. After y t substituting all values into Equation (5), the maximum allowable preload that bolt can withstand is obtained as 23.3 kN, and it is rounded up to 24 kN. By substituting this value into Equation (4), the maximum allowable torque on bolt is 43 N.m. For comparison purposes, six values of preload were utilized in this study. The preload values varied between 4 and 24 kN in the steps of 4 kN. To properly perform a bolt pretention analysis in ANSYS, the type of contact between bolted joint parts should be carefully selected. These contacts represent how each part moves relative to other parts. The coecient of friction, C , is assumed to be 0.2 for all frictional contact type [16]. Table 1 shows the details of the contacts used in this study. ANSYS 18.1 oers the bolt pretension feature as a boundary condition that applies to both static and transient structural analyses. This boundary condition is essential for obtaining the real behavior of bolt in joints since it oers ﬁve dierent types of loading on the bolt body. These types of loading are load, adjustment, lock, open, and increment. In the structural analysis, this study utilized two steps in the boundary condition. The ﬁrst step was to simulate the assembly of joints by applying a load on the bolt body as a preload. The second step was to simulate the real behavior of a joint under real conditions by applying a lock on the bolt [17]. From this analysis, the adjustment value (the change in the bolt length) was obtained. This is used in the ﬁrst step of the transient analysis. Acoustics 2019, 1 798 Table 1. Contacts set-up between the bolt and plate. Target Surface Source Surface Contact Type Clearance Hole The Bolt Body Frictional with C = 0.2 Top Surface of the Plate Bottom Surface of the Bolt Head Frictional with C = 0.2 The Bolt Body The Nut Inner Surface Bonded Bottom Surface of the Plate Top Surface of the Nut Frictional with C = 0.2 3.3. Transient Simulations The 3-D transient simulations were performed in two main steps: First, the adjustment value on the bolt body was applied, followed by locking the bolt to simulate the bolted joints behavior; then, applying a displacement load that acts perpendicular to the top surface of the plate at the twelve point sources. Before running the simulation, the basic parameters of a continuous ultrasound wave were deﬁned for simulating the real behavior of waves in a solid material. The ﬁrst parameter is the Rayleigh wave propagation speed, C , that mainly depends on the mechanical properties of the material, and it can be estimated using Equation (3). The second parameter is the wavelength () which is equal to the ratio of the ultrasonic wave speed (C ) to the frequency (ƒ) in a perfectly elastic material = C / f . R R By utilizing the material properties for steel and Equation (3), the Rayleigh wave speed (C ) was found to be 2917.58 m/s. The corresponding wavelength for this wave speed and the central frequency of interest (100 kHz) was obtained as 29.1 mm, which was rounded to 29 mm. In order to achieve an accurate solution of ultrasound wave propagation using the ﬁnite element analysis (FEA), an adequate temporal and spatial discretization is required. A precise time step, Dt, is needed to model the temporal behavior of the wave propagation and should be smaller than the critical time step of the studied phenomena. A proper element size should be selected for solving the wave propagation spatially. In general, the smaller the time step and element size, the higher is the solution accuracy. Although the solution accuracy improves with a smaller value of the time step and element size, unnecessarily small values necessitate signiﬁcantly increased computational time. For these reasons, it is generally recommended to use 8 to 20 elements per wavelength [13,18,19], and 10 to 20 points per cycle of the highest frequency for the time step [13,20]. In this study, the time step and the element size were evaluated using Equations (6) and (7). Dt = (6) 20 f max min l = (7) where f is the maximum frequency, l is the element size, and is the shortest wavelength. In the max e min case of this simulation, the wavelength was found to be 29 mm. Thus, the estimated values for time step and element size were obtained from Equations (6) and (7) as 0.5 s and 2.9 mm, respectively. In a ﬁnite element analysis, a surface acoustic wave can be produced in an elastic solid by either a localized impact on the surface or by a piezoelectric transducer. The loading condition of the elastic wave excitation in FEA can be achieved either by applying a force on the surface [21,22] or a nodal displacement [13,23]. In a real-life application, the surface acoustic wave is excited on the test specimen by using an angle beam transducer, a normal beam transducer, or a comb transducer. The angle beam transducer is the most ecient technique for producing the surface waves since it requires only one angle for all frequencies. The wedge material of an angle beam transducer with an incident angle of 58.42 produces surface waves which propagate with an angle of 90 on the surface of the test specimen(steel) [24]. Here, for excitation the SAW, a displacement load in Z-direction (perpendicular to the wave propagation direction) at 12 excitation source points was applied. These are located on the top left surface. The loading function was based on a 100-kHz 2-cycle Hanning- window tone burst that was obtained from Equation (8) [7]: Acoustics 2019, 1 799 Acoustics 2019, 2 FOR PEER REVIEW 6 " !# 1 2t 1 2πt z(t) = 1 cos sin(2ft), t 2 [0, T ] (8) z(t) = 1 − cos ∗ sin (2πft). t ∈ [0, T ] (8) 2 T Acoustics 2019, 2 FOR PEER REVIEW 6 2 T where ƒ is the modulation frequency, t is the time, T is the number of cycles divided by the frequency. Where ƒ is the modulation frequency, t is the tim He, TH is the number of cycles divided by the 1 2πt ( ) ( ) Substituting the time step of 0.5 zs, t the = fr 1equency − cos of 100 ∗ sinkHz, 2πft and . several t ∈ [0, Tcy ] cles of 2 into the Equation (8) (8), frequency. Substituting the time step of 0.5 µs, the frequency of 100 kHz , and several cycles of 2 into 2 T the time function of the input signal is obtained, as shown in Figure 2. Figure 2a illustrates that the the Equation (8), the time function of the input signal is obtained, as shown in Figure 2. Figure 2a Where ƒ is the modulation frequency, t is the time, TH is the number of cycles divided by the input signal completes two cycles in 20 s, and after this time, the displacement should be set to be illustrates that the input signal completes two cycles in 20 µs, and after this time, the displacement frequency. Substituting the time step of 0.5 µs, the frequency of 100 kHz, and several cycles of 2 into zero. The total duration time of the simulation, 60 s, was selected based on the initial simulation should be set to be zero. The total duration time of the simulation, 60 µs, was selected based on the the Equation (8), the time function of the input signal is obtained, as shown in Figure 2. Figure 2a results that showed the incident wave reaching the far end of the plate at the time of 60 s. A 42 ns initial simulation results that showed the incident wave reaching the far end of the plate at the time illustrates that the input signal completes two cycles in 20 µs, and after this time, the displacement (0.042 s) time step was taken between 20 and 60 s to record additional data of the wave amplitude; of 60 µs. A 42 ns (0.042 µs) time step was taken between 20 and 60 µs to record additional data of the should be set to be zero. The total duration time of the simulation, 60 µs, was selected based on the hence, an accurate comparison can be achieved. wave init amip alit l sim udu e; hence lation res , an ults t accur hat showed ate com tp hari e inci son c dent an b wae ve re achachin ieved. g the far end of the plate at the time of 60 µs. A 42 ns (0.042 µs) time step was taken between 20 and 60 µs to record additional data of the wave amplitude; hence, an accurate comparison can be achieved. (a) (b) (a) (b) Figure 2. (a) The time domain of the excitation signal; (b) the frequency spectrum of the excitation Figure 2. (a) The time domain of the excitation signal; (b) the frequency spectrum of the excitation Figure 2. (a) The time domain of the excitation signal; (b) the frequency spectrum of the excitation signal obtained by FFT. signal obtained by FFT. signal obtained by FFT. The mesh parameters were selected to characterize the surface wave propagation on the plate The mesh parameters were selected to characterize the surface wave propagation on the plate The mesh parameters were selected to characterize the surface wave propagation on the plate surface accurately. In all sets of simulations, the hexahedral element type and tetrahedron element type surface sur acc face u acc rateu ly ra . In a tely. In a ll seltl s s of ets simulat of simulat ions, the hex ions, the hexaahedr hedral elem al element type and ent type and tetrahedron e tetrahedron e lement lement were chosen for the plate and the bolt respectively. The mesh element size, 2.9 mm, was evaluated using type were chosen for the plate and the bolt respectively. The mesh element size, 2.9 mm, was type were chosen for the plate and the bolt respectively. The mesh element size, 2.9 mm, was Equation (7) to solve the wave propagation spatially. The mesh quality is a critical aspect in computing evaluated using Equation (7) to solve the wave propagation spatially. The mesh quality is a critical evaluated using Equation (7) to solve the wave propagation spatially. The mesh quality is a critical and solving aspect in co the mput FEA ing an modeldcorr solvin ectly g t , h and e FEA mode it shouldl be correct analyzed ly, and it in each shousimulation. ld be analyzTher ed in each e are thr ee aspect in computing and solving the FEA model correctly, and it should be analyzed in each simulation. There are three key mesh statistics can be checked in ANSYS, which are skewness, key mesh statistics can be checked in ANSYS, which are skewness, orthogonal quality, and element simulation. There are three key mesh statistics can be checked in ANSYS, which are skewness, orthogonal quality, and element quality. In this study, the average value of orthogonal quality, quality. In this study, the average value of orthogonal quality, element quality, and skewness were orthogonal quality, and element quality. In this study, the average value of orthogonal quality, element quality, and skewness were 0.9, 0.84, and 0.1, respectively. The excellent values should be 0.9, 0.84, and 0.1, respectively. The excellent values should be 0.95 for the average orthogonal quality, element quality, and skewness were 0.9, 0.84, and 0.1, respectively. The excellent values should be 0.95 for the average orthogonal quality, and 0.3 for element quality, and be less than 0.25 for the and 0.3 for element quality, and be less than 0.25 for the average skewness [11]. Figure 3a illustrates 0.95 for the average orthogonal quality, and 0.3 for element quality, and be less than 0.25 for the average skewness [11]. Figure 3a illustrates the geometry details and Figure 3b shows the top view the geometry details and Figure 3b shows the top view of the top surface of the plate including the average skewness [11]. Figure 3a illustrates the geometry details and Figure 3b shows the top view of the top surface of the plate including the location of the receivers, the excitation points, and the location of the receivers, the excitation points, and the reﬂected wave propagation path. of the top surf reflected w ac ae of the pla ve propagatio te i n p nacl th udi . ng the location of the receivers, the excitation points, and the reflected wave propagation path. (a) (b) Figure 3. (a) Schematic representation of the bolted joint geometry; (b) the top view of the geometry, (a) (b) including the receiver locations and the reﬂected wave propagation path. All dimensions are in mm. Acoustics 2019, 2 FOR PEER REVIEW 7 Figure 3. (a) Schematic representation of the bolted joint geometry; (b) the top view of the geometry, including the receiver locations and the reflected wave propagation path. All dimensions are in mm. Acoustics 2019, 1 800 4. Results and Discussion 4. Results and Discussion 4.1. Verification Studies 4.1. VThe goa eriﬁcation l oStudies f initial verification studies is to simulate the SAW propagation on a plate without a hole or a bolt and then verify the SAW propagation velocity and its central frequency. For verification The goal of initial veriﬁcation studies is to simulate the SAW propagation on a plate without purposes, three receivers along the direction of wave propagation (X-axis) were located on the plate a hole or a bolt and then verify the SAW propagation velocity and its central frequency. For veriﬁcation surface, and the displacement response of each receiver was recorded. The distance between each purposes, three receivers along the direction of wave propagation (X-axis) were located on the plate receiver (R1, R2, and R3 as illustrated in Figure 1) and the SAW source are 20.22, 40.44, and 60.667 mm, surface, and the displacement response of each receiver was recorded. The distance between each respectively. Figure 4a shows the response of each receiver (displacement in the Z direction) and the receiver (R , R , and R as illustrated in Figure 1) and the SAW source are 20.22, 40.44, and 60.667 mm, 1 2 3 time of zero-crossing for the three signals to verify the phase velocity based on a zero-crossing respectively. Figure 4a shows the response of each receiver (displacement in the Z direction) and algorithm using Equation (9) [25,26]. the time of zero-crossing for the three signals to verify the phase velocity based on a zero-crossing 𝑥 𝑥 algorithm using Equation (9) [25,26]. (9) x x 𝑡 i+1 −𝑡 i ( ) (9) t t im (i+1)m Where ∆𝑥 is the distance between two receivers and ∆𝑡 is the delay time calculated for the where Dx is the distance between two receivers and Dt is the delay time calculated for the zero-crossing zero-cross i ing point in the same phase for the two s im ignals, as illustrated in Figure 4. The results of point in the same phase for the two signals, as illustrated in Figure 4. The results of phase velocity that phase velocity that are obtained from Equation (10) and the theoretical phase velocity of the incident are obtained from Equation (10) and the theoretical phase velocity of the incident signal calculated signal calculated from Equation (3) are listed in Table 2. It can be observed from the Table 2 that the from Equation (3) are listed in Table 2. It can be observed from the Table 2 that the error between error between the theoretical and simulated phase velocity is reasonably small (less than 0.4% for all the theoretical and simulated phase velocity is reasonably small (less than 0.4% for all cases studied). cases studied). Furthermore, to verify the frequency of SAW, fast Fourier transform (FFT) of the Furthermore, to verify the frequency of SAW, fast Fourier transform (FFT) of the obtained signal was obtained signal was taken to convert the time domain of the signal into a frequency domain [27]. taken to convert the time domain of the signal into a frequency domain [27]. Figure 4b illustrates the Figure 4b illustrates the frequency domain of the signal response at R3. It is worth noting that the frequency domain of the signal response at R . It is worth noting that the central frequency is almost central frequency is almost 100 kHz, which is in 3 a good agreement with the central frequency of the 100 kHz, which is in a good agreement with the central frequency of the incident wave. incident wave. (a) (b) Figure 4. (a) Displacement responses in Z-axis for the three receivers R , R , and R as illustrated in 1 2 3 Figure 4. (a) Displacement responses in Z-axis for the three receivers R1, R2, and R3 as illustrated in Figure 1; (b) Frequency spectrum for R3 in Z-axis obtained by FFT. Figure 1; (b) Frequency spectrum for R3 in Z-axis obtained by FFT. Table 2. Phase Velocity Veriﬁcation. Table 2. Phase Velocity Verification. Dx C = Theoretical Value (m/s) Simulated Value (m/s) Error (%) phase Dt im ∆𝒙 Theoretical Value Simulated Value Error R1&R2 2917.58 2905.6 0.4 𝑪 = 𝒑𝒉𝒂𝒔𝒆 ∆𝒕 (m/s) (m/s) (%) R2&R3 2917.58 2926.48 0.3 R1&R3 2917.58 2916.01 0.05 −0.4 R1&R2 2917.58 2905.6 0.3 R2&R3 2917.58 2926.48 4.2. A Comparison Between No Hole and Fully-Loosened Bolt −0.05 R1&R3 2917.58 2916.01 As discussed in the previous section, two sets of simulation studies were conducted: (a) A plate with no hole or bolt; (b) a bolted plate having all components of the bolted joint. In the bolted plate simulation, the four contacts between the bolt and plate were removed to simulate the fully loosen bolt case. Figure 5 illustrates the displacement response in the Z-axis for both cases under the same 𝒊𝒎 Acoustics 2019, 2 FOR PEER REVIEW 8 4.2. A Comparison Between No hole and Fully-loosened Bolt As discussed in the previous section, two sets of simulation studies were conducted: (a) A plate with no hole or bolt; (b) a bolted plate having all components of the bolted joint. In the bolted plate Acoustics 2019, 1 801 simulation, the four contacts between the bolt and plate were removed to simulate the fully loosen bolt case. Figure 5 illustrates the displacement response in the Z-axis for both cases under the same boundary conditions. These responses were recorded using one receiver near the SAW source. It can boundary conditions. These responses were recorded using one receiver near the SAW source. It can be observed in Figure 5b that there is a noticeable dierence between the two amplitudes, particularly be observed in Figure 5b that there is a noticeable difference between the two amplitudes, particularly between the time of 40 and 60 s due to the reﬂected wave propagation from the hole toward the between the time of 40 and 60 µs due to the reflected wave propagation from the hole toward the source in case of the entirely-loosened bolt. Expectedly, this reﬂection is not observed for the plate source in case of the entirely-loosened bolt. Expectedly, this reflection is not observed for the plate without the hole in the design. without the hole in the design. (a) (b) Figure 5. (a) Displacement responses in Z-axis for no hole and fully-loosened bolt cases; (b) Zoomed- Figure 5. (a) Displacement responses in Z-axis for no hole and fully-loosened bolt cases; (b) Zoomed-in in view between the time of 30 and 60 µs for no hole and fully-loosened bolt cases. view between the time of 30 and 60 s for no hole and fully-loosened bolt cases. 4.3. Detection of Wave Reﬂection 4.3. Detection of Wave Reflection Based on the visual observation of incident wave which propagates in the positive x-direction Based on the visual observation of incident wave which propagates in the positive x-direction after the excitation, the major reﬂection, which propagates in the opposite direction of the original after the excitation, the major reflection, which propagates in the opposite direction of the original wave propagation (X direction), occurs when the incident wave completely passes the bolt head wave propagation (−X direction), occurs when the incident wave completely passes the bolt head approximately at the time stamp of 42 s. However, detecting the exact time of the major reﬂection approximately at the time stamp of 42 µs. However, detecting the exact time of the major reflection in the time history of the wave signal for any receiver requires placing one additional receiver -R4- in the time history of the wave signal for any receiver requires placing one additional receiver -R4- before the receiver -R5- with a very short distance between the two receivers as illustrated in Figure 3. before the receiver -R5- with a very short distance between the two receivers as illustrated in Figure The inclusion of the additional receiver in the analysis not only enables distinguishing between 3. The inclusion of the additional receiver in the analysis not only enables distinguishing between the the incident wave and the reﬂected wave from the bolt, but also enables implementing an accurate incident wave and the reflected wave from the bolt, but also enables implementing an accurate comparison when dierent preload values are applied to the bolt. Figure 6 shows the time history of comparison when different preload values are applied to the bolt. Figure 6 shows the time history of the wave at the two receivers. It can be noted that between time 15 and 46 s, the wave passes the R5 the wave at the two receivers. It can be noted that between time 15 and 46 µs, the wave passes the R5 after passing the R4. After 46 s the wave passes R5 earlier than R4, which indicates that the wave after passing the R4. After 46 µs the wave passes R5 earlier than R4, which indicates that the wave is is traveling in the opposite direction (due to reﬂection). It can also be observed that there is a slight traveling in the opposite direction (due to reflection). It can also be observed that there is a slight increase in the amplitude of the primary receiver at time 35 s, which might be due to the reﬂection increase in the amplitude of the primary receiver at time 35 µs, which might be due to the reflection from the bolt which adds with the incoming wave constructively. from the bolt which adds with the incoming wave constructively. Acoustics 2019, 1 802 Acoustics 2019, 2 FOR PEER REVIEW 9 Figure 6. Time history of the 4 kN signal that passes two neighboring receivers which are located Figure 6. Time history of the 4 kN signal that passes two neighboring receivers which are located before the bolt along the bolt center. before the bolt along the bolt center. 4.4. Particle Motion Along the Reﬂected Wave Propagation Path 4.4. Particle Motion Along the Reflected Wave Propagation Path The goal of the reﬂected wave particle motion investigation is to understand the wave reﬂection The goal of the reflected wave particle motion investigation is to understand the wave reflection behavior with varying the preload on the bolt. As previously discussed, the major reﬂection takes place behavior with varying the preload on the bolt. As previously discussed, the major reflection takes after time 42 s. Therefore, this study created a 41 mm-long path. This path is located between the bolt place after time 42 µs. Therefore, this study created a 41 mm-long path. This path is located between and the source, consisting of 50 segments as shown in Figure 3b. The direction of this path is in the the bolt and the source, consisting of 50 segments as shown in Figure 3b. The direction of this path is negative X-direction (reﬂected wave direction). The particle movement of the path in Z-direction for in the negative X-direction (reflected wave direction). The particle movement of the path in Z- the ﬁve preload values was recorded from time 42 to 50 s. Figure 7a shows the wave patterns reﬂected direction for the five preload values was recorded from time 42 to 50 µs. Figure 7a shows the wave from the bolt at the time stamp of 46 s. The displacement of the path at 46 s was chosen for a better patterns reflected from the bolt at the time stamp of 46 μs. The displacement of the path at 46 µs was illustration of how the reﬂection varies as the preload value increases, as shown in Figure 7b. It can chosen for a better illustration of how the reflection varies as the preload value increases, as shown be observed from this ﬁgure that as the preload increases, the width of the reﬂection wave increased. in Figure 7b. It can be observed from this figure that as the preload increases, the width of the It can also be observed that the reﬂection amplitude increases as the preload increases. There are four reflection wave increased. It can also be observed that the reflection amplitude increases as the zero-crossing points along this path. The last zero-crossing point was chosen for comparison as it preload increases. There are four zero-crossing points along this path. The last zero-crossing point was the nearest point to the source. The distance of zero-crossing from the source for each preload was chosen for comparison as it was the nearest point to the source. The distance of zero-crossing is listed in Table 3. It can be observed that the zero-crossing of the higher preload of 24 kN has the from the source for each preload is listed in Table 3. It can be observed that the zero-crossing of the shortest distance to the source, while the lower preload of 4 kN has the longest distance to the source. higher preload of 24 kN has the shortest distance to the source, while the lower preload of 4 kN has The dierence in the distance between the lower and higher preload is approximately 1.9 mm. Figure 7c the longest distance to the source. The difference in the distance between the lower and higher illustrates the average of the reﬂected wave amplitude between 42 and 50 s along this path. It can be preload is approximately 1.9 mm. Figure 7c illustrates the average of the reflected wave amplitude observed that the reﬂected wave associated with the higher preload has a higher amplitude average between 42 and 50 µs along this path. It can be observed that the reflected wave associated with the as the wave approaches the source. Table 3 also shows that the relationship between the reﬂection higher preload has a higher amplitude average as the wave approaches the source. Table 3 also shows location and the preload is nonlinear as has been predicted by our earlier study [6]. Furthermore, that the relationship between the reflection location and the preload is nonlinear as has been a slight dierence in the wavelength was observed in all ﬁve preload cases studied as shown in Table 3. predicted by our earlier study [6]. Furthermore, a slight difference in the wavelength was observed Estimating the location of the ﬁrst reﬂected wave has the potential to quantify the bolt tension. in all five preload cases studied as shown in Table 3. Estimating the location of the first reflected wave has the potential to quantify the bolt tension. Acoustics 2019, 1 803 Acoustics 2019, 2 FOR PEER REVIEW 10 (a) (b) (c) Figure 7. (a) The reflected waves from the bolt at time 46 µs; (b) The displacement in Z-direction along Figure 7. (a) The reﬂected waves from the bolt at time 46 s; (b) The displacement in Z-direction along the reflected wave path at time 46 µs before the bolt; (c) The average displacement in Z-direction along the reﬂected wave path at time 46 s before the bolt; (c) The average displacement in Z-direction along the path between time 42 and 50 µs. the path between time 42 and 50 s. Table 3. The distance of zero-crossing of each preload from the source and the corresponding wavelength. All dimensions in mm. Preload (kN) 4 8 12 16 20 24 Torque (N.m) 7.2 14.4 21.6 28.8 36 43.3 Distance from the source 17.4 16.7 16.4 16.2 15.9 15.5 Wavelength () 22.31 22.25 25.23 25.12 25.23 25.5 Acoustics 2019, 2 FOR PEER REVIEW 11 Table 3. The distance of zero-crossing of each preload from the source and the corresponding wavelength. All dimensions in mm. Preload (kN) 4 8 12 16 20 24 Torque (N.m) 7.2 14.4 21.6 28.8 36 43.3 Distance from the source 17.4 16.7 16.4 16.2 15.9 15.5 Acoustics 2019, 1 804 25.23 25.12 25.23 25.5 Wavelength (λ) 22.31 22.25 4.5. Time History of the Reflected Wave Near the Bolt 4.5. Time History of the Reﬂected Wave Near the Bolt A comparison of varying the preload on the SAW propagation and reflection wave A comparison of varying the preload on the SAW propagation and reﬂection wave characteristics characteristics involves three necessary steps: (1) Recording the time history of the signal for each involves three necessary steps: (1) Recording the time history of the signal for each preload at the preload at the receiver -R5-, as illustrated in Figure 3b; (2) locating the reflection zone which was receiver -R5-, as illustrated in Figure 3b; (2) locating the reﬂection zone which was found to be after found to be after 42 µs ; (3) measuring the difference using one preload case as a reference. Figure 8 42 s; (3) measuring the dierence using one preload case as a reference. Figure 8 shows the results of shows the results of the three steps followed to investigate the effect of preload on wave propagation the three steps followed to investigate the eect of preload on wave propagation and reﬂection. It can and reflection. It can be observed from Figure 8b that as the preload increases, the reflection be observed from Figure 8b that as the preload increases, the reﬂection amplitude increases. For this amplitude increases. For this reason, the 4 kN case was taken as a reference for measuring the reason, the 4 kN case was taken as a reference for measuring the dierence in the third step. It can difference in the third step. It can also be noted that the phase of reflection wave shifts to the left as also be noted that the phase of reﬂection wave shifts to the left as the preload increases, particularly at the preload increases, particularly at 46 µs. As a result, the peak amplitude of a higher preload occurs 46 s. As a result, the peak amplitude of a higher preload occurs earlier than the lower preload, as was earlier than the lower preload, as was also observed in our earlier study [6]. The most considerable also observed in our earlier study [6]. The most considerable dierence in time, 0.5 s, exists between difference in time, 0.5 µs, exists between the higher and the lower preload. Figure 8c shows the the higher and the lower preload. Figure 8c shows the amplitude dierence from 40 to 60 s, which is amplitude difference from 40 to 60 µs, which is obtained from Equation (10) obtained from Equation (10) Di f f erence = A A (10) i 4kN 𝑟𝑒𝑛𝑐𝑒𝐷𝑖𝑓𝑓𝑒 = 𝐴 − 𝐴 (10) where A is the amplitude of each preload, A the reference amplitude (4 kN). It can be noted in i iis the amplitude of each preload, A4kN 4kN the reference amplitude (4 kN). It can be noted in where A Figure 8c that the maximum dierence is found at time 46 s, where the phase shifting of reﬂected Figure 8c that the maximum difference is found at time 46 µs, where the phase shifting of reflected wave occurs. wave occurs. (a) (b) (c) Figure 8. (a) Time history of the propagating wave for all preload values; (b) Zoomed-in view of the reﬂection zone, located between 40 and 60 s; (c) The dierential representation of all preload cases with 4 kN case as the reference. Acoustics 2019, 2 FOR PEER REVIEW 12 Figure 8. (a) Time history of the propagating wave for all preload values; (b) Zoomed-in view of the reflection zone, located between 40 and 60 µs; (c) The differential representation of all preload cases with 4 kN case as the reference. 4.6. Simulation and Experimental Comparison In our previous experimental study, the authors employed a 5 MHz ultrasonic transducer which was connected to a judiciously designed ultrasonic wedge to create surface waves propagating on Acoustics 2019, 1 805 the steel test specimen [6]. In order to obtain a higher resolution ultrasonic image, a synthetic phased array system was utilized. This system is simply formed by precisely moving the transducer linearly 4.6. Simulation and Experimental Comparison with a known pitch. The synthetic phased array system simulated a total of 50 elements in an array In our previous experimental study, the authors employed a 5 MHz ultrasonic transducer which which was aimed to generate an acoustic beam. The bolted joints used in the experimental setup were was connected to a judiciously designed ultrasonic wedge to create surface waves propagating on the a 12.7 mm thick 1018 steel plate, a 6.4 mm grade 8 bolt and a 19 mm stainless steel washer. The steel test specimen [6]. In order to obtain a higher resolution ultrasonic image, a synthetic phased transducer was placed on the top surface of 1018 steel bolted joint from a distance of 29 mm from the array system was utilized. This system is simply formed by precisely moving the transducer linearly with a known pitch. The synthetic phased array system simulated a total of 50 elements in an array bolt. During the experiments, five bolt tension values, similar to those presented in this study, were which was aimed to generate an acoustic beam. The bolted joints used in the experimental setup applied to the bolt for comparison. The bolt tension values are indirectly related to the applied torque were a 12.7 mm thick 1018 steel plate, a 6.4 mm grade 8 bolt and a 19 mm stainless steel washer. as expressed in Equations 4 and 5. By analyzing the reconstructed ultrasonic images of SAWs for the The transducer was placed on the top surface of 1018 steel bolted joint from a distance of 29 mm from different torque val the bolt. During ues, it the wa experiments, s observed tha ﬁve boltttension as the values, applied t similar oto rqthose ue incre presented ased in , t this he dist study,ance of t were he first applied to the bolt for comparison. The bolt tension values are indirectly related to the applied torque reflected wave from the bolted joint received by the transducer decreased. It should be highlighted as expressed in Equations 4 and 5. By analyzing the reconstructed ultrasonic images of SAWs for the that the same finding was observed in this simulation study presented here, as discussed in section dierent torque values, it was observed that as the applied torque increased, the distance of the ﬁrst 4.4. Figure 9 illustrates the experimentally obtained distance of the first reflected wave from the reﬂected wave from the bolted joint received by the transducer decreased. It should be highlighted that the same ﬁnding was observed in this simulation study presented here, as discussed in Section 4.4. transducer as a function of the applied torque. As shown in Figure 9, the longest distance of 27.9 mm Figure 9 illustrates the experimentally obtained distance of the ﬁrst reﬂected wave from the transducer observed corresponds to the lowest torque value (3.61 N.m) from the transducer, while the shortest as a function of the applied torque. As shown in Figure 9, the longest distance of 27.9 mm observed distance of 27 mm is for a torque value of 10.31 N.m [6]. corresponds to the lowest torque value (3.61 N.m) from the transducer, while the shortest distance of 27 mm is for a torque value of 10.31 N.m [6]. 27.8 27.6 27.4 27.2 26.8 3.5 4.5 5.5 6.5 7.5 8.5 9.5 10.5 Applied Torque (N.m) Figure 9. The experimentally obtained distance of the ﬁrst reﬂected wave from the transducer as a function of the applied torque. Figure 9. The experimentally obtained distance of the first reflected wave from the transducer as a function of the applied torque. In order to make a fair comparison between the simulation and experiment results, three torque values were normalized to the maximum applied torque for each case. For instance, the three selected torque values from the experiment are 3.61, 7, and 10.39 N.m. Normalizing these values to the maximum In order to make a fair comparison between the simulation and experiment results, three torque torque applied (10.39 N.m), resulted in 34%, 67%, and 100%. The three equivalent torque values from values were normalized to the maximum applied torque for each case. For instance, the three selected the simulation that lead to the same percentages are 14.4, 28.8, and 43.3 N.m. The simulated and torque values from the experiment are 3.61, 7, and 10.39 N.m .Normalizing these values to the experimentally observed shift from the maximum allowable torque are listed in Table 4. The values in Table 4 are obtained by subtracting the positions of first reflected wave for each applied torque that maximum torque applied (10.39 N.m), resulted in 34%, 67%, and 100%. The three equivalent torque are shown in Table 3 and Figure 9 from the position of the first reflected wave associated with the values from the simulation that lead to the same percentages are 14.4, 28.8, and 43.3 N.m. The maximum torque. It can be observed that there is an excellent agreement between the simulation and the simulated and experimentally observed shift from the maximum allowable torque are listed in Table experimental results with just a few millimeter differences between the two studies. It should be pointed out that the discrepancy between the simulation and experiment results can be due to several factors, 4. The values in Table 4 are obtained by subtracting the positions of first reflected wave for each such as the difference in the wave frequency studied, the bolt size, the bolt material, or the plate thickness. applied torque that are shown in Table 3 and Figure 9 from the position of the first reflected wave associated with the maximum torque. It can be observed that there is an excellent agreement between the simulation and the experimental results with just a few millimeter differences between the two studies. It should be pointed out that the discrepancy between the simulation and experiment results can be due to several factors, such as the difference in the wave frequency studied, the bolt size, the bolt material, or the plate thickness. Distance from transducer (mm) Acoustics 2019, 1 806 Table 4. The simulated and experimentally shift from the maximum allowable torque. 34% 67% 100% Simulated shift from the max. allowable (mm) 1.2 0.7 0 Experimental shift from the max. allowable (mm) 0.9 0.8 0 5. Conclusions In this study, the eect of varying bolt tension on SAW was simulated via the commercially available ﬁnite element program, ANSYS 18.1. The study investigated the SAW reﬂection caused by the real contact area, which is located at the interface between the bolt head and plate surface. To ensure the validity of the FEA model developed, the phase velocity of the incident wave (SAW) and frequency were simulated and then compared with the theoretical values. The error in the phase velocity and frequency were 0.4% and negligible, respectively. The results show a clear trend of an increase in the reﬂected wave amplitude as the bolt tension are varied from 4 to 24 kN with a slight phase shifting of 0.5 s. The results further show that as the bolt tension increases, the distance of the ﬁrst reﬂected wave to the source becomes shorter. This was measured based on the nearest zero-crossing point to the source. The most signiﬁcant dierence in distance, 1.9 mm, was observed between 4 and 24 kN. Most importantly, the results also indicate that the bolt tension can be reliably quantiﬁed based on the location of the ﬁrst reﬂected wave. The simulation results show an excellent agreement with the experiment results. The results presented here have the potential to signiﬁcantly impact bolt tension quantifying and ensuring that all bolts are properly tightened for a large engineering structure or a machine. Author Contributions: Investigation, H.A. and R.G.; methodology, R.G.; project administration, R.G.; resources, R.G.; supervision, R.G.; writing—original draft, H.A.; writing—review & editing, R.G. Acknowledgments: The authors gratefully acknowledge the Research Computing Department at the University of South Florida for providing the computer resources which facilitate conducting the simulation in this paper. Conﬂicts of Interest: The authors declare no conﬂicts of interest. References 1. Sakai, T. Bolted Joint Engineering: Fundamentals and Applications; Beuth Verlag: Berlin, Germany, 2008. 2. Jhang, K.-Y.; Quan, H.-H.; Ha, J.; Kim, N.-Y. Estimation of clamping force in high-tension bolts through ultrasonic velocity measurement. Ultrasonics 2006, 44, e1339–e1342. [CrossRef] [PubMed] 3. Wang, T.; Song, G.; Liu, S.; Li, Y.; Xiao, H. Review of bolted connection monitoring. Int. J. Distrib. Sens. Netw. 2013, 9, 871213. [CrossRef] 4. Chaki, S.; Corneloup, G.; Lillamand, I.; Walaszek, H. Combination of longitudinal and transverse ultrasonic waves for in situ control of the tightening of bolts. J. Press. Vessel Technol. 2007, 129, 383–390. [CrossRef] 5. Wang, T.; Song, G.; Wang, Z.; Li, Y. Proof-of-concept study of monitoring bolt connection status using a piezoelectric based active sensing method. Smart Mater. Struct. 2013, 22, 087001. [CrossRef] 6. Martinez, J.; Sisman, A.; Onen, O.; Velasquez, D.; Guldiken, R. A synthetic phased array surface acoustic wave sensor for quantifying bolt tension. Sensors 2012, 12, 12265–12278. [CrossRef] 7. Giurgiutiu, V. Structural Health Monitoring with Piezoelectric Wafer Active Sensors: With Piezoelectric Wafer Active Sensors; Elsevier: Amsterdam, The Netherlands, 2007. 8. Rose, J.L. Ultrasonic Guided Waves in Solid Media; Cambridge University Press: Cambridge, UK, 2014. 9. White, R.M. Surface elastic waves. Proc. IEEE 1970, 58, 1238–1276. [CrossRef] 10. ANSYS®. Academic Research Mechanical Release 18.1, Static Structural, Engineering Data; ANSYS, Inc.: Canonsburg, PA, USA, 2019. 11. Padilla, S.; Tufekcioglu, E.; Guldiken, R. Simulation and veriﬁcation of polydimethylsiloxane (PDMS) channels on acoustic microﬂuidic devices. Microsyst. Technol. 2018, 24, 3503–3512. [CrossRef] 12. Duan, W.; Niu, X.; Gan, T.-H.; Kanfoud, J.; Chen, H.-P. A Numerical Study on the Excitation of Guided Waves in Rectangular Plates Using Multiple Point Sources. Metals 2017, 7, 552. [CrossRef] Acoustics 2019, 1 807 13. Moser, F.; Jacobs, L.J.; Qu, J. Modeling elastic wave propagation in waveguides with the ﬁnite element method. NDT E Int. 1999, 32, 225–234. [CrossRef] 14. Budynas, R.G.; Nisbett, J.K. Shigley’s Mechanical Engineering Design; McGraw-Hill: New York, NY, USA, 2008; Volume 8. 15. Fastenal Technical Reference Guide. Available online: https://www.fastenal.com/content/documents/ FastenalTechnicalReferenceGuide.pdf. (accessed on 13 March 2019). 16. Vandermaat, D. Factors Aecting Pre-Tension and Load Carrying Capacity in Rockbolts—A Review of Fastener Design. In Proceedings of the 18th Coal Operators’ Conference, Wollongong, Australia, 7–8 February 17. ANSYS®. Academic Research Mechanical, Products/Structures/Strength Analysis/Simulating Bolted Assemblies; ANSYS Inc.: Canonsburg, PA, USA, 2019; Available online: https://www.ansys.com/products/structures/ strength-analysis/simulating-bolted-assemblies (accessed on 20 March 2019). 18. Alleyne, D.; Cawley, P. A two-dimensional Fourier transform method for the measurement of propagating multimode signals. J. Acoust. Soc. Am. 1991, 89, 1159–1168. [CrossRef] 19. Yasuda, T.; Pang, B.; Nishino, H.; Yoshida, K. Dynamic behavior evaluation of martensitic transformation in Cu-Al-Ni shape memory alloy using acoustic emission simulation by FEM. Mater. Trans. 2011, 52, 397–405. [CrossRef] 20. Balvantín, A.J.; Diosdado-De-la-Peña, J.A.; Limon-Leyva, P.A.; Hernández-Rodríguez, E. Study of guided wave propagation on a plate between two solid bodies with imperfect contact conditions. Ultrasonics 2018, 83, 137–145. [CrossRef] [PubMed] 21. Fierro, G.P.M.; Ciampa, F.; Ginzburg, D.; Onder, E.; Meo, M. Nonlinear ultrasound modeling and validation of fatigue damage. J. Sound Vib. 2015, 343, 121–130. [CrossRef] 22. Shen, Y.; Giurgiutiu, V. Eective non-reﬂective boundary for Lamb waves: Theory, ﬁnite element implementation, and applications. Wave Motion 2015, 58, 22–41. [CrossRef] 23. Aiello, G.; Dilettoso, E.; Salerno, N. Finite element analysis of elastic transient ultrasonic wave propagation for NDT applications. In Proceedings of the 5th WSEAS/IASME International Conference on Systems Theory and Scientiﬁc Computation, Malta, 15–17 September 2005; World Scientiﬁc and Engineering Academy and Society (WSEAS): Stevens Point, WI, USA; pp. 114–119. 24. Rose, J.L. Ultrasonic Waves in Solid Media; Cambridge University Press: Cambridge, UK, 2004. 25. Draudviliene, L.; Aider, H.A.; Tumsys, O.; Mazeika, L. The Lamb waves phase velocity dispersion evaluation using an hybrid measurement technique. Compos. Struct. 2018, 184, 1156–1164. [CrossRef] 26. Mažeika, L.; Draudviliene, ˙ L. Analysis of the zero-crossing technique in relation to measurements of phase velocities of the Lamb waves. Ultragarsas Ultrasound 2010, 65, 7–12. 27. Lynch, K.M.; Marchuk, N.; Elwin, M.L. Chapter 22—Digital Signal Processing. In Embedded Computing and Mechatronics with the PIC32; Lynch, K.M., Marchuk, N., Elwin, M.L., Eds.; Newnes: Oxford, UK, 2016; pp. 341–374. © 2019 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

Acoustics – Multidisciplinary Digital Publishing Institute

**Published: ** Sep 27, 2019

**Keywords: **bolted joint; surface acoustic wave; bolt tension; finite element modeling; the real area of contact

Loading...

You can share this free article with as many people as you like with the url below! We hope you enjoy this feature!

Read and print from thousands of top scholarly journals.

System error. Please try again!

Already have an account? Log in

Bookmark this article. You can see your Bookmarks on your DeepDyve Library.

To save an article, **log in** first, or **sign up** for a DeepDyve account if you don’t already have one.

Copy and paste the desired citation format or use the link below to download a file formatted for EndNote

Access the full text.

Sign up today, get DeepDyve free for 14 days.

All DeepDyve websites use cookies to improve your online experience. They were placed on your computer when you launched this website. You can change your cookie settings through your browser.