[pt] Este trabalho visa simular situações que levam uma rocha a ruptura considerando o fluxo de fluidos, em uma tentativa de correlacionar os mecanismos envolvidos na ruptura ao processo conhecido como produção de areia em rochas reservatório, ou seja, a produção de sólidos durante a extração de hidrocarbonetos. Há algum tempo este assunto tem sido tratado do ponto de vista da mecânica das rochas dada a similaridade do seu comportamento com outros fenômenos estudados nesta área. Como se trata de um processo de extração de fluidos, a sua influência deve ser considerada no comportamento mecânico da rocha. Assim, simulações numéricas, representativas de situações favoráveis a este fenômeno, foram analisadas, utilizando-se, para tal, o programa comercial de elementos finitos ABAQUS, que permite a simulação do processo de acoplamento fluido mecânico. Dada a flexibilidade demonstrada pelo pacote numérico, possibilitou-se a implementação de um modelo elastoplástico e a verificação de sua eficiência na representação do comportamento da rocha reservatório neste processo. A implementação numérica baseou-se no conhecido modelo de Lade-Kim. Este modelo é apropriado à descrição do comportamento de materiais geomecânicos granulares. Trabalhos experimentais em paralelo procuram determinar os parâmetros do modelo de Lade-Kim. Estes parâmetros são usados na análise numérica. / [en] This work aims to simulate situations that take a rock the rupture considering the fluid flow, in an attempt to correlate the involved mechanisms in the rupture to the know process as sand production in rocks reservoir, that is, solid production during hydrocarbon extraction. Numerical simulations using the commercial finite element program ABAQUS, reproduce fluid-mechanical coupling, and allowing the implementation of an elastoplastic rock material behavior. More precisely, the author implemented a numerical routine based on the well-known Lade-Kim elastoplastic model, which can suitably describe the behavior of geomaterials subjected to external stress-strain fields.

The examiner of cloud computing systems in the last few years observes that there is a trend of the emergence of new Big Data frameworks every single year. Since Hadoop was developed in 2007, new frameworks followed it such as Spark, Storm, Heron, Apex, Flink, Samza, Kafka ... etc. Each framework is developed in a certain way to target and achieve certain objectives better than other frameworks do. However, there are few common functionalities and aspects that are shared between these frameworks. One vital aspect all these frameworks strive to achieve is better reliability and faster recovery time in case of failures. Despite all the advances in making datacenters dependable, failures actually still happen. This is particularly onerous for long-running “big data” applications, where partial failures can lead to significant losses and lengthy recomputations. This is also crucial for streaming systems where events are processed and monitored online in real time, and any delay in data delivery will cause a major inconvenience to the users.<div>Another observation is that some reliability implementations are redundant between different frameworks. Big data processing frameworks like Hadoop MapReduce include fault tolerance mechanisms, but these are commonly targeted at specific system/failure models, and are often redundant between frameworks. Encapsulating these implementations into one layer and making it shared between different applications will benefit more than one frame-work without the burden of re-implementing the same reliability approach in each single framework.<br></div><div>These observations motivated us to solve the problem by presenting two systems: Guardian and Warden. Guardian is tailored towards batch processing big data systems while Warden is targeted towards stream processing systems. Both systems are robust, RMS based, generic, multi-framework, flexible, customizable, low overhead systems that allow their users to run their applications with individually configurable fault tolerance granularity and degree, with only minor changes to their implementation.<br></div><div>Most reliability approaches carry out one rigid fault tolerance technique targeted towards one system at a time. It is more challenging to provide a reliability approach that is pluggable in multiple Big Data frameworks at a time and can achieve low overheads comparable with single targeted framework approaches, yet is flexible and customizable by its users to make it tailored towards their objectives. The genericity is attained by providing an interface that can be used in different applications from different frameworks in any part of the application code. The low overhead is achieved by providing faster application finish times with and without failures. The customizability is fulfilled by providing the users the options to choose between two fault tolerance guarantees (Crash Failures / Byzantine Failures) and, in case of streaming systems; it is combined with two delivery semantics (Exactly Once / At Most Once).<br></div><div>In other words, this thesis proposes the paradigm of dependable resources: big data processing frameworks are typically built on top of resource management systems (RMSs),and proposing fault tolerance support at the level of such an RMS yields generic fault tolerance mechanisms, which can be provided with low overhead by leveraging constraints on resources.<br></div><div>To the best of our knowledge, such approach was never tried on multiple big data batch processing and streaming frameworks before.<br></div><div>We demonstrate the benefits of Guardian by evaluating some batch processing frame-works such as Hadoop, Tez, Spark and Pig on a prototype of Guardian running on Amazon-EC2, improving completion time by around 68% in the presence of failures, while maintaining around 6% overhead. We’ve also built a prototype of Warden on the Flink and Samza (with Kafka) streaming frameworks. Our evaluations on Warden highlight the effectiveness of our approach in the presence of failures and without failures compared to other fault tolerance techniques (such as checkpointing)<br></div>

