Masonry macro-block analysis

The structural analysis involves the definition of the model and selection of the analysis type. The model should represent the stiffness, the mass and the loads of the structure. The structures can be represented using simplified models, such as the lumped mass models, and advanced models resorting the Finite Element Method (FEM) and Discrete Element Method (DEM). Depending on the characteristics of the structure, different types of analysis can be used such as limit analysis, linear and non-linear static analysis and linear and non-linear dynamic analysis. Unreinforced masonry structures present low tensile strength and the linear analyses seem to not be adequate for assessing their structural behaviour. On the other hand, the static and dynamic non-linear analyses are complex, since they involve large time computational requirements and advanced knowledge of the practitioner. The non-linear analysis requires advanced knowledge on the material properties, analysis tools and interpretation of results. The limit analysis with macro-blocks can be assumed as a more practical method in the estimation of maximum load capacity of structure. Furthermore, the limit analysis require a reduced number of parameters, which is an advantage for the assessment of ancient and historical masonry structures, due to the difficult in obtaining reliable data.

Numerical and experimental analysis of full scale arches reinforced with GFRP materials

In this contribution, original limit analysis numerical results are presented dealing with some reinforced masonry arches tested at the University of Minho-UMinho, PT. Twelve in-scale circular masonry arches were considered, reinforced in various ways at the intrados or at the extrados. GFRP reinforcements were applied either on undamaged or on previously damaged elements, in order to assess the role of external reinforcements even in repairing interventions. The experimental results were critically discussed at the light of limit analysis predictions, based on a 3D FE heterogeneous upper bound approach. Satisfactory agreement was found between experimental evidences and the numerical results, in terms of failure mechanisms and peak load.

The Islamic Crossed-Arch Domes in Cordoba: Geometry and Structural Analysis of the “Capilla de Villaviciosa”

Crossed-arch domes are a singular type of ribbed vaults. Their characteristic feature is that the ribs that form the vault are intertwined, forming polygons or stars, leaving an empty space in the centre. The earliest known vaults of this type are found in the Great Mosque of Córdoba, built ca. 960 a.C. The type spread through Spain, and the north of Africa in the 10th to the 16th Centuries, and was used by Guarini and Vittone in the 17th and 18th Centuries in Italy. However, it was used only in a few buildings. Though the literature about the structural behaviour of ribbed Gothic vaults is extensive, so far no structural analysis of crossed arch domes has been made. The purpose of this work is, first to show the way to attack such an analysis within the frame of Modern Limit Analysis of Masonry Structures (Heyman 1995), and then to apply the approach to study the stability of the dome of the Capilla de Villaviciosa. The work may give some clues to art and architectural historians to understand better the origin and development of Islamic dome architecture.

The Value Of The 2005 International Society Of Urological Pathology (isup) Modified Gleason Grading System As A Predictor Of Biochemical Recurrence After Radical Prostatectomy.

To compare time and risk to biochemical recurrence (BR) after radical prostatectomy of two chronologically different groups of patients using the standard and the modified Gleason system (MGS). Cohort 1 comprised biopsies of 197 patients graded according to the standard Gleason system (SGS) in the period 1997/2004, and cohort 2, 176 biopsies graded according to the modified system in the period 2005/2011. Time to BR was analyzed with the Kaplan-Meier product-limit analysis and prediction of shorter time to recurrence using univariate and multivariate Cox proportional hazards model. Patients in cohort 2 reflected time-related changes: striking increase in clinical stage T1c, systematic use of extended biopsies, and lower percentage of total length of cancer in millimeter in all cores. The MGS used in cohort 2 showed fewer biopsies with Gleason score ≤ 6 and more biopsies of the intermediate Gleason score 7. Time to BR using the Kaplan-Meier curves showed statistical significance using the MGS in cohort 2, but not the SGS in cohort 1. Only the MGS predicted shorter time to BR on univariate analysis and on multivariate analysis was an independent predictor. The results favor that the 2005 International Society of Urological Pathology modified system is a refinement of the Gleason grading and valuable for contemporary clinical practice.

