Rev. Int.

M
et. Num. C
alc. Dis. Ing.
Vol. 25, 2, 111-132 (2009)

Revista Internacional de
M
etodos Num
ericos para
C
alculo y Dise
no en Ingeniera

An
alisis de vibraci
on libre de una viga Timoshenko
escalonada, centrfugamente rigidizada, mediante el
m
etodo de cuadratura diferencial
Daniel H. Felix, Ra
ul E. Rossi
Instituto de Mec
anica Aplicada
Departamento de Ingeniera, Universidad Nacional del Sur
Alem 1253, 8000 Baha Blanca, Argentina

Diana V. Bambill
Instituto de Mec
anica Aplicada
Departamento de Ingeniera, Universidad Nacional del Sur
Miembro de la Carrera del Investigador Cientfico, CONICET
Alem 1253, 8000 Baha Blanca, Argentina
Tel.: 54-291-459 5100/01; Fax.: 54-291-459 5657
e-mail: [email protected]

Resumen
Este artculo trata sobre las vibraciones transversales de una viga rotante de secci
on no uniforme, utilizando
el metodo de cuadratura diferencial (MCD).
La formulaci
on se basa en la teora de Timoshenko, que considera las fuerzas de inercia debidas a la traslaci
on
lateral junto con las contribuciones de la deformaci
on por corte y de la inercia rotatoria. En este trabajo el
metodo de cuadratura diferencial se aplica al c
alculo de las frecuencias naturales. Los resultados obtenidos
son comparados con resultados publicados en la literatura abierta y con valores calculados por los autores
con un c
odigo de elementos finitos. El metodo propuesto ofrece un procedimiento directo y eficiente para
analizar con muy buena precisi
on las vibraciones naturales de vigas rotantes de secci
on no uniforme.

Palabras clave: Vibraciones libres, viga rotante, metodo de cuadratura diferencial, secci
on
no uniforme.

FREE VIBRATION ANALYSIS OF A CENTRIFUGALLY STIFFENED TIMOSHENKO

NON-UNIFORM BEAM USING THE DIFFERENTIAL QUADRATURE METHOD

Summary
This paper deals with the transverse vibration of non-uniform rotating beams, by means of the differential
quadrature method (DQM).
The formulation is based on the Timoshenko beam theory, which takes into account the inertia forces due
to the transverse translation together with the contributions of shear deformation and rotatory inertia. The
application of the differential quadrature method is demonstrated by investigating the natural frequencies.
The results are compared with published results given in the open literature and with values obtained by the
authors using a finite element code. The proposed method offers a direct and efficient procedure to analyze
the natural vibration of non-uniform rotating beams with very good accuracy.

Keywords: Free vibrations, rotating beam, differential quadrature method, non-uniform

beam.

