Next Article in Journal
Excitation of Wet Perovskite Films by Ultrasonic Vibration Improves the Device Performance
Next Article in Special Issue
Enhancing the Accuracy of Advanced High Temperature Mechanical Testing through Thermography
Previous Article in Journal
Evaluation of Low-Temperature Cracking Performance of Asphalt Pavements Using Acoustic Emission: A Review
Previous Article in Special Issue
Infrared Thermography for Inline Monitoring of Glass/Epoxy under Impact and Quasi-Static Bending
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

From Thermal Inspection to Updating a Numerical Model of a Race Bicycle: Comparison with Structural Dynamics Approach

Deptartment Electromechanics, University of Antwerp, Op3Mech research group, Groenenborgerlaan 171, B-2020 Antwerp, Belgium
*
Author to whom correspondence should be addressed.
Additional affiliation: Deptartment Electromechanics, Deptartment Mechanical Engineering, Acoustics & Vibration Research Group, Vrije Universiteit Brussel, Pleinlaan 2, B-1050 Brussels, Belgium.
Appl. Sci. 2018, 8(2), 307; https://doi.org/10.3390/app8020307
Submission received: 15 January 2018 / Revised: 9 February 2018 / Accepted: 14 February 2018 / Published: 21 February 2018

Abstract

:

Featured Application

The work described in this study provides the reader with a comparison between thermal model updating and well-known model updating in structural dynamics. This provides an idea of the applicability of thermal model updating and the advantages and disadvantages of well established techniques.

Abstract

Carbon fiber bicycle frames are complex-shaped structures and are prone to delaminations and difficult to inspect. The use of finite element model updating is common in structural dynamics but not so common in active thermography inspection. However, there are many advantages to using thermography when inspecting bicycle frames. These include the fact that the inspection can be performed in situ, can cover large areas, and is a quantitative method. In this paper, a numerical model of a bicycle frame will be updated and optimized by the surface temperature distribution captured with pulsed thermography. These results will be compared and benchmarked against frequency response function (FRF) measurement data as a reference. The chosen temperature decay measurements to be used as reference data will be of key importance. The goal of this manuscript is to compare both measurement results and model predictabilities after performing finite element model updating with respect to accuracy and speed.

1. Introduction

The application of composite materials in performance components used daily such as bicycle frames has increased due to its beneficial mechanical properties with respect to metals. Such materials are strong, lightweight, and fatigue- and corrosion-resistant when combined with the proper resins [1]. Composite materials are technically complex structures due to their anisotropic behavior, which makes quality control particularly challenging. Thus, structural health monitoring is essential in the production and lifetime stages of such composite structures. This work is concerned with carbon bicycle frames [2]. To detect and analyze defects and guarantee an extended lifetime, composite materials have been the subject of a tremendous amount of non-destructive inspection (NDI) research [3,4]. Non-destructive evaluation (NDE) techniques are broadly used to evaluate material properties and to inspect the integrity of structures without damaging the structure [5]. Infrared (IR) thermography has remarkable advantages with respect to conventional NDE techniques, such as ultrasound C-scanning, X-ray radiography, tap-testing, and visual testing [6]. The technology is contactless, safe, fast, and reliable [7]. The study of active thermography has become an important NDE technique for damage detection and material updating in metallic structural components [6,8,9] and in carbon fiber reinforced polymer (CFRP) and glass fiber reinforced polymer (GFRP) composites [10]. As the application of composite laminates such as CFRPs and GFRPs increases, the thermal conductance inside the components becomes anisotropic, and quantitative prediction becomes complex. The common form of active thermography makes use of heat pulses emitted to the surfaces and is called pulse thermography (PT) [11]. PT is used not only to identify defects but also to determine the thermal diffusivity of the material, inspect the structural integrity, and monitor the production process of composite structures.
High-performance carbon fiber bicycle frames are strong and lightweight. However, they are susceptible to damage caused by low-energy impact loading, which have many possible causes: small stones can be projected by nearby cars or the front wheel of the bicycle, the bicycle may fall and hit the ground, and defects can occur during manufacturing [12,13]. In addition, material thickness inconsistency can occur during manufacturing and cause false positive defect hits.
Finite element (FE) modeling is common in the design of large structures but generally remains in the design stage. Two of the largest drawbacks of the use of FE beyond the design stage are the approximated material property parameters and the high computational cost of building an accurate FE model of the full structure. Finite element model updating (FEMU) can deliver a solution to these two problems, whereby the model is updated with non-destructive testing data, mainly done with ultrasound vibration C-scanning [14,15,16]. The foundation of the technique has been described by Friswell and Mottershead in 1995 [14].
The main idea is to solve a set of initial numerical models with variable parameter values for the uncertain properties. By comparing the resulting response values (mostly done on vibration modes and/or eigenfrequencies) with experimental measurements, the variable parameters are modified iteratively until the response values correlate with the experimental verification data. Especially for complex-shaped objects and non-linear behavior, the methodology delivers faster and more accurate results with the inclusion of higher-order terms in the equations [14].
Unfortunately, the computational cost of an FEMU routine is even more expensive than that of a regular numerical simulation, as for each iteration at least one new numerical model needs to be computed. Therefore, a methodology based on the technology developed by Steenackers has been used to replace the FE model with an approximative regression model in order to reduce the calculation times [17]. When the relationship between the input parameters and the requested output information on the mesh nodes can be modeled with significant accuracy, a major part of the time-consuming FE computations can be bypassed during the model updating routine by implementing a so called “metamodel.” Some further computational intelligence techniques have been used within FEMU routines for structural dynamics and compared by Marwala in 2010 [15].
Some general challenges remain. Due to inaccuracy in the model and to measurement errors, there is a danger that wrong assumptions are made. Proper validation, physical understanding, and adequate modeling are essential to verify the value of the performed FEMU routines.
These updating methods are widely used in structural dynamics for modal vibration inspection [14,18,19,20,21,22], but their use is limited in heat transfer optimization problems related to non-linear convection and radiation boundaries [23,24]. FE model updating can be used to optimize the complex parameter set of active thermographic setups by correlating the numerical results with NDT data [25]. These methods are widely used in modal vibration inspection [20,21,22], but their usage is limited in thermal optimization problems [4,23,24] due to the difficulties in estimating the diffusivity.
In this paper, we compare ‘easy-to-perform’ FEMU on thermal non-destructive measurements, with respect to vibrational analysis, with regular finite element model updating. We therefore limit our experimental measurements to non-contact measurement techniques from a single view point. In practice, it is generally difficult to measure the same structure from different positions around the structure. The goal is to evaluate the influence of this limitation on the measurement results in order to evaluate local material thickness.

