Accelerated computational micromechanics

We present an approach to solving problems in micromechanics that is amenable to massively parallel calculations through the use of graphical processing units and other accelerators. The problems lead to nonlinear differential equations that are typically second order in space and first order in time. This combination of nonlinearity and nonlocality makes such problems difficult to solve in parallel. However, this combination is a result of collapsing nonlocal, but linear and universal physical laws (kinematic compatibility, balance laws), and nonlinear but local constitutive relations. We propose an operator-splitting scheme inspired by this structure. The governing equations are formulated as (incremental) variational problems, the differential constraints like compatibility are introduced using an augmented Lagrangian, and the resulting incremental variational principle is solved by the alternating direction method of multipliers. The resulting algorithm has a natural connection to physical principles, and also enables massively parallel implementation on structured grids. We present this method and use it to study two examples. The first concerns the long wavelength instability of finite elasticity, and allows us to verify the approach against previous numerical simulations. We also use this example to study convergence and parallel performance. The second example concerns microstructure evolution in liquid crystal elastomers and provides new insights into some counter-intuitive properties of these materials. We use this example to validate the model and the approach against experimental observations.

[1]  Matti Schneider,et al.  An efficient solution scheme for small-strain crystal-elasto-viscoplasticity in a dual framework , 2020 .

[2]  Sachin S. Gautam,et al.  GPU-warp based finite element matrices generation and assembly using coloring method , 2019, J. Comput. Des. Eng..

[3]  B. Bourdin,et al.  The Variational Approach to Fracture , 2008 .

[4]  Dennis Merkert,et al.  FFT-based homogenization on periodic anisotropic translation invariant spaces , 2017, 1701.04685.

[5]  Yang Chen,et al.  A FFT solver for variational phase-field modeling of brittle fracture , 2019, Computer Methods in Applied Mechanics and Engineering.

[6]  Radhi Abdelmoula,et al.  A damage model for crack prediction in brittle and quasi-brittle materials solved by the FFT method , 2012, International Journal of Fracture.

[7]  Stephen P. Boyd,et al.  Distributed Optimization and Statistical Learning via the Alternating Direction Method of Multipliers , 2011, Found. Trends Mach. Learn..

[8]  G. Bonnet,et al.  A polarization‐based FFT iterative scheme for computing the effective properties of elastic composites with arbitrary contrast , 2012 .

[9]  Roland Glowinski,et al.  Variational methods for the numerical solution of nonlinear elliptic problems , 2015, CBMS-NSF regional conference series in applied mathematics.

[10]  R. Ma,et al.  FFT-based homogenization of hypoelastic plasticity at finite strains , 2019, Computer Methods in Applied Mechanics and Engineering.

[11]  K. Washizu Variational Methods in Elasticity and Plasticity , 1982 .

[12]  Jaroslav Vondrejc,et al.  An FFT-based Galerkin method for homogenization of periodic media , 2013, Comput. Math. Appl..

[13]  J. Yvonnet Computational Homogenization of Heterogeneous Materials with Finite Elements , 2020, Solid Mechanics and Its Applications.

[14]  J. Segurado,et al.  Computational Homogenization of Polycrystals , 2018, 1804.02538.

[15]  M. Ortiz,et al.  The variational formulation of viscoplastic constitutive updates , 1999 .

[16]  P. Sheng,et al.  Introduction to Liquid Crystals , 1976 .

[17]  Y. Huo,et al.  Continuum mechanical modeling of liquid crystal elastomers as dissipative ordered solids , 2019, Journal of the Mechanics and Physics of Solids.

[18]  R. Glowinski,et al.  Augmented Lagrangian and Operator-Splitting Methods in Nonlinear Mechanics , 1987 .

[19]  Nicolas Triantafyllidis,et al.  Failure Surfaces for Finitely Strained Two-Phase Periodic Solids Under General In-Plane Loading , 2006 .

[20]  K. Bhattacharya,et al.  Probing the in-plane liquid-like behavior of liquid crystal elastomers , 2021, Science Advances.

[21]  Roland Glowinski,et al.  ADMM and Non-convex Variational Problems , 2016 .

[22]  A. Rollett,et al.  Validation of a numerical method based on Fast Fourier Transforms for heterogeneous thermoelastic materials by comparison with analytical solutions , 2014 .

[23]  D. Kochmann,et al.  Predicting the effective response of bulk polycrystalline ferroelectric ceramics via improved spectral phase field methods , 2017 .

[24]  A. DeSimone,et al.  Soft elastic response of stretched sheets of nematic elastomers: a numerical study , 2002 .

[25]  K. Bertoldi,et al.  Mechanics of deformation-triggered pattern transformations and superelastic behavior in periodic elastomeric structures , 2008 .

[26]  K. Bertoldi,et al.  Pattern transformation triggered by deformation. , 2007, Physical review letters.

[27]  H. Finkelmann,et al.  Nematic liquid single crystal elastomers , 1991 .

[28]  K. Bhattacharya,et al.  Supersoft elasticity in polydomain nematic elastomers. , 2009, Physical review letters.

