UNIVERSIDAD CENTRAL DEL ECUADOR

CARRERA DE INGENIERIA CIVIL

FACULTAD DE INGENIERÍA, CIENCIAS, FÍSICA Y MATEMÁTICA
ECUACIONES DIFERENCIALES

ASIGNATURA: Ecuaciones Diferenciales

PROFESOR: Ing. Pilaluisa Ramiro
PERÍODO ACADÉMICO: Abril 2015 – Agosto 2015

INFORME DE INVESTIGACIÓN

TÍTULO:

FLEXION DE UNA VIGA EN VOLADIZO PARA PEQUEÑAS FLEXIONES

FECHA DE ENTREGA: 17 de julio de 2015

MIEMBROS DEL GRUPO

1. ARMIJOS MACASLUIS ESTEBAN 100%

2. CACHAGO AIGAJE SEBASTIAN 100%

3. LLUMIGUANO DARIS EDISON 50%

4. SORIA CAIZA DARWIN SANTIAGO 50%

Ecuaciones Diferenciales Proyecto Final

Tabla de contenido
1 Resumen..................................................................................................................................2
2 Introducción.............................................................................................................................2
3 Materiales y Métodos..............................................................................................................4
4 Resultados................................................................................................................................4
5 Discusión.................................................................................................................................2
6 Conclusiones............................................................................................................................3
7 Recomendaciones....................................................................................................................3
8 Referencias..............................................................................................................................3
9 ANEXOS.................................................................................................................................4

Prof. Ing. José Ramiro Pilaluisa Q. M.Sc. Abril 2015 – Agosto 2015 1
Ecuaciones Diferenciales Proyecto Final

1 Resumen

En el presente trabajo se estudiara la deflexión de una viga en voladizo cuando aplicamos una carga,
ya sea, carga por peso propio de la viga que sería la carga muerta la cual se representa como una
carga uniformemente distribuida o la carga viva representada con una carga puntual, para el diseño
utilizaremos el software SAP 2000 para hacer la simulación de la aplicación de las cargas tanto
distribuidas uniformemente como la carga puntual, obteniendo así la deflexión de la viga en el
voladizo, y si obtenemos un valor muy alto rediseñaremos la viga para obtener una deflexión
pequeña.
Palabras Claves
Deflexión, voladizo, carga, viga.

Summary

In this paper the deflection of a beam will be study when we apply a load, either , load own weight
of the beam would be the dead load which is represented as a uniformly distributed load or live load
represented with a point load, to use the SAP 2000 design software for simulating the application
loads evenly distributed both as point load , thus obtaining the deflection of the cantilever beam ,
and if we get a very high value will redesign the beam to obtain a small deflection .
KeyWords
Deflection, cantiléver, load, beam.

2 Introducción

Una viga o una barra delgada son sólidos homogéneos e isótropos cuya longitud es grande
comparada con las dimensiones de su sección trasversal.

Figura 1: Flexión de una viga en voladizo de sección rectangular

Fuente: España, Estudio de la Flexión de una Viga de Material Elástico no Lineal, Departamento de
Ciencia y Tecnología de los Materiales
Universidad Miguel Hernández de Elche(Alicante), (2002).

Cuando una viga flexiona debido a las fuerzas exteriores que se aplican, existen algunas partes de la
viga que se acortan y hay otras zonas que se alargan. Pero hay una línea, denominada eje neutro,
Prof. Ing. José Ramiro Pilaluisa Q. M.Sc. Abril 2015 – Agosto 2015 2
Ecuaciones Diferenciales Proyecto Final

que no se acorta ni se alarga. Este eje se encuentra en el centro de gravedad de la sección trasversal.
Se usará una barra empotrada de un determinado material, de longitud L, de anchura a y de espesor
b. Se fijará uno de sus extremos y se aplicará una fuerza en su extremo libre. Mediremos el
desplazamiento del extremo libre y(L) o flecha en función de la fuerza aplicada F, comprobando su
relación de proporcionalidad, mientras que la flexión de la barra sea pequeña.
A continuación, calcularemos el desplazamiento de su extremo libre cuando se aplica una fuerza en
dicho extremo que produce una flexión considerable.
Este ejemplo, nos permite practicar con procedimientos numéricos aplicados al
 Cálculo de la raíz de una ecuación.
 Integral definida.

