基于全连接神经网络的地铁车站响应分析与地震强度指标优选

钟紫蓝; 倪博; 史跃波; 张成明; 申家旭; 杜修力

doi:10.11779/CJGE20221448

基于全连接神经网络的地铁车站响应分析与地震强度指标优选 English Version

北京工业大学城市与工程安全减灾教育部重点实验室, 北京 100124

基金项目:

国家自然科学基金面上项目 51978020

详细信息

作者简介:
钟紫蓝(1986—)，男，博士，副研究员，主要从事地下结构抗震研究。E-mail: zilanzhong@bjut.edu.cn

通讯作者:
申家旭, E-mail: shenjx@bjut.edu.cn

中图分类号: TU443
计量
- 文章访问数: 362
- HTML全文浏览量: 60
- PDF下载量: 95
出版历程
- 收稿日期: 2022-11-22
- 网络出版日期: 2024-03-14
- 刊出日期: 2024-02-29

Response analysis of subway station and optimization of seismic intensity measures based on fully connected neural network

Key Laboratory of Urban Security and Disaster Engineering of Ministry of Education, Beijing University of Technology, Beijing 100124, China

摘要

摘要: 为了降低随机地震响应分析的计算成本，将人工神经网络方法用于构建概率地震需求模型（PSDM），以预测地铁车站结构的地震响应，并对适用于地铁车站结构响应预测的地震强度指标（IM）进行了优选。首先选取了200条实测地震动，计算IM，并对典型的三层三跨地铁车站结构进行有限元建模，将IM与最大层间位移角作为输入与输出训练全连接神经网络模型（FCNN），得到了最大层间位移角的预测模型。最后基于训练后FCNN输入层到隐含层中的权重矩阵与传统方法对IM进行优选，得出了对最大层间位移角影响最大的IM。研究结果表明：训练后FCNN能以0.95的精度预测地铁车站最大层间位移角，且计算耗时仅为数值模拟的1/5；针对矩形地下结构最大层间位移角，速度型和速度反应谱型指标的影响明显高于其他类型指标，其中速度谱强度（VSI）对最大层间位移角的影响最大。
- 全连接神经网络 /
- 概率地震需求模型 /
- 地铁车站结构 /
- 地震强度指标
Abstract: In order to reduce the calculation cost of random seismic response analysis, the artificial neural network method is used to build a probabilistic seismic demand model (PSDM) to predict the seismic response of subway station structures, and the seismic intensity measure (IM) suitable for the prediction of subway station structural response is optimized. First, 200 measured ground motions are selected, IM is calculated, and the typical three-story and three-span subway station structure is modeled by the finite element method. Then, the IM and the maximum layer drift are used as the input and output to train the fully connected neural network (FCNN), and the prediction model for the maximum layer drift is obtained. Finally, the IM is optimized based on the weight matrix from the FCNN input layer to the hidden layer after training and the traditional methods, and the IM that has the greatest impact on the maximum layer drift is obtained. The results show that the FCNN after training can predict the maximum layer drift of subway station with an accuracy of 0.95, and the calculation efficiency is 18000 times higher than that of numerical simulation. For the maximum layer drift, the impact of velocity type and velocity response spectrum type indices is significantly higher than other types of indices, among which the velocity spectrum intensity (VSI) has the largest impact on the maximum layer drift.
- FCNN /
- probabilistic seismic demand model /
- subway station structure /
- seismic intensity measure

HTML全文

0. 引言

层状地基–薄板结构是一种工程上广泛应用的力学模型^[1]。随着高速公路、铁路、大运量航运的快速发展,其附属基础土建设施在高速荷载作用下的响应规律已与传统静力分析结论明显不符,因此移动荷载作用下该模型的动力响应已成为领域研究热点之一^[2]。早期受限于数学理论和电算水平,学者多在平面应变、各向同性、半无限空间地基的假设下研究该问题。但岩土体作为一种天然沉积的层状各向异性材料,其横观各向同性（transversely isotropic简称TI）及三维空间特性会对动力响应造成一定影响^[3-5]。近年来,数值计算水平的提升使得相应半解析研究逐渐成为可能,部分学者已在求解方法以及相关工程应用上获得了一定的启发性结论^[6-8]。

TI半无限体的Green函数最早由Rajapakse等^[9]通过引入3个势函数给出,其求解了内部荷载作用下的柱系响应结果,研究表明位移响应受激励频率、材料参数影响大,而应力则不受具体参数影响。对于层状地基,Barros等^[10]利用TRM得到层状弹性体作用移动线载后的稳态解答。薛松涛等^[11]通过研究斜SH波入射成层TI地基的传播特征建立直角坐标下的动力刚度矩阵,并指出横观各向同性可显著改变空间波动特性。艾智勇等^[12]从弹性动力学方程解法入手导出解析层元,进而获得直角坐标下的层状TI地基空间动力刚度矩阵。韩泽军等^[13]通过引入对偶变量获得精细积分形式的层状TI地基刚度矩阵,但存在数值计算程序复杂等具体问题。

Eskandari-Ghadi等^[14]证明可用2个位移势函数来完备表达TI弹性体不解耦的位移变量,当实际位移场旋度与材料对称轴垂直时,仅用一个势函数即可。在此基础上,Rahimain等^[15]研究了轴对称波动规律。Khojasteh等^[16]导出半无限空间、轴对称双层地基的Green函数。巴振宁等^[17]利用Ghadi解推导了TI层状地基的Green函数。但上述势函数需利用无穷级数构成的方程组确定多项式系数,求解过程繁琐且物理含义不强。

为简化分析,一般考虑简单荷载、层间完全连接下的结构响应规律。Achenbach^[18]研究了移动点载作用下半平面–薄板结构动力响应,指出振动水平受板、地基相对刚度、接触面粗糙性影响大。房营光^[19]解析了半平面–无限大板响应规律,计算表明板体接触应力随速度增加向地基转移。王博等^[20]的数值分析表明地基不同向弹性模量对板体变形影响较大。此外,王春玲等^[21]和Muho^[22]利用无穷级数逼近分别获得TI半空间上有限、无限大板的响应结果,但级数形式复杂且与荷载类型有关。

综上可见,现有研究多集中于半无限空间–薄板结构,而层状地基–薄板结构对应研究则较少,并且理论基于势函数分解等数学方法,程序复杂、物理含义弱。本文研究TI层状地基–薄板结构在典型移动荷载下的动力响应,意在寻求一种推导简洁、求解方便的算法,并在此基础上分析具体参数改变对动力响应的影响规律,为相关工程结构设计、施工提供参考。

1. 直角坐标系下的振动方程

简化抽象的层状地基–薄板体系如图1所示分为板体、层状体、半无限空间三部分。板体满足薄板假定,层状体及半无限空间为TI均质材料。层状体水平叠放,层间相互连接且几何、材料特性不同,顶面作用匀速移动矩形简谐荷载。

图 1 TI层状地基-薄板体系示意图

Figure 1. Transversely isotropic layered foundation-plate system

下载: 全尺寸图片幻灯片

为简化求解过程,在移动坐标系下建立力学模型,其建立原则、与固定坐标系间变换关系及模型参数在1.1节、1.2节、2.1节中详细说明。

1.1 坐标系统变换

参考坐标系建立在层状地基顶面,竖向坐标正向朝下且满足右手法则。当体系受到高速移动荷载作用时,基于X-Y-Z-O固定直角坐标系分析将使后续求解产生较大计算成本。在稳态振动的前提下,基于以速度 $c (\to)$ 运动的x-y-z-o移动直角坐标系分析将简化计算流程。

根据时间t内荷载移动位移,两坐标系间坐标变换关系为 $\vec{R} - \vec{c} t = \vec{r}$ ,若荷载仅沿x轴正向运动（略去 $c (\to)$ 箭头）,由位移场F_i变化规律与坐标系选择无关易得： $F_{X} (X, Y, Z, t) =$ $F_{x} (x, y, z, t) =$ $F_{x} (X - c t, y, z, t)$ 。则两系下位移分量的微分关系为

