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 559 records · Page 31

Geometrically Nonlinear Shell Analysis of Wrinkled Thin-Film Membranes with 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.

Tessler, Alexander↗

Large deflection behavior of quasi-isotropic laminates under low-velocity impact type point loading

Eight-ply quasi-isotropic circular composite plates of Thornel 300 graphite in Narmco 5208 epoxy resin (T300/5208) were analyzed to obtain the large deformation behavior under low-velocity impact type point loads. A simple plate-membrane coupling model was developed. The impact type point loads were replaced by equivalent quasi-static point loads. The plate-membrane coupling model was used to obtain the large deformation shapes for the thin circular composite laminates. The analyses indicated that the large deformation shapes of the composite plates under point loads vary with the center point displacements, and hence are different for different load levels. Quasi-isotropic plates were analyzed by replacing anisotropic bending stiffness components with the equivalent flexural stiffness for the isotropic plates. The plate-membrane coupling model was verified by conducting a series of tests on clamped circular quasi-isotropic laminates. Deflected shapes for the thin composite plates were experimentally obtained. These shapes agreed well with the analytically predicted shapes.

Kelkar, A.↗

Constitutive Cyclic Deformation Behavior in Single-crystal and Directionally Solidified SSME High-pressure Fuel Turbopump Airfoil Materials

The major goal of the project is to correlate mechanical properties with microstructural deformation behavior and to develop models for constitutive response under a variety of monotonic and cyclic loading cycles, temperatures, strain levels, strain rates, and environments. Two alloys are being studied as candidate SSME turbine blade materials. The first is PWA 1480, which is a single-crystal alloy whose nominal composition is reported. The second alloy being studied is D.S. Mar-M 246 + Hf. This is a directionally solidified material, and its nominal composition is also reported. The major areas of interest for the two materials will be slightly different. The single-crystal alloy lends itself well to fundamental deformation studies, since resolved shear stresses on slip planes are all known and only one grain is present. The D.S. material presents an excellent opportunity to study the effects of slightly misaligned grains on deformation behavior. The two materials will be studied by using approximately the same test matrix, so a good degree of direct comparison will also be possible.

Milligan, W. W.↗

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

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

Criteria for representing circular arc and sine wave spar webs by non-curved elements

The basic problem of how to simply represent a curved web of a spar in a finite element structural model was addressed. The ratio of flat web to curved web axial deformations and longitudinal rotations were calculated using NASTRAN models. Multiplying factors were developed from these calculations for various web thicknesses. These multiplying factors can be applied directly to the area and moment of inertia inputs of the finite element model. This allows the thermal stress relieving configurations of sine wave and circular arc webs to be simply accounted for in finite element structural models.

Jenkins, J. M.↗

Numerical Validation of Composite Panel Impact Tests

Researchers have met the challenge of modeling impact events involving composite targets for a variety of applications in a variety of ways. In this paper, the theory and implementation details of an orthotropic visco-elastic-visco-plastic material model with strain rate dependence are discussed. The model is driven by experimental data from quasi-static as well as high strain rate tests, and the data is used in defining the deformation, damage and failure sub-models. Validation data is generated by shooting a hollow 50 g Al-2024 projectile at different velocities against a flat panel target made of unidirectional composite material in a 16-ply [(0/90/45/-45)_(2)]_(S) layup. Explicit dynamic finite element analyses of four high speed tests involving one contained (projectile rebounded) and three uncontained (projectile penetrated) impact tests show that the developed material model and modeling techniques yield reasonable and acceptable predictions.

Orthotropic Composite↗

Structural Model Tuning Tool: User's Reference Manual

This report presents an efficient approach for tuning finite element models to match the measured ground vibration test and static test data. Frequencies, mode shapes, total weight, location of the center of gravity, and static deformation computed from the finite element model are matched to the measured data. The model tuning procedure used in this work is based on solving an optimization problem in which the errors for the considered metrics, between the finite element prediction and the measured data, are minimized. Analytical sensitivity values of performance indices are computed using the NASTRAN-generated sensitivity values together with the in-house computer codes, which allow for faster computational time and the use of gradient-based optimizers. The method is applied to the Aerostructures Test Wing IV model. The study shows that the military standard and the National Aeronautics and Space Administration standard for comparing analytical and experimental modal data are all satisfied. The final finite element model correlates well with the test data. The flutter speed decreases by 8.91 percent after model tuning compared with the original Aerostructures Test Wing IV design.

Chan-gi Pak↗

Designing for time-dependent material response in spacecraft structures

