Autor: Christian Picón M.

Análisis Matricial de Estructuras 11 de mayo de 2020

Tarea No3
Pregunta 1:
Determinar las funciones de forma y la Matriz de Rigidez para el elemento barra con 4 nodos en la configuración
libre – libre. Además, considerar una variación lineal de la sección de la barra de A(x = 0) = 3A a A(x=L) = A y
considere una interpolación para el desplazamiento en la barra de la forma:

u(x) = a + bx + cx2 + dx3

u1 u2 u3 u4

3A A
1 2 3 4

E, A(x)

L/3 L/3 L/3

x
Resolucion:
Aplicando las condiciones de borde del problema a la función del campo de desplazamiento se determinan las
constantes a, b, c y d en función de los desplazamientos ui.

Por lo tanto,
x = 0, u(x) = u1 ⇒ u1 = a
L bL cL2 dL3
x= , u(x) = u2 ⇒ u2 = u1 + + +
3 3 9 27
2L 2bL 4cL2 8 dL3
x= , u(x) = u3 ⇒ u3 = u1 + + +
3 3 9 27
x = L, u(x) = u4 ⇒ u4 = u1 + bL + cL2 + dL3

Resolviendo el sistema de ecuaciones de encima para b, c y d se tiene:

11 u1 − 18 u2 + 9 u3 − 2 u4
b=− (1)
2L
18 u1 − 45 u2 + 36 u3 − 9u4
c= (2)
2L2
9u1 − 27 u2 + 27 u3 − 9u4
d=− (3)
2L3

Los resultados de (1), (2) y (3) se reemplazan en la función del campo de desplazamientos obteniendo:

(11 u1 − 18 u2 + 9 u3 − 2 u4) x 9 (2 u1 − 5 u2 + 4 u3 − u4) x2 9 (u1 − 3 u2 + 3 u3 − u4) x3

u(x) = u1 − + −
2L 2L2 2L3

1/2
Autor: Christian Picón M. Análisis Matricial de Estructuras 11 de mayo de 2020

Agrupando en términos semejantes se obtiene la función de forma de cada grado de libertad:

9x2 9x3 45x2 27x3

   
11x 9x
u(x) = 1− + 2 − · u1 + − + · u2
2L L 2L3 L 2L2 2L3
18x2 27x3 9x2 9x3
   
9x x
+ − + − · u3 + − + · u4
2L L2 2 L3 L 2L2 2L3

9x2 9x3
 
11x
N1 (x) = 1 − + 2 −
2L L 2L3
45x2 27x3
 
9x
N2 (x) = − +
L 2L2 2L3
18x2 27x3
 
9x
N3 (x) = − + −
2L L2 2 L3
2 3
 
x 9x 9x
N4 (x) = − 2 +
L 2L 2L3
Cada función de norma evaluada en la distancia de su correspondiente grado de libertad da un valor de 1; y 0 si es
la distancia de otro grado de libertad.

Una vez encontradas las funciones de forma N (x), en base a la energía interna de deformación (4) se calcula la
matriz de rigidez aproximada de la barra.
Z L  T  
1 dN (x) dN (x)
U= AE dx (4)
2 0 dx dx

Donde: N (x)1×4 = [N1 (x) N2 (x) N3 (x) N4 (x)]

Derivando N (x)

dN (x) h 11 18x 27x2 9 45x 81x2 9 36x 81x2 1 9x 27x2

i
= − + − − + − 2L + − − +
dx 2L L2 2L3 L L2 2L3 L2 2L3 L L2 2L3

Resolviendo el término dentro de la integral de la ecuación (4) se obtiene la matriz de rigidez, que para este ejem-
plo la sección transversal es variable obedeciendo a la ecuación A(x) = 3A − 2A L x.

La ecuación (5) muestra la resolución para el primer término de la matriz de rigidez.

L
27 2 2
Z   
2A 11 18 41A
E 3A − x · − + 2x − 3
x dx = (5)
0 L 2L L 2L 4L

De igual manera, integrando con la misma ecuación de área y para cada término, se obtiene la matriz de rigidez
aproximada final:

41 A
− 129 A 33 A
− 13 A
 

 4L 10 L 10 L 20 L 

− 129 A 513 A
− 297 A 21 A 

 
 10 L 20 L 20 L 10 L 
[KL ] = E  
33 A
− 297 A 351 A 6A 
− L 



 10 L 20 L 20 L
− 13 A 21 A 6A 91 A
 
20 L 10 L − L 20 L

2/2