Visualización de campos escalares y vectoriales en elasticidad(A18)
 Este artículo está en versión beta. El autor de este artículo no lo ha terminado todavía, por favor no lo edites hasta que elimine este mensaje.
|
1 Introducción
Deseamos visualizar e interpretar campos escalares y vectoriales, para ello nos vamos a servir del estudio de una placa rectangular, basándonos en herramientas como el programa informático Matlab. Dicho programa nos proporcionará las imágenes que se mostrarán como ejemplo.
Consideramos una placa rectangular plana (en dimensión 2) que ocupa la región [math] [-1/2,1/2] \times [0,2][/math]. En ella vamos a suponer que tenemos definidas dos cantidades físicas: la temperatura [math]T(x,y,t)[/math], que depende de las dos variables espaciales [math](x,y)[/math] y el tiempo [math]t[/math], y los desplazamientos [math]\vec u(x,y,t)[/math]. De esta forma, si definimos [math]r_0(x,y)[/math] el vector de posición de los puntos de la placa en reposo, la posición de cada punto [math](x,y)[/math] de la placa en un instante de tiempo [math]t[/math] viene dada por: [math] \vec r (x,y,t)= \vec r_{0}(x,y)+\vec u(x,y,t). [/math] Vamos a suponer que sobre la placa se ha aplicado una fuerza que ha provocado una vibración de manera que los desplazamientos vienen dados por la onda: [math] \vec u(x,y,t) = \vec a \sin(\vec b \cdot \vec r_0-ct), [/math] donde [math]\vec a[/math] se conoce como amplitud, [math]\vec b[/math] es la fase que indica la dirección de propagación y [math]c/|\vec b|[/math] es la velocidad de propagación.
Si [math]\vec a [/math] es paralelo a [math]\vec b[/math] diremos que la onda es longitudinal mientras que si es perpendicular hablaremos de onda transversal.
En este trabajo vamos a centrarnos en las ondas transversales. Supondremos lo siguiente: [math] \vec a=\frac{\vec i}{10}, \qquad \vec b= \pi \vec j, \qquad t=0. [/math] En este caso, [math]\vec u(x,y)=\frac{\sin(\pi y)}{10}\vec i[/math].
2 Realización del Mallado
Comenzaremos dibujando el mallado de los puntos interiores del sólido. Mostramos el código matlab a continuación:
h=0.1;
x=-0.5:h:0.5; % Vector x con valores entre -0,5 y 0,5 con saltos de h
y=0:h:2; % Vector y con valores entre 0 y 2 con saltos de h
[xx,yy]=meshgrid(x,y); % Matrices de coordenadas
figure(1) % Preparamos la plantilla
mesh(xx,yy,0*xx) % Dibujo de la malla
axis([-2,2,-1,3]) %Se fijan y centran los ejes
view(2) % Nos permite verlo en 2 dimensiones
3 Obtención del gradiente de temperaturas
Continuamos obteniendo el campo escalar de temperaturas [math]T(\rho,\theta)=-\log(\rho+0.1)[/math]
h=0.1;
x=-0.5:h:0.5; % Vector x con valores entre -0,5 y 0,5 con saltos de h
y=0:h:2; % Vector y con valores entre 0 y 2 con saltos de h
[xx,yy]=meshgrid(x,y); % Matrices de coordenadas
figure(1) % Preparamos la plantilla
f=-log(0.1+sqrt(xx.^2+yy.^2)); % Campo escalar
surf(xx,yy,f) % Dibujo de la malla
axis([-2,2,-1,3]) % Región de dibujo
view(2) % Visualizado en 2 dimensiones
Se observa la variación de temperaturas en función de [math]\rho[/math] y [math]\theta[/math] como expresa el campo [math]T(\rho,\theta)=-\log(\rho+0.1)[/math] .