2. Test-Case Description

2.1. Thermal Experiments

A carbon fiber reinforced plastic (CFRP) bicycle frame was inspected in order to estimate the tube thickness and inspect for defects such as cracks and delaminations using an MWIR FLIR Phoenix thermal camera (3–5 μ m, 640 × 512 pixels, 25 mm lens) and an optical excitation source consisting of two 3.2 kJ Balcar flash excitation units. The bicycle was painted in black and yellow, and the tubular shapes and metallic inserts made evaluating the structure challenging. A representation of the test setup is shown in Figure 1.

2.2. Structural Dynamics Experiments

The same CFRP bicycle frame was evaluated by measuring frequency response functions (FRFs) and using state-of-the-art modal analysis techniques with a PSV 400 Polytec scanning laser doppler vibrometer system and an electrodynamic modal shaker. The bicycle frame was suspended freely, approximating free-free boundary conditions, and accelerations were measured on different positions over the structure, which is shown in Figure 2. From the measured FRFs, eigenfrequencies and corresponding mode shapes were estimated as reference data to be used for modal updating. In total, the setup contained 290 measurement points. Within all points, the deflection along 1 degree of freedom (DOF) was measured. With respect to the optimal excitation direction, we used the DOF perpendicular to the left side of the bicycle frame. FRF measurements were performed in free-free boundary conditions, using a swept sine excitation over a frequency band of 1–6000 Hz with a resolution of 0.5 Hz, and the modal parameters (modal frequencies, damping, and mode shapes) were estimated [26].

2.3. Numerical Simulation

The finite element model of the CFRP bicycle frame was built from a 3D scan measured with a Creaform 3D scanner with an accuracy of less than 0.1 mm using the structured light principle. The resulting .stl was cleaned and used to build a 2D structured mesh of CQUAD elements within the Siemens Simcenter 11 software (Siemens PLM, Plano, Texas, US), and an overall element size of 5 mm was used. The tubes were formed by an extruded 3D mesh from the 2D mesh basis. A representation of the 2D mesh is shown in Figure 3. The laminated properties of the bicycle frame were approximated using a five-layered laminate consisting of material properties shown in Table 1:
  • a thin uniform layer of paint 0.11 mm thick;
  • a CFRP layer 25% of the overall thickness with a global fiber orientation of 0 ;
  • a CFRP layer 25% of the overall thickness with a global fiber orientation of 90 ;
  • a CFRP layer 25% of the overall thickness with a global fiber orientation of 90 ;
  • a CFRP layer 25% of the overall thickness with a global fiber orientation of 0 .
