Visualizza wikitesto di Il problema di Saint Venant

{{Risorsa|tipo=lezione|materia1=Scienza delle costruzioni}}
Il '''problema di Saint Venant''' è un particolare problema dell'equilibrio elastico, proposto e risolto da [[w:Adhémar Jean Claude Barré de Saint-Venant|de Saint Venant]] alla metà del XIX secolo. 

==Posizione del problema==
{{Todo|Inserire immagine}}
Il problema di Saint Venant si riduce a valutare la [[I teoremi sull'elasticità lineare|soluzione dell'equilibrio elastico]] di un corpo avente alcune particolari caratteristiche geometriche, statiche e materiali.

Dal punto di vista geometrico il corpo in questione è un cilindro retto (quindi di sezione costante) a sezione trasversale qualsiasi di area <math>A</math>, "''sufficientemente allungato''", intendendo con quest'ultima affermazione che lo sviluppo longitudinale del cilindro sia molto maggiore della sua sezione, e cioè che la dimensione massima della sezione sia trascurabile rispetto all'intera lunghezza del corpo.

Il materiale di cui è composto il corpo è ipotizzato linearmente elastico, isotropo e omogeneo.

Staticamente si suppone che il corpo non sia soggetto a forze di volume (forza peso, forza centrifuga, etc), e abbia la sua superficie laterale scarica, e cioè priva di forze applicate. Le uniche forze applicate, dunque, insistono sulle due basi, le quali possono essere sottoposte ad una qualsiasi generica distribuzione di sforzi (anche differente tra le due basi), a patto che complessivamente le azioni agenti sul sistema siano equilibrate. Il corpo, inoltre, è supposto privo di vincoli, a parte un vincolo che si introduce esclusivamente per evitare traslazioni e rotazioni rigide del corpo stesso, e cioè tutti quei movimenti che lo interessano come se fosse un corpo rigido.

Il problema, dunque, è un ''problema nelle forze'', dal momento che su tutta la frontiera del corpo sono assegnate esclusivamente delle forze e non esistono porzioni di frontiera vincolate<ref>Al contrario, se fossero assegnati degli spostamenti lungo tutta la frontiera sarebbe stato un ''problema negli spostamenti'' (frontiera tutta vincolata). Il caso intermedio è definito ''problema misto''.</ref>.

Questo è il più importante dei problemi dell'equilibrio elastico in campo ingegneristico, dal momento che la quasi totalità delle strutture moderne condivide le caratteristiche geometriche del corpo studiato: in linea generale, infatti, la maggior parte delle strutture portanti moderne è costituita da travi e pilastri, e cioè strutture che effettivamente sono molto allungate e che molto raramente hanno sezioni variabili. La presenza di forze di volume nelle strutture reali (forza peso, azioni sismiche) e l'esistenza dei vincoli sarà tenuta in conto in seguito.

==Le equazioni alla base del problema==
Supponendo di considerare un sistema di riferimento <math>(O,x,y,z)</math> tale che <math>O</math> sia il baricentro della sezione di sinistra, <math>z</math> sia coincidente all'asse del cilindro e <math>x,y</math> siano due generiche direzioni appartenenti al piano della sezione trasversale, le equazioni che reggono il problema sono quelle di congruenza, equilibrio, costitutive e al contorno indicate di seguito:

Equazioni di equilibrio:
<math>\begin{cases}
\frac{\delta \sigma_{x}}{\delta x}+\frac{\delta \tau_{xy}}{\delta y}+\frac{\delta \tau_{xz}}{\delta z}=0 \\
\frac{\delta \tau_{xy}}{\delta x}+\frac{\delta \sigma_{y}}{\delta y}+\frac{\delta \tau_{yz}}{\delta z}=0 \\
\frac{\delta \tau_{xz}}{\delta x}+\frac{\delta \tau_{yz}}{\delta y}+\frac{\delta \sigma_{z}}{\delta z}=0
\end{cases}</math>

Equazioni di congruenza:
<math>\begin{cases}
\varepsilon_x=\frac{\delta v_x}{\delta x} \\
\varepsilon_{xy}=\frac{1}{2}\left(\frac{\delta v_x}{\delta y}+\frac{\delta v_y}{\delta x}\right) \\
\varepsilon_{xz}=\frac{1}{2}\left(\frac{\delta v_x}{\delta z}+\frac{\delta v_z}{\delta x}\right) \\
\varepsilon_y=\frac{\delta v_y}{\delta y} \\
\varepsilon_{yz}=\frac{1}{2}\left(\frac{\delta v_y}{\delta z}+\frac{\delta v_z}{\delta y}\right) \\
\varepsilon_z=\frac{\delta v_z}{\delta z}
\end{cases}</math>

