Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonlinear elasticity”

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 649 records · Page 36

A computer program for the geometrically nonlinear static and dynamic analysis of arbitrarily loaded shells of revolution, theory and users manual

A digital computer program known as SATANS (static and transient analysis, nonlinear, shells) for the geometrically nonlinear static and dynamic response of arbitrarily loaded shells of revolution is presented. Instructions for the preparation of the input data cards and other information necessary for the operation of the program are described in detail and two sample problems are included. The governing partial differential equations are based upon Sanders' nonlinear thin shell theory for the conditions of small strains and moderately small rotations. The governing equations are reduced to uncoupled sets of four linear, second order, partial differential equations in the meridional and time coordinates by expanding the dependent variables in a Fourier sine or cosine series in the circumferential coordinate and treating the nonlinear modal coupling terms as pseudo loads. The derivatives with respect to the meridional coordinate are approximated by central finite differences, and the displacement accelerations are approximated by the implicit Houbolt backward difference scheme with a constant time interval. The boundaries of the shell may be closed, free, fixed, or elastically restrained. The program is coded in the FORTRAN 4 language and is dimensioned to allow a maximum of 10 arbitrary Fourier harmonics and a maximum product of the total number of meridional stations and the total number of Fourier harmonics of 200. The program requires 155,000 bytes of core storage.

Ball, R. E.↗

Second-order microscopic nonlinear susceptibility in a centrosymmetric material: application to imaging valence electron motion

We report measurements of phase-matched nonlinear x-ray and optical sum-frequency generation from single-crystal silicon using sub-resonant 0.95 eV laser pulses and 9.5 keV hard x-ray pulses from the LCLS free-electron laser. The sum-frequency signal appears as energy and momentum sidebands to the elastic Bragg peak. It is proportional to the magnitude squared of the relevant temporal and spatial Fourier components of the optically induced microscopic charges/currents. We measure the first- and second-order sideband to the 220 Bragg peak and find that the efficiency is maximized when the applied field is along the reciprocal lattice vector. For an optical intensity of $\sim10^{12} \text{W}/\text{cm}^2$, we measure peak efficiencies of $3\times 10^{-7}$ and $3\times 10^{-10}$ for the first and second-order sideband respectively (relative to the elastic Bragg peak). The first-order sideband is consistent with induced microscopic currents along the applied electric field (consistent with an isotropic response). The second-order sideband depends nontrivially on the optical field orientation and is consistent with an anisotropic response originating from induced charges along the bonds with C$_{3v}$ site symmetry. The results agree well with first-principles Bloch-Floquet calculations.

FOS: Physical sciences↗

Controlled motion in an elastic world

The flexibility of the drives and structures of controlled motion systems are presented as an obstacle to be overcome in the design of high performance motion systems, particularly manipulator arms. The task and the measure of performance to be applied determine the technology appropriate to overcome this obstacle. Included in the technologies proposed are control algorithms (feedback and feed forward), passive damping enhancement, operational strategies, and structural design. Modeling of the distributed, nonlinear system is difficult, and alternative approaches are discussed. The author presents personal perspectives on the history, status, and future directions in this area.

Book, Wayne J.↗

Effect of Berry phase on nonlinear response of two-dimensional fermions

We develop a theory of nonlinear response to an electric field of two-dimensional (2D) fermions with topologically nontrivial wave functions characterized by the Berry phase Φ n = n π , n = 1 , 2 , ... . In particular, we find that owing to the suppression of backscattering at odd n , Hall field-induced resistance oscillations, which stem from elastic electron transitions between Hall field-tilted Landau levels, are qualitatively distinct from those at even n : Their amplitude decays with the electric field and their extrema are phase shifted by a quarter cycle. The theory unifies the cases of graphene ( n = 1 ) and graphite bilayer ( n = 2 ) with the case of conventional 2D electron gas ( n = 0 ) and suggests another method to probe backscattering in topological 2D systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Working toward a three-dimensional fatigue closure model for surface cracks

