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 757 records · Page 42

Recurrent neural network-based multiaxial plasticity model with regularization for physics-informed constraints

We report a recurrent neural network (RNN) based model is developed as a surrogate to predict nonlinear plastic response under multiaxial loading. The RNN-based model is trained and tested on stress versus strain curves generated using a numerical solution based on the classical radial return method. Besides simply learning the basic constitutive relationship, a novel approach is taken to enforce certain physical conditions. Specifically, regularization is employed to maintain non-negative plastic power density throughout the loading history thereby ensuring monotonically increasing plastic work and thermodynamic consistency. Enforcing physics in this manner permits coupling of the data-driven RNN approach with physics-based knowledge and laws. This has the effect of reducing the necessary amount of data and ensuring known physical laws are not violated. Since, once trained, the model need not perform the expensive task of solving nonlinear equations, its efficiency is orders of magnitude greater than its numerical counterpart. The RNN-based model has been trained on varied sets of data and the accuracy on test datasets validated. The developed model is general and robust and has widespread application such as in the simulation of metal forming, large scale plasticity, and part life prediction.

42 ENGINEERING↗

Solar Cruiser TRAC Boom Development

The TRAC boom studied is a 30 meter, high-aspect ratio and highly nonlinear structural element planned for the NASA Solar Cruiser spacecraft where it would serve as the solar sail’s skeletal system designed to deploy the packaged reflective sail membrane as well as support operational sail tension loads for flight. Key design drivers for this boom are understanding and characterizing load-deformation behavior and quantifying buckling performance–under flight conditions (zero gravity). Knowing that the full-scale solar sail system cannot be feasibly tested for flight conditions under Earth’s gravity, a building block approach is methodically followed to verify finite element analysis (FEA) predictions and to evolve the finite element model (FEM) accordingly, beginning with the lowest level complexity to match global elastic behavior to the highest level complexity to bound theoretical responses under conservative expected flight conditions. Although the testing campaign for the TRAC boom has not developed yet to produce a full-scale test data point that loads the boom to failure, the fidelity of simulations have evolved with success and warrant further interest and work on this high-aspect ratio, but highly stable, deployable structural boom.

Lee Nguyen↗

A quantitative assessment of the model form error of friction models across different interface representations for jointed structures

Hysteretic models are widely used to model frictional interactions in joints to recreate experimental behavior. However, it is unclear which models are best suited for fitting or predicting the responses of structures. Here, the present study evaluates 26 friction model/interface representation combinations to quantify the model form error. A Quasi-Static Modal Analysis approach (termed Rayleigh Quotient Nonlinear Modal Analysis) is adopted to calculate the nonlinear system response, and a Multi-Objective Optimization is solved to fit experimental data of the first mode of the Brake-Reuß Beam. Optimized parameters from the first mode are applied to the second and third bending modes to quantify the predictive ability of the models. Formulations for both tracing full hysteresis loops and recreating hysteresis loops from a single loading curve (Masing assumptions) are considered. Smoothly varying models applied to a five patch representation showed the highest flexibility (for fitting mode 1) and good predictive potential (for modes 2 and 3). For a second formulation, which uses 152 frictional elements to represent the interface, the physically motivated spring in series with a Coulomb slip model (elastic dry friction) has high error for fitting mode 1 and performs near the middle for predicting higher modes. For both interface representation, the best fit models are not the most physical, but rather the ones with the most parameters (as expected); however, the more physical models perform somewhat better for predicting the higher modes.

42 ENGINEERING↗

Oceanic lithosphere and asthenosphere - Thermal and mechanical structure

