更全的杂志信息网

成层饱和介质平面波斜入射问题的一维化时域方法1)

更新时间:2009-03-28

引言

地震波输入是目前大型结构抗震分析中亟待解决的问题之一.以往很多工程的抗震设计中,常假定地震波垂直入射,从而将场地反应分析简化为一维问题.而实际上,地震波从震源传播至地表的过程中,因传播路径、地质条件等复杂因素的影响,到达工程场地时并不是垂直入射的[1-2].现有研究表明,不同方向入射波产生的结构反应与垂直入射的结果有明显差别,尤其对于较大尺寸的工程结构[3-4],如大跨度桥梁、大坝、大型水电站、地下管道和地铁隧道[5]等.为了确保工程结构在地震作用下的安全性,应考虑地震波斜入射的影响.

近年来,地震波斜入射时场地的地震反应研究取得了很大进展[6].以传递矩阵法和刚度矩阵法等方法为代表的频域方法,基于解析求解,可以得到弹性场地或者非线性场地的等效线性化反应.目前,已在地震波斜入射时二维成层弹性介质自由场[7-8]以及非线性场地等效线性化[9-10]求解方面得到广泛应用.但频域方法无法反应地震过程中土体实际的非线性行为,且其计算精度往往与时域有限元方法不匹配,随着大型结构--地基系统抗震分析的发展,建立高精度和实用性强的时域分析方法已经迫在眉睫.廖河山等[11]用黏性吸收边界条件模拟底层半空间的吸收波效应,基于Snell定律将水平方向无限的二维空间问题转化为一维问题,建立了SH波倾斜入射时土层的非线性响应时域分析方法;李山有等[12]基于波动传播水平视波速不变且已知的特点,利用局部透射边界条件及显式有限差分法的内部节点位移计算公式,给出了地震波斜入射下水平成层半空间自由场的简化时域计算方法.刘晶波和王艳[13-15]仍采用黏性吸收边界条件来近似模拟截断半空间的吸波效应,在集中质量有限元法和时间中心差分相结合建立结点运动方程组的基础上,根据Snell定律,将水平方向相邻结点的运动用该结点相邻时刻的运动表示,从而将水平成层弹性半空间在SH波和PSV波斜入射下二维自由波场的计算问题简化为时域内的一维问题求解;赵密等[16]在刘晶波等[14]算法的基础上,提出一种模拟基岩半空间辐射阻尼的人工边界条件代替黏性边界条件,可获得更高的计算精度;卓卫东等[17]在刘晶波等[14]算法的基础上,建立P-SV波斜入射下有阻尼成层弹性半空间自由波场求解的一维化时域算法;Zhao等[18]提出一种精确人工边界条件模拟底部半空间的波动辐射和输入特性,并利用Snell定律将水平方向无限的二维空间问题转化为沿深度方向的一维空间问题,建立了一种计算平面波斜入射下成层半空间自由场地反应的一维化有限元时域解法.

