Diferencia entre revisiones de «Deformaciones de una columna.Grupo 28»
(→Gradiente de la temperatura) |
(→Gradiente de la temperatura) |
||
| Línea 47: | Línea 47: | ||
<center><math> \nabla T=\frac{2πx cos(2π\sqrt{x^2 + y^2})}{\sqrt{x^2 + y^2}}\vec{i}+\frac{2πy cos(2π\sqrt{x^2 + y^2})}{\sqrt{x^2 + y^2}}\vec{j}</math></center> | <center><math> \nabla T=\frac{2πx cos(2π\sqrt{x^2 + y^2})}{\sqrt{x^2 + y^2}}\vec{i}+\frac{2πy cos(2π\sqrt{x^2 + y^2})}{\sqrt{x^2 + y^2}}\vec{j}</math></center> | ||
| − | A continuación se muestra el código completo de matlab, del cual resultan las siguientes imágenes en las que podemos ver las curvas de nivel y donde se encuentra su máximo y mínimo, a partir de los colores de las gráficas y a partir de el propio matlab. Este nos índica al ejecutar el programa que el máximo es: | + | A continuación se muestra el código completo de matlab, del cual resultan las siguientes imágenes en las que podemos ver las curvas de nivel y donde se encuentra su máximo y mínimo, a partir de los colores de las gráficas y a partir de el propio matlab. Este nos índica al ejecutar el programa que el máximo es: 0.8 |
[[Archivo:Ctempcnivel.jpg|1120 × 840px|miniaturadeimagen|C.Temperatura C.Nivel]] | [[Archivo:Ctempcnivel.jpg|1120 × 840px|miniaturadeimagen|C.Temperatura C.Nivel]] | ||
| Línea 55: | Línea 55: | ||
h = 0.2; %MUESTREO | h = 0.2; %MUESTREO | ||
x =-1.0:h:1.0; %DOMINIO DE X | x =-1.0:h:1.0; %DOMINIO DE X | ||
| − | y =0.0:h: | + | y =0.0:h:10.0; %DOMINIO DE Y |
[X,Y]= meshgrid(x,y); % CREACIÓN DEL MALLADO | [X,Y]= meshgrid(x,y); % CREACIÓN DEL MALLADO | ||
| − | T= | + | T=sin(2.*pi.*(sqrt(X.^2+Y.^2))); %TEMPERATURA EN CADA PUNTO DEL MALLADO |
mesh(X,Y,0*X) %CREACIÓN DE LA MALLA | mesh(X,Y,0*X) %CREACIÓN DE LA MALLA | ||
subplot(1,3,1) %GRÁFICO SUPERIOR | subplot(1,3,1) %GRÁFICO SUPERIOR | ||
surf(X,Y,T); %DEGRADACIÓN DE COLORES | surf(X,Y,T); %DEGRADACIÓN DE COLORES | ||
| − | axis ([-1,1,-0.5, | + | axis ([-1,1,-0.5,10.5]) |
view(2) | view(2) | ||
colorbar %BARRA DE COLORES | colorbar %BARRA DE COLORES | ||
| Línea 70: | Línea 70: | ||
subplot(1,3,2) %GRÁFICO INFERIOR | subplot(1,3,2) %GRÁFICO INFERIOR | ||
contour(X,Y,T,11); | contour(X,Y,T,11); | ||
| − | axis ([-1.0,1.0,-0.5, | + | axis ([-1.0,1.0,-0.5,10.5]); |
%TÍTULO SEGUNDA GRÁFICA Y EJES | %TÍTULO SEGUNDA GRÁFICA Y EJES | ||
title('CURVAS DE NIVEL') | title('CURVAS DE NIVEL') | ||
| Línea 78: | Línea 78: | ||
Tmax=max(max(T)) %PUNTO MÁXIMO | Tmax=max(max(T)) %PUNTO MÁXIMO | ||
subplot(1,3,3) | subplot(1,3,3) | ||
| − | axis ([-1.25,1.25,-0.5, | + | axis ([-1.25,1.25,-0.5,10.5]) |
view(2) | view(2) | ||
%CURVAS DE NIVEL ENUMERADAS | %CURVAS DE NIVEL ENUMERADAS | ||
| Línea 93: | Línea 93: | ||
[[Archivo:definitivoo.png|300px|miniaturadeimagen|Figura 2]] | [[Archivo:definitivoo.png|300px|miniaturadeimagen|Figura 2]] | ||
{{matlab|codigo= | {{matlab|codigo= | ||
| − | + | clear;clc; | |
| − | clear | + | |
h = 2/10; %MUESTREO | h = 2/10; %MUESTREO | ||
x = [-1:h:1]; %DOMINIO DE X | x = [-1:h:1]; %DOMINIO DE X | ||
y = [0:h:12]; %DOMINIO DE Y | y = [0:h:12]; %DOMINIO DE Y | ||
[X,Y]= meshgrid(x,y); %CREACIÓN DEL MALLADO | [X,Y]= meshgrid(x,y); %CREACIÓN DEL MALLADO | ||
| − | T= | + | T=sin(2.*pi.*(sqrt(X.^2+Y.^2))); %TEMPERATURA EN CADA PUNTO DEL MALLADO |
contour(X,Y,T,11); %CURVAS DE NIVEL | contour(X,Y,T,11); %CURVAS DE NIVEL | ||
| − | dx=( | + | dx=2*pi*X.*cos(2.*pi.*sqrt(X.^2+Y.*2))./sqrt(X.^2+Y.*2); %PARCIAL DE X |
| − | dy=(2.*(Y | + | dy=2*pi*Y.*cos(2.*pi.*sqrt(X.^2+Y.*2))./sqrt(X.^2+Y.*2); %PARCIAL DE Y |
%TÍTULO Y EJES | %TÍTULO Y EJES | ||
title('Gradiente de temperatura'); | title('Gradiente de temperatura'); | ||
Revisión del 19:58 28 nov 2024
| Trabajo realizado por estudiantes | |
|---|---|
| Título | Deformaciones de una columna . Grupo 28 |
| Asignatura | Teoría de Campos |
| Curso | 2024-25 |
| Autores |
|
| Este artículo ha sido escrito por estudiantes como parte de su evaluación en la asignatura | |
1 Representación del mallado del sólido
Definimos la placa mediante dos vectores x e y que representan el intervalo en el cual está definida nuestra placa, tomando un muestreo h=1/10. Al crear la malla donde la placa está definida, al representarla con el comando surf, tenemos que tener en cuenta que es plana dando valor 0 a la altura de los puntos de la misma.
h=0.1; %muestreo
x=-1:h:1; %eje x de la placa
y=0:h:10; %eje y de la placa
[xx,yy]=meshgrid(x,y); %mallado de la placa
figure(1)
surf(xx,yy,0*xx); %representación de la placa
axis([-2,2,0,10]) %ejes del rectángulo
view(2)
2 Gradiente de la temperatura
El cálculo del gradiente de un campo escalar se expresa como la derivada parcial respecto de cada coordenada de dicho campo en función de la base ortonormal orientada positiva [math]\vec{i},\vec{j},\vec{k}[/math] (COORDENADAS CARTESIANAS). Es decir:
A continuación se muestra el código completo de matlab, del cual resultan las siguientes imágenes en las que podemos ver las curvas de nivel y donde se encuentra su máximo y mínimo, a partir de los colores de las gráficas y a partir de el propio matlab. Este nos índica al ejecutar el programa que el máximo es: 0.8
clc
clear
h = 0.2; %MUESTREO
x =-1.0:h:1.0; %DOMINIO DE X
y =0.0:h:10.0; %DOMINIO DE Y
[X,Y]= meshgrid(x,y); % CREACIÓN DEL MALLADO
T=sin(2.*pi.*(sqrt(X.^2+Y.^2))); %TEMPERATURA EN CADA PUNTO DEL MALLADO
mesh(X,Y,0*X) %CREACIÓN DE LA MALLA
subplot(1,3,1) %GRÁFICO SUPERIOR
surf(X,Y,T); %DEGRADACIÓN DE COLORES
axis ([-1,1,-0.5,10.5])
view(2)
colorbar %BARRA DE COLORES
%TÍTULO PRIMERA GRÁFICA Y EJES
title('CAMPO DE TEMPERATURAS')
xlabel('eje X')
ylabel('eje Y')
subplot(1,3,2) %GRÁFICO INFERIOR
contour(X,Y,T,11);
axis ([-1.0,1.0,-0.5,10.5]);
%TÍTULO SEGUNDA GRÁFICA Y EJES
title('CURVAS DE NIVEL')
xlabel('eje X')
ylabel('eje Y')
colorbar %BARRA DE COLORES
Tmax=max(max(T)) %PUNTO MÁXIMO
subplot(1,3,3)
axis ([-1.25,1.25,-0.5,10.5])
view(2)
%CURVAS DE NIVEL ENUMERADAS
[M,c]=contour(X,Y,T,[0:0.75:10],'ShowText','on')
%TÍTULO TERCERA GRÁFICA Y EJES
title('CURVAS DE NIVEL NUMERADAS')
xlabel('eje X')
ylabel('eje Y')
colorbar %BARRA DE COLORES
A continuación, se va a demostrar como el gradiente, previamente calculado, es perpendicular a las líneas de nivel anteriores.
Al no verse claro en la imagen original, con la cantidad de líneas de flujo, se opta por aplicar un zoom para una correcta visualización.
clear;clc;
h = 2/10; %MUESTREO
x = [-1:h:1]; %DOMINIO DE X
y = [0:h:12]; %DOMINIO DE Y
[X,Y]= meshgrid(x,y); %CREACIÓN DEL MALLADO
T=sin(2.*pi.*(sqrt(X.^2+Y.^2))); %TEMPERATURA EN CADA PUNTO DEL MALLADO
contour(X,Y,T,11); %CURVAS DE NIVEL
dx=2*pi*X.*cos(2.*pi.*sqrt(X.^2+Y.*2))./sqrt(X.^2+Y.*2); %PARCIAL DE X
dy=2*pi*Y.*cos(2.*pi.*sqrt(X.^2+Y.*2))./sqrt(X.^2+Y.*2); %PARCIAL DE Y
%TÍTULO Y EJES
title('Gradiente de temperatura');
xlabel('Eje X');
ylabel('Eje Y');
% Representación de la temperatura y las curvas de nivel
hold on
quiver(x,y,dx,dy);
axis equal
colorbar