Rev. Int. Met. Num. Calc. Dis. Ing.

Vol. 18, 1, 5977 (2002)

Revista Internacional de
Metodos Numericos para
Calculo y Diseno en Ingeniera

Pandeo elastico y diseno de estructuras espaciales

Antonio Aguero y Jose Ramon Atienza
Departamento de Mecanica de los Medios Continuos y Teora de Estructuras
ETS de Ingenieros de Caminos, Canales y Puertos
Universidad Politecnica de Valencia
Campus de Vera s/n, 46020 Valencia, Espana
Tel.: 34-96-387 76 75, Fax: 34-96-387 76 75
e-mail: [email protected]
e-mail: [email protected]

Resumen
El problema del pandeo elastico en estructuras espaciales con secciones doblemente simetricas bajo la accion
de unos esfuerzos primarios de axil y exion segun el eje fuerte es resuelto con el programa CRIT3D. Para
ello se calcula el efecto de los esfuerzos primarios sobre los terminos de la matriz de rigidez asociados a la
exo-torsion, planteando el equilibrio en una posicion ligeramente deformada. Con este programa se evalua
la precision de las cargas crticas obtenidas al emplear las formulas indicadas en el Eurocodigo 3. Por ultimo
se propone un metodo de diseno en el que los esfuerzos de calculo se obtienen planteando el equilibrio en
una posicion ligeramente deformada de un sistema con imperfecciones geometricas.

ELASTIC BUCKLING AND DESIGN OF SPACE FRAMES

Summary
The problem of elastic buckling of space frames with bi-symmetrical cross section taking into account the
eect of axial forces and bending moment about the major principal axis, is solved using the computer
program CRIT3D. First of all we calculate the eect of the primary eorts on the coecients of the stiness
matrix due to the exural-torsional eect, by using the equilibrium equations in a slightly deformed position.
With this program we can evaluate how accurate are the EC-3 formulas to estimate the critical loads. At
last we suggest a method to obtain the design eorts using the equilibrium equations in a slightly deformed
position considering the eect of geometrical imperfections.

INTRODUCCION

La gran utilizacion de vigas con secciones abiertas de pared delgada en la construccion

metalica ha llevado al estudio de su comportamiento estructural, que resulta ser muy
complejo por ser las vigas susceptibles de sufrir diversas formas de inestabilidad. El modo
de inestabilidad que se presente en cada caso viene condicionado por el tipo de carga y
las condiciones de apoyo. La Figura 1a muestra la inestabilidad por exion sobre el eje
fuerte o debil, la Figura 1b la inestabilidad por torsion y la Figura 1c la inestabilidad por
exo-torsion.
Con objeto de vericar el estado lmite ultimo de una estructura se tiene que comprobar
que resiste las posibles combinaciones de estados de carga a las que puede estar sometida.
Para realizar el calculo con exactitud hay que hacer un analisis no lineal geometrico y
mecanico, teniendo en cuenta que cada viga real tiene unas tensiones residuales e imperfec-
ciones tanto geometricas como del material. Este tipo de problemas no se puede abordar

Universitat
c Politecnica de Catalunya (Espana). ISSN: 02131315 Recibido: Febrero 2001
60 A. Aguero y J.R. Atienza

de forma analtica por lo que se hace necesario recurrir a otros metodos, como puede ser
el de los elementos nitos. Como el metodo anterior resulta excesivamente laborioso, las
normas actuales de diseno proponen disenar con los esfuerzos obtenidos de un calculo lineal
amplicandolos por unos coecientes mediante unas formulas con base analtica y experi-
mental. Para determinar estos coecientes se utiliza la esbeltez de la pieza, obteniendose la
misma para cada barra en funcion la carga crtica del sistema perfecto. Para la obtencion
de estos coecientes se han contrastado las cargas crticas con los ensayos realizados sobre
piezas reales concluyendo que la carga ultima es inferior a la del sistema perfecto equiva-
lente. Esto se aprecia al introducir en un sistema perfecto una imperfeccion geometrica, lo
que transforma los puntos de bifurcacion en partes lmite (Figura 2).

Figura 1. Modos de inestabilidad

Figura 2. Efectos de la magnitud de imperfeccion {} sobre la carga crtica {cr } en un

punto de bifurcacion inestable simetrico


Segun el EC-3 la esbeltez de una pieza sometida a axil es = A y
y sometida a exion

Ncr

