Thermoelastic Investigation of Carbon-Fiber-Reinforced Composites Using a Drop-Weight Impact Test

: Composite materials are becoming more popular in technological applications due to the significant weight savings and strength offered by these materials compared to metallic materials. In many of these practical situations, the structures suffer from drop-impact loads. Materials and structures significantly change their behavior when submitted to impact loading conditions compared to quasi-static loading. The present work is devoted to investigating the thermal process in carbon-fiber-reinforced polymers (CFRP) subjected to a drop test. A novel drop-weight impact test experiment is performed to evaluate parameters specific to 3D composite materials. A strain gauge rosette and infrared thermography are employed to record the kinematic and thermal fields on the composites’ surfaces. This technique is nondestructive and offers an extensive full-field investigation of a material’s response. The combination of strain and infrared thermography data allows a comprehensive analysis of thermoelastic effects in CFRP when subjected to impacts. The experimental results are validated using numerical analysis by developing a MATLAB ® code to analyze whether the coupled heat and wave equation phenomenon exists in a two-dimensional polar coordinate system by discretizing through a forward-time central-space (FTCS) ﬁnite-difference method (FDM). The results show the coupling has no signiﬁcant impact as the waves generated due to impact disappears in 0.015 s. In contrast, heat diffusion happens for over a one-second period. This study demonstrates that the heat equation alone governs the CFRP heat ﬂow process, and the thermoelastic effect is negligible for the speciﬁc drop-weight impact load.


Introduction
Materials are the fundamental elements of all natural and human-made structures. New materials usually emerge due to the need to improve structural efficiency and performance. One of the best examples of this interrelated process is a composite material. Composites consist of two or more different materials in various combinations that form a structural unit. Composites are generally used because they have required properties that cannot be achieved through the use of the constituent materials alone. The most common example is a fibrous composite of reinforcing fibers embedded in a binder or matrix material. In general, fibers are the principal load-bearing members, while the surrounding matrix keeps them in their allocated position, and their orientation acts as a load-transfer medium between them and protects them from environmental damage due to chemical corrosion and humidity. Carbon-fiber-reinforced polymers (CFRP) offer a combination of strength and modulus comparable to or better than many traditional metallic materials. Despite the many superior properties of composites over metals, composites are susceptible to damage caused by low-velocity impact that reduce their performance to a large extent [1].
Composites have become the leading structural material in many industries due to their superior properties, including their high strength-to-weight ratio [2]. However, some aspects of switching the traditional structural material from metal to composite have not been well-investigated. A significant amount of research has been conducted to study the mechanical behavior of composites prior to their elastic threshold in relation to their thermal response. However, researchers have made few attempts to analyze this phenomenon with the help of IR thermography.

Non-Destructive Testing (NDT) of Composites-Thermographic Evaluation Method
IR thermography offers ideal, cost-effective NDT solutions for a wide range of composites applications, including detection of impact damage, delamination, water entrapment, inclusions, density variations, and evaluation of the adhesive bonds [3][4][5]. This technique, which is fast and contact-free, offers wide-area coverage of flat or curved parts. The principle of infrared thermographic NDT systems is based on the detection of infrared energy emitted from an object, and detection devices range from a simple heat gun/IR camera combination to advanced systems capable of defect measurement and advanced materials characterization [6,7].
The IR thermographic method measures temperature changes on the surface of materials resulting from stress-generated or externally applied thermal fields. Thermography can be used to study thermal variations. Excitation sources can be applied externally using thermal or mechanical sources to generate coupled stress and thermal fields. These excitation sources can be in the form of mechanical loads, direct heat, and ultrasonic stress waves. [8][9][10][11].

Drop-Weight Impact Testing (DWIT)
A drop-weight impact test was designed to examine a material's behavior under the action of a falling impact weight ( Figure 1). It is used to assess the impact properties of composites [12]. The impact weight falls from a predefined height under the action of gravity and hits the material being studied. The specimen used for the test is often a flat plate cut from larger material. The test piece is supported at its edges and impacted centrally by a vertically falling weight. The impact performance of composites is concerned with absorbing energy in the system when the component is struck, either through deformation or damage development. A strain-time curve is obtained from the impact event by attaching a strain gauge rosette. The test results enable the investigation of the behavior of specimens during and after the impact test.

Thermoelastic and Thermoplastic Effects
The thermoelastic effect is a small reversible change in the temperature of an object when it undergoes elastic deformation. The thermoelastic effect depends upon the type of strain, which can be elastic or non-elastic. As described by Lord Kelvin in 1853 [13], a

