针对粒子增强材料的大规模数值模拟问题, 将局部Eshelby 矩阵的概念引入到本征应变边界积分方程计算模型中, 以解决粒子间的相互作用问题. 局部Eshelby 矩阵可以看作Eshelby 张量和等效夹杂物的概念在数值方面的一种拓广. 以全空间边界元子域法为参照, 利用计算模型对无限域中的若干椭球粒子进行了三维应力分析. 数值算例不仅验证了模型的正确性和方法的可行性, 也表现出较高的计算效率, 说明该计算模型和方法具有对粒子增强材料进行大规模数值分析的能力.
Aiming at large scale numerical simulation of particle reinforced materials, a concept of local Eshelby matrix is introduced into a computational model of the eigenstrain boundary integral equation to solve the problem of interactions among particles. The local Eshelby matrix can be considered as an extension of Eshelby tensor and an equivalent inclusion in a numerical form. Taking the sub-domain boundary element method as the control, three-dimensional stress analyses are carried out for some ellipsoidal particles in infinite media with the proposed computational model. Numerical examples verify correctness, feasibility and high efficiency of the present model with the corresponding solution procedure, showing potential of solving large scale numerical simulations for particle reinforced materials.
[1] Eshelby J D. The determination of the elastic field of an ellipsoidal inclusion and related problems [J]. Proceedings of the Royal Society of London Series A: Mathematical and Physical
Sciences, 1957, A241(1226): 376-396.
[2] Mura T, Shodja H M, Hirose Y. Inclusion problems (part 2) [J]. Applied Mechanics Review,1996, 49(10): S118-S127.
[3] Kiris A, Inan E. Eshelby tensors for a spherical inclusion in microelongated elastic fields [J]. International Journal of Engineering Science, 2005, 43: 49-58.
[4] Mercier S, Jacques N, Molinari A. Validation of an interaction law for the Eshelby inclusion problem in elasto-viscoplasticity [J]. International Journal of Solids and Structures, 2005, 42:
1923-1941.
[5] Shen L X, Yi S. An effective inclusion model for effective moluli of heterogeneous materials with ellipsoidal inhomogeneities [J]. International Journal of Solids and Structures, 2001, 38:
5789-5805.
[6] Pan E. Elastic or piezoelastic fields around a quantum dot: fully coupled or semicoupled model? [J] Journal of Applied Physics, 2002, 91(6): 3785-3796.
[7] Ma H, Deng H L. Nondestructive determination of welding residual stresses by boundary element method [J]. Advances in Engineering Software, 1998, 29: 89-95.
[8] Ma H, Guo Z, Dhanasekar M, et al. Efficient solution of multiple cracks in great number using eigen COD boundary integral equations with iteration procedure [J]. Engineering Analysis
with Boundary Elements, 2013, 37(3): 487-500.
[9] Kakavas P A, Kontoni D N. Numerical investigation of the stress field of particulate reinforced polymeric composites subjected to tension [J]. International Journal for Numerical Methods in Engineering, 2006, 65: 1145-1164.
[10] Lee J, Choi S, Mal A. Stress analysis of an unbounded elastic solid with orthotropic inclusions and voids using a new integral equation technique [J]. International Journal of Solids and
Structures, 2001, 38: 2789-2802.
[11] Dong C Y, Cheung Y K, Lo S H. A regularized domain integral formulation for inclusion problems of various shapes by equivalent inclusion method [J]. Computer Methods in Applied
Mechanics and Engineering, 2002, 191(31): 3411-3421.
[12] Nakasone Y, Nishiyama H, Nojiri T. Numerical equivalent inclusion method: a new computational method for analyzing stress fields in and around inclusions of various shapes [J].
Materials Science and Engineering, 2000, A285: 229-238.
[13] Greengard L F, Rokhlin V. A fast algorithm for particle simulations [J]. Journal of Computational Physics, 1987, 73: 325-48.
[14] Liu Y J, Nishimura N, Tanahashi T, et al. A fast boundary element method for the analysis of fiber-reinforced composites based on a rigid-inclusion model [J]. ASME Journal of Applied
Mechanics, 2005, 72: 115-128.
[15] Ma H, Yan C, Qin Q H. Eigenstrain formulation of boundary integral equations for modeling particle-reinforced composites [J]. Engineering Analysis with Boundary Elements, 2009, 33(3): 410-419.
[16] Ma H, Xia L W, Qin Q H. Computational model for short-fiber composites with eigen-strain formulation of boundary integral equations [J]. Applied Mathematics and Mechanics, 2008,
29(6): 757-767.
[17] Ma H, Fang J B, Qin Q H. Simulation of ellipsoidal particle-reinforced materials with eigenstrain formulation of 3D BIE [J]. Advances in Engineering Software, 2011, 42(10): 750-759.
[18] Chen Y Z. Boundary integral equation method for two dissimilar elastic inclusions in an infinite plate [J]. Engineering Analysis with Boundary Elements, 2012, 36(1): 137-146.