To study the influence on overall deformations of the time-dependent constitutive properties of fiber-reinforced polymeric matrix composite materials being considered for use in orbiting precision segmented reflectors, simple sandwich beam models are developed. The beam models include layers representing the face sheets, the core, and the adhesive bonding of the face sheets to the core. A three-layer model lumps the adhesive layers with the face sheets or core, while a five-layer model considers the adhesive layers explicitly. The deformation response of the three-layer and five-layer sandwich beam models to a midspan point load is studied. This elementary loading leads to a simple analysis, and it is easy to create this loading in the laboratory. Using the correspondence principle of viscoelasticity, the models representing the elastic behavior of the two beams are transformed into time-dependent models. Representative cases of time-dependent material behavior for the facesheet material, the core material, and the adhesive are used to evaluate the influence of these constituents being time-dependent on the deformations of the beam. As an example of the results presented, if it assumed that, as a worst case, the polymer-dominated shear properties of the core behave as a Maxwell fluid such that under constant shear stress the shear strain increases by a factor of 10 in 20 years, then it is shown that the beam deflection increases by a factor of 1.4 during that time. In addition to quantitative conclusions, several assumptions are discussed which simplify the analyses for use with more complicated material models. Finally, it is shown that the simpler three-layer model suffices in many situations.

Hyer, M. W.↗

A spectral formalism for computing three-dimensional deformations due to surface loads. 2: Present-day glacial isostatic adjustment

Using a spherically symmetric, self-gravitating, linear viscoelastic Earth model, we predict present-day three-dimensional surface deformation rates and baseline evolutions arising as a consequence of the late Pleistocene glacial cycles. In general, we use realistic models for the space-time geometry of the final late Pleistocene deglaciation event and incorporate a gravitationally self-consistent ocean meltwater redistribution. The predictions of horizontal velocity presented differ significantly, in both their amplitude and their spatial variation, from those presented in earlier analysis of others which adopted simplified models of both the late Pleistocene ice history and the Earth rheology. An important characteristic of our predicted velocity fields is that the melting of the Laurentide ice sheet over Canada is capable of contributing appreciably to the adjustment in Europe. The sensitivity of the predictions to variations in mantle rheology is investigated by considering a number of different Earth models, and by computing appropriate Frechet kernels. These calculations suggest that the sensitivity of the deformations to the Earth's rheology is significant and strongly dependent on the location of the site relative to the ancient ice sheet. The effects on the predictions of three-dimensional deformation rates of altering the ice history or adopting approximate models for the ocean meltwater redistribution have also been considered and found to be important (the former especially so). Finally, for a suite of Earth models we provide predictions of the velocity of a number of baselines in North America and Europe. We find that, in general, both radial and tangential motions contribute significantly to baseline length changes, and that these contributions are a strong function of the Earth model. We have, furthermore, found a set of Earth models which, together with the ICE-3G deglaciation chronology, produce predictions of baseline length changes that are consistent with very long baseline interferometry measurements of baselines within Europe.

Mitrovica, J. X.↗

Evolution of damage and plasticity in titanium-based, fiber-reinforced composites

The inelastic deformation mechanisms were evaluated for a model titanium-based, fiber-reinforced composite: a beta titanium alloy (Ti-15V-3Al-3Cr-3Sn) reinforced with SiC (SCS-6) fibers. The primary emphasis of this article is to illustrate the sequence in which damage and plasticity evolved for this system. The mechanical responses and the results of detailed microstructural evaluations for the 0(8), 90(8), and +/- 45(2s) line oriented laminates are provided. It is shown that the characteristics of the reaction zone around the fiber play a very important role in the way damage and plasticity evolve, particularly in the microyield regime of deformation, and must be included in any realistic constitutive model. Fiber-matrix debonding was a major damage mode for the off-axis systems. The tension test results are also compared with the predictions of a few constitutive models.

Majumdar, B. S.↗

Deformation and instability of underthrusting lithospheric plates

Models of the underthrusting lithosphere are constructed for the calculation of displacement and deflection. First, a mathematical theory is developed that rigorously demonstrates the elastic instability in the decending lithosphere. The theory states that lithospheric thrust beneath island arcs becomes unstable and suffers deflection as the compression increases. Thus, in the neighborhood of the edges where the lithospheric plate plunges into the asthenosphere and mesosphere its shape will be contorted. Next, the lateral displacement is calculated, and it is shown that, before contortion, the plate will thicken and contract at different positions with the variation in thickness following a parabolic profile. Finally, the depth distribution of the intermediate and deep focus earthquakes is explained in terms of plate buckling and contortion.

