TAREA 2

METODOS NUMERICOS Y DISEÑO DE ALGORITMOS I

MEC-130
Nombre: Erbin Javier

Apellido Paterno: Condori Apellido Materno: Espinoza

R.U.: 1780528

a b c d e f g
1 7 8 0 5 2 8
PREGUNTA 1:
El movimiento tridimensional de una partícula se define mediante el vector posición:

Donde r y t se expresan en pies y segundos. Determinar el valor diferente de cero más

pequeño de t para el cual el vector de posición y el vector de velocidad son perpendiculares
entre sí, justifique el resultado de manera clara.

Para determinar el valor de los parámetros A y B realizar el siguiente calculo:

A = (a + b + c + 3); B = (d + e + f + 3)

Utilice el método de la Bisección con seis decimales y una tolerancia 10-5.

SOLUCION: A=19 B=10

⃑ ^ ( 10. √t 2 +1 ) ^j+ ( t . sen ( t ) ) k^
r =( 19.t .cos ( t ) ) i+

⃑ ^
v=19 ( cos ( t )−t . sen ( t ) ) i+10
(√ ) t
2
t +1
^j+ ( sen ( t ) +t .cos ( t ) ) k^

Para que sean perpendiculares: ⃑ ⃑

r ∙ v=0

0=19 ( cos ( t )−t . sen ( t ) ) ( 19. t . cos ( t )) + 10

( √ )(
2
t
t +1
10. √ t 2+ 1 ) + ( sen ( t ) +t . cos ( t ) ) ( t . sen (t ) )

0=361 ( t . cos ( t ) −t . sen ( t ) cos ( t ) ) +100. t+ ( t . sen ( t ) +t . sen ( t ) cos ( t ))

2 2 2 2

0=361t . cos2 (t )−361 t 2 . sen (t ) cos ( t )+100. t+ t . sen 2 ( t )+ t 2 . sen ( t ) cos ( t )

0=361 cos2 (t )−361 t . sen ( t ) cos ( t )+ 100+sen 2 ( t ) +t . sen (t ) cos ( t )

2 2
0=361 cos (t )−360 t . sen ( t ) cos ( t ) +sen ( t ) +100
 Matriz de Evaluacion:

X F(X)
0 461
0,5 302,522026
1 42,4200326
1,5 64,6989484
2 435,793047
2,5 763,575117
3 604,715021
3,5 2,80084859
4 -457,527944
4,5 -216,81942
1ra ITERACIÓN:

a=3,5 b=4 tol=0,00001

a+ b 3,5+4
x r= = =3,75
2 2

|f ( xr )|=|f ( 3,75 )|=289,755627<0,00001 FALSO

f ( a ) . f (x r )=f ( 3,5 ) . f (3,75)=−811,561639<0 VERDADERO

f ( a ) . f ( xr )=f ( 3,5 ) . f ( 3,75 ) =−811,561639>0 FALSO

∴ a=3,5 b=3,75

METODO BISECCION

i Xa Xb Xr f(Xa) f(Xr) f(Xa).f(Xr) %Ea ERROR