A symmetric quasi-isotropic laminate was represented as the stacking recipient. The precise stacking and layer properties are unknown, so it was assumed that the represented stacking was sufficiently accurate as an input for the material parameter optimization process, where the global thermal diffusivity properties in the thermal case and the mode shapes in the structural case are updated based on reference measurements. A detailed picture of the 3D mesh with the different layers is shown in Figure 4. The mesh consisted of 28,623 elements with a structured material orientation according to the laminate properties. The metal inserts were neglected, as they were not positioned in the vicinity of the evaluated measurement areas.
In contrast to the mesh generation, numerical analysis was performed using Comsol Multiphysics 5.0 to facilitate interaction with the optimization routine using Matlab.
The excitation source was modeled as a combination of two point sources. This excitation source results in a non-homogeneous spatial distributed heat excitation. The thermal time constant results in, aside from a spatial non-homogeneity, a different distribution related to the type and history of the heat source, which needed to be considered in the numerical model. The used flash excitation units are visualized in Figure 5. It can be clearly seen that the excited radiation is delayed from the initiated trigger signal due to the thermal inertia. To correlate the numerical model accurately with the experimental data, this thermal inertia will be considered as a correction (weighting) of the input signal with the thermal response.

3. Methodology

Both the quality of the experimental data and the FE model data will contribute to the evolution of the optimization algorithm. The implemented optimization algorithm is based on the adaptive response surface method described in [24]. With respect to modal analysis updating, the original optimization routine is used, developed by Steenackers et al. [27]. In this section, we will start with a short description of the experimental data acquisition and post-processing of the thermal and modal data. The numerical model is briefly described in Section 2.3, which ends with a brief description of the used optimization routine for finite element model updating.

3.1. Thermal Experiments

The thermal experiments were performed using pulsed thermography (PT), schematically represented in Figure 6, via a flash excitation source and a thermal camera. This active-thermography NDE technique is based on an analysis of the response of the surface temperature of the structure to an emitted thermal wave initiated by an excitation pulse [28]. It is one of the most popular methods in active thermography due to its simplicity and fast operation.
The thermal data were stored for several time steps. A 3D matrix with the XY-direction containing the temperature of each pixel and the Z-direction capturing each time frame was delivered [4,29]. Local material properties, thickness, or experimental set-up parameters, such as excitation power and distance, could be defined [29]. The measurement accuracies were dependent on ambient temperature, emissivity, reflecting temperature, and humidity. These parameters were fixed when measurements were taken. An empirical rule of thumb was used to determine the minimally detectable diameter/depth ratio (should be equal to 2) [30].
Measurements were post-processed using pulsed phase thermography (PPT), which is a popular signal transformation technique for PT measurements. Considering its relation with thermal wave theory, as described in [7], the (dirac) pulse excitation of a PT measurement can be decomposed into a multitude of individual sinusoidal waves with different amplitudes and frequencies [31]. The benefit of a short dirac equivalent pulse is found in the flat frequency distribution of the discrete Fourier transformed equivalent spectrum. This method is a combination of the advantages of lock-in thermography with the speed of PT. As the dirac pulse is a combination of multiple sine waves, a fast examination can be performed for multiple depths through the structure. By computing the blind frequency f b and using Equation (1), it is possible to quantitatively estimate the defect depth.
z = C 1 α π · f b = C 1 · μ with α = κ / C p ρ .
f b is the blind frequency, and C 1 is a constant value typically equal to 1.82 [32].
The maximum penetration depth in a material is defined by the thermal diffusion length μ = α π · f b , which can be translated to the measured depth according to Equation (1). Defects are defined by their blind frequency f b . The lower the f b value is, the deeper into the material one can probe for defects.

3.2. Modal Analysis Theory

Multiple-degree-of-freedom (MDOF) systems can be expressed by the following equation [33]:
[ M x ¨ ( t ) + C x ˙ ( t ) + K x ( t ) = f ( t ) ]
with M , C , and K are the dynamic mass, damping, and stiffness matrices, and x ( t ) is the displacement of the MDOF system as a function of time.
Transforming Equation (2) to the Laplace domain yields the transfer function matrix H ( s ) between displacement and force vectors, X ( s ) = H ( s ) F ( s ) :
H ( s ) = [ M s 2 + C s + K ] 1 = N ( s ) d ( s )
where common-denominator polynomial d ( s ) = d e t ( M s 2 + C s + K ) , also known as the characteristic polynomial.
When the damping is small, the roots of the characteristic polynomial d(s) are complex conjugate pole pairs, λ m and λ m , m = 1, …, N m , where N m is the number of modes in the system. The transfer function H ( s ) can be rewritten in a pole-residue form, i.e., the so-called “modal” model:
H i k ( ω ) = = 1 N m R r j ω λ r + R r j ω λ r
with complex conjugate residues R r , R r , and poles:
λ , λ = σ ± j ω .
The pole is a complex number with numerical value of its real part ( σ ), the damping value, and its imaginary part, the damped natural (modal) frequency ( ω ).
Figure 7 shows a typical frequency response function measured on the bicycle frame. For a frequency range between 1 and 6000 Hz, 10 modes are selected. The modal parameters are estimated [26]. The estimated resonance frequencies of the performed measurements can be found in Table 2.