Thermoelastic and Thermoplastic Effects
The thermoelastic effect is a small reversible change in the temperature of an object when it undergoes elastic deformation. The thermoelastic effect depends upon the type of strain, which can be elastic or non-elastic. As described by Lord Kelvin in 1853 [13], a purely elastic tensile strain in an adiabatic configuration generates a thermoelastic effect that results in the cooling of isotropic specimens with positive thermal expansion coefficients. However, if inelastic deformation is also present, then the thermoplastic phenomenon occurs, resulting in heating due to the dissipation of the applied mechanical energy; this will normally overcome the thermoelastic contribution and result in a net specimen temperature increase. Mechanical deformation causes temperature changes, but they go unnoticed if the loading is slow, resulting in an isothermal process [14]. At the other extreme, under high-speed loading conditions, considerable viscous heating can occur, possibly leading to localized melting [15].
The thermoelastic response is produced by dilational strain. Tensile stress produces a drop in temperature, whereas compressive stress results in a temperature increase. Nowacki [16] developed the theory of thermoelastic diffusion. In this theory, the coupled thermoelastic model that is used implies infinite speeds for the propagation of thermoelastic waves. The overall effects arising from the coupling of fields of temperature, mass diffusion, and strain in an elastic cylinder have been discussed by Olesiak and Pyryev [17]. Sherief et al. [18][19][20][21] introduced the theory of generalized thermoelastic diffusion, which predicts finite speeds of propagation for thermoelastic and diffusive waves. The present study was motivated by the importance of studying the thermoelastic effect in CFRP.
When composites are subjected to dynamic loading under impact events, the zero or near-zero plasticity of these structures past their yielding points results in a susceptibility to delamination, matrix cracking, and fiber breakage [22].

Numerical Analysis
Different modeling techniques are used to model the impact dynamics of composite structures. In general, there are four categories considered: energy-balance models, which assume a quasi-static condition; spring-mass models, which simplify the dynamics of the structure; complete models, in which the full behavior of the structure is modeled; and infinite/finite plate models. The infinite plate model can be considered if the deformation wavefront does not reach the boundary of the plate. If the wave propagation reflects on the plate boundaries back to the point of impact, then the finite size of the test sample must be considered [23].

Methodology
In this study, the thermoelastic response of carbon-fiber-reinforced polymers (CFRP) subjected to high strain rates was investigated using a novel DWIT experimental setup. Fiber-reinforced composites are usually manufactured by stacking laminae of fiber sheets pre-impregnated with resin. The structure is therefore heterogeneous, and so the stress discontinuities between the fibers and matrix should be considered for the purpose of calculating the overall temperature field [24,25]. It would be expected that the fibers and matrix would experience different temperature fluctuations due to their dissimilar thermoelastic constants [14]. After the drop impact, viscous (frictional) energy is converted into thermal energy. This energy propagates from the point of impact toward the boundaries due to the temperature gradient. The goal of this study was to investigate the thermoelastic phenomenon in the CFRP material that produces spatial variations in the surface temperatures and allows correlation between measured time-averaged IR images, heat flow due to pure diffusion, and the strain waves on the surface.
The study was divided into two phases. The first phase was the experimental study, which involved performing the DWIT with a specialized rosette strain gauge, amplifiers, a high-speed oscilloscope, and an IR thermographic camera. This study yielded useful insight into the oscillatory elastic deformation (Video S1: Drop-Weight Impact Testing (DWIT) video) and internal friction phenomenon captured via IR thermography within the CFRP (Video S2: IR thermographic recording) when subjected to drop-weight impacts.
The second phase was the development of numerical models to describe the physical problem that was predicted from the experimental findings. The heat and wave equations in polar coordinate systems were coupled and solved using the forward-time central-space (FTCS) method and finite-difference method (FDM) using MATLAB ® software. The derived numerical analysis solutions were validated using experimental findings to optimize the models.

CFRP Sample
The CFRP test samples used in this study were sheets of DragonPlate ® (Allred and Associates Inc., Elbridge, NY, USA) [26]. Each sample was a solid quasi-isotropic carbon-fiber sheet 1 mm thick ( Figure 2) [27]. The CFRP sheets comprised five plies at Figure 2). The samples were composed of a tough and rigid carbon-reinforced epoxy matrix, with a textured finish on both sides. The properties of the samples are given in Table 1.