[29]  M. Geers,et al.  Finite strain FFT-based non-linear solvers made simple , 2016, 1603.08893.

[30]  Marko Knezevic,et al.  A high-performance computational framework for fast crystal plasticity simulations , 2014 .

[31]  J. Rice Inelastic constitutive relations for solids: An internal-variable theory and its application to metal plasticity , 1971 .

[32]  H. Finkelmann,et al.  Strain‐induced director reorientation in nematic liquid single crystal elastomers , 1995 .

[33]  Tom Goldstein,et al.  The Split Bregman Method for L1-Regularized Problems , 2009, SIAM J. Imaging Sci..

[34]  W. Ludwig,et al.  Imposing equilibrium on measured 3-D stress fields using Hodge decomposition and FFT-based optimization , 2021, 2105.01612.

[35]  H. Moulinec,et al.  A fast numerical method for computing the linear and nonlinear mechanical properties of composites , 1994 .

[36]  R. J. Atkin,et al.  An introduction to the theory of elasticity , 1981 .

[37]  Soft and nonsoft structural transitions in disordered nematic networks , 2000, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[38]  Laurent Capolungo,et al.  A FFT-based formulation for discrete dislocation dynamics in heterogeneous media , 2018, J. Comput. Phys..

[39]  Antonio DeSimone,et al.  Macroscopic Response of¶Nematic Elastomers via Relaxation of¶a Class of SO(3)-Invariant Energies , 2002 .

[40]  Wotao Yin,et al.  Bregman Iterative Algorithms for (cid:2) 1 -Minimization with Applications to Compressed Sensing ∗ , 2008 .

[41]  Roland Glowinski,et al.  Numerical Solution of Problems in Incompressible Finite Elasticity by Augmented Lagrangian Methods. I. Two-Dimensional and Axisymmetric Problems , 1982 .

[42]  Roland Glowinski,et al.  Numerical Solution of Problems in Incompressible Finite Elasticity by Augmented Lagrangian Methods II. Three-Dimensional Problems , 1984 .

[43]  Stephen Lee,et al.  Exascale Computing in the United States , 2019, Computing in Science & Engineering.

[44]  Ted Belytschko,et al.  Multiscale Methods in Computational Mechanics , 2010 .

[45]  L. Gélébart,et al.  Non-linear extension of FFT-based methods accelerated by conjugate gradients to evaluate the mechanical behavior of composite materials , 2013 .

[46]  Graeme W. Milton,et al.  A fast numerical scheme for computing the response of composites using grid refinement , 1999 .

[47]  Yongmei M Jin,et al.  Three-dimensional phase field model of proper martensitic transformation , 2001 .

[48]  Jacob Fish,et al.  Multiscale Methods: Bridging the Scales in Science and Engineering , 2009 .

[49]  J. Segurado,et al.  An algorithm for stress and mixed control in Galerkin‐based FFT homogenization , 2019, International Journal for Numerical Methods in Engineering.

[50]  H. Moulinec,et al.  Convergence of iterative methods based on Neumann series for composite materials: Theory and practice , 2017, 1711.05880.

[51]  R. Lebensohn,et al.  A fast Fourier transform-based mesoscale field dislocation mechanics study of grain size effects and reversible plasticity in polycrystals , 2020, Journal of the Mechanics and Physics of Solids.

[52]  Badel L. Mbanga,et al.  Modeling elastic instabilities in nematic elastomers. , 2009, Physical review. E, Statistical, nonlinear, and soft matter physics.

[53]  Alexander Mielke,et al.  A Mathematical Framework for Generalized Standard Materials in the Rate-Independent Case , 2006 .

[54]  Alexander Mielke,et al.  An Approach to Nonlinear Viscoelasticity via Metric Gradient Flows , 2013, SIAM J. Math. Anal..

[55]  Ricardo A. Lebensohn,et al.  Numerical implementation of non-local polycrystal plasticity using fast Fourier transforms , 2016 .

[56]  Yuki Ueda,et al.  Numerical computations of split Bregman method for fourth order total variation flow , 2019, J. Comput. Phys..

[57]  A Modified Split Bregman Algorithm for Computing Microstructures Through Young Measures , 2019, Multiscale Model. Simul..

[58]  M. Schneider,et al.  On Quasi‐Newton methods in fast Fourier transform‐based micromechanics , 2019, International Journal for Numerical Methods in Engineering.

[59]  Marko Knezevic,et al.  Three orders of magnitude improved efficiency with high‐performance spectral crystal plasticity on GPU platforms , 2014 .

[60]  Hervé Moulinec,et al.  A Computational Method Based on Augmented Lagrangians and Fast Fourier Transforms for Composites with High Contrast , 2000 .

[61]  Iam-Choon Khoo,et al.  Introduction to Liquid Crystals , 2006 .

[62]  K. Bhattacharya,et al.  Elasticity of polydomain liquid crystal elastomers , 2009, 0911.3513.

