Ogeâpkec"Eqorwvcekqpcn"Xqn"ZZXK.

"rr02429-2440
Ugtikq"C0"Gncumct."Gnxkq"C0"Rknqvvc."Igtoâp"C0"Vqttgu"*Gfu0+
Eôtfqdc."Ctigpvkpc."Qevwdtg"4229

FRECUENCIAS NATURALES DE VIGAS DE PARED DELGADA

DOBLEMENTE ASIMÉTRICAS CARGADAS AXIALMENTE
Fernanda M. de Borbóna, Daniel Ambrosinia
a
Área de Dinámica Experimental, Facultad de Ingeniería, Universidad Nacional de Cuyo. CONICET.
Centro Universitario - Parque Gral. San Martín - (5500) Mendoza. Fax 54 0261 4380120. e-mail:
[email protected] [email protected], web page:
http://fing.uncu.edu.ar/estructural/index.html

Palabras Clave: Vigas de pared delgada, Frecuencias naturales, Carga axial

Resumen. Las vigas de pared delgada y sección abierta son extensivamente usadas como componentes
estructurales en distintos tipos de estructuras en los campos de la Ingeniería Mecánica, Aeronáutica y
Civil. En el caso de vigas que tienen dos ejes de simetría en su sección transversal, la solución de
vibraciones libres es suficientemente conocida y estudiada. Inclusive cuando la sección transversal
tiene un eje de simetría, se ha estudiado ampliamente, en general usando la teoría de Bernoulli-Navier.
Sin embargo, es bastante limitado el número de estudios para vigas doblemente asimétricas en las
cuales existe un triple acoplamiento de vibraciones de flexión en dos direcciones y torsión.
En trabajos anteriores se ha desarrollado un modelo teórico-numérico que considera las vibraciones
flexotorsionales acopladas. En cuanto al modelo físico teórico-numérico, la formulación inicial basada
en las hipótesis de Vlasov, fue modificada permitiendo la inclusión de deformaciones angulares
inducidas por el esfuerzo de corte, sección variable en la dirección longitudinal e inercias rotacionales
en las ecuaciones de movimiento. Se utiliza un método de integración numérica conjuntamente con el
algoritmo de la transformada rápida de Fourier (FFT) para trabajar en el dominio de la frecuencia.
En este trabajo se incorpora el efecto de una carga axial a las ecuaciones flexotorsionales acopladas.
Se presenta además, una comparación numérica de las frecuencias naturales de vibración.

