Tesis Completa

INSTITUTO POLITCNICO NACIONAL
ESCUELA SUPERIOR DE INGENIERA Y ARQUITECTURA SECCIN DE ESTUDIOS DE POSGRADO E INVESTIGACIN
"REPERCUSIN DE LAS IMPERFECCIONES GEOMTRICAS EN CILINDROS DE PARED DELGADA SOMETIDOS A CARGA AXIAL APLICADA A SILOS"
T
PARA
E
OBTENER
S
EL
GRADO DE :
MAESTRO EN CIENCIAS CON ESPECIALIDAD EN ESTRUCTURAS P R E S E N T A :
ING. LUIS MANUEL ESPINOSA NIETO DIRECTOR DR. HCTOR A. SNCHEZ SNCHEZ
MEXICO D.F
2005
Esta tesis esta dedicada en especial a la empresa Silos y Camiones, S.A. de C.V., cuyo Director es un hombre entusiasta y convencido de que el desarrollo Tecnolgico es la base del progreso, Gracias Ing. Romualdo Tellera Beltrn, por su apoyo para que lograra culminar con xito mi maestra
A mi esposa Luz Edith Por todo el apoyo brindado y su gran comprensin
A mi hijo Jess Emmanuel Por entender mis ausencias Y siempre sonrer
A mis padres Por sus motivaciones y enseanzas
A mis maestros y amigos En especial al Dr. Hctor Snchez
Dr. Hctor Snchez Gracias por haberme entendido desde Mi inicio en la maestra, por las enseanzas En las aulas, por haberme ayudado a generar Una cultura de investigacin, y sobre todo Gracias porque se que siempre podr Contar con usted para solucionar juntos Los retos tecnolgicos que se nos presenten. Gracias por ser mi amigo. Dr. Joel Garca Desde que nos conocimos Siempre te brindaste como amigo En el aula siempre preocupado por ensear Y en la vida eres una gran persona Dr. Jos Luis Urrutia En la vida vamos conociendo mucha gente En verdad me hubiese gustado conocerlo Mucho antes, pero nunca es tarde para Estar cerca de tan gran investigador Nos seguiremos viendo y Gracias por sus recomendaciones
Dr. Jorge Luis Alamilla Gracias por tu gusto a mi tesis Y por tus contribuciones Mucho xito en las actividades que emprendas. Mtro. Ramss Rodrguez Rocha Gracias por tus recomendaciones que Fueron muy acertadas, xito en tu prximo doctorado. Gracias Instituto Tecnolgico de Pachuca, por haberme dado Las bases que me ayudaron a tener la visin y lograr mis objetivos. Gracias Instituto Politcnico Nacional, por acogerme En tu seno y formarme como investigador
REPERCUSIN DE LAS IMPERFECCIONES GEOMTRICAS EN CILINDROS DE PARED DELGADA SOMETIDOS A CARGA AXIAL APLICADA A SILOS
Ciudad de Mxico Diciembre del 2005
Repercusin de las imperfecciones geomtricas en cilindros de pared delgada sometidos a carga axial aplicada a silos.
CONTENIDO
Pg
Lista de figuras Lista de tablas Lista de fotografas Simbologa Resumen Abstract Captulo I
1 1.1 1.1.1 1.1.2 1.2 1.3 1.4 1.5 Generalidades y planteamiento de la investigacin Marco histrico Cargas actuantes sobre un silo Esfuerzos permisibles Estado de la tcnica Objetivo y justificacin Alcance Planteamiento propuesto
IV VII VIII IX X XI
2 2 5 6 9 11 12 12
Captulo II
2 2.1 2.1.1 2.1.2 2.1.3 2.1.4 2.1.5 2.2 2.2.1 2.2.2 Marco terico Compresin axial en cilindros de pared delgada Pandeo por bifurcacin con cada brusca de rigidez Pandeo terico de estructuras sometidas a compresin axial Efectos de imperfecciones iniciales en la estabilidad de cascarones cilndricos Formas de inestabilidad considerando la influencia de las imperfecciones geomtricas Influencia de los esfuerzos residuales en la estabilidad de cascarones Criterios tpicos de diseo Compresin axial en cilindros de pared delgada, [Baker. 1972] Eurocdigo 3 [Eurocode 3. 2003b] 15 15 16 17 24 26 28 31 31 34
Captulo III
3 3.1 3.1.1 3.1.2 3.1.3 3.1.4 3.1.5 3.1.6 Modelado y pandeo terico Probeta controlada Caractersticas geomtricas y mecnicas de los modelos Pandeo terico de la probeta a estudiar Pandeo axisimtrico Pandeo no axisimtrico Tcnica de Rish y Kovalesky Modelacin controlada (II) 42 42 42 44 47 49 50 50
Captulo IV
Pg 4 4.1 4.2 4.2.1 4.2.2 4.2.3 4.2.4 4.2.5 4.3 4.4 4.4.1 4.4.2 4.4.3 4.4.4 4.4.5 4.5 Experimentacin Probeta estudiada Medicin de los especimenes Medicin del espcimen 1 Medicin del espcimen 2 Medicin del espcimen 3 Medicin del especimen 4 Medicin del espcimen 5 Modelacin numrica Descripcin de las pruebas Espcimen 1 Espcimen 2 Espcimen 3 Espcimen 4 Espcimen 5 Comentarios 79 79 81 83 85 86 87 89 90 91 91 96 100 103 106 110 114 115 118 119 120 124 126 128
Discusin de resultados Conclusiones Recomendaciones Trabajos futuros Referencias Anexos 1 Medicin de imperfecciones Anexos 2 Matriz de transformacin Anexos 3 Video y archivos electrnicos
(III)
FIGURAS
Captulo I
Figura I-1 Clasificacin de silos Figura I-2 Apoyos tpicos de silos metlicos Figura I-3 Partes de un silo Figura I-4 Seccin de un faldn con perfiles tipo zeta (vista en planta) Figura I-5 La unin de dos o ms rolos forman un faldn o cuerpo del silo Pg 2 3 4 5 7
Captulo II
Figura II-1. Pandeo por bifurcacin con cada brusca de rigidez Figura II-2. Modos de pandeo Figura II-3. Modelo de panel curvo Figura II-4. Comportamiento del pandeo Figura II-5 Modos de pandeo tpicos para cilindros perfectos cargados axialmente Figura II-6 Esfuerzos de pandeo experimentales de cilindros isotrpicos comprimidos axialmente Figura II-7. Punto de bifurcacin y punto limite mediante un anlisis asinttico Figura II-8. Punto de bifurcacin y punto limite mediante un anlisis no-lineal Figura II-9 Medicin de imperfecciones, solo carga axial Figura II-10 Medicin de imperfecciones, cuando se incluyen efectos de cortante Figura II-11 Medicin imperfecciones, en soldaduras Figura II-12 a) Compresin uniforme , b) Compresin no uniforme 16 18 19 21 23 25 27 28 36 37 37 40
Captulo III
Figura III-1 Elemento Shell 63 Figura III-2 Punto de bifurcacin y punto limite Figura III-3 Cilindro perfecto, primer modo de inestabilidad Figura III-4 Curva de estabilidad neutra para diversas relaciones R/t Figura III-5 Curva de estabilidad neutra para valores de n y m Figura III-6. Esquema de probeta virtual Figura III-7 Comportamiento de la carga crtica normalizada en relacin con el factor de fabricacin Figura III-8. Comportamiento de la carga crtica normalizada en relacin con las imperfecciones Figura III-9 Amplificacin de la figura III-8 Figura III-10. Modos de circunferenciales, probeta 1 Figura III-11. Concentraciones de esfuerzos, probeta 1 Figura III-12 Modos de pandeo, probeta 1 Figura III-13 Distribucin de esfuerzos para la carga de bifurcacin de la probeta 1 Figura III-14 Anlisis no lineal y punto de bifurcacin para la probeta 1 Figura III-15 Modos de pandeo, probeta 2 Figura III-16 Distribucin de esfuerzos para la carga de bifurcacin de la probeta 2 Figura III-17 Anlisis no lineal y punto de bifurcacin para la probeta 2 Figura III-18 Modos de pandeo, probeta 3 Figura III-19 Distribucin de esfuerzos para la carga de bifurcacin de la probeta 3 (IV) 45 46 47 48 50 52 52 53 54 56 56 58 59 59 60 61 61 62 63
Pg Figura III-20 Anlisis no lineal y punto de bifurcacin para la probeta 3 Figura III-21 Modos de pandeo, probeta 4 Figura III-22 Distribucin de esfuerzos para la carga de bifurcacin de la probeta 4 Figura III-23 Anlisis no lineal y punto de bifurcacin para la probeta 4 Figura III-24 Modos de pandeo, probeta 5 Figura III-25 Distribucin de esfuerzos para la carga de bifurcacin de la probeta 5 Figura III-26 Anlisis no lineal y punto de bifurcacin para la probeta 5 Figura III-27 Modos de pandeo, probeta 6 Figura III-28 Distribucin de esfuerzos para la carga de bifurcacin de la probeta 6 Figura III-29 Anlisis no lineal y punto de bifurcacin para la probeta 6 Figura III-30 Modos de pandeo, probeta 7 Figura III-31 Distribucin de esfuerzos para la carga de bifurcacin de la probeta 7 Figura III-32 Anlisis no lineal y punto de bifurcacin para la probeta 7 Figura III-33 Modos de pandeo, probeta 8 Figura III-34 Distribucin de esfuerzos para la carga de bifurcacin de la probeta 8 Figura III-35 Anlisis no lineal y punto de bifurcacin para la probeta 8 Figura III-36 Modos de pandeo, probeta 9 Figura III-37 Distribucin de esfuerzos para la carga de bifurcacin de la probeta 9 Figura III-38 Anlisis no lineal y punto de bifurcacin para la probeta 9 Figura III-39 Modos de pandeo, probeta 10 Figura III-40 Distribucin de esfuerzos para la carga de bifurcacin de la probeta 10 Figura III-41 Anlisis no lineal y punto de bifurcacin para la probeta 10 63 64 65 65 66 67 67 68 69 69 70 71 71 72 73 73 74 75 75 76 77 77
Captulo IV
Figura IV-1 Espcimen estudiado Figura IV-2 Medicin de las imperfecciones geomtricas, espcimen 1 Figura IV-3 Mxima imperfeccin generada en el espcimen 1, paralelo 125 Figura IV-4 Mxima imperfeccin generada en el espcimen 1, meridiano 290 Figura IV-5 Medicin de las imperfecciones geomtricas, espcimen 2 Figura IV-6 Mxima imperfeccin generada en el espcimen 2, paralelo 250 Figura IV-7 Mxima imperfeccin generada en el espcimen 2, meridiano 180 Figura IV-8 Medicin de las imperfecciones geomtricas, espcimen 3 Figura IV-9 Mxima imperfeccin generada en el espcimen 3, paralelo 175 Figura IV-10 Mxima imperfeccin generada en el espcimen 3, meridiano 20 Figura IV-11 Medicin de las imperfecciones geomtricas, espcimen 4 Figura IV-12 Mxima imperfeccin generada en el espcimen 4, paralelo 50 Figura IV-13 Mxima imperfeccin generada en el espcimen 3, meridiano 290 Figura IV-14 Medicin de las imperfecciones geomtricas, espcimen 5 Figura IV-15 Mxima imperfeccin generada en el espcimen 5, paralelo 75 Figura IV-16 Mxima imperfeccin generada en el espcimen 5, meridiano 60 Figura IV-17 Correcta aplicacin de la carga Figura IV-18 Grfico que muestra la sensibilidad de no aplicar la carga en todos los nodos del modelo numrico, se normaliza respecto al dato de carga crtica numrico Figura IV-19 Modos de pandeo, espcimen 1 Figura IV-20 Distribucin de esfuerzos para la carga de bifurcacin del espcimen 1 Figura IV-21 Anlisis no lineal y punto de bifurcacin para el espcimen 1 Figura IV-22 Modos de pandeo, espcimen 2 Figura IV-23 Distribucin de esfuerzos para la carga de bifurcacin del espcimen 2 (V) 80 84 84 84 85 85 86 86 87 87 88 88 88 89 89 90 93 94 95 95 96 99 99
Pg Figura IV-24 Anlisis no lineal y punto de bifurcacin para el espcimen 2 Figura IV-25 Grfico de las cargas de inestabilidad del espcimen 3 Figura IV-26 Modos de pandeo, espcimen 3 Figura IV-27 Distribucin de esfuerzos para la carga de bifurcacin del espcimen 3 Figura IV-28 Anlisis no lineal y punto de bifurcacin para el espcimen 3 Figura IV-29 Modos de pandeo, espcimen 4 Figura IV-30 Distribucin de esfuerzos para la carga de bifurcacin del espcimen 4 Figura IV-31 Anlisis no lineal y punto de bifurcacin para el espcimen 4 Figura IV-32 Grfico de las cargas de inestabilidad de la probeta 5 Figura IV-33 Modos de pandeo, espcimen 5 Figura IV-34 Distribucin de esfuerzos para la carga de bifurcacin del espcimen 5 Figura IV-35 Anlisis no lineal y punto de bifurcacin para el espcimen 5 Figura IV-36 Comparacin de las cargas de inestabilidad de los 5 especimenes 100 101 102 103 103 105 105 106 108 109 109 110 112
(VI)
TABLAS
Captulo II
Tabla II-1 Clasificacin de cascarones cilndricos con base a la calidad de fabricacin Tabla II-2 Clasificacin de silos segn su clase Tabla II-3 Valores del parmetro de imperfeccin Uomax Pg 34 35 38
Captulo III
Tabla III-1 Valores de estabilidad neutra para diversas relaciones R/t Tabla III-2 Valores de estabilidad neutra para valores de n y m Tabla III-3 Generacin de imperfecciones y cargas encontradas 48 50 55
Captulo IV
Tabla IV-1 Resultados comparativos de varias metodologas Tabla IV-2 Cargas criticas obtenidas vs. carga clsica esperada Tabla IV-3 Carga critica de los especimenes usando el Eurocdigo 110 111 113
(VII)
FOTOGRAFAS
Captulo I
Fotografa I-1 Falla en un silo por pandeo local Pg 10
Captulo II
Fotografa II-1 Apariencia tpica de pandeo a compresin axial Fotografa II-2 Rolo o virola rolada (Cortesa Silos y Camiones, S.A. de C.V.) Fotografa II-3 Esfuerzos residuales por el proceso de soldadura Fotografa II-4 Hojas empaquetadas (Cortesa Silos y Camiones, S.A. de C.V.) Fotografa II-5 Silos atornillados (Cortesa Silos y Camiones, S.A. de C.V.) 26 29 30 30 31
Capitulo IV
Fotografa IV-1 Especimenes fabricados Fotografa IV-2 Anillos de rigidez Fotografa IV-3 Mesa de nivelacin Fotografa IV-4 Espcimen de prueba sobre la mesa de nivelacin Fotografa IV-5 Sistema de medicin Fotografa IV-6 Uso del indicador de cartula Fotografa IV-7 Espcimen 1 montado sobre la maquina universal Fotografa IV-8 Espcimen 1, modo de pandeo Fotografa IV-9 Deformacin permanente Fotografa IV-10 Espcimen 2 montado sobre la maquina universal Fotografa IV-11 Espcimen 2, modo de pandeo Fotografa IV-12 Espcimen 3, modo de pandeo Fotografa IV-13 Espcimen 4, modo de pandeo Fotografa IV-14 Espcimen 5, modo de pandeo 80 81 81 82 82 83 91 92 92 97 98 100 104 107
(VIII)
SIMBOLOGA
B D D E E Fs fy H, h Is K1, K2 Kc L lg lg lgw lgx m N n Ncr P P Pclsica Pcr Q R, r t Uo Uow Uox V w wok X, Y, Z z o V cr Distancia entre atiezadores Dimetro del cilindro o silo, Rigidez del cilindro a flexin Desplazamientos, distancia Modulo de elasticidad Excentricidad, Exponencial Esfuerzo admisible a compresin Esfuerzo de fluencia Altura de los especimenes o cilindros, distancia en general Segundo momento de rea de la combinacin atiezador cilindro Constantes para la determinacin del esfuerzo crtico a pandeo en cilindros c/atiezadores Coeficiente de pandeo Longitud, altura Carga critica normalizada Longitud del calibrador para medir imperfecciones Longitud del calibrador para medir imperfecciones en direccin circunferencial Longitud del calibrador para medir imperfecciones en soldaduras Longitud del calibrador para medir imperfecciones en direccin meridional Modos axiales Carga en general Modos circunferenciales Carga crtica Carga aplicada Carga de bifurcacin, Carga de pandeo experimental Carga de pandeo clsica o de Euler Carga crtica de pandeo Parmetro de calidad Radio de los cilindros, radio en general Espesor de los cilindros, espesor en general Profundidad de la imperfeccin normalizada en direccin circunferencial Profundidad de la imperfeccin normalizada para soldaduras Profundidad de la imperfeccin normalizada en direccin meridional Energa potencial Desplazamiento en general Amplitud de la imperfeccin Ejes rectangulares Parmetro de Batdorf Desplazamientos Factor de reduccin, para cilindros no presurizados Deformaciones Variacin de potencial Imperfeccin geomtrica normalizada Factor de correccin Factor de correccin por plasticidad Mdulo de Poisson, 0.3 para acero Esfuerzo crtico de pandeo a compresin axial Parmetro para esfuerzo no uniforme (IX)
Resumen
Debido a la constante necesidad de las industrias (plstico, petroqumica, construccin, etc) por ser ms competitivas, se han revolucionado los procesos para el manejo y almacenamiento de los materiales a granel. Por lo tanto se ha optado por almacenar este tipo de productos en silos; sin embargo, muchos son los silos que son diseados deficientemente y acompaados por lo tanto de variadas fallas estructurales, siendo la ms comn la falla por inestabilidad estructural (pandeo), efecto que se ve incrementado por la presencia de imperfecciones geomtricas iniciales. Esto es resultado del desconocimiento de muchos diseadores en la aplicacin de los modelos matemticos aplicables al diseo de silos, y sobre todo, al poco control que se tiene de las imperfecciones geomtricas existentes en este tipo de estructuras durante su construccin. En esta tesis se presenta brevemente la teora de estabilidad de cascarones cilndricos aplicada al comportamiento y diseo de silos metlicos, adems de pruebas experimentales, que comprenden la fabricacin y ensayo de cinco especmenes, a los cuales se les miden sus imperfecciones geomtricas iniciales y se les aplica una carga axial de compresin hasta la obtencin de su carga crtica de pandeo; para posteriormente comparar stos resultados con los obtenidos mediante tcnicas numricas, como es el mtodo del elemento finito, y la normatividad existente en silos. Este trabajo de tesis contribuye a proponer una metodologa prctica para el anlisis y diseo de silos metlicos; sustentados en pruebas experimentales que son verificadas mediante el modelado numrico de los especimenes estudiados en este trabajo, considerando la influencia de las imperfecciones geomtricas iniciales y basndose en la teora clsica sobre este tema.
(X)
Abstract
Due to a constant need of industries (plastics, petrochemical, construction, etc) to be more competitive, they have revolutionized handling and storage processes of bulk materials. Thus, it is convenient to store these materials in silos; however, many of these are deficiently designed therefore with potential structural failures, being buckling the most common structural type of failure. Geometric imperfections might increase the bucking effect. This problem is caused as a result of inadequate mathematical models applied for silos design moreover due to non-control of geometric imperfections during construction.
This work applies stability of cylindrical shells theory to design of steel silos. In addition, experimental tests are presented including the geometric imperfections effect on five specimens. These are studied under compression axial load to obtain its buckling critical load. Results are compared against those from numerical methods such as finite element models (FEM) as well as code guideline provisions for silos design. This thesis work contributes with a practical methodology for analysis and design of steel silos.
(XI)
Captulo I
Generalidades y planteamiento de la investigacin
1 1.1 1.1.1 1.1.2 1.2 1.3 1.4 Generalidades y planteamiento de la investigacin .......................................... 2 Marco histrico .................................................................................................. 2 Cargas actuantes sobre un silo ......................................................................... 5 Esfuerzos permisibles ....................................................................................... 6 Estado de la tcnica .......................................................................................... 9 Objetivo y justificacin..................................................................................... 11 Planteamiento propuesto ................................................................................ 12
Resumen
El presente captulo es una introduccin a la tesis desarrollada. Se inicia con el marco histrico, se definen las justificaciones y el objetivo que se pretende alcanzar en la realizacin de esta investigacin. Se presenta el estado del arte sobre este tema, y se comentan algunas metas a lograr al finalizar este trabajo de investigacin.
1 Generalidades y planteamiento de la investigacin

