Next Article in Journal
Multimodal Gold Nanostars as SERS Tags for Optically-Driven Doxorubicin Release Study in Cancer Cells
Next Article in Special Issue
Physical Properties of an Eco-Sustainable, Form-Stable Phase Change Material Included in Aerial-Lime-Based Mortar Intended for Different Climates
Previous Article in Journal
Dilatometric Analysis and Kinetics Research of Martensitic Transformation under a Temperature Gradient and Stress
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Experimental and Numerical Analysis of Multiple Low-Velocity Impact Damages in a Glass Fibered Composite Structure

by
Kaleeswaran Balasubramaniam
1,
Dominika Ziaja
2,
Michał Jurek
2,*,
Piotr Fiborek
1 and
Paweł Malinowski
1
1
Department of Mechanics of Intelligent Structures, Institute of Fluid Flow Machinery, Polish Academy of Sciences, Fiszera 14, 80-231 Gdansk, Poland
2
Department of Structural Mechanics, Rzeszow University of Technology, Al Powstańców Warszawy 12, 35-959 Rzeszow, Poland
*
Author to whom correspondence should be addressed.
Materials 2021, 14(23), 7268; https://doi.org/10.3390/ma14237268
Submission received: 26 October 2021 / Revised: 17 November 2021 / Accepted: 25 November 2021 / Published: 28 November 2021
(This article belongs to the Special Issue Innovative Materials for Construction II)

Abstract

:
Glass fiber-reinforced polymer structures (GFRPS) are widely used in civil and mechanical fields due to their light weight and corrosion resistance. However, these structures are prone to damage with very-low-energy impacts. The reliability of such structures is of prime importance before their installation and usage. This study aimed to identify, visualize, localize, and verify multiple barely visible impact damage (BVID) in a GFRPS using a combination of guided waves (GW)-based online structural health monitoring (SHM) and thermal strain-based nondestructive testing (NDT) approaches. Global NDT techniques like the use of a laser Doppler vibrometer (LDV) and digital image correlation (DIC) were used in the experimental analysis. The effectiveness of the experimental LDV-GW process was also checked numerically with the spectral element method (SEM). A threshold-based baseline free SHM approach to effectively localize the damages was proposed along with quick DIC verification of composite structure with thermal loading based on short-pulse heating as an excitation source. This study analyzed combined experimental- and numerical-based SHM-NDT methods in characterizing the multiple BVIDs located in a GFRPS.

1. Introduction

