Enunciado del trabajo 2 - Historial de revisiones

Carlos Castro en 11:30 9 dic 2013

2013-12-09T11:30:29Z

Carlos Castro en 11:29 9 dic 2013

2013-12-09T11:29:09Z

Carlos Castro en 08:33 25 nov 2013

2013-11-25T08:33:42Z

Carlos Castro en 08:24 25 nov 2013

2013-11-25T08:24:32Z

Carlos Castro en 08:42 22 nov 2013

2013-11-22T08:42:06Z

Carlos Castro en 08:37 22 nov 2013

2013-11-22T08:37:28Z

Carlos Castro en 15:14 21 nov 2013

2013-11-21T15:14:00Z

Carlos Castro: Página creada con «'''Visualización de campos escalares y vectoriales en elasticidad.''' Mallado de una placa plana rectangular Consideramos una...»

2013-11-21T15:00:53Z

Página creada con «'''Visualización de campos escalares y vectoriales en elasticidad.''' Mallado de una placa plana rectangular Consideramos una...»

Página nueva

'''Visualización de campos escalares y vectoriales en elasticidad.'''

[[Archivo:Mallado1.jpg|300px|thumb|right|Mallado de una placa plana rectangular]]

Consideramos una placa plana (en dimensión 2) que ocupa la región <math> [-1/2,1/2] \times [0,2]</math>. En ella vamos a suponer que tenemos definidas dos cantidades físicas: la temperatura <math>T(x,y,t)</math>, que depende de las dos variables espaciales <math>(x,y)</math> y el tiempo <math>t</math>, y los desplazamientos <math>\vec u(x,y,t)</math>. De esta forma, si definimos <math>r_0(x,y)</math> el vector de posición de los puntos de la placa en reposo, la posición de cada punto <math>(x,y)</math> de la placa en un instante de tiempo <math>t</math> viene dada por:
<math>
\vec r (x,y,t)= \vec r_{0}(x,y)+\vec u(x,y,t).
</math>
Vamos a suponer que sobre la placa se ha aplicado una fuerza que ha provocado una vibración de manera que los desplazamientos vienen dados por la onda:
<math>
\vec u(x,y,t) = \vec a \sin(\vec b \cdot \vec r_0-ct),
</math>
donde <math>\vec a</math> se conoce como amplitud, <math>\vec b</math> es la fase que indica la dirección de propagación y <math>c/|\vec b|</math> es la velocidad de propagación.

Si <math>\vec a </math> es paralelo a <math>\vec b</math> diremos que la onda es longitudinal mientras que si es perpendicular hablaremos de onda transversal.

En este trabajo vamos a centrarnos en las ondas longitudinales. Supondremos lo siguiente:
<math>
\vec a=\frac{\vec j}{10}, \qquad \vec b= \pi \vec j, \qquad t=0.
</math>
En este caso, <math>\vec u(x,y)=\frac{\sin(\pi y)}{10}\vec j</math>.

