-
Effects of transverse temperature field nonuniformity on stress in silicon sheet growth
NASA Technical Reports Server (NTRS)
Mataga, P. A.; Hutchinson, J. W.; Chalmers, B.; Bell, R. O.; Kalejs, J. P.
1987-01-01
Stress and strain rate distributions are calculated using finite element analysis for steady-state growth of thin silicon sheet temperature nonuniformities imposed in the transverse (sheet width) dimension. Significant reductions in residual stress are predicted to occur for the case where the sheet edge is cooled relative to its center provided plastic deformation with high creep rates is present.
-
Effect of material variation on the biomechanical behaviour of orthodontic fixed appliances: a finite element analysis
PubMed Central
Papageorgiou, Spyridon N.; Keilig, Ludger; Hasan, Istabrak; Jäger, Andreas
2016-01-01
Summary Introduction: Biomechanical analysis of orthodontic tooth movement is complex, as many different tissues and appliance components are involved. The aim of this finite element study was to assess the relative effect of material alteration of the various components of the orthodontic appliance on the biomechanical behaviour of tooth movement. Methods: A three-dimensional finite element solid model was constructed. The model consisted of a canine, a first, and a second premolar, including the surrounding tooth-supporting structures and fixed appliances. The materials of the orthodontic appliances were alternated between: (1) composite resin or resin-modified glass ionomer cement for the adhesive, (2) steel, titanium, ceramic, or plastic for the bracket, and (3) β-titanium or steel for the wire. After vertical activation of the first premolar by 0.5mm in occlusal direction, stress and strain calculations were performed at the periodontal ligament and the orthodontic appliance. Results: The finite element analysis indicated that strains developed at the periodontal ligament were mainly influenced by the orthodontic wire (up to +63 per cent), followed by the bracket (up to +44 per cent) and the adhesive (up to +4 per cent). As far as developed stresses at the orthodontic appliance are concerned, wire material had the greatest influence (up to +155 per cent), followed by bracket material (up to +148 per cent) and adhesive material (up to +8 per cent). Limitations: The results of this in silico study need to be validated by in vivo studies before they can be extrapolated to clinical practice. Conclusion: According to the results of this finite element study, all components of the orthodontic fixed appliance, including wire, bracket, and adhesive, seem to influence, to some extent, the biomechanics of tooth movement. PMID:26174769
-
VALIDATION OF ANSYS FINITE ELEMENT ANALYSIS SOFTWARE
DOE Office of Scientific and Technical Information (OSTI.GOV)
HAMM, E.R.
2003-06-27
This document provides a record of the verification and Validation of the ANSYS Version 7.0 software that is installed on selected CH2M HILL computers. The issues addressed include: Software verification, installation, validation, configuration management and error reporting. The ANSYS{reg_sign} computer program is a large scale multi-purpose finite element program which may be used for solving several classes of engineering analysis. The analysis capabilities of ANSYS Full Mechanical Version 7.0 installed on selected CH2M Hill Hanford Group (CH2M HILL) Intel processor based computers include the ability to solve static and dynamic structural analyses, steady-state and transient heat transfer problems, mode-frequency andmore » buckling eigenvalue problems, static or time-varying magnetic analyses and various types of field and coupled-field applications. The program contains many special features which allow nonlinearities or secondary effects to be included in the solution, such as plasticity, large strain, hyperelasticity, creep, swelling, large deflections, contact, stress stiffening, temperature dependency, material anisotropy, and thermal radiation. The ANSYS program has been in commercial use since 1970, and has been used extensively in the aerospace, automotive, construction, electronic, energy services, manufacturing, nuclear, plastics, oil and steel industries.« less
-
Drop Testing Representative Multi-Canister Overpacks
DOE Office of Scientific and Technical Information (OSTI.GOV)
Snow, Spencer D.; Morton, Dana K.
The objective of the work reported herein was to determine the ability of the Multi- Canister Overpack (MCO) canister design to maintain its containment boundary after an accidental drop event. Two test MCO canisters were assembled at Hanford, prepared for testing at the Idaho National Engineering and Environmental Laboratory (INEEL), drop tested at Sandia National Laboratories, and evaluated back at the INEEL. In addition to the actual testing efforts, finite element plastic analysis techniques were used to make both pre-test and post-test predictions of the test MCOs structural deformations. The completed effort has demonstrated that the canister design is capablemore » of maintaining a 50 psig pressure boundary after drop testing. Based on helium leak testing methods, one test MCO was determined to have a leakage rate not greater than 1x10 -5 std cc/sec (prior internal helium presence prevented a more rigorous test) and the remaining test MCO had a measured leakage rate less than 1x10 -7 std cc/sec (i.e., a leaktight containment) after the drop test. The effort has also demonstrated the capability of finite element methods using plastic analysis techniques to accurately predict the structural deformations of canisters subjected to an accidental drop event.« less
-
Porosity Defect Remodeling and Tensile Analysis of Cast Steel
PubMed Central
Sun, Linfeng; Liao, Ridong; Lu, Wei; Fu, Sibo
2016-01-01
Tensile properties on ASTM A216 WCB cast steel with centerline porosity defect were studied with radiographic mapping and finite element remodeling technique. Non-linear elastic and plastic behaviors dependent on porosity were mathematically described by relevant equation sets. According to the ASTM E8 tensile test standard, matrix and defect specimens were machined into two categories by two types of height. After applying radiographic inspection, defect morphologies were mapped to the mid-sections of the finite element models and the porosity fraction fields had been generated with interpolation method. ABAQUS input parameters were confirmed by trial simulations to the matrix specimen and comparison with experimental outcomes. Fine agreements of the result curves between simulations and experiments could be observed, and predicted positions of the tensile fracture were found to be in accordance with the tests. Chord modulus was used to obtain the equivalent elastic stiffness because of the non-linear features. The results showed that elongation was the most influenced term to the defect cast steel, compared with elastic stiffness and yield stress. Additional visual explanations on the tensile fracture caused by void propagation were also given by the result contours at different mechanical stages, including distributions of Mises stress and plastic strain. PMID:28787919
-
Stress and efficiency studies in EFG
NASA Technical Reports Server (NTRS)
1986-01-01
The goals of this program were: (1) to define minimum stress configurations for silicon sheet growth at high speeds; (2) to quantify dislocation electrical activity and their limits on minority carrier diffusion length in deformed silicon; and (3) to study reasons for degradation of lifetime with increases in doping level in edge-defined film-fed growth (EFG) materials. A finite element model was developed for calculating residual stress with plastic deformation. A finite element model was verified for EFG control variable relationships to temperature field of the sheet to permit prediction of profiles and stresses encountered in EFG systems. A residual stress measurement technique was developed for finite size EFG material blanks using shadow Moire interferometry. Transient creep response of silicon was investigated in the temperature range between 800 and 1400 C in strain and strain regimes of interest in stress analysis of sheet growth. Quantitative relationships were established between minority carrier diffusion length and dislocation densities using Electron Beam Induced Current (EBIC) measurement in FZ silicon deformed in four point bending tests.
-
Explicit 2-D Hydrodynamic FEM Program
DOE Office of Scientific and Technical Information (OSTI.GOV)
Lin, Jerry
1996-08-07
DYNA2D* is a vectorized, explicit, two-dimensional, axisymmetric and plane strain finite element program for analyzing the large deformation dynamic and hydrodynamic response of inelastic solids. DYNA2D* contains 13 material models and 9 equations of state (EOS) to cover a wide range of material behavior. The material models implemented in all machine versions are: elastic, orthotropic elastic, kinematic/isotropic elastic plasticity, thermoelastoplastic, soil and crushable foam, linear viscoelastic, rubber, high explosive burn, isotropic elastic-plastic, temperature-dependent elastic-plastic. The isotropic and temperature-dependent elastic-plastic models determine only the deviatoric stresses. Pressure is determined by one of 9 equations of state including linear polynomial, JWL highmore » explosive, Sack Tuesday high explosive, Gruneisen, ratio of polynomials, linear polynomial with energy deposition, ignition and growth of reaction in HE, tabulated compaction, and tabulated.« less
-
Buckling behavior of origami unit cell facets under compressive loads
NASA Astrophysics Data System (ADS)
Kshad, Mohamed Ali Emhmed; Naguib, Hani E.
2018-03-01
Origami structures as cores for sandwich structures are designed to withstand the compressive loads and to dissipate compressive energy. The deformation of the origami panels and the unit cell facets are the primary factors behind the compressive energy dissipation in origami structures. During the loading stage, the origami structures deform through the folding and unfolding process of the unit cell facets, and also through the plastic deformation of the facets. This work presents a numerical study of the buckling behavior of different origami unit cell elements under compressive loading. The studied origami configurations were Miura and Ron-Resch-like origami structures. Finite element package was used to model the origami structures. The study investigated the buckling behavior of the unit cell facets of two types of origami structures Miura origami and Ron-Resch-Like origami structures. The simulation was conducted using ANSYS finite element software, in which the model of the unit cell represented by shell elements, and the eigenvalues buckling solver was used to predict the theoretical buckling of the unit cell elements.
-
Variations of a global constraint factor in cracked bodies under tension and bending loads
NASA Technical Reports Server (NTRS)
Newman, J. C., Jr.; Crews, J. H., Jr.; Bigelow, C. A.; Dawicke, D. S.
1994-01-01
Elastic-plastic finite-element analyses were used to calculate stresses and displacements around a crack in finite-thickness plates for an elastic-perfectly plastic material. Middle- and edge-crack specimens were analyzed under tension and bending loads. Specimens were 1.25 to 20 mm thick with various widths and crack lengths. A global constraint factor alpha(sub g), an averaged normal-stress to flow-stress ratio over the plastic region, was defined to simulate three-dimensional (3D) effects in two-dimensional (2D) models. For crack lengths and uncracked ligament lengths greater than four times the thickness, the global constraint factor was found to be nearly a unique function of a normalized stress-intensity factor (related to plastic-zone size to thickness ratio) from small- to large-scale yielding conditions for various specimen types and thickness. For crack length-to-thickness ratios less than four, the global constraint factor was specimen type, crack length and thickness dependent. Using a 2D strip-yield model and the global constraint factors, plastic-zone sizes and crack-tip displacements agreed reasonably well with the 3D analyses. For a thin sheet aluminum alloy, the critical crack-tip-opening angle during stable tearing was found to be independent of specimen type and crack length for crack length-to-thickness ratios greater than 4.
-
Indentation-derived elastic modulus of multilayer thin films: Effect of unloading induced plasticity
DOE PAGES
Jamison, Ryan Dale; Shen, Yu -Lin
2015-08-13
Nanoindentation is useful for evaluating the mechanical properties, such as elastic modulus, of multilayer thin film materials. A fundamental assumption in the derivation of the elastic modulus from nanoindentation is that the unloading process is purely elastic. In this work, the validity of elastic assumption as it applies to multilayer thin films is studied using the finite element method. The elastic modulus and hardness from the model system are compared to experimental results to show validity of the model. Plastic strain is shown to increase in the multilayer system during the unloading process. Additionally, the indentation-derived modulus of a monolayermore » material shows no dependence on unloading plasticity while the modulus of the multilayer system is dependent on unloading-induced plasticity. Lastly, the cyclic behavior of the multilayer thin film is studied in relation to the influence of unloading-induced plasticity. Furthermore, it is found that several cycles are required to minimize unloading-induced plasticity.« less
-
Crack Growth Prediction Methodology for Multi-Site Damage: Layered Analysis and Growth During Plasticity
NASA Technical Reports Server (NTRS)
James, Mark Anthony
1999-01-01
A finite element program has been developed to perform quasi-static, elastic-plastic crack growth simulations. The model provides a general framework for mixed-mode I/II elastic-plastic fracture analysis using small strain assumptions and plane stress, plane strain, and axisymmetric finite elements. Cracks are modeled explicitly in the mesh. As the cracks propagate, automatic remeshing algorithms delete the mesh local to the crack tip, extend the crack, and build a new mesh around the new tip. State variable mapping algorithms transfer stresses and displacements from the old mesh to the new mesh. The von Mises material model is implemented in the context of a non-linear Newton solution scheme. The fracture criterion is the critical crack tip opening displacement, and crack direction is predicted by the maximum tensile stress criterion at the crack tip. The implementation can accommodate multiple curving and interacting cracks. An additional fracture algorithm based on nodal release can be used to simulate fracture along a horizontal plane of symmetry. A core of plane strain elements can be used with the nodal release algorithm to simulate the triaxial state of stress near the crack tip. Verification and validation studies compare analysis results with experimental data and published three-dimensional analysis results. Fracture predictions using nodal release for compact tension, middle-crack tension, and multi-site damage test specimens produced accurate results for residual strength and link-up loads. Curving crack predictions using remeshing/mapping were compared with experimental data for an Arcan mixed-mode specimen. Loading angles from 0 degrees to 90 degrees were analyzed. The maximum tensile stress criterion was able to predict the crack direction and path for all loading angles in which the material failed in tension. Residual strength was also accurately predicted for these cases.
-
Application of 1 D Finite Element Method in Combination with Laminar Solution Method for Pipe Network Analysis
NASA Astrophysics Data System (ADS)
Dudar, O. I.; Dudar, E. S.
2017-11-01
The features of application of the 1D dimensional finite element method (FEM) in combination with the laminar solutions method (LSM) for the calculation of underground ventilating networks are considered. In this case the processes of heat and mass transfer change the properties of a fluid (binary vapour-air mix). Under the action of gravitational forces it leads to such phenomena as natural draft, local circulation, etc. The FEM relations considering the action of gravity, the mass conservation law, the dependence of vapour-air mix properties on the thermodynamic parameters are derived so that it allows one to model the mentioned phenomena. The analogy of the elastic and plastic rod deformation processes to the processes of laminar and turbulent flow in a pipe is described. Owing to this analogy, the guaranteed convergence of the elastic solutions method for the materials of plastic type means the guaranteed convergence of the LSM for any regime of a turbulent flow in a rough pipe. By means of numerical experiments the convergence rate of the FEM - LSM is investigated. This convergence rate appeared much higher than the convergence rate of the Cross - Andriyashev method. Data of other authors on the convergence rate comparison for the finite element method, the Newton method and the method of gradient are provided. These data allow one to conclude that the FEM in combination with the LSM is one of the most effective methods of calculation of hydraulic and ventilating networks. The FEM - LSM has been used for creation of the research application programme package “MineClimate” allowing to calculate the microclimate parameters in the underground ventilating networks.
-
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
-
The effect of plate-scale rheology and plate interactions on intraplate seismicity
NASA Astrophysics Data System (ADS)
So, Byung-Dal; Capitanio, Fabio A.
2017-11-01
We use finite element modeling to investigate on the stress loading-unloading cycles and earthquakes occurrence in the plate interiors, resulting from the interactions of tectonic plates along their boundary. We model a visco-elasto-plastic plate embedding a single or multiple faults, while the tectonic stress is applied along the plate boundary by an external loading visco-elastic plate, reproducing the tectonic setting of two interacting lithospheres. Because the two plates deform viscously, the timescale of stress accumulation and release on the faults is self-consistently determined, from the boundary to the interiors, and seismic recurrence is an emerging feature. This approach overcomes the constraints on recurrence period imposed by stress (stress-drop) and velocity boundary conditions, while here it is unconstrained. We illustrate emerging macroscopic characteristics of this system, showing that the seismic recurrence period τ becomes shorter as Γ and Θ decreases, where Γ =ηI /ηL, the viscosity ratio of the viscosities of the internal fault-embedded to external loading plates, respectively, and Θ =σY /σL the stress ratio of the elastic limit of the fault to far-field loading stress. When the system embeds multiple, randomly distributed faults, stress transfer results in recurrence period deviations, however the time-averaged recurrence period of each fault show the same dependence on Γ and Θ, illustrating a characteristic collective behavior. The control of these parameters prevails even when initial pre-stress was randomly assigned in terms of the spatial arrangement and orientation on the internal plate, mimicking local fluctuations. Our study shows the relevance of macroscopic rheological properties of tectonic plates on the earthquake occurrence in plate interiors, as opposed to local factors, proposing a viable model for the seismic behavior of continent interiors in the context of large-scale, long-term deformation of interacting tectonic plates.
-
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.
-
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.
-
Numerical modeling of intraplate seismicity with a deformable loading plate
NASA Astrophysics Data System (ADS)
So, B. D.; Capitanio, F. A.
2017-12-01
We use finite element modeling to investigate on the stress loading-unloading cycles and earthquakes occurrence in the plate interiors, resulting from the interactions of tectonic plates along their boundary. We model a visco-elasto-plastic plate embedding a single or multiple faults, while the tectonic stress is applied along the plate boundary by an external loading visco-elastic plate, reproducing the tectonic setting of two interacting lithospheres. Because the two plates deform viscously, the timescale of stress accumulation and release on the faults is self-consistently determined, from the boundary to the interiors, and seismic recurrence is an emerging feature. This approach overcomes the constraints on recurrence period imposed by stress (stress-drop) and velocity boundary conditions, while here it is unconstrained. We illustrate emerging macroscopic characteristics of this system, showing that the seismic recurrence period τ becomes shorter as Γ and Θ decreases, where Γ = ηI/ηL the viscosity ratio of the viscosities of the internal fault-embedded to external loading plates, respectively, and Θ = σY/σL the stress ratio of the elastic limit of the fault to far-field loading stress. When the system embeds multiple, randomly distributed faults, stress transfer results in recurrence period deviations, however the time-averaged recurrence period of each fault show the same dependence on Γ and Θ, illustrating a characteristic collective behavior. The control of these parameters prevails even when initial pre-stress was randomly assigned in terms of the spatial arrangement and orientation on the internal plate, mimicking local fluctuations. Our study shows the relevance of macroscopic rheological properties of tectonic plates on the earthquake occurrence in plate interiors, as opposed to local factors, proposing a viable model for the seismic behavior of continent interiors in the context of large-scale, long-term deformation of interacting tectonic plates.
-
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.
-
Three dimensional modelling of earthquake rupture cycles on frictional faults
NASA Astrophysics Data System (ADS)
Simpson, Guy; May, Dave
2017-04-01
We are developing an efficient MPI-parallel numerical method to simulate earthquake sequences on preexisting faults embedding within a three dimensional viscoelastic half-space. We solve the velocity form of the elasto(visco)dynamic equations using a continuous Galerkin Finite Element Method on an unstructured pentahedral mesh, which thus permits local spatial refinement in the vicinity of the fault. Friction sliding is coupled to the viscoelastic solid via rate- and state-dependent friction laws using the split-node technique. Our coupled formulation employs a picard-type non-linear solver with a fully implicit, first order accurate time integrator that utilises an adaptive time step that efficiently evolves the system through multiple seismic cycles. The implementation leverages advanced parallel solvers, preconditioners and linear algebra from the Portable Extensible Toolkit for Scientific Computing (PETSc) library. The model can treat heterogeneous frictional properties and stress states on the fault and surrounding solid as well as non-planar fault geometries. Preliminary tests show that the model successfully reproduces dynamic rupture on a vertical strike-slip fault in a half-space governed by rate-state friction with the ageing law.
-
Structure and Composition of Isolated Core-Shell (In ,Ga )N /GaN Rods Based on Nanofocus X-Ray Diffraction and Scanning Transmission Electron Microscopy
NASA Astrophysics Data System (ADS)
Krause, Thilo; Hanke, Michael; Nicolai, Lars; Cheng, Zongzhe; Niehle, Michael; Trampert, Achim; Kahnt, Maik; Falkenberg, Gerald; Schroer, Christian G.; Hartmann, Jana; Zhou, Hao; Wehmann, Hergo-Heinrich; Waag, Andreas
2017-02-01
Nanofocus x-ray diffraction is used to investigate the structure and local strain field of an isolated (In ,Ga )N /GaN core-shell microrod. Because the high spatial resolution of the x-ray beam is only 80 ×90 nm2, we are able to investigate several distinct volumes on one individual side facet. Here, we find a drastic increase in thickness of the outer GaN shell along the rod height. Additionally, we performed high-angle annular dark-field scanning-transmission-electron-microscopy measurements on several rods from the same sample showing that (In,Ga)N double-quantum-well and GaN barrier thicknesses also increase strongly along the height. Moreover, plastic relaxation is observed in the top part of the rod. Based on the experimentally obtained structural parameters, we simulate the strain-induced deformation using the finite-element method, which serves as the input for subsequent kinematic scattering simulations. The simulations reveal a significant increase of elastic in-plane relaxation along the rod height. However, at a certain height, the occurrence of plastic relaxation yields a decrease of the elastic strain. Because of the experimentally obtained structural input for the finite-element simulations, we can exclude unknown structural influences on the strain distribution, and we are able to translate the elastic relaxation into an indium concentration which increases by a factor of 4 from the bottom to the height where plastic relaxation occurs.