A coupled thermomechanical subsolidus model of the oceanic lithosphere and asthenosphere is developed which includes vertical heat conduction, a temperature-dependent thermal conductivity, heat advection by a horizontal and vertical mass flow that depends on depth and age, contributions of viscous dissipation or shear heating, a linear or nonlinear deformation law relating shear stress and strain rate, as well as a temperature- and pressure-dependent viscosity. The model requires a constant horizontal velocity and temperature at the surface, but zero horizontal velocity and constant temperature at great depths. The depth- and age-dependent temperature, horizontal and vertical velocities, and viscosity structure of the lithosphere and asthenosphere are determined along with the age-dependent shear stress in those two zones. The ocean-floor topography, oceanic heat flow, and lithosphere thickness are deduced as functions of ocean-floor age; seismic velocity profiles which exhibit a marked low-velocity zone are constructed from the age-dependent geotherms and assumed values of the elastic parameters. It is found that simple boundary-layer cooling determines the thermal structure at young ages, while effects of viscous dissipation become more important at older ages.

Schubert, G.↗

Quasi-static extension of a tensile crack contained in a viscoelastic-plastic solid

Final stretch criterion of failure is applied to the problem of quasi-static extension of a crack embedded in an elastic-plastic or viscoelastic-plastic matrix. The slow growth under subcritical conditions in a rate-sensitive Tresca solid is shown to be a superposition of creep rupture and McClintock's ductile growth. This type of growth occurs at subcritical magnitude of the imposed K-factor and can be accounted for only through a recognition of inelastic properties of solids. In the subcritical range there is no unique value for K sub c independent of geometrical configuration and flaw size. Not only the produced states of stress and strain are dependent on the loading path, but also the material resistance to fracture turns out to be a function of the history of loading that precedes catastrophic failure. A nonlinear integro-differential equation of motion is derived for a crack progressing through a viscoelastic medium with some limited ability to plastic flow. Examples of numerical integration are given incorporating both monotonic and cyclic loading programs.

Wnuk, M. P.↗

Robust multivariable controller design for flexible spacecraft

Large, flexible spacecraft are typically characterized by a large number of significant elastic modes with very small inherent damping, low, closely spaced natural frequencies, and the lack of accurate knowledge of the structural parameters. Summarized here is some recent research on the design of robust controllers for such spacecraft, which will maintain stability, and possible performance, despite these problems. Two types of controllers are considered, the first being the linear-quadratic-Gaussian-(LQG)-type. The second type utilizes output feedback using collocated sensors and actuators. The problem of designing robust LQG-type controllers using the frequency domain loop transfer recovery (LTR) method is considered, and the method is applied to a large antenna model. Analytical results regarding the regions of stability for LQG-type controllers in the presence of actuator nonlinearities are also presented. The results obtained for the large antenna indicate that the LQG/LTR method is a promising approach for control system design for flexible spacecraft. For the second type of controllers (collocated controllers), it is proved that the stability is maintained in the presence of certain commonly encountered nonlinearities and first-order actuator dynamics. These results indicate that collocated controllers are good candidates for robust control in situations where model errors are large.

Joshi, Suresh M.↗

A Poromechanical Framework for Internal Interactions Induced by Solid Inclusions

The framework of poromechanics is generalized to simulate the multiscale behavior of porous media subjected to internal loadings stemming from the growth of solid inclusions. This generalization is designed to enable the study of anisotropic internal stress generation from solid growth within the pores, while recovering isotropic fluid-induced loading as a particular case. For this purpose, a mathematical strategy to define constitutive tensors in a thermodynamically consistent form is proposed, thus offering new opportunities for determining the poromechanical properties of a porous solid through advanced experimentation or micromechanical models. The framework is specialized by means of established elastic solutions for single pore–matrix interaction, as well as through homogenization schemes considering the interaction among congruent pores. In particular, the second Eshelby solution and the Tanaka–Mori–Benveniste homogenization scheme are used to derive a microporoelastic model. At an elemental scale, the model is tested under mixed control conditions by replicating different scenarios of geomaterial testing. In addition, the model characteristics are outlined with reference to inelastic microscopic loadings replicating chemo-mechanical forcing, such as expansive crystal formation. Through a series of parametric analyses, it is shown that the microstructure of the pores significantly influences the properties of porous media. Most notably, it is shown that the effects of a solid forming within the pores depend in a highly nonlinear fashion on the constitutive characteristics of the inhomogeneities and can therefore not be readily quantified or predicted without models capturing the diverse multiscale interactions among pores, inhomogeneities, and matrix.

