The storage of renewable hydrogen in salt caverns requires fast injection and production rates to cope with the imbalance between energy production and consumption. This raises concerns about the mechanical stability of salt caverns under such operational conditions. The use of appropriate constitutive models for salt mechanics is an important step in investigating this issue, therefore many constitutive models with several parameters have been presented in the literature. However, a robust calibration strategy to reliably determine which model and parameter set represents the given rock, based on stress–strain data sets, remains an unsolved challenge. In this paper, for the first time in the community, we present a multi-step strategy to determine a single parameter set based on many deformation data sets for salt rocks. Towards this end, we first develop a comprehensive constitutive model able to capture all relevant nonlinear deformation physics of transient, reverse, and steady-state creep. The determination of the single set of representative material parameters is then achieved by framing the calibration process as an optimization problem, for which the global Particle Swarm Optimization algorithm is employed. To allow for dynamic data integration, a multi-step calibration strategy is developed for a situation where experiments are included one at a time, as they become available. Additionally, due to the existing mild heterogeneity in the experimental rock samples, our optimization strategy is made flexible to allow for this slight heterogeneity. The devised optimization strategy, based on the multi-physics comprehensive constitutive modeling framework, results in a single set of representative material properties of all the deformation data sets. As a rigorous mathematical analysis for the presented method and the lack of relevant experimental data sets, we consider a wide range of synthetic experimental data sets, inspired by the existing sparse relevant data in the literature. The results of our performance analyses show that the proposed calibration strategy is robust. Moreover, the results become increasingly more accurate as more data sets become available.

Lightweight materials used for impact mitigation must be able to resist impact and absorb the maximum amount of energy from the impactor. Auxetic materials have the potential to achieve high resistance by drawing material into the impact zone and providing higher indentation and shear resistance. However, these materials must be artificially designed, and the large deformation dynamic effects of the created structures must be taken into consideration when deciding on a protection concept. Despite their promise, little attention has been given to understanding the working mechanisms of high-rate and finite deformation effects of architected auxetic lattice structures. This study compares the static and dynamic elastic properties of different auxetic structures with a honeycomb structure, a typical non-auxetic lattice, at equivalent mass and stiffness levels. In this study, we limit the investigation to elastic material behavior and do not consider contact between the beams of the lattices. It is demonstrated that the equivalent static and dynamic properties of individual lattices at an undeformed state are insufficient to explain the variations observed in impact situations. In particular, the initial Poisson's ratio does not determine the ability of a structure to resist impact. To gain a thorough comprehension of the overall behavior of these structures during localized, high rate compression, the evolution of the elastic tangent properties under compression and shear deformation was monitored, leading to a more profound understanding. Observations made in one configuration of stiffness and mass are replicated and analyzed in related configurations.

A yield function in the stress resultant space of geometrically exact beams based on the elastoplastic cross-sectional warping problem has been proposed by Herrnböck et al. (Comput Mech, 67(3):723–742, 2021). This plasticity framework has been extended with a hardening tensor to model the kinematic hardening effects in Herrnböck et al. (Comput Mech, 71(1):1–24, 2022). While this framework provides scaling for the yield surface in ideal plasticity, scaling in hardening plasticity has not yet been explored. This paper focuses on the numeric modelling of hardening beams and beam assemblies at different geometric scales. Discretization effects from the introduction of plasticity into the geometrically exact beam model are demonstrated. Furthermore, the effects of scaling are explored, and a method to mitigate undesirable effects in order to achieve a size-agnostic formulation is proposed. Consistent geometric scaling is demonstrated for two alternative scaling approaches of the yield function.

In quasi-brittle materials, such as reinforced concrete, localisation of initially diffuse cracking evolving in localised cracking patterns consists of a numerically challenging task. With conventional iterative methods, convergence of the numerical solution scheme to model crack localisation is often difficult to obtain. On the other hand, conventional total approaches, such as the Sequentially Linear Approach, although robust, fail to approximate properly the underlying material law. In the present work, a new model is introduced, designated the Total Iterative Approach, in which the internal damage variables are updated iteratively. It is found that this approach is robust, allows for the correct approximation of the material law and is a powerful tool for the analysis of softening behaviour. Some examples are presented to illustrate the performance of the model.

