Determinacin de Mdulos de Young - F.Arrufat, U. Novick, M.P. Frigerio y G.

Sardelli
Determinacin de Mdulos de Young

Arrufat, Francisco Toms [email protected]
Novick, Uriel Sebastin Tel: 4861-1542
Frigerio, Mara Paz [email protected]
Sardelli, Gastn [email protected]

Universidad Favaloro, Facultad de Ingeniera Bs. As. , Argentina- Julio 2001.

Utilizamos mtodos dinmicos y estticos para analizar los modos de vibracin
excitados en varas y tubos de diversos materiales. En el mtodo dinmico
estudiamos la dependencia de la frecuencia de vibracin con las longitudes de las
varas y tubos. En el mtodo esttico estudiamos la dependencia entre la
deformacin de las varas con la carga a la que son sometidas. Con la informacin,
buscamos obtener el mdulo de elasticidad de los materiales, encontramos que se
puede conseguir dicho valor con slo analizar el sonido emitido por los tubos al
golpearlos.


INTRODUCCIN

Todo material existente es elstico y se deforma en cierto grado. Al construir un edificio se debe buscar que el material
sea fuerte, pero a su vez flexible para no derrumbarse ante un terremoto; queremos que las alas de un avin puedan
flexionarse, pero que a su vez sean resistentes.
Para determinar las propiedades elsticas de un material dado, es necesario en general someter a ensayos una muestra
del material. Es claro que disponer de tcnicas de ensayos no destructivos es una gran ventaja en muchas aplicaciones
prcticas.
Definimos el mdulo de elasticidad E
1
de un material como la relacin entre el esfuerzo
A
F
y la deformacin
l
l
:
E

El objetivo de este trabajo es determinar el mdulo de elasticidad por distintos mtodos experimentales. Analizamos
tambin los diferentes modos de vibracin que surgen en los materiales al estar sometidos a las diferentes condiciones
de contorno.

I. EXPERIMENTO

El objetivo del presente experimento es, a travs del anlisis de varas y tubos de distintos materiales, variando las
condiciones de contorno en que se encuentran las mismas, determinar el mdulo de elasticidad de los materiales
causando el menor dao y deformacin a los mismos.
Para tal fin, optamos por dos mtodos de anlisis, uno dinmico y uno esttico.


MTODO DINMICO : Este mtodo es representado mediante una ecuacin diferencial de cuarto orden en la cual
varan los resultados segn las condiciones de contorno. Dicha ecuacin
2
es:

0
4
4
2
2

x
y
EI
t
y
A

(1)
donde A es el rea transversal a la onda de las varas, es la densidad del material, I es el momento de rea respecto de
una lnea neutra y E es el mdulo de elasticidad buscado.
Las experiencias que realizamos fueron las siguientes:

EXPERIMENTO CON AMBOS EXTREMOS LIBRES

En esta primera parte del experimento analizamos el sonido emitido por los distintos tubos al ser golpeados. Para ello
utilizamos tubos de bronce y de aluminio de distintos largos, un micrfono conectado a un digitalizador de seales
conectado a una PC.
Figura 1. Ilustracin del experimento realizado
con ambos extremos de la vara libres y sujetada
de un nodo tomando el sonido que emite al
golpearla con un micrfono.



El primer paso fue estimar la posicin del primer nodo desde
uno de los extremos de cada vara, considerando que los nodos
principales se encontraban ubicados aproximadamente al
22,4% de la longitud de las mismas. Esta medida se obtiene a
partir del anlisis de datos especificados en el apndice. (Es
fundamental que este paso sea efectuado con la mayor
precisin posible para optimizar al mximo la nitidez y
duracin del sonido, como as tambin cuidar de golpear el
tubo en un antinodo). Una vez hecho esto, procedimos a
golpear los distintos tubos cerca del micrfono, y registramos
los datos en la PC.
Una vez obtenidos los datos, y con la ayuda de la herramienta
FFT (Fast Fourier Transformer), efectuamos Transformada de
Fourier a los datos, obteniendo as los distintos modos que se
excitaban en los tubos y las correspondientes frecuencias.
Resolviendo la ecuacin (1) por el mtodo de separacin de
variables con las condiciones de contorno correspondientes
(*)