<p dir="ltr">The process of metal additive manufacturing is characterized by the layer-by-layer construction of components, where each individual layer may be subjected to distinct thermal variations, resulting in differences in cooling rates and thermal gradients. These variations can impact the microstructure and, subsequently, mechanical properties of the final product, especially as the height of the build increases. In the present investigation, an evaluation was undertaken to ascertain the impact of build height on the structural integrity of Ti-6Al-4V samples produced using the laser powder bed fusion (LPBF) technique. The study encompassed a comprehensive examination of microstructural features, the microhardness measurement, as well as an evaluation of defect characteristics including size, location, and distribution, with respect to the build height. Tensile and fatigue tests were conducted to elucidate the potential dependence of fatigue and tensile failures on the build height. Two groups of specimens were fabricated: the first, underwent continuous fabrication, while the second involved a pause at the half height, with the process resuming after a 24-hour interval. The results of this investigation unveiled a discernible influence of the height of the build on the structural integrity of components under cyclic loading. Most fatigue specimens were observed to exhibit failure in the upper portion of the gage section with respect to the build direction. Analyses of microstructure revealed a consistent grain morphology in alignment with the build direction, and a uniform distribution of hardness throughout the build height was noted. However, for the specimens in the first group, more process-induced defects were detected within the top half of the gage section in comparison to the bottom half, while there was no noticeable difference in the distribution of defects in the second group. The results suggest that in LPBF process, as the build height is increased, there is a higher likelihood of process-induced defect formation, ultimately resulting in a reduction in structural integrity at greater build heights.</p>

京都大学 / 新制・課程博士 / 博士(工学) / 甲第24847号 / 工博第5164号 / 新制||工||1986(附属図書館) / 京都大学大学院工学研究科社会基盤工学専攻 / (主査)教授 杉浦, 邦征, 教授 高橋, 良和, 教授 山本, 貴士, 教授 北根, 安雄 / 学位規則第4条第1項該当 / Doctor of Philosophy (Engineering) / Kyoto University / DGAM

Knots are usually regarded as defects when grading lumber. In order to evaluate a member under out-of-plane loading, shear strength is one of the major mechanical properties, specifically, rolling shear (RS) strength is one of the critical mechanical properties of Cross-Laminated Timber (CLT), which determines the flexural strength of CLT under short-span bending loads. Lower grade lumber with a higher percentage of knots is recommended to be utilized for the cross-layer laminations which are mainly responsible for resisting shear stresses. Firstly, shear tests were performed in order to evaluate the effect of knots on longitudinal shear strength using shear blocks. After that, the effect of knots on the RS strength of 3-ply southern yellow pine CLT were investigated by experimental tests and an analytical model. Center-point bending tests with a span-to-depth ratio of 6 and two-plate shear tests with a loading angle of 14° were conducted on six CLT configurations composed of different types of cross layer laminations: clear flatsawn lumber with/without pith, lumber with sound knots with/without pith, and lumber with decayed knots with/without pith. The shear analogy method was implemented to evaluate the RS strength values from the bending test results, which were also compared against the results from the two-plate shear tests. It was found that: (1) The shear blocks containing sound knots had higher shear strength than matched clear shear blocks, the shear blocks containing unsound knots had lower shear strength than the matched clear shear blocks. (2) CLT specimens with cross-layer laminations with either sound knots or decayed knots had higher RS strength. (3) In general, the shear analogy method underestimated the RS strength of CLT specimens containing knots and pith.

