We present an open-source software framework called PERMIX for multiscale modeling and simulation of fracture in solids. The framework is an object oriented open-source effort written primarily in ...Fortran 2003 standard with Fortran/C++ interfaces to a number of other libraries such as LAMMPS, ABAQUS, LS-DYNA and GMSH. Fracture on the continuum level is modeled by the extended finite element method (XFEM). Using several novel or state of the art methods, the piece software handles semi-concurrent multiscale methods as well as concurrent multiscale methods for fracture, coupling two continuum domains or atomistic domains to continuum domains, respectively. The efficiency of our open-source software is shown through several simulations including a 3D crack modeling in clay nanocomposites, a semi-concurrent FE-FE coupling, a 3D Arlequin multiscale example and an MD-XFEM coupling for dynamic crack propagation.
Multiscale computational modelling is challenging due to the high computational cost of direct numerical simulation by finite elements. To address this issue, concurrent multiscale methods use the ...solution of cheaper macroscale surrogates as boundary conditions to microscale sliding windows. The microscale problems remain a numerically challenging operation both in terms of implementation and cost. In this work we propose to replace the local microscale solution by an Encoder-Decoder Convolutional Neural Network that will generate fine-scale stress corrections to coarse predictions around unresolved microscale features, without prior parametrisation of local microscale problems. We deploy a Bayesian approach providing credible intervals to evaluate the uncertainty of the predictions, which is then used to investigate the merits of a selective learning framework. We will demonstrate the capability of the approach to predict equivalent stress fields in porous structures using linearised and finite strain elasticity theories.
We present a Nitche’s method to couple non-conforming two and three-dimensional non uniform rational b-splines (NURBS) patches in the context of isogeometric analysis. We present results for linear ...elastostatics in two and and three-dimensions. The method can deal with surface-surface or volume-volume coupling, and we show how it can be used to handle heterogeneities such as inclusions. We also present preliminary results on modal analysis. This simple coupling method has the potential to increase the applicability of NURBS-based isogeometric analysis for practical applications.
Spheroids encapsulated within alginate capsules are emerging as suitable in vitro tools to investigate the impact of mechanical forces on tumor growth since the internal tumor pressure can be ...retrieved from the deformation of the capsule. Here we focus on the particular case of Cellular Capsule Technology (CCT). We show in this contribution that a modeling approach accounting for the triphasic nature of the spheroid (extracellular matrix, tumor cells and interstitial fluid) offers a new perspective of analysis revealing that the pressure retrieved experimentally cannot be interpreted as a direct picture of the pressure sustained by the tumor cells and, as such, cannot therefore be used to quantify the critical pressure which induces stress-induced phenotype switch in tumor cells. The proposed multiphase reactive poro-mechanical model was cross-validated. Parameter sensitivity analyses on the digital twin revealed that the main parameters determining the encapsulated growth configuration are different from those driving growth in free condition, confirming that radically different phenomena are at play. Results reported in this contribution support the idea that multiphase reactive poro-mechanics is an exceptional theoretical framework to attain an in-depth understanding of CCT experiments, to confirm their hypotheses and to further improve their design.
The aim of this contribution is to explain in a straightforward manner how Bayesian inference can be used to identify material parameters of material models for solids. Bayesian approaches have ...already been used for this purpose, but most of the literature is not necessarily easy to understand for those new to the field. The reason for this is that most literature focuses either on complex statistical and machine learning concepts and/or on relatively complex mechanical models. In order to introduce the approach as gently as possible, we only focus on stress–strain measurements coming from uniaxial tensile tests and we only treat elastic and elastoplastic material models. Furthermore, the stress–strain measurements are created artificially in order to allow a one-to-one comparison between the true parameter values and the identified parameter distributions.
Anomalous proximity effects have been observed in adhesive systems ranging from proteins, bacteria, and gecko feet suspended over semiconductor surfaces to interfaces between graphene and different ...substrate materials. In the latter case, long-range forces are evidenced by measurements of non-vanishing stress that extends up to micrometer separations between graphene and the substrate. State-of-the-art models to describe adhesive properties are unable to explain these experimental observations, instead underestimating the measured stress distance range by 2-3 orders of magnitude. Here, we develop an analytical and numerical variational approach that combines continuum mechanics and elasticity with quantum many-body treatment of van der Waals dispersion interactions. A full relaxation of the coupled adsorbate/substrate geometry leads us to conclude that wavelike atomic deformation is largely responsible for the observed long-range proximity effect. The correct description of this seemingly general phenomenon for thin deformable membranes requires a direct coupling between quantum and continuum mechanics.
A numerical scheme is proposed for the detection of multiple cracks in three dimensional (3D) structures. The scheme is based on a variant of the extended finite element method (XFEM) and a hybrid ...optimizer solution. The proposed XFEM variant is particularly well-suited for the simulation of 3D fracture problems, and as such serves as an efficient solution to the so-called forward problem. A set of heuristic optimization algorithms are recombined into a multiscale optimization scheme. The introduced approach proves effective in tackling the complex inverse problem involved, where identification of multiple flaws is sought on the basis of sparse measurements collected near the structural boundary. The potential of the scheme is demonstrated through a set of numerical case studies of varying complexity.
Astrocytes with their specialised morphology are essential for brain homeostasis as metabolic mediators between blood vessels and neurons. In neurodegenerative diseases such as Alzheimer's disease ...(AD), astrocytes adopt reactive profiles with molecular and morphological changes that could lead to the impairment of their metabolic support and impact disease progression. However, the underlying mechanisms of how the metabolic function of human astrocytes is impaired by their morphological changes in AD are still elusive. To address this challenge, we developed and applied a metabolic multiscale modelling approach integrating the dynamics of metabolic energy pathways and physiological astrocyte morphologies acquired in human AD and age-matched control brain samples. The results demonstrate that the complex cell shape and intracellular organisation of energetic pathways determine the metabolic profile and support capacity of astrocytes in health and AD conditions. Thus, our mechanistic approach indicates the importance of spatial orchestration in metabolism and allows for the identification of protective mechanisms against disease-associated metabolic impairments.