从这些成果上看,地震波斜入射时场地地震反应的时域分析方法已经得到很大发展,但这些成果都是基于成层单相介质展开的.实际上地球表面的地震波在从震源通过各种途径传播时,通常要穿过饱和土层和非饱和土层最后到达地面,当将充满液体的土层进行动力分析时,按饱和多孔介质理论分析比按单相介质理论更为合理.自1956年Biot[19]建立饱和多孔介质波传播理论至今,人们对于饱和多孔介质层的波动问题已经有了系统深入的认识[20],但对于地震波斜入射下成层饱和多孔介质自由场地的时域分析研究还处在起步阶段.王子辉[21]在赵成刚等[22]建立的饱和多孔介质时域显式有限元方法的基础上,采用刘晶波等[14]提出的一维化求解方法,发展形成了基岩上覆单相固体介质与两相饱和介质互层分布的二维自由波场的一维化时域计算方法.该方法可用于地震波斜入射下成层饱和多孔介质自由场的计算,但其中,同样采用黏性吸收边界条件来近似模拟基岩的吸波效应,使得求解精度降低;且因为水平方向也要划分网格,使得有限元离散误差与二维有限元方法一致.为了克服这些缺点,本文建立了一种新的求解地震波斜入射下基岩上覆饱和多孔介质成层场地自由场分析的一维化时域计算方法.首先根据Snell定律将饱和多孔介质二维空间问题转化为一维时域问题,通过对深度方向的有限元离散,推导饱和多孔介质一维化有限元方程,然后基于Zhao等[18]建立了单相弹性介质精确人工边界条件,在考虑基岩与饱和多孔介质间透水或不透水边界条件以及不同饱和多孔介质交界面边界条件的基础上,形成基岩上覆成层饱和介质系统的整体有限元方程,最后采用中心差分法与Newmark平均加速度近似格式相结合的方法对时间进行离散,得到节点的动力时程的显式表达.通过与已有的矩阵传递法的解进行对比,验证了本文方法的有效性.

1 成层饱和多孔介质模型

分析模型如图1所示,基岩上覆成层饱和多孔介质,平面P(SV)波从基岩层以θ(β)角入射到成层饱和多孔介质层中.饱和多孔介质共N层,编号由上至下分别为 1,2,···,N,厚度分别为 h1h2,···,hN,总厚度为H.考虑固体颗粒及孔隙流体的压缩性以及流固两相的惯性耦合和黏性耦合,Biot[19]给出了流体饱和多孔介质的矢量波动方程

  

图1 成层饱和多孔介质场地模型Fig.1 The fiel model of layered saturated poroelastic medium

 

