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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
NASA NTRS2004

Explicit Finite Element Modeling of Multilayer Composite Fabric for Gas Turbine Engine Containment Systems: Ballistic Impact Testing - Part 2

Under the Federal Aviation Administration's Airworthiness Assurance Center of Excellence and the Aircraft Catastrophic Failure Prevention Program, National Aeronautics and Space Administration Glenn Research Center collaborated with Arizona State University, Honeywell Engines, Systems and Services, and SRI International to develop improved computational models for designing fabric-based engine containment systems. In the study described in this report, ballistic impact tests were conducted on layered dry fabric rings to provide impact response data for calibrating and verifying the improved numerical models. This report provides data on projectile velocity, impact and residual energy, and fabric deformation for a number of different test conditions.

ZYLON↗
NASA NTRS1990

Nonlinear analysis of AS4/PEEK thermoplastic composite laminate using a one parameter plasticity model

A one-parameter plasticity model was shown to adequately describe the orthotropic plastic deformation of AS4/PEEK (APC-2) unidirectional thermoplastic composite. This model was verified further for unidirectional and laminated composite panels with and without a hole. The nonlinear stress-strain relations were measured and compared with those predicted by the finite element analysis using the one-parameter elastic-plastic constitutive model. The results show that the one-parameter orthotropic plasticity model is suitable for the analysis of elastic-plastic deformation of AS4/PEEK composite laminates.

Sun, C. T.↗
NASA NTRSDate not supplied

Modeling Pore Closure in the Hot Isostatic Pressing of Additively Manufactured Inconel 718 Samples

Additive manufacturing provides opportunity for a new world of possibilities for metallic component design and fabrication; however, as-built part quality is notoriously variable, anisotropic, and degraded by porous defects. To improve quality, parts are commonly post processed using a variety of techniques including hot isostatic pressing (HIP). During HIP, high temperature and pressure is employed to improve microstructure and reduce pore volume by inducing visco-plastic deformation in the metallic substrate. Validated computational capabilities can provide predictions to guide the selection of processing condition, expected pore reduction, and the effects of varying entrapped gas in the pores. This discussion describes a micro-scale finite element structural model developed to predict pore closure due to HIP in as-built additively manufactured Inconel 718 samples produced by the laser powder bed fusion process. The model accounts for plastic deformation, creep, and internal pressure for a single pore. Pores are characterized in physical samples before and after HIP, and model predictions are compared to this measured data.

defect model prediction↗
NASA NTRS1972

Study of the dynamics of orbital assemblies including interactions with geometrical appendages. Unified flexible spacecraft load program (LOAD): Final report and user's and operation manual

The addition of a dynamic loads computation capability to the Unified Flexible Spacecraft Simulation (UFSS) program is discussed. The added capability provides a means for determining the internal member loads due to the time-varying external loading conditions experienced by an orbiting spacecraft/cluster. The flexible bodies are modeled as a system of joints or nodes which are interconnected by weightless finite element members. All masses are lumped at the joints. The orthogonal functions used to describe the spatial deformation of the bodies are normally taken to be the orthonormal cantilever modes produced by a standard structural dynamics program.

Ness, D. J.↗
NASA NTRS1985

Nonlinear transient responses of initially stressed composite plates

The nonlinear transient response of initially stressed composite plates is investigated using the finite element method. A nine-node isoparametric quadrilateral element is developed to model laminated plates under initial deformation and initial stress according to the Mindlin plate theory and von Karman large deflection assumptions. In the time integration, the Newmark constant acceleration method in conjunction with an efficient and accurate iteration scheme is used. Numerical results for deflections and bending moments for isotropic and laminated plates are obtained.

Chen, J. K.↗
NASA NTRS1993

Software to model AXAF image quality

