Showing 1–33 of 33 results for author: Auricchio, F

Numerical investigation of fluid-structure interaction in a pilot-operated microfluidic valve

Authors: Ahmed Aissa Berraies, E. Harald van Brummelen, Ferdinando Auricchio

Abstract: The present paper is concerned with numerical investigation of the performance of a pilot-operated control valve based on shape memory alloy actuation control. The valve under investigation can be integrated into miniaturized hydraulic systems and is developed to perform precise dispensing, mixing, or dosing tasks while being able to withstand relatively high pressure differences. The study evalua… ▽ More The present paper is concerned with numerical investigation of the performance of a pilot-operated control valve based on shape memory alloy actuation control. The valve under investigation can be integrated into miniaturized hydraulic systems and is developed to perform precise dispensing, mixing, or dosing tasks while being able to withstand relatively high pressure differences. The study evaluates the valve's response under the current ON/OFF and the desired proportional control regimes using numerical methods for fluid-structure interaction. The computational model replicates the operation of the valve, which requires an understanding of the complex interactions between the fluid flow with the pressurized valve and the contact with the valve seat during the opening and closing processes. In addition, the model leverages advanced numerical techniques to overcome several complexities arising mainly from the geometrical, material, and contact nonlinearities, and to mitigate the shortcomings of the partitioned fluid-structure interaction approach. Several simulations are conducted to examine the valve's structural and flow behavior under varying pressure conditions. Results indicate that the valve is adequate for ON/OFF actuation control but is susceptible to flow-induced vibrations during the proportional control regime that occurs due to the sharp pressure drop in the valve-seat gap and the ensuing Venturi effect, which counteract the opening of the main valve. The fluid-structure-interaction simulations provide insight into the mechanism underlying the flow-induced vibrations, which can serve to improve the design and enhance the performance of the valve in microfluidic applications. △ Less

Submitted 28 April, 2024; originally announced April 2024.

Comments: 41 pages, 18 figures, submitted to Journal of Fluids and Structures, published in SSRN

Analysis of a combined Filtered/phase-field approach to topology optimization in elasticit

Authors: Ferdinando Auricchio, Michele Marino, Idriss Mazari, Ulisse Stefanelli

Abstract: We advance a combined filtered/phase-field approach to topology optimization in the setting of linearized elasticity. Existence of minimizers is proved and rigorous parameter asymptotics are discussed by means of variational convergence techniques. Moreover, we investigate an abstract space discretization in the spirit of conformal finite elements. Eventually, stationarity is equivalently reformul… ▽ More We advance a combined filtered/phase-field approach to topology optimization in the setting of linearized elasticity. Existence of minimizers is proved and rigorous parameter asymptotics are discussed by means of variational convergence techniques. Moreover, we investigate an abstract space discretization in the spirit of conformal finite elements. Eventually, stationarity is equivalently reformulated in terms of a Lagrangian. △ Less

Submitted 23 January, 2024; originally announced January 2024.

Data-informed uncertainty quantification for laser-based powder bed fusion additive manufacturing

Authors: Mihaela Chiappetta, Chiara Piazzola, Lorenzo Tamellini, Alessandro Reali, Ferdinando Auricchio, Massimo Carraturo

Abstract: We present an efficient approach to quantify the uncertainties associated with the numerical simulations of the laser-based powder bed fusion of metals processes. Our study focuses on a thermomechanical model of an Inconel 625 cantilever beam, based on the AMBench2018-01 benchmark proposed by the National Institute of Standards and Technology (NIST). The proposed approach consists of a forward unc… ▽ More We present an efficient approach to quantify the uncertainties associated with the numerical simulations of the laser-based powder bed fusion of metals processes. Our study focuses on a thermomechanical model of an Inconel 625 cantilever beam, based on the AMBench2018-01 benchmark proposed by the National Institute of Standards and Technology (NIST). The proposed approach consists of a forward uncertainty quantification analysis of the residual strains of the cantilever beam given the uncertainty in some of the parameters of the numerical simulation, namely the powder convection coefficient and the activation temperature. The uncertainty on such parameters is modelled by a data-informed probability density function obtained by a Bayesian inversion procedure, based on the displacement experimental data provided by NIST. To overcome the computational challenges of both the Bayesian inversion and the forward uncertainty quantification analysis we employ a multi-fidelity surrogate modelling technique, specifically the multi-index stochastic collocation method. The proposed approach allows us to achieve a 33\% reduction in the uncertainties on the prediction of residual strains compared with what we would get basing the forward UQ analysis on a-priori ranges for the uncertain parameters, and in particular the mode of the probability density function of such quantities (i.e., its ``most likely value'', roughly speaking) results to be in good agreement with the experimental data provided by NIST, even though only displacement data were used for the Bayesian inversion procedure. △ Less

Submitted 20 March, 2024; v1 submitted 7 November, 2023; originally announced November 2023.

Comments: 28 pages, 14 figures

On the accuracy of interpolation based on single-layer artificial neural networks with a focus on defeating the Runge phenomenon

Authors: Ferdinando Auricchio, Maria Roberta Belardo, Gianluca Fabiani, Francesco Calabrò, Ariel F. Pascaner

Abstract: In the present paper, we consider one-hidden layer ANNs with a feedforward architecture, also referred to as shallow or two-layer networks, so that the structure is determined by the number and types of neurons. The determination of the parameters that define the function, called training, is done via the resolution of the approximation problem, so by imposing the interpolation through a set of sp… ▽ More In the present paper, we consider one-hidden layer ANNs with a feedforward architecture, also referred to as shallow or two-layer networks, so that the structure is determined by the number and types of neurons. The determination of the parameters that define the function, called training, is done via the resolution of the approximation problem, so by imposing the interpolation through a set of specific nodes. We present the case where the parameters are trained using a procedure that is referred to as Extreme Learning Machine (ELM) that leads to a linear interpolation problem. In such hypotheses, the existence of an ANN interpolating function is guaranteed. The focus is then on the accuracy of the interpolation outside of the given sampling interpolation nodes when they are the equispaced, the Chebychev, and the randomly selected ones. The study is motivated by the well-known bell-shaped Runge example, which makes it clear that the construction of a global interpolating polynomial is accurate only if trained on suitably chosen nodes, ad example the Chebychev ones. In order to evaluate the behavior when growing the number of interpolation nodes, we raise the number of neurons in our network and compare it with the interpolating polynomial. We test using Runge's function and other well-known examples with different regularities. As expected, the accuracy of the approximation with a global polynomial increases only if the Chebychev nodes are considered. Instead, the error for the ANN interpolating function always decays and in most cases we observe that the convergence follows what is observed in the polynomial case on Chebychev nodes, despite the set of nodes used for training. △ Less

