Diferencia entre revisiones de «Campos escalares en una placa plana»

De MateWiki
Saltar a: navegación, buscar
 
(No se muestran 26 ediciones intermedias de 2 usuarios)
Línea 9: Línea 9:
 
Representa la sección transversal de un sólido para el cual se desprecian las variaciones en la dirección ortogonal a la sección mencionada.
 
Representa la sección transversal de un sólido para el cual se desprecian las variaciones en la dirección ortogonal a la sección mencionada.
  
−
Suponemos que tenemos dos cantidades físicas definidas en nuestra sección: la temperatura T(x,y)=x<sup>2</sup>+(y-1)<sup>2</sup>, y los desplazamientos, producidos por una fuerza u(ρ,θ)=ρ/20 *sin(6(θ-π/4)) e<sub>ρ</sub> .
+
Suponemos que tenemos dos cantidades físicas definidas en nuestra sección: la temperatura T(x,y)=x<sup>2</sup>+(y-1)<sup>2</sup>, y los desplazamientos, producidos por una fuerza <math>\vec{u}</math>(ρ,θ)=ρ/20 *sin(6(θ-π/4)) <math>\vec{e}</math><sub>ρ</sub> .
  
−
El vector posición de la placa antes de la deformación es: r<sub>o</sub>(x,y)=x i +y j ,y el vector posición de la placa después de la deformación es: r<sub>d</sub>(x,y)=r<sub>o</sub>(x,y)+u(x,y) .
+
El vector posición de la placa antes de la deformación es: <math>\vec{r}</math><sub>o</sub>(x,y)=x <math>\vec{i}</math> +y <math>\vec{j}</math> ,y el vector posición de la placa después de la deformación es:  
 +
 
 +
<math>\vec{r}</math><sub>d</sub>(x,y) = <math>\vec{r}</math><sub>o</sub>(x,y)+<math>\vec{u}</math>(x,y) .
  
 
== .-Dibujo del mallado ==
 
== .-Dibujo del mallado ==
Línea 61: Línea 63:
 
% Aplicamos la función de la temperatura al mallado
 
% Aplicamos la función de la temperatura al mallado
 
surf(Mx,My,T);
 
surf(Mx,My,T);
−
%Establecemos el límite de los ejes
+
% Establecemos el límite de los ejes
 
axis([-3,3,-1,3]);
 
axis([-3,3,-1,3]);
 
view(2)
 
view(2)
Línea 76: Línea 78:
  
 
== .-Gradiente de la temperatura==
 
== .-Gradiente de la temperatura==
−
Vamos a calcular y a representar el gradiente de la temperatura. Se define el gradiente como: ∇T(x,y)=∂T/∂x i +∂T/∂y j .
+
Vamos a calcular y a representar el gradiente de la temperatura.  
−
 
+
Utilizando nuestra función de la temperatura, el gradiente es:  
−
Utilizando nuestra función de la temperatura, el gradiente es: ∇T=2x i + 2(y-1) j .
+
[[Archivo:Eje-3.png|miniaturadeimagen|centro]]
−
 
+
 
Con la gráfica, podemos observar que el gradiente es ortogonal a las curvas de nivel de la temperatura. Además, sabemos que el gradiente de cualquier función es la variación máxima entre varios puntos.
 
Con la gráfica, podemos observar que el gradiente es ortogonal a las curvas de nivel de la temperatura. Además, sabemos que el gradiente de cualquier función es la variación máxima entre varios puntos.
  
Línea 86: Línea 87:
 
clc
 
clc
 
clear
 
clear
−
%Definimos regiones
+
% Definimos regiones
 
h = 0.1;
 
h = 0.1;
 
x=[1:h:2]; 
 
x=[1:h:2]; 
 
y=[pi/4:h:3*pi/4];
 
y=[pi/4:h:3*pi/4];
−
%Creamos el mallado
+
% Creamos el mallado
 
[rr,tt]= meshgrid(x,y);
 
[rr,tt]= meshgrid(x,y);
 
Mx=rr.*cos(tt);
 
Mx=rr.*cos(tt);
 
My=rr.*sin(tt);
 