c
Universitat
Polit`
ecnica de Catalunya (Espa
na).

ISSN: 02131315

Recibido: Sept. 2008 Aceptado: Enero 2009

112

D.H. Felix, D.V. Bambill y R.E. Rossi

INTRODUCCION
En el campo de la ingeniera, hay una gran variedad de estructuras que pueden modelarse como vigas cantilever rotantes con vibracion transversal, tales como helices, aspas de
helicopteros, alabes de turbinas, brazos de robot, entre otros.
Existen diversas teoras y una gran variedad de tecnicas analticas y numericas para
estudiar el problema de la vibracion transversal de vigas rotantes. Este estudio presenta
el metodo de cuadratura diferencial como procedimiento de calculo de las frecuencias naturales de vigas cantilever que rotan alrededor de un eje, y se analiza la influencia de la
variacion de la seccion de la viga en forma escalonada, como asimismo de los dos parametros
mas importantes como son la velocidad de rotacion y la esbeltez de la viga rotante. La
formulacion del problema propuesto se basa en la teora de vigas de Timoshenko1,2 , que da
mejores resultados en vigas de poca esbeltez y especialmente para los modos de vibracion
correspondientes a las frecuencias superiores, en comparacion con la teora elemental de
Bernoulli-Euler que solo tiene en cuenta la deformacion por flexion.
Se considera que el tema es tanto de interes practico como academico. Muchos investigadores se han dedicado a su estudio, y en la literatura cientfica se encuentran numerosas
publicaciones sobre el tema. Entre los metodos utilizados para resolverlo el mas difundido
y utilizado es quizas el metodo de elementos finitos37 . Ademas pueden citarse entre otros:
el metodo de matriz de rigidez dinamica1,8 , el metodo de Ritz9,10 , el metodo de Galerkin11
y el de serie de potencias12 .
En a
nos recientes el metodo de cuadratura diferencial1319 , reaparece como una herramienta u
til para resolver ecuaciones diferenciales relacionadas con aplicaciones a problemas
de ingeniera. Este metodo permite modelar numericamente el problema del ttulo y obtener
valores de las frecuencias con un alto grado de precision y un gasto trivial en tiempo de
computacion, tanto en la introduccion de los datos como en la obtencion de los resultados que demanda unos pocos segundos20 . El algoritmo del MCD se desarrollo en lenguaje
Mathematica21 . Si bien en este trabajo se resuelven, a modo de ejemplo, vigas rotantes
escalonadas con uno, dos y tres tramos, el planteo es aplicable a una viga con cualquier
n
umero de tramos.
A fin de brindar una confirmacion de los buenos resultados logrados se los compara con
valores presentes en la literatura (Banerjee1 ), y con los obtenidos utilizando el metodo de
elementos finitos (MEF), cuya aplicacion se explica mas adelante. Cabe mencionar que al
utilizar un software profesional de uso general para problemas de Mecanica de Solidos, la
carga de los datos, muy numerosos aun para el caso de modelos con pocos elementos finitos,
resulta tediosa, con el consiguiente riesgo de cometer errores difcilmente detectables. Por
otra parte, al utilizar rutinas de calculo de uso general, los tiempos de resolucion son
considerablemente mayores que con el algoritmo desarrollado con MCD.
DEL MODELO DE VIGA ROTANTE
DESCRIPCION
En la Figura 1 se presenta un esquema de la viga cantilever rotante de seccion escalonada, dividida en tramos de seccion constante, que rota con velocidad angular constante
El subndice k designa al tramo
con eje de rotacion Y , perpendicular al eje de la viga X.
generico kesimo de la viga, cuya longitud es lk . La posicion de la seccion de cada tramo
0k . En particular la comas cercana al eje de rotacion queda definida por la coordenada X

ordenada X01 corresponde al radio del eje material de rotacion. En este estudio dicho radio
01 = 0).
se ha considerado nulo (X

Los ejes X e Y son los ejes del plano de vibracion transversal de la viga cantilever, y el
plano de rotacion de la viga resulta normal al eje Y .

An
alisis de vibraci
on libre de una viga Timoshenko escalonada, centrfugamente rigidizada

 
 

 







Figura 1. Modelo esquem

atico de una viga cantilever rotante compuesta por 3
tramos de secci
on constante


 

 





 






Figura 2. Modelo esquem

atico de la viga cantilever rotante, en el que se pone
en evidencia la fuerza centrfuga Fk+1 , que act
ua sobre el elemento
k-esimo: a) para k = 2, b) para k = 1

113

114

D.H. Felix, D.V. Bambill y R.E. Rossi

Al producirse la rotacion se genera una fuerza centrfuga, indicada con Fk , producida

por la masa del tramo k-esimo y de los sucesivos tramos hasta el extremo libre de la viga
rotante. En las Figuras 2a y 2b se muestran las fuerzas F3 y F2 en el caso particular de la
viga escalonada de tres tramos de la Figura 1.
En las Figuras 3a, 3b y 3c se muestran los distintos casos de viga elegidos para el analisis
numerico. La Figura 3a corresponde a la viga de un tramo y las Figuras 3b y 3c a las vigas
escalonadas de dos y tres tramos, respectivamente. Puede suponerse que los casos mostrados
en las Figuras 3b y 3c se generan al remover parte del material de la viga de un solo tramo,
buscando un efecto de alivianamiento y al mismo tiempo de mayor rigidez dinamica.

 

 












Figura 3. Modelos resueltos de viga cantilever rotante: a) de un tramo, b) de
secci
on escalonada de dos tramos, c) de secci
on escalonada de 3 tramos.

An
alisis de vibraci
on libre de una viga Timoshenko escalonada, centrfugamente rigidizada

115

ECUACIONES DIFERENCIALES GOBERNANTES

Utilizando la teora de Timoshenko para vigas, las dos ecuaciones de movimiento de un
elemento k-esimo de viga que rota alrededor de un eje con velocidad , se expresan:
k (
Q
x, t)
2w
k (
x, t)
= A k
(1a)
2
x

t
k (
w
k (
x, t) M
x, t)
2 k (
x, t)
k (
k (
Q
x, t) N
x)
(1b)
+
+ Ik 2 k (
x, t) = Ik
2
x

t
donde Ak es el area e Ik es el momento de inercia de la seccion transversal constante; es
la densidad del material.
Con w
k (
x, t) se indica al desplazamiento transversal y con k (
x, t) al giro de una seccion
de la viga ubicada en la posicion x
(0 x
lk ) del k-esimo tramo. En tanto que los
k (
esfuerzos internos en la seccion analizada son el esfuerzo de corte Q
x, t) y momento

flector Mk (
x, t), que expresados en funcion de w
k (
x, t) y k (
x, t), resultan:

w
k (
x, t)
w
k (
x, t)
k (
k (
Q
x, t) = N
x)
+ GAk
k (
x, t)
(2a)
x

k (
x, t)
k (
(2b)
M
x, t) = E Ik
x

donde G es el modulo de elasticidad transversal, E el de elasticidad longitudinal y el

factor de corte utilizado para definir el area de corte efectiva de la seccion.
k (
La fuerza centrfuga N
x) se expresa como:

1 2
k (
0 k lk + 1 l2 X
0k x
N
x) = 2 Ak X

+ Fk+1 , 0 x
lk
(3)
2 k
2
donde Fk+1 es la fuerza centrfuga que act
ua sobre el k-esimo tramo, en la seccion del
mismo mas alejada del eje de rotacion. Obviamente, para el u
ltimo tramo (k = u) es
Fk+1 = Fu+1 = 0.
Al reemplazar las ecuaciones constitutivas (2a, b) en las ecuaciones de movimiento (1a,
b) se tiene:
k (
dN
x) w
k (
x, t)
2 w
k (
x, t)
k (
+N
x)
2
d
x
x

w
k (
x, t) k (
x, t)
2w
k (
x, t)
+A k G

=
A
k
2
2
x

2 k (
x, t)
2 k (
x, t)
w
k (
x, t)
2
k (
x, t) + E Ik
+
I

(
x
,
t)
=
I
Ak G
k
k
k
2
2
x

(4a)

(4b)

Un desarrollo detallado que muestra la obtencion de las ecuaciones diferenciales gobernantes aplicando el principio de Hamilton, fue presentado por Banerjee 1 , quien tambien
puso de relieve la importancia de incluir en la ecuacion (1b) el termino en el que interviene el
cuadrado de la velocidad angular, para obtener resultados mas cercanos al comportamiento
dinamico real de la viga, especialmente en el caso de altas velocidades de rotacion.
x, t) se pueden escribir como:
Asumiendo una oscilacion armonica simple, w(
x, t) y (
k (
w(
x, t) = W
x).eit ;

x, t) =
k (
(
x).eit

116

D.H. Felix, D.V. Bambill y R.E. Rossi

donde es la frecuencia de vibracion natural del sistema para el modo normal correspondiente. Al ser sustituidas en las ecuaciones (4a, b), resultan:
2

k (x)
k (
d2 W
x)
d Wk (
x)
0k + x
k (
2 Ak X
d Wd

N
x
)

A
G

k
x
d
x2
d
x2
A k G

k (
d
x)
d
x

k (
= A k 2 W
x)
(5a)

k (
dWk (
x)
d2
x)
k (
k (
k (
x) E Ik
Ik 2
x) = Ik 2
x)
d
x
d
x2

(5b)

Introduciendo un formato adimensional, se pueden escribir las siguientes ecuaciones para

cada tramo k-esimo de viga:

d2 Wk

s21 d2 Wk
dk
dWk
s21

= 2 Wk
(X0k + x)
2 Nk
dx
dx2
2(1 + ) lk2
dx2
dx
lk

s21 d2 k
2 2 d Wk
2 2
2

+
s s
2
s s k = 2 k
2(1 + ) 1 k dx
2(1 + ) 1 k
lk dx2
2

(6a)
(6b)

donde es el coeficiente de Poisson, y ademas se han definido las siguientes magnitudes

adimensionales:
lk
: Longitud del tramo.
lk =
L
x

x = : Coordenada espacial sobre el eje de la viga, local de cada tramo; 0 x 1.

lk
0k
X
X0k = : Coordenada global del extremo mas cercano al eje Y , del tramo k-esimo.
lk
k (
W
x)
Wk =
: Amplitud de desplazamiento de la seccion.
lk
k (
k =
x) : Amplitud de giro de la seccion.
k (
N
x)
Nk =
: Esfuerzo normal o axil en la seccion.
EAk
L
lk
sk = p
: Factor de esbeltez del tramo k-esimo.
lk Ik /A k
r
A1 2
=
L : Coeficiente de velocidad de rotacion de la viga.
rE I1
A 1 2
=
L : Coeficiente de frecuencia natural de vibracion transversal.
E I1
En la viga compuesta por mas de un tramo (viga escalonada), deben considerarse las
ecuaciones de compatibilidad de desplazamientos y esfuerzos entre tramos adyacentes. En
el caso de los desplazamientos deben satisfacerse:

k
(k1)
+W
=0
(7a)
W
x
=0
x
=l (k1)

(k1)
k

+
=0
(7b)
x
=0
x
=l
(k1)

que en forma adimensional se escriben:

l(k1) W(k1) x=1 + lk Wk |x=0 = 0

(k1) x=1 + k |x=0 = 0

(8a)
(8b)

An
alisis de vibraci
on libre de una viga Timoshenko escalonada, centrfugamente rigidizada

117

De modo similar, se impone la compatibilidad de esfuerzos internos en las secciones que

limitan tramos contiguos:

(k1)
k
Q
+Q
=0
(9a)
x
=0
x
=l
(k1)

(k1)
M
x
=l

(k1)

k
+M
=0
x
=0

(9b)

cuya forma adimensional resulta:

(k1) Q(k1) x=1 + k Qk |x=0 = 0

(k1)
k
M(k1) x=1 +
Mk |x=0 = 0
l(k1)
lk

(10a)
(10b)

donde k = Ak /A 1 es la relacion entre el area de la seccion transversal del tramo kesimo y

k (
x)
el area de la seccion transversal del tramo 1 mas cercano al eje Y ; Qk = QEA
es el esfuerzo
k

lk
de corte y Mk =
Mk (
x) el momento flector, expresados en forma adimensional. Si se
EIk

indican los esfuerzos en funcion de los desplazamientos, las ecuaciones (10a, b) resultan:

dW(k1) (k1)

(k1)
N(k1) +

2(1
dx
+ )
2(1 + ) x=1

(11a)

dWk
k

=
0
+k
Nk +
2(1 + )
dx
2(1 + ) x=0

(k1) d(k1)
d

k
k

+
=0
(11b)

l(k1)
dx
lk dx x=0
x=1
Las condiciones de borde del primer tramo (k = 1, en x
= 0, borde empotrado):

1
W
=0
(12a)
x
=0

=0
(12b)
x
=0
y en forma adimensional:
W1 |x=0 = 0

(13a)

1 |x=0 = 0

(13b)

En tanto las condiciones de borde para el u

ltimo tramo (k = u, en x
= lu , borde libre),
son:

u = 0
Q
(14a)
x
=lu

u = 0
M
(14b)
x
=lu
en forma adimensional se tiene:
Qu |x=1 = 0

(15a)

Mu |x=1 = 0

(15b)

118

D.H. Felix, D.V. Bambill y R.E. Rossi

0.4

0.2

-0.2

0.2

0.4

0.6

0.8

Figura 4. Distribuci
on irregular de nodos sobre una viga escalonada de 3 tramos

y expresando las ecuaciones (15) en funcion de los corrimientos resultan:

dWu

=0
Nu +

u
2(1 + )
dx
2(1 + )
x=1

du
=0
dx x=1

(16a)
(16b)

EL METODO
DE CUADRATURA DIFERENCIAL
El metodo de cuadratura diferencial fue inicialmente propuesto en 1971 por Bellman13 ,
como una tecnica numerica destinada a resolver ecuaciones diferenciales.
Una excelente revision del tema, que abarca el perodo comprendido entre sus comienzos
y mediados de los 90, fue publicado en 1996 por Bert y Malik14 .
En las u
ltimas dos decadas la popularidad del MCD fue creciendo y actualmente estan
disponibles, en la extensa literatura existente, diversas variantes de la metodologa original.
En este trabajo se aplica a la resolucion del problema de vibraciones enunciado. Cada tramo
continuo de viga fue reemplazado por un mallado de puntos o nodos a lo largo de su eje.
Los puntos adoptados se distribuyen en forma irregular utilizando la expresion formulada
por ChebyshevGaussLobato18 :
1 cos [(i 1) / (n 1)]
xi =
, i = 1, 2, ..., n
(17)
2
donde n es el n
umero de nodos que se adopta en el tramo y xi la coordenada espacial
adimensional del nodo i. La Figura 4 muestra el mallado de puntos generado al utilizar la
expresion (17) en cada tramo de una viga escalonada de tres tramos.
Seg
un lo indicado por Bert y Malik14 , en cada punto i del mallado se aproximan las
derivadas de las funciones Wk y k mediante la aplicacion de las denominadas ecuaciones
an
alogas de cuadratura:

n
X
d(q) Wk
(q)
Ai j Wk j
(18a)
=
q
dx
j=1
xi

n
X
d(q) k
(q)
=
Ai j k j
(18b)

dx q
xi

j=1

An
alisis de vibraci
on libre de una viga Timoshenko escalonada, centrfugamente rigidizada

119

donde q indica el orden de derivacion; n es el n

umero de puntos del mallado en cada tramo;
i, j son subndices que se refieren a los puntos de la malla; y el subndice k identifica al
tramo de la viga. Con Wk j y k j se indican el desplazamiento y el giro, adimensionales, en
el punto j del tramo kesimo.
(q)
Los factores Ai j son los coeficientes de ponderacion relacionados con los valores de las
funciones en los n puntos del mallado. Para su calculo se siguieron las referencias 14,15,18 ,
como se describe a continuacion. Como primer paso, asociado a cada punto i del mallado se
calcula un polinomio de Lagrange cuyas races son las coordenadas de los restantes puntos
del tramo, seg
un la expresion:
Y

(x i ) =

n
Y

(x i x j )

(19)

j=1 j6=i
(1)

Seguidamente, en base a los polinomios de Lagrange, se determinan los coeficientes Ai j ,

correspondientes a las derivadas de primer orden, q = 1, seg
un las siguientes expresiones:
Q
(x i )
(1)
Q
A ij =
con i, j = 1, 2, ..., n para j 6= i
(20a)
(xi xj )
(xj )
(1)
Aii

n
X

(1)

A ij con i, j = 1, 2, ..., n para j = i

(20b)

j=1 j6=i

Finalmente, los coeficientes correspondientes a las derivadas de segundo orden y superiores, q > 1, se obtienen mediante expresiones recursivas de la forma:
(q1)

(q)

(q1)

Aij = Ai i

(1)

Aij
n
X

(q)

Aii =

Aij

xi xj
(q)

Aij

con i, j = 1, 2, ..., n para j 6= i

con

i, j = 1, 2, ..., n para j = i

(21a)
(21b)

j=1
con j 6= i

Ecuaciones an
alogas de cuadratura del sistema
Las ecuaciones diferenciales de movimiento (6a, b), junto con las de compatibilidad (8a,
b) y (11a, b), y las condiciones de borde (13a, b) y (16a, b), quedan reemplazadas por
un conjunto de ecuaciones lineales algebraicas. Tales expresiones se obtienen mediante la
aplicacion de las reglas de cuadratura diferencial dadas por las ecuaciones (18a, b).
De esta forma, las ecuaciones analogas de cuadratura correspondientes al sistema de
ecuaciones diferenciales (6a, b) resultan:
2 (X0k
+

2(1 +

+ x i)

n
P

(1)
A ij

j=1
n
2
X
s1
(1)
A ij k j
) lk2
j=1

s2 s2
2(1 + ) 1 k

n
X
j=1

(1)

A ij Wk j

Wk j

n
n

s21
s21 X (2)
(2)
A ij Wk j
A ij Wk j
2 Nk i
2(1 + ) lk2
lk
j=1
j=1

= 2 Wk i
n
s21 X (2)
A
lk2 j=1 ij

k j +

s2 s2 2
2(1 + ) 1 k

X
n

(22a)
(0)

A ij k j = 2 k i

j=1

(22b)

120

D.H. Felix, D.V. Bambill y R.E. Rossi

donde la fuerza centrfuga Nk i , de acuerdo a la ecuacion (3), queda expresada:

2
1
1 2
2 lk
Nk i = Nk |x=xi = 2 X0 k + X0k x i x i + Fk+1
2
2
s1
Las ecuaciones analogas de cuadratura correspondientes a las ecuaciones de compatibilidad de desplazamientos (8a, b) son:
l(k1)

n
X

(0)

A nj W(k1)j + lk

j=1

n
X

n
X

(0)

(23a)

A 1 j k j = 0

(23b)

A 1j Wk j = 0

j=1

(0)
A nj

(k1)j +

j=1

n
X

(0)

j=1

mientras que las correspondientes a las de compatibilidad de esfuerzos internos (11a, b)

son:

P
n
n
P
(0)
(1)

(k1)
A nj (k1) j
N(k1) n (k1) + 2(1+)
A nj W(k1) j 2(1+)
j=1

! j=1
(24a)

P
n
n
P
(1)
(0)

+k
Nk 1 + 2(1+)
A 1j Wk j 2(1+)
A 1j k j = 0
j=1

j=1

n
n
(k1) X
k X (1)
(1)

A 1j k j = 0
A nj (k1) j +
l(k1)
lk

(24b)

j=1

j=1

Las ecuaciones analogas de cuadratura de las condiciones de borde para el extremo

empotrado, ecuaciones (13a, b) resultan:
n
X

(0)

(25a)

(0)

(25b)

A 1j W1 j = 0

j=1
n
X

A 1j 1j = 0

j=1

y para el extremo libre, ecuaciones (16a, b) se tiene:

X
n
n
X

(1)
(0)
A n j Wu j
An j uj = 0
Nu |x=1 +
2(1 + )
2(1 + )
j=1

n
X

(26a)

j=1

(1)

A n j uj = 0

(26b)

j=1

El conjunto de ecuaciones analogas de cuadratura diferencial, ecuaciones (22) a (26),

conforma un sistema lineal, cuyos autovalores son los coeficientes de frecuencia de la viga
rotante a ser determinados.

An
alisis de vibraci
on libre de una viga Timoshenko escalonada, centrfugamente rigidizada

121

RESULTADOS NUMERICOS
Los resultados numericos que se presentan corresponden a las primeras frecuencias naturales de vibracion de una viga Timoshenko, empotrada en uno de sus extremos a un eje
que gira con velocidad angular constante.
Para elegir la cantidad conveniente de puntos (n) a adoptar en el mallado de cada tramo
de seccion constante se realizo un analisis de convergencia. Esto permitio determinar que es
suficiente utilizar un mallado con 15 puntos en el tramo, ya que adoptar mallas mas densas
conduce a valores que coinciden con cinco dgitos significativos, para todas las frecuencias
y velocidades de rotacion adoptadas en este estudio.
Por otra parte, con el fin de de corroborar el correcto dise
no y desarrollo del algoritmo
basado en el MCD, algunos casos tambien fueron resueltos mediante el metodo de elementos
finitos. Para ello se recurrio al software profesional ALGOR22 utilizando el elemento-viga
disponible de dos nodos, con seis grados de libertad en cada nodo, con el cual se contempla la
deformacion axial y su efecto rigidizante. Para la construccion de los modelos se utilizo una
misma malla rectilnea con mil elementos-viga de igual longitud, la cual se comprobo que
resultaba suficientemente precisa, mediante un analisis previo de convergencia. Los valores
obtenidos para el caso de la viga escalonada en dos y tres tramos son mostrados en las
Tablas VI y VII.
La Tabla I muestra los coeficientes de frecuencia fundamental, calculados con el MCD,
de la viga rotante de un tramo en toda su longitud L, en funcion del parametro de velocidad de rotacion . En el calculo se utilizaron las relaciones E/ G = 3,059 y s = 30, que
fueron tomadas de la referencia1 para realizarpla primera comprobacion de los resultados.
Los valores de frecuencia fundamental 1 = A/E I L2 1 coinciden en todas las cifras
significativas con los calculados por el metodo de la matriz de rigidez dinamica de la citada
referencia.

0
1
2
3
4
5

MCD
3.4798
3.6445
4.0971
4.7516
5.5314
6.3858

Ref. 1
3.4798
3.6445
4.0971
4.7516
5.5314
6.3858

p
Tabla I. Coeficientes de frecuencia fundamental 1 =
A/EI L2 1 de una
viga cantilever rotante de un tramo, en funci
on de , con E/G = 3,059;
= 0,30; s = 30

En las Tablas II a V se muestran los resultados obtenidos con el MCD de los coeficientes
i (con i = 1,2,3,4 ) para la viga de un tramo, en funcion del parametro de Timoshenko
r, definido como la inversa del factor de esbeltez s. La Tabla II muestra los coeficientes
de la frecuencia fundamental en funcion de la velocidad de rotacion , del parametro de
Timoshenko r, y del factor de corte . La incorporacion del valor r = 0,15, aunque sin
mucha aplicacion practica, permite mostrar la estabilidad numerica del metodo. Para cada
par de valores y r se presentan tres coeficientes 1 , el primero corresponde a una viga
de seccion rectangular ( = 0.8457), el segundo a una de seccion circular ( = 0.8864) y
el tercero, en letra italica, corresponde al de una viga rotante en que no se considera en el
calculo el efecto de la deformacion por corte.
Las Tablas III, IV y V, analogas a la Tabla II, muestran los coeficientes de las frecuencias
superiores 2 , 3 y 4 , respectivamente, para el mismo caso de viga de un tramo.

122

D.H. Felix, D.V. Bambill y R.E. Rossi

0
1
2
3
4
5
6
7
8
9
10

0.01
3,5127
3.5128
3.5152
3.6782
3.6784
3.6807
4.1336
4.1337
4.1362
4.7931
4.7932
4.7958
5.5800
5.5802
5.5832
6.4436
6.4438
6.4473
7.3533
7.3534
7.3578
8.2911
8.2914
8.2968
9.2468
9.2471
9.2536
10.214
10.214
10.222
11.189
11.189
11.198

0.02
3,5028
3.5032
3.5128
3.6681
3.6685
3.6780
4.1226
4.1231
4.1327
4.7806
4.7810
4.7913
5.5653
5.5658
5.5776
6.4260
6.4266
6.4407
7.3323
7.3330
7.3502
8.2662
8.2671
8.2881
9.2174
9.2185
9.2440
10.180
10.181
10.212
11.149
11.151
11.187

0.03
3,4866
3.4875
3.5087
3.6515
3.6524
3.6734
4.1046
4.1055
4.1268
4.7601
4.7610
4.7838
5.5413
5.5424
5.5684
6.3976
6.3989
6.4297
7.2985
7.3000
7.3374
8.2264
8.2283
8.2737
9.1710
9.1732
9.2280
10.126
10.129
10.194
11.088
11.091
11.168

r = 1/s
0.05
3.4364
3.4388
3.4958
3.6000
3.6024
3.6590
4.0492
4.0516
4.1083
4.6976
4.7001
4.7601
5.4688
5.4716
5.5389
6.3126
6.3157
6.3946
7.1992
7.2028
7.2968
8.1115
8.1158
8.2277
9.0399
9.0447
9.1768
9.9786
9.9840
10.138
10.924
10.930
11.107

0.08
3.3240
3.3295
3.4648
3.4853
3.4908
3.6242
3.9270
3.9323
4.0639
4.5620
4.5674
4.7028
5.3150
5.3206
5.4678
6.1371
6.1431
6.3095
7.0004
7.0069
7.1981
7.8893
7.8963
8.1159
8.7951
8.8025
9.0521
9.7125
9.7202
10.001
10.639
10.647
10.957

0.10
3.2309
3.2388
3.4368
3.3908
3.3985
3.5928
3.8277
3.8350
4.0237
4.4541
4.4613
4.6508
5.1955
5.2027
5.4030
6.0045
6.0119
6.2318
6.8544
6.8620
7.1077
7.7304
7.7381
8.0130
8.6238
8.6316
8.9372
9.5297
9.5374
9.8738
10.445
10.452
10.819

0.15
2.9655
2.9791
3.3444
3.1242
3.1372
3.4889
3.5545
3.5661
3.8895
4.1665
4.1766
4.4762
4.8863
4.8952
5.1839
5.6687
5.6765
5.9672
6.4878
6.4947
6.7980
7.3276
7.3339
7.6589
8.1775
8.1835
8.5393
9.0291
9.0354
9.4327
9.8745
9.8821
10.335

p
Tabla II. Coeficientes de frecuencia fundamental 1 =
A/EI L2 1 de una
viga cantilever rotante de un tramo, en funci
on de y r, ( = 0,30)

An
alisis de vibraci
on libre de una viga Timoshenko escalonada, centrfugamente rigidizada

0
1
2
3
4
5
6
7
8
9
10

0.01
21.891
21.895
21.999
22.038
22.042
22.145
22.472
22.476
22.578
23.178
23.182
23.282
24.131
24.135
24.233
25.303
25.307
25.403
26.665
26.669
26.763
28.189
28.193
28.284
29.847
29.851
29.942
31.619
31.623
31.713
33.485
33.489
33.578

0.02
21.478
21.495
21.893
21.625
21.642
22.038
22.060
22.077
22.468
22.767
22.783
23.167
23.721
23.737
24.111
24.893
24.909
25.273
26.254
26.268
26.624
27.773
27.787
28.135
29.425
29.440
29.782
31.189
31.202
31.541
33.043
33.057
33.394

0.03
20.841
20.875
21.720
20.989
21.023
21.863
21.426
21.460
22.288
22.137
22.169
22.979
23.093
23.125
23.913
24.267
24.298
25.062
25.626
25.656
26.397
27.141
27.170
27.892
28.787
28.815
29.520
30.540
30.567
31.259
32.381
32.408
33.091

r = 1/s
0.05
19.138
19.213
21.191
19.290
19.364
21.329
19.739
19.811
21.739
20.464
20.533
22.405
21.436
21.503
23.306
22.623
22.686
24.414
23.989
24.049
25.703
25.505
25.561
27.145
27.143
27.196
28.716
28.880
28.930
30.394
30.699
30.746
32.162

0.08
16.288
16.406
20.045
16.452
16.568
20.172
16.932
17.044
20.548
17.702
17.808
21.161
18.722
18.820
21.988
19.950
20.040
23.007
21.344
21.426
24.191
22.870
22.944
25.515
24.497
24.565
26.958
26.202
26.264
28.499
27.966
28.023
30.120

0.10
14.530
14.661
19.136
14.704
14.833
19.254
15.212
15.335
19.604
16.018
16.132
20.172
17.072
17.176
20.940
18.321
18.414
21.886
19.717
19.801
22.984
21.221
21.296
24.213
22.797
22.866
25.551
24.420
24.483
26.978
26.066
26.127
28.479

123

0.15
11.036
11.165
16.748
11.235
11.360
16.841
11.799
11.915
17.118
12.651
12.754
17.568
13.692
13.784
18.177
14.829
14.914
18.924
15.979
16.064
19.792
17.075
17.167
20.761
18.061
18.170
21.812
18.899
19.033
22.929
19.568
19.735
24.099

p
Tabla III. Coeficientes de la segunda frecuencia natural 2 = A/EI L2 2 de
una viga cantilever rotante de un tramo, en funci
on de y r, ( = 0,30)

124

D.H. Felix, D.V. Bambill y R.E. Rossi

0
1
2
3
4
5
6
7
8
9
10

0.01
60.755
60.783
61.460
60.900
60.928
61.604
61.333
61.362
62.033
62.048
62.077
62.742
63.034
63.062
63.719
64.277
64.304
64.951
65.760
65.787
66.423
67.466
67.493
68.117
69.376
69.402
70.015
71.472
71.497
72.098
73.734
73.759
74.349

0.02
58.197
58.297
60.765
58.345
58.445
60.907
58.785
58.884
61.330
59.511
59.609
62.029
60.511
60.608
62.992
61.770
61.864
64.207
63.269
63.362
65.658
64.990
65.080
67.327
66.913
67.001
69.198
69.019
69.104
71.251
71.287
71.371
73.470

0.03
54.638
54.821
59.657
54.790
54.973
59.796
55.244
55.426
60.209
55.991
56.170
60.892
57.018
57.194
61.834
58.308
58.479
63.021
59.840
60.006
64.438
61.593
61.755
66.069
63.546
63.702
67.896
65.677
65.828
69.901
67.967
68.113
72.067

r = 1/s
0.05
46.749
47.059
56.486
46.917
47.226
56.616
47.419
47.723
57.002
48.241
48.538
57.639
49.363
49.651
58.518
50.761
51.038
59.626
52.407
52.672
60.948
54.274
54.525
62.469
56.332
56.570
64.171
58.557
58.782
66.038
60.924
61.136
68.053

0.08
36.704
37.071
50.528
36.910
37.273
50.640
37.516
37.870
50.976
38.496
38.835
51.530
39.812
40.132
52.292
41.415
41.714
53.253
43.259
43.537
54.399
45.298
45.554
55.716
47.490
47.726
57.188
49.799
50.016
58.801
52.193
52.394
60.538

0.10
31.667
32.025
46.493
31.899
32.252
46.595
32.578
32.918
46.897
33.661
33.983
47.395
35.088
35.086
48.082
36.788
37.063
48.947
38.695
38.947
49.978
40.747
40.979
51.162
42.889
43.106
52.486
45.074
45.281
53.934
47.261
47.466
55.494

0.15
22.993
23.301
37.807
23.268
23.571
37.885
24.043
24.334
38.119
25.183
25.464
38.504
26.514
26.800
39.036
27.865
28.178
39.707
29.125
29.483
40.506
30.274
30.675
41.426
31.351
31.775
42.454
32.384
32.805
43.580
33.341
33.734
44.794

p
Tabla IV. Coeficientes de la tercera frecuencia natural 3 = A/EI L2 3 de
una viga cantilever rotante de un tramo, en funci
on de y r, ( = 0,30)

An
alisis de vibraci
on libre de una viga Timoshenko escalonada, centrfugamente rigidizada

0
1
2
3
4
5
6
7
8
9
10

0.01
117.56
117.66
120.05
117.71
117.81
120.20
118.17
118.26
120.64
118.91
119.01
121.37
119.95
120.05
122.39
121.27
121.37
123.68
122.86
122.95
125.25
124.71
124.80
127.06
126.80
126.89
129.12
129.12
129.22
131.41
131.67
131.76
133.92

0.02
109.18
109.49
117.59
109.33
109.65
117.73
109.80
110.11
118.16
110.58
110.89
118.88
111.65
111.96
119.88
113.02
113.32
121.14
114.66
114.96
122.67
116.56
116.85
124.44
118.71
119.00
126.45
121.10
121.38
128.69
123.69
123.97
131.14

0.03
98.779
99.293
113.81
98.947
99.459
113.95
99.447
99.956
114.36
100.27
100.78
115.05
101.42
101.91
116.01
102.86
103.35
117.23
104.60
105.08
118.70
106.61
107.07
120.71
108.86
109.32
122.35
111.35
111.79
124.51
114.06
114.48
126.86

r = 1/s
0.05
79.233
79.944
103.82
79.433
80.142
103.95
80.030
80.731
104.33
81.013
81.701
104.95
82.362
83.032
105.82
84.053
84.703
106.92
86.058
86.685
108.24
88.349
88.951
109.79
90.894
91.470
111.54
93.665
94.214
113.48
96.633
97.156
115.60

0.08
58.272
58.989
87.710
58.528
59.240
87.815
59.283
59.982
88.130
60.505
61.183
88.650
62.144
62.796
89.373
64.138
64.764
90.292
66.422
67.022
91.400
68.927
69.506
92.688
71.585
72.150
94.148
74.327
74.891
95.769
77.083
77.662
97.541

0.10
48.222
48.902
78.213
48.486
49.162
78.306
49.253
49.919
78.586
50.450
51.104
79.048
51.966
52.615
79.690
53.664
54.322
80.507
55.394
56.081
81.492
57.016
57.755
82.638
58.432
59.236
83.937
59.608
60.474
85.381
60.562
61.477
86.959

125

0.15
29.327
29.935
60.200
29.272
29.977
60.271
29.877
30.369
60.484
30.477
30.945
60.836
31.229
31.655
61.325
32.116
32.476
61.949
33.103
33.384
62.701
34.113
34.325
63.578
35.064
35.234
64.574
35.935
36.097
65.683
36.798
36.980
66.897

p
Tabla V. Coeficientes de la cuarta frecuencia natural 4 = A/EI L2 4 de
una viga cantilever rotante de un tramo, en funci
on de y r, ( = 0,30)

126

D.H. Felix, D.V. Bambill y R.E. Rossi

0
1
2
3
4
5
6
7
8
9
10

Metodo
MCD
MEF
MCD
MCD
MCD
MCD
MCD
MEF
MCD
MCD
MCD
MCD
MCD
MEF

1
3.9554
3.9554
4.1199
4.5701
5.2171
5.9836
6.8208
6.8301
7.7012
8.6092
9.5355
10.475
11.423
11.444

%
0.000
0.136
0.184

2
12.854
12.854
13.085
13.752
14.795
16.136
17.697
17.718
19.416
21.246
23.152
25.110
27.105
27.160

%
0.000
0.119
0.203

3
37.850
37.850
38.024
38.539
39.383
40.536
41.973
42.005
43.666
45.587
47.709
50.005
52.452
52.565

%
0.000
0.076
0.215

4
63.293
63.291
63.499
64.113
65.123
66.507
68.240
68.284
70.292
72.632
75.229
78.051
81.071
81.230

%
-0.003
0.064
0.196

p
Tabla VI. Coeficientes de las primeras 4 frecuencias i = A/EI L2 i de una
viga escalonada de 2 tramos, calculados con MCD y MEF, en funci
on de
(Fig. 3b), con E/G = 3,059; = 0,30; 2 = 0,6; l1 = 0,6; l2 = 0,4;
s1 = 30; s2 = 83,333

0
1
2
3
4
5
6
7
8
9
10

Metodo
MCD
MEF
MCD
MCD
MCD
MCD
MCD
MEF
MCD
MCD
MCD
MCD
MCD
MEF

1
4.3874
4.3873
4.5389
4.9571
5.5655
6.2953
7.1012
7.1110
7.9561
8.8438
9.7541
10.680
11.618
11.641

%
-0.002
0.138
0.198

2
11.861
11.861
12.110
12.824
13.924
15.313
16.904
16.928
18.632
20.451
22.330
24.250
26.198
26.253

%
0.000
0.142
0.210

3
29.468
29.466
29.709
30.419
31.559
33.074
34.902
34.938
36.981
39.258
41.687
44.236
46.878
46.991

%
-0.007
0.103
0.241

4
51.791
51.789
51.994
52.600
53.601
54.981
56.721
56.764
58.797
61.179
63.836
66.737
69.849
70.004

%
-0.004
0.076
0.222

p
Tabla VII. Coeficientes de las primeras 4 frecuencias i = A/EI L2 i de una
viga escalonada de 3 tramos, calculados con MCD y MEF, para diferentes valores del par
ametro de velocidad de rotaci
on (Fig. 3c), con
E/G = 3,059; = 0,30; 2 = 0,6; 3 = 0,4; l1 = 0,6; l2 = 0,2; s1 = 30;
s2 = 83,333; s3 = 187,5

p
A continuacion, en las Tablas VI y VII se presentan valores i = A1 /E I1 L 2 i
calculados para vigas rotantes de seccion no uniforme, con E/ G = 3,059. Las vigas
. anap
lizadas son las que se muestran en los modelos de las figuras 3b y 3c. r = I1 /A1 L . Las
tablas contienen los coeficientes de las cuatro primeras frecuencias naturales calculados por
los metodos de cuadratura diferencial y de elementos finitos.
La diferencia entre ambos grupos de resultados es mnima y se indica en la columna %.
Para la viga en reposo ( = 0) los valores son casi totalmente coincidentes. Para valores de
>0 aparecen peque
nas diferencias que se incrementan con el aumento de . En particular,

An
alisis de vibraci
on libre de una viga Timoshenko escalonada, centrfugamente rigidizada

127

para un parametro alto de la velocidad de rotacion, = 10, las diferencias mayores son del
orden del 0.241 % para 3 . Los autores atribuyen esta diferencia a la menor precision del
algoritmo de rigidizacion que utiliza el metodo de elementos finitos. No obstante, se destaca
el muy buen acuerdo de los resultados obtenidos por el MCD con los hallados con el MEF.
En las Tablas VIII a XI se muestran los coeficientes de las cuatro primeras frecuencias
calculados con el MCD para vigas de seccion escalonada de dos y tres tramos, en funcion de
r y de la velocidad . Las relaciones de longitud para los tramos de viga son las graficadas
en las Figuras 3b y 3c. De cada coeficiente se presentan dos valores, el primero corresponde
a una viga de seccion rectangular con = 0.8457 en la que se considera el efecto de corte
y el segundo valor, a una viga en la que este efecto ha sido despreciado en el calculo.

N
umero de tramos
2

3
2

3
2

10

0.01
4.2469
4,2501
4.7232
4,7272
6.9223
6,9265
7.2236
7,2284
11.542
11,552
11.749
11,759

0.02
4.2348
4,2478
4.7086
4,7244
6.9041
6,9209
7.2030
7,2221
11.504
11,542
11.708
11,747

r=
0.03
4.2149
4,2439
4.6845
4,7197
6.8745
6,9115
7.1696
7,2117
11.445
11,524
11.644
11,728

I1 /A1 /L
0.05
0.08
4.1534 4.0151
4,2314 4,2015
4.6102 4.4442
4,7049 4,6693
6.7854 6.5981
6,8816 6,8094
7.0688 6.8570
7,1785 7,0986
11.280 10.985
11,469 11,336
11.466 11.144
11,666 11,518

0.10
3.9003
4,1744
4.3072
4,6371
6.4546
6,7439
6.6946
7,0261
10.789
11,214
10.928
11,382

0.15
3.5716
4,0843
3.9200
4,5302
6.0934
6,5233
6.2873
6,7831
10.329
10,803
10.427
10,925

p
Tabla VIII. Coeficiente de frecuencia fundamental 1 = A/EI L2 1 de una viga
cantilever rotante de secci
on escalonada, de 2 y 3 escalones, calculado
con MCD, en funci
on de y r, con = 0,30; 2 = 0,6; 3 = 0,4;
s2 = 1,667/r; s3 = 2,5/r

N
umero de tramos
2

3
2

3
2

10

0.01
10.084
18,140
18.741
18,798
21.574
21,629
21.857
21,912
29.551
29,613
29.186
29,246

0.02
17.874
18,094
18.543
18,764
21.357
21,573
21.654
21,871
29.297
29,534
28.958
29,189

r=
0.03
17.543
18,019
18.228
18,708
21.016
21,480
21.335
21,804
28.907
29,404
28.607
29,095

I1 /A1 /L
0.05
0.08
16.603 14.851
17,784 17,243
17.321 15.594
18,533 18,125
20.071 18.392
21,189 20,517
20.436 18.799
21,591 21,096
27.884 26.221
28,996 28,048
27.671 26.124
28,800 28,111

0.10
13.661
16,782
14.394
17,770
17.309
19,941
17.719
20,665
25.196
27,227
25.176
27,509

p
Tabla IX. Coeficiente de la segunda frecuencia natural 2 = A/EI L2 2 de
una viga cantilever rotante de secci
on escalonada, de 2 y 3 escalones,
calculado con MCD, en funci
on de y r, con = 0,30; 2 = 0,6;
3 = 0,4; s2 = 1,667/r; s3 = 2,5/r

0.15
11.056
15,414
11.712
16,680
15.049
18,213
15.438
19,334
22.644
24,715
22.975
25,640

128

D.H. Felix, D.V. Bambill y R.E. Rossi

N
umero de tramos
2

3
2

3
2

10

0.01
49.116
49,557
45.939
46,292
52.565
52,985
49.101
49,441
61.784
62,167
57.582
57,899

0.02
47.534
49,182
44.689
46,027
51.014
52,577
47.872
49,154
60.270
61,672
56.378
57,552

r=
0.03
45.237
48,575
42.836
45,595
48.778
51,917
46.062
48,685
58.121
60,871
54.634
56,987

I1 /A1 /L
0.05
0.08
39.777 32.198
46,781 43,155
38.258 31.566
44,287 41,512
43.541 36.461
49,964 46,017
41.661 35.410
47,266 44,252
53.230 46.727
58,499 53,697
50.538 44.925
55,276 51,632

0.10
28.201
40,494
27.907
39,363
32.790
43,119
32.064
41,915
43.048
50,168
41.735
48,797

0.15
21.266
34,176
21.369
33,879
26.158
36,241
25.831
35,944
30.733
41,794
30.492
41,526

p
Tabla X. Coeficiente de la tercera frecuencia natural 3 =
A/EI L2 3 de
una viga cantilever rotante de secci
on escalonada, de 2 y 3 escalones,
calculado con MCD, en funci
on de y r, con = 0,30; 2 = 0,6;
3 = 0,4; s2 = 1,667/r; s3 = 2,5/r

N
umero de tramos
2

3
2

3
2

10

0.01
95.072
96,442
83.320
84,262
98.452
99,791
86.750
87,673
107.92
109,18
96.266
97,153

0.02
90.332
95,250
79.959
83,384
93.772
98,554
83.415
86,760
103.36
107,82
92.969
96,140

r=
0.03
83.910
93,355
75.279
81,980
87.476
96,587
78.802
85,298
97.332
105,65
88.484
94,518

I1 /A1 /L
0.05
0.08
70.283 53.857
87,953 77,873
64.875 51.506
77,917 70,081
74.306 58.779
90,978 80,507
68.687 56.038
81,065 72,894
85.113 70.870
99,457 87,887
78.974 67.546
89,817 80,721

0.10
45.830
71,088
44.722
64,610
51.141
73,458
49.743
67,183
61.484
80,095
61.004
74,345

0.15
31.141
56,630
32.486
52,398
33.066
58,436
34.664
54,416
37.639
63,499
40.449
60,041

p
Tabla XI. Coeficiente de la cuarta frecuencia natural 4 =
A/EI L2 4 de
una viga cantilever rotante de secci
on escalonada, de 2 y 3 escalones,
calculado con MCD, en funci
on de y r, con = 0,30; 2 = 0,6;
3 = 0,4; s2 = 1,667/r; s3 = 2,5/r

En la Figura 5 estan graficadas las tendencias de variaci

on de las cuatro primeras frecuencias de una viga de un tramo, al aumentar la velocidad de rotacion y variar el valor
del parametro de Timoshenko. Para todas las frecuencias, las vigas de mayor esbeltez son
las que tienen los coeficientes de frecuencia mas altos. Por otra parte, se observa que la
rigidizacion de la viga a causa del aumento de la velocidad es mucho mas importante para
la frecuencia fundamental que para las frecuencias superiores.
En la Figura 6 se muestran como varan los coeficientes de frecuencia 1 , 2 , 3 y 4
en funcion de la velocidad de rotacion , para vigas con r = 1/30, de seccion rectangular de
uno, dos y tres tramos. Las relaciones geometricas adimensionales adoptadas en cada caso

An
alisis de vibraci
on libre de una viga Timoshenko escalonada, centrfugamente rigidizada

:1

129

:2

12

35

r 0.05
10

30
r 0.01

25

r 0.08
r 0.10
6

r 0.05

20

r 0.01

r 0.10

r 0.08

15

K
0

10

(a)

10

10

(b)

:3

:4
140

70
r 0.01

r 0.01

120

60
r 0.05

100

50

r 0.05
80

r 0.08

40

r 0.08
60
r 0.10

r 0.10

K
0

10

K
0

(c)

(d)
p

Figura 5. Coeficientes de las cuatro primeras frecuencias i = A/EI L2 i de una viga

cantilever rotante, de secci
on transversal uniforme, en funci
on de , para diferentes
valores del par
ametro de Timoshenko r, con E/G = 3,059; = 0,30

analizado son:
- Viga de un tramo: l1 = 1; 1 = 1.
- Viga de dos tramos: l1 = 0,60; l2 = 0,40; 1 = 1; 2 = 0,6.
- Viga de tres tramos: l1 = 0,60; l2 = 0,20; l3 = 0,20; 1 = 1; 2 = 0,6; 3 = 0,4.
Como comentario puede remarcarse que los coeficientes de frecuencia fundamental de
las vigas mas livianas, seg
un se observa en la Figura 6, en todos los casos son mas altos para
la misma velocidad de rotacion (rigidizacion). Sin embargo, para las frecuencias superiores
2 , 3 y 4 se presenta la relacion inversa (flexibilizacion).
CONCLUSIONES
El metodo es aplicable a la resolucion de ecuaciones diferenciales como las que gobiernan el problema propuesto. Con el se logran muy buenos resultados para los valores de
frecuencia, tal como puede apreciarse de las comparaciones presentadas. Aun para el caso
en que la viga vare su seccion en forma escalonada, el modelado de la misma es sencillo. La
formulacion mediante las ecuaciones de cuadratura es relativamente simple de desarrollar,
facil de programar y calcular con una computadora de caractersticas estandar. Los tiempos
de carga de datos son extremadamente breves y es muy corto tambien el tiempo que insume
el calculo de las frecuencias.

130

D.H. Felix, D.V. Bambill y R.E. Rossi

!1

!2

12
32.5

10

30
27.5

3 tramos

1 tramo

25

22.5

1 tramo
20

2 tramos

3 tramos

17.5
2 tramos

10

(a)

10

(b)

!3

!4

70
65

110

60

100

1 tramo

1 tramo
55
90
2 tramos

50

2 tramos
80

45

3 tramos

3 tramos
0

10

(c)

10

(d)
p

Figura 6. Coeficientes de las primeras cuatro frecuencias i = A/EI L2 i de una viga

cantilever rotante, de secci
on escalonada, en funci
on de , con = 0,30; 2 = 0,6;
3 = 0,4; s1 = 30; s2 = 50; s3 = 75

El objetivo del trabajo es realizar un aporte, mostrando las ventajas que brinda este procedimiento y sin pretender generar una discusion sobre el metodo de calculo mas adecuado
a utilizar por el ingeniero de dise
no. Los autores esperan que otros investigadores se sientan
estimulados a desarrollar el metodo de cuadratura diferencial para resolver estructuras mas
complejas.
AGRADECIMIENTOS
Los autores desean agradecer a los revisores del trabajo que con sus sugerencias y
atinadas observaciones han permitido enriquecerlo.
Este estudio fue desarrollado en el Departamento de Ingeniera de la Universidad Nacional del Sur, patrocinado por el Programa de Investigaci
on y Desarrollo del CONICET,
por la Agencia Nacional de Promocion Cientfica y Tecnol
ogica, y por la Secretara General
de Ciencia y Tecnologa de la citada universidad.
REFERENCIAS
1 J.R. Banerjee, Dynamic stiffness formulation and free vibration analysis of centrifugally stiffened Timoshenko beam, Journal of Sound and Vibration, Vol. 247, pp. 97-115, (2001).
2 M.H. Seon, H. Benaroya y T. Wei, Dynamics of transversely vibrating beams using four engineering theories, Journal of Sound and Vibration, Vol. 225, pp. 935-988, (1999).

An
alisis de vibraci
on libre de una viga Timoshenko escalonada, centrfugamente rigidizada

131

3 R.E. Rossi, R.H. Gutierrez y P.A.A. Laura, Transverse vibrations of a Timoshenko beam of
non uniform cross section elastically restrained at one end and carrying a concentrated mass at
the other, J. Acoust Soc. Am, Vol. 89, pp. 2456-2458, (1991).
4 J.Chung y H.H. Yoo, Dynamic analysis of a rotating cantilever beam by using the finite element
method, Journal of Sound and Vibration, Vol. 249, pp. 147-164, (2002).
5 K.G.Vinod, S. Gopalakrishnan y R. Gangul, Free vibration and wave propagation analysis of
uniform and tapered rotating beams using spectrally formulated finite elements, International
Journal of Solids and Structures, Vol. 44, pp. 5875-5893, (2007).
6 R.E. Rossi, Introducci
on al an
alisis de Vibraciones con el Metodo de Elementos Finitos, Editorial EdiUNS, Universidad Nacional del Sur, Baha Blanca, Argentina. IBSN 978-987-1171-71-2,
(2007).
7 S.S. Rao y R.S. Gupta, Finite element vibration analysis of rotating Timoshenko beams,
Journal of Sound and Vibration, Vol. 242, pp. 103-124, (2001).
8 J.R. Banerjee, Free vibration of centrifugally stiffened uniform and tapered beams using the
dynamic stiffness method, Journal of Sound and Vibration, Vol. 233, pp. 857-875, (2000).
9 S.K. Sinha, Non-Linear dynamic response of a rotating radial Timoshenko beam with periodic
pulse loading at the free end, International Journal of Non-Linear Mechanics, Vol. 40, pp.
113-149, (2005).
10 S.Y. Lee y Y.H. Kuo, Bending frequency of a rotating Timoshenko beam with general elastically
restrained root, Journal of Sound and Vibrations, Vol. 162, pp. 243-250, (1993).
11 M.D. Al-Ansary, Flexural vibrations of rotating beams considering rotary inertia, Computers
& Structures, Vol. 69, pp. 321-328, (1998).
12 S.C. Lin y K.M. Hsiao, Vibration analysis of a rotating Timoshenko beam, Journal of Sound
and Vibration, Vol. 240, pp. 303-322, (2001).
13 R. Bellman y J. Casti, Differential quadrature and long-term integration, J. Math. Anal. App,
Vol. 34, pp. 235-238, (1971).
14 C.W. Bert y M. Malik, Differential quadrature method in computational mechanics: A review,
Applied Mechanics Review, Vol. 49, pp. 1-28, (1996).
15 G. Karami, P. Malekzadeh y S.A. Shahpari, A DQEM for vibration of shear deformable nonuniform beams with general boundary conditions, Engineering Structures, Vol. 25, pp. 1169-1178,
(2003).
16 S.T. Choi, J.D. Wu y Y.T. Chou, Dynamic analysis of a spinning Timoshenko beam by the
differential quadrature method, American Institute of Aeronautics and Astronautics, Vol. 38,
pp. 851-856, (2000).
17 G.R. Liu y T.Y. Wu, Vibration analysis of beams using the generalized differential quadrature
rule and domain decomposition, Journal of Sound and Vibration, Vol. 246, pp. 461-481, (2001).
18 C. Shu y W. Chen, On optimal selection of interior points for applying discretized boundary
conditions in DQ vibration analysis of beams and plates, Journal of Sound and Vibration, Vol.
222, pp. 239-257, (1999).
19 S.J. Escanes, D.V. Bambill y D.H. Felix, El metodo de cuadratura diferencial generalizado:
estudio de arcos circulares delgados vibrantes, Rev. Internac. Met. Num. C
alc. Dis. Ing., Vol.
24, pp. 93-100, (2008).
20 D.V. Bambill, D.H. Felix y R.E. Rossi, The differential quadrature method in the free vibration
analysis of uniform rotating beams, Enviado a publicacion.

132

D.H. Felix, D.V. Bambill y R.E. Rossi

21 S. Wolfram, Mathematica: A System for Doing Mathematics by Computer, Third Edition,

Addison-Wesley, (1996).
22 ALGOR Professional Software V21.3, Linear Stress and Dynamics Analysis, Algor Inc., Pittsburgh, PA, USA, (2007).