CAPÍTULO 5 Vigas R. Cueva P.

97

PRELIMINARES
La VIGA es una barra que recibe la carga mecánica denominada FLEXIÓN; en las figuras siguientes se
muestra su significado.
a) TIPOS DE FLEXIÓN

 Simple o pura:

 Con cortante:

b) SISTEMA DE COORDENADAS
La fibra neutra es una fibra longitudinal de la viga, la cual se curva únicamente, sin alargarse o
acortarse debido a la deformación de la viga; la fibra neutra está ubicada en el centroide de la
sección de la viga. El eje de coordenadas (X) coincide con la fibra neutra.

Un punto genérico de la viga, tendrá las siguientes coordenadas (antes de la deformación de la

viga):
( )
Después de la deformación de la viga, sus coordenadas serán:

( )

En la fibra neutra tendremos lo siguiente:

 R. Cueva P. 98
Resulta lo siguiente:
( )
En forma de función explícita:
( ) DEFORMADA DE LA VIGA

c) ANALISIS DE LA FLEXIÓN PURA

En la figura se muestra un segmento de la viga sometida a flexión pura; en la fibra longitudinal que
está separada una distancia (y) de la fibra neutra (F.N.).
En una sección genérica situada a una distancia (X) del origen de coordenadas (la sección gira
permaneciendo plana); tendremos lo siguiente:

La deformación en la fibra genérica (en color verde), será:

Por semejanza de triángulos:

( )

El desplazamiento resulta:

La deformación resultará:
(1)

El esfuerzo (tracción simple), debe ser:

(2)
 R. Cueva P. 99
La curvatura de la fibra neutra, está dada por la siguiente expresión:

⁄
*( ) +

Si las pendientes son pequeñas como ocurre dentro de la zona elástica de deformación de la viga;
lo anterior se reduce a:

(3)

Reemplazando (3) en (2):

(4)

Las condiciones de equilibrio en la sección de la viga, serán las siguientes:

 DE LOS MOMENTOS:
∫ (5)

 DE LAS FUERZAS NORMALES A LA SECCIÓN DE LA VIGA:

∫ (6)

De lo anterior, tendremos lo siguiente:

De (4) en (6):

∫ ∫

Según lo anterior, la fibra neutra ( ) debe pasar obligatoriamente por EL CENTROIDE DE LA

SECCIÓN DE LA VIGA.

De (4) en (5):

La integral es igual al momento de inercia de la sección respecto al centroide y aun eje ortogonal al
plano de la figura (en el papel); es decir ortogonal al PLANO DE LA FLEXIÓN de la viga.
De esta última expresión y de (4), tendremos lo siguiente:

(A)
 R. Cueva P. 100
De esta última expresión y de (3), tendremos lo siguiente:

(B)

d) FLEXIÓN CON CORTANTE

En la figura siguiente se muestra una viga con cargas transversales que están produciendo una
flexión con fuerza cortante; las cargas internas en el caso de flexión con cortante, son: el momento
flector y la fuerza cortante. En el caso de la flexión pura, la carga interna es el momento flector
únicamente.

(x)

El equilibrio de un elemento de longitud de la viga, debe ser satisfecho.

 R. Cueva P. 101

 DE LAS FUERZAS VERTICALES:

( )

(C)

 DE LOS MOMENTOS:
( )

(D)

Mediante (B), (C) y (D); se obtienen LAS ECUACIONES DE LA FLEXIÓN:

 EL ESFUERZO CORTANTE EN LA FLEXIÓN:

En el elemento de longitud de la viga, el equilibrio de las fuerzas horizontales, también debe
ser satisfecho:

El equilibrio de las fuerzas horizontales:

 R. Cueva P. 102

∫ ∫

En donde:

( )

( )

Reemplazando en lo anterior y ordenando términos, resulta:

Mediante (D); en la expresión anterior resulta lo siguiente:

e) LA ENERGIA ELASTICA EN UNA VIGA

En una viga sometida a flexión simple, se cumple lo siguiente:

 R. Cueva P. 103

En la figura se muestra el significado de la expresión anterior. En un elemento diferencial de

longitud de la viga , se tendrá:

En una viga cualquiera (flexión con cortante), el momento flector será variable a lo largo de (x); se
tendrá entonces:
∫
Mediante (B):

De lo cual se obtiene:

Reemplazando en la integral anterior, se obtiene lo siguiente:

∫ (E)

En función de la curvatura, según (B):

∫ ( ) (E’)

1 MODELADO DE UNA VIGA MEDIANTE ELEMENTOS FINITOS

a) GLOBAL:

