基于abaqus的高拱坝横缝配筋抗震措施效果分析.docxVIP

基于abaqus的高拱坝横缝配筋抗震措施效果分析 水库是一种有利于地震的水库。目前，该水库已被洪水淹没。但在强震作用下,拱坝中的横缝过度张开可能导致缝间止水破坏,造成库水渗漏和坝体应力的重分布,使拱向应力释放,梁向应力变大,严重影响坝体结构的整体性,甚至危害结构安全。因此,必须采取适当措施控制横缝开度,增强拱坝整体性。常用的措施有横缝插筋和布设阻尼器两种。跨缝钢筋的铺设工作量若较大将会影响灌浆时横缝开度、加大施工难度,影响工期,提高工程造价;而阻尼器的安置相对容易,仅需事先预留安装坑,对工期没有干扰。近年来,国内对高拱坝抗震措施研究已取得一定成果,但就横缝插筋和布设阻尼器分别应用于同一座高拱坝并对比这两种措施加固拱坝效果则较少见。鉴此,本文基于ABAQUS有限元软件,以大岗山高拱坝为例,对这两种措施进行了对比研究,为高拱坝抗震措施的选择提供一些有价值的参考。 1 基材料非线性伸缩缝模型 本文考虑了横缝开合的接触非线性,未考虑坝体、坝基材料的非线性;由于伸缩缝上设有键槽,模型中也不允许缝面间有相对的剪切滑移。在研究横缝动接触的非线性行为及两种抗震措施时,主要有以下3种模型。 1.1 缝合强度对坝体非线性反应的影响 拱坝横缝缝面的相对位移vi和缝应力qi之间的关系为非线性。据文献知,横缝抗拉强度对坝体非线性反应几乎无影响,因此在模拟拱坝横缝力学行为时忽略横缝的抗拉强度,qi与vi满足: 式中,ki为缝闭合时的刚度。 1.2 钢筋等效弹簧模型的建立 横缝配筋模型采用传统钢筋混凝土有限元理论中的整体式模型来考虑钢筋的宏观效果,即将横缝的钢筋面积弥散于钢筋所处的横缝缝面单元的节点上,采用点—点模型,通过在接触点对法向上增加一个与分布钢筋等效的弹簧值来表示钢筋作用。钢筋采用理想弹塑性模型。该等效弹簧值为: Fs=KnsΔL(2) 其中Kns=EsA0/l 式中,Kns为缝面单元上钢筋等效弥散刚度;Es为钢筋的弹模;A0为单元接触面上的钢筋截面面积总和;l为缝面两侧钢筋总的自由段长度;ΔL为横缝开度。 1.3 粘滞流体运动方程 阻尼器的两端分别固定在横缝两侧坝体混凝土单元面的节点上。粘滞阻尼器阻尼力的力学模型为: Fd=sign(v)C|v|ζ(v)C|v|ζ(3) 式中,Fd为阻尼力;sign(v)(v)为v的速度方向,阻尼器两端相对靠近时取正;v为阻尼器两端的相对速度;C为粘滞阻尼系数;ζ为粘滞流体的流动系数。采用法国的坚利特(Jareet)阻尼器,C=8.50×103kN/(m·s-1)ζ,流动系数ζ=0.08。 2 有限元计算分析 大岗山拱坝为双曲抛物线拱坝,坝高210.00 m,坝底高程925.00 m,顶拱中心线弧长622.42 m。坝体混凝土的静弹模为24 GPa,动弹模均取为静弹模的1.3倍,重度为24 kN/m3,泊松比为0.17。考虑地基为非均质弹性介质,按不同岩类选取相应的变形模量值和泊松比νr,见表1。设计地震反应谱采用《水工建筑物抗震设计规范》规定的标准谱,按重现期5 000年、超越概率为0.02确定的设计地震水平向最大峰值加速度为0.557 5 g,竖向取为水平向值的2/3。采用无质量地基,在截断边界三向均匀输入地震。按全部由梁承担的分缝自重来考虑坝体自重,用有限元法以附加质量来考虑动水压力的影响,同时还考虑了泥沙压力和温度荷载。 选取死水位1 120 m时的3种工况。①工况1。未采取任何加固措施。②工况2。在横缝上配置穿横缝钢筋;采用Ⅲ级钢筋,静态屈服强度为400 MPa,弹性模量为200 GPa,直径D为40 mm;据文献选取相对合理的配筋率布置,计算时横缝两侧的钢筋自由段长度共取8 m(每侧4 m);钢筋间距为10D、沿坝厚的上下游1/3部位配筋;钢筋布置位置见图1,图中阴影处表示配筋区。从厚度方向看,坝体有限元网格分为3层,在每条横缝上游侧和下游侧的单元配置钢筋,在中部坝体混凝土单元上不布设钢筋,配筋高度为31 m;上、下游配筋面积均为234.4 m2。③工况3。在全部28条横缝的缝面顶部各布置4台阻尼器,中间21条横缝在高程1 080 m处各布置2台阻尼器,布设方案见图2(图中黑点表示阻尼器)。 3 抗震加固设计 在假定坝体混凝土材料为线弹性、地基为无质量截断地基,并考虑横缝张开非线性条件下,当大岗山拱坝在死水位运行工况未配置钢筋、遭遇规范标准谱设计地震荷载袭击时,其最大横缝开度可达到18.67mm,坝体最大拉、压应力分别为6.66、13.40MPa,表明非线性地震反应显著,因此需采取抗震加固措施。图3、表2、图4分别为3种工况最大横缝开度分布、极值应力及最大横缝开度比较、下游坝面拱和梁向应力等值线。 3.1 端部钢筋变形 考虑地震时横缝反复开合,钢筋残余变形将逐次累积,地震后将阻止横缝闭合,有必要合