Figura 2: Perfil de una barra empotrada en un extremo y sometida a una fuerza que actua sobre el
extremo libre.

Fuente: Beléndez T., Neipp C., Beléndez A., Large and small defections of a cantilever
beam. (2002).

Supongamos que
 La barra tiene una longitud L mucho mayor que las dimensiones de su sección trasversal, y
que la deformación debida a su propio peso es despreciable.
 Que la sección de la barra no cambia cuando se dobla. Cuando el espesor de la barra es
pequeño comparado con el radio de curvatura, la sección trasversal cambia muy poco.
En estas condiciones es aplicable la ecuación de Euler-Bernoulli que relaciona el momento flector
M de la fuerza aplicada y el radio de curvatura ρ de la barra deformada

YI
M=
ρ
Donde: Y =modulo de elasticidad de Young (depende delmaterial )
I =momento de inercia de la secciontransversal de la viga

Tomando con respecto a una línea horizontal que pasa por el centro de gravedad de la sección:
El producto EI se llamarigidez a la flexion y es una constante .
ρ=es el radio de curvatura elastica de la seccion con ecuacion:
2 3

ds
ρ= =
[ ( )]
1+
dy
dx
2

dθ d2 y
d x2
2
dy
Para pequeñas pendientes ( )
dx
≈0

1 d2 y
=
ρ d x2

Prof. Ing. José Ramiro Pilaluisa Q. M.Sc. Abril 2015 – Agosto 2015 3
Ecuaciones Diferenciales Proyecto Final

Si despreciamos el peso de la propia barra, el momento de la fuerza F aplicada en el extremo libre,

respecto del punto P (x, y) es M =F ( x f −x ) ≈ F (L− x)
3
y=
FL
2Y . I (
x 2−
x
3L )
Que integramos dos veces con las siguientes condiciones iníciales x=0, y=0, dy/dx=0.
El desplazamiento y f del extremo libre x=L es proporcional a la fuerza F aplicada
3 3
L ab
y= F I=
3Y . I 12

Prof. Ing. José Ramiro Pilaluisa Q. M.Sc. Abril 2015 – Agosto 2015 4
Se considera que la aproximación de pequeñas flexiones: el desplazamiento y del extremo libre de
la barra, es proporcional a la fuerza F aplicada, produce resultados aceptables hasta un cierto valor
del parámetro a dimensional α<0.375, (véase en anexos el ejercicio) o bien, hasta un valor máximo
de la fuerza aplicada Fm=2Y·I·α/L2

Hipótesis

 Con este procedimiento obtendremos la deflexión de la viga.

 Podremos hacer una comparación con las deflexiones de las vigas apoyadas en sus
extremos.
 Análisis si la viga propuesta resiste o no las cargas.

Solución del problema

  Si la viga no resiste las cargas podemos hacer un nuevo dimensionamiento para que la
misma no se fleje.
  Aumentar la resistencia del hormigón.
  Aumentar el acero de refuerzo.

3 Materiales y Métodos

El material a emplearse para el diseño de la viga será hormigón armado, el mismo que consta de
hormigón simple y acero lo que le da el nombre de hormigón armado, el método utilizado será la
medición de la flexión con el programa SAP 2000 y el método de demostración será una maqueta
de un pórtico con voladizo para ver la flexión de forma que pueda ser apreciada a simple vista.

4 Resultados

 Sea L=30 cm=0.3 m, la longitud de la barra.

 Sea b=0.78 mm=0.00078 m, el espesor de la barra.

 La anchura a=0.03 m está fijada por el programa interactivo y no se puede cambiar.

 Elegimos como material, el Acero.

Figura 3: Barra delgada empotrada en un extremo y con una fuerza vertical en el extremo
libre.
 Fuente: España, Flexión de Una Barra Delgada Empotrada en un Extremo: Aproximación
para Pequeñas Pendientes, Departamento de Ciencia y Tecnología de los Materiales,
Universidad Miguel Hernández de Elche(Alicante), (2002).