<p> The use of AA 5754 Al-Mg alloy for automotive applications is limited by its rapid shear failure process, due to shear banding. This failure mechanism is further complicated by the presence of inhomogeneous plastic deformation, so-called Portevin-Le Chatelier (PLC) effect, during deformation. Therefore, the purpose of this study was primarily to investigate the impact of Portevin-Le Chatelier (PLC) banding towards shear banding in this commercial alloy. The second objective was to study the PLC banding as a function of prior deformation under positive strain rate sensitivity condition.</p> <p> The experimental work involved pre-straining experiments coupled with a non-contact strain measurement technique. Pre-straining experiments were carried out by deforming the sample at 223 K, at which the PLC effect is significantly suppressed, up to a prescribed amount of true strain prior to room temperature testing. A non-contact strain measurement technique, based on digital image correlation (DIC), was utilized in order to observe PLC band behavior during tensile tests at room temperature and subsequently to measure the amount of plastic strain carried within the band.</p> <p> The results showed the appearance of random nucleation deformation bands, associated with type B PLC banding, with short distance propagation during constant strain rate tensile test at room temperature. A change in the nature of PLC banding, marked by distinct band propagation, was observed once a critical amount of pre-strain is given. However, there is no evidence of a relationship between two existing phenomena, PLC banding and shear banding, in this alloy.</p> / Thesis / Master of Applied Science (MASc)

Amostras de um aço inoxidável martensítico AISI 410 temperado e revenido foram nitretadas a plasma em baixa temperatura usando o tratamento de nitretação plasma DC e a nitretação a plasma com tela ativa. Ambos os tratamentos foram realizados a 400 °C, utilizando mistura gasosa de 75 % de nitrogênio e 25 % de hidrogênio durante 20 horas e 400 Pa de pressão. As amostras de aço AISI 410 temperado e revenido foram caracterizadas antes e depois dos tratamentos termoquímicos, usando as técnicas de microscopia óptica, microscopia eletrônica de varredura, medidas de microdureza, difração de raios X e medidas de teor de nitrogênio em função da distância à superfície por espectrometria WDSX de raios X. A resistência à erosão por cavitação do aço AISI 410 nitretado DC e com tela ativa foi avaliada segundo a norma ASTM G32 (1998). Os ensaios de erosão, de erosão - corrosão e de esclerometria linear instrumentada segundo norma ASTM C1624 (2005) somente foram realizados no aço AISI 410 nitretado com tela ativa. Ensaios de nanoindentação instrumentada forma utilizados para medir a dureza (H) e o módulo de elasticidade reduzido (E*) e calcular as relações H/E* e H3/E*2 e a recuperação elástica (We), utilizando o método proposto por Oliver e Pharr. Ambos os tratamentos produziram camadas nitretadas de espessura homogênea constituídas por martensita expandida supersaturada em nitrogênio e nitretos de ferro com durezas superiores a 1200 HV, porém, a nitretação DC produziu maior quantidade de nitretos de ferro do que o tratamento de tela ativa. Os resultados de erosão por cavitação do aço nitretado DC mostraram que a precipitação de nitretos de ferro é prejudicial para a resistência à cavitação já que reduziu drasticamente o período de incubação e aumentou a taxa de perda de massa nos estágios iniciais do ensaio; entretanto, depois da remoção desses nitretos de ferro, a camada nitretada formada somente por martensita expandida resistiu bem ao dano por cavitação. Já no caso do aço nitretado com tela ativa, a resistência à erosão por cavitação aumentou 27 vezes quando comparada com o aço AISI 410 sem nitretar, fato atribuído à pequena fração volumétrica e ao menor tamanho dos nitretos de ferro presente na camada nitretada, às maiores relações H/E* e H3/E*2 e à alta recuperação elástica da martensita expandida. A remoção de massa ocorreu, principalmente, pela formação de crateras e de destacamento de material da superfície dos grãos por fratura frágil sem evidente deformação plástica. As perdas de massa acumulada mostradas pelo aço nitretado foram menores do que aquelas do aço AISI 410 nos ensaios de erosão e de erosão corrosão. O aço nitretado apresentou uma diminuição nas taxas de desgaste em ambos os ensaios de aproximadamente 50 % quando comparadas com o aço AISI 410. O mecanismo de remoção de material foi predominantemente dúctil, mesmo com o grande aumento na dureza. Os resultados de esclerometria linear instrumentada mostraram que a formação de martensita expandida possibilitou uma diminuição considerável do coeficiente de atrito em relação ao observado no caso do aço AISI 410 sem nitretar. O valor de carga crítica de falha foi de 14 N. O mecanismo de falha operante no aço nitretado foi trincamento por tensão. / Specimens of a quenched and tempered AISI 410 martensitic stainless were low temperature plasma nitrided using DC pulsed plasma treatment and the pulsed plasma active screen technic. Both treatments were carried out at 400 °C in a mixture of 75 % of nitrogen and 25 % of hydrogen during 20 hours and 400 Pa of pressure. Nitrided and non-nitrided AISI 410 specimens were characterized by optical and scanning electron microscopy, micro and nanohardness measurements, X ray diffraction and determination of the nitrogen content as a function of the depth using wavelength dispersive spectroscopy WDSX. Cavitation erosion tests were carried out according to ASTM G32 (1998) standard for both DC nitrided steel and active screen nitrided steel, whereas, the erosion, erosion - corrosion tests and scratch tests according to ASTM C1624 (2005) were conducted only for active screen nitrided steel. Nanoindentation tests were carried out in order to assess the hardness (H), the reduced elastic modulus (E*) the H/E* and H3/E*2 ratios and the elastic recovery (We) of the active screen nitrided steel according to the procedure proposed by Oliver and Pharr. Both nitrided treatments produced thick nitrided cases composed of nitrogen supersaturaded expanded martensite and iron nitrides, however, the DC treatment promoted the precipitation of large quantities of iron nitrides in comparison to the active screen technic. The cavitation erosion results of the DC nitrided steel showed that iron nitrides precipitation is harmful for the cavitation resistance as it drastically reduced the incubation period, despites this, after the removal of those iron nitrides, the nitrided case composed solely of expanded martensite resisted the cavitation damage. On the other hand, the active screen technic increased 27 times the cavitation erosion resistance of the AISI 410 steel. The increase in cavitation erosion resistance was attributed to minor quantities of smaller size iron nitrides, the higher H/E* and H3/E*2 ratios and to the higher elastic response of the expanded martensite. The material removal mainly comes from the formation of craters and from debris detachment from the grain surfaces due to brittle fracture, without plastic deformation. The active screen nitrided steel showed the lower cumulative mass losses in erosion and erosion - corrosion tests. The nitrogen addition decreased around 50 % the erosion rate in both tests. The active screen nitrided steel showed a ductile behavior despite the intense increase in hardness. The scratch tests showed that expanded martensite formation led to a significant decrease of the friction coefficient. The critical load was 14 N and the failure mechanism acting in the nitrided case was tensile cracking.

