Performance of Gabion Faced Reinforced Earth Retaining Walls

Gabion faced re.taining walls are essentially semi rigid structures that can generally accommodate large lateral and vertical movements without excessive structural distress. Because of this inherent feature, they offer technical and economical advantage over the conventional concrete gravity retaining walls. Although they can be constructed either as gravity type or reinforced soil type, this work mainly deals with gabion faced reinforced earth walls as they are more suitable to larger heights. The main focus of the present investigation was the development of a viable plane strain two dimensional non linear finite element analysis code which can predict the stress - strain behaviour of gabion faced retaining walls - both gravity type and reinforced soil type. The gabion facing, backfill soil, In - situ soil and foundation soil were modelled using 20 four noded isoparametric quadrilateral elements. The confinement provided by the gabion boxes was converted into an induced apparent cohesion as per the membrane correction theory proposed by Henkel and Gilbert (1952). The mesh reinforcement was modelled using 20 two noded linear truss elements. The interactions between the soil and the mesh reinforcement as well as the facing and backfill were modelled using 20 four noded zero thickness line interface elements (Desai et al., 1974) by incorporating the nonlinear hyperbolic formulation for the tangential shear stiffness. The well known hyperbolic formulation by Ouncan and Chang (1970) was used for modelling the non - linearity of the soil matrix. The failure of soil matrix, gabion facing and the interfaces were modelled using Mohr - Coulomb failure criterion. The construction stages were also modelled.Experimental investigations were conducted on small scale model walls (both in field as well as in laboratory) to suggest an alternative fill material for the gabion faced retaining walls. The same were also used to validate the finite element programme developed as a part of the study. The studies were conducted using different types of gabion fill materials. The variation was achieved by placing coarse aggregate and quarry dust in different proportions as layers one above the other or they were mixed together in the required proportions. The deformation of the wall face was measured and the behaviour of the walls with the variation of fill materials was analysed. It was seen that 25% of the fill material in gabions can be replaced by a soft material (any locally available material) without affecting the deformation behaviour to large extents. In circumstances where deformation can be allowed to some extents, even up to 50% replacement with soft material can be possible.The developed finite element code was validated using experimental test results and other published results. Encouraged by the close comparison between the theory and experiments, an extensive and systematic parametric study was conducted, in order to gain a closer understanding of the behaviour of the system. Geometric parameters as well as material parameters were varied to understand their effect on the behaviour of the walls. The final phase of the study consisted of developing a simplified method for the design of gabion faced retaining walls. The design was based on the limit state method considering both the stability and deformation criteria. The design parameters were selected for the system and converted to dimensionless parameters. Thus the procedure for fixing the dimensions of the wall was simplified by eliminating the conventional trial and error procedure. Handy design charts were developed which would prove as a hands - on - tool to the design engineers at site. Economic studies were also conducted to prove the cost effectiveness of the structures with respect to the conventional RCC gravity walls and cost prediction models and cost breakdown ratios were proposed. The studies as a whole are expected to contribute substantially to understand the actual behaviour of gabion faced retaining wall systems with particular reference to the lateral deformations.

Extensional rheology and stability of gas cell walls in bread doughs at elevated temperatures in relation to breadmaking performance

The rheological properties of gas cell walls in bread doughs are considered to be important in relation to their stability and gas retention during proof and baking. Large deformation rheological properties of gas cell walls were measured using biaxial extension for a number of doughs of varying breadmaking quality at constant strain rate and elevated temperatures of 25-60degreesC. Strain hardening and failure strain of cell walls both decreased with temperature, with cell walls in good breadmaking doughs remaining stable and retaining their strain hardening properties at higher temperatures (60degreesC), while the cell walls of poor breadmaking doughs became unstable at lower temperatures (45-50degreesC) and had lower strain hardening. Strain hardening measured at 50degreesC gave good correlations with baking volume, with the best correlations achieved between rheological measurements and baking tests that used similar mixing conditions. As predicted by the considered failure criterion, a strain hardening value of I defines a region below which gas cell walls become unstable, and discriminates well between the baking quality of a range of commercial flour blends of varying quality. This indicates that the stability of gas cell walls during baking is strongly related to strain hardening properties, and that extensional rheological measurements can be used as indicators of baking quality.

