Determination of sensitivity coefficients of the elastic effective properties for periodic fiber-reinforced composites using the response function method

Derivation and implementation of the homogenization method including determination of sensitivity gradients of the effective elasticity tensor using combined numerical-analytical approach are addressed in this paper. This is possible thanks to an application of the numerical response function together with the effective moduli method known from classical homogenization theory. Computational procedure is implemented using 4-noded quadrilateral plane strain finite elements (program MCCEFF) and the symbolic computations system MAPLE. The sensitivity coefficients are determined on the basis of partial derivatives of the homogenized elasticity tensor calculated using the response function method with respect to all composite components’ elastic characteristics. They are further separately subjected to normalization procedure for a final comparison with each other. Such an enriched homogenization procedure is tested on the periodic fiber-reinforced two component composites; the results of computational analysis are compared to the results of the central finite difference approach applied before. Computational methodology proposed here may be further successively applied not only in the context of homogenization method but also to extend various discrete computational techniques like boundary/finite element, finite difference and volumes together with various meshless methods.