Después de realizar la experiencia. La pendiente de la recta que relaciona la desviación del extremo
libre y (L) con la fuerza aplicada F en dicho extremo es:

m=3.683 [ ]
cm
N
=0.03683
m
N [ ]
 El momento de inercia I que se obtiene es:

1
I= ×0.03 × ( 0,00078 )2 =1.186× 10−12 [ m2 ]
12
 Dada la pendiente ( coeficiente de proporcionalidad de F) calculamos el módulo de Young
(Y)
( 0.3 )3
0.03683=
3 ×Y × ( 1.186 × 10−12)

Y =2.06 × 10
11
[ ]
N
m
2

5 Discusión

Este valor obtenido corresponde al modulo de Young del acero, y se lo puede comparar con los
datos teoricos que encontramos en los textos descritos en las referencias,entonces obtuvimos:
Y (experimental )=2.06 ×1011 2
N
m [ ]
Que esta correcto ya que el valor que encontramos en los textos está cercano a:
Y (teorico)=200× 10
9 N

m
2 [ ]
Lo que nos indica que la ecuaciones que se han plateado están correctamente determinadas.
En análisis de las mismas se encuentra ubicado en la parte de anexos del presente informe.
6 Conclusiones

 Partiendo de la hipótesis de que el material del que está fabricada una barra delgada
flexionada es elástico y se toma la barra como un elemento lineal, se ha presentado la
ecuación diferencial de la elástica de la barra de forma análoga a como se hace, esta
ecuación diferencial se ha simplificado considerando el caso de pequeños desplazamientos
de la directriz de la barra (pendientes pequeñas), pero sin hacer una aproximación tan
extrema a como se hace en los textos universitarios de Física y Mecánica.
 Un aspecto importante y de gran interés didáctico es la forma en como se ha aplicado la
condición de que la barra es inextensible (la longitud de su fibra neutra no cambia), aspecto
que no se considera en los libros de Física General.
 Se han obtenido las ecuaciones de la elástica y los desplazamientos vertical y horizontal del
extremo libre de la barra, comparándolos con los correspondientes al caso general de
pequeñas y fuertes pendientes.
 En este último caso se supone que el desplazamiento horizontal del extremo libre de la barra
es nulo.

7 Recomendaciones

 Se recomienda que se trabaje con mucha precisión y cuidado al momento de realizar las
lecturas experimentales ya que es este momento donde más se cometen erros que pueden
deberse a algunos factores como: la pericia del individuo que realiza la práctica, la calidad
de los equipos utilizados, las condiciones de trabajo, entre otras, así que es mejor buscar un
lugar y momento adecuado para realizar la práctica.
 Es recomendable y muy conveniente apoyarse con algunos software para realizar el análisis
de los datos leídos en la experimentación, ya que siempre podemos confiar en las
computadores ya que no cometen errores con en los que incurre el ser humano por esa
misma condición.
 Se recomienda empaparse de información sobre el tema presentado previo a la lectura de
este trabajo ya que siempre quedaran dudas sobre algunos temas o términos que no se
comprendan bien ya que la información es específica y puntual sobre el tema, y no se ha
adentrado en el contexto que se necesita para adentrarse en el estudio de deflexiones
pequeñas en vigas o barras.

8 Referencias

Beer F., Russell E., (2004). Mecánica de Materiales (3 ed.). Mexico: Mc-Graw-Hill.

Bedford A. y Fowler W., Mecánica para Ingeniería: Estática. (1987). U.S.A: Addison-Wesley
Iberoamericana.

Hibbeler R. C., (1998), Mecánica de Materiales. México: Prentice Hall.

Stiopin P., Gutierres A., (1976). Resistencia de Materiales (2 ed.). Moscú: Mir.

Pijaudier-Cabot. G., (2013). Damage Mechanics of Cementitious Materials and Structures. U.S.A:
John Wiley & Sons. Recuperado el 24 de julio de 2015 de:
http://bvirtual.uce.edu.ec:2054/lib/bgeneraluce/reader.action?docID=10660549
Khanna S., Jenkins C., (2005). Mechanics of Materials : A Modern Integration of Mechanics and
Materials in Structural Design. U.S.A: Academic Press. Recuperado el 24 de julio de 2015 de :
http://bvirtual.uce.edu.ec:2054/lib/bgeneraluce/reader.action?docID=10167019