Estudo da cabeceira da Sé do Porto

Dissertação de mestrado integrado em Engenharia Civil

The Lithium, Boron and Beryllium content of serpentinized peridotites from ODP Leg 209 (Sites 1272A and 1274A): Implications for lithium and boron budgets of oceanic lithosphere

Despite the key importance of altered oceanic mantle as a repository and carrier of light elements (B, Li, and Be) to depth, its inventory of these elements has hardly been explored and quantified. In order to constrain the systematics and budget of these elements we have studied samples of highly serpentinized (>50%) spinel harzburgite drilled at the Mid-Atlantic Ridge (Fifteen-Twenty Fracture zone, ODP Leg 209, Sites 1272A and 1274A). In-situ analysis by secondary ion mass spectrometry reveals that the B, Li and Be contents of mantle minerals (olivine, orthopyroxene, and clinopyroxene) remain unchanged during serpentinization. B and Li abundances largely correspond to those of unaltered mantle minerals whereas Be is close to the detection limit. The Li contents of clinopyroxene are slightly higher (0.44-2.8 mu g g(-1)) compared to unaltered mantle clinopyroxene, and olivine and clinopyroxene show an inverse Li partitioning compared to literature data. These findings along with textural observations and major element composition obtained from microprobe analysis suggest reaction of the peridotites with a mafic silicate melt before serpentinization. Serpentine minerals are enriched in B (most values between 10 and 100 mu g g(-1)), depleted in Li (most values below I mu g g(-1)) compared to the primary phases, with considerable variation within and between samples. Be is at the detection limit. Analysis of whole rock samples by prompt gamma activation shows that serpentinization tends to increase B (10.4-65.0 mu g g(-1)), H2O and Cl contents and to lower Li contents (0.07-3.37 mu g g(-1)) of peridotites, implying that-contrary to alteration of oceanic crust-B is fractionated from Li and that the B and Li inventory should depend essentially on rock-water ratios. Based on our results and on literature data, we calculate the inventory of B and Li contained in the oceanic lithosphere, and its partitioning between crust and mantle as a function of plate characteristics. We model four cases, an ODP Leg 209-type lithosphere with almost no igneous crust, and a Semail-type lithosphere with a thick igneous crust, both at I and 75 Ma, respectively. The results show that the Li contents of the oceanic lithosphere are highly variable (17-307 kg in a column of I m x I m x thickness of the lithosphere (kg/col)). They are controlled by the primary mantle phases and by altered crust, whereas the B contents (25-904 kg/col) depend entirely on serpentinization. In all cases, large quantities of B reside in the uppermost part of the plate and could hence be easily liberated during slab dehydration. The most prominent input of Li into subduction zones is to be expected from Semail-type lithosphere because most of the Li is stored at shallow levels in the plate. Subducting an ODP Leg 209-type lithosphere would mean only very little Li contribution from the slab. Serpentinized mantle thus plays an important role in B recycling in subduction zones, but it is of lesser importance for Li. (C) 2008 Elsevier Ltd. All rights reserved.

Modeling, Estimation, and Control of Robot-Soil Interactions

This thesis presents the development of hardware, theory, and experimental methods to enable a robotic manipulator arm to interact with soils and estimate soil properties from interaction forces. Unlike the majority of robotic systems interacting with soil, our objective is parameter estimation, not excavation. To this end, we design our manipulator with a flat plate for easy modeling of interactions. By using a flat plate, we take advantage of the wealth of research on the similar problem of earth pressure on retaining walls. There are a number of existing earth pressure models. These models typically provide estimates of force which are in uncertain relation to the true force. A recent technique, known as numerical limit analysis, provides upper and lower bounds on the true force. Predictions from the numerical limit analysis technique are shown to be in good agreement with other accepted models. Experimental methods for plate insertion, soil-tool interface friction estimation, and control of applied forces on the soil are presented. In addition, a novel graphical technique for inverting the soil models is developed, which is an improvement over standard nonlinear optimization. This graphical technique utilizes the uncertainties associated with each set of force measurements to obtain all possible parameters which could have produced the measured forces. The system is tested on three cohesionless soils, two in a loose state and one in a loose and dense state. The results are compared with friction angles obtained from direct shear tests. The results highlight a number of key points. Common assumptions are made in soil modeling. Most notably, the Mohr-Coulomb failure law and perfectly plastic behavior. In the direct shear tests, a marked dependence of friction angle on the normal stress at low stresses is found. This has ramifications for any study of friction done at low stresses. In addition, gradual failures are often observed for vertical tools and tools inclined away from the direction of motion. After accounting for the change in friction angle at low stresses, the results show good agreement with the direct shear values.