3.3. Finite Element Model Updating

Making use of the modal model described in [14], the frequency response function H i k ( ω ) between input DOF k and output DOF i at angular frequency ω can be written as
H i k ( ω ) = r = 1 N m Q r ψ i r ψ k r j ω λ r + Q r ψ i r ψ k r j ω λ r
where N m is the number of modes, Q r is the modal scaling factor, λ r is the pole number, and r, ψ i r and ψ k r are the mode shapes for mode r and ∗ for the complex conjugate.
We will use eigenfrequencies and mode shapes to correlate and update the model parameters to the real, optimized values.
In this section, the automated updating method [22] used to update the finite element frequencies based on two updating parameters will be briefly explained, using poles (resonance frequencies) and corresponding mode shapes as reference data. A non-linear least squares problem is solved in order to minimize an objective function consisting of the sum of the squares of the differences in corresponding resonance frequencies and mode shapes between the test modes and FE modes. In the objective function, the resonance frequencies and mode shapes are both measured but are weighted differently, taking their importance and value into account [34].
The parameter updating problem is stated as
minimize ( ϑ ) with ( ϑ ) = i = 1 n m W f i 2 f i ϑ F E f i E 2 + i = 1 n m W ϕ i 2 j = 1 n p ϕ j i ϑ F E ϕ j i E 2 subject to 0 ϑ 1
where f i and ϕ j i are the i t h natural frequency and associated mode shape at measurement point j, respectively. W f and W ϕ are the weight coefficients of the frequency and mode shape, respectively; n m denotes the number of modes and n p the number of measurement points. Subscripts F E and E denote the items from the FE model and experiments, respectively, ϑ denotes the normalized parameter vector, which allows for variation within the interval [0,1], and ( ϑ ) is a real valued scalar objective function.

3.4. Adaptive Response Surface Method

The optimization routine makes use of adaptive response surface optimization to optimize numerical models based on a multitude of measured data points. The routine is designed to handle multiple-output time series data for the thermal case and frequency data for the modal case [24,27]. The optimization procedure can be divided into the following steps:
  • Starting reference simulation points are selected, and an adapted object function is defined by taking into account the difference between the FE model values and the optimal target values.
  • The FE model is replaced by an interpolation-based meta-model of response surfaces to decrease optimization time, but this replacement remains an accurate approximation.
  • The optimization routine is run on the specific object function to find the minimum of the meta-model. An optimization routine is used according to a genetic algorithm computation to find the optimum speed of the complex numerical model.
  • The estimated parameter values are used as input parameters for an improved FE model that corrects the meta-model.
  • Only the points that are closer to the minimum are used to form the response surface (pan & zoom function).
The optimized response model is not built from a pre-defined number of design experiments but is adapted and refined during the optimization routine by the pan & zoom command [24]. The method automatically calculates the parameter influence sensitivity and sequentially resolves multiple sets with major influencing parameters and fixed parameter values for minor influencing parameters. Afterward, the second parameter set is optimized, and an influence check is then performed for all parameter sets until the results converge.

3.5. Optimized Parameter Choice

3.5.1. Thermal Case

The thermal excitation power, excitation duration, and the thermal camera frame rate are used as optimization parameters. The laminate thickness is used as an objection function, and there are physical boundary conditions as well. These limits are shown in Table 3.

3.5.2. Modal Case

The local thickness of each element is optimized using the modal frequencies (global information) estimated from the structural FRFs as well as the mode shapes (local information) measured using the LD-vibrometry measurements. It should be noted that the measurements only apply in one degree of freedom, as the flexure amplitude is only measured perpendicular to the left side of the bicycle frame.

3.6. Objective Function

The objective function locates the minimum of the response surface. The minimization is defined in Equation (8), where ( ϑ ) is defined as the objective function, and T is the response parameter measuring the relative temperature or the resonance frequencies. The response parameter is dependent on the optimized parameter ϑ . The objective function is tested if the boundary conditions are not violated.
( ϑ ) = T m o d e l ( ϑ ) T E x .