chemomechanics↗

The influence of Al2O3 on the structural properties of MgSiO3 akimotoite

Abstract Akimotoite, a MgSiO3 polymorph present in the lower transition zone within ultramafic portions of subducting slabs and potentially also in the ambient mantle, will partition some amount of Al, raising the question of how this will affect its crystal structure and properties. In this study, a series of samples along the MgSiO3-Al2O3 (akimotoite-corundum) solid solution have been investigated by means of single-crystal X-ray diffraction to examine their crystal chemistry. Results show a strong nonlinear behavior of the a- and c-axes as a function of Al content, which arises from fundamentally different accommodation mechanisms in the akimotoite and corundum structures. Furthermore, two Al2O3-bearing akimotoite samples were investigated at high pressure to determine the different compression mechanisms associated with Al substitution. Al2O3-bearing akimotoite becomes more compressible at least up to 20 mol% Al2O3, due likely to an increase in compressibility as the Al cation is incorporated into the SiO6 octahedron. This observation is in strong contrast to the stiffer corundum end-member having a KT = 250 GPa, which is larger than that of the akimotoite end-member [KT = 205(1) GPa]. These findings have implications for mineral physics models of elastic properties, which have in the past assumed linear mixing behavior between the MgSiO3 akimotoite and Al2O3 corundum end-members to calculate sound wave velocities for Al-bearing akimotoite at high pressure and temperature.

Geochemistry & Geophysics↗

Finite-strain large-deflection elastic-viscoplastic finite-element transient response analysis of structures

A method of analysis for thin structures that incorporates finite strain, elastic-plastic, strain hardening, time dependent material behavior implemented with respect to a fixed configuration and is consistently valid for finite strains and finite rotations is developed. The theory is formulated systematically in a body fixed system of convected coordinates with materially embedded vectors that deform in common with continuum. Tensors are considered as linear vector functions and use is made of the dyadic representation. The kinematics of a deformable continuum is treated in detail, carefully defining precisely all quantities necessary for the analysis. The finite strain theory developed gives much better predictions and agreement with experiment than does the traditional small strain theory, and at practically no additional cost. This represents a very significant advance in the capability for the reliable prediction of nonlinear transient structural responses, including the reliable prediction of strains large enough to produce ductile metal rupture.

Rodal, J. J. A.↗

Drops in Space: Super Oscillations and Surfactant Studies

An unprecedented microgravity observation of maximal shape oscillations of a surfactant-bearing water drop the size of a ping pong ball was observed during a mission of Space Shuttle Columbia as part of the second United States Microgravity Laboratory-USML-2 (STS-73, October 20-November 5, 1995). The observation was precipitated by the action of an intense sound field which produced a deforming force on the drop. When this deforming force was suddenly reduced, the drop executed nearly free and axisymmetric oscillations for several cycles, demonstrating a remarkable amplitude of nonlinear motion. Whether arising from the discussion of modes of oscillation of the atomic nucleus, or the explosion of stars, or how rain forms, the complex processes influencing the motion, fission, and coalescence of drops have fascinated scientists for centuries. Therefore, the axisymmetric oscillations of a maximally deformed liquid drop are noteworthy, not only for their scientific value but also for their aesthetic character. Scientists from Yale University, the Jet Propulsion Laboratory (JPL) and Vanderbilt University conducted liquid drop experiments in microgravity using the acoustic positioning/manipulation environment of the Drop Physics Module (DPM). The Yale/JPL group's objectives were to study the rheological properties of liquid drop surfaces on which are adsorbed surfactant molecules, and to infer surface properties such as surface tension, Gibb's elasticity, and surface dilatational viscosity by using a theory which relies on spherical symmetry to solve the momentum and mass transport equations.

Apfel, Robert E.↗

NASA Tech Briefs, December 2005

