Microsoft Word - numero_35_art_32 N. Oudni et alii, Frattura ed Integrità Strutturale, 35 (2016) 278-284; DOI: 10.3221/IGF-ESIS.35.32 278 Focussed on Crack Paths Implementation of a damage model in a finite element program for computation of structures under dynamic loading Nasserdine Oudni Department of civil engineering, faculty of technology, A. Mira University, Bejaia, Algeria Laboratory of modeling of materials and structures, University of Tizi-Ouzou, Algeria nasroudni@yahoo.fr Youcef Bouafia Laboratory of modeling of materials and structures, University of Tizi-Ouzou, Algeria youcef.bouafia2012@yahoo.com ABSTRACT. This work is a numerical simulation of nonlinear problems of the damage process and fracture of quasi-brittle materials especially concrete. In this study, we model the macroscopic behavior of concrete material, taking into account the phenomenon of damage. J. Mazars model whose principle is based on damage mechanics has been implemented in a finite element program written Fortran 90, it takes into account the dissymmetry of concrete behavior in tension and in compression, this model takes into account the cracking tensile and rupture in compression. It is a model that is commonly used for static and pseudo-static systems, but in this work, it was used in the dynamic case. KEYWORDS. Damage; Dynamic; Modeling; Displacement; Strain; Local approach; Time integration; explicit; quasi-brittle. INTRODUCTION he numerical program is employed in analyzing Koyna gravity dam and comparing the results with available results of many researchers find in literature. The choice of Koyna concrete gravity dam is due to the fact that many studies have been done on this structure which experienced a devastating earthquake in 1967. Since Koyna dam in the 1967 Koyna earthquake [1-4] is one of the few examples of cracking in a concrete dam, it has been selected to investigate its earthquake response. The earthquake records are in the form of accelerograms with horizontal and vertical components which are used as dynamic loads. The effect of the damping coefficient is obviously taken into account. This work proposes the implementation of a damage model based on the damage mechanics as part of the isotropic formulation proposed by J. Lemaitre in a finite element program. the J. Mazars model [5] whose principle is based on damage mechanics which is a theory describing the progressive reduction of the mechanical properties of a material due to initiation, growth and coalescence of microscopic cracks. These internal changes lead to the degradation of mechanical properties of the material. The model takes into account the asymmetry of concrete behavior and it considers the cracking in tension and compression failure. The method is based on the notion of equivalence of the strains. T N. Oudni et alii, Frattura ed Integrità Strutturale, 35 (2016) 278-284; DOI: 10.3221/IGF-ESIS.35.32 279 MODEL he influence of microcracking due to external loads is introduced via a single scalar damage variable d ranging from 0 for the undamaged material to 1 for completely damaged material. The stress-strain relation reads [6]:     0 0 ij ij kk ij 0 0 1 υ υ ε σ σ δ      E 1 d E 1 d        (1) 0E and 0 are the Young's modulus and the Poisson's ratio of the undamaged material; ijε and ijσ are the strain and stress components, and ijδ is the Kronecker symbol. The elastic (i.e., free) energy per unit mass of material is   0 1 1      2 ij ijkl kld C    (2) Where 0ijklC is the stiffness of the undamaged material. This energy is assumed to be the state potential. The damage energy release rate is 0 1   2 ij ijkl klY C d         (3) With the energy of dissipated energy: d d         (4) Since the dissipation of energy ought to be positive or zero, the damage rate is constrained to the same inequality because the damage energy release rate is always positive. DAMAGE EVOLUTION he evolution of damage is based on the amount of extension that the material is experiencing during the mechanical loading. An equivalent strain is defined as   3 2 1    i i       (5) Where .  is the Macauley bracket and i are the principal strains. The loading function of damage is  ,     f       (6) Where  is the threshold of damage growth. Initially, its value is 0 , which can be related to the peak stress tf of the material in uniaxial tension: 0 0   t f E   (7) In the course of loading  assumes the maximum value of the equivalent strain ever reached during the loading history. T T N. Oudni et alii, Frattura ed Integrità Strutturale, 35 (2016) 278-284; DOI: 10.3221/IGF-ESIS.35.32 280           ,  0       ,  0,  If f and f Then        0           0,       0 d h d With d else                 (8) The function  h  is detailed as follows: in order to capture the differences of mechanical responses of the material in tension and in compression, the damage variable is split into two parts:    t t c cd d d   (9) Where td and cd are the damage variables in tension and compression, respectively. They are combined with the weighting coefficients t and c , defined as functions of the principal values of the strains t ij and c ij due to positive and negative stresses:    1 1  1 ,      1  t t c cij ijkl kl ij ijkl kld C d C         (10) 3 3 2 2 1 1 ,     t c i i i i t c i i                            (11) Note that in these expressions, strains labeled with a single indicia are principal strains. In uniaxial tension 1t  and 0c  . In uniaxial compression 1c  and 0t  . Hence, td and td can be obtained separately from uniaxial tests. The evolution of damage is provided in an integrated form, as a function of the variable  :     0 0 1 1   exp t t t t A A d B            (12)     0 0 1 1   exp c c c c A A d B            . (13) Figure 1. Evolution of two parts of damage td and cd [5]. A direct tensile test or three point bend test can provide the parameters which are related to damage in tension ( 0 , tA ,   tB ). Note that Eq. 7 provides a first approximation of the initial threshold of damage, and the tensile strength of the material can be deduced from the compressive strength according to standard code formulas. The parameters ( cA , cB ) are fitted from the response of the material to uniaxial compression. N. Oudni et alii, Frattura ed Integrità Strutturale, 35 (2016) 278-284; DOI: 10.3221/IGF-ESIS.35.32 281 DYNAMIC EQUATIONS OF MOTION ynamic analysis of structures exhibiting non linear behavior is performed by using direct integration, to trace the response in the time domain. The nonlinear dynamic equilibrium equation can be written as   n n n nM Cu p fu     (14) where M and C are the global mass and damping matrices respectively, np is the global vector of internal resisting nodal forces, nf is the vector of consistent nodal forces for the applied body and surfaces traction forces grouped together, the body force term ( gMIu  ) due to seismic excitation, is included in the body forces which are taken into account in n f , nu is the global vector of nodal accelerations and nu is the global vector of nodal velocities [7]. Figure 2: Uniaxial response of the model in (a) tension and (b) compression [5] When the structure is subjected to seismic excitation, the external applied body forces is n gf MIu   (15) Where gu is ground acceleration and I is a vector indicating the direction of the earthquake excitation. Figure 3: Koyna accelerograms a) Transverse component b) Vertical component [8]. D N. Oudni et alii, Frattura ed Integrità Strutturale, 35 (2016) 278-284; DOI: 10.3221/IGF-ESIS.35.32 282 The geometry of a typical non-overflow monolith of the Koyna dam is illustrated in Fig. 4. The monolith is 103 m high and 71 m wide at its base. The upstream wall of the monolith is assumed to be straight and vertical, which is slightly different from the real configuration. The depth of the reservoir at the time of the earthquake is hw = 91.75 m. Following the work of other investigators, we consider a two-dimensional analysis of the non-overflow monolith assuming plane stress conditions. The finite element mesh used for the analysis is shown in Fig. 4. It consists of 760 first-order, reduced- integration, plane stress elements (CPS4R). Nodal definitions are referred to a global rectangular coordinate system centered at the lower left corner of the dam, with the vertical y-axis pointing in the upward direction and the horizontal x- axis pointing in the downstream direction. The transverse and vertical components of the ground accelerations recorded during the Koyna earthquake are shown in Fig. 3. (units of g = 9.81 m sec–2). Prior to the earthquake excitation, the dam is subjected to gravity loading due to its self-weight and to the hydrostatic pressure of the reservoir on the upstream wall [8].       Figure 4: Geometric properties of the KOYNA dam. Parameters used in Mazars’s model for concrete Young modulus E 31027 MPa Poisson’s ratio  0.2 Mass density  2643 kg/m3 Initial damage threshold 0D  1.5 10-4 At 1.0 Bt 30000 AC 1.4 BC 1545 Table1: Material properties.     RESULTS e note that the displacements are relatively low during the first two seconds because of low the amplitudes of the excitations. The displacements reach their maximum at 3.7 s and 7.5 s, 30 mm was recorded at 3.8 s, the maximum displacement value does not correspond to the maximum amplitude of the excitement that is recorded at 3.65 s. The nodal displacements decrease after 7.5 s.   W N. Oudni et alii, Frattura ed Integrità Strutturale, 35 (2016) 278-284; DOI: 10.3221/IGF-ESIS.35.32 283   Figure 5: Horizontal displacements at the crest of the dam We note that the displacements are relatively low during the first two seconds because of low the amplitudes of the excitations. The displacements reach their maximum at 3.7 s and 7.5 s, 30 mm was recorded at 3.8 s, the maximum displacement value does not correspond to the maximum amplitude of the excitement that is recorded at 3.65 s. The nodal displacements decrease after 7.5 s.   Figure 6: Damage evolution, tD and cD . Fig. 6. Shows the evolution of tD and cD , the tensile and compressive part of the damage variable D . We note that traction damage is more significant and more predominant than in compression, which leads to a brittle fracture of the material by traction. The total rupture of the material can be achieved by traction, in compression the damage is less significant, so the rupture will not occur under compression, Fig. 1b presents some ductility of the material. It is also seen from this figure (Fig. 6) that the tensile damage starts firstly and evolves quickly to a value close to 1 ( 1tD  ), then followed by the compression damage, but the temporary difference is only of the order of a thousandth of a second, so nearly or completely insignificant difference. In counterparty, compressive damage is less important, it increases until a value 0.15cD  a value that does not allow rupture in compression of the element. Damage were observed (see Fig. 7) in the dam after 3.53 s at the integration point (815) of the element 204, then from 3.54 s at integration point (811) of the element 203 and finally, from 3.54 s at Gauss point (807) of the element 202. It may be noted that the evolution of the damage is mainly concentrated in the time interval where the maximum values of positive and negative displacements occur. N. Oudni et alii, Frattura ed Integrità Strutturale, 35 (2016) 278-284; DOI: 10.3221/IGF-ESIS.35.32 284   Figure 7: Damage evolution in Gauss points. CONCLUSION e have taken the example of a concrete gravity dam subjected to seismic excitation in the form of accelerograms. The results are, first the response of gravity dam subjected to seismic loading as displacements history, secondly the history of the evolution of integration points damage at the cracked or damaged zone. The dynamic equilibrium equation is solved, an integration algorithm by the method of central differences for nonlinear systems, we have used a very small time step (0.0001 s). The choice of weight Koyna dam in India is the fact that many studies have been conducted on this structure to a known destructive earthquake in 1967. The records of the earthquake accelerograms as with horizontal and vertical components were used as dynamic loads.  REFERENCES [1] Calayir, Y., Karaton, M., A continuum damage concrete model for earthquake analysis of concrete gravity dam– reservoir systems, Soil Dynamics and Earthquake engineering, 25 (2005) 857-869. [2] Jianwen, P., Chuhann, Z., Xuyanjie, Feng, J., A comparative study of the different procedures for seismic cracking analysis of concrete dams, Soil Dynamics and Earthquake engineering, 31 (2011) 1594-1606. [3] Calayir, Y., Karaton, M., A continuum damage concrete model for earthquake analysis of concrete gravity dam– reservoir systems. Soil Dynamics and Earthquake Engineering, 25 (2005) 857–869. [4] Omidi, O., Valliappan, S., Lotfi, V., Seismic cracking of concrete gravity dams by plastic–damage model using different damping mechanisms. Finite Elements in Analysis and Design, 63 (2013) 80–97. [5] Mazars, J., Application de la mécanique de l’endommagement au comportement non linéaire et à la rupture du béton de structure, Doctoral thesis, University of Paris 6, France (1984). [6] Mazars, J., Pijaudier-Cabot, G., Continuum damage theory- application to concrete, J. Engrg. Mech. ASCE, 115 (1989) 345-365. [7] Owen, D.R.J., Hinton, E., Finite element in plasticity, Theory and Practice. Ed. Pineridge Press Limited, (1986). [8] ABAQUS 6.10. Example problems manual; volume I: Static and Dynamic Analyses. W << /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