An energy–momentum conserving temporal integration scheme is presented for a recently proposed geometrically exact shell formulation with drilling degrees of freedom. The scheme is based on a novel idea of defining mixed discrete derivatives for holonomic constraint functions with displacements and rotations. By defining general discrete derivative expressions with unknown terms, the mixed discrete

Phase-field modeling of fracture has gained popularity within the last decade due to the flexibility of the related computational framework in simulating three-dimensional arbitrarily complicated fracture processes. However, the numerical predictions are greatly affected by the presence of uncertainties in the mechanical properties of the material originating from unresolved heterogeneities and the

Grain boundary strengthening effect for polycrystalline copper is studied numerically using crystal plasticity in conjunction with cylindrical indentation simulations under the plane strain condition. In order to compare with an isotropic, heterogeneous continuum model a new constitutive relation is developed. This new nonlocal continuum model that encompasses the heterogeneity in yield strength based

In this work, we develop a dynamic hybrid local/nonlocal continuum model to study wave propagations in a linear elastic solid. The developed hybrid model couples, in the dynamic regime, a classical continuum mechanics model with a bond-based peridynamic model using the Morphing coupling method that introduced in a previous study (Lubineau et al., J Mech Phys Solids 60(6):1088–1102, 2012). The classical

This paper investigates the dynamic characteristics of a multilevel structure for the transportation and storage of spent nuclear fuel (SNF) from commercial power plants. The nuclear fuel is stored in slender rods that are grouped together into fuel assemblies (FA). In a sealed cylindrical container called “canister”, the FA are inserted into a honeycomb basket. The objective of this paper is to develop

Light structural metals have been extensively applied throughout the transportation sector in recent years with greater impetus to reduce vehicle weight and enhance energy efficiency. However, their use is restricted by the engineering challenge of fatigue crack growth and the need to understand, simulate, and predict crack propagation mechanisms with respect to materials’ microstructure. To address

A class of isothermal dissipative continua being often referred to as “standard dissipative” is considered. The initial boundary value problem describing the behavior of these continua can be conveniently formulated in terms of a Helmholtz free energy functional, a dissipation functional, a power functional, and, possibly, a Lagrangian multiplier functional. In order to obtain (approximate) solutions

Statistical Volume Elements (SVEs) are employed to evaluate homogenized mesoscopic ductile failure response of a carbon nanofiber reinforced composite under uniaxial tensile and compressive loadings. In the mesoscale analysis, after virtual reconstruction of the material microstructure, 2D finite element models are generated for each SVE using a non-iterative meshing algorithm named CISAMR, which fully

The paper deals with a numerical model to investigate the influence of stress state on damage and failure in the ductile steel X5CrNi18-10. The numerical analysis is based on an anisotropic continuum damage model taking into account yield and damage criteria as well as evolution equations for plastic and damage strain rate tensors. Results of numerical simulations of biaxial experiments with the X0-

The hierarchical deep-learning neural network (HiDeNN) is systematically developed through the construction of structured deep neural networks (DNNs) in a hierarchical manner, and a special case of HiDeNN for representing Finite Element Method (or HiDeNN-FEM in short) is established. In HiDeNN-FEM, weights and biases are functions of the nodal positions, hence the training process in HiDeNN-FEM includes

The paper presents an assumed strain formulation over polygonal meshes to accurately evaluate the strain fields in nonlocal damage models. An assume strained technique based on the Hu-Washizu variational principle is employed to generate a new strain approximation instead of direct derivation from the basis functions and the displacement fields. The underlying idea embedded in arbitrary finite polygons

Neural network constitutive models (NNCMs) for crystal structures are proposed based on computationally generated high-fidelity data. Stress, and tangent modulus data are generated under various strain states using empirical potentials and first-principles calculations. Strain–stress artificial neural network and strain-tangent modulus ANN are constructed. The symmetry conditions are considered for

This paper presents an adaptive strategy for phase-field simulations with transition to fracture. The phase-field equations are solved only in small subdomains around crack tips to determine propagation, while an extended finite element method (XFEM) discretization is used in the rest of the domain to represent sharp cracks, enabling to use a coarser discretization and therefore reducing the computational

We demonstrate a Bayesian method for the “real-time” characterization and forecasting of partially observed COVID-19 epidemic. Characterization is the estimation of infection spread parameters using daily counts of symptomatic patients. The method is designed to help guide medical resource allocation in the early epoch of the outbreak. The estimation problem is posed as one of Bayesian inference and