This draft final report describes the work performed under this delivery order from May 1992 through June 1993. The purpose of this contract was to enhance and develop an integrated optical performance modeling software for complex x-ray optical systems such as AXAF. The GRAZTRACE program developed by the MSFC Optical Systems Branch for modeling VETA-I was used as the starting baseline program. The original program was a large single file program and, therefore, could not be modified very efficiently. The original source code has been reorganized, and a 'Make Utility' has been written to update the original program. The new version of the source code consists of 36 small source files to make it easier for the code developer to manage and modify the program. A user library has also been built and a 'Makelib' utility has been furnished to update the library. With the user library, the users can easily access the GRAZTRACE source files and build a custom library. A user manual for the new version of GRAZTRACE has been compiled. The plotting capability for the 3-D point spread functions and contour plots has been provided in the GRAZTRACE using the graphics package DISPLAY. The Graphics emulator over the network has been set up for programming the graphics routine. The point spread function and the contour plot routines have also been modified to display the plot centroid, and to allow the user to specify the plot range, and the viewing angle options. A Command Mode version of GRAZTRACE has also been developed. More than 60 commands have been implemented in a Code-V like format. The functions covered in this version include data manipulation, performance evaluation, and inquiry and setting of internal parameters. The user manual for these commands has been formatted as in Code-V, showing the command syntax, synopsis, and options. An interactive on-line help system for the command mode has also been accomplished to allow the user to find valid commands, command syntax, and command function. A translation program has been written to convert FEA output from structural analysis to GRAZTRACE surface deformation file (.dfm file). The program can accept standard output files and list files from COSMOS/M and NASTRAN finite analysis programs. Some interactive options are also provided, such as Cartesian or cylindrical coordinate transformation, coordinate shift and scale, and axial length change. A computerized database for technical documents relating to the AXAF project has been established. Over 5000 technical documents have been entered into the master database. A user can now rapidly retrieve the desired documents relating to the AXAF project. The summary of the work performed under this contract is shown.

Ahmad, Anees↗
NASA NTRS2000

Onset of Convection Due to Surface Tension Variations in Multicomponent and Binary Fluid Layers

Under certain conditions, such as in thin liquid films or microgravity, surface tension variations along a free surface can induce convection. Convection onset due to surface tension variation is important to many terrestrial technological processes in addition to microgravity materials processing applications. Examples include coating, drying crystallization, solidification, liquid surface contamination, and containerless processing. In double-diffusive and multicomponent systems, the spatial variations of surface tension are associated with two or more stratifying agencies, respectively. For example, both temperature and species (concentration) gradients are associated with convection in the solidification of binary alloys or salt ponds. The direction of the two (or more) gradients has a profound effect on the nature of the flow at or slightly beyond the onset of convection. Our recent work at the NASA Lewis Research Center focused on characterizing surface-tension-induced onset of convection, often referred to as Marangoni-Benard convection. Exact solutions for the stationary neutral stability of multicomponent fluid layers with interfacial deformation were derived. These solutions also permit the computation of a boundary curve that separates the long and finite wavelength instabilities. Computing points along this boundary using the exact solution (when possible) is more efficient than the typical numerical approaches, such as finite difference or spectral methods. Above the curve, a long wavelength instability was predicted, suggesting that convection would occur principally through one large flow cell in the layer, whereas below the curve, finite wavelength instabilities occur which suggest multiple finite-sized circulation cells. For many common liquids with layer depths greater than 100 mm, finite wave instability is predicted under terrestrial conditions; however, with little exception, long wavelength instability is predicted in microgravity for the identical fluid systems.

Skarda, J. Raymond Lee↗
NASA NTRS1976

Vibration of a flexible spacecraft with momentum exchange controllers

Floating reference frames were investigated in order to allow first order vibration analysis in the presence of large system rotations. When the deformations of an elastic continuum are expanded in terms of the free-free modes of an unconstrained system, the rigid body modes are found to be fixed relative to the Tisserand frame, with respect to which the relative momentum is zero. The proof presented for this is based on the orthogonality condition for modes with distinct natural frequencies. This result also guarantees the independence of coordinates for all modes with nonzero natural frequencies. A Modified Tisserand Constraint is introduced in order to define a floating reference frame with similar properties for an elastic body which contains a spinning rotor. Finite element equations of motion are derived for a completely flexible spacecraft with momentum exchange controllers, using a Modified Tisserand Frame. The deformable systems covered in this application are assumed to undergo only small rotations, and therefore the rotor torques must formally be small, although in engineering applications it may be possible to relax this constraint. A modal analysis is performed for the system and the resulting set of equations is reduced in number by a truncation procedure for more efficient system simulation.

Canavin, J. R.↗
NASA NTRS1987

A simplified computer solution for the flexibility matrix of contacting teeth for spiral bevel gears

A computer code, FLEXM, was developed to calculate the flexibility matrices of contacting teeth for spiral bevel gears using a simplified analysis based on the elementary beam theory for the deformation of gear and shaft. The simplified theory requires a computer time at least one order of magnitude less than that needed for the complete finite element method analysis reported earlier by H. Chao, and it is much easier to apply for different gear and shaft geometries. Results were obtained for a set of spiral bevel gears. The teeth deflections due to torsion, bending moment, shearing strain and axial force were found to be in the order 10(-5), 10(-6), 10(-7), and 10(-8) respectively. Thus, the torsional deformation was the most predominant factor. In the analysis of dynamic load, response frequencies were found to be larger when the mass or moment of inertia was smaller or the stiffness was larger. The change in damping coefficient had little influence on the resonance frequency, but has a marked influence on the dynamic load at the resonant frequencies.

Hsu, C. Y.↗
NASA NTRS1994

HRMA calibration handbook: EKC gravity compensated XRCF models

This document, consisting of hardcopy printout of explanatory text, figures, and tables, represents one incarnation of the AXAF high resolution mirror assembly (HRMA) Calibration Handbook. However, as we have envisioned it, the handbook also consists of electronic versions of this hardcopy printout (in the form of postscript files), the individual scripts which produced the various figures and the associated input data, the model raytrace files, and all scripts, parameter files, and input data necessary to generate the raytraces. These data are all available electronically as either ASCII or FITS files. The handbook is intended to be a living document and will be updated as new information and/or fabrication data on the HRMA are obtained, or when the need for additional results are indicated. The SAO Mission Support Team (MST) is developing a high fidelity HRMA model, consisting of analytical and numerical calculations, computer software, and databases of fundamental physical constants, laboratory measurements, configuration data, finite element models, AXAF assembly data, and so on. This model serves as the basis for the simulations presented in the handbook. The 'core' of the model is the raytrace package OSAC, which we have substantially modified and now refer to as SAOsac. One major structural modification to the software has been to utilize the UNIX binary pipe data transport mechanism for passing rays between program modules. This change has made it possible to simulate rays which are distributed randomly over the entrance aperture of the telescope. It has also resulted in a highly efficient system for tracing large numbers of rays. In one application to date (the analysis of VETA-I ring focus data) we have employed 2 x 10(exp 7) rays, a substantial improvement over the limit of 1 x 10(exp 4) rays in the original OSAC module. A second major modification is the manner in which SAOsac incorporates low spatial frequency surface errors into the geometric raytrace. The original OSAC included the ability to use Legendre-Fourier polynomials to describe deviations from the basic optical prescription. To this we have added bicubic splines to address a deficiency in the handling of the sharper deformations in the areas of mirror support pads. SAO has developed software (TRANS-FIT) to translate the most common finite element analysis models into these forms for incorporation into the raytrace program.

Tananbaum, H. D.↗
NASA NTRS1981

Plastic deformation of engines and other nonlinear structures

Plastic Analysis Capability for Engines (BOPACE3D) in nonlinear stress-analysis program based on very general family of isoparametric finite elements. Although development of BOPACE3D has been heavily influenced by requirements for engine analysis (in particular Space Shuttle main engine), it is general program applicable to many nonlinear structures.

Vos, R. G.↗
NASA NTRS1994

Patterns of brittle deformation under extension on Venus

The development of fractures at regular length scales is a widespread feature of Venusian tectonics. Models of lithospheric deformation under extension based on non-Newtonian viscous flow and brittle-plastic flow develop localized failure at preferred wavelengths that depend on lithospheric thickness and stratification. The characteristic wavelengths seen in rift zones and tessera can therefore provide constraints on crustal and thermal structure. Analytic solutions were obtained for growth rates in infinitesimal perturbations imposed on a one-dimensional, layered rheology. Brittle layers were approximated by perfectly-plastic, uniform strength, overlying ductile layers exhibiting thermally-activated power-law creep. This study investigates the formation of faults under finite amounts of extension, employing a finite-element approach. Our model incorporates non-linear viscous rheology and a Coulomb failure envelope. An initial perturbation in crustal thickness gives rise to necking instabilities. A small amount of velocity weakening serves to localize deformation into planar regions of high strain rate. Such planes are analogous to normal faults seen in terrestrial rift zones. These 'faults' evolve to low angle under finite extension. Fault spacing, orientation and location, and the depth to the brittle-ductile transition, depend in a complex way on lateral variations in crustal thickness. In general, we find that multiple wavelengths of deformation can arise from the interaction of crustal and mantle lithosphere.

Neumann, G. A.↗
NASA NTRS1998

Analysis and Design of Variable Stiffness Composite Cylinders