Eqr{tkijv"B"4229"Cuqekcekôp"Ctigpvkpc"fg"Ogeâpkec"Eqorwvcekqpcn"
jvvr<11yyy0coecqpnkpg0qti0ct
2429
1 INTRODUCCIÓN
La determinación de frecuencias naturales y modos de vibración de vigas continuas no
amortiguadas, se discute en detalle en Pestel and Leckie(1963), donde también se describe el
cálculo de la respuesta dinámica a excitación armónica. Por otra parte, Ebner and
Billington(1968) usan integración numérica para estudiar las vibraciones estacionarias de
vigas de Timoshenko amortiguadas. En la literatura pueden encontrarse también numerosas
aplicaciones de vibración de vigas rectas y curvas, como así también arcos y estructuras de
cáscara. Por otra parte, la teoría formulada por Vlasov(1963) ha sido ampliamente usada en el
análisis dinámico de vigas de pared delgada y sección abierta, como es ejemplificado en los
estudios de Christiano and Culver(1969) y Yoo and Fehrenbach(1982) en los cuales existe
gran correlación entre trabajos experimentales y las predicciones teóricas de frecuencias
naturales y desplazamientos inducidos en vigas rectas y curvas sometidas a cargas móviles.
En vigas esbeltas, la modelación de la estructura con un gran número de elementos finitos,
elementos de cáscara u otros métodos más refinados, no conduce a ninguna mejora en la
correlación con resultados experimentales, porque las diferencias entre los modelos teóricos
normalmente son más pequeñas que los errores experimentales.
Sin embargo, aunque la teoría de Vlasov para vigas de sección abierta ya está firmemente
establecida, presenta algunas limitaciones tales como: a) Suponer que las deformaciones por
corte, al igual que en la teoría de Bernoulli, no contribuyen a la flexibilidad de la viga. Por lo
tanto, pueden encontrarse errores importantes en el caso de vigas poco esbeltas o en los
modos altos de vigas esbeltas (Timoshenko and Young, 1968). b) También se desprecia la
influencia de las inercias rotacionales en las resultantes de tensiones. c) Las ecuaciones de
cuarto orden de Vlasov son válidas solamente para vigas de sección uniforme. En Ambrosini
et al. (1995) y Ambrosini et al. (2000) se propuso una teoría modificada, la cual está basada
en la formulación de Vlasov pero que tiene en cuenta los tres efectos mencionados
anteriormente. Esta formulación, que utiliza el método de variables de estado en el dominio
de la frecuencia, conduce a un tratamiento numérico eficiente, el cual puede ser muy útil en
numerosas aplicaciones.
En la literatura pueden encontrarse otras teorías que también tienen en cuenta el
acoplamiento entre flexión y torsión, tales como las presentadas por Gere and Lin (1958) que
obtienen una ecuación simplificada para vigas uniformes de sección abierta y Muller (1983)
que formula una teoría general que incluye todos los efectos de acoplamiento entre las
ecuaciones de movimiento, pero que es de dificultosa aplicación práctica. Otras
contribuciones en este campo están restringidas a aplicaciones particulares. Por ejemplo,
Aggarwal and Cranch (1967) y Yaman (1997) trabajan con secciones tipo canal y Ali Hasan
and Barr (1974) con secciones ángulo.
Más recientemente, Tanaka and Bercin (1999) extendieron el modelo de Bishop et al
(1989) para estudiar acoplamiento triple de vigas uniformes usando Mathematica. La
ecuaciones diferenciales presentadas por los autores mencionados tienen una confusión del
sistema de coordenadas que fue clarificada por Arpaci and Bozdag (2002). En Ambrosini and
Danesi (2003) se presentó un estudio numérico comparando los resultados con los obtenidos
con Arpaci and Bozdag (2002). Luego L.Jun et al (2004) extendieron el modelo de Arpaci
and Bozdag (2002) considerando el efecto de la carga axial. Obtuvieron la matriz de
transferencia dinámica para una viga con sección arbitraria pero las deformaciones por corte
no fueron incluidas en la formulación.
La determinación con precisión de las frecuencias naturales de vibración y de las formas
modales, es de suma importancia en el diseño de vigas de pared delgada axialmente cargadas
sujetas a cargas dinámicas. Aspas de helicópteros, álabes de turbinas, ejes de automóviles son

2430
algunos ejemplos de este tipo de vigas. La utilización de materiales compuestos también
induce al acoplamiento de las vibraciones flexotorsionales. L. Jun et al (2004) presentaron un
modelo analítico para determinar la respuesta dinámica de vigas de pared delgada de
materiales compuestos bajo carga axial considerando las deformaciones por corte e inercias
rotacionales. Viola et al (2007) estudiaron el cambio en las vibraciones libres de vigas de
pared delgadas luego de producidas fisuras utilizando el método de la matriz de rigidez
dinámica.
En este trabajo se presenta una viga de pared delgada axialmente cargada. Se incorpora
este efecto en el modelo presentado en Ambrosini et al. (2000). Se presenta una comparación
numérica de las frecuencias naturales de vibración.

2 TEORÍA

2.1 Ecuaciones de movimiento

Siguiendo la convención adoptada por Vlasov, se adopta un sistema rectangular de
coordenadas globales izquierdo (x, y, z). (Ver Figura 1). Los desplazamientos asociados son
designados ξ,η y ζ

Figura 1: Definición de términos

En la figura, O representa el centroide de la sección transversal y A el centro de corte.