1 3,5 4 3,75 2,80084859 -289,755627 -811,561639 - -
2 3,5 3,75 3,625 2,80084859 -153,833922 -430,865524 3,44827586 FALSO
3 3,5 3,625 3,5625 2,80084859 -77,3784444 -216,725307 1,75438596 FALSO
4 3,5 3,5625 3,53125 2,80084859 -37,6610431 -105,482879 0,88495575 FALSO
3,51562
5 3,5 3,53125 2,80084859 -17,511429 -49,0468612 0,44444444 FALSO
5
3,50781
6 3,5 3,515625 2,80084859 -7,37415441 -20,65389 0,22271715 FALSO
3
3,50390
7 3,5 3,5078125 2,80084859 -2,29118528 -6,41726306 0,11148272 FALSO
6
3,50195
8 3,5 3,50390625 2,80084859 0,25372153 0,71063558 0,05577245 FALSO
3
3,50293
9 3,50195313 3,50390625 0,25372153 -1,01901228 -0,25854535 0,02787845 FALSO
0
1 3,50244
3,50195313 3,50292969 0,25372153 -0,38271512 -0,09710306 0,01394117 FALSO
0 1
1 3,50219
3,50195313 3,50244141 0,25372153 -0,06451419 -0,01636864 0,00697107 FALSO
1 7
1 3,50207
3,50195313 3,50219727 0,25372153 0,09459933 0,02400189 0,00348566 FALSO
2 5
1 3,50213
3,5020752 3,50219727 0,09459933 0,01504148 0,00142291 0,0017428 FALSO
3 6
1 3,50216
3,50213623 3,50219727 0,01504148 -0,02473662 -0,00037208 0,00087139 FALSO
4 7
1 3,50215
3,50213623 3,50216675 0,01504148 -0,00484764 -7,2916E-05 0,0004357 FALSO
5 1
1 3,50214
3,50213623 3,50215149 0,01504148 0,00509691 7,6665E-05 0,00021785 FALSO
6 4
1 3,50214
3,50214386 3,50215149 0,00509691 0,00012463 6,3523E-07 0,00010892 FALSO
7 8
1 3,50215
3,50214767 3,50215149 0,00012463 -0,00236151 -2,9431E-07 5,4462E-05 FALSO
8 0
1 3,50214
3,50214767 3,50214958 0,00012463 -0,00111844 -1,3939E-07 2,7231E-05 FALSO
9 9
2 3,50214
3,50214767 3,50214863 0,00012463 -0,0004969 -6,1929E-08 1,3616E-05 FALSO
0 8
2 3,50214
3,50214767 3,50214815 0,00012463 -0,00018614 -2,3198E-08 6,8078E-06 FALSO
1 8
2 3,50214
3,50214767 3,50214791 0,00012463 -3,0754E-05 -3,8328E-09 3,4039E-06 FALSO
2 8
2 3,50214
3,50214767 3,50214779 0,00012463 4,6938E-05 5,8499E-09 1,7019E-06 FALSO
3 8
2 3,50214
3,50214773 3,50214779 4,6938E-05 8,0923E-06 3,7984E-10 8,5097E-07 VERDADERO
4 8

x=3,502148

t =3,502148 [ s ]

Pregunta 2
Utilizando dos elementos finitos lineales y calcular los desplazamiento, reacciones de la Figura.
E = 210 (GPa); I = 4x10-4 (m4); k = 10 (kN/m)
Para determinar el valor de los parámetros A y B que representan las cargas realizar el
siguiente calculo:
A = (a + b + c + 5) [kN]; B = (a + b + c + 7) [kN/m]
SOLUCION:

A=21 B=23

Discretizando:

2 4 6
θ1 θ2 θ3

1 E1 2 E2 3

1 3 5
V1 V2 V3

Trabajando todo en KPa, m, KN, tenemos:

E = 210x106 (KPa); I = 4x10-4 (m4); k = 10 (KN/m)

L1 = 6m L2 = 6m

[ K ] = E .3I 1
L [−1 −11 ]
Para el elemento E1:

[ ]
12 6 ( 6 ) −12 6 ( 6 )
6 −4 2 2
210 ×10 ∗4 ×10 6 ( 6) 4 ( 6) −6 ( 6 ) 2 ( 6 )
[ K 1 ]= 3
6 −12 −6 ( 6 ) 12 −6 ( 6 )
2 2
6 ( 6 ) 2 ( 6 ) −6 ( 6 ) 4 ( 6 )

[ ]
4667 14000 −4667 14000 1
14000 56000 −14000 28000 2
[ K 1 ]=
−4667 −14000 4667 −14000 3
14000 28000 −14000 56000 4

Para el elemento E1:

 [ ]
