The main aim of this study is determination of the basic probabilistic characteristics of the effective stiffness for inelastic particulate composites with spherical reinforcement and an uncertain Gaussian volume fraction of the interphase defects. This is determined using a homogenization method with a cubic single-particle representative volume element (RVE) of such a composite and the finite element method solution. A reinforcing particle is spherical, located centrally in the RVE, surrounded by the thin interphase of constant thickness, and remains in an elastic reversible regime opposite to the matrix, which is hyper-elastic. The interphase defects are represented as semi-spherical voids, which are placed on the outer surface of this particle. The interphase is modeled as hyper-elastic and isotropic, whose effective stiffness is calculated by the spatial averaging of hyper-elastic parameters of the matrix and of the defects. A constitutive relation of the matrix is recovered experimentally by its uniaxial stretch. The 3D homogenization problem solution is based upon a numerical determination of strain energy density in the given RVE under specific uniaxial and biaxial stretches as well as under shear deformations. The analytical relation of the effective composite stiffness to the input uncertain parameter is recovered via the response function method, using a polynomial basis and an optimized order. Probabilistic calculations are completed using three concurrent approaches, namely the iterative stochastic finite element method (SFEM), Monte Carlo simulation and by the semi-analytical method. Previous papers consider the composite fully elastic, which limits the applicability of the resulting effective stiffness tensor computed therein. The current study voids this assumption and defines the composite as fully hyper-elastic, thus extending applicability of this tensor to strains up to 0.25. The most important research finding is that (1) the effective stiffness tensor is sensitive to random interface defects in its hyper-elastic range, (2) its resulting randomness is not close to Gaussian, (3) the semi-analytical method is not perfectly suited to stochastic calculations in this region of strains, as opposed to the linear elastic region, and (4) that the increase in random dispersion of defects volume fraction has a much higher effect on the stochastic characteristics of this stiffness tensor than fluctuation of the strain.
Read full abstract