Se considera una fuerza axial longitudinal aplicada en un punto arbitrario de la sección ex,
ey. Si P es una fuerza de compresión y Mx y My (Ver figura 2) son los momentos de flexión
producidos por la carga excéntrica, se obtiene para las tensiones normales σ

P My M
σ =− + x− x y (1)
FT Jy Jx

El desplazamiento de un punto M ubicado en la sección transversal queda determinado

por:

2431
 ξ M = ξ − ( y − a y )θ
(2)
η M = η + ( x − ax )θ

Donde ax y ay son las coordenadas del centro de corte y θ es el giro de la sección

transversal alrededor del polo.
Los esfuerzos normales σδds actuando en la sección transversal de una faja elemental, son
proyectados en las direcciones de los ejes coordenados x e y como resultado de las
deformaciones por flexión.

Figura 2: Proyecciones px y py

Llamando a las magnitudes de estas proyecciones px y py, y considerando que las tensiones
normales no varían a lo largo de la viga obtenemos de la figura

dz
p x dzds = σδds
ρx
(3)
dz
p y dzds = σδds
ρy

1 ∂ 2ξ M ∂γ mx
= φ x′ = −
ρx ∂z 2
∂z
(4)
1 ∂ 2η M ∂γ my
= φ y′ = −
ρy ∂z 2 ∂z

Donde ρx y ρy son los radios de curvatura y γmx y γmy son los valores medios de las
deformaciones debidas al esfuerzo de corte. Luego

2432
 ∂ 2 ξ ∂γ mx ∂ 2θ
p x = σδ ( 2 − ) − σδ ( y − a y ) 2
∂z ∂z ∂z
(5)
∂ η ∂γ my
2
∂ 2θ
p y = σδ ( 2 − ) + σδ ( x − a x ) 2
∂z ∂z ∂z
A partir de (5), se pueden determinar las cargas por unidad de longitud qx, qy y m

∂ 2 ξ ∂γ mx ∂ 2θ
qx = − σδds − σδ ( y − a y )ds
∂z 2 ∂z L ∂z 2 L
∂ 2η ∂γ my ∂ 2θ
qy = − σδds + σδ ( x − a x )ds
∂z 2 ∂z L ∂z 2 L
(6)
∂ 2 ξ ∂γ mx ∂ 2η ∂γ my
m=− − σδ ( y − a y )ds + − σδ ( x − a x )ds +
∂z 2 ∂z L ∂z 2 ∂z L

∂ 2θ
+ σδ [( x − a x ) 2 + ( y − a y ) 2 ]ds
∂z L
2

Sustituyendo (1) en las expresiones (6) y considerando que δds=dF;y que Sx=Sy=Jxy=0:

∂ 2ξ ∂γ mx ∂ 2θ
qx = −P − − P ( a y − e y )
∂z 2 ∂z ∂z 2
∂ 2η ∂γ my ∂ 2θ
q y = −P − + P ( a x − e x )
∂z 2 ∂z ∂z 2
(7)
∂ 2ξ ∂γ mx ∂ 2η ∂γ my
m = − P (a y − e y ) − + P ( a x − e x ) − +
∂z 2 ∂z ∂z 2 ∂z
∂ 2θ
+ P(−r 2 + 2β y e y − 2β x e x )
∂z 2

Donde r y βx y βy son características geométricas

Jx + Jy
r 2 = a x2 + a y2 + (8)
FT

Uy
βx = − ax
2J y
(9)
U
β y = x − ay
2J x

U x = y 3 dF + x 2 ydF
F F
(10)
U y = x 3 dF + y 2 xdF
F F

2433
 Los términos de la derecha de la expresión (10) corresponden a los momentos y productos
de inercia de tercer orden de la sección considerada.
Considerando el principio de D’Alembert se obtienen las ecuaciones de movimiento para
vibraciones flexotorsionales de vigas de pared delgada cargadas axialmente:

∂ 4ξ ∂ 3γ mx ∂ 2ξ ∂γ mx ∂ 2θ
EJ y − + P − + P ( a y − e y ) +
∂z 4 ∂z 3 ∂z 2 ∂z ∂z 2
(11a)
∂ ξ
2
∂θ2
∂ ξ ∂ γ mx
4 3
+ ρFT + ρFT a y 2 − ρJ y ( 2 2 − ) = qx
∂t 2
∂t ∂z ∂t ∂z∂t 2

∂ 4η ∂ γ my ∂ 2η ∂γ my
3
∂ 2θ
EJ x − +P − − P(a x − e x ) 2 +
∂z 4 ∂z 3 ∂z 2 ∂z ∂z
(11b)
∂η
2
∂θ2
∂η ∂ γ my
4 3

+ ρFT − ρFT a x 2 − ρJ x ( 2 2 − ) = qy
∂t 2
∂t ∂z ∂t ∂z∂t 2

∂ 4θ ∂ 2θ ∂ 2ξ ∂γ mx
EJ ϕ 4 − GJ d 2 + P(a y − e y ) − −
∂z ∂z ∂z 2 ∂z
∂ 2η ∂γ my ∂ 2θ
− P(a x − e x ) − − P ( − r 2
+ 2 β e
y y − 2 β e
x x ) (11c)
∂z 2 ∂z ∂z 2
∂ 2ξ ∂ 2η 2 ∂ θ
2
∂ 4θ
− ρFT (a y 2 − a x 2 + r ) − ρJ ϕ 2 2 = m
∂t ∂t ∂t 2 ∂z ∂t

En estas ecuaciones, FT es el área de la sección transversal, Jx y Jy son los momentos de

inercia de la sección transversal en relación a los ejes principales, Jϕ el momento de inercia
sectorial, Jd el módulo de torsión, ax y ay las coordenadas del centro de corte. ρ denota la
densidad de masa del material de la viga. E y G son los módulos de elasticidad longitudinal y
transversal respectivamente. Finalmente, γmx y γmy representan los valores medios de las
deformaciones por corte sobre una sección transversal z = constante.
El sistema representa un modelo general para vigas que tiene en cuenta el acoplamiento
triple entre las vibraciones flexionales y torsionales. Debe tenerse en cuenta que la ecuación
de vibración longitudinal relacionada al desplazamiento generalizado ζ (Figura 1) está
desacoplada del resto del sistema. En el caso que las vibraciones longitudinales sean de
interés pueden ser tratadas independientemente.

2.2 Método de variables de estado

Si se realiza la transformada de Fourier del sistema (11) se obtiene un sistema equivalente

de 12 ecuaciones diferenciales parciales con doce incógnitas en el dominio de la frecuencia.
( )
′′′ + P(ξ ′′ − γ mx ) + P (a y − e y )θ ′′ + ρFT ω 2ξ +
EJ y ξ IV − γ mx
(12a)
+ ρFT a y ω 2θ − ρJ y ω 2 (ξ ′′ − γ mx
′ ) = qx

2434
 ( )
′′′ + P(η ′′ − γ my ) − P(a x − e x )θ ′′ + ρFT ω 2η −
EJ x η IV − γ my
(12b)
− ρFT a x ω 2θ − ρJ x ω 2 (η ′′ − γ my
′ ) = qy

EJ ϕ θ IV − GJ d θ ′′ + P (a y − e y )(ξ ′′ − γ mx ) − P(a x − e x )(η ′′ − γ my ) −

(12c)
( ) ( )
− P − r 2 + 2 β y e y − 2 β x e x θ ′′ − − ρFT ω 2 a y ξ − a xη + r 2θ − ρJ ϕ ω 2θ ′′ = m

El esquema descrito anteriormente es conocido en la literatura como “método de variables

de estado”. Se seleccionan como componentes del vector v de estado seis componentes
geométricas y seis estáticas: Los desplazamientos ξ y η, las rotaciones de flexión φx y φy, los
esfuerzos de corte Qx y Qy, los momentos flectores Mx y My, la rotación torsional θ y su
derivada espacial θ´, el momento torsional total MT y el bimomento B.
v(z,ω) = {η, φy, Qy, Mx, ξ, φx, Qx, My, θ, θ’, MT, B}T (13)
Donde
MT = Hϕ + Hk (14)

