Dissertations / Theses on the topic 'Rock deformation Simulation methods'

Important orebody characteristics that determine viability of the mineral resource and ore reserve potential such as physical properties, mineralogical and geochemical compositions often vary substantially across an ore deposit. Geometallurgical models aim to capture the spatial relationships between mineral compositions, physical properties of rock and their interactions with mechanical and chemical processes during mining extraction and processing. This characterisation of physical and chemical properties of ores can in turn be used to inform mining and processing decisions that enable the extraction of the maximum value from the ore deposit most efficiently. During the construction of such spatial geometallurgical models, practitioners are presented with many challenges. These include modelling high-dimensional data of various types including categorical, continuous and compositional attributes and their uncertainties. Decisions on how to segregate samples data into spatially and statistically homogeneous groups to satisfy modelling assumptions such as stationarity are often a requirement. Secondary properties such as metallurgical test results are often few in number, acquired on larger scales than that of primary rock property data and non-additive in nature. In this thesis a data driven workflow that aims to address these challenges when constructing geometallurgical models of ore deposits is devised. Spatial machine learning techniques are used to derive geometallurgical categories, or classes, from multiscale, multiresolution, high dimensional rock properties. In supervised mode these methods are also used to predict geometallurgical classes at samples where rock property information is incomplete. Realisations of the layout of geometallurgical classes and the variabilities of associated rock properties are then mapped using geostatistical simulations and machine learning. The workflow is demonstrated using a case study at Orebody H; a complex stratabound Bedded Iron Ore deposit in Western Australia’s Pilbara. A detailed stochastic model of five compositions representing primary rock properties and geometallurgical responses in the form of lump and fine product iron ore quality specifications was constructed. The predicted product grade recoveries are realistic values that honour constraints of the predicted head grade compositions informed by more abundant and regularly spaced sampling than metallurgical tests. Finally, uncertainties are quantified to assess risk following a confidence interval based framework. This could be used to identify zones of high uncertainty where collection of additional data might help mitigate or minimise risks and in turn improve forecast production performances.

Characterization of rock masses and evaluation of their mechanical properties are important and challenging tasks in rock mechanics and rock engineering. Since in many cases rock quality designation (RQD) is the only rock mass classification index available, this paper outlines the key aspects on determination of RQD and evaluates the empirical methods based on RQD for determining the deformation modulus and unconfined compressive strength of rock masses. First, various methods for determining RQD are presented and the effects of different factors on determination of RQD are highlighted. Then, the empirical methods based on RQD for determining the deformation modulus and unconfined compressive strength of rock masses are briefly reviewed. Finally, the empirical methods based on RQD are used to determine the deformation modulus and unconfined compressive strength of rock masses at five different sites including 13 cases, and the results are compared with those obtained by other empirical methods based on rock mass classification indices such as rock mass rating (RMR), Q-system (Q) and geological strength index (GSI). It is shown that the empirical methods based on RQD tend to give deformation modulus values close to the lower bound (conservative) and unconfined compressive strength values in the middle of the corresponding values from different empirical methods based on RMR, Q and GSI. The empirical methods based on RQD provide a convenient way for estimating the mechanical properties of rock masses but, whenever possible, they should be used together with other empirical methods based on RMR, Q and GSI. (C) 2016 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting by Elsevier B.V.

Rock is a very heterogeneous material, containing structural weakness at all scales. These weaknesses include grain boundaries, pores, and cracks on the small scale, and joints, faults, and bedding planes on the large scale. Nonlinear rock deformation in the low-temperature, low-confinement regime is due primarily to the growth of cracks from these weaknesses and the coalescence of cracks to form macroscopic structural features. Another important aspect of rock deformation and failure is the statistical distribution of weaknesses in the initial microstructure. Borehole breakout is the process by which portions of a borehole wall fracture or spall when subjected to compressive stresses. Studies of borehole breakout in the past twenty years include experiments, field studies, and numerical modeling. With regards to the numerical modeling of borehole breakout, the rock surrounding the borehole is considered as a nonlinear continuum material in most of the previous approaches. Experiments and field studies, however, have shown that the heterogeneous and discontinuous nature of rock has a strong impact on the mechanics of borehole breakout. This dissertation describes a numerical model that has been developed to simulate the damage of rock and the corresponding non-linear stress-strain behavior, and also the progression of borehole breakout in heterogeneous and discontinuous rock by mixed mode crack growth, interaction, and coalescence. The rock is simulated as an elastic material containing a random distribution of cracks. As compressive load is applied, the initial cracks grow, interact, and coalesce to form macroscopic fractures. The numerical model was developed by making a series of modifications to the displacement discontinuity code of Crouch and Starfield (Crouch & Starfield, 1983). The most important modifications include modifying the boundary element for the calculation of stress intensity factors, adding Coulomb friction for closed portions of cracks, adding a crack generator, and adding an algorithm for crack coalescence. The numerical model is used to simulate the non-linear deformation and the progression of breakout in Westerly granite, and the results are realistic.