Approccio agli spostamenti per la valutazione della vulnerabilità sismica di telai in conglomerato cementizio armato.

Il presente lavoro ha lo scopo di presentare gli studi e i risultati ottenuti durante l’attività di ricerca svolta sul Displacement-based Assessment (DBA) dei telai in cemento armato. Dopo alcune considerazioni iniziali sul tema della vulnerabilità sismica e sui metodi di analisi e verifica, si procede alla descrizione teorica del metodo. Sono stati analizzati tre casi studio di telai piani, progettati per soli carichi verticali e secondo normative non più in vigore che non prevedevano l’applicazione della gerarchia delle resistenze. I telai considerati, destinati ad abitazione civile, hanno diversa altezza e numero di piani, e diverso numero di campate. Si è proceduto all’applicazione del metodo, alla valutazione della vulnerabilità sismica in base alla domanda in termini di spostamento costituita da uno spettro elastico previsto dall’EC8 e alla validazione dei risultati ottenuti mediante analisi non lineari statiche e dinamiche e mediante l’applicazione dei teoremi dell’Analisi limite dei telai, proposta come procedura alternativa per la determinazione del meccanismo anelastico e della capacità in termini di taglio alla base. In ultimo si è applicata la procedura DBA per la valutazione della vulnerabilità sismica di un edificio scolastico, realizzato tra il 1969 e il 1975 in un sito caratterizzato da una accelerazione di picco orizzontale pari a 0,24g e una probabilità di superamento del 10% in 75 anni.

Geochemistry of serpentinized peridotites from ODP Holes 207-1272A and 207-1274A

Despite the key importance of altered oceanic mantle as a repository and carrier of light elements (B, Li, and Be) to depth, its inventory of these elements has hardly been explored and quantified. In order to constrain the systematics and budget of these elements we have studied samples of highly serpentinized (>50%) spinel harzburgite drilled at the Mid-Atlantic Ridge (Fifteen-Twenty Fracture zone, ODP Leg 209, Sites 1272A and 1274A). In-situ analysis by secondary ion mass spectrometry reveals that the B, Li and Be contents of mantle minerals (olivine, orthopyroxene, and clinopyroxene) remain unchanged during serpentinization. B and Li abundances largely correspond to those of unaltered mantle minerals whereas Be is close to the detection limit. The Li contents of clinopyroxene are slightly higher (0.44-2.8 µg/g) compared to unaltered mantle clinopyroxene, and olivine and clinopyroxene show an inverse Li partitioning compared to literature data. These findings along with textural observations and major element composition obtained from microprobe analysis suggest reaction of the peridotites with a mafic silicate melt before serpentinization. Serpentine minerals are enriched in B (most values between 10 and 100 µg/g), depleted in Li (most values below 1 µg/g) compared to the primary phases, with considerable variation within and between samples. Be is at the detection limit. Analysis of whole rock samples by prompt gamma activation shows that serpentinization tends to increase B (10.4-65.0 µg/g), H2O and Cl contents and to lower Li contents (0.07-3.37 µg/g) of peridotites, implying that-contrary to alteration of oceanic crust-B is fractionated from Li and that the B and Li inventory should depend essentially on rock-water ratios. Based on our results and on literature data, we calculate the inventory of B and Li contained in the oceanic lithosphere, and its partitioning between crust and mantle as a function of plate characteristics. We model four cases, an ODP Leg 209-type lithosphere with almost no igneous crust, and a Semail-type lithosphere with a thick igneous crust, both at 1 and 75 Ma, respectively. The results show that the Li contents of the oceanic lithosphere are highly variable (17-307 kg in a column of 1 m * 1 m * thickness of the lithosphere (kg/col)). They are controlled by the primary mantle phases and by altered crust, whereas the B contents (25-904 kg/col) depend entirely on serpentinization. In all cases, large quantities of B reside in the uppermost part of the plate and could hence be easily liberated during slab dehydration. The most prominent input of Li into subduction zones is to be expected from Semail-type lithosphere because most of the Li is stored at shallow levels in the plate. Subducting an ODP Leg 209-type lithosphere would mean only very little Li contribution from the slab. Serpentinized mantle thus plays an important role in B recycling in subduction zones, but it is of lesser importance for Li.

