用时程分析法对坝进行动力分析.docVIP

PAGE PAGE 1 用时程分析法对大坝进行动力分析 摘要：时程分析法是由结构基本运动方程输入地震加速度记录进行积分，求得整个时间历程内结构地震作用效应的方法。时程分析法为国际通用的动力分析方法，我国《水工建筑物抗震设计规范》（SL—97）首次正式规定该方法为大坝动力分析的方法之一。 关键词：时程分析法大坝动力分析笔者在进行某国际招标项目的设计过程中，根据国际咨询工程师要求，采用时程分析法对大坝进行动力分析。本文对有关内容进行了论述，可供类似工程参考。1工程简介某国际工程以灌溉为主，兼顾防洪发电。工程主要由大坝、电站厂房、分水堰及6.6万公顷的灌区组成。大坝为碾压混凝土曲线重力坝，坝顶高程763.00m，最大坝高133m，坝顶全长231m。自左向右依次布置有左岸非溢流坝段、溢洪道、电站取水口和右岸非溢流坝段。坝体上游面直立，下游坝坡为1：0.6。2工程地质条件坝址区位于峡谷河段，河谷呈“V”字型，底宽约25~40m，两岸基岩裸露，岸坡陡立。河谷两岸岸坡略显不对称，总体上左岸岸坡较陡，右岸岸坡较缓。从河床（高程635.00m）到高程747.00m左右，两岸岸坡陡峻，左岸平均75#176;，右岸平均65#176;；从高程747.00m起向上地形坡度略缓些，但依然较陡，仅局部为45#176;或40#176;左右。坝址区主要由侏罗系灰岩、砂质页岩及第四系全新统冲洪积砂卵砾石组成。河床覆盖层厚度一般为1～5m，多处可见有基岩出露。两岸坝肩岩体较差，断层、裂隙、层间剪切带较发育，其相互组合对坝肩岩体的稳定性有一定的影响。尤其是与坝肩两岸小角度斜交的断层（或裂隙）和裂隙的追踪组合对坝肩稳定最为不利。如断层F2、F3、F13与层理或顺层剪切带组合，把坝肩切割成棱柱体、楔形体，对坝肩稳定十分不利。3计算模型考虑到坝体的规模以及坝址地区的地形地貌和坝基岩体岩性、结构面的分布特征，计算模型区域为：左右岸方向为1000m，上下游方向为800m，坝顶到模型底部为500m。计算单元采用八结点三维实体单元及接触单元，坝体沿坝厚度方向分8层实体混凝土单元，共4760个；围岩划分8807个实体单元；坝体与围岩间接触单元共684个。断层间接触元1316个。整个有限元计算模型单元划分示意见图1，坝体单元划分示意见图2。模型建立充分考虑坝基和坝肩岩体材料的实际分布情况，并模拟了对坝肩岩体稳定和大坝抗滑稳定明显不利的断层F2、F3和F13。同时模型中对大坝基础固结灌浆区域进行了模拟。图1整个计算模型单元图图2坝体网格剖分图坝体结构按弹性计算，围岩结构按弹塑性计算。在进行弹塑性有限元分析时，围岩采用理想弹塑性本构特性。计算中非线性应力屈服准则选用德鲁克-普拉格（Drucker-Prager）准则。αI1+=k式中：α=K=I1=s1+s2+s3J2=1/6[(s1-s2)2+(s2-s3)2+(s3-s1)2]当j0时，在主应力空间中德鲁克-普拉格准则的屈服面是莫尔-库仑内切圆锥；当j=0时，德鲁克-普拉格屈服准则就是米赛斯准则。德鲁克-普拉格屈服准则优点是考虑了s2的作用，适用于混凝土、岩石和土壤等颗粒状材料。4计算参数4.1岩体三维有限元计算模型中模拟了坝基岩体的分区，其分区见图3。坝基、坝肩岩体各区原始的物理力学参数见表1，断层、裂隙物理力学指标见表2。固结灌浆范围内岩体的计算参数见表3。图3坝基岩体分区图表1岩石材料参数表位置左岸坝基右岸分区编号①②③⑦⑧⑨⒀⒁④⑩⒂⑤⑥⑾⑿⒃弹性模量Eo(GPa)0.811.31.31.72.22.631.52.63.90.60.911.52.6泊松比0.34容重(t/m3)2.73表2断层破碎带主要物理力学指标表f′0.4~0.5c′(MPa)0.04弹性模量Eo(GPaF2、F3：1.2~1.8GPa；F13：0.6~2.4Gpa说明：680.00m高程以上取低值，以下（河床以下）取高值。组成物质比重F2、F3、F13：26.8kN/m3表3坝基坝肩岩体主要物理力学指标表部位强度指标坝基660以下660-680680-700700-730730以上左岸右岸左岸右岸左岸右岸左岸右岸左岸右岸岩石的弹性模量E(GPa)6664.43.753.83.753.42.52.62.54.2混凝土坝体碾压混凝土计算参数见表4。表4混凝土材料计算参数坝体碾压混凝土静弹性模量(GPa)19动弹性模量(GPa)28.5容重(kN/m3)24泊松比0.18线胀系数(1/℃)1#215;10-5导温系数(m2/h)0.0044.3坝体的允许应力和稳定安全系数4.3.1坝体的允许应力根据美国工程兵团的拱坝设计规范，坝体应力控制标准见表5。表5坝体应力控制标准工况允许压应力fc(MPa)允许拉应力ft(MPa)上游面下游面上游面下游面极端荷载