Paragrafo 1.7 tensore di deformazione U e direzioni principali (trasformazione degli elementi di linea)
L'individuazione delle direzioni principali del tensore di deformazione consente di ottenenere una rappresentazione “intrinseca” della deformazione a cui è soggetto un generico punto del corpo. L'aggettivazione “intrenseca” si giustifica perché, come vedremo, si ottiene una rappresentazione più semplice dello stato deformativo ma, soprattutto, indipendente dal sistema di riferimento utilizzato. La discussione che segue, in modo graduale, presenta l'idea che sta dietro alle direzioni principali facendo riferimento al tensore destro della deformazione, U, in una forma particolare: la forma diagonale. Con le dovute precisazioni, quanto si ottiene vale anche per il tensore sinistro della deformazione, V.Sottoparagrafo 1.7.1
Si supponga per U la seguente forma diagonale- \lambda > 1\text{,} il vettore si allunga;
- \lambda = 1\text{,} il vettore non modifica la lunghezza;
- \lambda < 1\text{,} il vettore si accorcia;
- \lambda = 0\text{,} non è ammissibile perché il vettore, per l'ammissibilità della trasformazione, non può annullare la propria lunghezza;
- \lambda < 0\text{,} non è ammissibile perché il vettore, per l'ammissibilità della trasformazione, non può cambiare verso.
Sottoparagrafo 1.7.2
Si supponga adesso per \tens{U} sempre una forma diagonale ma con coefficienti non tutti uguali fra loro. In particolare si assuma- Tutti i vettori appartenenti a direzioni giacenti nel piano individuato dai i primi due assi di riferimento vengono semplicemente scalati dal fattore \lambda_1\text{.}
- Vettori paralleli al terzo asse di riferimento vengono semplicemente scalati dal fattore \lambda_2\text{.}
- Tutti gli altri vettori vengono sia scalati sia cambiano direzione e il fattore di scala dipende sia \lambda_1 e sia da \lambda_2\text{.}
Sottoparagrafo 1.7.3
È ormai chiaro che la forma diagonale più generica del tensore \tens{U} sia data da- Solo gli assi di riferimento \vec{e}_{1}\text{,} \vec{e}_{2} e \vec{e}_{3} costituiscono le direzioni principali del tensore \tens{U}\text{.}
- Gli allungamenti principali \lambda_{1}\text{,} \lambda_{2} e \lambda_{3} sono tutti numeri reali positivi ed associati, rispettivamente, a \vec{e}_{1}\text{,} \vec{e}_{2} e \vec{e}_{3}\text{.} Elementi di linea appartanenti alle direzioni principali per effetto della trasformazione mantegono direzione e verso originari modificando solo la lunghezza sulla base di tali fattori di scala.
Sottoparagrafo 1.7.4 forma spettrale del tensore di deformazione
Osservazione 1.7.1.
Da quando si è cominciato a parlare di direzioni principali non si è mai fatto accenno alla simmetria di \tens{U}\text{.} La simmetria del tensore determina la mutua ortogonalità delle direzioni principali, qualunque direzione esse assumano.
Osservazione 1.7.2.
La forma spettrale (1.7.6) spiega come mai, rispetto alla base di riferimento standard, la forma matriciale del tensore \tens{U} perda quasi sempre la forma diagonale. Infatti, pur essendo sempre presenti le direzioni principali, queste in generale non coincidono con gli assi del sistema di riferimento standard. La perdita della forma diagonale viene evidenziata da una scrittura matriciale del forma spettrale del tensore
Dove per le componenti delle direzioni principali si è adottata la seguente notazione
Infine è interessante osservare come l'espressione (1.7.7) possa essere scritta nel modo seguente
che è la modalità con cui viene valutata l'espressione che assume il tensore \tens{U} quando si passa dal sistema di coordinate costituito dalle direzioni pricipali al sistema di coordinate standard \vec{e}_1\text{,} \vec{e}_2 e \vec{e}_3\text{.} La matrice utilizzata per la trasformazione viene costruita inserendo nelle colonne i vettori \vec{u}_{1}\text{,} \vec{u}_{2} e \vec{u}_{3}\text{.}
Assumendo per \tens{U} la seguente forma diagonale
si calcoli l'espressione che assume \tens{U} supponendo di applicare alle direzioni principali le seguenti rotazioni in successione
Le seguenti istruzioni MATLAB® consentono di calcolare quanto richiesto.
R1 = [cos(pi/4) -sin(pi/4) 0; sin(pi/4) cos(pi/4) 0; 0 0 1]; R2 = [cos(pi/4) 0 -sin(pi/4); 0 1 0; sin(pi/4) 0 cos(pi/4)]; u1 = R2*R1*[1; 0; 0]; u2 = R2*R1*[0; 1; 0]; u3 = R2*R1*[0; 0; 1]; lam1 = 2; lam2 = 5; lam3 = 7; U = lam1*u1*transpose(u1) + lam2*u2*transpose(u2) + ... lam3*u3*transpose(u3)
Si effettui anche il seguente calcolo e si rifletta sul risultato ottenuto.
u1*transpose(u1) + u2*transpose(u2) + u3*transpose(u3)
Al fine di verificare come il tensore \tens{U} agisca in modo differente sulle direzioni principali e sulle direzioni individuate dal riferimento standard, si consideri il seguente stato piano caratterizzato dalla seguente forma diagonale di \tens{U}
Si assuma per le direzioni principali un'inclinazione di \pi/4 e si applichi quindi il tensore della deformazione ai due quadrati, rappresentati nella seguente figura. Il quadrato di colore rosso ha lati allineati con il riferimento standard, il quadrato di colore blu ha invece lati allineati con le direzioni principali.