1.1 Marco histrico
Debido a la gran ventaja que representa el almacenamiento de material a granel en contenedores, actualmente el uso de silos metlicos se ha incrementado de manera considerable. Comnmente se acostumbraba a tener en grandes patios toda la materia prima (material a granel) ya sea en sacos o en la misma intemperie, esto ocasionaba que se necesitaran grandes reas de almacenamiento, lo cual se ha solucionado en gran parte con la tecnologa de ensilar el producto en grandes tanques cilndricos.
Como parte de la reduccin de costos en la recepcin de materia prima, un silo debe ser seguro y capaz de resistir las condiciones de carga a las que estar sujeto en su vida til. El tema silos, a pesar de ser un tema difundido y utilizado en la ingeniera civil, presenta un desconocimiento muy grande, prueba de ello son las deformaciones, grietas, fallas, hundimientos y explosiones aparecidas en muchas instalaciones en prcticamente todos los pases industrializados. Rara es la instalacin en la que se manejan productos granulares y pulverulentos, que no haya tenido problemas en sus silos, ya sean estructurales o de flujo de vaciado, con interrupciones, bvedas, etc. Se ha estudiado la bibliografa existente sobre el tema de silos desde su inicio en el ao 1882 hasta nuestros das, pudiendo de esta forma describir las teoras de los principales especialistas mundiales. En la actualidad hay una gran variedad de silos, pudindose clasificar como se muestra en la figura I-1.
SILOS
CONCRETO
ALUMINIO
A. INOXIDABLE
A. CARBN
SOLDADOS MODULARES
Figura I-1 Clasificacin de silos. -2-
Es obvio notar la ventaja de un silo metlico contra uno de concreto, siendo sta la posibilidad de reutilizacin del silo metlico en otro lugar. Un silo metlico se mantiene erecto por su sistema de soporte con el suelo y existen tres formas de mantener en pie un silo, siendo bsicamente las siguientes (figura I-2):
Mediante columnas, se colocan en el permetro del silo y son de perfiles metlicos comerciales. Mediante silletas, muy utilizadas para sentar o elevar un silo sobre estructuras metlicas o de concreto. Con faldn, esta es una continuacin del cuerpo del silo y cuyo xito esta en apoyarse uniformemente sobre la cimentacin.
Patas
Silletas
Faldn
Figura I-2 Apoyos tpicos de silos metlicos.
El faldn es un elemento estructural bastante utilizado en el diseo de tanques y silos, puesto que es el elemento que mantiene de pie al silo, y debido a que existe poca reglamentacin sobre el diseo de faldones y/o silos, es muy importante conocer su comportamiento bajo la accin de cargas axiales de compresin y/o de flexin. De manera general un faldn, tal y como se ve en la figura I-3, es un cilindro que apoya sobre una base de concreto (cimentacin) y que sirve de soporte al cuerpo del silo (donde se almacena material a granel). La tolva de un silo es el medio para fluidizar el -3-
material y trasmite el peso del material almacenado al faldn del silo cuando el material est en estado de no flujo (estacionario).
Figura I-3 Partes de un silo
El problema consiste en disear un faldn que soporte las descargas actuantes en el silo, generadas por las diferentes combinaciones de carga. Normalmente los silos con faldn se fabrican con lmina o placa de acero (inoxidable o al carbn), y/o en aluminio. Los espesores empleados en faldones de silos son variables, pudiendo ser desde lmina de 2.66 mm (calibre 12) hasta placa de 8 mm (5/16") o inclusive ms si el anlisis y el diseo lo requieren. Para reforzar an ms el faldn de un silo se pueden colocar elementos estructurales (columnas) en el permetro del silo, haciendo con esto un sistema estructural con atiezadores axiales, ms rgido y de mayor resistencia, ver figura I-4.
-4-
Faldon
b.- distancia entre atiezadores

Atiezadores
Figura I-4 Seccin de un faldn con perfiles tipo zeta (vista en planta)
Una falla comn en faldones metlicos, es el fenmeno de pandeo que da inicio a fallas locales y que pueden ser el principio de que un silo se colapse y se desplome. Debido a este efecto se hace necesario el estudio de estos elementos estructurales y sobre todo el uso de factores de seguridad bastante elevados (entre 2 y 3). El factor de seguridad empleado para faldones de silos metlicos es amplio debido a que es difcil controlar la redondez y verticalidad, adicionalmente a esto, en el sitio de instalacin de los silos se acostumbra hacer perforaciones al faldn para el paso de tuberas y otros elementos. Por esto, para llegar a un buen diseo del faldn (placa curva de pared delgada) se deben usar factores de seguridad apropiados de modo que exista una correcta seguridad estructural y que los costos de los materiales empleados sean adecuados.
1.1.1 Cargas actuantes sobre un silo

Un silo est sometido a mltiples condiciones de carga, debidas a presiones internas dentro del silo, ocasionadas por fenmenos de fluidizacin, por efectos de la forma en como se llene el silo (mecnica o neumticamente), asentamientos en la cimentacin, o bien por efectos accidentales, tales como efectos ssmicos o la fuerza del viento. Las principales condiciones de carga que se deben estudiar en un silo segn la experiencia que se ha tenido en su diseo y que actualmente es usada en algunas empresas de manufactura de silos es la siguiente: -5-
Presiones de llenado del silo Presiones de vaciado del silo Sobre presiones de vaciado del silo Peso del silo vaco ms peso nominal del material almacenado, ms carga viva esperada Peso del silo vaco, ms carga de viento, ms carga viva esperada Peso del silo vaco, ms capacidad de trabajo del silo, y ms carga ssmica
1.1.2 Esfuerzos permisibles

Para el estudio de las diferentes combinaciones de carga son vlidas todas las teoras de esfuerzos existentes aplicables a la combinacin pertinente, sin embargo para el anlisis ssmico, viento y carga gravitacional es importante conocer un fenmeno conocido como pandeo, el cual ocurre con frecuencia en elementos de pared delgada. A continuacin se describen dos mtodos prcticos empleados para determinar el esfuerzo terico resistente a pandeo de un cilindro sometido a carga axial.
1.1.2.1 Esfuerzo
de
pandeo
para
cascarones
cilndricos
sin
atiezadores longitudinales
La AWWA (American Water Works Association) ha establecido que para esfuerzos de compresin en anillos circulares de lminas o placas bajo cargas de viento o sismo combinadas con carga muerta, el esfuerzo admisible de compresin debe determinarse mediante la siguiente frmula:
lb t 2 t 2 f s = 15000 100 2 100 15000 2 in R 3 R 3
(I-1)
donde: fs = esfuerzo admisible de compresin, en libras por pulgada cuadrada. t = espesor de la lmina o placa usada, pulgadas. r = radio medio del rolo o anillo, pulgadas. -6-
Es importante mencionar que un cilindro se puede constituir de varias secciones cilndricas soldadas o atornilladas entre si, a cada seccin se le da el nombre de rolo o anillo. Es decir que una seccin cilndrica (rolo) puede dar como resultado un faldn o el cuerpo de almacenamiento del silo (figura I-5).
Figura I-5 La unin de dos o ms rolos forman un faldn o cuerpo del silo
1.1.2.2 Esfuerzo
de
pandeo
para
cascarones
cilndricos
con
atiezadores longitudinales
La siguiente expresin para el esfuerzo de pandeo para elementos cilndricos atiesados longitudinalmente y sometidos a compresin axial, fue desarrollada por Shang [Gaylord. 1984]:
f cr =
donde:
12 1 2 (b t )
k1 k 2 2 E
(I-2)
2
-7-
0.24b 2 b = distancia circunferencial entre atiezadores. rt (500 ) k 2 = 1 1.75 si p 500 k 2 = 1 si 500 k 1 k1 = 1.3 +
La distancia circunferencial entre atiezadores b se determina del cociente del permetro del cilindro entre el nmero de atiezadores determinado.
r , t = radio y espesor del cilindro, respectiva mente. Si el modulo de Poisson es 0.3 para acero, = 0.3 I 10.9 I s = 3 s = 2 bt 12 1 bt 3
Donde Is es el momento de inercia del atiesador, incluyendo el ancho efectivo del cilindro. El ancho efectivo de un atiesador en combinacin con el cilindro se puede calcular de manera aproximada en funcin del radio y espesor del cilindro como: bef = 0.76 rt
Para acero, se pueden sustituir el mdulo de elasticidad y el coeficiente de Poisson
E = 30000
f cr =
klb in 2
y = 0.3 reducindola a:
27000k1 k 2
(b t )2
Sustituyendo el valor de k1, el esfuerzo crtico fcr se simplifica en la ecuacin I-3 [lb/in2]:
35000 6480 f cr = k 2 (b t )2 + r t
(I-3)
-8-
1.1.2.3 Factor de seguridad

El factor de seguridad depende mucho del conocimiento que se tenga acerca del comportamiento de silos metlicos, ya que normalmente en la manufactura de silos es prcticamente imposible obtener un cilindro perfecto, debido a las imperfecciones geomtricas. Gary A. Chubb and Keith E. McGuire [ Chubb. 1996] en la revista Powder and Bulk Engineering, hacen mencin que dentro de la zona de almacenaje del producto, factores de seguridad de 2 y 3 son adecuados; esto debido a que la presin interna que produce el material almacenado estabiliza el efecto de pandeo dentro de esta zona, sin embargo si el silo va a ser soportado por un faldn con una tolva suspendida en la descarga (cono de descarga), con puerta, ventanas y otros ductos para paso de tubera y sabiendo que esta zona del faldn carece de presin interna que ayude a rigidizarlo, el factor de seguridad se ve incrementado y el uso de factores de seguridad de 3 y 4 se justifican.
1.2 Estado de la tcnica

En la actualidad existen diversas investigaciones realizadas en el extranjero acerca del tema de imperfecciones geomtricas en silos (cilindros de pared delgada), tal es el caso de las investigaciones de J. Michael Rotter y su grupo de investigadores, de J, Mark F.G. Holst, Chris R. Calladine, u otras investigaciones anteriores como Koiter, 1963; Yamaki, 1984. Todas estas investigaciones han tenido como conclusin importante que los cilindros de pared delgada sometidos a carga axial, son muy susceptibles de falla por efecto de inestabilidad, debido a las imperfecciones geomtricas generadas por los procesos de fabricacin con los que se manufacturan este tipo de elementos. A pesar de esto pocas son las normas o recomendaciones existentes de divulgacin internacional, tal es el caso de las normas DIN 18800.1990, ECCS.1988, ESDU.1974, [Baker. 1972], CRCJ.1971, NASA.1968. Probablemente uno de los mejores cdigos para placas metlicas de revolucin sean el ECCS, 1988 [Eurocode 3. 1988] y las normas alemanas DIN 18800.
-9-
El avance que se ha tenido en la realizacin de normas y reglamentos para este tipo de estructuras ha sido gracias a la realizacin de programas de computacin poderosos y exclusivos para este tipo de estudios.
Las estructuras denominadas silos son un caso especial para el cual no existen mayores normas en el mundo, los mejores estndares para silos comienzan con el Japanese Aluminum Silos Standard y el British Materials Handling Board / BSI Standard, Eurocode 3: Design of steel structures, part 4.1 Silos [Eurocode 3, 2003b]. La gua australiana AISC y textos como [Gaylord. 1984], Martens.1988, Hampe. 1991, proporcionan las mejores recomendaciones de diseo. Sin embargo, poco se ha estudiado para conocer los efectos generados por esfuerzos residuales en adicin con las imperfecciones geomtricas en cilindros de pared delgada. Sin duda alguna la falta de reglamentacin ocasiona que los fabricantes de silos construyan en muchas ocasiones con base a su experiencia o bien usando factores de seguridad elevados o incluso escasos, provocando con esto problemas estructurales de inestabilidad que pueden provocar el colapso de este tipo de estructuras, ver fotografa I-1.
Fotografa I-1 En esta fotografa se muestra una falla en un silo por pandeo local. (Cortesa Silos y Camiones, Pachuca Hgo. Mxico).
Por tanto, con base a lo anteriormente expuesto se propone en esta tesis evaluar los valores recomendados por el Eurocode 3 y compararlos con lo establecido por [Baker. - 10 -
1986]; para posteriormente realizar pruebas experimentales de especimenes metlicos y evaluar sus imperfecciones geomtricas para conocer su carga crtica y con esto tener un criterio de seleccin del factor de seguridad apropiado para el diseo de faldones en silos metlicos. En la actualidad en Mxico el diseo de faldones de silos metlicos se realiza mediante una reduccin emprica de la carga crtica terica, es decir reducindola hasta en un 30% del valor terico y no se consideran explcitamente la contribucin de las imperfecciones geomtricas en estos elementos estructurales.
1.3 Objetivo y justificacin

Se ha observado que an conociendo la teora de inestabilidad bsica de silos (generalmente usando comportamientos elsticos lineales) se generan fallas locales en diferentes partes de ste, como en el faldn y en el cuerpo de almacenamiento. Esto debido a que los factores de seguridad usados en el diseo de silos son variables y se dejan al criterio del diseador. Sin embargo, no se consideran los efectos causados por imperfecciones geomtricas en los silos y mucho menos los esfuerzos residuales causados por los procesos de fabricacin de este tipo de elementos estructurales.
El objetivo de esta investigacin se enfoca en la obtencin de resultados que ayudarn a tener un marco de referencia con respecto a la teora bsica, y poder sugerir valores de factor de seguridad y compararlos con los usados en la prctica profesional. La justificacin de la presente investigacin radica en que en Mxico hay poca o nula informacin que ayuden a determinar el factor de seguridad adecuado en el diseo de silos metlicos.
Para lograr el objetivo planteado se necesita estudiar como influyen las imperfecciones geomtricas (ocasionados por los procesos de fabricacin, tales como, rolado y al calor generado por los procesos de soldadura) en cilindros de pared delgada, conocidas como cascarones.
- 11 -
1.4 Alcance
Este trabajo busca estimar como afectan las imperfecciones geomtricas en el comportamiento de faldones en silos metlicos y dar un criterio de decisin sobre el uso adecuado del factor de seguridad apropiado en el diseo de silos metlicos.
Esta investigacin ayudar a que en el mbito nacional, en especial en las empresas dedicadas a la manufactura de silos, se comience con la realizacin de normas y reglamentos internos para el diseo y construccin de silos metlicos, buscando en un futuro que estas recomendaciones sean de legislacin nacional.
1.5 Planteamiento propuesto

En la actualidad se tiene contacto con J. Michael Rotter el cual ha compartido informacin sobre este tema, y con esto se ha logrado plantear el siguiente procedimiento para realizar la investigacin.
Seleccionar un cilindro de dimetro "D" , espesor "t" y longitud "L" Realizar un anlisis de inestabilidad, proponiendo imperfecciones controladas en modelos de computadora con la tcnica del elemento finito. Las imperfecciones sern establecidas con base a lo estipulado por el Eurocode 3. La razn de usar el Eurocdigo 3, es debido a que tiene clasificados los grados de imperfeccin en funcin del tipo de fabricacin que se est realizando y es uno de los pocos (quizs el nico) reglamentos que est considerando las imperfecciones geomtricas como parte de la manufactura y proporciona valores que estn en funcin del dimetro y el espesor del cascarn cilndrico. Comparar los resultados tericos obtenidos con la teora clsica de inestabilidad de cilindros sometidos a compresin axial, y con la teora de Baker. [Baker. 1972] Fabricacin de probetas de dimetro "D" , espesor "t" y longitud "L" Medicin de imperfecciones geomtricas en especimenes (probetas).
- 12 -
Realizar modelos de anlisis por computadora de las probetas fabricadas, introduciendo las imperfecciones geomtricas reales y realizar anlisis de estabilidad por pandeo. Realizar pruebas a compresin de las probetas en una maquina universal para encontrar la carga crtica real. Comparar los resultados y establecer conclusiones
Debido a la complejidad de realizar operaciones analticas del fenmeno de pandeo, se har uso de la herramienta de anlisis numrico basada en el mtodo del elemento finito, mediante el software ANSYS, usando elementos tipo shell, (quadrilteros o triangulares) que aceptan rotaciones y traslaciones en los tres ejes conocidos (X, Y y Z). Es importante seleccionar un elemento que permita los seis grados de libertad mencionados, debido a que no se puede conocer inicialmente en que direccin se generan las deformaciones en el cascarn, es por ello que no es conveniente la realizacin de un modelado tipo axisimetrico, ya que no reflejara la realidad de las imperfecciones en toda la superficie del modelo.
Se usa el software Ansys por tener varias tcnicas de solucin apropiadas para el fenmeno de pandeo, por contar con manuales, y recomendaciones personales por parte de la compaa que distribuye este software, adems de contar con varios cursos de capacitacin en el uso y aplicacin de este software.
- 13 -
Captulo II
Marco terico
2 2.1 2.1.1 2.1.2 2.1.3 2.1.4 2.1.5 2.2 2.2.1 2.2.2 Marco terico .................................................................................................. 15 Compresin axial en cilindros de pared delgada............................................. 15 Pandeo por bifurcacin con cada brusca de rigidez....................................... 16 Pandeo terico de estructuras sometidas a compresin axial. ....................... 17 Efecto de las imperfecciones en la estabilidad de cascarones cilndricos....... 24 Formas de inestabilidad considerando las imperfecciones geomtricas......... 26 Influencia de los esfuerzos residuales en la estabilidad de cascarones.......... 28 Criterios tpicos de diseo ............................................................................... 31 Compresin axial en cilindros de pared delgada, [Baker. 1972] ..................... 31 Eurocdigo 3 [Eurocode 3. 2003b] .................................................................. 33
Resumen
El objetivo particular de este captulo es presentar la teora clsica que existe respecto al comportamiento de cascarones cilndricos sometidos a carga axial. Se presentan las ecuaciones tericas que rigen el comportamiento de los cascarones cilndricos. Se muestra el proceso bsico de fabricacin de silos y se explica la procedencia de las imperfecciones geomtricas consideradas en este tipo de estructuras.
2 Marco terico 2.1 Compresin axial en cilindros de pared delgada

Miembros estructurales sujetos a compresin son susceptibles a fallar por inestabilidad estructural, fenmeno conocido como pandeo elastico de Euler. estructuras axisimtricas de pared delgada depende La estabilidad de de las
principalmente
caractersticas geomtricas de la estructura, de las caractersticas mecnicas y la no linealidad de los materiales, de las condiciones de frontera, de las imperfecciones geomtricas iniciales, alteraciones en el material y quizs de los esfuerzos residuales generados por los procesos de fabricacin para el formado de cilindros.
El fenmeno fsico de pandeo, es un problema de inestabilidad que se presenta en estructuras sometidas a acciones mecnicas o trmicas (que generen un estado de esfuerzos en compresin); si estas solicitaciones (cargas, gradientes, etc) rebasan un cierto valor denominado crtico, las estructuras describen un cambio repentino y brutal que se hace notar con la aparicin de ondulaciones y/ pliegues. Este cambio de forma est ligado generalmente a efectos geomtricos no - lineales. La nocin bsica de carga crtica ha sido propuesta a partir de los trabajos realizados por Euler, y el comportamiento corresponde a una estructura ideal. Segn este concepto, una estructura es estable hasta que alcanza un valor crtico. Una imperfeccin infinitamente pequea es suficiente para que la estructura pierda su equilibrio entorno de estabilidad. El pandeo depende principalmente de algunos parmetros tales como: la carga aplicada, la geometra, los esfuerzos esperados, las imperfecciones geomtricas iniciales, y las condiciones de frontera [Notas de mecnica estructural II. 2002].
La prdida de estabilidad est caracterizada por un punto de bifurcacin por un punto lmite. Por tanto, a travs de este tipo de pandeo se pueden identificar y clasificar los fenmenos de estabilidad, en formas de estabilidad.
- 15 -
2.1.1 Pandeo por bifurcacin con cada brusca de rigidez

Est tipo de pandeo corresponde, en ausencia de imperfecciones geomtricas iniciales, al caso de un cascarn cilndrico sometido a compresin axial, o al caso de un cascarn esfrico bajo presin externa. En el punto B de la figura II.1, la estructura (perfecta) pasa de su forma inicial (fundamental), a su forma adyacente (pandeada) que es estable para N < Ncr . La prdida de rigidez de la nueva forma (pandeada) ocasiona que la estructura no sea capaz de mantenerse en equilibrio estable.
N x
Ncr N'cr N''cr 0 B'
y B''
Figura II-1. Pandeo por bifurcacin con cada brusca de rigidez
En la figura II-1 la rama B-y, que representa la zona de equilibrio inestable del pandeo poscrtico, podra ser vecina o cercana a la rama O-B correspondiente a la forma inicial. La presencia de imperfecciones en la estructura en el proceso de carga, podra conducir a un pandeo prematuro, la estructura podra pandearse antes que la carga crtica esperar su valor crtico Ncr. Siguiendo las condiciones del proceso de carga, para este caso, el pandeo se puede producir para un valor inferior al crtico N'cr < Ncr [Snchez. 2003].
- 16 -
2.1.2 Pandeo terico de estructuras sometidas a compresin axial.

Acorde a la teora de pequeas deformaciones, la forma clsica para calcular el esfuerzo crtico de pandeo en cascarones cilndricos sometidos a compresin axial fue dada a conocer por Timoshenko [Notas de Mecnica estructural II. 2002], aunque no fue el nico, este esfuerzo es evaluado como se muestra en la ecuacin II-1:
cr =
3 1 2
Et R
(II-1)
donde: E = mdulo de Young t = espesor de las paredes del cascarn R = radio medio del cascarn
= relacin de Poisson
El pandeo de un cascarn sin imperfecciones, depende de las caractersticas mecnicas de material, de la relacin R/t y la relacin L/R. El efecto de la longitud (L) se puede estudiar utilizando el parmetro de Batdorf (ecuacin II-2), que es adimensional y depende de las caractersticas geomtricas y mecnicas del cascarn. El parmetro de Batdorf establece cuando un cilindro se considera corto o largo, ver seccin 2.2.1.
z=
L2 1 2 Rt
1 2
(II-2)
La ecuacin II-1 es aplicable para los modos axisimtricos, figura II-2 (a), y el pandeo no axisimtrico figura II-2 (b), que se caracteriza por la aparicin de ondulaciones horizontales y verticales. Sin embargo, es bien sabido que en cilindros reales la carga de pandeo ocurre por debajo de lo estipulado por la ecuacin II-1. Usando la teora de grandes desplazamientos y reconociendo que los cilindros se pandean en un modo diamantado, figura II-2 (c), Karman y Tsien fueron capaces de ofrecer una explicacin de la discrepancia entre la teora y la experimentacin. Ellos obtuvieron una curva de pandeo postcrtico, figura II-2 (d), (que muestra - 17 la presencia de grandes
deformaciones), obtuvieron configuraciones de equilibrio estable a causa de cargas considerablemente por debajo de la carga de pandeo, y ellos sugirieron que el cascarn puede saltar repentinamente desde el estado no pandeado a esas configuraciones de equilibrio de grandes deformaciones. En una extensin de la investigacin de Karman sobre grandes deformaciones, Donell y Wan muestran que imperfecciones iniciales reducen grandemente la mxima carga que un cilindro puede soportar. Una visin adicional dentro del comportamiento de pandeo postcrtico de cilindros fue obtenida cuando Yoshimura demostr el patrn diamantado con el que pandean los cascarones cilndricos, desarrollado en la superficie. La posibilidad de esfuerzos de membrana de configuracin libre bajo grandes desplazamientos, ocurre necesariamente para pequeas cargas para mantener el equilibrio en el intervalo de pandeo postcrtico [Chajes. 1969].
a) Modo axisimetrico
b) Modo no axisimetrico
Carga axial Acortamiento axial c) Modo diamantado d) Curva de postpandeo
Figura II-2. Modos de pandeo
Las placas curvas en compresin axial se comportan muy similarmente como cilindros cargados axialmente. Por consiguiente es posible ilustrar algunas de las caractersticas de pandeo de cilindros considerando un modelo de placa curva. El modelo representado en la figura II-3 consiste de los mismos elementos como los usados en el - 18 -
modelo para simular pandeo de placas [Chajes. 1969]. Este difiere del ltimo slo en la existencia de un nudo kink en el miembro transversal, representando la curvatura del cascarn en el estado descargado.
Para un desplazamiento lateral, d, la energa potencial total es definida como la diferencia entre la energa interna de deformacin (U) y el trabajo generado por las fuerzas externas (Wext), V = U Wext
U Wext =
4Cd 2 Pd 2 + K2 L L2
(II-3)
P P e
d h
P L
e P
(a)
(b)
Figura II-3. Modelo de panel curvo
En el cual , es el acortamiento axial de la barra transversal y est dado por
d d h L 2
(II-4)
- 19 -
Si la ecuacin II-4 se sustituye en la ecuacin II-3 y la variacin del potencial V resultante es igualado a cero ( V = 0 ), se obtiene:
P=
4C K 2 + d + 2h 2 3dh L 2L
(II-5)
La ecuacin II-5 es ms sencilla para interpretar si se introduce un parmetro de deformacin no dimensional:
d h
(II-6)
y si se define una relacin arbitraria entre la constante rotacional y la extensional del resorte:
C=
Kh 2 8
(II-7)
Por tanto sustituyendo las ecuacines II-6 y II-7 en la ecuacin II-5, nos conduce a:
P=
h2 K 3 + 2 3 2L
(II-8)
En un anlisis lineal, los esfuerzos de membrana resultantes y sus deformaciones por cambio en la geometra de la estructura, se desprecian. Esto permite eliminar los trminos que involucran y su influencia en la ecuacin II-8, y resulta en el criterio de pandeo lineal:
Pcr =
3h 2 K 2L
(II-9)
- 20 -
La carga de pandeo lineal del modelo tomado por la ecuacin II-9 es mostrada en la figura II-4 a travs de la curva 12. Este corresponde a la carga de pandeo del actual cilindro tomado por la ecuacin II-1.
Para un anlisis por pandeo donde los esfuerzos de membrana sean tomados en cuenta, el comportamiento para pandeo es descrito por la ecuacin II-8 y la curva 1-3-5 de la figura II-4. El modelo no lineal empieza a deformarse con la carga de pandeo lineal, Pcr, pero despus de esto la carga requerida para mantener el equilibrio disminuye rpidamente conforme la deformacin se incrementa. La energa de deformacin, la cual ha aumentado previa al pandeo, decrece posterior al pandeo y esto resulta en la cada de carga observada. Eventualmente, como las deformaciones continan incrementando en su magnitud, los esfuerzos de membrana que estaban originalmente a compresin cambian a tensin, y la carga se incrementa de nuevo. En un cilindro completo, el estado de distorsiones menores es caracterizado por un patrn de pandeo consistente de formas de diamantes, figura II-2(c), y uno de los cuales se comporta similar al modelo de la figura II-3(a).
1
2
Figura II-4. Comportamiento de pandeo.
- 21 -
Para comprender por qu los cilindros se pandean con cargas por debajo de la esperada, la teora elstico lineal requiere que sean consideradas las imperfecciones geometricas, as como en el comportamiento de grandes desplazamientos. Si el modelo es cargado excntricamente, como se muestra en la figura II-3 (b), el trmino
2 Ped es L
adicionado a la energa potencial y la formacin de la carga se relaciona como sigue:

1 1 + e h
h2 K P= 3 + 2 3 2L
(II-10)
La ecuacin II-10, en la cual
e = 0.1 , se muestra en la figura II-4 como la curva 0-6-7. h
De acuerdo a esta curva, el pandeo ocurre en el punto 6, a una carga por debajo de la tomada por la teora lineal. Como en el caso de una cpula , el pandeo toma lugar como resultado de un snap through (aparicin repentina del modo de pandeo) para una configuracin estable. Pequeas deformaciones en columnas y placas pueden despreciarse sin variaciones importantes en la carga de pandeo, no as en placas curvas donde los efectos de imperfecciones son muy influyentes [Chajes. 1969].
2.1.2.1 Tipos de pandeo terico de estructuras cilndricas

Bsicamente existen dos tipos de pandeo terico que se generan en estructuras cilndricas.
Pandeo simtrico de cascarones cilndricos sometidos a compresin axial uniforme. Si un cascarn cilndrico est sometido a compresin en la direccin de su eje longitudinal puede experimentar un pandeo simtrico, ver figura II-5a. El esfuerzo crtico de pandeo para el caso de configuracin axisimtrica est dado por:
m 2 2 N cr E L2 = D + 2 tL2 t R D m 2 2
cr =
(II-11)
- 22 -
El valor mnimo de esta expresin se conoce como el "esfuerzo crtico" de pandeo de cascarones cilndricos sometidos a compresin axial, conocida tambin como "esfuerzo clsico" y corresponde a la ecuacin II.1
a) Modo axisimtrico
b) Modo no axisimtrico
Figura II-5 Modos de pandeo tpicos para cilindros perfectos cargados axialmente.
Esfuerzo crtico de pandeo para el caso de configuracin no axisimtrica. En la realidad la configuracin deformada generada por el efecto de pandeo en estructuras axisimtricas de pared delgada experimentan configuraciones que distan mucho de aquellas cuya aplicacin de carga y comportamiento es perfectamente simtrica, presentndose por consecuencia configuraciones no axisimtricas por lo que el esfuerzo crtico en estructuras perfectas puede ser calculado mediante la ecuacin II.12, ver figura II-5b [Notas de mecnica estructural II. 2002]:
cr = k x
2D
L
2
; kx
(m =
+ 2 m2
(m +
4 2
12 z 2 m 2
2 2
(II-12)
- 23 -
donde:
L2 1 2 z= Rt
1 2
Et 2 nL Parmetro de Batdorf ; = ; D= Rigidez a flexin R 12 1 2
m = configuracin axial; n = configuracin circunferencial
2.1.3 Efecto de las imperfecciones en la estabilidad de cascarones cilndricos

En ciertas investigaciones se ha mostrado que las imperfecciones geomtricas iniciales son consideradas como la principal causa de la amplia dispersin de resultados experimentales (ver figura II.6), y una pobre correlacin entre las predicciones de la teora linealizada de pequeas deformaciones y los resultados experimentales. A pesar de saber esto, la incorporacin de la sensibilidad de imperfecciones dentro de la prctica de ingeniera no ha sido completamente satisfactoria. Los diseadores de estructuras de cascarones continan haciendo el diseo de la forma acostumbrada en donde usan un factor de reduccin emprico; esto resulta en muchas ocasiones antieconmico e incluso inseguro. Por lo que es importante encontrar un procedimiento racional que involucre estos efectos de imperfecciones y se lleven a reglamentos de dominio pblico [Sanchez. 2003].
Los trabajos de [Donnell. 1934], y [Donnell & Van. 1950] conducen a reflexionar que si se conoce la imperfeccin inicial, entonces la solucin terica producir la carga de pandeo de la estructura. De esta manera se podra pensar en un procedimiento para medir la imperfeccin de un cascarn terminado, y usar estos resultados en un mtodo analtico apropiado para predecir o estimar la carga de pandeo.
En la figura II-6 se muestran los resultados de varias pruebas realizadas a cascarones cilndricos. Se puede observar que se han ensayado cilindros con una amplia relacin radio espesor (r/t) y los resultados muestran una dispersin de resultados muy grande, que oscilan entre 0.8 y 0.1 del valor normalizado de la carga crtica. Es muy claro que existe una diferencia muy fuerte entre la carga clsica de pandeo contra la obtenida en laboratorios y ms aun con la que se puede presentar en cascarones fabricados reales. - 24 -
Sin duda alguna esta dispersin de resultados es debida a las imperfecciones geomtricas, imperfecciones en el material de fabricacin y quizs de los esfuerzos residuales en cascarones cilndricos.
En la fotografa II-1 se pueden apreciar cilindros con fallas estructurales por inestabilidad (pandeo). Estas fotografas ejemplifican un modo de pandeo diamantado. Es importante hacer notar que si no se considera el efecto de imperfecciones geomtricas, se pueden presentar fallas estructurales como las mostradas en dicha fotografa.
Esfuerzo de pandeo normalizado k/
Esfuerzo ideal: Esfuerzo crtico clsico
Experimentos
Relacin radio espesor r/t
Figura II-6 Esfuerzos de pandeo experimentales de cilindros isotrpicos comprimidos axialmente, Harris et al., 1957, [Rotter.1996]
- 25 -
Fotografa II-1 Configuracin deformada de pandeo a compresin axial, falla en servicio y prueba en laboratorio, Kndel and Schultz, 1992, [Rotter.1996]
2.1.4 Formas
de
inestabilidad
considerando
la
influencia
de
las
imperfecciones geomtricas
Se sabe que los cascarones de pared delgada lisos o atiesados muestran una muy favorable resistencia con relacin a su peso. Por lo que en la actualidad son muy usados en la ingeniera (civil, aeronutica, espacial) especialmente cuando su ligereza se convierte en un factor importante. Sin embargo el mayor problema al usar este tipo de estructuras es la sensibilidad de presentarse fallas de estabilidad por pandeo. Ahora bien cuando se hable de pandeo de cascarones de pared delgada, se debe distinguir entre colapso por punto limite y por bifurcacin, debido a que un anlisis por punto de bifurcacin es un caso terico que difcilmente se da en la prctica, mientras que un anlisis por punto limite considerando las imperfecciones geomtricas se acerca ms al comportamiento real de un cascarn cilndrico sometido a carga axial (ver resultados de esta tesis). Para obtener los niveles de carga crtica, se pueden realizar dos tipos de anlisis: un anlisis asinttico, y/o un anlisis no-lineal.
- 26 -
Al aplicar un anlisis (post-crtico) de un cascarn cargado axialmente, el desplazamiento radial inicial wb que puede experimentar la estructura debido a pandeo, ser nulo hasta llegar a obtener la carga de bifurcacin c = N N CL en el punto "B", ver figura II.7. Siguiendo la bifurcacin, la falla inicial de la estructura perfecta ser caracterizada por el crecimiento en la deformacin axisimtrica representada por la rama descendente BD, y por la carga decreciente .
Carga
= N/NCL
Punto de bifurcacin Pandeo Postcrtico de la estructura perfecta
c s
Carga por punto lmite
E F
Pandeo Postcrtico de la estructura con imperfecciones
Desplazamiento del pandeo modal w'b
Figura II-7. Punto de bifurcacin y punto limite mediante un anlisis asinttico.
De otra manera, si se emplea un anlisis no-lineal, el cascarn perfecto cargado axialmente se deformar siguiendo la rama ascendente OA hasta alcanzar la carga lmite L mxima, en el punto A, ver figura II.8. Sin embargo, el punto de bifurcacin "B" para este caso se encuentra en la trayectoria de la curva O y A. As, una vez que la carga de bifurcacin C ha sido alcanzada, la falla inicial por pandeo del cascarn perfecto est caracterizada por un crecimiento acelerado en la deformacin axisimtrica siguiendo la carga axial decreciente.
- 27 -
Carga
= N/NCL
Carga por punto lmite del cascarn perfecto
L
Carga de bifurcacin
A B C
Pandeo post crtico del cascarn perfecto
Carga por punto lmite
E F
Carga por punto lmite del cascarn con imperfecciones
Desplazamiento total w
Figura II-8. Punto de bifurcacin y punto limite mediante un anlisis no-lineal.
En el caso de cascarones reales, que presentan inevitablemente imperfecciones iniciales, mediante los dos mtodos que definen la respuesta estructural, se seguir una trayectoria representada por la rama fundamental OEF, con una falla que se presenta como un "snapthrough" en el punto E que representa la carga mxima[Sanchez. 2003].
Como se ver en el desarrollo de esta tesis el mtodo que se aproxima ms a la realidad, es un anlisis no lineal por punto limite, siempre y cuando sean consideradas las imperfecciones geomtricas del cascarn cilndrico.
2.1.5 Influencia de los esfuerzos residuales en la estabilidad de cascarones

Debido a la complejidad de separar las imperfecciones geomtricas y los esfuerzos residuales en cilindros de pared delgada no ha sido posible elaborar formulaciones tericas que agrupen estos dos fenmenos. Sin embargo con el desarrollo de software y hardware los mtodos numricos han tomado mucha fuerza ya que pueden resolver este tipo de problemas mediante mtodos como el del elemento finito.
- 28 -
Los esfuerzos residuales y las imperfecciones geomtricas en silos se generan debido a los procesos de fabricacin empleados en la manufactura de estas estructuras. Tal es el caso del proceso de rolado que genera esfuerzos iniciales para formar el cilindro, ver fotografa II.2. Otra contribucin importante es sin duda el proceso de soldadura al unir dos virolas o rolos para formar el silo, fotografa II.3. Esto se refiere exclusivamente a silos que son conocidos como soldados, sin embargo en silos modulares (atornillados) hay otros procesos que pueden generar esfuerzos adicionales como el proceso de punzonado y cejado (doblez).
Como se puede apreciar en la serie de fotografas siguientes, existen esfuerzos residuales generados por los diferentes procesos de fabricacin, as como imperfecciones iniciales generadas ya sea por estos procesos de fabricacin o por el manejo que se le da en taller a este tipo de estructuras.
Fotografa II-2 .Rolo o virola rolada (Cortesa Silos y Camiones S.A. de C.V.)
- 29 -
Fotografa II-3. Esfuerzos residuales generados por el proceso de soldadura (Cortesa Silos y Camiones S.A. de C.V.)
Fotografa II-4. Hojas empaquetadas las cuales estructurarn al silo atornillado (Cortesa Silos y Camiones S.A. de C.V.)
- 30 -
Fotografa II-5. Silos atornillados (cortesa Silos y Camiones S.A. de C.V.)
2.2 Criterios tpicos de diseo

En esta seccin se abordarn bsicamente dos metodologas que se usan en el diseo de silos metlicos, el mtodo de Baker [Baker. 1972] que es muy usado por las empresas que fabrican silos en los Estados Unidos de Amrica y dada la relacin que tiene Mxico con este pas, tambin ha sido muy aceptada por los fabricantes mexicanos. Sin embargo en la Unin Europea ha sido bien acogido el Eurocdigo 3, diseo de estructuras de acero en su parte 4.1 de Silos. En la actualidad el Eurocdigo es bien visto por los pases sudamericanos. Y es importante conocer su uso y aplicacin.
2.2.1 Compresin axial en cilindros de pared delgada, [Baker. 1972]

Para cilindros sujetos a compresin axial el fenmeno de pandeo est gobernado por la siguiente ecuacin: - 31 -
cr 2E t = Kc 12(1 2 ) L
(II-13)
El efecto de la plasticidad puede considerarse por el factor de correccin por plasticidad
, sin embargo dado que se esta trabajando en pandeo elstico, toma el valor de 1.
Para una mayor explicacin de este factor se recomienda consultar la seccin 10-7 de Baker [Baker. 1972].
Para cilindros cortos, el coeficiente de pandeo K c puede expresarse aproximadamente por la ecuacin II-14, siempre y cuando se cumpla la desigualdad z < 2 K c 0 2 3 . El efecto de la longitud del cilindro bajo el efecto de pandeo por compresin axial, puede estudiarse utilizando el parmetro de Batdorf, Z, que es un parmetro adimensional que depende de las caractersticas geomtricas del material del cascarn.
K c = K c0 +
12 2 z 2 4 K c0
(II-14)
K c 0 se expresa en funcin de las modos circunferenciales n y verticales m as como de

parmetros geomtricos como longitud (L) y radio (R) [Notas de mecnica estructural II. 2002]:
K co
(m =
+2 m2
(II-15)
Donde =
nL R
Y para cilindros moderadamente largos z > 2 K c 0 2 3 por:
Kc =
4 3
(II-16)
- 32 -
Donde : z= L2 1 2 Rt
Baker [Baker. 1972] sugiere usar K c 0 = 1 para cilindros con extremos simplemente apoyados, K c 0 = 4 para extremos empotrados. El factor ha sido incluido para tomar en cuenta la diferencia entre los resultados tericos y experimentales. El valor para el coeficiente para el pandeo clsico de cilindros sometidos a compresin axial es 1. Existe un gran nmero de pruebas experimentales cuyos resultados se pueden consultar en las referencias 10-3 y 10-4 de Baker [Baker. 1972].
Si se minimiza la ecuacin II-14, se encuentra la expresin para el esfuerzo de pandeo en cascarones cilndricos moderadamente largos sometidos a carga axial.
Et cr = Cc r
(II-17)
Donde:
Cc =
[3(1 )]
2
= 0.6 Siempre que z > 2 K c 0
2 3 ; ( = 0.3)
Datos experimentales sugieren que el valor de usado para calcular cr para cilindros en compresin con comportamiento isotrpico, puede tomarse como [Baker. 1972]:
= 1 0.901(1 e )
Donde:
(II-18)
1 R 16 t
- 33 -
2.2.2 Eurocdigo 3 [Eurocode 3. 2003b]

El Eurocdigo 3 Diseo de estructuras de acero Parte 4.1 Silos, tambin incorpora resultados sobre el tema de pandeo en cilindros axialmente cargados. Estas normas se han preocupado por establecer una relacin entre las imperfecciones encontradas en este tipo de estructuras y el esfuerzo crtico de pandeo. Estos estndares europeos proponen el uso de un factor de calidad de fabricacin, mediante el cual se calcula un valor para reducir el esfuerzo critico de pandeo y se toma en cuenta en el diseo de silos; que como anteriormente se ha establecido, son cascarones de pared delgada.
La valoracin de los esfuerzos reales por pandeo comienza con la medicin de la calidad en la fabricacin debido a que esto es una influencia importante en la capacidad del esfuerzo crtico de pandeo que soporta un cilindro a compresin axial. El Eurocdigo establece que el factor de calidad Q debe ser valorado como se muestra en la tabla II-1, con base en el valor seleccionado, se establece la amplitud de la imperfeccin wok, que es tomada como:
wok =
t r Q t
(II-19)
Tabla II-1 Clasificacin de cascarones cilndricos con base en la calidad de fabricacin. Tolerancia fabricacin construccin. Clase A Excelente 40 Solo permitido para clase 3 Clase B Clase C Alta Calidad Normal 25 16 Obligatorio clase 1 para en de Descripcin la Parmetro Calidad, Q de Restricciones segn su clase
Para la clase 3 (ver tabla II-2) el Eurocdigo establece que las fuerzas internas en un silo deben ser determinadas usando un mtodo numrico de anlisis, por ejemplo elementos finitos. Para los de clase 2 las fuerzas pueden encontrarse mediante teora - 34 -
de membranas u otro mtodo numrico vlido adems de otras restricciones. Para silos clase 1, la teora de membrana puede emplearse.
Cuando la junta en las paredes del cilindro sea a traslape, la imperfeccin wok debe incrementarse en un 50%.
Tabla II-2 Clasificacin de silos segn su clase. Clase Clase 3 Situaciones de diseo
Silos soportados al piso o silos soportados con un faldn con capacidad mayor a las 5000 toneladas Silos discretamente soportados con capacidad mayor a 1000 ton Silo con capacidad mayor a 200 toneladas en cualquiera de los casos siguientes:
a) Descarga excntrica b) Cargas locales por reparaciones (local match loading) c) Llenado excntrico (no simtrico) Clase 2
Todos los silo cubiertos por esta norma y no incluidos en otra clase
Clase 1
Silos con capacidades entre 10 y 100 toneladas
Silos con capacidades por debajo de las 10 toneladas no estn cubiertos por el Eurocdigo 3 Parte 4.1
2.2.2.1 Medicin de las imperfecciones geomtricas despus de la fabricacin

Despus de la fabricacin, el silo debe ser verificado para garantizar que las tolerancias de fabricacin correspondan a las de diseo.
La profundidad de cavidades iniciales debe medirse usando calibradores (gages) de longitud lg, esta longitud se debe medir en diferentes localizaciones como sigue: - 35 -
Cuando se tengan esfuerzos en compresin, incluyendo soldaduras, en direcciones meridional y circunferencial se tiene (ver figuras II-9 ,a y b):
l gx = 4 rt
(II-20)
Cuando ocurran esfuerzos en compresin o esfuerzos cortantes, la medicin de las imperfecciones en la direccin circunferencial debe ser hecha acorde a las figuras II-9, II-10, II-11 con un calibrador (gage) de longitud:
l g = 2.3 l 2 rt
( )
0.25
pero l g r
(II-21)
Donde l es la longitud axial del segmento de lmina estudiada, la cual puede tomarse como la longitud entre anillos atiesados.
Adicionalmente en los cruces de soldaduras la longitud en ambas direcciones es tomada como:
lgw = 25t o l gw = 25tmin pero con lgw 500mm

Donde tmin es el espesor ms delgado de la lmina o la soldadura.
r t
(II-22)
Wox
t
Wox
lgx lgx
a) Medicin de un meridiano
b) Primeras mediciones sobre un crculo
Figura II-9 Medicin de imperfecciones, slo carga axial - 36 -
t r
Wo
t
Wox
lgx lg
Soldadura
c) Primera medicin transversal de la soldadura
d) Segunda medicin de la circunferencia del crculo
Figura II-10 Medicin de imperfecciones, cuando se incluyen efectos de cortante.
Wox Wo Wow
Soldadura
lgw lgx o lg o lgw
Soldadura
Wow
e) Segunda medicin transversal de la soldadura con un calibrador especial (gage)
f) Medicin en una seccin circular de soldadura
Figura II-11 Medicin imperfecciones, en soldaduras.
2.2.2.2 Verificacin de la tolerancia despus de la fabricacin.

