Volume 29, Issue1 (January 2005)
Articles in the Current Issue:
Research Article
Homogenization framework for three-dimensional elastoplastic finite element analysis of a grouted pipe-roofing reinforcement method for tunnelling
NASA Astrophysics Data System (ADS)
Bae, G. J.; Shin, H. S.; Sicilia, C.; Choi, Y. G.; Lim, J. J.
2005-01-01
This paper deals with the grouted pipe-roofing reinforcement method that is used in the construction of tunnels through weak grounds. This system consists on installing, prior to the excavation of a length of tunnel, an array of pipes forming a kind of umbrella above the area to be excavated. In some cases, these pipes are later used to inject grout to strengthen the ground and connect the pipes.This system has proven to be very efficient in reducing tunnel convergence and water inflow when tunnelling through weak grounds. However, due to the geometrical and mechanical complexity of the problem, existing finite element frameworks are inappropriate to simulate tunnelling using this method.In this paper, a mathematical framework based on a homogenization technique to simulate grouted pipe-roofing reinforced ground and its implementation into a 3-D finite element programme that can consider stage construction situations are presented. The constitutive model developed allows considering the main design parameters of the problem and only requires geometrical and mechanical properties of the constituents. Additionally, the use of a homogenization approach implies that the generation of the finite element mesh can be easily produced and that re-meshing is not required as basic geometrical parameters such as the orientation of the pipes are changed.
Orientation Effects in Fault Reactivation in Geological CO2 Sequestration
NASA Astrophysics Data System (ADS)
Castelletto, N.; Ferronato, M.; Gambolati, G.; Janna, C.; Teatini, P.
2012-12-01
Geological CO2 sequestration remains one of the most promising option for reducing the greenhouse gases emission. The accurate simulation of the complex coupled physical processes occurring during the injection and the post-injection stage represents a key issue for investigating the feasibility and the safety of the sequestration. The fluid-dynamical and geochemical aspects related to sequestering CO2 underground have been widely debated in the scientific literature over more than one decade. Recently, the importance of geomechanical processes has been widely recognized. In the present modeling study, we focus on fault reactivation induced by injection, an essential aspect for the evaluation of CO2 sequestration projects that needs to be adequately investigated to avoid the generation of preferential leaking path for CO2 and the related risk of induced seismicity. We use a geomechanical model based on the structural equations of poroelasticity solved by the Finite Element (FE) - Interface Element (IE) approach. Standard FEs are used to represent a continuum, while IEs prove especially suited to assess the relative displacements of adjacent elements such as the opening and slippage of existing faults or the generation of new fractures [1]. The IEs allow for the modeling of fault mechanics using an elasto-plastic constitutive law based on the Mohr-Coulomb failure criterion. We analyze the reactivation of a single fault in a synthetic reservoir by varying the fault orientation and size, hydraulic conductivity of the faulted zone, initial vertical and horizontal stress state and Mohr-Coulomb parameters (i.e., friction angle and cohesion). References: [1] Ferronato, M., G. Gambolati, C. Janna, and P. Teatini (2008), Numerical modeling of regional faults in land subsidence prediction above gas/oil reservoirs, Int. J. Numer. Anal. Methods Geomech., 32, 633-657.
The rise and fall of axial highs during ridge jumps
NASA Astrophysics Data System (ADS)
Shah, Anjana K.; Buck, W. Roger
2006-08-01
We simulate jumps of ocean spreading centers with axial high topography using elastoplastic thin plate flexure models. Processes considered include ridge abandonment, the breaking of a stressed plate on the ridge flank, and renewed spreading at the site of this break. We compare model results to topography at the East Pacific Rise between 15°25'N and 16°N, where there is strong evidence of a recent ridge jump. At an apparently abandoned ridge, gravity data do not suggest buoyant support of topography. Model deflections during cooling and melt solidification stages of ridge abandonment are of small vertical amplitude because of plate strengthening, resulting in the preservation of a "frozen" fossil high. The present-day high is bounded by slopes with up to a 40% grade, a scenario very difficult to achieve flexurally given generally accepted constraints on lithospheric strength. We model these slopes by assuming that the height at which magma is accreted increases rapidly after the ridge jumps. This increase is attributed to high overburden pressure on melt that resided in an initially deep magma chamber, followed by a rapid increase in temperature and melt supply to the region shortly after spreading began. The high is widest at the segment center, suggesting that magmatic activity began near the center of the segment, propagated south and then north. The mantle Bouguer anomaly exhibits a "bull's-eye" pattern centered at the widest part of the high, but the depth of the axis is nearly constant along the length of the segment. We reconcile these observations by assigning different cross-axis widths to a low-density zone within the crust.
A theoretical study on pure bending of hexagonal close-packed metal sheet
NASA Astrophysics Data System (ADS)
Mehrabi, Hamed; Yang, Chunhui
2018-05-01
Hexagonal close-packed (HCP) metals have quite different mechanical behaviours in comparison to conventional cubic metals such as steels and aluminum alloys [1, 2]. They exhibit a significant tension-compression asymmetry in initial yielding and subsequent plastic hardening. The reason for this unique behaviour can be attributed to their limited symmetric crystal structure, which leads to twining deformation [3-5]. This unique behaviour strongly influences sheet metal forming of such metals, especially for roll forming, in which the bending is dominant. Hence, it is crucial to represent constitutive relations of HCP metals for accurate estimation of bending moment-curvature behaviours. In this paper, an analytical model for asymmetric elastoplastic pure bending with an application of Cazacu-Barlat asymmetric yield function [6] is presented. This yield function considers the asymmetrical tension-compression behaviour of HCP metals by using second and third invariants of the stress deviator tensor and a specified constant, which can be expressed in terms of uniaxial yield stresses in tension and compression. As a case study, the analytical model is applied to predict the moment-curvature behaviours of AZ31B magnesium alloy sheets under uniaxial loading condition. Furthermore, the analytical model is implemented as a user-defined material through the UMAT interface in Abaqus [7, 8] for conducting pure bending simulations. The results show that the analytical model can reasonably capture the asymmetric tension-compression behaviour of the magnesium alloy. The predicted moment-curvature behaviour has good agreement with the experimental results. Furthermore, numerical results show a better accuracy by the application of the Cazacu-Barlat yield function than those using the von-Mises yield function, which are more conservative than analytical results.
A Transversely Isotropic Thermo-mechanical Framework for Oil Shale
NASA Astrophysics Data System (ADS)
Semnani, S. J.; White, J. A.; Borja, R. I.
2014-12-01
The present study provides a thermo-mechanical framework for modeling the temperature dependent behavior of oil shale. As a result of heating, oil shale undergoes phase transformations, during which organic matter is converted to petroleum products, e.g. light oil, heavy oil, bitumen, and coke. The change in the constituents and microstructure of shale at high temperatures dramatically alters its mechanical behavior e.g. plastic deformations and strength, as demonstrated by triaxial tests conducted at multiple temperatures [1,2]. Accordingly, the present model formulates the effects of changes in the chemical constituents due to thermal loading. It is well known that due to the layered structure of shale its mechanical properties in the direction parallel to the bedding planes is significantly different from its properties in the perpendicular direction. Although isotropic models simplify the modeling process, they fail to accurately describe the mechanical behavior of these rocks. Therefore, many researchers have studied the anisotropic behavior of rocks, including shale [3]. The current study presents a framework to incorporate the effects of transverse isotropy within a thermo-mechanical formulation. The proposed constitutive model can be readily applied to existing finite element codes to predict the behavior of oil shale in applications such as in-situ retorting process and stability assessment in petroleum reservoirs. [1] Masri, M. et al."Experimental Study of the Thermomechanical Behavior of the Petroleum Reservoir." SPE Eastern Regional/AAPG Eastern Section Joint Meeting. Society of Petroleum Engineers, 2008. [2] Xu, B. et al. "Thermal impact on shale deformation/failure behaviors---laboratory studies." 45th US Rock Mechanics/Geomechanics Symposium. American Rock Mechanics Association, 2011. [3] Crook, AJL et al. "Development of an orthotropic 3D elastoplastic material model for shale." SPE/ISRM Rock Mechanics Conference. Society of Petroleum Engineers, 2002.
Gao, Chunxia; Rahaman, Mohamed N; Gao, Qiang; Teramoto, Akira; Abe, Koji
2013-07-01
The development of inorganic-organic hybrid scaffolds with controllable degradation and bioactive properties is receiving considerable interest for bone and tissue regeneration. The objective of this study was to create hybrid scaffolds of gelatin and bioactive glass (BG) with a controlled, three-dimensional (3D) architecture by a combined sol-gel and robotic deposition (robocasting) method and evaluate their mechanical response, bioactivity, and response to cells in vitro. Inks for robotic deposition of the scaffolds were prepared by dissolving gelatin in a sol-gel precursor solution of the bioactive glass (70SiO2 -25CaO-5P2 O5 ; mol%) and aging the solution to form a gel with the requisite viscosity. After drying and crosslinking, the gelatin-BG scaffolds, with a grid-like architecture (filament diameter ∼350 µm; pore width ∼550 µm), showed an elasto-plastic response, with a compressive strength of 5.1 ± 0.6 MPa, in the range of values for human trabecular bone (2-12 MPa). When immersed in phosphate-buffered saline, the crosslinked scaffolds rapidly absorbed water (∼440% of its dry weight after 2 h) and showed an elastic response at deformations up to ∼60%. Immersion of the scaffolds in a simulated body fluid resulted in the formation of a hydroxyapatite-like surface layer within 5 days, indicating their bioactivity in vitro. The scaffolds supported the proliferation, alkaline phosphatase activity, and mineralization of osteogenic MC3T3-E1 cells in vitro, showing their biocompatibility. Altogether, the results indicate that these gelatin-BG hybrid scaffolds with a controlled, 3D architecture of inter-connected pores have potential for use as implants for bone regeneration. Copyright © 2012 Wiley Periodicals, Inc.
Mathematical modeling of thermal stresses in basic oxygen furnace hood tubes
NASA Astrophysics Data System (ADS)
Samarasekera, I. V.
1985-06-01
The stress-strain history of Basic Oxygen Furnace hood tubes during thermal cycling has been computed using heat flow and stress analyses. The steady-state temperature distribution in a transverse section of the tube was computed at a location where gas temperature in the hood could be expected to be a maximum. Calculations were performed for peak gas temperatures in the range 1950 to 2480 °C (3500 to 4500 °F). The stress-strain history of an element of material located at the center of the tube hot face was traced for three consecutive cycles using elasto-plastic finite-element analysis. It has been shown that the state of stress in the element alternates between compression and tension as the tube successively heats and cools. Yielding and plastic flow occurs at the end of each half of a given cycle. It was postulated that owing to repctitive yielding, plastic strain energy accumulates causing failure of the tubes by fatigue in the low cycle region. Using fatigue theory a conservative estimate for tube life was arrived at. In-plant observations support this mechanism of failure, and the number of cycles within which tube cracking was observed compares reasonably with model predictions. Utilizing the heat flow and stress models it was recommended that tube life could be enhanced by changing the tube material to ARMCO 17-4 pH or AISI 405 steel or alternatively reconstructing hoods with AISI 316L tubes of reduced thickness. These recommendations were based on the criterion that low-cycle fatigue failure could be averted if the magnitude of the cyclic strain could be reduced or if macroscopic plastic flow could be prevented.
Mechano-Chemical Interactions at Cement-Geomaterial Interfaces in Repository and Borehole Scenarios
NASA Astrophysics Data System (ADS)
Mohagheghi, J. R.; Dewers, T. A.; Matteo, E. N.; Heath, J. E.; Jove Colon, C. F.; Fuller, T.
2017-12-01
A number of factors negatively affect wellbore integrity including interactions at boundaries between cement and surrounding geomaterial. These include mechanical and chemical mechanisms that can lead to wellbore failure. To examine these interactions, potential coupling, and pathways to failure, we discuss progress on an experimental and modeling study involving cement-clay and cement-salt interfaces. A sample shotcrete-bentonite interface from the FEBEX heater test at the Grimsel Test Site in Switzerland is examined using multi-beam scanning electron microscopy (mSEM) at 4 nm resolution over an area 10's of square millimeters. We examine changes in alteration as manifested by pore structural changes as a function of distance from the interface. A parallel effort examines time-dependent changes in interface structure in cement cores in a triaxial coreholder. Cores are exposed to conditions of 70oC, 14 MPa pressure, and small differential loads, with degradation being monitored by effluent pH, pulse-echo ultrasonics, and piston displacement (measuring sample shortening). We will measure the mechanical consequences of interface alteration using nano-indentation. Experimental results are being incorporated as a validation effort in a coupled reactive-transport mechanics model linking the Sandia ALBANY finite element code, the KAYENTA elasto-plastic constitutive model, with the reactive transport code PFLOTRAN. Plans call to apply the model to understanding the evolution of the FEBEX sample, as well as a cement-salt sample from the Waste Isolation Pilot Plant in Carlsbad, New Mexico. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology and Engineering Solutions of Sandia LLC, a wholly owned subsidiary of Honeywell International Inc. for the U.S. Department of Energy's National Nuclear Security Administration under contract DE-NA0003525. SAND 2017-8277 A
DOE Office of Scientific and Technical Information (OSTI.GOV)
Hilairet, Nadège; Wang, Yanbin; Sanehira, Takeshi
2012-03-15
Polycrystalline samples of San Carlos olivine were deformed at high-pressure (2.8-7.8 GPa), high-temperature (1153 to 1670 K), and strain rates between 7.10{sup -6} and 3.10{sup -5} s{sup -1}, using the D-DIA apparatus. Stress and strain were measured in situ using monochromatic X-rays diffraction and imaging, respectively. Based on the evolution of lattice strains with total bulk strain and texture development, we identified three deformation regimes, one at confining pressures below 3-4 GPa, one above 4 GPa, both below 1600 K, and one involving growth of diffracting domains associated with mechanical softening above {approx}1600 K. The softening is interpreted as enhancedmore » grain boundary migration and recovery. Below 1600 K, elasto-plastic self-consistent analysis suggests that below 3-4 GPa, deformation in olivine occurs with large contribution from the so-called 'a-slip' system [100](010). Above {approx}4 GPa, the contribution of the a-slip decreases relative to that of the 'c-slip' [001](010). This conclusion is further supported by texture refinements. Thus for polycrystalline olivine, the evolution in slip systems found by previous studies may be progressive, starting from as low as 3-4 GPa and up to 8 GPa. During such a gradual change, activation volumes measured on polycrystalline olivine cannot be linked to a particular slip system straightforwardly. The quest for 'the' activation volume of olivine at high pressure should cease at the expense of detailed work on the flow mechanisms implied. Such evolution in slip systems should also affect the interpretation of seismic anisotropy data in terms of upper mantle flow between 120 and 300 km depth.« less
NASA Astrophysics Data System (ADS)
Kiliclar, Yalin; Laurischkat, Roman; Vladimirov, Ivaylo N.; Reese, Stefanie
2011-08-01
The presented project deals with a robot based incremental sheet metal forming process, which is called roboforming and has been developed at the Chair of Production Systems. It is characterized by flexible shaping using a freely programmable path-synchronous movement of two industrial robots. The final shape is produced by the incremental infeed of the forming tool in depth direction and its movement along the part contour in lateral direction. However, the resulting geometries formed in roboforming deviate several millimeters from the reference geometry. This results from the compliance of the involved machine structures and the springback effects of the workpiece. The project aims to predict these deviations caused by resiliences and to carry out a compensative path planning based on this prediction. Therefore a planning tool is implemented which compensates the robots's compliance and the springback effects of the sheet metal. The forming process is simulated by means of a finite element analysis using a material model developed at the Institute of Applied Mechanics (IFAM). It is based on the multiplicative split of the deformation gradient in the context of hyperelasticity and combines nonlinear kinematic and isotropic hardening. Low-order finite elements used to simulate thin sheet structures, such as used for the experiments, have the major problem of locking, a nonphysical stiffening effect. For an efficient finite element analysis a special solid-shell finite element formulation based on reduced integration with hourglass stabilization has been developed. To circumvent different locking effects, the enhanced assumed strain (EAS) and the assumed natural strain (ANS) concepts are included in this formulation. Having such powerful tools available we obtain more accurate geometries.
Kelly, Nicola; McGarry, J Patrick
2012-05-01
The inelastic pressure dependent compressive behaviour of bovine trabecular bone is investigated through experimental and computational analysis. Two loading configurations are implemented, uniaxial and confined compression, providing two distinct loading paths in the von Mises-pressure stress plane. Experimental results reveal distinctive yielding followed by a constant nominal stress plateau for both uniaxial and confined compression. Computational simulation of the experimental tests using the Drucker-Prager and Mohr-Coulomb plasticity models fails to capture the confined compression behaviour of trabecular bone. The high pressure developed during confined compression does not result in plastic deformation using these formulations, and a near elastic response is computed. In contrast, the crushable foam plasticity models provide accurate simulation of the confined compression tests, with distinctive yield and plateau behaviour being predicted. The elliptical yield surfaces of the crushable foam formulations in the von Mises-pressure stress plane accurately characterise the plastic behaviour of trabecular bone. Results reveal that the hydrostatic yield stress is equal to the uniaxial yield stress for trabecular bone, demonstrating the importance of accurate characterisation and simulation of the pressure dependent plasticity. It is also demonstrated in this study that a commercially available trabecular bone analogue material, cellular rigid polyurethane foam, exhibits similar pressure dependent yield behaviour, despite having a lower stiffness and strength than trabecular bone. This study provides a novel insight into the pressure dependent yield behaviour of trabecular bone, demonstrating the inadequacy of uniaxial testing alone. For the first time, crushable foam plasticity formulations are implemented for trabecular bone. The enhanced understanding of the inelastic behaviour of trabecular bone established in this study will allow for more realistic simulation of orthopaedic device implantation and failure. Copyright © 2011 Elsevier Ltd. All rights reserved.
3D Finite Element Analysis of Particle-Reinforced Aluminum
NASA Technical Reports Server (NTRS)
Shen, H.; Lissenden, C. J.
2002-01-01
Deformation in particle-reinforced aluminum has been simulated using three distinct types of finite element model: a three-dimensional repeating unit cell, a three-dimensional multi-particle model, and two-dimensional multi-particle models. The repeating unit cell model represents a fictitious periodic cubic array of particles. The 3D multi-particle (3D-MP) model represents randomly placed and oriented particles. The 2D generalized plane strain multi-particle models were obtained from planar sections through the 3D-MP model. These models were used to study the tensile macroscopic stress-strain response and the associated stress and strain distributions in an elastoplastic matrix. The results indicate that the 2D model having a particle area fraction equal to the particle representative volume fraction of the 3D models predicted the same macroscopic stress-strain response as the 3D models. However, there are fluctuations in the particle area fraction in a representative volume element. As expected, predictions from 2D models having different particle area fractions do not agree with predictions from 3D models. More importantly, it was found that the microscopic stress and strain distributions from the 2D models do not agree with those from the 3D-MP model. Specifically, the plastic strain distribution predicted by the 2D model is banded along lines inclined at 45 deg from the loading axis while the 3D model prediction is not. Additionally, the triaxial stress and maximum principal stress distributions predicted by 2D and 3D models do not agree. Thus, it appears necessary to use a multi-particle 3D model to accurately predict material responses that depend on local effects, such as strain-to-failure, fracture toughness, and fatigue life.
Slice-push, formation of grooves and the scale effect in cutting.
Atkins, A G
2016-06-06
Three separate aspects of cutting are investigated which complement other papers on the mechanics of separation processes presented at this interdisciplinary Theo Murphy meeting. They apply in all types of cutting whether blades are sharp or blunt, and whether the material being cut is 'hard, stiff and strong' or 'soft, compliant and weak'. The first topic discusses why it is easier to cut when there is motion along (parallel to) the blade as well motion across (perpendicular to) the cutting edge, and the analysis is applied to optimization of blade geometries to produce minimum cutting forces and hence minimum damage to cut surfaces. The second topic concerns cutting with more than one edge with particular application to the formation of grooves in surfaces by hard pointed tools. The mechanics are investigated and applied to the topic of abrasive wear by hard particles. Traditional analyses say that abrasive wear resistance increases monotonically with the hardness of the workpiece, but we show that the fracture toughness of the surface material is also important, and that behaviour is determined by the toughness-to-hardness ratio rather than hardness alone. Scaling forms the third subject. As cutting is a branch of elasto-plastic fracture mechanics, cube-square energy scaling applies in which the important length scale is (ER/k (2)), where E is Young's modulus, R is the fracture toughness and k is the shear yield strength. Whether, in cutting, material is removed as ductile ribbons, as semi-ductile discontinuous chips, or by brittle 'knocking lumps out' is shown to depend on the depth of cut relative to this characteristic length parameter. Scaling in biology is called allometry and its relationship with engineering scaling is discussed. Some speculative predictions are made in relation to the action of teeth on food.
Experiment Evaluation of Bifurcation in Sands
NASA Technical Reports Server (NTRS)
Alshibi, Khalid A.; Sture, Stein
2000-01-01
The basic principles of bifurcation analysis have been established by several investigators, however several issues remain unresolved, specifically how do stress level, grain size distribution, and boundary conditions affect general bifurcation phenomena in pressure sensitive and dilatant materials. General geometrical and kinematics conditions for moving surfaces of discontinuity was derived and applied to problems of instability of solids. In 1962, the theoretical framework of bifurcation by studying the acceleration waves in elasto-plastic (J2) solids were presented. Bifurcation analysis for more specific forms of constitutive behavior was examined by studying localization in pressure-sensitive, dilatant materials, however, analyses were restricted to plane deformation states only. Bifurcation analyses were presented and applied to predict shear band formations in sand under plane strain condition. The properties of discontinuous bifurcation solutions for elastic-plastic solids under axisymmetric and plane strain loading conditions were studied. The study focused on theory, but also references and comparisons to experiments were made. The current paper includes a presentation of a summary of bifurcation analyses for biaxial and triaxial (axisymmetric) loading conditions. The Coulomb model is implemented using incremental piecewise scheme to predict the constitutive relations and shear band inclination angles. Then, a comprehensive evaluation of bifurcation phenomena is presented based on data from triaxial experiments performed under microgravity conditions aboard the Space Shuttle under very low effective confining pressure (0.05 to 1.30 kPa), in which very high peak friction angles (47 to 75 degrees) and dilatancy angles (30 to 31 degrees) were measured. The evaluation will be extended to include biaxial experiments performed on the same material under low (10 kPa) and moderate (100 kPa) confining pressures. A comparison between the behavior under biaxial and triaxial loading conditions will be presented, and related issues concerning influence of confining pressure will be discussed.
NASA Astrophysics Data System (ADS)
Ye, Shujun; Franceschini, Andrea; Zhang, Yan; Janna, Carlo; Gong, Xulong; Yu, Jun; Teatini, Pietro
2018-03-01
Initially observed in the semiarid basins of southwestern USA, earth fissures due to aquifer over-exploitation are presently threatening a large number of subsiding basins in various countries worldwide. Different mechanics have been proposed to explain this process, such as differential compaction, horizontal movements, and fault reactivation. Numerical modeling and prediction of this major geohazard caused by overuse of groundwater resources are challenging because of two main requirements: shifting from the classical continuous to discontinuous geomechanics and incorporating two-dimensional features (the earth fissures) into large three-dimensional (3-D) modeling domain (the subsiding basin). In this work, we proposed a novel modeling approach to simulate earth fissure generation and propagation in 3-D complex geological settings. A nested two-scale approach associated with an original nonlinear elastoplastic finite element/interface element simulator allows modeling the mechanics of earth discontinuities, in terms of both sliding and opening. The model is applied on a case study in Wuxi, China, where groundwater pumping between 1985 and 2004 has caused land subsidence larger than 2 m. The model outcomes highlight that the presence of a shallow (˜80 m deep) bedrock ridge crossing the Yangtze River delta is the key factor triggering the earth fissure development in this area. Bending of the alluvial deposits around the ridge tip and shear stress due to the uneven piezometric change and asymmetrical shape of the bedrock have caused the earth fissure to onset at the land surface and propagate downward to a maximum depth of about 20-30 m. Maximum sliding and opening are computed in the range of 10-40 cm, in agreement with the order of magnitude estimated in the field.
NASA Astrophysics Data System (ADS)
Lu, Yaqing; Hui, Hu; Gong, Jianguo
2018-05-01
Austenitic stainless steel is widely used in pressure vessels for the storage and transportation of liquid gases such as liquid nitrogen, liquid oxygen, and liquid hydrogen. Cryogenic pressure vessel manufacturing uses cold stretching technology, which relies heavily on welding joint performance, to construct lightweight and thin-walled vessels. Residual stress from welding is a primary factor in cases of austenitic stainless steel pressure vessel failure. In this paper, on the basis of Visual Environment 10.0 finite element simulation technology, the residual stress resulting from different welding strength matching coefficients (0.8, 1, 1.2, 1.4) for two S30408 plates welded with three-pass butt welds is calculated according to thermal elastoplastic theory. In addition, the stress field was calculated under a loading as high as 410 MPa and after the load was released. Path 1 was set to analyze stress along the welding line, and path 2 was set to analyze stress normal to the welding line. The welding strength matching coefficient strongly affected both the longitudinal residual stress (center of path 1) and the transverse residual stress (both ends of path 1) after the welding was completed. However, the coefficient had little effect on the longitudinal and transverse residual stress of path 2. Under the loading of 410 MPa, the longitudinal and transverse stress decreased and the stress distribution, with different welding strength matching coefficients, was less diverse. After the load was released, longitudinal and transverse stress distribution for both path 1 and path 2 decreased to a low level. Cold stretching could reduce the effect of residual stress to various degrees. Transverse strain along the stretching direction was also taken into consideration. The experimental results validated the reliability of the partial simulation.
Random lattice structures. Modelling, manufacture and FEA of their mechanical response
NASA Astrophysics Data System (ADS)
Maliaris, G.; Sarafis, I. T.; Lazaridis, T.; Varoutoglou, A.; Tsakataras, G.
2016-11-01
The implementation of lightweight structures in various applications, especially in Aerospace/ Automotive industries and Orthopaedics, has become a necessity due to their exceptional mechanical properties with respect to reduced weight. In this work we present a Voronoi tessellation based algorithm, which has been developed for modelling stochastic lattice structures. With the proposed algorithm, is possible to generate CAD geometry with controllable structural parameters, such as porosity, cell number and strut thickness. The digital structures were transformed into physical objects through the combination of 3D printing technics and investment casting. This process was applied to check the mechanical behaviour of generated digital models. Until now, the only way to materialize such structures into physical objects, was feasible through 3D printing methods such as Selective Laser Sintering/ Melting (SLS/ SLM). Investment casting possesses numerous advantages against SLS or SLA, with the major one being the material variety. On the other hand, several trials are required in order to calibrate the process parameters to have successful castings, which is the major drawback of investment casting. The manufactured specimens were subjected to compression tests, where their mechanical response was registered in the form of compressive load - displacement curves. Also, a finite element model was developed, using the specimens’ CAD data and compression test parameters. The FE assisted calculation of specimen plastic deformation is identical with the one of the physical object, which validates the conclusions drawn from the simulation results. As it was observed, strut contact is initiated when specimen deformation is approximately 5mm. Although FE calculated compressive force follows the same trend for the first 3mm of compression, then diverges because of the elasto-plastic FE model type definition and the occurred remeshing steps.
-
PLAN2D - A PROGRAM FOR ELASTO-PLASTIC ANALYSIS OF PLANAR FRAMES
NASA Technical Reports Server (NTRS)
Lawrence, C.
1994-01-01
PLAN2D is a FORTRAN computer program for the plastic analysis of planar rigid frame structures. Given a structure and loading pattern as input, PLAN2D calculates the ultimate load that the structure can sustain before collapse. Element moments and plastic hinge rotations are calculated for the ultimate load. The location of hinges required for a collapse mechanism to form are also determined. The program proceeds in an iterative series of linear elastic analyses. After each iteration the resulting elastic moments in each member are compared to the reserve plastic moment capacity of that member. The member or members that have moments closest to their reserve capacity will determine the minimum load factor and the site where the next hinge is to be inserted. Next, hinges are inserted and the structural stiffness matrix is reformulated. This cycle is repeated until the structure becomes unstable. At this point the ultimate collapse load is calculated by accumulating the minimum load factor from each previous iteration and multiplying them by the original input loads. PLAN2D is based on the program STAN, originally written by Dr. E.L. Wilson at U.C. Berkeley. PLAN2D has several limitations: 1) Although PLAN2D will detect unloading of hinges it does not contain the capability to remove hinges; 2) PLAN2D does not allow the user to input different positive and negative moment capacities and 3) PLAN2D does not consider the interaction between axial and plastic moment capacity. Axial yielding and buckling is ignored as is the reduction in moment capacity due to axial load. PLAN2D is written in FORTRAN and is machine independent. It has been tested on an IBM PC and a DEC MicroVAX. The program was developed in 1988.
-
Sub-discretized surface model with application to contact mechanics in multi-body simulation
DOE Office of Scientific and Technical Information (OSTI.GOV)
Johnson, S; Williams, J
2008-02-28
The mechanics of contact between rough and imperfectly spherical adhesive powder grains are often complicated by a variety of factors, including several which vary over sub-grain length scales. These include several traction factors that vary spatially over the surface of the individual grains, including high energy electron and acceptor sites (electrostatic), hydrophobic and hydrophilic sites (electrostatic and capillary), surface energy (general adhesion), geometry (van der Waals and mechanical), and elasto-plastic deformation (mechanical). For mechanical deformation and reaction, coupled motions, such as twisting with bending and sliding, as well as surface roughness add an asymmetry to the contact force which invalidatesmore » assumptions for popular models of contact, such as the Hertzian and its derivatives, for the non-adhesive case, and the JKR and DMT models for adhesive contacts. Though several contact laws have been offered to ameliorate these drawbacks, they are often constrained to particular loading paths (most often normal loading) and are relatively complicated for computational implementation. This paper offers a simple and general computational method for augmenting contact law predictions in multi-body simulations through characterization of the contact surfaces using a hierarchically-defined surface sub-discretization. For the case of adhesive contact between powder grains in low stress regimes, this technique can allow a variety of existing contact laws to be resolved across scales, allowing for moments and torques about the contact area as well as normal and tangential tractions to be resolved. This is especially useful for multi-body simulation applications where the modeler desires statistical distributions and calibration for parameters in contact laws commonly used for resolving near-surface contact mechanics. The approach is verified against analytical results for the case of rough, elastic spheres.« less
-
NASA Astrophysics Data System (ADS)
Tong, X.; Lavier, L.
2017-12-01
Cold and warm subduction zones usually have different seismicity and tectonic structure. Aseismic events like episodic tremor and slip (ETS) and slow slip event (SSE) are often observed in warm and young slabs which typically have less megathrust seismicity and smaller seismogenic area (e.g. southwest Japan). On the other hand, cold and old slabs (e.g. Northeast Japan) have more megathrust events and larger seismogenic area and few aseismic events. Recent studies have try to model the differences in seismic behaviors with different approaches, includes rheological heterogeneity (e.g. frictional vs. viscous), petrological heterogeneity (e.g. hydration-dehydration process and mineral phase changes), and the frictional heterogeneity (e.g. rate-and-state dependent friction). Following previous works, we proposed a new model in which the subduction channel has a temperature dependent material assembly which composed of an explicit mixture of basalt/eclogite and mantle peridotite. Our model also take into account rate and state dependent friction and pore fluid pressure. Depending on the temperature, the basalt and peridotite mixture can behave either as an elastoplastic frictional or a Maxwell viscoelastic material. To model the mixture numerically, we use DynEarthSol3D (DES3D). DES3D is a robust, adaptive, multi-dimensional, finite element method solver which has a composite Elasto-Visco-Plastic rheology. We vary the temperature profile, the ratio of basalt vs. peridotite, the rheology of the mantle peridotites and the loading rate of the subduction interface. Over multiple earthquake cycles, our two end member experiments show that megathrust earthquakes are dominate the seismicity for cold condition (e.g. Japan trench) while both coseismic and aseismic events account for the seismicity for warm condition (e.g. Nankai trench).
-
Asteroid Shapes Are Always Close To Fluid Equilibrium
NASA Astrophysics Data System (ADS)
Tanga, Paolo; Comito, C.; Hestroffer, D.; Richardson, D. C.
2010-10-01
The simple evidence that asteroid are composed by solid rocks suggests that their shape can be rather far from the theoretical equilibrium for rotating fluid bodies. The possible fragmented ("rubble-pile") nature of most of them has suggested interpretations based on elasto-plastic models (such as the Mohr-Coulomb theory) that take into account the static behavior of a granular structure. However, these approaches did not incorporate explicitly the possible evolution of shapes in time due to external factors such as crater forming impacts or tidal deformations. We revisited the theory of equilibrium shapes for fluids, quantitatively evaluating - by appropriate metrics - the distance of the observed shapes from fluid equilibrium. This distance turns out to be much smaller than previously expected. On the basis of this evidence, we simulated numerically the evolution of gravitational aggregates having a small degree of internal friction, consistent with the theoretical findings. Our results offer a global scenario for the evolution of asteroid shapes under the action of gradual stresses due to minor impacts, tidal forces and seismic shaking. We show that actual asteroid shapes are consistent with the evolution of aggregates tending towards minimum free energy states. We are able to explain the samples of observed shapes obtained by different techniques. Our findings strongly support a highly porous and fragmented nature for most asteroids, at least in an external layer. Reference: Tanga et al. 2009: ApJ Letters, 706, 1, L197-L202 Acknowledgments: PT and CC have been supported by the "Programme Nationale de Planetologie" of France; DCR acknowledges support by the NASA (grant no. NNX08AM39G issued through the Office of Space Science) and by the NSF (grant no. AST0708110).
-
NASA Technical Reports Server (NTRS)
Hoffarth, C.; Khaled, B.; Rajan, S. D.; Goldberg, R.; Carney, K.; DuBois, P.; Blankenhorn, Gunther
2016-01-01
An orthotropic elasto-plastic-damage three-dimensional model with tabulated input has been developed to analyze the impact response of composite materials. The theory has been implemented as MAT 213 into a tailored version of LS-DYNA being developed under a joint effort of the FAA and NASA and has the following features: (a) the theory addresses any composite architecture that can be experimentally characterized as an orthotropic material and includes rate and temperature sensitivities, (b) the formulation is applicable for solid as well as shell element implementations and utilizes input data in a tabulated form directly from processed experimental data, (c) deformation and damage mechanics are both accounted for within the material model, (d) failure criteria are established that are functions of strain and damage parameters, and mesh size dependence is included, and (e) the theory can be efficiently implemented into a commercial code for both sequential and parallel executions. The salient features of the theory as implemented in LS-DYNA are illustrated using a widely used composite - the T800S/3900-2B[P2352W-19] BMS8-276 Rev-H-Unitape fiber/resin unidirectional composite. First, the experimental tests to characterize the deformation, damage and failure parameters in the material behavior are discussed. Second, the MAT213 input model and implementation details are presented with particular attention given to procedures that have been incorporated to ensure that the yield surfaces in the rate and temperature dependent plasticity model are convex. Finally, the paper concludes with a validation test designed to test the stability, accuracy and efficiency of the implemented model.
-
NASA Astrophysics Data System (ADS)
Bistacchi, A.; Pisterna, R.; Romano, V.; Rust, D.; Tibaldi, A.
2009-04-01
The plumbing system that connects a sub-volcanic magma reservoir to the surface has been the object of field characterization and mechanical modelling efforts since the pioneering work by Anderson (1936), who produced a detailed account of the spectacular Cullin Cone-sheet Complex (Isle of Skye, UK) and a geometrical and mechanical model aimed at defining the depth to the magma chamber. Since this work, the definition of the stress state in the half space comprised between the magma reservoir and the surface (modelled either as a flat surface or a surface comprising a volcanic edifice) was considered the key point in reconstructing dike propagation paths from the magma chamber. In fact, this process is generally seen as the propagation in an elastic media of purely tensional joints (mode I or opening mode propagation), which follow trajectories perpendicular to the least compressive principal stress axis. Later works generally used different continuum mechanics methodologies (analytic, BEM, FEM) to solve the problem of a pressure source (the magma chamber, either a point source or a finite volume) in an elastic (in some cases heterogeneous) half space (bounded by a flat topography or topped by a "volcano"). All these models (with a few limited exceptions) disregard the effect of the regional stress field, which is caused by tectonic boundary forces and gravitational body load, and consider only the pressure source represented by the magma chamber (review in Gudmundsson, 2006). However, this is only a (sometimes subordinate) component of the total stress field. Grosfils (2007) first introduced the gravitational load (but not tectonic stresses) in an elastic model solved with FEM in a 2D axisymmetric half-space, showing that "failure to incorporate gravitational loading correctly" affect the calculated stress pattern and many of the predictions that can be drawn from the models. In this contribution we report on modelling results that include: 2D axisymmetric or true 3D geometry; gravitational body load; anisotropic tectonic stresses; different shapes and depths of the magma chamber; different overpressure levels in the magma chamber; different shapes of the topographic surface (e.g. flat, volcano, caldera); linear-elastic or elasto-plastic Drucker-Prager rheology. The latter point, which in our opinion constitutes a fundamental improvement in the model, has proven necessary because in a purely elastic model the stress state would rise at levels that cannot be sustained by geologic materials. Particularly around and above the magma chamber, yielding is expected, influencing the stress field in the remaining modelling domain. The non-linear problem has been solved with the commercial finite element package Comsol Multiphysics, using a parametric solver. At the same time, a field structural analysis of the classical Cuillin Cone-sheet Complex has been performed. This analysis has shown that four distinct families of cone sheets of different age do exist. Among these, the sheets with the higher dip angle range (80-65°) are confirmed as purely tensional joints, but those with a lower dip angle range (60-40°) are quite often (when suitable markers are available) associated with a measurable shear component. Combining these new field observations with mechanical modelling results, we propose a new interpretation for the Cuillin Cone Sheet Complex. The plumbing system was composed by both purely tensional joints and mesoscopic faults with a shear component, produced in response to the regional stress field perturbed by the magma chamber, and later passively re-used as magma emplacement conduits. Under this assumption, the observed geometry of the Cuillin Cone-sheet Complex is consistent with a relatively shallow magma chamber with a flattened laccolite shape. The shape of the palaeotopography, now completely eroded, has also been considered, but is more weakly constrained by modelling results. References: Anderson E.M., 1936. The dynamics of the formation of cone-sheets, ring-dykes and cauldron subsidences. Proc R Soc Edinburgh, 56, 128-157. Grosfils E.B., 2007. Magma reservoir failure on the terrestrial planets: Assessing the importance of gravitational loading in simple elastic models. Journal of Volcanology and Geothermal Research, 166 (2), 47-75. Gudmundsson A. , 2006. How local stresses control magma-chamber ruptures, dyke injections, and eruptions in composite volcanoes. Earth Science Reviews, 79 (1), 1-31.
-
DOE Office of Scientific and Technical Information (OSTI.GOV)
Chattopadhyay, J.; Dutta, B.K.; Kushwaha, H.S.
Leak-Before-Break (LBB) is being used to design the primary heat transport piping system of 500 MWe Indian Pressurized Heavy Water Reactors (IPHWR). The work is categorized in three directions to demonstrate three levels of safety against sudden catastrophic break. Level 1 is inherent in the design procedure of piping system as per ASME Sec.III with a well defined factor of safety. Level 2 consists of fatigue crack growth study of a postulated part-through flaw at the inside surface of pipes. Level 3 is stability analysis of a postulated leakage size flaw under the maximum credible loading condition. Developmental work relatedmore » to demonstration of level 2 and level 3 confidence is described in this paper. In a case study on fatigue crack growth on PHT straight pipes for level 2, negligible crack growth is predicted for the life of the reactor. For level 3 analysis, the R6 method has been adopted. A database to evaluate SIF of elbows with throughwall flaws under combined internal pressure and bending moment has been generated to provide one of the inputs for R6 method. The methodology of safety assessment of elbow using R6 method has been demonstrated for a typical pump discharge elbow. In this analysis, limit load of the cracked elbow has been determined by carrying out elasto-plastic finite element analysis. The limit load results compared well with those given by Miller. However, it requires further study to give a general form of limit load solution. On the experimental front, a set of small diameter pipe fracture experiments have been carried out at room temperature and 300{degrees}C. Two important observations of the experiments are - appreciable drop in maximum load at 300{degrees}C in case of SS pipes and out-of-plane crack growth in case of CS pipes. Experimental load deflection curves are finally compared with five J-estimation schemes predictions. A material database of PHT piping materials is also being generated for use in LBB analysis.« less
-
NASA Astrophysics Data System (ADS)
Burczynski, Grzegorz; Marcinowski, Jakub
2014-09-01
The paper deals with the numerical modelling of a complex, steel shell structure. The part under analysis is the upper segment of a steel pylon, which consists of several cylindrical shells and one conical segment. Particular parts of the structure are welded together. Geometrical and loading data calculations were performed for the particular material for both an ideally elastic case and an elasto-plastic case. The conclusion that the structural member analysed required strengthening were drawn on the basis of these results. The structural modification was proposed and additional calculations for this modified structure were also performed. Introduced additional shell elements locked the mechanism of plastic flow. The proposed modification can be treated as a possible strengthening concept. The whole analysis was performed by means of the ABAQUS system but some stages of calculations were also verified by the COSMOS/M system. Przedmiotem pracy jest numeryczne modelowanie pewnej bardzo złożonej, stalowej konstrukcji powłokowej. Analizowana szczegółowo czesc jest górnym fragmentem stalowego pylonu, na który składa sie kilka odcinków powłok cylindrycznych oraz jeden segment stożkowy. Te poszczególne fragmenty konstrukcji były ze soba połaczone spawaniem. Dla znanych parametrów materiałowych, geometrycznych i obciażeniowych wykonano obliczenia w zakresie idealnie spreżystym oraz w zakresie spreżystoplastycznym. Na podstawie tych obliczen wyciagnieto wniosek o koniecznosci wzmocnienia tej czesci pylonu. Zaproponowano istotna modyfikacje istniejacej konstrukcji i wykonano dla niej ponownie obliczenia. Wprowadzone dodatkowe elementy powłokowe zablokowały mechanizm plastycznego płyniecia. Zaproponowana modyfikacje można potraktowac jako jedna z możliwych koncepcji wzmocnienia konstrukcji. Wszystkie analizy numeryczne zostały wykonane za pomoca systemu ABAQUS. Pewne wybrane fragmenty obliczen były weryfikowane także z pomoca systemu COSMOS/M.
-
A generalized model for stability of trees under impact conditions
NASA Astrophysics Data System (ADS)
Dattola, Giuseppe; Crosta, Giovanni; Castellanza, Riccardo; di Prisco, Claudio; Canepa, Davide
2016-04-01
Stability of trees to external actions involve the combined effects of stem and tree root systems. A block impacting on the stem or an applied force pulling the stem can cause a tree instability involving stem bending or failure and tree root rotation. So different contributions are involved in the stability of the system. The rockfalls are common natural phenomena that can be unpredictable in terms of frequency and magnitude characteristics, and this makes difficult the estimate of potential hazard and risk for human lives and activities. In mountain areas a natural form of protection from rockfalls is provided by forest growing. The difficulties in the assessment of the real capability of this natural barrier by means of models is an open problem. Nevertheless, a large amount of experimental data are now available which provides support for the development of advanced theoretical framework and corresponding models. The aim of this contribution consists in presenting a model developed to predict the behavior of trees during a block impact. This model describes the tree stem by means of a linear elastic beam system consisting of two beams connected in series and with an equivalent geometry. The tree root system is described via an equivalent foundation, whose behavior is modelled through an elasto-plastic macro-element model. In order to calibrate the model parameters, simulations reproducing a series of winching tests, are performed. These numerical simulations confirm the capability of the model to predict the mechanical behavior of the stem-root system in terms of displacement vs force curves. Finally, numerical simulations of the impact of a boulder with a tree stem are carried out. These simulations, done under dynamic regime and with the model parameters obtained from the previous set of simulations, confirm the capability of the model to reproduce the effects on the stem-roots system generated by impulsive loads.
-
Consequences of the presence of a weak fault on the stress and strain within an active margin
NASA Astrophysics Data System (ADS)
Conin, M.; Henry, P.; Godard, V.; Bourlange, S.
2009-12-01
Accreting margins often display an outer thrust and fold belt and an inner forearc domain overlying the subduction plate. Assuming that this overlying material behaves as Coulomb material, the outer wedge and the inner wedge are classically approximated as a critical state and a stable state Coulomb wedge, respectively. Critical Coulomb wedge theory can account for the transition from wedge to forearc. However, it cannot be used to determine the state of stress in the transition zone, nor the consequences of a discontinuity within the margin. The presence of a discontinuity such as a splay fault having a low effective friction coefficient should affect the stress state within the wedge, at least locally around the splay fault. Moreover, the effective friction coefficient of the seismogenic zone is expected to vary during the seismic cycle, and this may influence the stability of the Coulomb wedges. We use the ADELI finite element code (Chery and Hassani, 2000) to model the quasi-static stress and strain of a decollement and splay fault system, within a two dimensional elasto-plastic wedge with Drucker-Prager rheology. The subduction plane, the basal decollement of the accretionary wedge and the splay fault are modeled with contact elements. The modeled margin comprises an inner and an outer domain with distinct tapers and basal friction coefficients. For a given splay fault geometry, we evaluate the friction coefficient threshold for splay fault activation as a function of the basal friction coefficients, and examine the consequences of motion along the splay fault on stress and strain within the wedge and on the surface slope at equilibrium. Friction coefficients are varied in time to mimic the consequence of the seismic cycle on the static stress state and strain distribution. Results show the possibility of coexistence of localized extensional regime above the splay fault within a regional compressional regime. Such coexistence is consistent with stress orientation estimation made from breakouts in the Nankai accretionary prim (Kinoshita et al, 2009).
-
Effects of Elastoplastic Material Properties on Shallow Fault Slip and Surface Displacement Fields
NASA Astrophysics Data System (ADS)
Nevitt, J. M.; Brooks, B. A.; Minson, S. E.; Lockner, D. A.; Moore, D. E.; Ericksen, T. L.; Hudnut, K. W.; Glennie, C. L.; Madugo, C. M.
2016-12-01
A primary control on the geodynamics of rifting is the thermal regime. To better understand the geodynamics of rifting in the northern Gulf of California we systematically measured heat-flow across the Wagner Basin, a tectonically active basin that lies near the southern terminus of the Cerro Prieto fault. Seismic reflection profiles show sediment in excess of 5 s two-way travel time, implying a sediment thickness of > 7 km. The heat flow profile is 40 km long, has a nominal measurement spacing of 1 km, and is collocated with a seismic reflection profile. Heat flow measurements were made with a 6.5-m violin-bow probe. Most measurements are of good quality in that the probe fully penetrated sediments and measurements were stable enough to invert for heat flow and thermal properties. We have estimated corrections for environmental perturbations due to changes in bottom water temperature and sedimentation. The mean and standard deviation of heat flow across the western, central, and eastern parts of the basin are 220±60, 99±14, 1058±519 mW m-2, respectively. Corrections for sedimentation would increase measured heat flow across the central part of basin by 40 to 60%. We interpret the relatively high heat flow and large variability on the western and eastern flanks in terms of upward fluid flow at depth below the seafloor, whereas the lower and more consistent values across the central part of the basin are suggestive of conductive heat transfer. Based on an observed fault depth of 1.75 km we estimated the maximum Darcy velocities through the western and eastern flanks as 3 and 10 cm yr-1, respectively. Heat flow across the central basin is consistent with gabbroic underplating at a depth of 15 km and suggests that continental rupture here has not gone to completion.
-
Hysteretic behavior using the explicit material point method
NASA Astrophysics Data System (ADS)
Sofianos, Christos D.; Koumousis, Vlasis K.
2018-05-01
The material point method (MPM) is an advancement of particle in cell method, in which Lagrangian bodies are discretized by a number of material points that hold all the properties and the state of the material. All internal variables, stress, strain, velocity, etc., which specify the current state, and are required to advance the solution, are stored in the material points. A background grid is employed to solve the governing equations by interpolating the material point data to the grid. The derived momentum conservation equations are solved at the grid nodes and information is transferred back to the material points and the background grid is reset, ready to handle the next iteration. In this work, the standard explicit MPM is extended to account for smooth elastoplastic material behavior with mixed isotropic and kinematic hardening and stiffness and strength degradation. The strains are decomposed into an elastic and an inelastic part according to the strain decomposition rule. To account for the different phases during elastic loading or unloading and smoothening the transition from the elastic to inelastic regime, two Heaviside-type functions are introduced. These act as switches and incorporate the yield function and the hardening laws to control the whole cyclic behavior. A single expression is thus established for the plastic multiplier for the whole range of stresses. This overpasses the need for a piecewise approach and a demanding bookkeeping mechanism especially when multilinear models are concerned that account for stiffness and strength degradation. The final form of the constitutive stress rate-strain rate relation incorporates the tangent modulus of elasticity, which now includes the Heaviside functions and gathers all the governing behavior, facilitating considerably the simulation of nonlinear response in the MPM framework. Numerical results are presented that validate the proposed formulation in the context of the MPM in comparison with finite element method and experimental results.
-
Blanchard, Romane; Morin, Claire; Malandrino, Andrea; Vella, Alain; Sant, Zdenka; Hellmich, Christian
2016-09-01
While in clinical settings, bone mineral density measured by computed tomography (CT) remains the key indicator for bone fracture risk, there is an ongoing quest for more engineering mechanics-based approaches for safety analyses of the skeleton. This calls for determination of suitable material properties from respective CT data, where the traditional approach consists of regression analyses between attenuation-related grey values and mechanical properties. We here present a physics-oriented approach, considering that elasticity and strength of bone tissue originate from the material microstructure and the mechanical properties of its elementary components. Firstly, we reconstruct the linear relation between the clinically accessible grey values making up a CT, and the X-ray attenuation coefficients quantifying the intensity losses from which the image is actually reconstructed. Therefore, we combine X-ray attenuation averaging at different length scales and over different tissues, with recently identified 'universal' composition characteristics of the latter. This gives access to both the normally non-disclosed X-ray energy employed in the CT-device and to in vivo patient-specific and location-specific bone composition variables, such as voxel-specific mass density, as well as collagen and mineral contents. The latter feed an experimentally validated multiscale elastoplastic model based on the hierarchical organization of bone. Corresponding elasticity maps across the organ enter a finite element simulation of a typical load case, and the resulting stress states are increased in a proportional fashion, so as to check the safety against ultimate material failure. In the young patient investigated, even normal physiological loading is probable to already imply plastic events associated with the hydrated mineral crystals in the bone ultrastructure, while the safety factor against failure is still as high as five. Copyright © 2016 John Wiley & Sons, Ltd. Copyright © 2016 John Wiley & Sons, Ltd.
-
NASA Astrophysics Data System (ADS)
Hatami, M. K.; Pardoen, T.; Lacroix, G.; Berke, P.; Jacques, P. J.; Massart, T. J.
2017-01-01
TRansformation Induced Plasticity (TRIP) is a very effective mechanism to increase the strain hardening capacity of multiphase steels containing a fraction of metastable austenite, leading to both high strength and large uniform elongation. Excellent performances have been reached in the past 20 years, with recent renewed interest through the development of the 3rd generation of high strength steels often involving a TRIP effect. The microstructure and composition optimization is complex due to the interplay of coupled effects on the transformation kinetics and work hardening such as phase stability, size of retained austenite grains, temperature and loading path. In particular, recent studies have shown that the TRIP effect can only be quantitatively captured for realistic microstructures if strain gradient plasticity effects are taken into account, although direct experimental validation of this claim is missing. Here, an original computational averaging scheme is developed for predicting the elastoplastic response of TRIP aided multiphase steels based on a strain gradient plasticity model. The microstructure is represented by an aggregate of many elementary unit cells involving each a fraction of retained austenite with a specified stability. The model parameters, involving the transformation kinetics, are identified based on experimental tensile tests performed at different temperatures. The model is further assessed towards original experiments, involving temperature changes during deformation. A classical size independent plasticity model is shown unable to capture the TRIP effect on the mechanical response. Conversely, the strain gradient formulation properly predicts substantial variations of the strain hardening with deformation and temperature, hence of the uniform elongation in good agreement with the experiments. A parametric study is performed to get more insight on the effect of the material length scale as well as to determine optimum transformation kinetics to reach the highest possible strength-ductility balance. It is shown that the uniform elongation can potentially be increased by 50% or more, paving the way towards future microstructure engineering efforts.
-
Three-dimensional, thermo-mechanical and dynamical analogue experiments of subduction: first results
NASA Astrophysics Data System (ADS)
Boutelier, D.; Oncken, O.
2008-12-01
We present a new analogue modeling technique developed to investigate the mechanics of the subduction process and the build-up of subduction orogenies. The model consists of a tank filled with water representing the asthenosphere and two lithospheric plates made of temperature-sensitive hydrocarbon compositional systems. These materials possess elasto-plastic properties allowing the scaling of thermal and mechanical processes. A conductive thermal gradient is imposed in the lithosphere prior to deformation. The temperature of the asthenosphere and model surface are imposed and controlled with an electric heater, two infrared ceramic heat emitters, two thermocouples and a thermo-regulator. This system allows an unobstructed view of the model surface, which is monitored using a stereoscopic particle image technique. This monitoring technique provides a precise quantification of the horizontal deformation and variations of elevation in the three-dimensional model. Convergence is imposed with a piston moving at a constant rate or pushing at a constant stress. The velocity is scaled using the dimensionless ratio of thermal conduction over advection. The experiments are first produced at a constant rate and the stress in the horizontal direction of the convergence is recorded. Then the experiment is reproduced with a constant stress boundary condition where the stress value is set to the averaged value obtained in the previous experiment. Therefore, an initial velocity allowing proper scaling of heat exchanges is obtained, but deformation in the model and spatial variations of parameters such as density or friction coefficient can produce variations of plate convergence velocity. This in turn impacts the strength of the model lithosphere because it changes the model thermal structure. In the first presented experiments the model lithosphere is one layer and the plate boundary is linear. The effects of variations of the subducting plate thickness, density and the lubrication of the interface between the plates are investigated.
-
Estimation of hysteretic losses for MgB2 tapes under the operating conditions of a generator
NASA Astrophysics Data System (ADS)
Vargas-Llanos, Carlos Roberto; Zermeño, Víctor M. R.; Sanz, Santiago; Trillaud, Frederic; Grilli, Francesco
2016-03-01
Hysteretic losses in the MgB2 wound superconducting coils of a 550 kW synchronous hybrid scaled generator were estimated as part of the European project SUPRAPOWER led by the Spanish Fundación Tecnalia Research & Innovation. Particular interest was given to the losses caused by the magnetic flux ripples in the rotor coils originating from the conventional stator during nominal operation. To compute these losses, a 2D finite element analysis was conducted and Maxwell’s equations written in the H-formulation were solved considering the nonlinear material properties of the conductor materials. The modeled tapes are made of multiple MgB2 filaments embedded in a Ni matrix and soldered to a high purity copper strip and insulated with Dacron braid. Three geometrical models of single tape cross sections of decreasing complexity were studied: (1) the first model reproduced closely the actual cross section obtained from tape micrographs. (2) The second model was obtained from the computed elasto-plastic deformation of a round Ni wire. (3) The third model was based on a simplified cross section with the superconducting filaments bundled in a single elliptical bulky structure. The last geometry allowed the validation of the modeling technique by comparing numerical losses with results from well-established analytical expressions. Additionally, the following cases of filament transpositions of the multi-filamentary tape were studied: no transposition, partial and full transposition; thereby improving understanding of the relevance of the tape fabrication process on the magnitude of the determination of ac losses. Finally, choosing the right level of geometrical detail, the following operational regimes of the machine and its impact on individual superconducting tape losses in the rotor were studied: bias-dc current, ramping current under ramping background field and magnetic flux ripples under dc background current and field.
-
A generalized self-consistent polycrystal model for the yield strength of nanocrystalline materials
NASA Astrophysics Data System (ADS)
Jiang, B.; Weng, G. J.
2004-05-01
Inspired by recent molecular dynamic simulations of nanocrystalline solids, a generalized self-consistent polycrystal model is proposed to study the transition of yield strength of polycrystalline metals as the grain size decreases from the traditional coarse grain to the nanometer scale. These atomic simulations revealed that a significant portion of atoms resides in the grain boundaries and the plastic flow of the grain-boundary region is responsible for the unique characteristics displayed by such materials. The proposed model takes each oriented grain and its immediate grain boundary to form a pair, which in turn is embedded in the infinite effective medium with a property representing the orientational average of all these pairs. We make use of the linear comparison composite to determine the nonlinear behavior of the nanocrystalline polycrystal through the concept of secant moduli. To this end an auxiliary problem of Christensen and Lo (J. Mech. Phys. Solids 27 (1979) 315) superimposed on the eigenstrain field of Luo and Weng (Mech. Mater. 6 (1987) 347) is first considered, and then the nonlinear elastoplastic polycrystal problem is addressed. The plastic flow of each grain is calculated from its crystallographic slips, but the plastic behavior of the grain-boundary phase is modeled as that of an amorphous material. The calculated yield stress for Cu is found to follow the classic Hall-Petch relation initially, but as the gain size decreases it begins to depart from it. The yield strength eventually attains a maximum at a critical grain size and then the Hall-Petch slope turns negative in the nano-range. It is also found that, when the Hall-Petch relation is observed, the plastic behavior of the polycrystal is governed by crystallographic slips in the grains, but when the slope is negative it is governed by the grain boundaries. During the transition both grains and grain boundaries contribute competitively.
-
Splay fault slip in a subduction margin, a new model of evolution
NASA Astrophysics Data System (ADS)
Conin, Marianne; Henry, Pierre; Godard, Vincent; Bourlange, Sylvain
2012-08-01
In subduction zones, major thrusts called splay faults are thought to slip coseismically during large earthquakes affecting the main plate interface. We propose an analytical condition for the activation of a splay fault based on force balance calculations and suggest thrusting along the splay fault is generally conditioned by the growth of the accretionary wedge, or by the erosion of the hanging wall. In theory, normal slip on the splay fault may occur when the décollement has a very low friction coefficient seaward. Such a low friction also implies an unstable extensional state within the outer wedge. Finite element elasto-plastic calculations with a geometry based on the Nankai Kumano section were performed and confirm that this analytical condition is a valid approximation. Furthermore, localized extension at a shallow level in the splay hanging wall is observed in models for a wide range of friction coefficients (from ∼0 to the value of internal friction coefficient of the rock, here equals to 0.4). The timing of slip established for the splay fault branch drilled on Nankai Kumano transect suggests a phase of concurrent splay and accretionary wedge growth ≈2 Ma to ≈1.5 Ma, followed by a locking of the splay ≈1.3 Ma. Active extension is observed in the hanging wall. This evolution can be explained by the activation of a deeper and weaker décollement, followed by an interruption of accretion. Activation of a splay as a normal fault, as hypothesized in the case of the Tohoku 2011 earthquake, can be achieved only if the friction coefficient on the décollement drops to near zero. We conclude that the tectonic stress state largely determines long-term variations of tightly related splay fault and outer décollement activity and thus influences where and how coseismic rupture ends, but that occurrence of normal slip on a splay fault requires coseismic friction reduction.
-
NASA Astrophysics Data System (ADS)
Wu, Guocai
This study systematically explores the mechanical behavior, damage tolerance and durability of fiber metal laminates, a promising candidate materials system for next generation aerospace structures. The experimental results indicated that GLARE laminates exhibited a bilinear deformation behavior under static in-plane loading. Both an analytical constitutive model based on a modified classical lamination theory which incorporates the elasto-plastic behavior of aluminum alloy and a numerical simulation based on finite element modeling are used to predict the nonlinear stress-strain response and deformation behavior of GLARE laminates. The blunt notched strength of GLARE laminates increased with decreasing specimen width and decreasing hole diameter. The notched strength of GLARE laminates was evaluated based on a modified point stress criterion. A computer simulation based on finite element method was performed to study stress concentration and distribution around the notch and verify the analytical and experimental results of notched strength. Good agreement is obtained between the model predictions and experimental results. Experimental results also indicate that GLARE laminates exhibited superior impact properties to those of monolithic 2024-T3 aluminum alloy at low velocity impact loading. The GLARE 5-2/1 laminate with 0°/90°/90°/0° fiber configuration exhibits a better impact resistance than the GLARE 4-3/2 laminate with 0°/90°/0° fiber orientation. The characteristic impact energies, the damage area, and the permanent deflection of laminates are used to evaluate the impact damage resistance. The post-impact residual tensile strength under various damage states ranging from the plastic dent, barely visible impact damage (BVID), clearly visible impact damage (CVID) up to the complete perforation was also measured and compared. The post-impact fatigue behavior under various stress levels and impact damage states was extensively explored. The damage initiation and progression, failure modes and crack propagation under different loading conditions were investigated and identified with microscopy, SEM, X-ray radiography, and by chemically removing outer aluminum layers.
-
Podczeck, Fridrun; Newton, J Michael; Fromme, Paul
2014-12-30
Flat, round tablets may have a breaking ("score") line. Pharmacopoeial tablet breaking load tests are diametral in their design, and industrially used breaking load testers often have automatic tablet feeding systems, which position the tablets between the loading platens of the machine with the breaking lines in random orientation to the applied load. The aim of this work was to ascertain the influence of the position of the breaking line in a diametral compression test using finite element methodology (FEM) and to compare the theoretical results with practical findings using commercially produced bevel-edged, scored tablets. Breaking line test positions at an angle of 0°, 22.5°, 45°, 67.5° and 90° relative to the loading plane were studied. FEM results obtained for fully elastic and elasto-plastic tablets were fairly similar, but they highlighted large differences in stress distributions depending on the position of the breaking line. The stress values at failure were predicted to be similar for tablets tested at an angle of 45° or above, whereas at lower test angles the predicted breaking loads were up to three times larger. The stress distributions suggested that not all breaking line angles would result in clean tensile failure. Practical results, however, did not confirm the differences in the predicted breaking loads, but they confirmed differences in the way tablets broke. The results suggest that it is not advisable to convert breaking loads obtained on scored tablets into tablet tensile strength values, and comparisons between different tablets or batches should carefully consider the orientation of the breaking line with respect to the loading plane, as the failure mechanisms appear to vary. Copyright © 2014 Elsevier B.V. All rights reserved.
-
Elbow stress indices using finite element analysis
NASA Astrophysics Data System (ADS)
Yu, Lixin
Section III of the ASME Boiler and Pressure Vessel Code (the Code) specifies rules for the design of nuclear power plant components. NB-3600 of the Code presents a simplified design method using stress indices---Scalar Coefficients used the modify straight pipe stress equations so that they can be applied to elbows, tees and other piping components. The stress indices of piping components are allowed to be determined both analytically and experimentally. This study concentrates on the determination of B2 stress indices for elbow components using finite element analysis (FEA). First, the previous theoretical, numerical and experimental investigations on elbow behavior were comprehensively reviewed, as was the philosophy behind the use of stress indices. The areas of further research was defined. Then, a comprehensive investigation was carried out to determine how the finite element method should be used to correctly simulate an elbow's structural behavior. This investigation included choice of element type, convergence of mesh density, use of boundary restraint and a reconciliation study between FEA and laboratory experiments or other theoretical formulations in both elastic and elasto-plastic domain. Results from different computer programs were also compared. Reasonably good reconciliation was obtained. Appendix II of the Code describes the experimental method to determine B2 stress indices based on load-deflection curves. This procedure was used to compute the B2 stress indices for various loading modes on one particular elbow configuration. The B2 stress indices thus determined were found to be about half of the value calculated from the Code equation. Then the effect on B2 stress indices of those factors such as internal pressure and flange attachments were studied. Finally, the investigation was extended to other configurations of elbow components. A parametric study was conducted on different elbow sizes and schedules. Regression analysis was then used to obtain a modified coefficient and exponent for the Code equation used to calculate B2 index for elbows.
-
Effect of texture on rheological properties: the case of ɛ-Fe (Invited)
NASA Astrophysics Data System (ADS)
Merkel, S.; Gruson, M.; Tomé, C. N.; Nishiyama, N.; Wang, Y.
2009-12-01
Lattice preferred orientations (LPO) are known to affect the physical properties of materials. However, in most high pressure deformation experiments, LPO are ignored when interpreting the measured stress-strain curves. In addition, stress measurements in those experiments are complicated by the effect of plastic deformation on the measured lattice strains(1). Here, we present a new interpretation of the results obtained on hcp-iron at up to 19 GPa and 600 K in the deformation-DIA(2). In those experiments, five independent stress-strain curves were obtained on axial shortening with a ductile behavior of the sample for all. Stress were studied using results of monochromatic X-ray diffraction and the elastic theory of lattice strains(3). However, measured stresses were inconsistent with a change of behavior after 4% axial strain, particularly for strains measured on the 0002 line. We use elasto-plastic self consistent modeling(1) to show that this change of behavior is due to the evolution of LPO in the sample. With compression, 10-10 planes in hcp-iron align parallel to the compression direction and this affects the rheological behavior of the sample, which can not be summarized in a simple average law. We will also discuss the implication of those results for the extraction of polycrystalline rheological properties for materials with non-random lattice preferred orientations and how this could affect our understanding of the Earth deep interior. 1- S. Merkel, C.N. Tomé, H.-R Wenk, A modeling analysis of the influence of plasticity on high pressure deformation of hcp-Co, Phys. Rev. B, 79, 064110 (2009) 2- N. Nishiyama, Y. Wang, M. L. Rivers, S. R. Sutton, D. Cookson, Rheology of e-iron up to 19 GPa and 600 K in the D-DIA, Geophys. Res. Lett., 34, L23304 (2007) 3- A. K. Singh, C. Balasingh, H. K. Mao, R. J. Hemley, J. Shu, Analysis of lattice strains measured under non-hydrostatic pressure, J. Appl. Phys., 83, 7567-7575 (1998)
-
Process and application of shock compression by nanosecond pulses of frequency-doubled Nd:YAG laser
NASA Astrophysics Data System (ADS)
Sano, Yuji; Kimura, Motohiko; Mukai, Naruhiko; Yoda, Masaki; Obata, Minoru; Ogisu, Tatsuki
2000-02-01
The authors have developed a new process of laser-induced shock compression to introduce a residual compressive stress on material surface, which is effective for prevention of stress corrosion cracking (SCC) and enhancement of fatigue strength of metal materials. The process developed is unique and beneficial. It requires no pre-conditioning for the surface, whereas the conventional process requires that the so-called sacrificial layer is made to protect the surface from damage. The new process can be freely applied to water- immersed components, since it uses water-penetrable green light of a frequency-doubled Nd:YAG laser. The process developed has the potential to open up new high-power laser applications in manufacturing and maintenance technologies. The laser-induced shock compression process (LSP) can be used to improve a residual stress field from tensile to compressive. In order to understand the physics and optimize the process, the propagation of a shock wave generated by the impulse of laser irradiation and the dynamic response of the material were analyzed by time-dependent elasto-plastic calculations with a finite element program using laser-induced plasma pressure as an external load. The analysis shows that a permanent strain and a residual compressive stress remain after the passage of the shock wave with amplitude exceeding the yield strength of the material. A practical system materializing the LSP was designed, manufactured, and tested to confirm the applicability to core components of light water reactors (LWRs). The system accesses the target component and remotely irradiates laser pulses to the heat affected zone (HAZ) along weld lines. Various functional tests were conducted using a full-scale mockup facility, in which remote maintenance work in a reactor vessel could be simulated. The results showed that the system remotely accessed the target weld lines and successfully introduced a residual compressive stress. After sufficient training for operational personnel, the system was applied to the core shroud of an existing nuclear power plant.
-
DOE Office of Scientific and Technical Information (OSTI.GOV)
Kim, Jihoon; Moridis, George J.
We investigate coupled flow and geomechanics in gas production from extremely low permeability reservoirs such as tight and shale gas reservoirs, using dynamic porosity and permeability during numerical simulation. In particular, we take the intrinsic permeability as a step function of the status of material failure, and the permeability is updated every time step. We consider gas reservoirs with the vertical and horizontal primary fractures, employing the single and dynamic double porosity (dual continuum) models. We modify the multiple porosity constitutive relations for modeling the double porous continua for flow and geomechanics. The numerical results indicate that production of gasmore » causes redistribution of the effective stress fields, increasing the effective shear stress and resulting in plasticity. Shear failure occurs not only near the fracture tips but also away from the primary fractures, which indicates generation of secondary fractures. These secondary fractures increase the permeability significantly, and change the flow pattern, which in turn causes a change in distribution of geomechanical variables. From various numerical tests, we find that shear failure is enhanced by a large pressure drop at the production well, high Biot's coefficient, low frictional and dilation angles. Smaller spacing between the horizontal wells also contributes to faster secondary fracturing. When the dynamic double porosity model is used, we observe a faster evolution of the enhanced permeability areas than that obtained from the single porosity model, mainly due to a higher permeability of the fractures in the double porosity model. These complicated physics for stress sensitive reservoirs cannot properly be captured by the uncoupled or flow-only simulation, and thus tightly coupled flow and geomechanical models are highly recommended to accurately describe the reservoir behavior during gas production in tight and shale gas reservoirs and to smartly design production scenarios.« less
-
NASA Astrophysics Data System (ADS)
Cerfontaine, B.; Charlier, R.; Collin, F.; Taiebat, M.
2017-10-01
Old mines or caverns may be used as reservoirs for fuel/gas storage or in the context of large-scale energy storage. In the first case, oil or gas is stored on annual basis. In the second case pressure due to water or compressed air varies on a daily basis or even faster. In both cases a cyclic loading on the cavern's/mine's walls must be considered for the design. The complexity of rockwork geometries or coupling with water flow requires finite element modelling and then a suitable constitutive law for the rock behaviour modelling. This paper presents and validates the formulation of a new constitutive law able to represent the inherently cyclic behaviour of rocks at low confinement. The main features of the behaviour evidenced by experiments in the literature depict a progressive degradation and strain of the material with the number of cycles. A constitutive law based on a boundary surface concept is developed. It represents the brittle failure of the material as well as its progressive degradation. Kinematic hardening of the yield surface allows the modelling of cycles. Isotropic softening on the cohesion variable leads to the progressive degradation of the rock strength. A limit surface is introduced and has a lower opening than the bounding surface. This surface describes the peak strength of the material and allows the modelling of a brittle behaviour. In addition a fatigue limit is introduced such that no cohesion degradation occurs if the stress state lies inside this surface. The model is validated against three different rock materials and types of experiments. Parameters of the constitutive laws are calibrated against uniaxial tests on Lorano marble, triaxial test on a sandstone and damage-controlled test on Lac du Bonnet granite. The model is shown to reproduce correctly experimental results, especially the evolution of strain with number of cycles.
-
NASA Astrophysics Data System (ADS)
Delage, Pierre; Karakostas, Foivos; Dhemaied, Amine; Belmokhtar, Malik; Lognonné, Philippe; Golombek, Matt; De Laure, Emmanuel; Hurst, Ken; Dupla, Jean-Claude; Kedar, Sharon; Cui, Yu Jun; Banerdt, Bruce
2017-10-01
In support of the InSight mission in which two instruments (the SEIS seismometer and the HP3 heat flow probe) will interact directly with the regolith on the surface of Mars, a series of mechanical tests were conducted on three different regolith simulants to better understand the observations of the physical and mechanical parameters that will be derived from InSight. The mechanical data obtained were also compared to data on terrestrial sands. The density of the regolith strongly influences its mechanical properties, as determined from the data on terrestrial sands. The elastoplastic compression volume changes were investigated through oedometer tests that also provided estimates of possible changes in density with depth. The results of direct shear tests provided values of friction angles that were compared with that of a terrestrial sand, and an extrapolation to lower density provided a friction angle compatible with that estimated from previous observations on the surface of Mars. The importance of the contracting/dilating shear volume changes of sands on the dynamic penetration of the mole was determined, with penetration facilitated by the ˜1.3 Mg/m3 density estimated at the landing site. Seismic velocities, measured by means of piezoelectric bender elements in triaxial specimens submitted to various isotropic confining stresses, show the importance of the confining stress, with lesser influence of density changes under compression. A power law relation of velocity as a function of confining stress with an exponent of 0.3 was identified from the tests, allowing an estimate of the surface seismic velocity of 150 m/s. The effect on the seismic velocity of a 10% proportion of rock in the regolith was also studied. These data will be compared with in situ data measured by InSight after landing.
-
Analysis of rainfall-induced slope instability using a field of local factor of safety
Lu, Ning; Şener-Kaya, Başak; Wayllace, Alexandra; Godt, Jonathan W.
2012-01-01
Slope-stability analyses are mostly conducted by identifying or assuming a potential failure surface and assessing the factor of safety (FS) of that surface. This approach of assigning a single FS to a potentially unstable slope provides little insight on where the failure initiates or the ultimate geometry and location of a landslide rupture surface. We describe a method to quantify a scalar field of FS based on the concept of the Coulomb stress and the shift in the state of stress toward failure that results from rainfall infiltration. The FS at each point within a hillslope is called the local factor of safety (LFS) and is defined as the ratio of the Coulomb stress at the current state of stress to the Coulomb stress of the potential failure state under the Mohr-Coulomb criterion. Comparative assessment with limit-equilibrium and hybrid finite element limit-equilibrium methods show that the proposed LFS is consistent with these approaches and yields additional insight into the geometry and location of the potential failure surface and how instability may initiate and evolve with changes in pore water conditions. Quantitative assessments applying the new LFS field method to slopes under infiltration conditions demonstrate that the LFS has the potential to overcome several major limitations in the classical FS methodologies such as the shape of the failure surface and the inherent underestimation of slope instability. Comparison with infinite-slope methods, including a recent extension to variably saturated conditions, shows further enhancement in assessing shallow landslide occurrence using the LFS methodology. Although we use only a linear elastic solution for the state of stress with no post-failure analysis that require more sophisticated elastoplastic or other theories, the LFS provides a new means to quantify the potential instability zones in hillslopes under variably saturated conditions using stress-field based methods.
-
Kim, Jihoon; Moridis, George J.
2014-12-01
We investigate coupled flow and geomechanics in gas production from extremely low permeability reservoirs such as tight and shale gas reservoirs, using dynamic porosity and permeability during numerical simulation. In particular, we take the intrinsic permeability as a step function of the status of material failure, and the permeability is updated every time step. We consider gas reservoirs with the vertical and horizontal primary fractures, employing the single and dynamic double porosity (dual continuum) models. We modify the multiple porosity constitutive relations for modeling the double porous continua for flow and geomechanics. The numerical results indicate that production of gasmore » causes redistribution of the effective stress fields, increasing the effective shear stress and resulting in plasticity. Shear failure occurs not only near the fracture tips but also away from the primary fractures, which indicates generation of secondary fractures. These secondary fractures increase the permeability significantly, and change the flow pattern, which in turn causes a change in distribution of geomechanical variables. From various numerical tests, we find that shear failure is enhanced by a large pressure drop at the production well, high Biot's coefficient, low frictional and dilation angles. Smaller spacing between the horizontal wells also contributes to faster secondary fracturing. When the dynamic double porosity model is used, we observe a faster evolution of the enhanced permeability areas than that obtained from the single porosity model, mainly due to a higher permeability of the fractures in the double porosity model. These complicated physics for stress sensitive reservoirs cannot properly be captured by the uncoupled or flow-only simulation, and thus tightly coupled flow and geomechanical models are highly recommended to accurately describe the reservoir behavior during gas production in tight and shale gas reservoirs and to smartly design production scenarios.« less
-
Development of BEM for ceramic composites
NASA Technical Reports Server (NTRS)
Henry, D. P.; Banerjee, P. K.; Dargush, G. F.; Hopkins, D. A.; Goldberg, R. K.
1993-01-01
BEST-CMS (boundary element solution technology - composite modeling system) is an advanced engineering system for the micro-analysis of fiber composite structures. BEST-CMS is based upon the boundary element program BEST3D which was developed for NASA by Pratt and Whitney Aircraft and the State University of New York at Buffalo under contract NAS3-23697. BEST-CMS presently has the capabilities for elastostatic analysis, steady-state and transient heat transfer analysis, steady-state and transient concurrent thermoelastic analysis, and elastoplastic and creep analysis. The fibers are assumed to be perfectly bonded to the composite matrix, or in the case of static or steady-state analysis, the fibers may be assumed to have spring connections, thermal resistance, and/or frictional sliding between the fibers and the composite matrix. The primary objective of this user's manual is to provide an overview of all BEST-CMS capabilities, along with detailed descriptions of the input data requirements. In the next chapter, a brief review of the theoretical background is presented for each analysis category. Then, chapter three discusses the key aspects of the numerical implementation, while chapter four provides a tutorial for the beginning BEST-CMS user. The heart of the manual, however, is in chapter five, where a complete description of all data input items is provided. Within this chapter, the individual entries are grouped on a functional basis for a more coherent presentation. Chapter six includes sample problems and should be of considerable assistance to the novice. Chapter seven includes capsules of a number of fiber-composite analysis problems that have been solved using BEST-CMS. This chapter is primarily descriptive in nature and is intended merely to illustrate the level of analysis that is possible within the present BEST-CMS system. Chapter eight contains a detail description of the BEST-CMS Neutral File which is helpful in writing an interface between BEST-CMS and any graphic post-processor program. Finally, all pertinent references are listed in chapter nine.
-
Features of tuned mass damper behavior under strong earthquakes
NASA Astrophysics Data System (ADS)
Nesterova, Olga; Uzdin, Alexander; Fedorova, Maria
2018-05-01
Plastic deformations, cracks and destruction of structure members appear in the constructions under strong earthquakes. Therefore constructions are characterized by a nonlinear deformation diagram. Two types of construction non-linearity are considered in the paper. The first type of nonlinearity is elastoplastic one. In this case, plastic deformations occur in the structural elements, and when the element is unloaded, its properties restores. Among such diagrams are the Prandtl diagram, the Prandtl diagram with hardening, the Ramberg-Osgood diagram and others. For systems with such nonlinearity there is an amplitude-frequency characteristic and resonance oscillation frequencies. In this case one can pick up the most dangerous accelerograms for the construction. The second type of nonlinearity is nonlinearity with degrading rigidity and dependence of behavior on the general loading history. The Kirikov-Amankulov model is one of such ones. Its behavior depends on the maximum displacement in the stress history. Such systems do not have gain frequency characteristic and resonance frequency. The period of oscillation of such system is increasing during the system loading, and the system eigen frequency decreases to zero at the time of collapse. In the cases under consideration, when investigating the system with MD behavior, the authors proposed new efficiency criteria. These include the work of plastic deformation forces for the first type of nonlinearity, which determines the possibility of progressive collapse or low cycle fatigue of the structure members. The period of system oscillations and the time to collapse of the structural support members are the criterion for systems with degrading rigidity. In the case of non-linear system behavior, the efficiency of MD application decreases, because the fundamental structure period is reduced because of structure damages and the MD will be rebound from the blanking regime. However, the MD using can significantly reduce the damageability of the protected object.
-
Estimation of a melting probe's penetration velocity range to reach icy moons' subsurface ocean
NASA Astrophysics Data System (ADS)
Erokhina, Olga; Chumachenko, Eugene
2014-05-01
In modern space science one of the actual branches is icy satellites explorations. The main interest is concentrated around Jovian's moons Europa and Ganymede, Saturn's moons Titan and Enceladus that are covered by thick icy layer according to "Voyager1", "Voyager2", "Galileo" and "Cassini" missions. There is a big possibility that under icy shell could be a deep ocean. Also conditions on these satellites allow speculating about possible habitability, and considering these moons from an astrobiological point of view. One of the possible tasks of planned missions is a subsurface study. For this goal it is necessary to design special equipment that could be suitable for planetary application. One of the possible means is to use a melting probe which operates by melting and moves by gravitational force. Such a probe should be relatively small, should not weight too much and should require not too much energy. In terrestrial case such kind of probe has been successfully used for glaciers study. And it is possible to extrapolate the usage of such probe to extraterrestrial application. One of the tasks is to estimate melting probe's penetration velocity. Although there are other unsolved problems such as analyzing how the probe will move in low gravity and low atmospheric pressure; knowing whether hole will be closed or not when probe penetrate thick enough; and considering what order could be a penetration velocity. This study explores two techniques of melting probe's movement. One of them based on elasto-plastic theory and so-called "solid water" theory, and other one takes phase changing into account. These two techniques allow estimating melting probe's velocity range and study whole process. Based on these technique several cases of melting probe movement were considered, melting probe's velocity range estimated, influence of different factors studied and discussed and an easy way to optimize parameters of the melting probe proposed.
-
NASA Technical Reports Server (NTRS)
Goldberg, Robert K.; Carney, Kelly S.; Dubois, Paul; Hoffarth, Canio; Khaled, Bilal; Rajan, Subramaniam; Blankenhorn, Gunther
2016-01-01
A material model which incorporates several key capabilities which have been identified by the aerospace community as lacking in the composite impact models currently available in LS-DYNA(Registered Trademark) is under development. In particular, the material model, which is being implemented as MAT 213 into a tailored version of LS-DYNA being jointly developed by the FAA and NASA, incorporates both plasticity and damage within the material model, utilizes experimentally based tabulated input to define the evolution of plasticity and damage as opposed to specifying discrete input parameters (such as modulus and strength), and is able to analyze the response of composites composed with a variety of fiber architectures. The plasticity portion of the orthotropic, three-dimensional, macroscopic composite constitutive model is based on an extension of the Tsai-Wu composite failure model into a generalized yield function with a non-associative flow rule. The capability to account for the rate and temperature dependent deformation response of composites has also been incorporated into the material model. For the damage model, a strain equivalent formulation is utilized to allow for the uncoupling of the deformation and damage analyses. In the damage model, a diagonal damage tensor is defined to account for the directionally dependent variation of damage. However, in composites it has been found that loading in one direction can lead to damage in multiple coordinate directions. To account for this phenomena, the terms in the damage matrix are semi-coupled such that the damage in a particular coordinate direction is a function of the stresses and plastic strains in all of the coordinate directions. The onset of material failure, and thus element deletion, is being developed to be a function of the stresses and plastic strains in the various coordinate directions. Systematic procedures are being developed to generate the required input parameters based on the results of experimental tests.
-
Micromechanical Modeling of Woven Metal Matrix Composites
NASA Technical Reports Server (NTRS)
Bednarcyk, Brett A.; Pindera, Marek-Jerzy
1997-01-01
This report presents the results of an extensive micromechanical modeling effort for woven metal matrix composites. The model is employed to predict the mechanical response of 8-harness (8H) satin weave carbon/copper (C/Cu) composites. Experimental mechanical results for this novel high thermal conductivity material were recently reported by Bednarcyk et al. along with preliminary model results. The micromechanics model developed herein is based on an embedded approach. A micromechanics model for the local (micro-scale) behavior of the woven composite, the original method of cells (Aboudi), is embedded in a global (macro-scale) micromechanics model (the three-dimensional generalized method of cells (GMC-3D) (Aboudi). This approach allows representation of true repeating unit cells for woven metal matrix composites via GMC-3D, and representation of local effects, such as matrix plasticity, yarn porosity, and imperfect fiber-matrix bonding. In addition, the equations of GMC-3D were reformulated to significantly reduce the number of unknown quantities that characterize the deformation fields at the microlevel in order to make possible the analysis of actual microstructures of woven composites. The resulting micromechanical model (WCGMC) provides an intermediate level of geometric representation, versatility, and computational efficiency with respect to previous analytical and numerical models for woven composites, but surpasses all previous modeling work by allowing the mechanical response of a woven metal matrix composite, with an elastoplastic matrix, to be examined for the first time. WCGMC is employed to examine the effects of composite microstructure, porosity, residual stresses, and imperfect fiber-matrix bonding on the predicted mechanical response of 8H satin C/Cu. The previously reported experimental results are summarized, and the model predictions are compared to monotonic and cyclic tensile and shear test data. By considering appropriate levels of porosity, residual stresses, and imperfect fiber-matrix debonding, reasonably good qualitative and quantitative correlation is achieved between model and experiment.
-
Analytical modeling of structure-soil systems for lunar bases
NASA Technical Reports Server (NTRS)
Macari-Pasqualino, Jose Emir
1989-01-01
The study of the behavior of granular materials in a reduced gravity environment and under low effective stresses became a subject of great interest in the mid 1960's when NASA's Surveyor missions to the Moon began the first extraterrestrial investigation and it was found that Lunar soils exhibited properties quite unlike those on Earth. This subject gained interest during the years of the Apollo missions and more recently due to NASA's plans for future exploration and colonization of Moon and Mars. It has since been clear that a good understanding of the mechanical properties of granular materials under reduced gravity and at low effective stress levels is of paramount importance for the design and construction of surface and buried structures on these bodies. In order to achieve such an understanding it is desirable to develop a set of constitutive equations that describes the response of such materials as they are subjected to tractions and displacements. This presentation examines issues associated with conducting experiments on highly nonlinear granular materials under high and low effective stresses. The friction and dilatancy properties which affect the behavior of granular soils with low cohesion values are assessed. In order to simulate the highly nonlinear strength and stress-strain behavior of soils at low as well as high effective stresses, a versatile isotropic, pressure sensitive, third stress invariant dependent, cone-cap elasto-plastic constitutive model was proposed. The integration of the constitutive relations is performed via a fully implicit Backward Euler technique known as the Closest Point Projection Method. The model was implemented into a finite element code in order to study nonlinear boundary value problems associated with homogeneous as well as nonhomogeneous deformations at low as well as high effective stresses. The effect of gravity (self-weight) on the stress-strain-strength response of these materials is evaluated. The calibration of the model is performed via three techniques: (1) physical identification, (2) optimized calibration at the constitutive level, and (3) optimized calibration at the finite element level (Inverse Identification). Activities are summarized in graphic and outline form.
-
Mid-ocean ridge jumps associated with hotspot magmatism
NASA Astrophysics Data System (ADS)
Mittelstaedt, Eric; Ito, Garrett; Behn, Mark D.
2008-02-01
Hotspot-ridge interaction produces a wide range of phenomena including excess crustal thickness, geochemical anomalies, off-axis volcanic ridges and ridge relocations or jumps. Ridges are recorded to have jumped toward many hotspots including, Iceland, Discovery, Galápagos, Kerguelen and Tristan de Cuhna. The causes of ridge jumps likely involve a number of interacting processes related to hotspots. One such process is reheating of the lithosphere as magma penetrates it to feed near-axis volcanism. We study this effect by using the hybrid, finite-element code, FLAC, to simulate two-dimensional (2-D, cross-section) viscous mantle flow, elasto-plastic deformation of the lithosphere and heat transport in a ridge setting near an off-axis hotspot. Heating due to magma transport through the lithosphere is implemented within a hotspot region of fixed width. To determine the conditions necessary to initiate a ridge jump, we vary four parameters: hotspot magmatic heating rate, spreading rate, seafloor age at the location of the hotspot and ridge migration rate. Our results indicate that the hotspot magmatic heating rate required to initiate a ridge jump increases non-linearly with increasing spreading rate and seafloor age. Models predict that magmatic heating, itself, is most likely to cause jumps at slow spreading rates such as at the Mid-Atlantic Ridge on Iceland. In contrast, despite the higher magma flux at the Galápagos hotspot, magmatic heating alone is probably insufficient to induce a ridge jump at the present-day due to the intermediate ridge spreading rate of the Galápagos Spreading Center. The time required to achieve a ridge jump, for fixed or migrating ridges, is found to be on the order of 10 5-10 6 years. Simulations that incorporate ridge migration predict that after a ridge jump occurs the hotspot and ridge migrate together for time periods that increase with magma flux. Model results also suggest a mechanism for ridge reorganizations not related to hotspots such as ridge jumps in back-arc settings and ridge segment propagation along the Mid-Atlantic Ridge.
-
Eberhardsteiner, Lukas; Hellmich, Christian; Scheiner, Stefan
2012-01-01
Extracellular bone material can be characterised as a nanocomposite where, in a liquid environment, nanometre-sized hydroxyapatite crystals precipitate within as well as between long fibre-like collagen fibrils (with diameters in the 100 nm range), as evidenced from neutron diffraction and transmission electron microscopy. Accordingly, these crystals are referred to as ‘interfibrillar mineral’ and ‘extrafibrillar mineral’, respectively. From a topological viewpoint, it is probable that the mineralisations start on the surfaces of the collagen fibrils (‘mineral-encrusted fibrils’), from where the crystals grow both into the fibril and into the extrafibrillar space. Since the mineral concentration depends on the pore spaces within the fibrils and between the fibrils (there is more space between them), the majority of the crystals (but clearly not all of them) typically lie in the extrafibrillar space. There, larger crystal agglomerations or clusters, spanning tens to hundreds of nanometers, develop in the course of mineralisation, and the micromechanics community has identified the pivotal role, which this extrafibrillar mineral plays for tissue elasticity. In such extrafibrillar crystal agglomerates, single crystals are stuck together, their surfaces being covered with very thin water layers. Recently, the latter have caught our interest regarding strength properties (Fritsch et al. 2009 J Theor Biol. 260(2): 230–252) – we have identified these water layers as weak interfaces in the extrafibrillar mineral of bone. Rate-independent gliding effects of crystals along the aforementioned interfaces, once an elastic threshold is surpassed, can be related to overall elastoplastic material behaviour of the hierarchical material ‘bone’. Extending this idea, the present paper is devoted to viscous gliding along these interfaces, expressing itself, at the macroscale, in the well-known experimentally evidenced phenomenon of bone viscoelasticity. In this context, a multiscale homogenisation scheme is extended to viscoelasticity, mineral-cluster-specific creep parameters are identified from three-point bending tests on hydrated bone samples, and the model is validated by statistically and physically independent experiments on partially dried samples. We expect this model to be relevant when it comes to prediction of time-dependent phenomena, e.g. in the context of bone remodelling. PMID:22563708
-
Eberhardsteiner, Lukas; Hellmich, Christian; Scheiner, Stefan
2014-01-01
Extracellular bone material can be characterised as a nanocomposite where, in a liquid environment, nanometre-sized hydroxyapatite crystals precipitate within as well as between long fibre-like collagen fibrils (with diameters in the 100 nm range), as evidenced from neutron diffraction and transmission electron microscopy. Accordingly, these crystals are referred to as 'interfibrillar mineral' and 'extrafibrillar mineral', respectively. From a topological viewpoint, it is probable that the mineralisations start on the surfaces of the collagen fibrils ('mineral-encrusted fibrils'), from where the crystals grow both into the fibril and into the extrafibrillar space. Since the mineral concentration depends on the pore spaces within the fibrils and between the fibrils (there is more space between them), the majority of the crystals (but clearly not all of them) typically lie in the extrafibrillar space. There, larger crystal agglomerations or clusters, spanning tens to hundreds of nanometers, develop in the course of mineralisation, and the micromechanics community has identified the pivotal role, which this extrafibrillar mineral plays for tissue elasticity. In such extrafibrillar crystal agglomerates, single crystals are stuck together, their surfaces being covered with very thin water layers. Recently, the latter have caught our interest regarding strength properties (Fritsch et al. 2009 J Theor Biol. 260(2): 230-252) - we have identified these water layers as weak interfaces in the extrafibrillar mineral of bone. Rate-independent gliding effects of crystals along the aforementioned interfaces, once an elastic threshold is surpassed, can be related to overall elastoplastic material behaviour of the hierarchical material 'bone'. Extending this idea, the present paper is devoted to viscous gliding along these interfaces, expressing itself, at the macroscale, in the well-known experimentally evidenced phenomenon of bone viscoelasticity. In this context, a multiscale homogenisation scheme is extended to viscoelasticity, mineral-cluster-specific creep parameters are identified from three-point bending tests on hydrated bone samples, and the model is validated by statistically and physically independent experiments on partially dried samples. We expect this model to be relevant when it comes to prediction of time-dependent phenomena, e.g. in the context of bone remodelling.
-
Subsurface and Surface Characterization using an Information Framework Model
NASA Astrophysics Data System (ADS)
Samuel-Ojo, Olusola
Groundwater plays a critical dual role as a reservoir of fresh water for human consumption and as a cause of the most severe problems when dealing with construction works below the water table. This is why it is critical to monitor groundwater recharge, distribution, and discharge on a continuous basis. The conventional method of monitoring groundwater employs a network of sparsely distributed monitoring wells and it is laborious, expensive, and intrusive. The problem of sparse data and undersampling reduces the accuracy of sampled survey data giving rise to poor interpretation. This dissertation addresses this problem by investigating groundwater-deformation response in order to augment the conventional method. A blend of three research methods was employed, namely design science research, geological methods, and geophysical methods, to examine whether persistent scatterer interferometry, a remote sensing technique, might augment conventional groundwater monitoring. Observation data (including phase information for displacement deformation from permanent scatterer interferometric synthetic aperture radar and depth to groundwater data) was obtained from the Water District, Santa Clara Valley, California. An information framework model was built and applied, and then evaluated. Data was preprocessed and decomposed into five components or parts: trend, seasonality, low frequency, high frequency and octave bandwidth. Digital elevation models of observed and predicted hydraulic head were produced, illustrating the piezometric or potentiometric surface. The potentiometric surface characterizes the regional aquifer of the valley showing areal variation of rate of percolation, velocity and permeability, and completely defines flow direction, advising characteristics and design levels. The findings show a geologic forcing phenomenon which explains in part the long-term deformation behavior of the valley, characterized by poroelastic, viscoelastic, elastoplastic and inelastic deformations under the influence of an underlying geologic southward plate motion within the theory of plate tectonics. It also explains the impact of a history of heavy pumpage of groundwater during the agricultural and urbanization era. Thus the persistent scatterer interferometry method offers an attractive, non-intrusive, cost-effective augmentation of the conventional method of monitoring groundwater for water resource development and stability of soil mass.
-
The elastic and inelastic behavior of woven graphite fabric reinforced polyimide composites
NASA Astrophysics Data System (ADS)
Searles, Kevin H.
In many aerospace and conventional engineering applications, load-bearing composite structures are designed with the intent of being subjected to uniaxial stresses that are predominantly tensile or compressive. However, it is likely that biaxial and possibly triaxial states of stress will exist throughout the in-service life of the structure or component. The existing paradigm suggests that unidirectional tape materials are superior under uniaxial conditions since the vast majority of fibers lie in-plane and can be aligned to the loading axis. This may be true, but not without detriment to impact performance, interlaminar strength, strain to failure and complexity of part geometry. In circumstances where a sufficient balance of these properties is required, composites based on woven fabric reinforcements become attractive choices. In this thesis, the micro- and mesoscale elastic behavior of composites based on 8HS woven graphite fabric architectures and polyimide matrices is studied analytically and numerically. An analytical model is proposed to predict the composite elastic constants and is verified using numerical strain energy methods of equivalence. The model shows good agreement with the experiments and numerical strain energy equivalence. Lamina stresses generated numerically from in-plane shear loading show substantial shear and transverse normal stress concentrations in the transverse undulated tow which potentially leads to intralaminar damage. The macroscale inelastic behavior of the same composites is also studied experimentally and numerically. On an experimental basis, the biaxial and modified biaxial Iosipescu test methods are employed to study the weaker-mode shear and biaxial failure properties at room and elevated temperatures. On a numerical basis, the macroscale inelastic shear behavior of the composites is studied. Structural nonlinearities and material nonlinearities are identified and resolved. In terms of specimen-to-fixture interactions, load eccentricities, geometric (large strains and rotations) nonlinearities and boundary contact (friction) nonlinearities are explored. In terms of material nonlinearities, anisotropic plasticity and progressive damage are explored. A progressive damage criterion is proposed which accounts for the elastic strain energy densities in three directions. Of the types of nonlinearities studied, the nonlinear shear stress-strain behavior of the composites is principally from progressive intralaminar damage. Structural nonlinearities and elastoplastic deformation appear to be inconsequential.
-
NASA Astrophysics Data System (ADS)
Salameh, Christelle; Bard, Pierre-Yves; Guillier, Bertrand; Harb, Jacques; Cornou, Cécile; Gérard, Jocelyne; Almakari, Michelle
2017-04-01
Post-seismic investigations repeatedly indicate that structures having frequencies close to foundation soil frequencies exhibit significantly heavier damages (Caracas 1967; Mexico 1985; Pujili, Ecuador 1996; L'Aquila 2009). However, observations of modal frequencies of soils and buildings in a region or within a current seismic risk analysis are not fully considered together, even when past earthquakes have demonstrated that coinciding soil and building frequencies leads to greater damage. The present paper thus focuses on a comprehensive numerical analysis to investigate the effect of coincidence between site and building frequencies. A total of 887 realistic soil profiles are coupled with a set of 141 single-degree-of-freedom elastoplastic oscillators, and their combined (nonlinear) response is computed for both linear and nonlinear soil behaviors, for a large number (60) of synthetic input signals with various PGA levels and frequency contents. The associated damage is quantified on the basis of the maximum displacement as compared to both yield and ultimate post-elastic displacements, according to the RISK-UE project recommendations (Lagomarsino and Giovinazzi in Bull Earthq Eng 4(4):415-443, 2006), and compared with the damage obtained in the case of a similar building located on rock. The correlation between this soil/rock damage increment and a number of simplified mechanical and loading parameters is then analyzed using a neural network approach. The results emphasize the key role played by the building/soil frequency ratio even when both soil and building behave nonlinearly; other important parameters are the PGA level, the soil/rock velocity contrast and the building ductility. A numerical investigation based on simulation of ambient noise for the whole set of 887 profiles also indicates that the amplitude of H/ V ratio may be considered as a satisfactory proxy for site amplification when applied to measurements at urban scale. A very easy implementation of this method, using ambient vibration measurements both at ground level and within buildings, is illustrated with an example application for the city of Beirut (Lebanon).[Figure not available: see fulltext.
-
NASA Astrophysics Data System (ADS)
Liu, Haitao
The objective of the present study is to investigate damage mechanisms and thermal residual stresses of composites, and to establish the frameworks to model the particle-reinforced metal matrix composites with particle-matrix interfacial debonding, particle cracking or thermal residual stresses. An evolutionary interfacial debonding model is proposed for the composites with spheroidal particles. The construction of the equivalent stiffness is based on the fact that when debonding occurs in a certain direction, the load-transfer ability will lose in that direction. By using this equivalent method, the interfacial debonding problem can be converted into a composite problem with perfectly bonded inclusions. Considering the interfacial debonding is a progressive process in which the debonding area increases in proportion to external loading, a progressive interfacial debonding model is proposed. In this model, the relation between external loading and the debonding area is established using a normal stress controlled debonding criterion. Furthermore, an equivalent orthotropic stiffness tensor is constructed based on the debonding areas. This model is able to study the composites with randomly distributed spherical particles. The double-inclusion theory is recalled to model the particle cracking problems. Cracks inside particles are treated as penny-shape particles with zero stiffness. The disturbed stress field due to the existence of a double-inclusion is expressed explicitly. Finally, a thermal mismatch eigenstrain is introduced to simulate the inconsistent expansions of the matrix and the particles due to the difference of the coefficients of thermal expansion. Micromechanical stress and strain fields are calculated due to the combination of applied external loads and the prescribed thermal mismatch eigenstrains. For all of the above models, ensemble-volume averaging procedures are employed to derive the effective yield function of the composites. Numerical simulations are performed to analyze the effects of various parameters and several good agreements between our model's predictions and experimental results are obtained. It should be mentioned that all of expressions in the frameworks are explicitly derived and these analytical results are easy to be adopted in other related investigations.
-
NASA Technical Reports Server (NTRS)
Banerjee, P. K.; Henry, D. P.; Hopkins, D. A.; Goldberg, R. K.
1997-01-01
BEST-CMS (Boundary Element Solution Technology - Composite Modeling System) is an advanced engineering system for the micro-analysis of fiber composite structures. BEST-CMS is based upon the boundary element program BEST3D which was developed for NASA by Pratt and Whitney Aircraft and the State University of New York at Buffalo under contract NAS3-23697. BEST-CMS presently has the capabilities for elastostatic analysis, steady-state and transient heat transfer analysis, steady-state and transient concurrent thermoelastic analysis and elastoplastic and creep analysis. The fibers are assumed to be perfectly bonded to the composite matrix, or in the case of static or steady-state analysis, the fibers may be assumed to have spring connections, thermal resistance, and/or frictional sliding between the fibers and the composite matrix. The primary objective of this User's Manual is to provide an overview of all BEST-CMS capabilities, along with detailed descriptions of the input data requirements. A brief review of the theoretical background is presented for each analysis category. Then, Chapter 3 discusses the key aspects of the numerical implementation, while Chapter 4 provides a tutorial for the beginning BEST-CMS user. The heart of the manual, however, is in Chapter 5, where a complete description of all data input items is provided. Within this chapter, the individual entries are grouped on a functional basis for a more coherent presentation. Chapter 6 includes sample problems and should be of considerable assistance to the novice. Chapter 7 includes capsules of a number of fiber-composite analysis problems that have been solved using BEST-CMS. This chapter is primarily descriptive in nature and is intended merely to illustrate the level of analysis that is possible within the present BEST-CMS system. Chapter 8 contains a detailed description of the BEST-CMS Neutral File which is helpful in writing an interface between BEST- CMS and any graphic post-processor program. Finally, all pertinent references are listed in Chapter 9.
-
NASA Astrophysics Data System (ADS)
Kotowski, A. J.; Behr, W. M.; Tong, X.; Lavier, L.
2017-12-01
The rheology of the deep subduction interface strongly influences the occurrence, recurrence, and migration of episodic tremor and slow slip (ETS) events. To better understand the environment of deep ETS, we characterize the length scales and types of rheological heterogeneities that decorate the deep interface using an exhumed subduction complex. The Cycladic Blueschist Unit on Syros, Greece, records Eocene subduction to 60 km, partial exhumation along the top of the slab, and final exhumation along Miocene detachment faults. The CBU reached 450-580˚C and 14-16 kbar, PT conditions similar to where ETS occurs in several modern subduction zones. Rheological heterogeneity is preserved in a range of rock types on Syros, with the most prominent type being brittle pods embedded within a viscous matrix. Prograde, blueschist-facies metabasalts show strong deformation fabrics characteristic of viscous flow; cm- to m-scale eclogitic lenses are embedded within them as massive, veined pods, foliated pods rotated with respect to the blueschist fabric, and attenuated, foliation-parallel lenses. Similar relationships are observed in blueschist-facies metasediments interpreted to have deformed during early exhumation. In these rocks, metabasalts form lenses ranging in size from m- to 10s of m and are distributed at the m-scale throughout the metasedimentary matrix. Several of the metamafic lenses, and the matrix rocks immediately adjacent to them, preserve multiple generations of dilational veins and shear fractures filled with quartz and high pressure minerals. These observations suggest that coupled brittle-viscous deformation under high fluid pressures may characterize the subduction interface in the deep tremor source region. To test this further, we modeled the behavior of an elasto-plastic pod in a viscous shear zone under high fluid pressures. Our models show that local stress concentrations around the pod are large enough to generate transient dilational shear at seismic strain rates. Scaling the model up to a typical source area for deep tremor suggests these heterogeneities may yield a seismic moment similar to those calculated for tremor bursts in modern subduction zones.
-
Asymptotic analysis of stability for prismatic solids under axial loads
NASA Astrophysics Data System (ADS)
Scherzinger, W.; Triantafyllidis, N.
1998-06-01
This work addresses the stability of axially loaded prismatic beams with any simply connected crosssection. The solids obey a general class of rate-independent constitutive laws, and can sustain finite strains in either compression or tension. The proposed method is based on multiple scale asymptotic analysis, and starts with the full Lagrangian formulation for the three-dimensional stability problem, where the boundary conditions are chosen to avoid the formation of boundary layers. The calculations proceed by taking the limit of the beam's slenderness parameter, ɛ (ɛ 2 ≡ area/length 2), going to zero, thus resulting in asymptotic expressions for the critical loads and modes. The analysis presents a consistent and unified treatment for both compressive (buckling) and tensile (necking) instabilities, and is carried out explicitly up to o( ɛ4) in each case. The present method circumvents the standard structural mechanics approach for the stability problem of beams which requires the choice of displacement and stress field approximations in order to construct a nonlinear beam theory. Moreover, this work provides a consistent way to calculate the effect of the beam's slenderness on the critical load and mode to any order of accuracy required. In contrast, engineering theories give accurately the lowest order terms ( O( ɛ2)—Euler load—in compression or O(1)—maximum load—in tension) but give only approximately the next higher order terms, with the exception of simple section geometries where exact stability results are available. The proposed method is used to calculate the critical loads and eigenmodes for bars of several different cross-sections (circular, square, cruciform and L-shaped). Elastic beams are considered in compression and elastoplastic beams are considered in tension. The O( ɛ2) and O( ɛ4) asymptotic results are compared to the exact finite element calculations for the corresponding three-dimensional prismatic solids. The O( ɛ4) results give significant improvement over the O( ɛ2) results, even for extremely stubby beams, and in particular for the case of cross-sections with commensurate dimensions.
-
NASA Astrophysics Data System (ADS)
Tabari, Mehdi Sherveen Ghofrani
Although true-triaxial test (TTT) of rocks is now more extensive worldwide, stress-induced heterogeneity is not accounted for and usually simplified anisotropic models are used. Data from a TTT on a cubic sample of Fontainebleau sandstone is used in this study to evaluate our velocity imaging methodology. An anisotropic P wave velocity tomography method was developed using a geometrical approach based on an ellipsoidal P wavefront surface. During the two non-damaging phases of the experiment, saturation of the rock sample with water resulted in inaccurate tomographic images; however, during the final elasto-plastic phase of the experiment comprising major AE activities, tomographic images demonstrated reasonable anomalies. Thus, the P-S1-S2 velocity survey was utilized to obtain an accurate and reliable velocity image of the sample during the two non-damaging phases. This was accomplished using a numerical investigation by FLAC3D on the non-uniform distribution of stress over the sample to estimate the compaction pseudo-boundary surfaces within the rock. Thus, the problem of breakdown in the expected symmetry of shear wave velocities was resolved. It was discovered that a homogeneous anisotropic core in the center of the sample is formed under the standard polyaxial setup where elastic parameters could be computed. Off-diagonal elastic tensor parameters were obtained by a combination of various velocity survey data and justified the ellipsoidal model as being the most appropriate and facilitated the calculation of Thomsen parameters. The ellipsoidal heterogeneous velocity model was also verified by AE event location of transducer shots through the cubic rock specimen especially at the final phase of the experiment consisting lower-velocity zones bearing partially saturated fractures. AE of the rock during the whole experiment recorded by the surrounding transducers were investigated by location methods developed for anisotropic heterogeneous medium. AE events occurred in the vicinity of the dilation pseudo-boundaries where, a relatively large velocity gradient was formed and along parallel fractures in the sigma1/sigma2 plane. This research facilitated the computation of anisotropic parameters for rock during polyaxial tests contributing to enhanced AE interpretation of fracture growth processes in the rock under laboratory true-triaxial stress conditions.
-
Experimental Techniques Verified for Determining Yield and Flow Surfaces
NASA Technical Reports Server (NTRS)
Lerch, Brad A.; Ellis, Rod; Lissenden, Cliff J.
1998-01-01
Structural components in aircraft engines are subjected to multiaxial loads when in service. For such components, life prediction methodologies are dependent on the accuracy of the constitutive models that determine the elastic and inelastic portions of a loading cycle. A threshold surface (such as a yield surface) is customarily used to differentiate between reversible and irreversible flow. For elastoplastic materials, a yield surface can be used to delimit the elastic region in a given stress space. The concept of a yield surface is central to the mathematical formulation of a classical plasticity theory, but at elevated temperatures, material response can be highly time dependent. Thus, viscoplastic theories have been developed to account for this time dependency. Since the key to many of these theories is experimental validation, the objective of this work (refs. 1 and 2) at the NASA Lewis Research Center was to verify that current laboratory techniques and equipment are sufficient to determine flow surfaces at elevated temperatures. By probing many times in the axial-torsional stress space, we could define the yield and flow surfaces. A small offset definition of yield (10 me) was used to delineate the boundary between reversible and irreversible behavior so that the material state remained essentially unchanged and multiple probes could be done on the same specimen. The strain was measured with an off-the-shelf multiaxial extensometer that could measure the axial and torsional strains over a wide range of temperatures. The accuracy and resolution of this extensometer was verified by comparing its data with strain gauge data at room temperature. The extensometer was found to have sufficient resolution for these experiments. In addition, the amount of crosstalk (i.e., the accumulation of apparent strain in one direction when strain in the other direction is applied) was found to be negligible. Tubular specimens were induction heated to determine the flow surfaces at elevated temperatures. The heating system induced a large amount of noise in the data. By reducing thermal fluctuations and using appropriate data averaging schemes, we could render the noise inconsequential. Thus, accurate and reproducible flow surfaces (see the figure) could be obtained.
-
NASA Astrophysics Data System (ADS)
Leith, Kerry; Kupp, Jan; Geisenhof, Benedikt; Krautblatter, Michael
2015-04-01
Bedrock stresses in alpine regions result from the combined effects of exhumation, tectonics, topography, inelastic strain (e.g. fault displacement and fracture formation), and external loading. Gravitational loading by glacial ice can significantly affect near-surface stress magnitudes, although the nature of this effect and it's impact on stress distributions and bedrock fracturing is strongly dependent on the stress history of the bedrock landscape. We assess the effects of recent (post-Little Ice Age , ~1850 AD) and future deglaciation on bedrock stresses in the region of the Zugspitzplatt, a glaciated plateau surrounded by 1500 m high bedrock walls in SE Germany. We address this by undertaking a 2-D elasto-plastic finite element method analysis of stress changes and fracture propagation due to repeated glacial - interglacial cycles. Our model is initialised with upper crustal stresses in equilibrium with bedrock strength and regional tectonics, and we then simulate two cycles of major Pleistocene glaciation and deglaciation in order to dissipate stress concentrations and incorporate path-dependent effects of glacial loading on the landscape. We then simulate a final glacial cycle, and remove 1 m of bedrock to approximate glacial erosion across the topography. Finally, ice levels are reduced in accordance with known late-glacial and recent ice retreat, allowing us to compare relative stress changes and predicted patterns of fracture propagation to observed fracture distributions on the Zugspitzplatt. Model results compare favourably to observed fracture patterns, and indicate the plateau is likely to be undergoing N-S extension as a result of deglaciation, with a strong reduction of horizontal stress magnitudes beneath the present-day Schneeferner glacier. As each glacial cycle has a similar effect on the plateau, it is likely that surficial stresses are slightly tensile, and each cycle of deglaciation produces additional sub-vertical tensile fractures, which are then exploited by the karst groundwater system. Here we show how stress histories and brittle deformation in near-surface stress models can provide a better understanding of long-term rock slope evolution and failure as well as karst co-evolution in Alpine Environments.
-
Fracture propagation in Indiana Limestone interpreted via linear softening cohesive fracture model
NASA Astrophysics Data System (ADS)
Rinehart, Alex J.; Bishop, Joseph E.; Dewers, Thomas
2015-04-01
We examine the use of a linear softening cohesive fracture model (LCFM) to predict single-trace fracture growth in short-rod (SR) and notched 3-point-bend (N3PB) test configurations in Indiana Limestone. The broad goal of this work is to (a) understand the underlying assumptions of LCFM and (b) use experimental similarities and deviations from the LCFM to understand the role of loading paths of tensile fracture propagation. Cohesive fracture models are being applied in prediction of structural and subsurface fracture propagation in geomaterials. They lump the inelastic processes occurring during fracture propagation into a thin zone between elastic subdomains. LCFM assumes that the cohesive zone initially deforms elastically to a maximum tensile stress (σmax) and then softens linearly from the crack opening width at σmax to zero stress at a critical crack opening width w1. Using commercial finite element software, we developed LCFMs for the SR and N3PB configurations. After fixing σmax with results from cylinder splitting tests and finding an initial Young's modulus (E) with unconfined compressive strength tests, we manually calibrate E and w1 in the SR model against an envelope of experimental data. We apply the calibrated LCFM parameters in the N3PB geometry and compare the model against an envelope of N3PB experiments. For accurate simulation of fracture propagation, simulated off-crack stresses are high enough to require inclusion of damage. Different elastic moduli are needed in tension and compression. We hypothesize that the timing and location of shear versus extensional micromechanical failures control the qualitative macroscopic force-versus-displacement response in different tests. For accurate prediction, the LCFM requires a constant style of failure, which the SR configuration maintains until very late in deformation. The N3PB configuration does not maintain this constancy. To be broadly applicable between geometries and failure styles, the LCFM would require additional physics, possibly including elastoplastic damage in the bulk material and more complicated cohesive softening models.
-
NASA Astrophysics Data System (ADS)
Lavier, L. L.; Bennett, R. A.; Anderson, M. L.; Matti, J. C.
2005-05-01
Recent displacement rate and geodetic data on the San Andreas, San Jacinto and eastern California shear zone suggest that changes in the geometry and/or the magnitude of the applied forces on the crust (e.g., a general or local change in fault strike relative to plate motion) can generate strain repartitioning within the crust on time scales of millions to thousands of years. The rates over which this repartitioning takes place in response to changing forces are controlled by the rheological evolution of the lithosphere. We investigate the implications of observed fault displacement histories for the rheology of the lithosphere using 2.5 D numerical experiments of deformation in an analogue system. The numerical technique used allows for the spontaneous formation of elastoplastic shear zones and flow in a Maxwell viscoelastic lower crust. The results show that when a strike slip fault is rotated to strike obliquely to the direction of relative plate motion it causes changes in bending and frictional stresses due to the formation of topography. To accommodate these changes, a conjugate system of oblique-striking strike slip faults develops. The total displacement is then slowly distributed over the new fault system on the time scale of mountain building (i.e. million of years). The rate of change is dependent on the strength of the lithosphere as well as the amount of obliquity applied on the initial strike-slip fault. In other numerical experiments we show that in a system of multiple strike-slip fault zones, displacement rate changes can occur over a time scale of about 100 kyr. This time scale corresponds to the Maxwell time at the brittle ductile transition (BDT). In such a system the lithospheric displacement is alternatively distributed (over 100 kyr) in clusters localized in lower crustal channels and over strike-slip fault zones. We show that the clustering time scale is controlled by the ratio of upper to lower crustal strength. This incomplete exercise shows how displacement rates data sets spanning thousands to millions of years can be used to constrain numerical experiments of lithospheric deformation and, in doing so, place new constraints on the rheology of the lithosphere.
-
Mechanics and Physics of Solids, Uncertainy, and the Archetype-Genome Exemplar
NASA Astrophysics Data System (ADS)
Greene, M. Steven
This dissertation argues that the mechanics and physics of solids rely on a fundamental exemplar: the apparent properties of a system depend on the building blocks that comprise it. Building blocks are referred to as archetypes and apparent system properties as the system genome. Three entities are of importance: the archetype properties, the conformation of archetypes, and the properties of interactions activated by that conformation. The combination of these entities into the system genome is called assembly. To show the utility of the archetype-genome exemplar, the dissertation presents the mathematical construction and computational implementation of a new theory for solid mechanics that is a continuum manifestation of the assembly process. The so-called archetype-blending continuum theory aligns the form of globally valid balance laws with physics evolving in a material's composite constitutive response so that, by rethinking conventional micromechanics, the theory accounts naturally for each piece of the genome assembly triplet: archetypes, interactions, and their conformation. With the pieces of the triplet isolated in the theory, materials genome design concepts that separately control microstructure and property may be gleaned from exploration of the constitutive parameter space. A suite of simulations that apply the new theory to polymer nanocomposite materials demonstrate the ability of the theory to predict a robust material genome that includes damping properties, modulus weakening, local strain amplification, and size effects. The dissertation also presents a theoretical assessment of the importance of uncertainty propagation in the archetype-genome exemplar. The findings from a set of computational experiments on instances of a general class of microstructured materials suggest that when overlap occurs between the size of the system geometry and the features of the conformation, material genomes become less certain. Increasing nonuniformity of boundary conditions and the size of random field correlation lengths exacerbate this conclusion. These criteria are combined into a scalar metric used to assess the impact of archetype-level uncertainties on the material genome for general scenarios in solid mechanics. Exemplary benchmark problems include bending in elastoplasticity and instability-induced pattern transition in porous elastomer. The contributions of this dissertation are threefold: (1) the mathematical construction of a new continuum theory for mechanics and physics of solids, (2) implementation of the theory, and (3) theoretical assessment of the scenarios in which material genomes deviate from determinism.
-
DOE Office of Scientific and Technical Information (OSTI.GOV)
Vilarrasa, Víctor; Rutqvist, Jonny; Blanco Martin, Laura
Expansive soils are suitable as backfill and buffer materials in engineered barrier systems to isolate heat-generating nuclear waste in deep geological formations. The canisters containing nuclear waste would be placed in tunnels excavated at a depth of several hundred meters. The expansive soil should provide enough swelling capacity to support the tunnel walls, thereby reducing the impact of the excavation-damaged zone on the long-term mechanical and flow-barrier performance. In addition to their swelling capacity, expansive soils are characterized by accumulating irreversible strain on suction cycles and by effects of microstructural swelling on water permeability that for backfill or buffer materialsmore » can significantly delay the time it takes to reach full saturation. In order to simulate these characteristics of expansive soils, a dual-structure constitutive model that includes two porosity levels is necessary. The authors present the formulation of a dual-structure model and describe its implementation into a coupled fluid flow and geomechanical numerical simulator. The authors use the Barcelona Basic Model (BBM), which is an elastoplastic constitutive model for unsaturated soils, to model the macrostructure, and it is assumed that the strains of the microstructure, which are volumetric and elastic, induce plastic strain to the macrostructure. The authors tested and demonstrated the capabilities of the implemented dual-structure model by modeling and reproducing observed behavior in two laboratory tests of expansive clay. As observed in the experiments, the simulations yielded nonreversible strain accumulation with suction cycles and a decreasing swelling capacity with increasing confining stress. Finally, the authors modeled, for the first time using a dual-structure model, the long-term (100,000 years) performance of a generic heat-generating nuclear waste repository with waste emplacement in horizontal tunnels backfilled with expansive clay and hosted in a clay rock formation. The thermo-hydro-mechanical results of the dual-structure model were compared with those of the standard single-structure BBM. The main difference between the simulation results from the two models is that the dual-structure model predicted a time to fully saturate the expansive clay barrier on the order of thousands of years, whereas the standard single-structure BBM yielded a time on the order of tens of years. These examples show that a dual-structure model, such as the one presented here, is necessary to properly model the thermo-hydro-mechanical behavior of expansive soils.« less
-
Generalization of the slip line field theory for temperature sensitive visco-plastic materials
NASA Astrophysics Data System (ADS)
Paesold, Martin; Peters, Max; Regenauer-Lieb, Klaus; Veveakis, Manolis; Bassom, Andrew
2015-04-01
Geological processes can be a combination of various effects such as heat production or consumption, chemical reactions or fluid flow. These individual effects are coupled to each other via feedbacks and the mathematical analysis becomes challenging due to these interdependencies. Here, we concentrate solely on thermo-mechanical coupling and a main result of this work is that the coupling can depend on material parameters and boundary conditions and the coupling is more or less pronounced depending on theses parameters. The transitions from weak to strong coupling can be studied in the context of a bifurcation analysis. classically, Material instabilities in solids are approached as material bifurcations of a rate-independent, isothermal, elasto-plastic solid. However, previous research has shown that temperature and deformation rate are important factors and are fully coupled with the mechanical deformation. Early experiments in steel revealed a distinct pattern of localized heat dissipation and plastic deformation known as heat lines. Further, earth materials, soils, rocks and ceramics are known to be greatly influenced by temperature with strain localization being strongly affected by thermal loading. In this work, we provide a theoretical framework for the evolution of plastic deformation for such coupled systems, with a two-pronged approach to the prediction of localized failure. First, slip line field theory is employed to predict the geometry of the failure patterns and second, failure criteria are derived from an energy bifurcation analysis. The bifurcation analysis is concerned with the local energy balance of a material and compares the effects of heat diffusion terms and heat production terms where the heat production is due to mechanical processes. Commonly, the heat is produced locally along the slip lines and if the heat production outweighs diffusion the material is locally weakened which eventually leads to failure. The effect of diffusion and heat production is captured by a dimensionless quantity, the Gruntfest number, and only if the Gruntfest number is larger than a critical value localized failure occurs. This critical Gruntfest number depends on boundary conditions such as temperature or pressure and hence this critical value gives rise to localization criteria. We find that the results of this approach agree with earlier contributions to the theory of plasticity but gives the advantage of a unified framework which might prove useful in numerical schemes for visco-plasticity.
-
NASA Astrophysics Data System (ADS)
Crosta, G.; Imposimato, S.; Roddeman, D.; Frattini, P.
2012-04-01
Fast moving landslides can be originated along slopes in mountainous terrains with natural and artificial lakes, or fjords at the slope foot. This landslides can reach extremely high speed and the impact with the immobile reservoir water can be influenced by the local topography and the landslide mass profile. The impact can generate large impulse waves and landslide tsunami. Initiation, propagation and runup are the three phases that need to be considered. The landslide evolution and the consequent wave can be controlled by the initial mass position (subaerial, partially or completely submerged), the landslide speed, the type of material, the subaerial and subaqueous slope geometry, the landslide depth and length at the impact, and the water depth. Extreme events have been caused by subaerial landslides: the 1963 Vajont rockslide (Italy), the 1958 Lituya Bay event (Alaska), the Tafjord and the Loen multiple events event (Norway), also from volcanic collapses (Hawaii and Canary islands). Various researchers completed a systematic experimental work on 2D and 3D wave generation and propagation (Kamphuis and Bowering, 1970; Huber, 1980; Müller, 1995; Huber and Hager, 1997; Fritz, 2002; Zweifel, 2004; Panizzo et al., 2005; Heller, 2007; Heller and Kinnear, 2010; Sælevik et al., 2009), using both rigid blocks and deformable granular" masses. Model data and results have been used to calibrate and validate numerical modelling tools (Harbitz, 1992; Jiang and LeBlond, 1993; Grilli et al., 2002; Grilli and Watts, 2005; Lynett and Liu, 2005; Tinti et al., 2006; Abadie et al., 2010) generally considering simplified rheologies (e.g. viscous rheologies) for subaerial subaqueous spreading. We use a FEM code (Roddeman, 2011; Crosta et al., 2006, 2009, 2010, 2011) adopting an Eulerian-Lagrangian approach to give accurate results for large deformations. We model both 2D and fully 3D events considering different settings. The material is considered as a fully deformable elasto-plastic continuum and water as nearly incompressible. In particular we modeled the Vajont rockslide both in 2D and 3D considering the landslide water interaction. More simulations have been performed to validate the model against 2D and 3D tank experiments considering different slope geometries and water depth.
-
Finite Element Modeling of In-Situ Stresses near Salt Bodies
NASA Astrophysics Data System (ADS)
Sanz, P.; Gray, G.; Albertz, M.
2011-12-01
The in-situ stress field is modified around salt bodies because salt rock has no ability to sustain shear stresses. A reliable prediction of stresses near salt is important for planning safe and economic drilling programs. A better understanding of in-situ stresses before drilling can be achieved using finite element models that account for the creeping salt behavior and the elastoplastic response of the surrounding sediments. Two different geomechanical modeling techniques can be distinguished: "dynamic" modeling and "static" modeling. "Dynamic" models, also known as forward models, simulate the development of structural processes in geologic time. This technique provides the evolution of stresses and so it is used to simulate the initiation and development of structural features, such as, faults, folds, fractures, and salt diapers. The original or initial configuration and the unknown final configuration of forward models are usually significantly different therefore geometric non-linearities need to be considered. These models may be difficult to constrain when different tectonic, deposition, and erosion events, and the timing among them, needs to be accounted for. While dynamic models provide insight into the stress evolution, in many cases is very challenging, if not impossible, to forward model a configuration to its known present-day geometry; particularly in the case of salt layers that evolve into highly irregular and complex geometries. Alternatively, "static" models use the present-day geometry and present-day far-field stresses to estimate the present-day in-situ stress field inside a domain. In this case, it is appropriate to use a small deformation approach because initial and final configurations should be very similar, and more important, because the equilibrium of stresses should be stated in the present-day initial configuration. The initial stresses and the applied boundary conditions are constrained by the geologic setting and available data. This modeling technique does not predict the evolution of structural elements or stresses with time; therefore it does not provide any insight into the formation of fractures that were previously developed under a different stress condition or the development of overpressure generated by a high sedimentation rate. This work provides a validation for predicting in-situ stresses near salt using "static" models. We compare synthetic examples using both modeling techniques and show that stresses near salt predicted with "static" models are comparable to the ones generated by "dynamic" models.
-
Modelling the graphite fracture mechanisms
DOE Office of Scientific and Technical Information (OSTI.GOV)
Jacquemoud, C.; Marie, S.; Nedelec, M.
2012-07-01
In order to define a design criterion for graphite components, it is important to identify the physical phenomena responsible for the graphite fracture, to include them in a more effective modelling. In a first step, a large panel of experiments have been realised in order to build up an important database; results of tensile tests, 3 and 4 point bending tests on smooth and notched specimens have been analysed and have demonstrated an important geometry related effects on the behavior up to fracture. Then, first simulations with an elastic or an elastoplastic bilinear constitutive law have not made it possiblemore » to simulate the experimental fracture stress variations with the specimen geometry, the fracture mechanisms of the graphite being at the microstructural scale. That is the reason why a specific F.E. model of the graphite structure has been developed in which every graphite grain has been meshed independently, the crack initiation along the basal plane of the particles as well as the crack propagation and coalescence have been modelled too. This specific model has been used to test two different approaches for fracture initiation: a critical stress criterion and two criteria of fracture mechanic type. They are all based on crystallographic considerations as a global critical stress criterion gave unsatisfactory results. The criteria of fracture mechanic type being extremely unstable and unable to represent the graphite global behaviour up to the final collapse, the critical stress criterion has been preferred to predict the results of the large range of available experiments, on both smooth and notched specimens. In so doing, the experimental observations have been correctly simulated: the geometry related effects on the experimental fracture stress dispersion, the specimen volume effects on the macroscopic fracture stress and the crack propagation at a constant stress intensity factor. In addition, the parameters of the criterion have been related to experimental observations: the local crack initiation stress of 8 MPa corresponds to the non-linearity apparition on the global behavior observed experimentally and the the maximal critical stress defined for the particle of 30 MPa is equivalent to the fracture stress of notched specimens. This innovative combination of crack modelling and a local crystallographic critical stress criterion made it possible to understand that cleavage initiation and propagation in the graphite microstructure was driven by a mean critical stress criterion. (authors)« less
-
NASA Astrophysics Data System (ADS)
Lollino, Piernicola; Andriani, Gioacchino Francesco
2017-07-01
The strength decay that occurs in the post-peak stage, under low confinement stress, represents a key factor of the stress-strain behaviour of rocks. However, for soft rocks this issue is generally underestimated or even neglected in the solution of boundary value problems, as for example those concerning the stability of underground cavities or rocky cliffs. In these cases, the constitutive models frequently used in limit equilibrium analyses or more sophisticated numerical calculations are, respectively, rigid-plastic or elastic-perfectly plastic. In particular, most of commercial continuum-based numerical codes propose a variety of constitutive models, including elasticity, elasto-plasticity, strain-softening and elasto-viscoplasticity, which are not exhaustive in simulating the progressive failure mechanisms affecting brittle rock materials, these being characterized by material detachment and crack opening and propagation. As a consequence, a numerical coupling with mechanical joint propagation is needed to cope with fracture mechanics. Therefore, continuum-based applications that treat the simulation of the failure processes of intact rock masses at low stress levels may need the adoption of numerical techniques capable of implementing fracture mechanics and rock brittleness concepts, as it is shown in this paper. This work is aimed at highlighting, for some applications of rock mechanics, the essential role of post-peak brittleness of soft rocks by means of the application of a hybrid finite-discrete element method. This method allows for a proper simulation of the brittle rock behaviour and the related mechanism of fracture propagation. In particular, the paper presents two ideal problems, represented by a shallow underground cave and a vertical cliff, for which the evolution of the stability conditions is investigated by comparing the solutions obtained implementing different brittle material responses with those resulting from the assumption of perfectly plastic behaviour. To this purpose, a series of petrophysical and mechanical tests were conducted on samples of soft calcarenite belonging to the Calcarenite di Gravina Fm. (Apulia, Southern Italy), focusing specific attention on the post-peak behaviour of the material under three types of loading (compression, indirect tension and shear). Typical geometrical features representative of real rock engineering problems observed in Southern Italy were assumed in the problems examined. The numerical results indicate the impact of soft rock brittleness in the assessment of stability and highlight the need for the adoption of innovative numerical techniques to analyse these types of problems properly.
-
DOE Office of Scientific and Technical Information (OSTI.GOV)
I. M. Robertson; A. Beaudoin; J. Lambros
2004-01-05
OAK-135 Development and validation of constitutive models for polycrystalline materials subjected to high strain rate loading over a range of temperatures are needed to predict the response of engineering materials to in-service type conditions (foreign object damage, high-strain rate forging, high-speed sheet forming, deformation behavior during forming, response to extreme conditions, etc.). To account accurately for the complex effects that can occur during extreme and variable loading conditions, requires significant and detailed computational and modeling efforts. These efforts must be closely coupled with precise and targeted experimental measurements that not only verify the predictions of the models, but also providemore » input about the fundamental processes responsible for the macroscopic response. Achieving this coupling between modeling and experimentation is the guiding principle of this program. Specifically, this program seeks to bridge the length scale between discrete dislocation interactions with grain boundaries and continuum models for polycrystalline plasticity. Achieving this goal requires incorporating these complex dislocation-interface interactions into the well-defined behavior of single crystals. Despite the widespread study of metal plasticity, this aspect is not well understood for simple loading conditions, let alone extreme ones. Our experimental approach includes determining the high-strain rate response as a function of strain and temperature with post-mortem characterization of the microstructure, quasi-static testing of pre-deformed material, and direct observation of the dislocation behavior during reloading by using the in situ transmission electron microscope deformation technique. These experiments will provide the basis for development and validation of physically-based constitutive models, which will include dislocation-grain boundary interactions for polycrystalline systems. One aspect of the program will involve the dire ct observation of specific mechanisms of micro-plasticity, as these will indicate the boundary value problem that should be addressed. This focus on the pre-yield region in the quasi-static effort (the elasto-plastic transition) is also a tractable one from an experimental and modeling viewpoint. In addition, our approach will minimize the need to fit model parameters to experimental data to obtain convergence. These are critical steps to reach the primary objective of simulating and modeling material performance under extreme loading conditions. In this annual report, we describe the progress made in the first year of this program.« less
-
Verification of the GIS-based Newmark method through 2D dynamic modelling of slope stability
NASA Astrophysics Data System (ADS)
Torgoev, A.; Havenith, H.-B.
2012-04-01
The goal of this work is to verify the simplified GIS-based Newmark displacement approach through 2D dynamic modelling of slope stability. The research is applied to a landslide-prone area in Central Asia, the Mailuu-Suu Valley, situated in the south of Kyrgyzstan. The comparison is carried out on the basis of 30 different profiles located in the target area, presenting different geological, tectonic and morphological settings. One part of the profiles were selected within landslide zones, the other part was selected in stable areas. Many of the landslides are complex slope failures involving falls, rotational sliding and/or planar sliding and flows. These input data were extracted from a 3D structural geological model built with the GOCAD software. Geophysical and geomechanical parameters were defined on the basis of results obtained by multiple surveys performed in the area over the past 15 years. These include geophysical investigation, seismological experiments and ambient noise measurements. Dynamic modelling of slope stability is performed with the UDEC version 4.01 software that is able to compute deformation of discrete elements. Inside these elements both elasto-plastic and purely elastic materials (similar to rigid blocks) were tested. Various parameter variations were tested to assess their influence on the final outputs. And even though no groundwater flow was included, the numerous simulations are very time-consuming (20 mins per model for 10 secs simulated shaking) - about 500 computation hours have been completed so far (more than 100 models). Preliminary results allow us to compare Newmark displacements computed using different GIS approaches (Jibson et al., 1998; Miles and Ho, 1999, among others) with the displacements computed using the original Newmark method (Newmark, 1965, here simulated seismograms were used) and displacements produced along joints by the corresponding 2D dynamical models. The generation of seismic amplification and its impact on peak-ground-acceleration, Arias Intensity and permanent slope movements (total and slip on joints) is assessed for numerous morphological-lithological settings (curvature, slope angle, surficial geology, various layer dips and orientations) throughout the target area. The final results of our studies should allow us to define the limitations of the simplified GIS-based Newmark displacement modelling; thus, the verified method would make landslide susceptibility and hazard mapping in seismically active regions more reliable.
-
Design and Experimental Verification of Deployable/Inflatable Ultra-Lightweight Structures
NASA Technical Reports Server (NTRS)
Pai, P. Frank
2004-01-01
Because launch cost of a space structural system is often proportional to the launch volume and mass and there is no significant gravity in space, NASA's space exploration programs and various science missions have stimulated extensive use of ultra-lightweight deployable/inflatable structures. These structures are named here as Highly Flexible Structures (HFSs) because they are designed to undergo large displacements, rotations, and/or buckling without plastic deformation under normal operation conditions. Except recent applications to space structural systems, HFSs have been used in many mechanical systems, civil structures, aerospace vehicles, home appliances, and medical devices to satisfy space limitations, provide special mechanisms, and/or reduce structural weight. The extensive use of HFSs in today's structural engineering reveals the need of a design and analysis software and a database system with design guidelines for practicing engineers to perform computer-aided design and rapid prototyping of HFSs. Also to prepare engineering students for future structural engineering requires a new and easy-to- understand method of presenting the complex mathematics of the modeling and analysis of HFSs. However, because of the high flexibility of HFSs, many unique challenging problems in the modeling, design and analysis of HFSs need to be studied. The current state of research on HFSs needs advances in the following areas: (1) modeling of large rotations using appropriate strain measures, (2) modeling of cross-section warpings of structures, (3) how to account for both large rotations and cross- section warpings in 2D (two-dimensional) and 1D structural theories, (4) modeling of thickness thinning of membranes due to inflation pressure, pretension, and temperature change, (5) prediction of inflated shapes and wrinkles of inflatable structures, (6) development of efficient numerical methods for nonlinear static and dynamic analyses, and (7) filling the gap between geometrically exact elastic analysis and elastoplastic analysis. The objectives of this research project were: (1) to study the modeling, design, and analysis of deployable/inflatable ultra-lightweight structures, (2) to perform numerical and experimental studies on the static and dynamic characteristics and deployability of HFSs, (3) to derive guidelines for designing HFSs, (4) to develop a MATLAB toolbox for the design, analysis, and dynamic animation of HFSs, and (5) to perform experiments and establish an adequate database of post-buckling characteristics of HFSs.
-
ISALE impact simulations in support of AIDA mission
NASA Astrophysics Data System (ADS)
Oklay, Nilda; Vincent, Jean-Baptiste; Michel, Patrick; Schwartz, Stephen
2016-07-01
Introduction: The Asteroid Impact Deflection Assessment (AIDA) mission is a joint project of ESA and NASA with two independent spacecraft. ESA's contribution is an observer satellite called Asteroid Impact Mission (AIM, [1]), and NASA's contribution is a projectile called Double Asteroid Redirection Test (DART, [2]). The target of the mission is a near-Earth binary asteroid system (65803) Didymos. The aim is to study the possibility of deflecting an asteroid by using a kinetic impactor, as well as to characterize the internal properties of the target and test various relevant technologies for other missions. The design is that the DART would impact the secondary of the binary system and AIM would characterize the target asteroid, observe the impact event and measure the changes in the relative orbit after the impact. Impact modeling will be used to interpret the results of the AIDA impact event. There are numerous impact simulation codes, which are planned to be used to understand the AIDA impact results. Therefore an international benchmarking program is ongoing for the comparison of the results of various codes on the defined test cases [3]. We will present the results of the test cases performed by iSALE hydrocode. Modeling: In this work we use the iSALE-2D shock physics code [4], which is based on the SALE hydrocode solution algorithm [5]. To simulate hypervelocity impact processes in solid materials SALE was modified to include an elastoplastic constitutive model, fragmentation models, various EOS, and multiple materials [6, 7]. More recent improvements include a modified strength model [8] and a porosity compaction model [4, 9]. References: [1] Michel P. et al., 2016, ASR, submitted [2] Cheng A. F. et al., (2016) PSS, 121, 27-35 [3] Stickle A. M. et al., (2016). 47th LPSC [4] Wünnemann,K. et al., (2006). Icarus, 180:514-527 [5] Amsden, A., et al., (1980) LANL Report, LA-8095:101p. [6] Melosh, H. J., et al., (1992). J. Geophys. Res., 97(E9):14735-14759 [7] Ivanov, B. A., et al., (1997) Int. J. Imp. Eng., 20:411-430; [8] Collins, G. S., et al., (2004). Met. & Planet. Sci., 39:217-231. [9] Collins, G., et al., (2011) Int. J. Imp. Eng., 38:434-439