其中,usuf分别为固相和液相位移;e=∇·us,ε = ∇ ·uf;ρ11= ρ1+ ρa22= ρ2+ ρa12= −ρa1=(1 − ns2=nρf;ρs为固相质量密度,ρf为液相质量密度,ρa为液固两相耦合质量密度;n为孔隙率;bn2/k(η为流体黏滞系数,k为渗透系数);A,N,R,Q为材料常数,可以分别用土骨架和孔隙流体的材料常数表示[23]

LTE-M在地铁隧道中的非稳态MIMO信道的几何建模······································陈旭敏 王大庆 潘韵天 郑国莘 (6,888)

 

饱和多孔介质骨架切向位移连续

SV波入射时

 

将式(7)代入式(4),得到

基岩视为弹性单相介质,其弹性常数为λjj,质量密度为ρj.

2 饱和多孔介质一维化显式有限元方法

一维化方法的关键是根据Snell定理将动力方程中对水平方向的偏导数转化为对时间的导数,将二维平面问题转化为一维时域问题.

2.1 饱和多孔介质波动方程的空间一维化

水泥稳定碎石基层所用水泥的初凝时间应不小于3h,终凝时间应为6~10h,且不得使用早强水泥、快硬水泥、受潮或安定期未到的水泥[2]。本文所使用的水泥为32.5级水泥,技术指标检测结果如表2所示。

我对野菜不大识得,别呦呦就教我认,那高有两尺,紫茎,叶子像薄荷,边上都是锯齿的,是马兰头;那叶子微长,每三片叶子攒生一处的,是苜蓿……

 

其中,D1D2D3为参数矩阵,表达式分别如下

 

L为微分算子矩阵

 

根据Snell定理,水平成层介质水平方向视波速相等,所以有

 

其中,cx是波沿x轴的视波速

文中首先进行实验一来确定最优的人脸人耳特征融合权重系数α,β。将权重系数α从0.1开始以0.1为步长逐步调整至0.9,对应的权重系数β也随之调整,β=1-α,统计每次系数对应的识别率,统计结果如图3所示。

 

式中和σf分别为固体骨架部分和流体部分承担的应力,为土骨架应变.

 

根据式(3),饱和多孔介质中的应力与位移的关系,可以表示为如下的矩阵关系

 

将式(7)代入,得

 

2.2 一维化方程的有限元离散

将每一层饱和多孔介质在z轴方向上离散为一定数量的二节点有限元单元,单元和节点编号由上而下进行编号.第j层饱和多孔介质的网格划分示意见图1.单元e的两节点在整体坐标系中的坐标为zizi+1,单元网格厚度为Δz=zi+1zi,固相和液相节点位移向量为

 

其中,N为形状函数矩阵

 

其中,

针对单元e,方程(9)的Galerkin弱式可以写成

②系统全面、突出重点原则。预算项目绩效指标应系统全面地反映项目绩效目标与工作内容,覆盖项目工作内容的所有重要方面和关键领域,并合理完整地设计关键指标。同时,抓住项目重点和本质,设计最具代表性的指标,与专业规划总体目标、项目年度任务数、计划数相对应,能清晰地反映投入和产出绩效之间的紧密对应关系。同时,项目关键绩效目标和重大经费支出必须有相关绩效指标进行反映,并对项目绩效目标的实现形成支撑、对应关系。在水利重大预算项目绩效指标管理工作中,一般占项目预算10%以上的单项支出,必须要求项目单位在产出指标中予以体现出来。

 

其中,feFe为单元载荷列阵,分别为

 

为单元集中质量矩阵,为单元阻尼矩阵,为单元刚度矩阵,表达式详见附录.

2.3 组装成整体有限元方程

在同一种饱和多孔介质内部,可以直接对单元有限元离散方程进行组装,对不同饱和多孔介质层之间,需要考虑介质层之间的边界条件.下面以两个单元(如图2)为例,说明不同饱和多孔介质层之间有限元离散方程的组装.

  

图2 单元节点的关系Fig.2 The relationship between element nodes

假定不同介质层间无相对运动,则不同饱和多孔介质层交界面的边界条件为:

饱和多孔介质骨架法向位移连续

 

其中,Ks,Kb分别为土颗粒和土骨架的体积模量;λ,µ是土骨架的Lame常数;Kf是孔隙流体的体积模量.

 

法向相对位移连续

 

饱和多孔介质骨架切向应力连续

 

法向总应力连续

 

孔隙水压力连续

页岩样品均采集自湖北某页岩气资源调查评价井页岩岩心;页岩标准物质GBW07107:地球物理地球化学勘查研究所、地质矿产部岩矿测试技术研究所研制。

方勋梅没有主动出击,程晓也懒得理她,只是每天,方勋梅刚刚钻出尼桑车,他的凯迪拉克也悄然驶入厂区;每天下班后,方勋梅的车刚一驶上公路,程晓的凯迪拉克就如一尾快鱼从她身边穿过,在暮色中倏地无影无踪。

 

不同饱和土层,下层饱和土层的各量用上角标“”以便区分.

为了便于应用边界条件,将饱和多孔介质单元的有限元方程式(16)展开为规范化形式

3)统一标准、科学评价。2016年基础性地理国情监测检查验收技术标准以普查时期为基础,2017年根据监测需要进行修改印发,全国执行统一标准,从国检中心到任务承担单位都积极开展质量培训,统一认识,明确要求,对成果质量依据标准进行客观、准确评价。

 

其中

 

上标(3−k)表示3−k阶导数.

 

其中,为该单元的孔隙率.将式(20)代入式(19),并将式(19b)中第2行和第4行分别对应加到式(19a)第2行和第4行上,式(19b)第2行和第4行等式左右两边除以单元的孔隙率ne得到

 
 

因为不同饱和多孔介质层间切向相对位移不存在连续性条件,进行单元有限元方程叠加时需要进行如下数学处理:在不增加约束条件的前提下在整体位移、速度和加速度列阵中各增加一个维度,并且保证等式的恒成立.处理后,不同饱和多孔介质层间相关联的两个单元叠后的结果可以表示为