La profundidad de las imperfecciones iniciales debe asegurarse en trminos de parmetros de imperfeccin Uox, Uo, Uow dados por [Eurocode 3. 2003a]: - 37 -
U ox =
wox l gx
U o =
wo l g
U ow =
wow l gw
(II-23)
El valor de los parmetros de imperfeccin Uox, Uo, Uow, deben satisfacer las condiciones siguientes [Eurocode 3. 2003a]:
U ox U o,max
U o U o,max
U ow U o,max
(II-24)
Donde Uo,max es el parmetro de tolerancia de imperfeccin para la clase de tolerancia de fabricacin seleccionada, el cual es tomado de la tabla II-3.
De acuerdo al Eurocdigo 3. Parte 4.1, la mxima medida de imperfeccin no debe exceder de:
wod = 0.6wok
(II-25)
Tabla II-3 Valores del parmetro de imperfeccin Uomax Tolerancia fabricacin construccin. Clase A Clase B Clase C Excelente Alta Calidad Normal 40 25 16 0.006 0.01 0.016 en de Descripcin la Parmetro de Valor de Uo,max Calidad, Q
2.2.2.3 Parmetro de reduccin

Las imperfecciones descritas en las secciones anteriores, tienen un efecto serio en el esfuerzo elstico de pandeo de cilindros delgados. La reduccin del esfuerzo es
representada por un factor de reduccin para pandeo elstico, el cual toma un valor entre 0 y 1. Se definen tres valores: o para compresin axial sin presin interna, pe - 38 -
cuando la presin interna incrementa la resistencia a pandeo, y pp para la condicin donde la presin interna causa una reduccin en la resistencia de pandeo por la induccin de un pandeo plstico (pata de elefante).
El factor de reduccin para cilindros no - presurizados y sometidos solamente a compresin axial, esta dada por:
o =
0.62 w 1 + 1.91 ok t
1.44
(II-26)
El parmetro para esfuerzo no uniforme para compresin circunferencialmente uniforme (figura II-12a) puede ser tomado como 1, y puede ser usado
conservadoramente para todos los casos. Cuando el esfuerzo de compresin axial no es uniforme alrededor de la circunferencia del cilindro (ver figura II-12b), el efecto debe ser considerado por el parmetro , el cual debe determinarse de la distribucin de esfuerzos elstico lineal de la distribucin de esfuerzos de compresin axiales actuantes. La distribucin de los esfuerzos de membrana a compresin axial alrededor de la circunferencia del nivel en cuestin debe transformarse como se indica en la seccin 5.3.2.4 prrafo 8 del Eurocdigo 3 [Eurocode 3. 2003b]. Este parmetro solo es vlido para compresin axial (uniforme o no) y el mximo valor que tomar ser igual a 1.
- 39 -
Cilindro
Cilindro
(a)
conservadoramente = 1 .
(b)
Figura II-12 a) Compresin uniforme = 1 , b) Compresin no uniforme < 1 ,
- 40 -
Captulo III
Modelado y pandeo terico
3 3.1 3.1.1 3.1.2 3.1.3 3.1.4 3.1.5 3.1.6 Modelado y pandeo terico ............................................................................. 42 Probeta terica ................................................................................................ 42 Caractersticas geomtricas y mecnicas de los modelos .............................. 42 Pandeo terico de la probeta a estudiar.......................................................... 44 Pandeo axisimtrico ........................................................................................ 47 Pandeo no axisimtrico ................................................................................... 49 Tcnica de Rish y Kovalesky .......................................................................... 50 Modelacin controlada .................................................................................... 50
Resumen
Para comprender el fenmeno de pandeo, se realizan algunos modelos numricos controlados para la obtencin de la carga crtica en cilindros de pared delgada, mediante diversas tcnicas previamente estudiadas en el captulo anterior. Adems en este captulo se realizan modelos de anlisis, y se proponen tcnicas numricas con el fin de establecer una metodologa de anlisis que ser utilizada en el siguiente captulo y que sirvan como base para encontrar la carga de inestabilidad de los especimenes fabricados.
3 Modelado y pandeo terico

Debido a la complejidad del anlisis del fenmeno de pandeo se hace necesaria la modelacin numrica por medio de la tcnica de los elementos finitos. Primeramente se analizar un modelo de dimetro "D", espesor "t" y longitud L generndole una imperfeccin geomtrica controlada. Este modelo servir como base de calibracin y para comenzar a establecer como influyen las imperfecciones geomtricas en la obtencin de la carga critica.
Para darle mayor sustento a estos resultados se probarn experimentalmente cinco especimenes (Captulo IV), los cuales se ensayarn en una mquina universal para realizarles una prueba a compresin que proporcione la historia de carga vs deformacin y as encontrar su carga crtica de pandeo experimental. A estos especimenes se les medirn sus imperfecciones geomtricas iniciales, que sern incluidas en un modelo numrico en computadora (Captulo IV), para posteriormente obtener el valor de la carga crtica, usando la tcnica de los elementos finitos (MEF), mediante el anlisis de estabilidad de Euler, adems de usar una tcnica de anlisis no lineal.
3.1 Probeta controlada

Es necesario partir de datos que se puedan comprobar, es por esta razn que se propone la generacin de un modelo de dimensiones e imperfecciones geomtricas controladas, para posteriormente realizar una comparacin entre los resultados tericos (Seccin 2.1.2), los clsicos de diseo (Seccin 2.2) y los experimentales (Captulo IV).
3.1.1 Caractersticas geomtricas y mecnicas de los modelos

Para la seleccin de las dimensiones de los modelos a estudiarse se toman consideraciones como: material que comercialmente se pueda conseguir, dimensiones comerciales, capacidad de la mquina de compresin a usarse. Estos parmetros son
- 42 -
importantes ya que servirn como base para las dimensiones de construccin de los especimenes de tipo experimental en el captulo IV.
Se puede conseguir lmina negra o galvanizada en espesores hasta de calibre 30 (0.3 mm), acero al carbn tipo ASTM 570 A36, la capacidad de la mquina universal es hasta de 100 toneladas y se puede experimentar con una probeta de dimensiones de hasta de 600 mm de altura; el dimetro puede ser variable, aunque la mquina universal no tiene mordazas para compresin (Se debern fabricar).
Una vez analizados, los detalles de aspecto y comerciales, se procede a seleccionar un cilindro con las siguientes caractersticas generales:
Las caractersticas geomtricas del modelo son las siguientes: Dimetro, D =: 30 cm; Altura efectiva, h o L: 30 cm; Espesor, t = 0.03 cm (cal 30)
L t = 0.299 0.3 r r
Valor tpico para silos, [Urrutia. 2002]
Las caractersticas mecnicas del material son: Material de construccin: Lmina galvanizada cal 30 (0.030 cm). Mdulo de elasticidad, E = 2.1E06 kg/cm2 Mdulo de Poisson, = 0.3 Esfuerzo de fluencia fy = 2530 kg/cm2
La relacin ancho espesor del cilindro es de:
R 15 = = 500 t 0.03
El parmetro de Batdorf para este modelo se determina con la ecuacin II-2:
- 43 -
z=
300 2 1 0 .3 2 (150)(0.3)
1 2
= 1908
Para cilindros moderadamente largos y con extremos empotrados (seccin 2.2.1), se tiene:
2 K c 0 2 3 = 11.4
Por lo que este modelo tiene una relacin radio espesor de 500 y es de longitud moderadamente larga z > 2 K c 0 2 3 .
3.1.2 Pandeo terico de la probeta a estudiar

Dado que se debe partir de datos tericos del modelo en estudio, se necesita evaluar el comportamiento de un cilindro perfecto de dimensiones antes descritas y evaluar la carga crtica clsica. As de la ecuacin II-1, se tiene lo siguiente:
(2.1 106 )(0.03) kg = 2541 2 2 15 cm 3 1 0 .3
cr =
Donde la carga crtica en un rea de 2rt est dada por: Pcr = 2rt cr = 7184.51 kg
Como es bien sabido, el pandeo se ve caracterizado por la generacin de modos de pandeo, por lo que resulta importante conocer cuales son los primeros modos de pandeo de nuestro cilindro y con esto hacer algunas deducciones iniciales.
Como anteriormente se estableci, se utilizar la tcnica numrica de los elementos finitos (MEF) con la ayuda de Ansys 8.1. Se va a trabajar con el elemento shell 63 (figura III-1). El elemento shell 63 tiene capacidades de flexin y membrana. Ambas para cargas en el plano y normales. Este elemento permite 6 grados de libertad en cada nodo, 3 traslaciones en x, y, z y 3 rotaciones alrededor de x, y, z. Tiene capacidades para anlisis de grandes desplazamientos [ Ansys. Manual R8.1]. - 44 -
El mtodo a usarse para la obtencin de los modos y cargas de pandeo ser mediante la extraccin de valores y vectores caractersticos que darn solucin a los modos de frecuencia y al anlisis de pandeo. Ansys maneja varias tcnicas para la solucin del problema de valores y vectores caractersticos como la de Subspace y la de Block Lanczos. Para esta tesis se usar la tcnica de extraccin de Lanczos por ser ms estable en su solucin y mostrar menos problemas de convergencia en la solucin del fenmeno estudiado. Es importante hacer notar que el mtodo de extraccin de Lanczos proporciona el punto de bifurcacin de la estructura estudiada, por lo que tambin ser importante realizar una anlisis no lineal y obtener as la carga de inestabilidad que sea ms cercana a la realidad.
Figura III-1 Elemento Shell 63
Como se puede apreciar en la figura III-2 el mtodo de extraccin de Lanczos, es un anlisis lineal y por ello debe ser usado conservadoramente, ya que el resultado que se obtiene es el punto de bifurcacin de la estructura pandeada. Es conveniente realizar un anlisis no lineal para obtener resultados ms adecuados [ Ansys. Manual R8.1]
En la figura III-2 se puede observar que para un cilindro a compresin con imperfecciones iniciales y bajo la aplicacin gradual de carga, ste llega hasta un punto de bifurcacin, que corresponde a punto donde la estructura puede pasar de su estado estable a su estado inestable. Se puede apreciar que a este punto de bifurcacin se llega de manera lineal, que es como el mtodo de Lanczos realiza la extraccin de los valores y vectores caractersticos, es decir el mtodo de Lanczos proporciona la - 45 -
solucin al problema de inestabilidad a travs de un anlisis lineal. Por lo anterior para el estudio de estructuras con imperfecciones geomtricas iniciales, es de suma importancia realizar una anlisis no lineal (figura III-2), y as encontrar el punto lmite, que se asemejar ms a un comportamiento real.
P
Punto de Bifurcacin Carga lmite
Eigenvalor (Lanczos)
Analisis no lineal
Figura III-2 Punto de bifurcacin y punto limite
Para la realizacin del anlisis no lineal, se usar la tcnica del Arc Length, que es un mtodo usado por Ansys para obtener soluciones numricamente estables para problemas con soluciones singulares (matriz de rigidez tiende a cero o es cero, matriz singular).
En la figura III-3 se muestra el modo de pandeo axisimtrico tpico para cascarones cilndricos bajo la accin de una fuerza de compresin axial, en esta figura se muestra el cascarn cilndrico perfecto analizado en Ansys. Se us una geometra con condiciones de frontera doblemente empotrada con elementos shell 63 y usando la tcnica de Lanczos para obtener los modos de pandeo
As los tres primeros modos de pandeo son los siguientes:
Modo 1 Pcr =7147.9067 kg, corresponde a un modo circunferencial n = 0 y un modo longitudinal m = 12 - 46 -
Modo 2 Pcr = 7147.9103 kg, corresponde a un modo circunferencial n = 1 y modo longitudinal m = 12 Modo 3 Pcr = 7147.5982 kg, corresponde a un modo circunferencial n = 2 y modo longitudinal m = 12
De este primer anlisis se puede observar que los modos de pandeo son muy cercanos entre si, por lo tanto se deduce que es una estructura sumamente inestable, y que precisamente se caracteriza porque los modos de inestabilidad son muy prximos entre si. [ Ansys. Manual R8.1].
MX
Z Y
Z Y
MN
< a) Estructura ideal, cargada b) axialmente De la configuracin c) Configuracin del tipo
inicial al modo pandeado, senoidal en la direccin axial m=1
Figura III-3 Cilindro perfecto, primer modo de inestabilidad (n = 0 modo cero circunferencial,)
3.1.3 Pandeo axisimtrico

Se presentan los resultados de diferentes geometras de cascarones cilndricos de pared delgada. Los valores de la tabla III-1 fueron obtenidos a partir de la ecuacin II-6 para el caso de pandeo axisimtrico con diferentes relaciones radio - espesor (R/t), manteniendo R = 15 cm y variando t desde 0.03 a 0.135, donde m representa las configuraciones axiales. De la tabla III-1 se observa que el valor mnimo de los
esfuerzos se presenta para la relacin R/t = 500 que indica la relacin ms alta y - 47 -
corresponde al espesor ms pequeo de los valores seleccionados, la carga crtica para este valor de esfuerzo es de 7187.5 kg.
Tabla III-1 Valores de estabilidad neutra para diversas relaciones R/t, con n=0 m R/t T (cm) 26 22 18 16 14 13 12 500.00 375.00 250.00 187.50 150.00 125.00 107.14 0.03 0.04 0.06 0.08 0.10 0.12 0.14 Esfuerzo critico Carga critica (kg) (kg/cm2) 2542.069 3391.591 5086.655 6779.806 8475.761 10168.277 11862.533 7187.531 12785.997 28764.354 51118.536 79882.161 115000.501 156522.436
100.00 90.00 R/t = 500
Esfuerzo crtico (kg/cm2) x 1000
80.00 70.00 60.00 50.00 40.00 30.00 20.00 10.00 0.00 0 5 10 15
R/t = 375 R/t = 250 R/t = 187.5 R/t = 150 R/t = 125 R/t = 107
n=0
20
25
30
35
40
45
50
m (modos axiales)
Figura III-4 Curva de estabilidad neutra para diversas relaciones R/t - 48 -
En la figura III-4 se muestran todas las configuraciones realizadas, reafirmndose que la carga crtica corresponde a una relacin r/t de 500 y al espesor ms pequeo, que para este caso es 0.30 mm (0.03 cm). En esta misma figura se puede apreciar que el esfuerzo de pandeo es de 2542 kg/cm2 , Sin embargo el modelo matemtico elegido (ecuacin II-11) obtiene la carga crtica terica, pero no soluciona el modo en el cual se obtiene dicha carga, que debera ser .en modos axiales bajos (entre 1 y mximo 6 o 7).
3.1.4 Pandeo no axisimtrico

En la tabla III-2 y la figura III-5 se muestran los valores correspondientes a los esfuerzos crticos de la estructura en estudio, para una configuracin no axisimtrica. Se hizo variar la configuracin axial m de 1 a 10. Se puede apreciar en la tabla III-2 que para cada configuracin axial m, existe su respectiva configuracin circunferencial n (1 a 30), la cual corresponde a la carga crtica buscada. En la grfica de la figura III-5, se puede apreciar que el valor mnimo que toma la curva para cada caso, corresponde al de la carga crtica (o esfuerzo crtico). Para la carga crtica buscada, los esfuerzos crticos de pandeo son siempre constantes, del orden de 2544 kg/cm2, independientemente de la configuracin axial en m elegida.
Tabla III-2 Valores de estabilidad neutra para valores de n y m R/t, R=30 cm, t= 0.03cm m 1 2 3 4 5 6 7 8 9 10 n 8 11 13 15 16 17 18 19 19 20 cr (kg/cm2) 2550.166 2545.01248 2541.97199 2548.1609 2542.07836 2542.93236 2542.06276 2543.21851 2544.95573 2542.79172 - 49 Pcr (kg) 7210.42449 7195.85326 7187.25648 7204.75522 7187.55724 7189.97185 7187.51314 7190.78093 7195.69282 7189.5742
10000 9000
m= 1 m= 4 m= 7 m = 10 m= 2 m= 5 m= 8 m= 3 m= 6 m= 9
Esfuerzo crtico (kg/cm2)
8000 7000 6000 5000 4000 3000 2000 1000 0 0 5 10 15
20
25
30
35
n (ondas circunferenciales)
Figura III-5 Curva de estabilidad neutra para valores de n y m
3.1.5 Tcnica de Rish y Kovalesky

Como se estudi en la seccin 2.1.6 [Baker. 1972] propone un factor de reduccin de la carga critica clsica en donde se ven involucradas las imperfecciones geomtricas y todos los posibles factores que disminuyen la capacidad de estabilidad de un cilindro de pared delgada. As el esfuerzo crtico de pandeo es de 813.67 kg/cm2, que corresponde a una carga crtica de 2335 kg. Es importante hacer notar que este tipo de elementos se ven caracterizados por el uso de factores de seguridad de hasta 3 cuando no hay algn medio rigidizante, por lo que la carga crtica admisible sera de 778.34 kg.
3.1.6 Modelacin controlada

Se hace necesario tener un conocimiento previo de cmo tiende a comportarse un cilindro con dimensiones D, t y longitud L al cual se le inducir una imperfeccin
- 50 -
geomtrica controlada y generada a partir de lo estipulado por el Eurocdigo 3, seccin 4.1 Silos.
Posteriormente se procede a generar una serie de imperfecciones geomtricas tomando como base el factor Q de calidad de fabricacin mencionado anteriormente en la seccin 2.2.2
Este anlisis se realiza usando el software Ansys 8.1 con elementos shell 63 y shell 181 (cuadrilteros y triangulares). Se hace el anlisis de pandeo usando el mtodo Lanczos de Ansys, con densidad sencilla y con doble densidad en el mallado, es por esto ltimo que para el elemento shell 181 se hace una diferencia de densidad de 120 y 240, siendo la densidad de 240 el mallado ms fino. Tambin se realiza un anlisis no lineal mediante el mtodo del Arc Length de Ansys.
Las imperfecciones inducidas al modelo analtico son como las mostradas en la figura III-6, donde se muestra una muesca wok que est en funcin del parmetro de calidad propuesto. Como se explica en la seccin 2.2.2 el parmetro de calidad Q puede tener tres valores, 16, 25 y 40. Para poder estudiar que sucede en valores inferiores y superiores de Q se realizarn imperfecciones con valores de Q desde 2 hasta 280. Las imperfecciones se calculan de acuerdo a la ecuacin II-19 y pueden revisarse en la tabla III-2.
En la figura III-7 se puede apreciar la sensibilidad de la imperfecciones geomtricas respecto al factor de calidad Q, es decir, si el valor de Q tiende a ser ms grande, se habla de que el cilindro tiene pocas imperfecciones, de manera que la relacin entre carga critica del eurocdigo vs carga crtica clsica tiende a ser asinttica a 0.6; para valores mayores de Q = 280, nunca se alcanza el valor normalizado de 1. Mientras que para valores inferiores a Q = 16 la carga crtica se ve muy afectada. Esto tambin se puede apreciar en la figura III-8, sin embargo aqu se aprecia que mientras la metodologa usada por el Eurocdigo sigue reduciendo la carga crtica; otras tcnicas de evaluacin del pandeo parecen tomar una resistencia poscrtica a medida que las imperfecciones son mayores. Esto es probablemente debido a que antes que se - 51 -
desarrolle un pandeo de tipo elstico, se desarrollan plastificaciones locales y una probable resistencia poscrtica del cilindro a compresin debido a endurecimiento del material por efectos de plasticidad, es decir existe un mejoramiento en la rigidez del elemento, sin embargo esto es un fenmeno que no se estudia en este trabajo.
D = 30 lgx;lgw = 6 30 wok
Figura III-6. Esquema de probeta virtual

1.20
1.00
Relacin Pcr/Pcr clas
0.80
0.60
0.40
Pcr/Pcr cl, Eurocodigo Pcr FEM/Pcr cl, Shell 63
0.20
Pcr FEM/Pcr cl, Shell 181, 120 Pcr FEM/Pcr clas, Shell 181, 240
0.00 0 50 100 150 200 250 300
Q, Factor de tolerancia en fabricacion

Figura III-7 Comportamiento de la carga crtica normalizada en relacin con el factor de fabricacin - 52 -
En el grfico de la figura III-7 se puede apreciar que para una Q menor de 16 hay incremento sbito en la carga crtica de pandeo obtenida. En este momento se puede decir que quizs sea por efectos no lineales, como endurecimiento del material o algn tipo de post resistencia. Esto ser discutido ms adelante y en las conclusiones.
Dado a lo anterior se deben usar las ecuaciones recomendadas del Eurocdigo para valores de Q entre 16 y 40, pero para valores menores de 16 se debe tener cuidado en su aplicacin directa, ya que pueden ser muy conservadoras. En el caso de tener imperfecciones geomtricas que correspondan a valores menores de Q igual a16, es recomendable realizar un anlisis no lineal con la medicin de las imperfecciones y as tener un mejor marco de comparacin.
1.20
1.00
Pcr/Pcr cl, Eurocodigo
Relacin Pcr/Pcrclas
0.80
Pcr FEM/Pcr cl, Shell 63 Pcr FEM/Pcr cl, Shell 181, 120 Pcr FEM/Pcr clas, Shell 181, 240
0.60
0.40
0.20
0.00 0.0000
2.0000
4.0000
6.0000
8.0000
10.0000
12.0000
Imperfecciones =Wok/ t
Figura III-8. Comportamiento de la carga critica normalizada en relacin con las imperfecciones
- 53 -
1.20
1.00
Pcr/Pcr cl, Eurocodigo Pcr FEM/Pcr cl, Shell 63
Relacin Pcr/Pcrclas
Pcr FEM/Pcr cl, Shell 181, 120

0.80
Pcr FEM/Pcr clas, Shell 181, 240
0.60
0.40
0.20
0.00 0.0000
0.1000
0.2000
0.3000
0.4000
0.5000
0.6000
0.7000
0.8000
0.9000
1.0000
Imperfecciones =Wok/ t
Figura III-9 Amplificacin de la figura III-8
En la tabla III-3 se puede apreciar que para valores de Q menores a 16, la carga crtica recomendada por el Eurocdigo se ve fuertemente reducida y esto es confirmado por los anlisis no lineales realizados (Arc Length). Tambin se puede apreciar que usando la tcnica de Lanczos, la carga por bifurcacin encontrada est muy por encima de los valores recomendados por el Eurocdigo y por los anlisis no lineales. Esto es debido a que el mtodo de Lanczos es una tcnica lineal para encontrar los valores caractersticos que producen inestabilidad en la estructura (ver seccin 3.1.2 y figura III2).
- 54 -
Tabla III-3 Generacin de imperfecciones y cargas encontradas