4. Results and Discussion

With the numerical model, it is possible to estimate the exposed heating flash, optimally define the heating power for the thermal case, and adjust the mass and structural material properties for the modal case. In this way, the tube thickness of the frame can be measured by the response, considering an anisotropic fiber distribution in the bicycle structure.
The numerical model used in the optimization routine is shown in Figure 8. The model consists of three different material mesh collectors: one for the cover paint with a thickness of 0.11 mm and two for the CFRP plies with a consecutive 90 degree angle between each ply and with the respective global orthotropic ply material properties. In the initial assumption, the full bicycle frame has a continuous thickness of four CFRP plies in a regular stacking of [ 0 / 90 ] s . The material properties and heat exchange properties are defined using the numerical updating process as described and compared with a preliminary static measurement.

4.1. Thermal Case

One time instance of the thermal results is shown in Figure 9, where an inter- and intralaminar crack through all plies including the surrounding delaminations is shown. The thickness of this area was difficult to compute due to the crack and type of delamination. The crack showed a thickness (crack depth) of 0 mm, which is below the lower limit. In addition, the delamination was nearly as thick as the minimal thickness. In reality, this thickness is basically the local thickness of the tube, which was further investigated.
The local thickness of different points was examined in comparison to certain benchmark points that could be physically measured. The measurements were performed with an accuracy of 20 μ m and are shown in Table 4. It can be seen that the estimated thickness at the three inspected points was accurate, and this was validated by physical measurements. This confirms the two values that cannot be verified without destroying the bicycle and confirms the assumption that a measurement inspection routine can be defined using advanced path planning in combination with thermal FE routines for complex-shaped composite materials.

4.2. Model Case

The measured FRF of the used mesh elements of the saddle tube is shown in Figure 7. The global resonance frequencies of the bicycle frame are listed in Table 2. The first natural bending resonance frequency occurred at 102.8 Hz. The excitation shaker was fixed in the middle of the frame at the saddle tube. After the model updating routine was performed, the updated material thicknesses of the evaluated points was changed according to the results in Table 4.
As a result of the updating procedure, an average error of 4.2% in relative difference was achieved. By taking the mode shape in the objective function into account (Equation (7)), three modes were switched in the updating process, demonstrating the usefulness of the mode shapes and not just the resonance frequencies in the objective function. If this value were not included, the wrong frequency value differences for these modes would have been minimized.

4.3. Comparison

By comparing the updated thickness values of the thermal and structural updating (Table 4), one comes to the conclusions that a thermal model update delivers more precise results than a structural model update. However, some remarks should be made:
  • The mode shapes are only measured in one degree of freedom instead of three, which results in missing mode shape data along the two other degrees of freedom. The main reason for this consists in the complexity of 3D non-contact vibration analysis. However, this problem does not occur when thermal data is used as a reference.
  • By performing model updating using vibrational data to evaluate local material thickness, it can be seen that the combination of local measurement and global responses biases the local thickness measurement.
  • Although estimating the thermal properties as structural properties is more complex, the thermal results appear to be better matched.
  • The thermal updating routine is far more time-consuming in contrast to modal updating. However, the thermal measurements on the structure are easier to perform.
Modal updating versus thermal updating appears to be comparatively complex when local parameters such as material thickness are evaluated. Considering complex-shaped structures with a high variation in thickness, measurement time during thermal inspection can be feasibly reduced by investing in extensive computation time. This is especially true when a series of similar products is inspected, as the thermal properties are less variable than structural properties during production. The main reason for this is the use of local temperature measurements in model updating. This was demonstrated by the use of local thermal material properties to evaluate local material thickness (Table 4). When the structure contains a wide variety of materials containing high-temperature materials with little conduction, different conclusions may be drawn. Further investigation is needed to make conclusions about those type of structures.

5. Conclusions

We can assume that the use of correlation and updating routines is of major importance in improving the predictive quality of numerical simulations. It has been demonstrated that techniques similar to numerical model updating can be used for both thermal updating and modal updating. However, an automatic updating routine based on mode shape correlation is thermally not available. The results demonstrate that similar accuracy can be achieved using thermal measurements for updating purposes. In addition, the link between temperature response and material thickness can be approximated using an analytic equation, which facilitates the definition of normalized parameter values in the objective function.

Acknowledgments

This research was funded by the University of Antwerp and the Institute for the Promotion of Innovation by Science and Technology in Flanders (VLAIO) via support for the TETRA project, “SINT: Smart Integration of Numerical modeling and Thermal inspection,” project number HBC.2017.0032. The researchers received funding from the Antwerp University IOF-council through project PSID-34924 entitled “Fast Broadband Lock-In Thermography for Fragile Structures Using System Identification.” The authors would like to end with special thanks to the MiViM research Chair, prof. Xavier Maldague, and express gratitude for the support of the full team in supporting some of the preliminary measurements of this research.

