Riferimenti Accademici
[1] C. B. Vreugdenhil, Numerical methods for shallow-water flow, volume 13, Springer Science & Business Media, 1994.
[2] M. A. Abdelrahman, On the shallow water equations, Z NATURFORSCH A 72 (2017) 873–879.
[3] M. Brocchini, N. Dodd, Nonlinear shallow water equation modelling for coastal engineering, J Water Port Coast 134 (2008) 104–120.
[4] P. L.-F. Liu, Model equations for wave propagations from deep to shallow water, Adv Coastal Ocean En: (Volume 1) (1995) 125–157.
[5] C. Altomare et al., Applicability of smoothed particle hydrodynamics for estimation of sea wave impact on coastal structures, Coast Eng. 96 (2015) 1–12.
[6] E. Onate et al., The particle finite element method—an overview, Int. J. Comp. Meth. 1 (2004) 267–307.
[7] S. R. Idelsohn, E. Onate, F. D. Pin, The particle finite element method: a powerful tool to solve incompressible flows with free-surfaces and breaking waves, IJNME 61 (2004) 964–989.
[8] M. Cremonesi et al., A state of the art review of the particle finite element method (pfem), Arch. Comp. Meth. Eng. 27 (2020) 1709–1735.
[9] M. Zhu, M. H. Scott, A pfem background mesh for simulating fluid and frame structure interaction, J. Str. Eng. 148 (2022) 04022051.
[10] O. Mavrouli et al., Damage analysis of masonry structures subjected to rockfalls, Landslides 14 (2017) 891–904.
[11] A. Montanino et al. Modelling with a meshfree approach the cornea-aqueous humor interaction during the air puff test. JMBBM, 77(2018), 205-216 .
[12] A. Montanino et al. Finite element formulation for compressible multiphase flows and its application to pyroclastic gravity currents. JCP, 451(2022), 110825.
[13] P.J. Ziółkowski et al. Adaptation of the arbitrary Lagrange–Euler approach to fluid–solid interaction on an example of high velocity flow over thin platelet. Continuum Mech. Therm., 33(6), pp. 2301-2314 2021
[14] P.J. Ziółkowski et al.. Fluid–solid interaction on a thin platelet with high-velocity flow: vibration modelling and experiment. Continuum Mech. Therm., pp. 1-27 (2022).
[15] E. Lorentz, S. Andrieux, Analysis of non-local models through energetic formulations, IJSS 40 (2003) 2905–2936.
[16] L. Ambrosio, V. M. Tortorelli, Approximation of functional depending on jumps by elliptic functional via t-convergence, Commun Pure Appl. Math 43 (1990) 999–1036.
[17] P. Sicsic, J.-J. Marigo, From gradient damage laws to Griffith’s theory of crack propagation, JElas 113 (2013) 55–74.
[18] J.-J. Marigo, Constitutive relations in plasticity, damage and fracture mechanics based on a work property, Nucl. Eng. Des. 114 (1989) 249–272.
[19] C. Comi, Computational modelling of gradient-enhanced damage in quasi-brittle materials, Mech. Cohes-fric Mat. 4 (1999) 17–36.
[20] K. Pham et al., Gradient damage models and their use to approximate brittle fracture, Int. Jour. Dam.Mech. 20 (2011) 618–652.
[21] C. Miehe et al., Phase field modeling of ductile fracture at finite strains: A variational gradient-extended plasticity-damage theory, Int. J.Plasticity 84 (2016) 1–32.
[22] F. dell’Isola et al., At the origins and in the vanguard of peridynamics, non-local and higher-gradient continuum mechanics: an underestimated and still topical contribution of Gabrio Piola, Math. Mech. Solids 20 (2015) 887–928.
[23] Z. P. Bazant, G. Pijaudier-Cabot, Nonlocal continuum damage, localization instability and convergence (1988).
[24] Z. P. Bazant, M. Jirasek, Nonlocal integral formulations of plasticity and damage: survey of progress, J. Eng. Mech. 128 (2002) 1119–1149.
[25] S. Forest, Micromorphic approach for gradient elasticity, viscoplasticity, and damage, J. Eng. Mech. 135 (2009) 117–131.
[26] R. Peerlings et al., A critical comparison of nonlocal and gradient-enhanced softening continua, IJSS 38 (2001) 7723–7746.
[27] U. Mühlich et al., A first-order strain gradient damage model for simulating quasi-brittle failure in porous elastic solids, Arch. Appl. Mech. 83 (2013) 955–967.
[28] L. Zybell et al., Constitutive equations for porous plane-strain gradient elasticity obtained by homogenization, Arch. Appl. Mech. 79 (2009) 359–375.
[29] C. Oliver-Leblond et al., A micro-mechanics based strain gradient damage model: formulation and solution for the torsion of a cylindrical bar, Eur. J. Mech.-A/Solids 56 (2016) 19–30.
[30] L. Placidi et al., Gedanken experiments for the determination of two-dimensional linear second gradient elasticity coefficients, ZAMP 66 (2015) 3699–3725.
[31] L. Placidi et al., Identification of two-dimensional pantographic structure via a linear D4 orthotropic second gradient elastic model, J.Eng.Math. 103 (2017) 1–21.
[32] Y. Rahali et al., Numerical identification of classical and nonclassical moduli of 3D woven textiles and analysis of scale effects, Comp. Str. 135 (2016) 122–139.
Risultati Pubblicati
I risultati principali sono stati pubblicati nell’articolo Diffusion of a damaging fluid through a dam-shaped bidimensional body for the estimation of its lifetime. Scrofani, A., Barchiesi, E., Chiaia, B., Misra, A., & Placidi, L. (2025).
La ricerca ha dimostrato con successo l’efficacia di un principio emivariazionale per simulare il decadimento meccanico indotto dalla diffusione di un fluido degradante. Sono state confrontati due modelli, uno rettangolare e uno trapezoidale, nella figura sottostante sono rappresentati il modello e le condizioni al contorno per le due forme di diga: rettangolare e trapezoidale.
Le analisi sono state condotte con il software COMSOL Multiphysics, il dominio delle due forme, con le diverse mesh è riportato in figura.
Il punto chiave dell’ approccio risiede nell’introduzione di un accoppiamento diretto tra la concentrazione del fluido e la variabile di danno . A differenza dei modelli classici, in questo framework il danno aumenta non solo per i carichi esterni, ma anche per la presenza chimica del fluido stesso; ciò permette di modellare il danneggiamento come una funzione monotonicamente crescente, compresa tra 0 e 1, che riduce progressivamente la rigidezza del materiale.
Per quanto riguarda l’influenza del coefficiente di diffusione sulla vita utile, l’analisi parametrica ha rivelato un comportamento non lineare e altamente realistico. In condizioni di bassa diffusività, il fluido resta confinato vicino alla superficie di contatto, portando la concentrazione a valori elevatissimi in un’area ristretta e causando un’evoluzione lenta del danno totale. Al contrario, in presenza di alta diffusività, il fluido penetra più facilmente distribuendosi su un’area maggiore; questo riduce il tasso di crescita della concentrazione locale e, paradossalmente, può rallentare l’evoluzione del danno globale dopo aver superato un valore minimo di vita utile. Il grafico della longevità in funzione di mostra quindi un punto di minimo, oltre il quale la permeabilità aiuta a distribuire l’agente degradante, rallentando il collasso strutturale. Inoltre, il parametro di accoppiamento concentrazione-danno si è rivelato critico, poiché anche piccole variazioni nel suo valore, ad esempio nel range da , producono cambiamenti significativi nella vita utile. Un modulo maggiore di questo fattore accelera drasticamente l’insorgenza del danno a parità di carico, evidenziando come la composizione chimica dell’interazione fluido-solido sia determinante per la durabilità.
Infine, il confronto tra la forma rettangolare e quella trapezoidale ha confermato scientificamente le intuizioni dell’ingegneria tradizionale. La forma trapezoidale garantisce che quasi l’intera struttura rimanga in compressione, mentre il rettangolo presenta ampie zone in trazione alla base, dove la traccia del tensore di deformazione risulta Poiché la soglia di attivazione del danno è dieci volte inferiore in trazione rispetto alla compressione, il modello rettangolare degrada molto più velocemente. Di conseguenza, la diga trapezoidale mostra una vita utile superiore del 15% rispetto a quella rettangolare, con un inizio dell’evoluzione del danno che avviene solo dopo oltre 50 anni, contro l’insorgenza quasi immediata rilevata nel caso rettangolare.
2) La ricerca ha portato alla luce i seguenti risultati:
1. Efficacia del Modello Emivariazionale e Accoppiamento
La ricerca ha dimostrato con successo l’efficacia di un principio emivariazionale per simulare il decadimento meccanico indotto dalla diffusione di un fluido degradante. Il punto chiave è l’introduzione di un accoppiamento diretto tra la concentrazione del fluido e la variabile di danno (𝜔). A differenza dei modelli classici, qui il danno aumenta non solo per i carichi esterni, ma anche per la presenza chimica del fluido, permettendo di modellare il danno come una funzione monotonicamente crescente (da 0 a 1) che riduce progressivamente la rigidezza del materiale.
2. Influenza del Coefficiente di Diffusione (𝐾𝐷𝐼𝐹) sulla Vita Utile
L’analisi parametrica ha rivelato un comportamento non lineare e altamente realistico della vita utile della struttura in relazione alla diffusività:
Bassa diffusività: Il fluido resta confinato vicino alla superficie di contatto. La concentrazione raggiunge valori elevatissimi in un’area ristretta, portando a un’evoluzione lenta del danno totale.
Alta diffusività: Il fluido penetra più facilmente, distribuendosi su un’area maggiore. Questo riduce il tasso di crescita della concentrazione locale e, paradossalmente, può rallentare l’evoluzione del danno globale dopo aver superato un valore minimo di vita utile.
Il grafico della vita utile in funzione di 𝐾𝐷𝐼𝐹 mostra quindi un punto di minimo, oltre il quale la permeabilità aiuta a “spalmare” l’agente degradante, rallentando il collasso.
3. Sensibilità al Fattore di Accoppiamento Concentrazione-Danno (𝐾𝑐𝜔)
Il parametro 𝐾𝑐𝜔 si è rivelato critico: anche piccole variazioni nel suo valore (da -4.5 a -5.5 m^2⁄s^2 ) producono cambiamenti significativi nella vita utile. Un modulo maggiore di questo fattore accelera drasticamente l’insorgenza del danno a parità di carico, evidenziando come la composizione chimica dell’interazione fluido-solido sia determinante per la durabilità.
4. Superiorità Geometrica della Diga Trapezoidale
Il confronto tra la forma rettangolare e quella trapezoidale (a parità di area e altezza) ha confermato scientificamente le intuizioni dell’ingegneria tradizionale:
Stato Tensionale: La forma trapezoidale garantisce che quasi l’intera struttura rimanga in compressione, mentre il rettangolo presenta ampie zone in trazione alla base (traccia del tensore di deformazione (𝐺>0).
Soglie di Danno: Poiché la soglia di attivazione del danno è 10 volte inferiore in trazione rispetto alla compressione, il modello rettangolare degrada molto più velocemente.
Vita Utile: La diga trapezoidale mostra una vita utile superiore del 15% rispetto a quella rettangolare, con l’inizio dell’evoluzione del danno che avviene dopo oltre 50 anni, contro l’insorgenza quasi immediata nel caso rettangolare.
Il sito è in costante aggiornamento: nuovi risultati, sviluppi e materiali saranno pubblicati qui man mano che la ricerca avanza