Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “elastic model”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8

Analytic model of the γ -surface deviation and influence on the stacking fault width between partial dislocations

The stacking fault width (R$_{eq}$) between partial dislocations within an FCC crystalline lattice characterizes the onset of numerous plastic flow mechanisms, as well as the relationship between material strength and grain size. Continuum models traditionally consider a complete unit of dislocation slip ( b ) along the $\langle{110}$$\rangle$ direction distributed between two discrete partial dislocations, each with a fixed partial Burgers vector ( b p ), which bound a stacking fault. Across the stacking fault, the vectorial slip is assumed to be constant, yielding a constant intrinsic stacking fault energy density, $\gamma$$_{isf}$ Here, we demonstrate that the vectorial displacement path taken in accomplishing a complete unit of slip (b) deviates from the expected displacement path containing the local minima, $\gamma$$_{isf}$ leading to a correction in the nominally derived stacking fault width. The magnitude of the correction depends on both the net orientation of the dislocation within the lattice, and also the degree of relaxation of each partial Burgers vector along the $\langle{\bar{1}12}$$\rangle$ direction. We derive a simple analytic model for the corrected stacking fault width, explicitly accounting for the deviation, by introducing a variable ξ (0 < ξ < 1), which governs the magnitude of each partial dislocation component along the $\langle{\bar{1}12}$$\rangle$ direction. Significantly, our model predicts a correction to R$_{eq}$ of ~$\mathscr{O}$ (|b|) for the nominally screw dislocation, and has no noticeable influence on nominally edge dislocations. We apply our model towards computing the stacking fault width of several FCC metals, and demonstrate excellent agreement both with our own numerical data as well as that obtained from ab initio and Molecular Statics (MS) methods within the literature. Lastly, the results from this study demonstrate that, upon judicious application, discrete linear elastic models are successful in reproducing elastic interactions as computed from higher fidelity models on the spatial scale of metallic dislocation cores.

36 MATERIALS SCIENCE↗

Hydrodynamic Impact of a System with a Single Elastic Mode II : Comparison of Experimental Force and Response with Theory

Hydrodynamic impact tests were made on an elastic model approximating a two-mass spring system to determine experimentally the effects of structural flexibility on the hydrodynamic loads encountered during seaplane landing impacts and to correlate the results with theory. A flexible seaplane was represented by a two-mass spring system consisting of a rigid prismatic float connected to a rigid upper mass by an elastic structure. The model had a ratio of sprung mass to hull mass of 0.6 and a natural frequency of 3.0 cycles per second. The tests were conducted in smooth water at fixed trims and included both high and low flight-path angles and a range of velocity. Theoretical and experimental comparisons indicated that the theoretical results agreed well with the experimental results.

Miller, Robert W↗

Modeling and Detection of Future Cyber-Enabled DSM Data Attacks

Demand-Side Management (DSM) is an essential tool to ensure power system reliability and stability. In future smart grids, certain portions of a customer’s load usage could be under the automatic control of a cyber-enabled DSM program, which selectively schedules loads as a function of electricity prices to improve power balance and grid stability. In this scenario, the security of DSM cyberinfrastructure will be critical as advanced metering infrastructure and communication systems are susceptible to cyber-attacks. Such attacks, in the form of false data injections, can manipulate customer load profiles and cause metering chaos and energy losses in the grid. The feedback mechanism between load management on the consumer side and dynamic price schemes employed by independent system operators can further exacerbate attacks. To study how this feedback mechanism may worsen attacks in future cyber-enabled DSM programs, we propose a novel mathematical framework for (i) modeling the nonlinear relationship between load management and real-time pricing, (ii) simulating residential load data and prices, (iii) creating cyber-attacks, and (iv) detecting said attacks. In this framework, we first develop time-series forecasts to model load demand and use them as inputs to an elasticity model for the price-demand relationship in the DSM loop. This work then investigates the behavior of such a feedback loop under intentional cyber-attacks. We simulate and examine load-price data under different DSM-participation levels with three types of random additive attacks: ramp, sudden, and point attacks. We conduct two investigations for the detection of DSM attacks. The first studies a supervised learning approach, with various classification models, and the second studies the performance of parametric and nonparametric change point detectors. Results conclude that higher amounts of DSM participation can exacerbate ramp and sudden attacks leading to better detection of such attacks, especially with supervised learning classifiers. We also find that nonparametric detection outperforms parametric for smaller user pools, and random point attacks are the hardest to detect with any method.

