地面大位移地下管线计算分析.pdfVIP

地面大位移时地下管线计算分析 徐凤萍 李鹏程 杨广智 侯忠良 (冶金部建筑研究总院抗震所 北京100088 提【要】论述了地下管线在由于地震土体液化产生地面大位移的情况下变形的计算方法。在弹塑性 理论分析的基础上，采用非线性有限元法编制了埋地管线在地面大位移作用下内力和变形计算的专 用程序.分析了地下管线的内力和变形的影响因素，推导了简化计算公式 规【川竺桩液化地面大位移地下， 弯曲变形 ·止奎‘ 一、引言 地震造成土体液化产生变形和位移，使埋在地下的管线产生应变和位移，会导致 管道的严重破坏。本文模拟大位移情况下管一土相互间的作用力，编制了非线性有限 元程序 ，对液化产生地面大位移的情况下埋地管线的的内力和变形进行 了计算和分 析. 二、埋地管胡的非线性分折 管埋于液化土时，土和管材具有非线性特性，管线的变形以弯曲变形为主，管的 抗弯刚度为El，抗剪刚度为GAla (其中参数a为考虑了横截面翘曲的修正系数). 有限元方程为:[K](u)=仃}.假定管线为一根长L、两端固定的Timoshenko梁，符 合如下假设: 1.当弯矩达到屈服弯矩时，管的整个横截面立刻变成塑性的; 2.变形以前垂直于中性轴的直线变形以后仍为直线，但却不一定再垂直于中性 轴; 3.任何一点(x,y)的横向位移w就等于中性轴上的横向位移; 4.管在xy平面内承载; 5.采用小挠度、小变形理论; 6.剪切力在单元 内不变，且剪应力与剪应变总是弹性关系. 在每个单元中，弯曲力矩均为常数:剪切力在单元内按线性变化，管达到极限弯 矩MO以前，在弯曲刚度El下作弹性变形，当达到极限弯矩从 时，弹塑性有限元分析 表明，整个单元就进人塑性状态，此时，弯曲刚度El就用弹塑性抗弯刚度(EI)eP来代 691 替，而剪切刚度GA/a则假定保持不变.当继续加大荷载时，假定材料具有线性应变 强化的特性，切线弯曲刚度用(EI):来表示. 材料开始屈服后，若进一步施加荷载，必然导致弯曲力矩的增加，同时伴随有曲 率的改变，假设将曲率分成弹性和塑性两部分，则有 def=(def)g+(def)fl;屈服 EI 后弯曲刚度EI可用EI(1一 )来代替m{ E1+IT rwl 在有限网格的任一点上，横向位移和转角的表达式为:I月=[N1{w} Loll 在整个有限网格内，任一点的曲率和剪应变可表示为: dB _ dco 、 _ 一，厂=B,P(; 气厂一a=ttsrp, 改尤 “ 以龙 本文采用的验算准则为:强度准则a【]=口:和延性准则，延性验算容许应变取 [e]=10%o. 三、程序介绍及计算步膝 本程序采用模块法编制，即用独立的子程序来实现有限单元法的各种主要运算. 有限元的计算过程均采用有限元离散化得到HO十f=0，式中:0为基本未知量矢量; f为外加的荷载矢量;H为组集的刚度矩阵。对于非线性问题，系数矩阵H依赖于未 知数沪或其导数.计算中采用了广义Newton-Raphson法(即切线刚度法)，在荷载增量 的范围内把问题线性化.由于M 为非线性函数，只能对矢量P作近似的估计，故令 P-.f等于残余值W(O)，计算时把YW(O)尽量减小. 求解过程采用切线刚度法.为了进一步减少对内存的要求，同时也大大减少方程 的消元和回代阶段所需的算术运算次数，对总刚度矩阵采用了一维存储的方法.程序 输人的主要参数包括结构的已知条件和各单元的性质、材料的物理性质、地面位移型 式、计算容许误差等. 各单元的加载情况:在液化区，土一管相互作用关系为弹一塑性曲线