Se pide:
# Dibujar un mallado que represente los puntos interiores del sólido (ver figura). Tomar como paso de muestreo <math>h=1/10</math> para las variables <math>x</math> e <math>y</math>.
# La temperatura del sólido viene dada por <math>T(\rho,\theta)=e^{-y}</math>. Dibujarla.
# Calcular <math>\nabla T</math> y dibujarlo como campo vectorial. Calcular las curvas de nivel de T y observar gráficamente que <math>\nabla T</math> es ortogonal a dichas curvas.
# Consideramos ahora el campo de vectores <math>\vec u=\frac{\sin(y)}{10}\vec i</math>. Dibujar el campo de vectores en los puntos del mallado del sólido.
# Si <math>\vec u</math> determina el desplazamiento que sufre cada punto del sólido, dibujar el sólido antes y después del desplazamiento. Dibujar ambos en la misma figura usando el comando subplot.
# Calcular la <math>\nabla \cdot \vec u</math> en todos los puntos del sólido y dibujarla. ¿Qué puntos tienen mayor divergencia? La divergencia es una medida del cambio de volumen local debido al desplazamiento.
# Calcular <math>|\nabla \times \vec u|</math> en todos los puntos del sólido y dibujarlo. ¿Qué puntos sufren un mayor rotacional?
# Definamos <math>\epsilon(\vec u)=(\nabla \vec u + \nabla \vec u^t)/2</math>, la parte simétrica del tensor gradiente de <math>\vec u</math>. En un medio elástico lineal, isótropo y homogéneo los desplazamientos permiten escribir el tensor de tensiones <math>\sigma_{ij}</math> a través de la fórmula:
<math>
\sigma_{ij}=\lambda \nabla \cdot \vec u \delta_{ij} + 2\mu \epsilon_{ij},
</math>
donde <math>\lambda</math> y <math>\mu</math> son los conocidos como coeficientes de Lamé que dependen de las propiedades elásticas de cada material. Tomando <math>\lambda=\mu=1</math>, dibujar las tensiones normales en la dirección que marca el eje <math>\vec i</math>, es decir <math>\vec i \cdot \sigma \cdot \vec i</math> y las tensiones normales en la dirección que marca el eje <math>\vec j</math>, es decir <math>\vec j \cdot \sigma \cdot \vec j</math>.
# Dibujar las tensiones tangenciales respecto al plano ortogonal a <math>\vec i</math>, es decir <math>|\sigma \cdot \vec i-(\vec i \cdot \sigma \cdot \vec i) \vec i|</math>. ¿Donde son mayores? Compararlas con los puntos de mayor deformación de la malla.
# Dibujar las tensiones tangenciales respecto al plano ortogonal a <math>\vec j</math>, es decir <math>|\sigma \cdot \vec j-(\vec j \cdot \sigma \cdot \vec j) \vec j|</math>. ¿Donde son mayores? Compararla con los puntos de mayor deformación de la malla.

[[Categoría:Teoría de Campos]]

@@ Línea 1: / Línea 1: @@
 '''Visualización de campos escalares y vectoriales en elasticidad.'''
-[[Archivo:Mallado1.jpg|300px|thumb|middle|Mallado de una placa plana rectangular]]
+[[Archivo:Mallado1.jpg|300px|thumb|right|Mallado de una placa plana rectangular]]
 Consideramos una placa plana (en dimensión 2) que ocupa la región <math> [-1/2,1/2] \times [0,2]</math>. En ella vamos a suponer que tenemos definidas dos cantidades físicas: la temperatura <math>T(x,y,t)</math>, que depende de las dos variables espaciales <math>(x,y)</math> y el tiempo <math>t</math>, y los desplazamientos <math>\vec u(x,y,t)</math>. De esta forma, si definimos <math>r_0(x,y)</math> el vector de posición de los puntos de la placa en reposo, la posición de cada punto <math>(x,y)</math> de la placa en un instante de tiempo <math>t</math> viene dada por:

@@ Línea 3: / Línea 3: @@
 [[Archivo:Mallado1.jpg|300px|thumb|right|Mallado de una placa plana rectangular]]
-Consideramos una placa plana (en dimensión 2) que ocupa la región <math> [-1/2,1/2] \times [0,4]</math>. En ella vamos a suponer que tenemos definidas dos cantidades físicas: la temperatura <math>T(x,y,t)</math>, que depende de las dos variables espaciales <math>(x,y)</math> y el tiempo <math>t</math>, y los desplazamientos <math>\vec u(x,y,t)</math>. De esta forma, si definimos <math>r_0(x,y)</math> el vector de posición de los puntos de la placa en reposo, la posición de cada punto <math>(x,y)</math> de la placa en un instante de tiempo <math>t</math> viene dada por:
+Consideramos una placa plana (en dimensión 2) que ocupa la región <math> [-1/2,1/2] \times [0,2]</math>. En ella vamos a suponer que tenemos definidas dos cantidades físicas: la temperatura <math>T(x,y,t)</math>, que depende de las dos variables espaciales <math>(x,y)</math> y el tiempo <math>t</math>, y los desplazamientos <math>\vec u(x,y,t)</math>. De esta forma, si definimos <math>r_0(x,y)</math> el vector de posición de los puntos de la placa en reposo, la posición de cada punto <math>(x,y)</math> de la placa en un instante de tiempo <math>t</math> viene dada por:
 <math>
 \vec r (x,y,t)= \vec r_{0}(x,y)+\vec u(x,y,t).

@@ Línea 29: / Línea 29: @@
 # Calcular la <math>\nabla \cdot \vec u</math> en todos los puntos del sólido y dibujarla. ¿Qué puntos tienen mayor divergencia? La divergencia es una medida del cambio de volumen local debido al desplazamiento.
 # Calcular <math>|\nabla \times \vec u|</math> en todos los puntos del sólido y dibujarlo. ¿Qué puntos sufren un mayor rotacional?