es = W y
Mcr
. A su vez, el EC-3 propone formas de estimar las cargas crticas tanto para
esfuerzos de axil como de exion. En el caso de vigas biapoyadas doblemente simetricas
sometidas a exion uniforme, con el giro de torsion impedido y alabeo libre en los extremos,
el momento crtico viene dado por (1). Para distintas condiciones de vinculacion {k, kw} y
distintos tipos de carga {C1 } se emplea la ecuacion (2).
 
2 EIz Iw L2 GIt
Mcr = + 2 (1)
L2 Iz EIz
Pandeo elastico y diseno de estructuras espaciales 61


  2
EIz 
2
 k Iw (kL) 2 GI
t
Mcr = C1 + (2)
(kL)2 kw Iz 2 EIz

La ecuacion (2) es la solucion exacta en los casos en los que k = kw y C1 = 1, para el

resto de los casos es una aproximacion al momento crtico. Mas adelante se presenta una
metodologa alternativa de diseno mediante la cual los esfuerzos de calculo se obtienen con-
siderando que tratamos con un sistema imperfecto y realizamos el analisis en una geometra
ligeramente deformada. Esta ultima hipotesis es correcta, ya que en el estado lmite de
servicio imponemos que los movimientos han de ser pequenos.
Con objeto de facilitar la obtencion de cargas crticas para estructuras complejas se
ha desarrollado un programa de ordenador Crit3d que puede llevar a cabo el analisis de
estructuras espaciales con esfuerzos primarios de exion segun eje fuerte y axil. Para ello se
tiene que obtener previamente la matriz de rigidez en una posicion ligeramente deformada,
que se deduce de las correspondientes ecuaciones de equilibrio y sus condiciones de contorno.

Figura 3. Movimientos y esfuerzos en ejes locales

Las ecuaciones de equilibrio en una posicion ligeramente deformada para una viga de
seccion doblemente simetrica inicialmente en equilibrio con unos desplazamientos {d} bajo
la accion de un axil y un momento de eje y. Se obtienen planteando el equilibrio en una
posicion perturbada con unos desplazamientos d = {u, v, w, x } del centro de esfuerzos
cortantes:7
2
A) Asociado al axil: EA xu2 = 0. Tendremos dos grados de libertad {ud, uf }.
4 2
B) Efecto del axil sobre la exion primaria: Ely xw4 + P xw2 = 0. Esta inuencia es tenida
en cuenta a traves de las funciones de estabilidad.4 Tendremos cuatro grados de libertad
{wd, y d, wf, y f }.
C) El efecto del axil y la exion primaria sobre la exion debil y torsion se puede
obtener mediante unas funciones de estabilidad (Apendice IV) si se desprecia la rigidez
al alabeo. Si no se desprecia, la matriz de rigidez explcita no se puede expresar
de forma tan compacta y es necesario recurrir a un programa de calculo simbolico
para su obtencion (ver siguiente apartado). Tendremos ocho grados de libertad
{vd, x d, z d, d, vf, x f, z f, f }. Siendo el alabeo.
62 A. Aguero y J.R. Atienza

CALCULO DE MATRIZ DE RIGIDEZ PARA LA FLEXION LATERAL CON

ESFUERZOS PRIMARIOS DE {AXIL Y FLEXION}

Las ecuaciones de equilibrio para la exion lateral, que se obtienen del desarrollo del
apartado anterior, son
4v 2 x 2v
EIz 4
+ My 2 + P 2 = 0 (3)
x x x

 2 x 2v 4 x
GIt + P ro2 + M y + EIa =0 (4)
x2 x2 x4
Resolviendo el sistema de ecuaciones diferenciales acopladas obtenemos los desplaza-
mientos v y x en funcion de ocho constantes indeterminadas {a, b, c, d, k, l, m, n}.
v = aex + bex + cex + dex + k + lx (5)

(EIz 2 P ) x (EIz 2 P ) x
x = (ae + bex ) + (ce + dex ) + m + nx (6)
My My
siendo {, , , } las soluciones de la ecuacion caracterstica de (3) y (4). Ecuacion
polinomica de cuarto grado en

(E 2 Iz Ia )4 + (EIz P ro2 EIz GIt + P EIa )2 + (P 2 ro2 My2 P GIt ) = 0 (7)

