Boundary element approach to the dynamic stiffness functions of circular foundations

A boundary element approach for time harmonic axisymmetric problems using the complete space point load fundamental solution is presented. The fundamental solution is integrated numerically along the azimuthal co-ordinate of each axisymmetric element. To increase the accuracy of the numerical integration a simple co-ordinate transformation is proposed. The approach is applied to the computation of the dynamic stiffness functions of rigid circular foundations on layered viscoelastic soils. Three different sites are considered: a uniform half-space, a soil layer on a half-space, and a soil consisting of four horizontal layers and a compliant half-space. The numerical results obtained by the proposed approach for surface circular foundations are very close to corresponding published results obtained by different procedures.

Dynamic impedances of bridge abutments

The soil-structure interaction at bridge abutments may introduce important changes in the dynamic properties of short to medium span bridges. The paper presents the results obtained, through the use of the Boundary Element Method (B.E.M.) technique in several typical situations, including semiinfinite and layered media. Both stiffness and damping properties are included.

Combining a Similar Coefficient Identification Algorithm with the Boundary Element Method

The boundary element method (BEM) has been applied successfully to many engineering problems during the last decades. Compared with domain type methods like the finite element method (FEM) or the finite difference method (FDM) the BEM can handle problems where the medium extends to infinity much easier than domain type methods as there is no need to develop special boundary conditions (quiet or absorbing boundaries) or infinite elements at the boundaries introduced to limit the domain studied. The determination of the dynamic stiffness of arbitrarily shaped footings is just one of these fields where the BEM has been the method of choice, especially in the 1980s. With the continuous development of computer technology and the available hardware equipment the size of the problems under study grew and, as the flop count for solving the resulting linear system of equations grows with the third power of the number of equations, there was a need for the development of iterative methods with better performance. In [1] the GMRES algorithm was presented which is now widely used for implementations of the collocation BEM. While the FEM results in sparsely populated coefficient matrices, the BEM leads, in general, to fully or densely populated ones, depending on the number of subregions, posing a serious memory problem even for todays computers. If the geometry of the problem permits the surface of the domain to be meshed with equally shaped elements a lot of the resulting coefficients will be calculated and stored repeatedly. The present paper shows how these unnecessary operations can be avoided reducing the calculation time as well as the storage requirement. To this end a similar coefficient identification algorithm (SCIA), has been developed and implemented in a program written in Fortran 90. The vertical dynamic stiffness of a single pile in layered soil has been chosen to test the performance of the implementation. The results obtained with the 3-d model may be compared with those obtained with an axisymmetric formulation which are considered to be the reference values as the mesh quality is much better. The entire 3D model comprises more than 35000 dofs being a soil region with 21168 dofs the biggest single region. Note that the memory necessary to store all coefficients of this single region is about 6.8 GB, an amount which is usually not available with personal computers. In the problem under study the interface zone between the two adjacent soil regions as well as the surface of the top layer may be meshed with equally sized elements. In this case the application of the SCIA leads to an important reduction in memory requirements. The maximum memory used during the calculation has been reduced to 1.2 GB. The application of the SCIA thus permits problems to be solved on personal computers which otherwise would require much more powerful hardware.

Vanadium-Doped In and Sn Sulphides: Photocatalysts able to use the whole visible light spectrum