In fractured geological formations, as a result of the in-situ stress changes, fractures can propagate or slide. This phenomenon can be efficiently modeled by the extended finite element method (XFEM) when there are only a few fractures present. However, geological reservoirs contain many fractures which can also cross and are densely populated. Therefore, the classical XFEM is too expensive to be applied for the simulation of propagating fractures in geological formations. To reduce the costs, typically, homogenization or upscaling is used. However, they result in inaccurate solutions, since no separation of scales exists in this process. To resolve this challenge, in this work, a multiscale XFEM (MS-XFEM) for propagating fractures is developed and presented. In each time step, given the current geometries of the fractures, local XFEM-based basis functions are constructed or adaptively updated. The adaptive update takes place in certain regions where fracture geometries are changed due to propagation. Using these basis functions, a very efficient FEM-based coarse-scale system is developed since it has no extra degrees of freedom (DOFs). Once the coarse-scale system is solved, its solution is prolonged to the fine-scale original resolution using the basis functions. This approximate fine-scale solution is then used to estimate the group of growing fractures tips and their growing angles. This allows for exploiting the locality of the propagation process fully while solving a global system. To control the error, an iterative procedure is also developed. Proof-of-concept test cases are presented to study the developed MS-XFEM algorithm. It is shown that MS-XFEM results are capable of predicting the propagating paths for complicated fracture patterns. As such, MS-XFEM casts a promising method for field-scale applications.

In this work, a recently proposed high-cycle fatigue cohesive zone model, which covers crack initiation and propagation with limited input parameters, is embedded in a robust and efficient numerical framework for simulating progressive failure in composite laminates under fatigue loading. The fatigue cohesive zone model is enhanced with an implicit time integration scheme of the fatigue damage variable which allows for larger cycle increments and more efficient analyses. The method is combined with an adaptive strategy for determining the cycle increment based on global convergence rates. Moreover, a consistent material tangent stiffness matrix has been derived by fully linearizing the underlying mixed-mode quasi-static model and the fatigue damage update. The enhanced fatigue cohesive zone model is used to describe matrix cracking and delamination in laminates. In order to allow for matrix cracks to initiate at arbitrary locations and to avoid complex and costly mesh generation, the phantom node version of the eXtended finite element method (XFEM) is employed. For the insertion of new crack segments, an XFEM fatigue crack insertion criterion is presented, which is consistent with the fatigue cohesive zone formulation. It is shown with numerical examples that the improved fatigue damage update enhances the accuracy, efficiency and robustness of the numerical simulations significantly. The numerical framework is applied to the simulation of progressive fatigue failure in an open-hole [±45]-laminate. It is demonstrated that the numerical model is capable of accurately and efficiently simulating the complete failure process from distributed damage to localized failure.

Corrosion in reinforced concrete is an important feature which can lead to increased deformation and cracking, as well as to premature failure. In the present work, macro-mechanical modelling of corrosion is performed, namely the degradation of bond–slip between concrete and steel. A mixed-mode damage model is adopted, in which the interaction between the bond–slip law and the stress acting in the neighbourhood of the concrete–steel bar interface is taken into account. Bond–slip degradation is modelled using an evolutionary bond–slip relationship, which depends on the level of corrosion. Different relevant loading cases are studied. Special attention is given to the evolution of corrosion in time, under constant load. This is done by adopting a Total Iterative Approach, in which the structure is reevaluated each time step, upon damage increase due to corrosion. Pullout tests are presented to illustrate the performance of the model. Bending tests are also performed to evaluate the influence of corrosion at structural level.

Accurate assessment of the overall mechanical behavior of masonry, composed of bricks and mortar joints, remains challenging due to its inhomogeneous and orthotropic nature. In this study, the feasibility of various mean-field homogenization schemes for the three-dimensional orthotropic elastic properties of masonry is comprehensively investigated. Three kinds of masonry patterns are considered, including the stack bonded pattern, the running bonded pattern and the double-leaf Flemish bonded pattern that has received limited attention so far. Special attention is paid to the homogenization schemes which have not been applied to the masonry case, such as Lielens’ interpolative double inclusions (D-I) and the interaction direct derivative (IDD) schemes. After a comparison between the well-known mean-field homogenization schemes, an improved micro-mechanical model is proposed by combining the advantages of the IDD and D-I models. The validation of the proposed model is conducted through a comparison against experimental data from literature and numerical results obtained via finite element analyses (FEA). The results show that the proposed model can accurately evaluate the orthotropic elastic properties of the three masonry typologies for a wide range of stiffness ratios between brick and mortar, ranging from 1 to 1000. The proposed model also shows better performance than the classical schemes especially when the stiffness ratios between brick and mortar are higher than 10, which is of major importance for the application of mean-field homogenization based multiscale methods to the nonlinear analysis of masonry. Furthermore, the presented homogenization method can be of interest for other anisotropic materials, e.g., laminate materials.

