针对含流体粒子固体问题, 将 Eshelby 的本征应变和等效夹杂理论引入边界积分方程中, 提出含大量流体粒子固体的本征应变边界积分方程的计算模型及迭代算法. 通过全空间中 多粒子近场群的定义, 在本征应变边界积分方程中引入局部 Eshelby 矩阵, 以解决粒子间相互作用的问题, 从而保证了迭代计算的收敛性. 通过数值算例的比较, 验证了计算 模型的正确性和计算方法的可行性, 其中对于全空间单个圆形粒子的情况采用了解析解 作为比较基准, 其他情况则采用子域边界积分方程法的结果作为比较基准. 采用代表性体积单元(representative volume element, RVE)模拟了含流体粒子固体材料的整体力学性能, 对含有 1 000 多个流体粒子的 RVE 在粒子规则分布和随机分布的条件下进行了计算, 验证了本计算模型和算法的可行性和高效率.

In view of the fact that elastic solids contain fluid-filled pores, Eshelby's idea of eigenstrain and equivalent inclusion has been incorporated into the boundary integral equations (BIE), and as a result, the computational model of eigenstrain boundary integral equations and the corresponding iterative solution procedures are presented in the paper for the numerical simulation of solids with fluid-filled pores in great number. In order to guarantee the convergence of iteration sufficiently, the local Eshelby matrix has been proposed and constructed from the BIE combined with Eshelby's idea. The feasibility and effectiveness of the proposed computational model are verified in comparison with the results of the analytical solution in the case of a single circular fluid-filled pore in full space and of the subdomain BIE in other cases. The overall mechanical properties of solids are computed using a representative volume element (RVE) with more than one thousand fluid-filled pores distributed either regularly or randomly with the proposed computational model, showing the feasibility and high efficiency of the present model and the solution procedures.