Using photocatalysis for energy applications depends, more than for environmental purposes or selective chemical synthesis, on converting as much of the solar spectrum as possible; the best photocatalyst, titania, is far from this. Many efforts are pursued to use better that spectrum in photocatalysis, by doping titania or using other materials (mainly oxides, nitrides and sulphides) to obtain a lower bandgap, even if this means decreasing the chemical potential of the electron-hole pairs. Here we introduce an alternative scheme, using an idea recently proposed for photovoltaics: the intermediate band (IB) materials. It consists in introducing in the gap of a semiconductor an intermediate level which, acting like a stepstone, allows an electron jumping from the valence band to the conduction band in two steps, each one absorbing one sub-bandgap photon. For this the IB must be partially filled, to allow both sub-bandgap transitions to proceed at comparable rates; must be made of delocalized states to minimize nonradiative recombination; and should not communicate electronically with the outer world. For photovoltaic use the optimum efficiency so achievable, over 1.5 times that given by a normal semiconductor, is obtained with an overall bandgap around 2.0 eV (which would be near-optimal also for water phtosplitting). Note that this scheme differs from the doping principle usually considered in photocatalysis, which just tries to decrease the bandgap; its aim is to keep the full bandgap chemical potential but using also lower energy photons. In the past we have proposed several IB materials based on extensively doping known semiconductors with light transition metals, checking first of all with quantum calculations that the desired IB structure results. Subsequently we have synthesized in powder form two of them: the thiospinel In2S3 and the layered compound SnS2 (having bandgaps of 2.0 and 2.2 eV respectively) where the octahedral cation is substituted at a â?10% level with vanadium, and we have verified that this substitution introduces in the absorption spectrum the sub-bandgap features predicted by the calculations. With these materials we have verified, using a simple reaction (formic acid oxidation), that the photocatalytic spectral response is indeed extended to longer wavelengths, being able to use even 700 nm photons, without largely degrading the response for above-bandgap photons (i.e. strong recombination is not induced) [3b, 4]. These materials are thus promising for efficient photoevolution of hydrogen from water; work on this is being pursued, the results of which will be presented.

Open-atmosphere structural depth profiling of multilayer samples of photovoltaic interest using laser-induced plasma spectrometry

The present work aims to assess Laser-Induced Plasma Spectrometry (LIPS) as a tool for the characterization of photovoltaic materials. Despite being a well-established technique with applications to many scientific and industrial fields, so far LIPS is little known to the photovoltaic scientific community. The technique allows the rapid characterization of layered samples without sample preparation, in open atmosphere and in real time. In this paper, we assess LIPS ability for the determination of elements that are difficult to analyze by other broadly used techniques, or for producing analytical information from very low-concentration elements. The results of the LIPS characterization of two different samples are presented: 1) a 90 nm, Al-doped ZnO layer deposited on a Si substrate by RF sputtering and 2) a Te-doped GaInP layer grown on GaAs by Metalorganic Vapor Phase Epitaxy. For both cases, the depth profile of the constituent and dopant elements is reported along with details of the experimental setup and the optimization of key parameters. It is remarkable that the longest time of analysis was ∼10 s, what, in conjunction with the other characteristics mentioned, makes of LIPS an appealing technique for rapid screening or quality control whether at the lab or at the production line.

New Intermediate band sulphide nanoparticles acting in the full visible light range spectra as an active photocatalyst

Nowadays one of the challenges of materials science is to find new technologies that will be able to make the most of renewable energies. An example of new proposals in this field are the intermediate-band (IB) materials, which promise higher efficiencies in photovoltaic applications (through the intermediate band solar cells), or in heterogeneous photocatalysis (using nanoparticles of them, for the light-induced degradation of pollutants or for the efficient photoevolution of hydrogen from water). An IB material consists in a semiconductor in which gap a new level is introduced [1], the intermediate band (IB), which should be partially filled by electrons and completely separated of the valence band (VB) and of the conduction band (CB). This scheme (figure 1) allows an electron from the VB to be promoted to the IB, and from the latter to the CB, upon absorption of photons with energy below the band gap Eg, so that energy can be absorbed in a wider range of the solar spectrum and a higher current can be obtained without sacrificing the photovoltage (or the chemical driving force) corresponding to the full bandgap Eg, thus increasing the overall efficiency. This concept, applied to photocatalysis, would allow using photons of a wider visible range while keeping the same redox capacity. It is important to note that this concept differs from the classic photocatalyst doping principle, which essentially tries just to decrease the bandgap. This new type of materials would keep the full bandgap potential but would use also lower energy photons. In our group several IB materials have been proposed, mainly for the photovoltaic application, based on extensively doping known semiconductors with transition metals [2], examining with DFT calculations their electronic structures. Here we refer to In2S3 and SnS2, which contain octahedral cations; when doped with Ti or V an IB is formed according to quantum calculations (see e.g. figure 2). We have used a solvotermal synthesis method to prepare in nanocrystalline form the In2S3 thiospinel and the layered compound SnS2 (which when undoped have bandgaps of 2.0 and 2.2 eV respectively) where the cation is substituted by vanadium at a ?10% level. This substitution has been studied, characterizing the materials by different physical and chemical techniques (TXRF, XRD, HR-TEM/EDS) (see e.g. figure 3) and verifying with UV spectrometry that this substitution introduces in the spectrum the sub-bandgap features predicted by the calculations (figure 4). For both sulphide type nanoparticles (doped and undoped) the photocatalytic activity was studied by following at room temperature the oxidation of formic acid in aqueous suspension, a simple reaction which is easily monitored by UV-Vis spectroscopy. The spectral response of the process is measured using a collection of band pass filters that allow only some wavelengths into the reaction system. Thanks to this method the spectral range in which the materials are active in the photodecomposition (which coincides with the band gap for the undoped samples) can be checked, proving that for the vanadium substituted samples this range is increased, making possible to cover all the visible light range. Furthermore it is checked that these new materials are more photocorrosion resistant than the toxic CdS witch is a well know compound frequently used in tests of visible light photocatalysis. These materials are thus promising not only for degradation of pollutants (or for photovoltaic cells) but also for efficient photoevolution of hydrogen from water; work in this direction is now being pursued.

