Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model 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 505 records · Page 28

Numerical simulations of unsteady, viscous, transonic flow over isolated and cascaded airfoils using a deforming grid

A compressible, unsteady, full Navier-Stokes, finite difference code was developed for modeling transonic flow through two-dimensional, oscillating cascades. The procedure introduces a deforming grid technique to capture the motion of the airfoils. Results using a deforming grid are presented for both isolated and cascaded airfoils. The load histories and unsteady pressure distributions are predicted for the NASA 64A010 isolated airfoil and compared with existing experimental data. Results show that the deforming grid technique can be used to successfully predict the unsteady flow properties around an oscillating airfoil. The deforming grid technique was extended for modeling unsteady flow in a cascade. The use of a deforming grid simplifies the specification of boundary conditions. Unsteady flow solutions similar to the isolated airfoil predictions are found for a NACA 0012 cascade with zero interblade phase angle and zero stagger. Experimental data for these cases are not available for code validation, but computational results are presented to show sample predictions from the code. Applications of the code to typical turbomachinery flow conditions will be presented in future work.

Huff, Dennis L.↗

Numerical simulations of unsteady, viscous, transonic flow over isolated and cascaded airfoils using a deforming grid

A compressible, unsteady, full Navier-Stokes, finite difference code was developed for modeling transonic flow through two-dimensional, oscillating cascades. The procedure introduces a deforming grid technique to capture the motion of the airfoils. Results using a deforming grid are presented for both isolated and cascaded airfoils. The load histories and unsteady pressure distributions are predicted for the NASA 64A010 isolated airfoil and compared with existing experimental data. Results show that the deforming grid technique can be used to successfully predict the unsteady flow properties around an oscillating airfoil. The deforming grid technique was extended for modeling unsteady flow in a cascade. The use of a deforming grid simplifies the specification of boundary conditions. Unsteady flow solutions similar to the isolated airfoil predictions are found for a NACA 0012 cascade with zero interblade phase angle and zero stagger. Experimental data for these cases are not available for code validation, but computational results are presented to show sample predictions from the code. Applications of the code to typical turbomachinery flow conditions will be presented in future work.

Huff, Dennis L.↗

Hierarchical tapered bar elements undergoing axial deformation

A method is described to model the dynamics of tapered axial bars of various cross sections based on the well-known Craig/Bampton component mode synthesis technique. This element is formed in terms of the static constraint modes and interface restrained normal modes. This is in contrast with the finite elements as implemented in NASTRAN where the interface restrained normal modes are neglected. These normal modes are in terms of Bessel functions. Restoration of a few of these modes leads to higher accuracy with fewer generalized coordinates. The proposed models are hierarchical so that all lower order element matrices are embedded in higher order element matrices. The advantages of this formulation compared to standard NASTRAN truss element formulation are demonstrated through simple numerical examples.

Ganesan, N.↗

Parga and Hecate Chasmata, Venus: Structure, volcanism and models of formation

Linear zones of deformation in the equatorial zone of Venus characterized by extension and volcanism have previously been identified in Pioneer Venus and Arecibo data. Two of these zones, Parga and Hecate Chasmata, are seen in Magellan data to be characterized by systems of fractures along which lie numerous coronae and corona-like features. The relationship between tectonic deformation, corona formation, and volcanism in Parga and Hecate Chasmata is examined using Magellan image and altimetry data. We examine several hypotheses for the origin of these zones, including extension and upwelling and delamination or subduction. In a companion abstract, we examine a section of Hecate Chasma in detail.

Stofan, Ellen R.↗

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

Mapping Process Model Results to Fracture Model Initial Conditions for Adhesive Bonding

