Premium
A computational scheme for linear and non‐linear composites with arbitrary phase contrast
Author(s) -
Michel J. C.,
Moulinec H.,
Suquet P.
Publication year - 2001
Publication title -
international journal for numerical methods in engineering
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.421
H-Index - 168
eISSN - 1097-0207
pISSN - 0029-5981
DOI - 10.1002/nme.275
Subject(s) - linear elasticity , rate of convergence , mathematics , contrast (vision) , fourier transform , lagrange multiplier , fourier series , elasticity (physics) , saddle point , saddle , mathematical analysis , finite element method , algorithm , mathematical optimization , geometry , materials science , computer science , composite material , structural engineering , computer network , channel (broadcasting) , artificial intelligence , engineering
A numerical method making use of fast Fourier transforms has been proposed in Moulinec and Suquet (1994, 1998) to investigate the effective properties of linear and non‐linear composites. This method is based on an iterative scheme the rate of convergence of which is proportional to the contrast between the phases. Composites with high contrast (typically above 10 4 ) or infinite contrast (those containing voids or rigid inclusions or highly non‐linear materials) cannot be handled by the method. This paper presents two modified schemes. The first one is an accelerated scheme for composites with high contrast which extends to elasticity a scheme initially proposed in Eyre and Milton (1999). Its rate of convergence varies as the square root of the contrast. The second scheme, adequate for composites with infinite contrast, is based on an augmented Lagrangian method. The resulting saddle‐point problem involves three steps. The first step consists of solving a linear elastic problem, using the fast Fourier transform method. In the second step, a non‐linear problem is solved at each individual point in the volume element. The third step consists of updating the Lagrange multiplier. Applications of this scheme to rigidly reinforced and to voided composites are shown. Copyright © 2001 John Wiley & Sons, Ltd.