Liu, H.↗

Analysis of large, non-isothermal elastic-visco-plastic deformations

The development of a general mathematical model and solutions of test problems to analyze large nonisothermal elasto-visco-plastic deformatisms of structures is discussed. Geometric and material type nonlinearities of higher order are present in the development of the mathematical model and in the developed solution methodology.

Riff, R.↗

A Modular Approach to Model Oscillating Control Surfaces Using Navier Stokes Equations

The use of active controls for rotorcraft is becoming more important for modern aerospace configurations. Efforts to reduce the vibrations of helicopter blades with use of active-controls are in progress. Modeling oscillating control surfaces using the linear aerodynamics theory is well established. However, higher-fidelity methods are needed to account for nonlinear effects, such as those that occur in transonic flow. The aeroelastic responses of a wing with an oscillating control surface, computed using the transonic small perturbation (TSP) theory, have been shown to cause important transonic flow effects such as a reversal of control surface effectiveness that occurs as the shock wave crosses the hinge line. In order to account for flow complexities such as blade-vortex interactions of rotor blades higher-fidelity methods based on the Navier-Stokes equations are used. Reference 6 presents a procedure that uses the Navier-Stokes equations with moving-sheared grids and demonstrates up to 8 degrees of control-surface amplitude, using a single grid. Later, this procedure was extended to accommodate larger amplitudes, based on sliding grid zones. The sheared grid method implemented in EulerlNavier-Stokes-based aeroelastic code ENS AERO was successfully applied to active control design by industry. Recently there are several papers that present results for oscillating control surface using Reynolds Averaged Navier-Stokes (RANS) equations. References 9 and 10 report 2-D cases by filling gaps with overset grids. Reference 9 compares integrated forces with the experiment at low oscillating frequencies whereas Ref. 10 reports parametric studies but with no validation. Reference II reports results for a 3D case by modeling the gap region with a deformed grid and compares force results with the experiment only at the mid-span of flap. In Ref. II grid is deformed to match the control surface deflections at the section where the measurements are made. However, there is no indication in Ref. II that the gaps are explicitly modeled as in Ref. 6. Computations using overset grids are reported in Ref. 12 for a case by adding moving control surface to an existing blade but with no validation either with an experiment or another computation.

Navier↗

Sliding control of pointing and tracking with operator spline estimation

It is shown how a variable structure control technique could be implemented to achieve precise pointing and good tracking of a deformable structure subject to fast slewing maneuvers. The correction torque that has to be applied to the structure is based on estimates of upper bounds on the model errors. For a rapid rotation of the deformable structure, the elastic response can be modeled by oscillators driven by angular acceleration, and where stiffness and damping coefficients are also angular velocity and acceleration dependent. By transforming this slew-driven elastic dynamics into bilinear form (be regarding the vector made up of the angular velocity, squared angular velocity and angular acceleration components, which appear in the coefficients as the input to the deformation dynamics), an operator spline can be constructed, that gives a low order estimate of the induced disturbance. Moreover, a worst case error bound between the estimated deformation and the unknown exact deformation is also generated, which can be used where required in the sliding control correction.

Dwyer, Thomas A. W., III↗

Lithospheric structure on Venus from tectonic modelling of compressional features

In previous studies, extensional models were used that incorporated realistic rheologies in order to constrain lithospheric structure. Lithospheric modelling is considered herein from the standpoint of compressional deformation. Features of presumed compressional tectonic origin are reviewed and a model for compressional folding based on lithospheric strength envelopes are presented that include the effects of both brittle and ductile yielding as well as finite elastic strength. Model predictions are then compared with the widths and spacings of observed tectonic features and it is concluded that the results are consistent with a thin crust overlying a relatively stronger mantle, with thermal gradients probably in the range of 10 to 15 deg/km.

Banerdt, W. B.↗

A kinematic model for the late Cenozoic development of southern California crust and upper mantle

A model is developed for the young and ongoing kinematic deformation of the southern California crust and upper mantle. The kinematic model qualitatively explains both the overall seismic structure of the upper mantle and much of the known geological history of the late Cenozoic as consequences of ongoing convection beneath southern California. In this model, the high-velocity upper-mantle anomaly of the Transverse ranges is created through the convergence and sinking of the entire thickness of subcrustal lihtosphere, and the low-velocity upper-mantle anomaly beneath the Salton Trough region is attributed to high temperatures and 1-4 percent partial melt related to adiabatic decompression during mantle upwelling.

Humphreys, Eugene D.↗