《用时程分析法对大坝进行动力分析》由会员分享，可在线阅读，更多相关《用时程分析法对大坝进行动力分析（16页珍藏版）》请在金锄头文库上搜索。
1、用时程分析法对大坝进行动力分析摘要：时程分析法是由结构基本运动方程输入地震加速度记录进行积分，求得整个时间历程内结构地震作用效应的方法。时程分析法为国际通用的动力分析方法，我国水工建筑物抗震设计规范（SL97）首次正式规定该方法为大坝动力分析的方法之一。用时程分析法对大坝进行动力分析摘要:时程分析法是由结构基本运动方程输入地震加速度记录进行积分,求得整个时间历程内结构地震作用效应的方法。时程分析法为国际通用的动力分析方法，我国水工建筑物抗震设计规范（SL-97）首次正式规定该方法为大坝动力分析的方法之一。用时程分析法对大坝进行动力分析摘要：时程分析法是由结构基本运动方程输入地震加速度记录进行积
2、分，求得整个时间历程内结构地震作用效应的方法。时程分析法为国际通用的动力分析方法，我国水工建筑物抗震设计规范（SL-97）首次正式规定该方法为大坝动力分析的方法之一。用时程分析法对大坝进行动力分析摘要：时程分析法是由结构基本运动方程输入地震加速度记录进行积分，求得整个时间历程内结构地震作用效应的方法。时程分析法为国际通用的动力分析方法，我国水工建筑物抗震设计规范（SL97）首次正式规定该方法为大坝动力分析的方法之一。用时程分析法对大坝进行动力分析摘要：时程分析法是由结构基本运动方程输入地震加速度记录进行积分，求得整个时间历程内结构地震作用效应的方法。时程分析法为国际通用的动力分析方法,我国水工
3、建筑物抗震设计规范（SL97）首次正式规定该方法为大坝动力分析的方法之一。用时程分析法对大坝进行动力分析摘要：时程分析法是由结构基本运动方程输入地震加速度记录进行积分，求得整个时间历程内结构地震作用效应的方法.时程分析法为国际通用的动力分析方法，我国水工建筑物抗震设计规范（SL97）首次正式规定该方法为大坝动力分析的方法之一。用时程分析法对大坝进行动力分析摘要:时程分析法是由结构基本运动方程输入地震加速度记录进行积分，求得整个时间历程内结构地震作用效应的方法。时程分析法为国际通用的动力分析方法，我国水工建筑物抗震设计规范（SL-97）首次正式规定该方法为大坝动力分析的方法之一。用时程分析法对大
4、坝进行动力分析摘要：时程分析法是由结构基本运动方程输入地震加速度记录进行积分，求得整个时间历程内结构地震作用效应的方法。时程分析法为国际通用的动力分析方法，我国水工建筑物抗震设计规范（SL97）首次正式规定该方法为大坝动力分析的方法之一.关键词：时程分析法大坝动力分析笔者在进行某国际招标项目的设计过程中，根据国际咨询工程师要求,采用时程分析法对大坝进行动力分析.本文对有关内容进行了论述，可供类似工程参考。1工程简介某国际工程以灌溉为主，兼顾防洪发电。工程主要由大坝、电站厂房、分水堰及6。6万公顷的灌区组成.大坝为碾压混凝土曲线重力坝，坝顶高程763.00m,最大坝高1331n，坝顶全长235。
5、自左向右依次布置有左岸非溢流坝段、溢洪道、电站取水口和右岸非溢流坝段。坝体上游面直立，下游坝坡为l:0o6o2工程地质条件坝址区位于峡谷河段，河谷呈“V”字型，底宽约2540m,两岸基岩裸露，岸坡陡立。河谷两岸岸坡略显不对称,总体上左岸岸坡较陡，右岸岸坡较缓。从河床（高程635。00m）到高程747。00m左右，两岸岸坡陡峻，左岸平均75°；,右岸平均65°：从高程747.00m起向上地形坡度略缓些，但依然较陡，仅局部为45°或40°；左右.坝址区主要山侏罗系灰岩、砂质页岩及第四系全新统冲洪积砂卵砾石组成。河床覆盖层厚度一般为15m,多处可见有基岩出露.两岸
6、坝肩岩体较差，断层、裂隙、层间剪切带较发育，其相互组合对坝肩岩体的稳定性有一定的影响.尤其是与坝肩两岸小角度斜交的断层（或裂隙）和裂隙的追踪组合对坝肩稳定最为不利。如断层F2、F3、F13与层理或顺层剪切带组合，把坝肩切割成棱柱体、楔形体，对坝肩稳定十分不利。3计算模型考虑到坝体的规模以及坝址地区的地形地貌和坝基岩体岩性、结构面的分布特征，计算模型区域为：左右岸方向为1000m,上下游方向为800nl，坝顶到模型底部为500m。计算单元采用八结点三维实体单元及接触单元，坝体沿坝原度方向分8层实体混凝土单元,共4760个；围岩划分8807个实体单元;坝体与围岩间接触单元共684个.断层间接触元1
7、316个.整个有限元计算模型单元划分示意见图1,坝体单元划分示意见图2.模型建立充分考虑坝基和坝肩岩体材料的实际分布情况，并模拟了对坝肩岩体稳定和大坝抗滑稳定明显不利的断层F2、F3和F13。同时模型中对大坝基础固结灌浆区域进行了模拟。图1整个计算模型单元图图2坝体网格剖分图坝体结构按弹性计算，围岩结构按弹塑性计算。在进行弹塑性有限元分析时，围岩采用理想弹塑性本构特性。计算中非线性应力屈服准则选用德鲁克-普拉格(DruckerPrager)准则.aIl+=k式中：a=K=H=sl+s2+s3J2=l/6(sl-s2)2+(s2s3)2+(s3si)2当jO时，在主应力空间中德鲁克-普拉格准则的
8、屈服面是莫尔一库仑内切圆锥；当F0时.，德鲁克一普拉格屈服准则就是米赛斯准则。德鲁克-普拉格屈服准则优点是考虑了s2的作用，适用于混凝土、岩石和土壤等颗粒状材料.4计算参数4.1岩体三维有限元计算模型中模拟了坝基岩体的分区，其分区见图3.坝基、坝肩岩体各区原始的物理力学参数见表1,断层、裂隙物理力学指标见表2。固结灌浆范围内岩体的计算参数见表3o图3坝基岩体分区图表1岩石材料参数表位置左岸坝基右岸分区编号(13)(11)(弹性模量Eo(GPa)0.811.31.31.72.22.631。52.63。90。60。911.52。6泊松比0.34容重(t/m3)2.73表2断层破碎带主要物理力学指标
9、表f0。40。5c(MPa)0。04弹性模量Eo(GPaF2、F3：1.21.8GPa;F13：0。62.4Gpa说明:680.00m高程以上取低值,以下(河床以下)取高值。组成物质比重F2、F3、F13：26.8kN/m3表3坝基坝肩岩体主要物理力学指标表部位强度指标坝基660以下660-680680700700-730730以上左岸右.岸左岸右岸左岸右岸左岸右岸左岸右岸岩石的弹性模量E(GPa)6664。43.753.83。753.42.52。62.54.2混凝土坝体碾压混凝土计算参数见表4。表4混凝土材料计算参数坝体碾压混凝土静弹性模量(GPa)19动弹性模量(GPa)28。5容重(kN
10、/m3)24泊松比0.18线胀系数(1/)1×105导温系数(m2/h)0。0044.3坝体的允许应力和稳定安全系数33。1坝体的允许应力根据美国工程兵团的拱坝设计规范，坝体应力控制标准见表5。表5坝体应力控制标准工况允许压应力fc(MPa)允许拉应力ft(MPa)上游面下游面上游面下游面极端荷载组合(MCE)13o316o72o082。64.3.2稳定安全系数根据美国工程兵团的拱坝设计规范，坝体、坝肩抗滑稳定安全系数见表6.表6坝体、坝肩抗滑稳定安全系数匚况坝体、坝肩极端荷载组合(MCE)1。15荷载及荷载组合计算工况为水库正常运用时遭遇最大可信地震(MCE),上游水位750.40
11、m,下游无水,作用在大坝上的荷载主要有自重、静水压力、地震动水压力、波浪压力、泥沙压力、地震动沙压力、扬压力以及温度荷载。特殊荷载的计算方法如下：1)温度荷载参照拱坝荷载组合的规定，在YCE计算时，温度荷载为温升荷载。以坝体封拱时的温度场为温度荷载计算的基准温度场，以10年后6月15日的坝体温度场为计算温升荷载的计算温度场，两个温度场的差值即为温升荷载。2)地震动沙压力水工建筑物抗震设计规范(SL97)规定，地震对泥沙压力的影响，一般不予考虑，但本工程坝前淤沙原度很大，动沙压力对坝体应力、变形和稳定有较大影响，必须予以考虑。由于目前对于动沙压力的作用机理仍处于理论研究阶段，还没有成熟的计算方法
12、.在美国小坝设计一书中提出“近似的假定为饱和泥沙的动力影响相当于水”，同时考虑到地震过程中淤沙的液化，故在动泥沙压力计算时参照动水压力的计算公式计算地震动沙压力。3)扬压力坝基部分帷幕后主排水孔处扬压力折减系数a=0.3：坝肩部分帷幕后主排水孔处扬压力折减系数a=0.4。4)地震荷载(MCE)最大可信地震(MCE)最大水平地震加速度为0。25g,其地震加速度时程曲线见图4.在计算时，考虑两个方向的地震同时作用。顺河流方向水平地震加速度为0。25g,铅直方向为0.125g。图4地震加速度时程曲线6计算结果及分析在坝体上选取若干具有代表性的特征点，分析各特征点的应力、位移和稳定时程曲线，找出坝体的
13、最不利时刻，然后对最不利时刻的坝体状态进行分析。为了使特征点具有代表性，在特征点的选取过程中，考虑了坝体上下游面、坝体内部、左右两岸和不同高程等因素，按高程选取了三个断面，分别为高程6301n(9个特征点)，高程7031n(9个特征点)和高程763(6个特征点)，共计24个特征点.6。1坝体应力、位移和坝基稳定安全系数时程曲线的分析6.1.1坝体应力坝体的应力状况是通过质点的应力时程曲线反映出来的，在此以拱冠梁上游面高程630nl处特征点D为例进行说明。特征点应力时程曲线见图57。图5特征点D第一主应力时程曲线(Pa)图6特征点D第三主应力时程曲线(Pa)图7特征点D正应力时程曲线(Pa)通过
14、图57,可以得到特征点D的第一主应力的最大值出现在8秒时刻,其值为0o6Mpa,最小值出现在7o8和8。2秒时刻，其值为一0.24Mpa；第三主应力的最小值出现在5.2秒时刻，其值为-1。25Mpa,最大值出现在8。2秒时刻，其值为一1.98Mpa；正应力的最小值出现在8秒时刻，其值为04Mpa,最大值出现在8。2秒时刻,其值为-1。4Mpa。通过对各特征点应力的综合分析可得：1）建基面（高程630m）在6.2s和8。0s时，坝体应力状况较差，第一主应力最大值0.7MPa,正应力最小值-0.04MPa,均发生在8s时刻。第三主应力最小值为一L25MPa,发生在8o2s时刻。2）高程703m高程
15、坝内应力状况,在8s左右坝体应力状况比较差8.0s应力状况最差。第一主应力力最大值为1。16MPa和正应力最小值为0.3MPa,为拉应力，均发生在8s时刻.3）坝顶（高程763m）处应力状况，在6。3s和8。4s间应力条件较差。第一主应力最大值为0。2MPa,发生在8s时刻。4）通过分析可以看出坝体在地震作用下，各特征点应力状况最差的时刻,大多出现在8s左右时刻，其中第一主应力最大值出现在8s时刻的几率最高,所以8s为应力的最不利时刻。6.1.2坝体位移计算结果与分析按照应力计算结果的统计方法对各特征点的位移计算结果进行统计分析可得：1）坝基（高程630m）处各特征点的顺河向最大位移均发生在8
16、.2秒时刻,其最大值-24。7mm;高程703m处各特征点的顺河向最大位移大部分发生在7。9秒时刻，其最大值一38mm；坝顶（高程763m）处各特征点的顺河向最大位移发生在7。9或8。0秒时刻，其最大值46mm。最大位移均发生在拱冠梁位置。2）铅直向位移最大值均发生8。2秒时刻，其最大值一46。5mm,坝体的铅直位移最小值几乎全发生在8s时，其最小值为14o7mmo3）由以上分析可得坝体在7。9、8.0和8。2秒坝体状态较差，而且在&0秒时坝体的铅直位移最小，对坝体的稳定最为不利。6olo3建基面抗滑稳定利用有限元分析的结果，分别取出在设计荷载作用下的坝体和地基内各点的应力，并利用抗剪点安全系数计算公式，计算各个点的安全系数，绘制坝基个特征点的抗剪稳定安全系数时程曲线，然后对坝体建基面的抗剪稳定状况进行评价。抗剪点安全系数计算公式如下：其中：T该点剪应力。一一该点正应力，