Los grados de libertad:

[ ]
 R. Cueva P. 104

Las cargas nodales:

[ ]

En donde:

b) LOCAL:

Los grados de libertad:

[ ]

Las cargas nodales:

[ ]

c) CONECTIVIDAD:

--- --- --- ---

2 DESPLAZAMIENTO DE LOS PUNTOS INTERNOS DEL ELEMENTO FINITO

El desplazamiento de un punto interno en el elemento finito, se obtiene mediante la
INTERPOLACIÓN DE LOS VALORES NODALES (los grados de libertad nodales).
a) LAS FUNCIONES DE HERMITÉ
En el elemento finito genérico, se crea un vector de grados de libertad ficticios UNITARIOS Y
POSITIVOS para establecer una proporción con respecto a los grados de libertad reales:

[ ]
 R. Cueva P. 105

Observando las deformadas que cada uno de estos grados de libertad nodales UNITARIOS Y
POSITIVOS, producen por separado en el elemento finito; las DEFORMADAS parciales resultantes
en el elemento finito, se pueden definir mediante un polinomio de ( ) de tercer grado (por lo
menos):

La forma que adquiere el elemento finito; debido a

Con las condiciones de contorno correspondientes en el vector desplazamiento:

[ ]

El polinomio que grafica la deformada del elemento finito resulta:

( )

Debido a

Cuyas condiciones de contorno son:

[ ]

El polinomio que grafica la deformada del elemento finito resulta:

( )
 R. Cueva P. 106
Debido a

Cuyas condiciones de contorno son:

[ ]

El polinomio que grafica la deformada del elemento finito resulta:

( )

Debido a

Cuyas condiciones de contorno son:

[ ]

El polinomio que grafica la deformada del elemento finito resulta:

( )

Por proporción con los valores reales de los grados de libertad:

1
 R. Cueva P. 107

Igualmente en las funciones de interpolación y ; SUPERPONIENDO los cuatro gráficos

parciales, tomando en cuenta además los valores reales de los grados de libertad, resulta:

El gráfico de la figura, tiene diferencias con la deformada de un elemento finito que modela una
viga (ver la primera figura de la página 103); sin embargo, al cambiar de signo al grado de
libertad y al gráfico de el diagrama de la deformada se corrige a la forma esperada.

Dando como resultado la DEFORMADA DEL ELEMENTO FINITO, respecto a su propia coordenada.

La expresión anterior corresponde al desarrollo del siguiente producto escalar:

b) LAS COORDENADAS GLOBALES EN FUNCIÓN DE LAS COORDENADAS LOCALES

En el capítulo 3 se dedujo la siguiente relación:

También se puede escribir:

El primer término es la coordenada del centro del elemento finito; el segundo término es el
jacobiano o factor de escala:
 R. Cueva P. 108
c) TRANSFORMACIÓN DE LA PENDIENTE
Mediante la regla de la cadena para la derivación:

De lo cual se obtiene:

d) DEFORMADA DEL ELEMENTO FINITO

Los grados de libertad locales coinciden con los globales porque el origen de coordenadas local está
siempre sobre el eje (X), únicamente cambia el factor escala de las pendientes; efectuando dicho
cambio, tendremos lo siguiente:

( )

Teniendo en cuenta este factor de escala de las pendientes, se podrá escribir:

( )

e) DEFORMADA DE LA VIGA
Según la conectividad del modelo, se podrá graficar la deformada de la viga con los grados de libertad
globales, utilizando el recuadro anterior para graficar dentro de cada elemento finito.

3 MATRICES DE RIGIDEZ LOCALES

La energía elástica en función de la curvatura de la viga según el recuadro (E’) (página 103), es lo
siguiente:

( ) ∫( )

En donde, la expresión numérica de la deformada, es:

Cuyo desarrollo, interpretando los grados de libertad como flechas y pendientes en la viga; es el
siguiente:

( ) ( ) ( ) ( )
( ) ( )
 R. Cueva P. 109
La primera derivada de la flecha respecto a la variable global, teniendo en cuenta el cambio de escala
con respecto a la variable local, es la siguiente:

La segunda derivada:

[ ]

Mediante la expresión numérica de la deformada:

( )

Elevando al cuadrado la segunda derivada de la flecha:

( ) ( ) *( )( )+

En donde:

* +

Resulta:

*( ) ( ) ( ) ( ) +

Reemplazando en la expresión de la energía (en la página 108):

( ) ∫ ( )

En donde:

Desarrollando:

( ) ( ) [( ) ( ) ( ) ( ) ]
 R. Cueva P. 110
El producto siguiente resulta:

( )

*( ) ( ) ( ) ( )

( ) ( ) ( ) ( )

( ) ( ) +

Ejecutando y acotando la siguiente integral:

∫ ( )

[
]

Reemplazando en la expresión de la energía elástica (el recuadro anterior); se obtiene lo siguiente:

( )[

La matriz de rigidez del elemento finito, será:

* +

[ ]
Reemplazando el recuadro anterior en esta última matriz; se obtiene:

( )
[ ]
 R. Cueva P. 111

Los grados de libertad de un elemento finito que modela una viga, son de traslación (mm) y de
rotación (adimensional); por esta razón, las unidades (ecuación dimensional) de las componentes de
la matriz de rigidez (recuadro anterior) no son iguales entre sí.

4 LA MATRIZ DE RIGIDEZ DEL MODELO

∑ { }

5 EL VECTOR DE CARGAS EXTERNAS

La fuerza distribuida en una viga es bastante común; ya sea como peso propio o como una carga
externa aplicada sobre la viga; sin embargo, las cargas externas sobre el modelo deben ser NODALES.
Se puede hacer el cambio de las fuerzas externas distribuidas a concentradas nodales, manteniendo
la misma ENERGIA ELÁSTICA para que el modelo sea equivalente.

La equivalencia energética entre ambos sistemas de cargas externas, es la siguiente: el trabajo

externo del primer sistema debe ser igual al trabajo externo del segundo sistema:

∫ [ ∫ ]

Ejecutando y acotando las cuatro integrales, se obtiene lo siguiente:

* +[ ]

De lo cual se deduce que el vector carga en el elemento finito, debe ser:

* +

En la figura siguiente se muestra los diagramas de momento flector (superpuestos) de los dos
sistemas de cargas externas, equivalentes entre sí: una flexión transversal con fuerza distribuida, y
una flexión pura con los momentos externos y en los extremos.
 R. Cueva P. 112

Las fuerzas de reacción en los apoyos deben ser iguales entre sí, por simetría y con el siguiente valor:

Las áreas bajo la parábola y del rectángulo, deben ser iguales entre sí, para que exista la equivalencia
de la energía; obsérvese sus dimensiones y compárese con el recuadro de la página anterior.

Las cargas concentradas, ya sea fuerzas o momentos, deben actuar directamente en los nodos del
modelo.

6 LA ECUACIÓN DE RIGIDEZ

7 LA FUERZA CORTANTE Y EL MOMENTO FLECTOR

Según las ecuaciones de la flexión (página 101); la segunda y la tercera ecuación, tenemos lo
siguiente:

 La relación de la cortante con la deformada de la viga:

 La relación del momento flector con la deformada de la viga:

Teniendo en cuenta que la expresión numérica de la deformada de la viga, dentro de cada elemento
finito, es:

Tendremos las siguientes expresiones numéricas de la fuerza cortante y del momento flector:

( ) ( )[ ]

( ) ( )[ ( ) ( ) ]
 R. Cueva P. 113

La expresión de la fuerza cortante es constante (independiente de ); la expresión del momento

flector es de primer grado en (distribución lineal), dentro del elemento finito. En las dos figuras
siguientes se muestran sus gráficos respectivos.

Significa que el momento flector máximo, estará en alguno de los dos nodos del elemento finito:

En el EJEMPLO SIGUIENTE se muestra una viga con una fuerza distribuida uniformemente, de sección
constante, el material constante y con dos apoyos simples.

Se ha modelado la viga con cuatro elementos finitos a título de ejemplo; en cada diagrama se observa
la forma analítica (con trazo fino) y la forma numérica (con trazo grueso), superpuestas. Al aumentar
el número de elementos finitos (afinado de la malla), ambas formas tienden a una mayor
coincidencia entre sí. En la figura siguiente se observa los diagramas de la fuerza cortante y del
momento flector.
 R. Cueva P. 114

8 LOS ESFUERZOS EN UNA VIGA

a) EL ESFUERZO DE FLEXIÓN (tracción y compresión):

En forma analítica:

En forma numérica:

( )[ ( ) ( ) ]

b) EL ESFUERZO CORTANTE (cizallante):

Proviene de la condición de equilibrio de las fuerzas en la dirección horizontal, de un elemento de
longitud de la viga, cuyo resultado del análisis es el siguiente:

∫
 R. Cueva P. 115

El valor máximo de que se encuentra a la altura del eje (X); en forma analítica está expresado de
la manera siguiente:

En forma numérica:

( )[ ]

En el caso de una sección rectangular, la distribución de esfuerzos cortantes, será:

( )* +