12 6 ( 6 ) −12 6 ( 6 )
6 −4 2
6 ( 6) 4 ( 6 ) −6 ( 6 ) 2 ( 6 )2
[ K 2 ] = 210 ×10 ∗4
3
×10
6 −12 −6 ( 6 ) 12 −6 ( 6 )
2 2
6 ( 6 ) 2 ( 6 ) −6 ( 6 ) 4 ( 6 )

[ ]
4667 14000 −4667 14000 3
[ K 2 ] = 14000 56000 −14000 28000 4
−4667 −14000 4667 −14000 5
14000 28000 −14000 56000 6

Calculamos los vectores elementales para las cargas:

 De la carga puntual:
a 2 1
s= = =
L 6 3

{ } { }{ }
3 2
2 s −3 s + 1 0,7407 −15,5 1
L ( s −2 s +s )
3 2
0,8889 −18,7 2
f =P . =−21. =
2
3 s −2 s 3
0,2593 −5,4 3
L ( s −s )
3 2 −0,4444 9,3 4

{ }{ }
 De la carga distribuida:
−3 wL −3∗23∗6
20 20

{ }
2
−w L −23∗62 −20,7 3
f= 30 = 30 = −27,6 4
−7 wL −7∗23∗6 −48,3 5
20 20 41,4 6
2 2
wL 23∗6
20 20

La matriz global:

[ ]
4667 14000 −4667 14000 0 0
14000 56000 −14000 28000 0 0
−4667 −14000 9334 0 −4667 14000
14000 28000 0 112000 −14000 28000
0 0 −4667 −14000 4667 −14000
0 0 14000 28000 −14000 56000

Sabemos que:

1 3
Entonces:

[ ]{ } { } { }
4667 14000 −4667 14000 0 0 0 −15,5 F 1
14000 56000 −14000 28000 0 0 θ 1 −18,7 0
−4667 −14000 9334 0 −4667 14000 . v 2 = −26,1 + 0
14000 28000 0 112000 −14000 28000 θ2 −18,3 0
0 0 −4667 −14000 4667 −14000 0 −48,3 F 3
0 0 14000 28000 −14000 56000 θ3 41,4 0

[ ]{ } { }
56000 −14000 28000 0 θ1 −18,7
−14000 9334 0 14000 v2 = −26,1
28000 0 112000 28000 θ2 −18.3
0 14000 28000 56000 θ 41,4
3

θ1=−4,562 ×10−3 v 2 =−0,0176 θ2=−3,530 × 10−4 θ3 =5,321×10−3

{ }{
−3
−4,562 ×10

[ 14000 −4667
0 −4667
14000 0 −0,0176
−14000 −14000 −3,530 ×10]−4 =

−3
−15,5+ F 1
−48,3+ F 3 }
5,321 ×10

F 1=28,30 F 3=60,88

v1 0
v2 −0,0176 [ m ]
v3 0
θ1 −4,562×10−3 [ rad ]
θ2 −3,530 ×10−4 [ rad ]
θ3 5,321× 10−3 [ rad ]
F1 28,30 [ kN ]
F2 0
F3 60,88 [ kN ]
M1 0
M2 0
M3 0
F resorte -176×10−3 [ kN ]

Pregunta 3
Una barra de sección transversal circular de radio 4cm y 80cm de longitud transporta calor en
dirección de su eje. El extremo izquierdo se mantiene a 100 ºC, el extremo derecho se
mantiene a 20 ºC. La conductividad térmica del material es 10W/ (m. ºC). La barra no está
aislada lateralmente, B=40W/ (ºC.m2) la temperatura del ambiente es 10ºC, la generación de
calor interna de calor es f = A * (10x + 20) (W/m).
Utilizando cuatro elementos finitos lineales calcular la distribución de temperatura.
Para determinar el valor del parámetro A realizar el siguiente cálculo con los dígitos de su R: U:
A = (a + b + c + d + e + f + g + 20)
SOLUCION:
A = 51

k e =10
[ W
m .ºC ]
Be =40
W
[
ºC . m
2 ]
he =0,2 [ m ] r e =0,04 [ m ] T ∞ =10ºC