De las ecuaciones (5), (6) y (7) obtenemos los valores de los desplazamientos v, x y
sus derivadas respectivas z y en el extremo dorsal x = 0 y en el frontal x = L.
Los esfuerzos asociados a estos desplazamientos son:
2v
Mz = EIz (8)
x2
x v 3 x
Mx = (GIt P ro2) My EIa 3 (9)
x x x
2 x
Bi = EIa (10)
x2
3v x v
Qy = EIz 3
My P (11)
x x x
Sustituyendo en (8) a (11) los desplazamientos se obtienen los valores de los momentos
{Mz , Mx }, cortante {Qy } y bimomento {Bi} en el extremo frontal x = L y en el extremo
dorsal x = 0, tras cambiar de signo a los del dorsal. Se obtienen tanto los esfuerzos como
los movimientos en funcion de las ocho constantes indeterminadas {a, b, c, d, k, l, m, n}.

vd a Qy d a



b


b

x d



M x d






c


c


z d




Mz d



d d Bid d

=
=

vf
Matriz 1
k Qy f
Matriz 2 k

xf
l




M xf


l











z f
m
Mz f
m

f n Bif n
Pandeo elastico y diseno de estructuras espaciales 63

Para obtener la matriz de rigidez no hay mas que expresar las constantes indetermi-
nadas {a, b, c, d, k, l, m, n} en funcion de los grados de libertad {vd, x d, z d, d, vf, x f, z f,
f }K= (Matriz 2) (Matriz 1)1 . Esta se obtiene de forma explcita con la ayuda de algun
programa de calculo simbolico (MapleV) (ver Apendice II).
En el caso de despreciar la rigidez al alabeo la expresion resultante para la matriz de
rigidez es mas simplicada y se puede ver en el Apendice IV.

PROCESO DE CALCULO DE CARGAS CRITICAS

Esquema general
La estructura se subdivide en elementos lineales de dos nodos. El criterio de subdivision
vendra condicionado por la forma geometrica y tipo de carga. La busqueda de la carga
crtica supone encontrar la magnitud de los esfuerzos primarios para los que existe mas
de una solucion que satisfaga las condiciones de equilibrio. Esto supone que la variacion
segunda del potencial deja de ser denida positiva, lo que se representa en la Figura 4
para un sistema de 2 GDL y es equivalente a encontrar la carga para la cual la matriz
de rigidez es singular presentando un autovalor nulo. Tras un analisis lineal en el que se
obtienen los esfuerzos primarios de axil y exion, se generan las matrices de rigidez de cada
elemento teniendo en cuenta las desconexiones existentes. Estas se pasan de ejes locales a los
globales, se ensamblan formando la matriz de rigidez global y se busca para que magnitud
de los esfuerzos la matriz de rigidez pasa a tener un autovalor nulo.

Figura 4. Representacion de la variacion del incremento de potencial para un sistema de 2GDL

Analisis lineal
Se realiza un calculo lineal de la estructura obteniendo de este los esfuerzos primarios
de axil y exion segun el eje fuerte.
64 A. Aguero y J.R. Atienza

Matriz de rigidez ejes locales

Teniendo en cuenta los esfuerzos primarios escalados por un factor se calcula la
matriz de rigidez en una geometra ligeramente deformada.

Paso de ejes locales a ejes globales

En el Crit3d se ha considerado para la obtencion de los ejes locales que el x coincide
con la directriz de la barra y el eje y local esta contenido en el plano eje x y, permitiendo
un giro gama sobre el eje x para admitir distintas orientaciones del eje fuerte. El calculo de
la matriz de rotacion esta en el Apendice 2.
En cuanto al paso de ejes locales a globales cabe hacer mencion a la variable alabeo.
Para el caso de que en el nudo se encuentren dos barras de igual seccion y alineadas habra
continuidad de alabeos. En el resto de los casos dependera de como esta pensado ejecutar
las uniones, pudiendo existir varios grados de libertad de alabeo por nudo.

Calculo de autovalores
Se busca el valor de para el cual se produce el primer autovalor nulo, estando por
tanto en un punto de bifurcacion. Para ello, en caso de no estar los sucientemente proximos
al autovalor nulo, volvemos a empezar modicando coherentemente el factor .
El modo de pandeo asociado a este autovalor sera su autovector.

EJEMPLOS
En los ejemplos 1 y 2 se analiza una variacion lineal del momento My = M (M
B M )/L x y se estudia lo que ocurre para B = 1 (momento uniforme), B = 0, 5, B = 0
y B = 0, 5.