Submitted 7 May, 2024; v1 submitted 21 August, 2023; originally announced August 2023.

MSC Class: 65D05

Impact of interaction forces in first order many-agent systems for swarm manufacturing

Authors: Ferdinando Auricchio, Massimo Carraturo, Giuseppe Toscani, Mattia Zanella

Abstract: We study the large time behavior of a system of interacting agents modeling the relaxation of a large swarm of robots, whose task is to uniformly cover a portion of the domain by communicating with each other in terms of their distance. To this end, we generalize a related result for a Fokker-Planck-type model with a nonlocal discontinuous drift and constant diffusion, recently introduced by three… ▽ More We study the large time behavior of a system of interacting agents modeling the relaxation of a large swarm of robots, whose task is to uniformly cover a portion of the domain by communicating with each other in terms of their distance. To this end, we generalize a related result for a Fokker-Planck-type model with a nonlocal discontinuous drift and constant diffusion, recently introduced by three of the authors, of which the steady distribution is explicitly computable. For this new nonlocal Fokker-Planck equation, existence, uniqueness and positivity of a global solution are proven, together with precise equilibration rates of the solution towards its quasi-stationary distribution. Numerical experiments are designed to verify the theoretical findings and explore possible extensions to more complex scenarios. △ Less

Submitted 9 July, 2023; originally announced July 2023.

Trends to equilibrium for a nonlocal Fokker-Planck equation

Authors: Ferdinando Auricchio, Giuseppe Toscani, Mattia Zanella

Abstract: We obtain equilibration rates for a one-dimensional nonlocal Fokker-Planck equation with time-dependent diffusion coefficient and drift, modeling the relaxation of a large swarm of robots, feeling each other in terms of their distance, towards the steady profile characterized by uniform spreading over a finite interval of the line. The result follows by combining entropy methods for quantifying th… ▽ More We obtain equilibration rates for a one-dimensional nonlocal Fokker-Planck equation with time-dependent diffusion coefficient and drift, modeling the relaxation of a large swarm of robots, feeling each other in terms of their distance, towards the steady profile characterized by uniform spreading over a finite interval of the line. The result follows by combining entropy methods for quantifying the decay of the solution towards its quasi-stationary distribution, with the properties of the quasi-stationary profile. △ Less

Submitted 5 June, 2023; v1 submitted 10 March, 2023; originally announced March 2023.

Sparse-grids uncertainty quantification of part-scale additive manufacturing processes

Authors: Mihaela Chiappetta, Chiara Piazzola, Massimo Carraturo, Lorenzo Tamellini, Alessandro Reali, Ferdinando Auricchio

Abstract: The present paper aims at applying uncertainty quantification methodologies to process simulations of powder bed fusion of metal. In particular, for a part-scale thermomechanical model of an Inconel 625 super-alloy beam, we study the uncertainties of three process parameters, namely the activation temperature, the powder convection coefficient and the gas convection coefficient. First, we perform… ▽ More The present paper aims at applying uncertainty quantification methodologies to process simulations of powder bed fusion of metal. In particular, for a part-scale thermomechanical model of an Inconel 625 super-alloy beam, we study the uncertainties of three process parameters, namely the activation temperature, the powder convection coefficient and the gas convection coefficient. First, we perform a variance-based global sensitivity analysis to study how each uncertain parameter contributes to the variability of the beam displacements. The results allow us to conclude that the gas convection coefficient has little impact and can therefore be fixed to a constant value for subsequent studies. Then, we conduct an inverse uncertainty quantification analysis, based on a Bayesian approach on synthetic displacements data, to quantify the uncertainties of the two remaining parameters, namely the activation temperature and the powder convection coefficient. Finally, we use the results of the inverse uncertainty quantification analysis to perform a data-informed forward uncertainty quantification analysis of the residual strains. Crucially, we make use of surrogate models based on sparse grids to keep to a minimum the computational burden of every step of the uncertainty quantification analysis. The proposed uncertainty quantification workflow allows us to substantially ease the typical trial-and-error approach used to calibrate power bed fusion part-scale models, and to greatly reduce uncertainties on the numerical prediction of the residual strains. In particular, we demonstrate the possibility of using displacement measurements to obtain a data-informed probability density function of the residual strains, a quantity much more complex to measure than displacements. △ Less

Submitted 19 April, 2023; v1 submitted 13 October, 2022; originally announced October 2022.

Fokker-Planck modeling of many-agent systems in swarm manufacturing: asymptotic analysis and numerical results

Authors: Ferdinando Auricchio, Giuseppe Toscani, Mattia Zanella

Abstract: In this paper we study a novel Fokker-Planck-type model that is designed to mimic manufacturing processes through the dynamics characterizing a large set of agents. In particular, we describe a many-agent system interacting with a target domain in such a way that each agent/particle is attracted by the center of mass of the target domain with the aim to uniformly cover this zone. To this end, we f… ▽ More In this paper we study a novel Fokker-Planck-type model that is designed to mimic manufacturing processes through the dynamics characterizing a large set of agents. In particular, we describe a many-agent system interacting with a target domain in such a way that each agent/particle is attracted by the center of mass of the target domain with the aim to uniformly cover this zone. To this end, we first introduce a mean-field model with discontinuous flux whose large time behavior is such that the steady state is globally continuous and uniform over a connected portion of the domain. We prove that a diffusion coefficient that guarantees that a given portion of mass enters in the target domain exists and that it is unique. Furthermore, convergence to equilibrium in 1D is provided through a reformulation of the initial problem involving a nonconstant diffusion function. The extension to 2D is explored numerically by means of recently introduced structure preserving methods for Fokker-Planck equations. △ Less

Submitted 7 December, 2022; v1 submitted 20 May, 2022; originally announced May 2022.

Well-posedness for a diffusion-reaction compartmental model simulating the spread of COVID-19

Authors: Ferdinando Auricchio, Pierluigi Colli, Gianni Gilardi, Alessandro Reali, Elisabetta Rocca