La rupture des versants rocheux est une source de risque majeur dans toutes les regions montagneuses. Les terrains metamorphiques anisotropes presentent tout l'eventail des mecanismes de deformation, lesquels sont analyses dans la premiere partie de ce travail a partir d'exemples europeens. L'etude pluridisciplinaire des ruines de sechilienne, presentee en seconde partie, revet un interet particulier du fait de la complexite du phenomene et par l'ampleur des moyens d'auscultation et de surveillance mis en uvre pendant une decennie (1988-1998). Nous avons menees trois approches complementaires : geologique, geomecanique, hydrogeologique. L'approche geologique a permise une description fine de la structure du versant a toutes les echelles. L'approche geomecanique a consiste a elaborer un modele structural simplifie et a tester numeriquement differentes hypotheses par la methodes des elements distincts. Ceci nous autorise a proposer un mecanisme en accord qualitatif avec les observations de terrain et les mesures de deplacement. Ce mecanisme de rupture interne, controle par les grandes familles de discontinuites preexistantes, induit des concentrations de contrainte et des efforts de traction pouvant provoquer un endommagement progressif et irreversible du massif. Enfin, l'approche hydrogeologique s'est attachee a analyser le debit, la temperature ainsi que la composition chimique et isotopique des eaux du massif, pour preciser leur origine et leur mode d'ecoulement. En couplant les trois approches, on a pu proposer un modele hydromecanique global expliquant l'influence de la pluviometrie sur le rythme de la deformation du versant. Des pistes sont avancees pour etudier la geometrie probable de la rupture. Par contre, la prevision de la date de la rupture reste hors d'atteinte. Enfin, l'observation du versant sud de la toura, a saint-christophe-en-oisans, suggere que le mecanisme de rupture interne pourrait etre plus repandu qu'on ne le croit.