The Ciborium or Lantern Tower of Valencia Cathedral: Geometry, Construction and Stability

Early 18th century treatise writer Tomas Vicente Tosca1 includes in his Tratado de la montea y cortes de Canteria [On Masonry Design and Stone Cutting], what is an important documentary source about the lantern of Valencia Cathedral. Tosca writes about this lantern as an example of vaulting over cross arches without the need of buttresses. A geometrical description is followed by an explanation of the structural behavior which manifests his deep understanding of the mechanics of masonry structures. He tries to demonstrate the absence of buttresses supporting his thesis on the appropriate distribution of loads which will reduce the "empujos" [horizontal thrusts] to the point of not requiring more than the thickness of the walls to stand (Tosca [1727] 1992, 227-230). The present article2 assesses T osca' s appreciation studying how loads and the thrusts they generate are transmitted through the different masonry elements that constitute this ciborium. In order to do so, we first present a geometrical analysis and make considerations regarding its materials and construction methods to, subsequently, analyze its stability adopting an equilibrium approach within the theoretical framework of the lower bound limit analysis.

The Ciborium or Lantern Tower of Valencia Cathedral: Geometry, Construction and Stability

Early 18th century treatise writer Tomas Vicente Tosca1 includes in his Tratado de la montea y cortes de Canteria [On Masonry Design and Stone Cutting], what is an important documentary source about the lantern of Valencia Cathedral. Tosca writes about this lantern as an example of vaulting over cross arches without the need of buttresses. A geometrical description is followed by an explanation of the structural behavior which manifests his deep understanding of the mechanics of masonry structures. He tries to demonstrate the absence of buttresses supporting his thesis on the appropriate distribution of loads which will reduce the "empujos" [horizontal thrusts] to the point of not requiring more than the thickness of the walls to stand (Tosca [1727] 1992, 227-230). The present article2 assesses T osca' s appreciation studying how loads and the thrusts they generate are transmitted through the different masonry elements that constitute this ciborium. In order to do so, we first present a geometrical analysis and make considerations regarding its materials and construction methods to, subsequently, analyze its stability adopting an equilibrium approach within the theoretical framework of the lower bound limit analysis.

Ethics and materials: Some Spanish Cases

The problem is general: modern architects and engineers are trying to understand historic structures using the wrong theoretical frame, the classic (elastic) thery of structures developed in the 19th Century for iron and stell, and in the 20th century for reinforced concrete, disguised with "modern" computer packages, mainly FEM, but also others. Masonry is an essentially different material, and the structural equations must be adapted accordingly. It is not a matter of "taste" or "opinion", and the consequences are before us. Since, say 1920s, historic monuments have suffered the aggression of generations of archietcts and engineers, trying to transform masonry in reinfored concrete or steel. The damage to the monuments and the expense has been, and is, enormous. However, as we have an adequate theory (modern limit analysis of masonry structures, Heyman 1966) which encompasses the "old theory" used successfully by the 18th and 19th Century practical engineers (from Perronet to Sejourné), it is a matter of "Ethics" not to use the wrong approach. It is also "contra natura" to modify the material masonry with indiscriminate injections, stitchings, etc. It is insane to consider, suddenly, that buildings which are Centuries or milennia old, are suddenly in danger of collapse. Maintenance is necessary but not the actual destruction of the constructive essence of the monument. A cocktail of "ignorance, fear and greed" is acting under the best of intentions.