Composites are light weight, have high strength, and are used in many automotive, civil, aerospace, etc. applications. The GFRPS laminates have greater strength along the direction of fiber [1] GFPRS are widely used in retrofitting structures in construction industries, and, nowadays, these glass fibers are combined to prepare reinforcement bars for structural strengthening [2]. They are highly preferred due to their high tensile strength, low density, and corrosion resistance [3]. Even though several advantages exist, they also have disadvantages like delamination defects, impact cracks [4], internal fiber matrix damage [5], etc. If such damage fails to be noticed, it could lead to severe structural disintegration [6]. BVID and impact cracks in the GFRPS are of most common that occurs even due to a low-energy impact force and are mostly not visible for internal visual identification.
The engineering industry offers many methods for analyzing the defects in the structures. The NDT analysis of glass fibers using a combination of X-ray tomography, ultrasonic imaging, and optical proliferometry was proposed by Miqoi et al. [7] and on carbon fibers using digital image correlation (DIC) and scanning electron microscopy by Sarasini et al. [8], in addition to experimental compression after impact tests by Talreja et al. [9]. SHM finds its way in detecting such impact damage by performing continuous health monitoring of structures 1. SHM reduces the laborious human hours of checking the whole structure for damage with the help of fewer piezoelectric lead zirconate titanate (PZT) to analyze the structure. These PZTs send and receive the signals and thus help in online monitoring.
GW are of recent trends in SHM applications due to their higher area coverage in a short interval of time. The GW analyzes the damage based on differences in the attenuation parameters [4] The GWs are classified into symmetric (S) and antisymmetric (A) waves based on their velocity, varying propagation patterns, etc. The slower GW mode is the A0 mode in lower frequencies, and the faster one is in the S0 mode [10,11]. The electrical input applied to the PZT results in the excitation of GW signals, which could be used for online monitoring of structures [1].
The impact damage on glass fibers was analyzed using wavelet packet transform [12]. Impact damage sensing based on fast Fourier transform peak analysis, nonlinear sideband peak variations, and time-of-flight analysis were performed on GFRP plates by Alnuaimi et al. [13]. The impact effects in the GFRP structure were studied using GW amplitude variation due to defect size, fiber breakage [14], and phase difference [15] with a circular PZT network. Zhang et al. used differences in mode conversion behaviour to analyze the glass fibered structure at different stages of impact [16]. LDV-based root mean square [17] and radially weighted root mean square (RRMS) formulations [18] were used by many in analyzing the entire area of composite structures. Numerical-based SEM codes were also used in studying impact damage [19]. A combination of experimental and numerical methods using the SEM was used in detecting the impacts in composite structures [20]. A comparison of BVID identification was performed [21] using the signal-to-noise ratio damage index method and the probabilistic damage algorithm.
The DIC method uses the analysis of changes in the random speckles pattern coated on the surface of the examined object, and the observation of displacements or strain field variation is possible. In many cases, these changes inform about damage appearance before the defects are visible [22]. An additional advantage of DIC is the possibility of simultaneous measurement of many points of the observed area; for this reason, this method was chosen as a verification technique in this study [23].
In many DIC applications, the specimens were tested in universal testing machines, where the test finishes with specimen destruction [24,25]. The main problem of applying DIC as a nondestructive method is the determination of the way the sample is loading such that the observed displacement/strain field is sufficiently large concerning the accuracy of the measurement and, at the same time, fully reversible. The short heating process was selected and proposed as an excitation in the presented approach to reduce the verification time. The effectiveness of the DIC technique for full-field high-temperature thermal deformation measurement has been demonstrated, e.g., in [26] where the authors examined the thermal properties of austin stainless steel. The thermal effects in composite specimens using DIC were examined by Foti et al. [22], but the authors examined the high-temperature fatigue. The thermal buckling of the circular composite plate was analyzed based on DIC measurement by Jin et al. [27] and Zhu et al. [28] who used DIC to observe the cooling process of turbine vane thermal barrier coatings after thermal shock. At the same time, Rajaram et al. [29] showed the utility of DIC to surface damage detection of aluminum specimens based on the coefficient of thermal expansion.
Most research works focus on the localization of single-impact damage during and after the impact [5,8,10,12,15]. There are fewer research works on analyzing and localizing the multiple impacts and on verification of the deformed areas after impact. In our research work, we focused not only on identification, visualization, and localization of the multiple linearly located impact damages but also studies were carried out to observe changes (verification) in the impact-deformed zones by applying a modified [23] quick thermal strain field technique.
A numerical method based on SEM code was developed and compared with the experimental LDV signals to analyze the damage. Since all three damages existed before, the numerical simulation was utilized to study each damage scenario separately. The RRMS method was applied to the signal data obtained from both LDV and SEM. The SHM-based baseline free elliptical threshold method (ETM) was proposed to show the effective online localization of the damages. Thus the proposed SHM-GW technique can be used as the first step of the quick analysis to locate the damages followed by NDT-DIC verification based on SHM results. The SHM combined with the NDT method can drastically reduce the time to analyze larger and complex structures.

2. Methodology, Material, and Approaches

2.1. Methodological Flow Chart

The flowchart (Figure 1) briefly explains the methodological process involved in analyzing the GFRPS using combined NDT and SHM methods. The process involves damage-identification, visualization, localization, and verification.

2.2. GFRP Material Studied

The 0.2 cm thick GFRPS has 12 layups stacked in [0/90]3s format, which are bonded together by HEXCEL-212Na adhesive. The material is procured from the G.ANGELONI group [30] (Venice, Italy). These materials are largely used in the mechanical and construction industries. The GFRPS is of dimensions 50 × 50 cm2 in length and breadth. Multiple BVIDs of Ø 0.9 cm (BVID1) and Ø 1.2 cm (BVID2 and BVID3) are made with an impact force of 20, 25, and 30J. This is done by dropping a steel ball attached to a steel beam from varying heights guided via a pipe. All the 4 PZTs (Ø 1 cm and 0.1 cm thin) are made from SONOX material P502 [31] and are attached to the GFRPS using cyanoacrylate glue. The PZTs (S1–S4) and BVIDs location coordinates are presented in Figure 2.

2.3. LDV Setup

A single head (1D) LDV was chosen, which analyzed the structures using out-of-plane wavefield components. Polytec PSV 400 was the LDV type used. LDV helps to visualize the GW propagation and also to acquire the GW full-wave fields. The GW propagation was done by exciting the PZT. The retro-reflective film was applied only to an area of 26 cm × 50 cm (blue dotted lines), as shown in Figure 2, to reduce the LDV calculation time. The LDV setup consisted of the signal generator, the power amplifier, and the oscilloscope. A 5-cycle sinusoidal tone burst generated with a Hanning window was applied to the PZT. Frequencies of 50 kHz, 100 kHz, and 200 kHz were used in this analysis. The LDV measurement at each grid point was time-averaged 5 times with an overall sampling period of 1024. A 16 Vpp voltage with a gain of 20 was applied to PZT via the generator. All four symmetrical positioned PZTs were excited via this process.

