外文翻译-基于离散单元法的砌石重力坝安全分析工具

外文翻译题目亭口水利枢纽斜墙坝设计及工程量计算专业水利水电工程班级工101学生指导教师2014年基于离散单元法的砌石重力坝安全分析工具摘要介绍一种基于离散单元法的砌石重力坝分析数值模型。大坝和岩基用3到4节点的基本块组成的块集合表示。复杂的块形状通过把基本块整合到宏模块来得到，允许模型应用在从等效连续到完全不连续分析的各种情形。开发了一个接触面公式，能够根据基本块之间建立的接触面，基于一种精确的边边方法表示宏模块之间的相互作用。描述了模型的主要数值方面，特别介绍了接触面的创建和更新步骤以及一个支持一种高效且能够得到明显结果的算法的数值设备。讨论了一个现有的砌石坝的安全评价的应用，包括结构的应力分析和涉及大坝岩石界面附近不同路径的滑动失效机制的评估。关键词砌石坝；离散元；应力分析；失效机制1简介结构分析必须使用适当的手段来实现它的最终目的。这些手段必须能够（1）模拟建筑物的几何和物理特征，尤其是不连续的和共有的特征；（2）用一套完整的方法来模拟荷载，能够考虑所涉及的不同现象间的相互影响；（3）评价非线性作用，特别是能够界定失效结构。砌石重力坝应该被理解为一个包含坝体本身，水库，岩石基础的系统。坝体和岩石是不均匀且不连续的介质。坝体和围岩的交界面也是不连续的，需要特别注意。不连续面控制着砌石坝的强度，因为它们是薄弱面，决定着主要失效机制。另外，大坝所受的各种不同的荷载需要一套完整的方法来处理，因为它们之间相互关联。这些独特的特点使得大多数的数值工具，无论商业的还是科学的，都不能完全适合地模拟浆砌石重力坝。在这种背景下，新的数值工具的发展就显得尤为重要。本文将描述一种为砌石重力坝的静力、动力及流体力学的分析量身定做的数值工具，离散单元法。离散单元法最初是为了处理岩石力学问题而提出的替代有限元法的方法。离散单元法的基本原理是将不连续介质看作受力作用不同的块的集合，因此不同于将其看作同一单元的有限元法。这些数值处理方法也被广泛地应用于砖石结构。在CUNDALL所开发的产品的基础上发展起来的二维编码通用离散元程序（UDEC），已经被用于包括混凝土坝的地基等方面的研究中，主要用于尝试通过岩石评估失效机理。TAONE等人做了一个在坝体与岩石交界面上的滑移离散单元法分析，充分地考虑了交界面的不规则几何形状及应力集中。离散单元法的代码通常通过将应变块离散成三角形均匀应变单元的内部网格来代表应变块。指令“离散有限单元法”通常被用于编码块单元的允许破坏量来模拟渐进破坏过程。本文所介绍的模型是基于离散单元法，满足三个设计要求，在一个完全由作者自己开发的新型软件工具里实施的。首先，本模型要将浆砌石坝和岩石基础用一套完整的方法做成块系统的部件；其次，软件工具需要提供一种可实现的方法来用同一网格表示等效的连续介质和块状模型；最后，该工具需要包含大坝工程分析的所有要素，比如水流量及接头处和加强部位的压力，比如主动或被动锚固，以及应用在静力分析和地震分析中涉及到的荷载的方法。所有这些部分通过一个兼容的数据结构相互作用。因此，本模型在一个更普遍的框架下结合了标准离散单元法的能力，该框架可以满足实际应用的兼顾刚性块和应变块，连续网格和离散单元的要求。此外，还采用了一个基于边界间相互作用的非传统接触面公式，能够得出一个更精确的交界面压力。新的系统与离散单元法一样有按照预先设定好的路径模拟一个连续体碎裂成块的能力，但是采用了一种基于联合刚度和适合砌石和岩石构成规律的表示接触面的方法，这点不同于后者，例如MUNJIZA的接触力公式。有关应力计算的方面将在下面的部分详细讨论，本文将介绍一个砌石坝的安全评价在实际中的应用的例子。大坝与岩石基础的水工分析以及最新开发的分析工具，将在另一篇论文中被介绍。图1砌石坝和岩石基础形成的不连续介质图2宏模块构成的连续和不连续模型2模型离散化和接触面本数值工具旨在模拟如图1所示的浆砌石坝与岩石基础所组成的系统。根据设计惯例以及大坝安全规范，结构的保守假设是二维假设，其中的实际原因是以往大坝都设计成重力坝，无法保证其拱效应，而且这样的计算模型更容易理解。结构离散化的基本要素是具有三或四个边界的块，可以是刚性的或可变形的，并且可以计划性地在同一模型中使用。块的选择应由结构，特征以及分析的目的来决定。在性能方面，由于刚性块只在单元的质心处建立运动方程，减少了模型自由度，所以刚性块的计算速度更快。刚性块的计算优势是只在显动态分析时有实质性作用，因为其静态的结果通常可以迅速得到。在大坝工程中，结构应力分析和地基应力分析通常是必须的，所以应变块是首选。每个应变块在这里都被假定为一个具有完全高斯积分的等参非线性有限元。一般形状的块可以通过将3和4节点的块整合成宏模块来创建。这是模拟如砌石坝和岩石基础这类不连续介质的重要特征。这样等效连续表示整个或部分系统，其中每个块正是有限元网格的一个单元。一个宏模块就是一个块的排列，相应块的顶点相互重合，组成一张连续的网络。在同一宏模块的不同块之间没有相对位移，所以没有接触压力。宏模块类似于有限元网格，但是有一个明确的结果，因为整体刚度矩阵的组合并不会发生。图2A展示了一个有点个单独模块组成的不连续模型。图2B表示的是一个类似但是连续的模型，它是由一个宏模块构成的。另一个与有限元网格相似的连续模型在图2C展示。图2D表示一个由两个宏模块构成的混合模型，在这两个宏模块之间有一个明显而重要的联接。宏模块有一个包含一系列块和一系列由一个主节点和多个从节点构成的宏节点的数据结构。宏节点与独立节点有着相同的自由度，所有数值操作只能在主节点上进行。在计算循环中，所有从节点的力必须集中在主节点上，在新坐标上计算完成后，从节点从各自的主节点中更新数据。尽管步骤较多，利用宏模块还是有减少接触面和自由读数目的优点。同一个模型可能有几个宏模块，且每个宏模块可以有不同材质的块。ADEMBASEDTOOLFORTHESAFETYANALYSISOFMASONRYGRAVITYDAMSABSTRCTANUMERICALMODELFORANALYSISOFMASONRYGRAVITYDAMSBASEDONTHEDISCRETEELEMENTMETHODISPRESENTEDTHEDAMANDTHEROCKFOUNDATIONAREREPRESENTEDASBLOCKASSEMBLIES,USINGELEMENTARY3AND4NODEBLOCKSCOMPLEXBLOCKSHAPESAREOBTAINEDBYASSEMBLINGTHEELEMENTARYBLOCKSINTOMACROBLOCKS,ALLOWINGTHEMODELTOBEAPPLIEDINVARIOUSSITUATIONSRANGINGFROMEQUIVALENTCONTINUUMTOFULLYDISCONTINUUMANALYSISACONTACTFORMULATIONWASDEVELOPED,WHICHREPRESENTSTHEINTERACTIONBETWEENMACROBLOCKSINTERMSOFCONTACTSESTABLISHEDBETWEENELEMENTARYBLOCKS,BASEDONANACCURATEEDGEEDGEAPPROACHTHEMAINNUMERICALASPECTSOFTHEMODELAREDESCRIBED,ADDRESSINGINPARTICULARTHECONTACTCREATIONANDUPDATEPROCEDURES,ANDTHENUMERICALDEVICESTHATSUPPORTANEFFICIENTEXPLICITSOLUTIONALGORITHMANAPPLICATIONTOTHESAFETYEVALUATIONOFANEXISTINGMASONRYDAMISDISCUSSED,INCLUDINGSTRESSANALYSISINTHESTRUCTURE,ANDTHEASSESSMENTOFSLIDINGFAILUREMECHANISMS,INVOLVINGDIFFERENTPATHSINTHEVICINITYOFTHEDAMROCKINTERFACE1INTRODUCTIONSTRUCTURALANALYSISMUSTUSEAPPROPRIATEMETHODSTOACHIEVEITSFINALPURPOSESTHESEMETHODSSHOULDBECAPABLEOFIMODELINGTHEGEOMETRICALANDPHYSICALCHARACTERISTICSOFTHESTRUCTURE,INPARTICULARTHEDISCONTINUITIESANDJOINTS,IIMODELLINGTHELOADSINANINTEGRATEDMANNER,TAKINGINTOACCOUNTTHEINTERACTIONBETWEENTHERELEVANTPHENOMENAINVOLVED,ANDIIIEVALUATINGTHENONLINEARBEHAVIOUR,PARTICULARLYALLOWINGTHEDEFINITIONOFFAILUREMECHANISMSMASONRYGRAVITYDAMSSHOULDBEUNDERSTOODASASYSTEMCOMPOSEDOFTHEDAMITSELF,THERESERVOIR,ANDTHEROCKMASSFOUNDATIONTHEDAMANDTHEROCKMASSAREHETEROGENEOUSANDDISCONTINUOUSMEDIATHEDAMROCKINTERFACEISALSOADISCONTINUITYWHICHREQUIRESPARTICULARATTENTIONTHEDISCONTINUITYSURFACESCONTROLTHEBEHAVIOUROFMASONRYDAMS,BECAUSETHEYAREWEAKNESSPLANESTHATDETERMINETHEMAINMECHANISMSOFFAILUREINADDITION,DAMSARESUBJECTTOAWIDEVARIETYOFLOADSREQUIRINGANINTEGRATEDAPPROACHSINCETHEYAREOFTENCORRELATEDTHESEPARTICULARFEATURESMAKETHEMAJORITYOFTHEAVAILABLENUMERICALTOOLS,BOTHCOMMERCIALANDSCIENTIFIC,NOTENTIRELYSUITABLEFORMODELINGMASONRYGRAVITYDAMSINTHISCONTEXT,THEDEVELOPMENTOFNEWANALYSISTOOLSISREQUIREDHERE,ATAILOREDNUMERICALIMPLEMENTATIONOFTHEDISCRETEELEMENTMETHODDEMFORSTATIC,DYNAMICANDHYDROMECHANICALANALYSISOFMASONRYGRAVITYDAMSISDESCRIBEDTHEDISCRETEELEMENTMETHODWASINITIALLYPROPOSEDASANALTERNATIVETOTHEFINITEELEMENTMETHODFEMTOADDRESSROCKMECHANICSPROBLEMS1DEMWASBASEDONTHEREPRESENTATIONOFTHEDISCONTINUOUSMEDIAASANASSEMBLYOFBLOCKSINMECHANICALINTERACTION,THUSDIFFERINGFROMTHESTANDARDFEMAPPROACHBASEDONJOINTELEMENTS2,3THESENUMERICALAPPROACHESHAVEALSOBEENWIDELYAPPLIEDTOMASONRYSTRUCTURESEG4THE2DCODEUDEC5,WHICHEVOLVEDFROMCUNDALLSPIONEERINGWORK,HASBEENUSEDINSEVERALSTUDIESINVOLVINGCONCRETEDAMFOUNDATIONS,MOSTLYINTENDINGTOASSESSFAILUREMECHANISMSTHROUGHTHEROCKMASSEG68TATONEETAL9PERFORMEDADEMANALYSISOFSLIDINGONTHEDAMROCKINTERFACE,CONSIDERINGADETAILEDREPRESENTATIONOFTHEIRREGULARGEOMETRYOFTHISSURFACEANDTHEENSUINGSTRESSCONCENTRATIONSDISCRETEELEMENTMETHODCODESUSUALLYREPRESENTDEFORMABLEBLOCKSBYDISCRETIZINGTHEMINTOANINTERNALMESHOFTRIANGULARUNIFORMSTRAINELEMENTSEG,5THEDESIGNATIONDISCRETEFINITEELEMENTMETHOD10,17ISOFTENAPPLIEDTOCODESTHATALLOWTHEBREAKAGEOFTHEBLOCKELEMENTSTOSIMULATEPROGRESSIVEFAILUREPROCESSESTHEMODELPRESENTEDINTHISPAPERISBASEDONDEMANDWASDEVISEDWITHTHREEMAINREQUIREMENTS,IMPLEMENTEDINANOVELSOFTWARETOOLFULLYDEVELOPEDBYTHEAUTHORSFIRSTLY,ITISINTENDEDTOMODELINANINTEGRATEDMANNERBOTHTHEMASONRYDAMANDTHEROCKFOUNDATIONASCOMPONENTSOFABLOCKYSYSTEMSECONDLY,THESOFTWARETOOLSHOULDPROVIDEAPRACTICALMEANSTOADDRESSBOTHEQUIVALENTCONTINUUMANDBLOCKYMODELS,USINGTHESAMEMESHFINALLY,THETOOLNEEDSTOINCLUDEALLTHEFEATURESREQUIREDINDAMENGINEERINGANALYSIS,SUCHASWATERFLOWANDPRESSURESINTHEJOINTS,REINFORCEMENTELEMENTS,SUCHASPASSIVEORACTIVEANCHORS,ANDTHEMEANSTOAPPLYTHELOADSINVOLVEDINSTATICANDSEISMICANALYSISALLTHESECOMPONENTSINTERACTTHROUGHACOMPATIBLEDATASTRUCTURETHEREFORE,THEPRESENTMODELCOMBINESTHESTANDARDDEMCAPABILITIESINAMOREGENERALFRAMEWORK,WHICHALLOWSCOMBININGRIGIDANDDEFORMABLEBLOCKS,CONTINUUMMESHESANDDISCRETECOMPONENTS,ASREQUIREDBYTHEAPPLICATIONMOREOVER,ANONTRADITIONALCONTACTFORMULATIONISADOPTED,BASEDONEDGEEDGEINTERACTION,WHICHPROVIDESAMOREACCURATESTRESSREPRESENTATIONINTHEINTERFACETHENEWCODESHARESWITHDEMTHECAPABILITYTOSIMULATEFRACTURINGOFACONTINUUMINTOBLOCKSTHROUGHPREDEFINEDPATHS,BUTADOPTSAREPRESENTATIONOFCONTACTBASEDONTHEJOINTSTIFFNESSESANDCONSTITUTIVELAWSAPPROPRIATEFORMASONRYANDROCK,NOTFOLLOWING,FOREXAMPLE,MUNJIZASFORMULATIONOFCONTACTFORCEPOTENTIALS17THEASPECTSRELATEDTOTHEMECHANICALCALCULATIONWILLBEDISCUSSEDINDETAILINTHEFOLLOWINGSECTIONS,ANDANEXAMPLEOFAPPLICATIONTOTHESAFETYASSESSMENTOFAMASONRYDAMINOPERATIONWILLBEPRESENTEDTHEHYDRAULICANALYSISOFTHEDAMANDROCKFOUNDATION,ALSOINCORPORATEDINTHENEWLYDEVELOPEDANALYSISTOOL,WASDESCRIBEDINADIFFERENTPAPER112MODELDISCRETIZATIONANDCONTACTSTHENUMERICALTOOLISINTENDEDTOMODELSYSTEMSCOMPOSEDOFAMASONRYDAMANDITSROCKFOUNDATION,ASSHOWNSCHEMATICALLYINFIG1TWODIMENSIONALANALYSISISCONSERVATIVELYASSUMEDFORTHESESTRUCTURES,FOLLOWINGCOMMONDESIGNPRACTICESANDDAMSAFETYCODESEG1214,FORPRACTICALREASONSHISTORICALLY,MASONRYDAMSWEREDESIGNEDASGRAVITYDAMSARCHINGEFFECTSCANNOTBEGUARANTEEDAND,THECOMPUTATIONALMODELISSIMPLERTOUNDERSTANDTHEFUNDAMENTALELEMENTOFDISCRETIZATIONOFTHESTRUCTUREISTHEBLOCKWITHTHREEORFOUREDGES,WHICHMAYBERIGIDORDEFORMABLE,ANDCANBEUSEDSIMULTANEOUSLYINTHESAMEMODELTHESTRUCTURE,CHARACTERISTICSANDOBJECTIVESOFTHEANALYSISSHOULDDICTATETHECHOICEOFBLOCKSINTERMSOFPERFORMANCE,THECALCULATIONISFASTERFORTHERIGIDBLOCKSBECAUSETHEEQUATIONOFMOTIONISESTABLISHEDONLYINTHECENTROIDOFTHEELEMENT,THUSREDUCINGTHEDEGREESOFFREEDOMOFTHEMODELTHECOMPUTATIONALADVANTAGEOFRIGIDBLOCKSISONLYRELEVANTINEXPLICITDYNAMICANALYSIS,SINCESTATICSOLUTIONSAREUSUALLYVERYFASTTOOBTAININDAMENGINEERING,STRESSANALYSISINTHESTRUCTUREANDFOUNDATIONISUSUALLYREQUIRED,SODEFORMABLEBLOCKSAREPREFERREDINCASEOFDEFORMABLEBLOCKS,EACHBLOCKISASSUMEDHEREASANISOPARAMETRICLINEARFINITEELEMENTWITHFULLGAUSSINTEGRATIONBLOCKSOFGENERALSHAPESMAYBECREATEDBYASSEMBLINGTHE3AND4NODEBLOCKSINTOMACROBLOCKSTHISISANIMPORTANTFEATURETOMODELDISCONTINUOUSMEDIA,SUCHASMASONRYDAMSANDROCKMASSFOUNDATIONSINTHISWAY,ITISPOSSIBLETOADOPTANEQUIVALENTCONTINUUMREPRESENTATIONOFTHEWHOLESYSTEM,ORPARTOFTHESYSTEM,INWHICHEACHBLOCKISJUSTANELEMENTOFTHEFEMMESHAMACROBLOCKISACOMBINATIONOFBLOCKS,FORMINGACONTINUOUSMESH,INWHICHTHEVERTICESARECOINCIDENTBETWEENTHEBLOCKSOFTHESAMEMACROBLOCK,RELATIVEMOVEMENTISNOTPERMITTED,SOTHEREARENOCONTACTFORCESTHEMACROBLOCKISSIMILARTOAFINITEELEMENTFEMESHBUTWITHANEXPLICITSOLUTIONBECAUSETHEASSEMBLAGEOFAGLOBALSTIFFNESSMATRIXDOESNOTTAKEPLACEFIG2ASHOWSADISCONTINUOUSMODELCOMPOSEDBYTWOINDIVIDUALBLOCKSASIMILARMODEL,BUTCONTINUOUS,COMPOSEDBYONEMACROBLOCK,ISPRESENTEDINFIG2BANOTHERCONTINUUMMODEL,SIMILARTOAFEMESHISSHOWEDINFIG2CAHYBRIDMODEL,COMPOSEDBYTWOMACROBLOCKS,ISPRESENTEDINFIG2D,WITHANEXPLICITJOINTCONSIDEREDBETWEENTHETWOMACROBLOCKSTHEMACROBLOCKHASADATASTRUCTURECONTAININGALISTOFBLOCKSANDALISTOFMACRONODES,WITHAMASTERNODEANDSEVERALSLAVENODESTHEMACRONODEHASTHESAMEDEGREESOFFREEDOMOFANYINDIVIDUALNODE,ANDALLNUMERICALOPERATIONSCANFOCUSONLYONTHEMASTERNODEDURINGTHECALCULATIONCYCLE,ALLFORCESFROMTHESLAVENODESMUSTBECONCENTRATEDINTHEMASTERNODE,ANDAFTERCALCULATIONOFNEWCOORDINATES,SLAVENODESAREUPDATEDFROMTHERESPECTIVEMASTERNODEDESPITETHESEPROCEDURES,THEUSEOFMACROBLOCKSHASTHEADVANTAGEOFREDUCINGTHENUMBEROFCONTACTSANDTHENUMBEROFDEGREESOFFREEDOMTHESAMEMODELMAYHAVESEVERALMACROBLOCKSANDEACHMACROBLOCKCANHAVEBLOCKSWITHDIFFERENTMATERIALSREFERENCES1CUNDALLPAACOMPUTERMODELFORSIMULATINGPROGRESSIVELARGESCALEMOVEMENTSINBLOCKYROCKSYSTEMSINROCKFRACTUREISRM,NANCY19712GOODMANRE,TAYLORRL,BREKKETLAMODELFORTHEMECHANICSOFJOINTEDROCKJSOILMECHFOUNDDIVASCE1968943637593WITTKEWROCKMECHANICSTHEORYANDAPPLICATIONSWITHCASEHISTORIESBERLINSPRINGERVERLAG19904LOURENPBCOMPUTATIONSOFHISTORICALMASONRYCONSTRUCTIONSPROGSTRUCTENGMATER20024301195ITASCA,UNIVERSALDISTINCTELEMENTCODEUDECVERSION50,MINNEAPOLIS20116LEMOSJVDISCRETEELEMENTANALYSISOFDAMFOUNDATIONSINSHARMAVM,SAXENAKR,WOODSRD,EDITORSDISTINCTELEMENTMODELLINGINGEOMECHANICS,BALKEMA,ROTTERDAM1999P891157BARLAG,BONINIM,CAMMARATAGSTRESSANDSEEPAGEANALYSESFORAGRAVITYDAMONAJOINTEDGRANITICROCKMASSIN1STINTERNATIONALUDEC/3DECSYMPOSIUM,BOCHUM2004P26388GIMENESE,FERNNDEZGHYDROMECHANICALANALYSISOFFLOWBEHAVIORINCONCRETEGRAVITYDAMFOUNDATIONSCANGEOTECHJ20064