本文首先对NDD体系下的网络入侵过程进行分析,指出其会因为脆弱性变换而具有非马尔可夫特性,进而导致入侵成功概率求解困难.为此,在随机抽样模型中引入时间概率密度函数,从脆弱性变换角度推导出NDD体系下入侵者的入侵成功概率计算公式.然后,分别给出单、多脆弱性变换条件下的入侵成功概率极限定理,并通过入侵成功概率计算公式极限变换的形式对定理进行证明,同时基于这两条定理给出单、多脆弱性变换情况下动态变换的最优变换空间计算方法.最后,通过仿真实验对文章的结论进行验证.

 

式(22)中

 
 

其中,下标t为节点号,U(k,l)为单元叠加后得到位移(速度或加速度)列阵.f1f2为单元叠加后得到的载荷列阵.分别表示t点隶属于上下两个单元时的土骨架水平方向位移和水平方向的相对位移,由边界条件可知

3 底部人工边界及输入

设波从基岩入射到饱和多孔介质层中,基岩视作单相弹性介质.Zhao等[18]提出了单相弹性介质精确动力人工边界条件

 

其中,上标 “j”表示基岩,为基岩应力为散射波速度矢量,为入射波速度矢量,式中

 

基岩层与第N层饱和多孔介质层交界面的边界条件为:

法向位移连续

 

切向位移连续

 

切向应力连续

 

法向总应力连续

 

当基岩透水时,孔隙水压力为零

 

当基岩不透水时,饱和多孔介质土骨架和液相法向位移相等

 

根据基岩层与第N层饱和多孔介质层交界面的边界条件式(24)以及人工边界条件式(23),装配全部有限单元,得到一维化后全场的有限元动力方程.

饱和多孔介质波动方程用矩阵形式表示为[22]

当基岩边界透水时

 

Ms,Msf,Mf,Mfs为整体质量矩阵,Cs,Csf,Cf,Cfs为整体阻尼矩阵,Ks,Ksf,Kf,Kfs为整体刚度矩阵,按式(22)的方式组装;u1u2为整体位移列阵

 

下标n为总节点数,t表示不同饱和多孔介质交界面节点;F1F2为整体载荷列阵

 

P波入射时

 

在Biot模型中,固相、液相部分的应力--应变关系分别为

 

当基岩边界不透水时

1.脑缺血模型的成功构建:(1)TTC染色:正常脑组织为红色,缺血脑组织为白色。见图1。(2)神经功能评分及死亡率:脑缺血模型家兔30只,神经功能评分1级4只,神经功能评分≥2级26只,其中4只神经功能评分4级的家兔于模型建立成功1 d后死亡。

 

以基岩透水为例,可以得到用前一时刻的位移和速度表示的当前时刻的位移

 

s为整体系数矩阵的总阶次;为饱和多孔介质层与基岩交界面面节点的孔隙水压力.

4 节点动力反应的显式表达

利用一维化的方法建立离散的有限元方程式(25)或式(27)之后,可采用中心差分法与Newmark平均加速度法相结合的显式积分方法求解各节点的动力时程[22].周正华等[24]和王进廷等[25]先后对比分析了几种常用显式积分方法的稳定性和精度,证明该方法具有二阶精度及良好的稳定性,并给出了其稳定性条件.李亮等[26]研究了采用该方法进行饱和两相介质波动问题分析时的稳定性影响因素,包括时间步长、空间步距和渗透系数取值等的作用规律.在此不再赘述.

其中

综上所述,在数字化时代背景下,出版编辑的内容服务随着时代的发展而不断进步。加之数字技术的日新月异,但凡可以承载数字内容的载体都可以改变数字编辑的从业方向,从最开始的网站、论坛到现在微博、微信、云服务、区块链、网络直播、人工智能等,不断变化的载体意味着编辑自身必须不断地学习,这样才能真正跟上时代的潮流,成为一名优秀的数字编辑。

 

