T6 MNI 2021b

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 2

Universidad Nacional de Ingeniería

Facultad de Ingeniería Civil 01-12-2021


Unidad de Posgrado

MÉTODOS NUMÉRICOS EN INGENIERÍA

Sexto Trabajo 2021-II

u   0 1 u 
1. Para resolver las ecuaciones diferenciales      
 
v    9 0 v 
¿Puede usarse el método de Euler? ¿Un método de Runge-Kutta de segundo orden? ¿Un método de
Runge-Kutta de cuarto orden? ¿O el método de Crank-Nicholson? Explique por qué basándose en una
descomposición modal. En cada caso afirmativo, ¿Cuál sería el máximo  x para que el procedimiento sea
estable?
Prepare una hoja de cálculo para resolver las ecuaciones (en su forma original, es decir, sin desacoplar) con
uno de los procedimientos antes mencionados y compare las soluciones con distintos  x .

2. Utilizando el método de Newton, resuelva el problema de valor frontera:

y  
1
8
32  2 x 3  y y   para 1 x  3 con y (1)  17 e y (3)  43 / 3

Con tal fin resuelva el problema de valor inicial:

1
u  
8
 
32  2 x 3  u u  para 1  x  3 con u (1,  )  17 y u (1, )  

Considere x  0.10 . Inicie el proceso con  0  14 (que sería la solución exacta)

u
Por otro lado, considerando v  resuelva el problema:


v   
1
u v  u v  para 1  x  3 con v(1,  )  0 y v (1,  )  1
8
u (3,  k )  43 / 3
Entonces, para corregir  itere con:  k 1   k 
v(3,  k )
Compare la solución obtenida con la solución exacta: y  x 2  16 / x

3. Para una viga en voladizo, de longitud L y de sección constante, empotrada en el extremo x  0 y


sometida a una carga q cualquiera, la energía potencial puede escribirse como:
L L
 (v )   1
0 2
EI (v) 2 dx   0
v q dx  mínimo

a. Determine la ecuación de Euler Lagrange del funcional


b. Considere el caso particular de una carga concentrada en el extremo: q   P  (L) de modo que
L
0
v q dx   P v( L) . Entonces, con la aproximación: v   1 x 3   2 x 2 determine  1 y  2 usando el

método de Rayleigh-Ritz y luego obtenga v(L) y v (L) .


 L  v  L  v 
Note que
 i
  0
EI v   dx 

  i 
 0

 
 i
 q dx  0


Px 2
c. Para el mismo caso de cargas, según la teoría de vigas: v( x)   ( L  x / 3) . Compare con los
2 EI
resultados antes obtenidos y comente los resultados.

d. Repita el paso b, pero considerando una fuerza uniformemente repartida q   w y compare con la
solución dada por la teoría de vigas:
w
v (( L  x) 4  4 L3 ( L  x)  3L4 )
24 EI
Los resultados solo coinciden en dos puntos. ¿A qué se debe que no coincidan en toda la longitud?

4. . El principio variacional que gobierna el pandeo de una columna como la que se muestra en la figura
siguiente es:
2
 d 2u   du 
2

  dx  12 k u ( L)   estacionar io
P
 (u )   EI  2  dx  
1 2
2 L
 dz  2 L  dz 

u
con las restricciones u (0)  0 y  0 . EI, k son constantes.
z z 0

a. Determine la ecuación de Euler-Lagrange del funcional y las


condiciones de borde implícitas al hacer estacionario el
funcional.

b. Con la aproximación: u  a1 z 2  a 2 z 3

 a1   a1 
obtenga ecuaciones de la forma: A    P B  
a 2  a 2 

A y B son matrices simétricas de orden 2.


c. Suponiendo valores numéricos para L, EI , k determine el menor valor de P  0 que satisface las
ecuaciones precedentes y la correspondiente forma de modo (es decir a1 , a 2 ).

El trabajo deberá ser enviado al aula virtual del curso el miércoles 15 de diciembre de 2021 (o antes). Como
excepción, se aceptarán también trabajos entregados hasta el 22 de diciembre, sin que ello afecte la
calificación. En caso se entregue el trabajo después de esa fecha, se reducirán tres puntos en la
calificación. No se aceptará este trabajo después del 29 de diciembre de 2021.

También podría gustarte