Comi, Claudia

M., Negri

Mariani, Stefano

Perego, Umberto

COMI, CLAUDIA

MARIANI, STEFANO

M. Negri

PEREGO, UMBERTO

Archivio istituzionale della ricerca - Politecnico di Milano

 1 UNA FORMULAZIONE VARIAZIONALE PER LA FRATTURA COESIVA IN UNA BARRA IN TRAZIONE   C. COMI1, S. MARIANI1, M. NEGRI2, U. PEREGO1  1 Dipartimento di Ingegneria Strutturale, Politecnico di Milano, Milano 2 Dipartimento di Matematica, Università di Pavia, Pavia    SOMMARIO In questo lavoro si studia il problema della formazione e della crescita di una frattura coesiva in una barra elastica. Si mostra come sia possibile ottenere la posizione di attivazione della cricca dalle condizioni di minimo di un funzionale di energia che contenga anche l’energia superficiale della cricca. Il danneggiamento progressivo dell’interfaccia coesiva è tenuto in conto grazie ad una procedura passo-passo. Si mostra che anche la soluzione del problema nel passo finito rende stazionario un funzionale definito ad ogni passo.  ABSTRACT In this work, the problem of inception and growth of a cohesive crack in an elastic bar is considered. The position where the crack actually forms is obtained from the minimality conditions of an energy functional which includes the surface energy. The progressive damage of the cohesive interface is taken into account by means of a step by step procedure. The finite step solution is also shown to make stationary a functional defined for each step.    1. INTRODUZIONE  Si considera una barra soggetta ad un carico distribuito assiale ed a spostamento imposto ad un estremo. Il materiale è assunto elastico lineare con un limite elastico variabile lungo la barra. Quando tale limite viene raggiunto si assume che si possa attivare una cricca coesiva con un legame softening tra tensione trasmessa e discontinuità di spostamento, fino a rottura completa. Nella letteratura recente sono presenti numerosi contributi intesi a fornire una caratterizzazione variazionale del problema della generazione di fratture in mezzi elastici [1-6]. Facendo riferimento all’impostazione di Braides et alii [2], inizialmente si assume per il sistema un comportamento perfettamente reversibile (il materiale ritorna integro a seguito dell’eventuale richiusura della cricca) e si dimostra che la soluzione del problema minimizza l’energia potenziale totale del sistema che comprende un termine di energia elastica, definito sul volume, ed un termine di superficie relativo alla cricca. La posizione di innesco della GIMC‘04XV Congresso Italiano di Meccanica Computazionale AIMETA Associazione Italiana di Meccanica Teorica e Applicata 2 frattura può anch’essa essere ottenuta dalle condizioni di ottimalità per valore nullo dell’apertura di fessura. Approssimando l’evoluzione reale con una sequenza di passi finiti reversibili, ed imponendo quindi la non reversibilità solo al termine del passo, si dimostra poi che la soluzione di ogni passo finito può essere ottenuta dalle condizioni di ottimalità di un opportuno funzionale. Un semplice esempio con carico assiale distribuito in modo uniforme e spostamento imposto ad un estremo crescente nel tempo illustra i risultati teorici conseguiti.   2. FRATTURA COESIVA DI UNA BARRA TESA  Si consideri una barra di lunghezza L=10 mm, di modulo elastico pari a E=1 MPa, uniforme su tutta la sua lunghezza, soggetta ad un carico distribuito uniforme in direzione del suo asse pari a b=0.2 E/L. Sia inoltre imposto al suo estremo uno spostamento ( )η t , monotonamente crescente nel tempo. Si assuma che la resistenza alla frattura della barra sia data da ( ) ( ) ( )0 0=p x a x g , variabile lungo la barra secondo la legge ( )2( ) 1 2= + −a x x L , con ( )0 2 0.1=p L  E. Il modello coesivo assunto prevede una diminuzione lineare della resistenza con l’apertura w della frattura (figura 1). Alla fine del passo che va dall’istante nt  all’istante 1+nt , la forza per unità di superficie 1+np  scambiata tra le facce di una cricca posta in z è governata dalle seguenti relazioni   ( ) ( ) ( )( ) ( )( ) ( ) ( ) ( )( )11 1 1 1 110 11 1 1 1,    for 0,0, [0, ]               for 00   0    0ξξ ξξξξ ξ ξ ξ++ + + + ++++ + + + = , = > 0 ∈ =− ≤ − ∆ = ∆ ≥nn n n n nnnn n n ngp p w z a z wp z p zp a z g p a z g  dove ξn  rappresenta il danno non reversibile accumulato sino all’inizio del passo corrente. E’ possibile definire un potenziale ( ),ξ% nG w  per l’energia superficiale della cricca nel passo    ( )( )( ) ( ) ( )2   for 2,   for w2nnnnn nn ngw wG wgG w Gξ ξξξ ξ ξξ ξ ≤=  − + ≥%  Nel caso considerato, in cui si ha una decrescita lineare della forza al crescere dell’apertura, la funzione ( ) = ∂ ∂%g w G w  per 0ξ =n  è data da (figura 1)  ( ) 011     for 0                        for w  − ≤  =    ≥critcritcritg w w wg w ww con 0.15=critw L  e 0 0.1=g  E.  Definendo il seguente funzionale energia per il passo corrente (u rappresenta lo spostamento in direzione assiale)  3           Figura 1– Legame tra tensione e discontinuità di spostamento.    ( ) ( ) ( ) ( )( )21, ,2η ξ ξ∈ = − +   ∑∫ ∫ %un nI I z SduU u E dx budx a z G w z zdx  dove uS  è il luogo delle possibili sedi di frattura lungo la barra, è possibile mostrare come la posizione di prima attivazione della cricca sia ottenibile minimizzando il funzionale rispetto alla posizione della frattura per 0ξ =  e come tutte le altre equazioni che governano il problema siano ottenibili come condizioni di stazionarietà del medesimo funzionale. Lo spostamento u può essere eliminato dal funzionale imponendo il rispetto dell’equilibrio, ottenendo così un nuovo funzionale U  che dipende solo dallo spostamento impresso η , da w e dalla posizione della cricca x :  ( ) ( ) ( )( ) ( )wGLxELbwLxwbLwExw   −++−+−+−=23222124222,, ηηηU  Imponendo la stazionarietà rispetto a w e x  ed imponendo w = 0, è possibile ricavare il valore 0 0.099η = L  dello spostamento impresso che porta alla prima attivazione della cricca e la posizione 0.49=x L  della frattura. La figura 2a mostra l’andamento del funzionale ( )xw,,0ηU  in corrispondenza della prima attivazione della cricca. E’ evidente l’avvallamento in corrispondenza del valore 0.49=x L  che minimizza il funzionale. Dal punto di vista meccanico, la posizione di prima attivazione della cricca corrisponde al punto x  in cui per la prima volta lo sforzo raggiunge il valore limite ( ) ( ) ( )0 0=p x a x g , come mostrato in figura 2b. Al crescere dello spostamento impresso, è sempre energeticamente più conveniente aprire la fessura coesiva già formatasi piuttosto che deformare ulteriormente la barra. Questo si evince chiaramente dal grafico in figura 3a, dove è mostrato l’andamento del funzionale ( )xw,,ηU in soluzione per valori crescenti di η. Nella figura 3a, la curva rossa rappresenta l’energia della barra integra, mentre la curva blu, l’energia della barra fessurata. E’ opportuno notare che per valori di 0η η> , la soluzione si può ottenere dalle condizioni di stazionarietà di U  solo per x  fissato (punto bianco in figura 3b). Il minimo assoluto di U prevederebbe la presenza di una cricca in una posizione diversa da 0.49=x L  e non rappresenta quindi una soluzione ammissibile del problema.  Ringraziamenti – Questo lavoro è stato svolto nell’ambito del progetto MIUR- PRIN03 su  Interfacial damage failure in structural systems: applications to civil engineering and emerging research fields.   ξnw wcritp p0 ( )ξξ nngapn  4 0246810 00.511.5-0.100.10.2x w U               (a)       (b) Figura 2 – (a) energia per 0ηη = ; (b) posizione della frattura.                Figura 3 – (a) energia in funzione dello spostamento imposto; (b) curve di livello dell’energia per L14.0=η .  BIBLIOGRAFIA  [1] G. Francfort, J.J. Marigo, Revisiting brittle fracture as an energy minimization problem, J. Mech. Phys. Solids, vol. 46, pp. 1319–1342, (1998). [2] A. Braides, G. Dal Maso, A. Garroni, Variational formulation of softening phenomena in fracture mechanics: The one-dimensional case, Arch. Ration. Mech. Anal., vol. 146, pp. 23–58, (1999).  [3] G. Maier, C. Comi, Energy properties of solutions to quasi-brittle fracture mechanics problems with piecewise linear cohesive crack models, in: Continuous damage and fracture, A. Benallal ed., Elsevier (2000). [4] G. Del Piero, L. Truskinovsky, Macro- and micro-cracking in one-dimensional elasticity, Internat. J. Solids Structures, vol. 38, pp. 1135–1148, (2001). [5] M. Negri, A finite element approximation of the Griffith’s model in fracture mechanics, Numer. Math., vol. 95, pp. 653–687, (2003).  [6] M. Angelillo, E. Babilio, A. Fortunato, A computational approach to fracture of brittle solids based on energy minimization, Preprint Università di Salerno, (2004). 4.6 4.7 4.8 4.9 5 5.1 5.211.11.21.31.41.5w x η U 0.25 0.5 0.75 1 1.25 1.5-0.1-0.08-0.06-0.04-0.022 4 6 8 100.20.40.60.811.2x p0(x) σ (x) 

Formulazione variazionale per la frattura coesiva in una barra in trazione

https://re.public.polimi.it/bitstream/11311/273320/1/2004_Comi_Mariani_Negri_Perego_Gimc_XV_Genova.pdf

Formulazione variazionale per la frattura coesiva in una barra in trazione

Abstract

Similar works

Full text

Available Versions

Archivio istituzionale della ricerca - Politecnico di Milano