Amostras de um aço inoxidável martensítico AISI 410 temperado e revenido foram nitretadas a plasma em baixa temperatura usando o tratamento de nitretação plasma DC e a nitretação a plasma com tela ativa. Ambos os tratamentos foram realizados a 400 °C, utilizando mistura gasosa de 75 % de nitrogênio e 25 % de hidrogênio durante 20 horas e 400 Pa de pressão. As amostras de aço AISI 410 temperado e revenido foram caracterizadas antes e depois dos tratamentos termoquímicos, usando as técnicas de microscopia óptica, microscopia eletrônica de varredura, medidas de microdureza, difração de raios X e medidas de teor de nitrogênio em função da distância à superfície por espectrometria WDSX de raios X. A resistência à erosão por cavitação do aço AISI 410 nitretado DC e com tela ativa foi avaliada segundo a norma ASTM G32 (1998). Os ensaios de erosão, de erosão - corrosão e de esclerometria linear instrumentada segundo norma ASTM C1624 (2005) somente foram realizados no aço AISI 410 nitretado com tela ativa. Ensaios de nanoindentação instrumentada forma utilizados para medir a dureza (H) e o módulo de elasticidade reduzido (E*) e calcular as relações H/E* e H3/E*2 e a recuperação elástica (We), utilizando o método proposto por Oliver e Pharr. Ambos os tratamentos produziram camadas nitretadas de espessura homogênea constituídas por martensita expandida supersaturada em nitrogênio e nitretos de ferro com durezas superiores a 1200 HV, porém, a nitretação DC produziu maior quantidade de nitretos de ferro do que o tratamento de tela ativa. Os resultados de erosão por cavitação do aço nitretado DC mostraram que a precipitação de nitretos de ferro é prejudicial para a resistência à cavitação já que reduziu drasticamente o período de incubação e aumentou a taxa de perda de massa nos estágios iniciais do ensaio; entretanto, depois da remoção desses nitretos de ferro, a camada nitretada formada somente por martensita expandida resistiu bem ao dano por cavitação. Já no caso do aço nitretado com tela ativa, a resistência à erosão por cavitação aumentou 27 vezes quando comparada com o aço AISI 410 sem nitretar, fato atribuído à pequena fração volumétrica e ao menor tamanho dos nitretos de ferro presente na camada nitretada, às maiores relações H/E* e H3/E*2 e à alta recuperação elástica da martensita expandida. A remoção de massa ocorreu, principalmente, pela formação de crateras e de destacamento de material da superfície dos grãos por fratura frágil sem evidente deformação plástica. As perdas de massa acumulada mostradas pelo aço nitretado foram menores do que aquelas do aço AISI 410 nos ensaios de erosão e de erosão corrosão. O aço nitretado apresentou uma diminuição nas taxas de desgaste em ambos os ensaios de aproximadamente 50 % quando comparadas com o aço AISI 410. O mecanismo de remoção de material foi predominantemente dúctil, mesmo com o grande aumento na dureza. Os resultados de esclerometria linear instrumentada mostraram que a formação de martensita expandida possibilitou uma diminuição considerável do coeficiente de atrito em relação ao observado no caso do aço AISI 410 sem nitretar. O valor de carga crítica de falha foi de 14 N. O mecanismo de falha operante no aço nitretado foi trincamento por tensão. / Specimens of a quenched and tempered AISI 410 martensitic stainless were low temperature plasma nitrided using DC pulsed plasma treatment and the pulsed plasma active screen technic. Both treatments were carried out at 400 °C in a mixture of 75 % of nitrogen and 25 % of hydrogen during 20 hours and 400 Pa of pressure. Nitrided and non-nitrided AISI 410 specimens were characterized by optical and scanning electron microscopy, micro and nanohardness measurements, X ray diffraction and determination of the nitrogen content as a function of the depth using wavelength dispersive spectroscopy WDSX. Cavitation erosion tests were carried out according to ASTM G32 (1998) standard for both DC nitrided steel and active screen nitrided steel, whereas, the erosion, erosion - corrosion tests and scratch tests according to ASTM C1624 (2005) were conducted only for active screen nitrided steel. Nanoindentation tests were carried out in order to assess the hardness (H), the reduced elastic modulus (E*) the H/E* and H3/E*2 ratios and the elastic recovery (We) of the active screen nitrided steel according to the procedure proposed by Oliver and Pharr. Both nitrided treatments produced thick nitrided cases composed of nitrogen supersaturaded expanded martensite and iron nitrides, however, the DC treatment promoted the precipitation of large quantities of iron nitrides in comparison to the active screen technic. The cavitation erosion results of the DC nitrided steel showed that iron nitrides precipitation is harmful for the cavitation resistance as it drastically reduced the incubation period, despites this, after the removal of those iron nitrides, the nitrided case composed solely of expanded martensite resisted the cavitation damage. On the other hand, the active screen technic increased 27 times the cavitation erosion resistance of the AISI 410 steel. The increase in cavitation erosion resistance was attributed to minor quantities of smaller size iron nitrides, the higher H/E* and H3/E*2 ratios and to the higher elastic response of the expanded martensite. The material removal mainly comes from the formation of craters and from debris detachment from the grain surfaces due to brittle fracture, without plastic deformation. The active screen nitrided steel showed the lower cumulative mass losses in erosion and erosion - corrosion tests. The nitrogen addition decreased around 50 % the erosion rate in both tests. The active screen nitrided steel showed a ductile behavior despite the intense increase in hardness. The scratch tests showed that expanded martensite formation led to a significant decrease of the friction coefficient. The critical load was 14 N and the failure mechanism acting in the nitrided case was tensile cracking.