Corrosion is a phenomenon observed in structural components in corrosive environments such as pipelines, bridges, aircrafts, turbines, etc. The computational model of corrosion should enjoy two features: a) accurately considering the electrochemistry of corrosion and b) properly dealing with the moving interface between solid and electrolyte. There are several approaches to model corrosion such as using FEM with mesh refinement algorithms, combining FEM and level-set method, employing finite volume methods, adopting peridynamic formulation, and utilizing phase field models. Because of its accuracy, lower computational cost, and robust dealing with multiple pit merging, the model which combines FEM with level-set method is selected to be more extensively assessed in this paper. Part I focuses on demonstrating the model's capabilities of simulating pitting corrosion through a set of numerical examples which include numerical solution verification, experimental validation, and uncertainty quantification of model parameters and properties.

Structural components in corrosive environments such as pipelines, bridges, aircrafts, and turbines are imposed to stress corrosion. A stress corrosion model for pit growth should a) accurately consider the electrochemistry of the corrosion process, b) properly deal with the moving interface between solid and electrolyte, and c) effectively incorporates the synergism between corrosion and mechanical field at the interface. In Part II, the influence of mechanical loading is added to the approach described in Part I. Part II investigates the model's capabilities of simulating stress corrosion via a set of numerical examples of corrosion pitting which include experimental validation and uncertainty quantification of model parameters and properties.

Lattice models have been used to simulate mass transport to predict durability of cementitious materials. In particular, the use of dual lattice meshes allows for the coupling of fracture and transport processes, which commonly occur at the same time in these materials. Literature has shown good agreement between simulations and experimental results. Nevertheless, work regarding relevant computational aspects of the numerical model are scarce. In this study, a Voronoi-discretized lattice model is used to simulate unsaturated moisture transport in cement-base materials through the Richards equation. First, investigations regarding the choice of elemental volume approximation, time-stepping procedure and quadrature are evaluated. After validation of the approximations, simulated moisture transport in sound concrete was compared to experiments and mesh and time step sensitivity were discussed. A new approach to model capillary absorption of water in cracked concrete was also proposed and its advantages with respect to existing approaches are discussed by comparing to experimental measurements. The results confirm that the model can accurately predict the transport processes for the earlier stage of capillary absorption. Furthermore, moisture ingress in cracked concrete is simulated for different crack configurations and the use of different approaches is suggested accordingly. Finally, guidelines regarding the approximations used for optimization of the computations are presented.

Altering the state of the stress of the subsurface reservoirs can lead fractures to slip and extend their lengths (i.e., to propagate). This process can even be engineered, in many applications, e.g., enhanced geothermal systems. As such, accurate and efficient simulation of the mechanical deformation of the subsurface geological reservoirs, allowing for fracture propagation, is at the core of many geoscientific operational designs. Subsurface reservoirs entail many fractures at multiple scales. Implementation of 3D complex grids on these complex fractured systems, for mechanical deformation analyses, is extremely challenging. An alternative approach can be developed by using extended finite element methods (XFEM). XFEM allows for capturing the fractures effects on a conveniently-generated structured matrix mesh. The cracks are introduced by extra degrees of freedom (DOFs) on the nodes of the matrix rock mesh. For geoscientific applications, however, XFEM results in too many DOFs which are beyond the scope of simulators. Additionally, for propagating fractures, these DOFs need to be updated in response to the dynamic extension of the fractures in the domain. The propagation process not only adds to the sensitivity of the outputs to the accuracy of the estimated stress field, but also increases the size of the linear systems. In addition to these, matrix rocks are often highly heterogeneous, at high resolutions. In this work, we present a novel multiscale procedure for propagating fractures in heterogeneous geological reservoirs. For the first time in the community, we present the highly fractured systems at coarser resolutions via XFEM-based basis functions, which also account for the propagating effects. Fractures are allowed to extend their scale and the enriched basis functions are locally updated. Using these bases, the coarse scale system is obtained in which no extra DOFs due to fractures exist. This significantly reduces the computational complexity. As a significant step forward compared with our recently-published journal paper [Xu, Hajibeygi, Sluys, Journal of Computational Physics, 2021], in this conference contribution we allow the fractures to propagate. Specially, we introduce a local-global-based approach, in which fracture propagation is treated only at local stage; while the stress and deformation are modelled at global scale. In the search of convenient implementation, the procedure is presented algebraically. Through several test cases, we demonstrate the applicability of the method for complex fractured media. Specially we demonstrate that propagation can be modeled at local scale, while accurate stress and deformation fields are obtained at global scale.