El valor máximo, ocurrirá en ( ):

( )

El coeficiente del corte por cizallamiento, para la sección rectangular, será:

 R. Cueva P. 116

En una sección circular de radio , la distribución del esfuerzo cortante, será:

( )* +
El coeficiente del corte por cizallamiento, para la sección circular, será:

9 EJEMPLO DE APLICACIÓN
Determinar las pendientes de la viga en los apoyos simples, la flecha en el punto medio de la zona
con la carga distribuida, y la flecha máxima y su ubicación. El material y la sección, son los siguientes:

SOLUCIÓN
a) MODELADO Y GRADOS DE LIBERTAD
Mediante dos elementos finitos, con seis grados de libertad (como corresponde); tal como se
muestra en la figura siguiente:

b) LA FUERZA DISTRIBUIDA CONVERTIDA A CONCENTRADAS NODALES

 R. Cueva P. 117
En donde:

N-mm

N-mm

c) LAS MATRICES DE RIGIDEZ LOCALES

El factor numérico:

( )( )
( )

( )
( )[ ]

3 4 5 6

( )
( )[ ]

d) LA MATRIZ DE RIGIDEZ DEL MODELO

( )

[ ]

e) EL VECTOR CARGA EN EL MODELO

 R. Cueva P. 118

Resulta lo siguiente:

[ ( ) ( ) ]

f) LAS CONDICIONES DE CONTORNO PARA LOS GRADOS DE LIBERTAD

[ ]

g) LA ECUACIÓN DE RIGIDEZ

Resolviendo el subsistema:

* + ( )* + [ ]

Resulta:

h) LAS REACCIONES EN LOS APOYOS

El vector carga se puede escribir:

( )

En la ecuación de rigidez:

i) LA DEFORMADA

Dentro de cada elemento finito:

 R. Cueva P. 119
Cuyo desarrollo es:

[( ) ( ) ( )

( ) ]

La flecha en el punto medio del elemento finito (2); en sus coordenadas particulares, ( ):

( )
[ ]

Según la conectividad (de la segunda figura de la página 116), se tendrá lo siguiente:

Reemplazando en la expresión anterior:

j) LA PENDIENTE
Con respecto a las coordenadas locales:

[ ( ) ( ) ( ) ( ) ]

La flecha máxima en el elemento finito (1):

( ) {

La flecha máxima en el elemento finito (2):

[ ( ) ( ) ] {

k) EL ESFUERZO
En cada elemento finito:

( )[ ( ) ( ) ]

En el elemento finito (1):

( )
( ) [( ) ] {
 R. Cueva P. 120
En el elemento finito (2):

( )
( ) [( ) ( ) ] {

10 FLEXIÓN EN EL ESPACIO
El área de la sección de una viga, cuya forma es una característica propia; en las figuras siguientes se
muestran dos (a título de ejemplo). La sección de la viga se encuentra en el plano formado por las
coordenadas (Y, Z) y la longitud de la viga en la coordenada (X).
El momento de inercia del área de la sección, es una característica que tiene que ver con el
comportamiento de la viga en la flexión; el valor de este momento de inercia depende del eje que
tenga, este eje a la vez, debe ser ortogonal al plano de flexión de la viga.
El momento de inercia tiene una naturaleza DIÁDICA; en el plano (Y, Z), tendremos lo siguiente:

[ ]

En las coordenadas principales (que tiene toda diádica):

[ ]

Los ejes de simetría (si la sección tiene simetría) coinciden con las coordenadas principales, porque
debido a la simetría, los productos de inercia de la sección se anulan.

Observemos el esfuerzo normal debido a la flexión, en las dos situaciones; tomando en cuenta que el
plano de la flexión es vertical: en el plano (Y, X).

En la figura de la derecha, NO EXISTE COINCIDENCIA de las coordenadas (Y, Z) con los ejes principales
(1, 2):

( ) ( )
 R. Cueva P. 121

La curvatura de la viga, en ambos planos:

( )
( )

( )
( )
En la figura de la izquierda, SI EXISTE COINCIDENCIA de las coordenadas (Y, Z) con los ejes principales
(1, 2):

Observando la curvatura de la viga, en ambos planos:

( )

( )

En este segundo caso (figura de la izquierda), las flexiones en ambos planos SON INDEPENDIENTES
ENTRE SÍ; lo cual se debe a la ausencia del producto de inercia de la sección de la viga ( ).

Cuando una carga de flexión no se encuentra en una dirección principal del momento de inercia de la
sección de la viga, entonces dicha carga se debe descomponer en las direcciones principales.

(Fin del capítulo 5)