97 MATHEMATICS AND COMPUTING↗

Finite element modeling of stress in the Nazca plate - Driving forces and plate boundary earthquakes

The state of stress within the Nazca plate due to plate driving forces and large plate boundary earthquakes has been analyzed by applying a finite element method using the wave front solution technique to models of the intraplate stress field in a single plate using a refined grid. Although only static elastic models have been explicitly calculated, certain limiting cases of an elastic plate over a viscous asthenosphere were also treated. A state of nearly east-west compression inferred from the source mechanism of thrust earthquakes in the interior of the plate requires ridge pushing forces. The net pulling force on the oceanic plate by the subducted slab has a maximum value comparable to pushing forces. The estimated horizontal deviatoric stress in intraplate regions, based on potential forces associated with the ridge, is on the order of a few hundred bars. The intraplate stress field in the region of the 1960 earthquake may change by a few tens of bars at most once the asthenosphere has relaxed, with changes on the order of one bar occurring at greater distances into the plate. The changes in the intraplate stress field are probably not noticeable unless the lithosphere is near failure.

Richardson, R. M.↗

High-performance modeling acoustic and elastic waves using the parallel Dichotomy Algorithm

A high-performance parallel algorithm is proposed for modeling the propagation of acoustic and elastic waves in inhomogeneous media. An initial boundary-value problem is replaced by a series of boundary-value problems for a constant elliptic operator and different right-hand sides via the integral Laguerre transform. It is proposed to solve difference equations by the conjugate gradient method for acoustic equations and by the GMRES(k) method for modeling elastic waves. A preconditioning operator was the Laplace operator that is inverted using the variable separation method. The novelty of the proposed algorithm is using the Dichotomy Algorithm , which was designed for solving a series of tridiagonal systems of linear equations, in the context of the preconditioning operator inversion. Via considering analytical solutions, it is shown that modeling wave processes for long instants of time requires high-resolution meshes. The proposed parallel fine-mesh algorithm enabled to solve real application seismic problems in acceptable time and with high accuracy. By solving model problems, it is demonstrated that the considered parallel algorithm possesses high performance and efficiency over a wide range of the number of processors (from 2 to 8192).

97 MATHEMATICS AND COMPUTING↗

Vector elastic deconvolution migration with dual wavefield decomposition

Here, we report elastic-wave imaging using multicomponent data can provide more useful subsurface information than acoustic-wave imaging, but it is usually algorithmically challenging. We have developed a vector elastic deconvolution migration method for high-resolution imaging of subsurface structures in isotropic and anisotropic elastic media. Our new method uses a vector deconvolution imaging condition based on dual wavefield decomposition, including an explicit directional wavefield separation using the Hilbert transform and a P/S vector wavefield decomposition using the low-rank decomposition method. Using three elastic models, we numerically determine that our new method produces notably higher resolution and more amplitude-balanced elastic images compared with a crosscorrelation-based vector elastic reverse time migration method.

58 GEOSCIENCES↗

Single crystal plasticity by modeling dislocation density rate behavior

The goal of this work is to formulate a constitutive model for the deformation of metals over a wide range of strain rates. Damage and failure of materials frequently occurs at a variety of deformation rates within the same sample. The present state of the art in single crystal constitutive models relies on thermally-activated models which are believed to become less reliable for problems exceeding strain rates of 10{sup 4} s{sup -1}. This talk presents work in which we extend the applicability of the single crystal model to the strain rate region where dislocation drag is believed to dominate. The elastic model includes effects from volumetric change and pressure sensitive moduli. The plastic model transitions from the low-rate thermally-activated regime to the high-rate drag dominated regime. The direct use of dislocation density as a state parameter gives a measurable physical mechanism to strain hardening. Dislocation densities are separated according to type and given a systematic set of interactions rates adaptable by type. The form of the constitutive model is motivated by previously published dislocation dynamics work which articulated important behaviors unique to high-rate response in fcc systems. The proposed material model incorporates thermal coupling. The hardening model tracks the varying dislocation population with respect to each slip plane and computes the slip resistance based on those values. Comparisons can be made between the responses of single crystals and polycrystals at a variety of strain rates. The material model is fit to copper.

