重力式码头等价粘弹性动力反应分析.pdfVIP

重力式码头等价粘弹性动力反应分析 吴敏敏，张建民，张 嘎 (清华大学水利水电工程系水沙科学与水利水电工程国家重点实验室，北京，100084) 摘 要:本文采用等价粘弹性模型模拟地基土在地震中的应力应变规律，并应用商业有限元程序Mac‘及其子程序 HYPELA二次开发对重力式码头的地震过程进行了数值棋拟.首先验证了有限元分析中单元尺寸和截断边界设里距离 对计算结果的影响，在此基础上参照神户港沉箱重力式码头典型断面建立模型，进行了地震反应分析，并来用一残余 变形公式对震后码头的变形规律进行预测 ‘同时，针对不同级别码头分别建立有限元模型进行动力反应分析，对现行 《规范》中码头顶部相对底部加速度放大系数提出了修改建议. 关键词:重力式码头;有限元法;等价粘弹性;动力反应分析 1 前言 随着我国社会经济的迅猛发展以及对外贸易活动的日趋频繁，水运吞吐量快速增长，港口码头建 设的规模越来越大，其中许多大型码头结构已经超出了现行设计规范的适用范围，缺少相应的设计理 论及经验积累。确保这些大型码头结构的安全性、特别是强震区的抗震安全性问题己经成为需要认真 研究的热点和重点课题。重力式码头在我国港口工程中应用广泛，但其动力设计仍采用拟静力方法， 随着我国港口建设的飞速发展和码头规模的不断扩大，一方面根据传统的拟静力法进行抗震设计存在 诸多的不足之处，另一方面数值分析方法及计算机技术的发展，为码头动力有限元分析提供了强大的 支撑。本文采用等价粘弹性模型对重力式码头进行了动力反应分析，探讨了码头在地震过程中的加速 度响应及变形规律。 2 等价粘弹性模型 本文采用Hardin和Dmevieh提出的等价粘弹性模型1I]描述土的动力特性。该模型对实际滞回圈用 倾角和面积相等的椭圆代替，并由此确定粘弹性体的两个基本参数:剪切模量G和阻尼比几。假定滞 回圈顶点轨迹为双曲线，模型主要计算公式如下: _ . (。生、 k,Y} J‘x‘=anal- I 几 = 月峪，~ — (1) lP, G_ 1十k.7, 一 1+k,Y} 公式中，P,为大气压，0为有效球应力，参数k、k2和凡曰可通过共振柱或动三轴的模量试验确定: Y}为参考剪应变，由下式确定: l :=}Ya)}}/f;〔}z (2) 公式中，(Yde为有效剪应变，且 (Yd),=0.65(Ya)_ (3) 523 第，届全国土动力学学术会议论文集 北京 2006-11-3-6 公式中，(Y)二为该时段的最大剪应变· 由于粘弹性模型不能模拟动力作用条件下土的残余变形，因此，本文采用沈珠江院士提出的在等 价粘弹性范围内考虑残余变形的经验公式，残余体积应变△#V.R和残余剪切应变△YR分别如下: △、，:一。Or,)b)`-P一(、合 △、一。(Y(d)g),SI总 (4) 公式中，ClC2C3C;和CS为模型的输入参数，可由常规动三轴液化试验测定;又为剪应力动用 水平，AN为该时段的等效振动周次，N,为累计的等效振动周次。 基于上述模型，本文基于Marc有限元软件进行了二次开发，建立了码头的动力分析数值方法。 3 单元尺寸和边界范围的验证分析 3.，单元尺寸验证分析 单元尺寸的大小对动力计算精度有较大影响。如果单元网格尺寸过大，振动的高频部分无法通过。 一般单元尺寸越小，计算精度越高，但由此会导致自由度数增加，使得计算量增大，非线性计算的收 敛变得更为困难。为此，本文通过如下算例对合理的