Abstract
A method of multi-scale computations is developed for analyzing the quasi-static mechanical behavior of elastoplastic heterogeneous media with interfacial debonding. After posing the two-scale boundary value problem, we developed the numerical algorithm for microscopic self-equlibrated problem involving discontinuous deformations. This algorithm, which accommodates the kinematical constraints on the interfaces by the penalty method, enables us to analyze the elastic-plastic behavior with debonding without any numerical instability. The representative numerical examples well illustrate the remarkable influence of microstructures with debonding interfaces on the macroscopic elastic-plastic material behavior; in particular the macroscopic yield stress.
