Microsoft Word - numero_38_art_48 G. S. Serovaev et alii, Frattura ed Integrità Strutturale, 38 (2016) 392-398; DOI: 10.3221/IGF-ESIS.38.48 392 Numerical study of the response of dynamic parameters to defects in composite structures Grigorii S. Serovaev, Valery P. Matveenko Perm National Research Polytechnic University, Department of Applied Mathematics and Mechanics, 29 Komsomolsky prospekt, Perm, 614990, Russia serovaev@icmm.ru, mvp@icmm.ru ABSTRACT. The purpose of this work is to study three different models of delamination in composite plate (the free mode model, the constrained model and the contact model) and applicability of this models to the vibrational method of damage detection based on frequency shifts. The results of numerical simulation have shown that the free mode model leads to abrupt changes in natural frequencies due to non-physical condition of mutual penetration of adjacent volumes in the defect zone. The constrained and the contact models yield qualitative agreement in shifts of natural frequencies with a change of the defect size. All models have shown the necessity of analyzing shifts of high frequencies to detect small size delamination. KEYWORDS. Natural frequencies; Vibrational methods of damage detection; Delamination. Citation: Serovaev, G.S., Matveenko, V. P., Numerical study of the response of dynamic parameters to defects in composite structures, Frattura ed Integrità Strutturale, 38 (2016) 392-398. Received: 28.08.2016 Accepted: 15.09.2016 Published: 01.10.2016 Copyright: © 2016 This is an open access article under the terms of the CC-BY 4.0, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. INTRODUCTION he estimation of the mechanical state of a structure during its operation based on detection of internal damages has gained a great importance for such rapidly developing fields of production as mechanical engineering, aerospace industry, etc., whose recent trends are toward an increased use of new types of materials. The most widely spread type of defects in composite constructions is delamination, which leads to a separation of the object in the defect zone into several parts, which begin to respond to loads independently of each other [1] while the stiffness of each part is significantly lower than the stiffness of the whole object. The appearance of delamination may be the cause of the deprivation of bearing capacity of the structure, its breakage and decommissioning which indicates the relevance of the problem of damage detection at an early stage of development. The application of visual methods of damage detection is difficult due to internal location of delamination. It is well known that the appearance of defects in a particular place of the structure influences the local stiffness in this zone, which leads to a shift of natural frequencies. This specific feature of structure response is behind the global vibrational methods of damage detection. These methods use different dynamic parameters (eigenfrequencies, mode shapes, curvatures of mode shapes, etc.) to evaluate the state of the investigated structure [2-5]. T G. S. Serovaev et alii, Frattura ed Integrità Strutturale, 38 (2016) 392-398; DOI: 10.3221/IGF-ESIS.38.48 393 The methods of numerical simulation allow us to investigate the processes occurring in structures during application of different damage detection methods and receive the answers to many questions without performing expensive experiments. Moreover, it is due to numerical methods that we can track the evolution of different parameters with a change of the defect size. Studies of objects with delamination type defects have received considerable attention in scientific literature. Analytical description of the free mode model of delamination in a beam, which allows the mutual penetration of volumes in the zone of delamination, is given in [6]. This disadvantage is eliminated in the constrained model where areas in the delamination zone have equal vertical displacements [7]. Only few low frequencies are considered in these studies. The presence of multiple delaminations is discussed in [8]. Experimental investigations of vibrations of delaminated structures are performed in [9, 10]. Excitation of vibrations and measurement of the signal must be implemented using actuators and special measurement devices (sensors). The capabilities of using piezoelements for registration of high frequency vibrations are described in [11]. The Authors of this research capture vibrations of the plate up to 40 kHz with the help of a piezoelectric sensor which proves the possibility of using such devices for measuring high frequency vibrations. The method of electromechanical impedance (EMI), which is one of the vibrational methods of damage detection applied to delaminated composite beam, is given in [12]. The need to place the actuator close to the defect and hence the requirement of arrangement of a dense grid of piezoelectric devices or prior knowledge of the location of defect is the main disadvantage of such a kind of damage detection method. On the other hand natural frequencies give integral characteristic of the object of research. In [13-15] vibrations of delaminated structures with different geometries such as beam, cylindrical and conical shells are studied. A numerical study of the dynamic parameter response to defects of different sizes is carried out in the framework of three models of delamination in a composite structure, the advantages and drawbacks of each model are estimated in the process of simulation. First, we consider a free mode model of delamination, in which the adjacent volumes in the zone of the defect are not coupled with each other and therefore they are assumed to be mutually penetrable. In the second, the so-called constrained model, the coincident nodes in the zone of delamination are coupled by one component of the displacement, while other components of the displacement remain independent. The last one is the model that takes into account the contact forces. The free mode and constrained models allow us to perform a modal analysis for computing the eigenfrequencies of the structure. With contact forces taken into account the problem becomes nonlinear. The algorithm for calculating the spectrum of eigenfrequencies in the nonlinear problem includes the following steps: setting of impact load, transient analysis, measurements of signal at a certain point of the structure with a subsequent conversion of the received response of the structure from the amplitude-time to amplitude-frequency domain with the help of the Fourier Transform NUMERICAL MODELING numerical study was performed on a square plate of size L = 0.15х0.15 m, made of a layered composite material. The thickness of one layer is h1c = 0.0003 m. There are a total of n = 15 layers across the thickness of the plate therefore the total thickness is h1c*n = 0.0045 m (Fig. 1). The square shape delamination of size Ld m is located between 6 and 7th layers. The center of delamination coincides with the center of the plate. The composite material is modeled as a homogeneous body with the following orthotropic effective mechanical properties (reinforcing direction is not taken into account). Ex = 24 GPa, Ey = 18 GPa, Ez = 6 GPa, Gxy = 4 GPa, Gyz = 3 GPa, Gxz = 3 GPa, vxy = 0.15, vyz = 0.18, vxz = 0.42, ρ = 1800 kg/m3. The presence of composite layers is taken into account only by geometrical dimensions. Finite elements with quadratic approximation providing more accurate results in terms of out of plane strains were used in the simulation. The principle of virtual work is used for the mathematical formulation of the problem i i i i i ij ij V S V u F u dV p u dS dV t 2 2 (             (1) where iu - components of the displacement vector, ji ij j i uu x x 1 ( ) 2       - components of strain tensor, ij ijkl klC   - components of stress tensor, ijklC - stiffness tensor,  - material density, iF - components of body forces vector, ip - components of surface tractions vector. A G. S. Serovaev et alii, Frattura ed Integrità Strutturale, 38 (2016) 392-398; DOI: 10.3221/IGF-ESIS.38.48 394 Figure 1: Geometrical representation of the plate. The boundary condition corresponds to the cantilever fixing of the plate on the border x = L The free mode model allows the mutual penetration of volumes (Fig. 2), hence surfaces S1 and S2 are free from stresses and boundary conditions for this models have the following form: zz zx zyS S S S S S, , ,1 2 1 2 1 2 0     For the constrained model, the coincident nodes associated with the surfaces S1 and S2 have equal displacements in the z axis direction while other components of the displacement remain independent therefore the boundary conditions can be written as follows zx zyS S S S z zS S U U , ,1 2 1 2 1 2 0    A numerical solution of the dynamic problem of vibration of a plate was found by the finite element method using the commercial package ANSYS. The finite-element formulation in matrix form can be written as follows        M u K u F t( )  (2) where  M - mass matrix of the system,  K - stiffness matrix of the system,  F t( ) - time dependent load function,  u - the vector of nodal displacements,  u - the vector of nodal accelerations. In the absence of external influences, the problem is reduced to the typical problem of finding eigenvalues and eigenvectors i t i iu x t e x( , ) ( )    (3) where  - natural frequency, i - eigenvector or mode shape. The finite element formulation of the problem has the following form        K M u2 0 0   (4) where  u0 - the vector of nodal values of natural vibration modes. G. S. Serovaev et alii, Frattura ed Integrità Strutturale, 38 (2016) 392-398; DOI: 10.3221/IGF-ESIS.38.48 395 Figure 2: Models of delamination. RESULTS he frequency range from 0-20 kHz was analyzed. The size of delamination was changed from 0 to 0,05 m with a step of 0.002 m. The change of eigenfrequencies, depending on the size of delamination in the range of 0-20 kHz, is shown in Fig. 3 and 4 (the upper plots show the results for the free mode model of delamination, the lower – for constrained model). The vertical axis represents the defect size, the horizontal axis corresponds to the frequency of oscillations in Hertz. The reaction of the natural frequencies to increase in the size of the defect for the free mode model is not observed up to 4 kHz. At higher frequencies a sudden drop in the natural frequencies for the defect of the large size is observed. Figure 3: Shifts of natural frequencies in the range of 0-10 kHz. In the range of frequencies from 10-20 kHz a similar abrupt pattern of natural frequency shifts is observed (Fig. 4). Restriction applied to the component of displacement in the delamination zone and thereby precluding mutual penetration of the adjacent volumes in this area significantly affects the results. This model is characterized by a more gradual and smaller change in magnitude of natural frequencies with increasing the size of the defect. For both models a clear tendency towards greater sensitivity of high frequency vibrations to the defects of a smaller size is observed, which is explained by the fact that the reaction of frequencies to the appearance of a defect depends not only on its size but also on its location. If a defect is located in the area of small or zero strains of mode shape, the frequency will remain unchanged. The mode shapes of high frequencies have a large number of zones with non-zero strains (bends), providing greater sensitivity of these frequencies to the defect. That's what makes these frequencies more appropriate for detection of small sized defects. T G. S. Serovaev et alii, Frattura ed Integrità Strutturale, 38 (2016) 392-398; DOI: 10.3221/IGF-ESIS.38.48 396 Figure 4: Shifts of natural frequencies in the range of 10-20 kHz. The third variant of the model of delamination is the most accurate and complete, as it uses contact elements to model the interaction of adjacent surfaces in the delamination zone. Thus, several options for the contact status are possible (Fig. 5): 1) open contact where the nodes, between which the contact is defined, are divided; 2) closed contact, in this case, the nodes are in contact with each other and their relative behavior is determined by the contact stiffness, by default, equal to the modulus of elasticity of the material multiplied by the size of the element adjacent to the contact surface. In this problem the node to node type of contact (element conta178) is used because mode superposition transient analysis, which supports only this type of nonlinearity, was performed. A significant advantage of this type of calculation is its performance, compared to the full method of transient analysis. Figure 5: State of the contact. In contrast to the free mode model and constrained model where natural frequencies could be calculated by the modal analysis, the spectrum of frequencies in transient analyses is calculated using the Fourier Transform of a signal recorded at a certain point of the plate. Since accelerometers are usually used as sensors to measure the acceleration in the experiment, in the numerical model the accelerations are also calculated at the point. Increasing the size of the centrally located delamination, the graphs similar to ones in Fig. 3 and 4 are received which reflect the change of the natural frequencies of the spectrum when you change the size of the defect. This type of analysis requires a more profound approach, as a reflection of the resonant frequencies on the spectrum significantly depends on the type of the input signal, the points of impact and measurement. If these points are located on the nodal line of mode shape, this frequency will not be excited and will not affect the resulting spectrum. The algorithm of calculation of the nonlinear dynamic problem with contact interaction taken into account consists of the following steps (Fig. 6). Impact load is applied to a certain area on the surface of the plate (step I) G. S. Serovaev et alii, Frattura ed Integrità Strutturale, 38 (2016) 392-398; DOI: 10.3221/IGF-ESIS.38.48 397 m oF t Sin f t Sin f t( ) ( 2 ) ( 2 )  where om f f n2  – modulation frequency, of – central frequency. On the time interval T = 0.4 s. mode superposition transient analysis is performed. The result of this analysis is the component of acceleration az registered at the node at the specified time period (step II). The received signal includes a set of steady-state oscillations of different frequencies, which is obtained by decomposition of the signal using the Fourier Transform (step III). Figure 6: Algorithm of transient analysis Resonance peaks on the resulting graph correspond to the eigenfrequencies. Internal dissipation was not considered in the current study, therefore, the analysis of the amplitudes of the resonance peaks was not performed. By implementing this algorithm, and increasing the size of the defect at each step, it is possible to monitor changes in natural frequencies. Figure 7: Shifts of natural frequencies in the contact model of delamination. The general nature of the frequency response to an increase of the damage size is similar to the models described above. We can see the lack of influence of damage on low frequencies and high frequency sensitivity to the defect of a small size. G. S. Serovaev et alii, Frattura ed Integrità Strutturale, 38 (2016) 392-398; DOI: 10.3221/IGF-ESIS.38.48 398 The sharpness of the frequency spectrum picture will significantly depend on relative location of impact and the measurement and type of impact load. CONCLUSIONS he comparison of the three considered models of delamination with respect to their applicability to the solution of the problem of analysis of changes of natural frequencies of a composite plate shows that in spite of the indisputable simplicity of the free mode model, the results show an excessive drop in eigenfrequencies. The constrained and the contact models yield qualitative agreement in shifts of natural frequencies with a change of the defect size. The calculation of the constrained model is much easier and faster. However, the algorithm of calculation of the model with the contact repeats the steps performed during the experiment and more fully reflects special features typical for real structures. The results suggest that for detection of defects at an early stage of their development it is necessary to record the change in the spectrum in the high frequency range (in this case greater than 4 kHz). A numerical model of the studied structure allows to analyse frequency response to the occurrence of the defect and to determine the most sensitive of them to a specific type of defect. ACKNOWLEDGMENTS he study was performed in Perm National Research Polytechnic University with the support of the Russian Science Foundation (project №15-19-00243) REFERENCES [1] Zhang, Z., Shankar, K., Ray, T., Morozov, E.V., Tahtali, M., Vibration-based inverse algorithms for detection of delamination in composites, Composite Structures, 102 (2013) 226-236. [2] Stepinski, T., Uhl, T., Staszewski, W., Advanced Structural Damage Detection: From Theory to Engineering Applications, John Wiley & Sons, (2013). [3] Adams, D.E., Health Monitoring of Structural Materials and Components, John Wiley & Sons, (2007). [4] Morassi, A., Vestroni, F., Dynamic methods for damage detection in structures. Springer Wien New York, (2008). [5] Zou, Y., Tong, L., Steven, G.P., Vibration-based model-dependent damage (delamination) identification and health monitoring for composite structures – a review. Journal of Sound and Vibration, 230(2) (2000) 357-378. [6] Wang, J.T.S., Liu, Y.Y., Gibby, J.B., Vibrations of split beams. Journal of Sound and Vibration, 84 (1982) 491-502. [7] Mujumdar, P.M., Suryanarayan, S., Flexural vibrations of beams with delaminations, Journal of Sound and Vibration, 125(3) (1988) 441-461. [8] Lee, J., Free vibration analysis of delaminated composite beams. Computers and Structures, 74 (2000) 121-129. [9] Shen, M.H., J.E. Grady. Free vibrations of delaminated beams. AIAA Journal. 30(5) (1992) 1361-1370. [10] Valdes Diaz, S.H., Soutis, C., Delamination detection in composite laminates from variations of their modal characteristics. Journal of Sound and vibration, 228(1) (1999) 1-9. [11] Chuang, K.-C., Liou, H.-C., Ma, C.-C., Investigation of Polyvinylidene Fluoride(PVDF) Films in Identifying High- Frequency Vibration Modes of Flexible Plates. IEEE Transactions on Ultrasonic, Ferroelectrics, and Frequency Control, 61(6) (2014) 1047-1058. [12] Yan, W., Wang, J., Chen, W.Q., Delamination assessment of a laminated composite beam using distributed piezoelectric sensor/actuator, Smart. Mater. Struct., 20 (2011) 1-14. [13] Rout., M., Baishya, N., Effects of delamination on the vibration characteristics of composite beams, Noize and Vibration Worldwide, 24-29 (2010). [14] Muc, A., Stawiarski, A., Identification of damages in composite multilayered cylindrical panels with delaminations. Composite Structures. 94 (2012) 1871-1879. [15] Dey, S., Karmakar, A., Free vibration analyses of multiple delaminated angle-ply composite conical shells – A finite element approach, Composite Structures. 94 (2012) 2188-2196. T T << /ASCII85EncodePages false /AllowTransparency false /AutoPositionEPSFiles true /AutoRotatePages /None /Binding /Left /CalGrayProfile (Dot Gain 20%) /CalRGBProfile (sRGB IEC61966-2.1) /CalCMYKProfile (U.S. Web Coated \050SWOP\051 v2) /sRGBProfile (sRGB IEC61966-2.1) /CannotEmbedFontPolicy /Error /CompatibilityLevel 1.4 /CompressObjects /Tags /CompressPages true /ConvertImagesToIndexed true /PassThroughJPEGImages true /CreateJobTicket false /DefaultRenderingIntent /Default /DetectBlends true /DetectCurves 0.0000 /ColorConversionStrategy /CMYK /DoThumbnails false /EmbedAllFonts true /EmbedOpenType false /ParseICCProfilesInComments true /EmbedJobOptions true /DSCReportingLevel 0 /EmitDSCWarnings false /EndPage -1 /ImageMemory 1048576 /LockDistillerParams false /MaxSubsetPct 100 /Optimize true /OPM 1 /ParseDSCComments true /ParseDSCCommentsForDocInfo true /PreserveCopyPage true /PreserveDICMYKValues true /PreserveEPSInfo true /PreserveFlatness true /PreserveHalftoneInfo false /PreserveOPIComments true /PreserveOverprintSettings true /StartPage 1 /SubsetFonts true /TransferFunctionInfo /Apply /UCRandBGInfo /Preserve /UsePrologue false /ColorSettingsFile () /AlwaysEmbed [ true ] /NeverEmbed [ true ] /AntiAliasColorImages false /CropColorImages true /ColorImageMinResolution 300 /ColorImageMinResolutionPolicy /OK /DownsampleColorImages true /ColorImageDownsampleType /Bicubic /ColorImageResolution 300 /ColorImageDepth -1 /ColorImageMinDownsampleDepth 1 /ColorImageDownsampleThreshold 1.50000 /EncodeColorImages true /ColorImageFilter /DCTEncode /AutoFilterColorImages true /ColorImageAutoFilterStrategy /JPEG /ColorACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /ColorImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000ColorACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000ColorImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasGrayImages false /CropGrayImages true /GrayImageMinResolution 300 /GrayImageMinResolutionPolicy /OK /DownsampleGrayImages true /GrayImageDownsampleType /Bicubic /GrayImageResolution 300 /GrayImageDepth -1 /GrayImageMinDownsampleDepth 2 /GrayImageDownsampleThreshold 1.50000 /EncodeGrayImages true /GrayImageFilter /DCTEncode /AutoFilterGrayImages true /GrayImageAutoFilterStrategy /JPEG /GrayACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /GrayImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000GrayACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000GrayImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasMonoImages false /CropMonoImages true /MonoImageMinResolution 1200 /MonoImageMinResolutionPolicy /OK /DownsampleMonoImages true /MonoImageDownsampleType /Bicubic /MonoImageResolution 1200 /MonoImageDepth -1 /MonoImageDownsampleThreshold 1.50000 /EncodeMonoImages true /MonoImageFilter /CCITTFaxEncode /MonoImageDict << /K -1 >> /AllowPSXObjects false /CheckCompliance [ /None ] /PDFX1aCheck false /PDFX3Check false /PDFXCompliantPDFOnly false /PDFXNoTrimBoxError true /PDFXTrimBoxToMediaBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXSetBleedBoxToMediaBox true /PDFXBleedBoxToTrimBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXOutputIntentProfile () /PDFXOutputConditionIdentifier () /PDFXOutputCondition () /PDFXRegistryName () /PDFXTrapped /False /CreateJDFFile false /Description << /ARA /BGR /CHS /CHT /CZE /DAN /DEU /ESP /ETI /FRA /GRE /HEB /HRV (Za stvaranje Adobe PDF dokumenata najpogodnijih za visokokvalitetni ispis prije tiskanja koristite ove postavke. Stvoreni PDF dokumenti mogu se otvoriti Acrobat i Adobe Reader 5.0 i kasnijim verzijama.) /HUN /ITA /JPN /KOR /LTH /LVI /NLD (Gebruik deze instellingen om Adobe PDF-documenten te maken die zijn geoptimaliseerd voor prepress-afdrukken van hoge kwaliteit. De gemaakte PDF-documenten kunnen worden geopend met Acrobat en Adobe Reader 5.0 en hoger.) /NOR /POL /PTB /RUM /RUS /SKY /SLV /SUO /SVE /TUR /UKR /ENU (Use these settings to create Adobe PDF documents best suited for high-quality prepress printing. Created PDF documents can be opened with Acrobat and Adobe Reader 5.0 and later.) >> /Namespace [ (Adobe) (Common) (1.0) ] /OtherNamespaces [ << /AsReaderSpreads false /CropImagesToFrames true /ErrorControl /WarnAndContinue /FlattenerIgnoreSpreadOverrides false /IncludeGuidesGrids false /IncludeNonPrinting false /IncludeSlug false /Namespace [ (Adobe) (InDesign) (4.0) ] /OmitPlacedBitmaps false /OmitPlacedEPS false /OmitPlacedPDF false /SimulateOverprint /Legacy >> << /AddBleedMarks false /AddColorBars false /AddCropMarks false /AddPageInfo false /AddRegMarks false /ConvertColors /ConvertToCMYK /DestinationProfileName () /DestinationProfileSelector /DocumentCMYK /Downsample16BitImages true /FlattenerPreset << /PresetSelector /MediumResolution >> /FormElements false /GenerateStructure false /IncludeBookmarks false /IncludeHyperlinks false /IncludeInteractive false /IncludeLayers false /IncludeProfiles false /MultimediaHandling /UseObjectSettings /Namespace [ (Adobe) (CreativeSuite) (2.0) ] /PDFXOutputIntentProfileSelector /DocumentCMYK /PreserveEditing true /UntaggedCMYKHandling /LeaveUntagged /UntaggedRGBHandling /UseDocumentProfile /UseDocumentBleed false >> ] >> setdistillerparams << /HWResolution [2400 2400] /PageSize [612.000 792.000] >> setpagedevice