Author Contributions

J.P. and B.R. conceived and designed the experiments; J.P., G.S., B.R., and S.V. performed the experiments; G.S. and S.V. analyzed the data; B.R. and G.S. contributed materials and analysis tools; J.P., B.R., and G.S. wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest. The founding sponsors had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Chung, D.D.L. Carbon Fiber Composites; Butterworth-Heinemann: Oxford, UK, 1994; Available online: https://www.sciencedirect.com/science/article/pii/B9780080500737500105 (accessed on 2 December 2012). [CrossRef]
  2. Meola, C.; Boccardi, S.; Carlomagno, G.M.; Boffa, N.D.; Monaco, E.; Ricci, F. Nondestructive evaluation of carbon fibre reinforced composites with infrared thermography and ultrasonics. Compos. Struct. 2015, 134, 845–853. [Google Scholar] [CrossRef]
  3. Czichos, H. (Ed.) Handbook of Technical Diagnostics; Springer: Berlin, Germany, 2013; pp. 3–9. [Google Scholar]
  4. Peeters, J. Optimised Failure Detection and Advanced Non-Destructive Inspection Using Active IR Thermography and Numerical Modelling. Ph.D. Thesis, University of Antwerp, Antwerp, Belgium, 2017. [Google Scholar]
  5. Hellier, C. Handbook of Nondestructive Evaluation; McGraw-Hill Education: New York, NY, USA, 2001. [Google Scholar]
  6. Vollmer, M.; Möllmann, K. Infrared Thermal Imaging: Fundamentals, Research and Applications; Wiley-VCH: Berlin, Germany, 2010. [Google Scholar]
  7. Maldague, X. Theory and Practice of Infrared Thermography for Nondestructive Testing; Wiley: New York, NY, USA, 2001; p. 684. [Google Scholar]
  8. Chaudhuri, P.; Santra, P.; Yoele, S.; Prakash, A.; Reddy, D.; Lachhvani, L.; Govindarajan, J.; Saxena, Y. Non-destructive evaluation of brazed joints between cooling tube and heat sink by IR thermography and its verification using FE analysis. NDT E Int. 2006, 39, 88–95. [Google Scholar] [CrossRef]
  9. Jeong, W.; Earls, C.; Philpot, W.; Zehnder, A. Inverse thermographic characterization of optically unresolvable through cracks in thin metal plates. Mech. Syst. Signal Process. 2012, 27, 634–650. [Google Scholar] [CrossRef]
  10. Susa, M.; Ibarra-Castanedo, C. Pulse thermography applied on a complex structure sample: Comparison and analysis of numerical and experimental results. In Proceedings of the IV Conferencia Panamericana de END, Buenos Aires, Argentina, 22–26 October 2007; p. 12. [Google Scholar]
  11. Lau, S.; Almond, D.; Milne, J. A quantitative analysis of pulsed video thermography. NDT E Int. 1991, 24, 195–202. [Google Scholar] [CrossRef]
  12. Usamentiaga, R.; Ibarra-castanedo, C.; Klein, M.; Maldague, X. Nondestructive Evaluation of Carbon Fiber Bicycle Frames Using Infrared Thermography. Sensors 2017, 17, 2679. [Google Scholar] [CrossRef] [PubMed]
  13. Bowkett, M.; Thanapalan, K. Comparative analysis of failure detection methods of composites materials’ systems. Syst. Sci. Control Eng. 2017, 5, 168–177. [Google Scholar] [CrossRef]
  14. Friswell, M.I.; Mottershead, J.E. Finite Element Model Updating in Structural Dynamics; Springer: Dordrecht, The Netherlands, 1995; Volume 38. [Google Scholar]
  15. Marwala, T. Finite-Element-Model Updating Using Computational Intelligence Techniques; Springer: Berlin/Heidelberg, Germany, 2010; p. 263. [Google Scholar]
  16. Sun, R.; Chen, G.; He, H.; Zhang, B. The impact force identification of composite stiffened panels under material uncertainty. Finite Elem. Anal. Des. 2014, 81, 38–47. [Google Scholar] [CrossRef]
  17. Steenackers, G. Finite Element Model Updating and Optimization of Mechanical Systems Making Use of Regressive Techniques. Ph.D. Thesis, Vrije Universiteit Brussel, Brussels, Belgium, 2008. [Google Scholar]
  18. Abd El-Ghany, K.M.; Farag, M.M. Expert system to automate the finite element analysis for non-destructive testing. NDT E Int. 2000, 33, 409–415. [Google Scholar] [CrossRef]
  19. Foti, D.; Diaferio, M.; Giannoccaro, N.I.; Mongelli, M. Ambient vibration testing, dynamic identification and model updating of a historic tower. NDT E Int. 2012, 47, 88–95. [Google Scholar] [CrossRef]
  20. Meo, M.; Zumpano, G. Damage assessment on plate-like structures using a global-local optimization approach. Optim. Eng. 2008, 9, 161–177. [Google Scholar] [CrossRef]
  21. Mottershead, J.; Friswell, M. Model Updating In Structural Dynamics: A Survey. J. Sound Vib. 1993, 167, 347–375. [Google Scholar] [CrossRef]
  22. Steenackers, G.; Guillaume, P. Finite element model updating taking into account the uncertainty on the modal parameters estimates. J. Sound Vib. 2006, 296, 919–934. [Google Scholar] [CrossRef]
  23. Cheng, H.C.; Tsai, Y.H.; Chen, K.N.; Fang, J. Thermal placement optimization of multichip modules using a sequential metamodeling-based optimization approach. Appl. Therm. Eng. 2010, 30, 2632–2642. [Google Scholar] [CrossRef]
  24. Peeters, J.; Arroud, G.; Ribbens, B.; Dirckx, J.; Steenackers, G. Updating a finite element model to the real experimental setup by thermographic measurements and adaptive regression optimization. Mech. Syst. Signal Process. 2015, 64–65, 428–440. [Google Scholar] [CrossRef]
  25. Peeters, J.; Ibarra-Castanedo, C.; Khodayar, F.; Mokhtari, Y.; Sfarra, S.; Zhang, H.; Maldague, X.; Dirckx, J.; Steenackers, G. Optimised dynamic line scan thermographic detection of CFRP inserts using FE updating and POD analysis. NDT E Int. 2018, 93, 141–149. [Google Scholar] [CrossRef]
  26. Heylen, W.; Lammens, S.; Sas, P. Modal Analysis Theory and Testing; K.U. Leuven: Leuven, Belgium, 1998. [Google Scholar]
  27. Steenackers, G.; Presezniak, F.; Guillaume, P. Development of an adaptive response surface method for optimization of computation-intensive models. Comput. Ind. Eng. 2009, 57, 847–855. [Google Scholar] [CrossRef]
  28. Maldague, X.; Galmiche, F.; Ziadi, A. Advances in pulsed phase thermography. Infrared Phys. Technol. 2002, 1, 175–181. [Google Scholar] [CrossRef]
  29. Sfarra, S.; Ibarra-Castanedo, C.; Lambiase, F.; Paoletti, D.; Di Ilio, A.; Maldague, X. From the experimental simulation to integrated non-destructive analysis by means of optical and infrared techniques: Results compared. Meas. Sci. Technol. 2012, 23, 115601. [Google Scholar] [CrossRef]
  30. Vavilov, V.; Taylor, R. Theoretical and practical aspects of the thermal nondestructive testing of bonded structures. Res. Tech. NDT Test. 1982, 5, 240–279. [Google Scholar]
  31. Maldague, X.; Marinetti, S. Pulse phase infrared thermography. J. Appl. Phys. 1996, 79, 2694. [Google Scholar] [CrossRef]
  32. Thomas, R.L.; Pouch, J.J.; Wong, Y.H.; Favro, L.D.; Kuo, P.K.; Rosencwaig, A. Subsurface flaw detection in metals by photoacoustic microscopya. J. Appl. Phys. 1980, 51, 1152–1156. [Google Scholar] [CrossRef]
  33. Guillaume, P. Modal Analysis. In Control Systems, Robotics and Automation—VOL V; Unbehauen, H., Ed.; Hindawi Publishing Corporation: Cairo, Egypte, 2009; p. 558. Available online: http://mech.vub.ac.be/avrg/publications/ModalAnalysis.pdf (accessed on 1 January 2010).
  34. Hao, H.; Xia, Y. Vibration-based damage detection of structures by genetic algorithm. J. Comput. Civ. Eng. 2002, 16, 222–229. [Google Scholar] [CrossRef]
