加筋土挡土墙动力特性分析

  • 格式:pdf
  • 大小:118.87 KB
  • 文档页数:5

第17卷 第2期2004年4月中 国 公 路 学 报ChinaJournalofHighwayandTransportVol.17 No.2Apr.2004文章编号:100127372(2004)0220028204

收稿日期:2003208203作者简介:李海深(19642),男,湖南宁乡人,湖南大学副教授,工学博士研究生.E2mail:xiangtanlhs

@yahoo.com.cn加筋土挡土墙动力特性分析李海深1,杨果林1,2,邹银生1(1.湖南大学土木工程学院,湖南长沙 410082;2.中南大学土木建筑学院,湖南长沙 410075)摘 要:通过对常用数值计算方法的分析、对比及评价,运用有限元计算方法,建立了加筋土弹塑性本构模型,编制了加筋土挡土墙在动荷载作用下的通用数值计算程序。筋材与其上下表层一定厚度的填土层共同工作,视为筋—土单元,成功地解决了筋材与土介质的接触问题,较好地解决了墙面板与筋材连接处两个接触面单元的特殊问题,不仅使计算方便,而且也符合实际情况。关键词:道路工程;加筋土挡土墙;有限元分析;动力特性;弹塑性本构模型中图分类号:U417.115 文献标识码:AAnalysisofdynamiccharacterofreinforcedearthretainingwallLIHai2shen1,YANGGuo2lin1,2,ZOUYin2sheng1(1.SchoolofCivilEngineering,HunanUniversity,Changsha410082,China;2.SchoolofCivilEngineeringandArchitecture,CentralSouthUniversity,Changsha410075,China)Abstract:Byanalyzing,comparingandappraisingnumericalanalysismethodsincommonused,acom2monanalyticalprogrammehasbeenworkedoutwithelastic2plasticstress2strainrelationshipandfiniteelementmethodisusedtoanalyzereinforcedearthretainingwallunderdynamicload.Asatisfyingso2lutionhasbeenputforwardtothecontactboundary,whichelementofreinforcement2soilwillbeseenatotal.Theespecialproblemofthejointsofcontactfaceelementbetweenplateandreinforcementissolved.Itnotonlycalculatesconvenientbutalsotallieswiththefact.Keywords:roadengineering;reinforcedearthretainingwall;finiteelementanalysis;dynamiccharac2ter;elastic2plasticstress2strainrelationship0引 言求解加筋土挡土结构弹塑性动力问题,实际上可归结为按初始条件和边界条件求解偏微分方程MU・・+CU・t+KUt=Rt的初值—边值问题。当只有初始条件而没有边界条件就成为初值问题,反之则为边值问题。区域内的偏微分方程称为基本方程,初始条件是表示初始状态的条件,边界条件是表示边界约束情况的条件。对于工程上提出的问题能采用解析法按照边值条件求解偏微分方程的仅限于极少数情况。所以,一般只能采用近似方法求解。随着计算机的广泛应用,数值解法逐渐地成为解边值问题的一种有效的方法。数值解法分为区域型和边界型两大类。区域型数值解法主要是有限元法(FiniteElementMethod,FEM)和有限差分法(FiniteDifferenceMethod,FDM)。边界型数值解法主要是边界元法(Bound2aryElementMethod,BEM)。采用差分法时,将所考虑的区域织成网格,用差分近似微分,把差分方程变成微分方程。通过数学上的近似,把求解微分方程的问题变换成求解关于结点未知量的代数方程的问题。采用有限元法时,将所考虑的区域分割成有限小区域,称为有限单元。这些单元仅在有限结点上