ANSYS
SHELL 93
Pcr,
SHELL 63
SHELL 181
Eurocdigo
RELACIN Pcr/Pcrclasica Carga de Bifurcacin No lineal
Punto limite Arc Length kg 354.86 Pcr FEM (Carga de bifurcacin) 240 120 120 3242 2151 2713
Pcr/Pcr cl Pcr/Pcr cl Pcr/Pcr cl
t cm 0.03 0.03 0.03 0.03 0.03 0.03 0.03 0.03 0.03 0.03 0.03 0.03 0.03 0.03 0.03 0.03 0.03 0.03 0.03 0.03 0.03 0.03 0.03 0.03 0.03 0.03 0.03 0.03 0.03 0.03 0.03 0.03 0.03 0.03 0.03 0.03 0.03 0.03 0.03 0.03 0.03 0.03 0.03 0.03
Wok
cm
Area cm2 2.827 2.827 2.827 2.827 2.827 2.827 2.827 2.827 2.827 2.827 2.827 2.827 2.827 2.827 2.827 2.827 2.827 2.827 2.827 2.827 2.827 2.827 2.827 2.827 2.827 2.827 2.827 2.827 2.827 2.827 2.827 2.827 2.827 2.827 2.827 2.827 2.827 2.827 2.827 2.827 2.827 2.827 2.827 2.827
Pcr clsico
Eurocdigo
kg 7,184.51 7,184.51 7,184.51 7,184.51 7,184.51 7,184.51 7,184.51 7,184.51 7,184.51 7,184.51 7,184.51 7,184.51 7,184.51 7,184.51 7,184.51 7,184.51 7,184.51 7,184.51 7,184.51 7,184.51 7,184.51 7,184.51 7,184.51 7,184.51 7,184.51 7,184.51 7,184.51 7,184.51 7,184.51 7,184.51 7,184.51 7,184.51 7,184.51 7,184.51 7,184.51 7,184.51 7,184.51 7,184.51 7,184.51 7,184.51 7,184.51 7,184.51 7,184.51 7,184.51
kg 70.96 187.41 325.17 474.30 628.67 784.19 938.05 1,088.34 1,233.78 1,373.59 1,507.28 1,634.64 1,755.63 1,870.33 1,978.90 2,081.57 2,178.58 2,270.23 2,356.79 2,438.55 2,784.90 3,048.75 3,253.10 3,414.17 3,543.26 3,648.31 3,734.99 3,807.40 3,868.57 3,920.77 3,965.70 4,004.69 4,038.78 4,068.77 4,095.33 4,118.97 4,140.12 4,159.13 4,176.30 4,191.86 4,206.01 4,218.93 4,230.76 4,241.62
120 3020
240 2496 0.01 0.03 0.05 0.07 0.09 0.11 0.13 0.15 0.17 0.19 0.21 0.23 0.24 0.26 0.28 0.29 0.30 0.32 0.33 0.34 0.39 0.42 0.45 0.48 0.49 0.51 0.52 0.53 0.54 0.55 0.55 0.56 0.56 0.57 0.57 0.57 0.58 0.58 0.58 0.58 0.59 0.59 0.59 0.59
Punto limite 0.05
SHELL 63 0.30
2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 36 38 40 50 60 70 80 90 100 110 120 130 140 150 160 170 180 190 200 210 220 230 240 250 260 270 280
0.33541 0.16771 0.11180 0.08385 0.06708 0.05590 0.04792 0.04193 0.03727 0.03354 0.03049 0.02795 0.02580 0.02396 0.02236 0.02096 0.01973 0.01863 0.01765 0.01677 0.01342 0.01118 0.00958 0.00839 0.00745 0.00671 0.00610 0.00559 0.00516 0.00479 0.00447 0.00419 0.00395 0.00373 0.00353 0.00335 0.00319 0.00305 0.00292 0.00280 0.00268 0.00258 0.00248 0.00240
0.00988 0.02609 0.04526 0.06602 0.08750 0.10915 0.13057 0.15148 0.17173 0.19119 0.20980 0.22752 0.24436 0.26033 0.27544 0.28973 0.30323 0.31599 0.32804 0.33942 0.38763 0.42435 0.45279 0.47521 0.49318 0.50780 0.51987 0.52995 0.53846 0.54573 0.55198 0.55741 0.56215 0.56633 0.57002 0.57331 0.57626 0.57890 0.58129 0.58346 0.58543 0.58723 0.58887 0.59038
2.68328 2.68328 2.68328 2.68328 2.68328 2.68328 2.68328 2.68328 2.68328 2.68328 2.68328 2.68328 2.68328 2.68328 2.68328 2.68328 2.68328 2.68328 2.68328 2.68328 2.68328 2.68328 2.68328 2.68328 2.68328 2.68328 2.68328 2.68328 2.68328 2.68328 2.68328 2.68328 2.68328 2.68328 2.68328 2.68328 2.68328 2.68328 2.68328 2.68328 2.68328 2.68328 2.68328 2.68328
SHELL 181 120 240 0.38 0.35
908.35
1511
1481
1577
1492
1452
0.13
0.22
0.21
0.20
1,079.61
2030
1927
2001
2011
1896
0.15
0.28
0.28
0.26
1,186.23
2596
2424
2495
2630
2409
0.17
0.35
0.37
0.34
1,295.90
3310
3079
3179
3418
3066
0.18
0.44
0.48
0.43
1,386.00
4044
3734
3848
4251
3764
0.19
0.54
0.59
0.52
1,423.16
5934
5667
6406
5590
0.20
0.79
0.89
0.78
1,434.38
6414
6124
6860
6073
0.20
0.85
0.95
0.85
1,554.21
7038
6507
7306
6802
0.22
0.91
1.02
0.95
1,483.12
7260
6722
7246
6747
0.21
0.94
1.01
0.94
- 55 -
Para la probeta nmero 1 se obtuvo una configuracin no axisimtrica, como la mostrada en la figura III-10, se puede observar que difcilmente se alcanzara el punto de bifurcacin, debido a que antes que eso pasara, la estructura se plastificara debido a la aparicin de fallas de tipo local, a causa de las grandes concentraciones de esfuerzos generados (figura III-11).
Figura III-10. Modos de circunferenciales, probeta 1, n=13, m=1
Figura III-11. Concentraciones de esfuerzos, probeta 1
- 56 -
A continuacin se muestra una serie de imgenes y figuras de las 10 diferentes probetas analizadas. Primero se muestra el anlisis hecho para encontrar la carga de bifurcacin (Lanczos), y sus modos de pandeo as como los estados de esfuerzos desarrollados. Al final se proporciona una grfica que muestra la historia de carga vs deformacin hasta llegar al punto lmite, se aprovecha esta grfica para comparar la carga de bifurcacin y la de pandeo clsica obtenidas.
- 57 -
3.1.6.1 Probeta 1
Para la realizacin de este anlisis se usaron los siguientes parmetros:
Dimetro: 300 mm (30 cm) Tipo de elemento: shell 181 Tamao de elemento: 0.7854 Nmero de nodos: 240 Factor de Calidad Q = 2 Imperfeccin generada: Wok = 0.33541 cm Carga de bifurcacin (Lanczos) = 2496 kg Carga crtica del Eurocdigo: 70.96 kg Relacin de la imperfeccin vs espesor Wok/t= 11.18 Para el mtodo del Arc Length (Punto lmite) Carga de referencia 3020 kg Factor de carga de inestabilidad = .105615 Carga introducida 28 * 120 nodos = 3360 kg. Nodo analizado 1945 Carga de inestabilidad = 354.8664 kg.
a) Modo de pandeo axial, m = 1 b) Modo circunferencial, n = 13 Figura III-12 Modos de pandeo, probeta 1 - 58 -
Figura III-13 Esfuerzos de Von Mises para la carga de bifurcacin de la probeta 1, kg/cm2
0.400 P = 2496 kg Pclasica = 7184.51 kg
P/Pcl
0.200
P = 354.85 kg
0.000 0.000
0.005
0.010
0.015
0.020
0.025
0.030
0.035
0.040
0.045
0.050
Desplazamientos (cm)
Carga no lineal Punto de Bifurcacin
Figura III-14 Anlisis no lineal y punto de bifurcacin para la probeta 1
- 59 -
3.1.6.2 Probeta 2
Para la realizacin de este anlisis se usaron los siguientes parmetros
Dimetro : 300 mm (30 cm) Tipo de elemento: shell 181 Tamao de elemento: 0.7854 Nmero de nodos : 240 Factor de Calidad Q = 10 Imperfeccin generada: Wok = 0.06798 cm Carga de bifurcacin (Lanczos) = 1452 kg Carga crtica del Eurocdigo: 628.67 kg Relacin de la imperfeccin vs espesor Wok/t= 2.266 Para el mtodo del Arc Length (Punto lmite) Carga de referencia 1511 kg Factor de carga de inestabilidad = 0.540685 Carga introducida 14 * 120 nodos = 1680 kg. Nodo analizado 2122 Carga de inestabilidad = 908.3508 kg.
a) Modo de pandeo axial, m = 1 b) Modo circunferencial, n= 9 Figura III-15 Modos de pandeo, probeta 2
- 60 -
0.400
Pclasica = 7184.51 kg
P/Pcl
0.200
P = 1452 kg
P = 908.33 kg
0.000 0.000000
0.001000
0.002000
0.003000
0.004000
0.005000
0.006000
0.007000
0.008000
- 61 -
3.1.6.3 Probeta 3
a) Modo de pandeo axial, m = 1
b) Modo circunferencial, n = 12
Figura III-18 Modos de pandeo, probeta 3 - 62 -
0.400
P = 1896 kg
P/Pcl
0.200 P = 1079 kg
0.000 0.0000
0.0005
0.0010
0.0015
0.0020
0.0025
0.0030
0.0035
0.0040
0.0045
0.0050
- 63 -
3.1.6.4 Probeta 4
a) Modo de pandeo, m = 1
Figura III-21 Modos de pandeo, probeta 4
- 64 -
0.400
P = 2409 kg
P/Pcl
0.200
P = 1186.23 kg
0.000 0.00000
0.00100
0.00200
0.00300
0.00400
0.00500
- 65 -
3.1.6.5 Probeta 5
a) Modo de pandeo, m = 1
Figura III-24 Modos de pandeo, probeta 5
- 66 -
P = 3066 kg 0.400
P/Pcl
P = 1296 kg 0.200
0.000 0.0000
0.0040
0.0080
0.0120
- 67 -
3.1.6.6 Probeta 6
Dimetro : 300 mm (30 cm) Tipo de elemento: shell 181 Tamao de elemento: 0.7854 Nmero de nodos : 240 Factor de Calidad :Q = 40 Imperfeccin generada: Wok = 0.01677 cm Carga de bifurcacin (Lanczos) = 3764 kg Carga crtica del Eurocdigo: 2438.55 kg Relacin de la imperfeccin vs espesor Wok/t= 0.559 Para el mtodo del Arc Length (Punto lmite) Carga de referencia 4044 kg Factor de carga de inestabilidad = 0.303947 Carga introducida 38 * 120 nodos = 4560 kg. Nodo analizado 702 Carga de inestabilidad = 1385.99832 kg.
a) Modo de pandeo, m > 1
P = 3764 kg
0.400
P/Pcl
P = 1386 kg 0.200
0.000 0.0000
0.0010
0.0020
0.0030
0.0040
0.0050
0.0060
0.0070
- 69 -
3.1.6.7 Probeta 7
0.800
P = 5590 kg
0.600
P/Pcl
0.400 P = 1423 kg 0.200
0.000 0.0000
0.0001 Desplazam ientos (cm ) Carga no lineal Punto de Bifurcacin
0.0002
- 71 -
3.1.6.8 Probeta 8
Dimetro: 300 mm (30 cm) Tipo de elemento: shell 181 Tamao de elemento: 0.7854 cm Nmero de nodos : 240 Factor de Calidad Q = 150 Imperfeccin generada: Wok = 0.00447 cm Carga de bifurcacin (Lanczos) = 6073 kg Carga crtica del Eurocdigo: 3965.70 kg Relacin de la imperfeccin vs espesor Wok/t= 0.149 Para el mtodo del Arc Length (Punto lmite) Carga de referencia 6414 kg Factor de carga de inestabilidad = 0.202597 Carga introducida 59 * 120 nodos = 7080 kg. Nodo analizado 2056 Carga de inestabilidad = 1439.152 kg.
0.800
P = 6073 kg
0.600 P/Pcl
0.400 P = 1439 kg 0.200
0.000 0.0000
0.0010
0.0020 Desplazam ientos (cm )
0.0030
Carga no lineal
Punto de Bifurcacin
- 73 -
3.1.6.9 Probeta 9
1.000 P = 6802 kg 0.800
P/Pcl
0.600 0.400 P = 1554 kg 0.200
0.000 0.0000
0.0010
0.0020
0.0030
0.0040
0.0050
0.0060
0.0070
Desplazam ientos (cm ) Carga no lineal Punto de Bifurcacin
- 75 -
3.1.6.10
Probeta 10
Figura III-40 Distribucin de esfuerzos para la carga de bifurcacin de la probeta 10, kg/cm2
1.000 P = 6747 kg 0.800
P/Pcl
0.600 0.400 P = 1483 kg 0.200
0.000 0.0000
0.0050
0.0100 Desplazam ientos (cm ) Carga no lineal
0.0150
0.0200
Punto de Bifurcacin
- 77 -
Captulo IV
Experimentacin
4 4.1 4.2 4.2.1 4.2.2 4.2.3 4.2.4 4.2.5 4.3 4.4 4.4.1 4.4.2 4.4.3 4.4.4 4.4.5 4.5 Experimentacin.............................................................................................. 79 Probeta estudiada ........................................................................................... 79 Medicin de los especimenes ......................................................................... 81 Medicin del espcimen 1............................................................................... 83 Medicin del espcimen 2............................................................................... 85 Medicin del espcimen 3............................................................................... 86 Medicin del especimen 4............................................................................... 87 Medicin del espcimen 5............................................................................... 89 Modelacin numrica ...................................................................................... 90 Descripcin de las pruebas ............................................................................. 91 Espcimen 1 ................................................................................................... 91 Espcimen 2 ................................................................................................... 96 Espcimen 3 ................................................................................................. 100 Espcimen 4 ................................................................................................. 103 Espcimen 5 ................................................................................................. 106 Comentarios .................................................................................................. 110
Resumen
Una vez conocidas las tcnicas que se pueden usar para encontrar la carga crtica de pandeo, y haber elegido las dimensiones de fabricacin del especimen, se procede a realizar la fase de experimentacin. De este especimen se obtienen las imperfecciones geomtricas, que son estudiadas a travs de la realizacin de un modelo numrico de anlisis con ayuda de la computadora, y realizar as el anlisis de estabilidad para encontrar la carga critica de pandeo y compararla con la experimental.
4 Experimentacin
Una vez conocidos los parmetros a controlar y el intervalo de valores de la carga crtica para las probetas antes estudiadas, se procede a fabricar 5 especimenes con las dimensiones que a continuacin se detallan, para posteriormente realizar la medicin en estas y as poder realizar un modelo de anlisis que sirva para establecer una comparacin con los especimenes fabricados que sern ensayados a compresin axial en una mquina universal en el laboratorio.
4.1 Probeta estudiada

Se propone un cilindro con las dimensiones mostradas en la figura IV-1, en cuyos extremos se colocaran un par de anillos rigidizantes y se pueda simular (en lo posible) un cilindro con extremos rgidos (ver figura IV-1 y fotografa IV-1).
Las caractersticas geomtricas del modelo son las siguientes: Dimetro, D =: 30 cm; Altura efectiva, h: 30 cm; Espesor, t = 0.03 cm (cal 30) Longitud de la muestra con las soleras: 40 cm Altura de los anillos: 5.08 cm Solera: Acero ASTM A36 de 5.08 cm (2") y espesor de 0.0476 cm (3/16")
Mientras que las caractersticas mecnicas del material son: Material de construccin: Lmina galvanizada cal 30 (0.030 cm). Mdulo de elasticidad, E= 2.1E06 kg/cm2 Mdulo de Poisson, = 0.3 Esfuerzo de fluencia fy = 2530 kg/cm2 Proceso de fabricacin: Se rola la lmina y los anillos, para poder cerrar el cuerpo de lamina cal 30, se proceder a realizar un engargolado.
- 79 -
Planta
3D Elevacin Figura IV-1 Espcimen estudiado
Se fabricaron cinco especmenes similares a los que se ven en la fotografa IV-1.
Fotografa IV-1 Especmenes fabricados
- 80 -
Fotografa IV-2 Se observan las soleras de los extremos para dar rigidez y restringir el giro en las condiciones de frontera
4.2 Medicin de los especmenes

Para llevar a cabo los anlisis comparativos de las pruebas experimentales a realizarse es, como se sabe, de vital importancia conocer las caractersticas geomtricas de los especimenes de prueba y sobre todo las imperfecciones geomtricas iniciales de estos especimenes, las cuales son generadas por los procesos de fabricacin llevados a efecto. As para poder medir fsicamente los especimenes, primeramente se habilit una mesa de nivelacin, que a su vez fuera giratoria para poder realizar mediciones en las paredes de los especimenes cada 10, ver fotografas IV-3 y IV-4.
Fotografa IV-3 Mesa de nivelacin y rotacin habilitada - 81 -
Los especimenes fuern mapeados mediante el rallado de una cudricula que se aprecia en la fotografa IV-4. Para realizar las mediciones, se habilit un indicador de cartula (trusquin) en un calibrador para medir alturas (fotografia IV-5 y IV-6), de esta manera se midieron las imperfecciones geomtricas existentes a lo largo de los meridianos de los especimenes, y conforme se fuera rotando la mesa, se tendra un mapeo de todo el espcimen en rotaciones a cada 10.
Fotografa IV-4 Espcimen de prueba sobre la mesa de nivelacin, observe el rallado del espcimen a cada 10
Fotografa IV-5 Sistema de medicin del espcimen habilitado, el calibrador de alturas solo tiene desplazamientos verticales - 82 -
Fotografa IV-6 Uso del indicador de cartula sobre el calibrador de alturas, el indicador de cartula mide las imperfecciones existentes.
Las figuras IV-5 y IV-6 muestran el dispositivo de medicin habilitado para conocer las imperfecciones geomtricas iniciales de los especimenes de prueba, se observa el proceso de medicin que se sigui para conocer las imperfecciones geomtricas iniciales a travs de los meridianos y paralelos, tomando medidas a cada 25 mm a lo largo de los meridianos y a cada 10 a lo largo de los paralelos.
4.2.1 Medicin del espcimen 1

La mxima imperfeccin de este espcimen se da en el paralelo 125 y el meridiano 290 y corresponde a un valor -0.853 mm, el signo negativo indica que la imperfeccin se gener en direccin radial hacia el interior del cilindro. En la figura IV-2 se puede observar que el espcimen 1 tiene imperfecciones geomtricas en toda su superficie y son graduales o de comportamiento senoidal (ver figura IV-3). Mientras que a lo largo de la altura (meridianos), la mxima imperfeccin tiene un comportamiento gradual, tambin de tipo senoidal (ver figura IV-4).
- 83 -
0.600 0.400 Imperfecciones(mm) 0.200 0.000 -0.200 -0.400 -0.600 -0.800 -1.000 -0 50 100 150 200 250 300 0 Altura, paralelos (mm) 80 40 120 Meridianos, grados 200 160 240 320 280
Figura IV-2 Superficie de las imperfecciones geomtricas, espcimen 1

Espcimen 1 (paralelo 125)
0.800
Imperfecciones (mm)
0.600 0.400 0.200 0.000 -0.200 0 -0.400 -0.600 -0.800 -1.000 50 100 150 200 250 300 350 400
Grados
Figura IV-3 Mxima imperfeccin generada en el espcimen 1, paralelo 125
0.000
Espcimen 1 (meridiano 290)

50 100 150 200 250 300
Imperfecciones (mm)
-0.100 0 -0.200 -0.300 -0.400 -0.500 -0.600 -0.700 -0.800 -0.900
Altura (mm)
Figura IV-4 Mxima imperfeccin generada en el espcimen 1, meridiano 290 - 84 -

La mxima imperfeccin de este espcimen se da en el paralelo 250 y el meridiano 180 y corresponde a un valor 2.498 mm, aproximadamente 8 veces el espesor. Esta imperfeccin corresponde quizs a algn tipo de imperfeccin generada por un golpe de tipo puntual, esto se puede deducir debido a su amplitud correspondiente, tal y como se muestra en la figura IV-5. A pesar de esto el espcimen 2 es el que presenta mayor uniformidad en su superficie y de esto se pudiera intuir que es el espcimen que se retrasar ms en la generacin de su carga de inestabilidad. En las figuras IV-6 y IV-7 se puede observar que las imperfecciones pudieran corresponder a un comportamiento puntual, como se presenta en las funciones de tipo impulso.
2.500 2.000 Imperfecciones(mm) 1.500 1.000 320 0.500 240 0.000 -0.500 -1.000 0 50 100 150 200 250 300 0 Altura, paralelos (mm) 80 40 160 120 200 Meridianos, grados 280
Figura IV-5 Medicin de las imperfecciones geomtricas, espcimen 2

3.000
Imperfecciones (mm)
2.500 2.000 1.500 1.000 0.500 0.000 -0.500 0 50 100 150 200 250 300 350 400
Grados
Figura IV-6 Mxima imperfeccin generada en el espcimen 2, paralelo 250 - 85 -
3.000
Imperfecciones (mm)
2.500 2.000 1.500 1.000 0.500 0.000 -0.500 0 50 100 150 200 250 300
Altura (mm)
Figura IV-7 Mxima imperfeccin generada en el espcimen 2, meridiano 180