Hk = GJdθ’ (15)
El sistema es
∂v
= Av (16)
∂z
Donde A es la matriz del sistema.

Las componentes del vector de estado v son designadas “variables de estado”. En el

dominio de la frecuencia las variables de estado dependen de la frecuencia ω y la variable
longitudinal z. Por simplicidad, se usa la misma notación para las variables de estado en el
dominio del tiempo y de la frecuencia, puesto que el dominio de trabajo puede identificarse
con la indicación de los argumentos de la función. Por ejemplo, η(z, t) y η(z,ω) se refieren a
los desplazamientos según el eje y en el dominio del tiempo y su transformada de Fourier,
respectivamente. Es importante notar que la presente formulación constituye una teoría
general de vigas aplicable tanto a vigas de pared delgada como a vigas sólidas.
En (17) k x′ y k ′y son los coeficientes de corte de Cowper y

Bθ’ = ρJϕω2 - GJd (17)

R=-r2+2βyey - 2βxex (18)

2435
 1
0 1 0 0 0 0 0 0 0 0 0
k ′y FG
1
0 0 0 − 0 0 0 0 0 0 0 0
EJ x
P − P(ac − ex )
− ρFω 2 0 0 0 0 0 0 ρFω 2 a x 0 0
EJ x EJ ϕ

0 ρJ xω 2 1 0 0 0 0 0 0 0 0 0

1
0 0 0 0 0 1 0 0 0 0 0
k x′ FG
1
0 0 0 0 0 0 0 0 0 0 0
EJ y
A=
−P P (a y − e y )
0 0 0 0 − ρFω 2 0 0 − ρFω 2 a y 0 0
EJ Y EJ ϕ

0 0 0 0 0 − ρJ yω 2
−1 0 0 0 0 0

0 0 0 0 0 0 0 0 0 1 0 0 (19)
1
0 0 0 0 0 0 0 0 0 0 0 −
EJ ϕ
− P(ac − ex ) − P (a y − e y ) PR
ρFω 2 ax 0 0 − ρFω 2 a y 0 0 − ρFω 2 r 2 0 0 −
EJ X EJ Y EJ ϕ

0 0 0 0 0 0 0 0 0 Bθ ′ 1 0

2.3 Procedimiento numérico y condiciones de contorno

El sistema (16) puede integrarse fácilmente usando procedimientos numéricos estándar,
tales como el método de cuarto orden Runge-Kutta, el algoritmo predector-corrector u otros
métodos. Para resolver el problema, tanto para frecuencias naturales como para la respuesta
dinámica, debe transformarse el mismo de un problema de valores de borde en dos puntos a
un problema de valores iniciales, como es demostrado por Ebner and Billington (1968). El
procedimiento es normalmente aplicado en el método de matrices de transferencia (Pestel and
Leckie, 1963). Las frecuencias naturales son determinadas por el método de Thomson.
En este trabajo se consideran las clásicas condiciones de contorno: empotrada, libre o
simplemente apoyada.
Borde empotrado
ξ = η = 0, φx = φy = θ = 0, θ´ = 0
Borde libre
Qy = Qx = 0, Mx = My = MT = 0, B = 0
Borde articulado
ξ = η = 0, Mx = My = 0, θ = 0, B = 0

3 RESULTADOS NUMÉRICOS
3.1 Modelo desarrollado
Se incorporó al programa computacional DYBEAM (Ambrosini et al, 1995, 2000) la
solución obtenida en la sección 2.
Se analiza una viga doblemente asimétrica de aluminio, con sección transversal presentada
en Figura 3.

2436
 40 mm

e = 2 mm 20 mm

100 mm

Figura 3: Sección transversal de viga de aluminio