Topics covered include: Video Mosaicking for Inspection of Gas Pipelines; Shuttle-Data-Tape XML Translator; Highly Reliable, High-Speed, Unidirectional Serial Data Links; Data-Analysis System for Entry, Descent, and Landing; Hybrid UV Imager Containing Face-Up AlGaN/GaN Photodiodes; Multiple Embedded Processors for Fault-Tolerant Computing; Hybrid Power Management; Magnetometer Based on Optoelectronic Microwave Oscillator; Program Predicts Time Courses of Human/ Computer Interactions; Chimera Grid Tools; Astronomer's Proposal Tool; Conservative Patch Algorithm and Mesh Sequencing for PAB3D; Fitting Nonlinear Curves by Use of Optimization Techniques; Tool for Viewing Faults Under Terrain; Automated Synthesis of Long Communication Delays for Testing; Solving Nonlinear Euler Equations With Arbitrary Accuracy; Self-Organizing-Map Program for Analyzing Multivariate Data; Tool for Sizing Analysis of the Advanced Life Support System; Control Software for a High-Performance Telerobot; Java Radar Analysis Tool; Architecture for Verifiable Software; Tool for Ranking Research Options; Enhanced, Partially Redundant Emergency Notification System; Close-Call Action Log Form; Task Description Language; Improved Small-Particle Powders for Plasma Spraying; Bonding-Compatible Corrosion Inhibitor for Rinsing Metals; Wipes, Coatings, and Patches for Detecting Hydrazines; Rotating Vessels for Growing Protein Crystals; Oscillating-Linear-Drive Vacuum Compressor for CO2; Mechanically Biased, Hinged Pairs of Piezoelectric Benders; Apparatus for Precise Indium-Bump Bonding of Microchips; Radiation Dosimetry via Automated Fluorescence Microscopy; Multistage Magnetic Separator of Cells and Proteins; Elastic-Tether Suits for Artificial Gravity and Exercise; Multichannel Brain-Signal-Amplifying and Digitizing System; Ester-Based Electrolytes for Low-Temperature Li-Ion Cells; Hygrometer for Detecting Water in Partially Enclosed Volumes; Radio-Frequency Plasma Cleaning of a Penning Malmberg Trap; Reduction of Flap Side Edge Noise - the Blowing Flap; and Preventing Accidental Ignition of Upper-Stage Rocket Motors.

Source record↗

Nonlinear analysis for the response and failure of compression-loaded angle-ply laminates with a hole

The objective was to determine the effect of nonlinear material behavior on the response and failure of unnotched and notched angle-ply laminates under uniaxial compressive loading. The endochronic theory was chosen as the constitutive theory to model the AS4/3502 graphite-epoxy material system. Three-dimensional finite element analysis incorporating the endochronic theory was used to determine the stresses and strains in the laminates. An incremental/iterative initial strain algorithm was used in the finite element program. To increase computational efficiency, a 180 deg rotational symmetry relationship was utilized and the finite element program was vectorized to run on a supercomputer. Laminate response was compared to experimentation revealing excellent agreement for both the unnotched and notched angle-ply laminates. Predicted stresses in the region of the hole were examined and are presented, comparing linear elastic analysis to the inelastic endochronic theory analysis. A failure analysis of the unnotched and notched laminates was performed using the quadratic tensor polynomial. Predicted fracture loads compared well with experimentation for the unnotched laminates, but were very conservative in comparison with experiments for the notched laminates.

Mathison, Steven R.↗

A non-isothermal breakage-damage model for plastic-bonded granular materials incorporating temperature, pressure, and rate dependencies