Rheology and the breadmaking process

The applications of rheology to the main processes encountered during breadmaking (mixing, sheeting, fermentation and baking) are reviewed. The most commonly used rheological test methods and their relationships to product functionality are reviewed. It is shown that the most commonly used method for rheological testing of doughs, shear oscillation dynamic rheology, is generally used under deformation conditions inappropriate for breadmaking and shows little relationship with end-use performance. The frequency range used in conventional shear oscillation tests is limited to the plateau region, which is insensitive to changes in the HMW glutenin polymers thought to be responsible for variations in baking quality. The appropriate deformation conditions can be accessed either by long-time creep or relaxation measurements, or by large deformation extensional measurements at low strain rates and elevated temperatures. Molecular size and structure of the gluten polymers that make up the major structural components of wheat are related to their rheological properties via modern polymer rheology concepts. Interactions between polymer chain entanglements and branching are seen to be the key mechanisms determining the rheology of HMW polymers. Recent work confirms the observation that the dynamic shear plateau modulus is essentially independent of variations in MW of glutens amongst wheat varieties of varying baking performance and also that it is not the size of the soluble glutenin polymers, but the secondary structural and rheological properties of the insoluble polymer fraction that are mainly responsible for variations in baking performance. Extensional strain hardening has been shown to be a sensitive indicator of entanglements and long-chain branching in HMW polymers, and is well related to baking performance of bread doughs. The Considere failure criterion for instability in extension of polymers defines a region below which bubble walls become unstable, and predicts that when strain hardening falls below a value of around 1, bubble walls are no longer stable and coalesce rapidly, resulting in loss of gas retention and lower volume and texture. Strain hardening in doughs has been shown to reach this value at increasingly higher temperatures for better breadmaking varieties and is directly related to bubble stability and baking performance. (C) 2003 Elsevier Ltd. All rights reserved.

The physics of baking: rheological and polymer molecular structure-function relationships in breadmaking

Molecular size and structure of the gluten polymers that make up the major structural components of wheat are related to their rheological properties via modem polymer rheology concepts. Interactions between polymer chain entanglements and branching are seen to be the key mechanisms determining the rheology of HMW polymers. Recent work confirms the observation that dynamic shear plateau modulus is essentially independent of variations in MW amongst wheat varieties of varying baking performance and is not related to variations in baking performance, and that it is not the size of the soluble glutenin polymers, but the structural and rheological properties of the insoluble polymer fraction that are mainly responsible for variations in baking performance. The rheological properties of gas cell walls in bread doughs are considered to be important in relation to their stability and gas retention during proof and baking, in particular their extensional strain hardening properties. Large deformation rheological properties of gas cell walls were measured using biaxial extension for a number of doughs of varying breadmaking quality at constant strain rate and elevated temperatures in the range 25-60 degrees C. Strain hardening and failure strain of cell walls were both seen to decrease with temperature, with cell walls in good breadmaking doughs remaining stable and retaining their strain hardening properties to higher temperatures (60 degrees C), whilst the cell walls of poor breadmaking doughs became unstable at lower temperatures (45-50 degrees C) and had lower strain hardening. Strain hardening measured at 50 degrees C gave good correlations with baking volume, with the best correlations achieved between those rheological measurements and baking tests which used similar mixing conditions. As predicted by the Considere failure criterion, a strain hardening value of I defines a region below which gas cell walls become unstable, and discriminates well between the baking quality of a range of commercial flour blends of varying quality. This indicates that the stability of gas cell walls during baking is strongly related to their strain hardening properties, and that extensional rheological measurements can be used as predictors of baking quality. (C) 2004 Elsevier B.V. All rights reserved.