2.4. SEM Numerical Approach

The simulation was conducted with the algorithm developed in Matlab software (version 2020a) based on the implementation found in [32]. The GFRP structure was modelled according to the Mindlin–Reissner first-order shear deformation theory [33], and the laminate theory [34] was used to determine the effective engineering constants shown in Table 1. The details of the determination of the effective constants are given in Appendix A. The model of PZT transducers was omitted in the analysis to simplify calculations. Instead, the out-of-plane forces were applied to the nodes in the place of the actuator. The damage was modulated in two ways. First, Young’s modulus of the structure was reduced by 50% in the circular of the damaged area. Second, three cracks were introduced using node separation [35] at the location marked by the dashed line in Figure 3. The LDV experimental setup values, as mentioned earlier, were used in studying the numerical GFRP models.
The four model cases (Figure 3) were made numerically as the structure obtained already had three BVIDs in the initial state itself. Thus, SEM models added an advantage to visualize how the GW propagation would be with the individual BVIDs.

2.5. DIC Setup

As a verification procedure, 3D DIC was proposed; so, the measuring setup consisted of two USB3.0 Baumer 12.3 Mpx cameras (Baumer GmbH, Radeberg, Germany) with VS-1220HV lenses, combined in a Q400 system by DANTEC DYNAMICS GmbH. The strain field changes resulted from heating the GFRP-plate using a halogen lamp HEDLER HF 65 (HEDLER GmbH, Hesse, Germany); the exposition time was 2 min. The size of the observed area was a compromise between the highest possible accuracy and the widest observed area. The area observed by each camera was about 340 mm in width and 249 mm high, which considering the camera resolution of 4096 × 3000 px gave the spatial resolution equal to 83 μm/px. The calibration was made using the Al-15-BMB_9×9 calibration target. During the evaluation of pictures, the facet size was established to 17 pixels. In the presented examination, the mean accuracy of the displacement measurement was: in the horizontal direction, 0.00006 mm; in the vertical direction, 0.00007 mm; and, perpendicular to the observed plane, 0.002 mm. The whole specimen was not observed. The measuring stand is shown in Figure 4a, and the scheme of the experiment is shown in Figure 4b, top view.

3. Damage Analysis Methods

The damage identification and visualization were done using RRMS for both LDV and SEM data (Section 3.1), localized by the developed SHM-based ETM method (Section 3.2) and verified by the proposed DIC thermal speckle-based variation method (Section 3.3).

3.1. RRMS-Based Damage Analysis

The visualization method using RRMS [18] for four PZTs excitations were analyzed. The RMS provides the energy distribution of the signals. However, if more energy gets accumulated in a certain region, it shows a higher value of RMS. It results in blurring of the damage zone sometimes. The RRMS function Equation (1) was added to counteract such cases, which enhanced the results obtained.
R R M S ( x , y ) = i = 1 N ( k i ( x , y , t ) ) 2 N × [ R W F ] a
where k i ith sample amplitude, and N—the number of registered time samples.
RWF is the radially weighted factored value for each scanning point, and “a” is the power of the weight factor. In this function, the excitation points were the origin point (x0, y0), and the radial distance was measured for all the points (xi, yin), as shown in Equation (2).
( R W F ) i = ( x 0 x i ) 2 + ( y 0 y i ) 2
With the help of the obtained Equation (2), a combination of LDV and SEM models were analyzed as mentioned in Table 2.

3.2. Threshold-Based Elliptical Method in Damage Analysis

The localization of the damages was done using the SHM-based strategy using the ETM. It works on the methodology of convergence and interaction of ellipses. The method requires the time of arrival (TOA) of the wave peaks in the damage calculation [36,37]. The time of arrival (theoretical) (TOAT) was obtained between GW paths to the arbitrary grid points, the sensor, and the actuator. After applying the Hilbert transformation (HT) envelope to the signals, amplitude peaks were obtained. The obtained TOAT was then tried to match up with the HT peak values (TOAE).
A higher value in the mesh grid was obtained when the TOAT matched with an HT value. A fine mesh of 0.1 mm spacing was used in the calculation. Since most reflections occurred at the damage spot or near the damage spot, that region was highlighted in the grid. A threshold value (90%) was applied to the algorithm to separate the higher values in the grid during the run. The threshold value was selected based on the experience and the references from the threshold-based probabilistic study [38]. The damage localization approach is explained in a flowchart, as shown in Figure 5.
(G1, G2) are arbitrary grid point coordinates; (A1, A2, S1, S2) are actuator and sensor coordinates; VA, VS represents the group velocities obtained from the actuator to damage and the sensor to damage, respectively; AT represents the amplification factor; TOAE represents the time of arrival obtained from the experiment; D (x, y) represents the damage index; AS represents the total number of actuator sensor paths studied; τ is the decay factor value of 5 μs, and N represents threshold value used in the calculation.