Resonant absorption in a plasma step profile

Resonant absorption of p-polarized light shined on a plane-layered plasma with a step profile, is discussed as a function of wavelength (or critical density n,) of the light: for simplicity the incidence angle is assumed small. If n, lies within or above the step, the absorption A is given by Ginzburg’s result modified by strong reflections at the foot and top of the step. The absorption above is total for particular values of nc and U. For n, crossing the top of the density step the absorption is not monotonical: it exhibits a minimum that vanishes for zero radius of curvature U there and zero collision frequency 1’ (A - Iln VI-’). The results are applied to the profile produced by irradiating a solid target with a high-intensity pulse that steepens the plasma by radiation pressure.

Estudio teórico de propiedades termodinámicas y optoelectronicas de nuevos materiales fotovoltaicos de banda intermedia = Theoretical study of thermodynamic and optoelectronic properties of new intermediate band photovoltaic materials

Los materiales de banda intermedia han atraido la atención de la comunidad científica en el campo de la energía solar fotovoltaica en los últimos años. Sin embargo, con el objetivo de entender los fundamentos de las células solares de banda intermedia, se debe llevar a cabo un estudio profundo de la características de los materiales. Esto se puede hacer mediante un modelo teórico usando Primeros Principios. A partir de este enfoque se pueden obtener resultados tales como la estructura electrónica y propiedades ópticas, entre otras, de los semiconductores fuertemente dopados y sus precursores. Con el fin de desentrañar las estructuras de estos sistemas electrónicos, esta tesis presenta un estudio termodinámico y optoelectrónico de varios materiales fotovoltaicos. Específicamente se caracterizaron los materiales avanzados de banda intermedia y sus precursores. El estudio se hizo en términos de caracterización teórica de la estructura electrónica, la energética del sistema, entre otros. Además la estabilidad se obtuvo usando configuraciones adaptadas a la simetría del sistema y basado en la combinatoria. Las configuraciones de los sitios ocupados por defectos permiten obtener información sobre un espacio de configuraciones donde las posiciones de los dopantes sustituidos se basan en la simetría del sólido cristalino. El resultado puede ser tratado usando elementos de termodinámica estadística y da información de la estabilidad de todo el espacio simétrico. Además se estudiaron otras características importantes de los semiconductores de base. En concreto, el análisis de las interacciones de van der Waals fueron incluidas en el semiconductor en capas SnS2, y el grado de inversión en el caso de las espinelas [M]In2S4. En este trabajo además realizamos una descripción teórica exhaustiva del sistema CdTe:Bi. Este material de banda-intermedia muestra características que son distintas a las de los otros materiales estudiados. También se analizó el Zn como agente modulador de la posición de las sub-bandas prohibidas en el material de banda-intermedia CuGaS2:Ti. Analizándose además la viabilidad termodinámica de la formación de este compuesto. Finalmente, también se describió el GaN:Cr como material de banda intermedia, en la estructura zinc-blenda y en wurtztite, usando configuraciones de sitios ocupados de acuerdo a la simetría del sistema cristalino del semiconductor de base. Todos los resultados, siempre que fue posible, fueron comparados con los resultados experimentales. ABSTRACT The intermediate-band materials have attracted the attention of the scientific community in the field of the photovoltaics in recent years. Nevertheless, in order to understand the intermediate-band solar cell fundamentals, a profound study of the characteristics of the materials is required. This can be done using theoretical modelling from first-principles. The electronic structure and optical properties of heavily doped semiconductors and their precursor semiconductors are, among others, results that can be obtained from this approach. In order to unravel the structures of these crystalline systems, this thesis presents a thermodynamic and optoelectronic study of several photovoltaic materials. Specifically advanced intermediate-band materials and their precursor semiconductors were characterized. The study was made in terms of theoretical characterization of the electronic structure, energetics among others. The stability was obtained using site-occupancy-disorder configurations adapted to the symmetry of the system and based on combinatorics. The site-occupancy-disorder method allows the formation of a configurational space of substitutional dopant positions based on the symmetry of the crystalline solid. The result, that can be treated using statistical thermodynamics, gives information of the stability of the whole space of symmetry of the crystalline lattice. Furthermore, certain other important characteristics of host semiconductors were studied. Specifically, the van der Waal interactions were included in the SnS2 layered semiconductor, and the inversion degree in cases of [M]In2S4 spinels. In this work we also carried out an exhaustive theoretical description of the CdTe:Bi system. This intermediate-band material shows characteristics that are distinct from those of the other studied intermediate-band materials. In addition, Zn was analysed as a modulator of the positions of the sub-band gaps in the CuGaS2:Ti intermediate-band material. The thermodynamic feasibility of the formation of this compound was also carried out. Finally GaN:Cr intermediate-band material was also described both in the zinc-blende and the wurtztite type structures, using the symmetry-adapted-space of configurations. All results, whenever possible, were compared with experimental results.