-# Definamos <math>\epsilon(\vec u)=(\nabla \vec u + \nabla \vec u^t)/2</math>, la parte simétrica del tensor gradiente de <math>\vec u</math>. En un medio elástico lineal, isótropo y homogéneo los desplazamientos permiten escribir el tensor de tensiones <math>\sigma_{ij}</math> a través de la fórmula:
+# Definamos <math>\epsilon(\vec u)=(\nabla \vec u + \nabla \vec u^t)/2</math>, la parte simétrica del tensor gradiente de <math>\vec u</math>, que se denomina [http://en.wikipedia.org/wiki/Infinitesimal_strain_theory tensor de deformaciones]. En un medio elástico lineal, isótropo y homogéneo los desplazamientos permiten escribir el tensor de tensiones <math>\sigma_{ij}</math> a través de la fórmula:
 <math>
 \sigma_{ij}=\lambda \nabla \cdot \vec u \delta_{ij} + 2\mu \epsilon_{ij},

@@ Línea 33: / Línea 33: @@
 \sigma_{ij}=\lambda \nabla \cdot \vec u \delta_{ij} + 2\mu \epsilon_{ij},
 </math>
-donde <math>\lambda</math> y <math>\mu</math> son los conocidos como coeficientes de Lamé que dependen de las propiedades elásticas de cada material. Tomando <math>\lambda=\mu=1</math>, dibujar las tensiones normales en la dirección que marca el eje <math>\vec i</math>, es decir <math>\vec i \cdot \sigma \cdot \vec i</math> y las tensiones normales en la dirección que marca el eje <math>\vec j</math>, es decir <math>\vec j \cdot \sigma \cdot \vec j</math>.
+donde <math>\lambda</math> y <math>\mu</math> son los conocidos como [http://en.wikipedia.org/wiki/Lam%C3%A9_parameters coeficientes de Lamé] que dependen de las propiedades elásticas de cada material. Tomando <math>\lambda=\mu=1</math>, dibujar las tensiones normales en la dirección que marca el eje <math>\vec i</math>, es decir <math>\vec i \cdot \sigma \cdot \vec i</math> y las tensiones normales en la dirección que marca el eje <math>\vec j</math>, es decir <math>\vec j \cdot \sigma \cdot \vec j</math>.
 # Dibujar las tensiones tangenciales respecto al plano ortogonal a <math>\vec i</math>, es decir <math>|\sigma \cdot \vec i-(\vec i \cdot \sigma \cdot \vec i) \vec i|</math>. ¿Donde son mayores? Compararlas con los puntos de mayor deformación de la malla.
 # Dibujar las tensiones tangenciales respecto al plano ortogonal a <math>\vec j</math>, es decir <math>|\sigma \cdot \vec j-(\vec j \cdot \sigma \cdot \vec j) \vec j|</math>. ¿Donde son mayores? Compararla con los puntos de mayor deformación de la malla.
 [[Categoría:Teoría de Campos]]

@@ Línea 25: / Línea 25: @@
 # La temperatura del sólido viene dada por <math>T(\rho,\theta)=e^{-y}</math>. Dibujarla.
 # Calcular <math>\nabla T</math> y dibujarlo como campo vectorial. Calcular las curvas de nivel de T y observar gráficamente que <math>\nabla T</math> es ortogonal a dichas curvas.
-# Consideramos ahora el campo de vectores <math>\vec u=\frac{\sin(y)}{10}\vec i</math>. Dibujar el campo de vectores en los puntos del mallado del sólido.
+# Consideramos ahora el campo de vectores <math>\vec u=\frac{\sin(y)}{10}\vec j</math>. Dibujar el campo de vectores en los puntos del mallado del sólido.
 # Si <math>\vec u</math> determina el desplazamiento que sufre cada punto del sólido, dibujar el sólido antes y después del desplazamiento. Dibujar ambos en la misma figura usando el comando subplot.
 # Calcular la <math>\nabla \cdot \vec u</math> en todos los puntos del sólido y dibujarla. ¿Qué puntos tienen mayor divergencia? La divergencia es una medida del cambio de volumen local debido al desplazamiento.