[pt] Este trabalho apresenta resultados de um estudo realizado visando avaliar o mecanismo de ruptura associado ao escorregamento em solo não saturado ocorrido na Vista Chinesa (RJ) em fevereiro de 1988. São apresentados os resultados de um programa que envolveu investigações de campo e de laboratório, através do qual foram definidos os parâmetros necessários à realização de diversas análises. Foram utilizados os programas de computador SWMS-2D (análises de problemas de fluxo permanente/transiente em meios saturados / não saturados) e PC-SLOPE (análise de estabilidade de taludes utilizando de métodos de equilíbrio limite), com os quais procurou-se reproduzir, da maneira mais realista possível, os fenômenos observados na encosta em estudo. / [en] This thesys presents the results of one study aiming to evaluate failure mechanisms associated with the unsaturated soil lanslides which ocurred in Vista Chinesa (RJ) in february, 1988. There are showed the field and laboratory investigation program, which was the basis to define the parameters required for numerical simulations. In these simulations the SWMS-2D(analysis of permanent/transient flow problems in saturated/unsaturated media) and PC-SLOPE (stability analysis using the limit equilibrium methods) programs were used, in na attempt to reproduce as close as possiible the instability phenomena in the slope under investigation. / [es] Este trabajo presenta resultados de un estudio realizado para evaluar el mecanismo de ruptura asociado al deslizamiento de suelo no saturado, hecho que ocurrió en la zona de la Vista Chinesa (Rio de JaneiroJ) en Febrero de 1988. Se presentan los resultados de un programa que abarco investigaciones de campo y de laboratorio, a través del cual fueron definidos los parámetros necesarios a la realización de diversos análisis. Se utilizaron los programas de computador SWMS-2D (análisis de problemas de flujo permanente/transiente en medios saturados / no saturados) y PC-SLOPE (análisis de estabilidad de taludes utilizando métodos de equilibrio límite), con los cuales se intentó reproducir, de la manera más realista posible, los fenómenos observados en el lugar de estudio.