Quantitative evaluation of patient-specific conforming hexahedral meshes of abdominal aortic aneurysms and intraluminal thrombus generated from MRI

A novel method for generating patient-specific high quality conforming hexahedral meshes is presented. The meshes are directly obtained from the segmentation of patient magnetic resonance (MR) images of abdominal aortic aneu-rysms (AAA). The MRI permits distinguishing between struc-tures of interest in soft tissue. Being so, the contours of the lumen, the aortic wall and the intraluminal thrombus (ILT) are available and thus the meshes represent the actual anato-my of the patient?s aneurysm, including the layered morpholo-gies of these structures. Most AAAs are located in the lower part of the aorta and the upper section of the iliac arteries, where the inherent tortuosity of the anatomy and the presence of the ILT makes the generation of high-quality elements at the bifurcation is a challenging task. In this work we propose a novel approach for building quadrilateral meshes for each surface of the sectioned geometry, and generating conforming hexahedral meshes by combining the quadrilateral meshes. Conforming hexahedral meshes are created for the wall and the ILT. The resulting elements are evaluated on four patients? datasets using the Scaled Jacobian metric. Hexahedral meshes of 25,000 elements with 94.8% of elements well-suited for FE analysis are generated.

Optimum high temperature strength of two-dimensional nanocomposites

High-temperature nanoindentation was used to reveal nano-layer size effects on the hardness of two-dimensional metallic nanocomposites. We report the existence of a critical layer thickness at which strength achieves optimal thermal stability. Transmission electron microscopy and theoretical bicrystal calculations show that this optimum arises due to a transition from thermally activated glide within the layers to dislocation transmission across the layers. We demonstrate experimentally that the atomic-scale properties of the interfaces profoundly affect this critical transition. The strong implications are that interfaces can be tuned to achieve an optimum in high temperature strength in layered nanocomposite structures.

GMonE: a complete approach to cloud monitoring

