Computational micromechanics and homogenization require the solution of the mechanical equilibrium of a periodic cell that comprises a (generally complex) microstructure. Techniques that apply the Fast Fourier Transform have attracted much attention as they outperform other methods in terms of speed and memory footprint. Moreover, the Fast Fourier Transform is a natural companion of pixel-based digital images which often serve as input. In its original form, one of the biggest challenges for the method is the treatment of (geometrically) non-linear problems, partially due to the need for a uniform linear reference problem. In a geometrically linear setting, the problem has recently been treated in a variational form resulting in an unconditionally stable scheme that combines Newton iterations with an iterative linear solver, and therefore exhibits robust and quadratic convergence behavior. Through this approach, well-known key ingredients were recovered in terms of discretization, numerical quadrature, consistent linearization of the material model, and the iterative solution of the resulting linear system. As a result, the extension to finite strains, using arbitrary constitutive models, is at hand. Because of the application of the Fast Fourier Transform, the implementation is substantially easier than that of other (Finite Element) methods. Both claims are demonstrated in this paper and substantiated with a simple code in Python of just 59 lines (without comments). The aim is to render the method transparent and accessible, whereby researchers that are new to this method should be able to implement it efficiently. The potential of this method is demonstrated using two examples, each with a different material model

No abstract

Sliding at a quasi-statically loaded frictional interface can occur via macroscopic slip events, which nucleate locally before propagating as rupture fronts very similar to fracture. We introduce a novel microscopic model of a frictional interface that includes asperity-level disorder, elastic interaction between local slip events, and inertia. For a perfectly flat and homogeneously loaded interface, we find that slip is nucleated by avalanches of asperity detachments of extension larger than a critical radius Ac governed by a Griffith criterion. We find that after slip, the density of asperities at a local distance to yielding xσ presents a pseudo-gap P (xσ) ∼ (xσ) θ , where θ is a non-universal exponent that depends on the statistics of the disorder. This result makes a link between friction and the plasticity of amorphous materials where a pseudo-gap is also present. For friction, we find that a consequence is that stick-slip is an extremely slowly decaying finite size effect, while the slip nucleation radius Ac diverges as a θ-dependent power law of the system size. We discuss how these predictions can be tested experimentally. Significance statementUnderstanding how slip at a frictional interface initiates is important for a range of problems including earthquake prediction and precision engineering. The force needed to start sliding a solid object over a flat surface is classically described by a 'static friction coefficient': a constant established by measurements. It was recently questioned if such constant exists, as it was shown to be poorly reproducible. We provide a model supporting that it is stochastic even for very large system sizes: sliding is nucleated when, by chance, an avalanche of microscopic detachments reaches a critical radius, beyond which slip becomes unstable and propagates along the interface. It leads to testable predictions on key observables characterising the stability of the interface.

Quasi-localised modes appear in the vibrational spectrum of amorphous solids at low-frequency. Though never formalised, these modes are believed to have a close relationship with other important local excitations, including shear transformations and two-level systems. We provide a theory for their frequency density, DL(ω) ∼ ω α , that establishes this link for systems at zero temperature under quasi-static loading. It predicts two regimes depending on the density of shear transformations P (x) ∼ x θ (with x the additional stress needed to trigger a shear transformation). If θ > 1/4, α = 4 and a finite fraction of quasi-localised modes form shear transformations, whose amplitudes vanish at low frequencies. If θ < 1/4, α = 3 + 4θ and all quasi-localised modes form shear transformations with a finite amplitude at vanishing frequencies. We confirm our predictions numerically. arXiv:1806.01561v2 [cond-mat.soft]