The physics of baking: rheological and polymer molecular structure-function relationships in breadmaking

Molecular size and structure of the gluten polymers that make up the major structural components of wheat are related to their rheological properties via modern polymer rheology concepts. Interactions between polymer chain entanglements and branching are seen to be the key mechanisms determining the rheology of HMW polymers. Recent work confirms the observation that dynamic shear plateau modulus is essentially independent of variations in MW amongst wheat varieties of varying baking performance and is not related to variations in baking performance, and that it is not the size of the soluble glutenin polymers, but the structural and rheological properties of the insoluble polymer fraction that are mainly responsible for variations in baking performance. The rheological properties of gas cell walls in bread doughs are considered to be important in relation to their stability and gas retention during proof and baking, in particular their extensional strain hardening properties. Large deformation rheological properties of gas cell walls were measured using biaxial extension for a number of doughs of varying breadmaking quality at constant strain rate and elevated temperatures in the range 25oC to 60oC. Strain hardening and failure strain of cell walls were both seen to decrease with temperature, with cell walls in good breadmaking doughs remaining stable and retaining their strain hardening properties to higher temperatures (60oC), whilst the cell walls of poor breadmaking doughs became unstable at lower temperatures (45oC to 50oC) and had lower strain hardening. Strain hardening measured at 50oC gave good correlations with baking volume, with the best correlations achieved between those rheological measurements and baking tests which used similar mixing conditions. As predicted by the Considere failure criterion, a strain hardening value of 1 defines a region below which gas cell walls become unstable, and discriminates well between the baking quality of a range of commercial flour blends of varying quality. This indicates that the stability of gas cell walls during baking is strongly related to their strain hardening properties, and that extensional rheological measurements can be used as predictors of baking quality.

Estudo comparativo de pás para aerogeradores de grande porte fabricadas em materiais compósitos reforçadas com fibra de carbono ou fibra de vidro

The research and development of wind turbine blades are essential to keep pace with worldwide growth in the renewable energy sector. Although currently blades are typically produced using glass fiber reinforced composite materials, the tendency for larger size blades, particularly for offshore applications, has increased the interest on carbon fiber reinforced composites because of the potential for increased stiffness and weight reduction. In this study a model of blade designed for large generators (5 MW) was studied on a small scale. A numerical simulation was performed to determine the aerodynamic loading using a Computational Fluid Dynamics (CFD) software. Two blades were then designed and manufactured using epoxy matrix composites: one reinforced with glass fibers and the other with carbon fibers. For the structural calculations, maximum stress failure criterion was adopted. The blades were manufactured by Vacuum Assisted Resin Transfer Molding (VARTM), typical for this type of component. A weight comparison of the two blades was performed and the weight of the carbon fiber blade was approximately 45% of the weight of the fiberglass reinforced blade. Static bending tests were carried out on the blades for various percentages of the design load and deflections measurements were compared with the values obtained from finite element simulations. A good agreement was observed between the measured and calculated deflections. In summary, the results of this study confirm that the low density combined with high mechanical properties of carbon fibers are particularly attractive for the production of large size wind turbine blades

Analysis of failed commercially pure titanium dental implants: A scanning electron microscopy and energy-dispersive spectrometer X-ray study