L'applicazione del tensore viene effettuata facendo crescere il parametro \lambda\text{,} da cui dipendono i due autovalori del tensore, dal valore 1.1 al valore 1.5\text{.}
Si riportano le istruzioni MATLAB® utilizzabili per il calcolo e la visualizzazione delle deformate dei due quadrati a partire dalle configurazioni iniziali.
Q1 = [1 -1 -1 1 1; 1 1 -1 -1 1]; Q2 = [1 0 -1 0 1; 0 1 0 -1 0]; x = Q1(1,:); y = Q1(2,:); pbaspect([1 1 1]) line(x,y,'Color','red','LineWidth',2) x = Q2(1,:); y = Q2(2,:); line(x,y,'Color','blue','LineWidth',2) u1 = [1/sqrt(2); 1/sqrt(2)]; u2 = [-1/sqrt(2); 1/sqrt(2)]; syms U(lam) U(lam) = [u1 u2]*diag([1.2*lam lam])*transpose([u1 u2]); for v = 1.1:0.1:1.5 UQ1 = U(v)*Q1; x = UQ1(1,:); y = UQ1(2,:); line(x,y,'Color','red','LineWidth',2) UQ2 = U(v)*Q2; x = UQ2(1,:); y = UQ2(2,:); line(x,y,'Color','blue','LineWidth',2) end
Sottoparagrafo 1.7.5 autovalori ed autovettori
Noti gli allungamenti principali \lambda_i e le direzioni principali ad essi associate \vec{u}_i (i=1 \dots 3), la forma spettrale (1.7.6) consente di valutare facilmente il tensore \tens{U} rispetto al sistema di riferimento utilizzato. Lo scenario tipico è, comunque, ribaltato. Ovvero dato il tensore \tens{U} si vogliono valutare gli allungamenti e le direzioni principali. Risolvere tale quesito vuol dire calcolare gli autovalori e gli autovettori della matrice associata al tensore \tens{U}\text{.} Nella presente sezione non verrà presentata la procedura da utilizzare nella soluzione del problema agli autovalori perché questa costituisce un argomento che sicuramente è già stato oggetto di studio in corsi precedenti e per il quale la letteratura disponibile, anche in rete, è vasta. In questa sede si è ritenuto più importante mettere in evidenza il significato meccanico di autovalori ed autovettori, significato ampiamente discusso nelle sezioni precedenti. Inoltre per il calcolo degli autovalori ed autovettori si utilizzerà MATLAB® come viene illustrato nel seguente esempio.A partire dall'espressione di \tens{U} ottenuta nell'esempio precedente, ovvero
si calcolino con MATLAB® gli autovalori ed autovettori.
U = ... [5.2500 -1.0607 -1.7500; ... -1.0607 3.5000 -1.0607; ... -1.7500 -1.0607 5.2500]; [V, D] = eig(U)
Intuizione 1.7.8. il polinomio caratteristico.
La valutazione degli autovalori ed autovettori associati al tensore \tens{U} passa attraverso la soluzione della seguente equazione cubica
dove compare il polinomio caratteristico il quale dipendende da tre quantità che prendono il nome di invarianti del tensore \tens{U}\text{.} Il termine invariante sta ad indicare che tali valori non cambiano al variare del sistema di riferimento rispetto al quale viene rappresentato \tens{U}\text{.} Si riportano le espressioni da utilizzare per il calcolo degli invarianti nel caso il tensore \tens{U} sia espresso alle direzioni principali e nel caso generico.
-
Primo invariante\begin{equation} \func{I_1}{\tens{U}} = \lambda_1 + \lambda_2 + \lambda_3 = \trace{\tens{U}}\,,\tag{1.7.8} \end{equation}dove compare la traccia del tensore ovvero la somma dei coefficienti appartenenti alla diagonale principale.
-
Secondo invariante\begin{equation} \func{I_2}{\tens{U}} = \lambda_1\, \lambda_2 + \lambda_2\, \lambda_3 + \lambda_3\, \lambda_1 = \frac{1}{2}\left[(\trace{\tens{U}})^2 - \trace{\tens{U}^2}\right]\,.\tag{1.7.9} \end{equation}
-
Terzo invariante\begin{equation} \func{I_3}{\tens{U}} = \lambda_1\, \lambda_2\, \lambda_3 = \det{\tens{U}}\,.\tag{1.7.10} \end{equation}
Si verifichi che gli invarianti calcolati su \tens{U} nella forma diagonale
e nella forma
non cambiano.
Si riportano le istruzioni MATLAB® utilizzabili per la verifica.
lam1 = 2; lam2 = 5; lam3 = 7; I1 = lam1 + lam2 + lam3 I2 = lam1*lam2 + lam2*lam3 + lam3*lam1 I3 = lam1*lam2*lam3 U = ... [5.2500 -1.0607 -1.7500; ... -1.0607 3.5000 -1.0607; ... -1.7500 -1.0607 5.2500]; I1 = trace(U) I2 = 1/2*(trace(U)^2 - trace(U*U)) I3 = det(U)