Adobe is a traditional masonry made of sundried earthen bricks and mud mortar. Despite a millennial history of buildings of architectural value, adobe still connotes a so called ‘not engineered’ construction type. Namely, the material and structural properties of adobe are still not entirely addressed, resulting in an equally uncertain normative framework for adobe buildings design. However, over the last ten years, a large research program has been conducted in the Netherlands to qualify the material and structural properties of this sustainable building technology. In this paper, a critical analysis of the current normative body for the material characterization of adobe is addressed. Guidelines, prescriptions and requirements related to test methods, materials selection and properties contained in the available building codes for adobe around the world are assessed. A critical normative review is performed using the most recent literature produced on adobe, with particular regards to the results of experimental tests and numerical simulations performed by the authors. On the basis of these findings, some issues have been identified in relation to the knowledge currently condensed in the norms for adobe. A series of programmatic guidelines is aimed at orienting future research on adobe as well as fostering the process of updating its current normative body.

Calcium phosphate cements (CPCs) have been widely used during the past decades as biocompatible bone substitution in maxillofacial, oral and orthopedic surgery. CPCs are injectable and are chemically resemblant to the mineral phase of native bone. Nevertheless, their low fracture toughness and high brittleness reduce their clinical applicability to weakly loaded bones. Reinforcement of CPC matrix with polymeric fibers can overcome these mechanical drawbacks and significantly enhance their toughness and strength. Such fiber-reinforced calcium phosphate cements (FRCPCs) have the potential to act as advanced bone substitute in load-bearing anatomical sites. This work achieves integrated experimental and numerical characterization of the mechanical properties of FRCPCs under bending and tensile loading. To this end, a 3-D numerical gradient enhanced damage model combined with a dimensionally-reduced fiber model are employed to develop a computational model for material characterization and to simulate the failure process of fiber-reinforced CPC matrix based on experimental data. In addition, an advanced interfacial constitutive law, derived from micromechanical pull-out tests, is used to represent the interaction between the polymeric fiber and CPC matrix. The presented computational model is successfully validated with the experimental results and offers a firm basis for further investigations on the development of numerical and experimental analysis of fiber-reinforced bone cements.

This study presents calibration and validation of a cohesive extended finite element model for fatigue crack propagation in ductile materials. The approach relies on a separation between plasticity around the crack tip and fatigue crack growth at the crack tip such that the influence of plasticity on fatigue driving forces is predicted. This implies that characterization of crack growth requires effective Paris parameters. It is shown that the calibrated model can capture fatigue crack growth behaviour in ductile materials for in-phase and out-of-phase biaxial fatigue loading as well as in-phase biaxial loading with an overload.

Deformable fractured porous media appear in many geoscience applications. While the extended finite element method (XFEM) has been successfully developed within the computational mechanics community for accurate modeling of deformation, its application in geoscientific applications is not straightforward. This is mainly due to the fact that subsurface formations are heterogeneous and span large length scales with many fractures at different scales. To resolve this limitation, in this work, we propose a novel multiscale formulation for XFEM, based on locally computed enriched basis functions. The local multiscale basis functions capture heterogeneity of th e porous rock properties, and discontinuities introduced by the fractures. In order to preserve accuracy of these basis functions, reduced-dimensional boundary conditions are set as localization condition. Using these multiscale bases, a multiscale coarse-scale system is then governed algebraically and solved. The coarse scale system entails no enrichment due to the fractures. Such formulation allows for significant computational cost reduction, at the same time, it preserves the accuracy of the discrete displacement vector space. The coarse-scale solution is finally interpolated back to the fine scale system, using the same multiscale basis functions. The proposed multiscale XFEM (MS-XFEM) is also integrated within a two-stage algebraic iterative solver, through which error reduction to any desired level can be achieved. Several proof-of-concept numerical tests are presented to assess the performance of the developed method. It is shown that the MS-XFEM is accurate, when compared with the fine-scale reference XFEM solutions. At the same time, it is significantly more efficient than the XFEM on fine-scale resolution, as it significantly reduces the size of the linear systems. As such, it develops a promising scalable XFEM method for large-scale heavily fractured porous media.

A numerical model for corrosion pit propagation under mechanical loading is presented. The level set method is used for corrosion front tracking and also enables the domain to be split into a solid and a pit domain. In the pit the diffusion of atoms originating from the dissolution process occurring at the pit front is simulated. The model is capable of automatically capturing lacy cover formation due to the inclusion of activation control, diffusion control, and passivation. In the solid static equilibrium is solved to obtain strains and stresses. A parameter, dependent on the signs of the plastic strain increment and the back stress, is introduced to define the influence of plasticity on the corrosion rate. The model is used to study pit growth under electrochemical and mechanical loading. Under activation control combined with an elastic material response, pits propagate faster under constant loading than under cyclic loading. When plastic deformation occurs, cyclic loading can significantly increase the pit growth rate. Increasing the cyclic load frequency results in faster propagation due to kinematic hardening. Under diffusion control, mechanical loading does not influence the pit growth rate, given that the salt layer leading to diffusion control remains intact.