Figure 1. Photograph of the thermography measurement setup.
Figure 1. Photograph of the thermography measurement setup.
Applsci 08 00307 g001
Figure 2. Photograph of the modal analysis measurement setup.
Figure 2. Photograph of the modal analysis measurement setup.
Applsci 08 00307 g002
Figure 3. Overview of the 2D mesh of the bicycle frame.
Figure 3. Overview of the 2D mesh of the bicycle frame.
Applsci 08 00307 g003
Figure 4. Close-up of the 3D mesh with the different plies schematically represented by different colors.
Figure 4. Close-up of the 3D mesh with the different plies schematically represented by different colors.
Applsci 08 00307 g004
Figure 5. A Xenon flash excitation source output of 6.4 kJ.
Figure 5. A Xenon flash excitation source output of 6.4 kJ.
Applsci 08 00307 g005
Figure 6. Schematic view of a pulsed thermography inspection, modified with permission from MiViM. 1: Flash lamps, 2: specimen, 3; camera, 4: processing computer.
Figure 6. Schematic view of a pulsed thermography inspection, modified with permission from MiViM. 1: Flash lamps, 2: specimen, 3; camera, 4: processing computer.
Applsci 08 00307 g006
Figure 7. Frequency response function of the expermimental LDV measurement.
Figure 7. Frequency response function of the expermimental LDV measurement.
Applsci 08 00307 g007
Figure 8. Numerical finite element model result after 1 s. Temperatures in C.
Figure 8. Numerical finite element model result after 1 s. Temperatures in C.
Applsci 08 00307 g008
Figure 9. Thermal response of the crack after PCT (Bin 4).
Figure 9. Thermal response of the crack after PCT (Bin 4).
Applsci 08 00307 g009
Table 1. Material properties used within the numerical model.
Table 1. Material properties used within the numerical model.
ValueUnits
Young’s modulus E165GPa
Young’s modulus E2, E314.9GPa
Poisson’s ratio NU12, NU230.058
Poisson’s ratio NU130.024
Shear modulus G12, G133.5MPa
Shear modulus G231.4MPa
Mass density1570Kg/m 3
Thermal conductivity60(W/m·K)
Thermal conductivity4.2(W/m·K)
Specific heat capacity1000(J/(Kg·K))
Coefficient of convection10W/(m 2 ·K)
Emissivity0.85
Table 2. Optimization results.
Table 2. Optimization results.
NoTest (Hz)FE_init (HZ)Relative Difference (%)FE_opt (Hz)Relative Difference (%)
1102.81393.581 8.98 98.252 4.35
2138.438107.621 22.26 131.874 4.74
3206.563190.744 7.66 200.493 2.94
4242.500233.786 3.59 243.1520.27
5278.125282.7031.65295.0336.08
6367.500310.999 15.37 342.514 6.80
7382.188325.293 14.89 359.383 5.97
8454.375525.07615.56477.1425.01
9569.688377.917 33.66 484.327 14.98
10593.125397.058 33.06 512.527 13.59
Table 3. Thermal optimization parameter limits.
Table 3. Thermal optimization parameter limits.
Lower LimitUpper Limit
Excitation power (W)100500
Excitation duration (s)2.020
Frame rate (Hz)5.050
Table 4. Optimization results.
Table 4. Optimization results.
Experimental Geometric Measured Value (mm)Thermal Model Computation (mm)Error (%)Structural Model Computation (mm)Error (%)
Saddle tube (top)2.642.53 4 , 2 2.681,5
Head tube (top)4.054.050,03.19 21 , 2
Lower tube (middle)2.052.092,01.59 22 , 4
Top Tube (middle)x2.28x2.65x
Rear Tube (middle)x3.55x2.21x

Share and Cite

MDPI and ACS Style

Steenackers, G.; Peeters, J.; Verspeek, S.; Ribbens, B. From Thermal Inspection to Updating a Numerical Model of a Race Bicycle: Comparison with Structural Dynamics Approach. Appl. Sci. 2018, 8, 307. https://doi.org/10.3390/app8020307

AMA Style

Steenackers G, Peeters J, Verspeek S, Ribbens B. From Thermal Inspection to Updating a Numerical Model of a Race Bicycle: Comparison with Structural Dynamics Approach. Applied Sciences. 2018; 8(2):307. https://doi.org/10.3390/app8020307

Chicago/Turabian Style

Steenackers, Gunther, Jeroen Peeters, Simon Verspeek, and Bart Ribbens. 2018. "From Thermal Inspection to Updating a Numerical Model of a Race Bicycle: Comparison with Structural Dynamics Approach" Applied Sciences 8, no. 2: 307. https://doi.org/10.3390/app8020307

APA Style

Steenackers, G., Peeters, J., Verspeek, S., & Ribbens, B. (2018). From Thermal Inspection to Updating a Numerical Model of a Race Bicycle: Comparison with Structural Dynamics Approach. Applied Sciences, 8(2), 307. https://doi.org/10.3390/app8020307

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