La mxima imperfeccin de este espcimen se da en el paralelo 175 y el meridiano 20 y corresponde a un valor 1.007 mm. Como se puede apreciar en las figuras IV-8, IV-9 y IV-10 las imperfecciones presentadas son graduales y del tipo senoidal.
1.200 1.000 Imperfecciones(mm) 0.800 0.600 0.400 0.200 0.000 -0.200 -0.400 0 50 100 150 200 250 300 0 Altura, paralelos (mm) 80 40 120 200 160 Meridianos, grados 240 320 280
Figura IV-8 Medicin de las imperfecciones geomtricas, espcimen 3 - 86 -

1.200
Imperfecciones (mm)
1.000 0.800 0.600 0.400 0.200 0.000 -0.200 0 -0.400 -0.600 50 100 150 200 250 300 350 400
Grados
1.200
Imperfecciones (mm)
1.000 0.800 0.600 0.400 0.200 0.000 -0.200 0 50 100 150 200 250 300
Altura (mm)

La mxima imperfeccin de este espcimen se da en el paralelo 50 y el meridiano 290 y corresponde a un valor 1.122 mm. Este espcimen presenta una combinacin de imperfecciones, la mxima que corresponde a una falla puntual (tipo impulso) y el resto que tienen un comportamiento gradual del tipo senoidal, tal y como se muestra en las figuras IV-11, IV-12 y IV-13. En la figura IV-13 se puede apreciar que la imperfeccin es transitoria, es decir primero va del interior al exterior del cilindro, para posteriormente cambiar de direccin, hacia el interior del especimen. - 87 -
0.600 0.400 Imperfecciones(mm) 0.200 0.000 -0.200 -0.400 -0.600 -0.800 -1.000 -1.200 0 50 100 150 200 250 300 0 Altura, paralelos (mm) 80 40 120 160 240 200 Meridianos, grados 280 320

0.600
Imperfecciones (mm)
0.400 0.200 0.000 -0.200 0 -0.400 -0.600 -0.800 -1.000 -1.200 50 100 150 200 250 300 350 400
Grados
0.400
Imperfecciones (mm)
0.200 0.000 -0.200 -0.400 -0.600 -0.800 -1.000 -1.200 0 50 100 150 200 250 300
Altura (mm)
Figura IV-13 Mxima imperfeccin generada en el espcimen 3, meridiano 290 - 88 -

La mxima imperfeccin de este espcimen se da en el paralelo 75 y el meridiano 60 y corresponde a un valor 0.398 mm, ligeramente mayor al espesor. Es el especimen que presenta las imperfecciones geomtricas menores y son del tipo cosenoidal, como se puede apreciar en la figuras IV-14 a la IV-16.
0.400 0.300 Imperfecciones(cm) 0.200 0.100 0.000 -0.100 -0.200 -0.300 -0.400 0 50 100 150 200 250 300 0 Altura, paralelos (mm) 80 40 120 200 160 Meridianos, grados 280 240 320

0.500
Imperfecciones (mm)
0.400 0.300 0.200 0.100 0.000 -0.100 0 -0.200 -0.300 50 100 150 200 250 300 350 400
Grados
- 89 -
0.400
Imperfecciones (mm)
0.300 0.200 0.100 0.000 -0.100 0 -0.200 -0.300 -0.400 -0.500 50 100 150 200 250 300
Altura (mm)
4.3 Modelacin numrica

Se hace uso del mtodo del elemento finito mediante la aplicacin del software ANSYS 8.1; primeramente se realizar un anlisis de pandeo para encontrar la carga crtica del espcimen, posteriormente se realizar un anlisis no lineal mediante la tcnica del ArcLength de ANSYS.
El anlisis de estabilidad por pandeo se realizar con la tcnica que se ha venido manejando, es decir el mtodo Lanczos. Sin embargo, como se podr observar en el resumen de resultados, se obtienen cargas crticas bastante altas respecto a las cargas reales de inestabilidad que se presentan en los ensayos experimentales, esto debido a que con el mtodo Lanczos se est encontrando el punto de bifurcacin de la estructura, que es un comportamiento terico lineal, tal y como se explic en la seccin 3.1.2. Por esta razn, se hace indispensable realizar un anlisis no lineal, (para este caso se usar la tcnica del Arclength) que representa un anlisis por punto limite de un cascarn cilndrico con imperfecciones geomtricas. Primeramente se encontrar la carga crtica de bifurcacin del espcimen como valor de referencia, para posteriormente introducir esta misma carga para realizar el anlisis no lineal. El mtodo del Arclength tiene la caracterstica de poder encontrar el punto de inestabilidad por pandeo de la estructura estudiada, mediante la aplicacin gradual de incrementos de carga, hasta que el proceso numrico detecta una singularidad en la matriz de rigidez del sistema, lo que se origina cuando en algn elemento del espcimen se tiene una - 90 -
rigidez cero o muy prxima a cero, ocasionando con esto un punto de inestabilidad, siendo el punto buscado. Al finalizar el proceso numrico, el software lo que proporciona es un factor de carga en el cual se gener la singularidad en la matriz de rigidez del sistema y para traducirlo a carga de inestabilidad por pandeo slo hay que multiplicar dicho factor por la carga inicial introducida.
4.4 Descripcin de las pruebas

Se probaron 5 especimenes con las caractersticas geomtricas mostradas en la figura IV-1 (fotografa IV-1). Para cada espcimen, su descripcin es presentada a continuacin.
4.4.1 Espcimen 1
El cascarn cilndrico del espcimen 1 fue sometido a una prueba de compresin axial en una prensa universal y se tuvo que adaptar un par de mordazas (fotografa IV-7) para trasmitir en lo posible una carga perpendicular distribuida en el permetro de la probeta. Se pudo apreciar que el punto de inestabilidad por pandeo se alcanz con una carga de 350 kg, como se puede ver en la fotografa IV-8. El tipo de pandeo corresponde a una falla de tipo local y se vi caracterizada por la aparicin repentina del modo de pandeo snap through. La forma de inestabilidad presentada corresponde a la que se presenta por punto lmite y que es explicada en la seccin 2.1.4, adems se presenta en un modo no axisimtrico.
Fotografa IV-7 Espcimen 1 montado sobre la mquina universal con las mordazas de prueba - 91 -
Fotografa IV-8 Espcimen 1 modo de pandeo
La carga de inestabilidad por pandeo es del orden del 4.8% de la clsica y del 26% respecto del anlisis no lineal (punto limite), ver figura IV-21, que es ms pequea de lo que se esperaba en el modelo numrico, ya que la aplicacin de la carga se vi afectada por una deformacin en las mordazas fabricadas, es decir la transferencia de carga no fue perfectamente perpendicular al espcimen. Dado que se observ que las mordazas no se apoyaban perfectamente sobre el permetro de la probeta, se aprovech esta experiencia para encontrar la carga critica por pandeo del espcimen bajo estas condiciones, para posteriormente retirar la carga aplicada, encontrando que la deformacin generada por la aplicacin de la carga en el cascarn cilndrico, regresaba histerticamente, es decir que a pesar que la carga se retirar al 100% se mantena una deformacin permanente en el punto de aparicin de la falla, ver fotografa IV-9. Lo anterior a pesar de que la carga de inestabilidad por pandeo de 350 kg encontrada corresponde a un esfuerzo actuante de tipo elstico.
Fotografa IV-9 Espcimen 1, deformacin permanente, despus de retirar la carga al 100% - 92 -
En el video anexo a este documento se puede apreciar con mayor detalle lo expresado anteriormente. Se puede observar que queda una huella que representa la zona donde comenz a generarse el primer modo de pandeo, por lo que despus de volver a aplicar carga, se puede ver que en esa zona de induccin se vuelve a generar el modo de pandeo, pero con la aplicacin de una carga menor y del orden de 200 kg.
La deformacin de las mordazas en el espcimen 1 fue ocasionada, debido a que el calor generado por el proceso de soldadura en la fabricacin de dichas mordazas no fue debidamente controlado. Esta deformacin de las mordazas, gener que la carga no se trasmitiera uniformemente en el permetro del cilindro, concluyendo as que s la carga comienza a tener excentricidad en su punto de aplicacin, la carga crtica se ve fuertemente afectada. Para llegar a la conclusin anterior, fue necesario realizar un anlisis numrico mediante ANSYS 8.1 del espcimen. En este anlisis se compar el efecto de la excentricidad en la aplicacin de la carga respecto a la carga crtica por pandeo del espcimen. Como se puede apreciar en la figura IV-17, la transferencia de la carga debe realizarse a travs de todo el permetro del espcimen, que para el caso numrico corresponde a la totalidad de nodos en el extremo superior (108 nodos). El anlisis consisti en ir aplicando carga a los nodos del extremo superior del espcimen numrico, desde un solo nodo, hasta el 100% de los nodos.
Z Y
Figura IV-17 Correcta aplicacin de la carga
- 93 -
Como se puede ver en la figura IV-18, el comportamiento en la aparicin del punto de bifurcacin y/o carga de inestabilidad es muy sensible a la aplicacin de la carga. Y es muy notorio si nos posicionamos en la aplicacin del 100% de los nodos, donde esperaramos una reduccin en el valor de la carga crtica clsica vs carga critica desarrollada (carga de bifurcacin, en este caso) del orden del 30 al 40% debido a las imperfecciones geomtricas existentes en el espcimen, sin embargo si vamos disminuyendo la cantidad de nodos donde se aplica la carga (carga excntrica), se puede apreciar reduccin de carga hasta del orden del 96.8%.
1.000 0.900 0.800 0.700 0.600
Pcr/Pcl
0.500 0.400 0.300 0.200 0.100 0.000 0% 10% 20% 30% 40% 50% 60% 70% 80% 90% 100%
% No. de nodos
Figura IV-18 Grfico que muestra la sensibilidad de no aplicar la carga en todos los nodos del modelo numrico, se normaliza respecto al dato de carga crtica numrico
La figura IV-18 muestra la sensibilidad de la carga crtica de pandeo por las condiciones de frontera, es decir si el espcimen no logra estar perfectamente apoyado en el 100% de su permetro, la carga crtica se ve afectada, tal y como se puede ver en el grfico IV-18, por ejemplo si por alguna razn el permetro del especimen solo logra apoyarse en el 70% de su totalidad, solo se logra una carga crtica de 0.33 respecto a la carga crtica clsica
- 94 -
4.4.1.1 Modelacin numrica, espcimen 1

Tal y como se aprecia en la fotografa IV-8, se esperara un pandeo del tipo local, y que es replicado en la figura IV-19. Adems que la distribucin de esfuerzos debera ser del tipo elstico lineal, el cual es verificado en la figura IV-20.
3300 (AVG) 0 .92984 .92984 13:04
MN
MX
MX
Y Z
Y
MN
.206631 .103316 .309947
.413262 .516578
.619894 .723209
.82652
a) Modo de pandeo axial
b) Modo circunferencial
Figura IV-19 Modos de pandeo, espcimen 1
.363366 (AVG) =.020425 =235.087 =577.062
13:32:38
MX
Y Z 235.087 273.084 311.081 349.078 X 463.07 425.073 501.067 539.064 577.062
387.076
Figura IV-20 Distribucin de esfuerzos para la carga de bifurcacin del espcimen 1
En el grfico de la figura IV-21 se puede observar que se alcanza el punto lmite con una carga de 1311.64 kg que corresponde al 18% del valor de la carga critica clsica, - 95 -
mientras que la carga de pandeo por experimentacin se desarrolla en un 4.8% de la clsica. La carga de 3299.63 kg, corresponde a la carga de bifurcacin del espcimen y slo es vlida como un valor de referencia, ya que en la prctica difcilmente se obtendra esta carga. La distribucin de esfuerzos, tal y como lo indica la figura IV-20 son del tipo elstico, es decir por debajo del punto de fluencia del material (2530 kg/cm2).
1.000
P = 7184.51 kg
0.800
0.600
P/Pcl
P = 3299.63 kg 0.400
0.200
P = 1311.64 kg P = 350 kg
0.000 0.000
0.001
0.002
0.003
0.004
0.005
0.006
0.007
0.008
0.009
Desplazamientos
Anlisis no lineal Carga critica clsica Punto de Bifurcacin Carga experimental
Figura IV-21 Anlisis no lineal y punto de bifurcacin para el espcimen 1
4.4.2 Espcimen 2
La probeta 2 fue sometida a una prueba de compresin axial en una maquina universal como el ensaye anterior, se construyo un nuevo par de mordazas. Las dimensiones de las nuevas mordazas son de 400 X 400 X 25.4 mm. La finalidad de estas nuevas mordazas es evitar deformaciones en el momento de la aplicacin de la carga y que pudieran tener influencia en las condiciones de frontera del espcimen y que afecten en los resultados obtenidos. En la fotografa IV-10 se puede observar el espcimen 2 con el nuevo diseo de mordazas.
- 96 -
Fotografa IV-10 Espcimen 2 montado sobre la mquina universal con las nuevas mordazas de prueba
Al colocar esta probeta en la mquina universal se tuvieron que adicionar unas pequeas lainas con la finalidad de lograr un apoyo uniforme sobre la superficie de la mordaza, y aplicar perpendicularmente la carga, evitando el mismo efecto que en el espcimen anterior.
Con las nuevas mordazas no se observ ningn tipo de deformacin al momento de aplicar la carga de prueba, y por consiguiente la carga vertical fue transmitida en su totalidad sobre el permetro del espcimen, logrando as el efecto deseado y adems evitar alguna excentricidad. El valor de la carga de inestabilidad por pandeo del espcimen fue de 1922 Kg.
Notse que el valor de la carga de inestabilidad por pandeo de esta prueba es mayor en comparacin del primer espcimen. Esto se logr debido a que las condiciones de frontera del espcimen son diferentes a las de la prueba anterior. Tambin se logr que la carga vertical aplicada fuera completamente perpendicular al permetro de la probeta.
El tipo de pandeo presentado tambin corresponde a uno de tipo no axisimtrico por punto lmite y se caracteriza por la aparicin repentina del modo de pandeo local, como se puede ver en la fotografa IV-11. - 97 -
Fotografa IV-11 Espcimen 2, modo de pandeo

Se puede apreciar en la figura IV-22 que los modos de pandeo alcanzados son muy parecidos geomtricamente a los que se presentan en las pruebas experimentales.
En el grfico de la figura IV-24 se puede observar que se alcanza el punto lmite con una carga de 1426.31 kg que corresponde al 19.85% del valor de la carga critica clsica, mientras que la carga de pandeo por experimentacin se desarrolla en un 26.75% de la clsica. La carga de 4521.5 kg, corresponde a la carga de bifurcacin del espcimen y slo es vlida como un valor de referencia, ya que en la prctica difcilmente se obtendra esta carga. La distribucin de esfuerzos, tal y como lo indica la figura IV-23 son del tipo elstico, es decir por debajo del punto de fluencia del material (2530 kg/cm2).
- 98 -
OCT 19 11:1 (AVG) 3 3
G)
MX
MX
Y
MN
(AVG) =.037011 =231.587 =598.169
MN MX
Y XZ
231.587 272.318
313.05 353.781
394.512 435.244
475.975 516.706
557.438 598
- 99 -
1.000
P = 7184.51 kg
0.800
0.600
P = 4521.5 kg
P/Pcl
0.400 P = 1922 kg 0.200 P = 1426.31 kg 0.000 0.000
0.005
0.010
0.015
0.020
0.025
0.030
Desplazamientos
4.4.3 Espcimen 3
El espcimen 3, fue sometido a una prueba de compresin axial en una mquina universal. La carga de inestabilidad de 524 kg as como el modo de pandeo presentado se puede ver en la fotografa IV-12. El modo de pandeo para este caso se presenta muy cerca de los anillos rigidizantes del espcimen. Sin embargo se tiene cierta duda que las lainas colocadas no hayan sido suficientes para dar un apoyo axial uniforme al espcimen, dado que exista mucha variacin de paralelismo entre las caras del cilindro.
Fotografa IV-12 Especimen 3, modo de pandeo - 100 -
La razn de no lograr paralelismo en las caras se dio porque todos los especimenes fueron fabricados de forma independiente, por lo que dichos especmenes no fueron 100% idnticos. Por lo que al colocar lainas de nivelacin y dado que haba poco paralelismo en las caras del cilindro, no se logro un buen apoyo.
En esta prueba se decidi realizar varios ciclos para poder obtener diferentes cargas de pandeo, y adems, para ver el comportamiento del espcimen. El procedimiento usado para la repeticin de la carga, consisti en retirar la carga a 0 una vez que el cilindro alcanzaba su carga de inestabilidad por pandeo, este procedimiento fue repetido hasta 10 veces.
Las dems cargas de inestabilidad obtenidas despus de la primera carga de inestabilidad se reducen notoriamente como se puede ver en el grfico de la figura IV25. En el momento de retirar la carga las deformaciones generadas no regresan por completo a su posicin inicial, generando una deformacin permanente en esa zona de la probeta.
550 500
524
Ciclos de carga de inestabilidad

Carga (kg) 450 400 350 300 0 2 4 6 Nmero de ciclos 8 10 12
380 369 359 348 350 353 359 359 367
374
Figura IV-25 Grfico de las cargas de inestabilidad del espcimen 3
- 101 -

En la figura IV-26 se puede constatar que el modo de pandeo obtenido es muy similar al desarrollado en la prueba fsica.
En el grfico de la figura IV-28 se puede observar que se alcanza el punto lmite con una carga de 772.55 kg que corresponde al 10.75% del valor de la carga critica clsica, mientras que la carga de pandeo por experimentacin se desarrolla en un 7.29% de la clsica. La carga de 2180.3 kg, corresponde a la carga de bifurcacin del especimen y slo es vlida como un valor de referencia, ya que en la prctica difcilmente se obtendra esta carga. La distribucin de esfuerzos, tal y como lo indica la figura IV-27 son del tipo elstico, es decir por debajo del punto de fluencia del material (2530 kg/cm2).
MX
MX
MN
Y Z
a) Modo de pandeo b) Modo circunferencial Figura IV-26 Modos de pandeo, espcimen 3
- 102 -
B =40 ME=.311011 QV (AVG) X =.166136 N =54.586 X =833.943
MN MX
OCT 28 2004 09:49:50
Y XZ
54.586 141.181
227.777 314.372
400.967 487.562
574.158 660.753
747.348 833.943

1.000 P = 7184.51 kg
0.800
0.600
P/Pcl
0.400 P = 2180.3 kg 0.200 P = 772.55 kg P = 524 kg 0.000 0.000 0.020 0.040 0.060 0.080 0.100 0.120 0.140
Desplazamientos
4.4.4 Espcimen 4
El especimen 4 fue sometido al igual que los anteriores a una prueba de compresin axial en una mquina universal. Al igual que la probeta 2 y la probeta 3, se le tuvieron que colocar lainas para que la carga se aplicara uniformemente. Las lainas se colocaron - 103 -
a travs de la mitad de la circunferencia, debido a que exista mucha variacin de paralelismo entre las caras del espcimen. Las mencionadas lainas se dispusieron perimetralmente entre el espcimen y las mordazas para trasmitir perpendicular y axialmente la carga aplicada. La carga de inestabilidad encontrada fue de 411 kg y el modo de pandeo multimodal se puede apreciar en la fotografa IV-13.
Fotografa IV-13 Espcimen 4, modo de pandeo multimodal

En las figura IV-29 se puede verificar que el modo de pandeo obtenido es muy similar al desarrollado en la prueba fsica, siendo este del tipo diamantado.
En la curva de la figura IV-31 se puede observar que se alcanza el punto lmite con una carga de 1297.70 kg que corresponde al 18.10% del valor de la carga critica clsica, mientras que la carga de pandeo por experimentacin se desarrolla en un 5.7% de la clsica. La carga de 2641.64 kg, corresponde a la carga de bifurcacin del espcimen y slo es vlida como un valor de referencia, ya que en la prctica difcilmente se obtendra esta carga. La distribucin de esfuerzos, tal y como lo indica la figura IV-30 son del tipo elstico, es decir por debajo del punto de fluencia del material (2530 kg/cm2).
- 104 -
MX
MN
Y Z
MX
a) Modo de pandeo axia
P=1 =37 E=.496034 V (AVG) =.035071 =287.433 =569.233
NOV 4 2004 17:31:35
MX
Y
MN X
Z 475.3 443.989 506.611 537.922 569.233
287.433 318.744
350.055 381.366
412.677
- 105 -
1.000
P = 7184.51 kg
0.800
0.600
P/Pcl
0.400 P = 2641.64 kg P = 1297.70 kg
0.200
P = 411 kg 0.000 -0.001 0.000 0.001 0.002 0.003 0.004 0.005 0.006 0.007 0.008 0.009
Desplazamientos
En la figura IV-31 se puede apreciar en la lnea que identifica el anlisis no lineal, que existe un incremento gradual en la deformacin presentada conforme aumenta la aplicacin de carga. Este incremento es gradual hasta alrededor del 18% de la carga crtica de pandeo clsica, despus de esto se puede apreciar un brinco repentino de la configuracin deformada a otra (snap through), y se puede apreciar que el snap through esta muy bien indicado, ya que se pasa de una deformacin hacia dentro del especimen, hacia otro estado de deformacin fuera del espcimen.
4.4.5 Especimen 5
El espcimen 5 fue sometido a una prueba de compresin axial. En esta prueba tambin se opt por realizar algunos ciclos de carga para ver el comportamiento de la probeta. Al igual que a los especmenes 2, 3 y 4, se le tuvieron que colocar lainas para que la carga se aplicara uniformemente.
El modo de pandeo encontrado para la carga de inestabilidad por pandeo de 1800 kg es el mostrado en la fotografa IV-14.
- 106 -
Fotografa IV-14 Especimen 5, modo de pandeo
Cabe mencionar que al momento de colocar la probeta en la mquina universal se observ que este espcimen era el que presentaba menor perpendicularidad con las mordazas, en consecuencia se le tuvieron que colocar ms lainas que a los especmenes anteriores. Por lo que se procuro en lo posible dejar la menor cantidad de separacin entre las mordazas y el cilindro.
En la figura IV-32 se muestra que la primera carga de inestabilidad por pandeo se presenta a los 1800 kg de carga. Se observa que entre la primera y la segunda carga de inestabilidad hay una variacin, esto se debe a que la probeta sufre una deformacin permanente despus de presentarse la primera carga de inestabilidad por pandeo. Este espcimen muestra un comportamiento diferente al espcimen 3, ya que se puede apreciar un aumento en la obtencin de la carga de pandeo en los ciclos de carga 3, 4 y 5 (figura IV-32). Este fenmeno puede estar ligado a la capacidad de pandeo poscrtica del cilindro y evitar as un colapso prematuro, es decir es un mecanismo de sobre vivencia del espcimen. En el ciclo de carga nmero 6 la aparicin del modo de pandeo no se registro en la misma zona que los ciclos 1 al 5, esto quizs se debe a que se gener una plastificacin local y un endurecimiento del material en los primeros cinco modos de pandeo.
- 107 -
1900 1700 Carga (kg) 1500 1300 1100 900 700 500 0
1800
Ciclos de carga de inestabilidad