Equazioni costitutive:
<math>\begin{cases}
\varepsilon_{x}=\frac{1}{E} \left[ \sigma_{x} - \nu (\sigma_{y}+\sigma_{z}) \right] \\
\varepsilon_{y}=\frac{1}{E} \left[ \sigma_{y} - \nu (\sigma_{x}+\sigma_{z}) \right] \\
\varepsilon_{z}=\frac{1}{E} \left[ \sigma_{z} - \nu (\sigma_{x}+\sigma_{y}) \right] \\
\gamma_{xy}=\frac{1}{G}\tau_{xy} \\
\gamma_{xz}=\frac{1}{G}\tau_{xz} \\
\gamma_{yz}=\frac{1}{G}\tau_{yz}
\end{cases}</math>

Condizioni al contorno sulla superficie laterale (con <math>n</math> direzione normale alla superficie laterale):
<math>\begin{cases}
\sigma_xn_x+\tau_{xy}n_y=0 \\
\tau_{yx}n_x+\sigma_yn_y=0 \\
\tau_{zx}n_x+\tau_{zy}n_y=0
\end{cases}</math><ref>Si fa notare che tali posizioni corrispondono effettivamente ad imporre che le componenti del vettore di tensione lungo i tre assi siano nulle. Infatti la componente secondo l'asse <math>x</math> del vettore di tensione è espresso da:

<math>\sigma_xn_x+\tau_{xy}n_y+\tau_{xz}n_z</math>

ed ugualmente per le altre due direzioni. Tuttavia, poiché il cilindro è retto per ipotesi, in ogni punto della superficie laterale la normale <math>n</math> è anche normale all'asse del cilindro, per cui la sua componente lungo tale direzione è nulla (<math>n_z=0</math>), per cui scompaiono tutti i termini che lo contengono. Porre nulle le componenti del vettore di tensione secondo le tre direzioni ortogonali, infine, corrisponde ad imporre che sia nullo il vettore stesso.</ref>

Condizioni al contorno sulle due basi:
<math>\begin{cases}
\tau_{zx}=p_x \\
\tau_{zy}=p_y \\
\sigma_z=p_z
\end{cases}
(z=l)</math>

<math>\begin{cases}
\tau_{zx}=-p_x' \\
\tau_{zy}=-p_y' \\
\sigma_z=-p_z'
\end{cases}
(z=0)</math><ref>Dal momento che sulle basi la normale alla superficie ha la stessa direzione dell'asse <math>z</math>, le componenti del vettore di tensione si riducono ad essere quelle espresse nelle condizioni precedenti. A differenza del caso della superficie laterale, tuttavia, le tensioni non si devono annullare ma devono avere un valore uguale all'azione esterna. Si fa notare che si considerano tutte le componenti del vettore di tensione: ad un primo impatto, infatti, si potrebbe ritenere necessario considerare solo la <math>\sigma_z</math>, ma ad un'osservazione più attenta dovrebbe rivelarsi ovvio che è necessario considerare anche la possibilità che sulle basi vengano impresse azioni tangenziali</ref>

==Il metodo risolutivo di Saint Venant==
Il problema, posto in questi termini, presenta enormi difficoltà analitiche per la sua risoluzione. Lo stesso de Saint Venant non giunse alla soluzione del problema risolvendo direttamente questi sistemi di equazioni differenziali, ma per mezzo di un metodo ''seminverso'': tale metodo consiste nell'imporre a priori alcune caratteristiche della soluzione cercata per poi ricercarne le altre per mezzo delle equazioni a disposizione. La scelta di tali caratteristiche, naturalmente, deve essere tale da presentare una certa validità fisica.

Le condizioni che de Saint Venant impose essere soddisfatte a priori in ogni punto del corpo furono le seguenti:

<math>\sigma_x=\sigma_y=\tau_{xy}=0</math>