The standard Yee FDTD algorithm is widely used in computational electromagnetics because of its simplicity and divergence free nature. A generalization of this classical scheme to 3D unstructured co-volume meshes is adopted, based on the use of a Delaunay primal mesh and its high quality Voronoi dual. This circumvents the problem of accuracy losses, which are normally associated with the use of a staircased

The original article was published with errors in some sentences. The correct sentences are provided in this correction.

Fully coupled global equations are proposed for enhancing the performance of Finite Element analysis of unsaturated soils. The governing equation describing mechanical equilibrium is formulated in terms of net stress, and in the mass conservation equation the contribution of this net stress in determining the change of degree of saturation is also included. The novelty of this paper is the development

In this work we present a novel computational method for embedding arbitrary curved one-dimensional (1D) fibers into three-dimensional (3D) solid volumes, as e.g. in fiber-reinforced materials. The fibers are explicitly modeled with highly efficient 1D geometrically exact beam finite elements, based on various types of geometrically nonlinear beam theories. The surrounding solid volume is modeled with

Thermo-mechanical finite element (FE) models predict the thermal behavior of machine tools and the associated mechanical deviations. However, one disadvantage is their high computational expense, linked to the evaluation of the large systems of differential equations. Therefore, projection-based model order reduction (MOR) methods are required in order to create efficient surrogate models. This paper

The surface-to-surface master–master contact treatment is a technique that addresses pointwise contact between bodies with no prior election of slave points, as in master–slave case. For a given configuration of contact-candidate surfaces, one needs to find the material points associated with a pointwise contact interaction. This is the local contact problem (LCP). The methodology can be applied together

Diffusion-type problems in (nearly) unbounded domains play important roles in various fields of fluid dynamics, biology, and materials science. The aim of this paper is to construct accurate absorbing boundary conditions (ABCs) suitable for classical (local) as well as nonlocal peridynamic (PD) diffusion models. The main focus of the present study is on the PD diffusion formulation. The majority of

The processes in which cardiac cells are reorganized for tissue regeneration is still unclear. It is a complicated process that is orchestrated by many factors such as mechanical, chemical, thermal, and/or electrical cues. Studying and optimizing these conditions in-vitro is complicated and time costly. In such cases, in-silico numerical simulations can offer a reliable solution to predict and optimize

The paper deals with the gradient based shape optimization of the biaxial X0-specimen, which has been introduced and examined in various papers, under producibility restrictions and the related experimental verification. The original, engineering based design of the X0-specimen has been applied successfully to different loading conditions persisting the question if relevant stress states could be reached

A recently introduced NURBS mesh generation method for complex-geometry Isogeometric Analysis (IGA) is applied to building a high-quality mesh for a gas turbine. The compressible flow in the turbine is computed using the IGA and a stabilized method with improved discontinuity-capturing, weakly-enforced no-slip boundary-condition, and sliding-interface operators. The IGA results are compared with the

This paper presents a Lagrangian approach to simulating multibody dynamics in a tensegrity framework with an ability to tackle holonomic constraint violations in an energy-preserving scheme. Governing equations are described using non-minimum coordinates to simplify descriptions of the structure’s kinematics. To minimize constraint drift arising from this redundant system, the direct correction method

Magnetostrictive materials that couple mechanical and magnetic domains have been widely explored for use in sensors and actuators. These materials often exhibit a nonlinear material response under applied magnetic fields, which limits the use of linear constitutive models. Furthermore, the nonlinear constitutive relations tend to be implicit in nature. Hence, a finite element scheme that can handle

Phase-field models for brittle fracture in anisotropic materials result in a fourth-order partial differential equation for the damage evolution. This necessitates a \(\mathcal {C}^1\) continuity of the basis functions. Here, Powell-Sabin B-splines, which are based on triangles, are used for the approximation of the field variables as well as for the the description of the geometry. The use of triangles

In the original publication of the article, the sentence “If we assume that everyone in the beginning is susceptible (S(t = 0) = 1).

Fractional diffusion equations have been widely used to accurately describe anomalous solute transport in complex media. This paper proposes a local meshless method named the generalized finite difference method (GFDM), to solve a class of multidimensional space fractional diffusion equations (SFDEs) in a finite domain. In the GFDM, the spatial derivative terms are expressed as linear combinations