The first reliable elastic fracture mechanics solutions for a surface crack in a plate were obtained by Newman and Raju. The authors, both from the Mechanics of Materials Branch at NASA-Langley, used a highly detailed finite element solution requiring substantial computational resources. Computers have since become more powerful and available; however, many important related problems remain computationally expensive. The problem of three-dimensional fatigue crack growth taking into account plasticity-induced crack closure is one such problem. It is the goal of this research to provide an efficient method to account for three-dimensional crack closure in fatigue. Newman developed a two-dimensional plasticity-induced crack closure model for center cracked specimens. This model requires iterations to determine both the contact solution at each growth step and the extent of the plastic zone at the crack tip. A three-dimensional version of this model would require obtaining these nonlinear variables all along the crack front. This model must be efficient so that repeated calculations can be performed during crack growth simulations. The highly versatile line spring model (LSM) with contact, fatigue, and plasticity will form the basis of the closure model. There are several required additions to past work to address the three-dimensional crack closure problem. Initially, these additions will include (1) an improved LSM to more accurately obtain the crack opening displacement, stress intensity factors, and elastic T-stress near the ends of the surface crack; (2) a method to determine the extent of the plastic zone all along the crack front; (3) a method to determine the contact zone given a perfectly plastic layer of material on the crack surfaces; (4) a method to determine the magnitude of the compressive contact stress; and (5) a way to implement the degree of constraint along the curved crack front. During the summer ASEE program an enhanced LSM was developed. A method similar to that of 'strip synthesis' first introduced by Fujimoto was used. Briefly, the crack opening displacements of 'slices' of the surface crack in a direction parallel to the plate surface are considered in addition to the standard LSM approach that makes use of springs obtained from slices perpendicular to the plate surface. This enhancement is necessary so that an accurate three-dimensional representation of quantities such as contact zone size, plastic zone size, stress intensity factors, T-stress, and crack opening displacement can be determined. By combining results of previous investigations with the LSM, the problem of three-dimensional crack closure will be addressed. In addition to a closure model, the enhanced LSM can be used for many other problems including interacting surface cracks and fatigue crack growth of a through crack with a curved crack front.

Joseph, Paul F.↗

A fast Fourier transform-based solver for elastic micropolar composites

This work presents a spectral micromechanical formulation for obtaining the full-field and homogenized response of elastic micropolar composites. The algorithm relies on a coupled set of convolution integral equations for the micropolar strains, where periodic Green’s operators associated with a linear homogeneous reference medium are convolved with functions of the Cauchy and couple stress fields that encode the material’s heterogeneity, as well as any potential material nonlinearity. Such convolution integral equations take an algebraic form in the reciprocal Fourier space that can be solved iteratively. In this vein, the fast Fourier transform (FFT) algorithm is leveraged to accelerate the numerical solution, resulting in a mesh-free formulation in which the periodic unit cell representing the heterogeneous material can be discretized by a regular grid of pixels in two dimensions (or voxels in three dimensions). For verification, the numerical solutions obtained with the micropolar FFT solver are compared with analytical solutions for a matrix with a dilute circular inclusion subjected to plane strain loading. The developed computational framework is then used to study length-scale effects and effective (micropolar) moduli of composites with various topological configurations.

97 MATHEMATICS AND COMPUTING↗

The lattice parameter – composition relationship of the body centered cubic uranium-niobium alloys

In this paper, in-situ time-of-flight neutron diffraction experiments were performed on uranium-niobium alloy with 6 wt% Nb to study the lattice parameter-composition relationship for the body centered cubic (bcc) alloys. Based on lattice-parameter measurements of Nb-supersaturated γs phase over the range of 450 – 790 °C, an improved Vegard's-type relationship was established for determination of Nb concentrations at elevated temperatures. Neutron diffraction data were also collected as a function of aging time over the 450 - 600 °C range when γs phase decomposed isothermally into the mixture of orthorhombic α-U and Nb-rich bcc phase, as well as at room-temperature on the alloys aged ex-situ at 500 °C up to five years. From these measurements, the average Nb concentrations in the bcc phase were determined based on the Rietveld refinements of weight fraction and mass conservation relations (lever rule). Over the 15at% – 78at% range of Nb concentrations that correspond to different experimental aging times, the lattice parameters at constant temperatures exhibit a nonlinear S-shaped variation with Nb concentration, and the associated excess volumes of mixing can be described by a subregular solution model of the Redlich–Kister type of polynomial. Over the full range of composition the S-shaped deviation from Vegard's law can be modeled using a combination of an elastic continuum model and a perturbation to the radii of the solute atoms in the solvent, suggesting that electronic interactions between solute and solvent atoms could play an important role in the compositional dependence of lattice parameter for the γ-phase U-Nb alloys. While Vegard's law is a straightforward and reasonably good approximation for the bcc solid solutions in the U-Nb system, the Nb concentrations determined from the weight fraction refinements of diffraction data provide internally consistent, mass-conserving estimates of solute redistribution for the monotectoid reaction.