36 MATERIALS SCIENCE↗

Turbofan Volume Dynamics Model for Investigations of Aero-Propulso-Servo-Elastic Effects in a Supersonic Commercial Transport

A turbofan simulation has been developed for use in aero-propulso-servo-elastic coupling studies, on supersonic vehicles. A one-dimensional lumped volume approach is used whereby each component (fan, high-pressure compressor, combustor, etc.) is represented as a single volume using characteristic performance maps and conservation equations for continuity, momentum and energy. The simulation is developed in the MATLAB/SIMULINK (The MathWorks, Inc.) environment in order to facilitate controls development, and ease of integration with a future aero-servo-elastic vehicle model being developed at NASA Langley. The complete simulation demonstrated steady state results that closely match a proposed engine suitable for a supersonic business jet at the cruise condition. Preliminary investigation of the transient simulation revealed expected trends for fuel flow disturbances as well as upstream pressure disturbances. A framework for system identification enables development of linear models for controller design. Utilizing this framework, a transfer function modeling an upstream pressure disturbance s impacts on the engine speed is developed as an illustrative case of the system identification. This work will eventually enable an overall vehicle aero-propulso-servo-elastic model

Connolly, Joseph W.↗

SEASAT observations of lithospheric flexure seaward of trenches

Lithospheric flexure seaward of deep ocean trenches in SEASAT altimeter observations of the marine geoid. In fact, mechanical models of lithospheric flexure can be tested directly on the SEASAT altimeter data. A simple elastic model was used for the oceanic lithosphere and, after least squares adjustments, estimates of model parameters were recovered including Outer Rise (OR) amplitude, OR wavelength, and effective lithospheric thickness. Effective lithospheric thickness was recovered for five regions: the Mariana, the Kuril, the Philippine, the Aleutian and the Middle America OR. These results support the suggestion of Bodine et al. (1981) that effective thickness, T, increased with age of lithosphere in approximate accord with the relation T approximately equals x age to the 1/2 power where C approximately equals 4 km x my to the -1/2 power. Altimetric results agree more closely with this relation than do published results based on bathymetric data. The close agreement with the thickness-age relation suggests that there is no longer any need to assume that significant horizontal compression acts across the Kuril, Marianas and Izo-Bonin trenches. This thickness-age relation implies that flexural strength of the oceanic lithosphere is temperature controlled.

Mcadoo, D. C.↗

Seasat observations of lithospheric flexure seaward of trenches

Lithospheric flexure seaward of deep ocean trenches is evident in Seasat altimeter observations of the marine geoid. In fact, mechanical models of lithospheric flexure can be tested directly on the Seasat altimeter data. A simple elastic model has been used for the oceanic lithosphere and, after least squares adjustments, estimates have been recovered of model parameters including outer rise (OR) amplitude, OR wavelength, and effective lithospheric thickness. Effective lithospheric thicknesses have been recovered for six regions: the Mariana, the Kuril, the Philippine, the Aleutian, the Izu-Bonin, and the Middle America OR's. These results support the proposition that effective thickness Te increases with age of lithosphere in approximate accord with the relation Te approximately C x square root of age where C approximately 4 km/square root of (m.y.). In fact, altimetric results agree more closely with this relation than do published results based on bathymetric data. The close agreement with the thickness-age relation suggests that there is no longer any need to assume that significant horizontal compression acts across the Kuril, Marianas, and Izu-Bonin trenches. This thickness-age relation implies that flexural strength of the oceanic lithosphere is temperature controlled.