$\begin{array}{l} \frac{\partial^{n} F_{R}}{\partial R^{n}} = \frac{\partial^{n} F_{r}}{\partial r^{n}}, \\ \frac{\partial^{n} F_{X}}{\partial t^{n}} = {(\frac{\partial F_{x}}{\partial x} - c \frac{\partial F_{x}}{\partial t})}^{(n)}, \end{array}}$

(1)

式中,F_i为位移场以下标区分对应坐标系,R为固定坐标系中任意点坐标分量如X,r为移动坐标系中任意点坐标分量如x。

因稳态响应变量可用荷载圆频率 $ω$ 表示,故各函数有其分离时间变量形式即F(x,y,z,t)= F(x,y,z)e^iωt,则式（1）中时间导数项变为

$\begin{array}{l} \frac{\partial F_{X}}{\partial t} = [i ω F_{x} (x, y, z) - c \frac{\partial F_{x} (x, y, z)}{\partial x}] e^{i ω t}, \\ \frac{\partial^{2} F_{X}}{\partial t^{2}} = [c^{2} \frac{\partial^{2} F_{x} (x, y, z)}{\partial x^{2}} - 2 c ω \frac{\partial F_{x} (x, y, z)}{\partial x} i - \\ ω^{2} F_{x} (x, y, z) \begin{array}{l} \end{array}] e^{i ω t} 。 \end{array}}$

(2)

移动坐标系下求解结果可通过 $\vec{R} = \vec{r} + c t$ 映射到常用固定坐标系对应坐标点处。

1.2 模型基本方程及变量解耦

固定坐标系下各位移变量表示的横观各向同性体物理方程组为

$(\begin{matrix} σ_{11} \\ σ_{22} \\ σ_{33} \\ σ_{13} \\ σ_{23} \\ σ_{12} \end{matrix}) = (\begin{matrix} C_{11} & C_{12} & C_{13} & 0 & 0 & 0 \\ C_{11} & C_{13} & 0 & 0 & 0 \\ C_{33} & 0 & 0 & 0 \\ C_{44} & 0 & 0 \\ Sym . & C_{44} & 0 \\ C_{66} \end{matrix}) (\begin{matrix} \frac{\partial u}{\partial X} \\ \frac{\partial v}{\partial Y} \\ \frac{\partial w}{\partial Z} \\ \frac{\partial w}{\partial X} + \frac{\partial u}{\partial Z} \\ \frac{\partial w}{\partial Y} + \frac{\partial v}{\partial Z} \\ \frac{\partial u}{\partial Y} + \frac{\partial v}{\partial X} \end{matrix})$ 。

(3)

式中 $u, v, w$ 为任意点在固定坐标系下的位移分量（省略下标R）；X,Y,Z为固定系坐标； $σ_{11}$ - $σ_{12}$ 等表示应力分量；中部方阵代表弹性矩阵,各弹性常数需满足正定性要求^[23]：

$\begin{array}{l} C_{11}, C_{33}, C_{44}, C_{12} > 0, \\ C_{11} C_{33} - C_{13}^{2} - C_{33} C_{66} > 0 。 \end{array}}$

(4)

为引入质量密度ρ及时间变量t,将式（3）代入平衡方程得到忽略体力的Lame-Navier运动方程组：

(5)

薄板结构采用Kirchhoff薄板模型,其在固定坐标系下的动力控制方程为

$D_{q} \nabla^{4} γ (X, Y, t) + ρ_{b} h_{b} \frac{\partial^{2} γ (X, Y, t)}{\partial t^{2}} = Q (X, Y, t) - P (X, Y, t),$

(6)

式中, $D_{q} = E_{p} h_{b}^{3} / [12 (1 - μ_{p}^{2})]$ 为薄板抗弯刚度,E_p,μ_p,ρ_b,h_b分别为薄板弹性模量、泊松比、质量密度、板厚, $γ$ 为待求板体挠度,Q为外部竖向荷载,P为待求板间竖向接触力。

为便于求解,式（5）各方程需解耦,解耦原理与波的传播特性有关,引入下列变换^[24]：

$\begin{array}{l} [\begin{matrix} u^{'} \\ v^{'} \end{matrix}] = T [\begin{matrix} u \\ v \end{matrix}], \\ [\begin{matrix} σ_{c} \\ σ_{d} \end{matrix}] = T [\begin{matrix} σ_{13} \\ σ_{23} \end{matrix}], \\ T = [\begin{matrix} \frac{\partial}{\partial x} & \frac{\partial}{\partial y} \\ \frac{\partial}{\partial y} & \frac{- \partial}{\partial x} \end{matrix}], \end{array}}$

(7)

式中,T为变换矩阵, $u^{'}, v^{'}, σ_{c}, σ_{d}$ 分别为变换后的位移及应力分量。

联立式（1）,（3）,（7）可得解耦后的部分物理方程：

$\begin{array}{l} σ_{c} = C_{44} \nabla^{2} w + C_{44} \frac{\partial u^{'}}{\partial z}, \\ σ_{d} = C_{44} \frac{\partial v^{'}}{\partial z}, \\ σ_{z} = σ_{33} = C_{13} u^{'} + C_{33} \frac{\partial w}{\partial z} 。 \end{array}}$

(8)

联立式（1）,（2）,（5）,（7）得到解耦后的运动方程：

(9)

实际应用中,多采用工程常数描述模型材料特性,其与弹性常数转换按下式计算：

$\begin{array}{l} E_{h} = - \frac{8 C_{13}^{2} (C_{11} - C_{12}) (C_{11} C_{33} - C_{13}^{2})}{{(2 C_{13}^{2} - C_{11} C_{33} + C_{12} C_{33})}^{2}}, \\ E_{v} = \frac{2 (C_{11} C_{33} - C_{13}^{2})}{C_{11} + C_{12}}, \\ μ_{vh} = \frac{2 C_{13}^{2} - C_{11} C_{33} + C_{12} C_{33}}{2 C_{11} C_{13} + 2 C_{12} C_{13}}, \\ μ_{h} = \frac{2 C_{13}^{2} + C_{11} C_{33} - C_{12} C_{33}}{2 C_{13}^{2} - C_{11} C_{13} + C_{12} C_{13}}, \\ G_{v} = C_{44}, \end{array}}$

(10)

式中, $E_{v}$ , $E_{h}$ 为竖、水平向弹性模量, $μ_{h}$ , $μ_{vh}$ 为平行、垂直材料TI面的泊松比, $G_{v}$ 为垂直TI面的剪切模量。

对板体动力方程式（6）按移动系简谐变量化简得

$\begin{matrix} ρ_{b} h_{b} (\frac{\partial^{2} γ (x, y)}{\partial x^{2}} - 2 c ω \frac{\partial γ (x, y)}{\partial x} i - ω^{2} γ (x, y)) \\ = Q (x, y) - P (x, y) - D_{q} \nabla^{4} γ (x, y) \end{matrix}$ 。

(11)

1.3 控制方程求解

上述控制方程较复杂,本文通过双重傅里叶变换将偏微分方程组转化为常微分方程组,从而获得频率域下的位移解答。若有三元函数f(x,y,z),其双重傅里叶正逆变换规定如下：

$\bar{\bar{f}} (k_{x}, k_{y}, z) = \int_{- \infty}^{+ \infty} \int_{- \infty}^{+ \infty} f (x, y, z) e^{- i (k_{x} x + k_{y} y)} d x d y$ ,

(12)

$f (x, y, z) = \frac{1}{4 π^{2}} \int_{- \infty}^{+ \infty} \int_{- \infty}^{+ \infty} \bar{\bar{f}} (k_{x}, k_{y}, z) e^{i (k_{x} x + k_{y} y)} d k_{x} d k_{y}$ 。