Wedges and plate-bandes: mechanical theories after De la Hire

Plate-bandes are straight masonry arches (they are called, also, flat arches or lintel arches). Ideally they have the surfaces of extrados and intrados plane and horizontal. The stones or bricks have radial joints converging usually in one centre. The voussoirs have the form of wedges and in French they are called "claveaux". A plate-bande is, in fact, a lintel made of several stones and the proportions of lintels and plate-bandes are similar. Proportions of plate-bandes, that is the relationship between the thickness t and the span s (t/s)varies, typically between 1/4–1/3 in thick plate-bandes, and is less than 1/20 in the most slender ones. A ratio of circa 1/8 was usual in the 18th Century and follows a simple geometrical rule: the centre form with the intrados an equilateral triangle and the plate-bande should contain an arc of circle. The joints are usually plane, but in some cases present a «rebated» or «stepped» form. Plate-bandes exert an inclined thrust as any masonry arch. This thrust is usually very high and it requires either massive buttresses, or to be built in the middle of thick walls. Master builders and architects have tried since antiquity to calculate the abutment necessary for any arch. A modern architect or engineer will measure the arch thrust in units of force, kN or tons. Traditionally, the thrust has been measured as the size of the buttresses to resist it safely. Old structural rules, then, addressed the design problem establishing a relationship between the span and the depth of the buttress. These were empirical rules, particular for every type of arch or structure in every epoch. Thus, the typical gothic buttress is 1/4 of the vault span, but a Renaissance or baroque barrel vault will need more than 1/3 of the span. A plate-bande would require more than one half of the span; this is precisely the rule cited by the French engineer Gautier, who tried unsuccessfully to justify it by static reasons. They were used, typically, to form the lintels of windows or doors (1-2 m, typically); in Antiquity they were used, also, though rarely, at the gates of city walls or in niches (ca. 2 m, reaching 5.2 m). Plate-bandes may show particular problems: it is not unusual that some sliding of the voussoirs can be observed, particularly in thick plate-bandes. The stepped joints on Fig. 1, left, were used to avoid this problem. There are other «hidden» methods, like iron cramps or the use of stone wedges, etc. In seismic zones these devices were usual. Another problem relates to the deformation; a slight yielding of the abutments, or even the compression of the mortar joints, may lead to some cracking and the descent of the central keystone. Even a tiny descent will convert the original straight line of the intrados in a broken line with a visible «kink» or angle in the middle. Of course, both problems should be avoided. Finally, the wedge form of the voussoirs lead to acute angles in the stones and this can produce partial fractures; this occurs usually at the inferior border of the springers at the abutments. It follows, that to build a successful plate-bande is not an easy matter. Also, the structural study of plate-bandes is far from simple, and mechanics and geometry are related in a particular way. In the present paper we will concentrate on the structural aspects and their constructive consequences, with a historical approach. We will outline the development of structural analysis of plate-bandes from ca. 1700 until today. This brief history has a more than purely academic interest. Different approaches and theories pointed to particular problem, and though the solution given may have been incorrect, the question posed was often pertinent. The paper ends with the application of modern Limit Analysis of Masonry Structures, developed mainly by professor Heyman in the last fifty years. The work aims, also, to give some clues for the actual architect and engineer involved in the analysis or restoration of masonry buildings.

Bóvedas y contrarresto del Colegio de Nuestra Señora de la Antigua de Monforte de Lemos : geometría, construcción y mecánica

