p波倾斜入射对河谷两岸地震动影响分析.docxVIP

p波倾斜入射对河谷两岸地震动影响分析 地震发生时，地震波以地震波的形式传播到地面。地震波作为传播路径和地形条件的一部分，随着传播路径的不同，地表各点的振动不可避免。然而，在当前的抗强试验中，只有欧洲标准考虑了地震波的空间变化，而其他标准采用了相同的地震波输入。也就是说，假设结构与周围接触面上的不同振动完全相同。这一假设适用于小型基础设施，但对于大体积结构的情况，如果采用这一假设，显然是非常粗糙的。根据中国网络智能1的最新数据，地震波沿地面的水平传播速度有限，波形有空间变化。虽然在50米范围内，地震波也存在明显差异。 地震动的非一致输入体现为地震动的时空变化, 主要有行波效应和局部场地效应.行波效应是引起地震动时空变化的重要原因, 从物理机制来讲, 地震波的斜入射是引起行波效应的主要因素.当震源离场地较近时, 地震波并非垂直向上入射, 是以一定角度倾斜向上传播.研究表明, 当震源距离场地较近时, 基岩场地地震波入射角平均在30°左右.杜修力等对地震波斜入射条件下地下结构的地震反应进行分析发现:地震波倾斜入射时, 地下结构的动力反应与地震波垂直入射时有较为明显的差异, 地震波斜入射条件下结构的动力反应要大于地震波垂直入射时的动力反应.另外, 地下结构的地震反应对入射波作用方向十分敏感, 地震波入射方向变化时, 地下结构的变形和各部位应力均会发生明显变化. 局部场地条件对地震动有很强的放大与缩小效应, 受不规则地形影响, 地震波在河谷地形中传播与平坦场地相比有很大的不同.河谷地震动分布在小尺度空间上并不一致, 呈现明显的差异性.台湾翡翠拱坝坝址的实测地震记录表明:河谷不同高程的地震动幅值有很大差异, 而且河谷两岸同一高程处的地震动加速度也存在很大差异. 河谷场地是实际工程中常见的重要场地, 两岸的地震动效应与河谷形状、地震波入射角度有关.本文主要研究不同形状河谷在地震波倾斜入射条件下的动力响应情况, 并采用黏弹性人工边界结合等效荷载的方法, 实现地震波斜入射的有限元输入. 1 斜入射时地震波的有限输入 将地震波动问题转化为波源问题, 采用黏弹性边界模拟吸收散射波场, 通过在边界节点施加等效荷载的方式来实现斜入射地震波的波动输入. 1.1 人工边界模拟 采用有限单元法进行动力分析时, 通常从半无限空间连续介质中截取有限域进行分析计算, 这会导致传播的地震动在边界节点处发生反射, 人为夸大计算区域的动力响应.因此, 在模拟半无限地基时, 必须设置合适的人工边界条件来模拟无限域介质带来的影响.本文选用黏弹性边界来模拟, 即在人工边界位置设置由线性弹簧与黏滞阻尼器并联的弹簧-阻尼元件, 确保计算区域产生的散射波穿过人工边界时不发生反射. 二维人工边界的弹簧刚度系数和黏滞阻尼系数按下式计算: 法向: 切向: 1.2 入射波场位移场 设置黏弹性人工边界后, 由于黏弹性边界只能模拟吸收外行散射波, 对于外源激振问题, 入射地震波需要以一定的方式输入计算区域内.本文采用刘晶波等提出的等效荷载波动输入方法, 把输入地震动转化为等效荷载来模拟波动入射过程. 边界节点的等效荷载计算公式如下所示: 式中:Fl(t) 为边界节点l处的等效荷载;Al为人工边界面上l节点的影响面积;σ0(x, y, t) 为连续介质中自由场产生的应力;K、C分别为黏弹性人工边界的法向、切向弹簧刚度和黏滞阻尼系数;u (x, y, t) 、6) u (x, y, t) 分别为入射波场位移和速度. 由波动理论可知, 地震波在自由表面处会发生波形转换, 即当P波倾斜入射到自由表面时, 会衍生出另外的P波和SV波.假定平面P波从左侧人工边界以α角度入射 (如图1所示) , 位移时程为up(t) .考虑地震波在传播过程中的时间延迟, 左侧人工边界处的内行场由入射角为α的入射P波up(t-Δt1) 、反射角为α的反射P波A1up(t-Δt2) 和反射角为β的反射S波A2up(t-Δt3) 构成;底部人工边界的内行场只有入射角为α的入射P波up(t-Δt4) , 内行位移场可表示如下. 左侧人工边界: 底部人工边界: 根据平面波传播时的应力状态和应力状态变换公式, 可以得到平面P波以α角度从左侧边界入射时人边界节点l处的应力计算公式. 左侧人工边界: 底部人工边界: 式中, Δt1、Δt2、Δt3、Δt4分别为左侧人工边界直接入射P波、地表反射P波和S波、底边界直接入射P波相对于入射波初始时刻波振面的延迟时间, 可以通过波的传播距离除以相关波速获得.相应的速度场可以通过对内行位移场差分或求导计算获得;反射P波和S波的幅值和反射角可根据弹性半空间表面反射规律确定. 1.3 入射地震波时程计算模型 以二维半空间问题验证上述斜入射方法的模拟精度, 取介质密度2 400kg/m3, 弹性模量10G