Ejemplo 1
En este caso las condiciones de vinculacion son {vd = 0, x d = 0, vf = 0, x f = 0}.
{Mz d = 0, Bid = 0, Mz f = 0, Bif = 0} {EC3 k = kw = 1} (Figura 5 y Figura 13).

Figura 5. {vd = 0, x d = 0, vf = 0, x f = 0}. {Mz d = 0, Bid = 0, Mz f = 0, Bif = 0}

Pandeo elastico y diseno de estructuras espaciales 65

Si B = 1, se superponen los resultados obtenidos con el programa y el Eurocodigo y ambos

dan la solucion analtica. Si B = 0, 5 y 0, las soluciones no son identicas, pero se parecen.
Si B = 0, 5, existe una diferencia apreciable entre ambos resultados, los resultados del
programa pueden ser erroneos porque se requiere aumentar la discretizacion. Al aumentar
el numero de elementos (Figura 6) vemos, que la solucion obtenida sigue quedando por
debajo de la del Eurocodigo, por lo que estamos del lado de la inseguridad utilizando el
Eurocodigo

Figura 6. {vd = 0, x d = 0, d = 0, vf = 0, x f = 0, f = 0}. {Mz d = 0, Mz f = 0}

En la Figura 7 se puede ver la conguracion de pandeo para el caso B = 0, 5.

Figura 7. {vd = 0, x d = 0, vf = 0, x f = 0}. {Mz d = 0, Bid = 0, Mz f = 0, Bif = 0}

Ejemplo 2
En este caso las condiciones de vinculacion son {vd = 0, x d = 0, d = 0, vf =
0, x f = 0, f = 0}. {Mz d = 0, Mz f = 0}. {EC3 k = 1, kw = 0, 5} (Figura 8 y
Figura 14).
66 A. Aguero y J.R. Atienza

Figura 8. {vd = 0, x d = 0, d = 0, vf = 0, x f = 0, f = 0}. {Mz d = 0, Mz f = 0}

Se aprecia que para el caso B = 1, momento uniforme la solucion del Eurocodigo es

inferior a la exacta dada por el programa. Esto es debido a que k = kw, quedando en este
caso el Eurocodigo del lado de la seguridad.
En la Figura 9 se puede ver la conguracion de pandeo para el caso B = 0, 5.

Figura 9. {vd = 0, x d = 0, d = 0, vf = 0, x f = 0, f = 0}. {Mz d = 0, Mz f = 0}

Ejemplo 3
En este caso las condiciones de vinculacion son {vd = 0, x d = 0, z d = 0, vf =
0, x f = 0, z f = 0}. {Bid = 0, Bif = 0}. {EC3 k = 0, 5; kw = 1}. Se utiliza un
unico elemento, ya que el momento es constante y por tanto obtenemos la solucion analtica.
Se aprecia en las Figuras 10 y 15 el error cometido por el Eurocodigo, ya que k = kw, la
solucion del Eurocodigo queda del lado de la seguridad.
Pandeo elastico y diseno de estructuras espaciales 67

Figura 10. {vd = 0, x d = 0, z d = 0, vf = 0, x f = 0, z f = 0}. {Bid = 0, Bif = 0}