My=rr.*sin(tt);
−
%Función temperatura
+
% Función temperatura
 
T = Mx.^2 + (My-1).^2;
 
T = Mx.^2 + (My-1).^2;
 
contour(Mx,My,T,30);
 
contour(Mx,My,T,30);
−
%Gradiente de T
+
% Gradiente de T
 
dx = 2*Mx;
 
dx = 2*Mx;
 
dy = 2*(My-1);
 
dy = 2*(My-1);
−
%Título y ejes
+
% Título y ejes
 
title('Gradiente de temperatura');
 
title('Gradiente de temperatura');
 
xlabel('EJE X');
 
xlabel('EJE X');
 
ylabel('EJE Y');
 
ylabel('EJE Y');
−
%Representación de la temperatura y las curvas de nivel
+
% Representación de la temperatura y las curvas de nivel
 
hold on
 
hold on
 
quiver(Mx,My,dx,dy);
 
quiver(Mx,My,dx,dy);
Línea 112: Línea 113:
  
 
== .-Campo de vectores de deslizamiento==
 
== .-Campo de vectores de deslizamiento==
−
Representamos el campo de vectores de deslizamiento en la placa plana. El campo dado en el enunciado está en coordenadas cilíndricas y hacemos un cambio a coordenadas cartesianas, ya que conocemos que el vector e<sub>ρ</sub> = cosθ i + sinθ j.
+
Representamos el campo de vectores de deslizamiento en la placa plana. El campo dado en el enunciado está en coordenadas cilíndricas y hacemos un cambio a coordenadas cartesianas, ya que conocemos que el vector <math>\vec{e}</math><sub>ρ</sub> = cosθ <math>\vec{i}</math> + sinθ <math>\vec{j}</math>.
  
−
El campo de vectores de deslizamiento en coordenadas cartesianas es: u(ρ,θ,z)= ρ/20 *sin(6(θ-π/4))cosθ i + ρ/20 *sin(6(θ-π/4))sinθ j.
+
El campo de vectores de deslizamiento en coordenadas cartesianas es:  
 +
[[Archivo:Eje-4.png|miniaturadeimagen|centro]]
  
 
Utilizaremos el comando quiver() para representar el campo de vectores.
 
Utilizaremos el comando quiver() para representar el campo de vectores.
Línea 136: Línea 138:
 
mx=(rr.*sin(6.*(tt-pi/4)).*cos(tt))/20;
 
mx=(rr.*sin(6.*(tt-pi/4)).*cos(tt))/20;
 
my=(rr.*sin(6.*(tt-pi/4)).*sin(tt))/20;
 
my=(rr.*sin(6.*(tt-pi/4)).*sin(tt))/20;
−
%campo vectorial en 2D
+
% Campo vectorial en 2D
 
quiver(Mx,My,mx,my);
 
quiver(Mx,My,mx,my);
 
axis equal
 
axis equal
Línea 143: Línea 145:
  
 
== .-Deslizamiento del sólido==
 
== .-Deslizamiento del sólido==
−
A continuación, vamos a representar el sólido antes y después del desplazamiento dado por el campo de vectores u (en t=0). Realizaremos tres gráficas: situación inicial, situación final y comparación. En nuestro campo, podemos observar que sólo varía en la dirección del vector unitario e<sub>ρ</sub>.
+
A continuación, vamos a representar el sólido antes y después del desplazamiento dado por el campo de vectores <math>\vec{u}</math> (en t=0). Realizaremos tres gráficas: situación inicial, situación final y comparación. En nuestro campo, podemos observar que sólo varía en la dirección del vector unitario <math>\vec{e}</math><sub>ρ</sub>.
  
−
[[Archivo:Ej-5.png|405px|thumb|right|Variación del deslizamiento]]
+
[[Archivo:Ej-5.png|2000px|thumb|right|Variación del deslizamiento]]
 
