Source author record

Laura De Lorenzis

Laura De Lorenzis appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

9works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

9 published item(s)

preprint2025arXiv

Variational phase-field modeling of cohesive fracture with flexibly tunable strength surface

Variational phase-field models of brittle fracture are powerful tools for studying Griffith-type crack propagation in complex scenarios. However, as approximations of Griffith's theory-which does not incorporate a strength criterion-these models lack flexibility in prescribing material-specific strength surfaces. Consequently, they struggle to accurately capture crack nucleation under multiaxial stress conditions. In this paper, inspired by Alessi et al. (2014), we propose a variational phase-field model that approximates cohesive fracture. The model accommodates an arbitrary (convex) strength surface, independent of the regularization length scale, and allows for flexible tuning of the cohesive response. Our formulation results in sharp cohesive cracks and naturally enforces a sharp non-interpenetration condition, thereby eliminating the need for additional energy decomposition strategies. It inherently satisfies stress softening and produces "crack-like" residual stresses by construction. To ensure strain hardening, the ratio of the regularization length to the material's cohesive length must be sufficiently small; however, if crack nucleation is desired, this ratio must also be large enough to make the homogeneous damaged state unstable. We investigate the model in one and three dimensions, establishing first- and second-order stability results. The theoretical findings are validated through numerical simulations using the finite element method, employing standard discretization and solution techniques.

preprint2025arXiv

Variational phase-field modeling of fracture and fatigue in shape memory alloys

We propose a novel variational phase-field model for fracture and fatigue in pseudoelastic shape memory alloys (SMAs). The model, developed in a one-dimensional setting, builds upon the Auricchio-Petrini constitutive formulation for SMAs and couples damage evolution with phase transformation. We study analytically and numerically the homogeneous and localization responses of a bar under both monotonic and cyclic loading, and we investigate various macroscopic behaviors by tuning the constitutive parameters. A key feature of the model is the introduction of a transformation strain limit, beyond which the material is fully martensitic and behaves elastically. This leads to a distinctive behavior in which the region of localized damage widens, yielding a delay of fracture. The capability of the model to predict the fatigue performance is demonstrated by simulating the uniaxial response of Ni-Ti multi-wire samples under different loading conditions. The results show promising agreement with experimental fatigue life data, enabling the discrimination between safe and critical loading scenarios.

preprint2022arXiv

Automated discovery of generalized standard material models with EUCLID

We extend the scope of our approach for unsupervised automated discovery of material laws (EUCLID) to the case of a material belonging to an unknown class of behavior. To this end, we leverage the theory of generalized standard materials, which encompasses a plethora of important constitutive classes. We show that, based only on full-field kinematic measurements and net reaction forces, EUCLID is able to automatically discover the two scalar thermodynamic potentials, namely, the Helmholtz free energy and the dissipation potential, which completely define the behavior of generalized standard materials. The a priori enforced constraint of convexity on these potentials guarantees by construction stability and thermodynamic consistency of the discovered model; balance of linear momentum acts as a fundamental constraint to replace the availability of stress-strain labeled pairs; sparsity promoting regularization enables the automatic selection of a small subset from a possibly large number of candidate model features and thus leads to a parsimonious, i.e., simple and interpretable, model. Importantly, since model features go hand in hand with the correspondingly active internal variables, sparse regression automatically induces a parsimonious selection of the few internal variables needed for an accurate but simple description of the material behavior. A fully automatic procedure leads to the selection of the hyperparameter controlling the weight of the sparsity promoting regularization term, in order to strike a user-defined balance between model accuracy and simplicity. By testing the method on synthetic data including artificial noise, we demonstrate that EUCLID is able to automatically discover the true hidden material model from a large catalog of constitutive classes, including elasticity, viscoelasticity, elastoplasticity, viscoplasticity, isotropic and kinematic hardening.

preprint2022arXiv

The isogeometric collocated contact surface approach

We propose a frictionless contact formulation for isogeometric analysis, which combines a collocated formulation for the contact surfaces with a standard Galerkin treatment of the bulk. We denote it as isogeometric Collocated Contact Surface (CCS) formulation. The approach is based on a simple pointwise enforcement of the contact constraints, performed in this study with the penalty method. Unlike pointwise (node-to-surface or point-to-surface) contact algorithms in the Galerkin framework, the CCS formulation passes the contact patch test to machine precision by naturally exploiting the favorable properties of isogeometric collocation. Compared with approaches where the discretization of both bulk and contact surfaces is based on collocation, the CCS approach does not need enhancements to remove oscillations for highly non-uniform meshes. With respect to integral contact approaches, the CCS algorithm is less expensive, easier to code and can be added to a pre-existing isogeometric analysis code with minimal effort. Numerical examples in both small and large deformations are investigated to compare the CCS approach with some available contact formulations and to demonstrate its accuracy.