Les séismes le long de grandes failles crustales représentent un danger énorme pour de nombreuses populations. Le mécanique de ces failles est influencé par des zones endommagées qui entourent le coeur de faille. La fracturation dans ces zones contrôle chaque étape du cycle sismique. En effet, cette zone contrôle la mécanique de la rupture sismique, elle est un conduit pour les fluides, réagit chimiquement sous l'effet de fluides réactifs, et facilite la déformation pendant les périodes post- et inter-sismiques. Dans cette thèse de doctorat, des expériences de laboratoire ont été réalisées pour mieux comprendre 1) la façon dont l'endommagement est généré pendant le chargement transitoire co-sismique, 2) comment l'endommagement permet de mieux contraindre le chargement co-sismique le long de grandes failles, et iii) comment les fractures peuvent se cicatriser au fil du temps et contrôler l'évolution de la perméabilité et de la résistance mécanique de la faille.L'introduction de la thèse propose une revue critique de la littérature sur la génération de dommages co-sismiques et en particulier sur la formation des roches pulvérisées. Le potentiel de ces roches comme marqueur des déformations co-sismiques est discuté. Bien que ces roches pulvérisées soient prometteuses pour ces aspects, plusieurs questions restent ouvertes.L'une de ces questions concerne les conditions de chargement transitoire nécessaires pour atteindre la pulvérisation. Le seuil de taux de deformation pour atteindre la pulvérisation peut être réduit par des endommagemments progressifs, au cours de ruptures sismiques successives. Des barres de Hopkinson ont été utilisées pour effectuer des chargements dynamique successifs d'une roche cristalline (monzonite). Les résultats montrent que le seuil pour atteindre la pulvérisation est réduit d'au moins 50% lorsque des chargements successives sont imposés. Cette thèse discute aussi pourquoi les roches pulvérisées sont presque toujours observées dans des roches cristallines et peu dans des roches sédimentaires poreuses. Pour comprendre cette observation, des expériences à haute vitesse de déformation ont été effectuées sur des grès de Rothbach. Les résultats montrent que la pulvérisation des grains eux mêmes ne se produit pas dans les grès. L'endommagement reste se produit principalement à une échelle supérieure à celle grains, et des bandes de compaction sont observées. La compétition entre l'endommagement inter- et intra-granulaire est expliquée par les paramètres microstructuraux en combinant deux modèles micromécaniques classiques. Les microstructures observées dans les grès peuvent se former dans le régime quasi-statiques et aussi dans le régime dynamique. Par conséquent, il est recommandée d'être prudent lors de l'interprétation du mécanisme de deformation dans les roches sédimentaires proches de la surface. La dernière question abordée durant la thèse est la cicatrisation post-sismique de fractures co-sismiques. Des expériences ont été réalisées pour cicatriser des fissures par précipitation de calcite. Le but est l'étude du couplage entre l'augmentation de résistance mécanique de la roche fissurée et l'évolution de la perméabilité. Les échantillons fracturées ont été soumis à des conditions de pression et températures similaires de la croûte supérieure et à une percolation d'un fluide sursaturé en calcite pendant plusieurs mois. Ce couplage non-existe dans les premières étapes de la cicatrisation. Il est révélé par l'imagerie par tomographie aux rayons X que le scellement naissant des fractures se produit dans les porosités situées en aval de barrières d'écoulement, et donc dans des régions qui ne touchent pas les principales voies d'écoulement du fluide. Le découplage entre l'augmentation de résistance de la roche et la perméabilité suggère que les zones d'endommagement peu profondes dans les failles actives peuvent rester des conduits actifs pour les fluides plusieurs années après un séisme
Earthquakes along large crustal scale faults are a huge hazard threatening large populations. The behavior of such faults is influenced by the fault damage zone that surrounds the fault core. Fracture damage in such fault damage zones influences each stage of the seismic cycle. The damage zone influences rupture mechanics, behaves as a fluid conduit to release pressurized fluids at depth or to give access to reactive fluids to alter the fault core, and facilitates strain during post- and interseismic periods. Also, it acts as an energy sink for earthquake energy. Here, laboratory experiments were performed to come to a better understanding of how this fracture damage is formed during coseismic transient loading, what this fracture damage can tell us about the earthquake rupture conditions along large faults, and how fracture damage is annihilated over time.First, coseismic damage generation, and specifically the formation of pulverized fault damage zone rock, is reviewed. The potential of these pulverized rocks as a coseismic marker for rupture mechanisms is discussed. Although these rocks are promising in that aspect, several open questions remain.One of these open questions is if the transient loading conditions needed for pulverization can be reduced by progressively damaging during many seismic events. The successive high strain rate loadings performed on quartz monzonites using a split Hopkinson pressure bar reveal that indeed the pulverization strain rate threshold is reduced by at least 50%.Another open question is why pulverized rocks are almost always observed in crystalline lithologies and not in more porous rock, even when crystalline and porous rocks are juxtaposed by a fault. To study this observation, high strain rate experiments were performed on porous Rothbach sandstone. The results show that pervasive pulverization below the grain scale, such as observed in crystalline rock, does not occur in the sandstone samples for the explored strain rate range (60-150 s-1). Damage is mainly occurs at a scale superior to that of the scale of the grains, with intragranular deformation occurring only in weaker regions where compaction bands are formed. The competition between inter- and intragranular damage during dynamic loading is explained with the geometric parameters of the rock in combination with two classic micromechanical models: the Hertzian contact model and the pore-emanated crack model. In conclusion, the observed microstructures can form in both quasi-static and dynamic loading regimes. Therefore caution is advised when interpreting the mechanism responsible for near-fault damage in sedimentary rock near the surface. Moreover, the results suggest that different responses of different lithologies to transient loading are responsible for sub-surface damage zone asymmetry.Finally, post-seismic annihilation of coseismic damage by calcite assisted fracture sealing has been studied in experiments, so that the coupling between strengthening and permeability of the fracture network could be studied. A sample-scale fracture network was introduced in quartz monzonite samples, followed exposure to upper crustal conditions and percolation of a fluid saturated with calcite for several months. A large recovery of up to 50% of the initial P-wave velocity drop has been observed after the sealing experiment. In contrast, the permeability remained more or less constant for the duration of the experiment. This lack of coupling between strengthening and permeability in the first stages of sealing is explained by X-ray computed micro tomography. Incipient sealing in the fracture spaces occurs downstream of flow barriers, thus in regions that do not affect the main fluid flow pathways. The decoupling of strength recovery and permeability suggests that shallow fault damage zones can remain fluid conduits for years after a seismic event, leading to significant transformations of the core and the damage zone of faults with time

Thin and slender structures are widely occurring both in nature and in human creations. Clever geometries of thin structures can produce strong constructions while using a minimal amount of material. Computer modeling and analysis of thin and slender structures has its own set of problems stemming from assumptions made when deriving the equations modeling their behavior from the theory of continuum mechanics. In this thesis we consider two kinds of thin elastic structures; threads and plates. Real-time simulation of threads are of interest in various types of virtual simulations such as surgery simulation for instance. In the first paper of this thesis we develop a thread model for use in interactive applications. By viewing the thread as a continuum rather than a truly one dimensional object existing in three dimensional space we derive a thread model that naturally handles both bending, torsion and inertial effects. We apply a corotational framework to simulate large deformation in real-time. On the fly adaptive resolution is used to minimize corotational artifacts. Plates are flat elastic structures only allowing deflection in the normal direction. In the second paper in this thesis we propose a family of finite elements for approximating solutions to the Kirchhoff-Love plate equation using a continuous piecewise linear deflection field. We reconstruct a discontinuous piecewise quadratic deflection field which is applied in a discontinuous Galerkin method. Given a criterion on the reconstruction operator we prove a priori estimates in energy and L2 norms. Numerical results for the method using three possible reconstructions are presented.