3.3. Thermal Speckle-Based Verification Method

The examined GFRPS was coated by the typical DIC speckle pattern (registered by two cameras as shown in Figure 6), so the damages were invisible. The specimen was observed under the heating process realized by the light source. Since the damage location was initially localized with SHM, the lamps were focused in the localized area. The reference picture was made at stable room temperature, and then the whole heating took two minutes. During this time, each camera registered two new photos every one second (sampling frequency equal to 2 Hz). For the analysis of picture sequences, the commercial ISTRA 4D software (version 4.6) was used.

4. Discussion on the Results Obtained

4.1. RRMS Models-Visualization of Damages

RRMS studies were carried out with the LDV and SEM data to visualize the damage. A constant power/weight factor value of a = 0.5 was used for all the calculations. Table 2’s case results are shown in Figure 7 and Figure 8. Initially, cases 4 and 5 (Figure 7) were considered for the visualization study as they represent the real damage scenarios.
After analyzing all the frequencies, it was found that both LDV_RRMS and SEM_RRMS 200 kHz plots identified all the BVIDs more prominently with rounded red circles (Figure 7). Table 2’s other cases (1–3) were also checked only with 200 kHz. Thus, only 200 kHz was chosen as the identified frequency for other RRMS and SHM results (Section 4.2). The SEM_RRMS plots, as shown in Figure 8, identified the inner numerically made BVIDs (1, 2, and 3) separately and proved that the methodology can be used for single damages.
After checking the results from S1, similar BVID identification and visualization checks were done for other symmetrically located sensors, as shown in Figure 9, using the RRMS-based weightage formulation. All sensors identified the damages clearly both from the experimental and the numerical data.

4.2. ETM-Localization of Damages

The amplification and threshold factor-added elliptical-based algorithm was applied to Table 2 in damage localization. The damage points were identified based on elliptical intersections. Equidistant nodal points (P1–P10) were taken from the LDV area scan results with excitation from S1 (Figure 10a) for the process. Since S2 and S3 are in a straight line to the detected damages, they were not considered for taking the SHM study. After analyzing S1 and S4, the S1 RRMS maps detected the damage better than its symmetrical counterpart S4. A frequency range of 200 kHz was considered for this study, as mentioned earlier. The full wavefield data of the cases showed that the A0 GW mode identified the damage cases clearly, which was chosen to study further. The A0 mode velocity values obtained for the LDV by scanning points at 90° and 0° were 2231 m/s, 1768 m/s, 2135 m/s, and 1677 m/s for the SEM calculations. The velocity profile was obtained by fitting the velocity values to an elliptical function and determining the 360° velocity profile.
The damage localization plots are shown in Figure 10b–f. The maps were plotted for all five cases (4 SEM and 1 LDV), and higher-energy regions were isolated with threshold values. The region shows that higher energy was nearer to the created BVIDs (circled in black). The error estimation is shown in Table 3 and Table 4. The values obtained were closer to the BVID regions and are in an error marginal range of >1 cm difference maximum.

4.3. Thermal Speckle-Based Verification of Damages

The collected data were carefully analyzed, considering the displacements in three directions and strain fields. In Figure 11a–c, the horizontal displacement of points in the first, 120th, and 240th steps (after two minutes of heating) were shown, respectively. In the figures, the horizontal displacement field was quasi even in each step and corresponds to the boundary conditions (clamping on the right side). In such a short period of heating, there were no visible differences between point displacements in damaged and undamaged regions.
However, the estimated standard deviation of the displacement in the horizontal direction for every data point (shown in Figure 11d) allows to differentiate the failure areas (damages); the standard deviation for damages was greater than for healthy regions in the neighbourhood. Similar observations were made for the other directions.

5. Conclusions

In this study, a combined SHM-GW and NDT-DIC method was proposed in studying BVIDs presented in the GFRPS. A frequency range of 200 kHz was analyzed for the GW-based studies. SEM numerical models were compared with the experimental data.
SEM helped to model BVID separately and together (experimental data) to validate the SHM methodology.
The RRMS study provided results about the damage locations even near the higher-energy zones.
S1 provided good analyses of the results, and scanning points were taken to analyze the damage zones.
The threshold-based proposed ETM algorithm predicted the location of the damages.
The error range was less than 1%, as shown in the difference (cm) of Table 4 in all the analyzed cases.
DIC-based diagrams of variation in displacements differentiated the damage region and thus verified it.
Combined online monitoring and verification of the structure was proposed.
Therefore, it is promising to use the proposed SHM-GW and DIC-NDT application in monitoring mechanical, civil structures to verify structural integrity. Future research studies involve implementing similar threshold-based localization and verification approaches in analyzing different structures for certification and reliability purposes.

