Strumenti Utente

Strumenti Sito


wikipaom2019:lez_2019-03-21lab

Differenze

Queste sono le differenze tra la revisione selezionata e la versione attuale della pagina.

Link a questa pagina di confronto

Entrambe le parti precedenti la revisione Revisione precedente
Prossima revisione
Revisione precedente
wikipaom2019:lez_2019-03-21lab [2019/05/20 12:44]
249692
wikipaom2019:lez_2019-03-21lab [2019/07/11 09:13] (versione attuale)
ebertocchi [Soluzione di Michell]
Linea 1: Linea 1:
  
 +A cura di Antonio Loriso, Francesco della torca, Francesco Marino 
 +{{ :wikipaom2019:lezione_21-03-19_.pdf |}}
 +{{ :wikipaom2019:lezione_21-03-19_1_1_.docx |}}
 +
 +A cura di Fabio Veruschi, Gianmarco Rigon
 +{{ :wikipaom2019:lezione_21-03-2019_part2.pdf |}}
 +
 +
 +----
 +
 +====== Area a cura del Docente ======
 +
 +===== Soluzione di Michell =====
 +
 +testo di riferimento: J.R. Barber, Elasticity, da [[restricted:materiale_didattico_copyright|Materiale corsi di NON libera distribuzione]], .
 +
 +{{ :wikipaom2019:estratto_michell_barber_elasticity_p119_p130.pdf |estratto termini soluzione di Michell}}
 +
 +
 +Argomenti di riferimento:
 +  * stato piano di tensione e stato piano di deformazione;
 +  * costante di Kolosov per tensione e deformazione piana, definizione delle comp. entro piano di deformazione p.43 (59 pdf);
 +    * DP: $\kappa=\left(3-4\nu\right)$; TP: $\kappa=\left(\frac{3-\nu}{1+\nu}\right)$;
 +    * $\epsilon_x=\left(\frac{\kappa+1}{8\mu} \right) \sigma_x - \left(\frac{3-\kappa}{8\mu} \right) \sigma_y$
 +    * $\epsilon_y=\left(\frac{\kappa+1}{8\mu} \right) \sigma_y - \left(\frac{3-\kappa}{8\mu} \right) \sigma_x$
 +    * $\gamma_{xy}=\frac{\tau_{xy}}{\mu}$
 +    * $\kappa$ e modulo di taglio $\mu$ definiscono completamente il legame costitutivo per un materiale isotropo in stati piani.
 +  * componenti fuori piano di tensione e deformazione: 
 +    * TP: $\epsilon_z=-\frac{\nu}{1-\nu}\left(\epsilon_x + \epsilon_y \right)=-\frac{\nu}{E}\left(\sigma_x+\sigma_y\right)$
 +    * DP: $\sigma_z = \nu \left( \sigma_x+\sigma_y \right)$
 +  * deformazione piana generalizzata come sovrapposizione ad uno stato di deformazione piana di una soluzione $\epsilon_z=\bar{\epsilon}-\frac{1}{\rho_y}x+\frac{1}{\rho_x}y$ costruita in compensazione delle risultanti di sforzo normale e momento flettente;
 +  * equazioni di equilibrio in stati piani, in eventuale presenza di azioni distribuite $q_x$ e $p_y$:
 +    * eq. tx: $\frac{\partial \sigma_x}{\partial x} + \frac{\partial \tau_{xy}}{\partial y}+p_x=0$ 
 +    * eq. ty: $\frac{\partial \tau_{xy}}{\partial x} + \frac{\partial \sigma_y}{\partial y}+p_y=0$ 
 +  * rotazione del sistema di riferimento per le componenti di tensione, p. 9 (27 pdf), da utilizzarsi nel passaggio da tensioni in coordinate cartesiane a tensioni in coordinate polari 
 +
 +{{:wikipaom2017:componenti_di_tensione_da_cartesiano_a_polare_v000.png?600|}}
 +
 +{{:wikipaom2017:componenti_di_tensione_da_cartesiano_a_polare_v000.pdf|sorgente ipe}}
 +  * equazioni di equilibrio in coordinate polari, basato su di un settore di corona circolare di ampiezza radiale $dr$ e ampiezza angolare $d\theta$:
 +    * eq. tr. radiale: $\frac{\partial \sigma_r}{\partial r}+\frac{1}{r}\frac{\partial \tau_{r\theta}}{\partial \theta}+\frac{\sigma_r-\sigma_\theta}{r}+p_r=0$
 +    * eq. tr. circonf.:$\frac{1}{r}\frac{\partial \sigma_\theta}{\partial \theta}+\frac{\partial \tau_{r\theta}}{\partial r}+\frac{2\tau_{r\theta}}{r}+p_\theta=0$
 +  * relazioni spostamento-deformazione in coordinate polari:
 +    * $\epsilon_r=\frac{\partial u}{\partial r}$ 
 +    * $\epsilon_\theta=\frac{1}{r}\frac{\partial v}{\partial \theta}+\frac{u}{r}$
 +    * $\gamma_{r\theta}=\frac{\partial v}{\partial r}+\frac{1}{r}\frac{\partial u}{\partial \theta} -\frac{v}{r}$
 +  * operatore laplaciano e bilaplaciano, p. 49 (65 pdf) in coord. cartesiane, p. 111 (125 pdf) in coordinate polari;
 +    * coord. cart.: $\nabla^2=\left(\frac{\partial^2}{\partial x^2}+\frac{\partial^2}{\partial y^2}\right)$
 +    * coord. polari:$\nabla^2=\left(\frac{\partial^2}{\partial r^2}+\frac{1}{r}\frac{\partial}{\partial r}+\frac{1}{r^2}\frac{\partial^2}{\partial \theta^2}\right)$
 +    * bilaplaciano: $\nabla^4 \phi=\nabla^2 \nabla^2 \phi$
 +  * Airy stress function $\phi$, p. 46 (62 pdf); 
 +    * l'associata definizione (4.6) di $\sigma_{xx},\sigma_{yy},\sigma_{xy}=\tau_{xy}$ soddisfa nativamente (=sulla base del teorema di Schwarz) le equazioni di equilibrio in assenza di azione distribuita (4.1), in particolare:
 +      * $\sigma_x=\frac{\partial^2\phi}{\partial y^2}$
 +      * $\sigma_y=\frac{\partial^2\phi}{\partial x^2}$
 +      * $\tau_{xy}=-\frac{\partial^2\phi}{\partial x\partial y}$
 +    * l'equazione $\nabla^4 \phi=0$ caratterizzante la funzione di Airy deriva dall'equazione di compatibilità, una volta sostituita in essa il legame elastico omogeneo isotropo, p. 48 (64 pdf);
 +  * note sull'equazione di compatibilità:
 +    * uno stato di deformazione è compatibile se è definibile in termini di un campo di spostamenti monodromo,differenziabile e a derivate parziali continue;
 +    * uno stato di deformazione è compatibile se lo spostamento relativo tra due punti A e B entro il solido elastico è definibile per accumulo (integrazione) dei contributi deformativi su di un percorso A->B, e se tale integrale è indipendente dal percorso, per piccole variazioni del percorso stesso.
 +    * uno stato di deformazione è compatibile se non genera dislocazioni nel corpo elastico, ove non ne preesistessero.
 +  * derivazione di $\sigma_{rr},\sigma_{\theta\theta},\tau_{r\theta}$ da $\phi$, p. 110 (124);
 +    *$\sigma_r=\frac{1}{r}\frac{\partial \phi}{\partial r} + \frac{1}{r^2}\frac{\partial^2 \phi}{\partial \theta^2}$
 +    *$\sigma_\theta=\frac{\partial^2\phi}{\partial r^2}$
 +    *$\tau_{r\theta}=\frac{1}{r^2}\frac{\partial\phi}{\partial\theta}-\frac{1}{r}\frac{\partial^2\phi}{\partial r \partial \theta}$ 
 +  * termini della soluzione di Michell, componenti di tensione p. 119 (133), inseriti nel foglio maxima di seguito;
 +    * gli stati tensionali, deformativi e di spostamento descritti dai termini della soluzione di Michell (o funzioni di Airy in generale) sono stati "naturali" o "propri" del materiale elastico nel senso che possono sussistere in assenza di azioni esterne applicate ai punti interni del dominio
 +  * termini della soluzione di Michell, spostamenti p. 130 (144).
 +===== Lastra forata =====
 +
 +==== Stralci di codice da inserire ====
 +
 +Conversione delle componenti di tensione (stati piani) da sistema cartesiano a polare, cfr. pp. 8-9 Barber.
 +
 +  define(
 +      srr_from_xy(sxx,syy,sxy,t),
 +      sxx*c^2 + syy*s^2+2*sxy*s*c
 +  ), [c = cos(t) , s = sin(t)];
 +
 +  define(
 +      srt_from_xy(sxx,syy,sxy,t),
 +      sxy*(c^2-s^2)+(syy-sxx)*s*c
 +  ), [c = cos(t) , s = sin(t)];
 +
 +  define(
 +      stt_from_xy(sxx,syy,sxy,t),
 +      syy*c^2 + sxx*s^2-2*sxy*s*c
 +  ), [c = cos(t) , s = sin(t)];
 +
 +Tabelle termini Michell utilizzati
 +<code>
 +philist : [ 
 +    r^2 , 
 +    log(r) ,
 +    t, 
 +    r^(-n+2)*cos(n*t), 
 +    r^n*cos(n*t) , 
 +    r^(-n)*cos(n*t)
 +],n=2;
 +</code>
 +
 +
 +<code>
 +twomu_ur_list : [ 
 +    (kappa-1)*r , 
 +    -1/r ,
 +    0, 
 +    (kappa+n-1)*r^(-n+1)*cos(n*t), 
 +            -n *r^( n-1)*cos(n*t) , 
 +             n *r^(-n-1)*cos(n*t)
 +],n=2;
 +</code>
 +
 +<code>
 +twomu_ut_list : [ 
 +    0 , 
 +    0 ,
 +    -1/r, 
 +   -(kappa-n+1)*r^(-n+1)*sin(n*t), 
 +             n *r^( n-1)*sin(n*t) , 
 +             n *r^(-n-1)*sin(n*t)
 +],n=2;
 +</code>
 + 
 +Derivazione dei termini di tensione dai relativi della funzione di Airy
 +<code>
 +srr : 1/r * diff( phi , r , 1) + 1/r^2*diff( phi , t , 2 );
 +stt : diff( phi , r , 2 );
 +srt : 1/r^2 * diff(phi, t,1) - 1/r * diff(phi,r,1,t,1);
 +</code>
 +
 +==== file maxima lato cattedra ====
 +
 +{{ :wikipaom2019:foro_in_lastra_infinita_a2019_v001.wxmx | concentrazione di tensioni in lastra forata, foglio a fine lezione}}
 +
 +~~DISCUSSION~~