Orientador: Sérgio Tonini Button
Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecânica
Made available in DSpace on 2018-08-19T02:21:16Z (GMT). No. of bitstreams: 1 Alves_JoseAugustoCamargo_M.pdf: 12204746 bytes, checksum: 78b1a8db87093cdc58baaa854932096e (MD5) Previous issue date: 2011
Resumo: O processo sheet hydroform, ou hidroconformação de chapas, é realizado por meio de uma prensa composta por uma membrana de borracha, preenchido com um fluido hidráulico cuja função é atuar como uma matriz, exercendo esforços sobre uma chapa de metal (blank), que é então empurrada contra um punção rígido, fazendo-o adquirir o formato deste. Embora este processo seja amplamente utilizado para a produção de pequenos lotes de peças metálicas de formato complexo e de espessura reduzida, ele exige habilidades de quem o define, pois se por um lado pode ser visto como simples por empregar apenas um molde maciço, por outro, a ausência de um sistema macho-fêmea capaz de garantir um completo travamento do blank pode ser encarado como um problema por permitir movimentos indesejados do material, que muitas vezes resultam na formação de rugas ou outros defeitos. Baseando-se nesses conceitos, o propósito deste trabalho foi o de mapear, por meio de simulações e ensaios práticos, a conformabilidade de abas convexas da liga de alumínio AA2024-0 de quatro espessuras quando submetidas a diferentes combinações de raio de curvatura e comprimento de aba. Deste modo, foi possível definir quais combinações destes parâmetros possibilitam a obtenção de componentes conformados adequadamente, isentos de rugas e outros defeitos macroscópicos, e ainda, quais condições levam à formação de irregularidades na aba conformada acima do limite preestabelecido, exigindo o emprego de prensa-chapas especiais, também conhecidos como dams, capazes de evitar a ocorrência de tais desvios. Com base nos resultados obtidos pode-se constatar que a ocorrência de rugas está associada principalmente à altura da aba conformada e não se altera significativamente quando a espessura do blank é modificada. Além disso, foi possível notar que raios de curvaturas maiores proporcionam menores valores de deformação compressiva na região conformada, permitindo obter abas mais altas e sem rugas
Abstract: Hydroform, or sheet metal fluid forming, is performed using a fluid cell press, in which the hydraulic fluid acts on the metallic blank pushing it against the male tool, acquiring its geometry. It is widely employed to manufacture small batches of complex and low thick components. If by one point of view it can be seen as simple, involving just a single rigid block as tool, by the other hand the absence of a rigid punch in certain cases can be a limitation, since it may allow the blank to move incorrectly during the process, causing wrinkles or other macro defects. Based on this limitation, the aim of this study was to define, using computational simulations and practical tests, the shrink flange formability limit of four different thickness aluminum alloy sheets when submitted to different combinations of curvature radius and flange length. As result, it could be seen which combinations can lead the material to be formed properly and which may cause failures, requiring special blank holders, known as dams, to avoid this problems. Based on the results, it can be verified that wrinkles nucleation is mainly associated with flange height and it does not change significantly when using blanks with different thicknesses. Furthermore, it could be noted that bigger curvature radius implies in smaller compressive strain on formed region, allowing to obtain higher flanges without wrinkles in these conditions
Mestrado
Materiais e Processos de Fabricação
Mestre em Engenharia Mecânica