Background: the failure of osseointegration in oral rehabilitation has gained importance in current literature and in clinical practice. The integration of titanium dental implants in alveolar bone has been partly ascribed to the biocompatibility of the implant surface oxide layer. The aim of this investigation was to analyze the surface topography and composition of failed titanium dental implants in order to determine possible causes of failure.Methods: Twenty-one commercially pure titanium (cpTi) implants were retrieved from 16 patients (mean age of 50.33 +/- 11.81 years). Fourteen implants were retrieved before loading (early failures), six after loading (late failures), and one because of mandibular canal damage. The failure criterion was lack of osseointegration characterized as dental implant mobility. Two unused implants were used as a control group. All implant surfaces were examined by scanning electron microscopy (SEM) and energy-dispersive spectrometer x-ray (EDS) to element analysis. Evaluations were performed on several locations of the same implant.Results: SEM showed that the surface of all retrieved implants consisted of different degrees of organic residues, appearing mainly as dark stains. The surface topography presented as grooves and ridges along the machined surface similar to control group. Overall, foreign elements such as carbon, oxygen, sodium, calcium, silicon, and aluminum were detected in failed implants. The implants from control group presented no macroscopic contamination and clear signs of titanium.Conclusion: These preliminary results do not suggest any material-related cause for implant failures, although different element composition was assessed between failed implants and control implants.

An implicit BEM formulation to model strong discontinuities in solids

A boundary element method (BEM) formulation to predict the behavior of solids exhibiting displacement (strong) discontinuity is presented. In this formulation, the effects of the displacement jump of a discontinuity interface embedded in an internal cell are reproduced by an equivalent strain field over the cell. To compute the stresses, this equivalent strain field is assumed as the inelastic part of the total strain. As a consequence, the non-linear BEM integral equations that result from the proposed approach are similar to those of the implicit BEM based on initial strains. Since discontinuity interfaces can be introduced inside the cell independently on the cell boundaries, the proposed BEM formulation, combined with a tracking scheme to trace the discontinuity path during the analysis, allows for arbitrary discontinuity propagation using a fixed mesh. A simple technique to track the crack path is outlined. This technique is based on the construction of a polygonal line formed by segments inside the cells, in which the assumed failure criterion is reached. Two experimental concrete fracture tests were analyzed to assess the performance of the proposed formulation.

A finite element approach for predicting the ultimate rotation capacity of RC beams

This paper presents a numerical approach to model the complex failure mechanisms that define the ultimate rotational capacity of reinforced concrete beams. The behavior in tension and compression is described by a constitutive damage model derived from a combination of two specific damage models [1]. The nonlinear behavior of the compressed region is treated by the compressive damage model based on the Drucker-Prager criterion written in terms of the effective stresses. The tensile damage model employs a failure criterion based on the strain energy associated with the positive part the effective stress tensor. This model is used to describe the behavior of very thin bands of strain localization, which are embedded in finite elements to represent multiple cracks that occur in the tensioned region [2]. The softening law establishes dissipation energy compatible with the fracture energy of the concrete. The reinforcing steel bars are modeled by truss elements with elastic-perfect plastic behavior. It is shown that the resulting approach is able to predict the different stages of the collapse mechanism of beams with distinct sizes and reinforcement ratios. The tensile damage model and the finite element embedded crack approach are able to describe the stiffness reduction due to concrete cracking in the tensile zone. The truss elements are able to reproduce the effects of steel yielding and, finally, the compressive damage model is able to describe the non-linear behavior of the compressive zone until the complete collapse of the beam due to crushing of concrete. The proposed approach is able to predict well the plastic rotation capacity of tested beams [3], including size-scale effects.

Um estudo da estabilidade mecânica de poços a partir de perfis geofísicos em formações pouco consolidadas

