Diferencia entre revisiones de «Usuario:Ve.cedillo»

De MateWiki
Saltar a: navegación, buscar
(→‎Efecto de los desplazamientos)
(→‎Representación campo de vectores)
Línea 178: Línea 178:
 
set(w,'maxheadsize',0.33)        %cambiamos formato flechas
 
set(w,'maxheadsize',0.33)        %cambiamos formato flechas
 
}}
 
}}
 +
Según nuestro gráfico, el campo de vectores es mayor en la zona central del semianillo, donde las flechas son más largas y hay más próximas.
  
 
==Efecto de los desplazamientos==
 
==Efecto de los desplazamientos==

Revisión del 06:35 4 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.









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.

Semianillo circular plano
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]


Distribución y curvas de nivel de la temperatura
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 
                                  %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)                 %Líneas 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 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);
T=exp(-x.^2+y.^2-1);             %Campo escalar de temperatura
figure
hold on
contour(x,y,T,40)                %Líneas 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


Campo vectorial del gradiente
Detalle de la ortogonalidad del gradiente respecto a las curvas de nivel
















4 Campo vectorial en los puntos del sólido

Conociendo el campo de deformaciones de la placa: 

                                                  u(ρ,θ) = sinθf(ρ)\(\vec g_{ρ}\) + 0\(\vec g_{\theta}\) + 0\(\vec g_{z}\)

Tenemos dos condiciones que se deben cumplir:

1º Que los puntos situados en ρ = 1 no sufren desplazamiento, lo que significa que: u(ρ,θ) = sin(θ)f(1)\(\vec g_{ρ}\) = 0 ⇔ f(1) = 0


2º |∇ · u| = sinθ(2ρ − 1)/4


Ahora calculamos la divergencia del campo: ∇ · u = 1/ρ{∂/∂ρ (ρ · sin(θ)f(ρ))} = (sin(θ)(ρf(ρ)))'/ρ = sinθ(2ρ − 1)/4 (ρf(ρ))' = (2ρ^2-ρ)/4



5 Representación campo de vectores

A continuación, se representa el campo de desplazamientos calculado en el apartado anterior. Para ello se empleará la base ortonormal cartesiana por lo que se pasa el \(\vec g_{\rho}\) a [math] cosθi + senθ j [/math], obteniendo: [math] u(ρ,θ)=senθ. cosθ. (2ρ-1)/4 i + sen(θ)^2 .(2ρ-1)/4[/math]

izquierda
h= 0.1;                           %paso de muestreo
                                  %usamos coordenadas polares.
rr= 1:h:2;
tt= 0:h:pi;
[RR,TT]= meshgrid(rr,tt);         %mallado
x=RR.*cos(TT);
y=RR.*sin(TT);                     
axis ([-3,3,-1,3])
a= sin(TT).*cos(TT).*(2.*RR-1)/4;
b=(sin(TT)).^2.*(2.*RR-1)/4;
w=quiver(x,y,a,b);
axis ([-3,3,-1,3])
set(w,'maxheadsize',0.33)        %cambiamos formato flechas

Según nuestro gráfico, el campo de vectores es mayor en la zona central del semianillo, donde las flechas son más largas y hay más próximas.

6 Efecto de los desplazamientos

Para poder observar cómo afecta el campo de desplazamientos a nuestro sólido, lo representaremos antes y después de su efecto.


Efecto de los desplazamientos.
h= 0.1;                          %Paso de muestreo                                 
rr= 1:h:2;                       %Usamos coordenadas polares
tt= 0:h:pi;
[RR,TT]= meshgrid(rr,tt);        %Mallado
x=RR.*cos(TT);
y=RR.*sin(TT);
clf
%Sólido antes de los desplazamientos
subplot(1,2,1)                  
i=mesh(x,y,0*x);
view(2)
set(i,'EdgeColor','g');
axis ([-3,3,-1,3])
title('Placa no desplazada','Fontsize',25); 

%Sólido después de los desplazamientos
subplot(1,2,2)                     
A= sin(TT).*cos(TT).*(2.*RR-1)/4;
B=(sin(TT)).^2.*(2.*RR-1)/4; 
X=x+A;
Y=y+B;
mesh(X,Y,0*X)
view(2)
axis ([-3,3,-1,3])
title('Placa desplazada','Fontsize',25);

%Comparación de ambas representaciones. 
figure                                 
mesh(X,Y,0*X)
axis ([-3,3,-1,3])
view(2)
title('Desplazamiento de la placa','Fontsize',16);
hold on
i = mesh(x,y,0*x);
set(i,'EdgeColor','g');


Como podemos ver, el desplazamiento sigue la dirección del campo de vectores hallado anteriormente, deformándose más en la zona central del semianillo.

Comparación del antes y después del efecto.

7 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.

Matriz de componentes de sigma

Pese a que los desplazamientos se realizan en el plano, las tensiones no tienen por qué ser planas, y se puede representar en las tres direcciones del espacio. Puesto que estamos empleando coordenadas cilíndricas, se representan las tensiones en las direcciones de la base ortogonal cilíndrica {\(\vec g_{ρ}\), \(\vec g_{\theta}\), \(\vec g_{z}\)}. Para ello, multiplicamos escalarmente la matriz sigma por el vector \(\vec g_{i}\) por ambos lados: tensión normales en dirección i = \(\vec g_{i}\) * σ * \(\vec g_{i}\)