(13)

式中,f为变换域下与f对应的函数变量,k_x,k_y为变换域下的坐标变量即波数,定义 $k^{2} = k_{x}^{2} + k_{y}^{2}$ 。

联立式（8）,（9）,（12）得常微分形式对应方程：

$\begin{array}{l} {\bar{\bar{σ}}}_{c} = - C_{44} k^{2} \bar{\bar{w}} + C_{44} \frac{d \bar{\bar{u'}}}{d z}, \\ {\bar{\bar{σ}}}_{d} = C_{44} \frac{d \bar{\bar{v'}}}{d z}, \\ {\bar{\bar{σ}}}_{z} = C_{13} \bar{\bar{u}}' + C_{33} \frac{d \bar{\bar{w}}}{d z}, \end{array}}$

(14)

$\begin{array}{l} - C_{11} k^{2} \bar{\bar{u}}' - (C_{13} + C_{44}) k^{2} \frac{d \bar{\bar{w}}}{d z} + C_{44} \frac{d^{2} \bar{\bar{u'}}}{d z^{2}} + \\ ρ {(ω - c k_{x})}^{2} \bar{\bar{u}}' = 0, \\ - C_{44} k^{2} \bar{\bar{w}} + C_{33} \frac{d^{2} \bar{\bar{w}}}{d z^{2}} + (C_{13} + C_{44}) \frac{d \bar{\bar{u'}}}{d z} + \\ ρ {(ω - c k_{x})}^{2} \bar{\bar{w}} = 0, \\ - (\frac{C_{11} - C_{12}}{2}) k^{2} \bar{\bar{v}}' + C_{44} \frac{d^{2} \bar{\bar{v}}'}{d z^{2}} + ρ {(ω - c k_{x})}^{2} \bar{\bar{v}}' = 0 。 \end{array}}$

(15)

为便求解,引入广义位移变量 $\bar{\bar{U}}', \bar{\bar{V}}'$ 及其组成的状态变量 ${\bar{\bar{W}}}_{i}$ 将位移变量统一：

$\begin{array}{l} \bar{\bar{U}}' (k_{x}, k_{y}, z) = {(\bar{\bar{u}}', k^{2} \bar{\bar{w}})}^{T}, \\ \bar{\bar{V}}' (k_{x}, k_{y}, z) = \bar{\bar{v}}', \end{array}}$

(16)

$\begin{array}{l} {\bar{\bar{W}}}_{1} (k_{x}, k_{y}, z) = {(\bar{\bar{U}}', \frac{d \bar{\bar{U}}'}{d z})}^{T}, \\ {\bar{\bar{W}}}_{2} (k_{x}, k_{y}, z) = {(\bar{\bar{V}}', \frac{d \bar{\bar{V}}'}{d z})}^{T} 。 \end{array}}$

(17)

联立式（15）,（17）可得两组矩阵微分方程（A,B系数矩阵中元素的具体形式见附录）：

$\begin{array}{l} {\frac{d {\bar{\bar{W}}}_{1} (k_{x}, k_{y}, z)}{d z}}_{4 \times 1} = A {(k_{x}, k_{y})}_{4 \times 4} {\bar{\bar{W}}}_{1} {(k_{x}, k_{y}, z)}_{4 \times 1}, \\ {\frac{d {\bar{\bar{W}}}_{2} (k_{x}, k_{y}, z)}{d z}}_{2 \times 1} = B {(k_{x}, k_{y})}_{2 \times 2} {\bar{\bar{W}}}_{2} {(k_{x}, k_{y}, z)}_{2 \times 1} 。 \end{array}}$

(18)

根据矩阵分析理论,上述两组矩阵方程存在exp(A(k_x,k_y)z)C(k_x,k_y)形式的一般解满足要求,令 $\bar{\bar{W_{1}}} (k_{x}, k_{y}, z) = \exp (A z) \bar{\bar{W_{1}}} (k_{x}, k_{y}, 0)$ ,根据汉密尔顿-凯莱定理, ${\bar{\bar{W}}}_{1} (k_{x}, k_{y}, z)$ 可以传递矩阵T₁的形式表达：

$\begin{matrix} \bar{\bar{W_{1}}} {(k_{x}, k_{y}, z)}_{4 \times 1} = \sum_{i = 0}^{3} b_{i} (z) A^{i}_{4 \times 4} {\bar{\bar{W}}}_{1} {(k, 0)}_{4 \times 1} \\ = T_{1} {(k, z)}_{4 \times 4} \bar{\bar{W_{1}}} {(k, 0)}_{4 \times 1} \end{matrix}$ 。

(19)

具体先计算矩阵A的特征多项式：

$\det (λ I - A) = \prod_{i = 1}^{4} (λ - λ_{i}) = 0$ 。

(20)

其解经计算：

$\begin{array}{l} λ_{1, 2}^{2} = a + b, \\ λ_{3, 4}^{2} = a - b 。 \end{array}}$

(21)

式中,a,b见附录。以4个特征值互异为例,可列得其含参多项式 $r (λ_{i}, z)$ ：

$r (λ_{i}, z) = \exp (λ_{i} z) = \sum_{j = 1}^{4} b_{j} (z) λ_{i}^{j}$ 。

(22)

联立含参多项式解得b_j,进而由A得T₁表达式：

$r (A, z) = \exp (A z) = \sum_{i = 1}^{4} b_{i} (z) A^{i} = T_{1} {(k, z)}_{4 \times 4}$ 。

(23)

则可求得状态变量 ${\bar{\bar{W}}}_{1}$ 并分割传递矩阵T₁：

$\bar{\bar{W_{1}}} {(k, z)}_{4 \times 1} = T_{1} (k, z) \bar{\bar{W_{1}}} (k, 0) = [\begin{matrix} T_{11} & T_{12} \\ T_{13} & T_{14} \end{matrix}] \bar{\bar{W_{1}}} {(k, 0)}_{4 \times 1}$ 。

(24)

进而得到 $\bar{\bar{U}}' (k, z)$ 与顶面状态变量 ${\bar{\bar{W}}}_{1} (k, 0)$ 关系：

$\bar{\bar{U}}' (k, z) = (\begin{matrix} T_{11} & T_{12} \end{matrix}) {\bar{\bar{W}}}_{1} (k, 0)$ 。

(25)

同理可推得

$\bar{\bar{V}}' (k, z) = (\begin{matrix} T_{11} & T_{12} \end{matrix}) {\bar{\bar{W}}}_{2} (k, 0)$ 。

(26)

为推导刚度矩阵,定义广义应力变量 ${\bar{\bar{V}}}_{1} 、 {\bar{\bar{V}}}_{2}$ 如下。联立式（14）,（16）求取广义应力与状态变量间关系：

${\bar{\bar{V}}}_{1} (k, z) = {[\bar{\bar{σ_{c}}}, \bar{\bar{σ_{z}}}]}^{T}, {\bar{\bar{V}}}_{2} (k, z) = [\bar{\bar{σ_{d}}}]$ ,

(27)

$\begin{array}{l} {\bar{\bar{V}}}_{1} {(k, z)}_{2 \times 1} = L_{1}_{2 \times 4} {\bar{\bar{W}}}_{1} {(k, z)}_{4 \times 1}, \\ {\bar{\bar{V}}}_{2} {(k, z)}_{1 \times 1} = L_{2}_{1 \times 2} {\bar{\bar{W}}}_{2} {(k, z)}_{2 \times 1}, \end{array}}$

(28)

式中,矩阵L₁,L₂具体形式见附录。联立式（25）,（28）并代入z=z_i可推得单层结构的刚度矩阵K₁：