An investigation of the possible performance improvements of thin circular cylindrical shells through the use of the variable stiffness concept is presented. The variable stiffness concept implies that the stiffness parameters change spatially throughout the structure. This situation is achieved mainly through the use of curvilinear fibers within a fiber-reinforced composite laminate, though the possibility of thickness variations and discrete stiffening elements is also allowed. These three mechanisms are incorporated into the constitutive laws for thin shells through the use of Classical Lamination Theory. The existence of stiffness variation within the structure warrants a formulation of the static equilibrium equations from the most basic principles. The governing equations include sufficient detail to correctly model several types of nonlinearity, including the formation of a nonlinear shell boundary layer as well as the Brazier effect due to nonlinear bending of long cylinders. Stress analysis and initial buckling estimates are formulated for a general variable stiffness cylinder. Results and comparisons for several simplifications of these highly complex governing equations are presented so that the ensuing numerical solutions are considered reliable and efficient enough for in-depth optimization studies. Four distinct cases of loading and stiffness variation are chosen to investigate possible areas of improvement that the variable stiffness concept may offer over traditional constant stiffness and/or stiffened structures. The initial investigation deals with the simplest solution for cylindrical shells in which all quantities are constant around the circumference of the cylinder. This axisymmetric case includes a stiffness variation exclusively in the axial direction, and the only pertinent loading scenarios include constant loads of axial compression, pressure, and torsion. The results for these cases indicate that little improvement over traditional laminates exists through the use of curvilinear fibers, mainly due to the presence of a weak link area within the stiffness variation that limits the ultimate load that the structure can withstand. Rigorous optimization studies reveal that even though slight increases in the critical loads can be produced for designs with an arbitrary variation of the fiber orientation angle, the improvements are not significant when compared to traditional design techniques that utilize ring stiffeners and frames. The second problem that is studied involves arbitrary loading of a cylinder with a stiffness variation that changes only in the circumferential direction. The end effects of the cylinder are ignored, so that the problem takes the form of an analysis of a cross-section for a short cylinder segment. Various load cases including axial compression, pressure, torsion, bending, and transverse shear forces are investigated. It is found that the most significant improvements in load-carrying capability exist for cases which involve loads that also vary around the circumference of the shell, namely bending and shear forces. The stiffness variation of the optimal designs contribute to the increased performance in two ways: lowering the stresses in the critical areas through redistribution of the stresses; and providing a relatively stiff region that alters the buckling behavior of the structure. These results lead to an in-depth optimization study involving weight optimization of a fuselage structure subjected to typical design constraints. Comparisons of the curvilinear fiber format to traditional stiffened structures constructed of isotropic and composite materials are included. It is found that standard variable stiffness designs are quite comparable in terms of weight and load-carrying capability yet offer the added advantage of tailorability of distinct regions of the structure that experience drastically different loading conditions. The last two problems presented in this work involve the nonlinear phenomenon of long tubes under bending. Though this scenario is not as applicable to fuselage structures as the previous problems, the mechanisms that produce the nonlinear effect are ideally suited to be controlled by the variable stiffness concept. This is due to the fact that the dominating influence for long cylinders under bending is the ovalization of the cross-section, which is governed mainly by the stiffness parameters of the cylindrical shell. Possible improvement of the critical buckling moments for these structures is investigated using either a circumferential or axial stiffness variation. For the circumferential case involving infinite length cylinders, it is found that slight improvements can be observed by designing structures that resist the cross-sectional deformation yet do not detract from the buckling resistance at the critical location. The results also indicate that buckling behavior is extremely dependent on cylinder length. This effect is most easily seen in the solution of finite length cylinders under bending that contain an axial stiffness variation. For these structures, the only mechanism that exhibits improved response are those that effectively shorten the length of the cylinder, thus reducing the cross-sectional deformation due to the forced restraint at the ends. It was found that the use of curvilinear fibers was not able to achieve this effect in sufficient degree to resist the deformation, but that ring stiffeners produced the desired response admirably. Thus, it is shown that the variable stiffness concept is most effective at improving the bending response of long cylinders through the use of a circumferential stiffness variation.

Tatting, Brian F.↗
NASA NTRS1989

Stress and free vibration analyses of multilayered composite plates

