外文翻译-准脆性材料在动荷载作用下结构的数值计算：以混凝土重力坝为例

外文翻译题目南沟门常态混凝土重力坝及溢流坝段设计专业水利水电工程2015年准脆性材料在动荷载作用下结构的数值计算以混凝土重力坝为例摘要本文使用有限元数值计算方法对动态加载后的混凝土结构进行计算,混凝土结构进行动态加载通常情况下类似于自然地震的效果，其后果往往是毁灭性的。如今在动态的一般方程中,已经建立了基于损伤力学且考虑了混凝土材料的性质的模型。从现象学的角度看,选择行为模型考虑造成的损害的微裂纹。本文采用1967年KOYNA地震的记录，通过数值应用程序,对二维KOYNA重力坝进行抗震分析。对混凝土重力坝的地震反应破坏的影响进行了研究。关键词法律行为，混凝土，损伤，动态加载，开裂，模型。1简介准脆性材料的损伤与断裂过程是属于非线性问题，特别是对混凝土材料性质的模拟，是研究混凝土材料在土木工程领域中应用的脆弱性研究（LEMAITREANDMAZARS,1982LEMAITREANDCHABOCHE,1985）。模拟过程中的现象往往是复杂的，即使在简单的情况下，问题的解析解也是困难的。至少，要知道材料的性质，模拟的结构需进行静态和或动态负荷。混凝土被认为是设计中最常用的材料，其可能会出现在强动载荷的结构中。无论是意外情况或人为的，它都属于一类比较复杂且具有非线性性质的非均质材料。大量存在的文献中的实验测试表明该材料的性质为“准脆性的材料的拉伸强度明显低于压缩，在高速度对混凝土的冲击试验（KLEPACZKO和BRARA,2001HERVA和GATUINGT,2002HENTZETAL,2004）或循环加载试验（DUB1994LABORDERIE1991）中表现很明显。这些实验的目的是预测宏观行为模型的材料损伤原理。在这些模型中，JMAZARS模型的原理是基于损伤力学理论，描述由于起始材料的力学性能逐渐减少，从而产生显微裂纹的生长和聚结这些内在的变化导致材料力学性能的退化。该模型考虑了混凝土的不对称性质和裂纹在拉伸和压缩时破坏的程度。这是一个通常用于静态或伪静态系统模型。它对强度和安全监管技术要求非常严格，混凝土结构或钢筋混凝土结构的大小的动态载荷是属于地震类型的载荷。一些研究已经进行了在这个方向来考虑结构承受地震荷载时对结构损伤的影响（CALAYIR和KARATON2005OMIDIETAL2012DAVENNEETAL2003）结构设计与混凝土坝的抗震课题属于动态加载的类型。动态横向和纵向的加速度输入，在1967年11月KOYNA地震资料中，我们将考虑了每平方米的表面结构的静力效应应用到大坝墙，但不影响流体动力效果的阻尼系数是在计算中考虑。该材料是基于由JMAZARS1984在当地研发的各向同性模型它是一个在拉伸和压缩时损伤不同的模型，可以准确地描述该行为，它也允许二次损坏，其压缩刚度会减少。这种模式已经在用90种FORTRAN语言，旨在与非线性材料模型专门处理有限元程序得到落实。该程序的功能是为了解决动态的一般方程，从而使我们能够计算位移的节点，应变，应力和损坏的高斯积分点。时间积分方案称为显性积分，采用动态的平衡方程的数值解与小的时间间隔。尤其是当我们采用差分法它是基于速度和加速度与位移与时间的已知值的差分的商数的近似。解决方法的一个条件是稳定的，必须具有非常小的时间间隔。另一种方法是有限元方法，该方法是基于连续的描述物质，它可以被应用到数值模拟准脆性材料的降解处理的线性和非线性的问题。2动力平衡方程对于运动物体的动态平衡，在时间站的运动方程中，给出下面的NT表达式（1）（OWEN和HINTON1986）NFPUCM（1）其中M和C是总体质量和阻尼矩阵，是内部节点力的总体矢量，是NP所有的节点力向量组合在一起，身体力项（MUG由于地震，包括在有了考虑到身体的力量，N是节点的加速度和总体向量是节点速度NF矢量的总和。运动方程（1）可以通过使用中心差分方程如下重写（2）。121122NNNNUCTMUFPTCTMU（2）为了简便求解方程（2），质量矩阵和阻尼矩阵M、C做对角矩阵变换。在这些假设下，方程（2）可以改写为一个标量方程（3）（OWEN和HINTON1986；PAULTRE2005）。IINIIINININICTMUCTUFPTU221121（3）方程式（3）的稳定性与（T）的时间间隔大小程度有关。分析其内部抵抗力可由下式（4）给出，DBPITIN（4）其中和是有效应力和整体应变的位移矩阵。I3阈值函数损伤加载函数F（，D），根据损伤变量值而引入。这种损伤阈值函数定义其过程是可逆的，只要当F（，D）0时，损伤就不会增加。所以，损伤是在一个给定的状态下，其加载的函数是（5）式。5加载函数的变量为相关加载的历史系数,为等效应变,主DKI要取决于压力给定的表达式6,比如，如果，则和如I00果，则损伤D生长，当等效应变达到一个阈值KD0I0初始。如果6,0,DKTHENDF破坏由MAZARS分为两个部分。CTD（7）剪切实验的参数它通常被认为是常数105。和损伤变TDC量分别表示张力和压缩。破坏的过程以一个集成的形式表示,作为变量的函数。0EXP1,0,DBAADCTCTDCT（8）4KOYNA大坝的损伤反应KOYNA混凝土重力坝,高103米,宽702米，如图1所示,位于KOYNA河印度半岛的西部。1967年,该地区发生65级地震与，其基础面坝轴线水平方向和垂直方向的最大加速度分别为049G和034G。KOYNA地震的时间历时记录图2所示。观测发现大坝上方的上游和下游出现严重破坏的水平裂缝。本文对大坝进行非线性动态分析，使用将混凝土材料当做各向同性的材料的模型进行观察损伤的方法。我们注意到大坝蓄水过程影响是不考虑。大坝部分的网格图如图2B所示。选择分析节点313N，此点的时间历时图表代表在坝顶垂直和水平运动。三个元素集成点807P,811P807P选择代表的时间历时损伤图。A几何性质B网格状的大坝图1在此研究中使用的材料性质如下表1中。表1材料特性弹性模量泊松比密度抗拉强度D0ATBTACBCE31027MPA20643GMKTFMPA241051020000141545本研究中使用的关键的粘性阻尼系数可以分别采取3、5和7。积分时间间隔是0001秒。动态加载历时是通过横向和纵向的加速度图时间记录或记录KOYNA地震图2所示。本文总结了地震造成的影响如位移、应变、应力和损伤,对重力坝结构图表的结果进行了讨论和比较。A垂直分量B水平分量图2KOYNA加速度图图3是位于坝顶节点313的垂直位移和水平位移的历时曲线。（A）垂直位移B水平位移图3坝顶节点313的位移历时曲线5讨论的结果在本文中,我们讨论了影响混凝土材料安全和破坏重力坝稳定的地震荷载。选用的阻尼比是5,在文中选用了现实例子进行计算。图3时间历时图显示了位于坝顶的的节点313）处垂直和水平位移情况。我们发现到前两秒期间位移相对较低,主要由于低振幅的影响。四秒到五秒的位移达到最大,462秒记录为20毫米,最大位移值与在365秒时的最大振幅点不对应。5S后的节点位移减少。在运动振幅处没有损失，损伤相对较低。损伤记录参见图4。在353秒后的大坝的集成点815达到最大,然后在354秒集成点811达到最大。可以看出,损伤的演变主要是集中在时间间隔正负位移的最大值处。结论相比在我们的研究中获得的结果CALAYIR和KARATON2005JIANWENETAL2011年相对一致。图4在各集成点处累积损伤记录815、811和807点6结论本研究的目的首先是研究当重力坝受到地震荷载时的位移,应变和应力随时间的变化过程其次是研究几个代表集成点处损伤演化的历史位移图，其损坏程度和起始结构有关。选择KOYNA混凝土重力坝的原因是因为,在这个大坝结构在1967年经历了一场毁灭性的地震后，对于这个大坝的许多研究已经完成。地震记录的水平和垂直加速度图形式和结构用作动态载荷时的性质相同。在这项研究中显然考虑了阻尼系数的影响，而不考虑水的动力作用。参考文献1CALAYIR,Y、KARATONM,2005年混凝土重力坝蓄水系统中混凝土连续损伤模型的地震分析土动力学和地震工程25200585720052DAVENNEL，RAGUENEAUF，玛泽J，IBRAHIMBEGOVICA，2003年有限元分析地震工程中的高效的方法。计算机和结构81（2003）1223年至1239年。3DUBJF，1994年混凝土结构建模，博士论文的简化粘易损行为，巴黎6S4HENTZS,DAUDEVILLEL,DONZFV,，动态加载在高应变离散元建模。计算机和结构82（2004）2509至24年5HERVG,GATUINGTF，损坏水泥板和混凝土受一个理想化的动态加载的数值模拟。LMT，ENSCACHAN。6JIANWENPCHUHANNZYANJIEX,FENGJ不同程序对1594年至1606年混凝土大坝土动力学与地震工程地震开裂分析比较研究。7KACHANOVLM，蠕变条件下破裂过程的时间，IZV。AKAD。NAUK。OTD。TEKH。NAUK。8（1958）2631。8KLEPACZKORJ,BRARAA2001年的实验方法混凝土剥落通过动态拉伸试验。影响工程的国际杂志。第25卷，第4期，2001年4月，页387409。9LABORDERIEC,1991在易损材料1991单方面现象建模和应用到混凝土结构，博士论文的分析，巴黎大学610LEMAITREJ,MAZARSJ1982理论的1982年的应用损坏结构混凝土的非线性行为和破裂。理论和计算249的方法。11LEMAITREJ,CHABOCHEJL,1985力学的固体材料,DUNOD,PARIS12MAZARSJ,1984机械损坏的混凝土结构的非线性行为和破损应用博士论文，大学皮埃尔与玛丽居里CNRS13OMIDIO,VALLIAPPANS,LOTFIV,，采用不同的阻尼机制由塑料损伤模型的混凝土重力坝地震2012年破裂。在分析和设计有限元632013809714OWENDRJ,HINTONE,1986年，理论与实践可塑性的有限要素。15PAULTREP，2005年土木工程中的结构动力学。爱马仕科学出版物。16PEYROTIBOUCHARDPO,BAYFBERNARDF,GARCIADIAZE2007年数值方面来模拟准脆性材料的力学行为。计算材料科学40（2007）327340。NUMERICALCOMPUTATIONOFSTRUCTURESWITHQUASIBRITTLEMATERIELUNDERDYNAMICLOADINGCASESTUDYOFCONCRETEGRAVITYDAMABSTRACTTHEPRESENTPAPERDEALSWITHANUMERICALCOMPUTATIONUSINGFINITEELEMENTMETHODFORCONCRETESTRUCTURESSUBJECTEDTODYNAMICLOADS,USUALLYOFSEISMICNATURETHECONSEQUENCESOFWHICHAREOFTENDEVASTATINGINTHEGENERALEQUATIONOFDYNAMIC,ALAWOFCONCRETEMATERIALBEHAVIORTAKINGINTOACCOUNTMODELSBASEDONTHEDAMAGEMECHANICSHASBEENIMPLEMENTEDFROMAPHENOMENOLOGICALVIEWPOINT,THECHOSENBEHAVIORMODELTAKESINTOACCOUNTTHEDAMAGECAUSEDBYTHEOPENINGOFMICROCRACKSFORNUMERICALAPPLICATION,TWODIMENSIONALSEISMICANALYSISOFKOYNAGRAVITYDAMISPERFORMEDBYUSINGTHE1967KOYNAEARTHQUAKERECORDSTHEEFFECTSOFDAMAGEONEARTHQUAKERESPONSEOFTHECONCRETEGRAVITYDAMARESTUDIEDKEYWORDSLAWOFBEHAVIOR,CONCRETE,DAMAGE,DYNAMICLOADING,CRACKING,MODEL1INTRODUCTIONTHESIMULATIONOFNONLINEARPROBLEMSOFDAMAGEANDFRACTUREPROCESSESFORQUASIBRITTLEMATERIALS,PARTICULARLYCONCRETE,ISTHESUBJECTOFRESEARCHSTUDIESINTHEFIELDOFCIVILENGINEERINGINTHEGOALTOSERVEFORVULNERABILITYSTUDIESLEMAITREANDMAZARS,1982LEMAITREANDCHABOCHE,1985THEPHENOMENATAKINGPARTINTHESEPROCESSESAREOFTENCOMPLEXAND,EVENINSIMPLECASES,THEANALYTICALRESOLUTIONOFTHEPROBLEMSTURNSOUTTOBE,FORTHELEAST,DIFFICULTTOSIMULATECOMPLEXSTRUCTURESSUBJECTEDTOSTATICAND/ORDYNAMICLOADS,ITISNECESSARYTOKNOWTHEBEHAVIOROFTHEMATERIALCONCRETEISCONSIDEREDASTHEMOSTUSEDMATERIALINTHEDESIGNOFSTRUCTURESWHICHMAYOCCASIONALLYBEEXPOSEDTOINTENSEDYNAMICLOADING,WHETHEROFACCIDENTALORINTENTIONALNATUREITBELONGSTOACLASSOFHETEROGENEOUSMATERIALSWITHQUITECOMPLEXNONLINEARBEHAVIORITISAQUASIBRITTLEMATERIALWITHATENSILESTRENGTHSIGNIFICANTLYLOWCOMPAREDTOTHECOMPRESSIONONENUMEROUSEXPERIMENTALTESTSSHOWINGTHEASPECTSOFTHEBEHAVIOROFTHISMATERIALEXISTSINTHELITERATURE,SUCHASTHEIMPACTTESTINGONCONCRETEATHIGHSPEEDKLEPACZKOANDBRARA,2001HERVANDGATUINGT,2002HENTZETAL,2004ORCYCLICLOADINGTENSIONCOMPRESSIONEXPERIMENTALTESTSDUB1994LABORDERIE1991THESEWORKSAIMSTOMODELTHEMACROSCOPICBEHAVIORANDPREDICTDAMAGEOFTHEMATERIALAMONGTHESEMODELS,THEJMAZARSMODELWHOSEPRINCIPLEISBASEDONDAMAGEMECHANICSWHICHISATHEORYDESCRIBINGTHEPROGRESSIVEREDUCTIONOFTHEMECHANICALPROPERTIESOFAMATERIALDUETOINITIATION,GROWTHANDCOALESCENCEOFMICROSCOPICCRACKSTHESEINTERNALCHANGESLEADTOTHEDEGRADATIONOFMECHANICALPROPERTIESOFTHEMATERIALTHEMODELTAKESINTOACCOUNTTHEASYMMETRYOFCONCRETEBEHAVIORANDITCONSIDERSTHECRACKINGINTENSIONANDCOMPRESSIONFAILURETHISISAMODELTHATISGENERALLYUSEDFORSTATICORPSEUDOSTATICSYSTEMSTHEREGULATORYTECHNICALDOCUMENTSFORREASONSOFSTRENGTHANDSAFETYREQUIRERIGOROUSLYTHATTHESTRUCTURESOFCONCRETEORREINFORCEDCONCRETETOBESIZEDUNDERDYNAMICLOADSESPECIALLYSEISMICTYPEFORDYNAMICINPUT,THETRANSVERSEANDVERTICALACCELERATIONCOMPONENTSOF11DECEMBER1967KOYNAEARTHQUAKEARESELECTEDWEWILLCONSIDERTHEHYDROSTATICEFFECTSOFTHEBASINPERSQUAREMETEROFSURFACEAPPLIEDTOTHEDAMWALL,BUTHYDRODYNAMICEFFECTSARENOTTAKENINTOCONSIDERATIONTHEEFFECTOFTHEDAMPINGCOEFFICIENTISTAKENINTOACCOUNTINTHECALCULATIONSTHEMATERIALISBASEDONTHEISOTROPICMODELDEVELOPEDBYJMAZARS1984INLOCALAPPROACHITISANEXPLICITDAMAGEMODELTHATCANACCURATELYDESCRIBETHEBEHAVIORINTENSIONANDCOMPRESSIONITALSOALLOWSREPRODUCINGTHEDAMAGEOFTENSION,COMPRESSIONANDDECREASEOFSTIFFNESSTHISMODELHASBEENIMPLEMENTEDINFINITEELEMENTPROGRAMWRITTENINFORTRAN90LANGUAGE,DESIGNEDTODEALSPECIFICALLYWITHNONLINEARMATERIALMODELINGTHEPROGRAMISDESIGNEDTOSOLVETHEGENERALEQUATIONOFDYNAMIC,CONSEQUENTLYALLOWINGUSTOCALCULATETHEDISPLACEMENTSATNODES,STRAINS,STRESSESANDDAMAGEATTHEGAUSSINTEGRATIONPOINTSTIMEINTEGRATIONSCHEMESCALLEDEXPLICITAREADOPTEDINNUMERICALSOLUTIONOFDYNAMICEQUILIBRIUMEQUATIONSFORDAMSTRUCTUREWITHSMALLTIMESTEPINPARTICULARWEADOPTACENTRALDIFFERENCEMETHODITISBASEDONAPPROXIMATIONOFVELOCITIESANDACCELERATIONSWITHQUOTIENTSOFFINITEDIFFERENCEOFKNOWNVALUESOFDISPLACEMENTSWITHREGULARTIMESTEPTHISMETHODISKNOWNCONDITIONALLYSTABLE,WITHAVERYSMALLTIMEINTERVALWEADOPTINTHISWORKTHEFINITEELEMENTMETHOD,THISMETHODISBASEDONACONTINUOUSDESCRIPTIONOFMATTER,ITCANBEAPPLIEDTONUMERICALSIMULATIONOFTHEDEGRADATIONOFQUASIBRITTLEMATERIALSFORTREATMENTOFLINEARANDNONLINEARPROBLEMS2DYNAMICEQUILIBRIUMEQUATIONSFORDYNAMICEQUILIBRIUMOFABODYINMOTION,THEEQUATIONOFMOTIONATTIMESTATIONTNISGIVENASTHEFOLLOWINGEXPRESSION1OWENANDHINTON19861NFPUCMWHEREMANDCARETHEGLOBALMASSANDDAMPINGMATRICESRESPECTIVELY,PNISTHEGLOBALVECTOROFINTERNALRESISTINGNODALFORCES,FNISTHEVECTOROFCONSISTENTNODALFORCESFORTHEAPPLIEDBODYANDSURFACESTRACTIONFORCESGROUPEDTOGETHER,THEBODYFORCETERM（MUGDUETOSEISMICEXCITATION,ISINCLUDEDINTHEBODYFORCESWHICHARETAKENINTOACCOUNTINFN,NISTHEGLOBALVECTOROFNODALACCELERATIONSANDISTHEGLOBALVECTOROFNODALVELOCITIESEQUATIONOFMOTION1CANBEREWRITTENBYUSINGCENTRALDIFFERENCEEQUATIONASFOLLOWING2121122NNNNUCTMUFPTCTMU（2）INORDERTOSOLVETHEEQUATION2EXPLICITLY,THEMASSMATRIXMANDTHEDAMPINGMATRIXCARETRANSFORMEDINDIAGONALMATRICESUNDERTHESEASSUMPTIONS,EQUATION2CANBEREWRITTENASASCALAREQUATION3OWENANDHINTON1986PAULTRE2005（3）THESTABILITYOFEQUATION3ISLINKEDTOTAVERYSMALLTIMEINTERVALTHEINTERNALRESISTINGFORCESUSINGINTHENUMERICALPROGRAMISGIVENBYTHEFOLLOWINGEXPRESSION4,（4）DBPITINWHEREANDBIARETHEEFFECTIVESTRESSANDTHEGLOBALSTRAINDISPLACEMENTMATRIX3THRESHOLDFUNCTIONADAMAGELOADINGFUNCTIONF,D,DEPENDINGONDAMAGEVARIABLE,ISINTRODUCEDTHISDAMAGETHRESHOLDFUNCTIONDEFINESTHEDOMAINWHERETHEBEHAVIORISREVERSIBLE,ASLONGASF,D0,THEDAMAGEDOESNOTINCREASESO,FORAGIVENSTATEOFDAMAGE,THELOADINGFUNCTIONIS（5）WHEREKDISTHEVARIABLERELATEDTOTHEHISTORYOFTHEDAMAGECALLEDEQUIVALENTSTRAINANDDEPENDSONMAINCALLEDEQUIVALENTSTRAINANDDEPENDSONMAINSTRAINSIGIVENBYTHEEXPRESSION6,SUCHAS0IFIIIFI0DAMAGEDGROWSWHENTHEEQUIVALENTSTRAINREACHESATHRESHOLDKDINITIALIZEDATD06DAMAGEDEFINEDBYMAZARSISSPLITINTOTWOPARTS7THEPARAMETERISREPRESENTATIVEOFSHEAREXPERIMENTSITISUSUALLYCONSIDEREDASCONSTANT105PIEROTETAL2007DTANDDCAREDAMAGEVARIABLESINTENSIONANDCOMPRESSIONRESPECTIVELYTHEEVOLUTIONOFDAMAGEISPROVIDEDINANINTEGRATEDFORM,ASAFUNCTIONOFTHEVARIABLE84DAMAGERESPONSEOFKOYNADAMTHEKOYNACONCRETEGRAVITYDAM,103MINHEIGHTAND702MINWIDTHSHOWNINFIG1,ISLOCATEDONTHEKOYNARIVERINTHEWESTOFTHEINDIANPENINSULAIN1967,A65MAGNITUDEEARTHQUAKESHOOKTHEREGIONWITHMAXIMUMACCELERATIONMEASUREDATTHEFOUNDATIONGALLERYOF049AND034GINHORIZONTALDIRECTIONNORMALTOTHEDAMAXISANDINTHEVERTICALONE,RESPECTIVELYTHETIMEHISTORIESRECORDSOFTHEKOYNAEARTHQUAKEARESHOWNINFIG2SEVEREDAMAGEWASFOUNDINTHEFORMOFHORIZONTALCRACKINGOBSERVEDONBOTHTHEUPSTREAMANDDOWNSTREAMFACESOFTHEUPPERPARTOFDAMMONOLITHSINTHISPAPER,THENONLINEARDYNAMICANALYSISOFKOYNADAMISPERFORMEDUSINGTHECONCRETEMODELWITHISOTROPICDAMAGEINLOCALAPPROACHWENOTETHATTHEDAMRESERVOIRINTERACTIONISNOTCONSIDEREDTHEMESHOFTHEDAMSECTIONISSHOWNINFIG2BNODEN313ISSELECTEDTOREPRESENTTHETIMEHISTORYGRAPHSOFVERTICALANDHORIZONTALMOVEMENTATTHEDAMCRESTTHETHREEELEMENTINTEGRATIONPOINTSP807,P811,P815ARESELECTEDTOREPRESENTTHETIMEHISTORYGRAPHOFDAMAGETHEMATERIALPROPERTIESUSEDINTHISSTUDYAREASFOLLOWINTHETABLE1TABLE1MATERIALPROPERTIESTHECRITICALVISCOUSDAMPINGCOEFFICIENTUSEDINTHISSTUDYCANTAKETHREEVALUESESTIMATEDAT3,5AND7INTEGRATIONTIMESTEPISTAKENASEQUALTO0001STHEDYNAMICEXCITATIONISBYACCELEROGRAMSWHOSEHORIZONTALANDVERTICALCOMPONENTSARETHETIMEHISTORIESORRECORDSOFTHEKOYNAEARTHQUAKESHOWNINFIG2THEEFFECTSCAUSEDBYTHEEARTHQUAKESUCHASDISPLACEMENTS,STRAINS,STRESSESANDDAMAGETOTHEGRAVITYDAMSTRUCTUREARESUMMARIZEDINGRAPHSTHATAREDISCUSSEDANDCOMPAREDTORESULTSOBTAINEDINTHELITERATUREFIGURE3SHOWTHETIMEHISTORYGRAPHSOFTHEVERTICALANDHORIZONTALDISPLACEMENTSOFTHENODALPOINT313LOCATEDATTHECRESTOFTHEDAM5DISCUSSIONOFRESULTSINTHISPART,WEWILLREPORTTHEEFFECTSOFTHECONCRETEDAMAGEONTHESEISMICRESPONSEOFTHEGRAVITYDAMTHEVALUEOFTHEDAMPINGRATIOIS5FORTHESOLUTIONSOFTHISCALCULATIONEXAMPLEFIG3SHOWTHETIMEHISTORYGRAPHSOFTHEVERTICALANDHORIZONTALDISPLACEMENTSOFTHENODALPOINT313LOCATEDATTHECRESTOFTHEDAMWENOTETHATTHEDISPLACEMENTSARERELATIVELYLOWDURINGTHEFIRSTTWOSECONDSBECAUSEOFLOWTHEAMPLITUDESOFTHEEXCITATIONSTHEDISPLACEMENTSREACHTHEIRMAXIMUMAT3SAND47S,20MMWASRECORDEDAT462S,THEMAXIMUMDISPLACEMENTVALUEDOESNOTCORRESPONDTOTHEMAXIMUMAMPLITUDEOFTHEEXCITEMENTTHATISRECORDEDAT365STHENODALDISPLACEMENTSDECREASEAFTER5STHEREISNODAMAGEDURINGMOVEMENTWHERETHEAMPLITUDEISRELATIVELYLOWDAMAGEWEREOBSERVEDSEEFIG4INTHEDAMAFTER353SATTHEINTEGRATIONPOINT815OFTHEELEMENT204,THENFROM354SATINTEGRATIONPOINT811OFTHEELEMENT203ANDFINALLY,FROM354SATGAUSSPOINT807OFTHEELEMENT202ITMAYBENOTEDTHATTHEEVOLUTIONOFTHEDAMAGEISMAINLYCONCENTRATEDINTHETIMEINTERVALWHERETHEMAXIMUMVALUESOFPOSITIVEANDNEGATIVEDISPLACEMENTSOCCURTHERESULTSOBTAINEDINOURSTUDYCOMPAREDTOTHOSEOFREFERENCESCALAYIRANDKARATON2005JIANWENETAL2011ARERELATIVELYSATISFACTORY6CONCLUSIONTHEAIMOFTHISSTUDYIS,FIRSTLYTOHAVETHERESPONSEOFGRAVITYDAMSUBJECTEDTOSEISMICLOADINGASATIMEHISTORYOFTHEDISPLACEMENTS,STRAINSANDSTRESSESSECONDLYTOREPRESENTTHETIMEHISTORYOFDAMAGEEVOLUTIONINTHEINTEGRATIONPOINTSANDTODEDUCETHEAREASLIKELYTOBEDAMAGEDFIRSTLYANDHOWDOTHEYDEVELOPINTHESTRUCTURETHECHOICEOFTHEMETHODISRELATEDTOTHEVERYSMALLTIMEINTERVALCONSIDEREDTHECHOICEOFKOYNACONCRETEGRAVITYDAMISDUETOTHEFACTTHATMANYSTUDIESHAVEBEENDONEONTHISSTRUCTUREWHICHEXPERIENCEDADEVASTATINGEARTHQUAKEIN1967THEEARTHQUAKERECORDSAREINTHEFORMOFACCELEROGRAMSWITHHORIZONTALANDVERTICALCOMPONENTSWHICHAREUSEDASDYNAMICLOADSTHEEFFECTOFTHEDAMPINGCOEFFICIENTISOBVIOUSLYTAKENINTOACCOUNTHYDRODYNAMICEFFECTSARENOTCONSIDEREDINTHISSTUDYREFERENCES1CALAYIR,Y,KARATON,M,2005ACONTINUUMDAMAGECONCRETEMODELFOREARTHQUAKEANALYSISOFCONCRETEGRAVITYDAMRESERVOIRSYSTEMSSOILDYNAMICSANDEARTHQUAKEENGINEERING2520058578692DAVENNEL,RAGUENEAUF,MAZARSJ,IBRAHIMBEGOVICA,2003EFFICIENTAPPROACHESTOFINITEELEMENTANALYSISINEARTHQUAKEENGINEERINGCOMPUTERSANDSTRUCTURES812003122312393DUBJF,1994MODLISATIONSIMPLIFIEETCOMPORTEMENTVISCOENDOMMAGEABLEDESSTRUCTURESENBTON,PHDTHESIS,UNIVERSITYOFPARIS64HENTZS,DAUDEVILLEL,DONZFV,2004DISCRETEELEMENTMODELLINGOFCONCRETESUBMITTEDTODYNAMICLOADINGATHIGHSTRAINRATESCOMPUTERSANDSTRUCTURES822004250925245HERVG,GATUINGTF,2