pandeo por flexotorsion - aero.ing.unlp.edu.ar · se calculó el momento torsor que absorbe un...
Post on 29-Jun-2018
221 Views
Preview:
TRANSCRIPT
Página 1 de 21
Facultad de Ingeniería
Universidad Nacional de La Plata
ESTRUCTURAS IV
FLEXOTORSION
Autores: Ing. Juan Pablo Durruty
Ing. Marcos D. Actis Ing. Alejandro J. Patanella
Página 2 de 21
INTRODUCCION: En los perfiles abiertos de paredes delgadas con dos ejes de simetría (centro de corte coincidente con el centro de gravedad), el pandeo puede darse según algunos de estos ejes o por torsión. En el caso de tener un sólo eje de simetría o ninguno, la pieza puede fallar por alguno de estos tres casos o por una combinación de flexión con torsión, lo que se denomina Pandeo por Flexotorsión. ECUACIÓN GENERAL DE PANDEO: Cuando se estudia la flexión simple se obtiene que:
∂∂Mz
Qxy= −
∂∂
2
2
Mz
qxy= −
y el Mx está dado por:
M y EJY
zEJ
Yz
EJ q
x x x
x y
= − ′′ = −
∴ =
∂∂
∂∂
2
2
4
4
(1)
análogamente para el otro eje:
∂∂
4
4
Xz
EJ qy x= (2)
Por otra parte, al estudiar el problema de estabilidad, se obtiene: según el eje y:
EJY
zP Yx cr
∂∂
2
2 = − (3)
Página 3 de 21
según el eje x:
EJX
zP Xy cr
∂∂
2
2 = − (4)
Derivando (3) y (4):
EJY
zP
Yzx cr
∂∂
∂∂
4
4
2
2= − y EJX
zP
Xzy cr
∂∂
∂∂
4
4
2
2= −
Reemplazando en (1) y (2):
q PY
zy cr= −∂∂
2
2 (5) y q PX
zx cr= −∂∂
2
2 (6)
Luego para encontrar la ecuación general de pandeo, provocaremos una perturbación infinitesimal en las direcciones x e y, compuestas con una rotación θ. (figura 1) M’ α M(xM,yM) Y θ β θ B” β G’ B’ F D G” B(xB,yB) ds β Y YG” t α G X X Figura 1 XG” Como se observa en la figura 1, se dan a la sección, los desplazamientos X e Y, y un giro θ que se verifica alrededor del centro de corte de la misma (M); pero como éste también se ha desplazado el giro es alrededor de su nueva posición (M’). Las nuevas coordenadas del centro de gravedad (G) serán:
X X G D X G GG ′′ = + ′ ≅ + ′ ′′senα
Página 4 de 21
pero:
′ ′′ ≅ ′ ′ =
→ = +∴ = +
G G G M GMX X GMX X Y
G
G M
θ θθ α
θ
"
"
sen
análogamente:
Y Y G D Y G GY Y GMY Y X
G
G
G M
"
"
"
coscos
= − ′′ ≅ − ′ ′′
= −= −
αθ αθ
Como la carga está ubicada en el centro de gravedad, al desplazarse éste, se generan los momentos:
( )M PY P Y Xx G M= = −" θ y ( )M PX P X Yy G M= = +" θ reemplazando en (3) y (4):
( )
( )
EJY
zP Y X M
EJX
zP X Y M
x M x
y M y
∂∂
θ
∂∂
θ
2
2
2
2
= − − = −
= − + = −
Considerando un punto cualquiera sobre la línea media de la sección, como por ejemplo el punto B (xB, yB), su desplazamiento será similar al de G:
Y Y Y B F Y B BB B B" sen= + + ′ ≅ + ′ ′′ β
( )
′ ′′ ≅ ′ ′ =
= + +
= + − −
B B B M BMY Y Y BM
Y Y Y X X
B B
B B M B
θ θθ β
θ
"
"
sen
(7)
(8)
(9) (10)
(11)
Página 5 de 21
( )
X X X B F
X X X Y Y
B B
B B M B
" cos= + + ′
= + + −
β
θ
Reemplazando (11) y (12) en (5) y (6):
( )[ ] ( )[ ]
( )[ ] ( )[ ]
q PY
zP
zY Y X X P
zY X X
q PX
zP
zX X Y Y P
zX Y Y
y cr cr B M B cr M B
x cr cr B M B cr M B
= − = − + + − = − − −
= − = − + + − = − + −
∂∂
∂∂
θ∂∂
θ
∂∂
∂∂
θ∂∂
θ
2
2
2
2
2
2
2
2
2
2
2
2
Tomando Momento torsor con respecto al centro de corte y siendo:
( )[ ]
( )[ ]
∂ ∂∂∂
θ ∂ σ
∂ ∂∂∂
θ ∂ σ
q Pz
Y X X P tds
q Pz
X Y Y P tds
y cr M B cr
x cr M B cr
= − − =
= + − =
2
2
2
2
y
y
( ) ( )[ ]( )
( )[ ]( )
∂ σ∂∂
θ
σ∂∂
θ
T z tdsz
Y X X X X
tdsz
X Y Y Y Y
M B M B
M B M B
= − − − −
− + − −
2
2
2
2
Integrando a lo largo de toda la línea media de la sección:
( ) ( ) ( )
( ) ( )
T z tY
zX X ds t
zX X ds
tX
zY Y ds t
zY Y ds
sM B M B
s
sM B M B
s
= − − − −
− − − −
∫ ∫
∫ ∫
σ∂∂
σ∂ θ∂
σ∂∂
σ∂ θ∂
2
2
2
2
2
2
2
2
2
2
(12)
Página 6 de 21
( ) ( )
( )
( )
( )
T z XY
ztds
Yz
tX dsz
X tds
zX tX ds
zt X ds
Xz
Y tds
Xz
tY dsz
Y tdsz
Y tY ds
zt Y ds
Ms
Bs
Ms
M BS s
B Ms
Bs
Ms
M BS
sB
= − − +
− − +
− + −
−
∫ ∫ ∫
∫ ∫ ∫
∫ ∫ ∫
∫
σ∂∂
σ∂∂
σ∂ θ∂
σ∂ θ∂
σ∂ θ∂
σ∂∂
σ∂∂
σ∂ θ∂
σ∂ θ∂
σ∂ θ∂
2
2
2
2
2
22
2
2
2
2
22
2
2
2
2
2
22
2
2
22
+ 2
+ 2
( ) ( )
( )
T z XY
zA
zX A
zJ
YX
zA
zY A
zJ
M M y
M M x
= − − −
− − −
σ∂∂
σ∂ θ∂
σ∂ θ∂
σ∂∂
σ∂ θ∂
σ∂ θ∂
2
2
2
22
2
2
2
2
2
22
2
2
( ) ( )T z P XY
zY
Xz
PA z
X A J Y A JM M M y M x= −⎛
⎝⎜
⎞
⎠⎟ − + + +
∂∂
∂∂
∂ θ∂
2
2
2
2
2
22 2
Llamando: ( )J X Y A J JM M y x02 2= + + +
( )T z P XY
zY
Xz
PA z
JM M= −⎛
⎝⎜
⎞
⎠⎟−
∂∂
∂∂
∂ θ∂
2
2
2
2
2
2 0 (13)
T(z) es un momento torsor distribuido en la longitud de la pieza. En la “referencia 1” se calculó el momento torsor que absorbe un perfil abierto de paredes delgadas con alabeo restringido, cuando el momento aplicado es constante a lo largo de la pieza.
T C Ez
GJzM t= −
∂ θ∂
∂θ∂
3
3
Por lo tanto para obtener la expresión correspondiente al caso de un momento distribuido se hace:
Página 7 de 21
T(z) T T + dT dz dz dz Figura 2
( )T z dz TdTdz
dz TdTdz
C Eddz
GJddzM t= + − = = −
4
4
2
2
θ θ (14)
Reemplazando (13) en (14):
P XY
zY
Xz
PA z
J C Eddz
GJddzM M M t
∂∂
∂∂
∂ θ∂
θ θ2
2
2
2
2
2 0
4
4
2
2−⎛
⎝⎜
⎞
⎠⎟− = −
C Eddz
GJ JPA
ddz
PXY
zPY
XzM t M M
4
4 0
2
2
2
2
2
2 0θ θ ∂
∂∂∂
− −⎛⎝⎜
⎞⎠⎟ − + = (15)
Resolvemos la ecuación general para el caso que presenta las siguientes condiciones de borde: en z = 0 y z = L : X = Y = 0; θ = 0 y además Mx = My = 0
o sea: d Ydz
d Xdz
ddz
2
2
2
2
2
2 0= = =θ
en z = 0 y z = L
Proponiendo la siguiente solución:
X Az
LY A
zL
Az
L=
⎛⎝⎜
⎞⎠⎟ =
⎛⎝⎜
⎞⎠⎟ =
⎛⎝⎜
⎞⎠⎟1 2 3sen sen sen
π πθ
π
Reemplazando en (9) y (10):
Página 8 de 21
−⎛⎝⎜
⎞⎠⎟
⎛⎝⎜
⎞⎠⎟ = −
⎛⎝⎜
⎞⎠⎟−
⎛⎝⎜
⎞⎠⎟
−⎛⎝⎜
⎞⎠⎟
⎡
⎣⎢
⎤
⎦⎥ − =
EJL
Az
LPA
zL
PX Az
L
P EJL
A PX A
x M
x M
π π π π
π
2
2 2 3
2
2 3 0
sen sen sen
−⎛⎝⎜
⎞⎠⎟
⎛⎝⎜
⎞⎠⎟ = −
⎛⎝⎜
⎞⎠⎟−
⎛⎝⎜
⎞⎠⎟
−⎛⎝⎜
⎞⎠⎟
⎡
⎣⎢
⎤
⎦⎥ − =
EJL
Az
LPA
zL
PY Az
L
P EJL
A PY A
y M
y M
π π π π
π
2
1 1 3
2
1 3 0
sen sen sen
de la (15):
C EL
Az
LGJ J
PA L
Az
L
PXL
Az
LPY
LA
zL
M t
M M
π π π π
π π π π
⎛⎝⎜
⎞⎠⎟
⎛⎝⎜
⎞⎠⎟+ −
⎛⎝⎜
⎞⎠⎟⎛⎝⎜
⎞⎠⎟
⎛⎝⎜
⎞⎠⎟+
+⎛⎝⎜
⎞⎠⎟
⎛⎝⎜
⎞⎠⎟−
⎛⎝⎜
⎞⎠⎟
⎛⎝⎜
⎞⎠⎟ =
4
3 0
2
3
2
2
2
1 0
sen sen
sen sen
PY A PX A C EL
GJPA
J AM M M t1 2
2
0 3 0− −⎛⎝⎜
⎞⎠⎟ + −
⎡
⎣⎢
⎤
⎦⎥ =
π (18)
Se obtuvo entonces un sistema de 3 ecuaciones [(16), (17) y (18)] con tres incógnitas (A1, A2 y A3); para encontrar una solución única el determinante deberá ser nulo:
0
0 0
2
2
0
2
P EJL
PX
P EJL
PY
PY PXPA
J C EL
GJ
x M
y M
M M M t
−⎛⎝⎜
⎞⎠⎟
⎡
⎣⎢
⎤
⎦⎥ −
−⎛⎝⎜
⎞⎠⎟
⎡
⎣⎢
⎤
⎦⎥
− −⎛⎝⎜
⎞⎠⎟ −
⎡
⎣⎢
⎤
⎦⎥
=
π
π
π
(16)
(17)
Página 9 de 21
( )
( )
P EJL
PY P EJL
PA
J C EL
GJ
PX P EJL
x M y M t
M y
−⎛⎝⎜
⎞⎠⎟
⎡
⎣⎢
⎤
⎦⎥ − −
⎛⎝⎜
⎞⎠⎟
⎡
⎣⎢
⎤
⎦⎥ −
⎛⎝⎜
⎞⎠⎟ −
⎡
⎣⎢
⎤
⎦⎥
⎧⎨⎩
⎫⎬⎭+
+ −⎛⎝⎜
⎞⎠⎟
⎡
⎣⎢
⎤
⎦⎥ =
π π π
π
22
2
0
2
22
0
Considerando ahora el caso en que el centro de corte coincida con el centro de gravedad, de las ecuaciones (16), (17) y (18) se obtiene:
PL
EJ
PL
EJ
PAJ
C EL
GJ
crx x
cry y
cr M t
=⎛⎝⎜
⎞⎠⎟
=⎛⎝⎜
⎞⎠⎟
=⎛⎝⎜
⎞⎠⎟ +
⎡
⎣⎢
⎤
⎦⎥
π
π
πθ
2
2
0
2
Estas ecuaciones dan las cargas críticas de pandeo por flexión en ambos ejes y por torsión. La menor de estas tres será la carga crítica de la columna y determinará entonces la forma de pandeo de la misma. Considerando las (20), (21) y (22); y desarrollando (19) en P, se obtiene la ecuación general de pandeo por flexotorsión:
( )
( )[ ]P
J JJ
P P P P P Y P XAJ
P P P P P P P P P
x y
ocrx cry cr crx M cry M
cr crx cry crx cry crx cry cr
3 2 2 2
0
0
+⎡
⎣⎢⎤
⎦⎥− + + − +
⎡
⎣⎢⎤
⎦⎥+
+ + + − =
θ
θ θ (23)
Para el caso de tener un solo eje de simetría se debe resolver el determinante reducido, suponiendo YM = 0, se obtiene:
( )
[ ] [ ]
P EJL
PA
J C EL
GJ PX
P P P PJA
P X
x M t M
crx cr M
−⎛⎝⎜
⎞⎠⎟
⎡
⎣⎢
⎤
⎦⎥ −
⎛⎝⎜
⎞⎠⎟ −
⎡
⎣⎢
⎤
⎦⎥ + =
− − + =
π π
θ
2
0
22
0 2 2
0
0
(19)
(20) (21) (22)
Página 10 de 21
( )P XAJ
P P P P PM crx cr crx cr2 2
01 0−⎛
⎝⎜
⎞
⎠⎟− + + =θ θ (24)
Para el caso de tener dos ejes de simetría, las cargas se calculan por las ecuaciones (20), (21) y (22). La carga crítica será entonces la menor de las tres. EJEMPLOS DE PANDEO POR FLEXOTORSION: 1) Sección con dos ejes de simetría: Datos: E = 75000 N/mm2 t G = 21000 N/mm2 σp = 100 N/mm2 h L = 5000 mm t = 2 mm h = 80 mm b = 40 mm A = 320 mm2 b Momentos de inercia: Jx = 2.(40.23/12 + 40.2.402) + 2.803/12 = 341386 mm4 Jy = 2.2.403/12 = 21333 mm4 Condiciones de borde: En ambos ejes la condición es articulado - articulado, por lo tanto, C = 1
Página 11 de 21
λπσp
p
C E= =
2
86
ρx = 32,66 mm λx = 153 > 86 ⇒ ρy = 8,165 mm λy = 612 > 86 Por lo tanto en ambos ejes, el pandeo será elástico CALCULO DE LAS CARGAS CRITICAS: σcrx = C.π2.E / λx
2 = 31,62 N/mm2 Pcrx = σcrx . A = 10118 N σcry = C.π2.E / λy
2 = 1,976 N/mm2 Pcry = σcry . A = 632,42 N Coeficiente de torsión: CM = tf . h2 . b3 = 34133333 mm6, tf = t 24 El centro de corte coincide con el baricéntrico ⇒ x0 = y0 = 0 Pcrθ = A(G Jt + E π2 CM) , donde: J0 = Jx + Jy + A(x0
2 + y02) = 362720 mm4
J0 L2
Jt = t3 (2b + h) = 426,66 mm4 ⇒ Pcrθ = 8796,2 N 3 En secciones con dos ejes de simetría, la carga crítica por flexotorsión es la menor de las tres anteriormente calculadas, ⇒ 2) Sección con un eje de simetría: Datos: E = 70000 N/mm2 t G = 27000 N/mm2 σp = 100 N/mm2 h L = 800 mm t = 1 mm (tf = tw = t) h = 20 mm b = 15 mm A = 50 mm2 b CALCULO DEL CENTRO DE GRAVEDAD: La ordenada del baricentro se encuentra en el eje de simetría.
Pcr = 632,42 N
Página 12 de 21
La abcisa es: xg = 2.15.7,5 / 50 = 4,5 (medida a partir del tabique vertical) CALCULO DE MOMENTOS DE INERCIA, COEFICIENTE DE TORSION Y UBICACION DEL CENTRO DE CORTE: Sistema C.G. Jx = 2.(15.13/12+15.1.102)+1.203/12 = Jx = 3669,16 mm4 c.c G Jy = 2.(1.153/12+15.1.32)+20.13/12+20.1.4,52 = x x Jy = 1239,16 mm4 e e = 3.b2.tf = 6,136 mm x0 6.b.tf+h.tw x0 = xg + e = 10,636 mm, y0 = 0 CM = tf.b3.h2 . (3.b.tf+2.h.tw) = 86931,8 mm6
12 (6.b.tf+h.tw) ρx = 8,566 mm λx = 93,4 ⇒ ρy = 4,978 mm λy = 160,7 Diagrama de área sectorial principal: w0 = 0 1 2 w1 = -10.e = 61,36 P - c.c w2 = -61,36 + 10.b = 88,64 0 Se verifica que el punto donde w =0, en las alas está a una distancia “e” de 1. CALCULO DE LAS CARGAS CRITICAS: Condiciones de borde: Según el eje X (de simetría), la condición es articulado - articulado ⇒ Cx = 1 Según el eje Y, la barra esta empotrada-empotrada ⇒ Cy = 4 λpx
2 = Cx . π2 . E / σp ⇒ λpx = 83,12 < λx λpy
2 = Cy . π2 . E / σp ⇒ λpy = 166,23 > λy Por lo tanto en el eje X el pandeo será elástico y en el eje Y será inelástico.
Página 13 de 21
σcrx = Cx . π2 . E / λx
2 = 79,2 N/mm2 ⇒ Pcrx = 3960 N σcry = σp = 100 N/mm2 ⇒ Pcry = 5000 N De datos y cálculos anteriores: J0 = 10564,54 mm4 Jt = 50/3 mm4 ⇒ Pcrθ = 2574 N G = 27000 N/mm2 Para secciones con un solo eje de simetría, la carga crítica la da la ecuación (24), en la que si reemplazamos los valores obtenidos de Pcrx y Pcrθ, obtenemos la siguiente ecuación cuadrática: Pcr
2 - 14063,62. Pcr + 21939247,06 = 0 de donde surgen los valores: Pcr1 = 12276,53 N Pcr2 = 1787 N La carga crítica será la menor de todas las calculadas anteriormente: Pcr = 1787 N
Página 14 de 21
3) Cálculo de la carga crítica para los ejemplos anteriores, por el programa de elementos finitos MSC Nastran. Se cargaron los perfiles con tres cargas puntuales de valor P = 1N, en los puntos donde el área sectorial es cero. Estos son los puntos donde no hay alabeo de la sección. e Y X Analizando las barras a pandeo, el programa da un valor (Eigenvalue) el cual se multiplica por la carga aplicada. En nuestro caso: EV . 3 . 1N = Pcr Ejemplo 1) Restricciones en ambos extremos: (T = traslaciones, R = rotaciones) Tx = 0 Ty = 0 Tz = libre Rx = libre Ry = libre Rz = 0 EV = 212,131 ⇒ Pcr = 636,4 N
Ejemplo 2)
Página 15 de 21
Restricciones en ambos extremos: (T = traslaciones, R = rotaciones) Tx = 0 Ty = 0 Tz = libre Rx = libre Ry = 0 Rz = 0 Ry = 0 considerando empotrado en el eje Y EV = 659,913 ⇒ Pcr = 1980 N
La discrepancia de valores en este caso, puede deberse a que si bien se restringieron los nodos con Ry = 0, en el resultado del análisis se aprecia un giro del plano de la sección en el eje Y.
Página 16 de 21
VERIFICACION EXPERIMENTAL Se ensayaron cinco perfiles de chapa doblada de aluminio 5052, de geometría similar a la del ejemplo 2. Para obtener las condiciones de borde planteadas en dicho ejemplo, se construyeron un par de mordazas con las características que se describen en las siguientes figuras: 120 20 Medidas en mm 60 15 20 80 20 30 10 30 10 Como se aprecia en las figuras, mediante estas mordazas se permite el giro en un eje (condición articulado - articulado) mientras que en el otro, se encuentra impedido (empotrado - empotrado). Además en la plataforma de sujeción del perfil se encuentran dispuestos tres elementos que toman al perfil en los puntos donde el diagrama de área sectorial es nulo (condición de alabeo impedido), como se muestra en el siguiente diagrama:
Página 17 de 21
Los ensayos se realizaron en un pórtico sujetando una de las mordazas en la viga horizontal superior, mientras que sobre la otra se aplicó la carga por medio de un pistón hidráulico. La fuerza aplicada fue medida mediante un manómetro.
Página 18 de 21
Características de los perfiles ensayados: Sección: t = 1 mm h = 20 mm b = 15 mm Resultados obtenidos:
Perfil Longitud (m) Carga Crítica Kg (N) 2 0,8 117 (1146,6) 3 0,8 128 (1254,4) 5 0,8 112 (1097,6) 1 1,0 96 (940,8) 4 1,0 96 (940,8)
Página 19 de 21
CONCLUSIONES Las diferencias entre los valores calculados y los obtenidos en los ensayos se deben
a problemas como centrado de carga, doblado de los perfiles (diferencias con la
sección transversal teórica, y condiciones de borde reales no exactamente iguales a
las teóricas.
No obstante, se pudo apreciar el fenómeno de flexotorsión ya que los perfiles
fallaron en la forma prevista por los cálculos.
Debe señalarse que en las fotos anteriores se aprecia en los perfiles, fallas del tipo
de pandeo local, esto se debe a que luego de haber fallado por flexotorsión, se les
continuó aplicando carga hasta que se verificaron las fallas mencionadas.
Página 20 de 21
REFERENCIAS: 1- CALCULO DE TENSIONES NORMALES EN PERFILES DELGADOS ABIERTOS CON ALABEO IMPEDIDO, SOMETIDOS A TORSION, Ing. Pisoni, Ing. Actis. BIBLIOGRAFIA: 1- STATIK UND STABILITAT DER BAUKONSTRUKTIONEN - Dr. Ing. Christian Petersen. 2- RESISTENCIA DE MATERIALES - V. I. Foedosiev. 3- DISEÑO DE ESTRUCTURAS DE ACERO - Bresler, Lin y Scalzi.
Página 21 de 21
TABLA I : Ubicación del centro de corte y coeficiente de alabeo para varios perfiles delgados
b tf tw o h CM = tf h2 b3 h/2 24 b1 tf tw o e h e = h b1
3 CM = tf h2 b13 b2
3 b1
3 + b23 12 b1
3 + b23
b2 b e tf o h e = 3 b2 tf CM = tf b3 h2 3 b tf + 2 h tw 6 b tf + h tw 12 6 b tf + h tw tw b tw tf CM = b3 h2 [ 2 tf (b2 + b h + h2) + 3 tw b h] o h 12 (2 b + h2) h/2 a t e = 2 a sen(α) - α cos(α) α - sen(α) cos(α) α o α CM = 2 t a5 α3 - 6 [sen(α) - α cos(α)]2 3 α - sen(α) cos(α) e
TABLA I : Continuación
a
Página 22 de 21
t CM = A3 A = área de la sección b e 144 o t b2 CM = t3 (b1
3 + b23)
b1 o 36 t2 d CM = t13 b3 + t23 d3 c 144 36 t1 o b b1 e = b b1
2 (3 b - 2 b1) c √ 2 [2b3 - (b - b1)3] b e o t w o e = 3 (b2 - b1
2) h (w/t) h + 6 (b + b1) e t b1 < b
b1 b
top related