$\begin{array}{l} {(\begin{matrix} - {\bar{\bar{V}}}_{1} (k, 0) \\ {\bar{\bar{V}}}_{1} (k, z_{i}) \end{matrix})}_{4 \times 1} = K_{14 \times 4} {(\begin{array}{l} \bar{\bar{U}}' (k, 0) \\ \bar{\bar{U}}' (k, z_{i}) \end{array})}_{4 \times 1} \\ = (\begin{array}{l} - L_{1} I (k) \\ L_{1} T_{1} (k, z_{i}) \end{array}) {(\begin{matrix} I & 0 \\ T_{11} (k, z_{i}) & T_{12} (k, z_{i}) \end{matrix})}^{- 1} (\begin{array}{l} \bar{\bar{U}}' (k, 0) \\ \bar{\bar{U}}' (k, z_{i}) \end{array}) \end{array}$ 。

(29)

同理可以推得v′位移变量对应的刚度矩阵：

${(\begin{matrix} - {\bar{\bar{V}}}_{2} (k, 0) \\ {\bar{\bar{V}}}_{2} (k, z) \end{matrix})}_{2 \times 1} = K_{22 \times 2} {(\begin{array}{l} \bar{\bar{V}}' (k, 0) \\ \bar{\bar{V}}' (k, z) \end{array})}_{2 \times 1}$ 。

(30)

上述刚度矩阵与层厚、材料属性及波数有关。对于下覆半无限体,此法存在伪逆不满足边界条件的困难,对此笔者参考Gao等^[25]所用方法推导对应刚度矩阵。以 $\bar{\bar{U}}'$ 为例,联立式（14）～（16）并整理得