Abstract: This paper is concerned with the well-posedness of a diffusion-reaction system for a Susceptible-Exposed-Infected-Recovered (SEIR) mathematical model. This model is written in terms of four nonlinear partial differential equations with nonlinear diffusions, depending on the total amount of the SEIR populations. The model aims at describing the spatio-temporal spread of the COVID-19 pandemic and is… ▽ More This paper is concerned with the well-posedness of a diffusion-reaction system for a Susceptible-Exposed-Infected-Recovered (SEIR) mathematical model. This model is written in terms of four nonlinear partial differential equations with nonlinear diffusions, depending on the total amount of the SEIR populations. The model aims at describing the spatio-temporal spread of the COVID-19 pandemic and is a variation of the one recently introduced, discussed and tested in [A. Viguerie et al, Diffusion-reaction compartmental models formulated in a continuum mechanics framework: application to COVID-19, mathematical analysis, and numerical study, Comput. Mech. 66 (2020) 1131-1152]. Here, we deal with the mathematical analysis of the resulting Cauchy-Neumann problem: the existence of solutions is proved in a rather general setting and a suitable time discretization procedure is employed. It is worth mentioning that the uniform boundedness of the discrete solution is shown by carefully exploiting the structure of the system. Uniform estimates and passage to the limit with respect to the time step allow to complete the existence proof. Then, two uniqueness theorems are offered, one in the case of a constant diffusion coefficient and the other for more regular data, in combination with a regularity result for the solutions. △ Less

Submitted 1 March, 2023; v1 submitted 21 March, 2022; originally announced March 2022.

Comments: Key words: COVID-19, compartmental model, partial differential equations, diffusion-reaction system, initial-boundary value problem, existence of solutions, uniqueness

MSC Class: 35Q92; 46N60; 92D30; 35K57; 35K55; 35K20

Residual-based error estimation and adaptivity for stabilized immersed isogeometric analysis using truncated hierarchical B-splines

Authors: Sai C. Divi, Pieter H. van Zuijlen, Tuong Hoang, Frits de Prenter, Ferdinando Auricchio, Alessandro Reali, E. Harald van Brummelen, Clemens V. Verhoosel

Abstract: We propose an adaptive mesh refinement strategy for immersed isogeometric analysis, with application to steady heat conduction and viscous flow problems. The proposed strategy is based on residual-based error estimation, which has been tailored to the immersed setting by the incorporation of appropriately scaled stabilization and boundary terms. Element-wise error indicators are elaborated for the… ▽ More We propose an adaptive mesh refinement strategy for immersed isogeometric analysis, with application to steady heat conduction and viscous flow problems. The proposed strategy is based on residual-based error estimation, which has been tailored to the immersed setting by the incorporation of appropriately scaled stabilization and boundary terms. Element-wise error indicators are elaborated for the Laplace and Stokes problems, and a THB-spline-based local mesh refinement strategy is proposed. The error estimation .and adaptivity procedure is applied to a series of benchmark problems, demonstrating the suitability of the technique for a range of smooth and non-smooth problems. The adaptivity strategy is also integrated in a scan-based analysis workflow, capable of generating reliable, error-controlled, results from scan data, without the need for extensive user interactions or interventions. △ Less

Submitted 17 February, 2022; originally announced February 2022.

Comments: Submitted to Journal of Mechanics

Two-level Method Part-scale Thermal Analysis of Laser Powder Bed Fusion Additive Manufacturing

Authors: Massimo Carraturo, Alex Viguerie, Alessandro Reali, Ferdinando Auricchio

Abstract: Numerical simulations of a complete laser powder bed fusion (LPBF) additive manufacturing (AM) process are extremely challenging or even impossible to achieve without a radical model reduction of the complex physical phenomena occurring during the process. However, even when we adopt reduced model with simplified physics, the complex geometries of parts usually produced by LPBF AM processes make t… ▽ More Numerical simulations of a complete laser powder bed fusion (LPBF) additive manufacturing (AM) process are extremely challenging or even impossible to achieve without a radical model reduction of the complex physical phenomena occurring during the process. However, even when we adopt reduced model with simplified physics, the complex geometries of parts usually produced by LPBF AM processes make this kind of analysis computationally expensive. In fact, small geometrical features - which might be generated when the part is design following the principal of the so-called design for AM, for instance, by means of topology optimization procedures - often require complex conformal meshes. Immersed boundary methods seem to offer a valid alternative to deal with this kind of complexity. The two-level method lies within this family of numerical methods and presents a very flexible tool to deal with multi-scale problems. In this contribution, we apply the recently introduced two-level method to part-scale thermal analysis of LPBF manufactured components, first validating the proposed part-scale model with respect to experimental measurements from the literature and then applying the presented numerical framework to simulate a complete LPBF process of a topologically optimized structure, showing the capability of the method to easily deal with complex geometrical features. △ Less

Submitted 17 December, 2021; originally announced December 2021.

A theoretical and numerical analysis of a Dirichlet-Neumann domain decomposition method for diffusion problems in heterogeneous media

Authors: Alex Viguerie, Silvia Bertoluzza, Alessandro Veneziani, Ferdinando Auricchio

Abstract: Problems with localized nonhomogeneous material properties present well-known challenges for numerical simulations. In particular, such problems may feature large differences in length scales, causing difficulties with meshing and preconditioning. These difficulties are increased if the region of localized dynamics changes in time. Overlapping domain decomposition methods, which split the problem… ▽ More Problems with localized nonhomogeneous material properties present well-known challenges for numerical simulations. In particular, such problems may feature large differences in length scales, causing difficulties with meshing and preconditioning. These difficulties are increased if the region of localized dynamics changes in time. Overlapping domain decomposition methods, which split the problem at the continuous level, show promise due to their ease of implementation and computational efficiency. Accordingly, the present work aims to further develop the mathematical theory of such methods at both the continuous and discrete levels. For the continuous formulation of the problem, we provide a full convergence analysis. For the discrete problem, we show how the described method may be interpreted as a Gauss-Seidel scheme or as a Neumann series approximation, establishing a convergence criterion in terms of the spectral radius of the system. We then provide a spectral scaling argument and provide numerical evidence for its justification. △ Less

Submitted 22 November, 2021; originally announced November 2021.

A spatiotemporal two-level method for high-fidelity thermal analysis of laserpowder bed fusion

Authors: Alex Viguerie, Massimo Carraturo, Alessandro Reali, Ferdinando Auricchio