A systematic approach is proposed for mapping adhesive bonding process outcomes to initial conditions for progressive damage analysis (PDA). Herein, two mapping procedures are developed: 1) direct mapping for residual stresses, strains, and deformations obtained from a process model; and 2) functional mapping for quantities such as bondline thickness, degree of cure, and porosity that are assumed to have a functional relationship with fracture toughness but are not discretely modeled in the PDA. The spatial distribution of functionally mapped quantities may be determined by process modeling or inspection. The functional maps between the quantity of interest and fracture toughness may be obtained from independent lower-scale numerical simulations or empirical test campaigns. The proposed mapping procedure links the adhesive bonding process to the structural performance. Therefore, by accounting for the effects of bondline nonuniformities in structural analysis, analytical predictive accuracy can be improved, and off-nominal conditions can be evaluated. The mapping procedures are demonstrated and verified with example problems.

Andrew C Bergan↗

Mapping Process Model Results to Fracture Model Initial Conditions for Adhesive Bonding

A systematic approach is proposed for mapping adhesive bonding process outcomes to initial conditions for progressive damage analysis (PDA). Herein, two mapping procedures are developed: 1) direct mapping for residual stresses, strains, and deformations obtained from a process model; and 2) functional mapping for quantities such as bondline thickness, degree of cure, and porosity that are assumed to have a functional relationship with fracture toughness but are not discretely modeled in the PDA. The spatial distribution of functionally mapped quantities may be determined by process modeling or inspection. The functional maps between the quantity of interest and fracture toughness may be obtained from independent lower-scale numerical simulations or empirical test campaigns. The proposed mapping procedure links the adhesive bonding process to the structural performance. Therefore, by accounting for the effects of bondline nonuniformities in structural analysis, analytical predictive accuracy can be improved, and off-nominal conditions can be evaluated. The mapping procedures are demonstrated and verified with example problems.

Andrew Bergan↗

Igneous intrusion models for floor fracturing in lunar craters

Lunar floor-fractured craters are primarily located near the maria and frequently contain ponded mare units and dark mantling deposits. Fracturing is confined to the crater interior, often producing a moat-like feature near the floor edge, and crater depth is commonly reduced by uplift of the crater floor. Although viscous relaxation of crater topography can produce such uplift, the close association of modification with surface volcanism supports a model linking floor fracture to crater-centered igneous intrusions. The consequences of two intrusion models for the lunar interior are quantitatively explored. The first model is based on terrestrial laccoliths and describes a shallow intrusion beneath the crater. The second model is based on cone sheet complexes where surface deformation results from a deeper magma chamber. Both models, their fit to observed crater modifications and possible implications for local volcanism are described.

Wichman, R. W.↗

Asteroid Fragmentation Approaches for Modeling Atmospheric Energy Deposition

During asteroid entry, energy is deposited in the atmosphere through thermal ablation and momentum-loss due to aerodynamic drag. Analytic models of asteroid entry and breakup physics are used to compute the energy deposition, which can then be compared against measured light curves and used to estimate ground damage due to airburst events. This work assesses and compares energy deposition results from four existing approaches to asteroid breakup modeling, and presents a new model that combines key elements of those approaches. The existing approaches considered include a liquid drop or "pancake" model where the object is treated as a single deforming body, and a set of discrete fragment models where the object breaks progressively into individual fragments. The new model incorporates both independent fragments and aggregate debris clouds to represent a broader range of fragmentation behaviors and reproduce more detailed light curve features. All five models are used to estimate the energy deposition rate versus altitude for the Chelyabinsk meteor impact, and results are compared with an observationally derived energy deposition curve. Comparisons show that four of the five approaches are able to match the overall observed energy deposition profile, but the features of the combined model are needed to better replicate both the primary and secondary peaks of the Chelyabinsk curve.

Modeling↗

A computational procedure for multibody systems including flexible beam dynamics

A computational procedure suitable for the solution of equations of motions for flexible multibody systems has been developed. A fully nonlinear continuum approach capable of accounting for both finite rotations and large deformations has been used to model a flexible beam component. The beam kinematics are referred directly to an inertial reference frame such that the degrees of freedom embody both the rigid and flexible deformation motions. As such, the beam inertia expression is identical to that of rigid body dynamics. The nonlinear coupling between gross body motion and elastic deformation is contained in the internal force expression. Numerical solution procedures for the integration of spatial kinematic systems can be directily applied to the generalized coordinates of both the rigid and flexible components. An accurate computation of the internal force term which is invariant to rigid motions is incorporated into the general solution procedure.