1712 1500 1535 1374
700 750 750 709 747
10
12
Nmero de ciclos
Figura IV-32 Grfico de las cargas de inestabilidad de la probeta 5

En la figura IV-33 se puede constatar que el modo de pandeo obtenido es muy similar al desarrollado en la prueba fsica siendo este del tipo diamantado.
En el grfico de la figura IV-35 se puede observar que se alcanza el punto lmite con una carga de 1389.26 kg que corresponde al 19.34% del valor de la carga critica clsica, mientras que la carga de pandeo experimental se desarrolla en un 25.1% de la clsica. La carga de 3721 kg, corresponde a la carga de bifurcacin del especimen y slo es vlida como un valor de referencia, ya que en estructuras reales difcilmente se obtendra esta carga. La distribucin de esfuerzos, tal y como lo indica la figura IV-34 son del tipo elstico, es decir por debajo del punto de fluencia del material (2530 kg/cm2), ver figura IV-34.
- 108 -
MX
MN
Y Z
Y
MX
SEQV (AVG) DMX =.032503 SMN =216.483 SMX =543.554
MX MN
Y Z 216.483 252.824 289.165 325.507 361.848 398.189 X 434.53 470.871 507.212 543
- 109 -
1.000
P = 7184.51 kg
0.800
0.600
P/Pcl
P = 3721 kg 0.400 P = 1800 kg 0.200 P = 1389.26 kg 0.000 0.000
0.001
0.001
0.002
0.002
0.003
0.003
0.004
Desplazamientos
4.5 Comentarios
En la tabla IV-1 se muestra un resumen de los resultados obtenidos, y en la tabla IV-2 se presentan los valores de las cargas obtenidos respecto a la carga crtica clsica. Los valores presentados en las tablas IV-1 y IV-2 son los correspondientes a los especmenes del 1 al 5, con sus imperfecciones geomtricas iniciales, en estas tablas se proporciona el valor de la carga critica clsica por pandeo de esos cilindros sin considerar las imperfecciones geomtricas (7184.51 kg) como valor de referencia y con el cual se han venido presentando todos los grficos normalizados de esta tesis.
Tabla IV-1 Resultados comparativos de varias metodologas

Especimen 1 2 3 4 5 Pcr, clasica (kg) 7184.51 7184.51 7184.51 7184.51 7184.51 Carga Metodo Carga de bifurcacin Carga por punto experimental Eurocdigo Baker & MEF (Ansys), metodo limite (kg) (kg) Kovalesky Lanczos 3299.6259 1311.6384 350 466.57 2312 4521.4858 1426.30524 1922 130.43 2312 2180.268 772.54884 524 231.45 2312 2641.6535 1297.698192 411 324.24 2312 3721.0018 1389.261456 1800 1154.95 2312
- 110 -
Tabla IV-2 Cargas criticas obtenidas vs carga clsica esperada

Especi men 1 2 3 4 5 Metodo del Arclength Carga de Carga Carga de bifurcacin Factor de Carga inestabilidad experimental MEF (ANSYS), metodo carga de introducida (kg) por punto (kg) LYNCZOS inestabilidad 3299.626 0.357 3672 0.183 0.049 4521.486 0.281 5076 0.199 0.268 2180.268 0.311 2484 0.108 0.073 2641.654 0.496 2616 0.181 0.057 3721.002 0.339 4104 0.193 0.251 Metodo Baker Metodo Baker & Kovalesky & Kovalesky/3 0.322 0.322 0.322 0.322 0.322 0.107 0.107 0.107 0.107 0.107
Pcr, clasica (kg) 7184.510 7184.510 7184.510 7184.510 7184.510
Eurocdigo 0.065 0.018 0.032 0.045 0.161
Es importante hacer una comparacin adicional contra lo establecido en el Eurocdigo y el mtodo propuesto de Baker & Kovalesky y as tener un criterio adicional sobre el comportamiento de estos especimenes. Para el mtodo de Baker & Kovalesky existe cierta indiferencia en la magnitud de la imperfeccin geomtrica y nicamente las considera de manera global, la carga se calcula segn lo establecido en las ecuaciones II-17 y II-18 del captulo 2:
:=
1 16
R t
:= 1 0.901 1 e cr := E
)
t cr = 817.825 Esfuerzo crtico
2 3 ( 1 )
Pcr := cr A
Pcr = 2.312 10
3 kg, carga crtica
En la prctica para disear un cilindro sometido a carga axial y usando el mtodo de Baker, es muy comn el uso de factores de seguridad de 3, es decir que segn el diseo clsico y usando la metodologa de Baker la carga mxima a la que se deberan someter los especimenes en estudio seria de 2312 kg / 3 = 770.6 kg, que representa hasta un 10.72% de la carga crtica de pandeo clsica (7184.51 kg).
As de las curvas de la figura IV-36, se puede observar que la metodologa expuesta por Baker&Kovalesky modificada con un factor de seguridad de 3, parece dar buenos resultados para el caso donde se controlan las condiciones de frontera, tales como aplicacin de la carga y de empotramiento, ya que los especimenes 2 y 5 fueron en los que mejor se pudo controlar estos aspectos. Sin embargo en la metodologa propuesta por el Eurocdigo, para valores menores de Q=16 (ver tabla IV-3) las cargas crticas - 111 -
propuestas son muy conservadoras y pueden ser poco econmicas. Por lo que la aplicacin del Eurocdigo se reduce a valores de factor de calidad Q entre 16 y 40, como previamente se haba establecido en el captulo 3.
0.400
0.300
Pcr/Pcl
0.200
0.100
0.000 0 1 2 3 especim enes Carga de inestabilidad por punto lmite(kg) Carga experimental (kg) Eurocdigo Metodo Baker & Kovalesky Metodo Baker&Kovalesky/3 4 5 6
Figura IV-36 Comparacin de las cargas de inestabilidad de los 5 especimenes
En la tabla IV-3 se puede ver que los especimenes del 1 al 4 corresponden a factores de calidad Q menores de 16 (entre 4 y 8), mientras que el especimen 5 tiene un factor de calidad Q de 16.86.
Dado que el factor de calidad Q del especimen 2 es de 2.69 (espcimen con mayor grado de imperfeccin) se esperara una carga crtica del orden del 6% como los especmenes 1, 3 y 4, sin embargo no fue as, se obtuvo una carga del 26.8% respecto de la clsica. - 112 -
El parmetro de imperfeccin (U) se relaciona con el factor de calidad Q, tal y como se explico en la seccin 2.2.2.2.
Tabla IV-3 Carga critica de los especimenes usando el Eurocdigo

