Revisión del 04:01 3 dic 2014

¡Bienvenido a MateWiki! Esperamos que contribuyas mucho y bien. Probablemente quieras leer las páginas de ayuda. Nuevamente, bienvenido y ¡diviértete! Carlos Castro (discusión) 15:40 27 nov 2014 (CET)

1 ENUNCIADO

Consideramos una placa plana (en dimension 2) que ocupa la region comprendida entre las parábolas P1 : 18y-􀀀81x^2-􀀀1 = 0 y P2 : 2y+x^2-􀀀1 = 0. Para representarla usaremos un sistema de coordenadas adaptado a la geometría que nos dan:

• x=uv
• y=(1/2)(u^2-v^2)

con u y v definidas en (u,v) 2 [1/3,1]x[􀀀-1,1].En ella vamos a suponer que tenemos definidas dos cantidades fisicas: la temperatura T(u,v), que depende de las dos coordenadas curvilineas (u,v), y los desplazamientos u(x,y) producidos por la accion de una fuerza determinada. De esta forma, si definimos r0(u,v) el vector de posicion de los puntos de la placa antes de la deformacion, la posicion de cada punto (u,v) de la placa despues de la deformacion viene dada por:

• r(u,v) = r0(u,v) + u(u,v)

Vamos a suponer que la fuerza aplicada sobre la placa ha provocado un desplazamiento de los puntos de la misma dado por el vector de desplazamientos u(u,v) = a (b . r0) donde a y b son vectores dados. En este trabajo supondremos lo siguiente: a=g_u/|g_u| y b=(-4).(g_u/|g_u|)

1.1 Representación gráfica de la placa

En primer lugar, se obtiene el gráfico que representa la placa que consideraremos a lo largo de todo este trabajo, que corresponde a la región comprendida entre dos parábolas dadas(P1 : 18y-􀀀81x^2-􀀀1 = 0 y P2 : 2y+x2^􀀀-1 = 0).Este programa realiza el mallado correspondiente a (u,v) en el intervalo [-1,1]x[-1,1], y con un paso de 1/20 (=h).Para ello utilizamos el siguiente programa:

h=(1/20);
u=(1/3):h:1;
v=-1:h:1;
[uu,vv]=meshgrid(u,v);
xx=uu.*vv;
yy=(1/2).*((uu.^2)-(vv.^2));
figure(1)
subplot(1,2,1);
mesh(xx,yy,0*xx);
axis([-1,1,-1,1]);
view(2);
subplot(1,2,2);
mesh(xx,yy,0*xx);
axis([-1,1,-1,1]);

2 Lineas coordenadas y vectores de la base natural

Definimos las líneas coordenadas como superficies que se obtienen al variar una componente y fijar las restantes. En este apartado dibujaremos dichas líneas y los vectores de la base natural, que cambiarán de dirección en cada punto de la placa, y esto sucede ya que la base natural en estas coordenadas no es constante.

(Con el siguiente programa obtenemos dos gráficos, en el primero de ellos se aprecia g_u aplicado en cada punto de la malla, y en el segundo aparece aplicado g_v )

h=(1/20);
u=(1/3):h:1;
v=-1:h:1;
[uu,vv]=meshgrid(u,v);
xx=uu.*vv;
yy=(1/2).*((uu.^2)-(vv.^2));
guu=vv; guv=uu; gvu=uu; gvv=-vv;
figure(2)
subplot(1,2,1);
quiver(xx,yy,guu,guv);
axis([-1.6,1.6,-0.8,0.8]);
subplot(1,2,2);
quiver(xx,yy,gvu,gvv);
axis([-1.6,1.6,-0.8,0.8]);

Aplicando ambas componentes de la base natural sobre cada punto de la malla, todo en un mismo gráfico, se obtiene:

figure (3)
hold on
quiver(xx,yy,guu,guv);
quiver(xx,yy,gvu,gvv);
axis([-1.6,1.6,-0.8,0.8]);
hold off

3 Influencia de un foco de calor

4 Influencia de un foco de calor