Las propiedades mecánicas y geométricas usadas para aplicar la teoría presentada en este
trabajo son:
l = 2.00 m
FT = 3.2 10-4 m2
Jx = 3.024 10-8 m4 Jy = 4.595 10-7 m4
Jϕ = 2.838 10-11 m6 Jd = 4.270 10-10 m4
ax = 2.67 10-2 m ay = -1.30 10-2 m
E = 4.50 10 N/m ρ = 2650 kg/m3
10 2

ν = 0.25
P = 2000 N
La carga axial se aplicará en el centroide de la sección.

3.2 Modelo de Elementos Finitos

Con el objetivo de realizar una comparación con un modelo de elementos finitos, se utilizó
el programa SAP2000 (2001) para modelar la viga presentada anteriormente. Se utilizaron 40
elementos tipo Shell en la dirección longitudinal, considerando los efectos flexionales y
membranales. Las propiedades materiales fueron las mismas que las usadas en 3.1. A modo
de ilustración se presenta en figura 4 uno de los modos de vibración obtenidos en uno de los
ejemplos analizados.

Figura 4: Modo 2. Caso 1 Empotrado – Libre. SAP 2000

2437
3.3 Comparación de resultados
Se realiza en este punto la comparación entre los resultados numéricos. Debido a la
asimetría que presentan las vigas, todos los modos de vibración son modos acoplados de
flexión y torsión. Se analizan dos casos: Caso 1 Empotrado – Libre y Caso 2: Articulado –
Articulado. Los resultados obtenidos se presentan en Tablas 1 y 2.
Se observa que la máxima diferencia en la obtención de las frecuencias naturales
considerando la viga descargada y cargada es del orden del 25% con el programa DYBEAM.
La influencia de la carga axial se reduce en los modos superiores de vibración. Estos
resultados coinciden cualitativamente con los resultados obtenidos por L.Jun et al (2004).
Similares resultados son obtenidos con SAP 2000, lo que valida el modelo desarrollado.

Caso 1: Empotrado – Libre

DYBEAM SAP 2000
Modo f(Hz) f(Hz)
P=0 P=2000 P=0 P=2000
1 5.29 4.59 5.35 3.76
2 12.16 12.80 11.66 8.19
3 22.18 22.00 23.76 21.43
4 28.53 31.80 30.49 35.30
5 51.39 53.52 50.22 52.80
Tabla 1: Caso 1 Empotrado - Libre.

Caso 2: Articulado – Articulado

DYBEAM SAP 2000
Modo f(Hz) f(Hz)
P=0 P=2000 P=0 P=2000
1 13.48 18.14 14.09 17.20
2 25.28 28.04 24.45 32.09
3 43.85 50.05 45.65 48.79
4 61.26 62.55 69.39 59.28
5 85.52 88.50 84.33 65.89
Tabla 2: Caso 2 Articulado - Articulado.

En general, existe buena correlación entre las frecuencias obtenidas con DYBEAM y SAP
2000, tanto para el caso de viga descargada como con carga axial.

4 CONCLUSIONES
Se presentan en este trabajo las ecuaciones de movimiento de vigas de pared delgada y
sección abierta bajo carga axial, en el marco del método de variables de estado en el dominio
de la frecuencia. Las ecuaciones tienen en cuenta la influencia de la flexibilidad por corte e
inercias rotacionales, efectos éstos que son despreciados en la teoría original de Vlasov. Las
ecuaciones permiten el análisis de problemas prácticos usando integración numérica
conjuntamente con técnicas aplicadas rutinariamente en análisis de matrices de transferencia.
Además, la teoría propuesta puede usarse para vigas sólidas en las cuales exista acoplamiento
de las vibraciones flexionales y torsionales.

2438
 Se presenta una comparación numérica de las frecuencias naturales de vigas de pared
delgada y sección abierta doblemente asimétrica. Puede observarse en los resultados
obtenidos que el modelo presentado tiene una precisión comparable al obtenido con el método
de elementos finitos usando elementos de placa, que resulta de mayor costo computacional.
Por último, se aprecia el considerable efecto de la carga axial en los primeros modos de
vibración para todas las condiciones de borde consideradas.