36 MATERIALS SCIENCE↗

Minimum weight design of ring and stringer stiffeners for axially compressed cylindrical shells with and without internal pressure

Results of analytical study to determine desirable ring and stringer stiffener parameters and proportions for axially compressed stiffened isotropic cylinders with and without internal pressure are presented. This investigation examines the panel and general instability buckling modes of a stiffened cylindrical shell and from this determines desirable stiffener parameters and proportions. Classical buckling equations are used which retain the important effects of the stiffeners. The results determined by using the simpler classical buckling equations are then spot checked and verified using buckling equations which considered discrete ring stiffeners and nonlinear prebuckling deformations. For both rings and stringers, T-shaped stiffeners are preferable and the effects to stiffener shape are much more pronounced at low or zero values of the internal pressure parameter. Simple analytical expressions are developed and presented which express the stiffener area parameter, the ratio of stiffener area and elastics to shell wall area and elastic modulus, in terms of the cylinder geometry and internal pressure parameter.

Block, D. L.↗

MHOST: An efficient finite element program for inelastic analysis of solids and structures

An efficient finite element program for 3-D inelastic analysis of gas turbine hot section components was constructed and validated. A novel mixed iterative solution strategy is derived from the augmented Hu-Washizu variational principle in order to nodally interpolate coordinates, displacements, deformation, strains, stresses and material properties. A series of increasingly sophisticated material models incorporated in MHOST include elasticity, secant plasticity, infinitesimal and finite deformation plasticity, creep and unified viscoplastic constitutive model proposed by Walker. A library of high performance elements is built into this computer program utilizing the concepts of selective reduced integrations and independent strain interpolations. A family of efficient solution algorithms is implemented in MHOST for linear and nonlinear equation solution including the classical Newton-Raphson, modified, quasi and secant Newton methods with optional line search and the conjugate gradient method.

Nakazawa, S.↗

Fully Coupled Aero-Thermochemical-Elastic Simulations of an Eroding Graphite Nozzle

A multiphysics simulation capability has been developed that incorporates mutual interactions between aerodynamics, structural response from aero/thermal loading, ablation/pyrolysis, heating, and surface-to-surface radiation to perform high-fidelity, fully coupled aerothermoelastic ablation simulations, which to date had been unattainable. The multiphysics framework couples CHAR (a 3-D implicit charring ablator solver), Loci/CHEM (a computational fluid dynamics solver for high-speed chemically reacting flows), and Abaqus (a nonlinear structural dynamics solver) to create a fully coupled aerothermoelastic charring ablative solver. The solvers are tightly coupled in a fully integrated fashion to resolve the effects of the ablation pyrolysis and charring process and chemistry products upon the flow field, the changes in surface geometry due to recession upon the flow field, and thermal-structural analysis of the body from the induced aerodynamic heating from the flow field. The multiphysics framework was successfully demonstrated on a solid rocket motor graphite nozzle erosion application. Comparisons were made with available experimental data that measured the throat erosion during the motor firing. The erosion data is well characterized, as the test rig was equipped with a windowed nozzle section for real-time X-ray radiography diagnostics of the instantaneous throat variations for deducing the instantaneous erosion rates. The nozzle initially undergoes a nozzle contraction due to thermal expansion before ablation effects are able to widen the throat. A series of parameters studies were conducted using the coupled simulation capability to determine the sensitivity of the nozzle erosion to different parameters. The parameter studies included the shape of the nozzle throat (flat versus rounded), the material properties, the effect of the choice of turbulence model, and the inclusion or exclusion of the mechanical thermal expansion. Overall, the predicted results match the experiment very well, and the predictions were able to bound the data within acceptable limits.

