岩土力学 ›› 2026, Vol. 47 ›› Issue (10): 3662-3670.doi: 10.16285/j.rsm.2025.1068CSTR: 32223.14.j.rsm.2025.1068

• 数值分析 • 上一篇    

时域边界元法求解多孔弹性介质瞬态流固耦合问题

赵峰1,罗俊成1,曾绮华2,雷卫东2,李宏军3,林坚1   

  1. 1.广州地铁设计研究院股份有限公司,广东 广州 510410;2.哈尔滨工业大学(深圳) 智能土木与海洋工程学院,广东 深圳 518055; 3.河北农业大学 河北省乡村建筑结构安全与低碳建造重点实验室,河北 保定 071001
  • 收稿日期:2025-10-08 接受日期:2026-05-27 出版日期:2026-10-10 发布日期:2026-10-09
  • 通讯作者: 李宏军,男,1988年生,硕士,讲师,主要从事计算力学、岩土工程数值方法研究。E-mail: chjlhj@hebau.edu.cn
  • 作者简介:赵峰,男,1987年生,硕士,高级工程师,主要从事地下空间结构设计与研究工作。E-mail: zhaofeng3@dtsjy.com
  • 基金资助:
    国家自然科学基金(No. 12471497);广州地铁设计研究院股份有限公司资助项目(No.KY-2022-014)。

Time domain boundary element method for solving transient fluid-structure coupling problems in porous elastic media

ZHAO Feng1, LUO Jun-cheng1, ZENG Qi-hua2, LEI Wei-dong2, LI Hong-jun3, LIN Jian1   

  1. 1. Guangzhou Metro Design & Research Institute Co., Ltd., Guangzhou, Guangdong 510410, China; 2. School of Intelligent Civil and Ocean Engineering, Harbin Institute of Technology Shenzhen, Shenzhen, Guangdong 518055, China; 3. Hebei Key Laboratory of Structural Safety and Low-Carbon Construction for Rural Buildings, Hebei Agricultural University, Baoding, Hebei 071001, China
  • Received:2025-10-08 Accepted:2026-05-27 Online:2026-10-10 Published:2026-10-09
  • Supported by:
    This work was supported by the National Natural Science Foundation of China (12471497) and the Project of Guangzhou Metro Design and Research Institute Co., Ltd. (KY-2022-014).

摘要: 采用时域边界元法研究多孔弹性介质瞬态流固耦合问题,自动满足无穷远处边界条件,因此,只需要离散边界,特别适合具有无穷边界的问题。首先,基于多孔弹性介质偏微分控制方程推导出瞬态流固耦合问题的边界积分方程。其次,对边界积分方程进行数值处理,重点求解了含有强异性的单元影响系数的Hadamard主值,避免了间接数值解法的误差累积,使奇异性处理精度较高。最后,采用二阶形函数确定了半无限单元影响系数,组装得到边界积分方程的代数方程组形式。采用两个典型算例进行验证,将本研究所提出方法的计算结果与解析解进行了对比。结果表明,所提出的时域边界元法适用于求解多孔弹性介质瞬态流固耦合问题。

关键词: 多孔弹性介质, 瞬态流固耦合, 时域边界元法, 奇异性处理, 无穷域

Abstract: In this paper, the time-domain boundary element method is employed to investigate transient fluid-structure interaction problems in porous elastic media. The proposed method automatically satisfies the far-field boundary condition at infinity, thus only necessitating boundary discretization, which renders it particularly well-suited for addressing problems characterized by infinite boundaries. Initially, the boundary integral equation for the transient fluid-structure interaction problem is derived from the governing partial differential equations that describe porous elastic media. Subsequently, numerical discretization is applied to the boundary integral equation, with a particular focus on accurately computing the Hadamard principal value of the element influence coefficients that exhibit strong singularities. This strategy effectively circumvents error accumulation in indirect numerical schemes and improves computational precision. Following this, second-order shape functions are adopted to determine the influence coefficients of semi-infinite elements, and the algebraic system corresponding to the boundary integral equation is constructed through assembly. Two representative numerical examples are used for validation purposes, wherein the numerical results generated by the proposed method are compared with analytical solutions. The results substantiate that the developed time-domain boundary element method is both applicable and reliable for solving transient fluid-structure interaction problems in porous elastic media.

Key words: porous elastic medium, transient fluid-structure coupling, time domain boundary element method, treatment of singularities, infinite domain

中图分类号: TU 433
[1] 袁 帅,钟宏志. 无穷域问题的弱形式求积元分析[J]. , 2016, 37(4): 1187-1194.
[2] 周 晔,郑荣跃,刘干斌. 列车荷载下上覆弹性层横观各向同性饱和地基的动力响应[J]. , 2011, 32(2): 604-610.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
No Suggested Reading articles found!