Agradecimientos
Se agradece el apoyo financiero del CONICET y la Universidad Nacional de Cuyo.

REFERENCIAS
Aggarwal, H., and Cranch, E. “A Theory of Torsional and Coupled Bending Torsional Waves
in Thin-Walled Open Section Beams.” Journal of Applied Mechanics, ASME, 34, 337-343,
(1967).
Ali Hasan, S., and Barr, A. “Linear Vibration of Thin-Walled Beams of Equal Angle-
Section.” Journal of Sound and Vibration, 32, 3-23, (1974).
Ambrosini, R.D., Riera, J.D., and Danesi, R.F. “Dynamic Analysis of Thin-Walled and
Variable Open Section Beams with Shear Flexibility”, International Journal for Numerical
Methods in Engineering, 38(17), 2867-2885, (1995).
Ambrosini, R.D., Riera, J.D., and Danesi, R.F. “A modified Vlasov theory for dynamic
analysis of thin-Walled and variable open section beams”, Engineering Structures, 22, 890-
900, (2000).
Ambrosini D., and Danesi R., “Free vibration of nonsymmetrical thin walled beams”.
Mecánica Computacional Vol. XXII, XIII Congreso sobre Métodos Numéricos y sus
Aplicaciones, ENIEF 2003, pp. 666-677. (2003).
Arpaci, A., and Bozdag E. “On free vibration analysis of thin-walled beams with
nonsymmetrical open cross-sections” Computers and Structures, 80, 691-695, (2002).
Bishop, R.E.D., Cannon S.M. and Miao, S. “On Coupled bending and Torsional vibration of
uniform beams.” Journal of Sound and Vibration, 131(3), 457-64, (1989).
Christiano, P., and Culver, C. “Horizontally Curved Bridges Subject to Moving Load.”
Journal of Structural Division, ASCE, New York, 1615-43 (1969).
Ebner, A., and Billington, D. “Steady State Vibrations of Damped Timoshenko Beams.”
Journal of the Structural Division, ASCE, New York, 737-60, (1968).
Gere, J., and Lin, Y. “Coupled Vibrations of Thin-Walled Beams of Open Cross Section.”
Journal of Applied Mechanics, ASME, 25, 373-378, (1958).
Jun, L. et al. “Coupled Bending an torsional vibration of nonsymmetrical axially loades thin-
walled Bernoulli-Euler beam.” Mechanics Research Communications,31,697-711, (2004).
Jun, L.et al. “Bending-torsional coupled dynamic response of axially loaded composite
Timoshenko thin-walled beam with close cross section.” Composite Structures,64,23-35,
(2004).
Muller, P. “Torsional-Flexural Waves in Thin-Walled Open Beams.” Journal of Sound and
Vibration, 87(1), 115-141, (1983).
Pestel, E.C., and Leckie, F.A. Matrix Methods in Elastomechanics, McGraw-Hill, New York.
(1963).
SAP2000 Integrated Finite Element Analysis and Design of Structures, v. 7.42. Computers
and Structures Inc. (2001).
Tanaka, M., and Bercin A.N. “Free vibration solution for uniform beams of nonsymmetrical
cross section using Mathematica” Computers and Structures, 71, 1-8, (1999).

2439
Timoshenko, S., and Young, D. Vibration Problems in Engineering, 3rd. Ed., Van Nostrand,
Princeton, NJ. (1968)
Viola et al. Journal of Sound and Vibration,304, 124-153, (2007)
Vlasov, V. Thin-walled Elastic Beams, Israel Program for Scientific Translations, 2ed.
Jerusalem. (1963)
Yaman, Y. “Vibrations of Open-section Channels: A Coupled Flexural and Torsional Wave
Analysis.” Journal of Sound and Vibration, 204(1), 131-158, (1997).
Yoo, C., and Fehrenbach, J. “Natural Frequencies of Curved Girders.” Journal of the
Engineering Mechanics Division, ASCE, New York, 339-54, (1982).

2440