Author Contributions

Conceptualization, K.B. and D.Z.; methodology, K.B., D.Z., P.F., P.M. and M.J.; software, K.B., P.M., D.Z. and P.F.; validation, K.B., D.Z. and P.F.; investigation and formal analysis, K.B., D.Z. and P.F.; data curation, K.B. and D.Z.; writing original draft preparation, K.B., D.Z. and P.F.; writing review and editing, K.B., D.Z., P.F., M.J. and P.M.; visualization, K.B., D.Z. and P.F.; supervision, P.M. and M.J.; project administration, P.M. and M.J.; funding acquisition P.M., D.Z., M.J. and P.F. All authors read and agreed to the published version of the manuscript.

Funding

The research work was carried out within project Nos. 2016/22/E/ST8/00068, 2019/03/X/ST8/01596, 2018/31/N/ST8/02865 and 2019/35/N/ST8/01086 financed by the National Science Centre, Poland.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the first, corresponding authors.

Acknowledgments

The authors also acknowledge Task-CI for the use of computational resources. This publication is based upon work from COST Action CA18123 (ODIN), supported by COST (European Cooperation in Science and Technology-Brussels, Belgium).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

The elastic constants of the unidirectional fibre composite in the form of the Hahn functional are gathered in Table A1.
Table A1. Elastic constants.
Table A1. Elastic constants.
Elastic ConstantPPfPmη
E11E11E11fE11m1
ν12ν12ν12fνm1
G121/G121/G12f1/Gmη6
G231/G231/G23f1/Gmη4
KT1/KT1/Kf1/KmηK
The Hahn function is given in the form:
P = ( P f V f + η P m V m ) ( V f + η V m )
where V f , V m are the volume fraction of the fibers and matrix, respectively. The bulk modulus (KT) for fibers and the matrix is equal to K f = E f / 2 ( 1 υ f ) , K m = E m / 2 ( 1 υ m ) , respectively, and the expressions for η are given as follow:
η 6 = 1 + G m / G 12 f 2 η 4 = 3 4 υ m + G m / G 23 f 4 ( 1 υ m ) , η K = 1 + G m / K f 2 ( 1 υ m ) .
The transverse moduli and υ 23 of the composite are given:
E 22 = E 33 = 4 K T G 23 K T + m G 23 , υ 23 = υ 12 f V f + υ m ( 1 V f ) [ 1 + υ m υ 12 E m / E 11 1 υ m 2 υ m υ 12 E m / E 11 ]
where m = 1 + 4 K T υ 12 2 E 11 .
According to the rules, the mixture mass density of the composite is equal ρ = ρ f V f + ρ m V m .