The inherent complexity of modern cloud infrastructures has created the need for innovative monitoring approaches, as state-of-the-art solutions used for other large-scale environments do not address specific cloud features. Although cloud monitoring is nowadays an active research field, a comprehensive study covering all its aspects has not been presented yet. This paper provides a deep insight into cloud monitoring. It proposes a unified cloud monitoring taxonomy, based on which it defines a layered cloud monitoring architecture. To illustrate it, we have implemented GMonE, a general-purpose cloud monitoring tool which covers all aspects of cloud monitoring by specifically addressing the needs of modern cloud infrastructures. Furthermore, we have evaluated the performance, scalability and overhead of GMonE with Yahoo Cloud Serving Benchmark (YCSB), by using the OpenNebula cloud middleware on the Grid’5000 experimental testbed. The results of this evaluation demonstrate the benefits of our approach, surpassing the monitoring performance and capabilities of cloud monitoring alternatives such as those present in state-of-the-art systems such as Amazon EC2 and OpenNebula.

Opportunities and challenges in the use of inorganic fullerene-like nanoparticles to produce advanced polymer nanocomposites

Polymer/inorganic nanoparticle nanocomposites have garnered considerable academic and industrial interest over recent decades in the development of advanced materials for a wide range of applications. In this respect, the dispersion of so-called inorganic fullerene-like (IF) nanoparticles, e.g., tungsten disulfide (IF-WS2) or molybdenum disulfide (IF-MoS2), into polymeric matrices is emerging as a new strategy. The surprising properties of these layered metal dichalcogenides such as high impact resistance and superior tribological behavior, attributed to their nanoscale size and hollow quasi-spherical shape, open up a wide variety of opportunities for applications of these inorganic compounds. The present work presents a detailed overview on research in the area of IF-based polymer nanocomposites, with special emphasis on the use of IF-WS2 nanoparticles as environmentally friendly reinforcing fillers. The incorporation of IF particles has been shown to be efficient for improving thermal, mechanical and tribological properties of various thermoplastic polymers, such as polypropylene, nylon-6, poly(phenylene sulfide), poly(ether ether ketone), where nanocomposites were fabricated by simple melt-processing routes without the need for modifiers or surfactants. This new family of nanocomposites exhibits similar or enhanced performance when compared with nanocomposites that incorporate carbon nanotubes, carbon nanofibers or nanoclays, but are substantially more cost-effective, efficient and environmentally satisfactory. Most recently, innovative approaches have been described that exploit synergistic effects to produce new materials with enhanced properties, including the combined use of micro- and nanoparticles such as IF-WS2/nucleating agent or IF-WS2/carbon fiber, as well as dual nanoparticle systems such as SWCNT/IF-WS2 where each nanoparticle has different characteristics. The structure–property relationships of these nanocomposites are discussed and potential applications proposed ranging from medicine to the aerospace, automotive and electronics industries.

Eigenvector centrality of nodes in multiplex networks

We extend the concept of eigenvector centrality to multiplex networks, and introduce several alternative parameters that quantify the importance of nodes in a multi-layered networked system, including the definition of vectorial-type centralities. In addition, we rigorously show that, under reasonable conditions, such centrality measures exist and are unique. Computer experiments and simulations demonstrate that the proposed measures provide substantially different results when applied to the same multiplex structure, and highlight the non-trivial relationships between the different measures of centrality introduced.

Operational Aspects of Full Reduction in Lambda Calculi