L’accent mis actuellement par l’industrie pétrolière sur l’augmentation de l’efficacité des réservoirs, ainsi que sur l’intérêt grandissant pour l’exploitation d’autres sources d’énergie enfouies profondément sous terre a suscité un regain d’intérêt pour la mécanique de la fracturation des roches en général et la fracturation hydraulique en particulier. La fracturation hydraulique, appelée de manière informelle “fracturation”, est un processus qui consiste généralement à injecter de l’eau, sous haute pression dans une formation rocheuse via le puits. Ce processus vise à créer de nouvelles fractures dans la roche et à augmenter la taille, l’étendue et la connectivité des fractures existantes. Des avancées récentes dans la modélisation et la simulation de fractures hydrauliques ont eu lieu, au confluent de facteurs qui incluent une activité accrue, une tendance vers une complexité accrue et une compréhension approfondie du modèle mathématique sous-jacent et de ses défis intrinsèques. Cependant, certaines des caractéristiques très importantes de ce processus ont été négligées. Parmi les caractéristiques négligées, on peut citer l’incapacité de la grande majorité des modèles existants de s’attaquer à la fois à la propagation de fractures hydrauliques dans la roche intacte, à l’inititation de nouvelles fractures ainsi qu’à la réactivation des fractures existantes. Une autre caractéristique qui a été ignorée est sa dimension intrinsèque en trois dimensions, négligée par la plupart des modèles actuallement proposés. Parmi tous les différents types de méthodes numériques développées pour évaluer le mécanisme du phénomène de fracturation, très peu sont capables de représenter la totalité des mecanismes mis en jeu. Dans la présente thèse, l’initiation et la propagation de fissures induites par les fluides dans des roches isotropes transversales sont simulées à l’aide d’un modèle hydromécanique (HM) couplé basé la méthode XFEM (eXtended Finite Element Method) et un modèle de zones cohésives. Le HM-XFEM développé dans cette thèse est une extension des modèles précédemment développés dans l’équipe hydro-géomécanique multi échelle de GeoRessources. L’accent a été porté plus particulièrement sur la prise en compte de l’anisotropie du milieu et sur son influence sur le chemin de propagation. Ce dernier est défini à partir du le concept d’angle de bifurcation introduit auparavant dans la littérature. En complément des développements réalisés dans le modèle HM-XFEM, effort a été fait pour mieux comprendre l’initiation de la fissure en utilisant la méthode des éléments discrets (DEM) à l’aide du logiciel open source YADE Open DEM. La nature différente des deux méthodes, DEM étant une méthode discontinue et XFEM, une méthode continue, révèle les potentiels des deux méthodes et permet de comparer correctement la méthode qui convient le mieux au problème à résoudre, compte tenu des objectifs de la conception
Current emphasis in petroleom industry toward increasing the reservoirs efficiency, along with the interest in exploitation of other sources of energy buried deep underground created a renewed interest in rock fracture mechanics in general and hydraulic fracturing specifically. Hydraulic fracturing, informally referred to as “fracking,” is an oil and gas well development process that typically involves injecting water, under high pressure into a bedrock formation via the well. This process is intended to create new fractures in the rock as well as increase the size, extent, and connectivity of existing fractures. However some of the very important features of this process have been overlooked. Among these neglected features one can name of inability of the vast majority of existing models to tackle at once the propagation of hydraulic fractures in fractured rocks-masses where a competing dipole mechanism exists between fracturing of the intact rock and re-activation of exiting fracture networks. Another feature that has been ignored is its intrinsically three dimensionality which is neglected by most models. Among all different types of numerical methods that have been developed in order to assess the mechanism of fracturing phenomenon very few, if any, can handle the entire complexity of such process. In the present thesis, fluid-driven crack initiation and propagation in transverse isotropic rocks is simulated using a coupled model comprising of eXtended Finite Element Method (XFEM) and cohesive zone models. The HM XFEM developed in this thesis is an extension to previous models developed introduced in multiscale hydrogeomechanics team of GeoRessources. An emphasis is put on considering the anisotropic nature of the medium and on studying its influence on the propagation path. This latter is investigated by the concept of bifurcation angle previously introduced in literature. In complementary efforts was made to have a better understanding of crack initiation in transversely isotropic media, we also used the discrete element method (DEM) in order to gain insights into the mechanisms at stake. Both methods exhibit their advantages and disadvantages in modeling fracturing phenomenon. The different nature of two methods, DEM being a discontinuous and XFEM being a continuous method, reveals potentials of both methods and renders a good comparison of which method suits the problem in hand the best, considering the the objectives of the design

Pour juger de la pertinence du concept de bande de cisaillement applique a la rupture des roches aussi bien dans le domaine fragile que ductile, deux methodes de suivi de la deformation lors d'essais homogenes ont ete mises en uvre, a savoir, la tomodensitometrie et la stereophotogrammetrie de faux relief. Appliquees a des essais triaxiaux et biaxiaux realises avec ou sans confinement, ces deux methodes permettent, chacune a leur echelle, de cerner la transition d'une deformation diffuse a un mode localise avec apparition de fissures. Avec la stereocomparaison, l'ouverture de fissures, leurs interactions et leur contribution a la deformation globale peuvent entre autres etre quantifiees. Le geomateriau cohesif qui a ete l'objet d'une telle etude est la marne de beaucaire, une roche tendre saturee finement grenue dont la pression de consolidation est de 5 mpa. La rupture de type fragile avec propagation brutale en mode de cisaillement d'une macro-fissure que l'on observe en compression simple et sous tres faible confinement (0. 05 ou 0. 1 mpa) a pu etre analysee de meme que la localisation en bandes de cisaillement contractantes obtenue sous 0. 25 mpa ou la densification du materiau observee sous plus fort confinement (1 mpa). Dans la mesure ou la notion de bande de cisaillement a ete, a partir de nos resultats experimentaux, jugee adequate pour decrire le comportement a la rupture de la marne de beaucaire dans les domaines de contraintes consideres, une analyse de bifurcation en mode localise visant a predire le seuil de localisation et l'orientation de la bande de cisaillement, a ete conduite avec une loi incrementale non lineaire particulierement souple, la loi cloe. Pour mener a bien une telle analyse sur un cas concret: le trou de forage, les parametres de cloe pour la marne de beaucaire ont d'abord ete identifies sur des chemins triaxiaux et biaxiaux elementaires

Description d'une naturelle, l'argile de pornic, suivant les concepts d'etat limite et d'etat critique. On realise une serie d'essais triaxiaux axisymetriques sur argile consolidee et on met en evidence la courbe d'etat limite du materiau, intersection de la surface d'etat limite et de la surface dite "d'etat surconsolide". La representation des chemins contraintes-deformations dans un diagramme original ou dans le diagramme normalise permet de cartacteriser leur reciprocite et de suivre leur evolution au passage de la courbe d'etat limite. L'etude de l'influence du temps ou de la vitesse de deformation demontre que les enveloppes de fluage se contractent dans l'espace des contraibntes quand le temps augmente et qu'une unique fait defaut pour caracteriser reciproquement les chemins de contraintes draines et non-draines en fonction des vitesses de sollicitations

