NDT.net • Aug 2005 • Vol. 10 No.8 
Finite element method applied to the three dimension thermal characterization of a delaminationA. Obbadi, S. Belattar*Laboratory of Energetic and processing of the signal U.F.R Instrumentation and processing of the signal Faculty of Sciences, 24000 El Jadida, Morocco *Corresponding Author Contact: Email: belattars@ucd.ac.ma, Internet: AbstractIn this work we present the application of a 3D numerical method of heat conduction for the determination and characterization of subsurface defects in material. The numerical model is based on the finite elements method. The object of this study is to investigate, the influence of position of defect in the structure and the thickness of this later on the thermal response. Some configurations are simulated in which we have taken defects in form of cylinders containing air placed at different positions. In this case twelve artificial cylindrical defects were taken and are machined in a steel metal specimen. A heat source is applied to the input surface of metal. The thermal interaction of defects on the material surface temperature was analyzed and its influence on the current contrast determined with the aim to recognize the possible determination of defect geometry and its position.1. IntroductionThermal nondestructive testing (TNDT) is a method of inspection and of analysis of homogeneity of various structures by simple acquisition of information at the possible accesses of considered system. It plays a fundamental role in all domains where the requirements of safety are important, like nuclear industry, aeronautical, automobile, railway... The techniques used are varied (thermography [1], [2], [3], ultrasounds, eddy currents, xrays, radiography, sweating, visual inspections...).. The TNDT is usually used to seek defects in metallurgical parts [2] or composites [3]. In this search of quality thermics can play a very important part. Many thermal nondestructive testing methods were used. The most interesting are the active methods [5], [6], [7] which consist in subjecting material to be controlled with various thermal excitations, then to observe their relaxation. The heat source can be thermal [1],[2],[3], [4], [5], [6], [7] or optics. The temperature can be measured by contact (micro manufactured resistors, thermocouples, sensors ) or without contact (infrared thermography, acoustic detection...). The static methods are more sensitive to the losses by radiation than the dynamic methods. The optical methods are very sensitive to the surface quality of the sample. Among all these existing control methods, we chose to study thermal method "TNDT" [4], [5], [6], [7] resting on detection of a possible disturbance of the thermal field which appears when in the structure to be inspected is induced a permanent or transitory heating. The disturbance is due to a local or total deterioration of the thermal properties relating to the nominal distribution. This deterioration can be the fruit of a physicochemical modification of one of the components (effect of ageing), or the consequence of the appearance of a foreign body: inclusion [4] introduced accidentally, separation [1] or delamination [2], [4] associated to appearance of a chink.2. Description of the wallIn order to illustrate the application of the TNDT method we present the results of the non destructive testing of a standard sample (steel metal) (fig.1) of thickness e = 20 mm, length L = 156mm, and width l=146mm, containing 12 equidistant defects in honeycomb. The defects have a cylindrical form of diameter d = 9, 18 and 36 mm, height h = 15, 17, 18 and 19 mm, located at the position l_{1} = 1, 2, 3, and 5 mm from the entry face (fig.2). Lines A_{2}A
3. Mathematical description of the model, boundary and initial conditionsTo solve the following thermal equation [5], [6], [7]:The ratio is called thermal diffusivity. We call upon the numerical method of the finite elements [6], [8]. The analytical resolution is indeed impossible being given the geometry of the problem. The method consists in using an approximation by finite elements of the unknown functions T to discretize the variational form of the equation (1) and to transform it into system of algebraic equations of the form: With:
F a vector of N components T the vector of the temperatures to be calculated
4. Resolution of the equations4.1. Description of the gridWe present the whole wall to emphasize the grid density around the defect (fig. 4). We chose a grid made up of triangular elements. Its density increases when one is around the defects (fig. 4). The latter is simply simulated by an absence of matter.
5. Results of simulationsIn order to illustrate the previous theoretical considerations, we present the computation results of the thermal response in the case of an isotropic material, steel in this case, characterized by K = 58W/m.k (thermal conductivity) = 7800kg/m^{3} (density) and C = 460J/kg.k (specific heat) containing delamination characterized by K = 0,0272W/m.k (thermal conductivity) r = 1,057kg/m^{3} (density) and C = 717,8J/kg.k (specific heat). The diameters d of the defects are equal to 9, 18, and 36mm and the positions of the front face l_{1} are equal to 1, 2, 3, and 5mm. After resolution of the considered problem, it is possible to plot the temperature distribution on all or a part of the wall at a given moment, as well as the temporal evolution of the temperature in a given point (fig. 5). In the presence of the defect the heat flow has tendency to propagate by avoiding the defect as showed in figures 5 and 6. This phenomenon explains the rise of the temperature at the place of the defect, represented by a hot thermal patch fig. 5 and 6. The maximum of this temperature gives an estimate on the required resolution of the non destructive testing equipment. At the exit of the defect the flow lines tend to be uniform. This could be information on the form and the position of the defect in material.
5.1. Influence of the defect parametersLet us interest in the influence of paramount parameters, namely the defect diameter d, and defect position l_{1}. In this study, we consider the case of delamination with the following thermophysical characteristics: K = 0,0272W/m.k (thermal conductivity) = 1,057kg/m^{3} (density) and C = 717,8J/kg.k (specific heat), and located at depths l_{1} of 1, 2, 3, and 5mm of the front face.
5.1.1. Influence of defect diameter
5.1.2. Influence of defect position by report to the front face
6. ConclusionIn this work, we studied the case of a material containing a cylindrical delamination subjected to a thermal stress. We carried out a systematic study of the temperature evolution according to the intrinsic thermophysical characteristics of the defect. This enabled us to conclude that, on the assumption of the less deep defects, diameter and position play a determining role in the thermal response of delaminated material. By studying the influence of the defect parameters on the measurable magnitude, we showed that it is theoretically possible to detect any defect, with the proviso of applying a sufficient energy of excitation and that of the defect, it is different in practice. Indeed, this model relates to a resistive defect in a rigorously plane plate. However, if one introduces a light initial curve (what is practically always the case in reality), one realizes that the heat gradient which exists between the heated face and the back face is at the origin of a total deformation of the plate which can completely occult the deformation at the location of defect. This phenomenon which is very often observed while a measurement makes that it very difficult to detect the defects whose diameter is very low. All calculations were carried out in the case of an isotropic material, but the taking into account of the anisotropy would be possible with the proviso of using a computer code which allows it.References

© NDT.net  Top 