用前一时刻的位移和速度以及当前时刻的位移表示的当前时刻的速度

 

由前一时刻的速度和加速度以及当前时刻的速度表示的此时刻的加速度

 

对于应力,按照传统有限元法,单元为常应力单元.为了提高应力结果的精度,可以采用如下两种方法计算.其一,采用单元离散方程(21)求解;其二,由式(12)求解,具体表达式如下

 

对于可以根据本构关系联立式(31)求得

 

为验证本文方法,采用如图3所示的2个不同场地模型计算,将计算结果与传递矩阵法结合傅里叶变换的计算结果进行对比.其中图3(a)为基岩上覆单一饱和多孔介质层,图3(b)为基岩上覆两层不

5 验证

  

图3 饱和多孔介质覆盖层场地简化模型Fig.3 The simplifie model of saturated soil overburden site

同饱和多孔介质层.基岩层和饱和介质层参数如表1和表2所示.考虑基岩与饱和介质层交界面为透水边界.

协同理论引入中国后,学界进行了深入研究和探讨。张刚(1997年)根据协同创新的路径不同,将协同分为内部协同和外部协同:内部协同是内部要素之间的互动与协同;外部协同是研究主体和外部要素之间的互动与协同。陈光(2005年)根据研究主体与相关主体之间的关系,将外部协同细分为横向协同和纵向协同:横向协同是指同一大类产业中细分产业主体间的协同;纵向协同是指同一功能链不同环节上主体之间的协同[8]。

 

表1 基岩的材料参数Table 1 Material parameters of bedrock

  

Parameters Value Units ρj 2385 kg/m3 λj 15.6 GPaµj 15.6 GPa

 

表2 饱和多孔介质层的材料参数Table 2 Material parameters of saturated poroelastic medium

  

Parameters Layer① Layer② Layer③ Units ρs 2650 2700 2700 kg/m3 ρf 1000 1000 1000 kg/m3 ρa 0 0 0 kg/m3 λ 15.6×103 22.0 26.2 MPaµ 15.6×103 22.0 26.2 MPa η 1.0×10−3 1.0×10−3 1.0×10−3 Pa·s kp 1.0×10−10 1.0×10−10 1.0×10−10 m2 n 0.1 0.60 0.27 —Ks 36.0 36.0 36.0 GPa Kf 2.0 2.0 2.0 GPa

5.1 地震波输入

入射波的位移时程采用狄拉克函数的有限差分近似,即

 

其中为单位脉冲的作用时间,A为单位脉冲峰值,H(τ)为Heaviside函数.计算中取T=0.5s,A=1m,对应的位移时程曲线如图4所示.

5.2 基岩上覆单一饱和多孔介质层

  

图4 入射波位移时程Fig.4 Time history of incident wave displacement

  

图5 模型一P波60°入射时半空间表面O点固相位移时程Fig.5 Time history of solid phase displacement for point O on the half-space surface excited by P waves with incident angle 60°for model one

  

图6 模型一SV波30°入射时半空间表面O点固相位移时程Fig.6 Time history of solid phase displacement for point O on the half-space surface excited by SV waves with incident angle 30°for model one

在图3(a)所示的模型中,取饱和多孔介质覆盖层厚度h=100m,饱和多孔介质取表2中第①组参数,网格厚度取1m,时间步长取10−4s.图5和图6分别给出采用本文方法得到的P波60°入射和SV波30°入射时半空间表面O点的固体骨架位移时程.由于本算例中所取饱和多孔介质参数特殊,孔隙率较小,固体骨架刚度与基岩刚度相同,且大于孔隙流体体积模量,因此此时的饱和多孔介质可以看作和基岩材料相近的单相介质[27],基岩上覆单一饱和多孔介质层模型相应地也可以看成是均匀半空间场地,图中的精确解是根据波传播理论得到的半空间场地地表位移的解析解[28-29].从图中可以看出,本文方法的计算结果与解析解完全吻合.一方面说明本文方法在模拟基岩上覆单一饱和多孔介质层场地时的有效性和精度,另一方面也再次证明,当孔隙率较小,固体骨架刚度大于孔隙流体体积模量时,饱和多孔介质可以作为单相介质处理.