[pt] O método de exploração Sublevel Caving é um dos métodos de extração massiva mais usados na indústria mundial de exploração subterrânea, sendo considerado pela indústria de mineração, num futuro próximo dentre os substitutos naturais das atuais minas a céu aberto. Uma operação Sublevel Caving requer que o maciço rochoso circundante ao minério rompa continuamente e se movimente para dentro do espaço criado pela extração do minério. Análises existentes na literatura consideram apenas configurações parciais do processo Sublevel Caving sem considerar o processo evolutivo da extração, e o dano induzido ao maciço rochoso decorrente deste processo. Esta dissertação desenvolve uma modelagem numérica utilizando o método dos elementos discretos para simular o mecanismo de ruptura e a subsidência causada pelo método de exploração Sublevel Caving, analisando o referido efeito e suas consequências na evolução do mecanismo de ruptura e subsidência no Sublevel Caving. O software comercial Particle Flow Code (PFC2D) foi selecionado para esta modelagem devido à capacidade de simular, em um evento de excesso de tensão, o fraturamento do maciço rochoso e sua desintegração, desta forma, originam-se o fluxo do material e os deslocamentos em grande escala, os quais são considerados fenômenos físicos dominantes que regem a formação da subsidência e fraturamento num processo Sublevel Caving. Os resultados obtidos nesse estudo mostraram-se satisfatórios, reproduzindo adequadamente a superfície de subsidência induzida por Sublevel Caving, conseguindo-se uma simulação física realista da sua evolução desde o início do fraturamento até à subsidência final. / [en] The Sublevel Caving Method is one of the most massive extraction methods used in underground mining industry worldwide and is considered by the mining industry as one of the natural replacements of the current open cut mines in the near future. A Sublevel Caving operation requires that the rockmass surrounding the orebody continually fails and moves into the void created by ore extraction. This dissertation develops a modeling using the discrete element method to simulate the failure mechanism and subsidence caused by Sublevel Caving method. Analyses reported in the literature consider only partial configurations of process Sublevel Caving, without taking into consideration the excavation evolution process, and damage induced to the rock mass resulting from this process. This dissertation analyzes this effect and its consequences on the evolution of failure mechanism and subsidence in Sublevel Caving. Particle Flow Code (PFC2D) was selected for modeling because of its ability to simulate, if the event of excess stress, fracturing and disintegration of the rock mass and large-scale deformation and material flow, to be simulated, which are believed to be the dominant physical phenomena governing the formation of subsidence and fracturing of Sublevel Caving. The results obtained in this study were satisfactory, reproducing properly the surface subsidence induced by Sublevel Caving, allowing physically realistic simulation of its evolution since the beginning of the fracturing to final subsidence.