Beléndez T., Neipp C., Beléndez A., Flexión de una barra delgada empotrada en un extremo:
Aproximación para pequeñas pendientes. Revista Brasileira de Ensino de Física. 24 (4) Dezembro
2002, págs, 399- 407. Recuperado el 24 de julio de 2015 de: http://www.scielo.br/scielo.php?
pid=S1806-11172002000400006&script=sci_arttext

Guerrero, N. L., Análisis de vigas tubulares de acero estructural propensas a la falla por pandeo
local sometidas a flexión biaxial. Rev. Fac. Ing. UCV, Dic 2006, vol.21, no.4, p.17-21. ISSN 0798-
4065. Recuperado el 24 de julio de 2015 de: http://www.scielo.org.ve/scielo.php?pid=S0798-
40652006000400002&script=sci_arttext

Kimiaeifar A., Tolou N., Barari A., & Herder J.L., (01 August 2013), Large deflection analysis of
cantilever beam under end point and distributed loads. Recuperado el 24 july 2015 de:
http://bvirtual.uce.edu.ec:2053/doi/full/10.1080/02533839.2013.814991#abstract

Ma M., McNatt T., Hays B., Hunter S., (03 January 2013), Elastic lateral distortional buckling
analysis of cantilever I-beams. Recuperado el 24 julio 2015 de:
http://bvirtual.uce.edu.ec:2053/doi/full/10.1080/17445302.2012.747282#abstract

9 ANEXOS

ESTUDIO DE LA FLEXION DE UNA VIGA EN VOLADIZO:

Anexo 1: Perfil de una barra que se flexiona debido a la acción de una carga aplicada en su
extremo libre

Consideremos una barra delgada de longitud L en posición horizontal, empotrada por un extremo y
sometida a una fuera vertical F en el extremo libre. Determinaremos la forma de la barra y las
coordenadas (xf, yf) del extremo libre para grandes flexiones de la barra.
Supongamos que:

 La barra tiene una longitud L mucho mayor que las dimensiones de su sección trasversal, y
que la deformación debida a su propio peso es despreciable.

 Que la sección de la barra no cambia cuando se dobla. Cuando el espesor de la barra es

pequeño comparado con el radio de curvatura, la sección trasversal cambia muy poco.

En estas condiciones es aplicable la ecuación de Euler-Bernoulli que relaciona el momento flector

M de la fuerza aplicada y el radio de curvatura ρ de la barra deformada

YI
M=
ρ

Donde Y es el módulo de Young del material e I es el momento de inercia de la sección trasversal

respecto del eje neutro. El radio de curvatura:

ds
ρ=
dφ

dφ M
=
ds Y . I

Anexo 2: Definición de parámetros presentes cuando una barra se flexiona.

El momento flector M de la fuerza F aplicada en el extremo libre de la barra respecto del punto P
(x, y) es M=F(xf-x)

dφ F
= (xf −x)
ds Y . I

Derivando con respecto a s, y teniendo en cuanta que cosφ=dx/ds:

d2 φ F
+ cosφ=0
d s2 Y . I
Para determinar φ(s) se resuelve la ecuación diferencial con las siguientes condiciones iníciales

φ ( L )=φ 0

dφ
|s=L=0
ds

Para obtener una solución de la ecuación diferencial, multiplicamos por dφ/ds la ecuación
diferencial

dφ d 2 φ F
+ senφ=0
ds d s 2 Y . I
2

(( )
d 1 dφ
ds 2 ds
2
+
F
Y .I
senφ =0 )
1 dφ F
( )
2 ds
+
Y .I
senφ=cte

La constante de integración la determinamos a partir de las condiciones iníciales especificadas

anteriormente:

dφ 2 2 F
( )
ds
=
Y .I
(sen φo−senφ)

ds=
√ Y .I
2F
dφ
√ sen φ o−senφ

La Longitud L de la barra y las coordenadas x e y de cada uno de los puntos de la misma se

obtienen:
φ0

L=∫ ds
0
φ0