El Colegio de Nuestra Señora de la Antigua de Monforte de Lemos (1593-1619) constituye el principal ejemplo del clasicismo herreriano de Galicia. Esta tesis analiza los aspectos técnicos de la construcción de sus bóvedas de cantería, en el marco disciplinar de la Historia de la Construcción. El Colegio contiene varios tipos de bóvedas entre las que destacan su cúpula trasdosada y su escalera monumental. Además, se conservan los documentos notariales históricos que describen las obras y una de las monteas más importantes de Galicia. Las bóvedas han sufrido movimientos y algunas presentan importantes deformaciones, que han suscitado preocupación sobre su estabilidad. El estudio se organiza en tres epígrafes: geometría, construcción y mecánica, tres aspectos estrechamente interrelacionados en este tipo de estructuras. La parte de estudio geométrico parte de un levantamiento en el que se emplea una estación total sin reflectante, una técnica que permite una medición muy precisa de la forma real de las bóvedas. Como complemento se ha recurrido a técnicas de escáner fotográfico. Con los datos de las mediciones se han elaborado los dibujos de las bóvedas, y analizado sus medidas e irregularidades. Se ha incluido un análisis metrológico por métodos inductivos, que ha puesto de manifiesto la unidad utilizada en la construcción. El estudio se completa con la recopilación de las reglas de dimensionamiento que figuran en los tratados, que se comparan con las dimensiones de las bóvedas analizadas. En el estudio constructivo se registran los despieces y la configuración constructiva y se contrastan, por un lado, con las soluciones que aparecen en los principales tratados de los siglos XVI, XVII y XVIII, y por otro, con las descripciones que se encuentran en los documentos históricos del Colegio. Se ha prestado especial atención al estudio de estos documentos, que han permitido conocer algunos aspectos originales del proceso de construcción de estas bóvedas. En una de las paredes del claustro se ha identificado una montea inédita cuyo levantamiento y análisis se incorpora en esta parte. En el estudio mecánico se analiza la estabilidad dentro del marco teórico del Moderno Análisis Límite de Estructuras de Fábricas. La comprensión del comportamiento estructural, unida a la precisión de los levantamientos, ha permitido explicar los movimientos de las fábricas e inferir la geometría original de algunas de las bóvedas, no siempre evidente en los casos de fuertes deformaciones. Los estudios realizados hasta ahora en bóvedas de los siglos XVI y XVII permiten vislumbrar un panorama de soluciones constructivas más rico que el que recogen los tratados históricos. Las nuevas técnicas de medición sin contacto abren nuevas posibilidades para analizar estas estructuras. Los casos estudiados todavía son escasos y se necesita un conjunto más amplio para poder realizar un análisis general. Esta tesis aporta el estudio de uno de esos casos con el fin de contribuir al mejor conocimiento de la construcción en España en los siglos XVI y XVII. ABSTRACT The Colegio de Nuestra Señora de la Antigua (1593-1619), in Monforte de Lemos, Spain, is the best example in Galicia of the classicism influenced by the Monastery of El Escorial and his architect, Juan de Herrera. This thesis analyzes the technical aspects of the construction of its masonry vaults within the discipline of Construction History. The Colegio exhibits various types of vaults, among which the extradosed ashlar dome and the grand staircase are particularly worth mentioning. In addition, the legal documents containing thorough specifications of the work, as well as one of the construction drawings at actual size —one of the best examples in Galicia—, have both been preserved to this date. Some of the vaults have undergone important deformations that have raised concerns about their stability. The study is organized in three sections: geometry, construction, and mechanics, the three clearly interrelated in this type of structures. The geometrical study starts out with a metric survey using a reflectorless total station, a technique that allows a very precise measurement of the actual shape of the vault. This technique was complemented with the use of a photo-based 3D scanner. The resulting measurements were used to draw the vaults, and analyze their dimensions and irregularities. An inductive metrological analysis, which was able to reveal the exact metric unit utilized during the construction, is included in this section. The section is completed with the gathering of the dimensional rules appearing in the various historical treatises, which are compared with the dimensions of the actual vaults. The construction study deals with the quartering of the structure, whose results are compared, on the one hand, with the approaches appearing in the main treatises of the 16th- 18th centuries, and on the other, with the descriptions in the historical documents of the Colegio itself. Special attention has been paid to the study of these documents, which have revealed some original aspects in the construction process of the vaults. In one of the cloister’s walls, a hitherto unheard-of construction drawing at actual size was found. Its survey and analysis is included in this section. The mechanical study analyzes the Colegio's vaults stability within the theoretical frame of the Modern Limit Analysis of Masonry Structures. The interpretation of the structural behavior, coupled with the precision of the surveys, has allowed a deep understanding of the masonry’s movements and the original geometry of some of the vaults —not always evident in the case of strong deformations— has been inferred. The studies dealing with XVI and XVII century vaults conducted up to this date show a landscape of constructive solutions far richer than the one suggested in the historical treatises. The new contactless measuring techniques offer exciting possibilities for the analysis of these structures. The cases studied are still few and more would be needed before a general analysis can be attempted. This thesis constitutes one such additional case and its goal is to improve our understanding of Spanish construction in the 16th-17th centuries.