5.3 基岩上覆两层不同饱和多孔介质层

  

图7 模型二P波60°入射时半空间表面O点固相位移时程Fig.7 Time history of solid phase displacement for point O on the half-space surface excited by P waves with incident angle 60°for model two

  

图8 模型二SV波30°入射时半空间表面O点固相位移时程Fig.8 Time history of solid phase displacement for point O on the half-space surface excited by SV waves with incident angle 30°for model two

计算模型如图3(b)所示,两层不同饱和介质层的参数分别取表2中第②、③组参数,两层饱和介质层厚度h1=h2=50m,网格厚度取1m,时间步长取10−4s.图7给出了P波60°入射时半空间表面O点固体骨架的水平和竖向位移.图8为SV波30°入射时O点固体骨架的水平和竖向位移.图中精确解是根据文献[30]建立的入射平面谐波在基岩上覆成层饱和多孔地基中传播的解析解答,结合快速傅里叶变换得到的结果.从图中可以看出,本文方法得到的数值解与精确解吻合较好,说明本文方法在处理成层饱和多孔介质波动问题时同样具有良好的精度.

6 结语

本文以Biot流体饱和介质动力方程为基础,根据Snell定律,并结合不同介质交界面的边界条件以及单相弹性介质精确人工边界条件,建立了一种求解地震波斜入射下基岩上覆饱和多孔介质成层场地自由场分析的一维化时域计算方法.该方法可以方便地求出平面波斜入射下饱和多孔介质成层场地任意时刻的自由场,解决了饱和多孔介质成层场地地震波输入问题.在典型场地模型下,通过与已有解析解的对比,表明本文方法有效,且具有良好的精度.

参考文献

1 Jin X,Liao ZP.Statistical research on S-wave incident angle.Earthquake Research in China,1994,8(1)∶121-131

2 Takahiro S,Takeshi U,Ryoichi T,et al.Estimation of earthquake motion incident angle at rock site.Proc of 12th World Conference Earthquake Engineering,NZ National Society for Earthquake Engineering.Auckland,New Zealand,2000∶956

3 Wolf JP,Obernhuber P.Effects of horizontally traveling waves in soil-structure interaction.Nuclear Engineering and Design,1979,57(2)∶221-244

4 李山有,廖振鹏,周正华.大型结构地震反应数值模拟中的波动输入.地震工程与工程振动,2001,6(2)∶1-5(Li Shanyou,Liao Zhenpeng,Zhou Zhenghua.Wave motion input in numerical simulation of seismic response for large-scale structure.Earthquake Engineering and Engineering Vibration,2001,6(2)∶1-5(in Chinese))

5 林皋.地下结构地震响应的计算模型.力学学报,2017,49(3)∶528-542(Lin Gao.A computational model for seismic response analysis of underground stauctures.Chinese Journal of Theoretical and Applied Mechanics,2017,49(3)∶528-542(in Chinese))

6 章小龙,李小军,陈国兴等.黏弹性人工边界等效载荷计算的改进方法.力学学报,2016,48(5)∶1126-1135(Zhang Xiaolong,Li Xiaojun,Chen Guoxing,Zhou Zhenghua.An improved method of the calculation of equivalent nodal forces in viscous-elastic artificia boundary.Chinese Journal of Theoretical and Applied Mechanics,2016,48(5)∶1126-1135(in Chinese))

7 Thomson WT.Transmission of elastic waves through a stratifie solid medium.Journal of Applied Physics,1950,21(1)∶89-93