preprint2021arXiv

Floating Isogeometric Analysis

We propose Floating Isogeometric Analysis (FLIGA), which extends the concepts of IGA to Lagrangian extreme deformation analysis. The method is based on a novel tensor-product construction of B-Splines for the update of the basis functions in one direction of the parametric space. With basis functions 'floating' deformation-dependently in this direction, mesh distortion is overcome for problems in which extreme deformations occur predominantly along the associated (possibly curved) physical axis. In doing so, we preserve the numerical advantages of splines over many meshless basis functions, while avoiding remeshing. We employ material point integration for numerical quadrature attributing a Lagrangian character to our technique. The paper introduces the method and reviews the fundamental properties of the FLIGA basis functions, including a numerical patch test. The performance of FLIGA is then numerically investigated on the benchmark of Newtonian and viscoelastic Taylor-Couette flow. Finally, we simulate a viscoelastic extrusion-based additive manufacturing process, which served as the original motivation for the new approach.

preprint2020arXiv

Data-driven fracture mechanics

We present a new data-driven paradigm for variational brittle fracture mechanics. The fracture-related material modeling assumptions are removed and the governing equations stemming from variational principles are combined with a set of discrete data points, leading to a model-free data-driven method of solution. The solution at a given load step is identified as the point within the data set that best satisfies either the Kuhn-Tucker conditions stemming from the variational fracture problem or global minimization of a suitable energy functional, leading to data-driven counterparts of both the local and the global minimization approaches of variational fracture mechanics. Both formulations are tested on different test configurations with and without noise and for Griffith and R-curve type fracture behavior.

preprint2015arXiv

Computational homogenization of fibrous piezoelectric materials

Flexible piezoelectric devices made of polymeric materials are widely used for micro- and nano-electro-mechanical systems. In particular, numerous recent applications concern energy harvesting. Due to the importance of computational modeling to understand the influence that microscale geometry and constitutive variables exert on the macroscopic behavior, a numerical approach is developed here for multiscale and multiphysics modeling of thin piezoelectric sheets made of aligned arrays of polymeric nanofibers, manufactured by electrospinning. At the microscale, the representative volume element consists in piezoelectric polymeric nanofibers, assumed to feature a piezoelastic behavior and subjected to electromechanical contact constraints. The latter are incorporated into the virtual work equations by formulating suitable electric, mechanical and coupling potentials and the constraints are enforced by using the penalty method. From the solution of the micro-scale boundary value problem, a suitable scale transition procedure leads to identifying the performance of a macroscopic thin piezoelectric shell element.

preprint2015arXiv

Cooperativity in the enhanced piezoelectric response of polymer nanowires

We provide a detailed insight into piezoelectric energy generation from arrays of polymer nanofibers. For sake of comparison, we firstly measure individual poly(vinylidenefluoride-co-trifluoroethylene) (P(VDF-TrFe)) fibers at well-defined levels of compressive stress. Under an applied load of 2 mN, single nanostructures generate a voltage of 0.45 mV. We show that under the same load conditions, fibers in dense arrays exhibit a voltage output higher by about two orders of magnitude. Numerical modelling studies demonstrate that the enhancement of the piezoelectric response is a general phenomenon associated to the electromechanical interaction among adjacent fibers, namely a cooperative effect depending on specific geometrical parameters. This establishes new design rules for next piezoelectric nano-generators and sensors.

preprint2014arXiv

A multiscale-multiphysics strategy for numerical modeling of thin piezoelectric sheets

Flexible piezoelectric devices made of polymeric materials are widely used for micro- and nano-electro-mechanical systems. In particular, numerous recent applications concern energy harvesting. Due to the importance of computational modeling to understand the influence that microscale geometry and constitutive variables exert on the macroscopic behavior, a numerical approach is developed here for multiscale and multiphysics modeling of piezoelectric materials made of aligned arrays of polymeric nanofibers. At the microscale, the representative volume element consists in piezoelectric polymeric nanofibers, assumed to feature a linear piezoelastic constitutive behavior and subjected to electromechanical contact constraints using the penalty method. To avoid the drawbacks associated with the non-smooth discretization of the master surface, a contact smoothing approach based on Bézier patches is extended to the multiphysics framework providing an improved continuity of the parameterization. The contact element contributions to the virtual work equations are included through suitable electric, mechanical and coupling potentials. From the solution of the micro-scale boundary value problem, a suitable scale transition procedure leads to the formulation of a macroscopic thin piezoelectric shell element.