L=
√
Y .I
∫
dφ
2 F 0 √ sen φ0−senφ
dx=ds . cosφ
φ
x=
√
Y .I
∫ cosφdφ
2 F 0 √ sen φ0−senφ
dy=ds . senφ
=
2Y . L
F √
( √ sen φ0− √ sen φ0 −senφ )

φ
y=
√
Y.I
∫ s enφdφ
2 F 0 √ sen φ0 −senφ
Dada la fuerza F aplicada en el extremo libre de la barra y conocida la longitud L de la barra, se
resuelve la primera ecuación para calcular el ángulo φ0, que forma la recta tangente a la barra en su
extremo libre con la parte negativa del eje horizontal X

Una vez que se conoce este ángulo φ0, se calcula la abscisa x dando valores al ángulo φ en el
intervalo (0, φ0)

El cálculo de la ordenada y es más complicado, ya que para cada valor del ángulo φ hay que hallar
una integral definida en el intervalo (0, φ) empleando procedimientos numéricos.

Cálculo numérico

Las ecuaciones anteriores las podemos expresar

φ0
dφ F L2
2 √ ∝=∫ ; ∝=
0 √ sen φ0−senφ
2Y .I
x 1
L √∝ √
= ( sen φ0−√ sen φ0−senφ )
φ
y 1 senφdφ
= ∫
L 2 √∝ 0 √ sen φ0−senφ

Donde α es un parámetro a dimensional que engloba las características geométricas de la barra, del
material del que está hecha, y de la fuerza aplicada en su extremo libre

Cálculo de φ0.

Empezamos con la primera ecuación que nos determina el ángulo φ0 que forma la recta tangente a la
barra en su extremo libre con la parte negativa del eje horizontal X, tal como se ve en la figura:

Anexo 3: Ubicación de del ángulo que forma la tangente a la barra en su extremo libre

Requiere dos pasos:

 1. Hallar la integral
φ0
dφ
∫
0 √ sen φ 0−senφ
2. Calcular la raíz de la ecuación

φ
f (¿¿ 0)=0
¿

La integral se puede expresar en términos de la suma de dos integrales elípticas de primera

especie, haciendo cambios de variable. El primer cambio es:

π
θ=φ+
2
φ0
dφ
∫
0 √ sen φ 0−senφ
π π
φ0 + φ 0+
2 2
dθ 1 dθ
¿ ∫ = ∫
π
2
√cosθ−cos ( φ0 −π /2 ) √2 π
2
√ sen ( φ0 /2+ π /4 )−sen2 ( θ/2 )
2

El segundo cambio de variable es

φ π
sen ∅=
sen ( θ/2 )
k (
; k=sen 0 +
2 4 )
2 kcos ∅d ∅
dθ=
√1−k 2 sen2 ∅
π
φ0 2
dφ d∅
∫ = √ 2∫ =¿
0 √ sen φ 0−senφ 0∅ √ 1−k sen ∅
2 2

¿√2 ∫
d∅
(
∅ √ 1−k sen ∅
2 2
2
−
0
∫ d∅
0 √ 1−k sen ∅
2 2
; sen
∅0

∅ 0 =
sen ( π /4 )
k )
Finalmente, calculamos la raíz de la ecuación

k , π /2
Y (¿)−Y (k ,∅ 0)
¿
√2 ¿
Cálculo de las coordenadas (x/L, y/L) de cada punto de la barra deformada
El cálculo de x/L no reviste dificultad alguna. Conocido φ0, se calcula x/L para cada ángulo φ en el
intervalo (0, φ0). La posición xf del extremo libre es

xf 1
=
L √∝
√ sen φ0

El cálculo de y/L es más problemático. Conocido φ0, se determina la ordenada y/L para cada ángulo
φ en el intervalo (0, φ0) calculando la integral definida,
φ
y 1 senφdφ
= ∫
L 2 √∝ 0 √ sen φ0−senφ

por el procedimiento numérico de Simpson

Cuando φ→φ0 el denominador de la integral tiende a cero. El ordenador no calcula correctamente la

ordenada yf/L del extremo libre de la barra cuando φ=φ0. Para solucionar este inconveniente,
empleamos el procedimiento de interpolación que se muestra en la figura.