References

  1. Giurgiutiu, V. Structural Health Monitoring; Academic Press: London, UK, 2008. [Google Scholar]
  2. Li, J.; Lu, Y.; Guan, R.; Qu, W. Guided waves for debonding identification in CFRP-reinforced concrete beams. Constr. Build. Mater. 2017, 131, 388–399. [Google Scholar] [CrossRef]
  3. Stickel, J.; Nagarajan, M. Glass fiber-reinforced composites: From formulation to application. Int. J. Appl. Glass Sci. 2012, 3, 122–136. [Google Scholar] [CrossRef]
  4. Balasubramaniam, K.; Fiborek, P.; Malinowski, P. Lamb waves based assessment of impact damage in multilayered CFRP plate. In Proceedings of the Health Monitoring of Structural and Biological Systems XV, Online, 22–26 March 2021; SPIE Library: Belingham, WA, USA, 2021; p. 11593. [Google Scholar]
  5. Bogenfeld, R.; Kreikemeier, J.; Wille, T. An analytical scaling approach for low-velocity impact on composite structures. Compos. Struct. 2018, 187, 71–84. [Google Scholar] [CrossRef]
  6. Reis, P.; Ferreira, J.; Richardson, M. Fatigue damage characterization by NDT in polypropylene/glass fibre composites. Appl. Compos. Mater. 2010, 18, 409–419. [Google Scholar] [CrossRef]
  7. Miqoi, N.; Pomarede, P.; Meraghni, F.; Declercq, N.; Guillaumat, L.; Le Coz, G.; Delalande, S. Detection and evaluation of barely visible impact damage in woven glass fabric reinforced polyamide 6.6/6 composite using ultrasonic imaging, X-ray tomography and optical profilometry. Int. J. Damage Mech. 2020, 30, 323–348. [Google Scholar] [CrossRef]
  8. Sarasini, F.; Tirillò, J.; D’Altilia, S.; Valente, T.; Santulli, C.; Touchard, F.; Chocinski-Arnault, L.; Mellier, D.; Lampani, L.; Gaudenzi, P. Damage tolerance of carbon/flax hybrid composites subjected to low velocity impact. Compos. Part B Eng. 2016, 91, 144–153. [Google Scholar] [CrossRef]
  9. Talreja, R.; Phan, N. Assessment of damage tolerance approaches for composite aircraft with focus on barely visible impact damage. Compos. Struct. 2019, 219, 1–7. [Google Scholar] [CrossRef]
  10. Zima, B.; Kędra, R. Numerical study of concrete mesostructure effect on lamb wave propagation. Materials 2020, 13, 2570. [Google Scholar] [CrossRef]
  11. Tiwari, K.; Raisutis, R. Identification and characterization of defects in glass fiber reinforced plastic by refining the guided lamb waves. Materials 2018, 11, 1173. [Google Scholar] [CrossRef] [Green Version]
  12. Yan, G.; Lu, X.; Tang, J. Guided wave-based monitoring of evolution of fatigue damage in glass fiber/epoxy composites. Appl. Sci. 2019, 9, 1394. [Google Scholar] [CrossRef] [Green Version]
  13. Alnuaimi, H.; Amjad, U.; Russo, P.; Lopresto, V.; Kundu, T. Monitoring damage in composite plates from crack initiation to macro-crack propagation combining linear and nonlinear ultrasonic techniques. Struct. Health Monit. 2020, 20, 139–150. [Google Scholar] [CrossRef]
  14. Santos, M.; Santos, J.; Amaro, A.; Neto, M. Low velocity impact damage evaluation in fiber glass composite plates using PZT sensors. Compos. Part B Eng. 2013, 55, 269–276. [Google Scholar] [CrossRef]
  15. Ren, B.; Lissenden, C. Ultrasonic guided wave inspection of adhesive bonds between composite laminates. Int. J. Adhes. Adhes. 2013, 45, 59–68. [Google Scholar] [CrossRef]
  16. Zhang, C.; Zhang, Z.; Ji, H.; Qiu, J.; Tao, C. Mode conversion behavior of guided wave in glass fiber reinforced polymer with fatigue damage accumulation. Compos. Sci. Technol. 2020, 192, 108073. [Google Scholar] [CrossRef]
  17. Lee, C.; Zhang, A.; Yu, B.; Park, S. Comparison study between RMS and Edge detection image processing algorithms for a pulsed laser UWPI (ultrasonic wave propagation imaging)-based NDT technique. Sensors 2017, 17, 1224. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Kudela, P.; Wandowski, T.; Malinowski, P.; Ostachowicz, W. Application of scanning laser Doppler vibrometry for delamination detection in composite structures. Opt. Lasers Eng. 2017, 99, 46–57. [Google Scholar] [CrossRef]
  19. Kudela, P.; Żak, A.; Krawczuk, M.; Ostachowicz, W. Modelling of wave propagation in composite plates using the time domain spectral element method. J. Sound Vib. 2007, 302, 728–745. [Google Scholar] [CrossRef]
  20. Balasubramaniam, K.; Sikdar, S.; Fiborek, P.; Malinowski, P. Ultrasonic guided wave signal based nondestructive testing of a bonded composite structure using piezoelectric transducers. Signals 2021, 2, 13–24. [Google Scholar] [CrossRef]
  21. Dziendzikowski, M.; Kurnyta, A.; Dragan, K.; Klysz, S.; Leski, A. In situ barely visible impact damage detection and localization for composite structures using surface mounted and embedded PZT transducers: A comparative study. Mech. Syst. Signal Process. 2016, 78, 91–106. [Google Scholar] [CrossRef]
  22. Foti, F.; Pannier, Y.; Mellier, D.; Gigliotti, M. Damage characterization during high temperature fatigue of off-axis woven organic matrix composites for aircraft applications. IOP Conf. Ser. Mater. Sci. Eng. 2018, 406, 012055. [Google Scholar] [CrossRef]
  23. Ziaja, D.; Turoń, B.; Miller, B. Detection of anomaly in a pretensioned bolted beam-to-column connection node using digital image correlation and neural networks. Appl. Sci. 2020, 10, 2400. [Google Scholar] [CrossRef] [Green Version]
  24. Harrell, T.; Thomsen, O.; Dulieu-Barton, J.; Madsen, S. Damage in CFRP composites subjected to simulated lighting strikes-Assessment of thermal and mechanical responses. Compos. Part B Eng. 2019, 176, 107298. [Google Scholar] [CrossRef] [Green Version]
  25. Ha, N.; Le, V.; Goo, N. Investigation of fracture properties of a piezoelectric stack actuator using the digital image correlation technique. Int. J. Fatigue 2017, 101, 106–111. [Google Scholar] [CrossRef]
  26. Ha, N.; Le, V.; Goo, N.; Kim, J. Thermal strain measurement of austin stainless steel (SS304) during a heating-cooling process. Int. J. Aeronaut. Space Sci. 2017, 18, 206–214. [Google Scholar] [CrossRef]
  27. Jin, T.; Ha, N.; Le, V.; Goo, N.; Jeon, H. Thermal buckling measurement of a laminated composite plate under a uniform temperature distribution using the digital image correlation method. Compos. Struct. 2015, 123, 420–429. [Google Scholar] [CrossRef]
  28. Zhu, W.; Zhang, C.; Yang, L.; Zhou, Y.; Liu, Z. Real-time detection of damage evolution and fracture of EB-PVD thermal barrier coatings under thermal shock: An acoustic emission combined with digital image correlation method. Surf. Coat. Technol. 2020, 399, 126151. [Google Scholar] [CrossRef]
  29. Rajaram, S.; Cuadra, J.; Saralaya, R.; Bartoli, I.; Kontsos, A. In situ CTE measurements and damage detection using optical metrology. Measur. Sci. Technol. 2015, 27, 025202. [Google Scholar] [CrossRef]
  30. G-ANGELONI. Available online: https://www.g-angeloni.com/ (accessed on 10 November 2021).
  31. Ceramtec. Available online: https://www.ceramtec.com/files/ms_piezoceramic-soft-materials_en_de.pdf337 (accessed on 14 February 2020).
  32. Kudela, P.; Moll, J.; Fiborek, P. Parallel spectral element method for guided wave based structural health monitoring. Smart Mater. Struct. 2020, 29, 095010. [Google Scholar] [CrossRef]
  33. Auricchio, F.; Sacco, E. Refined first-order shear deformation theory models for composite laminates. J. Appl. Mech. 2003, 70, 381–390. [Google Scholar] [CrossRef]
  34. Vinson, J.R.; Sierakowski, R.L. The Behavior of Structures Composed of Composite Materials; Springer: Cham, Switzerland, 2006; Volume 105. [Google Scholar]
  35. Balasubramaniam, K.; Fiborek, P.; Ziaja, D.; Jurek, M.; Sawczak, M.; Soman, R.; Malinowski, P. Global and local area inspection methods in damage detection of carbon fiber composite structures. Measurement 2021, 33, 110336. [Google Scholar] [CrossRef]
  36. Fendzi, C.; Mechbal, N.; Rébillat, M.; Guskov, M.; Coffignal, G. A general Bayesian framework for ellipse-based and hyperbola-based damage localization in anisotropic composite plates. J. Intell. Mater. Syst. Struct. 2015, 27, 350–374. [Google Scholar] [CrossRef] [Green Version]
  37. Michaels, J.; Michaels, T. Guided wave signal processing and image fusion for in situ damage localization in plates. Wave Motion 2007, 44, 482–492. [Google Scholar] [CrossRef]
  38. Memmolo, V.; Maio, L.; Boffa, N.; Monaco, E.; Ricci, F. Damage detection tomography based on guided waves in composite structures using a distributed sensor network. Opt. Eng. 2015, 55, 011007. [Google Scholar] [CrossRef]