Blades, E. L.↗

Nonlinear structural and life analyses of a combustor liner

Three-dimensional, nonlinear, finite element structural analyses were performed for a simulated aircraft combustor liner specimen in order to assess the capability of nonlinear analyses using classical inelastic material models to represent the thermoplastic-creep response of the component. In addition, the computed stress-strain history at the critical location was input into life prediction methods in order to evaluate the ability of these procedures to predict crack initiation life. It is concluded that: (1) elastic analysis is adequate for obtaining strain range and critical location; (2) inelastic analyses did not accurately represent cyclic behavior of materials; and (3) none of the crack initiation life prediction methods were satisfactory.

Kaufman, A.↗

Elastoviscoplastic snap-through behavior of shallow arches subjected to thermomechanical loads

The problem of snap-through buckling of clamped shallow arches under thermomechanical loads is investigated. The analysis is based on nonlinear kinematic relations and nonlinear rate-dependent unified constitutive equations. A finite element approach is employed to predict the, in general, inelastic buckling behavior. The construction material is alloy B1900 + Hf, which is commonly utilized in high-temperature environments. The effect of several parameters is assessed. These parameters include the rise parameter and temperature. Comparison between elastic and elastoviscoplastic responses is also presented.

Simitses, George J.↗

Internal strain analysis of ceramics using scanning laser acoustic microscopy

Quantitative studies of material behavior characteristics are essential for predicting the functionality of a material under its operating conditions. A nonintrusive methodology for measuring the in situ strain of small dimeter (to 11 microns) ceramic fibers under uniaxial tensile loading and the local internal strains of ceramics and ceramic composites under flexural loading is introduced. The strain measurements and experimentally observed mechanical behavior are analyzed in terms of the microstructural development and fracture behavior of each test specimen evaluated. Measurement and analysis of Nicalon silicon carbide (SiC) fiber (15 microns diameter) indicate that the mean elastic modulus of the individual fiber is 185.3 GPa. Deviations observed in the experimentally determined elastic modulus values between specimens were attributed to microstructural variations which occur during processing. Corresponding variations in the fracture surface morphology were also observed. The observed local mechanical behavior of a lithium alumino-silicate (LAS) glass ceramic, a LAS/SiC monofilament composite, and a calcium alumino-silicate (CAS)/SiC fully reinforced composite exhibits nonlinearities and apparent hysteresis due to the subcritical mechanical loading. Local hysteresis in the LAS matrices coincided with the occurrence of multiple fracture initiation sites, localized microcracking, and secondary cracking. The observed microcracking phenomenon was attributed to stress relaxation of residual stresses developed during processing, and local interaction of the crack front with the microstructure. The relaxation strain and stress predicted on apparent mechanical hysteresis effects were defined and correlated with the magnitude of the measured fracture stress for each specimen studied. This quantitative correlation indicated a repeatable measure of the stress at which matrix microcracking occurred for stress relief of each material system. Stress relaxation occurred prior to the onset of steady state cracking conditions. The relaxation stress occurred at 18.5 percent of the fracture stress in LAS and 11.0 percent of the yield stress in CAS/SiC. The relaxation stress ratio was dependent upon the dominant fracture mode of the LAS/SiC specimens. Relaxation stress ratios greater than 0.30 were observed for specimens which fractured due to shear at the fiber matrix interface; specimens which fracture due to tensile cracking had relaxation stress ratios less than 0.30. The stress relaxation ratio appeared to be a specific characteristic of the glass ceramic material. The measured stress relaxation for LAS indicated a measure of the inherent residual stresses in the material due to processing and suggested localized toughening mechanisms for brittle material structures.

Kent, Renee M.↗

Thermo-elasto-viscoplastic analysis of problems in extension and shear

