Limit analysis of composite materials with anisotropic microstructures: A homogenization approach

A nonlinear mathematical programming approach together with the finite element method and homogenization technique is developed to implement kinematic limit analysis for a microstructure and the macroscopic strength of a composite with anisotropic constituents can be directly calculated. By means of the homogenization theory, the classical kinematic theorem of limit analysis is generalized to incorporate the microstructure – Representative Volume Element (RVE) chosen from a periodic composite/heterogeneous material. Then, using an associated plastic flow rule, a general yield function is directly introduced into limit analysis and a purely-kinematic formulation is obtained. Based on the mathematical programming technique, the finite element model of microstructure is finally formulated as a nonlinear programming problem subject to only one equality constraint, which is solved by a direct iterative algorithm. The calculation is entirely based on a purely-kinematical velocity field without calculation of stress fields. Meanwhile, only one equality constraint is introduced into the nonlinear programming problem. So the computational cost is very modest. Both anisotropy and pressure-dependence of material yielding behavior are considered in the general form of kinematic limit analysis. The developed method provides a direct approach for determining the macroscopic strength domain of anisotropic composites and can serve as a powerful tool for microstructure design of composites.