Figure 1. The methodological flow chart in damage analysis.
Figure 1. The methodological flow chart in damage analysis.
Materials 14 07268 g001
Figure 2. Schematic LDV experimental setup used with PZT and BVID coordinates.
Figure 2. Schematic LDV experimental setup used with PZT and BVID coordinates.
Materials 14 07268 g002
Figure 3. SEM_GFRP model cases studied.
Figure 3. SEM_GFRP model cases studied.
Materials 14 07268 g003
Figure 4. (a) DIC measuring stand; (b) scheme of the experimental setup.
Figure 4. (a) DIC measuring stand; (b) scheme of the experimental setup.
Materials 14 07268 g004
Figure 5. Damage analysis based on ETM.
Figure 5. Damage analysis based on ETM.
Materials 14 07268 g005
Figure 6. The exemplary view, which was registered by camera no. 1 (left) and camera no. 2 (right). The defects were invisible.
Figure 6. The exemplary view, which was registered by camera no. 1 (left) and camera no. 2 (right). The defects were invisible.
Materials 14 07268 g006
Figure 7. Cases 4 and 5: LDV_RRMS and SEM_RRMS plots for S1 excitation.
Figure 7. Cases 4 and 5: LDV_RRMS and SEM_RRMS plots for S1 excitation.
Materials 14 07268 g007
Figure 8. Cases 1, 2, and 3: SEM_RRMS plots for 200 kHz excitation.
Figure 8. Cases 1, 2, and 3: SEM_RRMS plots for 200 kHz excitation.
Materials 14 07268 g008
Figure 9. LDV_RRMS and SEM_RRMS plots for 200 kHz excitation of all S1–S4.
Figure 9. LDV_RRMS and SEM_RRMS plots for 200 kHz excitation of all S1–S4.
Materials 14 07268 g009aMaterials 14 07268 g009b
Figure 10. (a) Distribution of scanning nodal points taken; ETM for (b) case 1, (c) case 2, (d) case 3, (e) case 4, and (f) case 5.
Figure 10. (a) Distribution of scanning nodal points taken; ETM for (b) case 1, (c) case 2, (d) case 3, (e) case 4, and (f) case 5.
Materials 14 07268 g010aMaterials 14 07268 g010b
Figure 11. Horizontal displacement of points obtained by DIC measurements: (a) the first step (0.5 s of heating), (b) the 120 steps (after one minute of heating), (c) the 240 steps (after two minutes of heating), and (d) the estimated standard deviation of the displacement in the horizontal direction (mm). In (ac) the color was kept: from 0.35 mm (violet) to 0 mm (red).
Figure 11. Horizontal displacement of points obtained by DIC measurements: (a) the first step (0.5 s of heating), (b) the 120 steps (after one minute of heating), (c) the 240 steps (after two minutes of heating), and (d) the estimated standard deviation of the displacement in the horizontal direction (mm). In (ac) the color was kept: from 0.35 mm (violet) to 0 mm (red).
Materials 14 07268 g011
Table 1. Engineering constants for GFRP.
Table 1. Engineering constants for GFRP.
MaterialE11E22, E33G12G23ν12ν23ρV
GPaGPaGPaGPa--kg/m3%
Glass fiber75.075.030.7030.700.220.22260050
Epoxy3.403.401.301.300.350.35125050
GFRP39.207.753.412.820.290.371925-
Where E represents the Young’s modulus, G represents the shear modulus, ν represents the Poisson’s ratio, ρ is the density, and V is the volume fraction.
Table 2. RRMS cases studied.
Table 2. RRMS cases studied.
Case12345
TypeSEMSEMSEMSEMLDV
ModelBVID 1BVID 2BVID 3BVID 1,2,3BVID 1,2,3
Table 3. Error estimation for cases 1, 2, and 3 (center coordinates).
Table 3. Error estimation for cases 1, 2, and 3 (center coordinates).
SEM CasesCase 1Case 2Case 3
Actual location (cm)(38,14)(38,17.5)(38,24)
Estimated location (cm)(38.3,13.4)(38.5,17.4)(37.4,24.1)
Difference (cm)0.670.500.60
Table 4. Error estimation for cases 4 and 5 (center coordinates).
Table 4. Error estimation for cases 4 and 5 (center coordinates).
SEM and LDV CasesCase 4Case 5
BVID1BVID2BVID3BVID1BVID2BVID3
Actual location (cm)(38,14)(38,17.5)(38,24)(38,14)(38,17.5)(38,24)
Estimated location (cm)(38.7,13)(38.6,17)(37.7,24)(38.3,13.9)(38.3,17.3)(37.2,24)
Difference (cm)0.740.600.760.320.340.89
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Balasubramaniam, K.; Ziaja, D.; Jurek, M.; Fiborek, P.; Malinowski, P. Experimental and Numerical Analysis of Multiple Low-Velocity Impact Damages in a Glass Fibered Composite Structure. Materials 2021, 14, 7268. https://doi.org/10.3390/ma14237268

AMA Style

Balasubramaniam K, Ziaja D, Jurek M, Fiborek P, Malinowski P. Experimental and Numerical Analysis of Multiple Low-Velocity Impact Damages in a Glass Fibered Composite Structure. Materials. 2021; 14(23):7268. https://doi.org/10.3390/ma14237268

Chicago/Turabian Style

Balasubramaniam, Kaleeswaran, Dominika Ziaja, Michał Jurek, Piotr Fiborek, and Paweł Malinowski. 2021. "Experimental and Numerical Analysis of Multiple Low-Velocity Impact Damages in a Glass Fibered Composite Structure" Materials 14, no. 23: 7268. https://doi.org/10.3390/ma14237268

APA Style

Balasubramaniam, K., Ziaja, D., Jurek, M., Fiborek, P., & Malinowski, P. (2021). Experimental and Numerical Analysis of Multiple Low-Velocity Impact Damages in a Glass Fibered Composite Structure. Materials, 14(23), 7268. https://doi.org/10.3390/ma14237268

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop