Search NASA⌕ Search

SEARCH · Search NASA

Results for “finite deformation”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
NASA NTRS2011

A Deployable 4 Meter 180 to 680 GHz Antenna for the Scanning Microwave Limb Sounder

The Scanning Microwave Limb Sounder (SMLS) is a space-borne heterodyne radiometer which will measure pressure, temperature and atmospheric constituents from thermal emission between 180 and 680 GHz. SMLS, planned for the Global Atmospheric Composition Mission of the NRC Decadal Survey, uses a novel toric Cassegrain antenna to perform both elevation and azimuth scanning. These provide better horizontal and temporal resolution and coverage than were possible with elevation-only scanning at typical Low-Earth orbit spacing in the two previous MLS satellite instruments. Development of the SMLS antenna was the focus of a 2006 Small Business Innovative Research (SBIR) program whose phase II culminated in the fabrication and thermal stability testing of a composite demonstration model of the SMLS primary reflector. This reflector has the full 4m height and 1/3 the width planned for flight. An Instrument Incubator Program (IIP) titled "A deployable 4 Meter 180 to 680 GHz antenna for the Scanning Microwave Limb Sounder" continues development of the SMLS antenna with the study of 5 topics: 1) detailed mathematical modeling of the antenna patterns from which we simulate geophysical parameter retrievals in order to establish FOV performance requirements; 2) thorough correlation of finite element model predictions with measurements made on the SBIR reflector. We will again measure deformations of this reflector, under more flight-like thermal gradients, using higher precision metrology techniques available in a new large-aperture facility at JPL; 3) fabrication of a full-width primary reflector whose asbuilt surface figure will better meet the figure requirements of SMLS than did the SBIR reflector; 4) integration of the primary with other reflectors, and with residual front ends built in a 2007 IIP, in a breadboard antenna; and finally 5) RF testing of the breadboard on a Near Field Range at JPL. We report on significant progress in 3 areas of the current IIP: development of the mathematical model to predict SMLS antenna patterns and their application in a preliminary set of geophysical retrievals; the correlation between surface deformation predicted by finite element models and measurement in the 2009 isothermal stability tests of the SBIR, with implications for the thermal gradient testing to be performed at JPL; and aspects of the conceptual design of the full-width primary reflector to be fabricated and tested in the 2nd and 3rd years of the IIP.

heterodyne radiometer↗
NASA NTRS2009

Modeling and Simulation of Variable Mass, Flexible Structures

The advent of the new Ares I launch vehicle has highlighted the need for advanced dynamic analysis tools for variable mass, flexible structures. This system is composed of interconnected flexible stages or components undergoing rapid mass depletion through the consumption of solid or liquid propellant. In addition to large rigid body configuration changes, the system simultaneously experiences elastic deformations. In most applications, the elastic deformations are compatible with linear strain-displacement relationships and are typically modeled using the assumed modes technique. The deformation of the system is approximated through the linear combination of the products of spatial shape functions and generalized time coordinates. Spatial shape functions are traditionally composed of normal mode shapes of the system or even constraint modes and static deformations derived from finite element models of the system. Equations of motion for systems undergoing coupled large rigid body motion and elastic deformation have previously been derived through a number of techniques [1]. However, in these derivations, the mode shapes or spatial shape functions of the system components were considered constant. But with the Ares I vehicle, the structural characteristics of the system are changing with the mass of the system. Previous approaches to solving this problem involve periodic updates to the spatial shape functions or interpolation between shape functions based on system mass or elapsed mission time. These solutions often introduce misleading or even unstable numerical transients into the system. Plus, interpolation on a shape function is not intuitive. This paper presents an approach in which the shape functions are held constant and operate on the changing mass and stiffness matrices of the vehicle components. Each vehicle stage or component finite element model is broken into dry structure and propellant models. A library of propellant models is used to describe the distribution of mass in the fuel tank or Solid Rocket Booster (SRB) case for various propellant levels. Based on the mass consumed by the liquid engine or SRB, the appropriate propellant model is coupled with the dry structure model for the stage. Then using vehicle configuration data, the integrated vehicle model is assembled and operated on by the constant system shape functions. The system mode shapes and frequencies can then be computed from the resulting generalized mass and stiffness matrices for that mass configuration. The rigid body mass properties of the vehicle are derived from the integrated vehicle model. The coupling terms between the vehicle rigid body motion and elastic deformation are also updated from the constant system shape functions and the integrated vehicle model. This approach was first used to analyze variable mass spinning beams and then prototyped into a generic dynamics simulation engine. The resulting code was tested against Crew Launch Vehicle (CLV-)class problems worked in the TREETOPS simulation package and by Wilson [2]. The Ares I System Integration Laboratory (SIL) is currently being developed at the Marshall Space Flight Center (MSFC) to test vehicle avionics hardware and software in a hardware-in-the-loop (HWIL) environment and certify that the integrated system is prepared for flight. The Ares I SIL utilizes the Ares Real-Time Environment for Modeling, Integration, and Simulation (ARTEMIS) tool to simulate the launch vehicle and stimulate avionics hardware. Due to the presence of vehicle control system filters and the thrust oscillation suppression system, which are tuned to the structural characteristics of the vehicle, ARTEMIS must incorporate accurate structural models of the Ares I launch vehicle. The ARTEMIS core dynamics simulation models the highly coupled nature of the vehicle flexible body dynamics, propellant slosh, and vehicle nozzle inertia effects combined with mass and flexible body properties that vary significant with time during the flight. All forces that act on the vehicle during flight must be simulated, including deflected engine thrust force, spatially distributed aerodynamic forces, gravity, and reaction control jet thrust forces. These forces are used to excite an integrated flexible vehicle, slosh, and nozzle dynamics model for the vehicle stack that simulates large rigid body translations and rotations along with small elastic deformations. Highly effective matrix math operations on a distributed, threaded high-performance simulation node allow ARTEMIS to retain up to 30 modes of flex for real-time simulation. Stage elements that separate from the stack during flight are propagated as independent rigid six degrees of freedom (6DOF) bodies. This paper will present the formulation of the resulting equations of motion, solutions to example problems, and describe the resulting dynamics simulation engine within ARTEMIS.

Tobbe, Patrick A.↗
NASA NTRS1983

Analyses of large quasistatic deformations of inelastic bodies by a new hybrid-stress finite element algorithm - Applications

A new hybrid-stress finite element algorithm suitable for analyzing large quasistatic deformations of inelastic solids is presented and its feasibility and performance are demonstrated with examples. The algorithm provides extremely accurate bifurcation analysis which is stable with respect to variation in the finite element mesh, so long as the same type of element is used in every mesh. When the mesh element is varied, the result changes in a predictable manner. The method does not necessarily lead to an upper or lower bound for the critical load. An explicit forward gradient scheme is used to improve stability and is shown to be useful also for elongation-dominated deformations. The application of the method to the onset of necking in plane extension and to deformation and stress in plane extension of an elasticoviscous fluid with an array of cylindrical voids is given in detail.

Reed, K. W.↗
NASA NTRS2015

Face Sheet/Core Disbond Growth in Honeycomb Sandwich Panels Subjected to Ground-Air-Ground Pressurization and In-Plane Loading

Typical damage modes in light honeycomb sandwich structures include face sheet/core disbonding and core fracture, both of which can pose a threat to the structural integrity of a component. These damage modes are of particular interest to aviation certification authorities since several in-service occurrences, such as rudder structural failure and other control surface malfunctions, have been attributed to face sheet/core disbonding. Extensive studies have shown that face sheet/core disbonding and core fracture can lead to damage propagation caused by internal pressure changes in the core. The increasing use of composite sandwich construction in aircraft applications makes it vitally important to understand the effect of ground-air-ground (GAG) cycles and conditions such as maneuver and gust loads on face sheet/core disbonding. The objective of the present study was to use a fracture mechanics based approach developed earlier to evaluate the loading at the disbond front caused by ground-air-ground pressurization and in-plane loading. A honeycomb sandwich panel containing a circular disbond at one face sheet/core interface was modeled with three-dimensional (3D) solid finite elements. The disbond was modeled as a discrete discontinuity and the strain energy release rate along the disbond front was computed using the Virtual Crack Closure Technique (VCCT). Special attention was paid to the pressure-deformation coupling which can decrease the pressure load within the disbonded sandwich section significantly when the structure is highly deformed. The commercial finite element analysis software, Abaqus/Standard, was used for the analyses. The recursive pressure-deformation coupling problem was solved by representing the entrapped air in the honeycomb cells as filled cavities in Abaqus/Standard. The results show that disbond size, face sheet thickness and core thickness are important parameters that determine crack tip loading at the disbond front. Further, the pressure-deformation coupling was found to have an important load decreasing effect [6]. In this paper, a detailed problem description is provided first. Second, the analysis methodology is presented. The fracture mechanics approach used is described and the specifics of the finite element model, including the fluid-filled cavities, are introduced. Third, the initial model verification and validation are discussed. Fourth, the findings from a closely related earlier study [6] are summarized. These findings provided the basis for the current investigation. Fifth, an aircraft ascent scenario from 0 to 12192 m (0 to 40000 ft) is considered and the resulting crack tip loading at the disbond front is determined. In-plane loading to simulate maneuvers and gust conditions are also considered. Sixth, the results are shown for a curved panel, which was used to simulate potential fuselage applications. Finally, a brief summary of observations is presented and recommendations for improvement are provided.

Chen, Zhi M.↗
NASA NTRS1999

Non-Linear Finite Element Modeling of THUNDER Piezoelectric Actuators

A NASTRAN non-linear finite element model has been developed for predicting the dome heights of THUNDER (THin Layer UNimorph Ferroelectric DrivER) piezoelectric actuators. To analytically validate the finite element model, a comparison was made with a non-linear plate solution using Von Karmen's approximation. A 500 volt input was used to examine the actuator deformation. The NASTRAN finite element model was also compared with experimental results. Four groups of specimens were fabricated and tested. Four different input voltages, which included 120, 160, 200, and 240 Vp-p with a 0 volts offset, were used for this comparison.

Taleghani, Barmac K.↗
NASA NTRS1993

Current capabilities for simulating the extreme distortion of thin structures subjected to severe impacts

The explicit transient dynamics technology in use today for simulating the impact and subsequent transient dynamic response of a structure has its origins in the 'hydrocodes' dating back to the late 1940's. The growth in capability in explicit transient dynamics technology parallels the growth in speed and size of digital computers. Computer software for simulating the explicit transient dynamic response of a structure is characterized by algorithms that use a large number of small steps. In explicit transient dynamics software there is a significant emphasis on speed and simplicity. The finite element technology used to generate the spatial discretization of a structure is based on a compromise between completeness of the representation for the physical processes modelled and speed in execution. That is, since it is expected in every calculation that the deformation will be finite and the material will be strained beyond the elastic range, the geometry and the associated gradient operators must be reconstructed, as well as complex stress-strain models evaluated at every time step. As a result, finite elements derived for explicit transient dynamics software use the simplest and barest constructions possible for computational efficiency while retaining an essential representation of the physical behavior. The best example of this technology is the four-node bending quadrilateral derived by Belytschko, Lin and Tsay. Today, the speed, memory capacity and availability of computer hardware allows a number of the previously used algorithms to be 'improved.' That is, it is possible with today's computing hardware to modify many of the standard algorithms to improve their representation of the physical process at the expense of added complexity and computational effort. The purpose is to review a number of these algorithms and identify the improvements possible. In many instances, both the older, faster version of the algorithm and the improved and somewhat slower version of the algorithm are found implemented together in software. Specifically, the following seven algorithmic items are examined: the invariant time derivatives of stress used in material models expressed in rate form; incremental objectivity and strain used in the numerical integration of the material models; the use of one-point element integration versus mean quadrature; shell elements used to represent the behavior of thin structural components; beam elements based on stress-resultant plasticity versus cross-section integration; the fidelity of elastic-plastic material models in their representation of ductile metals; and the use of Courant subcycling to reduce computational effort.

Key, Samuel W.↗
NASA NTRS2017

Curved Displacement Transfer Functions for Geometric Nonlinear Large Deformation Structure Shape Predictions

For shape predictions of structures under large geometrically nonlinear deformations, Curved Displacement Transfer Functions were formulated based on a curved displacement, traced by a material point from the undeformed position to deformed position. The embedded beam (depth-wise cross section of a structure along a surface strain-sensing line) was discretized into multiple small domains, with domain junctures matching the strain-sensing stations. Thus, the surface strain distribution could be described with a piecewise linear or a piecewise nonlinear function. The discretization approach enabled piecewise integrations of the embedded-beam curvature equations to yield the Curved Displacement Transfer Functions, expressed in terms of embedded beam geometrical parameters and surface strains. By entering the surface strain data into the Displacement Transfer Functions, deflections along each embedded beam can be calculated at multiple points for mapping the overall structural deformed shapes. Finite-element linear and nonlinear analyses of a tapered cantilever tubular beam were performed to generate linear and nonlinear surface strains and the associated deflections to be used for validation. The shape prediction accuracies were then determined by comparing the theoretical deflections with the finiteelement- generated deflections. The results show that the newly developed Curved Displacement Transfer Functions are very accurate for shape predictions of structures under large geometrically nonlinear deformations.

displacement theory↗
NASA NTRS2016

Full-Field Reconstruction of Structural Deformations and Loads from Measured Strain Data on a Wing Using the Inverse Finite Element Method

A study was undertaken to investigate the measurement of wing deformation and internal loads using measured strain data. Future aerospace vehicle research depends on the ability to accurately measure the deformation and internal loads during ground testing and in flight. The approach uses the inverse Finite Element Method (iFEM). The iFEM is a robust, computationally efficient method that is well suited for real-time measurement of real-time structural deformation and loads. The method has been validated in previous work, but has yet to be applied to a large-scale test article. This work is in preparation for an upcoming loads test of a half-span test wing in the Flight Loads Laboratory at the National Aeronautics and Space Administration Armstrong Flight Research Center (Edwards, California). The method has been implemented into an efficient MATLAB® (The MathWorks, Inc., Natick, Massachusetts) code for testing different sensor configurations. This report discusses formulation and implementation along with the preliminary results from a representative aerospace structure. The end goal is to investigate the modeling and sensor placement approach so that the best practices can be applied to future aerospace projects.

Finite element↗
NASA NTRS1989

Prediction of thermal-stress and deformations due to phase change in solidifying objects via flux/stress based finite element representations

The paper presents numerical simulations for the prediction of thermal-stress and deformation fields resulting from phase change in solidifying bodies employing new finite element representations. The formulations herein demonstrated provide different perspectives and physical interpretation for the modeling/analysis of thermo-mechanical problems and possess several inherent advantages. In comparison to traditional approaches for solving similar problems, the paper employs new flux/stress based representations to enhance the overall effectiveness. Comparative numerical applications validate applicability of the formulations for predicting the temperature induced deformations and stresses resulting from effects due to phase change.

Tamma, K. K.↗
NASA NTRS1991

Postbuckling response and failure prediction of flat rectangular graphite-epoxy plates loaded in axial compression

The objectives of the study are to assess the capability of a first-order shear deformable degenerated shell finite element theory to predict the postbuckling response and failure modes of various graphite-epoxy panels loaded in axial compression. Two panels without holes and one with a hole are studied, and the resulting responses and failure modes correlated well with the experimental results. A progressive damage failure mechanism is applied in the nonlinear analysis, which proved successful in predictions of failure location, mode, and load.

Engelstad, S. P.↗
NASA NTRS1999

Plantar soft tissue thickness during ground contact in walking

A technique is introduced for the measurement of plantar soft tissue thickness during barefoot walking. Subjects stepped into an adjustable Plexiglas frame which ensured that the required bony landmarks were appropriately positioned relative to a linear ultrasound probe connected to a conventional 7.5 MHz ultrasound scanner. Clear images of the metatarsal condyles or other foot bones were obtained throughout ground contact. Subsequent analysis of the video taped images using a motion analysis system allowed the tissue displacement to be calculated as a function of time. The tissue underneath the second metatarsal head was shown to undergo an average maximum compression of 45.7% during the late stages of ground contact during first step gait in a group of five normal subjects with a mean unloaded tissue thickness of 15.2 mm. The technique has a number of applications, including use in the validation of deformation predicted by finite element models of the soft tissue of the foot, and the study of alterations in the cushioning properties of the heel by devices which constrain the displacement of the heel pad.

Non-NASA Center↗
NASA NTRS2007

Measurement of Creep Properties of Ultra-High-Temperature Materials by a Novel Non-Contact Technique

A non-contact technique for measuring the creep properties of materials has been developed and validated as part of a collaboration among the University of Massachusetts, NASA Marshall Space Flight Center Electrostatic Levitation Facility (ESL), and the University of Tennessee. This novel method has several advantages over conventional creep testing. The sample is deformed by the centripetal acceleration from the rapid rotation, and the deformed shapes are analyzed to determine the strain. Since there is no contact with grips, there is no theoretical maximum temperature and no concern about chemical compatibility. Materials may be tested at the service temperature even for extreme environments such as rocket nozzles, or above the service temperature for accelerated testing of materials for applications such as jet engines or turbopumps for liquid-fueled engines. The creep measurements have been demonstrated to 2400 C with niobium, while the test facility, the NASA MSFC ESL, has processed materials up to 3400 C. Furthermore, the ESL creep method employs a distribution of stress to determine the stress exponent from a single test, versus the many tests required by conventional methods. Determination of the stress exponent from the ESL creep tests requires very precise measurement of the surface shape of the deformed sample for comparison to deformations predicted by finite element models for different stress exponents. An error analysis shows that the stress exponent can be determined to about 1% accuracy with the current methods and apparatus. The creep properties of single-crystal niobium at 1985 C showed excellent agreement with conventional tests performed according to ASTM Standard E-139. Tests on other metals, ceramics, and composites relevant to rocket propulsion and turbine engines are underway.

Hyers, Robert W.↗
NASA NTRSDate not supplied

Modeling Uprighting of the Orion Crew Module using a Finite Element Method with Wave Coupling and Structural Deformability

The Crew Module Uprighting System (CMUS) for the Orion Crew Module (CM) includes five bags that inflate after splashdown to rotate and upright the CM as necessary. Wave conditions can be severe enough to damage the bags and tethers, which could make the recovery operation hazardous. The CMUS therefore has design limits that ultimately constrain where the CM can land based upon sea conditions. This manuscript describes a full fluid-structure interaction computational model being used for the CMUS. It is a finite element model that uses explicit time integration to directly simulate uprighting, a computational approach usually reserved for short duration, highly transient phenomena and novelly applied here. The goal of this current effort is to determine the feasibility of this approach for CMUS uprighting, and to estimate the loads the CMUS bags experience in different wave conditions. The results from this work will be used to supplement other experimental and computational data in support of NASA’s Artemis program.

John Puryear↗
NASA NTRS2018

Curvilinear Displacement Transfer Functions for Deformed Shape Predictions of Curved Structures Using Distributed Surface Strains

Curvilinear Displacement Transfer Functions were formulated for deformed shape predictions of different curved structures using surface strains. The embedded curved beam (depth-wise cross section of a curved structure along a surface strain-sensing line) was discretized into multiple small domains, with domain junctures matching the strain-sensing stations. Thus, the surface strain distribution can be described with a piecewise linear or a piecewise nonlinear function. The discrete approach enabled piecewise integrations of a curvature-strain differential equation for the embedded curved beam to yield closed-form Curvilinear Displacement Transfer Functions, which are written in terms of embedded curved-beam geometrical parameters and surface strains. By inputting the surface strain data, the Curvilinear Displacement Transfer Functions can transform surface strains into deflections along each embedded curved beam for mapping out the overall structural deformed shapes. The finite-element method was used to analytically generate the surface strains of the curved beams. The deformed shape prediction accuracies were then determined by comparing the theoretical deflections with the finite-element-generated deflections, which were used as yardsticks. By introducing the correction factors in simple mathematical forms, the Curvilinear Displacement Transfer Functions can be quite accurate for shape predictions of different curved-beam structures ranging from limit case of straight beam up to semicircular curved beam.