Abstract: Numerical simulation of the laser powder bed fusion (LPBF) procedure for additive manufacturing (AM) is difficult due to the presence of multiple scales in both time and space, ranging from the part scale (order of millimeters/seconds) to the powder scale (order of microns/milliseconds). This difficulty is compounded by the fact that the regions of small-scale behavior are not fixed, but change in… ▽ More Numerical simulation of the laser powder bed fusion (LPBF) procedure for additive manufacturing (AM) is difficult due to the presence of multiple scales in both time and space, ranging from the part scale (order of millimeters/seconds) to the powder scale (order of microns/milliseconds). This difficulty is compounded by the fact that the regions of small-scale behavior are not fixed, but change in time as the geometry is produced. While much work in recent years has been focused on resolving the problem of multiple scales in space, there has been less work done on multiscale approaches for the temporal discretization of LPBF problems. In the present contribution, we extend on a previously introduced two-level method in space by combining it with a multiscale time integration method. The unique transfer of information through the transmission conditions allows for interaction between the space and time scales while reducing computational costs. At the same time, all of the advantages of the two-level method in space (namely its geometrical flexibility and the ease in which one may deploy structured, uniform meshes) remain intact. Adopting the proposed multiscale time integration scheme, we observe a computational speed-up by a factor $\times 2.44$ compared to the same two-level approach with uniform time integration, when simulating a laser source traveling on a bare plate of nickel-based superalloy material following an alternating scan path of fifty laser tracks. △ Less

Submitted 25 October, 2021; originally announced October 2021.

Topology-preserving Scan-based Immersed Isogeometric Analysis

Authors: S. C. Divi, C. V. Verhoosel, F. Auricchio, A. Reali, E. H. van Brummelen

Abstract: To exploit the advantageous properties of isogeometric analysis (IGA) in a scan-based setting, it is important to extract a smooth geometric domain from the scan data (e.g., voxel data). IGA-suitable domains can be constructed by convoluting the grayscale data using B-splines. A negative side-effect of this convolution technique is, however, that it can induce topological changes in the process of… ▽ More To exploit the advantageous properties of isogeometric analysis (IGA) in a scan-based setting, it is important to extract a smooth geometric domain from the scan data (e.g., voxel data). IGA-suitable domains can be constructed by convoluting the grayscale data using B-splines. A negative side-effect of this convolution technique is, however, that it can induce topological changes in the process of smoothing when features with a size similar to that of the voxels are encountered. This manuscript presents an enhanced B-spline-based segmentation procedure using a refinement strategy based on truncated hierarchical (TH)B-splines. A Fourier analysis is presented to explain the effectiveness of local grayscale function refinement in repairing topological anomalies. A moving-window-based topological anomaly detection algorithm is proposed to identify regions in which the grayscale function refinements must be performed. The criterion to identify topological anomalies is based on the Euler characteristic, giving it the capability to distinguish between topological and shape changes. The proposed topology-preserving THB-spline image segmentation strategy is studied using a range of test cases. These tests pertain to both the segmentation procedure itself, and its application in an immersed IGA setting. △ Less

Submitted 23 June, 2021; originally announced June 2021.

Comments: Submitted to CMAME

Mixed variational formulations for structural topology optimization based on the phase-field approach

Authors: Michele Marino, Ferdinando Auricchio, Alessandro Reali, Elisabetta Rocca, Ulisse Stefanelli

