INTRODUCCIÓN AL MÉTODO DE LOS ELEMENTOS

FINITOS

TRABAJO PRÁCTICO 2
FUNCIONES DE FORMA

Estudiante

FREDY ANDRÉS MERCADO NAVARRO

Pasaporte: 98´773.532
Maestría en Simulación Numérica y Control
Cuatrimestre: II-2011
21 de Noviembre de 2011

Universidad de Buenos Aires

Ciudad Autónoma de Buenos Aires
Argentina
IMEF Facultad de Ingeniería - UBA

PROBLEMA 1

qué sucede con el campo de las deformaciones cuando
 =
 + /4.

En este problema se demostrará que existe una zona de exclusión para el nodo del medio y

Parte a)

Primero planteamos las coordenadas dentro del elemento como:



 =  ℎ

 = ℎ

 + ℎ

 + ℎ




Determinamos la expresión para el operador Jacobiano como:




 ℎ
 ℎ
 ℎ
 ℎ

=
 =
 +
 +









ℎ = − 1 −
, ℎ = 1 +
, ℎ = 1 −
 
2 2
ℎ
 1 ℎ
 1 ℎ

=
− , =
+ , = −2


2 
2 



 1 1
= 
− 
 + 
+ 
 + −2



2 2

Para efectos de nuestro problema específico:

 = 0,
 = ,
 =


La expresión 

/
representa al operador Jacobiano de la transformación de coordenadas
naturales a cartesianas. Se desea mostrar valores de r para los cuales la ubicación de
 no es
válida como herramienta para el análisis por elementos finitos. Iniciamos igualando la
pendiente a cero:



 1
= 0 = 
+   + −2



2
 1

 = 1 + 
2 2

Dado que las integrales serán desarrolladas en el intervalo −1 ≤

≤ 1, evaluamos
 en los
dos extremos:

 3


= −1 = ,

= 1 =
4 4

Page 2 of 20
IMEF Facultad de Ingeniería - UBA

Podemos concluir, que estos dos valores son los límites de la ubicación de
 . Cualquier valor
por debajo de L/4 implicaría que para un X dado existen dos valores de r, y de igual manera
para cualquier valor sobre 3L/4. Esto se comprende mucho más claramente observando las
figuras siguientes:

r vs x x3=L/4 r vs x x3=3L/4
1 1

0.8 0.8

0.6 0.6
x

x
0.4 0.4

0.2 0.2

0 0
-2 -1 0 1 -1 0 1 2
r r

Parte b)

Para hallar la respuesta se planteará y desarrollará la matriz de Deformaciones-

Desplazamientos. Expresando la coordenada en términos de las funciones de forma tenemos:


=  ℎ
 = ℎ
 + ℎ
 + ℎ





 =
 ,
 =
 + ,
 =
 +
4

Para
 =
 +  y tomando
 = 0, tenemos:


•



= − 1 −
 0 +  1 +
  + 1 −
 !  
2 2 4

 

= 1 +
 + 1 −
  ,
= 1 +

2 4 4

Evaluando la matriz de Deformaciones-Desplazamientos, tenemos:

ℎ ℎ ℎ 1 1

" = #$% . ' ( = #$% . '− +
  +
 −2