Ko, William L.↗
NASA NTRS1993

Wrinkle ridges, reverse faulting, and the depth penetration of lithospheric stress in lunae planum, Mars

Tectonic features on a planetary surface are commonly used as constraints on models to determine the state of stress at the time the features formed. Quantitative global stress models applied to understand the formation of the Tharsis province on Mars constrained by observed tectonics have calculated stresses at the surface of a thin elastic shell and have neglected the role of vertical structure in influencing the predicted pattern of surface deformation. Wrinkle ridges in the Lunae Planum region of Mars form a conentric pattern of regularly spaced features in the eastern and southeastern part of Tharsis; they are formed due to compressional stresses related to the response of the Martian lithosphere to the Tharsis bulge. As observed in the exposures of valley walls in areas such as the Kasei Valles, the surface plains unit is underlain by an unconsolidated impact-generated megaregolith that grades with depth into structurally competent lithospheric basement. The ridges have alternatively been hypothesized to reflect deformation restricted to the surface plains unit ('thin skinned deformation') and deformation that includes the surface unit, megaregolith and basement lithosphere ('thick skinned deformation'). We have adopted a finite element approach to quantify the nature of deformation associated with the development of wrinkle ridges in a vertically stratified elastic lithosphere. We used the program TECTON, which contains a slippery node capability that allowed us to explicitly take into account the presence of reverse faults believed to be associated with the ridges. In this study we focused on the strain field in the vicinity of a single ridge when slip occurs along the fault. We considered two initial model geometries. In the first, the reverse fault was assumed to be in the surface plains unit, and in the second the initial fault was located in lithospheric basement, immediately beneath the weak megaregolith. We are interested in the conditions underwhich strain in the surface layer and basement either penetrates or fails to penetrate through the megaregolith. We thus address the conditions required for an initial basement fault to propagate through the megaregolith to the surface, as well as the effect of the megareolith on the strain tensor in the vicinity of a fault that nucleates in the surface plains unit.

Zuber, M. T.↗
NASA NTRS1999

Strain Rate Dependent Modeling of Polymer Matrix Composites

A research program is in progress to develop strain rate dependent deformation and failure models for the analysis of polymer matrix composites subject to high strain rate impact loads. Strain rate dependent inelastic constitutive equations have been developed to model the polymer matrix, and have been incorporated into a micromechanics approach to analyze polymer matrix composites. The Hashin failure criterion has been implemented within the micromechanics results to predict ply failure strengths. The deformation model has been implemented within LS-DYNA, a commercially available transient dynamic finite element code. The deformation response and ply failure stresses for the representative polymer matrix composite AS4/PEEK have been predicted for a variety of fiber orientations and strain rates. The predicted results compare favorably to experimentally obtained values.

Goldberg, Robert K.↗
NASA NTRS1988

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.↗
NASA NTRS1987

Advances in computational structural and solid mechanics; Proceedings of the First World Congress on Computational Mechanics, University of Texas, Austin, Sept. 22-26, 1986

Papers are presented on the nonlinear dynamic analysis of quasi-symmetric anisotropic structures, the analysis of large deformations of membrane shells by the generalized finite difference method, the vibration of shear-deformable laminated plate structures by the finite strip method, and Coon's surface method for the formulation of finite elements for plates and shells. Other topics include an evaluation of higher-order modal methods for calculating the transient structural response, a nonlinear dynamic analysis of frame structures, and a nonlinear analysis of a truss by energy minimization. Also considered are an incremental Galerkin method for plates and stiffened plates, a comparison of the postbuckling behavior of plates and shells, and a general solution of bending in a cylindrical shell.

Noor, Ahmed K.↗