$\begin{matrix} K_{22} \frac{d^{2} \bar{\bar{U'}}}{d z^{2}} + (K_{21} + K_{12}) \frac{d \bar{\bar{U'}}}{d z} - \\ (K_{11} - ρ {(ω - c k_{x})}^{2} I) \bar{\bar{U'}} = 0 \end{matrix}$ ,

(31)

$V_{1} = K_{22} \frac{d^{2} \bar{\bar{U'}}}{d z^{2}} + K_{21} \bar{\bar{U'}}$ ,

(32)

式中,K_ij为系数矩阵。

整理式（31）,（32）得

$\frac{d}{d z} {(\begin{matrix} \bar{\bar{U'}} & \bar{\bar{V_{1}}} \end{matrix})}^{T} = H_{4 \times 4} {(\begin{matrix} \bar{\bar{U'}} & \bar{\bar{V_{1}}} \end{matrix})}^{T}$ ,

(33)

式中,H系数矩阵有如下形式：

$\begin{array}{l} H = [\begin{matrix} - K_{22}^{- 1} K_{21} & K_{22}^{- 1} \\ K_{11} - ρ {(ω - c k_{x})}^{2} I' + K_{12} K_{22}^{- 1} K_{21} & - K_{12} K_{22}^{- 1} \end{matrix}], \\ I' = [\begin{matrix} 1 & 0 \\ 0 & k^{- 2} \end{matrix}] 。 \end{array}}$

(34)

求H矩阵的右特征向量：

$H Θ = Θ Λ$ ,

(35)

式中,Θ,Λ由特征向量及特征值组成：

$\begin{array}{l} Θ_{4 \times 4} = [\begin{matrix} Θ_{11} & Θ_{12} \\ Θ_{21} & Θ_{22} \end{matrix}], \\ Λ_{4 \times 4} = d i a g (λ_{i}, \dots, - λ_{i}, \dots) 。 \end{array}}$

(36)

易推得

$(\begin{matrix} \bar{\bar{U'}} \\ \bar{\bar{V_{1}}} \end{matrix}) = [\begin{matrix} Θ_{11} & Θ_{12} \\ Θ_{21} & Θ_{22} \end{matrix}] {[c_{i} \exp (λ_{i} z), \dots, c_{i} \exp (- λ_{i} z)]}^{T}$ 。

(37)

因半空间体在边界无穷处广义位移存在极限,故取 $c_{i} \exp (- λ_{i})$ 所对应Θ矩阵列元素计算：

$\bar{\bar{V_{1}}} = (Θ_{22} Θ_{12}^{- 1}) \bar{\bar{U'}} = K_{3} \bar{\bar{U'}}$ 。

(38)

同理可以推得 $\bar{\bar{V'}}, \bar{\bar{V_{2}}}$ 对应刚度矩阵 $K_{4}$ 。

对薄板动力控制方程式（6）同样进行傅里叶变换：

$D_{q} k^{4} \bar{\bar{γ}} - \bar{\bar{Q}} (x, y) + \bar{\bar{P}} (x, y) = ρ_{b} h_{b} {(ω - c k_{x})}^{2} \bar{\bar{γ}} (x, y)$ 。

(39)

2. 体系控制方程

综上利用数学技巧推得层状体和半空间体对应刚度矩阵。以下根据边界条件推导下覆半空间的层状横观各向同性体总体刚度方程。

2.1 边界条件

对本问题,各层间相接面完全连接,其广义应力、位移间协调关系为

$\begin{array}{l} \bar{\bar{V_{1}}} (k, H_{i}^{1}) = \bar{\bar{V_{1}}} (k, H_{i + 1}^{0}), \\ \bar{\bar{V_{2}}} (k, H_{i}^{1}) = \bar{\bar{V_{2}}} (k, H_{i + 1}^{0}), \\ \bar{\bar{U'}} (k, H_{i}^{1}) = \bar{\bar{U'}} (k, H_{i + 1}^{0}), \\ \bar{\bar{V'}} (k, H_{i}^{1}) = \bar{\bar{V'}} (k, H_{i + 1}^{0}), \end{array}}$

(40)

式中,H_i上标0/1分别代表上/下层界面。

层状结构与板体接触条件为

$\begin{array}{l} {\bar{\bar{V}}}_{1} (k, H_{1}^{0}) = {(0 \bar{\bar{P}} (k, 0))}^{T}, \\ {\bar{\bar{V}}}_{2} (k, H_{1}^{0}) = 0, \end{array}}$

(41)

$\begin{array}{l} \bar{\bar{U'}} (k, H_{1}^{0}) = {(\begin{matrix} \bar{\bar{u'}} (k, 0) & k^{2} \bar{\bar{γ}} (k, 0) \end{matrix})}^{T}, \\ \bar{\bar{V'}} (k, H_{1}^{0}) = \bar{\bar{v'}} (k, 0), \end{array}}$

(42)

式中, $\bar{\bar{P}} (k, 0)$ 为二重正变换下的层间接触应力, $\bar{\bar{γ}} (k, 0)$ 为板体挠度变换值。

板上作用有移动简谐矩阵荷载,其表达式如下：

$\begin{matrix} Q (x, y, t) = \frac{q e^{i ω t}}{4 l_{1} l_{2}} [(H (x + l_{1}) - H (x - l_{1})) \cdot \\ (H (y + l_{2}) - H (y - l_{2}))] \end{matrix}$ ,

(43)

式中,q为矩形荷载作用区域合力,l₁,l₂为矩形区域尺寸,H为单位阶跃函数。

2.2 整体方程建立

依据上述关系,联立式（29）,（30）,（38）,（40）～（42）即得TI层状地基–薄板结构的整体动力刚度方程式（44）,（45）,当广义应力取单位值时即可求得动力Green函数。位移分量为复数,其模值可作振动幅值,幅角则表示相位情况。

此外当层数n较大时,刚度矩阵求逆成本剧增,高效算法的编制难度会大于模型求解本身。卢正等^[26]及杨广庆^[27]对路基的动力学研究表明,力学参数差异不大的层状体可用Odemark当量理论换算成等效层结构,且理论计算结果可满足工程设计需求。

综上建立由薄板、第一层TI地基、第二层等效TI地基、半无限TI空间组成的典型力学模型,进而分析其动力响应规律。

$\begin{matrix} {(\begin{matrix} 0 \\ ⋮ \\ 0 \end{matrix})}_{2 n + 2} = {(\begin{matrix} K_{21}^{1} & K_{22}^{1} & \dots & 0 \\ K_{23}^{1} & K_{24}^{1} + K_{21}^{2} & ⋮ \\ ⋮ & ⋱ & K_{24}^{n - 1} + K_{21}^{n} & K_{22}^{n} \\ 0 & \dots & K_{23}^{n} & K_{24}^{n} + K_{21}^{n + 1} \end{matrix})}_{{(2 n + 2)}^{2}} \cdot \\ {(\begin{matrix} \bar{\bar{v'}} (k, 0) \\ ⋮ \\ \bar{\bar{V}}' (k, H_{n}^{1}) \end{matrix})}_{2 n + 2} \end{matrix}$ ,

(44)

$\begin{matrix} {(\begin{matrix} 0 \\ \bar{\bar{P}} (k, 0) \\ \begin{array}{l} ⋮ \\ 0 \end{array} \end{matrix})}_{2 n + 2} \\ = {(\begin{matrix} K_{11}^{1} & K_{12}^{1} & \dots & 0 \\ K_{13}^{1} & K_{14}^{1} + K_{11}^{2} & \dots & ⋮ \\ ⋮ & \dots & K_{14}^{n - 1} + K_{11}^{n} & K_{12}^{n} \\ 0 & \dots & K_{13}^{n} & K_{14}^{n} + K_{11}^{n + 1} \end{matrix})}_{{(2 n + 2)}^{2}} \cdot \\ {(\begin{array}{l} \bar{\bar{u'}} (k, 0) \\ k^{2} \bar{\bar{γ}} (k, 0) \\ ⋮ \\ \bar{\bar{U}}' (k, H_{n}^{1}) \end{array})}_{2 n + 2} \end{matrix}$ 。

(45)

3. 数值计算及参数分析

联立式（7）,（39）,（43）先后解得 $\bar{\bar{P}}$ 、变换位移,经式（13）逆变换最终得位移振幅值。为避免积分奇异性并考虑土体黏弹性,对材料参数引入阻尼系数 $ξ$ ,形如 $E^{'}$ =E(1+ξi)。表面位移w可通过式（45）单独求解,工程应用中常将该值作为基本控制指标,下文以竖向位移为主讨论。

利用符号计算求解上述方程较为繁琐,拟采用数值计算求解空间域位移。具体对波数域位移变量作离散快速傅里叶逆变换,经试算在|k_x|,|k_y|<10、2048×2048积分域下所得结果可满足精度要求。本文使用Mathematica进行公式推导,数值计算程序则采用MATLAB编制。

3.1 验证

因TI层状地基-薄板结构动力响应研究相对较少,令各层参数一致（板厚h_b≈0）将模型退化为半空间体,重复已有算例参数计算说明上述方法准确性。

算例一研究移动点载下半空间体响应规律,BA等^[17]所建模型参数为 $μ_{h}$ = $μ_{vh}$ =0.25,E_h=5×10⁹ N/m²,E_v=2E_h,G_v=0.6E_h,ρ=2000 kg/m³, $c = 0.5 \sqrt{G_{v} / ρ}$ ,其解出（0,0,10）处无量纲竖向位移u′_z=G_vzu_z /q分布规律如图2所示（已将时间坐标转化为移动系下x坐标）。

图 2 移动恒点载作用TI半空间顶面竖向位移计算结果

Figure 2. Calculated results of top vertical displacement of TI half space under moving dead point loads

下载: 全尺寸图片幻灯片

从图2可见,上文算法所得结果与原结果吻合度较高,说明导出模型具备一定的合理性。

算例二研究移动矩形荷载下半空间-薄板结构响应规律,Muho^[22]所建模型参数为 $μ_{h}$ = $μ_{vh}$ =0.35,E_h=0.2×10⁹ N/m²,E_v=10E_h,G_v=3E_h,ρ=2100 kg/m³,E_p=30×10⁹ N/m²,μ_p=0.2,ρ_b=2300 kg/m³,h_b=0.2 m,l₁=l₂=0.15 m,q=80 kN,c=40 m/s,其解出板上中心点无量纲挠度γ′=γE_p/(qh_b)时程曲线如图3所示。

图 3 移动矩形荷载作用下TI半空间-薄板结构板体位移时程曲线计算结果

Figure 3. Calculated results of time-history curve of displacement of TI half-space-thin plate structure under moving rectangular loads

下载: 全尺寸图片幻灯片

由图3可见,本文模型结果与既有结果基本一致。两算例对比表明本文推导模型及数值算法具备一定的准确性,可用于移动荷载作用下TI层状地基-薄板结构的动力响应研究。

3.2 地基各向异性的影响

利用材料水平、竖直弹性模量比n表示各向异性程度,n取值一般在0.5～3.0,则对比n为0.5,1.0（各向同性处理）,1.5时单层地基的顶面竖向位移。板体及荷载参数：E_p=30×10⁹ Pa,μ_p=0.15,h_b=0.25 m,ρ_b=2400 kg/m³,c=200 m/s,ω=40 rad/s,l₁=0.25 m,l₂=0.25 m,q=100000 Pa。

n=0.5时模型各部分力学参数如表1所示,其余工况通过固定E_v,改变E_h调整n值,如n=1.5时,第一层弹模E_v=10×10⁹ Pa、E_h=15×10⁹ Pa,下文同。

表 1 层状地基各部分基本材料参数

Table 1. Basic material parameters of each part of layered foundation

参数/位置	第一层	半无限空间
E_v/(10⁹ Pa)	10	20
E_h/(10⁹ Pa)	5	10
G_v/(10⁹ Pa)	3	3
μ_h	0.25	0.25
μ_vh	0.25	0.25
z/m	10	∞
ρ/(kg·m^-3)	2000	2000

下载: 导出CSV

| 显示表格

各工况x方向板体位移幅值波形如图4所示,为便于比较曲线最值,图例及图左侧附有文字、楔形标记说明最值大小及位置,下文同。由图可知板体中心位移幅值在n=0.5时最大、n=1.5时最小。板体水平面波形受到地基各向异性影响,n为1.0,1.5时波形曲线呈对称分布,n=0.5时波形曲线向荷载移动方向偏移,虽出现零振幅点,但荷载影响范围显著增大。可见此类计算工况下,地基各向异性水平影响顶部面波传递规律,荷载作用中心板体位移幅值变化最大,相比各向同性情形出现一定幅度增减。

图 4 不同地基水平、竖直弹性模量比n下板体位移幅值在x方向的波形（y=0）

Figure 4. Waveshapes of position displacement amplitude of plate under different foundation levels and vertical elastic modulus ratios in x direction (y=0)

下载: 全尺寸图片幻灯片

3.3 荷载移动速率的影响

为研究荷载移动速率对动力响应的影响,取c为0,100,200,300,400 m/s 5种情况对比说明荷载移动速度对顶面竖向位移幅值的影响规律。模型取为单层地基-薄板结构,其余参数均与3.2节一致。

各工况x方向板体位移幅值波形如图5所示,可见移动速度对顶面竖向位移影响较大。当c从0逐渐增长至400 m/s时,波形出现不对称,零振幅点向移动方向偏移。增至300 m/s时,正向波形呈放大趋势,增至400 m/s时,正向波形迅速缩小。此外中心幅值随移动速率变化先增大后减小,300 m/s时幅值最大。可见此类计算工况下,临界荷载移动速率在300～400 m/s,其变化直接影响荷载作用区域面波具体分布形式。

图 5 不同荷载移动速率下板体位移幅值在x方向的波形（y=0）

Figure 5. Waveshapes of plate position displacement amplitude in x direction at different load movement rates (y=0)

下载: 全尺寸图片幻灯片

3.4 荷载加载频率的影响

荷载移动会影响水平向波形分布,为单独分析荷载加载频率对动力响应的影响,取c=0 m/s,ω为0,10,20,30,40 rad/s对比计算。模型取为单层地基-薄板结构,其余参数均与3.2节一致。

各工况x方向板体位移幅值波形如图6所示。由图可知：板体中心位移幅值受加载频率影响较小,不同频率主要影响水平面位移波形分布。当ω=0 rad/s时,位移沿水平轴的分布范围最广,荷载作用范围（拐点）延伸到了半径20 m之外。增大加载频率,荷载作用范围缩小,ω=40 rad/s时,波形集中分布在10 m之内。可见荷载加载频率主要影响顶面板体位移波形的集中程度或动荷载作用范围大小。

图 6 不同荷载加载频率下板体位移幅值在x方向的波形（y=0）

Figure 6. Waveshapes of plate position displacement amplitude in x direction under different loading frequencies (y=0)

下载: 全尺寸图片幻灯片

3.5 层状结构的影响

为研究地基各层结构参数改变对板体位移幅值的影响,对双层地基-板体结构按表2调整形成4种工况,由此分析各层地基对动力响应的影响。各工况材料、荷载参数在前基础上调整。

表 2 各计算工况下单层地基弹性模量比

Table 2. Elastic modulus ratios of single-layer foundation under various calculation conditions

层数/工况参数	水平、竖直弹性模量之比n				层厚/m
层数/工况参数	情形一	情形二	情形三	情形四	层厚/m
第一层	0.5	0.75	0.5	0.5	5
第二层	0.5	0.5	0.75	0.5	5
半无限空间	0.5	0.5	0.5	0.75	∞

下载: 导出CSV

| 显示表格

各工况计算结果如图7所示,由图可知：不同层n值改变或提升某层力学性能对顶面竖向位移影响不同。增大第二层或半空间n值对中心位移影响较小,而增加第一层将使位移大幅减小。此外单独增加各层n值会使荷载作用范围相比一般情形增加,且移动方向增长幅度更大。可见第一层各向异性水平变化将对顶面竖向位移产生较大影响。

图 7 某层地基水平、竖向弹性模量比改变时板体位移幅值在x方向的波形（y=0）

Figure 7. Waveshapes of plate position displacement amplitude in x direction when the horizontal and vertical elastic modulus ratio of a certain layer of foundation is changed (y=0)

下载: 全尺寸图片幻灯片

4. 结论

本文在已有各向异性体弹性动力学理论基础上,针对横观各向同性层状地基–薄板（无限大）这一典型力学模型,推导出二重积分形式的稳态响应解,并通过数值计算得到移动简谐荷载作用下的半解析解答。进而研究材料、荷载参数对板体位移的影响规律。

（1）对比横观各向同性、各向同性地基下计算结果,板体位移响应最值出现一定幅度增减。现有框架下的计算、设计结果可能与实际情况存在一定偏差。

（2）移动荷载相较静载使板体位移幅值的水平分布出现左右不对称,不对称程度与具体移动速率大小有关。存在临界速度使得荷载作用区域幅值最大。

（3）荷载加载频率对板体位移响应影响不如移动速率大,其主要影响荷载作用下的响应范围。

（4）分层各向异性水平改变对顶面竖向位移的影响程度不同,第一层地基相比其余层对降低竖向位移响应有较大作用,可利用此特性改善层状体系结构在动力作用下的响应特征。

附录：

（1）A矩阵元素形式

$\begin{array}{l} A_{13} = 1, A_{24} = 1, A_{31} = [C_{11} k^{2} - ρ {(ω - c k_{x})}^{2}] / C_{44}, \\ A_{34} = (C_{13} + C_{44}) / C_{44}, \\ A_{42} = [C_{44} k^{2} - ρ {(ω - c k_{x})}^{2}] / C_{44}, \\ A_{43} = - k^{2} (C_{13} + C_{44}) / C_{33}, \\ A_{11} = A_{12} = A_{14} = A_{21} = A_{22} = 0, \\ A_{23} = A_{32} = A_{33} = A_{41} = A_{44} = 0 。 \end{array}$

（2）B矩阵的元素形式

$\begin{array}{l} B_{11} = B_{22} = 0, B_{12} = 1, \\ B_{21} = [k^{2} (C_{11} - C_{12}) / 2 - ρ {(ω - c k_{x})}^{2}] / C_{44} 。 \end{array}$

（3）L₁,L₂矩阵的元素形式

$L_{1} = [\begin{matrix} 0 & - C_{44} & C_{44} & 0 \\ C_{13} & 0 & 0 & C_{33} / k^{2} \end{matrix}], L_{2} = [\begin{matrix} 0 & C_{44} \end{matrix}]$ 。

（4）a,b变量形式

图 1 传统PSDM

Figure 1. Traditional PSDM

下载: 全尺寸图片幻灯片

图 2 全连接神经网络结构

Figure 2. Structure of FCNN

下载: 全尺寸图片幻灯片

图 3 动剪切模量比和动阻尼比曲线

Figure 3. G/G_max-γ and λ-γ cureves of soil layers

下载: 全尺寸图片幻灯片

图 4 地铁车站结构横断面示意图

Figure 4. Schematic diagnam of cross-section of subway station structure

下载: 全尺寸图片幻灯片

图 5 钢筋和混凝土的本构模型

Figure 5. Constitutive models for concrete and reinforcement

下载: 全尺寸图片幻灯片

图 6 土-结构相互作用有限元模型

Figure 6. Finite element model for soil-structure model

下载: 全尺寸图片幻灯片

图 7 地震动加速度反应谱曲线

Figure 7. Spectral curves of seismic acceleration response

下载: 全尺寸图片幻灯片

图 8 损失函数变化趋势

Figure 8. Tendency of loss function

下载: 全尺寸图片幻灯片

图 9 FCNN对于训练集预测结果

Figure 9. Predicted results of FCNN for training set

下载: 全尺寸图片幻灯片

图 10 基于FCNN和对数线性回归结果的残差分布

Figure 10. Residual distribution based on FCNN and log-linear regression results

下载: 全尺寸图片幻灯片

图 11 FCNN模型的测试结果

Figure 11. Predicted results of FCNN model for testing set

下载: 全尺寸图片幻灯片

图 12 基于FCNN权值绝对值和的IM排序

Figure 12. Ranking of IM based on FCNN

下载: 全尺寸图片幻灯片

图 13 IM指标与θ对数线性回归图

Figure 13. IM and θ Log linear regression

下载: 全尺寸图片幻灯片

图 14 基于最大层间位移角的各IM指标α值

Figure 14. Values of α based on maximum layer drift

下载: 全尺寸图片幻灯片

图 15 基于最大层间位移角的各IM指标q值

Figure 15. Values of q based on maximum layer drift

下载: 全尺寸图片幻灯片

图 16 基于最大层间位移角的各IM指标ζ值

Figure 16. Values of ζ based on maximum layer drift

下载: 全尺寸图片幻灯片

表 1 采用的地震强度指标

Table 1 Selected intensity measures

类型	编号	名称	定义
加速度型	1	峰值加速度(PGA)	PGA = max(\|a(t)\|)
	2	Arias强度(AI)	${\text{AI}} = \frac{{\rm{\mathsf{π}}}}{{{\text{2g}}}}\int_0^{{t_{\text{f}}}} a {(t)^2}{\text{d}}t$
	3	持续最大加速度(SMA)	加速度时程中第三个最大值
	4	均方根加速度(a_rms)	${a_{{\text{rms}}}} = \sqrt {\int_{{t_{\text{5}}}}^{{t_{{\text{95}}}}} {a{{{\text{(}}t{\text{)}}}^{\text{2}}}} {\text{d}}t/{t_{\text{d}}}}$
	5	特征强度(I_C)	I_C = a_rms^1.5t_d^0.5
速度型	6	累计绝对速度(CAV)	${\text{CAV}} = \int_0^{{t_{\text{f}}}} {\left\| {a(t)} \right\|{\text{d}}t}$
	7	峰值速度(PGV)	PGV = max(\|v(t)\|)
	8	均方根速度(v_rms)	${v_{{\text{rms}}}} = \sqrt {\int_{{t_5}}^{{t_{95}}} {v{{(t)}^2}{\text{d}}t/{t_{\text{d}}}} }$
	9	能量密度指标(SED)	${\text{SED}} = \int_0^{{t_{\text{f}}}} {v{{(t)}^2}{\text{d}}t}$
	10	复合指标(I_F)	I_F = PGV × t_d^0.25
	11	持续最大速度(SMV)	速度时程中第三个最大值
位移型	12	峰值位移(PGD)	PGD = max(\|d(t)\|)
位移型	13	均方根位移(d_rms)	${d_{{\text{rms}}}} = \sqrt {\int_{{t_5}}^{{t_{95}}} {d{{(t)}^2}{\text{d}}t/{t_{\text{d}}}} }$
与反应谱单值相关	14	伪加速度谱峰值(PSA)	PSA=max(PS_a)
与反应谱单值相关	15	伪速度谱峰值(PSV)	PSV=max(PS_V)
与反应谱时间段积分相关	16	加速度谱强度(ASI)	${\text{ASI}} = \int_{0.1}^{0.5} {P{S_{\text{a}}}{\text{d}}T}$
	17	速度谱强度(VSI)	${\text{VSI}} = \int_{0.1}^{2.5} {{S_{\text{v}}}{\text{d}}T}$
	18	Housner强度(HI)	${\text{HI}} = \int_{0.1}^{2.5} {P{S_{\text{v}}}{\text{d}}T}$
	19	有效峰值加速度(EPA)	${\text{EPA}} = {\text{mean(}}\int_{{\text{0}}{\text{.1}}}^{{\text{0}}{\text{.5}}} {{S_{\text{a}}}{\text{d}}T/2.5} )$

下载: 导出CSV

表 2 Ⅱ类场地土层物理参数

Table 2 Physical parameters of Type Ⅱ site soil

土层	土类	土层厚度/m	密度/(kg·m^-3)	剪切波速/(m·s^-1)	剪切模量/MPa
1	人工填土	5.0	1750	180	57
2	粉质黏土	10.0	1900	250	118.8
3	细中砂	10.0	2000	300	180
4	细粉砂	15.0	2000	320	205
5	卵石	20.0	2280	500	525

下载: 导出CSV

表 3 计算效率对比

Table 3 Comparison of efficiency

耗时项目	数值模拟	FCNN
数据处理	—	10 min
模型训练^*	—	2 h
响应计算	20 min×30=10 h	2 s
合计	10 h	2.2 h
^*注：FCNN训练耗时包括搭建至调至最优模型耗时。

下载: 导出CSV

参考文献(36)

[1]	CORNELL C A, JALAYER F, HAMBURGER R O, et al. Probabilistic basis for 2000 SAC federal emergency management agency steel moment frame guidelines[J]. Journal of Structural Engineering, 2002, 128(4): 526-533 doi: 10.1061/(ASCE)0733-9445(2002)128:4(526)
[2]	钟紫蓝, 申轶尧, 郝亚茹, 等. 基于IDA方法的两层三跨地铁地下结构地震易损性分析[J]. 岩土工程学报, 2020, 42(5): 916-924. doi: 10.11779/CJGE202005014 ZHONG Zilan, SHEN Yiyao, HAO Yaru, et al. Seismic fragility analysis of two-story and three-span metro station structures based on IDA method[J]. Chinese Journal of Geotechnical Engineering, 2020, 42(5): 916-924. (in Chinese) doi: 10.11779/CJGE202005014
[3]	ARGYROUDIS S A, PITILAKIS K D. Seismic fragility curves of shallow tunnels in alluvial deposits[J]. Soil Dynamics and Earthquake Engineering, 2012, 35: 1-12. doi: 10.1016/j.soildyn.2011.11.004
[4]	LIU T, CHEN Z Y, YUAN Y, et al. Fragility analysis of a subway station structure by incremental dynamic analysis[J]. Advances in Structural Engineering, 2016, 20(7): 1111-1124.
[5]	RAMAMOORTHY S K, GARDONI P, BRACCI J M. Probabilistic demand models and fragility curves for reinforced concrete frames[J]. Journal of Structural Engineering, 2006, 132(10): 1563-1572. doi: 10.1061/(ASCE)0733-9445(2006)132:10(1563)
[6]	SIMON J, BRACCI J M, GARDONI P. Seismic response and fragility of deteriorated reinforced concrete bridges[J]. Journal of Structural Engineering, 2010, 136(10): 1273-1281. doi: 10.1061/(ASCE)ST.1943-541X.0000220
[7]	HE Z M, XU H, GARDONI P, et al. Seismic demand and capacity models, and fragility estimate for underground structures spatially varying soil properties[J]. Tunnelling and Underground Space Technology incorporating Trenchless Technology Research, 2022(119): 104231.
[8]	XU H, GARDONI P. Probabilistic capacity and seismic demand models and fragility estimates for reinforced concrete buildings based on three-dimensional analyses[J]. Engineering Structures, 2016(112): 200-214.
[9]	ZHANG W, SHOKRABADI M, BOZORGNIA Y, et al. A methodology for fragility analysis of buried water pipes considering coupled horizontal and vertical ground motions[J]. Computers and Geotechnics, 2020, 126(4): 103709.
[10]	ZHENG K F, CHEN L B, ZHUANG W L, et al. Bridge vulnerability analysis based on probabilistic seismic demand models[J]. Engineering Mechanics, 2013, 30(5): 165-171.
[11]	于晓辉, 吕大刚, 王光远. 关于概率地震需求模型的讨论[J]. 工程力学, 2013, 30(8): 172-179. https://www.cnki.com.cn/Article/CJFDTOTAL-GCLX201308027.htm YU Xizohui, LÜ Dagang, WANG Guangyuan. Discussions on probabilistic seismic demand models[J]. Engineering Mechanics, 2013, 30(8): 172-179. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-GCLX201308027.htm
[12]	YANG C T, XIE L L, LI A Q, et al. Ground motion intensity measures for seismically isolated RC tall buildings[J]. Soil Dynamics and Earthquake Engineering, 2019, 125: 105727. doi: 10.1016/j.soildyn.2019.105727
[13]	PADGETT J E, NIELSON B G, DESROCHES R. Selection of optimal intensity measures in probabilistic seismic demand models of highway bridge portfolios[J]. Earthquake Engineering and Structural Dynamics, 2008, 37(5): 711-725. doi: 10.1002/eqe.782
[14]	ZHANG C M, ZHAO M, ZHONG Z L, DU X L. Seismic intensity measures and fragility analysis for subway stations subjected to near-fault ground motions with velocity pulses[J]. Journal of Earthquake Engineering, 2022, 26(16): 8724-8750. doi: 10.1080/13632469.2021.1994056
[15]	DING Y J, CHEN J, SHEN J X. Prediction of spectral accelerations of aftershock ground motion with deep learning method[J]. Soil Dynamics and Earthquake Engineering, 2021, 150(25): 106951.
[16]	DING Y J, CHEN J, SHEN J X. Conditional generative adversarial network model for simulating intensity measures of aftershocks[J]. Soil Dynamics and Earthquake Engineering, 2020, 139(17): 106281.
[17]	HUANG P F, CHEN Z Y. Deep learning for nonlinear seismic responses prediction of subway station[J]. Engineering Structures, 2021, 244(13): 112735.
[18]	XU Z K, CHEN J, SHEN J X, XIANG M J. Recursive long short-term memory network for predicting nonlinear structural seismic response[J]. Engineering Structures, 2022, 250(26): 113406.
[19]	TORKY AHMED A, SUSUMU O. Deep learning techniques for predicting nonlinear multi-component seismic responses of structural buildings[J]. Computers and Structures, 2021, 252(15): 106570.
[20]	HUANG P F, CHEN Z Y. Fragility analysis for subway station using artificial neural network[J]. Journal of Earthquake Engineering, 2021(3): 1-21.
[21]	RIDDELL R. On ground motion intensity indices[J]. Earthquake Spectra, 2007, 23(1): 147-173. doi: 10.1193/1.2424748
[22]	PENZIEN J. Seismically induced racking of tunnel linings[J]. Earthquake Engineering and Structural Dynamics, 2000, 29(5): 683-691. doi: 10.1002/(SICI)1096-9845(200005)29:5<683::AID-EQE932>3.0.CO;2-1
[23]	CHEN Z Y, LIU Z Q. Stochastic seismic lateral deformation of a multi-story subway station structure based on the probability density evolution method[J]. Tunnelling and Underground Space Technology, 2019, 103114.
[24]	城市轨道交通结构抗震设计规范[S]. 北京: 中国建筑工业出版社, 2014. Code for Seismic Design of Urban Rail Transit Structures[S]. Beijing: China Architecture & Building Press, 2014. (in Chinese)
[25]	陈国兴, 刘雪珠, 朱定华, 等. 江苏长江以南地区新近沉积土动剪切试验研究[J]. 地下空间与工程学报, 2007, 3(4): 745-750. https://www.cnki.com.cn/Article/CJFDTOTAL-BASE200704034.htm CHEN Guoxing, LIU Xuezhu, ZHU Dinghua, et al. A study on dynamic shear of recently deposited soils for southern Jiangsu Province along Yangtze River[J]. Chinese Journal of Underground Space and Engineering, 2007, 3(4): 745-750. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-BASE200704034.htm
[26]	杜修力, 许紫刚, 许成顺, 等. 基于等效线性化的土-地下结构整体动力时程分析方法研究[J]. 岩土工程学报, 2018, 40(12): 2155-2163. doi: 10.11779/CJGE201812001 DU Xiuli, XU Zigang, XU Chengshun, et al. Time-history analysis method for soil-underground structure system based on equivalent linear method[J]. Chinese Journal of Geotechnical Engineering, 2018, 40(12): 2155-2163. (in Chinese) doi: 10.11779/CJGE201812001
[27]	曲哲, 叶列平. 基于有效累积滞回耗能的钢筋混凝土构件承载力退化模型[J]. 工程力学, 2011, 28(6): 45-51. https://www.cnki.com.cn/Article/CJFDTOTAL-GCLX201106009.htm QU Zhe, YE Lieping. Strength deterioration model based on effective hysteretic energy dissipation for rc members under cyclic loading[J]. Engineering Mechanics, 2011, 28(6): 45-51. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-GCLX201106009.htm
[28]	TSINIDIS G, PITILAKIS K, TRIKALIOTIA D. Numerical simulation of round robin numerical test on tunnels using a simplified kinematic hardening model[J]. Acta Geotechnica, 2013, 9(4): 641-659.
[29]	YAROTSKY D. Error bounds for approximations with deep ReLU networks[J]. Neural Networks, 2017, 94: 103-114. doi: 10.1016/j.neunet.2017.07.002
[30]	KINGMA D P, BA J L. Adam: a method for stochastic optimization[J]. Computer Research Repository, 2014(1412): 74-89.
[31]	HE K M, ZHANG X Y, REN S Q, SUN J. Delving deep into rectifiers: surpassing human-level performance on imagenet classification[C]//2015 IEEE International Conference on Computer Vision (ICCV), Santiago, Chile: Institute of Electrical and Electronics Engineers, 2015.
[32]	ZAREMBA W, SUTSKEVER I, VINYALS O. Recurrent neural network regularization[J]. Computer Research Repository, 2014(1409): 68-76.
[33]	ABADI M, BARHAM P, CHEN J M, et al. TensorFlow: a system for large-scale machine learning[J]. Computer Research Repository, 2016(1605): 1-18.
[34]	SREEJAYA K P, BASU J, RAGHUKANTH S, SRINAGESH D. Prediction of ground motion intensity measures using an artificial neural network[J]. Pure and Applied Geophysics, 2021, 178: 2025-2058. doi: 10.1007/s00024-021-02752-9
[35]	张成明, 钟紫蓝, 甄立斌, 等. 适用于圆形隧道损伤评价的地震动强度指标研究[J]. 工程力学, 2021, 38(1): 100-108. https://www.cnki.com.cn/Article/CJFDTOTAL-GCLX202101011.htm ZHANG Chengming, ZHONG Zilan, ZHEN Libin, et al. Seismic intensity measures for the damage evaluation of circular tunnels[J]. Engineering Mechanics, 2021, 38(1): 100-108. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-GCLX202101011.htm
[36]	LUCO N, CORNELL C A. Structure-specific scalar intensity measures for near-source and ordinary earthquake ground motions[J]. Earthquake Spectra, 2007, 23(2): 357-392. doi: 10.1193/1.2723158

施引文献

资源附件(1)

其他相关附件
- PDF格式
  12-202402-G22-1448审稿意见与作者答复点击下载(1141KB)

图(16) / 表(3)

计量

文章访问数: 362
HTML全文浏览量: 60
PDF下载量: 95
被引次数: 0

0. 引言
1. 直角坐标系下的振动方程
1.1 坐标系统变换
1.2 模型基本方程及变量解耦
1.3 控制方程求解
2. 体系控制方程
2.1 边界条件
2.2 整体方程建立
3. 数值计算及参数分析
3.1 验证
3.2 地基各向异性的影响
3.3 荷载移动速率的影响
3.4 荷载加载频率的影响
3.5 层状结构的影响
4. 结论
附录：

0. 引言
1. 直角坐标系下的振动方程
1.1 坐标系统变换
1.2 模型基本方程及变量解耦
1.3 控制方程求解
2. 体系控制方程
2.1 边界条件
2.2 整体方程建立
3. 数值计算及参数分析
3.1 验证
3.2 地基各向异性的影响
3.3 荷载移动速率的影响
3.4 荷载加载频率的影响
3.5 层状结构的影响
4. 结论
附录：

参考文献(36)

施引文献

资源附件(1)

基于全连接神经网络的地铁车站响应分析与地震强度指标优选 English Version

作者简介: 钟紫蓝(1986—)，男，博士，副研究员，主要从事地下结构抗震研究。E-mail: zilanzhong@bjut.edu.cn

通讯作者: 申家旭, E-mail: shenjx@bjut.edu.cn

计量

出版历程

Response analysis of subway station and optimization of seismic intensity measures based on fully connected neural network

0. 引言

1. 直角坐标系下的振动方程

1.1 坐标系统变换

1.2 模型基本方程及变量解耦

1.3 控制方程求解

2. 体系控制方程

2.1 边界条件

2.2 整体方程建立

3. 数值计算及参数分析

3.1 验证

3.2 地基各向异性的影响

3.3 荷载移动速率的影响

3.4 荷载加载频率的影响

3.5 层状结构的影响

4. 结论

附录：

其他相关附件

PDF格式

计量

出版历程

目录

0. 引言

1. 直角坐标系下的振动方程

1.1 坐标系统变换

1.2 模型基本方程及变量解耦

1.3 控制方程求解

2. 体系控制方程

2.1 边界条件

2.2 整体方程建立

3. 数值计算及参数分析

3.1 验证

3.2 地基各向异性的影响

3.3 荷载移动速率的影响

3.4 荷载加载频率的影响

3.5 层状结构的影响

4. 结论

附录：

作者简介:
钟紫蓝(1986—)，男，博士，副研究员，主要从事地下结构抗震研究。E-mail: zilanzhong@bjut.edu.cn

通讯作者:
申家旭, E-mail: shenjx@bjut.edu.cn