Ejemplo 4
Se estudia la conguracion de pandeo para una estructura formada por dos arcos
empotrados en los extremos y sometidos a un axil uniforme; luz 10 m, echa vertical 5
m, separacion entre apoyos 5 m estando formados los mismos por: Caso 1 Perl abierto
(IPN-300) Ncr = 615KN (Figura 11); Caso 2 Perl cerrado (#200 150) Ncr = 2059KN
(Figura 12). Calculos realizados con el CRIT 3D.

Figura 11. Modo de pandeo. Arcos perl abierto IPN300

Figura 12. Modo de Pandeo. Arcos perl cerrado rectangular 200 150
68 A. Aguero y J.R. Atienza

Obtencion de esfuerzos de calculo para la comprobacion del estado lmite ultimo

siendo:
K = matriz de rigidez lineal
KN L (d) = fuerzas internas formulacion nolineal exacta
KN = matriz de rigidez tangente en una posicion ligeramente deformada
de un sistema perfecto
KNI d + F p = fuerzas internas en una posicion ligeramente deformada
de un sistema imperfecto

Grado de dicultad
bajo medio muy alto
Pandeo elastico y diseno de estructuras espaciales 69

ESTUDIO DE SISTEMAS IMPERFECTOS

El objeto del estudio de los sistemas imperfectos es el obtener los esfuerzos de calculo
que sirven para la comprobacion del estado lmite ultimo, realizando los calculos de un
sistema imperfecto en una geometra ligeramente deformada. La imperfeccion se considera
unicamente geometrica.

Definicion de las imperfecciones

La imperfeccion de la estructura sera el modo de pandeo asociado a la primera
carga crtica de la estructura escalada convenientemente teniendo en cuenta las toleran-
cias geometricas constructivas. La carga crtica del sistema perfecto se alcanza para unos
esfuerzos primarios en cada barra: Pcr = cr P0 ; Mcr = cr My0 , siendo {P0 , My0 } los esfuer-
zos de un calculo lineal.
En cada barra el modo de pandeo lleva asociados un {ccr}. {cr , cr } de-
nen los exponentes y unas constantes que habra que escalar {fcr , gcr , hcr , icr },
{acr , bcr , ccr , dcr , ecr , kcr , lcr } y que pueden obtenerse a partir de los movimientos de
los nudos en ejes locales haciendo uso de la [ Matriz 1cr ]. Se obtienen las siguientes imper-
fecciones a lo largo de la barra

Pcr
wimp = fcr ecr x + gcr ecr x + hcr + icr x; siendo cr =
EIy

vimp = acr ecr x + bcr ecr x + ccr ecr x + dcr ecr x + kcr + lcr x

(EIz 2cr Pcr ) 2

(EIz cr Pcr )
ximp = (acrecr x +bcr ecr x )+ (ccr ecr x +dcr ecr x )+mcr +ncr x
Mycr Mycr

siendo {cr , cr , cr , cr } las soluciones de la ecuacion caracterstica. Ecuacion

polinomica de cuarto grado en

(E 2 Iz Ia )4 + (EIz Pcr ro2 EIz GIt + Pcr EIa )2 + (Pcr2 r02 Mycr
2
Pcr GIt ) = 0

Obtencion de las fuerzas internas de un sistema imperfecto

Ecuaciones de equilibrio de un sistema imperfecto en una posicion ligeramente defor-
mada:
2
A) Asociado al axil: EA xu2 = 0. Tendremos dos grados de libertad {ud, uf }.
4 2 2
B) Efecto del axil sobre la exion primaria: EIy xw4 + P xw2 = P x
wimp
2

 
2w 3w w wimp
My = EIy 2 Qz = EIy 3 P +
x x x x

Solucion homogenea:

P
wh = f ex + gex + h + ix, siendo =
EIy
70 A. Aguero y J.R. Atienza

wdh

f

z
Q dh

f

y dh g My dh g
= =
wf h
Matriz 1 h

Qz f h
Matriz 2 h

y f h i My f h i
Solucion particular:
P
wp = (fcr ecr x + gcr ecr x ); {Qz dp, My dp, Qz f p, My f p}
EIy 2cr + P

C) Efecto del axil y la exion primaria sobre la exion lateral

4v 2 x 2v 2 ximp 2 vimp
EIz + M y + P = M y P
x4 x2 x2 x2 x2

2 x 2v 4 x 2
2 ximp 2 vimp
(GIt + P ro2 ) + M y + EIa = P r 0 M y
x2 x2 x4 x2 x2

   
2v 3v x ximp v vimp
Mz = EIz 2 Qy = EIz 3 My + P +
x x x x x x

   
2 x 3 x x x ximp v vimp
Bi = EIa 2 Mx = EIa 3 + GIt P ro2 + My +
x x x x x x x
Solucion homogenea:

vh = aex + bex + cex + dex + k + lx

(EIz 2 P ) (EIz 2 P ) x
xh = (aex + bex ) + (ce + dex ) + m + nx
My My

siendo {, , , } las soluciones de la ecuacion caracterstica. Ecuacion polinomica

de cuarto grado en

(E 2 Iz Ia )4 + (EIz P ro2 EIz GIt + P EIa )2 + (P 2 r02 My2 P GIt ) = 0

vdh a Qy dh a

x dh
b
Mx dh

b



















z dh

c


M z dh

c

dh d Bidh d
=
Matriz 1 k =
Matriz 2 k

vf h
Qy f h








x f h
l

M x f h
l






z f h
m M
z f h
m
f h n Bif h n
Pandeo elastico y diseno de estructuras espaciales 71

Solucion particular
vp1 = Amp1(acr ecr x + bcr ecr x ) + Amp2(ccrecr x + dcr ecr x )