Imperfecciones del especimen mx. mx. imperfeccin plg. imperfeccin cm. 0.0336 0.0983 0.0590 0.0442 0.0157 0.0853 0.2498 0.1498 0.1122 0.0398
probeta 1 2 3 4 5
Pcr w ok 466.5728 130.4382 231.4502 324.2479 1154.9571
Uox
0.032 0.093 0.056 0.042 0.015
<16 7.86 <16 2.69 <16 4.48 <16 5.98 16 16.86
0.400
0.300
Pcr/Pcl
0.200
0.100
0.000 0 1 2 3 4 5 6 Carga de inestabilidad por punto lmite(kg) especim enes Carga experimental (kg) Eurocdigo Metodo Baker & Kovalesky Metodo Baker&Kovalesky/3 Carga de inestabilidad por punto lmite(kg)(media) Carga experimental (kg) (media) Eurocdigo (media)
Figura IV-37 Comparacin de las cargas de inestabilidad de los 5 especimenes
Para tener unas conclusiones ms contundentes, se determina el valor medio de los mtodos elegidos, por lo que en el grfico IV-37 se puede apreciar que la metodologa de Baker&Kovalesky es la que mejor aproximacin tiene con respecto a las cargas experimentales obtenidas. - 113 -
Discusin de resultados
El no lograr una buena distribucin de carga aplicada en el permetro de un silo o un cascarn cilndrico, ocasiona una disminucin en el valor de la carga crtica de inestabilidad por pandeo esperada. Esto se pudo comprobar en los especimenes ensayados 1, 3 y 4, en los cuales no se logr un buen apoyo axial, provocando con esto una disminucin en la carga de inestabilidad. En el espcimen nmero 3 se logro un mediano apoyo perimetral, sin embargo se obtuvo un valor bajo de carga, que se verifico mediante un anlisis no lineal (punto limite) cuyo valor es cercano al experimental, considerando las imperfecciones geomtricas existentes que presentaba. Lo anteriormente descrito da pauta a recomendar que en silos con faldn apoyados sobre una cimentacin, se debe proveer de un sellador (grout) para asegurar que el permetro del silo se apoye perfectamente y con ello disminuir la posibilidad de falla por inestabilidad en el faldn.
En los especimenes 2 y 5 se logr un apoyo total en su permetro, obtenindose valores de carga de alrededor del 25% respecto al valor de la carga de pandeo terica, sin duda alguna debido a las imperfecciones geomtricas existentes en estos especimenes.
Como se puede observar en el grfico de la figura IV-36 en combinacin con la tabla IV3, un anlisis no lineal (punto lmite) da buen resultado para diseo, siempre y cuando se controlen las condiciones de apoyo del cascarn, es decir apoyo perimetral uniforme.
Debido a que es difcil realizar mediciones en los faldones de los silos de sus imperfecciones geomtricas. Se debe partir de un diseo donde se supone la existencia de dichas imperfecciones, aplicando la metodologa de Baker&Kovalesky afectada por un factor de seguridad de 3.
Si se tiene la seguridad que los silos fabricados van a tener factores de calidad entre 16 y 40, la aplicacin del Eurocdigo es aplicable.
114
Conclusiones
Los valores encontrados en esta tesis quedan restringidos en su aplicacin para relaciones L/r = 2 y una relacin r/t = 500, sin embargo la metodologa empleada es de aplicacin general y puede aplicarse para las diferentes relaciones L/r y r/t necesarias.
Las cargas de inestabilidad experimentales por pandeo encontradas en todos los especmenes fueron del tipo elstico, y as lo demuestran los ensayos hechos a los especimenes 1, 3 y 5, donde al retirar la carga se pudo apreciar que los modos de pandeo inducidos regresaban a su estado inicial, dejando slo una pequea deformacin permanente.
Es importante hacer notar que al retirar la carga de un espcimen pandeado, y al regresar a su estado fundamental sin pandear, ningn cilindro regreso al 100% a su estado inicial, siempre permaneci una pequea deformacin que daba pauta a la generacin de un nuevo modo de pandeo, si se volva a aplicar carga al espcimen. Es decir que a los especmenes despus de aplicarle varios ciclos de carga (aplicando carga y retirando a un valor de cero), la carga de pandeo tiende a tomar valores diferentes respecto de la primera carga de pandeo, pero el modo de pandeo se puede volver a registrar donde queda la deformacin permanente del primer modo de pandeo. Los primeros ciclos de carga pueden generar el mismo modo de pandeo, sin embargo ciclos de carga superiores pueden dar como resultado otros modos de pandeo, que no necesariamente corresponden al primero. Esto da paso a mencionar una hiptesis que pudiera establecerse de la siguiente manera, "El aplicar y retirar cclicamente la carga que provoca pandeo en un cascarn cargado axisimetricamente, parece provocar una resistencia poscrtica, ya que en los especimenes ensayados se puede apreciar que despus de N ciclos se produce un modo diferente de pandeo, que ya no es el inicial, probablemente debido a endurecimiento en el material de fabricacin".
El diseo de faldones para silos se ve caracterizado por el uso de factores de seguridad de 3 y usan la metodologa de Baker&Kovalesky, sin embargo factores de seguridad de 115
2.5 pueden dar buenos resultados, siempre y cuando se controlen las condiciones de frontera del silo.
El uso del eurocdigo parece estar limitado a factores de calidad (Q) mayores de 16 y menores o iguales a 40. Usar factores de calidad menores de 16 en las formulas del eurocdigo disminuyen drsticamente las cargas de inestabilidad esperadas, sin embargo los anlisis no lineales presentados en este trabajo muestran que con imperfecciones muy grandes (factores de calidad menores a 16) se producen ciertas plastificaciones locales en las cercanas de la imperfeccines, provocando con esto un posible endurecimiento y/o redistribucin en la rigidez del especimen que da cmo resultado una resistencia poscrtica, sin embargo esto sera algo que habra que estudiarse ms a detalle.
El uso de factores de calidad (Q) mayores de 40 provocan un comportamiento asinttico que se va aproximando a la carga critica clsica, sin embargo los anlisis no lineales de inestabilidad hechos en esta tesis muestran que se generan cargas de inestabilidad muy por debajo de esta metodologa, para factores Q mayores de 40.
Para el uso del eurocdigo, es de suma importancia el control en las condiciones de frontera, es decir que se debe controlar que un cilindro sometido a compresin axial est perfectamente apoyado sobre su base o cimentacin en todo su permetro, en caso de no lograrlo, factores de calidad Q de 5 e incluso 6, son justificables, tal y como se aprecia en las los especmenes 1, 3 y 5. Mientras que si se usa la tcnica de Rish y no se tiene un apoyo total sobre la cimentacin, factores de seguridad de 5 y 6 pueden ser aplicables con un elevado costo en los materiales empleados.
Es necesario que en Mxico se establezcan normas de aplicacin nacional que regulen el diseo de este tipo de elementos estructurales. En Mxico es muy comn usar factores de seguridad de 3 en el diseo de faldones de silos metlicos, mientras que en el cuerpo de almacenamiento factores de 2.5 y 2 son usados. Despus de analizar el comportamiento de faldones de silos, as como la influencia de las imperfecciones geomtricas en estos elementos estructurales se puede concluir que factores de 116
seguridad de 3 son justificados para el diseo de faldones metlicos. Pero es muy importante verificar que el faldn asiente perfectamente sobre la cimentacin o base sobre la que se apoye. Parece un poco confuso e incluso complicado el medir y controlar las imperfecciones geomtricas en un silo para su diseo, ya que stas dependen de los procesos de fabricacin que intervienen en la manufactura de silos. Sin embargo, se pueden recomendar tolerancias de fabricacin que tiendan a minimizar las imperfecciones geomtricas en los silos, esto depender mucho del costo que implique manejar tolerancias de fabricacin ms justas.
En la prctica no es fcil eliminar las imperfecciones geomtricas en los faldones de silos, pero si se puede controlar su magnitud; por lo que en funcin de esto se puede tener un mtodo de diseo controlado.
La contribucin ms importante de esta tesis es la aplicacin industrial que tiene, ya que se estn usando los factores de seguridad recomendados para el diseo de silos metlicos.
Es importante realizar una metodologa de medicin en silos de tal forma que se puedan clasificar los silos por su grado de calidad y que a su vez sea un parmetro de control de calidad.
117
Recomendaciones
Es necesario que en Mxico se establezcan normas de legislacin nacional que regulen el diseo de este tipo de elementos estructurales. Desarrollar metodologas que ayuden a medir y controlar las imperfecciones geomtricas en silos metlicos.
En la instalacin de un silo metlico, asegurarse que el permetro apoye perfectamente en la cimentacin, mediante la aplicacin de un aglutinante (grouting)
Desarrollar una metodologa especial para silos del tipo atornillado, ya que lo anteriormente expuesto en esta tesis es aplicable a silos soldados, aunque se puede aplicar a silos atornillados. En silos atornillados se tienen mayores imperfecciones geomtricas que las que se dan en un silo soldado.
118
Trabajo futuros
Se pueden visualizar algunos trabajos a futuro que pueden contribuir en el proceso de mejorar el diseo de silos metlicos, como por ejemplo: Caracterizacin de materiales granulares y su influencia en el diseo de silos. Realizar una metodologa prctica para medir y controlar las imperfecciones geomtricas en los procesos de manufactura de silos. Verificar cual es la influencia en la determinacin de la carga de inestabilidad en silos debida al paso de tuberas Como contribuye el material almacenado en el cuerpo del silo a mejorar su resistencia al pandeo? Realizar anlisis combinados de faldones con elementos rigidizantes verticales y su influencia en la obtencin de la carga de inestabilidad Influencia de costuras de traslape verticales en silos atornillados
119
REFERENCIAS
ARTICULOS
[Chubb. 1996] Chubb Gary A., McGuire Keith E. 1996. Designing your storage tank to resist calculated loads. Proc. Powder and bulk Engineering. Volume 10 number 11, November 1996 pag 37. ISSN 0897-6627.
[Chajes. 1969] Chajes Alexander. 1969. Stability behavior illustrated by simple models. Journal of the American Society of Civil Engineers. ST 6 pp 1153
[Eurocode 3. 2003a] European Standard Intermediate Draft Confidential. 2003. Eurocode 3: Design of steel structures. Part 1.6: General Rules: Strength and stability of shell structures. Draft prEN1993-1-6, June 2003. CEN/TC250/SC3/
[Eurocode 3. 2003b] European Standard. Second preliminary Draft Confidential. 2003. Eurocode 3: Design of steel structures. Part 4.1: Silos. Draft prEN1993-4-1, June 2003. CEN/TC250/SC3/
[Rotter. 1996] Rotter, J.M. 1996. Buckling and collapse in internally pressurized axially compressed silo cylinders with measured axisymetric imperfections: Imperfection, residual stresses and local collapse. Proc. International Workshop on Imperfections in Metal Silos: Measurement, Characterization and Strength Analysis, CA-Silo, Lyon, France, 19 April pp 119-140
[Rotter. 1998a] Rotter, J.M. 1998a. Shell structures: The new European standard and current research needs. Thin Walled Structures, Vol. 31, pp 3-23
- 120 -
[Rotter. 1998b] Rotter, J.M. .1998b. Develop of proposed European Design rules for buckling of axially compressed cylinders. Proc. Advances in Structural Engineering, Vol 1, No. 4, July, pp 273-286.
[Rotter. 1999] Rotter, J.M., Holst Mark F.G. y Calladine C.R. 1999. Imperfections and Buckling in cylindrical shells with consistent residual stresses. Proc. Journal of Constructional Steel Research 54 (2000). Elsevier Science Ltd PII:S0143-974X(99)00047-4.
www.elsevier.com/locate/jcsr Agosto 2003
[Rotter. 2003] Rotter, J.M. 2003. Buckling of cylindrical shell under axial compression. Chapter 2 the book "Buckling of Thin Metal Structures" edited by JG Teng and JM Rotter, to be published by Spon, London in about December 2003.
[Sanchez. 2003] Sanchez S.H. Minchaca M.E., Romero A.J. 2003. Pandeo de estructuras axisimtricas con imperfecciones geomtricas iniciales. Proc. XII Congreso de Ingeniera Estructural, Sociedad Mexicana de Ingeniera Estructural, A.C. vol 159 pp 837-848
[Urrutia. 2002] Nachtigall Ingolf, Gebbeken Norbert, Urrutia Galicia J.L.. 2002. On the analysis of vertical circular cylindrical tanks under earthquake excitation at its base. Elsevier. August 2002.
- 121 -
BIBLIOGRAFA
[AISC. 1989] AISC. 1989. Manual of Steel Construction, Allowable Stress Desigs, 9th ed. American Institute of Steel construction Inc
[API Standard 620. 1978] API Standard 620. 1978. Recommended Rules for Design and Construction of Large Welded Low - Pressure Tanks, 6th ed. American Petroleum Institute, Washington, D.C.
[AWWA D100-79. 1979] AWWA D100-79. 1979. Standard for Welded Steel Tank for Water Storage, American Water Works Association, Denver, Colorado.
[Ansys. Manual R8.1] Ansys tutorial para la versin 8.1.
[Baker. 1972] Baker E.H., Kovalesky L., Rish F.L. 1986. Structural Analysis of Shells. Mc Graw Hill. ISBN 0-89874-118-1
[Bresler. 1970] Bresler, Lin y Scalzi. 1970. Diseo de Estructuras de Acero. Editorial Limusa.
[Calladine. 1988] Calladine C.R. 1988 Theory of structures, Ed. CambrIdge University Press, ISBN 0-52136945-2. 763 P.
[Fung. 1972] Fung Y.C., Sechler E.E. 1972. Thin Shell Structures, Theory, Experiment and Design. Prentice Hall, Inc. ISBN 0-13-918193-8
- 122 -
[Gaylord. 1984] Gaylord E.H. and Gaylord, C.N. 1984. Design of steel Bins for Storage of Bulks Solids, Prentice Hall, Inc Englewood Cliffs, N.J.
[Gaylord. 1990] Gaylord E.H. and Gaylord, C.N. 1990. Structural Engineering Handbook. Third edition, Mac Graw Hill. ISBN 0-07-023188-5
[Goldenveizer. 1963] Goldenveizer A.L. 1963. Cascarones Elsticos Delgados. Compaa Editorial Continental, S.A.
[Notas de Mecnica estructural II. 2002] Sanchez S.H. 2002. Notas de Mecnica estructural II. Curso de Posgrado. Escuela Superior de Ingenieria y Arquitectura del Instituto Politcnico Nacional. Mxico, DF. 2002. Semestre escolar 2002/2
[Ravenet. 1977] Ravenet Catalan Juan. 1977. Silos, Teora, Investigacin, Construccin. Tcnicos Asociados, S.A. ISBN 84-7146-088-2 Editores
INTERNET
[http 1. 2003] http://www.columbiantank.com/capacity/dry_capacity_charts.pdf 2003 6 de noviembre del
- 123 -
Anexo 1
Medicin de imperfecciones.
Las imperfecciones se midieron dividiendo los cilindros a cada 10 grados, respetando la longitud de medicin establecida por el Eurocdigo.
Por cuestiones de espacio, no se muestran las imperfecciones de cada meridiano, sin embargo se pueden consultar en las hojas de excel que se agregan al trabajo en el CD del anexo 3 de videos. Las imperfecciones aqu mostradas estn corregidas por la matriz de transformacin que se muestra en el anexo 2 y se presentan en mm.
Especimen 1
ANGULO ALTURA 1' 2' 3' 4' 5' 6' 7' 8' 9' 10' 11' 12' 13' MXIMO MNIMO 0 25 50 75 100 125 150 175 200 225 250 275 300 0.580 -0.853 0 yi 10 yi 20 yi 30 yi 40 yi 50 yi 60 yi 70 yi 80 yi 90 yi 100 yi 110 yi 120 yi 130 yi 140 yi 150 yi 160 yi 170 yi 180 yi
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 -0.038 0.008 0.002 0.004 0.063 0.059 -0.047 -0.148 -0.061 0.023 0.047 -0.178 -0.066 -0.239 0.078 -0.004 -0.032 0.025 0.025 -0.051 -0.008 -0.021 -0.017 0.178 0.093 -0.144 -0.246 -0.047 0.097 0.068 -0.229 -0.106 -0.174 0.131 -0.034 -0.089 0.076 0.051 -0.063 -0.025 -0.044 0.038 0.241 0.127 -0.241 -0.292 -0.083 0.171 0.038 -0.152 -0.171 -0.083 0.184 -0.063 -0.146 0.051 0.381 -0.076 -0.042 -0.119 0.068 0.279 0.135 -0.288 -0.339 -0.093 0.246 0.034 -0.254 -0.186 -0.068 0.186 -0.068 -0.178 -0.025 0.076 -0.089 -0.008 -0.142 0.580 0.292 0.144 -0.284 -0.360 -0.104 0.345 -0.021 -0.254 -0.176 -0.002 0.188 -0.097 -0.210 -0.076 0.051 -0.102 0.000 -0.190 0.051 0.330 0.076 -0.279 -0.381 -0.089 0.419 -0.076 -0.254 -0.190 -0.216 0.190 -0.127 -0.190 -0.025 0.102 -0.089 0.008 -0.188 0.030 0.368 0.034 -0.275 -0.377 -0.074 0.392 -0.055 -0.229 -0.180 0.078 0.167 -0.106 -0.197 0.000 0.127 -0.051 0.017 -0.161 -0.017 0.330 0.042 -0.296 -0.296 -0.034 0.262 -0.008 -0.203 -0.169 0.093 0.119 -0.135 -0.178 0.025 0.102 -0.038 0.000 -0.133 -0.038 0.317 0.051 -0.292 -0.241 0.006 0.184 0.013 -0.152 -0.133 0.083 0.044 -0.140 -0.159 0.000 0.051 -0.025 -0.042 -0.080 -0.059 0.254 0.034 -0.212 -0.161 0.021 0.106 0.008 -0.102 -0.123 0.072 0.047 -0.119 -0.140 -0.051 -0.025 0.241 -0.034 -0.053 -0.055 0.114 0.017 -0.131 -0.106 -0.015 0.028 0.004 -0.076 -0.163 0.036 -0.028 -0.072 -0.095 -0.076 -0.051 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.241 0.017 0.002 0.580 0.368 0.144 0.000 0.000 0.021 0.419 0.068 0.000 0.000 0.093 0.190 0.000 0.000 0.076 0.381 -0.102 -0.042 -0.190 -0.059 0.000 0.000 -0.296 -0.381 -0.104 0.000 -0.076 -0.254 -0.190 -0.239 -0.028 -0.140 -0.210 -0.076 -0.051
Especimen 2 ANGULO
ALTURA 1' 2' 3' 4' 5' 6' 7' 8' 9' 10' 11' 12' 13' MXIMO MNIMO 0 25 50 75 100 125 150 175 200 225 250 275 300 2.498 -0.832
0 yi 0.000
10 yi 0.000
20 yi
30 yi
40 yi
50 yi
60 yi 0.000
70 yi
80 yi
90 yi
100 yi
110 yi
120 yi
130 yi
140 yi
150 yi
160 yi
170 yi
180 yi
0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.036 -0.017 0.021 -0.034 0.032 -0.051 0.028 -0.034 0.063 -0.051 0.023 -0.068 0.008 -0.059
0.070 -0.091 0.104 0.038 -0.353 0.004 0.159 -0.273 0.337 0.267 -0.832 0.013 0.203 -0.313 0.389 0.330 -0.804 0.017 0.197 -0.303 0.417 0.343 -0.701 0.047 0.140 -0.292 0.419 0.305 -0.622 0.076 0.133 -0.282 0.370 0.267 -0.519 0.080 0.127 -0.271 0.296 0.229 -0.364 0.085 0.095 -0.260 0.248 0.216 -0.286 0.114 0.114 -0.199 0.148 0.127 -0.182 0.144 0.057 -0.112 0.074 0.038 -0.078 0.123 0.000 0.203 0.000 0.000 0.000 0.000 0.000 0.000 0.419 0.343 0.000 0.144
0.063 -0.002 0.066 0.163 0.135 -0.002 0.180 0.091 0.129 0.000 0.066 0.165 -0.083 0.070 0.108 0.152 -0.083 0.083 0.019 0.159 -0.025 0.095 0.216 -0.112 0.049 0.231 0.169 -0.112 -0.040 -0.028 0.112 0.000 0.099 0.203 -0.114 0.013 0.114 0.152 0.114 -0.038 -0.038 0.063 0.000 0.089 0.216 -0.142 0.028 0.074 0.135 -0.116 -0.036 -0.049 0.040 -0.025 0.078 0.203 0.102 0.038 0.000 0.216 0.008 0.068 0.059 0.068 -0.119 -0.034 -0.059 0.296 -0.025 0.068
0.089 -0.182 0.258 0.152 -0.656 -0.017 0.127 -0.030 0.080 0.123 0.169 -0.055 0.131 0.055 0.157 -0.025 0.106
0.203 -0.110 0.034 0.119 0.161 -0.110 0.008 -0.017 0.135 -0.025 0.110 -0.008 -0.017
0.165 -0.070 0.057 0.019 0.102 -0.070 -0.032 -0.044 -0.006 -0.025 0.057 -0.006 -0.025 0.030 0.097 0.004 0.008 -0.047 -0.004 -0.055 -0.004 -0.025 0.021 -0.021 2.498 0.053 0.036 -0.011 0.017 -0.023 -0.028 -0.040 -0.028 -0.025 0.011 -0.036 -0.008 0.000 0.000 0.000 0.000 0.000 0.053 0.097 0.231 0.169 0.114 0.000 0.000 0.000 0.000 0.000 0.180 0.091 0.296 0.000 0.110 0.000 0.000 0.063 2.498
0.000 -0.313 0.000 0.000 -0.832 -0.017 0.000 -0.142 0.000 -0.011 0.000 -0.119 -0.040 -0.059 -0.028 -0.025 0.000 -0.036 -0.068
124
Especimen 3 ANGULO
0 yi 0.000
10 yi 0.000
20 yi
30 yi
40 yi
50 yi
60 yi 0.000 0.070 0.159
70 yi
80 yi
90 yi
100 yi
110 yi
120 yi
130 yi
140 yi
150 yi
160 yi
170 yi
180 yi
0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.059 0.059 -0.114 -0.085 0.127 0.076 0.013 -0.171 -0.102 0.140 0.093 0.017 -0.203 -0.093 0.152 0.110 0.021 -0.235 0.195 0.165 0.127 0.025 -0.241 0.000 0.127 0.119 0.004 -0.222 0.034 0.114 0.110 -0.017 -0.203 0.042 0.178 0.127 -0.013 -0.159 0.076 0.140 0.093 -0.008 -0.089 0.059 0.076 0.085 -0.004 -0.019 0.042 0.038 0.000 0.000 0.000 0.000 0.000 0.127 0.059 0.000 0.195 0.178
0.000 0.000 0.008 -0.089 0.025 -0.146 0.017 -0.152 0.008 -0.133 0.025 -0.089 0.042 -0.044 0.161 0.000 0.076 0.019 0.017 0.013 0.034 0.032 0.000 0.000 0.161 0.032
-0.068 -0.076 0.047 0.226 0.195 0.036 -0.059 -0.203 0.093 0.402 0.415 0.148 -0.076 -0.279 0.140 0.603 0.609 0.286 -0.093 -0.330 0.186 0.703 0.753 0.398 -0.085 -0.356 0.233 0.777 0.872 0.485 -0.076 -0.381 0.279 0.825 0.940 0.521 -0.068 -0.356 0.275 0.747 1.007 0.557 -0.059 -0.330 0.296 0.847 0.897 0.567 -0.051 -0.305 0.292 0.590 0.736 0.502 -0.068 -0.254 0.288 0.461 0.576 0.411 -0.034 -0.127 0.182 0.307 0.364 0.269 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.296 0.847 1.007 0.567
0.013 0.823 0.114 0.059 0.047 -0.008 0.030 -0.032 -0.042 0.063 -0.008 -0.083 0.038 0.565 0.140 0.076 0.089
0.114 -0.025 0.681 0.152 0.042 0.068 0.229 -0.025 0.474 0.102 0.110 0.110 0.324 -0.013 0.358 0.063 0.144 0.131 0.394 0.432 0.400 0.292 0.210 0.000 0.438 0.000 0.267 0.051 0.102 0.178 0.000 0.161 0.000 0.119 0.144 0.013 0.121 0.013 0.076 0.140 0.025 0.131 0.025 0.085 0.135 0.013 0.091 0.038 0.068 0.106 0.000 0.000 0.000 0.000 0.000 0.038 0.823 0.152 0.144 0.178 0.438 -0.013 0.201 0.038 0.135 0.148
-0.093 -0.381 0.000 0.000 0.000 0.000
0.000 -0.025 0.000 0.000 0.000 0.000 -0.008 -0.017 -0.241 -0.102 0.000 -0.008 -0.152
Especimen 4 ANGULO
0 yi 0.000 0.237 0.220 0.178 0.161 0.119 0.076 0.110 0.042 0.025 0.008 0.017 0.000 0.237 0.000
10 yi 0.000 0.203 0.254 0.254 0.279 0.254 0.229 0.229 0.178 0.127 0.127 0.051 0.000 0.279 0.000
20 yi
30 yi
40 yi
50 yi
60 yi 0.000 0.042
70 yi
80 yi
90 yi
100 yi
110 yi
120 yi
130 yi
140 yi
150 yi
160 yi
170 yi
180 yi
0.000 0.000 0.000 0.000 0.102 0.028 0.030 0.049 0.102 -0.021 -0.017 0.072 0.127 -0.070 -0.038 0.095 0.152 -0.068 -0.085 0.093 0.203 -0.091 -0.131 0.142 0.203 -0.114 -0.102 0.114 0.178 -0.112 -0.148 0.112 0.152 -0.110 -0.144 0.110 0.102 -0.108 -0.140 0.108 0.102 -0.055 -0.085 0.080 0.051 -0.053 -0.055 0.028 0.000 0.000 0.000 0.000 0.203 0.028 0.030 0.142 0.000 -0.114 -0.148 0.000
0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000
0.000 0.000
0.019 -0.260 0.008 0.047 0.019 -0.008 0.040 0.091 0.013 -0.138 0.123 -0.506 0.000 -0.032 0.019 -0.038 -0.159 0.165 -0.044 0.017 -0.017 0.008 -0.051 -0.144 0.186 -0.017 0.008 -0.028 -0.002 -0.063 -0.155 0.157 -0.040 0.000 0.013 -0.013 -0.102 -0.140 0.152 -0.063 0.042 0.002 -0.049 -0.038 -0.125 0.148 -0.061 0.008 0.017 -0.059 -0.025 -0.135 0.119 -0.059 0.025 -0.019 -0.095 0.013 -0.146 0.089 -0.032 0.093 -0.004 -0.080 0.000 -0.131 0.110 -0.004 0.000 0.000 0.000 0.000 0.000 0.135 0.040 0.091 0.013 0.000 0.000 0.000 0.186 0.049
0.034 -0.013 -0.038 0.017 0.017 0.063 -0.017 0.004 0.055 0.000 -0.148 0.144 -0.021 0.051 -0.044 -0.044 -0.025 0.013 0.032 0.042 -0.051 -0.076 -0.042 0.034 0.051 0.034 -0.057 -0.057 -0.034 0.004 0.070 0.076 -0.063 -0.038 -0.025 0.000 0.063 0.042 -0.070 -0.044 -0.017 -0.004 0.083 0.034 -0.076 -0.051 -0.008 0.017 0.076 0.051 -0.057 -0.032 0.000 -0.013 0.070 0.042 -0.038 -0.013 0.008 0.034 0.013 0.000 0.076 0.000 0.000 0.000 0.000 0.000 0.019 0.032 0.017 0.047 0.083
0.034 -0.044 0.032 -0.008 0.030 -0.121 0.135 -0.015 -0.040 0.013 -0.040 0.080 0.049
0.000 -0.076 -0.260 -0.042 -0.013 -0.121 -0.017 -0.032 -0.095 -0.102 -0.159 0.000 -0.506
Especimen 5 ANGULO
1' 2' 3' 4' 5' 6' 7' 8' 9' 10' 11' 12' 13' MXIMO MNIMO 0 25 50 75 100 125 150 175 200 225 250 275 300 0.387 -0.398
0 yi 0.000
10 yi 0.000
20 yi
30 yi
40 yi
50 yi
60 yi 0.000
70 yi
80 yi
90 yi
100 yi
110 yi
120 yi
130 yi
140 yi
150 yi
160 yi
170 yi
180 yi
0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.159 0.119 0.036 -0.013 0.047 0.171 0.305 0.032 -0.114 0.089
0.000 0.000 0.110 -0.013 0.152 -0.089 0.161 -0.076 0.195 -0.114 0.178 -0.152 0.135 -0.216 0.093 -0.203 0.102 -0.165 0.085 -0.102 0.068 0.013 0.000 0.000 0.195 0.013
-0.002 -0.112 -0.072 -0.089 0.066 -0.142 0.239 -0.133 0.239 -0.260 0.207 0.165 0.021 -0.097 -0.042 0.000 0.106 -0.207 0.275 -0.140 0.377 -0.216 0.085 0.127 -0.032 -0.083 -0.063 0.063 0.095 -0.197 0.184 -0.197 0.387 -0.171 -0.013 0.038
0.190 0.237 0.021 -0.076 0.068 -0.364 -0.102
-0.034 -0.042 -0.059 0.051 0.135 -0.186 0.017 -0.254 0.373 -0.102 -0.059 -0.025 0.203 0.271 -0.034 -0.127 0.110 -0.036 -0.028 -0.030 -0.013 0.125 -0.201 -0.125 -0.286 0.282 -0.057 -0.080 -0.089 0.184 0.237 -0.074 -0.089 0.157 -0.063 -0.013 0.000 -0.051 0.165 -0.190 -0.267 -0.317 0.241 -0.038 -0.102 -0.127 0.190 0.127 -0.063 -0.051 0.152 -0.066 0.002 -0.093 0.017 -0.146 0.006 -0.174 0.047 0.000 0.021 0.000 0.047 0.004 -0.063 0.180 -0.332 -0.358 -0.324 0.176 -0.019 -0.072 -0.165 0.146 0.068 -0.078 -0.013 0.174 0.059 -0.102 0.195 -0.195 -0.398 -0.356 0.161 0.025 -0.042 -0.178 0.102 0.034 -0.093 0.000 0.220 0.063 -0.114 0.210 -0.159 -0.387 -0.362 0.146 -0.006 -0.038 -0.165 0.032 0.025 -0.108 0.013 0.165 0.042 -0.102 0.174 -0.148 -0.377 -0.368 0.131 0.013 -0.008 -0.127 0.013 -0.034 0.157 0.025 0.135 0.000 0.000 0.000 0.000 0.063 0.063 0.210 0.000 0.000 0.275 0.000 0.000 0.000 0.000 0.000 0.000 0.387 0.032 0.207 0.165 0.000 0.000 0.000 0.000 0.000 0.203 0.305 0.157 0.025 0.220
-0.099 0.011 -0.004 -0.089 0.138 -0.061 -0.239 -0.222 0.091 0.032 -0.004 -0.114 -0.006 -0.042 -0.061 0.013 0.080
-0.174 -0.112 -0.072 -0.114 0.000 -0.332 -0.398 -0.368 0.000 -0.260 -0.102 -0.178 -0.006 -0.042 -0.108 -0.127 0.000 -0.364 -0.216
125
Anexo 2
Matriz de transformacin
La matriz de transformacin empleada para mostrar las imperfecciones en los cilindros de manera absoluta es la siguiente. 0 0 U (X ,Y ) = X Y 0 0 cos( ) sin ( ) 0 0 sin ( ) cos( ) T 1( ) = 0 0 cos( ) sin ( ) 0 0 sin ( ) cos( )
Donde representa el ngulo que forman el sistema de referencia local contra el sistema de referencia global (ver figura A2-1), U(X,Y) representan las coordenadas de la imperfeccin geomtrica en la interseccin de un meridiano contra un paralelo, respecto al sistema de referencia local.
Por lo que la nueva ubicacin respecto al sistema de referencia global queda como sigue:
[XY ] = [T1( )]1 [U ( X , Y )]

0 0 [XY ] = X cos( ) + Y sin ( ) X sin ( ) + Y cos( )
126
Figura A2-1 Transformacin del sistema de referencia local a un sistema de referencia global
127
Anexo 3
Videos y archivos electrnicos
En el CD anexo se presentan los videos de las pruebas realizadas, as como las herramientas usadas para el desarrollo de esta tesis.
128

Tesis Completa

Cargado por

Copyright:

Formatos disponibles

Tesis Completa

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Tesis Completa

Cargado por

Copyright:

Formatos disponibles

INSTITUTO POLITCNICO NACIONAL

ESCUELA SUPERIOR DE INGENIERA Y ARQUITECTURA SECCIN DE ESTUDIOS DE POSGRADO E INVESTIGACIN

MAESTRO EN CIENCIAS CON ESPECIALIDAD EN ESTRUCTURAS P R E S E N T A :

A mi esposa Luz Edith Por todo el apoyo brindado y su gran comprensin

A mi hijo Jess Emmanuel Por entender mis ausencias Y siempre sonrer

A mis padres Por sus motivaciones y enseanzas

A mis maestros y amigos En especial al Dr. Hctor Snchez

Ciudad de Mxico Diciembre del 2005

1 Generalidades y planteamiento de la investigacin

Figura I-1 Clasificacin de silos. -2-

Figura I-2 Apoyos tpicos de silos metlicos.

Figura I-3 Partes de un silo

b.- distancia entre atiezadores

1.1.1 Cargas actuantes sobre un silo

1.1.2 Esfuerzos permisibles

Para acero, se pueden sustituir el mdulo de elasticidad y el coeficiente de Poisson

1.1.2.3 Factor de seguridad

1.2 Estado de la tcnica

1.3 Objetivo y justificacin

1.5 Planteamiento propuesto

2 Marco terico 2.1 Compresin axial en cilindros de pared delgada

2.1.1 Pandeo por bifurcacin con cada brusca de rigidez

Ncr N'cr N''cr 0 B'

Figura II-1. Pandeo por bifurcacin con cada brusca de rigidez

2.1.2 Pandeo terico de estructuras sometidas a compresin axial.

Carga axial Acortamiento axial c) Modo diamantado d) Curva de postpandeo

Figura II-2. Modos de pandeo

Figura II-3. Modelo de panel curvo

En el cual , es el acortamiento axial de la barra transversal y est dado por

La ecuacin II-5 es ms sencilla para interpretar si se introduce un parmetro de deformacin no dimensional:

Figura II-4. Comportamiento de pandeo.

adicionado a la energa potencial y la formacin de la carga se relaciona como sigue:

La ecuacin II-10, en la cual

e = 0.1 , se muestra en la figura II-4 como la curva 0-6-7. h

2.1.2.1 Tipos de pandeo terico de estructuras cilndricas

Et 2 nL Parmetro de Batdorf ; = ; D= Rigidez a flexin R 12 1 2

m = configuracin axial; n = configuracin circunferencial

2.1.3 Efecto de las imperfecciones en la estabilidad de cascarones cilndricos

Esfuerzo de pandeo normalizado k/

Esfuerzo ideal: Esfuerzo crtico clsico

Relacin radio espesor r/t

Carga por punto lmite

Pandeo Postcrtico de la estructura con imperfecciones

Desplazamiento del pandeo modal w'b

Figura II-7. Punto de bifurcacin y punto limite mediante un anlisis asinttico.

Carga por punto lmite

Carga por punto lmite del cascarn con imperfecciones

Figura II-8. Punto de bifurcacin y punto limite mediante un anlisis no-lineal.

2.1.5 Influencia de los esfuerzos residuales en la estabilidad de cascarones

Fotografa II-5. Silos atornillados (cortesa Silos y Camiones S.A. de C.V.)

2.2 Criterios tpicos de diseo

2.2.1 Compresin axial en cilindros de pared delgada, [Baker. 1972]

El efecto de la plasticidad puede considerarse por el factor de correccin por plasticidad

K c 0 se expresa en funcin de las modos circunferenciales n y verticales m as como de

Y para cilindros moderadamente largos z > 2 K c 0 2 3 por:

= 0.6 Siempre que z > 2 K c 0

2.2.2 Eurocdigo 3 [Eurocode 3. 2003b]

Silos con capacidades entre 10 y 100 toneladas

2.2.2.1 Medicin de las imperfecciones geomtricas despus de la fabricacin