The paper proposes a mathematical model aimed at calculating the effective elastic moduli of a micro-inhomogeneous two-component isotropic composite material, which components are connected randomly depending on the level of their relative volumetric contents. A stochastic equation is formulated for the connectivity parameter of the constituent components, according to which, with an increase in the volumetric content of the filler, individual inclusions build the structures of the matrix mixture in the form of interpenetrating frameworks, and then turn into a new binding matrix with individual inclusions from the material of the rest of the old matrix. The algorithm for the numerical solution of this stochastic differential equation is constructed in accordance with the Euler-Maruyama method. For each implementation of this algorithm, the corresponding stochastic trajectories are constructed for the random connectivity function of the constituent components of the composite material. A variant of the method aimed at calculating the mathematical expectation of a random connectivity function of the constituent components has been developed and the corresponding differential equation has been obtained for it. It is shown that the numerical solution of this equation and the average value of the production factor function calculated for all realizations of stochastic trajectories give close numerical values. New macroscopic constitutive relations are found for microinhomogeneous materials with a variable microstructure and their effective elastic moduli are calculated. It is noted that the formulas for these effective elastic moduli generalize the known results for isotropic composite materials. The values of the effective elastic moduli, constructed according to the expressions obtained in the paper, lie within the Khashin-Shtrikman range for the lower and upper bounds of the effective elastic moduli of the composite materials. The numerical analysis of the developed models showed a good agreement with the known experimental data.
Numerical solution of Painlevand#233; II differential equation