Esta tesis estudia la reducción plena (‘full reduction’ en inglés) en distintos cálculos lambda. 1 En esencia, la reducción plena consiste en evaluar los cuerpos de las funciones en los lenguajes de programación funcional con ligaduras. Se toma el cálculo lambda clásico (i.e., puro y sin tipos) como el sistema formal que modela el paradigma de programación funcional. La reducción plena es una técnica fundamental cuando se considera a los programas como datos, por ejemplo para la optimización de programas mediante evaluación parcial, o cuando algún atributo del programa se representa a su vez por un programa, como el tipo en los demostradores automáticos de teoremas actuales. Muchas semánticas operacionales que realizan reducción plena tienen naturaleza híbrida. Se introduce formalmente la noción de naturaleza híbrida, que constituye el hilo conductor de todo el trabajo. En el cálculo lambda la naturaleza híbrida se manifiesta como una ‘distinción de fase’ en el tratamiento de las abstracciones, ya sean consideradas desde fuera o desde dentro de si mismas. Esta distinción de fase conlleva una estructura en capas en la que una semántica híbrida depende de una o más semánticas subsidiarias. Desde el punto de vista de los lenguajes de programación, la tesis muestra como derivar, mediante técnicas de transformación de programas, implementaciones de semánticas operacionales que reducen plenamente a partir de sus especificaciones. Las técnicas de transformación de programas consisten en transformaciones sintácticas que preservan la equivalencia semántica de los programas. Se ajustan las técnicas de transformación de programas existentes para trabajar con implementaciones de semánticas híbridas. Además, se muestra el impacto que tiene la reducción plena en las implementaciones que utilizan entornos. Los entornos son un ingrediente fundamental en las implementaciones realistas de una máquina abstracta. Desde el punto de vista de los sistemas formales, la tesis desvela una teoría novedosa para el cálculo lambda con paso por valor (‘call-by-value lambda calculus’ en inglés) que es consistente con la reducción plena. Dicha teoría induce una noción de equivalencia observacional que distingue más puntos que las teorías existentes para dicho cálculo. Esta contribución ayuda a establecer una ‘teoría estándar’ en el cálculo lambda con paso por valor que es análoga a la ‘teoría estándar’ del cálculo lambda clásico propugnada por Barendregt. Se presentan resultados de teoría de la demostración, y se sugiere como abordar el estudio de teoría de modelos. ABSTRACT This thesis studies full reduction in lambda calculi. In a nutshell, full reduction consists in evaluating the body of the functions in a functional programming language with binders. The classical (i.e., pure untyped) lambda calculus is set as the formal system that models the functional paradigm. Full reduction is a prominent technique when programs are treated as data objects, for instance when performing optimisations by partial evaluation, or when some attribute of the program is represented by a program itself, like the type in modern proof assistants. A notable feature of many full-reducing operational semantics is its hybrid nature, which is introduced and which constitutes the guiding theme of the thesis. In the lambda calculus, the hybrid nature amounts to a ‘phase distinction’ in the treatment of abstractions when considered either from outside or from inside themselves. This distinction entails a layered structure in which a hybrid semantics depends on one or more subsidiary semantics. From a programming languages standpoint, the thesis shows how to derive implementations of full-reducing operational semantics from their specifications, by using program transformations techniques. The program transformation techniques are syntactical transformations which preserve the semantic equivalence of programs. The existing program transformation techniques are adjusted to work with implementations of hybrid semantics. The thesis also shows how full reduction impacts the implementations that use the environment technique. The environment technique is a key ingredient of real-world implementations of abstract machines which helps to circumvent the issue with binders. From a formal systems standpoint, the thesis discloses a novel consistent theory for the call-by-value variant of the lambda calculus which accounts for full reduction. This novel theory entails a notion of observational equivalence which distinguishes more points than other existing theories for the call-by-value lambda calculus. This contribution helps to establish a ‘standard theory’ in that calculus which constitutes the analogous of the ‘standard theory’ advocated by Barendregt in the classical lambda calculus. Some prooftheoretical results are presented, and insights on the model-theoretical study are given.

La estabilidad del frente de excavación de un túnel en terrenos heterogéneos o con un criterio de rotura no-lineal