Plastic-bonded granular materials (PBM) are widely used in industrial sectors, including building construction, abrasive applications, and defense applications such as plastic-bonded explosives. The mechanical behavior of PBM is highly nonlinear, irreversible, rate dependent, and temperature sensitive governed by various micromechanical attributions such as grain crushing and binder damage. This paper presents a thermodynamically consistent, microstructure-informed constitutive model to capture these characteristic behaviors of PBM. Key features of the model include a breakage internal variable to upscale the grain-scale information to the continuum level and to predict grain size evolution under mechanical loading. In addition, a damage internal state variable is introduced to account for the damage, deterioration, and debonding of the binder matrix upon loading. Temperature is taken as a fundamental external state variable to handle non-isothermal loading paths. The proposed model is able to capture with good accuracy several important aspects of the mechanical properties of PBM, such as pressure-dependent elasticity, pressure-dependent yield strength, brittle-to-ductile transition, temperature dependency, and rate dependency in the post-yielding regime. Furthermore, the model is validated against multiple published datasets obtained from confined and unconfined compression tests, covering various PBM compositions, confining pressures, temperatures, and strain rates.

Breakage↗

Demonstration of Probabilistic Sensitivity Analyses Tools on the Structural Response of a Representative Inflatable Space Structure

This work provides an initial step toward demonstrating a probabilistic numerical simulation capability to support trade studies and the development of a certification plan for inflatable space habitats. This study concentrates on interpreting the results from probabilistic analysis and numerical simulation tools to identify parameter sensitivities for a novel inflatable airlock concept, specifically the Non‐Axisymmetric Inflatable Pressure Structure (NAIPS) that was designed and tested under NASA's Minimalistic Advanced Softgoods Hatch (MASH) Program. A brief overview of the finite element model is provided along with the probabilistic sensitivity analysis approach. The sensitivity studies required a model that was numerically stable and efficient enough that hundreds of simulations could be completed in the allotted time. Therefore, the existing full model was simplified by: extracting a quarter symmetry section of the dome; focusing on a single inflation pressure; and replacing the non‐linear material stress‐strain curves with linear, isotropic materials defined by elastic moduli. Responses of interest include the sensitivity of various structural component loads to material properties, cord lengths, inflation pressure and friction. Multiple sensitivity studies were completed and three are reported here. The first study focused on utilizing wide input parameter ranges to provide an opportunity to assess numerical robustness. The next two studies narrowed the parameter ranges to enable focus on understanding uncertainty at a fixed operating condition. The completion of the sensitivity studies improved understanding of the interdependence of multiple inputs on the responses. In addition, numerical stability of the simulations over wide parameter ranges, shows the feasibility of incorporating uncertainty‐based methods in the design and certification of inflatable space habitats. With the experience and trust gained, it is anticipated that these same methods will be applied to nonlinear, orthotropic models in the future.

Lyle, Karen H.↗

A general-purpose approach to computer-aided dynamic analysis of a flexible helicopter

A general purpose mathematical formulation is described for dynamic analysis of a helicopter consisting of flexible and/or rigid bodies that undergo large translations and rotations. Rigid body and elastic sets of generalized coordinates are used. The rigid body coordinates define the location and the orientation of a body coordinate frame (global frame) with respect to an inertial frame. The elastic coordinates are introduced using a finite element approach in order to model flexible components. The compatibility conditions between two adjacent elements in a flexible body are imposed using a Boolean matrix, whereas the compatibility conditions between two adjacent bodies are imposed using the Lagrange multiplier approach. Since the form of the constraint equations depends upon the type of kinematic joint and involves only the generalized coordinates of the two participating elements, then a library of constraint elements can be developed to impose the kinematic constraint in an automated fashion. For the body constraints, the Lagrange multipliers yield the reaction forces and torques of the bodies at the joints. The virtual work approach is used to derive the equations of motion, which are a system of differential and algebraic equations that are highly nonlinear. The formulation presented is general and is compared with hard-wired formulations commonly used in helicopter analysis.

Agrawal, Om P.↗

Buckling Testing and Analysis of Honeycomb Sandwich Panel Arc Segments of a Full-Scale Fairing Barrel: Comparison of In- and Out-of-Autoclave Facesheet Configurations