El problema del rozamiento en el análisis de estructuras de fábrica mediante modelos de sólidos rígidos

La evaluación de la seguridad de estructuras antiguas de fábrica es un problema abierto.El material es heterogéneo y anisótropo, el estado previo de tensiones difícil de conocer y las condiciones de contorno inciertas. A comienzos de los años 50 se demostró que el análisis límite era aplicable a este tipo de estructuras, considerándose desde entonces como una herramienta adecuada. En los casos en los que no se produce deslizamiento la aplicación de los teoremas del análisis límite estándar constituye una herramienta formidable por su simplicidad y robustez. No es necesario conocer el estado real de tensiones. Basta con encontrar cualquier solución de equilibrio, y que satisfaga las condiciones de límite del material, en la seguridad de que su carga será igual o inferior a la carga real de inicio de colapso. Además esta carga de inicio de colapso es única (teorema de la unicidad) y se puede obtener como el óptimo de uno cualquiera entre un par de programas matemáticos convexos duales. Sin embargo, cuando puedan existir mecanismos de inicio de colapso que impliquen deslizamientos, cualquier solución debe satisfacer tanto las restricciones estáticas como las cinemáticas, así como un tipo especial de restricciones disyuntivas que ligan las anteriores y que pueden plantearse como de complementariedad. En este último caso no está asegurada la existencia de una solución única, por lo que es necesaria la búsqueda de otros métodos para tratar la incertidumbre asociada a su multiplicidad. En los últimos años, la investigación se ha centrado en la búsqueda de un mínimo absoluto por debajo del cual el colapso sea imposible. Este método es fácil de plantear desde el punto de vista matemático, pero intratable computacionalmente, debido a las restricciones de complementariedad 0 y z 0 que no son ni convexas ni suaves. El problema de decisión resultante es de complejidad computacional No determinista Polinomial (NP)- completo y el problema de optimización global NP-difícil. A pesar de ello, obtener una solución (sin garantía de exito) es un problema asequible. La presente tesis propone resolver el problema mediante Programación Lineal Secuencial, aprovechando las especiales características de las restricciones de complementariedad, que escritas en forma bilineal son del tipo y z = 0; y 0; z 0 , y aprovechando que el error de complementariedad (en forma bilineal) es una función de penalización exacta. Pero cuando se trata de encontrar la peor solución, el problema de optimización global equivalente es intratable (NP-difícil). Además, en tanto no se demuestre la existencia de un principio de máximo o mínimo, existe la duda de que el esfuerzo empleado en aproximar este mínimo esté justificado. En el capítulo 5, se propone hallar la distribución de frecuencias del factor de carga, para todas las soluciones de inicio de colapso posibles, sobre un sencillo ejemplo. Para ello, se realiza un muestreo de soluciones mediante el método de Monte Carlo, utilizando como contraste un método exacto de computación de politopos. El objetivo final es plantear hasta que punto está justificada la busqueda del mínimo absoluto y proponer un método alternativo de evaluación de la seguridad basado en probabilidades. Las distribuciones de frecuencias, de los factores de carga correspondientes a las soluciones de inicio de colapso obtenidas para el caso estudiado, muestran que tanto el valor máximo como el mínimo de los factores de carga son muy infrecuentes, y tanto más, cuanto más perfecto y contínuo es el contacto. Los resultados obtenidos confirman el interés de desarrollar nuevos métodos probabilistas. En el capítulo 6, se propone un método de este tipo basado en la obtención de múltiples soluciones, desde puntos de partida aleatorios y calificando los resultados mediante la Estadística de Orden. El propósito es determinar la probabilidad de inicio de colapso para cada solución.El método se aplica (de acuerdo a la reducción de expectativas propuesta por la Optimización Ordinal) para obtener una solución que se encuentre en un porcentaje determinado de las peores. Finalmente, en el capítulo 7, se proponen métodos híbridos, incorporando metaheurísticas, para los casos en que la búsqueda del mínimo global esté justificada. Abstract Safety assessment of the historic masonry structures is an open problem. The material is heterogeneous and anisotropic, the previous state of stress is hard to know and the boundary conditions are uncertain. In the early 50's it was proven that limit analysis was applicable to this kind of structures, being considered a suitable tool since then. In cases where no slip occurs, the application of the standard limit analysis theorems constitutes an excellent tool due to its simplicity and robustness. It is enough find any equilibrium solution which satisfy the limit constraints of the material. As we are certain that this load will be equal to or less than the actual load of the onset of collapse, it is not necessary to know the actual stresses state. Furthermore this load for the onset of collapse is unique (uniqueness theorem), and it can be obtained as the optimal from any of two mathematical convex duals programs However, if the mechanisms of the onset of collapse involve sliding, any solution must satisfy both static and kinematic constraints, and also a special kind of disjunctive constraints linking the previous ones, which can be formulated as complementarity constraints. In the latter case, it is not guaranted the existence of a single solution, so it is necessary to look for other ways to treat the uncertainty associated with its multiplicity. In recent years, research has been focused on finding an absolute minimum below which collapse is impossible. This method is easy to set from a mathematical point of view, but computationally intractable. This is due to the complementarity constraints 0 y z 0 , which are neither convex nor smooth. The computational complexity of the resulting decision problem is "Not-deterministic Polynomialcomplete" (NP-complete), and the corresponding global optimization problem is NP-hard. However, obtaining a solution (success is not guaranteed) is an affordable problem. This thesis proposes solve that problem through Successive Linear Programming: taking advantage of the special characteristics of complementarity constraints, which written in bilinear form are y z = 0; y 0; z 0 ; and taking advantage of the fact that the complementarity error (bilinear form) is an exact penalty function. But when it comes to finding the worst solution, the (equivalent) global optimization problem is intractable (NP-hard). Furthermore, until a minimum or maximum principle is not demonstrated, it is questionable that the effort expended in approximating this minimum is justified. XIV In chapter 5, it is proposed find the frequency distribution of the load factor, for all possible solutions of the onset of collapse, on a simple example. For this purpose, a Monte Carlo sampling of solutions is performed using a contrast method "exact computation of polytopes". The ultimate goal is to determine to which extent the search of the global minimum is justified, and to propose an alternative approach to safety assessment based on probabilities. The frequency distributions for the case study show that both the maximum and the minimum load factors are very infrequent, especially when the contact gets more perfect and more continuous. The results indicates the interest of developing new probabilistic methods. In Chapter 6, is proposed a method based on multiple solutions obtained from random starting points, and qualifying the results through Order Statistics. The purpose is to determine the probability for each solution of the onset of collapse. The method is applied (according to expectations reduction given by the Ordinal Optimization) to obtain a solution that is in a certain percentage of the worst. Finally, in Chapter 7, hybrid methods incorporating metaheuristics are proposed for cases in which the search for the global minimum is justified.