Experimental Setup to Perform a Drop Test
The drop-weight impact test was performed using a spherical steel ball 20 mm in diameter ( Figure 3a). The ball was wrapped with black tape to decrease steel reflectance for better thermal images( Figure 3b) [29,30]. A rosette strain gauge in a stacked configuration (at 0°, 45°, and 90° angles) was attached to the back of the CFRP plate ( Figure 4). The acquisition of the differential strain vector from the baseline was achieved using an op-amp circuit and a high-speed oscilloscope at a 100 MHz sampling frequency rate. The CFRP plate was clamped in a circular flange to prevent the plate from moving during the DWIT (Figure 5a

Experimental Setup to Perform a Drop Test
The drop-weight impact test was performed using a spherical steel ball 20 mm in diameter ( Figure 3a). The ball was wrapped with black tape to decrease steel reflectance for better thermal images( Figure 3b) [29,30]. A rosette strain gauge in a stacked configuration (at 0 • , 45 • , and 90 • angles) was attached to the back of the CFRP plate ( Figure 4). The acquisition of the differential strain vector from the baseline was achieved using an op-amp circuit and a high-speed oscilloscope at a 100 MHz sampling frequency rate. The CFRP plate was clamped in a circular flange to prevent the plate from moving during the DWIT (Figure 5a,b).
The test specimens were monitored using high-speed IR thermography throughout the course of the experiment using an FLIR ® T1030sc thermal camera [31], and were analyzed using Researcher IR Max software [32] to investigate the heat generation due to friction and diffusion associated with the thermoelastic effects. The physical process to collect IR thermographic images when excitation occurs through a mechanical impact can be explained in three steps (diagrammatic representation of the experiment is shown in Figure 6): for better thermal images( Figure 3b) [29,30]. A rosette strain gauge in a stacked configuration (at 0°, 45°, and 90° angles) was attached to the back of the CFRP plate ( Figure 4). The acquisition of the differential strain vector from the baseline was achieved using an op-amp circuit and a high-speed oscilloscope at a 100 MHz sampling frequency rate. The CFRP plate was clamped in a circular flange to prevent the plate from moving during the DWIT (Figure 5a,b).   The test specimens were monitored using high-speed IR thermography throughout the course of the experiment using an FLIR ® T1030sc thermal camera [31], and were analyzed using Researcher IR Max software [32] to investigate the heat generation due to friction and diffusion associated with the thermoelastic effects. The physical process to  The test specimens were monitored using high-speed IR thermography throughout the course of the experiment using an FLIR ® T1030sc thermal camera [31], and were analyzed using Researcher IR Max software [32] to investigate the heat generation due to friction and diffusion associated with the thermoelastic effects. The physical process to collect IR thermographic images when excitation occurs through a mechanical impact can be explained in three steps (diagrammatic representation of the experiment is shown in  (3) An image of the CFRP plate's surface temperature was monitored using an IR camera, which collected time-averaged IR radiations emitted from the surface of the sample. The frame rate of the IR camera was approximately 125 Hz.

Finite Difference Method (MATLAB ® )
The FDM is a numerical method for solving differential equations such as the twodimensional wave [34,35] and heat equation [36][37][38], as given in Equations (1) and (4). This method approximates the differentials by discretizing the dependent variables (strain and temperature) in the independent variable domains (space and time) [39][40][41]. Each discretized value of the dependent variable is referred to as a nodal value. In this case, 2D wave and heat equations in the cylindrical coordinate system are discretized using FTCS FDM. The discretized equations are given in Equations (3) and (5).
The two-dimensional wave equation in the cylindrical coordinate system is shown in Equation (1): where is the strain, is the radius, is the angle in degrees, and is the speed of a sound in a medium as shown in Equation (2) [42].
The two-dimensional discretized wave equation in the cylindrical coordinate system is given in Equation (3).
The two-dimensional heat conduction equation in a cylindrical coordinate system with no energy-generation term is given in Equation (4).
The two-dimensional discretized heat equation in a cylindrical coordinate system is shown in Equation (5): (1) The surface of the component was initially at a uniform ambient temperature. The excitation was accomplished using a lightweight indenter; in this study, a 33 g weight steel ball was used.
(2) The frictional energy was converted into thermal energy at the point of impact; heat diffused toward the cooler edges. The CFRP steel friction coefficient was 0.33 [33].
(3) An image of the CFRP plate's surface temperature was monitored using an IR camera, which collected time-averaged IR radiations emitted from the surface of the sample. The frame rate of the IR camera was approximately 125 Hz.

Finite Difference Method (MATLAB ® )
The FDM is a numerical method for solving differential equations such as the twodimensional wave [34,35] and heat equation [36][37][38], as given in Equations (1) and (4). This method approximates the differentials by discretizing the dependent variables (strain and temperature) in the independent variable domains (space and time) [39][40][41]. Each discretized value of the dependent variable is referred to as a nodal value. In this case, 2D wave and heat equations in the cylindrical coordinate system are discretized using FTCS FDM. The discretized equations are given in Equations (3) and (5).
The two-dimensional wave equation in the cylindrical coordinate system is shown in Equation (1): where u is the strain, r is the radius, θ is the angle in degrees, and c is the speed of a sound in a medium as shown in Equation (2) [42].
The two-dimensional discretized wave equation in the cylindrical coordinate system is given in Equation (3).
7 of 17 The two-dimensional heat conduction equation in a cylindrical coordinate system with no energy-generation term is given in Equation (4).
The two-dimensional discretized heat equation in a cylindrical coordinate system is shown in Equation (5): where T is the temperature; t, r, and θ refer to time and positions; ∆t is the timestep size in (s); and ∆r and ∆θ are the differences in the spatial positions of the temperature nodes. α = 0.12 mm 2 /s in-plane thermal conductivity of CFRPs [43,44].
Overall temperature changes in the CFRP plate can be obtained by solving the coupled heat and wave equation [45][46][47] as given in Equation (6): where ∆T conduction = T t+1 i,j T t i,j ∆t , which can be obtained from the heat equation in a cylindrical coordinate system as given in Equation (5). ∆T thermoelastic = ε CTE , ε is the thermoelastic strain obtained from wave equation in a cylindrical coordinate system as given in Equation (3), and CTE = −4.5·10 −5 • C −1 [33,42] is the coefficient of thermal expansion of the CFRP.

Mathematical Model
The mathematical model of the experiment is given by governing Equations (7) and (8) and the initial-boundary values as given in Equations (9) and (10): where T is the temperature and is the trace of the stress tensor. The parameters c, κ, α are the wave speed, thermal conductivity, and thermal expansion factor, respectively. To reach an analytical solution, we introduce τ = T − T 0 and linearize T∂ t ≈ T 0 ∂ t . Then we solve Equation (11) with the boundary conditions given in Equations (12) and (13): To solve the linearized version of the problem with axisymmetric initial and boundary data, we introduce the eigenvalues and eigenfunctions of the axisymmetric Laplace operator on the circular CFRP sheet with radius R, D R . The eigenfunctions for m = 1, · · · , ∞ are given in Equation (14): where J 0 is the 0th Bessel function of the first kind, and ρ m is the mth zero of J 0 . To simplify the notation, we introduceρ m = ρ m R , the eigenvalue corresponding to the eigenfunctions u m is -ρ 2 m . The solution is given in Equations (15) and (16) The coefficients are then computed as follows in Equation (17) through (19): where J 1 is the first Bessel function of the first kind.

Numerical Analysis
The two-dimensional discretized wave equation in cylindrical coordinates was initially solved using the FTCS FDM, a numerical method in MATLAB ® . The results were compared with the experiments, as discussed in Section 4.2.1. Afterward, the coupled heat and wave equation in the two-dimensional polar coordinate system was solved. As an initial condition, a constant temperature of 21.9 • C was specified throughout the domain, except at the source. It is vital for the stability and accuracy of the FDM to choose the correct timestep value. In this work, the Courant-Friedrichs-Lewy (CFL) condition [48] was used to decide the timestep size. The CFL condition is given in Equation (20): 2α∆t ≤ min (∆r) 2 , (∆θ) 2 (20) where α is the velocity magnitude, ∆t = dt is the timestep, and ∆r = dr and ∆θ = dr are the radius and angle between mesh elements. For solving the wave equation for a total of 3000 timesteps, the following were taken: dt = 5; % microseconds; dr = 1; % mm; dtheta = 1; % degrees.

Results and Discussion
The experiment's strain data for a fully clamped CFRP plate (with 21,000 times amplification using an op-amp circuit and a high-speed oscilloscope) at a sampling frequency of 100 MHz after signal processing showed a strain-wave frequency of 205 Hz. At the time of impact, the highest value of strain recorded was 2.5 microstrain, but was due to the ball's momentum and was not representative of the natural response of the CFRP plate. After the ball was removed, the plate vibrated naturally, and a value of 0.5 microstrain was recorded (Figure 7). This value was used to compute the thermoelastic effect while solving the coupled wave and heat equation by the FDM method in MATLAB ® . The strain waves experienced significant damping, and the amplitude dropped to 0.15 in 0.012 s. The dynamic strain reflected the behavior of the CFRP under drop impact for in-plane measurement.

Thermographic Experimental Results
When the steel ball hit the CFRP plate in DWIT, the kinetic energy of the ball was converted into frictional energy that caused localized heating of the specimen at the point of impact. The heat generated due to friction was recorded with an IR thermal camera. Upon impact, some of the energy was also consumed in generating sinusoidal vibrations in the plate, which was being recorded with the help of a rosette strain gauge attached at the center on the backside of the plate. The strain gauge data allowed the interpretation of the thermoelastic effect in the CFRP when it was subjected to impacts.
During the impact event, the infrared camera viewed the CFRP plate from the top, and recorded the thermal signatures that resulted from both the thermoelastic and heatdiffusion phenomenon. The recorded sequences of images were post-processed using FLIR Research IR software. As a primary step, a sequence of ΔT images was created by subtracting the temperature from the ambient temperature recorded before the impact for every image of the sequence recorded during the impact event. More specifically, ΔT was obtained from the relationship given in Equation (21): where ( , , 0) = 21.9 ℃, and t is the time instant at which each image is recorded starting from impact. Some ∆T images, at different time instants, are shown in Figure 8. The CFRP plate was impacted from t = 0.216 to 0.258 s. Heat was visible (bright region) at the point of impact due to the friction; a maximum temperature increase of 2.6 ℃ was observed.
As the ball moved away from the plate, the point of impact (hotspot) is visible on the CFRP plate. This point of impact slowly dissipates the thermal energy due to the thermal gradient between the center point and plate's edges. The hotspot can be seen disappearing in the frame at t = 0.649 s.

Thermographic Experimental Results
When the steel ball hit the CFRP plate in DWIT, the kinetic energy of the ball was converted into frictional energy that caused localized heating of the specimen at the point of impact. The heat generated due to friction was recorded with an IR thermal camera. Upon impact, some of the energy was also consumed in generating sinusoidal vibrations in the plate, which was being recorded with the help of a rosette strain gauge attached at the center on the backside of the plate. The strain gauge data allowed the interpretation of the thermoelastic effect in the CFRP when it was subjected to impacts.
During the impact event, the infrared camera viewed the CFRP plate from the top, and recorded the thermal signatures that resulted from both the thermoelastic and heatdiffusion phenomenon. The recorded sequences of images were post-processed using FLIR Research IR software. As a primary step, a sequence of ∆T images was created by subtracting the temperature from the ambient temperature recorded before the impact for every image of the sequence recorded during the impact event. More specifically, ∆T was obtained from the relationship given in Equation (21): where T(r, θ, 0) = 21.9 • C, and t is the time instant at which each image is recorded starting from impact. Some ∆T images, at different time instants, are shown in Figure 8. The CFRP plate was impacted from t = 0.216 to 0.258 s. Heat was visible (bright region) at the point of impact due to the friction; a maximum temperature increase of 2.6 • C was observed. As the ball moved away from the plate, the point of impact (hotspot) is visible on the CFRP plate. This point of impact slowly dissipates the thermal energy due to the thermal gradient between the center point and plate's edges. The hotspot can be seen disappearing in the frame at t = 0.649 s. Figure 9 shows the temperature versus time plot of the DWIT as captured through a thermographic imaging technique. The plot indicates a sharp rise in temperature as the ball strikes the CFRP plate; the temperature increased from 21.9 to 24.5 • C. After the impact, the plate began dissipating energy through diffusion and thermoelastic vibrations, and the plate returned to a constant temperature after t = 0.875 s.  Figure 9 shows the temperature versus time plot of the DWIT as captured through a thermographic imaging technique. The plot indicates a sharp rise in temperature as the ball strikes the CFRP plate; the temperature increased from 21.9 to 24.5 ℃. After the impact, the plate began dissipating energy through diffusion and thermoelastic vibrations, and the plate returned to a constant temperature after t = 0.875 s.

Wave Equation Results
The discretized wave equation was solved using MATLAB ® (time = 3000 timesteps, radius = 87 mm). The initial conditions given in the code were the quadratic profile and a

Wave Equation Results
The discretized wave equation was solved using MATLAB ® (time = 3000 timesteps, radius = 87 mm). The initial conditions given in the code were the quadratic profile and a maximum deformation of 0.5 V (recorded by high-frequency oscilloscope from experiments; this value is directly proportional to the strain). The damping term was also added, with a coefficient of damping = 0.006. After plotting the experimental and MATLAB ® simulation results together (Figure 10), the wave equation predicted a profile and frequency of 205 Hz, as obtained from the experiments (code is attached in the Appendix). Figure 9. The temperature vs. time plot of the DWIT captured using thermography.

Wave Equation Results
The discretized wave equation was solved using MATLAB ® (time = 3000 timestep radius = 87 mm). The initial conditions given in the code were the quadratic profile and maximum deformation of 0.5 V (recorded by high-frequency oscilloscope from expe ments; this value is directly proportional to the strain). The damping term was also adde with a coefficient of damping = 0.006. After plotting the experimental and MATLAB ® sim ulation results together (Figure 10), the wave equation predicted a profile and frequen of 205 Hz, as obtained from the experiments (code is attached in the Appendix).

Coupled Heat and Wave Equation Results
The discretized coupled-wave and heat equation was solved using the FTCS FDM numerical method in a two-dimensional cylindrical coordinate system. The resu showed the coupling had no significant impact, as the waves generated due to impa disappeared in 0.015 s. In contrast, heat diffusion occurred for nearly 1.2 s. This prov the heat equation alone is governing the heat-flow process in the CFRP, and the thermo lastic effect is negligible. The coupled wave and heat equations produced the same resu

Coupled Heat and Wave Equation Results
The discretized coupled-wave and heat equation was solved using the FTCS FDM, a numerical method in a two-dimensional cylindrical coordinate system. The results showed the coupling had no significant impact, as the waves generated due to impact disappeared in 0.015 s. In contrast, heat diffusion occurred for nearly 1.2 s. This proves the heat equation alone is governing the heat-flow process in the CFRP, and the thermoelastic effect is negligible. The coupled wave and heat equations produced the same results as the heat equation alone. The heat equation was solved for a total of 300 timesteps, and one timestep of 25 microseconds was chosen, and space steps of 1 mm and 1 • were defined in the code. The IR thermographic frames captured during the experiment at specific time instants were compared to numerical simulations shown as two-dimensional plots ( Figure 11). The experimental and simulation frames show heat dissipation from t = 0.258 to 0.892 s. For a better comparison, the experiment and simulation results were plotted on the same grid ( Figure 12). The experiment and simulation results were similar, as both showed a similar temperature drop, from 24.5 to 22 • C in 1.6 s. Afterward, the CFRP plate was in equilibrium with the surroundings, as the room temperature was 21.

Summary
The DWIT and nondestructive thermographic study results illustrate the convenience of using IR thermography to monitor impact tests. Useful information about a material's behavior can be derived quickly, at almost the same time as the test, without the

Summary
The DWIT and nondestructive thermographic study results illustrate the convenience of using IR thermography to monitor impact tests. Useful information about a material's behavior can be derived quickly, at almost the same time as the test, without the cumbersome need for post-impact processing.
The temperature versus time images and plots obtained in our study showed that heat is generated in the CFRP upon impact, and then dissipates over time. The coupled wave and heat equation was solved using a numerical method to understand the mechanism for heat generation. The results showed that the thermoelastic phenomenon is high-speed, dissipating in approximately 0.015 s, while heat diffusion occurred for nearly 1.2 s. The results for the heat-diffusion equation, which matched the experimental results captured using IR thermography, validated the findings. Our study found that the thermoelasticity (elastic waves) has a negligible effect, and that heat diffusion is responsible for heat dissipation in the CFRP. This indicates that the heat-diffusion equation entirely governed the heat-conduction mechanism of the CFRP upon impact.
Overall, temperature changes in the CFRP plate upon impact began with a rise in temperature due to friction, followed by the appearance of a hotspot and heat flow towards the edges. The flow of heat obeyed the heat-diffusion equation and mathematical model in a two-dimensional cylindrical coordinate system, as the results agreed with the experimental findings.

•
The thermoelastic phenomenon generated upon impact in a quasi-isotropic CFRP sheet is high-speed and completes in almost 0.015 s.

•
Heat diffusion from the hotspot, which developed due to friction between the steel ball and the CFRP surface, occurs for nearly 1.2 s. • No significant effect of the coupling heat and wave equation was observed, as temperature change due to thermoelasticity (elastic waves) was minimal compared to temperature change due to heat conduction.

•
We found that thermoelasticity has a negligible effect, and heat diffusion alone is responsible for the heat dissipation in the CFRP.