a) dirección \(\vec g_{\rho}\): (1,0,0)*σ*(1,0,0) = [math](7ρ/6-1/2+1/(12ρ^2 ))sinθ [/math]

b) dirección \(\vec g_{\theta}\)/ρ: (0,1,0)*σ*(0,1,0) = [math] (5ρ/6-1/2-1/(12ρ^2 ))sinθ[/math]

c) dirección \(\vec g_{z}\): (0,0,1)*σ*(0,0,1) = [math](ρ/2-1/4)sinθ [/math]

Representamos las tensiones primero en 2D y, posteriormente en 3D en las dos direcciones en las que existen tensiones.

Representación de las tensiones en las direcciones rho y sigma
h=0.1;
r=1:h:2;
t=0:h:pi;
[rr,tt]=meshgrid(r,t); 
clf
subplot(1,3,1);
xx=rr.*cos(tt);
yy=rr.*sin(tt);
a=(7*rr./6-1/2+1./(12*rr.^2)).*sin(tt); %Elemento (1,1) de la matriz sigma
b=(5*rr/6-1/2-1./(12*rr.^2)).*sin(tt);
surf(xx,yy,a)
axis equal
view(2)
colorbar
subplot(1,3,2);   %Elemento (2,2) de la matriz sigma
surf(xx,yy,b)
axis equal
view(2)
colorbar
subplot(1,3,3)
c=(rr/2-1/4).*sin(rr);            %Elemento (3,3) de la matriz sigma
surf(xx,yy,c)
axis equal
view(2)
colorbar



izquierda
h=0.1;
  r=1:h:2;
  t=0:h:pi;
 [rr,tt]=meshgrid(r,t); 
 clf
 subplot(3,1,1);
 xx=rr.*cos(tt);
 yy=rr.*sin(tt);
 a=(7*rr./6-1/2+1./(12*rr.^2)).*sin(tt); %Elemento (1,1) de la matriz sigma
 b=(5*rr/6-1/2-1./(12*rr.^2)).*sin(tt);
 c=(rr/2-1/4).*sin(rr); 
 surf(xx,yy,a)
 colorbar
 subplot(3,1,2);   %Elemento (2,2) de la matriz sigma
 surf(xx,yy,b)
 colorbar
 subplot(3,1,3);   %Elemento (3,3) de la matriz sigma
 surf(xx,yy,c)
 colorbar


8 10 Tensiones tangenciales

Una vez calculadas las tensiones normales podemos pasar a analizar las tensiones tangenciales, estas vienen dadas por σ.*\(\vec g_{\rho}\)+(\(\vec g_{\rho}\)*σ*\(\vec g_{\rho}\))\(\vec g_{\rho}\), y como \(\vec g_{\rho}\) tiene como coordenadas cilíndricas (1,0,0), esto da [math](ρ/6-1/8-1/(24ρ^2 ))cosθ[/math].

Tensiones tangenciales en la dirección de sigma
h=0.1;
r=1:h:2;
t=0:h:pi;
[rr,tt]=meshgrid(r,t);
xx=rr.*cos(tt);
yy=rr.*sin(tt);
a=(rr/6-1/8-1./(24*rr.^2 )).*cos(tt);
subplot(2,1,1)
surf(xx,yy,a)
colorbar
subplot(2,1,2)
surf(xx,yy,a)
colorbar
view(2)





9 Tensión de Von Mises

Conociendo que la tensión de Von Mises se define como:

Fórmula de Von Mises

siendo \(\ σ_{i}\) con i=1,2,3 los autovalores de σ (también son conocidos como tensiones principales)

La tensión de Von Mises es una magnitud escalar que es un indicador para saber cuando un material inicia un comportamiento plástico y no elástico puro

Tensión de Von Mises
h=0.1;
r=[1:h:2];
t=[0:h:pi];
[rr,tt]=meshgrid(r,t);
%Creación de la matriz sigma
s=zeros(3,3);
%Creamos la matriz de von mises 
vm=zeros(32,11); %Tiene un tamaño 32x11 porque es el tamaño del mallado
for i=1:32*11
    rrr=rr(i)';
    ttt=tt(i)';
    s(1,1)=(7*rrr./6-1/2+1./(12*rrr.^2)).*sin(ttt);
    s(1,2)=(rrr/6-1/8-1./(24*rrr.^2 )).*cos(ttt);
    s(2,1)=s(1,2);
    s(2,2)=(5*rrr/6-1/2-1./(12*rrr.^2)).*sin(ttt);
    s(3,3)=(rrr/2-1/4).*sin(rrr);
    [v,d]=eig(s);                                      %Autovalores de la matriz sigma
    vm(i)=sqrt(((d(1,1)-d(2,2))^2+(d(2,2)-d(3,3))^2+(d(3,3)-d(1,1))^2)/2);  %Fórmula de Von Mises
end
%Cambiamos a cartesianas
xx=rr.*cos(tt);
yy=rr.*sin(tt);
surf(xx,yy,vm)
colorbar
view(2)


Como se puede apreciar en la gráfica de la tensión, el punto donde ésta alcanza el mayor valor es en las inmediaciones del punto (0,2)

10 Masa total de la placa

La densidad de la placa viene dada en coordenadas cartesianas
[math] d(x, y) = 1 + x*y*log(1 + x + y^2) [/math]

Obtenido de «https://mat.caminos.upm.es/w/index.php?title=Usuario:Ve.cedillo&oldid=45136»