Pe =
2π
25
[m ] A e=
π
625 [ ]
[ m2 ] q ( x )=51 ( 10 x+ 20 ) W
m

Discretizando:

1 e1 2 e2 3 e3 4 e4 5

0,2 m 0,2 m 0,2 m 0,2 m

Modelo del elemento:

Matriz elemental para e1:

( ]) ( [ ]) [ ]
π 2π 22 π 2π
∗10 ∗40∗0,2
k 1=
625
0,2
1 −1 + 25
−1 1 [
6
2 1 = 75
1 2 2π
75 1
22 π 2
75 75
Matriz elemental para e2:

( ]) ( ]) [ ]
π 2π 22 π 2π
∗10 ∗40∗0,2
k 2=
625
0,2
1 −1 + 25
−1 1 [6
2 1 = 75
1 2 2π [ 75 2
22 π 3
75 75
Matriz elemental para e3:

( ]) ( ]) [ ]
π 2π 22 π 2π
∗10 ∗40∗0,2
k 3=
625
0,2
1 −1
−1 1
+
25
[6
2 1
1 2
=
75
2π [ 75 3
22 π 4
75 75
Matriz elemental para e4:
 ( ]) ( [ ]) [ ]
π 2π 22 π 2π
∗10 ∗40∗0,2
k 4=
625
0,2
1
[−1
−1 1
+
25
6
2 1
1 2
= 75
2π
75 4
22 π 5
75 75

Vector elemental para e1:

[ ]
2π 16 π
∗40∗10∗0,2
f 1=
25
2
1
1
=
[]5 1
16 π 2
5
Vector elemental para e2:

[][ ]
2π 16 π
∗40∗10∗0,2
25 1 5 2
f 2= =
2 1 16 π 3
5
Vector elemental para e3:

[][ ]
2π 16 π
∗40∗10∗0,2
25 1 5 3
f 3= =
2 1 16 π 4
5
Vector elemental para e4:

[][ ]
2π 16 π
∗40∗10∗0,2
25 1= 5 4
f 4=
2 1 16 π 5
5

Distribución de calor interno variable para e1: x 1=0 x 2=0,2

r 1=¿

Distribución de calor interno variable para e2: x 1=0,2 x2 =0,4

r 1=¿

Distribución de calor interno variable para e3: x 1=0,4 x 2=0,6

r 1=¿

Distribución de calor interno variable para e4: x 1=0,6 x 2=0,8

r 1=¿
 [ ]
MATRIZ GENERAL:

22 π 2π
0
75 75 0 0
2π 44 π 2π 0 0
75 75 75 2π
0
2π 44 π 75
0
75 75
44 π 2π
2π
0 0 75 75
75
2π 22 π
0 0 0
75 75
Uniendo todo tenemos:

[ ] []
22 π 2π 16 π
0
75 75 0 0 5

{ }[ ] { }
2π 44 π 2π 0 0 100 32 π
105,4 Q1
75 75 75 2π T2 5
0 224,4 0
2π 44 π 75 32 π
0 ∗ T 3 = 244,8 + + 0
75 75 5
T4 256,2 0
44 π 2π 32 π
2π 20 139,4 Q5
0 0 75 75 5
75
2π 22 π 16 π
0 0 0
75 75 5

[ ][ ]
44 π 2π 32 π 8 π
0 224,4+
−

{}
75 75 5 3
T2
2π 44 π 2π 32 π
∗ T 3 = 244,8+
75 75 75 5
2π 44 π T 4 32 π 8 π
0 256,2+ −
75 75 5 3

T 2=122,12T 3=131,84 T 4 =139,38

T1 100 [ ºC ] Q1 −13,07 [ J ]
T2 122,12 [ ºC ] Q2 0
T3 131,84 [ ºC ] Q3 0
T4 139,38 [ ºC ] Q4 0
T5 20 [ ºC ] Q5 -119,34[ J ]