© 1994-2008 China Academic Journal Electronic Publishing House. All rights reserved. http://www.cnki.net相连接,根据变分原理把微分方程变换成变分方程,它是通过物理上的近似,把求解微分方程的问题变换成求解关于结点未知量的代数方程的问题。边界元法是继有限差分法、有限元法之后发展起来的又一数值计算方法。采用边界元法求解时,根据积分定理,将区域内微分方程变换成边界上的积分方程,然后,将边界分割成有限大小的边界元,把边界积分方程离散成代数方程,把求解微分方程的问题变换成求解关于结点未知量的代数方程的问题。它具有降一维的特性,所占的计算机时小、计算时间省,与有限元相耦合能较好地解决工程实际问题。由于边界元法本身适用于无限域和半无限域,所以这一方法已在各个工程领域获得了广泛应用。边界单元法只有在边界上剖分单元,通过基本解把域内未知量化为边界未知量来求解,这就使自由度数目大大减少,而且由于基本解本身的奇异性特点,使得边界元法在解决奇异问题时精度较高,另外,基本解可以根据实际问题的特点适当选择,以达到最大限度地节约之功效,甚至可以避免直接处理无限边界问题。再者,边界元法的降维作用使得问题简化许多,特别是三维问题。边界元法特别适用于大区域、无限域和断裂、耦合问题[1]。结合的目的在于利用边界元法处理无限域边界所具有的优点以及有限元法求解非线性问题的灵活性,用以消除有限元法上的“边界效用”及边界元法域内剖分的不便,在减小解题规模方面也有明显的效果,这种方法最先被用于弹塑性分析,以后又发展到用于求解动力问题。文献[2]中,用有限元—边界元耦合法进行动态响应分析,将求解区域分割为二,分别采用有限元法和边界元法,在其交界面上通过迭代法满足界面条件,然后利用WILSON2θ法求解动力微分方程,此法兼有有限元法和边界元法的特点,具有效率高、输入数据少、节省机时等优点。文献[3]中从分区势能和分区混合能原理出发,讨论了弹性、弹塑性力学问题的有限元—边界元对称耦合法,它能够用于常规有限元和杂交有限元等不同的有限元模型和边界元的耦合分析。无限域单元或称无限元是有限元中专门模拟无限域边界的特殊单元,可视为另一种耦合方法。无限元的特点是采用一种特殊的形函数及位移插值函数使其能反映在无限元处的边界条件。无限元应用的主要优点在于:①有效地解决了有限元分析中的“边界效应”及人工边界的缺点,在动力问题中尤为突出;②提高了求解精度及计算效率,对三维问题尤为显著;③显著减小了解题规模,为微机应用提供了十分有利的条件。另外,马立明等人[4],从运动平衡方程的形式入手,应用变分原理,在由空间和时间张成的四维广义空间上,如同建立常规的空间有限元模型一样,对时间域进行离散,建立时间有限元模型。导出了动力分析的单步及两步时间元法。在对两种时间元法的稳定性进行分析的基础上,构造出相应的无条件稳定格式。有限元法、有限差分法和边界元法这些有力的数值技术建立在连续性假设的基础上,离散元法主要处理物体间具有不连续性问题的数值方法,重点是求解多个物体间的接触和冲击问题。一个物体与另一个物体是通过边界接触联系的,边界接触可以随时间变化,一个单元与另一个单元在任何处接触没有限制,可以是结点与结点的接触,也可以是结点与边界接触,接触单元之间产生的接触力遵循不同的接触原则。综上所述,当今有限元法被公认为是一种用数值方法求解工程中遇到的各种问题的最有效、最通用的方法。有限元法目前已成为求解具有已知边界和初始条件或两者之一的偏微分方程组的一种通用的数值解法。鉴于此,笔者用有限元法求解加筋土挡土结构的动力问题。1有限元动力分析方法有限元系统的运动平衡方程Md・・+Cd・+Kd=R(1)式中:M、C、K分别为质量矩阵、阻尼矩阵和刚度矩阵;R是外荷载向量;d、d・和d・分别为有限元集合体的位移、速度和加速度向量。式(1)是通过考虑在时刻t的静力平衡而推导出来的,即式(1)可写成F1(t)+FD(t)+FE(t)=R(t)(2)式中:F1(t)为惯性力,F1(t)=Md・;FD(t)=Cd・为阻尼力;FE(t)为弹性力,FE(t)=Kd。它们均与时

间t有关,因此,在动力分析中,原则上可认为是考虑与加速度有关的惯性力和与速度有关的阻尼力作用在时刻t的静力平衡。在数学上,式(1)是一个二阶线性微分方程组,原则上可用于求解常系数微分方程组的标准过程来求得方程的解。但是,如果矩阵的阶数很高,除非特别利用系数矩阵M、C和K的特殊性质,否则,采用求解一般微分方程组的过程可能要付出很高的代92第2期 李海深,等:加筋土挡土墙动力特性分析