{{Cassetto|Significato fisico di queste ipotesi| 
{{Todo|Inserire immagine}}
Per comprendere cosa significa fisicamente questa imposizione matematica si consideri una generica fibra elementare del corpo parallela all'asse <math>z</math> del corpo stesso. Noto che le generiche componenti di tensioni <math>\sigma_{n\nu}</math> possono essere espresse in funzione delle componenti di tensione secondo tre direzioni normali <math>1,2,3</math> nel modo seguente:

<math>\sigma_{n\nu}=\sum_{i,k=1}^3 \sigma_{ik}n_i\nu_k</math>

Esplicitando le componenti della tensione e applicandolo al nostro caso (in cui gli assi sono <math>x,y,z</math>) si ottiene:

<math>\begin{cases}
\sigma_n=\sigma_xn_x^2+2\tau_{xy}n_xn_y+2\tau_{xz}n_xn_z+\sigma_yn_y^2+\tau_{yz}n_yn_z+\sigma_zn_z^2 \\
\tau_{n\nu}=\sigma_xn_x\nu_x+\tau_{xy}n_x\nu_y+\tau_{xz}n_x\nu_z+\tau_{yx}n_y\nu_x+\sigma_yn_y\nu_y+\tau_{yz}n_y\nu_z+\tau_{zx}n_z\nu_x+\tau_{zy}n_z\nu_y+\sigma_zn_z\nu_z \\
\tau_{nz}=\sigma_xn_xz_x+\tau_{xy}n_xz_y+\tau_{xz}n_xz_z+\tau_{yx}n_yz_x+\sigma_yn_yz_y+\tau_{yz}n_yz_z+\tau_{zx}n_zz_x+\tau_{zy}n_zz_y+\sigma_zn_zz_z
\end{cases}</math>

Poiché il cilindro è retto, i coseni direttori delle tre direzioni sono <math>n(n_x,n_y,0)\;\nu(\nu_x,\nu_y,0)\;z(0,0,1)</math><ref>cfr. nota 2</ref> per cui le equazioni precedenti si trasformano nelle seguenti:

<math>\begin{cases}
\sigma_n=\sigma_xn_x^2+2\tau_{xy}n_xn_y+\sigma_yn_y^2 \\
\tau_{n\nu}=\sigma_xn_x\nu_x+\tau_{xy}n_x\nu_y+\tau_{yx}n_y\nu_x+\sigma_yn_y\nu_y \\
\tau_{nz}=\tau_{xz}n_x+\tau_{yz}n_y
\end{cases}</math>

Sostituendo le espressioni precedentemente ipotizzate si ottiene:

<math>\begin{cases}
\sigma_n=0 \\
\tau_{n\nu}=0 \\
\tau_{nz}=\tau_{xz}n_x+\tau_{yz}n_y
\end{cases}</math>

L'unica componente di tensione non nulla, cioè, è la <math>\tau_{nz}</math>.

Questo equivale a dire, in pratica, che le fibre del corpo non si scambiano tra loro azioni normali, dal momento che la trattazione precedente è stata fatta con il massimo della generalità. Le uniche azioni che le fibre longitudinali del corpo possono offrirsi mutuamente sono esclusivamente azioni tangenziali dirette secondo la direzione parallela all'asse del corpo stesso. Le azioni mutue tra le varie fibre del corpo, cioè, si riducono ad essere solo azioni di attrito.
}}

In virtù di tale ipotesi il tensore di tensione assume la forma seguente:

<math>[T]=\begin{vmatrix}
0 & 0 & \tau_{xz} \\
0 & 0 & \tau_{yz} \\
\tau_{zx} & \tau_{zy} & \sigma_z
\end{vmatrix}</math><ref>Lo stato di tensione descritto da questo tensore si dimostra essere al più biassiale. Naturalmente il tensore può degenerare in uno stato di tensione monoassiale (se <math>\tau_{xz}=\tau_{yz}=0</math></ref>

Con tale ipotesi le equazioni di equilibrio e costitutive precedenti si riducono a:

<math>\begin{cases}
\frac{\delta \tau_{xz}}{\delta z}=0 \\
\frac{\delta \tau_{yz}}{\delta z}=0 \\
\frac{\delta \tau_{xz}}{\delta x}+\frac{\delta \tau_{yz}}{\delta y}+\frac{\delta \sigma_{z}}{\delta z}=0
\end{cases}</math>

<math>\begin{cases}
\epsilon_{x}=-\frac{1}{E}  \nu \sigma_{z} \\
\epsilon_{y}=-\frac{1}{E} \nu \sigma_{z} \\
\epsilon_{z}=\frac{1}{E} \sigma_{z} \\
\gamma_{xy}=0 \\
\gamma_{xz}=\frac{1}{G}\tau_{xz} \\
\gamma_{yz}=\frac{1}{G}\tau_{yz}
\end{cases}</math><ref>In base a quanto scritto si deduce che due fibre parallele agli assi <math>x,y</math> si mantengono ortogonali tra loro dopo la deformazione (dal momento che <math>\gamma_{xy}=0</math> rappresenta lo scorrimento mutuo tra queste due direzioni). Dal momento che non sono state fatte specificazioni di sorta sulle direzioni in analisi eccetto di appartenere al piano della sezione trasversale, si può affermare che due direzioni qualsiasi appartenenti al piano della sezione trasversale originariamente ortogonali si mantengono tali anche dopo la deformazione</ref>

Le condizioni al contorno lungo la superficie laterale si riducono a:

<math>\begin{cases}
0=0 \\
0=0 \\
\tau_{zx}n_x+\tau_{zy}n_y=0
\end{cases}</math><ref>Tale espressione, che è l'unica non identicamente nulla in questa serie di condizioni, equivale a dire che la tensione lungo la frontiera del corpo ha direzione tangente alla frontiera stessa. In forma vettoriale, infatti, l'espressione equivale al prodotto scalare tra il vettore della tensione <math>t</math> e la direzione perpendicolare al contorno <math>n</math>, il quale è nullo solo quando i due vettori considerati sono ortogonali tra loro</ref>

Le equazioni al contorno lungo le basi e quelle di congruenza si mantengono, invece, inalterate:

<math>\begin{cases}
\tau_{zx}=p_x \\
\tau_{zy}=p_y \\
\sigma_z=p_z
\end{cases}
(z=l)</math>

<math>\begin{cases}
\tau_{zx}=-p_x' \\
\tau_{zy}=-p_y' \\
\sigma_z=-p_z'
\end{cases}
(z=0)</math>

<math>\begin{cases}
\epsilon_x=\frac{\delta v_x}{\delta x} \\
\epsilon_{xy}=\frac{1}{2}\left(\frac{\delta v_x}{\delta y}+\frac{\delta v_y}{\delta x}\right) \\
\epsilon_{xz}=\frac{1}{2}\left(\frac{\delta v_x}{\delta z}+\frac{\delta v_z}{\delta x}\right) \\
\epsilon_y=\frac{\delta v_y}{\delta y} \\
\epsilon_{yz}=\frac{1}{2}\left(\frac{\delta v_y}{\delta z}+\frac{\delta v_z}{\delta y}\right) \\
\epsilon_z=\frac{\delta v_z}{\delta z}
\end{cases}</math>

Volendo esprimere tutte le equazioni in funzione solo delle componenti di spostamento <math>v_i</math>, basta partire dall'ultima serie di equazioni riportate e sostituire via via tutti i termini. 

Dalle equazioni costitutive si ottiene prima di tutto che, dovendo essere <math>\epsilon_x=\epsilon_y</math>, deve essere:

<math>\frac{\delta v_x}{\delta x}=\frac{\delta v_y}{\delta y}=-\nu\frac{\delta v_z}{\delta z}</math>

poi, dovendo essere <math>\gamma_{xy}=0</math>:

<math>\frac{\delta v_x}{\delta y}=-\frac{\delta v_y}{\delta x}</math>

e infine:

<math>\begin{cases}
\sigma_{z}=-\frac{E}{\nu}\frac{\delta v_x}{\delta x}=-\frac{E}{\nu}\frac{\delta v_y}{\delta y}=E\frac{\delta v_z}{\delta z} \\
\tau_{xz}=G\left(\frac{\delta v_x}{\delta z}+\frac{\delta v_z}{\delta x}\right) \\
\tau_{yz}=G \left(\frac{\delta v_y}{\delta z}+\frac{\delta v_z}{\delta y}\right)
\end{cases}</math>

Dalle equazioni di equilibrio si deduce:

<math>\begin{cases}
\frac{\delta}{\delta z} \left(\frac{\delta v_x}{\delta z}+\frac{\delta v_z}{\delta x}\right)=0 \\
\frac{\delta}{\delta z}\left(\frac{\delta v_y}{\delta z}+\frac{\delta v_z}{\delta y}\right)=0 \\
\frac{\delta^2 v_z}{\delta x^2}+\frac{\delta^2v_z}{\delta y^2}+2\frac{\delta^2v_z}{\delta z^2}=0
\end{cases}</math>

Le condizioni ai limiti sulla superficie laterale impongono:

<math>\left(\frac{\delta v_x}{\delta z}+\frac{\delta v_z}{\delta x}\right)n_x+\left(\frac{\delta v_y}{\delta z}+\frac{\delta v_z}{\delta y}\right)n_y=0</math>

Sulle basi si ha:

<math>\begin{cases}
G\left(\frac{\delta v_x}{\delta z}+\frac{\delta v_z}{\delta x}\right)=p_x \\
G \left(\frac{\delta v_y}{\delta z}+\frac{\delta v_z}{\delta y}\right)=p_y \\
E\frac{\delta v_z}{\delta z}=p_z
\end{cases}
(z=l)</math>

<math>\begin{cases}
-G\left(\frac{\delta v_x}{\delta z}+\frac{\delta v_z}{\delta x}\right)=p_x' \\
-G \left(\frac{\delta v_y}{\delta z}+\frac{\delta v_z}{\delta y}\right)=p_y' \\
-E\frac{\delta v_z}{\delta z}=p_z'
\end{cases}
(z=0)</math>

== Note ==

<references/>
== Altri progetti ==
{{Interprogetto}}