The problems of extension and shear behavior of structural elements made of carbon steel and subjected to large thermomechanical loads are investigated. The analysis is based on nonlinear geometric and constitutive relations, and is expressed in a rate form. The material constitutive equations are capable of reproducing all nonisothermal, elasto-viscoplastic characteristics. The results of the test problems show that: (1) the formulation can accommodate very large strains and rotations; (2) the model incorporates the simplification associated with rate-insensitive elastic response without losing the ability to model a rate-temperature dependent yield strength and plasticity; and (3) the formulation does not display oscillatory behavior in the stresses for the simple shear problem.

Riff, R.↗

Thermo-elasto-viscoplastic analysis of problems in extension and shear

The problems of extension and shear behavior of structural elements made of carbon steel and subjected to large thermomechanical loads are investigated. The analysis is based on nonlinear geometric and constitutive relations, and is expressed in a rate form. The material constitutive equations are capable of reproducing all nonisothermal, elasto-viscoplastic characteristics. The results of the test problems show that: (1) the formulation can accommodate very large strains and rotations; (2) the model incorporates the simplification associated with rate-insensitive elastic response without losing the ability to model a rate-temperature dependent yield strength and plasticity; and (3) the formulation does not display oscillatory behavior in the stresses for the simple shear problem.

Riff, R.↗

Spacecraft dynamics characterization and control system failure detection, volume 1

The work under this grant has been directed to two aspects of the control of flexible spacecraft: (1) the modeling of deployed or erected structures including nonlinear joint characteristics; and (2) the detection and isolation of failures of the components of control systems for large space structures. The motivation for the first of these research tasks is the fact that very large assemblies in space will have to be built or deployed in situ. A likely scenario is, in fact, a combination of these wherein modules which are folded for transportation into orbit are erected to their final configuration and then jointed with other such erected modules to form the full assembly. Any such erectable modules will have joints. It remains to be seen whether or not joints designed for operational assemblies will have nonlinear properties, but it seems prudent to develop a methodology for dealing with that possibility. The motivation for the second of these research tasks is the fact that we foresee large assemblies in space which will require active control to damp vibrations and/or hold a desired shape. Lightweight structures will be very flexible, with many elastic modes having very low frequencies. In order to control these modes well, the control system will likely require many sensors and many actuators, probably distributed over much of the structure. The combination of a large number of control system components with long operational periods virtually guarantees that these systems will suffer control system component failures during operation. The control system must be designed to tolerate failures of some sensors and actuators, and still be able to continue to perform its function.

Source record↗

Stimulated resonant inelastic X-ray scattering in a solid

Abstract When materials are exposed to X-ray pulses with sufficiently high intensity, various nonlinear effects can occur. The most fundamental one consists of stimulated electronic decays after resonant absorption of X-rays. Such stimulated decays enhance the number of emitted photons and the emission direction is confined to that of the stimulating incident photons which clone themselves in the process. Here we report the observation of stimulated resonant elastic (REXS) and inelastic (RIXS) X-ray scattering near the cobalt L 3 edge in solid Co/Pd multilayer samples. We observe an enhancement of order 10 6 of the stimulated over the conventional spontaneous RIXS signal into the small acceptance angle of the RIXS spectrometer. We also find that in solids both stimulated REXS and RIXS spectra contain contributions from inelastic electron scattering processes, even for ultrashort 5 fs pulses. Our results reveal the potential and caveats of the development of stimulated RIXS in condensed matter.

36 MATERIALS SCIENCE↗

Debye mechanism of giant microwave absorption in superconductors

Here, we discuss a novel mechanism of microwave absorption in superconductors, which is similar to the Debye absorption mechanism in molecular gases. The contribution of this mechanism to the ac conductivity is proportional to the inelastic quasiparticle relaxation time $τ_{\text{in}}$ rather than the elastic one $τ_{\text{el}}$ and therefore it can be much larger than the conventional one. The Debye contribution to the linear conductivity arises only in the presence of a dc supercurrent in the system and its magnitude depends strongly on the orientation of the microwave field relative to the supercurrent. The Debye contribution to the nonlinear conductivity exists even in the absence of dc supercurrent, and it is proportional to $τ_{\text{in}}$. Therefore the nonlinear threshold is anomalously low. Microwave absorption measurements may provide direct information about $τ_{\text{in}}$ in superconductors.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