Mcadoo, D. C.↗

Earthquake deformation cycle in the Andean Back Arc, western Argentina

It is pointed out that the 1977 Caucete earthquake was an east-west compressive event which occurred at a depth of about 17 km in the Andean back arc approximately 450 km east of the Chile trench. The present study is primarily concerned with the mechanics of fault behavior prior to and after the 1977 earthquake, as evidenced by preseismic and postseismic vertical deformation. Simple fault models for each phase of deformation are presented. The models of preseismic slip, coseismic slip, and postseismic slip consist of uniform dislocations in an elastic half-space following the method of Chinnery (1961). The basically elastic models presented demonstrate that deformations in the vicinity of the 1977 Caucete, Argentina, earthquake are consistent with activity on some portion of the coseismic fault. This is evidence that these deformations are directly related to the processes of strain accumulation and release associated with the event.

Reilinger, R.↗

Finite Element Modeling Used to Study Stress Distribution on the Foot

A method to study the stress distribution inside the forefoot during walking was developed at the Cleveland Clinic Foundation by a researcher from the NASA Glenn Research Center. In this method, a semiautomated process was outlined to create a three-dimensional, patient-specific, finite element model (FEM) of the forefoot using magnetic resonance images (MRI). The images were processed in Matlab using the k-nearest neighbor (k-NN) classification algorithm and Sobel edge detection to separate the different tissue types: bone, skin, fat, and muscle. This information was used to create curves and surfaces that were exported to an FEM preprocessor known as Truegrid. In Truegrid, eight-noded or brick elements were created by using surface mapping. The FEM was processed and postprocessed in Abaqus. Material properties of the models were obtained from past experiments such as fat pad confined compression, skin axial and biaxial tests, muscle in vivo compressive tests, and reference literature (bone properties). Nonlinear (hyperelastic) material models were used for the skin (epidermis and dermis), fat, and muscles; and a linear elastic model was used for the bones. Muscle activation during walking yielded uncertainties in the muscle material model since contracted muscles are stiffer than relaxed muscles. These uncertainties were resolved by performing a sensitivity analysis of the muscle material properties. The original properties were multiplied by arbitrary factors of 2, 3, 0.5, and 0.33. The strain and stress distributions, as well as the locations of peak values, were similar in all cases. The peak contact pressure P obtained for each case varied with respect to the applied factor f as follows:

Morales, Nelson↗

Coupling 2D atomistic information to 3D kink-pair enthalpy models of screw dislocations in bcc metals

The kink-pair activation enthalpy is a fundamental parameter in the theory of plasticity of body-centered cubic (bcc) metals. It controls the thermally activated motion of screw dislocation at low and intermediate temperatures. While direct atomistic calculations of kink pairs on screw dislocations have reached a high degree of accuracy, they can only be practically performed using semiempirical interatomic force fields, as electronic structure methods have not yet reached the level of efficiency needed to capture the system sizes required to model kink-pair structures. In this context, an alternative approach based on standard three-dimensional elastic models, which are efficient but lack atomic-level information, coupled to a substrate potential that represents the underlying lattice, has been widely applied over the past few years. This class of methods, known as ‘line-on-substrate’ (LOS) models, uses the substrate potential to calculate the lattice contribution to the kink-pair energies. In this work, we introduce the stress dependence of the substrate potential into LOS models to evaluate its impact on kink-pair energies. In addition, we include asymmetric dislocation core energies as an extra descriptor of the dislocation character. This asymmetry is also elevated to the continuum level by adding core energies to the general LOS formulation and used to explain potential energy differences known to exist between left and right kinks in bcc metals. Finally, more importantly, by matching the total LOS energies to previously calculated atomistic energies of kink-pair configurations, we issue a rule to establish the value of the so-called core width in nonsingular elasticity theories and reduce its arbitrariness as a mathematical construct.

36 MATERIALS SCIENCE↗

Elastic Properties of the Solid Electrolyte Li7La3Zr2O12 (LLZO)