Four honeycomb sandwich panels, representing 1/16th arc segments of a 10-m diameter barrel section of the Heavy Lift Launch Vehicle, were manufactured and tested under the NASA Composites for Exploration and the NASA Constellation Ares V programs. Two configurations were chosen for the panels: 6-ply facesheets with 1.125 in. honeycomb core and 8-ply facesheets with 1.0 in. honeycomb core. Additionally, two separate carbon fiber/epoxy material systems were chosen for the facesheets: in-autoclave IM7/977-3 and out-of-autoclave T40-800b/5320-1. Smaller 3 ft. by 5 ft. panels were cut from the 1/16th barrel sections and tested under compressive loading. Furthermore, linear eigenvalue and geometrically nonlinear finite element analyses were performed to predict the compressive response of each 3 ft. by 5 ft. panel. To improve the robustness of the geometrically nonlinear finite element model, measured surface imperfections were included in the geometry of the model. Both the linear and nonlinear models yielded good qualitative and quantitative predictions. Additionally, it was correctly predicted that the panel would fail in buckling prior to failing in strength. Furthermore, several imperfection studies were performed to investigate the influence of geometric imperfections, fiber angle misalignments, and three-dimensional effects on the compressive response of the panel.

epoxy matrix composites↗

NUMERICAL CALCULATION OF LOSSES OF TRAPPED VORTICES UNDER STRONG RF MEISSNER CURRENT AND DC SUPERHEATING FIELD IN TYPE II SUPERCONDUCTORS

Research on the vortex dynamics and enhancing of superheating field in superconductors has attracted much attention in accelerator physics community to develop next-generation high-performance accelerator cavities. However, the extreme dynamics of curvilinear elastic vortices driven by very strong currents close to the depairing limit or superheating field of a superconductor with a nanostructured surface has not been well understood. We calculated the superheating field Hsh and critical momentum kc characterizing the wavelength of the instability ?m of the Meissner state to flux penetration by solving numerically the Ginzburg?Landau equations. A bulk superconductor, superconductor with the inhomogeneous surface disorder (S-S), and multilayered surface (S-I-S) have been thoroughly investigated in this work. Our result showed that S-S and S-I-S structures can enhance the superheating field well above their clean limit. In this work extensive numerical simulation of the power dissipated by an oscillating vortex segment driven by the surface ac Meissner currents was performed. Our simulations take into account the nonlinear vortex line tension, vortex mass, Bardeen?Stephen viscous vortex drag applicable at low fields, and nonlinear Larkin-Ovchinnikov (LO) viscous drag coefficient ?(v) at high fields and pinning force. We showed that the LO decrease of ?(v) with the vortex velocity v could radically change the field dependence of the surface resistance Ri(H) caused by trapped vortices. At low frequencies Ri(H) exhibits a conventional increase with H. However, as frequency increases, the surface resistance becomes a nonmonotonic function of H which decreases with H at higher fields irrespective of the pinning distribution. Overheating can mask the descending field dependence of Ri(H) as frequency increases. Our numerical simulations also show that the LO effect can cause a vortex bending instability at high field amplitudes and frequencies, giving rise to the formation of dynamic kinks along with the vortex when a vortex is pinned strongly to one end. Nonlinear losses of trapped vortices in thick films under high-amplitude RF fields as functions of frequency, mean free path and pinning characteristics have been calculated.

Pathirana, Walive↗

Analysis of plane-plastic stress problems with axial symmetry in strain-hardening range

A simple method is developed for solving plane-plastic-stress problems with axial symmetry in the strain-hardening range which is based on the deformation theory of plasticity employing the finite-strain concept. The equations defining the problems are first reduced to two simultaneous nonlinear differential equations involving two dependent variables: (a) the octahedral shear strain, and (b) a parameter indicating the ratio of principal stresses. By multiplying the load and dividing the radius by an arbitrary constant, it is possible to solve these problems without iteration for any value of the modified load. The constant is determined by the boundary condition. This method is applied to a circular membrane under pressure, a rotating disk without and with a central hole, and an infinite plate with a circular hole. Two materials, inconel x and 16-25-6, the octahedral shear stress-strain relations of which do not follow the power law, are used. Distributions of octahedral shear strain, as well as of principal stresses and strains, are obtained. These results are compared with the results of the same problems in the elastic range.

Wu, M H Lee↗