Anexo 4: Intervalos que forman las coordenadas de las deformaciones de la barra

 Calculamos las coordenadas (x/L, y/L) para el ángulo φ=φ0-Δφ, siendo Δφ un ángulo
pequeño.

 Calculamos la abscisa xf/L para el ángulo φ0.

La ordenada yf/L se obtiene resolviendo el triángulo rectángulo de la figura

y f y xf x
(
= + − tan φ0
L L L L )
Aproximación de pequeñas flexiones

Para pequeñas flexiones cuando el ángulo φ0 es pequeño. Sustituimos senφ≈φ y escribimos la

ecuación que calcula φ0.
 φ0
dφ
∫ =2 √ α
0 √ φ0−φ
El resultado es φ0 =α

Las coordenadas (x, y) de cada punto de la barra se aproximan a

x
L
=1− 1−
φ
α √
Para el extremos libre de la barra, cuando φ= φ0=α, xf=L, lo que implica que en la aproximación de
pequeñas flexiones, no hay desplazamiento horizontal del extremo libre de la barra.

La ordenada y la podemos aproximar

φ
y 1 φdφ
= ∫
L 2 √ α 0 √ α −φ

Integrando por partes y después de hacer algunas simplificaciones obtenemos la siguiente expresión

y 2 φ
= α− +α
L 3 2 ( ( )√ ) 1−
φ
α

Las coordenadas x e y, las hemos expresado en función del parámetro φ, eliminando el parámetro
obtenemos la función y=f(x) que describe la flexión de la barra cuando se aplica una fuerza F en su
extremo libre.

x2 1 x 3
y
L (
=α 2 −
L 3 L3 )
3
y=
FL
2Y ∗I (
x 2−
1 x
3 L3 )
Para el extremos libre de la barra, cuando

φ=φ0 =α ; x=L
yf 2 1 L
3
= α y f= F
L 3 3 Y .I

Límite de la aproximación de pequeñas flexiones

En la figura, se muestra la desviación y/L del extremo libre de la barra en función del parámetro a
dimensional α.
  En color rojo, los resultados del cálculo, empleando procedimientos numéricos, descrito en
el apartado anterior

 En color negro, la recta y/L=2α/3, aproximación de pequeñas flexiones

Anexo 5: Grafica de la recta del límite de aproximación de pequeñas flexiones

Podemos considerar, que la aproximación lineal produce resultados aceptables hasta un cierto valor
límite del parámetro αm o bien, hasta un cierto valor máximo de la fuerza aplicada Fm en el extremo
libre de la barra

2Y . I . α m
Fm =
L2

Ejemplo:

Sea una regla de acero de longitud L=30 cm, sección rectangular a=3.04 cm, y b=0.078 cm. El
módulo de Young es Y=2.06·1011 N/m2

El momento de inercia I vale

a b3
I= =1.20 ×10−12 [ m4 ]
12

Cuando aplicamos en el extremo libre de la barra una fuerza tal que α=0.25, es decir

F L2
α= F=1.38 [ N ]
2Y . I

Aplicando la aproximación de pequeñas flexiones

yf 2
= α
L 3
2
1 0.3
yf = × ×1.38=0.05 [ m ] =5.0 [ cm ]
3 2.06 ×1011∗1.2× 10−12
En la aproximación de pequeñas flexiones xf≈L, no hay desviación apreciable en sentido horizontal
y la desviación en sentido vertical yf es proporcional a la fuerza F aplicada en el extremo libre.

Cuando aplicamos en el extremo libre de la barra una fuerza tal que α=1.25, es decir

F L2 2 α .2Y . I
α= ; F= 2
; F =6.88 [ N ]
2Y . I L

Aplicando la aproximación de pequeñas flexiones

1 3
∗L
3
yf = ∗F
Y .I
1
∗0,33
3
yf = ∗6.88
32.06 × 1011∗1.2 ×10−12
y f =0.25 [ m ]=25 [ cm ]

En la aproximación de pequeñas flexiones deja de ser válida ya que hay una desviación apreciable
en sentido horizontal y la desviación en sentido vertical yf ya no es proporcional al a la fuerza F
aplicada en el extremo libre.