{{matlab|codigo=
 
{{matlab|codigo=
 
clear
 
clear
Línea 160: Línea 162:
 
view(2)
 
view(2)
 
title('SITUACIÓN INICIAL')  
 
title('SITUACIÓN INICIAL')  
−
%situación final
+
% Situación final
 
subplot(2,2,2)
 
subplot(2,2,2)
 
mx=(rr.*sin(6.*(tt-pi/4)).*cos(tt))/20;
 
mx=(rr.*sin(6.*(tt-pi/4)).*cos(tt))/20;
Línea 182: Línea 184:
  
 
== .-Divergencia del campo vectorial ==
 
== .-Divergencia del campo vectorial ==
−
 
+
La divergencia es una medida del cambio de volumen local debido al deslizamiento.
 +
Calculamos la divergencia del campo vectorial <math>\vec{u}</math>, y además, calculamos analíticamente los puntos en los que la divergencia de u es máxima, mínima y nula.
 +
[[Archivo:Eje-6.png|miniaturadeimagen|centro]]
  
 
[[Archivo:Ej-6.png|405px|thumb|right|Divergencia del campo vectorial]]
 
[[Archivo:Ej-6.png|405px|thumb|right|Divergencia del campo vectorial]]
Línea 209: Línea 213:
  
 
== .-Rotacional del campo ==
 
== .-Rotacional del campo ==
 +
Calculamos el rotacional de campo <math>\vec{u}</math> y deducimos que puntos sufren un mayor rotacional:
 +
[[Archivo:Eje-7.png|miniaturadeimagen|centro]]
  
−
 
+
[[Archivo:Ej-7.png|500px|thumb|left|Rotacional del campo]]
−
[[Archivo:Ej-7.png|405px|thumb|right|Rotacional del campo]]
+
 
[[Archivo:Ej-7.2.png|405px|thumb|right|Rotacional del campo]]
 
[[Archivo:Ej-7.2.png|405px|thumb|right|Rotacional del campo]]
 
{{matlab|codigo=
 
{{matlab|codigo=
Línea 234: Línea 239:
 
shading interp
 
shading interp
 
}}
 
}}
 +
 +
== .-Tensor de tensiones ==
 +
La placa se ve sometida a numerosas tensiones normales.En un medio elástico lineal, isótropo y homogéneo los desplazamientos permiten escribir el tensor de tensiones σ<sub>ij</sub> a través de la fórmula: σ=λ∇ * (<math>\vec{u}</math>)1 + 2μϵ , donde:
 +
* 1 es el tensor identidad en el conjunto de vectores libres del espacio ℝ<sup>3</sup>.
 +
* λ, µ son los conocidos como «coeficientes de Lamé», que dependen de las propiedades elásticas de cada material.
 +
 +
Tomando λ=μ=1 debido a las propiedades de nuestro material, calculamos las tensiones normales en la dirección del eje <math>\vec{e}</math><sub>ρ</sub>,<math>\vec{e}</math><sub>θ</sub>, y <math>\vec{e}</math><sub>z</sub>.
 +
 +
=== Tensor de deformaciones y de tensiones ===
 +
Calculamos el tensor de deformaciones y de tensiones:
 +
[[Archivo:Eje-8.png|miniaturadeimagen|centro]]
 +
[[Archivo:Eje-81.png|miniaturadeimagen|centro]]
 +
 +
=== Tensor de tensiones en la dirección <math>\vec{e}</math><sub>ρ</sub> ===
 +
 +
[[Archivo:Ej-8.png|2000px|thumb|left]]
 +
[[Archivo:Ej-8.2.png|400px|thumb|right]]
 +
{{matlab|codigo=
 +
%tensiones normales al eje e-ro
 +
clear
 +
clc
 +
x=1:0.1:2;
 +
y=pi/4:0.1:(3*pi)/4;
 +
[rr,tt]=meshgrid(x,y);
 +
Mx=rr.*cos(tt);
 +
My=rr.*sin(tt);
 +
%elemento (1,1) de la matriz tensión
 +
s1=(sin(6.*(tt-pi/4)).*(1/5));
 +
surf(Mx,My,s1);
 +
shading interp
 +
view(2)
 +
axis([-3,3,-1,3]);
 +
title('TENSIONES NORMALES AL EJE E DE RO')
 +
xlabel('EJE X');
 +
ylabel('EJE Y');
 +
zlabel('EJE Z');
 +
colorbar
 +
}}
 +
 +
=== Tensor de tensiones en la dirección <math>\vec{e}</math><sub>θ</sub> ===
 +
[[Archivo:Ej-8.3.png|2000px|thumb|left]]
 +
[[Archivo:Ej-8.4.png|350px|thumb|right]]
 +
{{matlab|codigo=
 +
%tensiones normales al eje e-teta
 +
clear
 +
clc
 +
x=1:0.1:2;
 +
y=pi/4:0.1:(3*pi)/4;
 +
[rr,tt]=meshgrid(x,y);
 +
Mx=rr.*cos(tt);
 +
My=rr.*sin(tt);
 +
%elemento (1,1) de la matriz tensión
 +
s2=(sin(6.*(tt-pi/4)).*(1/5));
 +
surf(Mx,My,s2);
 +
shading interp
 +
view(2)
 +
axis([-3,3,-1,3]);
 +
title('TENSIONES NORMALES AL EJE E DE TETA')
 +
xlabel('EJE X');
 +
ylabel('EJE Y');
 +
zlabel('EJE Z');+
 +
colorbar
 +
}}
 +
 +
=== Tensor de tensiones en la dirección <math>\vec{e}</math><sub>z</sub> ===
 +
[[Archivo:Ej-8.6.png|2000px|thumb|left]]
 +
[[Archivo:Ej-8.5.png|350px|thumb|right]]
 +
{{matlab|codigo=
 +
%tensiones normales al eje e-z
 +
clear
 +
clc
 +
x=1:0.1:2;
 +
y=pi/4:0.1:(3*pi)/4;
 +
[rr,tt]=meshgrid(x,y);
 +
Mx=rr.*cos(tt);
 +
My=rr.*sin(tt);
 +
%elemento (1,1) de la matriz tensión
 +
s3=(sin(6.*(tt-pi/4)).*(1/10));
 +
surf(Mx,My,s3);
 +
shading interp
 +
view(2)
 +
axis([-3,3,-1,3]);
 +
title('TENSIONES NORMALES AL EJE E DE Z')
 +
xlabel('EJE X');
 +
ylabel('EJE Y');
 +
zlabel('EJE Z');
 +
colorbar
 +
}}
 +
 +
 +
== .-Tensiones tangenciales ==
 +
Calculamos las tensiones tangenciales respecto al plano ortogonal <math>\vec{e}</math><sub>ρ</sub> :
 +
 +
[[Archivo:Ejercicio_9.png|miniaturadeimagen|centro]]
 +
 +
[[Archivo:Ej-9.2.png|1000px|thumb|left|En 3D]]
 +
[[Archivo:Ej-9.png|350px|thumb|right|En 2D]]
 +
{{matlab|codigo=
 +
clear
 +
clc
 +
%Definición de regiones
 +
x=1:0.1:2;
 +
y=pi/4:0.1:3*pi/4;
 +
[rr,tt]=meshgrid(x,y);
 +
Mx=rr.*cos(tt);
 +
My=rr.*sin(tt);
 +
%Tensión tangencial
 +
tg=abs(6*cos(6*(tt-pi/4))/20);
 +
surf(Mx,My,tg);
 +
shading interp
 +
axis([-3,3,-1,3]);
 +
title('TENSIONES TANGENCIALES RESPECTO AL PLANO ORTOGONAL')
 +
xlabel('EJE X');
 +
ylabel('EJE Y');
 +
zlabel('EJE Z');
 +
colorbar
 +
}}
 +
 +
[[Categoría:Teoría de Campos]]
 +
[[Categoría:TC22/23]]

Revisión actual del 14:51 10 dic 2022

Trabajo realizado por estudiantes
Título Campos escalares en una placa plana. Grupo 2-A
Asignatura Teoría de Campos
Curso 2022-23
Autores María González García-Nieto

Cristina Rubio Yanes

Marta García-Moris Fontcuberta

María Jiménez Estríngana

Este artículo ha sido escrito por estudiantes como parte de su evaluación en la asignatura


Consideramos una placa plana que ocupa un cuarto de un anillo circular centrado en el origen y comprendido entre los radios 1 y 2, en el plano y≥|x|, con los ejes en el cuadrado [-3,3]x[-1,3]. Representa la sección transversal de un sólido para el cual se desprecian las variaciones en la dirección ortogonal a la sección mencionada.

Suponemos que tenemos dos cantidades físicas definidas en nuestra sección: la temperatura T(x,y)=x2+(y-1)2, y los desplazamientos, producidos por una fuerza [math]\vec{u}[/math](ρ,θ)=ρ/20 *sin(6(θ-π/4)) [math]\vec{e}[/math]ρ .

El vector posición de la placa antes de la deformación es: [math]\vec{r}[/math]o(x,y)=x [math]\vec{i}[/math] +y [math]\vec{j}[/math] ,y el vector posición de la placa después de la deformación es:

[math]\vec{r}[/math]d(x,y) = [math]\vec{r}[/math]o(x,y)+[math]\vec{u}[/math](x,y) .

1 .-Dibujo del mallado

Dibujamos el mallado con forma de un cuarto de anillo circular, que representa los puntos interiores del sólido:

Mallado
%Limpeza de programas anteriores 
clear
clc
%Mallado interior de la placa plana
x=1:0.05:2;
y=pi/4:0.05:3/4*pi;
[rr,tt]=meshgrid(x,y);
%Cambio de coordenadas cilíndricas a cartesianas
Mx=rr.*cos(tt);
My=rr.*sin(tt);
mesh(Mx,My,My*0);
%Establecemos los límites
axis([-3,3,-1,3]);
view(2)
%Escribimos el título 
title('MALLADO DE LA PLACA')
%Nombramos los ejes
xlabel('EJE X')
ylabel('EJE Y')


2 .-Dibujo de las curvas de nivel de la temperatura y punto en el que es máxima

Dibujamos las curvas de nivel de la temperatura dada al principio del problema. En función de la temperatura, se puede observar una variación de colores con una gráfica. La zona más fría son los colores oscuros y la zona más caliente son los colores claros. En la gráfica, podemos encontrar el punto donde la temperatura es máxima: [-1.4,1.4] y [1.4,1.4].

Campo de temperaturas y curvas de nivel
clear
clc
% Definimos las regiones
x=1:0.05:2;
y=pi/4:0.05:3*pi/4;
% Creamos el mallado
[rr,tt]=meshgrid(x,y);
Mx=rr.*cos(tt); 
My=rr.*sin(tt);
mesh(Mx,My,My*0);
% Función temperatura
T=Mx.^2+(My-1).^2;
subplot(1,2,1)
% Aplicamos la función de la temperatura al mallado
surf(Mx,My,T);
% Establecemos el límite de los ejes
axis([-3,3,-1,3]);
view(2)
% Escribimos el título y ejes
title('CAMPO DE TEMPERATURAS')
subplot(1,2,2)
% Representación de la temperatura y las curvas de nivel
contour(Mx,My,T);
axis([-3,3,-1,3]);
title('CURVAS DE NIVEL')
% Barra de indicación de colores
colorbar


3 .-Gradiente de la temperatura

Vamos a calcular y a representar el gradiente de la temperatura. Utilizando nuestra función de la temperatura, el gradiente es:

centro

Con la gráfica, podemos observar que el gradiente es ortogonal a las curvas de nivel de la temperatura. Además, sabemos que el gradiente de cualquier función es la variación máxima entre varios puntos.

Gradiente de la temperatura
clc
clear
% Definimos regiones
h = 0.1;
x=[1:h:2];   	
y=[pi/4:h:3*pi/4];
% Creamos el mallado
[rr,tt]= meshgrid(x,y);
Mx=rr.*cos(tt);
My=rr.*sin(tt);
% Función temperatura
T = Mx.^2 + (My-1).^2;
contour(Mx,My,T,30);
% Gradiente de T
dx = 2*Mx;
dy = 2*(My-1);
% 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(Mx,My,dx,dy);
axis equal
colorbar


4 .-Campo de vectores de deslizamiento

Representamos el campo de vectores de deslizamiento en la placa plana. El campo dado en el enunciado está en coordenadas cilíndricas y hacemos un cambio a coordenadas cartesianas, ya que conocemos que el vector [math]\vec{e}[/math]ρ = cosθ [math]\vec{i}[/math] + sinθ [math]\vec{j}[/math].

El campo de vectores de deslizamiento en coordenadas cartesianas es:

centro

Utilizaremos el comando quiver() para representar el campo de vectores.

Campo de vectores
clear
clc
x=1:0.1:2;
y=pi/4:0.1:3/4*pi;
[rr,tt]=meshgrid(x,y);
Mx=rr.*cos(tt); 
My=rr.*sin(tt);
mesh(Mx,My,My*0);
axis([-3,3,-1,3]);
view(2)
title('CAMPO DE VECTORES')
xlabel ('EJE X')
ylabel ('EJE Y')
hold on
mx=(rr.*sin(6.*(tt-pi/4)).*cos(tt))/20;
my=(rr.*sin(6.*(tt-pi/4)).*sin(tt))/20;
% Campo vectorial en 2D
quiver(Mx,My,mx,my);
axis equal
hold off


5 .-Deslizamiento del sólido

A continuación, vamos a representar el sólido antes y después del desplazamiento dado por el campo de vectores [math]\vec{u}[/math] (en t=0). Realizaremos tres gráficas: situación inicial, situación final y comparación. En nuestro campo, podemos observar que sólo varía en la dirección del vector unitario [math]\vec{e}[/math]ρ.

Variación del deslizamiento
clear
clc
x=1:0.1:2;
y=pi/4:0.1:3/4*pi;
[rr,tt]=meshgrid(x,y);
Mx=rr.*cos(tt); 
My=rr.*sin(tt);
% Situación inicial
subplot(2,2,1)
surf(Mx,My,My*0);
axis([-3,3,-1,3]);
view(2)
title('SITUACIÓN INICIAL') 
% Situación final
subplot(2,2,2)
mx=(rr.*sin(6.*(tt-pi/4)).*cos(tt))/20;
my=(rr.*sin(6.*(tt-pi/4)).*sin(tt))/20;
rx=Mx+mx;
ry=My+my;
surf(rx,ry,ry*0);
axis([-3,3,-1,3]);
view(2)
title('SITUACIÓN FINAL')
% comparación
subplot(2,2,3)
plot3(Mx,My,My*0);
hold on
plot3(rx,ry,ry*0);
axis([-3,3,-1,3]);
view(2)
title('COMPARACIÓN')
hold off


6 .-Divergencia del campo vectorial

La divergencia es una medida del cambio de volumen local debido al deslizamiento. Calculamos la divergencia del campo vectorial [math]\vec{u}[/math], y además, calculamos analíticamente los puntos en los que la divergencia de u es máxima, mínima y nula.

centro
Divergencia del campo vectorial
clear
clc
x=1:0.1:2; 
y=pi/4:0.1:3/4*pi; 
[rr,tt]=meshgrid(x,y); 
Mx=rr.*cos(tt); 
My=rr.*sin(tt);
%representamos la divergencia 
div=(sin(6.*(tt-pi/4)))/10; 
surf(Mx,My,div);
shading flat 
axis([-3,3,-1,3]);
axis equal
view(2) 
title('DIVERGENCIA')
xlabel('EJE X'); 
ylabel('EJE Y');
colorbar
%para mostrar un degradado de colores
shading interp


7 .-Rotacional del campo

Calculamos el rotacional de campo [math]\vec{u}[/math] y deducimos que puntos sufren un mayor rotacional:

centro
Rotacional del campo
Rotacional del campo
clear clc
x=1:0.1:2;
y=pi/4:0.1:3/4*pi; 
[rr,tt]=meshgrid(x,y); 
Mx=rr.*cos(tt);
My=rr.*sin(tt);
%módulo del rotacional 
rot=(cos(6.*(tt-pi/4)).*(3/10)); 
%dibujamos el rotacional
surf(Mx,My,rot);
shading flat;
axis([-3,3,-1,3]);
axis equal
view(2)
title('MÓDULO ROTACIONAL') 
xlabel('EJE X');
ylabel('EJE Y'); 
colorbar 
shading interp


8 .-Tensor de tensiones

La placa se ve sometida a numerosas tensiones normales.En un medio elástico lineal, isótropo y homogéneo los desplazamientos permiten escribir el tensor de tensiones σij a través de la fórmula: σ=λ∇ * ([math]\vec{u}[/math])1 + 2μϵ , donde:

  • 1 es el tensor identidad en el conjunto de vectores libres del espacio ℝ3.
  • λ, µ son los conocidos como «coeficientes de Lamé», que dependen de las propiedades elásticas de cada material.

Tomando λ=μ=1 debido a las propiedades de nuestro material, calculamos las tensiones normales en la dirección del eje [math]\vec{e}[/math]ρ,[math]\vec{e}[/math]θ, y [math]\vec{e}[/math]z.

8.1 Tensor de deformaciones y de tensiones

Calculamos el tensor de deformaciones y de tensiones:

centro
centro

8.2 Tensor de tensiones en la dirección [math]\vec{e}[/math]ρ

left
right
%tensiones normales al eje e-ro
clear
clc
x=1:0.1:2;
y=pi/4:0.1:(3*pi)/4;
[rr,tt]=meshgrid(x,y);
Mx=rr.*cos(tt); 
My=rr.*sin(tt);
%elemento (1,1) de la matriz tensión
s1=(sin(6.*(tt-pi/4)).*(1/5));
surf(Mx,My,s1);
shading interp
view(2)
axis([-3,3,-1,3]);
title('TENSIONES NORMALES AL EJE E DE RO')
xlabel('EJE X');
ylabel('EJE Y');
zlabel('EJE Z');
colorbar


8.3 Tensor de tensiones en la dirección [math]\vec{e}[/math]θ

left
right
%tensiones normales al eje e-teta
clear
clc
x=1:0.1:2;
y=pi/4:0.1:(3*pi)/4;
[rr,tt]=meshgrid(x,y);
Mx=rr.*cos(tt); 
My=rr.*sin(tt);
%elemento (1,1) de la matriz tensión
s2=(sin(6.*(tt-pi/4)).*(1/5));
surf(Mx,My,s2);
shading interp
view(2)
axis([-3,3,-1,3]);
title('TENSIONES NORMALES AL EJE E DE TETA')
xlabel('EJE X');
ylabel('EJE Y');
zlabel('EJE Z');+
colorbar


8.4 Tensor de tensiones en la dirección [math]\vec{e}[/math]z

left
right
%tensiones normales al eje e-z
clear
clc
x=1:0.1:2;
y=pi/4:0.1:(3*pi)/4;
[rr,tt]=meshgrid(x,y);
Mx=rr.*cos(tt); 
My=rr.*sin(tt);
%elemento (1,1) de la matriz tensión
s3=(sin(6.*(tt-pi/4)).*(1/10));
surf(Mx,My,s3);
shading interp
view(2)
axis([-3,3,-1,3]);
title('TENSIONES NORMALES AL EJE E DE Z')
xlabel('EJE X');
ylabel('EJE Y');
zlabel('EJE Z');
colorbar


9 .-Tensiones tangenciales

Calculamos las tensiones tangenciales respecto al plano ortogonal [math]\vec{e}[/math]ρ :

centro
En 3D
En 2D
clear
clc
%Definición de regiones
x=1:0.1:2;
y=pi/4:0.1:3*pi/4;
[rr,tt]=meshgrid(x,y);
Mx=rr.*cos(tt);
My=rr.*sin(tt);
%Tensión tangencial
tg=abs(6*cos(6*(tt-pi/4))/20);
surf(Mx,My,tg);
shading interp
axis([-3,3,-1,3]);
title('TENSIONES TANGENCIALES RESPECTO AL PLANO ORTOGONAL')
xlabel('EJE X');
ylabel('EJE Y');
zlabel('EJE Z');
colorbar
Obtenido de «https://mat.caminos.upm.es/w/index.php?title=Campos_escalares_en_una_placa_plana&oldid=54010»