In this paper we alleviate limitations of the conventional embedded reinforcement formulation for applications with dense fiber contents. We demonstrate that by condensing the fiber degrees of freedom during the assembly stage, the condition number of the resultant stiffness matrix is effectively reduced and the use of iterative solvers is facilitated even without preconditioning. Numerical benchmarks

A large number of massive crystal-plasticity-finite-element (CPFE) simulations are performed and post-processed to reveal the effects of element type and mesh resolution on accuracy of predicted mechanical fields over explicit grain structures. A CPFE model coupled with Abaqus/Standard is used to simulate simple-tension and simple-shear deformations to facilitate such quantitative mesh sensitivity

In this paper layered shells subjected to static loading are considered. The displacements of the Reissner–Mindlin theory are enriched by a an additional part. These so-called fluctuation displacements include warping displacements and thickness changes. They lead to additional terms for the material deformation gradient and the Green–Lagrangian strain tensor. Within a nonlinear multi-field variational

In this work, a new thermo-mechanical formulation for the conventional and localizing gradient damage method is proposed. The proposed formulation is based on the generalized micromorphic theory, which accounts for the underlying fracture processes at the micro-level. The thermal and mechanical effects on the fracture response are incorporated in the formulation through three primary variables. These

The work is devoted to a coupling method for the finite element method (FEM) and the distance potential discrete element method. In this work, a well-defined distance potential function is developed. Meanwhile, a holonomic and precise algorithm for contact interaction is established, accounting for the influence of the tangential contact force. In addition, the measurement of deformation behaviors

We consider the numerical simulation of the left ventricle of the human heart by a hyperelastic fiber reinforced transversely isotropic model. This is an important model problem for the understanding of the mechanical properties of the human heart but its calculation is very time consuming because the lack of fast, scalable method that is also robust with respect to the model parameters. In this paper

We show that the logarithmic (Hencky) strain and its derivatives can be approximated, in a straightforward manner and with a high accuracy, using Padé approximants of the tensor (matrix) logarithm. Accuracy and computational efficiency of the Padé approximants are favourably compared to an alternative approximation method employing the truncated Taylor series. As an application, Hencky-type hyperelasticity

The article “The computational framework for continuum-kinematics-inspired peridynamics”.

Non-intrusive methods have been used since two decades to derive reduced-order models for geometrically nonlinear structures, with a particular emphasis on the so-called STiffness Evaluation Procedure (STEP), relying on the static application of prescribed displacements in a finite-element context. We show that a particularly slow convergence of the modal expansion is observed when applying the method

Selective laser melting (SLM) has gained large interest due to advanced manufacturing possibilities. However, the growing potential also necessitates reliable predictions of structures in particular regarding their long-term behaviour. The constitutive and structural response is thereby challenging to reproduce, due to the complex material behaviour. This motivates the aims of this contribution: To

We present computational flow analysis of a vertical-axis wind turbine (VAWT) that has been proposed to also serve as a tsunami shelter. In addition to the three-blade rotor, the turbine has four support columns at the periphery. The columns support the turbine rotor and the shelter. Computational challenges encountered in flow analysis of wind turbines in general include accurate representation of

X-ray tomography techniques give researchers the full access to material inner structures. With such ample information, employing numerical simulation on real material images becomes more and more common. Material behavior, especially, of heterogeneous materials, e.g. polycrystalline, composite, can be observed at the microscopic scale with the computed tomography (CT) techniques. In this work, an

A key strategy to prevent a local outbreak during the COVID-19 pandemic is to restrict incoming travel. Once a region has successfully contained the disease, it becomes critical to decide when and how to reopen the borders. Here we explore the impact of border reopening for the example of Newfoundland and Labrador, a Canadian province that has enjoyed no new cases since late April, 2020. We combine

The COVID-19 pandemic has captivated scientific activity since its early days. Particular attention has been dedicated to the identification of underlying dynamics and prediction of future trend. In this work, a switching Kalman filter formalism is applied on dynamics learning and forecasting of the daily new cases of COVID-19. The main feature of this dynamical system is its ability to switch between

A multiscale clustering-based topology optimization for thermo-elastic lattice structures is studied based on the extended multiscale finite element method (EMsFEM). The strain energy of thermo-elastic lattice structures is chosen as the objective function. The microstructural configuration and the macrostructural distribution of the thermo-elastic lattice material are designed through topology optimization

In this paper, we derive an analytical expression for reproducing kernel particle method (RKPM) shape functions. Based on this, we propose a necessary and sufficient stability condition for general RKPM in arbitrary function space, and illustrate with degenerate cases. By selecting proper basis vectors and the support of the kernel functions, we demonstrate that the RKPM framework allows generating