A análise da estabilidade mecânica de um poço pode ser feita a partir do cálculo de parâmetros elásticos da formação utilizando a densidade do meio e as velocidades de propagação das ondas compressional e cisalhante na formação rochosa, os quais podem ser obtidos de perfis geofísicos do poço. Em formações sedimentares pouco consolidadas as ferramentas de perfilagem sônica convencionais (monopolares) não conseguem registrar com acuidade a velocidade da onda cisalhante pois a primeira chegada dessa onda é camuflada pela chegada de outras ondas que podem ser mais rápidas que a onda cisalhante num poço perfurado neste tipo de formação. Medidas das velocidades sônicas são feitas em laboratório em amostras da formação, sob condições semelhantes às condições in situ, servindo como ajuste das velocidades registradas no poço pela ferramenta de perfilagem sônica. Para a análise de estabilidade da formação, perfis auxiliares são necessários como o perfil de porosidade, saturação de fluidos e perfis de composição mineralógica da formação rochosa. Exige-se ainda dados de testes de avaliação da formação e de condições do reservatório, mas que são comuns em poços de petróleo, como o teste de formação e os testes de pressurização do poço, tais como o teste de micro-fraturamento hidráulico ou o teste de absorção. A avaliação das tensões principais efetivas que atuam distante do poço, e que não são afetadas pela sua presença, é feita através da associação de um modelo de deformação elástica apropriado e o resultado do teste de pressurização disponível para o poço em estudo. Utilizando resultados clássicos da teoria da elasticidade geral pode-se calcular o campo de tensões modificado na vizinhança da parede do poço devido ao efeito da própria presença do poço ali perfurado e da diferença de pressão existente entre o interior do poço e a formação rochosa. A determinação das propriedades mecânicas da formação a partir das velocidades sônicas e a avaliação do campo de tensões assumindo um modelo elástico de deformação, supõem o meio rochoso no qual as ondas se propagam como elástico, homogêneo e isotrópico. Esta suposição representa a principal aproximação assumida pela metodologia descrita neste trabalho. De posse das propriedades mecânicas da formação e do campo de tensões que age na vizinhança do poço resta definir o critério segundo o qual a rocha sofre instabilidade mecânica quando submetida aquele campo de tensões. Isto permite determinar se, nas condições avaliadas do poço e da formação, haverá quebra da parede do poço por excesso de tensão e, se houver, qual a sua extensão. Assim o problema é como analisar o comportamento mecânico de um poço em uma formação pouco consolidada a partir de perfis geofísicos os quais podem ter problemas no registro das propriedades físicas do meio em formações deste tipo. A metodologia proposta é aplicada a dois intervalos de profundidade pertencentes a dois poços onde arenitos e folhelhos se intercalam e nos quais todos os dados necessários estão disponíveis. Os resultados obtidos mostram que, exceto quando outros mecanismos de quebra da parede do poço agem na formação, a metodologia proposta consegue com sucesso detectar zonas de ocorrência de instabilidade mecânica do poço provocadas por um campo de tensões que excede a resistência mecânica da formação.

Predição de vida em fadiga de compósitos estruturais baseado em análises dinâmico-mecânicas

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Constitutive models for biodegradable thermoplastic ropes for ligament repair

The concept behind a biodegradable ligament device is to temporarily replace the biomechanical functions of the ruptured ligament, while it progressively regenerates its capacities. However, there is a lack of methods to predict the mechanical behaviour evolution of the biodegradable devices during degradation, which is an important aspect of the project. In this work, a hyper elastic constitutive model will be used to predict the mechanical behaviour of a biodegradable rope made of aliphatic polyesters. A numerical approach using ABAQUS is presented, where the material parameters of the model proposal are automatically updated in correspondence to the degradation time, by means of a script in PYTHON. In this method we also use a User Material subroutine (UMAT) to apply a failure criterion base on the strength that decreases according to a first order differential equation. The parameterization of the material model proposal for different degradation times were achieved by fitting the theoretical curves with the experimental data of tensile tests on fibres. To model all the rope behaviour we had considered one step of homogenisation considering the fibres architectures in an elementary volume. (C) 2012 Elsevier Ltd. All rights reserved.

Crustal heterogeneities and complexities of fault processes

