Design sensitivity analysis for the homogenized elasticity tensor of a polymer filled with rubber particles

The main purpose of this work is the computational simulation of the sensitivity coefficients of the homogenized tensor for a polymer filled with rubber particles with respect to the material parameters of the constituents. The Representative Volume Element (RVE) of this composite contains a single spherical particle, and the composite components are treated as homogeneous isotropic media, resulting in an isotropic effective homogenized material. The sensitivity analysis presented in this paper is performed via the provided semi-analytical technique using the commercial FEM code ABAQUS and the symbolic computation package MAPLE. The analytical method applied for comparison uses the additional algebraic formulas derived for the homogenized tensor for a medium filled with spherical inclusions, while the FEM-based technique employs the polynomial response functions recovered from the Weighted Least-Squares Method. The homogenization technique consists of equating the strain energies for the real composite and the artificial isotropic material characterized by the effective elasticity tensor. The homogenization problem is solved using ABAQUS by the application of uniform deformations on specific outer surfaces of the composite RVE and the use of tetrahedral finite elements C3D4. The energy approach will allow for the future application of more realistic constitutive models of rubber-filled polymers such as that of Mullins and for RVEs of larger size that contain an agglomeration of rubber particles.