Dreiaxiale Druckprüfungen können als Einstufenversuche, als Mehrstufenversuche oder als Versuche mit kontinuierlichen Bruchzuständen ausgeführt werden. Bei der Anwendung der Mehrstufentechnik ergeben sich insbesondere Fragestellungen hinsichtlich der richtigen Wahl des Umschaltpunktes und des optimalen Verlaufs des Spannungspfades zwischen den einzelnen Versuchsstufen. Fraglich beim Versuch mit kontinuierlichen Bruchzuständen bleibt, ob im Versuchsverlauf tatsächlich Spannungszustände erfasst werden, welche die Höchstfestigkeit des untersuchten Materials repräsentieren. Die Dissertation greift diese Fragestellungen auf, ermöglicht den Einstieg in die beschriebene Thematik und schafft die Voraussetzungen, die zur Lösung der aufgeführten Problemstellungen notwendig sind. Auf der Grundlage einer umfangreichen Datenbasis gesteinsmechanischer und petrophysikalischer Kennwerte wurde ein numerisches Modell entwickelt, welches das Spannungs-Verformungs-, Festigkeits- und Bruchverhalten eines Sandsteins im direkten Zug- und im einaxialen Druckversuch sowie in dreiaxialen Druckprüfungen zufriedenstellend wiedergibt. Das Festigkeitsverhalten des entwickelten Modells wurde in Mehrstufentests mit unterschiedlichen Spannungspfaden analysiert und mit den entsprechenden Laborbefunden verglichen.

碩士
國立屏東科技大學
水土保持系所
102
This study on the slate creep behavior in Gaoping River was performed by the discrete element method. In this paper, Based on the rockmass creep patterns by Chigira (1992) and Richard (2000) and the field observation, the essential boundary and material condition and the mechanism of the creep are studied by UDEC simulation. The bending fold of a creep rockmass occurs frequently in the lower angle jointed slope while the numerical simulation obtain the similar result. Furthermore, the creep pattern converted as under the different materials and terrain conditions. The results shows when the joints inclination angle less than 45 degrees, caused by gravity, the rock mass yields the bending fold type of creep; On the other hand, for the one’s inclination angle is greater than 45 degrees, the rock mass yields the topple on the surface of slope. In the weak rocks with strong joints condition, the rock mass will yields a bending fold; next, in the strong rock with weak joints condition, the rock mass yields a block sliding; oppositely, in the weak rock with weak joints condition, the rock mass will slide and bending along weak plane; last, the anti-dip slope with steep angle in the strong rock with weak joints, The rock mass experiences the topple and falling down to slope in block. River erosion effect is also studied by cutting the foot of river bank in several slopes in numerical model. The results show the multiple small folds occurred from the toe to top of slope under the 30 degrees of river bank. Moreover, the discrete block sliding on surface of slope, and the squeezing at toe of slope are occurring under the 45 degrees of river bank.

As a good and competent surgical simulator, it should provide surgeons with visual, tactile and behavioral illusion of reality. In literature, methods for object deformation range from non-physically based models to physically based models. Early works of non-physically based models focused on pure geometrical models that were originally employed in computer-aided design. These methods could be used to produce vivid deformable effects in computer animation. However, the soft tissue simulation in surgical applications requires more realistic models based on physical properties of human tissues. As a result, the mass-spring model and the finite element model have become the most popular representations for deformable organs in surgical simulation. Our research focuses on the real-time soft tissue deformable model based on the finite element method for surgical application.
Extended from the hybrid condensed finite element model, an interactive hybrid condensed model with hardware acceleration by the graphics processing unit (GPU) is proposed. Two methods are developed in order to map the data onto the GPU in accordance with the application data structure. The performance of the primary calculation task in the solver is enhanced. Furthermore, an improved scheme is presented to conduct the newly applied forces induced by dragging or poking operations in the non-operational region.
In the thesis, new approaches to establish a physically based model for soft tissue deformation and cutting in virtual-reality-based simulators are proposed. A deformable model, called the hybrid condensed finite element model, based on the volumetric finite element method is presented. By this method, three-dimensional organs can be represented as tetrahedral meshes, divided into two regions: the operational region and the non-operational one. Different methods treat the regions with different properties in order to balance the computational time and the level of the simulation realism. The condensation technique is applied to only involve the calculation of the surface nodes in the non-operational region while the fully calculation of the volumetric deformation is processed in the operational part. This model guarantees the smooth simulation of cutting operation with the exact cutting path when users manipulate a virtual scalpel. Moreover, we discuss the relevant aspects on what affect the efficiency of implementing the finite element method, as well as the issues considered for choosing the effective solving method to our problem. Three numerical methods have been examined in our model.
Surgical simulator, which benefits from virtual reality techniques, presents a realistic and feasible approach to train inexperienced surgeons within a safe environment. It plays more and more important role in medical field and also changes the world of surgical training. Especially, the minimally invasive microsurgery, which offers patients various attractive advantages over the traditional surgery, has been widely used in otolaryngology, gastroenterology, gynecology and neurology in the last two decades.
Through the combination of these approaches, a physically based model which allows users to freely perform the soft tissue cutting and detecting, such as poking or dragging operations, with soft tissue deformation is achieved in real-time.
Wu Wen.
"August 2006."
Adviser: Pheng Ann Heng.
Source: Dissertation Abstracts International, Volume: 68-03, Section: B, page: 1745.
Thesis (Ph.D.)--Chinese University of Hong Kong, 2006.
Includes bibliographical references (p. 112-127).
Electronic reproduction. Hong Kong : Chinese University of Hong Kong, [2012] System requirements: Adobe Acrobat Reader. Available via World Wide Web.
Electronic reproduction. [Ann Arbor, MI] : ProQuest Information and Learning, [200-] System requirements: Adobe Acrobat Reader. Available via World Wide Web.
Abstracts in English and Chinese.
School code: 1307.