Abstract: We propose a variational principle combining a phase-field functional for structural topology optimization with a mixed (three-field) Hu-Washizu functional, then including directly in the formulation equilibrium, constitutive, and compatibility equations. The resulting mixed variational functional is then specialized to derive a classical topology optimization formulation (where the amount of mate… ▽ More We propose a variational principle combining a phase-field functional for structural topology optimization with a mixed (three-field) Hu-Washizu functional, then including directly in the formulation equilibrium, constitutive, and compatibility equations. The resulting mixed variational functional is then specialized to derive a classical topology optimization formulation (where the amount of material to be distributed is an \emph{a priori} assigned quantity acting as a global constraint for the problem) as well as a novel topology optimization formulation (where the amount of material to be distributed is minimized, hence with no pre imposed constraint for the problem). Both formulations are numerically solved by implementing a mixed finite element scheme, with the second approach avoiding the introduction of a global constraint, hence respecting the convenient local nature of the finite element discretization. Furthermore, within the proposed approach it is possible to obtain guidelines for settings proper values of phase-field-related simulation parameters and, thanks to the combined phase-field and Hu-Washizu rationale, a monolithic algorithm solution scheme can be easily adopted. An insightful and extensive numerical investigation results in a detailed convergence study and a discussion on the obtained final designs. The numerical results clearly highlight differences between the two formulations as well as advantages related to the monolithic solution strategy; numerical investigations address both two-dimensional and three-dimensional applications. △ Less

Submitted 1 March, 2021; originally announced March 2021.

Comments: 35 pages, 14 figures

Bending behavior of additively manufactured lattice structures: numerical characterization and experimental validation

Authors: Nina Korshunova, Gianluca Alaimo, Seyyed Bahram Hosseini, Massimo Carraturo, Alessandro Reali, Jarkko Niiranen, Ferdinando Auricchio, Ernst Rank, Stefan Kollmannsberger

Abstract: Selective Laser Melting (SLM) technology has undergone significant development in the past years providing unique flexibility for the fabrication of complex metamaterials such as octet-truss lattices. However, the microstructure of the final parts can exhibit significant variations due to the high complexity of the manufacturing process. Consequently, the mechanical behavior of these lattices is s… ▽ More Selective Laser Melting (SLM) technology has undergone significant development in the past years providing unique flexibility for the fabrication of complex metamaterials such as octet-truss lattices. However, the microstructure of the final parts can exhibit significant variations due to the high complexity of the manufacturing process. Consequently, the mechanical behavior of these lattices is strongly dependent on the process-induced defects, raising the importance on the incorporation of as-manufactured geometries into the computational structural analysis. This, in turn, challenges the traditional mesh-conforming methods making the computational costs prohibitively large. In the present work, an immersed image-to-analysis framework is applied to efficiently evaluate the bending behavior of AM lattices. To this end, we employ the Finite Cell Method (FCM) to perform a three-dimensional numerical analysis of the three-point bending test of a lattice structure and compare the as-designed to as-manufactured effective properties. Furthermore, we undertake a comprehensive study on the applicability of dimensionally reduced beam models to the prediction of the bending behavior of lattice beams and validate classical and strain gradient beam theories applied in combination with the FCM. The numerical findings suggest that the SLM octet-truss lattices exhibit size effects, thus, requiring a flexible framework to incorporate high-order continuum theories. △ Less

Submitted 22 January, 2021; originally announced January 2021.

Studies on Coupled Flight Dynamics and Aeroelasticity of a Prandtlplane Configuration

Authors: Rocco Bombardieri, Rauno Cavallaro, Rodrigo Castellanos, Francesco Auricchio

Abstract: The PrandtlPlane aircraft has been recently considered as a possible candidate to foster the ambition of a greener aviation. Despite the relevant amount of research carried out in the last years, several aspects of this novel configuration still need further insight to pave the way to future applications. Among them, the coupled flight-dynamic and aeroelastic response is addressed in this work by… ▽ More The PrandtlPlane aircraft has been recently considered as a possible candidate to foster the ambition of a greener aviation. Despite the relevant amount of research carried out in the last years, several aspects of this novel configuration still need further insight to pave the way to future applications. Among them, the coupled flight-dynamic and aeroelastic response is addressed in this work by means of a dedicated in-house framework. For the evaluation of the aerodynamic forces, an enhanced Doublet Lattice Method, able to take into account terms typically neglected by classic formulations, is employed. First, flight-dynamic aspects are considered, showing how effects of interaction between Short Period or Dutch Roll with elastic modes remarkably deteriorate the flying qualities. Then, focus is on the aeroelastic stability of the aircraft. As observed also in previous literature efforts on this configuration, flutter onset is considerably different when considering the aircraft being free in the air or fixed in space. Thanks to the adopted formulation it is shown how, for this PrandtlPlane, the aerodynamic coupling of elastic and rigid modes has a beneficial effect on flutter onset. However, the different modal properties, consequence of the diverse boundary conditions, when switching from fixed-in-space to free-flying aircraft, also play a relevant role in determining the occurrence of flutter. Whereas for the longitudinal case both effects contribute increasing flutter speed, for the lateral-directional case the variation in modal properties has a detrimental and dominating effect, leading to a flutter speedwell within the flight envelope.Finally,the work discusses the contributions of the additional terms modeled by the enhanced Doublet Lattice Method, showing how they induce a considerable effect when modeling the flight dynamics of the flexible aircraft. △ Less

Submitted 11 January, 2021; originally announced January 2021.

Image-based numerical characterization and experimental validation of tensile behavior of octet-truss lattice structures

Authors: Nina Korshunova, Gianluca Alaimo, Seyyed Bahram Hosseini, Massimo Carraturo, Alessandro Reali, Jarkko Niiranen, Ferdinando Auricchio, Ernst Rank, Stefan Kollmannsberger

Abstract: The production of lightweight metal lattice structures has received much attention due to the recent developments in additive manufacturing (AM). The design flexibility comes, however, with the complexity of the underlying physics. In fact, metal additive manufacturing introduces process-induced geometrical defects that mainly result in deviations of the effective geometry from the nominal one. Th… ▽ More The production of lightweight metal lattice structures has received much attention due to the recent developments in additive manufacturing (AM). The design flexibility comes, however, with the complexity of the underlying physics. In fact, metal additive manufacturing introduces process-induced geometrical defects that mainly result in deviations of the effective geometry from the nominal one. This change in the final printed shape is the primary cause of the gap between the as-designed and as-manufactured mechanical behavior of AM products. Thus, the possibility to incorporate the precise manufactured geometries into the computational analysis is crucial for the quality and performance assessment of the final parts. Computed tomography (CT) is an accurate method for the acquisition of the manufactured shape. However, it is often not feasible to integrate the CT-based geometrical information into the traditional computational analysis due to the complexity of the meshing procedure for such high-resolution geometrical models and the prohibitive numerical costs. In this work, an embedded numerical framework is applied to efficiently simulate and compare the mechanical behavior of as-designed to as-manufactured octet-truss lattice structures. The parts are produced using laser powder bed fusion (LPBF). Employing an immersed boundary method, namely the Finite Cell Method (FCM), we perform direct numerical simulations (DNS) and first-order numerical homogenization analysis of a tensile test for a 3D printed octet-truss structure. Numerical results based on CT scan (as-manufactured geometry) show an excellent agreement with experimental measurements, whereas both DNS and first-order numerical homogenization performed directly on the 3D virtual model (as-designed geometry) of the structure show a significant deviation from experimental data. △ Less

Submitted 11 March, 2021; v1 submitted 14 December, 2020; originally announced December 2020.

Simulating the spread of COVID-19 via spatially-resolved susceptible-exposed-infected-recovered-deceased (SEIRD) model with heterogeneous diffusion

Authors: Alex Viguerie, Guillermo Lorenzo, Ferdinando Auricchio, Davide Baroli, Thomas J. R. Hughes, Alessia Patton, Alessandro Reali, Thomas E. Yankeelov, Alessandro Veneziani

Abstract: We present an early version of a Susceptible-Exposed-Infected-Recovered-Deceased (SEIRD) mathematical model based on partial differential equations coupled with a heterogeneous diffusion model. The model describes the spatio-temporal spread of the COVID-19 pandemic, and aims to capture dynamics also based on human habits and geographical features. To test the model, we compare the outputs generate… ▽ More We present an early version of a Susceptible-Exposed-Infected-Recovered-Deceased (SEIRD) mathematical model based on partial differential equations coupled with a heterogeneous diffusion model. The model describes the spatio-temporal spread of the COVID-19 pandemic, and aims to capture dynamics also based on human habits and geographical features. To test the model, we compare the outputs generated by a finite-element solver with measured data over the Italian region of Lombardy, which has been heavily impacted by this crisis between February and April 2020. Our results show a strong qualitative agreement between the simulated forecast of the spatio-temporal COVID-19 spread in Lombardy and epidemiological data collected at the municipality level. Additional simulations exploring alternative scenarios for the relaxation of lockdown restrictions suggest that reopening strategies should account for local population densities and the specific dynamics of the contagion. Thus, we argue that data-driven simulations of our model could ultimately inform health authorities to design effective pandemic-arresting measures and anticipate the geographical allocation of crucial medical resources. △ Less

Submitted 11 May, 2020; originally announced May 2020.

Comments: 7 pages, 3 figures

ACM Class: G.1.8; G.1.10

Automatic differentiation for solid mechanics

Authors: Andrea Vigliotti, Ferdinando Auricchio

Abstract: Automatic differentiation (AD) is an ensemble of techniques that allow to evaluate accurate numerical derivatives of a mathematical function expressed in a computer programming language. In this paper we use AD for stating and solving solid mechanics problems. Given a finite element discretization of the domain, we evaluate the free energy of the solid as the integral of its strain energy density,… ▽ More Automatic differentiation (AD) is an ensemble of techniques that allow to evaluate accurate numerical derivatives of a mathematical function expressed in a computer programming language. In this paper we use AD for stating and solving solid mechanics problems. Given a finite element discretization of the domain, we evaluate the free energy of the solid as the integral of its strain energy density, and we make use of AD for directly obtaining the residual force vector and the tangent stiffness matrix of the problem, as the gradient and the Hessian of the free energy respectively. The result is a remarkable simplification in the statement and the solution of complex problems involving non trivial constraints systems and both geometrical and material non linearities. Together with the continuum mechanics theoretical basis, and with a description of the specific AD technique adopted, the paper illustrates the solution of a number of solid mechanics problems, with the aim of presenting a convenient numerical implementation approach, made easily available by recent programming languages, to the solid mechanics community. △ Less

Submitted 21 January, 2020; originally announced January 2020.

Comments: 30 pages, 9 figures, 2 appendices, accepted on Archives of Computational Methods in Engineering

Error-estimate-based Adaptive Integration For Immersed Isogeometric Analysis

Authors: Sai C Divi, Clemens V Verhoosel, Ferdinando Auricchio, Alessandro Reali, E Harald van Brummelen

Abstract: The Finite Cell Method (FCM) together with Isogeometric analysis (IGA) has been applied successfully in various problems in solid mechanics, in image-based analysis, fluid-structure interaction and in many other applications. A challenging aspect of the isogeometric finite cell method is the integration of cut cells. In particular in three-dimensional simulations the computational effort associate… ▽ More The Finite Cell Method (FCM) together with Isogeometric analysis (IGA) has been applied successfully in various problems in solid mechanics, in image-based analysis, fluid-structure interaction and in many other applications. A challenging aspect of the isogeometric finite cell method is the integration of cut cells. In particular in three-dimensional simulations the computational effort associated with integration can be the critical component of a simulation. A myriad of integration strategies has been proposed over the past years to ameliorate the difficulties associated with integration, but a general optimal integration framework that suits a broad class of engineering problems is not yet available. In this contribution we provide a thorough investigation of the accuracy and computational effort of the octree integration scheme. We quantify the contribution of the integration error using the theoretical basis provided by Strang's first lemma. Based on this study we propose an error-estimate-based adaptive integration procedure for immersed isogeometric analysis. Additionally, we present a detailed numerical investigation of the proposed optimal integration algorithm and its application to immersed isogeometric analysis using two- and three-dimensional linear elasticity problems. △ Less

Submitted 26 November, 2019; originally announced November 2019.

Comments: To CAMWA

A Fat boundary-type method for localized nonhomogeneous material problems

Authors: Alex Viguerie, Silvia Bertoluzza, Ferdinando Auricchio

Abstract: Problems with localized nonhomogeneous material properties arise frequently in many applications and are a well-known source of difficulty in numerical simulations. In certain applications (including additive manufacturing), the physics of the problem may be considerably more complicated in relatively small portions of the domain, requiring a significantly finer local mesh compared to elsewhere in… ▽ More Problems with localized nonhomogeneous material properties arise frequently in many applications and are a well-known source of difficulty in numerical simulations. In certain applications (including additive manufacturing), the physics of the problem may be considerably more complicated in relatively small portions of the domain, requiring a significantly finer local mesh compared to elsewhere in the domain. This can make the use of a uniform mesh numerically unfeasible. While nonuniform meshes can be employed, they may be challenging to generate (particularly for regions with complex boundaries) and more difficult to precondition. The problem becomes even more prohibitive when the region requiring a finer-level mesh changes in time, requiring the introduction of refinement and derefinement techniques. To address the aforementioned challenges, we employ a technique related to the Fat boundary method as a possible alternative. We analyze the proposed methodology, from a mathematical point of view and validate our findings on two-dimensional numerical tests. △ Less

Submitted 21 October, 2019; originally announced October 2019.

Structural multiscale topology optimization with stress constraint for additive manufacturing

Authors: Ferdinando Auricchio, Elena Bonetti, Massimo Carraturo, Dietmar Hömberg, Alessandro Reali, Elisabetta Rocca

Abstract: In this paper a phase-field approach for structural topology optimization for a 3D-printing process which includes stress constraint and potentially multiple materials or multiscales is analyzed. First order necessary optimality conditions are rigorously derived and a numerical algorithm which implements the method is presented. A sensitivity study with respect to some parameters is conducted for… ▽ More In this paper a phase-field approach for structural topology optimization for a 3D-printing process which includes stress constraint and potentially multiple materials or multiscales is analyzed. First order necessary optimality conditions are rigorously derived and a numerical algorithm which implements the method is presented. A sensitivity study with respect to some parameters is conducted for a two-dimensional cantilever beam problem. Finally, a possible workflow to obtain a 3D-printed object from the numerical solutions is described and the final structure is printed using a fused deposition modeling (FDM) 3D printer. △ Less

Submitted 15 July, 2019; originally announced July 2019.

MSC Class: 74P05; 49M05; 74B99

A novel layered topology of auxetic materials based on the tetrachiral honeycomb microstructure

Authors: Ferdinando Auricchio, Andrea Bacigalupo, Luigi Gambarotta, Marco Lepidi, Simone Morganti, Francesca Vadala

Abstract: Microstructured honeycomb materials may exhibit exotic, extreme and tailorable mechanical properties, suited for innovative technological applications in a variety of modern engineering fields. The paper is focused on analysing the directional auxeticity of tetrachiral materials, through analytical, numerical and experimental methods. Theoretical predictions about the global elastic properties hav… ▽ More Microstructured honeycomb materials may exhibit exotic, extreme and tailorable mechanical properties, suited for innovative technological applications in a variety of modern engineering fields. The paper is focused on analysing the directional auxeticity of tetrachiral materials, through analytical, numerical and experimental methods. Theoretical predictions about the global elastic properties have been successfully validated by performing tensile laboratory tests on tetrachiral samples, realized with high precision 3D printing technologies. Inspired by the kinematic behaviour of the tetrachiral material, a newly-design bi-layered topology, referred to as bi-tetrachiral material, has been theoretically conceived and mechanically modelled. The novel topology virtuously exploits the mutual collaboration between two tetrachiral layers with opposite chiralities. The bi-tetrachiral material has been verified to outperform the tetrachiral material in terms of global Young modulus and, as major achievement, to exhibit a remarkable auxetic behaviour. Specifically, experimental results, confirmed by parametric analytical and computational analyses, have highlighted the effective possibility to attain strongly negative Poisson ratios, identified as a peculiar global elastic property of the novel bi-layered topology. △ Less

Submitted 26 March, 2019; originally announced March 2019.

Accurate prediction of melt pool shapes in laser powder bed fusion by the non-linear temperature equation including phase changes - isotropic versus anisotropic conductivity

Authors: Stefan Kollmannsberger, Massimo Carraturo, Alessandro Reali, Ferdinando Auricchio

Abstract: In this contribution, we validate a physical model based on a transient temperature equation (including latent heat) w.r.t. the experimental set AMB2018-02 provided within the additive manufacturing benchmark series, established at the National Institute of Standards and Technology, USA. We aim at predicting the following quantities of interest: width, depth, and length of the melt pool by numeric… ▽ More In this contribution, we validate a physical model based on a transient temperature equation (including latent heat) w.r.t. the experimental set AMB2018-02 provided within the additive manufacturing benchmark series, established at the National Institute of Standards and Technology, USA. We aim at predicting the following quantities of interest: width, depth, and length of the melt pool by numerical simulation and report also on the obtainable numerical results of the cooling rate. We first assume the laser to posses a double ellipsoidal shape and demonstrate that a well calibrated, purely thermal model based on isotropic thermal conductivity is able to predict all the quantities of interest, up to a deviation of maximum 7.3\% from the experimentally measured values. However, it is interesting to observe that if we directly introduce, whenever available, the measured laser profile in the model (instead of the double ellipsoidal shape) the investigated model returns a deviation of 19.3\% from the experimental values. This motivates a model update by introducing anisotropic conductivity, which is intended to be a simplistic model for heat material convection inside the melt pool. Such an anisotropic model enables the prediction of all quantities of interest mentioned above with a maximum deviation from the experimental values of 6.5\%. We note that, although more predictive, the anisotropic model induces only a marginal increase in computational complexity. △ Less

Submitted 21 March, 2019; originally announced March 2019.

Journal ref: Integrating Materials and Manufacturing Innovation 2019, p 1-11

Modeling the non-trivial behavior of anisotropic beams: a simple Timoshenko beam with enhanced stress recovery and constitutive relations

Authors: Giuseppe Balduzzi, Simone Morganti, Josef Fussl, Mehdi Aminbaghai, Alessandro Reali, Ferdinando Auricchio

Abstract: This paper analyzes the non-trivial influence of the material anisotropy on the structural behavior of an anisotropic multilayer planar beam. Indeed, analytical results available in literature are limited to homogeneous beams and several aspects has not been addressed yet, impeding an in-depth understanding of the mechanical response of anisotropic structural elements. This paper proposes an effec… ▽ More This paper analyzes the non-trivial influence of the material anisotropy on the structural behavior of an anisotropic multilayer planar beam. Indeed, analytical results available in literature are limited to homogeneous beams and several aspects has not been addressed yet, impeding an in-depth understanding of the mechanical response of anisotropic structural elements. This paper proposes an effective recovery of stress distribution and an energetically consistent evaluation of constitutive relations to be used within a planar Timoshenko beam model. The resulting structural-analysis tool highlights the following peculiarities of anisotropic beams: (i) the axial stress explicitly depends on transversal internal force, which can weigh up to 30% on the maximal magnitude of axial stress, and (ii) the anisotropy influences the beam displacements more than standard shear deformation and even for extremely slender beams. A rigorous comparison with analytical and accurate 2D Finite Element solutions confirms the accuracy of the proposed approach that leads to errors exceptionally greater than 5%. △ Less

Submitted 1 March, 2019; originally announced March 2019.

Graded-material Design based on Phase-field and Topology Optimization

Authors: Massimo Carraturo, Elisabetta Rocca, Elena Bonetti, Dietmar Hömberg, Alessandro Reali, Ferdinando Auricchio

Abstract: In the present work we introduce a novel graded-material design based on phase-field and topology optimization. The main novelty of this work comes from the introduction of an additional phase-field variable in the classical single-material phase-field topology optimization algorithm. This new variable is used to grade the material properties in a continuous fashion. Two different numerical exampl… ▽ More In the present work we introduce a novel graded-material design based on phase-field and topology optimization. The main novelty of this work comes from the introduction of an additional phase-field variable in the classical single-material phase-field topology optimization algorithm. This new variable is used to grade the material properties in a continuous fashion. Two different numerical examples are discussed, in both of them we perform sensitivity studies to asses the effects of different model parameters onto the resulting structure. From the presented results we can observe that the proposed algorithm adds additional freedom in the design, exploiting the higher flexibility coming from additive manufacturing technology. △ Less

Submitted 31 May, 2019; v1 submitted 17 November, 2018; originally announced November 2018.

Skeleton-stabilized ImmersoGeometric Analysis for incompressible viscous flow problems

Authors: Tuong Hoang, Clemens V. Verhoosel, Chao-Zhong Qin, Ferdinando Auricchio, Alessandro Reali, E. Harald van Brummelen

Abstract: A Skeleton-stabilized ImmersoGeometric Analysis technique is proposed for incompressible viscous flow problems with moderate Reynolds number. The proposed formulation fits within the framework of the finite cell method, where essential boundary conditions are imposed weakly using a Nitsche-type method. The key idea of the proposed formulation is to stabilize the jumps of high-order derivatives of… ▽ More A Skeleton-stabilized ImmersoGeometric Analysis technique is proposed for incompressible viscous flow problems with moderate Reynolds number. The proposed formulation fits within the framework of the finite cell method, where essential boundary conditions are imposed weakly using a Nitsche-type method. The key idea of the proposed formulation is to stabilize the jumps of high-order derivatives of variables over the skeleton of the background mesh. The formulation allows the use of identical finite-dimensional spaces for the approximation of the pressure and velocity fields in immersed domains. The stability issues observed for inf-sup stable discretizations of immersed incompressible flow problems are avoided with this formulation. For B-spline basis functions of degree $k$ with highest regularity, only the derivative of order $k$ has to be controlled, which requires specification of only a single stabilization parameter for the pressure field. The Stokes and Navier-Stokes equations are studied numerically in two and three dimensions using various immersed test cases. Oscillation-free solutions and high-order optimal convergence rates can be obtained. The formulation is shown to be stable even in limit cases where almost every elements of the physical domain is cut, and hence it does not require the existence of interior cells. In terms of the sparsity pattern, the algebraic system has a considerably smaller stencil than counterpart approaches based on Lagrange basis functions. This important property makes the proposed skeleton-stabilized technique computationally practical. To demonstrate the stability and robustness of the method, we perform a simulation of fluid flow through a porous medium, of which the geometry is directly extracted from 3D $μ{CT}$ scan data. △ Less

Submitted 19 July, 2018; originally announced July 2018.

Computational methods in cardiovascular mechanics

Authors: Ferdinando Auricchio, Michele Conti, Adrian Lefieux, Simone Morganti, Alessandro Reali, Gianluigi Rozza, Alessandro Veneziani

Abstract: The introduction of computational models in cardiovascular sciences has been progressively bringing new and unique tools for the investigation of the physiopathology. Together with the dramatic improvement of imaging and measuring devices on one side, and of computational architectures on the other one, mathematical and numerical models have provided a new, clearly noninvasive, approach for unders… ▽ More The introduction of computational models in cardiovascular sciences has been progressively bringing new and unique tools for the investigation of the physiopathology. Together with the dramatic improvement of imaging and measuring devices on one side, and of computational architectures on the other one, mathematical and numerical models have provided a new, clearly noninvasive, approach for understanding not only basic mechanisms but also patient-specific conditions, and for supporting the design and the development of new therapeutic options. The terminology in silico is, nowadays, commonly accepted for indicating this new source of knowledge added to traditional in vitro and in vivo investigations. The advantages of in silico methodologies are basically the low cost in terms of infrastructures and facilities, the reduced invasiveness and, in general, the intrinsic predictive capabilities based on the use of mathematical models. The disadvantages are generally identified in the distance between the real cases and their virtual counterpart required by the conceptual modeling that can be detrimental for the reliability of numerical simulations. △ Less

Submitted 8 March, 2018; originally announced March 2018.

Comments: 54 pages, Book Chapter

Skeleton-stabilized IsoGeometric Analysis: High-regularity Interior-Penalty methods for incompressible viscous flow problems

Authors: Tuong Hoang, Clemens V. Verhoosel, Ferdinando Auricchio, E. Harald van Brummelen, Alessandro Reali

Abstract: A Skeleton-stabilized IsoGeometric Analysis (SIGA) technique is proposed for incompressible viscous flow problems with moderate Reynolds number. The proposed method allows utilizing identical finite dimensional spaces (with arbitrary B-splines/NURBS order and regularity) for the approximation of the pressure and velocity components. The key idea is to stabilize the jumps of high-order derivatives… ▽ More A Skeleton-stabilized IsoGeometric Analysis (SIGA) technique is proposed for incompressible viscous flow problems with moderate Reynolds number. The proposed method allows utilizing identical finite dimensional spaces (with arbitrary B-splines/NURBS order and regularity) for the approximation of the pressure and velocity components. The key idea is to stabilize the jumps of high-order derivatives of variables over the skeleton of the mesh. For B-splines/NURBS basis functions of degree $k$ with $C^α$-regularity ($0 \leq α< k$), only the derivative of order $α+1$ has to be controlled. This stabilization technique thus can be viewed as a high-regularity generalization of the (Continuous) Interior-Penalty Finite Element Method. Numerical experiments are performed for the Stokes and Navier-Stokes equations in two and three dimensions. Oscillation-free solutions and optimal convergence rates are obtained. In terms of the sparsity pattern of the algebraic system, we demonstrate that the block matrix associated with the stabilization term has a considerably smaller bandwidth when using B-splines than when using Lagrange basis functions, even in the case of $C^0$-continuity. This important property makes the proposed isogeometric framework practical from a computational effort point of view. △ Less

Submitted 19 July, 2018; v1 submitted 13 November, 2017; originally announced November 2017.

Journal ref: Comput. Methods Appl. Mech. Engrg. (2018)

A cost-effective isogeometric approach for composite plates based on a stress recovery procedure

Authors: John-Eric Dufour, Pablo Antolin, Giancarlo Sangalli, Ferdinando Auricchio, Alessandro Reali

Abstract: This paper introduces a cost-effective strategy to simulate the behavior of laminated plates by means of isogeometric 3D solid elements. Exploiting the high continuity of spline functions and their properties, a proper out-of-plane stress state is recovered from a coarse displacement solution using a post-processing step based on the enforcement of equilibrium in strong form. Appealing results are… ▽ More This paper introduces a cost-effective strategy to simulate the behavior of laminated plates by means of isogeometric 3D solid elements. Exploiting the high continuity of spline functions and their properties, a proper out-of-plane stress state is recovered from a coarse displacement solution using a post-processing step based on the enforcement of equilibrium in strong form. Appealing results are obtained and the method is shown to be particularly Peffective on slender composite stacks with a large number of layers. △ Less

Submitted 21 April, 2017; originally announced April 2017.

A new 3D macroscopic model for shape memory alloys describing martensite reorientation

Authors: Ferdinando Auricchio, Elena Bonetti

Abstract: In this paper we introduce a 3D phenomenological model for shape memory behavior, accounting for: martensite reorientation, asymmetric response of the material to tension/compression, different kinetics between forward and reverse phase transformation. We combine two modeling approaches using scalar and tensorial internal variables. Indeed, we use volume proportions of different configurations of… ▽ More In this paper we introduce a 3D phenomenological model for shape memory behavior, accounting for: martensite reorientation, asymmetric response of the material to tension/compression, different kinetics between forward and reverse phase transformation. We combine two modeling approaches using scalar and tensorial internal variables. Indeed, we use volume proportions of different configurations of the crystal lattice (austenite and two variants of martensite) as scalar internal variables and the preferred direction of stress-induced martensite as tensorial internal variable. Then, we derive evolution equations by a generalization of the principle of virtual powers, including microforces and micromovements responsible for phase transformation. In addition, we prescribe an evolution law for phase proportions ensuring different evolution laws during forward and reverse transformation of the oriented martensite. △ Less

Submitted 3 October, 2012; originally announced October 2012.

A rate-independent model for the isothermal quasi-static evolution of shape-memory materials

Authors: Ferdinando Auricchio, Alexander Mielke, Ulisse Stefanelli

Abstract: This note addresses a three-dimensional model for isothermal stress-induced transformation in shape-memory polycrystalline materials. We treat the problem within the framework of the energetic formulation of rate-independent processes and investigate existence and continuous dependence issues at both the constitutive relation and quasi-static evolution level. Moreover, we focus on time and space… ▽ More This note addresses a three-dimensional model for isothermal stress-induced transformation in shape-memory polycrystalline materials. We treat the problem within the framework of the energetic formulation of rate-independent processes and investigate existence and continuous dependence issues at both the constitutive relation and quasi-static evolution level. Moreover, we focus on time and space approximation as well as on regularization and parameter asymptotics. △ Less

Submitted 31 August, 2007; originally announced August 2007.

Comments: 33 pages, 3 figures

MSC Class: 74C05; 49J40