8 Kausel E,Röesset JM.Stiffness matrices for layered soils.Bulletin of the Seismological Society of America,1981,71(6)∶1743-1761

9 尤红兵,赵凤新,荣棉水.地震波斜入射时水平层状场地的非线性地震反应.岩土工程学报,2009,31(2)∶234-240(You Hongbing,Zhao Fengxin,Rong Mianshui.Nonlinear seismic response of horizontal layered site due to inclined wave.Chinese Journal of Geotechnical Engineering,2009,31(2)∶234-240(in Chinese))

10 王笃国,赵成刚.地震波斜入射时二维成层介质自由场求解的等效线性化方法.岩土工程学报,2016,38(3)∶554-561(Wang Duguo,Zhao Chenggang.Two-dimensional equivalent linear seismic analysis of free fiel in layered half-space due to oblique incidence.Chinese Journal of Geotechnical Engineering,2016,38(3)∶554-561(in Chinese))

11 廖河山,陈清军.SH波倾斜入射时土层的非线性响应分析.同济大学学报 (自然科学版),1994,22(4)∶517-522(Liao Heshan,Chen Qingjun.Nonlinear responses of layered soils to obliquely incident SH waves.Journal of Tongji University(Natural Science),1994,22(4)∶517-522(in Chinese))

12 李山有,王学良,周正华.地震波斜入射情形下水平成层半空间自由场的时域计算.吉林大学学报(地球科学版),2003,33(3)∶372-375(Li Shanyou,Wang Xueliang,Zhou Zhenghua.The time-step numerical simulation of free fiel motion of layered half-space for inclined seismic waves.Journal of Jilin University(Earth Science Edition),2003,33(3)∶372-375(in Chinese))

13 刘晶波,王艳.成层半空间出平面自由波场的一维化时域算法.力学学报,2006,38(2)∶219-225(Liu Jingbo,Wang Yan.A 1-D timedomain method for 2-D wave motion in elastic layered half-space by antiplane wave oblique incidence.Chinese Journal of Theoretical and Applied Mechanics,2006,38(2)∶219-225(in Chinese))

14 刘晶波,王艳.成层介质中平面内自由波场的一维化时域算法.工程力学,2007,24(7)∶16-22(Liu Jingbo,Wang Yan.1D timedomain method for in-plane wave motion of free fielin layered media.Engineering Mechanics,2007,24(7)∶16-22(in Chinese))

15 Liu Jingbo,Wang Yan.A 1D time-domain method for in-plane wave motions in a layered half-space.Acta Mechanica Sinica,2007,23(6)∶673-680

16 赵密,杜修力,刘晶波等.P-SV波斜入射时成层半空间自由场的时域算法.地震工程学报,2013,35(1)∶84-90(Zhao Mi,Du Xiuli,Liu Jingbo,et al.Time-domain method for free fiel in layered half space under P-SV Waves of oblique incidence.China Earthquake Engineering Journal,2013,35(1)∶84-90(in Chinese))

17 卓卫东,高智能,谷音.P-SV波斜入射时有阻尼成层介质自由波场的一维化时域算法.水利与建筑工程学报,2016,14(6)∶18-24,34(Zhuo Weidong,Gao Zhineng,Gu Yin.A 1D time-domain method for in-plane wave motion of free fiel in layered media with damping under obliquely incident P-SV waves.Journal of Water Resources and Architectural Engineering,2016,14(6)∶18-24,34(in Chinese))

18 Zhao M,Yin HQ,Du XL,et al.1D finit element artificia boundary method for layered half space site response from obliquely incident earthquake.Earthquakes and Structures,2015,9(1)∶173-194

19 Biot MA.Theory of propagation of elastic wave in fluid-saturate porous soil.The Journal of the Acoustical Society of America,1956,28(2)∶168-178