This paper presents a two-phase computational procedure for the accurate prediction of vibration frequencies, stresses, and deformations in simply supported bidirectional multilayered composite plates. The range of applicability of this procedure in terms of the plate's parameters is determined, and the potential of the proposed procedure for use in conjunction with large-scale finite element modeling of composite structures is discussed. It is shown that the use of the shear correction factors calculated by this procedure extends the range of the validity of the first-order shear deformation theory, used in the first phase of the procedure, to fairly thick multilayered plates with the h/L1 ratio of the order of 0.3.

Noor, Ahmed K.↗
NASA NTRS2018

Design of a Modular Monolithic Implicit Solver for Multi-Physics Applications

The design of a modular multi-physics high-order space-time finite-element framework is presented together with its extension to allow monolithic coupling of different physics. One of the main objectives of the framework is to perform efficient high- fidelity simulations of capsule/parachute systems. This problem requires simulating multiple physics including, but not limited to, the compressible Navier-Stokes equations, the dynamics of a moving body with mesh deformations and adaptation, the linear shell equations, non-re effective boundary conditions and wall modeling. The solver is based on high-order space-time - finite element methods. Continuous, discontinuous and C1-discontinuous Galerkin methods are implemented, allowing one to discretize various physical models. Tangent and adjoint sensitivity analysis are also targeted in order to conduct gradient-based optimization, error estimation, mesh adaptation, and flow control, adding another layer of complexity to the framework. The decisions made to tackle these challenges are presented. The discussion focuses first on the "single-physics" solver and later on its extension to the monolithic coupling of different physics. The implementation of different physics modules, relevant to the capsule/parachute system, are also presented. Finally, examples of coupled computations are presented, paving the way to the simulation of the full capsule/parachute system.

Carton De Wiart, Corentin↗
NASA NTRS1988

Time-accurate unsteady aerodynamic and aeroelastic calculations for wings using Euler equations

A time-accurate approach to simultaneously solve the Euler flow equations and modal structural equations of motion is presented for computing aeroelastic responses of wings. The Euler flow eauations are solved by a time-accurate finite difference scheme with dynamic grids. The coupled aeroelastic equations of motion are solved using the linear acceleration method. The aeroelastic configuration adaptive dynamic grids are time accurately generated using the aeroelastically deformed shape of the wing. The unsteady flow calculations are validated wih experiment, both for a semi-infinite wing and a wall-mounted cantilever rectangular wings. Aeroelastic responses are computed for a rectangular wing using the modal data generated by the finite-element method. The robustness of the present approach in computing unsteady flows and aeroelastic responses that are beyond the capability of earlier approaches using the potential equations are demonstrated.

Guruswamy, Guru P.↗
NASA NTRS1984

Analysis of a transport fuselage section drop test

Transport fuselage section drop tests provided useful information about the crash behavior of metal aircraft in preparation for a full-scale Boeing 720 controlled impact demonstration (CID). The fuselage sections have also provided an operational test environment for the data acquisition system designed for the CID test, and data for analysis and correlation with the DYCAST nonlinear finite-element program. The correlation of the DYCAST section model predictions was quite good for the total fuselage crushing deflection (22 to 24 inches predicted versus 24 to 26 inches measured), floor deformation, and accelerations for the floor and fuselage. The DYCAST seat and occupant model was adequate to approximate dynamic loading to the floor, but a more sophisticated model would be required for good correlation with dummy accelerations. Although a full-section model using only finite elements for the subfloor was desirable, constraints of time and computer resources limited the finite-element subfloor model to a two-frame model. Results from the two-frame model indicate that DYCAST can provide excellent correlation with experimental crash behavior of fuselage structure with a minimum of empirical force-deflection data representing structure in the analytical model.

E. L. Fasanella↗
NASA NTRS2010

Finite Element Model for Failure Study of Two-Dimensional Triaxially Braided Composite

A new three-dimensional finite element model of two-dimensional triaxially braided composites is presented in this paper. This meso-scale modeling technique is used to examine and predict the deformation and damage observed in tests of straight sided specimens. A unit cell based approach is used to take into account the braiding architecture as well as the mechanical properties of the fiber tows, the matrix and the fiber tow-matrix interface. A 0 deg / plus or minus 60 deg. braiding configuration has been investigated by conducting static finite element analyses. Failure initiation and progressive degradation has been simulated in the fiber tows by use of the Hashin failure criteria and a damage evolution law. The fiber tow-matrix interface was modeled by using a cohesive zone approach to capture any fiber-matrix debonding. By comparing the analytical results to those obtained experimentally, the applicability of the developed model was assessed and the failure process was investigated.

Li, Xuetao↗