-
FFT-based surrogate modeling of auxetic metamaterials with real-time prediction of effective elastic properties and swift inverse design
Authors:
Hooman Danesh,
Daniele Di Lorenzo,
Francisco Chinesta,
Stefanie Reese,
Tim Brepols
Abstract:
Auxetic structures, known for their negative Poisson's ratio, exhibit effective elastic properties heavily influenced by their underlying structural geometry and base material properties. While periodic homogenization of auxetic unit cells can be used to investigate these properties, it is computationally expensive and limits design space exploration and inverse analysis. In this paper, surrogate…
▽ More
Auxetic structures, known for their negative Poisson's ratio, exhibit effective elastic properties heavily influenced by their underlying structural geometry and base material properties. While periodic homogenization of auxetic unit cells can be used to investigate these properties, it is computationally expensive and limits design space exploration and inverse analysis. In this paper, surrogate models are developed for the real-time prediction of the effective elastic properties of auxetic unit cells with orthogonal voids of different shapes. The unit cells feature orthogonal voids in four distinct shapes, including rectangular, diamond, oval, and peanut-shaped voids, each characterized by specific void diameters. The generated surrogate models accept geometric parameters and the elastic properties of the base material as inputs to predict the effective elastic constants in real-time. This rapid evaluation enables a practical inverse analysis framework for obtaining the optimal design parameters that yield the desired effective response. The fast Fourier transform (FFT)-based homogenization approach is adopted to efficiently generate data for developing the surrogate models, bypassing concerns about periodic mesh generation and boundary conditions typically associated with the finite element method (FEM). The performance of the generated surrogate models is rigorously examined through a train/test split methodology, a parametric study, and an inverse problem. Finally, a graphical user interface (GUI) is developed, offering real-time prediction of the effective tangent stiffness and performing inverse analysis to determine optimal geometric parameters.
△ Less
Submitted 24 August, 2024;
originally announced August 2024.
-
An anisotropic, brittle damage model for finite strains with a generic damage tensor regularization
Authors:
Tim van der Velden,
Stefanie Reese,
Hagen Holthusen,
Tim Brepols
Abstract:
This paper establishes a universal framework for the nonlocal modeling of anisotropic damage at finite strains. By the combination of two recent works, the new framework allows for the flexible incorporation of different established hyperelastic finite strain material formulations into anisotropic damage whilst ensuring mesh-independent results by employing a generic set of micromorphic gradient-e…
▽ More
This paper establishes a universal framework for the nonlocal modeling of anisotropic damage at finite strains. By the combination of two recent works, the new framework allows for the flexible incorporation of different established hyperelastic finite strain material formulations into anisotropic damage whilst ensuring mesh-independent results by employing a generic set of micromorphic gradient-extensions. First, the anisotropic damage model, generally satisfying the damage growth criterion, is investigated for the specific choice of a Neo-Hookean material on a single element. Next, the model is applied with different gradient-extensions in structural simulations of an asymmetrically notched specimen to identify an efficient choice in the form of a volumetric-deviatoric regularization. Thereafter, the universal framework, which is without loss of generality here specified for a Neo-Hookean material with a volumetric-deviatoric gradient-extension, successfully serves for the complex simulation of a pressure loaded rotor blade.
After acceptance of the manuscript, we make the codes of the material subroutines accessible to the public at https://doi.org/10.5281/zenodo.11171630.
△ Less
Submitted 12 August, 2024;
originally announced August 2024.
-
Component based model order reduction with mortar tied contact for nonlinear quasi-static mechanical problems
Authors:
Stephan Ritzert,
Jannick Kehls,
Stefanie Reese,
Tim Brepols
Abstract:
In this work, we present a model order reduction technique for nonlinear structures assembled from components.The reduced order model is constructed by reducing the substructures with proper orthogonal decomposition and connecting them by a mortar-tied contact formulation. The snapshots for the substructure projection matrices are computed on the substructure level by the proper orthogonal decompo…
▽ More
In this work, we present a model order reduction technique for nonlinear structures assembled from components.The reduced order model is constructed by reducing the substructures with proper orthogonal decomposition and connecting them by a mortar-tied contact formulation. The snapshots for the substructure projection matrices are computed on the substructure level by the proper orthogonal decomposition (POD) method. The snapshots are computed using a random sampling procedure based on a parametrization of boundary conditions. To reduce the computational effort of the snapshot computation full-order simulations of the substructures are only computed when the error of the reduced solution is above a threshold. In numerical examples, we show the accuracy and efficiency of the method for nonlinear problems involving material and geometric nonlinearity as well as non-matching meshes. We are able to predict solutions of systems that we did not compute in our snapshots.
△ Less
Submitted 2 July, 2024;
originally announced July 2024.
-
A finite element-based physics-informed operator learning framework for spatiotemporal partial differential equations on arbitrary domains
Authors:
Yusuke Yamazaki,
Ali Harandi,
Mayu Muramatsu,
Alexandre Viardin,
Markus Apel,
Tim Brepols,
Stefanie Reese,
Shahed Rezaei
Abstract:
We propose a novel finite element-based physics-informed operator learning framework that allows for predicting spatiotemporal dynamics governed by partial differential equations (PDEs). The proposed framework employs a loss function inspired by the finite element method (FEM) with the implicit Euler time integration scheme. A transient thermal conduction problem is considered to benchmark the per…
▽ More
We propose a novel finite element-based physics-informed operator learning framework that allows for predicting spatiotemporal dynamics governed by partial differential equations (PDEs). The proposed framework employs a loss function inspired by the finite element method (FEM) with the implicit Euler time integration scheme. A transient thermal conduction problem is considered to benchmark the performance. The proposed operator learning framework takes a temperature field at the current time step as input and predicts a temperature field at the next time step. The Galerkin discretized weak formulation of the heat equation is employed to incorporate physics into the loss function, which is coined finite operator learning (FOL). Upon training, the networks successfully predict the temperature evolution over time for any initial temperature field at high accuracy compared to the FEM solution. The framework is also confirmed to be applicable to a heterogeneous thermal conductivity and arbitrary geometry. The advantages of FOL can be summarized as follows: First, the training is performed in an unsupervised manner, avoiding the need for a large data set prepared from costly simulations or experiments. Instead, random temperature patterns generated by the Gaussian random process and the Fourier series, combined with constant temperature fields, are used as training data to cover possible temperature cases. Second, shape functions and backward difference approximation are exploited for the domain discretization, resulting in a purely algebraic equation. This enhances training efficiency, as one avoids time-consuming automatic differentiation when optimizing weights and biases while accepting possible discretization errors. Finally, thanks to the interpolation power of FEM, any arbitrary geometry can be handled with FOL, which is crucial to addressing various engineering application scenarios.
△ Less
Submitted 6 August, 2024; v1 submitted 20 May, 2024;
originally announced May 2024.
-
A finite operator learning technique for mapping the elastic properties of microstructures to their mechanical deformations
Authors:
Shahed Rezaei,
Reza Najian Asl,
Shirko Faroughi,
Mahdi Asgharzadeh,
Ali Harandi,
Rasoul Najafi Koopas,
Gottfried Laschet,
Stefanie Reese,
Markus Apel
Abstract:
To obtain fast solutions for governing physical equations in solid mechanics, we introduce a method that integrates the core ideas of the finite element method with physics-informed neural networks and concept of neural operators. This approach generalizes and enhances each method, learning the parametric solution for mechanical problems without relying on data from other resources (e.g. other num…
▽ More
To obtain fast solutions for governing physical equations in solid mechanics, we introduce a method that integrates the core ideas of the finite element method with physics-informed neural networks and concept of neural operators. This approach generalizes and enhances each method, learning the parametric solution for mechanical problems without relying on data from other resources (e.g. other numerical solvers). We propose directly utilizing the available discretized weak form in finite element packages to construct the loss functions algebraically, thereby demonstrating the ability to find solutions even in the presence of sharp discontinuities. Our focus is on micromechanics as an example, where knowledge of deformation and stress fields for a given heterogeneous microstructure is crucial for further design applications. The primary parameter under investigation is the Young's modulus distribution within the heterogeneous solid system. Our investigations reveal that physics-based training yields higher accuracy compared to purely data-driven approaches for unseen microstructures. Additionally, we offer two methods to directly improve the process of obtaining high-resolution solutions, avoiding the need to use basic interpolation techniques. First is based on an autoencoder approach to enhance the efficiency for calculation on high resolution grid point. Next, Fourier-based parametrization is utilized to address complex 2D and 3D problems in micromechanics. The latter idea aims to represent complex microstructures efficiently using Fourier coefficients. Comparisons with other well-known operator learning algorithms, further emphasize the advantages of the newly proposed method.
△ Less
Submitted 3 June, 2024; v1 submitted 28 March, 2024;
originally announced April 2024.
-
Towards Understanding Counseling Conversations: Domain Knowledge and Large Language Models
Authors:
Younghun Lee,
Dan Goldwasser,
Laura Schwab Reese
Abstract:
Understanding the dynamics of counseling conversations is an important task, yet it is a challenging NLP problem regardless of the recent advance of Transformer-based pre-trained language models. This paper proposes a systematic approach to examine the efficacy of domain knowledge and large language models (LLMs) in better representing conversations between a crisis counselor and a help seeker. We…
▽ More
Understanding the dynamics of counseling conversations is an important task, yet it is a challenging NLP problem regardless of the recent advance of Transformer-based pre-trained language models. This paper proposes a systematic approach to examine the efficacy of domain knowledge and large language models (LLMs) in better representing conversations between a crisis counselor and a help seeker. We empirically show that state-of-the-art language models such as Transformer-based models and GPT models fail to predict the conversation outcome. To provide richer context to conversations, we incorporate human-annotated domain knowledge and LLM-generated features; simple integration of domain knowledge and LLM features improves the model performance by approximately 15%. We argue that both domain knowledge and LLM-generated features can be exploited to better characterize counseling conversations when they are used as an additional context to conversations.
△ Less
Submitted 21 February, 2024;
originally announced February 2024.
-
Mechanical Comparison of Arrangement Strategies for Topological Interlocking Assemblies
Authors:
Tom Goertzen,
Domen Macek,
Lukas Schnelle,
Meike Weiß,
Stefanie Reese,
Hagen Holthusen,
Alice C. Niemeyer
Abstract:
Topological Interlocking assemblies are arrangements of blocks kinematically constrained by a fixed frame, such that all rigid body motions of each block are constrained only by its permanent contact with other blocks and the frame. In the literature several blocks are introduced that can be arranged into different interlocking assemblies. In this study we investigate the influence of arrangement…
▽ More
Topological Interlocking assemblies are arrangements of blocks kinematically constrained by a fixed frame, such that all rigid body motions of each block are constrained only by its permanent contact with other blocks and the frame. In the literature several blocks are introduced that can be arranged into different interlocking assemblies. In this study we investigate the influence of arrangement on the overall structural behaviour of the resulting interlocking assemblies. This is performed using the Versatile Block, as it can be arranged in three different doubly periodic ways given by wallpaper symmetries. Our focus lies on the load transfer mechanisms from the assembly onto the frame. For fast a priori evaluation of the assemblies we introduce a combinatorial model called Interlocking Flows. To investigate our assemblies from a mechanical point of view we conduct several finite element studies. These reveal a strong influence of arrangement on the structural behaviour, for instance, an impact on both the point and amount of maximum deflection. The results of the finite element analysis are in very good agreement with the predictions of the Interlocking Flow model. Our source code, data and examples are available under https://doi.org/10.5281/zenodo.10246034.
△ Less
Submitted 4 December, 2023;
originally announced December 2023.
-
Discrete Empirical Interpolation Method for nonlinear softening problems involving damage and plasticity
Authors:
Steffen Kastian,
Jannick Kehls,
Tim Brepols,
Stefanie Reese
Abstract:
Accurate simulations are essential for engineering applications, and intricate continuum mechanical material models are constructed to achieve this goal. However, the increasing complexity of the material models and geometrical properties lead to a significant increase in computational effort. Model order reduction aims to implement efficient methods for accelerating the simulation process while p…
▽ More
Accurate simulations are essential for engineering applications, and intricate continuum mechanical material models are constructed to achieve this goal. However, the increasing complexity of the material models and geometrical properties lead to a significant increase in computational effort. Model order reduction aims to implement efficient methods for accelerating the simulation process while preserving a high degree of accuracy. Numerous studies have already demonstrated the benefits of this method for linear elastic material modeling. However, in the present work, we investigate a two-surface gradient-extended damage-plasticity model. We conducted complex simulations with this model, demonstrating both damage behavior and softening. The POD-based discrete empirical interpolation method (DEIM) is introduced and implemented. To accomplish simulations with DEIM and softening behaviour, we propose the implementation of a reduced form of the arc-length method. Existing research on calculating models with both damage and softening behavior using the DEIM and arc-length method is limited. To validate the methods, two numerical examples are thoroughly investigated in this study: a plate with a hole and an asymmetrically notched specimen. The results show that the proposed methods can create a reduced order model with high accuracy and a significant speedup of the simulation. For both examples, the analysis is conducted in three steps: first, plasticity without damage is examined, followed by damage without plasticity, and finally, the combination of plasticity and damage is investigated.
△ Less
Submitted 29 November, 2023;
originally announced November 2023.
-
A comparative study of micromorphic gradient-extensions for anisotropic damage at finite strains
Authors:
Tim van der Velden,
Tim Brepols,
Stefanie Reese,
Hagen Holthusen
Abstract:
Modern inelastic material model formulations rely on the use of tensor-valued internal variables. When inelastic phenomena include softening, simulations of the former are prone to localization. Thus, an accurate regularization of the tensor-valued internal variables is essential to obtain physically correct results. Here, we focus on the regularization of anisotropic damage at finite strains. Thu…
▽ More
Modern inelastic material model formulations rely on the use of tensor-valued internal variables. When inelastic phenomena include softening, simulations of the former are prone to localization. Thus, an accurate regularization of the tensor-valued internal variables is essential to obtain physically correct results. Here, we focus on the regularization of anisotropic damage at finite strains. Thus, a flexible anisotropic damage model with isotropic, kinematic, and distortional hardening is equipped with three gradient-extensions using a full and two reduced regularizations of the damage tensor. Theoretical and numerical comparisons of the three gradient-extensions yield excellent agreement between the full and the reduced regularization based on a volumetric-deviatoric regularization using only two nonlocal degrees of freedom.
△ Less
Submitted 27 November, 2023;
originally announced November 2023.
-
Theory and implementation of inelastic Constitutive Artificial Neural Networks
Authors:
Hagen Holthusen,
Lukas Lamm,
Tim Brepols,
Stefanie Reese,
Ellen Kuhl
Abstract:
Nature has always been our inspiration in the research, design and development of materials and has driven us to gain a deep understanding of the mechanisms that characterize anisotropy and inelastic behavior. All this knowledge has been accumulated in the principles of thermodynamics. Deduced from these principles, the multiplicative decomposition combined with pseudo potentials are powerful and…
▽ More
Nature has always been our inspiration in the research, design and development of materials and has driven us to gain a deep understanding of the mechanisms that characterize anisotropy and inelastic behavior. All this knowledge has been accumulated in the principles of thermodynamics. Deduced from these principles, the multiplicative decomposition combined with pseudo potentials are powerful and universal concepts. Simultaneously, the tremendous increase in computational performance enabled us to investigate and rethink our history-dependent material models to make the most of our predictions. Today, we have reached a point where materials and their models are becoming increasingly sophisticated. This raises the question: How do we find the best model that includes all inelastic effects to explain our complex data? Constitutive Artificial Neural Networks (CANN) may answer this question. Here, we extend the CANNs to inelastic materials (iCANN). Rigorous considerations of objectivity, rigid motion of the reference configuration, multiplicative decomposition and its inherent non-uniqueness, restrictions of energy and pseudo potential, and consistent evolution guide us towards the architecture of the iCANN satisfying thermodynamics per design. We combine feed-forward networks of the free energy and pseudo potential with a recurrent neural network approach to take time dependencies into account. We demonstrate that the iCANN is capable of autonomously discovering models for artificially generated data, the response of polymers for cyclic loading and the relaxation behavior of muscle data. As the design of the network is not limited to visco-elasticity, our vision is that the iCANN will reveal to us new ways to find the various inelastic phenomena hidden in the data and to understand their interaction. Our source code, data, and examples are available at doi.org/10.5281/zenodo.10066805
△ Less
Submitted 10 November, 2023;
originally announced November 2023.
-
Computational modeling of in-stent restenosis: Pharmacokinetic and pharmacodynamic evaluation
Authors:
Kiran Manjunatha,
Nicole Schaaps,
Marek Behr,
Felix Vogt,
Stefanie Reese
Abstract:
Persistence of the pathology of in-stent restenosis even with the advent of drug-eluting stents warrants the development of highly resolved in silico models. These computational models assist in gaining insights into the transient biochemical and cellular mechanisms involved and thereby optimize the stent implantation parameters. Within this work, an already established fully-coupled Lagrangian fi…
▽ More
Persistence of the pathology of in-stent restenosis even with the advent of drug-eluting stents warrants the development of highly resolved in silico models. These computational models assist in gaining insights into the transient biochemical and cellular mechanisms involved and thereby optimize the stent implantation parameters. Within this work, an already established fully-coupled Lagrangian finite element framework for modeling the restenotic growth is enhanced with the incorporation of endothelium-mediated effects and pharmacological influences of rapamycin-based drugs embedded in the polymeric layers of the current generation drug-eluting stents. The continuum mechanical description of growth is further justified in the context of thermodynamic consistency. Qualitative inferences are drawn from the model developed herein regarding the efficacy of the level of drug embedment within the struts as well as the release profiles adopted. The framework is then intended to serve as a tool for clinicians to tune the interventional procedures patient-specifically.
△ Less
Submitted 26 July, 2023; v1 submitted 24 July, 2023;
originally announced July 2023.
-
Mechanical modeling of the maturation process for tissue-engineered implants: application to biohybrid heart valves
Authors:
Mahmoud Sesa,
Hagen Holthusen,
Lukas Lamm,
Christian Böhm,
Tim Brepols,
Stefan Jockenhövel,
Stefanie Reese
Abstract:
The development of tissue-engineered cardiovascular implants can improve the lives of large segments of our society who suffer from cardiovascular diseases. Regenerative tissues are fabricated using a process called tissue maturation. Furthermore, it is highly challenging to produce cardiovascular regenerative implants with sufficient mechanical strength to withstand the loading conditions within…
▽ More
The development of tissue-engineered cardiovascular implants can improve the lives of large segments of our society who suffer from cardiovascular diseases. Regenerative tissues are fabricated using a process called tissue maturation. Furthermore, it is highly challenging to produce cardiovascular regenerative implants with sufficient mechanical strength to withstand the loading conditions within the human body. Therefore, biohybrid implants for which the regenerative tissue is reinforced by standard reinforcement material (e.g. textile or 3d printed scaffold) can be an interesting solution. In silico models can significantly contribute to characterizing, designing, and optimizing biohybrid implants. The first step towards this goal is to develop a computational model for the maturation process of tissue-engineered implants. This paper focuses on the mechanical modeling of textile-reinforced tissue-engineered cardiovascular implants. First, we propose an energy-based approach to compute the collagen evolution during the maturation process. Then, we apply the concept of structural tensors to model the anisotropic behavior of the extracellular matrix and the textile scaffold. Next, the newly developed material model is embedded into a special solid-shell finite element formulation with reduced integration. Finally, we use our framework to compute two structural problems: a pressurized shell construct and a tubular-shaped heart valve. The results show the ability of the model to predict collagen growth in response to the boundary conditions applied during the maturation process. Consequently, we can predict the implant's mechanical response, such as the deformation and stresses of the implant.
△ Less
Submitted 23 July, 2023;
originally announced July 2023.
-
Mixed formulation of physics-informed neural networks for thermo-mechanically coupled systems and heterogeneous domains
Authors:
Ali Harandi,
Ahmad Moeineddin,
Michael Kaliske,
Stefanie Reese,
Shahed Rezaei
Abstract:
Physics-informed neural networks (PINNs) are a new tool for solving boundary value problems by defining loss functions of neural networks based on governing equations, boundary conditions, and initial conditions. Recent investigations have shown that when designing loss functions for many engineering problems, using first-order derivatives and combining equations from both strong and weak forms ca…
▽ More
Physics-informed neural networks (PINNs) are a new tool for solving boundary value problems by defining loss functions of neural networks based on governing equations, boundary conditions, and initial conditions. Recent investigations have shown that when designing loss functions for many engineering problems, using first-order derivatives and combining equations from both strong and weak forms can lead to much better accuracy, especially when there are heterogeneity and variable jumps in the domain. This new approach is called the mixed formulation for PINNs, which takes ideas from the mixed finite element method. In this method, the PDE is reformulated as a system of equations where the primary unknowns are the fluxes or gradients of the solution, and the secondary unknowns are the solution itself. In this work, we propose applying the mixed formulation to solve multi-physical problems, specifically a stationary thermo-mechanically coupled system of equations. Additionally, we discuss both sequential and fully coupled unsupervised training and compare their accuracy and computational cost. To improve the accuracy of the network, we incorporate hard boundary constraints to ensure valid predictions. We then investigate how different optimizers and architectures affect accuracy and efficiency. Finally, we introduce a simple approach for parametric learning that is similar to transfer learning. This approach combines data and physics to address the limitations of PINNs regarding computational cost and improves the network's ability to predict the response of the system for unseen cases. The outcomes of this work will be useful for many other engineering applications where deep learning is employed on multiple coupled systems of equations for fast and reliable computations.
△ Less
Submitted 6 September, 2023; v1 submitted 9 February, 2023;
originally announced February 2023.
-
Model-Free Data-Driven Inference in Computational Mechanics
Authors:
Erik Prume,
Stefanie Reese,
Michael Ortiz
Abstract:
We extend the model-free Data-Driven computing paradigm to solids and structures that are stochastic due to intrinsic randomness in the material behavior. The behavior of such materials is characterized by a likelihood measure instead of a constitutive relation. We specifically assume that the material likelihood measure is known only through an empirical point-data set in material or phase space.…
▽ More
We extend the model-free Data-Driven computing paradigm to solids and structures that are stochastic due to intrinsic randomness in the material behavior. The behavior of such materials is characterized by a likelihood measure instead of a constitutive relation. We specifically assume that the material likelihood measure is known only through an empirical point-data set in material or phase space. The state of the solid or structure is additionally subject to compatibility and equilibrium constraints. The problem is then to infer the likelihood of a given structural outcome of interest. In this work, we present a Data-Driven method of inference that determines likelihoods of outcomes from the empirical material data and that requires no material or prior modeling. In particular, the computation of expectations is reduced to explicit sums over local material data sets and to quadratures over admissible states, i. e., states satisfying compatibility and equilibrium. The complexity of the material data-set sums is linear in the number of data points and in the number of members in the structure. Efficient population annealing procedures and fast search algorithms for accelerating the calculations are presented. The scope, cost and convergence properties of the method are assessed with the aid selected applications and benchmark tests.
△ Less
Submitted 28 October, 2022; v1 submitted 13 July, 2022;
originally announced July 2022.
-
A mixed formulation for physics-informed neural networks as a potential solver for engineering problems in heterogeneous domains: comparison with finite element method
Authors:
Shahed Rezaei,
Ali Harandi,
Ahmad Moeineddin,
Bai-Xiang Xu,
Stefanie Reese
Abstract:
Physics-informed neural networks (PINNs) are capable of finding the solution for a given boundary value problem. We employ several ideas from the finite element method (FEM) to enhance the performance of existing PINNs in engineering problems. The main contribution of the current work is to promote using the spatial gradient of the primary variable as an output from separated neural networks. Late…
▽ More
Physics-informed neural networks (PINNs) are capable of finding the solution for a given boundary value problem. We employ several ideas from the finite element method (FEM) to enhance the performance of existing PINNs in engineering problems. The main contribution of the current work is to promote using the spatial gradient of the primary variable as an output from separated neural networks. Later on, the strong form which has a higher order of derivatives is applied to the spatial gradients of the primary variable as the physical constraint. In addition, the so-called energy form of the problem is applied to the primary variable as an additional constraint for training. The proposed approach only required up to first-order derivatives to construct the physical loss functions. We discuss why this point is beneficial through various comparisons between different models. The mixed formulation-based PINNs and FE methods share some similarities. While the former minimizes the PDE and its energy form at given collocation points utilizing a complex nonlinear interpolation through a neural network, the latter does the same at element nodes with the help of shape functions. We focus on heterogeneous solids to show the capability of deep learning for predicting the solution in a complex environment under different boundary conditions. The performance of the proposed PINN model is checked against the solution from FEM on two prototype problems: elasticity and the Poisson equation (steady-state diffusion problem). We concluded that by properly designing the network architecture in PINN, the deep learning model has the potential to solve the unknowns in a heterogeneous domain without any available initial data from other sources. Finally, discussions are provided on the combination of PINN and FEM for a fast and accurate design of composite materials in future developments.
△ Less
Submitted 27 June, 2022;
originally announced June 2022.
-
A multiphysics modeling approach for in-stent restenosis: Theoretical aspects and finite element implementation
Authors:
Kiran Manjunatha,
Marek Behr,
Felix Vogt,
Stefanie Reese
Abstract:
Development of in silico models are intrinsic in understanding disease progression in soft biological tissues. Within this work, we propose a fully-coupled Lagrangian finite element framework which replicates the process of in-stent restenosis observed post stent implantation in a coronary artery. Coupled advection-reaction-diffusion reactions are set up that track the evolution of the concentrati…
▽ More
Development of in silico models are intrinsic in understanding disease progression in soft biological tissues. Within this work, we propose a fully-coupled Lagrangian finite element framework which replicates the process of in-stent restenosis observed post stent implantation in a coronary artery. Coupled advection-reaction-diffusion reactions are set up that track the evolution of the concentrations of the platelet-derived growth factor, the transforming growth factor-$β$, the extracellular matrix, and the density of the smooth muscle cells. A continuum mechanical description of growth incorporating the evolution of arterial wall constituents is developed, and a suitable finite element implementation discussed. Qualitative validation of the computational model are presented by emulating a stented artery. Patient-specific data can be integrated into the model to predict the risk of restenosis and thereby assist in tuning of stent implantation parameters to mitigate the risk.
△ Less
Submitted 7 April, 2022; v1 submitted 5 April, 2022;
originally announced April 2022.
-
Modeling moving boundary value problems in electrochemical machining
Authors:
Tim van der Velden,
Stephan Ritzert,
Stefanie Reese,
Johanna Waimann
Abstract:
This work presents a new approach to efficiently model the cathode in the moving boundary value problem of electrochemical machining. Until recently, the process simulation with finite elements had the drawback of remeshing required by the changing surface geometries. This disadvantage was overcome by a novel model formulation for the anodic dissolution that utilizes effective material parameters…
▽ More
This work presents a new approach to efficiently model the cathode in the moving boundary value problem of electrochemical machining. Until recently, the process simulation with finite elements had the drawback of remeshing required by the changing surface geometries. This disadvantage was overcome by a novel model formulation for the anodic dissolution that utilizes effective material parameters as well as the dissolution level as an internal variable and, thereby, does not require remeshing. Now, we extend this concept to model arbitrarily shaped and moving cathodes. Two methodologies are investigated to describe the time varying position of the cathode. In the first approach, we change the electric conductivity of elements within the cathode and, in a second approach, we apply Dirichlet boundary conditions on the nodes of corresponding elements. For both methods, elements on the cathode's surface are treated with effective material parameters. This procedure allows for the efficient simulation of industrially relevant, complex geometries without mesh adaptation. The model's performance is validated by means of analytical, numerical and experimental results from the literature. The short computation times make the approach interesting for industrial applications.
△ Less
Submitted 25 February, 2022;
originally announced February 2022.
-
An anisotropic cohesive fracture model: advantages and limitations of length-scale insensitive phase-field damage models
Authors:
Shahed Rezaei,
Ali Harandi,
Tim Brepols,
Stefanie Reese
Abstract:
The goal of the current work is to explore direction-dependent damage initiation and propagation within an arbitrary anisotropic solid. In particular, we aim at developing anisotropic cohesive phase-field (PF) damage models by extending the idea introduced in \cite{REZAEI2021a} for direction-dependent fracture energy and also anisotropic PF damage models based on structural tensors. The cohesive P…
▽ More
The goal of the current work is to explore direction-dependent damage initiation and propagation within an arbitrary anisotropic solid. In particular, we aim at developing anisotropic cohesive phase-field (PF) damage models by extending the idea introduced in \cite{REZAEI2021a} for direction-dependent fracture energy and also anisotropic PF damage models based on structural tensors. The cohesive PF damage formulation used in the current contribution is motivated by the works of \cite{LORENTZ201120, wu2018, GEELEN2019}. The results of the latter models are shown to be insensitive with respect to the length scale parameter for the isotropic case. This is because they manage to formulate the fracture energy as a function of diffuse displacement jumps in the localized damaged zone. In the present paper, we discuss numerical examples and details on finite element implementations where the fracture energy, as well as the material strength, are introduced as an arbitrary function of the crack direction. Using the current formulation for anisotropic cohesive fracture, the obtained results are almost insensitive with respect to the length scale parameter. The latter is achieved by including the direction-dependent strength of the material in addition to its fracture energy. Utilizing the current formulation, one can increase the mesh size which reduces the computational time significantly without any severe change in the predicted crack path and overall obtained load-displacement curves. We also argue that these models still lack to capture mode-dependent fracture properties. Open issues and possible remedies for future developments are finally discussed as well.
△ Less
Submitted 28 August, 2021;
originally announced August 2021.
-
A macroscopic approach for stress driven anisotropic growth in bioengineered soft tissues
Authors:
L. Lamm,
H. Holthusen,
T. Brepols,
S. Jockenhövel,
S. Reese
Abstract:
The simulation of growth processes within soft biological tissues is of utmost importance for many applications in the medical sector. Within this contribution we propose a new macroscopic approach fro modelling stress-driven volumetric growth occurring in soft tissues. Instead of using the standard approach of a-priori defining the structure of the growth tensor, we postulate the existance of a g…
▽ More
The simulation of growth processes within soft biological tissues is of utmost importance for many applications in the medical sector. Within this contribution we propose a new macroscopic approach fro modelling stress-driven volumetric growth occurring in soft tissues. Instead of using the standard approach of a-priori defining the structure of the growth tensor, we postulate the existance of a general growth potential. Such a potential describes all eligable homeostatic stress states that can ultimately be reached as a result of the growth process. Making use of well established methods from visco-plasticity, the evolution of the growth related right Cauchy-Green tensor is subsequently defined as a time dependent associative evolution law with respect to the introduced potential. This approach naturally leads to a formulation that is able to cover both, isotropic and anisotropic growth related changes in geometry. It furthermore allows the model to flexibly adapt to changing boundary and loading conditions. Besides the theoretical development, we also describe the algorithmic implementation and furthermore compare the newly derived model with a standard formulation of isotropic growth.
△ Less
Submitted 8 July, 2021;
originally announced July 2021.
-
A geometrically adapted reduced set of frequencies for a FFT-based microstructure simulation
Authors:
Christian Gierden,
Johanna Waimann,
Bob Svendsen,
Stefanie Reese
Abstract:
We present a modified model order reduction (MOR) technique for the FFT-based simulation of composite microstructures. It utilizes the earlier introduced MOR technique (Kochmann et al. [2019]), which is based on solving the Lippmann-Schwinger equation in Fourier space by a reduced set of frequencies. Crucial for the accuracy of this MOR technique is on the one hand the amount of used frequencies a…
▽ More
We present a modified model order reduction (MOR) technique for the FFT-based simulation of composite microstructures. It utilizes the earlier introduced MOR technique (Kochmann et al. [2019]), which is based on solving the Lippmann-Schwinger equation in Fourier space by a reduced set of frequencies. Crucial for the accuracy of this MOR technique is on the one hand the amount of used frequencies and on the other hand the choice of frequencies used within the simulation. Kochmann et al. [2019] defined the reduced set of frequencies by using a fixed sampling pattern, which is most general but leads to poor microstructural results when considering only a few frequencies. Consequently, a reconstruction algorithm based on the TV1-algorithm [Candes et al., 2006] was used in a post-processing step to generate highly resolved micromechanical fields. The present work deals with a modified sampling pattern generation for this MOR technique. Based on the idea, that the micromechanical material response strongly depends on the phase-wise material behavior, we propose the usage of sampling patterns adapted to the spatial arrangement of the individual phases. This leads to significantly improved microscopic and overall results. Hence, the time-consuming reconstruction in the post-processing step that was necessary in the earlier work is no longer required. To show the adaptability and robustness of this new choice of sampling patterns, several two dimensional examples are investigated. In addition, also the 3D extension of the algorithm is presented.
△ Less
Submitted 18 March, 2021;
originally announced March 2021.
-
A novel approach for the efficient modeling of material dissolution in electrochemical machining
Authors:
Tim van der Velden,
Bob Rommes,
Andreas Klink,
Stefanie Reese,
Johanna Waimann
Abstract:
This work presents a novel approach to efficiently model anodic dissolution in electrochemical machining. Earlier modeling approaches employ a strict space discretization of the anodic surface that is associated with a remeshing procedure at every time step. Besides that, the presented model is formulated by means of effective material parameters. Thereby, it allows to use a constant mesh for the…
▽ More
This work presents a novel approach to efficiently model anodic dissolution in electrochemical machining. Earlier modeling approaches employ a strict space discretization of the anodic surface that is associated with a remeshing procedure at every time step. Besides that, the presented model is formulated by means of effective material parameters. Thereby, it allows to use a constant mesh for the entire simulation and, thus, decreases the computational costs. Based on Faraday's law of electrolysis, an effective dissolution level is introduced, which describes the ratio of a dissolved volume and its corresponding reference volume. This inner variable allows the modeling of the complex dissolution process without the necessity of computationally expensive remeshing by controlling the effective material parameters. Additionally, full coupling of the thermoelectric problem is considered and its linearization and numerical implementation are presented. The model shows good agreement with analytical and experimental validation examples by yielding realistic results. Furthermore, simulations of a pulsed electrochemical machining process yield a process signature of the surface roughness related to the specific accumulated electric charge. The numerical examples confirm the simulation's computational efficiency and accurate modeling qualities.
△ Less
Submitted 15 March, 2021;
originally announced March 2021.
-
Finite element modelling of in-stent restenosis
Authors:
Kiran Manjunatha,
Marek Behr,
Felix Vogt,
Stefanie Reese
Abstract:
From the perspective of coronary heart disease, the development of stents has come significantly far in reducing the associated mortality rate, drug-eluting stents being the epitome of innovative and effective solutions. Within this work, the intricate process of in-stent restenosis is modelled considering one of the significant growth factors and its effect on constituents of the arterial wall. A…
▽ More
From the perspective of coronary heart disease, the development of stents has come significantly far in reducing the associated mortality rate, drug-eluting stents being the epitome of innovative and effective solutions. Within this work, the intricate process of in-stent restenosis is modelled considering one of the significant growth factors and its effect on constituents of the arterial wall. A multiphysical modelling approach is adopted in this regard. Experimental investigations from the literature have been used to hypothesize the governing equations and the corresponding parameters. A staggered solution strategy is utilised to capture the transport phenomena as well as the growth and remodeling that follows stent implantation. The model herein developed serves as a tool to predict in-stent restenosis depending on the endothelial injury sustained and the protuberance of stents into the lumen of the arteries. Keywords: in-stent restenosis, smooth mucsle cells, platelet-derived growth factor, extracellular matrix, growth
△ Less
Submitted 14 December, 2020; v1 submitted 5 December, 2020;
originally announced December 2020.
-
Efficient Data Structures for Model-free Data-Driven Computational Mechanics
Authors:
Robert Eggersmann,
Laurent Stainier,
Michael Ortiz,
Stefanie Reese
Abstract:
The data-driven computing paradigm initially introduced by Kirchdoerfer and Ortiz (2016) enables finite element computations in solid mechanics to be performed directly from material data sets, without an explicit material model. From a computational effort point of view, the most challenging task is the projection of admissible states at material points onto their closest states in the material d…
▽ More
The data-driven computing paradigm initially introduced by Kirchdoerfer and Ortiz (2016) enables finite element computations in solid mechanics to be performed directly from material data sets, without an explicit material model. From a computational effort point of view, the most challenging task is the projection of admissible states at material points onto their closest states in the material data set. In this study, we compare and develop several possible data structures for solving the nearest-neighbor problem. We show that approximate nearest-neighbor (ANN) algorithms can accelerate material data searches by several orders of magnitude relative to exact searching algorithms. The approximations are suggested by--and adapted to--the structure of the data-driven iterative solver and result in no significant loss of solution accuracy. We assess the performance of the ANN algorithm with respect to material data set size with the aid of a 3D elasticity test case. We show that computations on a single processor with up to one billion material data points are feasible within a few seconds execution time with a speedup of more than 106 with respect to exact k-d trees.
△ Less
Submitted 1 December, 2020;
originally announced December 2020.
-
Model-free Data-Driven Computational Mechanics Enhanced by Tensor Voting
Authors:
Robert Eggersmann,
Laurent Stainier,
Michael Ortiz,
Stefanie Reese
Abstract:
The data-driven computing paradigm initially introduced by Kirchdoerfer & Ortiz (2016) is extended by incorporating locally linear tangent spaces into the data set. These tangent spaces are constructed by means of the tensor voting method introduced by Mordohai & Medioni (2010) which improves the learning of the underlying structure of a data set. Tensor voting is an instance-based machine learnin…
▽ More
The data-driven computing paradigm initially introduced by Kirchdoerfer & Ortiz (2016) is extended by incorporating locally linear tangent spaces into the data set. These tangent spaces are constructed by means of the tensor voting method introduced by Mordohai & Medioni (2010) which improves the learning of the underlying structure of a data set. Tensor voting is an instance-based machine learning technique which accumulates votes from the nearest neighbors to build up second-order tensors encoding tangents and normals to the underlying data structure. The here proposed second-order data-driven paradigm is a plug-in method for distance-minimizing as well as entropy-maximizing data-driven schemes. Like its predecessor, the resulting method aims to minimize a suitably defined free energy over phase space subject to compatibility and equilibrium constraints. The method's implementation is straightforward and numerically efficient since the data structure analysis is performed in an offline step. Selected numerical examples are presented that establish the higher-order convergence properties of the data-driven solvers enhanced by tensor voting for ideal and noisy data sets.
△ Less
Submitted 7 April, 2020; v1 submitted 6 April, 2020;
originally announced April 2020.
-
Model-Free Data-Driven Inelasticity
Authors:
Robert Eggersmann,
Trenton Kirchdoerfer,
Stefanie Reese,
Laurent Stainier,
Michael Ortiz
Abstract:
We extend the Data-Driven formulation of problems in elasticity of Kirchdoerfer and Ortiz (2016) to inelasticity. This extension differs fundamentally from Data-Driven problems in elasticity in that the material data set evolves in time as a consequence of the history dependence of the material. We investigate three representational paradigms for the evolving material data sets: i) materials with…
▽ More
We extend the Data-Driven formulation of problems in elasticity of Kirchdoerfer and Ortiz (2016) to inelasticity. This extension differs fundamentally from Data-Driven problems in elasticity in that the material data set evolves in time as a consequence of the history dependence of the material. We investigate three representational paradigms for the evolving material data sets: i) materials with memory, i.e., conditioning the material data set to the past history of deformation; ii) differential materials, i.e., conditioning the material data set to short histories of stress and strain; and iii) history variables, i.e., conditioning the material data set to ad hoc variables encoding partial information about the history of stress and strain. We also consider combinations of the three paradigms thereof and investigate their ability to represent the evolving data sets of different classes of inelastic materials, including viscoelasticity, viscoplasticity and plasticity. We present selected numerical examples that demonstrate the range and scope of Data-Driven inelasticity and the numerical performance of implementations thereof.
△ Less
Submitted 10 September, 2018; v1 submitted 31 August, 2018;
originally announced August 2018.