线性加速度显式积分格式的海岸两相介质场地地震反应分析.pdfVIP

第38卷第2期烟台大学学报（自然科学与工程版Vol.38No.2

2025年4月JournalofYantaiUniversity(NaturalScienceandErngineeringEdition)Apr.2025

文章编号：1004-8820(2025)02-0125-09doi:10.13951/ki.39-1213/n.240903

线性加速度显式积分格式的海岸两相介质场地

地震反应分析

陈学良，张波？，李铁飞！

（1.中国地震局地球物理研究所，北京100081；2.河北建筑工程学院土木工程学院，河北张家口

075000）

摘要：为了建立覆水海岸两相介质地震反应分析的全显式有限元法，以饱和两相介质土

力学模型为基础，对近场波动问题的时域显式数值算法进行研究。首先，采用伽辽金加权

余量法对一维压缩波波动的动力方程实现空间离散；然后，分别应用中心差分和Newmark

线性加速度法求解位移和速度，而加速度是由速度积分法求得，进而得到该动力方程自起

步的时域逐步积分公式；最后，引入透射人工边界，得到饱和两相介质近场波动问题求解

的全显式时域有限元法。基于时域递推公式的积分逼近算子矩阵，研究了该算法的稳定

性条件。将该算法的数值解与相应的解析解进行比较,验证了算法的正确性和有效性

最后，将此方法应用于覆水海岸两相介质场地土层的地震反应分析中，研究了海水对其地

震响应的影响。该方法具有计算量小、效率高、精度高的优势，可为结构动力方程和波动

方程的求解提供借鉴。

关键词：饱和两相介质；土力学模型；显式有限元；透射人工边界；线性加速度法

中图分类号：0357.3；TU435文献标志码：A

饱和两相介质中地震波传播及其工程场地在地震波作用下的动力反应，受到国内外广大学者的高度关

注,该课题在土动力学、岩土工程、地震工程和地球物理等领域具有重要意义[1-3]。BIOT建立了描述饱和两

相介质波动的动力方程，研究了P,、P2和SV三种体波的传播特性，奠定了饱和多孔介质波传播理论研究的

基础[4-5]。此后,许多学者对饱和两相介质开展了相关研究，提出了各种不同的理论模型，如Zienkiewicz模

型[6]、槽路模型[7]、土力学模型[8]以及基于混合物理论的数学模型[9]。尤其是我国学者门福录提出的土力

学模型[7],因其参数具有明确的物理意义且易于测定,适用于模拟颗粒接触的两相介质。因此，本文以土力

学模型为基础，研究海岸饱和两相介质场地土层的波动特性。

饱和两相介质波动方程中固相与液相耦合、时间与空间耦合，力学性质比较复杂，难以求得闭合条件下

数值法求解主要以有限元为主，其

的解析解，这使得采用数值法求解两相介质的动力响应成为必然[10-1。娄

动力响应计算的准确性关键在于时域数值积分格式。根据求解过程中是否需要解耦联的方程组，时域数值

积分法分为隐式法【12】、隐-显式法[13]和显式法[14-20]。其中显式法因其计算量小、速度快的优点而备受关注，

研究者们多年来致力于探索不同的显式时域积分格式，以求高效解决饱和两相介质的动力响应问题。赵成

收稿日期：2024-09-18

基金项目：基本科研业务专项（DQJB24Z01，DQJB24R31）；国家重点研发计划（2023YFC3007305，2019YFC1509403）；国家

自然科学基金资助项目51978633）；河北省教育厅科学研究项目（QN2025419）。

通信作者：陈学良（chenxueliang007@126.com），研究员，博士生导师，博士，烟台大学结构工程专业1994级校友，主要研

究方向为工程地震、强地震动数值模拟和复杂场地非线性地震响应研究

126烟台自然科学与程版）