(



2 2


   + 
1 2
#=' ( = ' +
( = ' (, #$% = ) *=' (

2 2 2  + 
 + 

Para completar el análisis, expresamos

en términos de
:




= 1 +
 ,
= 2+ − 1
4 

Desarrollando la matriz de Deformaciones-Desplazamientos tenemos:

Page 3 of 20
IMEF Facultad de Ingeniería - UBA

2 1 1
"=' ( '− +
  +
 −2
(
 + 
2 2
2 3 1 2 1 1 2 1 4
" = ' −   −   − (
 2√ √
 2√ √
√ √


De este modo se concluye que en
= 0, el elemento posee una singularidad en el campo de

las deformaciones del orden 1/√
.

PROBLEMA 2

En este problema se determinarán las funciones de forma del siguiente elemento utilizando
polinomios de Lagrange.

Para hallar las funciones de forma definimos primero a las funciones de forma como:

ℎ
 = -,. /.
 + -, /
 + -, /
 + -,/
, 0 = 1 … 4 34546.

−
. 
−
  … 8
−
7$98
−
7:9 …
−
; 
/7
 = , < = 0 … 3.
8
7 −
. 98
7 −
 9 … 8
7 −
7$ 98
7 −
7:9 … 8
7 −
; 9

/7
  = =7 , >
43?@A?
5?/BC, =7 = 1 60 0 = <, =7 = 0 60 0 ≠ <

Dado que poseemos cuatro nodos, nuestro polinomio será de grado 3:

−
 
−
 
−
 
−
. 
−
 
−
 
/.
 = , /
 =

. −
 
. −
 
. −
 
 −
. 
 −
 
 −
 

−
. 
−
 
−
 
−
. 
−
 
−
 
/
 = , /
 =

 −
. 
 −
 
 −
 
 −
. 
 −
 
 −
 

Se cumplirá para todos los puntos:

1 1

. = −1,
 = − ,
 = ,
 = 1
3 3

• Para hallar ℎ

ℎ
 = -. /.


Nótese que la expresión anterior resulta de cumplir con la condición ℎ = 1 en

= −1, y cero
en los demás puntos.

Page 4 of 20
IMEF Facultad de Ingeniería - UBA

-. = -
= −1 = 1, - = -
= −1/3 = 0

- = -
= 1/3 = 0, - = -
= 1 = 0

1 1
E
− F− 3GH F
− 3G
− 1

−
 
−
 
−
 
/.
 = =

. −
 
. −
 
. −
  1 1
E−1 − F− 3GH F−1 − 3G −1 − 1

1 1
F
+ 3G F
− 3G
− 1 1
/.
 = =− 9
 − 9
 −
+ 1
2 4 16
F− 3G F− 3G −2

De este modo, se concluye que:

1
ℎ
 = − 9
 − 9
 −
+ 1
16

• Para hallar ℎ

ℎ
 = - /


- = -
= −1/3 = 1, -.,, = 0

1

−
. 
−
 
−
  8
− −19 F
− 3G
− 1
/
 = =

 −
. 
 −
 
 −
  1 1 1 1
E− 3 − −1H F− 3 − 3G F− 3 − 1G

1

+ 1 F
− G
− 1 1
/.
 = 3 = 27
 − 9
 − 27
+ 9
2 2 4 16
F3G F− 3G F− 3G

De este modo, se concluye que:

1
ℎ
 = 27
 − 9
 − 27
+ 9
16

• Para hallar ℎ

ℎ
 = - /


- = -
= 1/3 = 1, -.,, = 0

1
8
− −19 E
− F− 3GH
− 1

−
. 
−
 
−
 
/
 = =

 −
. 
 −
 
 −
  1 1 1 1
E3 − −1H E3 − F− 3GH F3 − 1G

Page 5 of 20
IMEF Facultad de Ingeniería - UBA

1

+ 1 F
+ G
− 1 1
/.
 = 3 = − 27
 + 9
 − 27
− 9
4 2 2 16
F3G F3G F− 3G

De este modo, se concluye que:

1
ℎ
 = − 27
 + 9
 − 27
− 9
16

• Para hallar ℎ

ℎ
 = - /


- = -
= 1 = 1, -.,, = 0

1 1
8
− −19 E
− F− 3GH F
− 3G

−
. 
−
 
−
 
/
 = =

 −
. 
 −
 
 −
  1 1
81 − −19 E1 − F− 3GH F1 − 3G

1 1

+ 1 F
+ G F
− G
/.
 = 3 3 = 1 9
 + 9
 −
− 1
4 2 16
2 F G F G
3 3

De este modo, se concluye que:

1
ℎ
 = 9
 + 9
 −
− 1
16

PROBLEMA 3

A continuación se realizará el procedimiento para colapsar dos nodos en un elemento

cuadrangular, para formar un nuevo elemento de forma triangular. Las funciones de forma
para el elemento de cuatro nodos son:

1
ℎ = 1 +
 1 + 6
4

Page 6 of 20
IMEF Facultad de Ingeniería - UBA

1
ℎ = 1 −
 1 + 6
4
1
ℎ = 1 −
 1 − 6
4
1
ℎ = 1 +
 1 − 6
4

A continuación se suman las funciones de forma de los nodos 3 y 4 para dar origen a una
nueva función de forma correspondiente al nodo resultante:

1 1
ℎ + ℎ = ℎ∆ = 1 −
 1 − 6 + 1 +
 1 − 6
4 4
1
ℎ∆ = 1 − 6 1 −
+ 1 +

4
1
ℎ∆ = 1 − 6
2

De este modo, las tres funciones de forma resultantes para el nuevo elemento triangular son:

1
ℎ = 1 +
 1 + 6
4
1
ℎ = 1 −
 1 + 6
4
1
ℎ∆ = 1 − 6
2

PROBLEMA 4

las coordenadas cartesianas para los siguientes elementos en los puntos

= 0, 6 = 0 y
Se evaluará la matriz Jacobiana, el determinante del Jacobiano y el gradiente con respecto de

= ,6 = .
 
√ √

La matriz Jacobiana de define a partir de la siguiente expresión:

Page 7 of 20
IMEF Facultad de Ingeniería - UBA

 
O 
R U
M N = M 



N . Q
T
 
O Q  T
6 6 6 POS


O
 
=V , V = M 

N

, 6 
, O 
O
6 6

Para los casos a) y b)

Ambos elementos poseen el mismo número de nodos dispuestos sobre los mismos lados del
elemento, y en consecuencia, el mismo número y formulación de funciones de forma.

 O

O ℎ ℎ ℎ ℎ ℎW ℎX R
 O U
Q O T
V = M 

N = M 





N . Q
 T

O ℎ ℎ ℎ ℎ ℎW ℎX Q
 O T
6 6 6 6 6 6 6 6 Q
W OW T
P
X OX S

Para Y% Para YZ

1 1 1 1 1
ℎ = 1 +
 1 + 6 − ℎW − ℎX ℎ = 1 −
 1 + 6 − ℎW
4 2 2 4 2

ℎ 1 1 1 ℎ 1 1
= ' 1 + 6  +
 − 1 − 6  ( = 
−  1 + 6

2 2 2 
2 2
ℎ 1 1 1 ℎ 1
= ' 1 +
  + 6 − 1 −
 ( =

− 1
6 2 2 2 6 4

Para Y[ Para Y\

1 1 1
ℎ = 1 −
 1 − 6 ℎ = 1 +
 1 − 6 − ℎX
4 4 2

ℎ 1 ℎ 1
= − 1 − 6 = 6 6 − 1

4 
4
ℎ 1 ℎ 1
= − 1 −
 = 6 6 − 1 1 +

6 4 6 4

Para Y] Para Y^

1 1
ℎW = 1 −
  1 + 6 ℎX = 1 − 6   1 +

2 2

ℎW ℎX 1
= −
1 + 6 = 1 − 6  


2
ℎW 1 ℎX
= 1 −
  = −6 1 +

6 2 6

Page 8 of 20
IMEF Facultad de Ingeniería - UBA

Evaluando las expresiones anteriores con Matlab ®, los resultados son los siguientes:

Para el caso a)

La matriz de coordenadas cartesianas de los nodos es:

 O 2 2
R
 O U R 0 2U
Q
O T Q 0 T
Q  T Q 0T
Q
 O T = Q 2 0T
Q
W OW T Q1.5 2T
P
X OX S P 2 1S

• En
= 0, 6 = 0


O 
6
R U
V = M 

N = ' 1 0
(, det V = 1, V$


= Q 


T 1 0

O 3/4 1 Q 6 T = '−3/4 1(
6 6 PO OS

En F
= ,6 = G
 
√ √
•


O 311
= 0.5447 −2.2204? − 16 311
V = M 

N = M 571 N, det V = = 0.5447

O 338 571
− = −0.1478 1
6 6 2287


6 571
R U = 1.8360 4.0768? − 16


T 311
V$ = Q 
6 T = M N
Q 769
= 0.2713 1
PO OS 2834

Page 9 of 20
IMEF Facultad de Ingeniería - UBA

Para el caso b)

La matriz de coordenadas cartesianas de los nodos es:

 O 3 2
R
 O U R 1 2U
Q
O T Q 0 T
Q  T Q 0T
Q
 O T = Q 2 0T
Q
W OW T Q2.5 2T
P
X OX S P2.5 1S

• En
= 0, 6 = 0


O 
6
R U
V = M 

N = ' 1 0
(, det V = 1, V$


= Q 


T 1 0

O 5/4 1 Q 6 T = '−5/4 1(
6 6 PO OS

En F
= ,6 = G
 
√ √
•


O 311
= 0.5447 −2.2204? − 16
V = M 

N = M571 N

O 311
= 0.3522 1
6 6 883
311
det V = = 0.5447
571

6 571
R U = 1.8360 4.0768? − 16


T 311
V$ = Q 
6 T = M N
Q 571
− = −0.6467 1
PO OS 883

Page 10 of 20
IMEF Facultad de Ingeniería - UBA

Para el Caso c)


O ℎ ℎ ℎ ℎ
 O

V = M 

N = M 



N . M
 O

O ℎ ℎ ℎ ℎ
 O N
6 6 6 6 6 6
 O

 O 50 60

 O 15 50
M
O N = M20 N
 20

 O 40 30

Para Y% Para YZ

1 1
ℎ = 1 +
 1 + 6 ℎ = 1 −
 1 + 6
4 4

ℎ 1 ℎ 1
= 1 + 6 = − 1 + 6

4 
4
ℎ 1 ℎ 1
= 1 +
 = 1 −

6 4 6 4

Para Y[ Para Y\

1 1
ℎ = 1 −
 1 − 6 ℎ = 1 +
 1 − 6
4 4

ℎ 1 ℎ 1
= − 1 − 6 = 1 − 6

4 
4
ℎ 1 ℎ 1
= − 1 −
 = − 1 +

6 4 6 4

• En
= 0, 6 = 0


O 55
= 13.75 5
V = M

N = M 4 N, det V = 200

O 5
= 1.25 15
6 6 4

Page 11 of 20
IMEF Facultad de Ingeniería - UBA


6 3 1
R U = 0.0750 − = −0.0250


T 40 40
V$ = Q 
6 T = M N
Q 1 11
− = −0.0063 = 0.0688
PO OS 160 160

En F
= ,6 = G
 
√ √
•


O 4122
= 15.9151 5 18397
V = M

N = M 259 N, det V = = 221.6506

O 2131 83
= 3.4151 15
6 6 624

6 103 103
R U = 0.0677 − = 0.0226


T 1522 4566
V$ = Q 
6 T = M N
Q 52 96
− = −0.0154 = 0.0718
PO OS 3375 1337

PROBLEMA 5

Se evaluará el vector de fuerza volumétrico de cuando el elemento de la figura siguiente está

sometido a una fuerza por unidad de volumen de fg h = 10?i + 5?j .

En forma general, para un conjunto de elementos finitos, dicho vector está definido como:

de =  k l hm fg h 5n h
h o p

Dado que solo se considera un elemento, solo se tendrá en cuenta una integral sobre el
volumen:

de = k l hm fg h 5n h
o p

Se presentan las siguientes definiciones:

ℎ 0 ℎ 0 ℎ 0 rig 10
lm = ' (, fg h = q g s =   , 5n h = 5?B V B 5
56
0 ℎ 0 ℎ 0 ℎ rj 5

Page 12 of 20
IMEF Facultad de Ingeniería - UBA

ℎ 0
R0 ℎ U 
O ℎ ℎ ℎ
Q T
 O 0 0
 O
ℎ 0T
l=Q  , t
 O u = t3 0u , V = M 

N = M 


N . t
O u
Q0 ℎ T 
O ℎ ℎ ℎ



 O 0 4 O
Qℎ 0T

6 6 6 6 6
P0 ℎ S

Reemplazamos:

ℎ rg rg
ℎ 0 R  ig U R i U
R0 ℎ U Qrj T
Qℎ rj T
g

w v Q T w v Q Qr g T
ℎ 0 T rig ℎ rig T
de = k k Q  . q g s 5?B #  B 56 5
= k k Q 5?B # B 56 5
= Q g T
i
0 ℎ T rj ℎ r T r
Q j T
w v Q
g
Q gT
w v  j
Qℎ 0T ℎ r Qr g T
P0 ℎ S Q  i T Q i T
Pℎ rj S Prj S
g g

Para determinar las funciones de interpolación, partimos de un elemento cuadrangular de

cuatro nodos. Las funciones de forma del elemento cuadrilateral son:

1 1
ℎ = 1 −
 1 − 6, ℎ = 1 +
 1 − 6
4 4
1 1
ℎ∗ = 1 +
 1 + 6, ℎ = 1 −
 1 + 6
4 4

Sumamos las funciones 3 y 4 para colapsar dos nodos en uno:

1
ℎ = ℎ∗ + ℎ = 1 + 6
2

De este modo, las funciones de forma definitivas para el elemento triangular son:

Para Y% Para YZ

1 1
ℎ = 1 −
 1 − 6 ℎ = 1 +
 1 − 6
4 4

ℎ 1 ℎ 1
= − 1 − 6 = 1 − 6

4 
4
ℎ 1 ℎ 1
= − 1 −
 = − 1 +

6 4 6 4

Para Y[

1
ℎ = 1 + 6
2

ℎ
=0


ℎ 1
=
6 2

Page 13 of 20
IMEF Facultad de Ingeniería - UBA

naturales
y 6.
Se concluye hasta este punto, que el operador Jacobiano es función de las coordenadas


O 1 1
− 1 − 6 1 − 6 0 0 0
V = M 

N = M 4 4 N.t u

O 1 1 1 3 0
− 1 −
 − 1 +
 0 4
6 6 4 4 2

O 3
1 − 6 0
V = M 

N = M 4 N

O 3
− 1 +
 2
6 6 4

Se debe evaluar entonces el operador Jacobiano en cada punto

, 6:

−1, −1 , 1, −1 , −1,1

3
V = t2 0u , det V  = 3
0 2
3
0
V = M 2 N, det V  = 3
3
− 2
2
0 0
V =  , det V  = 0
0 2

Ahora podemos calcular las componentes de las fuerzas concentradas en cada nodo, que
surgen como producto de la fuerza que actúa sobre el volumen del elemento:

rg
ℎ 0 R i U
R0 ℎ U Qrj T
g

w v Q T Qr g T
ℎ 0 T rig
de = k k Q  . q g s 5?B V  B 56 5
= Q g T
i
w v Q 0 ℎ T rj Qrj T
Qℎ 0T Qr g T
P0 ℎ S Q i T
Prj S
g

Page 14 of 20
IMEF Facultad de Ingeniería - UBA

de
  1 −
 1 − 6
Rk k q s 5?B V  B 56 5
0 U
Q $ $ 4 T
Q   1 −
 1 − 6 T
Q 0 k k q s 5?B V  B 56 5
T
4
Q   $ $
T
Qk k q 1 +
 1 − 6s 5?B V  B 56 5
0 T
Q 4 
T rig
= Q $ $  1 T . qr g s

+
 1 − 6
Q 0 k k q s 5?B V  B 56 5
T
j

Q $ $ 4 T
 1 + 6
Q 
T
Q k k q s 5?B V  B 56 5
0 T
2
Q T
$ $
  1 + 6
Q 0 k k q s 5?B V  B 56 5
T
P $ $ 2 S

Los valores de las integrales sobre cada función de forma son:

  1
k k ' 1 −
 1 − 6( 56 5
= 1
$ $ 4

  1
k k 1 +
 1 − 6 56 5
= 1
$ $ 4

  1
k k 1 + 6 56 5
= 2
$ $ 2

Finalmente, las componentes de las fuerzas nodales son:

rg
3 0 R i U 30
R0 U Q rj
g
T R15U
3
Q T Qr g T Q T
3 0T 10 30
de = Q .  = Q g T = Q T
i
Q0 3T 5 Qrj T Q15T
Q0 0T Qr g T Q 0 T
P0 0S Q i T P0S
Prj S
g

PROBLEMA 6

Se evaluará el vector de fuerza superficial dy para dos elementos.

Parte a)

El primero de ellos está descrito por la siguiente figura:

Page 15 of 20
IMEF Facultad de Ingeniería - UBA

dy = klz p fz p 5{ h
m

Para conformar la matriz de interpolación para los nodos de la superficie (línea) se requiere

influencia la carga uniformemente distribuída o fuerza -j .

evaluar las funciones de forma del elemento sobre la coordenada natural sobre la cual ejerce

1 1 1 1
ℎ = ' 1 +
 1 + 6( = 1 +
, ℎ = ' 1 −
 1 + 6( = 1 −

4 v 2 4 v 2

1 1
1 +
 0 1 −
 0 0 0 0 0
lz = M2 2 N
1 1
0 1 +
 0 1 −
 0 0 0 0
2 2
1
R 1 +
 0 U
Q2 T
1
Q 0 1 +
T
Q 2 T
Q1 1 −
 0 T
lz = Q2 T
m

Q 1 T
0 1 −

Q 2 T
Q 0 0 T
Q 0 0 T
Q 0 0 T
P 0 0 S

Definimos ahora el vector de fuerzas por unidad de area. En este caso, por unidad de longitud.

riz 0
fz = q z s =  
rj 5

transformación o det #z .
Luego, hallamos la corrección para el diferencial de longitud al introducir el Jacobiano de la

/
z z

 O 
5{ = 5/ = det # 5
, det # = q  +   s



Page 16 of 20
IMEF Facultad de Ingeniería - UBA



 −
 O O − O
= , =

2 
2

O
= 1, =0



det V z  = 1

Finalmente, evaluamos el vector de fuerzas superficiales:

dy = klz p fz p 5{ h
m

1 +
 0
R U
0 1 +

Q T
1 −
 0 T
 1Q
dy = k Q 0 1 −
T . 0 5

2
$ Q 0 0 T 5
Q 0 0 T
Q 0 0 T
P 0 0 S
0
U R ri U
z
R 
0 Qk 1 +
5
T Qrj T R0U
z
R 1 +
U 5
Q $ T Qr z T
Q 0 T Q 0 T Q i Q0T
 5Q T z T Q T
1 −
T 5Q  T rj 5
dy = k Q 5
= Qk 1 −
5
T = Q z T = Q T
$ 2 Q 0 T 2 $ Q ri T Q0T
Q T
Q 0 T Q 0 T Qrj T Q T
z 0
Q 0 T Q 0 T Q T Q0T
P 0 S rj
z
P
Q 0 T Q zT 0S
P r
S P j S
0

La fuerza resultante es de 10 unidades de fuerza en la dirección O. A cada nodo le corresponde

la mitad de la fuerza resultante debido a que la carga es uniformemente distribuída y es
simétrica.

Page 17 of 20
IMEF Facultad de Ingeniería - UBA

Parte b)

Se evaluará nuevamente el vector de fuerzas superficiales dy .

Nuevamente, evaluamos las funciones de forma sobre la línea sobre la cual está aplicada la
carga superficial.

1 1 1 1
ℎ = ' 1 +
 1 + 6( = 1 +
, ℎ = ' 1 −
 1 + 6( = 1 −

4 v 2 4 v 2

1 1
ℎ = ' 1 −
 1 − 6( = 0, ℎ = ' 1 +
 1 − 6( =0
4 v 4 v

1 1
1 +
 0 1 −
 0 0 0 0 0
| = M2 2 N
1 1
0 1 +
 0 1 −
 0 0 0 0
2 2
1
R 1 +
 0 U
Q2 T
1
Q 0 1 +
T
Q 2 T
Q1 1 −
 0 T
| = Q2
m T
Q 1 T
0 1 −

Q 2 T
Q 0 0 T
Q 0 0 T
Q 0 0 T
P 0 0 S

Definimos el vector de carga superficial como:

riz −- . 6?3 } −2.2361

r z = q zs = ' ; (= 
rj -; . cos } 4.4721

Page 18 of 20
IMEF Facultad de Ingeniería - UBA

Donde } es el ángulo que forma la vertical con -; , y también el mismo ángulo que forma la
horizontal con la línea entre los nodos 2 y 1. En este caso, el vector r z posee unidades de
fuerza sobre unidad de longitud.

1
tan } =
2

transformación o det #z .
Hallamos también la corrección para el diferencial de longitud al introducir el Jacobiano de la

/
z z

 O 
5{ = 5/ = det V 5
, det V = q  +   s



 O 2 3

 O 0 2
M N=M N

 O 0 0

 O 2 0



 −
 O O − O

= , =

2 
2

2 − 0 O 3 − 2 1
= = 1, = =

2 
2 2

1 / 5 / √5
det V z  = '1 + ( = ' ( =
4 4 2

Procedemos a calcular el vector de fuerza superficial:

dy = klz p fz p 5{ h
m

1 +
 0
R U
0 1 +

Q T
1 −
 0 T
 1Q
dy = k Q 0 1 −
T . −2.2361 √5 5

2
$ Q 0 0 T 4.4721 2
Q 0 0 T
Q 0 0 T
P 0 0 S

Se calculan las integrales:


k 1 +
5
= 2
$


k 1 −
5
= 2
$

El vector resultante posee componentes en
e O, como se esperaba.

Page 19 of 20
IMEF Facultad de Ingeniería - UBA

rz
R i U
2 0 rj
z −2.5
R0 Q T R U
2U Qri T Q 5 T
Q2 0 T z
Qr z T Q−2.5T
√5 Q0 T
2T −2.2361 5 T
dy = Q . =Q zT=Q
j
4 Q0 0T 4.4721 r
Q i T Q 0 T
Q0 0T Qrjz T Q 0 T
Q0 0T Q zT Q 0 T
P0 r
0S Q j T P 0 S
Prj S
z

Page 20 of 20