Come prevedere la deformazione dei pezzi sottili in pochi secondi?

Quando una fresa taglia un pezzo di metallo a parete sottile, il lavoro può diventare una sfida: la parte si deforma e perde precisione. Tradizionalmente, gli ingegneri usano il metodo degli elementi finiti (FEM) per simulare queste deformazioni, ma il calcolo richiede ore, troppo per decisioni rapide in fabbrica.

Un team di ricercatori ha ideato un’alternativa più veloce, basata sul modello voxel octree. Immagina di suddividere il pezzo in piccoli cubetti (voxel) organizzati in una struttura ad albero a otto rami. Ogni livello dell’albero rappresenta una risoluzione diversa: i voxel più grandi descrivono la forma generale, quelli più piccoli catturano i dettagli critici.

Gli scienziati hanno calcolato la matrice di rigidità per ogni voxel e hanno inserito le forze di taglio microscopiche generate dalla fresa. Con questi dati, il software stima la deformazione in tempo reale. I test mostrano che le previsioni coincidono quasi perfettamente con le misurazioni reali, ma il tempo di calcolo è ridotto del 90% rispetto al FEM e del 13% rispetto al modello voxel tradizionale.

Questa velocità permette agli operatori di ottimizzare subito i parametri di taglio, evitando errori costosi e migliorando la qualità dei prodotti finiti. In pratica, è come avere una sfera di cristallo digitale che indica subito come il pezzo reagirà al taglio.

Il modello voxel octree sfrutta la gerarchia spaziale per ridurre drasticamente il numero di gradi di libertà rispetto al FEM tradizionale. In pratica, il dominio del pezzo viene suddiviso in sei livelli di voxel, dove il livello più alto contiene pochi voxel grossolani e il livello più basso ne contiene migliaia di dimensioni micrometriche. Per ciascun voxel si definisce una matrice di rigidità locale K_i, derivata dalle proprietà elastiche del materiale (modulo di Young, coefficiente di Poisson) e dalla geometria del voxel stesso.

Le forze di taglio micro‑milling F_m sono calcolate mediante modelli di contatto tribologico, tenendo conto della velocità di taglio, dell’angolo di attacco e della rugosità dell’utensile. La matrice globale di rigidità K_tot si ottiene assemblando le K_i secondo le regole di continuità meccanica, ma grazie alla struttura ad albero è possibile eseguire l’assemblaggio in modo ricorsivo, aggiornando solo i rami interessati dal contatto.

Il problema di equilibrio K_tot·u = F_m viene risolto con metodi iterativi precondizionati, che convergono rapidamente grazie alla sparsità della matrice. Il risultato è il vettore spostamento u, che fornisce la deformazione locale in ogni voxel. Confrontando u con le misurazioni laser‑tracker, si osserva una correlazione R^2 > 0,98, confermando l’accuratezza del modello. Inoltre, la riduzione del tempo computazionale (da ore a pochi minuti) rende il metodo adatto a simulazioni in fase di progettazione, consentendo iterazioni rapide per l’ottimizzazione dei parametri di taglio, come velocità di avanzamento, profondità di passata e geometria dell’utensile.