El frente de un túnel puede colapsar si la presión aplicada sobre el es inferior a un valor limite denominado presión “critica” o “de colapso”. En este trabajo se desarrolla y presenta un mecanismo de rotura rotacional generado punto a punto para el cálculo de la presión de colapso del frente de túneles excavados en terrenos estratificados o en materiales que siguen un criterio de rotura nolineal. La solución propuesta es una solución de contorno superior en el marco del Análisis Límite y supone una generalización del mecanismo de rotura mas reciente existente en la bibliografía. La presencia de un terreno estratificado o con un criterio de rotura no-lineal implica una variabilidad espacial de las propiedades resistentes. Debido a esto, se generaliza el mecanismo desarrollado por Mollon et al. (2011b) para suelos, de tal forma que se puedan considerar valores locales del ángulo de rozamiento y de la cohesión. Además, la estratificación del terreno permite una rotura parcial del frente, por lo que se implementa esta posibilidad en el mecanismo, siendo la primera solución que emplea un mecanismo de rotura que se ajusta a la estratigrafía del terreno. Por otro lado, la presencia de un material con un criterio de rotura no-lineal exige introducir en el modelo, como variable de estudio, el estado tensional en el frente, el cual se somete al mismo proceso de optimización que las variables geométricas del mecanismo. Se emplea un modelo numérico 3D para validar las predicciones del mecanismo de Análisis Limite, demostrando que proporciona, con un esfuerzo computacional significativamente reducido, buenas predicciones de la presión critica, del tipo de rotura (global o parcial) en terrenos estratificados y de la geometría de fallo. El mecanismo validado se utiliza para realizar diferentes estudios paramétricos sobre la influencia de la estratigrafía en la presión de colapso. Igualmente, se emplea para elaborar cuadros de diseño de la presión de colapso para túneles ejecutados con tuneladora en macizos rocosos de mala calidad y para analizar la influencia en la estabilidad del frente del método constructivo. Asimismo, se lleva a cabo un estudio de fiabilidad de la estabilidad del frente de un túnel excavado en un macizo rocoso altamente fracturado. A partir de el se analiza como afectan las diferentes hipótesis acerca de los tipos de distribución y de las estructuras de correlación a los resultados de fiabilidad. Se investiga también la sensibilidad de los índices de fiabilidad a los cambios en las variables aleatorias, identificando las mas relevantes para el diseño. Por ultimo, se lleva a cabo un estudio experimental mediante un modelo de laboratorio a escala reducida. El modelo representa medio túnel, lo cual permite registrar el movimiento del material mediante una técnica de correlación de imágenes fotográficas. El ensayo se realiza con una arena seca y se controla por deformaciones mediante un pistón que simula el frente. Los resultados obtenidos se comparan con las estimaciones de la solución de Análisis Límite, obteniéndose un ajuste razonable, de acuerdo a la literatura, tanto en la geometría de rotura como en la presión de colapso. A tunnel face may collapse if the applied support pressure is lower than a limit value called the ‘critical’ or ‘collapse’ pressure. In this work, an advanced rotational failure mechanism generated ‘‘point-by-point” is developed to compute the collapse pressure for tunnel faces in layered (or stratified) grounds or in materials that follow a non-linear failure criterion. The proposed solution is an upper bound solution in the framework of limit analysis which extends the most advanced face failure mechanism in the literature. The excavation of the tunnel in a layered ground or in materials with a non-linear failure criterion may lead to a spatial variability of the strength properties. Because of this, the rotational mechanism recently proposed by Mollon et al. (2011b) for Mohr-Coulomb soils is generalized so that it can consider local values of the friction angle and of the cohesion. For layered soils, the mechanism needs to be extended to consider the possibility for partial collapse. The proposed methodology is the first solution with a partial collapse mechanism that can fit to the stratification. Similarly, the use of a nonlinear failure criterion introduces the need to introduce new parameters in the optimization problem to consider the distribution of normal stresses along the failure surface. A 3D numerical model is employed to validate the predictions of the limit analysis mechanism, demonstrating that it provides, with a significantly reduced computational effort, good predictions of critical pressure, of the type of collapse (global or partial) in layered soils, and of its geometry. The mechanism is then employed to conduct parametric studies of the influence of several geometrical and mechanical parameters on face stability of tunnels in layered soils. Similarly, the methodology has been further employed to develop simple design charts that provide the face collapse pressure of tunnels driven by TBM in low quality rock masses and to study the influence of the construction method. Finally, a reliability analysis of the stability of a tunnel face driven in a highly fractured rock mass is performed. The objective is to analyze how different assumptions about distributions types and correlation structures affect the reliability results. In addition, the sensitivity of the reliability index to changes in the random variables is studied, identifying the most relevant variables for engineering design. Finally, an experimental study is carried out using a small-scale laboratory model. The problem is modeled in half, cutting through the tunnel axis vertically, so that displacements of soil particles can be recorded by a digital image correlation technique. The tests were performed with dry sand and displacements are controlled by a piston that supports the soil. The results of the model are compared with the predictions of the Limit Analysis mechanism. A reasonable agreement, according to literature, is obtained between the shapes of the failure surfaces and between the collapse pressures observed in the model tests and computed with the analytical solution.