In this work we study the relation between crustal heterogeneities and complexities in fault processes. The first kind of heterogeneity considered involves the concept of asperity. The presence of an asperity in the hypocentral region of the M = 6.5 earthquake of June 17-th, 2000 in the South Iceland Seismic Zone was invoked to explain the change of seismicity pattern before and after the mainshock: in particular, the spatial distribution of foreshock epicentres trends NW while the strike of the main fault is N 7◦ E and aftershocks trend accordingly; the foreshock depths were typically deeper than average aftershock depths. A model is devised which simulates the presence of an asperity in terms of a spherical inclusion, within a softer elastic medium in a transform domain with a deviatoric stress field imposed at remote distances (compressive NE − SW, tensile NW − SE). An isotropic compressive stress component is induced outside the asperity, in the direction of the compressive stress axis, and a tensile component in the direction of the tensile axis; as a consequence, fluid flow is inhibited in the compressive quadrants while it is favoured in tensile quadrants. Within the asperity the isotropic stress vanishes but the deviatoric stress increases substantially, without any significant change in the principal stress directions. Hydrofracture processes in the tensile quadrants and viscoelastic relaxation at depth may contribute to lower the effective rigidity of the medium surrounding the asperity. According to the present model, foreshocks may be interpreted as induced, close to the brittle-ductile transition, by high pressure fluids migrating upwards within the tensile quadrants; this process increases the deviatoric stress within the asperity which eventually fails, becoming the hypocenter of the mainshock, on the optimally oriented fault plane. In the second part of our work we study the complexities induced in fault processes by the layered structure of the crust. In the first model proposed we study the case in which fault bending takes place in a shallow layer. The problem can be addressed in terms of a deep vertical planar crack, interacting with a shallower inclined planar crack. An asymptotic study of the singular behaviour of the dislocation density at the interface reveals that the density distribution has an algebraic singularity at the interface of degree ω between -1 and 0, depending on the dip angle of the upper crack section and on the rigidity contrast between the two media. From the welded boundary condition at the interface between medium 1 and 2, a stress drop discontinuity condition is obtained which can be fulfilled if the stress drop in the upper medium is lower than required for a planar trough-going surface: as a corollary, a vertically dipping strike-slip fault at depth may cross the interface with a sedimentary layer, provided that the shallower section is suitably inclined (fault "refraction"); this results has important implications for our understanding of the complexity of the fault system in the SISZ; in particular, we may understand the observed offset of secondary surface fractures with respect to the strike direction of the seismic fault. The results of this model also suggest that further fractures can develop in the opposite quadrant and so a second model describing fault branching in the upper layer is proposed. As the previous model, this model can be applied only when the stress drop in the shallow layer is lower than the value prescribed for a vertical planar crack surface. Alternative solutions must be considered if the stress drop in the upper layer is higher than in the other layer, which may be the case when anelastic processes relax deviatoric stress in layer 2. In such a case one through-going crack cannot fulfil the welded boundary conditions and unwelding of the interface may take place. We have solved this problem within the theory of fracture mechanics, employing the boundary element method. The fault terminates against the interface in a T-shaped configuration, whose segments interact among each other: the lateral extent of the unwelded surface can be computed in terms of the main fault parameters and the stress field resulting in the shallower layer can be modelled. A wide stripe of high and nearly uniform shear stress develops above the unwelded surface, whose width is controlled by the lateral extension of unwelding. Secondary shear fractures may then open within this stripe, according to the Coulomb failure criterion, and the depth of open fractures opening in mixed mode may be computed and compared with the well studied fault complexities observed in the field. In absence of the T-shaped decollement structure, stress concentration above the seismic fault would be difficult to reconcile with observations, being much higher and narrower.

Sviluppo di una metodologia per la modellazione geologica-geomeccanica ed il monitoraggio di reservoir adibiti a siti di stoccaggio di gas naturale