The oxide known as LLZO, with nominal composition Li 7 La 3 Zr 2 O 12 , is a promising solid electrolyte for Li-based batteries due to its high Li-ion conductivity and chemical stability with respect to lithium. Solid electrolytes may also enable the use of metallic Li anodes by serving as a physical barrier that suppresses dendrite initiation and propagation during cycling. Prior linear elasticity models of the Li electrode/solid electrolyte interface suggest that the stability of this interface is highly dependent on the elastic properties of the solid separator. For example, dendritic suppression is predicted to be enhanced as the electrolyte s shear modulus increases. In the present study a combination of first-principles calculations, acoustic impulse excitation measurements, and nanoindentation experiments are used to determine the elastic constants and moduli for highconductivity LLZO compositions based on Al and Ta doping. The calculated and measured isotropic shear moduli are in good agreement and fall within the range of 56-61 GPa. These values are an order of magnitude larger than that for Li metal and far exceed the minimum value ( 8.5 GPa) believed to be necessary to suppress dendrite initiation. These data suggest that LLZO exhibits sufficient stiffness to warrant additional development as a solid electrolyte for Li batteries.

36 MATERIALS SCIENCE↗

Elastic Green’s Function in Anisotropic Bimaterials Considering Interfacial Elasticity

Here, the two-dimensional elastic Green’s function is calculated for a general anisotropic elastic bimaterial containing a line dislocation and a concentrated force while accounting for the interfacial structure by means of a generalized interfacial elasticity paradigm. The introduction of the interface elasticity model gives rise to boundary conditions that are effectively equivalent to those of a weakly bounded interface. The equations of elastic equilibrium are solved by complex variable techniques and the method of analytical continuation. The solution is decomposed into the sum of the Green’s function corresponding to the perfectly bonded interface and a perturbation term corresponding to the complex coupling nature between the interface structure and a line dislocation/concentrated force. Such construct can be implemented into the boundary integral equations and the boundary element method for analysis of nano-layered structures and epitaxial systems where the interface structure plays an important role.

36 MATERIALS SCIENCE↗

An analytical elastic plastic contact model with strain hardening and frictional effects for normal and oblique impacts

Impact between metallic surfaces is a phenomenon that is ubiquitous in the design and analysis of mechanical systems. We found that to model this phenomenon, a new formulation for frictional elastic–plastic contact between two surfaces is developed. The formulation is developed to consider both frictional, oblique contact (of which normal, frictionless contact is a limiting case) and strain hardening effects. The constitutive model for normal contact is developed as two contiguous loading domains: the elastic regime and a transitionary region in which the plastic response of the materials develops and the elastic response abates. For unloading, the constitutive model is based on an elastic process. Moreover, the normal contact model is assumed to only couple one-way with the frictional/tangential contact model, which results in the normal contact model being independent of the frictional effects. Frictional, tangential contact is modeled using a microslip model that is developed to consider the pressure distribution that develops from the elastic–plastic normal contact. This model is validated through comparisons with experimental results reported in the literature, and is demonstrated to be significantly more accurate than 10 other normal contact models and three other tangential contact models found in the literature.

36 MATERIALS SCIENCE↗

Comparison of Multiscale Method of Cells-Based Models for Predicting Elastic Properties of Filament Wound C/C-SiC

Three different multiscale models, based on the method of cells (generalized and high fidelity) micromechanics models were developed and used to predict the elastic properties of C/C-SiC composites. In particular, the following multiscale modeling strategies were employed: Concurrent multiscale modeling of all phases using the generalized method of cells, synergistic (two-way coupling in space) multiscale modeling with the generalized method of cells, and hierarchical (one-way coupling in space) multiscale modeling with the high fidelity generalized method of cells. The three models are validated against data from a hierarchical multiscale finite element model in the literature for a repeating unit cell of C/C-SiC. Furthermore, the multiscale models are used in conjunction with classical lamination theory to predict the stiffness of C/C-SiC plates manufactured via a wet filament winding and liquid silicon infiltration process recently developed by the German Aerospace Institute.

multiscale modeling↗