Currently, various needs from industry, science and national defense strategy demand materials with cutting-edge ultra-high temperature performances. Typical applications of ultra-high temperature materials (UHTMs) are supersonic airplanes, gas turbines and rocket nozzles which usually require continuous service of critical components at temperatures higher than 1600°C. Creep resistance is a critical criterion in designing materials for these applications. Traditional creep characterization methods, however, due to limitations on cost, accuracy and most importantly temperature capability, gradually emerge as a bottleneck. Since 2004, a group of researchers in the University of Massachusetts, Amherst have been designing a new high temperature characterization scheme that can break through the limits of traditional methods. Their method is based on non-contact creep tests conducted with Electrostatic levitation (ESL) facilities in NASA Marshall Space Flight Center in Huntsville Alabama. The tested sample is levitated in electric field and is heated as well as rotated with specially positioned laser beam. After certain amount of time, the sample deforms under centripetal forces. By comparison of the shape of the deformed sample with results from finite element simulation, creep behavior of the tested material can be characterized. Based on the same theory, this thesis presents a computational creep characterization method based on non-contact method. A finite element model was built to simulate non-contact creep behavior and results were compared to ESL experiments to determine the creep characteristic. This method was validated both theoretically and numerically and then applied to creep characterization of a promising ultra-high temperature composite from General electric (GE).

Wireline formation testers are widely used to measure in-situ fluid pressure, to retrieve reservoir fluid samples, and to estimate formation mobility. However, formation-tester measurements are invariably influenced by mud-filtrate invasion due to drilling overbalance pressure, thereby affecting the acquisition of uncontaminated fluid samples and the estimation of in-situ petrophysical properties. Moreover, in cases of stress-sensitive formations, rock mechanical deformation may take place due to the combined effects of in-situ stress, wellbore stress imposed by mud overbalance, and wellbore pressure exerted by the formation tester itself. The latter deformation causes near-borehole perturbations of porosity and permeability that are evidenced by pressure transients measured during build-up and shut-in stages of formation testing, especially when using dual-packer pressure probes. If unaccounted for, such perturbations can also bias the estimation of in-situ fluid and petrophysical properties. Conversely, the detection and quantification of elastic mechanical deformation effects on measured pressure transients can be used to infer the underlying rock elastic and petrophysical properties of the stressed formation. The purpose of this dissertation is twofold: (a) to quantify the relative effects of mud-filtrate invasion and geomechanical deformation on pressure-transient measurements acquired with dual-packer formation testers, with special emphasis on the appraisal of near-borehole porosity and permeability enhancement due to elastic mechanical deformation, and (b) to develop a new method to estimate elastic and petrophysical properties of rock formations from dual-packer pressure transients acquired in mechanically deformable rocks. Numerical simulations of mud-filtrate invasion are performed with an axialsymmetric two-phase (water-oil) method that enforces the specific boundary and source conditions of a wellbore that penetrates horizontal layers. Simulations are performed in a cylindrical system of coordinates using finite differences together with an implicit-pressure, explicit-saturation time-marching approach that also incorporates the dynamic conditions of immiscible mudcake growth due to filtration of solids at the wellbore. Laboratory experiments are conducted to further study pressure transients due to formation testing in the presence of invasion with water-base mud. Experiments include the effects of both mud circulation and mudcake on pressure-transient measurements and are performed on a variety of rock-core samples. Measurements are successfully validated with both the developed simulator and a commercial simulator, thereby lending credence to the assumed model of dynamic solid filtration. The developed mud-filtrate fluid-flow simulator is coupled with a finite-element code that assumes 2D axial-symmetric linear elasticity to quantify geomechanical deformation. Coupling of mechanical deformation with variations of porosity and permeability assumes a staggered-in-time, iteratively coupled volumetric model. We assume a dual-packer formation tester to quantify elastic deformation effects in stress-sensitive formations as a preamble to estimating in-situ elastic and petrophysical properties. It is shown that near-wellbore spatial variations of porosity and permeability due to mechanical deformation can bias the corresponding pressure-transient measurements acquired with the dual-packer formation-tester. The degree of biasing depends on the rigidity of the stressed formation. Finally, we develop a method to estimate in-situ petrophysical and elastic rock properties from pressure-transient measurements acquired with formation-testers in mechanically deformable rocks. Petrophysical and elastic properties will change in both time and space depending on the time evolution of the conditions that influence mechanical deformation. We use a commercial reservoir simulator to calculate pressure transients due to fluid pumpout in the presence of both invasion and mechanical deformation. A pre-stressed initial condition due to mud overbalance is assumed with incremental deformation due to surface force applied by the packers or probes, and active flow imposed by the formation-tester. In so doing, we consider pressure data sets acquired with both flow and observation probes during draw-down and build-up periods. For cases where a-priori information can be sufficiently constrained, our estimation method provides reliable and accurate estimates of petrophysical and elastic properties in the presence of moderate levels of random noise.
text

