Page 93 - 2023年第54卷第9期
P. 93

如图 4所示,在三维场景中除规定地震波入射与 z轴正方向的夹角为 θ (入射角),还引入了传播
              平面与 x轴的夹角 α (方位角)。















                                                  图 4 三维场地中地震波的传播
                  假设地震波从左下角点起振,设任意节点 M(x,y,z)在图 4(a)中 r方向上距离起振点为 r,当 P
              波和 SV波入射时,位移表达方式如下:
                                             u(x,y,z,t) =u(z,t - r?c)cos α
                                            { u(x,y,z,t) =u(z,t - r?c)sin α                            (11)
                                              x
                                                             0r
                                                                     ij r
                                              y
                                                                     ij r
                                                             0r
                                             u(x,y,z,t) =u(z,t - r?c)
                                                                     ij r
                                                             0z
                                              z
                  当出平面 SH波入射时,振动方向垂直于图 4(a)中的灰色平面,位移表达方式如下:
                                            u(x,y,z,t) =u(z,t - r?c)sin α
                                           { u(x,y,z,t) =- u(z,t - r?c)cos α                           (12)
                                                             0r
                                                                     ij r
                                             x
                                                                      ij r
                                                              0r
                                             y
                                            u(x,y,z,t) =u(z,t - r?c)
                                             z               0z     ij r
              式中 c为地震波沿 r方向传播的水平视波速。
                    r
                  图 4(b)为底面上任意点在 r方向投影后的新坐标,经过几何计算后,投影点 M 的坐标计算如下:
                                                                                          1
                                                                   2
                                              { x= (x + ytan α )?(1 + tanα )  2                        (13)
                                                1
                                               y= tan α ·(x + ytan α )?(1 + tanα )
                                                1
                  当距离 y轴的 r方向长度不是 Δ r的整倍数时,需对其进行线性插值,参考式(10)即可。
              3 算例验证及方法讨论
              3.1 自由波场位移验证 图 5为三维的成层场地有限域模型;地基材料参数如表 1所示。假设地震波
                                                                                                          =
              起振点为地基左下角,入射脉冲波的位移时程如图 6所示。分别分析垂直入射 P波、斜入射 P波( θ P
                                                                       = 60°)、垂直入射 SH波和斜入射 SH波
              45° ,α P  = 40°)、垂直入射 SV波、斜入射 SV波( θ SV        = 30°,α SV
                            = 20° )情况下,地基表面中点的位移响应情况。图 7分别为 6类情况下成层场地表面中
              ( θ SH  = 50° ,α SH















                            图 5 三维场地模型及观测点                               图 6 输入地震动位移时程


                                                                                                   1
                                                                                              —   1 0 3 —
   88   89   90   91   92   93   94   95   96   97   98