Diferencia entre revisiones de «Deformaciones de un anillo circular»
De MateWiki
(→Lineas de nivel de la temperatura) |
|||
| Línea 8: | Línea 8: | ||
=Mallado= | =Mallado= | ||
| − | Realizaremos en Matlab un mallado de los puntos de la placa sobre la que vamos a trabajar | + | Realizaremos en Matlab un mallado de los puntos de la placa sobre la que vamos a trabajar. Dibujaremos la curva en los ejes <math>(x,y) ∈ [-3,3]\times [-1,3]</math> con un paso de muestro <math>h=1/10</math>. dado que nuestro solido es medio anillo circular debemos proceder usando coordenadas cilindricas. |
{{matlab|codigo= | {{matlab|codigo= | ||
h=0.1; %se define h para dividir los intervalos | h=0.1; %se define h para dividir los intervalos | ||
Revisión del 23:20 3 dic 2021
| Trabajo realizado por estudiantes | |
|---|---|
| Título | Deformaciones de un anillo circular. Grupo C-19 |
| Asignatura | Teoría de Campos |
| Curso | 2021-22 |
| Autores | Sandra Poza Diez
Eduardo Martinez Marinez Jaime Santi Alonso |
| Este artículo ha sido escrito por estudiantes como parte de su evaluación en la asignatura | |
En este articulo
1 Mallado
Realizaremos en Matlab un mallado de los puntos de la placa sobre la que vamos a trabajar. Dibujaremos la curva en los ejes [math](x,y) ∈ [-3,3]\times [-1,3][/math] con un paso de muestro [math]h=1/10[/math]. dado que nuestro solido es medio anillo circular debemos proceder usando coordenadas cilindricas.
h=0.1; %se define h para dividir los intervalos
u=1:h:2; %intervalo rho y teta
v=0:h:pi;
[uu,vv]=meshgrid(u,v); %malla de u x v
xx=uu.*cos(vv); %parametrización con los puntos de la malla
yy=uu.*sin(vv);
zz=0.*uu;
figure(1)
mesh(xx,yy,0*xx); %dibujamos la grafica
axis([-3,3,-1,3]); %región de la grafica
2 Lineas de nivel de la temperatura
La temperatura nos viene dada por el campo [math]T=[/math] que viene en coordenadas cartesianas, como nosostros estamos trabajando en coordenadas cilindricas debemos hacer un cambio, resultando nuestro campo Como se puede apreciar en la gráfica, el punto en el que la temperatura es máxima será (-2,0)
T=log(((xx-3).^2)+2); %Campo escalar temperatura
figure(2)
hold on
contour(xx,yy,T,125) %contour para las superficies de nivel;
colorbar
hold off