Infrastructure is developing at an extremely fast pace which includes construction of metros, underground storage places, railway bridges, caverns and tunnels. Very often these structures are found in or on the rock masses. Rock masses are seldom found in nature without joints or discontinuities. Jointed rocks are characterized by the presence of inherent discontinuities of varied sizes with different orientations and intensities, which can have significant effect on their mechanical response. Constructions involving jointed rocks often become challenging jobs for Civil Engineers as the instability of slopes or excavations in these jointed rocks poses serious concerns, sometimes leading to the failure of structures built on them. Experimental investigations on jointed rock masses are not always feasible and pose formidable problems to the engineers. Apart from the technical difficulties of extracting undisturbed rock samples, it is very expensive and time consuming to conduct the experiments on jointed rock masses of huge dimensions. The most popular methods of evaluating the rock mass behaviour are the Numerical methods. In this thesis, numerical modelling of jointed rock masses is carried out using computer program FLAC (Fast Lagrangian Analysis of Continua). The objective of the present study is to study the effect of various joint parameters on the response of jointed rock masses in static as well as seismic shaking conditions. This is achieved through systematic series of numerical simulations of jointed rocks in triaxial compression, in underground openings and in large rock slopes. This thesis is an attempt to study the individual effect of different joint parameters on the rock mass behaviour and to integrate these results to provide useful insight into the behaviour of jointed rock mass under various joint conditions. In practice, it is almost impossible to explore all of the joint systems or to investigate all their mechanical characteristics and implementing them explicitly in the model. In these cases, the use of the equivalent continuum model to simulate the behaviour of jointed rock masses could be valuable. Hence this approach is mainly used in this thesis. Some numerical simulations with explicitly modelled joints are also presented for comparison with the continuum modelling. The applicability of Artificial Neural Networks for the prediction of stress-strain response of jointed rocks is also explored. Static, pseudo-static and dynamic analyses of a large rock slope in Himalayas is carried out and parametric seismic analysis of rock slope is carried out with varying input shaking, material damping and shear strength parameters. Results from the numerical studies showed that joint inclination is the most influencing parameter for the jointed rock mass behaviour. Rock masses exhibit lowest strength at critical angle of joint inclination and the deformations around excavations will be highest when the joints are inclined at an angle close to the critical angle. However at very high confining pressures, the influence of joint inclination gets subdued. Under seismic base shaking conditions, the deformations of rock masses largely depend on the acceleration response with time, frequency content and duration rather than the peak amplitude or the magnitude of earthquake. All these aspects are discussed in the light of results from numerical studies presented in this thesis.

Dreiaxiale Druckprüfungen können als Einstufenversuche, als Mehrstufenversuche oder als Versuche mit kontinuierlichen Bruchzuständen ausgeführt werden. Bei der Anwendung der Mehrstufentechnik ergeben sich insbesondere Fragestellungen hinsichtlich der richtigen Wahl des Umschaltpunktes und des optimalen Verlaufs des Spannungspfades zwischen den einzelnen Versuchsstufen. Fraglich beim Versuch mit kontinuierlichen Bruchzuständen bleibt, ob im Versuchsverlauf tatsächlich Spannungszustände erfasst werden, welche die Höchstfestigkeit des untersuchten Materials repräsentieren. Die Dissertation greift diese Fragestellungen auf, ermöglicht den Einstieg in die beschriebene Thematik und schafft die Voraussetzungen, die zur Lösung der aufgeführten Problemstellungen notwendig sind. Auf der Grundlage einer umfangreichen Datenbasis gesteinsmechanischer und petrophysikalischer Kennwerte wurde ein numerisches Modell entwickelt, welches das Spannungs-Verformungs-, Festigkeits- und Bruchverhalten eines Sandsteins im direkten Zug- und im einaxialen Druckversuch sowie in dreiaxialen Druckprüfungen zufriedenstellend wiedergibt. Das Festigkeitsverhalten des entwickelten Modells wurde in Mehrstufentests mit unterschiedlichen Spannungspfaden analysiert und mit den entsprechenden Laborbefunden verglichen.