20 Schanz M.Poroelastodynamics∶linear models,analytical solutions,and numerical methods.Applied Mechanics Reviews,2009,62(3)∶030803-1-15

21 王子辉.饱和两相与单相土互层场地中地铁车站地震反应分析.[博士论文].北京∶北京交通大学,2008(Wang Zihui.Seismic response analysis of subway station in saturated phase and single phase soil.[PhD Thesis].Beijing∶Beijing Jiaotong University,2008(in Chinese))

22 赵成刚,王进廷,史培新等.流体饱和两相多孔介质动力反应分析的显式有限元法.岩土工程学报,2001,23(2)∶178-182(Zhao Chenggang,Wang Jinting,Shi Peixin,et al.Dynamic analysis of fluid-saturate porousmediabyusingexplicitfinit elementmethod.ChineseJournalofGeotechnicalEngineering,2001,23(2)∶178-182(in Chinese))

23 Biot MA,Willis DG.The elastic coefficients of the theory of consolidation.Journal of Applied Mechanics,1957,15(2)∶594-601

24 周正华,周扣华.有阻尼振动方程常用显式积分格式稳定性分析.地震工程与工程振动,2001,21(3)∶22-28(Zhou Zhenghua,Zhou Kouhua.Stability analysis of an integral method for damped vibration equation.Earthquake Engineering and Engineering Vibration,2001,21(3)∶22-28(in Chinese))

25 王进廷,张楚汉,金峰.有阻尼动力方程显式积分方法的精度研究.工程力学,2006,23(3)∶1-5(Wang Jinting,Zhang Chuhan,Jin Feng.Ontheaccuracyofseveralexplicitintegrationschemesfordynamicequationwithdamping.EngineeringMechanics,2006,23(3)∶1-5(in Chinese))

26 李亮,杜俢力,李立云等.两相介质波动问题显式有限元方法稳定性研究.西北地震学报,2011,33(3)∶218-222,227(Li Liang,Du Xiuli,Li Liyun,et al.Study on stability of explicit finit element method for wave motion of fluid-saturate porous media.Northwestern Seismological Journal,2011,33(3)∶218-222,227(in Chinese))

27 Lin CH,Lee VW,Trifunac MD.The reflectio of plane waves in a poroelastic half-space saturated with inviscid fluid Soil Dynamics and Earthquake Engineering,2005.25(13)∶205-223

28 傅淑芳,刘宝诚.地震学教程.北京∶地震出版社,1991∶73-84(Fu Shufang,Liu Baocheng.Seismology Tutorial.Beijing∶Seismic Press,1991∶73-84(in Chinese))

29 王艳.非一致地震动场数值方法研究及在结构动力分析中的应用.[博士论文].北京∶清华大学,2007(Wang Yan.Research on the numerical method for asynchronous seismic wave motions and its application in dynamic analysis of atructures.[PhD Thesis].Beijing∶Tsinghua University,2007(in Chinese))

30 李伟华,赵成刚.地下水位变化对地震地面运动的影响.地震学报,2015,37(3)∶482-492(Li Weihua,Zhao Chenggang.Effects of the groundwater level variation on earthquake ground motions.Acta Seismologica Sinica,2015,37(3)∶482-492(in Chinese))

 
李伟华,夏佩林,张奎,赵成刚
《力学学报》 2018年第02期
《力学学报》2018年第02期文献

服务严谨可靠 7×14小时在线支持 支持宝特邀商家 不满意退款

本站非杂志社官网,上千家国家级期刊、省级期刊、北大核心、南大核心、专业的职称论文发表网站。
职称论文发表、杂志论文发表、期刊征稿、期刊投稿,论文发表指导正规机构。是您首选最可靠,最快速的期刊论文发表网站。
免责声明:本网站部分资源、信息来源于网络,完全免费共享,仅供学习和研究使用,版权和著作权归原作者所有
如有不愿意被转载的情况,请通知我们删除已转载的信息