Diferencia entre revisiones de «Usuario:Ve.cedillo»

De MateWiki
Saltar a: navegación, buscar
(Página blanqueada)
 
(No se muestran 173 ediciones intermedias de 2 usuarios)
Línea 1: Línea 1:
{{ TrabajoED | Deformaciones de un semianillo circular en 2-D. Grupo 8-A | [[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:TC20/21|2020-21]] | Victoria Elena Cedillo García, María Cristina Pérez-Pozuelo López, Ángela Masiel Zaragoza Paredes }}
 
  
 
En el siguiente artículo procederemos a analizar las deformaciones que sufre la placa bajo la acción de la temperatura y las fuerzas.
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
== 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.
 
 
[[Archivo:A8-1.jpg|350px|miniaturadeimagen|izquierda|Semianillo circular plano]]
 
 
{{matlab|codigo=
 
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])
 
}}
 
 
 
== 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>
 
 
 
[[Archivo:A8-2.jpg|450px|miniatura|derecha|Distribución y curvas de nivel de la temperatura]]
 
{{matlab|codigo=
 
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.
 
 
== 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>
 
[[Archivo:A8-33.jpg|350px|miniaturadeimagen|izquierda|Campo vectorial del gradiente ]]
 
[[Archivo:A8-44.png|350px|miniaturadeimagen|izquierda|Detalle de la ortogonalidad del gradiente respecto a las curvas de nivel ]]
 
 
 
 
 
 
 
 
{{matlab|codigo=
 
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
 
}}
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
[[Categoría:Teoría de Campos]]
 
[[Categoría:TC20/21]]
 

Revisión actual del 20:02 4 dic 2020