The article “Non-localised contact between beams with circular and elliptical cross-sections”, written by “Marco Magliulo, Jakub Lengiewicz, Andreas Zilian and Lars A. A. Beex”, was originally published Online First without open access. After publication in volume 65, issue 5, page 1247–1266 the authors decided to opt for Open Choice and to make the article an open access publication. Therefore, the

Marine gas hydrate is an important energy source while its extraction may induce environmental problems such as subsea landslide, which is usually challengeable for numerical simulation due to the marine environment with high pressure and the existence of gas hydrate. Smoothed particle hydrodynamics (SPH) is a Lagrangian particle method which is attractive in modeling problems with large deformations

The virtual element method is a lively field of research, in which considerable progress has been made during the last decade and applied to many problems in physics and engineering. The method allows ansatz function of arbitrary polynomial degree. However, one of the prerequisite of the formulation is that the element edges have to be straight. In the literature there are several new formulations

The outbreak of COVID-19 in 2020 has led to a surge in interest in the research of the mathematical modeling of epidemics. Many of the introduced models are so-called compartmental models, in which the total quantities characterizing a certain system may be decomposed into two (or more) species that are distributed into two (or more) homogeneous units called compartments. We propose herein a formulation

We extend the classical SIR model of infectious disease spread to account for time dependence in the parameters, which also include diffusivities. The temporal dependence accounts for the changing characteristics of testing, quarantine and treatment protocols, while diffusivity incorporates a mobile population. This model has been applied to data on the evolution of the COVID-19 pandemic in the US

A multiscale crack band model is proposed to alleviate spurious mesh size dependence in the solution for a strain softening constitutive model. A continuum damage mechanics-based material model for the constituent phases of a carbon fiber reinforced polymer (CFRP) based composite material is applied to the Eigenstrain based Reduced-order Homogenization framework for multiscale modeling. This paper

Phase-field methods for fracture have been integrated with plasticity for better describing constitutive behaviours. In most of the previous phase-field models, however, the length-scale parameter must be interpreted as a material property in order to match the material strength in experiments. This study presents a phase-field model for fracture coupled with plasticity for quasi-brittle materials

A summary is given of the mechanical characteristics of virus contaminants and the transmission via droplets and aerosols. The ordinary and partial differential equations describing the physics of these processes with high fidelity are presented, as well as appropriate numerical schemes to solve them. Several examples taken from recent evaluations of the built environment are shown, as well as the

The scaled boundary finite element method (SBFEM) has recently been employed as an efficient tool to model three-dimensional structures, in particular when the geometry is provided as a voxel-based image. To this end, an octree decomposition of the computational domain is deployed, and each cubic cell is treated as an SBFE subdomain. The surfaces of each subdomain are discretized in the finite element

The COVID-19 pandemic has led to an unprecedented world-wide effort to gather data, model, and understand the viral spread. Entire societies and economies are desperate to recover and get back to normality. However, to this end accurate models are of essence that capture both the viral spread and the courses of disease in space and time at reasonable resolution. Here, we combine a spatially resolved

The presence of uncertainty in material properties and geometry of a structure is ubiquitous. The design of robust engineering structures, therefore, needs to incorporate uncertainty in the optimization process. Stochastic gradient descent (SGD) method can alleviate the cost of optimization under uncertainty, which includes statistical moments of quantities of interest in the objective and constraints

The increase in readily available computational power raises the possibility that direct agent-based modeling can play a key role in the analysis of epidemiological population dynamics. Specifically, the objective of this work is to develop a robust agent-based computational framework to investigate the emergent structure of Susceptible-Infected-Removed/Recovered (SIR)-type populations and variants

We consider a mixture-theoretic continuum model of the spread of COVID-19 in Texas. The model consists of multiple coupled partial differential reaction–diffusion equations governing the evolution of susceptible, exposed, infectious, recovered, and deceased fractions of the total population in a given region. We consider the problem of model calibration, validation, and prediction following a Bayesian

Peridynamics (PD) is a non-local continuum formulation. The original version of PD was restricted to bond-based interactions. Bond-based PD is geometrically exact and its kinematics are similar to classical continuum mechanics (CCM). However, it cannot capture the Poisson effect correctly. This shortcoming was addressed via state-based PD, but the kinematics are not accurately preserved. Continuum