Diferencia entre revisiones de «Usuario:Ve.cedillo»
| Línea 149: | Línea 149: | ||
<br /> | <br /> | ||
σ(x,y)=λ∇.u.1 + 2με<br /> | σ(x,y)=λ∇.u.1 + 2με<br /> | ||
| − | Siendo ε la parte simétrica del | + | Siendo ε la parte simétrica del ∇\(\vec u\) y λ y μ los coeficientes de Lamé que varían en función del material y que en este caso son ambos de valor 1. Las tensiones pueden producirse en cualquier dirección del espacio, pero puesto que, como se observa en la matriz el término asociado a \(\vec g_{z}\) es 0 no existen tensiones en dicha dirección. |
Revisión del 10:33 2 dic 2020
| Trabajo realizado por estudiantes | |
|---|---|
| Título | Deformaciones de un semianillo circular en 2-D. Grupo 8-A |
| Asignatura | Teoría de Campos |
| Curso | 2020-21 |
| Autores | Victoria Elena Cedillo García, María Cristina Pérez-Pozuelo López, Ángela Masiel Zaragoza Paredes |
| Este artículo ha sido escrito por estudiantes como parte de su evaluación en la asignatura | |
En el siguiente artículo procederemos a analizar las deformaciones que sufre la placa bajo la acción de la temperatura y las fuerzas.
Contenido
1 Visualización de la placa
Para comenzar se representará los puntos interiores de una placa plana con forma de semianillo circular con radios 1 y 2.
h= 0.1; %Paso de muestreo
%Usamos coordenadas cilíndricas
r= 1:h:2;
tetha= 0:h:pi;
[rr,tt]= meshgrid(r,tetha); %Mallado
%Parametrizamos en cartesianas
x=rr.*cos(tt);
y=rr.*sin(tt);
clf %borramos las posibles gráficas que hubiera
mesh(x,y,0*x); %Visualización de la placa
view(2)
axis ([-3,3,-1,3])
2 Distribución de temperaturas del sólido
En el siguiente apartado se representará la temperatura del sólido mediante tintas hipsométricas y se deducirá el punto en el que esta es máxima. La distribución de temperaturas en cada punto del sólido viene dada por el campo escalar: [math]\ (T(x,y)=exp({-x^2+y^2-1})[/math]
h= 0.1; %paso de muestreo
%usamos coordenadas polares
r= 1:h:2;
tetha= 0:h:pi;
[rr,tt]= meshgrid(r,tetha); %mallado
%parametrizamos en cartesianas
x=rr.*cos(tt);
y=rr.*sin(tt);
clf
%campo temperatura (en cartesianas)
T=exp(-x.^2+y.^2-1); %campo escalar de temperatura
subplot(1,2,1) %Dividimos la pantalla en dos
surf(x,y,T) %representamos el campo escalar de temperaturas
view(2)
axis ([-3,3,-1,3])
colorbar %mostramos la escala
subplot(1,2,2) %escribimos en la segunda pantalla
contour(x,y,T,60) %lineas de nivel
colorbar %mostramos las escala
axis ([-3,3,-1,3])Como se puede observar en la primera gráfica, la temperatura alcanza su apogeo en el punto (0,2), donde el color es más cálido. Esto también se puede apreciar en la segunda gráfica ya que hay mayor proximidad entre las curvas.
3 Estudio del gradiente de temperaturas
A continuación, se estudiará la dirección de máximo crecimiento de la temperatura, es decir, el gradiente, comprobando que este sea ortogonal a las curvas de nivel del apartado anterior. [math]grad T = -2x.exp(-x^2+y^2-1)\vec i\ + 2y.exp(-x^2+y^2-1)\vec j\ [/math]
h= 0.1; %paso de muestreo
%usamos coordenadas polares
r= 1:h:2;
tetha= 0:h:pi;
[rr,tt]= meshgrid(r,tetha); %mallado
%parametrizamos en cartesianas
x=rr.*cos(tt);
y=rr.*sin(tt);
T=exp(-x.^2+y.^2-1); %campo escalar de temperatura
figure
hold on
contour(x,y,T,40) %lineas de nivel
colorbar % mostramos las escala
axis ([-3,3,-1,3])
%gradiente temperatura (en cartesianas)
tx=inline('-2*x.*exp(-x.^2+y.^2-1)','x','y'); %derivada parcial respecto de x
ty=inline('2*y.*exp(-x.^2+y.^2-1)','x','y'); %derivada parcial respecto de y
TX=tx(x,y);
TY=ty(x,y);
h= quiver(x,y,TX,TY); %representamos el campo vectorial
axis ([-3,3,-1,3])
set(h,'maxheadsize',0.33) %cambiamos formato flechas
4 9 Tensiones normales
A continuación, analizaremos las tensiones a las que se ve sometida la placa, estas vienen definidas, en un medio elástico lineal, isótropo y homogéneo, por el tensor:
σ(x,y)=λ∇.u.1 + 2με
Siendo ε la parte simétrica del ∇\(\vec u\) y λ y μ los coeficientes de Lamé que varían en función del material y que en este caso son ambos de valor 1. Las tensiones pueden producirse en cualquier dirección del espacio, pero puesto que, como se observa en la matriz el término asociado a \(\vec g_{z}\) es 0 no existen tensiones en dicha dirección.