© 1994-2008 China Academic Journal Electronic Publishing House. All rights reserved. http://www.cnki.net价。在实际有限元分析中有两种基本的求解方法:直接积分法和振型叠加法。1.1直接积分法在直接积分法中对式(1)逐步地进行数值积分,进行数值积分前没有把方程进行另一种形式的变换,实质上,直接积分是基于下面的两种想法,第一个想法是只在相隔Δt的一些离散的时间区间上而不是试图在任意时刻t满足式(1)。即包含有惯性力和阻尼力作用的(静力)平衡是在求解区间上的一些离散时刻点上获得的,因此,似乎在静力分析中使用过的所有求解方法,在直接积分法中可有效地使用。第二个想法是假定位移、速度和加速度在每一时间区间Δt内变化,位移、速度和加速度在每一个时间区间内变化,决定着求解的精度、稳定性和求解过程的费用。假设分别用d0、d・0、d・0来表示初始时刻(t=0)的位移、速度和加速度向量为已知且要求出式(1)从t=0到t=T的解,在求解时把时间全程划分为n个相等的时间区间Δt(Δt=T/n),所用的积分格式是在时刻0,Δt,2Δt,3Δt,…,t,t+Δt,…,T上解方程的近似解。由于计算下一个时刻的解的算法要考虑到前面各个时刻的解,因此假定在时刻0,Δt,2Δt,3Δt,…,t的解为已知,来推导出求解时刻t+Δt的解的算法。计算时刻t+Δt的解对于计算自此以后的时刻Δt上的解是具有代表意义的,这样可建立用来计算所有离散时间点上解的一般算法。1.2振型叠加法直接积分法所需的运算次数直接正比于分析中的时间步数,因此,一般来说,当要求较短的时间(n个时间步数)的响应时,可以预料,使用直接积分法是有效的,但是,如果积分必须对许多时间步进行,则先把平衡方程式(1)变换,使之能以较少的费用进行逐步求解就可能会更有效。由于所需的运算次数直接正比刚度矩阵的半带宽mk,因而mk的减小会按比例地降低逐步解题的费用。振型叠加法和直接积分法之间的惟一区别是振型叠加法在时间积分之前进行了基的变换,即从有限元坐标基变换为广义特征向量问题Kφ=ω2Mφ的特征向量基。从数学观点上看,由于n个特征向量所张成的空间和有限元的结点位移所张成的空间是相同的,所以两种分析所得到的解必定是相同的,因而选择直接积分还是选择振型叠加仅需根据效果来考虑。2隐式一显式瞬变动态分析在非线性瞬变动态应力分析中,可用振型叠加法,但实际上,通常都采用时间步进法,这种直接积分法大体上可分为显式法和隐式法两类[5]。采用非常通用且易于执行的显式中心差分法,在计算每一时间步时,由于没有必要进行形式上的矩阵分解,因此只需要较小的计算量。可惜这种方法的稳定是有条件的,而且常常只有当时间步取得很小时才能稳定。采用隐式方法,需要对矩阵进行因式分解,然而可以选择一种无条件稳定的隐式算法,其时间步长只取决于精度的要求。本程序中可任意选择下列算法之一:隐式解法、显式解法、隐式—显式组合解法。在结构有限元分析时,结构中总有一部分区域由于变形大,而通过减小单元尺寸加密单元网格方法来精确模拟其变形和应力分布,其余区域部分则由于变形和应力梯度相对较小而使用较大尺寸单元和粗疏网格剖分形式。在这种情况下,通常一部分区域作隐式积分计算,而后一部分区域作显式积分计算,这样可充分应用显式、隐式算法的各自优点。BELYSTCHKO、HUGHES分别基于结点分离法和单元分离法提出了两种不同的显式—隐式混合算法。在此基础上CHIANG通过设计方程组并行求解器来实现显式—隐式混合并行的功能。冠哲君等人[6]基于区域并行法,提出了一种新的较前者具有更高并行粒度的显式—隐式混合并行法,给出了其算法执行过程。在物理直观意义上解释了该算法的设计思想,其本质是一种单元弱耦合的混合积分算法,并在网络机群上得以实现。隐式、显式时间积分方法结合起来,这时有限元网格包括两组单元:隐式组和显式组,用上标I和E分别表示隐式组和显式组。隐式—显式法中,为满足方程Man+1+PI(dn+1,Vn+1)+PE(󰁪dn+1,󰁴Vn+1)=fn+1式中:fn+1=fIn+1+fEn+1;M=MI+ME,并假定ME为对角线矩阵,就需要在每个时间步进行迭代。本程序参考文献[5]的程序,在这些程序基础上作了补充、修改、改进工作,使之适用于加筋土挡土结构动力分析,兼顾通用性,采用4、8和9结点等参四边形单元处理平面应力、平面应变和轴对称问题,并用总的LAGRANGE公式来处理几何非线性问题,假定材料特性是弹塑性材料的模型。03中 国 公 路 学 报 2004年