{"id":24423,"date":"2026-01-24T04:11:32","date_gmt":"2026-01-24T03:11:32","guid":{"rendered":"https:\/\/test.cfdfeaservice.it\/index.php\/2026\/01\/24\/la-modellazione-dei-contatti-nelle-analisi-agli-elementi-finiti\/"},"modified":"2026-01-24T04:11:32","modified_gmt":"2026-01-24T03:11:32","slug":"la-modellazione-dei-contatti-nelle-analisi-agli-elementi-finiti","status":"publish","type":"post","link":"https:\/\/test.cfdfeaservice.it\/index.php\/2026\/01\/24\/la-modellazione-dei-contatti-nelle-analisi-agli-elementi-finiti\/","title":{"rendered":"La modellazione dei contatti nelle analisi agli elementi finiti"},"content":{"rendered":"
\n
\"La<\/div>\n
\n

Un cosiddetto \u201cproblema di contatto\u201d \u00e8 classificabile come una non linearit\u00e0 al contorno, in cui sia lo spostamento, sia la forza di contatto sono incognite del problema.<\/strong><\/p>\n

A differenza di altri tipi di simulazione in cui alcuni nodi della mesh del componente sono sottoposti ad un carico noto, nella simulazione di contatto sia la posizione del contatto (e quindi il set di nodi), sia il valore del carico sono sconosciute.<\/p>\n

Con un\u2019analisi dei contatti \u00e8 possibile determinare:<\/p>\n

\u2022\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 Se due o pi\u00f9 corpi sono in contatto<\/p>\n

\u2022\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 La regione di contatto<\/p>\n

\u2022\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 La forza di contatto o pressione all\u2019interfaccia<\/p>\n

\u2022\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 Se all\u2019interfaccia \u00e8 presente un movimento relativo dopo il contatto<\/p>\n

La procedura FEA per un problema di contatto prevede due fasi principali:<\/p>\n

\u2022\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 Valutazione della sussistenza di un contatto<\/p>\n

\u2022\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 Nel caso un contatto sia presente, determinazione della forza di contatto corrispondente<\/p>\n

Trovare l\u2019equilibrio del sistema implica trovare un campo di spostamento che minimizzi l\u2019energia potenziale. Questa minimizzazione \u00e8 limitata dalle condizioni di contatto, il che significa che la soluzione deve soddisfare ulteriori restrizioni geometriche e\/o meccaniche.<\/p>\n

Numericamente, questo problema pu\u00f2 essere risolto utilizzando il metodo della penalizzazione, che approssima i vincoli attraverso una rigidezza fittizia, oppure mediante il metodo del moltiplicatore di Lagrange, che introduce variabili ausiliarie per imporre i vincoli in modo esatto.<\/p>\n

I problemi di contatto sono onerosi dal punto di vista computazionale per varie ragioni.<\/p>\n

La prima \u00e8 legata al fatto che la reale area di contatto \u00e8 sconosciuta: la vera zona di contatto non \u00e8 nota a priori, ma deve essere determinata come parte stessa della soluzione. Sebbene spesso sia possibile identificare una zona massima in cui il contatto avverr\u00e0 e, dunque, limitare la ricerca a quella zona, la accuratezza con cui l\u2019area viene determinata dipende dalla convergenza iterativa.<\/p>\n

La seconda causa \u00e8 legata ad una risposta discontinua della forza: le forze di contatto, infatti, mostrano cambiamenti bruschi (ad esempio, da zero ad un valore finito quando le superfici si toccano). Questa discontinuit\u00e0 assomiglia ai vincoli di incomprimibilit\u00e0, che richiedono metodi specializzati (moltiplicatori di Lagrange o metodi di penalit\u00e0)<\/p>\n

La terza va individuata nella sensibilit\u00e0 discreto al contorno: nei modelli numerici, l\u2019area di contatto evolve da nodo a nodo e non in modo uniforme come nella teoria del continuo. Ci\u00f2 rende i risultati altamente sensibili alla discretizzazione (dimensione della mesh).<\/p>\n

\"\"<\/a>
Figura 1: Illustrazione dell\u2019evoluzione dei confini di contatto.<\/em><\/figcaption><\/figure>\n

1. Approcci di risoluzione<\/strong><\/h2>\n

Si consideri una trave a sbalzo sottoposta ad un carico distribuito che entri in contatto con una superficie rigida. Il sistema \u00e8 caratterizzato da un carico distribuito q = 1 kN\/m agente sulla trave di lunghezza L = 1m con rigidezza flessionale EI = 105Nm2<\/sup>. La deflessione massima consentita \u00e8 pari a \u03b4=1mm.<\/p>\n

Ci sono tre approcci principali per risolvere il problema: soluzione per tentativi (trial and error), moltiplicatore di Lagrange o metodo della penalizzazione.<\/p>\n

1.1. Soluzione per tentativi<\/em><\/strong><\/h3>\n

La deflessione della punta di una trave a sbalzo sottoposta ad un carico distribuito \u00e8 data da:<\/p>\n

\"\"<\/a><\/figure>\n
\"\"<\/a><\/figure>\n

Quando la deflessione della punta VN<\/sub>(L) supera la distanza consentita \u03b4=1mm, si verifica il contatto con la superficie rigida. Ci\u00f2 richiede l\u2019introduzione di una forza di reazione verso l\u2019alto \u03bb agente sulla punta della trave in modo da impedire la compenetrazione. Il contributo alla deflessione di questa forza di reazione \u00e8 descritto da:<\/p>\n

\"\"<\/a><\/figure>\n
\"\"<\/a><\/figure>\n

Di conseguenza, la forza \u03bb necessaria ad impedire la compenetrazione pu\u00f2 essere calcolata imponendo impostando Vtip<\/sub>=\u03b4.<\/p>\n

\"\"<\/a><\/figure>\n

Da cui \u03bb = 75N.<\/p>\n

\u00c8 possibile trattare sia la forza di contatto, sia la distanza come sconosciute<\/p>\n

\"\"<\/a><\/figure>\n

e aggiungere vincoli. In particolare, sono tre gli stati distinti che caratterizzano il comportamento del sistema. Innanzitutto, quando VN<\/sub>(L)<\u03b4, non si verifica alcun contatto per cui \u03bb=0. Il contatto attivo (\u03bb>0) richiede la condizione Vtip<\/sub>=\u03b4. La terza condizione \u00e8 rappresentata dal caso fisicamente inammissibile di \u03bb<0 che implicherebbe una reazione di trazione al vincolo esterno. La funzione gioco g=Vtip<\/sub>\u2212\u03b4 deve soddisfare g\u22640 per impedire la penetrazione, con valori positivi di g che rappresentano valori non fisici .<\/p>\n

1.2. Moltiplicatore di Lagrange<\/em><\/strong><\/h3>\n

Per l\u2019esempio precedente, la condizione di contatto prevedeva:<\/p>\n

\u2022\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 Nessuna penetrazione (g \u2264 0)<\/p>\n

\u2022\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 Forza di contatto positiva (\u03bb \u2265 0)<\/p>\n

\u2022\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 Condizione di consistenza (\u03bbg = 0)<\/p>\n

Il moltiplicatore di Lagrange pu\u00f2 essere descritto come una strategia matematica per trovare i massimi e i minimi locali di una funzione soggetta a vincoli. \u03bbg risulta infatti<\/p>\n

\"\"<\/a><\/figure>\n

Quando \u03bb=0N \u2192 g=0.00025>0 \u2192 viene violata la condizione di contatto.<\/p>\n

Quando \u03bb=75N \u2192 g = 0 \u2192 la condizione di contatto risulta soddisfatta.<\/p>\n

1. 3. Metodo di penalizzazione<\/em><\/strong><\/h3>\n

Nel metodo della penalizzazione, \u00e8 ammessa una piccola penetrazione tra le parti. La forza di contatto risulta proporzionale ad essa. La funzione di penetrazione \u00e8 definita come:<\/p>\n

\"\"<\/a><\/figure>\n

che risulta pari a 0 quando g\u22640 (nessun contatto) e pari a g quando g>0 (contatto). La forza di contatto \u00e8 definita come<\/p>\n

\"\"<\/a><\/figure>\n

dove KN<\/sub> \u00e8 il parametro di penalizzazione (una sorta di rigidezza).<\/p>\n

Da \"Titolo:\u00a0si ha<\/p>\n

\"\"<\/a><\/figure>\n

Sembra dunque che il divario dipenda dal parametro della penalit\u00e0.<\/p>\n

Tabella 1: risultati in funzione del parametro di penalizzazione scelto<\/em><\/p>\n

\n\n\n\n\n\n\n
KN<\/sub><\/td>\nPenetrazione (-g) [m]<\/td>\nForza di contatto (\u03bb) [N]<\/td>\n<\/tr>\n
\"Titolo:<\/td>\n\"Titolo:<\/td>\n\"Titolo:<\/td>\n<\/tr>\n
\"Titolo:<\/td>\n\"Titolo:<\/td>\n\"Titolo:<\/td>\n<\/tr>\n
\"Titolo:<\/td>\n\"Titolo:<\/td>\n\"Titolo:<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<\/figure>\n

I risultati numerici mostrano il comportamento caratteristico del contatto (forza e compenetrazione) al variare del parametro di penalizzazione: all\u2019aumentare del parametro di penalizzazione KN<\/sub>, la penetrazione fisica \u2212g diminuisce, seguendo la relazione inversa \u2212g\u2248 \u03bb\/KN<\/sub>, mentre la forza di contatto \u03bb converge verso la soluzione esatta. Questo illustra la necessita di trovare un compromesso ideale per l\u2019impiego dei metodi di penalizzazione: valori elevati di KN<\/sub> applicano la regola di \u201cnon penetrazione\u201d in modo pi\u00f9 rigoroso fornendo dunque soluzioni accurate. Tuttavia, valori eccessivi rendono i calcoli pi\u00f9 complessi. Al contrario, valori eccessivamente piccoli mantengono i calcoli pi\u00f9 semplici ma consentono una maggiore penetrazione (effetto non fisico).<\/p>\n

\"\"<\/a>
Figura 2: Relazione discontinua forza-spostamento.<\/em><\/figcaption><\/figure>\n

2. Attrito<\/strong><\/h2>\n

In alcuni casi, \u00e8 anche interessante capire se esista un movimento relativo tra i corpi a contatto. Per questo motivo, dopo aver applicato una forza q, la stessa trave a sbalzo dell\u2019esempio precedente viene caricata assialmente anche da una forza concentrata P=100N applicata all\u2019estremit\u00e0 di destra. Si presume che il coefficiente di attrito sia \u03bc=0.5.<\/p>\n

2.1 Soluzione per tentativi ed errori<\/em><\/strong><\/h3>\n

In caso di assenza di attrito, l\u2019applicazione della forza P induce uno spostamento assiale dell\u2019apice della trave.<\/p>\n

\"\"<\/a><\/figure>\n

Il vincolo ulteriore dato dall\u2019attrito pu\u00f2 essere espresso come:<\/p>\n

\u2022\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 Condizione di adesione se t\u2212\u03bc\u03bb<0, Utip<\/sub>=0<\/p>\n

\u2022\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 Condizione di scorrimento se t\u2212\u03bc\u03bb=0, Utip<\/sub>>0<\/p>\n

\u2022\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 Condizione di consistenza Utip<\/sub>(t-\u03bc\u03bb)=0<\/p>\n

in cui t \u00e8 la forza di attrito tangenziale.<\/p>\n

Per prima cosa valutiamo la condizione di adesione<\/p>\n

\"\"<\/a><\/figure>\n

Avere t=100N viola la condizione di adesione t\u2212\u03bc\u03bb<0. Da questo \u00e8 possibile dedurre che vi sar\u00e0 dello scorrimento.<\/p>\n

Sapendo che \u03bb=75N, t potrebbe essere calcolato come t=\u03bc\u03bb=37.5N.<\/p>\n

\"\"<\/a><\/figure>\n

Questa soluzione soddisfa la condizione di scorrimento Utip<\/sub>>0 e, pertanto, \u00e8 una soluzione valida.<\/p>\n

2.2 Moltiplicatore di Lagrange<\/em><\/strong><\/h3>\n

In questo caso viene utilizzata la condizione di consistenza Utip<\/sub>(t\u2212\u03bc\u03bb)=0. In caso di scorrimento,<\/p>\n

Utip<\/sub> pu\u00f2 essere scelto come moltiplicatore Lagrangiano e (t\u2212\u03bc\u03bb) come vincolo per cui<\/p>\n

\"\"<\/a><\/figure>\n

Quando Utip<\/sub>=0, t=P e t\u2212\u03bc\u03bb=62.5>0. Ci\u00f2 viola la condizione di adesione. Quando Utip<\/sub>=(P\u2212\u03bc\u03bb)L=0.625 mm>0, t=\u03bc\u03bb la condizione di scorrimento \u00e8 soddisfatta (soluzione valida).<\/p>\n

2.3 Metodo di penalizzazione<\/em><\/strong><\/h3>\n

Si penalizza quando t\u2212\u03bc\u03bb>0.<\/p>\n

\"\"<\/a><\/figure>\n

Lo scorrimento (traslazione assiale relativa) viene calcolato come<\/p>\n

\"\"<\/a><\/figure>\n

con KT<\/sub> parametro di penalizzazione per lo scorrimento tangenziale.<\/p>\n

Quando t\u2212\u03bc\u2212\u03bb<0 \u2192 adesione, Utip<\/sub>=0 (nessuna penalizzazione). Quando t\u2212\u03bc\u2212\u03bb>0 \u2192 scorrimento, si penalizza per rimanere vicino a t=\u03bc\u03bb . In caso di scorrimento, la forza di attrito pu\u00f2 essere derivata da<\/p>\n

\"\"<\/a><\/figure>\n
\"\"<\/a><\/figure>\n

Lo spostamento della punta risultata<\/p>\n

\"\"<\/a><\/figure>\n

Tabella 1: risultati in funzione del parametro di penalizzazione scelto<\/em><\/p>\n

\n\n\n\n\n\n\n
KT<\/sub><\/td>\nSpostamento della punta (Utip<\/sub>) [m]<\/td>\nForza tangenaziale (t) [N]<\/td>\n<\/tr>\n
\"Titolo:<\/td>\n\"Titolo:<\/td>\n\"Titolo:<\/td>\n<\/tr>\n
\"Titolo:<\/td>\n\"Titolo:<\/td>\n\"Titolo:<\/td>\n<\/tr>\n
\"Titolo:<\/td>\n\"Titolo:<\/td>\n\"Titolo:<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<\/figure>\n
\"\"<\/a>
Figura 3: Schema di carico dell\u2019esempio considerato.<\/em><\/figcaption><\/figure>\n

3. Formulazione generale per problemi di contatto<\/strong><\/h2>\n

Gli esempi di contatto puntuali descritti nella sezione precedente sono limitati da un punto di vista pratico, poich\u00e9 la maggior parte dei contatti nelle applicazioni ingegneristiche prevede un contatto distribuito lungo curve o superfici. Nella figura, \u0393C<\/sub> denota la potenziale zona di contatto sul corpo deformabile (slave) che pu\u00f2 interagire con una superficie rigida (master) a seguito dell\u2019applicazione di un carico.<\/p>\n

L\u2019analisi del contatto diventa significativamente pi\u00f9 impegnativa quando si passa da punti discreti a superfici continue. Innanzitutto, l\u2019area di contatto effettiva si evolve durante il carico e deve essere determinata come parte della soluzione. L\u2019utente deve specificare un limite di contatto \u0393C<\/sub> (potenziale area di contatto) sufficientemente grande da comprendere tutte le possibili zone di interazione. Di conseguenza, i problemi di contatto richiedono la determinazione simultanea delle posizioni di contatto precise xc<\/sub>(\u03bec<\/sub>) sulla superficie master e delle corrispondenti forze di contatto (sia pressioni di contatto che eventuali forze di attrito). Infine, le condizioni di contatto introducono vincoli geometrici non lineari che in genere richiedono strategie di soluzione iterative.<\/p>\n

Per ogni punto candidato x \u2208 \u0393C<\/sub>, si deve trovare la sua proiezione ortogonale<\/p>\n

xc<\/sub>(\u03bec<\/sub>) sulla superficie master. Ci\u00f2 si riduce alla risoluzione dell\u2019equazione:<\/p>\n

\"\"<\/a><\/figure>\n

dove il vettore tangente \"Titolo:\u00a0caratterizza l\u2019origine della superficie locale a \u03beC<\/sub>. La funzione gioco normale:<\/p>\n

\"\"<\/a><\/figure>\n

permette di determinare se si verifica un contatto (gn<\/sub>\u22640) o se rimane del gioco (gn<\/sub>>0). Il versore normale \"Titolo:\u00a0completa la descrizione della cinematica dei contatti.<\/p>\n

Quando il contatto \u00e8 stabilito, il moto tangenziale relativo viene quantificato attraverso la funzione di scorrimento:<\/p>\n

\"\"<\/a><\/figure>\n

tracciando il movimento lungo la superficie master della configurazione precedente.<\/p>\n

ESEMPIO: Superficie di contatto parabolica<\/em><\/strong><\/p>\n

Si consideri la proiezione di un punto x=[3, 1]T<\/sup> sulla superficie rigida parabolica y=x2<\/sup><\/p>\n

Step di soluzione:<\/p>\n

1.\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 Parametrizzazione della superficie<\/strong>: la superficie master \u00e8 descritta parametricamente come xc<\/sub>(\u03be)=[\u03be,\u03be2]T<\/sup>, con vettore tangente xc,\u03be<\/sub>=[1,2\u03be]T<\/sup> .<\/p>\n

2.\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 Calcolo della proiezione<\/strong>: risoluzione della condizione di coerenza del contatto da cui \u00a0\u00a0 l\u2019equazione cubica:<\/p>\n

\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 3+\ud835\udf09\u22122\ud835\udf093<\/sup>=0<\/p>\n

\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 La radice fisicamente significativa \u03bec<\/sub>\u22481.29 che individua il punto di contatto in xc<\/sub>= \u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 [1.29,1.66]T<\/sup>.<\/p>\n

3.\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 Valutazione del gioco<\/strong>: il calcolo del gioco normale permette di discernere tra \u00a0\u00a0\u00a0 condizione di penetrazione e non, mentre la funzione di scorrimento traccerebbe qualsiasi movimento tangenziale successivo.<\/p>\n

\"\"<\/a>
Figura 4: Schema di carico dell\u2019esempio considerato (attrito).<\/em><\/figcaption><\/figure>\n

Conclusioni<\/strong><\/h2>\n

Nelle simulazioni numeriche agli elementi finiti, la risoluzione dei contatti \u00e8 una delle sfide maggiori. Questa viene affrontata con molteplici approcci. Il pi\u00f9 comune \u00e8 quello della penalizzazione che per\u00f2 richiede una accurata scelta dei parametri per garantire risultati accurati ed un onere computazionale ragionevole.<\/p>\n

\"\"<\/a>
Figura 5: Scenario di contatto generalizzato tra superfici deformabili (slave) e rigide (master).<\/em><\/figcaption><\/figure>\n
\"\"<\/a>
Figura 6: Proiezione ortogonale del punto sulla superficie master parabolica.<\/em><\/figcaption><\/figure>\n<\/div>\n

L’articolo La modellazione dei contatti nelle analisi agli elementi finiti<\/a> sembra essere il primo su Il Progettista Industriale<\/a>.<\/p>\n<\/div>\n

Vai alla fonte.<\/a><\/p>\n

Autore: Emanuela Bianchi<\/p>\n

Powered by WPeMatico<\/a><\/small><\/p>\n

_________________________________<\/strong><\/p>\n

CFD FEA Service SRL<\/strong> è una società di servizi che offre consulenza<\/em> e formazione<\/em> in ambito ingegneria<\/strong> e IT<\/strong>. Se questo post\/prodotto ti è piaciuto ti invitiamo a:<\/p>\n