xp1 = Amp3(acrecr x + bcr ecr x + Amp4(ccr ecr x + dcr ecr x )

(EIz 2cr Pcr ) 2

(EIz cr Pcr )
vp2 = Amp5 (acrecr x +bcr ecr x )+Amp6 (ccrecr x +dcr ecr x )
Mycr Mycr

(EIz 2cr Pcr ) 2

(EIz cr Pcr )
xp2 = Amp7 (acr ecr x +bcr ecr x )+Amp8 (ccr ecr x+dcr ecr x )
Mycr Mycr

    
EIz 2cr + P My Amp1 Amp5 P My
=
My EIa 2cr + (GIt + P ro2) Amp3 Amp7 My P ro2

 2
   
EIz cr +P My Amp2 Amp6 P My
2 =
My EIa cr + (GIt + P ro2) Amp4 Amp8 My P ro2

vp = vp1 + vp2 ; xp = xp1 + xp2 {Qy dp, Mx dp, Mz dp, Bidp, Qy df, Mx df, Mz df, Bidf }
resultando que las fuerzas internas estan compuestas por un termino constante debido a la
solucion particular mas un termino lineal en los desplazamientos
Fint = Kdh + F p

Ensamblaje y equilibrio
Tras expresar las fuerzas internas y desplazamientos en ejes globales se procede al
ensamblaje y se impone el equilibrio de fuerzas en los nudos igualando la resultante de las
fuerzas internas a las fuerzas exteriores.
 
FintG = (KG dhG + F pG ) = Fext
e e

resultando el sistema de ecuaciones

 
(KG dhG ) = Fext (F pG )
e e

Analisis de la estructura imperfecta

Para cada combinacion de acciones posibles, tras realizar un calculo lineal, obtendremos
unos esfuerzos primarios {P, My }. Con estos resolveremos el sistema indicado en el apartado
anterior, tras lo cual podemos obtener los esfuerzos en cada elemento y comprobar si son
admisibles.
72 A. Aguero y J.R. Atienza

CONCLUSIONES

Se ha presentado un programa de calculo de cargas crticas en estructuras espaciales con

secciones doblemente simetricas bajo la accion de unos esfuerzos primarios de axil y exion
segun el eje fuerte y una metodologa de obtencion los esfuerzos de diseno considerando
imperfecciones geometricas.
De los casos estudiados se puede concluir que las formulas del Eurocodigo 3 quedan
en algunos casos del lado de la inseguridad en una magnitud superior al 10 % (Figura 10),
dandose los mayores errores cuando k = kw. En las Figuras 13, 14 y 15 se representa el
error relativo de los casos estudiados.
Una limitacion tanto para el calculo de cargas crticas como de sistemas imperfectos
es que se ha considerado la inuencia de unos esfuerzos primarios constantes. En caso de
ser estos esfuerzos variables se puede o bien tomar un esfuerzo primario equivalente segun
el tipo de variacion o aumentar el numero de elementos.

Figura 13. Error {vd = 0, x d = 0, vf = 0, x f = 0}. {Mz d = 0, Bid = 0, Mz f = 0, Bif = 0}

Figura 14. Error {vd = 0, x d = 0, d = 0, vf = 0, x f = 0, f = 0}. {Mz d = 0, Mz f = 0}

Pandeo elastico y diseno de estructuras espaciales 73

Figura 15. Error. Caso biempotrada

REFERENCIAS

1 Eurocodigo 3, Proyecto de estructuras de acero, Parte 1-1: Reglas generales y reglas

para edicacion, UNE-ENV, (1993-1-1).

2 J.D. Renton, Stability of space frames by computer analisis, Journal of the Structural
Division, Proceedings of the American Society of Civil Engineers, Vol. 88, N ST4,
Agosto, (1962).

3 W.F. Chen y T. Atsuta, Theory of beam columns, Vol. 2: Space behaviour and
design, McGraw-Hill, New York, (1977).

4 R. Livesley y K. Chandler, Stability functions for structural frameworks, Manchester

University Press, (1956).

5 Don O. Brush y B.O. Almroth, Buckling of bars, plates, and shells, McGraw-Hill,
(1975).

6 J.M.T. Thompson y G.W. Hunt, Elastic instability phenomena, John Wiley and
Sons, (1984).

7 R. Kappus, Drillknicken zentrisch gedruckter Stabe mit oenem Prol im elastischen

Bereich, Luftfahrt-Forschung, (1937).
74 A. Aguero y J.R. Atienza

APENDICE I. NOTACION

u, v, w = desplazamientos segun eje x, y, z resp.

x , y , z = giros segun eje x, y, z resp.
= alabeo
= alabeo unitario
Mx , My , Mz = momentos segun eje x, y, z resp.
P = axil (positivo para compresion)
Qy , Qz = cortantes segun eje y, z resp.
Bi = bimomento
E = modulo de elasticidad
G = modulo de elasticidad transversal
Iy , Iz = momentos de inercia de eje y, z resp.
It = modulo de torsion
Ia o Iw = modulo de alabeo
W = modulo resistente
ud, vd, wd, xd , yd , zd , d = desplazamientos cara dorsal
uf, vf, wf, xf , yf , zf , f = desplazamientos cara frontal
P d, Qy d, Qz d, Mx d, My d, Mz d = esfuerzos cara dorsal
P f, Qy f, Qz f, Mx f, My f, Mz f = esfuerzos cara frontal
L = longitud de la barra
A = area de la seccion transversal
ro = radio de giro polar respecto del centro
de esfuerzos cortantes
y = lmite elastico

APENDICE II. CALCULO DE LA MATRIZ DE RIGIDEZ

> #
> #CAMPO DESPLAZAMIENTOS
> #
> v:=a*exp(x*C1)+b*exp(-x*C1)+c*exp(x*C2)+d*exp(-
x*C2)+k+l*x;dv:=diff(v,x):ddv:=diff(v,x,x):dddv:=diff(v,x,x,x):
> aa:=(-E*Iz*C1^2-P)*a/My:bb:=(-E*Iz*C1^2-P)*b/My:cc:=(-E*Iz*C2^2
-P)*c/My:dd:=(-E*Iz*C2^2-P)*d/My:
> g:=aa*exp(x*C1)+bb*exp(-x*C1)+cc*exp(x*C2)+dd*exp(-
x*C2)+m+n*x;dg:=diff(g,x):ddg:=diff(g,x,x):dddg:=diff(g,x,x,x):
>#
> #CALCULO ESFUERZOS
> #
> Mz:=E*Iz*ddv:Mx:=(G*It-P*ro^2)*dg-E*Ia*dddg-My*dv:Bi:=E*Ia*ddg:Qy:
=simplify(-E*Iz*dddv-My*dg-P*dv):
>x:=0:vd:=simplify(v):dvd:=simplify(dv):gd:=simplify(g):dgd:
=simplify(dg):
> x:=L:vf:=simplify(v):dvf:=simplify(dv):gf:=simplify(g):dgf:
=simplify(dg):
> x:=L:Mzf:=simplify(Mz):Mxf:=simplify(Mx):Bif:=simplify(Bi):Qyf:
=simplify(Qy):
> x:=0:Mzd:=simplify(-Mz);Mxd:=simplify(-Mx);Bid:=simplify(-Bi);Qyd:
=simplify(-Qy);
> with(linalg):
Pandeo elastico y diseno de estructuras espaciales 75

> K:=matrix(8,8):RIG:=matrix(8,8):
> BBB:= array(identity, 1..8,1..8):
> #
> # CALCULO MATRIZ 1
> #
> for II from 1 by 1 while II<9 do
>a:=BBB[1,II]:b:=BBB[2,II]:c:=BBB[3,II]:d:=BBB[4,II]:k:=BBB[5,II]:l:
=BBB[6,II]:m:=BBB[7,II]:n:=BBB[8,II]:
>K[1,II]:=vd:K[2,II]:=dvd:K[3,II]:=gd:K[4,II]:=dgd:K[5,II]:=vf:K[6,II]:
=dvf:K[7,II]:=gf:K[8,II]:=dgf:
> od:
> #
> # CALCULO INVERSA DE MATRIZ 1
> #
> SOL:=linsolve(K,BBB):
> #
> # CALCULO MATRIZ DE RIGIDEZ
> #
> for II from 1 by 1 while II<9 do
>a:=SOL[1,II]:b:=SOL[2,II]:c:=SOL[3,II]:d:=SOL[4,II]:k:=SOL[5,II]:l:
=SOL[6,II]:m:=SOL[7,II]:n:=SOL[8,II]:
>RIG[1,II]:=simplify(Qyd):RIG[2,II]:=simplify(Mxd):RIG[3,II]:
=simplify(Myd):RIG[4,II]:=simplify(Bid):RIG[5,II]:=simplify(Qyf):RIG[6,II]:
=simplify(Mxf):RIG[7,II]:=simplify(Myf):RIG[8,II]:=simplify(Bif):
> od:
> fortran(RIG,precision=double,filename=MATRIG,optimized):
> #
> # CALCULO C1 Y C2 . ECUACION CARACTERISTICA
> #
> A:=matrix([[E*Iz*alpha^2+P,My],[My,-G*It+P*ro^2+E*Ia*alpha^2]]):
> solve(det(A)=0,alpha);

APENDICE III. CALCULO DE LA MATRIZ DE ROTACION

># CALCULO DE MATRIZ DE ROTACION

> n1:=vector([DX,DY,DZ]):n1:=evalm(n1/XL(KB)):
> n3:=vector([-DX*DZ,-DZ*DY,DX^2+DY^2]):b:=(n3[1]^2+n3[2]^2+n3[3]^2)^(1/
> 2):n3[1]:=(n3[1]-e)/(b+e):n3[2]:=(n3[2])/(b+e):n3[3]:=(n3[3])/(b+e):
> n2:=vector([-DY,DX,0]):b:=(n2[1]^2+n2[2]^2+n2[3]^2)^(1/2);n2[1]:=(n2[1
> ])/(b+e)*DZ/abs(DZ+e):n2[3]:=0.0;n2[2]:=(n2[2]+e)/(b+e);
> T:=matrix([[n1[1],n2[1],n3[1]],[n1[2],n2[2],n3[2]],[n1[3],n2[3],n3[3]]
> ]);
># GIRO GAMA DE EJE X LOCAL
> T1:=matrix([[1,0,0],[0,cos(gama),-sin(gama)],[0,sin(gama),cos(gama)]]);
># MATRIZ DE ROTACION GENERAL
> Tr:=matrix(3,3,0):Tr:=evalm(T&*T1);
> fortran(Tr,precision=double,filename=rotacion,optimized);
76 A. Aguero y J.R. Atienza

APENDICE IV. MATRIZ DE RIGIDEZ Ia = 0

Se presentan las matrices de rigidez para el caso en el que la rigidez al alabeo sea
despreciable y para esfuerzos primarios de axil, exion, exion mas axil con objeto de poder
compararlas, siendo las funciones de estabilidad comunes
[ cos() sin()
s=
[2 + sin() + 2 cos()]

[ sin()]
c=
[sin() cos()]

t = (1 + c)s

FLEXION + AXIL

EIz
L3 MLy EIz
L2
t EIz
L3
My EIz
L2
t

Qy d
(GIt P ro2 )
L
ro2 )
vd

M d 0
My
(GIt P 0 x d
x
L L L



Mz d EIz
s EI
L2
z
t 0 EIz
sc d
= L L
z

Q f
y L3
EIz
MLy EI z
t vf
L2



x f
x
M f (GIt P ro2 )
0


Mz f L
EIz
z f
L
s

siendo


L EIz (GIt P ro2)(My2 + P GIt P 2 ro2)
=
EIz (GIt P ro2 )

P L2 2 [2 2 cos()]
= +
EIz 2 + sin() + 2 cos()]

FLEXION
EI
z
2t MLy EIz
t EI z
2t My EIz
t

Q d L3 L2 L3 L L2 vd
M d 0 d

y


GIt My
GI
x
t
0
x L L L


Mz d EIz
s EI z
t 0 EIz
sc z d
= L L2 L

Qy f EIz
2t MLy EI t
L2
z vf
M f
L3
f

x

GIt
0 x


L
Mz f EIz z f
L
s

siendo
My 
= Mcr = EIz GIt
Mcr L
Pandeo elastico y diseno de estructuras espaciales 77

AXIL

EIz 2t
L3 m 0 EIz
t EIz 2t
0 EIz
t
L2 L3 m L2

Q d
vd
0
y 2 2


M d
GIt P ro
0 My
GIt P ro
x d
x
L L L


Mz d
EIz
s EI
L2
z
t 0 EIz
sc d
= L
EIz 2t
L
z

y
Q f L3 m 0 EI z
t vf
L2

x f

xf GIt P ro2
0


M L
0Mz f EIz z f
L
s

siendo

P EIz 2
= Pcr =
Pcr L2

sin()
m=
2[cos() 1]