se llega a la siguiente expresin para cada modo:

A
EI
L
kL
f
n
n
2
2
2
) (


(2)

donde kL satisface Cos(kL)Cosh(kL)=1. Esta ecuacin tiene
varias soluciones, la cual da lugar a los distintos modos.
(kL)
1
=4,73004; (kL)
2
=7,8532; (kL)
3
=10,99561; etc.

EXPERIMENTO CON VARA EMPOTRADA

En esta parte del experimento utilizamos slo una vara de bronce. Aqu necesitamos un fotointerruptor conectado a una
PC y una prensa para empotrar la vara, como as tambin el programa para el procesamiento de los datos.
La disposicin de los elementos ha de ser la siguiente:

Figura 2. Ilustracin del experimento realizado con la vara empotrada tomando con un fotointerruptor las vibraciones
de la misma al moverla de su posicin de equilibrio.

Variando la longitud L del extremo libre de la vara, e imprimindole una ligera fuerza inicial para ponerla a oscilar, se
recogen los datos del periodo de dicha oscilacin para las distintas longitudes de L.
Resolviendo la ecuacin (1) con las condiciones de contorno correspondientes
(*)
se llega a la siguiente ecuacin para el
modo fundamental:

EI
L
f
2
2
2
597 . 0


(3)



*
La resolucin de la ecuacin (1) con las condiciones de contorno para llegar a (2) y (3) pueden verse en el apndice.
Determinacin de Mdulos de Young - F.Arrufat, U. Novick, M.P. Frigerio y G. Sardelli
MTODO ESTTICO : Este mtodo se representa mediante otra ecuacin diferencial, la cual describe el
desplazamiento en el eje y de la vara cuando se le suspenden diversos pesos. La ecuacin correspondiente
3
es:

2
2
) (
dx
y d
EI x M

(4)
siendo M(x) el momento flexor, I el momento de rea respecto de una lnea neutra y E el mdulo de elasticidad buscado.


EXPERIMENTO CON VARA EMPOTRADA Y DIVERSOS PESOS

Para esta parte del experimento, utilizamos la misma vara que en la etapa anterior. Esta vez, precisamos de la prensa
antes utilizada, una regla milimetrada y diversos pesos (P
o
) bien especificados (Tuvimos cuidado de no utilizar pesos
excesivos, ya que le imprimen una deformacin permanente a la vara, lo cual es no deseable).
Dejando libre una porcin de longitud L constante, procedimos a colgar los pesos a diferentes distancias l desde donde
empotramos la vara, midiendo el desplazamiento en el extremo libre de la vara con los distintos pesos.

Figura 3. Ilustracin del experimento realizado con
la vara empotrada agregando diversos pesos a la
misma.

La ecuacin (4) puede ser resuelta aplicando la transformacin
de Laplace
(*)
para llegar al siguiente resultado:

,
`

.
|

6 2
) (
2
x l
EI
x P
x y
o
0 < x < l

,
`

.
|

6 2
) (
2
l x
EI
l P
x y
o
l < x < L



(5)




II. RESULTADOS

Para cada experimento en particular analizamos el comportamiento de las varas y los tubos, y la dependencia que
pudieran las frecuencias con la distancia donde empotramos la vara y los pesos aplicados con la deformacin
consiguiente.
Conseguimos mediante grficos el mdulo de elasticidad, realizando una regresin lineal con los datos experimentales,
guindonos con las ecuaciones averiguadas de manera que la pendiente de la recta nos proporcionara el valor
correspondiente con su respectivo error.


MTODO DINMICO

ANLISIS DE FRECUENCIAS Y MODOS DE VIBRACIN

En la primer parte del experimento, determinamos las frecuencias correspondientes al modo fundamental para tubos de
bronce libres de distintos largos. El primer grfico realizado (Figura 4) se hizo con el fin de comprobar la
proporcionalidad entre la frecuencia fundamental con 1/L
2
, como seala la ecuacin (2).


*
Para ver el procedimiento de resolucin de la ecuacin (5) se puede consultar: Pipes, Louis A., Matemticas Aplicadas
para Ingenieros y Fsicos, Cap. 9, Seccin 4 (Ed. McGraw Hill, 1963).
y = 68,268x - 2,3643
R
2
= 0,9998
0
200
400
600
800
1000
1200
0 5 10 15 20
1/L [1/m]
f

[
h
z
]

Figura 4. Grfico de las frecuencias fundamentales en
funcin de 1/L
2
para los tubos de bronce.


Hicimos la misma comprobacin para los tubos de aluminio como puede observarse en la Figura 5.

y = 91,343x + 6,7212
R
2
= 0,9999
0
500
1000
1500
2000
2500
0 5 10 15 20 25 30
1/L [1/m]
f

[
h
z
]

Figura 5. Grfico de las frecuencias fundamentales en
funcin de 1/L
2
para los tubos de aluminio.

Para la vara de bronce utilizada en la segunda experiencia del mtodo dinmico, recogimos los datos del perodo de
oscilacin y realizamos el mismo grfico de la frecuencia en funcin de 1/L
2
realizado para la primera experiencia.

y = 3,3246x + 0,7351
R
2
= 0,9909
4
4,5
5
5,5
6
6,5
7
7,5
8
1,3 1,5 1,7 1,9 2,1
1/L [1/m]
f

[
h
z
]

Figura 6. Grfico de la frecuencia fundamental de la vara en
funcin de 1/L
2
.





Determinacin de Mdulos de Young - F.Arrufat, U. Novick, M.P. Frigerio y G. Sardelli
DETERMINACIN DEL MDULO DE ELASTICIDAD

Con el fin de conseguir el valor del mdulo de elasticidad, podemos rescribir la ecuacin (2) como:

CE f
n

2
con
4 2
4
4
) (
L A
kL I
C
n



(6)
La pendiente ser igual a E.

y = 9,6E+10x + 6,0E+04
R
2
= 1,0E+00
0,E+00
1,E+06
2,E+06
3,E+06
4,E+06
5,E+06
6,E+06
7,E+06
8,E+06
9,E+06
1,E+07
0 0,00002 0,00004 0,00006 0,00008 0,0001
C[m/kg]
f

[
h
z

]

y = 6,3E+10x + 4,2E+04
R
2
= 1,0E+00
0,E+00
1,E+06
2,E+06
3,E+06
4,E+06
5,E+06
6,E+06
7,E+06
8,E+06
9,E+06
1,E+07
0 0,00005 0,0001 0,00015
C[m/kg]
f

[
h
z

]

Figura 7. Grfico para la determinacin del mdulo de Figura 8. Grfico para la determinacin del mdulo de
elasticidad de los tubos de bronce. elasticidad de los tubos de aluminio.

Una ventaja de este ltimo procedimiento es que comprueba la validez general de la ecuacin (2) y nos permiti utilizar
todas las frecuencias correspondientes a los diversos modos excitados para obtener E.
A partir de los grficos, obtuvimos los valores del mdulo de elasticidad para ambos materiales:

E
bronce
= (9,6 0,2)10
10
Pa
E
aluminio
= (6,32 0,03)10
10
Pa

En el segundo experimento del mtodo dinmico, para determinar el mdulo de elasticidad de la vara, rescribimos la
ecuacin (3) como:

DE f
2
con
4
2
4
4
597 . 0
L
I
D

(7)

Analizamos el grfico (Figura 9) que muestra a E como la pendiente de la recta.

y = 8,1259E+10x + 5,2074E+00
R
2
= 9,9011E-01
20
25
30
35
40
45
50
55
60
2,5E-10 3,5E-10 4,5E-10 5,5E-10 6,5E-10
D[m/kg]
f

[
h
z

]

Figura 9. Grfico para la determinacin del mdulo de elasticidad
de la vara mediante regresin lineal.

El valor de E obtenido fue: E = (8,1 0,6)10
10
Pa


MTODO ESTTICO

ANLISIS DE LA DEFORMACIN EN FUNCIN DE LOS PESOS

Como sealan las ecuaciones (5), la deformacin en un punto es directamente proporcional al peso aplicado
transversalmente, manteniendo constante la posicin donde ste se aplica. Para corroborar esta relacin realizamos un
grfico ilustrativo.

y = 0,0077x - 0,0003
R
2
= 0,9828
0
0,002
0,004
0,006
0,008
0,01
0,012
0,014
0,016
0,018
0,02
0 0,5 1 1,5 2 2,5 3
Po [N]
y

[
m
]

Figura 10. Grfico de la deformacin en funcin del peso.

DETERMINACIN DEL MDULO DE ELASTICIDAD

Para determinar el mdulo de elasticidad de la misma vara utilizada en el experimento anterior, graficamos

,
`

.
|

6 2
2
l x
I
l P
G
o o
en funcin de la deformacin medida cerca del extremo libre en un punto x
o
en la Figura 11.
y = 9E+10x + 5E+07
R
2
= 0,9821
0,0E+00
5,0E+08
1,0E+09
1,5E+09
2,0E+09
2,5E+09
3,0E+09
3,5E+09
4,0E+09
0 0,01 0,02 0,03 0,04
y[m]
G
[
N
/
m
]

Figura 11. Grfico para determinar el mdulo de elasticidad de
la vara que puede apreciarse como la pendiente de la recta.

El mdulo de elasticidad obtenido en este caso fue: E = (9,3 0,5)10
10
Pa


III. CONCLUSIN

En este experimento se estudi, por una parte, los modos de vibracin de barras con distintas condiciones de contorno, y
tambin buscamos la determinacin del mdulo de elasticidad de un material causndole el menor dao posible.
En cuanto a los modos de vibracin, observamos por medio de las figuras 4, 5 y 6, que en todos los mtodos dinmicos,
la frecuencia de vibracin de los modos excitados depende de forma inversa con la longitud elevada al cuadrado de las
varas y los tubos. A partir de esto, observamos que cunto ms larga es la vara o el tubo, mayor es la cantidad de modos
que se excitan en los mismos.
Para los mdulos de elasticidad, construimos la siguiente tabla con los valores de E obtenidos para los distintos
materiales, y los valores que proporcionan los libros de texto.

Determinacin de Mdulos de Young - F.Arrufat, U. Novick, M.P. Frigerio y G. Sardelli
MATERIAL CASO 1 CASO 2 CASO 3 MEJOR VALOR TABLAS
BRONCE (9,6 0,2)10
10
Pa (8,1 0,6)10
10
Pa (9,3 0,5)10
10
Pa (9,4 0,2)10
10
Pa 10,8.10
10
Pa
ALUMINIO (6,32 0,03)10
10
Pa - - (6,32 0,03)10
10
Pa 7,0.10
10
Pa
Tabla 1. Cuadro comparativo de los valores de E obtenidos experimentalmente para los diversos materiales y los
valores conocidos.

De aqu observamos que los valores de los mdulos de elasticidad obtenidos, varan, para el bronce, en un rango de 8 a
10x10
10
Pa, lo cual concuerda con otras experiencias consultadas y es prxima al valor de las tablas, y en cuanto a los
tubos de aluminio, el valor obtenido fue prximo al valor que se encuentra en tablas. Estas variaciones son aceptables
ya que se deben a que, a pesar de tratarse de un mismo material, las varas y los tubos analizados poseen propiedades
propias de elaboracin.
Haciendo una comparacin de los resultados obtenidos por medio de los sonidos emitidos por los tubos versus los
obtenidos por medio de los otros mtodos, vemos que podemos obtener un valor del mdulo de elasticidad que se
encuentra en el rango general de aceptacin de estos valores para cada material. Como se ve, al poder obtener el valor
del mdulo de elasticidad con slo golpear el material y analizar los sonidos emitidos, poseemos una muy buena forma
de obtener este valor, sin causar prcticamente dao alguno al material. Esto sera un ejemplo de un ensayo no
destructivo.

BIBLIOGRAFA
1. Sears, Semansky y Young, Fsica Universitaria (Ed. Fondo Educativo Interamericano, 1986), Cap.10, Pg. 255.
2. D.L.R. Oliver, Hollow-Tube Chimes, Phys.Teach. 36, 209 (Abril 1998).
3. Pipes, Louis A., Matemticas Aplicadas para Ingenieros y Fsicos (Ed. McGraw Hill, 1963), Cap.9, Seccin 4.
4. Lafita Babio, Felipe y Mata Corts, Hilario, Introduccin a la teora de vibraciones mecnicas (Ed. Labor,
Barcelona).
5. www.fisicarecreativa.com Martnez, Patricia y Azuaga, Marcelo, Medicin del mdulo de elasticidad de Young.

APNDICE

En esta seccin pueden encontrarse los desarrollos tericos realizados para los distintos experimentos, los cuales son
tiles para analizar las ondas en cada situacin y para determinar el mdulo de elasticidad de los materiales:


EXPERIMENTO CON AMBOS EXTREMOS LIBRES

La ecuacin diferencial que describe las vibraciones de un tubo libre en ambos extremos
2
es la siguiente:

0
4
4
2
2

x
y
EI
t
y
A

(8)
donde A es el rea transversal a la onda de las varas, es la densidad del material, I es el momento de rea respecto de
una lnea neutra y E es el mdulo de elasticidad del material.
La ecuacin (8) puede resolverse por el mtodo de separacin de variables, considerando y(x,t)=X(x)T(t). La solucin de
la ecuacin es:

t i
e x X t x y
t
) ( ) , (
(9)
Como la vara est libre en los extremos, no se doblan ya que no hay fuerzas restauradoras, por lo tanto, las derivadas
segunda y tercera de X son nulas.

0
3
3
2
2

dx
X d
dx
X d


(10)
Aplicando las condiciones (10) a X(x) se obtiene:

[ ]

'

'

+
]
]
]

+ + ) ( ) (
) ( ) (
) ( ) (
) ( ) ( ) ( kx Senh kx Sen
kL Senh kL Sen
kL Cos kL Cosh
kx Cosh kx Cos A x X

(11)
donde kL satisface para cada modo:

1 ) ( ) ( kL Cosh kL Cos
(12)
DE la ecuacin (12) se deduce tambin los puntos de los nodos para los distintos modos, estando al 22,4% del extremo
para el primer modo. Por lo tanto, la ecuacin para la frecuencia angular es:

A
EI
k
2


(13)
De (13) sale que:

A
EI
L
kL
f
n
n 2
2
2
) (


(14)
A partir de (12) se est en condiciones de conocer el mdulo de elasticidad simplemente reemplazando los valores de
frecuencias obtenidos para cada modo con el kL correspondiente. Sabiendo que

A
dA y I
2
se da:

4
2 2
int ext
R R
A
I +


(15)
para los tubos con radio interno R
int
y radio externo R
ext
.


EXPERIMENTO CON VARA EMPOTRADA

La ecuacin diferencial que corresponde al movimiento de vibracin libre de la vara
4
es:

0
2
2
4
4
2

t
x
y
x
c

(16)
Siendo

EI
c
2
una constante del sistema, donde E es el mdulo de elasticidad, I el momento de rea respecto de
una lnea neutra y es la masa de la vara por unidad de longitud.

Este caso, al igual que el anterior, puede resolverse con el mtodo de separacin de variables tomando x(y,t)=Y(y)T(t) y
ajustando las condiciones de contorno para este caso:

0 ) , ( ) , (
0 ) , 0 ( ) , 0 (
3
3
2
2

t L
y
x
t L
y
x
t
y
x
t x



(17)
La solucin ser de la forma:

) ( ) ( ) , ( + t Sen y Y t y x
(18)
A partir de (16) y con las condiciones dadas en (17) se obtiene que:

'

]
]
]
]
]
]
]

+
+ +
c
Cosh y
c
Cos
L
c
Senh L
c
Sen
L
c
Cosh L
c
Cos
y
c
Senh y
c
Sen A y Y




) (



(19)
donde satisface:

L
c
Cosh
L
c
Cos

1


(20)
Lo cual da una infinidad de frecuencias naturales del sistema:

c
L
n
n 2
2
2
2
1

,
`

.
|


(21)
A partir de (21) se puede encontrar E reemplazando el valor de la frecuencia con su correspondiente n. En este caso, I
puede calcularse como:

4
4
R I



(22)
donde R es el radio de la vara maciza.