Downer, J. D.↗

Effect of propellant deformation on ignition and combustion processes in solid propellant cracks

A comprehensive theoretical model was formulated to study the development of convective burning in a solid propellant crack which continually deforms due to burning and pressure loading. In the theoretical model, the effect of interrelated structural deformation and combustion processes was taken into account by considering (1) transient, one dimensional mass, momentum, and energy conservation equations in the gas phase; (2) a transient, one dimensional heat conduction equation in the solid phase; and (3) quasi-static deformation of the two dimensional, linear viscoelastic propellant crack caused by pressure loading. Partial closures may generate substantial local pressure peaks along the crack, implying a strong coupling between chamber pressurization, crack combustion, and propellant deformation, especially when the cracks are narrow and the chamber pressurization rates high. The maximum pressure in the crack cavity is generally higher than that in the chamber. The initial flame-spreading process is not affected by propellant deformation.

Kumar, M.↗

Toward Effective Shell Modeling of Wrinkled Thin-Film Membranes Exhibiting Stress Concentrations

Geometrically nonlinear shell finite element analysis has recently been applied to solar-sail membrane problems in order to model the out-of-plane deformations due to structural wrinkling. Whereas certain problems lend themselves to achieving converged nonlinear solutions that compare favorably with experimental observations, solutions to tensioned membranes exhibiting high stress concentrations have been difficult to obtain even with the best nonlinear finite element codes and advanced shell element technology. In this paper, two numerical studies are presented that pave the way to improving the modeling of this class of nonlinear problems. The studies address the issues of mesh refinement and stress-concentration alleviation, and the effects of these modeling strategies on the ability to attain converged nonlinear deformations due to wrinkling. The numerical studies demonstrate that excessive mesh refinement in the regions of stress concentration may be disadvantageous to achieving wrinkled equilibrium states, causing the nonlinear solution to lock in the membrane response mode, while totally discarding the very low-energy bending response that is necessary to cause wrinkling deformation patterns. An element-level, strain-energy density criterion is suggested for facilitating automated, adaptive mesh refinements specifically aimed at the modeling of thin-film membranes undergoing wrinkling deformations.

Tessler, Alexander↗

Simulating Thermal Cycling and Isothermal Deformation Response of Polycrystalline NiTi

A microstructure-based FEM model that couples crystal plasticity, crystallographic descriptions of the B2-B19' martensitic phase transformation, and anisotropic elasticity is used to simulate thermal cycling and isothermal deformation in polycrystalline NiTi (49.9at% Ni). The model inputs include anisotropic elastic properties, polycrystalline texture, DSC data, and a subset of isothermal deformation and load-biased thermal cycling data. A key experimental trend is captured.namely, the transformation strain during thermal cycling is predicted to reach a peak with increasing bias stress, due to the onset of plasticity at larger bias stress. Plasticity induces internal stress that affects both thermal cycling and isothermal deformation responses. Affected thermal cycling features include hysteretic width, two-way shape memory effect, and evolution of texture with increasing bias stress. Affected isothermal deformation features include increased hardening during loading and retained martensite after unloading. These trends are not captured by microstructural models that lack plasticity, nor are they all captured in a robust manner by phenomenological approaches. Despite this advance in microstructural modeling, quantitative differences exist, such as underprediction of open loop strain during thermal cycling.

Manchiraju, Sivom↗

A crack-closure model for predicting fatigue crack growth under aircraft spectrum loading

The present paper is concerned with the development and application of an analytical model of cyclic crack growth that includes the effects of crack closure. The model was based on a concept like the Dugdale model, but was modified to leave plastically deformed material in the wake of the advancing crack tip. The model was used to correlate crack growth rates under constant-amplitude loading and to predict crack growth under aircraft spectrum loading on 2219-T851 aluminum alloy plate material. The predicted crack growth lives agreed well with experimental data. The ratio of predicted-to-experimental lives ranged from 0.66 to 1.48.

Newman, J. C., Jr.↗

Finite element implementation of state variable-based viscoplasticity models

The implementation of state variable-based viscoplasticity models is made in a general purpose finite element code for structural applications of metals deformed at elevated temperatures. Two constitutive models, Walker's and Robinson's models, are studied in conjunction with two implicit integration methods: the trapezoidal rule with Newton-Raphson iterations and an asymptotic integration algorithm. A comparison is made between the two integration methods, and the latter method appears to be computationally more appealing in terms of numerical accuracy and CPU time. However, in order to make the asymptotic algorithm robust, it is necessary to include a self adaptive scheme with subincremental step control and error checking of the Jacobian matrix at the integration points. Three examples are given to illustrate the numerical aspects of the integration methods tested.

Iskovitz, I.↗

Finite element implementation of state variable-based viscoplasticity models

The implementation of state variable-based viscoplasticity models is made in a general purpose finite element code for structural applications of metals deformed at elevated temperatures. Two constitutive models, Walker's and Robinson's models, are studied in conjunction with two implicit integration methods: the trapezoidal rule with Newton-Raphson iterations and an asymptotic integration algorithm. A comparison is made between the two integration methods, and the latter method appears to be computationally more appealing in terms of numerical accuracy and CPU time. However, in order to make the asymptotic algorithm robust, it is necessary to include a self adaptive scheme with subincremental step control and error checking of the Jacobian matrix at the integration points. Three examples are given to illustrate the numerical aspects of the integration methods tested.

Chang, T. Y. P.↗

A Study of Failure in Small Pressurized Cylindrical Shells Containing a Crack

The deformation in the vicinity of axial cracks in thin pressurized cylinders is examined using small experimental The deformation in the vicinity of axial cracks in thin pressurized cylinders is examined using small experimental models. The loading applied was either symmetric or unsymmetric about the crack plane, the latter being caused by structural constraints such as stringers. The objective was two fold - one, to provide the experimental results which will allow computer modeling techniques to be evaluated for deformations that are significantly different from that experienced by flat plates, and the other to examine the deformations and conditions associated with the onset of crack kinking which often precedes crack curving. The stresses which control crack growth in a cylindrical geometry depend on conditions introduced by the axial bulging, which is an integral part of this type of failure. For the symmetric geometry, both the hoop and radial strain just ahead off the crack, r = a, were measured and these results compared with those obtained from a variety of structural analysis codes, in particular STAGS [1], ABAQUS and ANSYS. In addition to these measurements, the pressures at the onset of stable and unstable crack growth were obtained and the corresponding crack deformations measured as the pressures were increased to failure. For the unsymmetric cases, measurements were taken of the crack kinking angle, and the displacements in the vicinity of the crack. In general, the strains ahead of the crack showed good agreement between the three computer codes and between the codes and the experiments. In the case of crack behavior, it was determined that modeling stable tearing with a crack-tip opening displacement fracture criterion could be successfully combined with the finite-element analysis techniques as used in structural analysis codes. The analytic results obtained in this study were very compatible with the experimental observations of crack growth. Measured crack kinking angles also showed good agreement with theories based on the maximum principle stress criterion.

Barwell, Craig A.↗

The shallow structure of the Martian lithosphere in the vicinity of the ridged plains

On Mars, the most prominent assemblage of regularly spaced plain ridges occurs in the plains of the Coprates and Lunae Planum regions. In this paper, a number of tectonic models is developed and quantitatively evaluated, that relate the regular spacing of these ridges to the shallow internal structure of Mars. It is shown that models with both rigid and deformable megaregolith-basement interface conditions can provide solutions which can explain the ridge spacing within the constraint provided by the estimated thickness of the smooth plains materials. Implications of the models for the structure and state of stress in the Tharsis region at the time of ridge formation are discussed.

Zuber, M. T.↗