-
Intrinsic mixed-dimensional beam-shell-solid couplings in linear Cosserat continua via tangential differential calculus
Authors:
Adam Sky,
Jack S. Hale,
Andreas Zilian,
Stéphane P. A. Bordas,
Patrizio Neff
Abstract:
We present an approach to the coupling of mixed-dimensional continua by employing the mathematically enriched linear Cosserat micropolar model. The kinematical reduction of the model to lower dimensional domains leaves its fundamental degrees of freedom intact. Consequently, the degrees of freedom intrinsically agree even at the interface with a domain of a different dimensionality. Thus, this app…
▽ More
We present an approach to the coupling of mixed-dimensional continua by employing the mathematically enriched linear Cosserat micropolar model. The kinematical reduction of the model to lower dimensional domains leaves its fundamental degrees of freedom intact. Consequently, the degrees of freedom intrinsically agree even at the interface with a domain of a different dimensionality. Thus, this approach circumvents the need for intermediate finite elements or mortar methods. We introduce the derivations of all models of various dimensions using tangential differential calculus. The coupling itself is then achieved by defining a mixed-dimensional action functional with consistent Sobolev trace operators. Finally, we present numerical examples involving a three-dimensional silicone-rubber block reinforced with a curved graphite shell on its lower surface, a three-dimensional silver block reinforced with a graphite plate and beams, and lastly, intersecting silver shells reinforced with graphite beams.
△ Less
Submitted 1 September, 2024; v1 submitted 17 July, 2024;
originally announced July 2024.
-
Gaussian process regression + deep neural network autoencoder for probabilistic surrogate modeling in nonlinear mechanics of solids
Authors:
Saurabh Deshpande,
Hussein Rappel,
Mark Hobbs,
Stéphane P. A. Bordas,
Jakub Lengiewicz
Abstract:
Many real-world applications demand accurate and fast predictions, as well as reliable uncertainty estimates. However, quantifying uncertainty on high-dimensional predictions is still a severely under-invested problem, especially when input-output relationships are non-linear. To handle this problem, the present work introduces an innovative approach that combines autoencoder deep neural networks…
▽ More
Many real-world applications demand accurate and fast predictions, as well as reliable uncertainty estimates. However, quantifying uncertainty on high-dimensional predictions is still a severely under-invested problem, especially when input-output relationships are non-linear. To handle this problem, the present work introduces an innovative approach that combines autoencoder deep neural networks with the probabilistic regression capabilities of Gaussian processes. The autoencoder provides a low-dimensional representation of the solution space, while the Gaussian process is a Bayesian method that provides a probabilistic mapping between the low-dimensional inputs and outputs. We validate the proposed framework for its application to surrogate modeling of non-linear finite element simulations. Our findings highlight that the proposed framework is computationally efficient as well as accurate in predicting non-linear deformations of solid bodies subjected to external forces, all the while providing insightful uncertainty assessments.
△ Less
Submitted 15 July, 2024;
originally announced July 2024.
-
Graph Representation Learning Strategies for Omics Data: A Case Study on Parkinson's Disease
Authors:
Elisa Gómez de Lope,
Saurabh Deshpande,
Ramón Viñas Torné,
Pietro Liò,
Enrico Glaab,
Stéphane P. A. Bordas
Abstract:
Omics data analysis is crucial for studying complex diseases, but its high dimensionality and heterogeneity challenge classical statistical and machine learning methods. Graph neural networks have emerged as promising alternatives, yet the optimal strategies for their design and optimization in real-world biomedical challenges remain unclear. This study evaluates various graph representation learn…
▽ More
Omics data analysis is crucial for studying complex diseases, but its high dimensionality and heterogeneity challenge classical statistical and machine learning methods. Graph neural networks have emerged as promising alternatives, yet the optimal strategies for their design and optimization in real-world biomedical challenges remain unclear. This study evaluates various graph representation learning models for case-control classification using high-throughput biological data from Parkinson's disease and control samples. We compare topologies derived from sample similarity networks and molecular interaction networks, including protein-protein and metabolite-metabolite interactions (PPI, MMI). Graph Convolutional Network (GCNs), Chebyshev spectral graph convolution (ChebyNet), and Graph Attention Network (GAT), are evaluated alongside advanced architectures like graph transformers, the graph U-net, and simpler models like multilayer perceptron (MLP).
These models are systematically applied to transcriptomics and metabolomics data independently. Our comparative analysis highlights the benefits and limitations of various architectures in extracting patterns from omics data, paving the way for more accurate and interpretable models in biomedical research.
△ Less
Submitted 20 June, 2024;
originally announced June 2024.
-
Discovering deposition process regimes: leveraging unsupervised learning for process insights, surrogate modeling, and sensitivity analysis
Authors:
Geremy Loachamín Suntaxi,
Paris Papavasileiou,
Eleni D. Koronaki,
Dimitrios G. Giovanis,
Georgios Gakis,
Ioannis G. Aviziotis,
Martin Kathrein,
Gabriele Pozzetti,
Christoph Czettl,
Stéphane P. A. Bordas,
Andreas G. Boudouvis
Abstract:
This work introduces a comprehensive approach utilizing data-driven methods to elucidate the deposition process regimes in Chemical Vapor Deposition (CVD) reactors and the interplay of physical mechanism that dominate in each one of them. Through this work, we address three key objectives. Firstly, our methodology relies on process outcomes, derived by a detailed CFD model, to identify clusters of…
▽ More
This work introduces a comprehensive approach utilizing data-driven methods to elucidate the deposition process regimes in Chemical Vapor Deposition (CVD) reactors and the interplay of physical mechanism that dominate in each one of them. Through this work, we address three key objectives. Firstly, our methodology relies on process outcomes, derived by a detailed CFD model, to identify clusters of "outcomes" corresponding to distinct process regimes, wherein the relative influence of input variables undergoes notable shifts. This phenomenon is experimentally validated through Arrhenius plot analysis, affirming the efficacy of our approach. Secondly, we demonstrate the development of an efficient surrogate model, based on Polynomial Chaos Expansion (PCE), that maintains accuracy, facilitating streamlined computational analyses. Finally, as a result of PCE, sensitivity analysis is made possible by means of Sobol' indices, that quantify the impact of process inputs across identified regimes. The insights gained from our analysis contribute to the formulation of hypotheses regarding phenomena occurring beyond the transition regime. Notably, the significance of temperature even in the diffusion-limited regime, as evidenced by the Arrhenius plot, suggests activation of gas phase reactions at elevated temperatures. Importantly, our proposed methods yield insights that align with experimental observations and theoretical principles, aiding decision-making in process design and optimization. By circumventing the need for costly and time-consuming experiments, our approach offers a pragmatic pathway towards enhanced process efficiency. Moreover, this study underscores the potential of data-driven computational methods for innovating reactor design paradigms.
△ Less
Submitted 24 May, 2024;
originally announced May 2024.
-
Integrating supervised and unsupervised learning approaches to unveil critical process inputs
Authors:
Paris Papavasileiou,
Dimitrios G. Giovanis,
Gabriele Pozzetti,
Martin Kathrein,
Christoph Czettl,
Ioannis G. Kevrekidis,
Andreas G. Boudouvis,
Stéphane P. A. Bordas,
Eleni D. Koronaki
Abstract:
This study introduces a machine learning framework tailored to large-scale industrial processes characterized by a plethora of numerical and categorical inputs. The framework aims to (i) discern critical parameters influencing the output and (ii) generate accurate out-of-sample qualitative and quantitative predictions of production outcomes. Specifically, we address the pivotal question of the sig…
▽ More
This study introduces a machine learning framework tailored to large-scale industrial processes characterized by a plethora of numerical and categorical inputs. The framework aims to (i) discern critical parameters influencing the output and (ii) generate accurate out-of-sample qualitative and quantitative predictions of production outcomes. Specifically, we address the pivotal question of the significance of each input in shaping the process outcome, using an industrial Chemical Vapor Deposition (CVD) process as an example. The initial objective involves merging subject matter expertise and clustering techniques exclusively on the process output, here, coating thickness measurements at various positions in the reactor. This approach identifies groups of production runs that share similar qualitative characteristics, such as film mean thickness and standard deviation. In particular, the differences of the outcomes represented by the different clusters can be attributed to differences in specific inputs, indicating that these inputs are critical for the production outcome. Leveraging this insight, we subsequently implement supervised classification and regression methods using the identified critical process inputs. The proposed methodology proves to be valuable in scenarios with a multitude of inputs and insufficient data for the direct application of deep learning techniques, providing meaningful insights into the underlying processes.
△ Less
Submitted 13 May, 2024;
originally announced May 2024.
-
Quantum-informed simulations for mechanics of materials: DFTB+MBD framework
Authors:
Zhaoxiang Shen,
Raúl I. Sosa,
Stéphane P. A. Bordas,
Alexandre Tkatchenko,
Jakub Lengiewicz
Abstract:
The macroscopic behaviors of materials are determined by interactions that occur at multiple lengths and time scales. Depending on the application, describing, predicting, and understanding these behaviors require models that rely on insights from electronic and atomic scales. In such cases, classical simplified approximations at those scales are insufficient, and quantum-based modeling is require…
▽ More
The macroscopic behaviors of materials are determined by interactions that occur at multiple lengths and time scales. Depending on the application, describing, predicting, and understanding these behaviors require models that rely on insights from electronic and atomic scales. In such cases, classical simplified approximations at those scales are insufficient, and quantum-based modeling is required. In this paper, we study how quantum effects can modify the mechanical properties of systems relevant to materials engineering. We base our study on a high-fidelity modeling framework that combines two computationally efficient models rooted in quantum first principles: Density Functional Tight Binding (DFTB) and many-body dispersion (MBD). The MBD model is applied to accurately describe non-covalent van der Waals interactions. Through various benchmark applications, we demonstrate the capabilities of this framework and the limitations of simplified modeling. We provide an open-source repository containing all codes, datasets, and examples presented in this work. This repository serves as a practical toolkit that we hope will support the development of future research in effective large-scale and multiscale modeling with quantum-mechanical fidelity.
△ Less
Submitted 5 April, 2024;
originally announced April 2024.
-
Breast simulation pipeline: from medical imaging to patient-specific simulations
Authors:
Arnaud Mazier,
Stéphane P. A. Bordas
Abstract:
Breast-conserving surgery is the most acceptable operation for breast cancer removal from an invasive and psychological point of view. Before the surgical procedure, a preoperative MRI is performed in the prone configuration, while the surgery is achieved in the supine position. This leads to a considerable movement of the breast, including the tumor, between the two poses, complicating the surgeo…
▽ More
Breast-conserving surgery is the most acceptable operation for breast cancer removal from an invasive and psychological point of view. Before the surgical procedure, a preoperative MRI is performed in the prone configuration, while the surgery is achieved in the supine position. This leads to a considerable movement of the breast, including the tumor, between the two poses, complicating the surgeon's task. In this work, a simulation pipeline allowing the computation of patient-specific geometry and the prediction of personalized breast material properties was put forward. Through image segmentation, a finite element model including the subject-specific geometry is established. By first computing an undeformed state of the breast, the geometrico-material model is calibrated by surface acquisition in the intra-operative stance. Using an elastic corotational formulation, the patient-specific mechanical properties of the breast and skin were identified to obtain the best estimates of the supine configuration. The final results are a Mean Absolute Error of 4.00mm for the mechanical parameters E_breast = 0.32 kPa and E_skin = 22.72 kPa, congruent with the current state-of-the-art. The Covariance Matrix Adaptation Evolution Strategy optimizer converges on average between 5 to 30 min depending on the initial parameters, reaching a simulation speed of 20s. To our knowledge, our model offers one of the best compromises between accuracy and speed. Satisfactory results were obtained for the estimation of breast deformation from preoperative to intra-operative configuration. Furthermore, we have demonstrated the clinical feasibility of such applications using a simulation framework that aims at the smallest disturbance of the actual surgical pipeline.
△ Less
Submitted 6 March, 2023;
originally announced March 2023.
-
Level set topology optimization of metamaterial-based heat manipulators using isogeometric analysis
Authors:
Chintan Jansari,
Stéphane P. A. Bordas,
Elena Atroshchenko
Abstract:
We exploit level set topology optimization to find the optimal material distribution for metamaterial-based heat manipulators. The level set function, geometry, and solution field are parameterized using the non-uniform rational B-spline (NURBS) basis functions in order to take advantage of easy control of smoothness and continuity. In addition, NURBS approximations can produce conic geometries ex…
▽ More
We exploit level set topology optimization to find the optimal material distribution for metamaterial-based heat manipulators. The level set function, geometry, and solution field are parameterized using the non-uniform rational B-spline (NURBS) basis functions in order to take advantage of easy control of smoothness and continuity. In addition, NURBS approximations can produce conic geometries exactly and provide higher efficiency for higher-order elements. The values of the level set function at the control points (called expansion coefficients) are utilized as design variables. For optimization, we use an advanced mathematical programming technique, Sequential Quadratic Programming (SQP). Taking into account a large number of design variables and the small number of constraints associated with our optimization problem, the adjoint method is utilized to calculate the required sensitivities with respect to the design variables. The efficiency and robustness of the proposed method are demonstrated by solving three numerical examples. We have also shown that the current method can handle different geometries and types of objective functions. In addition, regularization techniques such as Tikhonov regularization and volume regularization have been explored to reduce unnecessary complexity and increase the manufacturability of optimized topologies.
△ Less
Submitted 6 March, 2023;
originally announced March 2023.
-
Convolution, aggregation and attention based deep neural networks for accelerating simulations in mechanics
Authors:
Saurabh Deshpande,
Raúl I. Sosa,
Stéphane P. A. Bordas,
Jakub Lengiewicz
Abstract:
Deep learning surrogate models are being increasingly used in accelerating scientific simulations as a replacement for costly conventional numerical techniques. However, their use remains a significant challenge when dealing with real-world complex examples. In this work, we demonstrate three types of neural network architectures for efficient learning of highly non-linear deformations of solid bo…
▽ More
Deep learning surrogate models are being increasingly used in accelerating scientific simulations as a replacement for costly conventional numerical techniques. However, their use remains a significant challenge when dealing with real-world complex examples. In this work, we demonstrate three types of neural network architectures for efficient learning of highly non-linear deformations of solid bodies. The first two architectures are based on the recently proposed CNN U-NET and MAgNET (graph U-NET) frameworks which have shown promising performance for learning on mesh-based data. The third architecture is Perceiver IO, a very recent architecture that belongs to the family of attention-based neural networks--a class that has revolutionised diverse engineering fields and is still unexplored in computational mechanics. We study and compare the performance of all three networks on two benchmark examples, and show their capabilities to accurately predict the non-linear mechanical responses of soft bodies.
△ Less
Submitted 24 March, 2023; v1 submitted 1 December, 2022;
originally announced December 2022.
-
Oncology and mechanics: landmark studies and promising clinical applications
Authors:
Stéphane Urcun,
Guillermo Lorenzo,
Davide Baroli,
Pierre-Yves Rohan,
Giuseppe Sciumè,
Wafa Skalli,
Vincent Lubrano,
Stéphane P. A. Bordas
Abstract:
Clinical management of cancer has continuously evolved for several decades. Biochemical, molecular and genomics approaches have brought and still bring numerous insights into cancerous diseases. It is now accepted that some phenomena, allowed by favorable biological conditions, emerge via mechanical signaling at the cellular scale and via mechanical forces at the macroscale. Mechanical phenomena i…
▽ More
Clinical management of cancer has continuously evolved for several decades. Biochemical, molecular and genomics approaches have brought and still bring numerous insights into cancerous diseases. It is now accepted that some phenomena, allowed by favorable biological conditions, emerge via mechanical signaling at the cellular scale and via mechanical forces at the macroscale. Mechanical phenomena in cancer have been studied in-depth over the last decades, and their clinical applications are starting to be understood. If numerous models and experimental setups have been proposed, only a few have led to clinical applications. The objective of this contribution is to propose to review a large scope of mechanical findings which have consequences on the clinical management of cancer. This review is mainly addressed to doctoral candidates in mechanics and applied mathematics who are faced with the challenge of the mechanics-based modeling of cancer with the aim of clinical applications. We show that the collaboration of the biological and mechanical approaches has led to promising advances in terms of modeling, experimental design and therapeutic targets. Additionally, a specific focus is brought on imaging-informed mechanics-based models, which we believe can further the development of new therapeutic targets and the advent of personalized medicine. We study in detail several successful workflows on patient-specific targeted therapies based on mechanistic modeling.
△ Less
Submitted 16 November, 2022;
originally announced November 2022.
-
MAgNET: A Graph U-Net Architecture for Mesh-Based Simulations
Authors:
Saurabh Deshpande,
Stéphane P. A. Bordas,
Jakub Lengiewicz
Abstract:
In many cutting-edge applications, high-fidelity computational models prove to be too slow for practical use and are therefore replaced by much faster surrogate models. Recently, deep learning techniques have increasingly been utilized to accelerate such predictions. To enable learning on large-dimensional and complex data, specific neural network architectures have been developed, including convo…
▽ More
In many cutting-edge applications, high-fidelity computational models prove to be too slow for practical use and are therefore replaced by much faster surrogate models. Recently, deep learning techniques have increasingly been utilized to accelerate such predictions. To enable learning on large-dimensional and complex data, specific neural network architectures have been developed, including convolutional and graph neural networks. In this work, we present a novel encoder-decoder geometric deep learning framework called MAgNET, which extends the well-known convolutional neural networks to accommodate arbitrary graph-structured data. MAgNET consists of innovative Multichannel Aggregation (MAg) layers and graph pooling/unpooling layers, forming a graph U-Net architecture that is analogous to convolutional U-Nets. We demonstrate the predictive capabilities of MAgNET in surrogate modeling for non-linear finite element simulations in the mechanics of solids.
△ Less
Submitted 2 April, 2024; v1 submitted 1 November, 2022;
originally announced November 2022.
-
SOniCS: Develop intuition on biomechanical systems through interactive error controlled simulations
Authors:
Arnaud Mazier,
Sidaty El Hadramy,
Jean-Nicolas Brunet,
Jack S. Hale,
Stéphane Cotin,
Stéphane P. A. Bordas
Abstract:
This new approach allows the user to experiment with model choices easily and quickly without requiring in-depth expertise, as constitutive models can be modified by one line of code only. This ease in building new models makes SOniCS ideal to develop surrogate, reduced order models and to train machine learning algorithms for uncertainty quantification or to enable patient-specific simulations. S…
▽ More
This new approach allows the user to experiment with model choices easily and quickly without requiring in-depth expertise, as constitutive models can be modified by one line of code only. This ease in building new models makes SOniCS ideal to develop surrogate, reduced order models and to train machine learning algorithms for uncertainty quantification or to enable patient-specific simulations. SOniCS is thus not only a tool that facilitates the development of surgical training simulations but also, and perhaps more importantly, paves the way to increase the intuition of users or otherwise non-intuitive behaviors of (bio)mechanical systems. The plugin uses new developments of the FEniCSx project enabling automatic generation with FFCx of finite element tensors such as the local residual vector and Jacobian matrix. We validate our approach with numerical simulations such as manufactured solutions, cantilever beams, and benchmarks provided by FEBio. We reach machine precision accuracy and demonstrate the use of the plugin for a real-time haptic simulation involving a surgical tool controlled by the user in contact with a hyperelastic liver. We include complete examples showing the use of our plugin for simulations involving Saint Venant-Kirchhoff, Neo-Hookean, Mooney-Rivlin, and Holzapfel Ogden anisotropic models as supplementary material.
△ Less
Submitted 24 August, 2022;
originally announced August 2022.
-
Design of metamaterial-based heat manipulators by isogeometric shape optimization
Authors:
Chintan Jansari,
Stéphane P. A. Bordas,
Elena Atroshchenko
Abstract:
There has been a growing interest in controlled heat flux manipulation to increase the efficiency of thermal apparatus. Heat manipulators control and manipulate heat flow. A key to the effective performance of these heat manipulators is their thermal design. Such designs can be achieved by a periodic assembly of unit cells (known as metamaterials or meta-structure), whose geometry and material pro…
▽ More
There has been a growing interest in controlled heat flux manipulation to increase the efficiency of thermal apparatus. Heat manipulators control and manipulate heat flow. A key to the effective performance of these heat manipulators is their thermal design. Such designs can be achieved by a periodic assembly of unit cells (known as metamaterials or meta-structure), whose geometry and material properties can be optimized for a specific objective. In this work, we focus on thermal metamaterial-based heat manipulators such as thermal concentrator (which concentrates the heat flux in a specified region of the domain). The main scope of the current work is to optimize the shape of the heat manipulators using Particle Swarm Optimization (PSO) method. The geometry is defined using NURBS basis functions due to the higher smoothness and continuity and the thermal boundary value problem is solved using Isogeometric Analysis (IGA). Often, nodes as design variables (as in Lagrange finite element method) generate the serrate shapes of boundaries which need to be smoothened later. For the NURBS-based boundary with the control points as design variables, the required smoothness can be predefined through knot vectors and smoothening in the post-processing can be avoided. The optimized shape generated by PSO is compared with the other shape exploited in the literature. The effects of the number of design variables, the thermal conductivity of the materials used, as well as some of the geometry parameters on the optimum shapes are also demonstrated.
△ Less
Submitted 18 May, 2022;
originally announced May 2022.
-
A graph-based probabilistic geometric deep learning framework with online enforcement of physical constraints to predict the criticality of defects in porous materials
Authors:
Vasilis Krokos,
Stéphane P. A. Bordas,
Pierre Kerfriden
Abstract:
Stress prediction in porous materials and structures is challenging due to the high computational cost associated with direct numerical simulations. Convolutional Neural Network (CNN) based architectures have recently been proposed as surrogates to approximate and extrapolate the solution of such multiscale simulations. These methodologies are usually limited to 2D problems due to the high computa…
▽ More
Stress prediction in porous materials and structures is challenging due to the high computational cost associated with direct numerical simulations. Convolutional Neural Network (CNN) based architectures have recently been proposed as surrogates to approximate and extrapolate the solution of such multiscale simulations. These methodologies are usually limited to 2D problems due to the high computational cost of 3D voxel based CNNs. We propose a novel geometric learning approach based on a Graph Neural Network (GNN) that efficiently deals with three-dimensional problems by performing convolutions over 2D surfaces only. Following our previous developments using pixel-based CNN, we train the GNN to automatically add local fine-scale stress corrections to an inexpensively computed coarse stress prediction in the porous structure of interest. Our method is Bayesian and generates densities of stress fields, from which credible intervals may be extracted. As a second scientific contribution, we propose to improve the extrapolation ability of our network by deploying a strategy of online physics-based corrections. Specifically, we condition the posterior predictions of our probabilistic predictions to satisfy partial equilibrium at the microscale, at the inference stage. This is done using an Ensemble Kalman algorithm, to ensure tractability of the Bayesian conditioning operation. We show that this innovative methodology allows us to alleviate the effect of undesirable biases observed in the outputs of the uncorrected GNN, and improves the accuracy of the predictions in general.
△ Less
Submitted 5 November, 2023; v1 submitted 13 May, 2022;
originally announced May 2022.
-
An adapted deflated conjugate gradient solver for robust extended/generalised finite element solutions of large scale, 3D crack propagation problems
Authors:
Konstantinos Agathos,
Tim Dodwell,
Eleni Chatzi,
Stephane P. A. Bordas
Abstract:
An adapted deflation preconditioner is employed to accelerate the solution of linear systems resulting from the discretization of fracture mechanics problems with well-conditioned extended/generalized finite elements. The deflation space typically used for linear elasticity problems is enriched with additional vectors, accounting for the enrichment functions used, thus effectively removing low fre…
▽ More
An adapted deflation preconditioner is employed to accelerate the solution of linear systems resulting from the discretization of fracture mechanics problems with well-conditioned extended/generalized finite elements. The deflation space typically used for linear elasticity problems is enriched with additional vectors, accounting for the enrichment functions used, thus effectively removing low frequency components of the error. To further improve performance, deflation is combined, in a multiplicative way, with a block-Jacobi preconditioner, which removes high frequency components of the error as well as linear dependencies introduced by enrichment. The resulting scheme is tested on a series of non-planar crack propagation problems and compared to alternative linear solvers in terms of performance.
△ Less
Submitted 17 November, 2021;
originally announced November 2021.
-
Probabilistic Deep Learning for Real-Time Large Deformation Simulations
Authors:
Saurabh Deshpande,
Jakub Lengiewicz,
Stéphane P. A. Bordas
Abstract:
For many novel applications, such as patient-specific computer-aided surgery, conventional solution techniques of the underlying nonlinear problems are usually computationally too expensive and are lacking information about how certain can we be about their predictions. In the present work, we propose a highly efficient deep-learning surrogate framework that is able to accurately predict the respo…
▽ More
For many novel applications, such as patient-specific computer-aided surgery, conventional solution techniques of the underlying nonlinear problems are usually computationally too expensive and are lacking information about how certain can we be about their predictions. In the present work, we propose a highly efficient deep-learning surrogate framework that is able to accurately predict the response of bodies undergoing large deformations in real-time. The surrogate model has a convolutional neural network architecture, called U-Net, which is trained with force-displacement data obtained with the finite element method. We propose deterministic and probabilistic versions of the framework. The probabilistic framework utilizes the Variational Bayes Inference approach and is able to capture all the uncertainties present in the data as well as in the deep-learning model. Based on several benchmark examples, we show the predictive capabilities of the framework and discuss its possible limitations
△ Less
Submitted 4 July, 2022; v1 submitted 2 November, 2021;
originally announced November 2021.
-
Neural-network acceleration of projection-based model-order-reduction for finite plasticity: Application to RVEs
Authors:
S. Vijayaraghavan,
L. Wu,
L. Noels,
S. P. A. Bordas,
S. Natarajan,
L. A. A. Beex
Abstract:
Compared to conventional projection-based model-order-reduction, its neural-network acceleration has the advantage that the online simulations are equation-free, meaning that no system of equations needs to be solved iteratively. Consequently, no stiffness matrix needs to be constructed and the stress update needs to be computed only once per increment. In this contribution, a recurrent neural net…
▽ More
Compared to conventional projection-based model-order-reduction, its neural-network acceleration has the advantage that the online simulations are equation-free, meaning that no system of equations needs to be solved iteratively. Consequently, no stiffness matrix needs to be constructed and the stress update needs to be computed only once per increment. In this contribution, a recurrent neural network is developed to accelerate a projection-based model-order-reduction of the elastoplastic mechanical behaviour of an RVE. In contrast to a neural network that merely emulates the relation between the macroscopic deformation (path) and the macroscopic stress, the neural network acceleration of projection-based model-order-reduction preserves all microstructural information, at the price of computing this information once per increment.
△ Less
Submitted 16 September, 2021;
originally announced September 2021.
-
Machine learning in the social and health sciences
Authors:
Anja K. Leist,
Matthias Klee,
Jung Hyun Kim,
David H. Rehkopf,
Stéphane P. A. Bordas,
Graciela Muniz-Terrera,
Sara Wade
Abstract:
The uptake of machine learning (ML) approaches in the social and health sciences has been rather slow, and research using ML for social and health research questions remains fragmented. This may be due to the separate development of research in the computational/data versus social and health sciences as well as a lack of accessible overviews and adequate training in ML techniques for non data scie…
▽ More
The uptake of machine learning (ML) approaches in the social and health sciences has been rather slow, and research using ML for social and health research questions remains fragmented. This may be due to the separate development of research in the computational/data versus social and health sciences as well as a lack of accessible overviews and adequate training in ML techniques for non data science researchers. This paper provides a meta-mapping of research questions in the social and health sciences to appropriate ML approaches, by incorporating the necessary requirements to statistical analysis in these disciplines. We map the established classification into description, prediction, and causal inference to common research goals, such as estimating prevalence of adverse health or social outcomes, predicting the risk of an event, and identifying risk factors or causes of adverse outcomes. This meta-mapping aims at overcoming disciplinary barriers and starting a fluid dialogue between researchers from the social and health sciences and methodologically trained researchers. Such mapping may also help to fully exploit the benefits of ML while considering domain-specific aspects relevant to the social and health sciences, and hopefully contribute to the acceleration of the uptake of ML applications to advance both basic and applied social and health sciences research.
△ Less
Submitted 20 June, 2021;
originally announced June 2021.
-
Circumferential Crack Modeling of Thin Cylindrical Shells in Modal Deformation
Authors:
Ali Alijani,
Olga Barrera,
Stephane P. A. Bordas
Abstract:
An innovative technique, called conversion, is introduced to model circumferential cracks in thin cylindrical shells. The semi-analytical finite element method is applied to investigate the modal deformation of the cylinder. An element including the crack is divided into three sub-elements with four nodes in which the stiffness matrix is enriched. The crack characteristics are included in the fini…
▽ More
An innovative technique, called conversion, is introduced to model circumferential cracks in thin cylindrical shells. The semi-analytical finite element method is applied to investigate the modal deformation of the cylinder. An element including the crack is divided into three sub-elements with four nodes in which the stiffness matrix is enriched. The crack characteristics are included in the finite element method relations through conversion matrices and a rotational spring corresponding to the crack. Conversion matrices obtained by applying continuity conditions at the crack tip are used to transform displacements of the middle nodes to those of the main nodes. Moreover, another technique, called spring set, is represented based on a set of springs to model the crack as a separated element. Components of the stiffness matrix related to the separated element are incorporated while the geometric boundary conditions at the crack tip are satisfied. The effects of the circumferential mode number, the crack depth and the length of the cylinder on the critical buckling load are investigated. Experimental tests, ABAQUS modeling and results from literature are used to verify and validate the results and derived relations. In addition, the crack effect on the natural frequency is examined using the vibration analysis based on the conversion technique.
△ Less
Submitted 15 May, 2021;
originally announced May 2021.
-
A rigged model of the breast for preoperative surgical planning
Authors:
Arnaud Mazier,
Sophie Ribes,
Benjamin Gilles,
Stéphane P. A Bordas,
.
Abstract:
In breast surgical practice, drawing is part of the preoperative planning procedure and is essential for a successful operation. In this study, we design a pipeline to assist surgeons with patient-specific breast surgical drawings. We use a deformable torso model containing the surgical patterns to match any breast surface scan. To be compatible with surgical timing, we build an articulated model…
▽ More
In breast surgical practice, drawing is part of the preoperative planning procedure and is essential for a successful operation. In this study, we design a pipeline to assist surgeons with patient-specific breast surgical drawings. We use a deformable torso model containing the surgical patterns to match any breast surface scan. To be compatible with surgical timing, we build an articulated model through a skinning process coupled with shape deformers to enhance a fast registration process. On one hand, the scalable bones of the skinning account for pose and morphological variations of the patients. On the other hand, pre-designed artistic blendshapes create a linear space for guaranteeing anatomical variations. Then, we apply meaningful constraints to the model to find a trade-off between precision and speed. The experiments were conducted on 7 patients, in 2 different poses (prone and supine) with a breast size ranging from 36A and 42C (US/UK bra sizing). The acquisitions were obtained using the depth camera Structure Sensor, and the breast scans were acquired in less than 1 minute. The result is a registration method converging within a few seconds (3 maximum), reaching a Mean Absolute Error of 2.3 mm for mesh registration and 8.0 mm for breast anatomical landmarks. Compared to the existing literature, our model can be personalized and does not require any database. Finally, our registered model can be used to transfer surgical reference patterns onto any patient in any position.
△ Less
Submitted 31 May, 2022; v1 submitted 3 May, 2021;
originally announced May 2021.
-
Inverse deformation analysis: an experimental and numerical assessment using the FEniCS Project
Authors:
Arnaud Mazier,
Alexandre Bilger,
Antonio E. Forte,
Igor Peterlik,
Jack S. Hale,
Stéphane P. A. Bordas,
.,
Institute of Computational Engineering,
Department of Engineering,
University of Luxembourg,
Esch-sur-Alzette,
Luxembourg.,
Harvard University,
Cambridge,
USA.,
Department of Electronics,
Information,
Bioengineering,
Politecnico di Milano,
Milan,
Italy.,
Institute of Computer Science,
Masaryk University,
Czech Republic.,
Institute of Research
, et al. (3 additional authors not shown)
Abstract:
In this paper, we develop a framework for solving inverse deformation problems using the FEniCS Project finite element software. We validate our approach with experimental imaging data acquired from a soft silicone beam under gravity. In contrast with inverse iterative algorithms that require multiple solutions of a standard elasticity problem, the proposed method can compute the undeformed config…
▽ More
In this paper, we develop a framework for solving inverse deformation problems using the FEniCS Project finite element software. We validate our approach with experimental imaging data acquired from a soft silicone beam under gravity. In contrast with inverse iterative algorithms that require multiple solutions of a standard elasticity problem, the proposed method can compute the undeformed configuration by solving only one modified elasticity problem. This modified problem has a complexity comparable to the standard one. The framework is implemented within an open-source pipeline enabling the direct and inverse deformation simulation directly from imaging data. We use the high-level Unified Form Language (UFL) of the FEniCS Project to express the finite element model in variational form and to automatically derive the consistent Jacobian. Consequently, the design of the pipeline is flexible: for example, it allows the modification of the constitutive models by changing a single line of code. We include a complete working example showing the inverse deformation of a beam deformed by gravity as supplementary material.
△ Less
Submitted 26 February, 2021;
originally announced February 2021.
-
Hierarchical a posteriori error estimation of Bank-Weiser type in the FEniCS Project
Authors:
Raphaël Bulle,
Jack S. Hale,
Alexei Lozinski,
Stéphane P. A. Bordas,
Franz Chouly
Abstract:
In the seminal paper of Bank and Weiser [Math. Comp., 44 (1985), pp.283-301] a new a posteriori estimator was introduced. This estimator requires the solution of a local Neumann problem on every cell of the finite element mesh. Despite the promise of Bank-Weiser type estimators, namely locality, computational efficiency, and asymptotic sharpness, they have seen little use in practical computationa…
▽ More
In the seminal paper of Bank and Weiser [Math. Comp., 44 (1985), pp.283-301] a new a posteriori estimator was introduced. This estimator requires the solution of a local Neumann problem on every cell of the finite element mesh. Despite the promise of Bank-Weiser type estimators, namely locality, computational efficiency, and asymptotic sharpness, they have seen little use in practical computational problems. The focus of this contribution is to describe a novel implementation of hierarchical estimators of the Bank-Weiser type in a modern high-level finite element software with automatic code generation capabilities. We show how to use the estimator to drive (goal-oriented) adaptive mesh refinement and to mixed approximations of the nearly-incompressible elasticity problems. We provide comparisons with various other used estimators. An open-source implementation based on the FEniCS Project finite element software is provided as supplementary material.
△ Less
Submitted 7 February, 2022; v1 submitted 8 February, 2021;
originally announced February 2021.
-
A Bayesian multiscale CNN framework to predict local stress fields in structures with microscale features
Authors:
Vasilis Krokos,
Viet Bui Xuan,
Stéphane P. A. Bordas,
Philippe Young,
Pierre Kerfriden
Abstract:
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…
▽ More
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.
△ Less
Submitted 19 January, 2022; v1 submitted 17 December, 2020;
originally announced December 2020.
-
Distributed prediction of unsafe reconfiguration scenarios of modular robotic Programmable Matter
Authors:
Benoît Piranda,
Paweł Chodkiewicz,
Paweł Hołobut,
Stéphane P. A. Bordas,
Julien Bourgeois,
Jakub Lengiewicz
Abstract:
We present a distributed framework for predicting whether a planned reconfiguration step of a modular robot will mechanically overload the structure, causing it to break or lose stability under its own weight. The algorithm is executed by the modular robot itself and based on a distributed iterative solution of mechanical equilibrium equations derived from a simplified model of the robot. The mode…
▽ More
We present a distributed framework for predicting whether a planned reconfiguration step of a modular robot will mechanically overload the structure, causing it to break or lose stability under its own weight. The algorithm is executed by the modular robot itself and based on a distributed iterative solution of mechanical equilibrium equations derived from a simplified model of the robot. The model treats inter-modular connections as beams and assumes no-sliding contact between the modules and the ground. We also provide a procedure for simplified instability detection. The algorithm is verified in the Programmable Matter simulator VisibleSim, and in real-life experiments on the modular robotic system Blinky Blocks.
△ Less
Submitted 3 May, 2021; v1 submitted 19 June, 2020;
originally announced June 2020.
-
Quantifying discretization errors for soft-tissue simulation in computer assisted surgery: a preliminary study
Authors:
Michel Duprez,
Stéphane P. A. Bordas,
Marek Bucki,
Huu Phuoc Bui,
Franz Chouly,
Vanessa Lleras,
Claudio Lobos,
Alexei Lozinski,
Pierre-Yves Rohan,
Satyendra Tomar
Abstract:
Errors in biomechanics simulations arise from modeling and discretization. Modeling errors are due to the choice of the mathematical model whilst discretization errors measure the impact of the choice of the numerical method on the accuracy of the approximated solution to this specific mathematical model. A major source of discretization errors is mesh generation from medical images, that remains…
▽ More
Errors in biomechanics simulations arise from modeling and discretization. Modeling errors are due to the choice of the mathematical model whilst discretization errors measure the impact of the choice of the numerical method on the accuracy of the approximated solution to this specific mathematical model. A major source of discretization errors is mesh generation from medical images, that remains one of the major bottlenecks in the development of reliable, accurate, automatic and efficient personalized, clinically-relevant Finite Element (FE) models in biomechanics. The impact of mesh quality and density on the accuracy of the FE solution can be quantified with \emph{a posteriori} error estimates. Yet, to our knowledge, the relevance of such error estimates for practical biomechanics problems has seldom been addressed, see [25]. In this contribution, we propose an implementation of some a posteriori error estimates to quantify the discretization errors and to optimize the mesh. More precisely, we focus on error estimation for a user-defined quantity of interest with the Dual Weighted Residual (DWR) technique. We test its applicability and relevance in two situations, corresponding to computations for a tongue and an artery, using a simplified setting, i.e., plane linearized elasticity with contractility of the soft-tissue modeled as a pre-stress. Our results demonstrate the feasibility of such methodology to estimate the actual solution errors and to reduce them economically through mesh refinement.
△ Less
Submitted 18 June, 2018;
originally announced June 2018.
-
Corotational Cut Finite Element Method for real-time surgical simulation: application to needle insertion simulation
Authors:
Huu Phuoc Bui,
Satyendra Tomar,
Stéphane P. A. Bordas
Abstract:
This paper describes the use of the corotational cut Finite Element Method (FEM) for real-time surgical simulation. Users only need to provide a background mesh which is not necessarily conforming to the boundaries/interfaces of the simulated object. The details of the surface, which can be directly obtained from binary images, are taken into account by a multilevel embedding algorithm applied to…
▽ More
This paper describes the use of the corotational cut Finite Element Method (FEM) for real-time surgical simulation. Users only need to provide a background mesh which is not necessarily conforming to the boundaries/interfaces of the simulated object. The details of the surface, which can be directly obtained from binary images, are taken into account by a multilevel embedding algorithm applied to elements of the background mesh that cut by the surface. Boundary conditions can be implicitly imposed on the surface using Lagrange multipliers. The implementation is verified by convergence studies with optimal rates. The algorithm is applied to various needle insertion simulations (e.g. for biopsy or brachytherapy) into brain and liver to verify the reliability of method, and numerical results show that the present method can make the discretisation independent from geometric description, and can avoid the complexity of mesh generation of complex geometries while retaining the accuracy of the standard FEM. Using the proposed approach is very suitable for real-time and patient specific simulations as it improves the simulation accuracy by taking into account automatically and properly the simulated geometry.
△ Less
Submitted 8 December, 2017;
originally announced December 2017.
-
Isogeometric analysis of thin Reissner-Mindlin plates and shells: locking phenomena and B-bar method
Authors:
Qingyuan Hu,
Yang Xia,
Sundararajan Natarajan,
Andreas Zilian,
Ping Hu,
Stéphane P. A. Bordas
Abstract:
We propose a local type of B-bar formulation, addressing locking in degenerated Reissner-Mindlin plate and shell formulations in the context of isogeometric analysis. Parasitic strain components are projected onto the physical space locally, i.e. at the element level, using a least-squares approach. The formulation allows the flexible utilization of basis functions of different order as the projec…
▽ More
We propose a local type of B-bar formulation, addressing locking in degenerated Reissner-Mindlin plate and shell formulations in the context of isogeometric analysis. Parasitic strain components are projected onto the physical space locally, i.e. at the element level, using a least-squares approach. The formulation allows the flexible utilization of basis functions of different order as the projection bases. The present formulation is much cheaper computationally than the classical $\bar{B}$ method. We show the numerical consistency of the scheme through numerical examples, moreover they show that the proposed formulation alleviates locking and yields good accuracy even for slenderness ratios of $1 \times 10^5$, and has the ability to capture deformations of thin shells using relatively coarse meshes. In addition it can be opined that the proposed method is less sensitive to locking and mesh distortion.
△ Less
Submitted 22 February, 2018; v1 submitted 1 September, 2017;
originally announced September 2017.
-
Constructing IGA-suitable planar parameterization from complex CAD boundary by domain partition and global/local optimization
Authors:
Gang Xu,
Ming Li,
Bernard Mourrain,
Timon Rabczuk,
Jinlan Xu,
Stephane P. A. Bordas
Abstract:
In this paper, we propose a general framework for constructing IGA-suitable planar B-spline parameterizations from given complex CAD boundaries consisting of a set of B-spline curves. Instead of forming the computational domain by a simple boundary, planar domains with high genus and more complex boundary curves are considered. Firstly, some pre-processing operations including Bézier extraction an…
▽ More
In this paper, we propose a general framework for constructing IGA-suitable planar B-spline parameterizations from given complex CAD boundaries consisting of a set of B-spline curves. Instead of forming the computational domain by a simple boundary, planar domains with high genus and more complex boundary curves are considered. Firstly, some pre-processing operations including Bézier extraction and subdivision are performed on each boundary curve in order to generate a high-quality planar parameterization; then a robust planar domain partition framework is proposed to construct high-quality patch-meshing results with few singularities from the discrete boundary formed by connecting the end points of the resulting boundary segments. After the topology information generation of quadrilateral decomposition, the optimal placement of interior Bézier curves corresponding to the interior edges of the quadrangulation is constructed by a global optimization method to achieve a patch-partition with high quality. Finally, after the imposition of C1=G1-continuity constraints on the interface of neighboring Bézier patches with respect to each quad in the quadrangulation, the high-quality Bézier patch parameterization is obtained by a C1-constrained local optimization method to achieve uniform and orthogonal iso-parametric structures while keeping the continuity conditions between patches. The efficiency and robustness of the proposed method are demonstrated by several examples which are compared to results obtained by the skeleton-based parameterization approach.
△ Less
Submitted 8 July, 2017; v1 submitted 3 July, 2017;
originally announced July 2017.
-
Controlling the Error on Target Motion through Real-time Mesh Adaptation: Applications to Deep Brain Stimulation
Authors:
Huu Phuoc Bui,
Satyendra Tomar,
Hadrien Courtecuisse,
Michel Audette,
Stéphane Cotin,
Stéphane P. A. Bordas
Abstract:
We present an error-controlled mesh refinement procedure for needle insertion simulation and apply it to the simulation of electrode implantation for deep brain stimulation, including brain shift. Our approach enables to control the error in the computation of the displacement and stress fields around the needle tip and needle shaft by suitably refining the mesh, whilst maintaining a coarser mesh…
▽ More
We present an error-controlled mesh refinement procedure for needle insertion simulation and apply it to the simulation of electrode implantation for deep brain stimulation, including brain shift. Our approach enables to control the error in the computation of the displacement and stress fields around the needle tip and needle shaft by suitably refining the mesh, whilst maintaining a coarser mesh in other parts of the domain. We demonstrate through academic and practical examples that our approach increases the accuracy of the displacement and stress fields around the needle without increasing the computational expense. This enables real-time simulations. The proposed methodology has direct implications to increase the accuracy and control the computational expense of the simulation of percutaneous procedures such as biopsy, brachytherapy, regional anesthesia, or cryotherapy and can be essential to the development of robotic guidance.
△ Less
Submitted 30 September, 2017; v1 submitted 25 April, 2017;
originally announced April 2017.
-
Bayesian inference for the stochastic identification of elastoplastic material parameters: Introduction, misconceptions and insights
Authors:
Hussein Rappel,
Lars A. A. Beex,
Jack S. Hale,
Stephane P. A. Bordas
Abstract:
We discuss Bayesian inference (BI) for the probabilistic identification of material parameters. This contribution aims to shed light on the use of BI for the identification of elastoplastic material parameters. For this purpose a single spring is considered, for which the stress-strain curves are artificially created. Besides offering a didactic introduction to BI, this paper proposes an approach…
▽ More
We discuss Bayesian inference (BI) for the probabilistic identification of material parameters. This contribution aims to shed light on the use of BI for the identification of elastoplastic material parameters. For this purpose a single spring is considered, for which the stress-strain curves are artificially created. Besides offering a didactic introduction to BI, this paper proposes an approach to incorporate statistical errors both in the measured stresses, and in the measured strains. It is assumed that the uncertainty is only due to measurement errors and the material is homogeneous. Furthermore, a number of possible misconceptions on BI are highlighted based on the purely elastic case.
△ Less
Submitted 15 January, 2017; v1 submitted 8 June, 2016;
originally announced June 2016.
-
Isogeometric analysis for functionally graded microplates based on modified couple stress theory
Authors:
Hoang X. Nguyen,
Tuan N. Nguyen,
M. Abdel-Wahab,
S. P. A. Bordas,
H. Nguyen-Xuan,
Thuc P. Vo
Abstract:
Analysis of static bending, free vibration and buckling behaviours of functionally graded microplates is investigated in this study. The main idea is to use the isogeometric analysis in associated with novel four-variable refined plate theory and quasi-3D theory. More importantly, the modified couple stress theory with only one material length scale parameter is employed to effectively capture the…
▽ More
Analysis of static bending, free vibration and buckling behaviours of functionally graded microplates is investigated in this study. The main idea is to use the isogeometric analysis in associated with novel four-variable refined plate theory and quasi-3D theory. More importantly, the modified couple stress theory with only one material length scale parameter is employed to effectively capture the size-dependent effects within the microplates. Meanwhile, the quasi-3D theory which is constructed from a novel seventh-order shear deformation refined plate theory with four unknowns is able to consider both shear deformations and thickness stretching effect without requiring shear correction factors. The NURBS-based isogeometric analysis is integrated to exactly describe the geometry and approximately calculate the unknown fields with higher-order derivative and continuity requirements. The convergence and verification show the validity and efficiency of this proposed computational approach in comparison with those existing in the literature. It is further applied to study the static bending, free vibration and buckling responses of rectangular and circular functionally graded microplates with various types of boundary conditions. A number of investigations are also conducted to illustrate the effects of the material length scale, material index, and length-to-thickness ratios on the responses of the microplates.
△ Less
Submitted 2 April, 2016;
originally announced April 2016.
-
Isogeometric analysis: an overview and computer implementation aspects
Authors:
Vinh Phu Nguyen,
Stéphane P. A. Bordas,
Timon Rabczuk
Abstract:
Isogeometric analysis (IGA) represents a recently developed technology in computational mechanics that offers the possibility of integrating methods for analysis and Computer Aided Design (CAD) into a single, unified process. The implications to practical engineering design scenarios are profound, since the time taken from design to analysis is greatly reduced, leading to dramatic gains in efficie…
▽ More
Isogeometric analysis (IGA) represents a recently developed technology in computational mechanics that offers the possibility of integrating methods for analysis and Computer Aided Design (CAD) into a single, unified process. The implications to practical engineering design scenarios are profound, since the time taken from design to analysis is greatly reduced, leading to dramatic gains in efficiency. The tight coupling of CAD and analysis within IGA requires knowledge from both fields and it is one of the goals of the present paper to outline much of the commonly used notation. In this manuscript, through a clear and simple Matlab implementation, we present an introduction to IGA applied to the Finite Element (FE) method and related computer implementation aspects. Furthermore, implemen- tation of the extended IGA which incorporates enrichment functions through the partition of unity method (PUM) is also presented, where several examples for both two-dimensional and three-dimensional fracture are illustrated. The open source Matlab code which accompanies the present paper can be applied to one, two and three-dimensional problems for linear elasticity, linear elastic fracture mechanics, structural mechanics (beams/plates/shells including large displacements and rotations) and Poisson problems with or without enrichment. The Bezier extraction concept that allows FE analysis to be performed efficiently on T-spline geometries is also incorporated. The article includes a summary of recent trends and developments within the field of IGA.
△ Less
Submitted 27 September, 2013; v1 submitted 9 May, 2012;
originally announced May 2012.
-
Efficient recovery-based error estimation for the smoothed finite element method for smooth and singular linear elasticity
Authors:
Octavio A. González-Estrada,
Sundararajan Natarajan,
Juan José Ródenas,
Hung Nguyen-Xuan,
Stéphane P. A. Bordas
Abstract:
An error control technique aimed to assess the quality of smoothed finite element approximations is presented in this paper. Finite element techniques based on strain smoothing appeared in 2007 were shown to provide significant advantages compared to conventional finite element approximations. In particular, a widely cited strength of such methods is improved accuracy for the same computational co…
▽ More
An error control technique aimed to assess the quality of smoothed finite element approximations is presented in this paper. Finite element techniques based on strain smoothing appeared in 2007 were shown to provide significant advantages compared to conventional finite element approximations. In particular, a widely cited strength of such methods is improved accuracy for the same computational cost. Yet, few attempts have been made to directly assess the quality of the results obtained during the simulation by evaluating an estimate of the discretization error. Here we propose a recovery type error estimator based on an enhanced recovery technique. The salient features of the recovery are: enforcement of local equilibrium and, for singular problems a "smooth+singular" decomposition of the recovered stress. We evaluate the proposed estimator on a number of test cases from linear elastic structural mechanics and obtain precise error estimations whose effectivities, both at local and global levels, are improved compared to recovery procedures not implementing these features.
△ Less
Submitted 6 March, 2012;
originally announced March 2012.