Diferencia entre revisiones de «Deformaciones de una columna.Grupo 28»

De MateWiki
Saltar a: navegación, buscar

Deprecated: The each() function is deprecated. This message will be suppressed on further calls in /home/mat/public_html/w/includes/diff/DairikiDiff.php on line 434
(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:12.0; %DOMINIO DE Y
+
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=3*log(1+(X+1).^2)+log(1+(Y-2).^2); %TEMPERATURA EN CADA PUNTO 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
 
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,12.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,12.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,12.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=
clc
+
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=3*log(1+(X+1).^2)+log(1+(Y-2).^2); %TEMPERATURA EN CADA PUNTO 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
 
contour(X,Y,T,11); %CURVAS DE NIVEL
dx=(6.*(X+1))./(1+(X+1).^2); %PARCIAL DE X
+
dx=2*pi*X.*cos(2.*pi.*sqrt(X.^2+Y.*2))./sqrt(X.^2+Y.*2); %PARCIAL DE X
dy=(2.*(Y-2))./(1+(Y-2).^2); %PARCIAL DE 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
  • Mateo Navarro Díaz
  • María Victoria González Junco
  • Fernando Benítez Pérez
  • Rodrigo Prado Fornos
  • Claudia Elimar Manrique
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)
Mallado de la placa








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:

[math] \nabla T=\frac{\partial{T}}{\partial{x}}\vec{i}+\frac{\partial{T}}{\partial{y}}\vec{j}[/math]
[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]

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

C.Temperatura C.Nivel
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.

Figura 2
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