[63]  Claude Fressengeas,et al.  A numerical spectral approach for solving elasto-static field dislocation and g-disclination mechanics , 2014 .

[64]  Adnan Eghtesad,et al.  Spectral database constitutive representation within a spectral micromechanical solver for computationally efficient polycrystal plasticity modelling , 2018 .

[65]  Hervé Moulinec,et al.  Comparison of three accelerated FFT‐based schemes for computing the mechanical response of composite materials , 2014 .

[66]  E. Terentjev,et al.  Semisoft elastic response of nematic elastomers to complex deformations. , 2008, Physical review. E, Statistical, nonlinear, and soft matter physics.

[67]  Luc Dormieux,et al.  FFT-based methods for the mechanics of composites: A general variational framework , 2010 .

[68]  Jaroslav Vondrejc,et al.  A comparative study on low-memory iterative solvers for FFT-based homogenization of periodic media , 2015, J. Comput. Phys..

[69]  L. Dormieux,et al.  Combining Galerkin approximation techniques with the principle of Hashin and Shtrikman to derive a new FFT-based numerical method for the homogenization of composites , 2012 .

[70]  K. Urayama,et al.  Polydomain−Monodomain Transition of Randomly Disordered Nematic Elastomers with Different Cross-Linking Histories , 2009 .

[71]  Elastic effects in disordered nematic networks. , 1999, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[72]  M. Schneider,et al.  Efficient fixed point and Newton–Krylov solvers for FFT-based homogenization of elasticity at large deformations , 2014 .

[73]  Renald Brenner,et al.  Computational approach for composite materials with coupled constitutive laws , 2010 .

[74]  N. Triantafyllidis,et al.  Homogenization of nonlinearly elastic materials, microscopic bifurcation and macroscopic loss of rank-one convexity , 1993 .

[75]  K. Bhattacharya,et al.  Microstructure-enabled control of wrinkling in nematic elastomer sheets , 2016, 1611.08621.

[76]  B. He,et al.  Alternating Direction Method with Self-Adaptive Penalty Parameters for Monotone Variational Inequalities , 2000 .

[77]  Graeme W. Milton,et al.  Substitution of subspace collections with nonorthogonal subspaces to accelerate Fast Fourier Transform methods applied to conducting composites , 2020, ArXiv.

[78]  M. Geers,et al.  A finite element perspective on nonlinear FFT‐based micromechanical simulations , 2016, 1601.05970.

[79]  Jan Novák,et al.  Accelerating a FFT-based solver for numerical homogenization of periodic media by conjugate gradients , 2010, J. Comput. Phys..

[80]  Wolfgang Paul,et al.  GPU accelerated Monte Carlo simulation of the 2D and 3D Ising model , 2009, J. Comput. Phys..

[81]  Hervé Moulinec,et al.  A numerical method for computing the overall response of nonlinear composites with complex microstructure , 1998, ArXiv.

[82]  T. Nguyen,et al.  Viscoelasticity of the polydomain-monodomain transition in main-chain liquid crystal elastomers , 2016 .

[83]  J. Schröder A numerical two-scale homogenization scheme: the FE 2 -method , 2014 .

[84]  K. Bhattacharya,et al.  Effective Behavior of Nematic Elastomer Membranes , 2015, Archive for Rational Mechanics and Analysis.

[85]  K. Bhattacharya,et al.  Interplay of martensitic phase transformation and plastic slip in polycrystals , 2013 .

[86]  E. Cerreta,et al.  Modeling void growth in polycrystalline materials , 2013 .

[87]  Philip Eisenlohr,et al.  An elasto-viscoplastic formulation based on fast Fourier transforms for the prediction of micromechanical fields in polycrystalline materials , 2012 .

[88]  M. Schneider An FFT‐based method for computing weighted minimal surfaces in microstructures with applications to the computational homogenization of brittle fracture , 2019, International Journal for Numerical Methods in Engineering.

[89]  P. Bladon,et al.  Transitions and instabilities in liquid crystal elastomers. , 1993, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[90]  Javier Segurado,et al.  DBFFT: A displacement based FFT approach for non-linear homogenization of the mechanical behavior , 2019, International Journal of Engineering Science.

[91]  A. DeSimone,et al.  Ogden-type energies for nematic elastomers , 2012 .

[92]  Dimitri P. Bertsekas,et al.  On the Douglas—Rachford splitting method and the proximal point algorithm for maximal monotone operators , 1992, Math. Program..

[93]  M. Schneider On the Barzilai‐Borwein basic scheme in FFT‐based computational homogenization , 2019, International Journal for Numerical Methods in Engineering.

[94]  Kaushik Bhattacharya,et al.  A model problem concerning recoverable strains of shape-memory polycrystals , 2005, Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[95]  Diego Rossinelli,et al.  High performance CPU/GPU multiresolution Poisson solver , 2013, PARCO.

[96]  Martin Diehl,et al.  Numerically robust spectral methods for crystal plasticity simulations of heterogeneous materials , 2013 .

[97]  E. M. Terentjev,et al.  Liquid Crystal Elastomers , 2003 .