La valutazione dei rischi associati all’operatività dei sistemi di stoccaggio, quali la sismicità indotta e la subsidenza, è requisito basilare per una loro corretta gestione e progettazione, e passa attraverso la definizione dell’influenza sullo stato tensionale delle variazioni di pressione di poro nel sottosuolo. Principale scopo di questo progetto è lo sviluppo di una metodologia in grado di quantificare le deformazioni dei reservoir in funzione della pressione di poro, di tarare i modelli utilizzati con casi studio che presentino dati di monitoraggio reali, tali da consentire un confronto con le previsioni di modello. In questa tesi, la teoria delle inomogeneità è stata utilizzata, tramite un approccio semianalitico, per definire le variazioni dei campi elastici derivanti dalle operazioni di prelievo e immissione di fluidi in serbatoi geologici. Estensione, forma e magnitudo delle variazioni di stress indotte sono state valutate tramite il concetto di variazione dello sforzo critico secondo il criterio di rottura di Coulomb, tramite un’analisi numerica agli elementi finiti. La metodologia sviluppata è stata applicata e tarata su due reservoir sfruttati e riconvertiti a sistemi di stoccaggio che presentano dataset, geologia, petrofisica, e condizioni operative differenti. Sono state calcolate le variazioni dei campi elastici e la subsidenza; è stata mappata la variazione di sforzo critico di Coulomb per entrambi i casi. I risultati ottenuti mostrano buon accordo con le osservazioni dei monitoraggi, suggerendo la bontà della metodologia e indicando la scarsa probabilità di sismicità indotta. Questo progetto ha consentito la creazione di una piattaforma metodologica di rapido ed efficace utilizzo, per stimare l’influenza dei sistemi di stoccaggio di gas sullo stato tensionale della crosta terrestre; in fase di stoccaggio, permette di monitorare le deformazioni e gli sforzi indotti; in fase di progettazione, consente di valutare le strategie operative per monitorare e mitigare i rischi geologici associati a questi sistemi.

Slope stability analysis of the Pacaya Volcano, Guatemala, using limit equilibrium and finite element method

The Pacaya volcanic complex is part of the Central American volcanic arc, which is associated with the subduction of the Cocos tectonic plate under the Caribbean plate. Located 30 km south of Guatemala City, Pacaya is situated on the southern rim of the Amatitlan Caldera. It is the largest post-caldera volcano, and has been one of Central America’s most active volcanoes over the last 500 years. Between 400 and 2000 years B.P, the Pacaya volcano had experienced a huge collapse, which resulted in the formation of horseshoe-shaped scarp that is still visible. In the recent years, several smaller collapses have been associated with the activity of the volcano (in 1961 and 2010) affecting its northwestern flanks, which are likely to be induced by the local and regional stress changes. The similar orientation of dry and volcanic fissures and the distribution of new vents would likely explain the reactivation of the pre-existing stress configuration responsible for the old-collapse. This paper presents the first stability analysis of the Pacaya volcanic flank. The inputs for the geological and geotechnical models were defined based on the stratigraphical, lithological, structural data, and material properties obtained from field survey and lab tests. According to the mechanical characteristics, three lithotechnical units were defined: Lava, Lava-Breccia and Breccia-Lava. The Hoek and Brown’s failure criterion was applied for each lithotechnical unit and the rock mass friction angle, apparent cohesion, and strength and deformation characteristics were computed in a specified stress range. Further, the stability of the volcano was evaluated by two-dimensional analysis performed by Limit Equilibrium (LEM, ROCSCIENCE) and Finite Element Method (FEM, PHASE 2 7.0). The stability analysis mainly focused on the modern Pacaya volcano built inside the collapse amphitheatre of “Old Pacaya”. The volcanic instability was assessed based on the variability of safety factor using deterministic, sensitivity, and probabilistic analysis considering the gravitational instability and the effects of external forces such as magma pressure and seismicity as potential triggering mechanisms of lateral collapse. The preliminary results from the analysis provide two insights: first, the least stable sector is on the south-western flank of the volcano; second, the lowest safety factor value suggests that the edifice is stable under gravity alone, and the external triggering mechanism can represent a likely destabilizing factor.