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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

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

Dynamics of spacecraft with deploying flexible appendages

A computational formulation for the dynamic analysis of spacecraft with deploying appendages is presented. The appendage model is based on a geometrically nonlinear beam formulation which accurately accounts for large rotational and large deformation motions. A moving finite element reference grid is incorporated within the nonlinear beam formulation to model the deployment motion. Hamilton's Law is used to formulate the general equations of motion, and a transient integration solution procedure is derived from a space-time finite element discretization of the Hamiltonian variational statement. Computational results of the methodology are presented for a classical gravity gradient stabilized satellite configuration.

Downer, Janice D.↗
NASA NTRS1992

The gabbro-eclogite phase transition and the elevation of mountain belts on Venus

Among the four mountain belts surrounding Lakshmi Planum, Maxwell Montes is the highest and stands up to 11 km above the mean planetary radius and 7 km above Lakshmi Planum. The bulk composition and radioactive heat production of the crust on Venus, where measured, are similar to those of terrestrial tholeiitic basalt. Because the thickness of the low-density crust may be limited by the gabbro-garnet granulite-eclogite phase transitions, the 7-11 km maximum elevation of Maxwell Montes is difficult to understand except in the unlikely situation that the crust contains a large volume of magma. A possible explanation is that the base of the crust is not in phase equilibrium. It has been suggested that under completely dry conditions, the gabbro-eclogite phase transition takes place by solid-state diffusion and may require a geologically significant time to run to completion. Solid-state diffusion is a strongly temperature-dependent process. In this paper we solve the thermal evolution of the mountain belt to attempt to constrain the depth of the gabbro-eclogite transition and thus to assess this hypothesis quantitatively. The one-dimensional heat equation is solved numerically by a finite difference approximation. The deformation of the horizontally shortening crustal and mantle portions of the thermal boundary layer is assumed to occur by pure shear, and therefore the vertical velocity is given by the product of the horizontal strain rate and depth.

Namiki, Noriyuki↗
NASA NTRS1990

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. The flexible beams are modeled using a fully nonlinear theory which accounts for both finite rotations and large deformations. The present formulation incorporates physical measures of conjugate Cauchy stress and covariant strain increments. As a consequence, the beam model can easily be interfaced with real-time strain measurements and feedback control systems. A distinct feature of the present work is the computational preservation of total energy for undamped systems; this is obtained via an objective strain increment/stress update procedure combined with an energy-conserving time integration algorithm which contains an accurate update of angular orientations. The procedure is demonstrated via several example problems.

Downer, J. D.↗
NASA NTRS2002

3D Finite Element Analysis of Particle-Reinforced Aluminum

Deformation in particle-reinforced aluminum has been simulated using three distinct types of finite element model: a three-dimensional repeating unit cell, a three-dimensional multi-particle model, and two-dimensional multi-particle models. The repeating unit cell model represents a fictitious periodic cubic array of particles. The 3D multi-particle (3D-MP) model represents randomly placed and oriented particles. The 2D generalized plane strain multi-particle models were obtained from planar sections through the 3D-MP model. These models were used to study the tensile macroscopic stress-strain response and the associated stress and strain distributions in an elastoplastic matrix. The results indicate that the 2D model having a particle area fraction equal to the particle representative volume fraction of the 3D models predicted the same macroscopic stress-strain response as the 3D models. However, there are fluctuations in the particle area fraction in a representative volume element. As expected, predictions from 2D models having different particle area fractions do not agree with predictions from 3D models. More importantly, it was found that the microscopic stress and strain distributions from the 2D models do not agree with those from the 3D-MP model. Specifically, the plastic strain distribution predicted by the 2D model is banded along lines inclined at 45 deg from the loading axis while the 3D model prediction is not. Additionally, the triaxial stress and maximum principal stress distributions predicted by 2D and 3D models do not agree. Thus, it appears necessary to use a multi-particle 3D model to accurately predict material responses that depend on local effects, such as strain-to-failure, fracture toughness, and fatigue life.

Shen, H.↗
NASA NTRS2009

Advances in Non-contact Measurement of Creep Properties

Our team has developed a novel approach to measuring creep at extremely high temperatures using electrostatic levitation (ESL). This method has been demonstrated on niobium up to 2300 C, while ESL has melted tungsten (3400 C). High-precision machined spheres of the sample are levitated in the NASA MSFC ESL, a national user facility, and heated with a laser. The laser is aligned off-center so that the absorbed photons transfer their momentum to the sample, causing it to rotate at up to 250,000+ RPM. The rapid rotation loads the sample through centripetal acceleration, causing it to deform. The deformation of the sample is captured on high-speed video, which is analyzed by machine-vision software from the University of Massachusetts. The deformations are compared to finite element models to determine the constitutive constants in the creep relation. Furthermore, the noncontact method exploits stress gradients within the sample to determine the stress exponent in a single test. This method was validated in collaboration with the University of Tennessee for niobium at 1985 C, with agreement within the uncertainty of the conventional measurements. A similar method is being employed on Ultra-High-Temperature ZrB2- SiC composites, which may see application in rocket nozzles and sharp leading edges for hypersonic vehicles.

Hyers, Robert W.↗
NASA NTRS2010

Improved Creep Measurements for Ultra-High Temperature Materials

Our team has developed a novel approach to measuring creep at extremely high temperatures using electrostatic levitation (ESL). This method has been demonstrated on niobium up to 2300 C, while ESL has melted tungsten (3400 C). This method has been extended to lower temperatures and higher stresses and applied to new materials, including a niobium-based superalloy, MASC. High-precision machined spheres of the sample are levitated in the NASA MSFC ESL, a national user facility and heated with a laser. The samples are rotated with an induction motor at up to 30,000 revolutions per second. The rapid rotation loads the sample through centripetal acceleration, producing a shear stress of about 60 MPa at the center, causing the sample to deform. The deformation of the sample is captured on high-speed video, which is analyzed by machine-vision software from the University of Massachusetts. The deformations are compared to finite element models to determine the constitutive constants in the creep relation. Furthermore, the non-contact method exploits stress gradients within the sample to determine the stress exponent in a single test.

Hyers, Robert W.↗
NASA NTRS1982

Dynamic responses of graphite/epoxy laminated beam to impact of elastic spheres

Wave propagation in 90/45/90/-45/902s and 0/45/0/-45/02s laminates of a graphite/epoxy composite due to impact of a steel ball was investigated experimentally and also by using a high order beam finite element. Dynamic strain responses at several locations were obtained using strain gages. The finite element program which incorporated statically determined contact laws was employed to calculate the contact force history as well as the target beam dynamic deformation. The comparison of the finite element solutions with the experimental data indicated that the static contact laws for loading and unloading (developed under this grant) are adequate for the dynamic impact analysis. It was found that for the 0/45/0/-45/02s laminate which has a much larger longitudinal bending rigidity, the use of beam finite elements is not suitable and plate finite element should be used instead.

Sun, C. T.↗
NASA NTRS1997

Application of Interface Technology in Nonlinear Analysis of a Stitched/RFI Composite Wing Stub Box

A recently developed interface technology was successfully employed in the geometrically nonlinear analysis of a full-scale stitched/RFI composite wing box loaded in bending. The technology allows mismatched finite element models to be joined in a variationally consistent manner and reduces the modeling complexity by eliminating transition meshing. In the analysis, local finite element models of nonlinearly deformed wide bays of the wing box are refined without the need for transition meshing to the surrounding coarse mesh. The COMET-AR finite element code, which has the interface technology capability, was used to perform the analyses. The COMET-AR analysis is compared to both a NASTRAN analysis and to experimental data. The interface technology solution is shown to be in good agreement with both. The viability of interface technology for coupled global/local analysis of large scale aircraft structures is demonstrated.

Wang, John T.↗
NASA NTRS1977

Hamilton's principle - Finite-element methods and flexible body dynamics

A variational formulation is given for the equations of motion for an unconstrained elastic body, and Hamilton's principle is used to derive the equations of motion and deformation of the body. Finite-element approximations are developed for these dynamical equations with respect to a body axis system satisfying the mean axis conditions. The free-body influence matrix for the body then is developed in terms of the finite-element model parameters.

Cavin, R. K., III↗
NASA NTRS1990

A computational procedure for the dynamics of flexible beams within multibody systems

The dynamic analysis of three dimensional elastic beams which experience large rotational and large deformational motions are examined. The beam motion is modeled using an inertial reference for the translational displacements and a body-fixed reference for the rotational quantities. Finite strain rod theories are then defined in conjunction with the beam kinematic description which accounts for the effects of stretching, bending, torsion, and transverse shear deformations. A convected coordinate representation of the Cauchy stress tensor and a conjugate strain definition is introduced to model the beam deformation. To treat the beam dynamics, a two-stage modification of the central difference algorithm is presented to integrate the translational coordinates and the angular velocity vector. The angular orientation is then obtained from the application of an implicit integration algorithm to the Euler parameter/angular velocity kinematical relation. The combined developments of the objective internal force computation with the dynamic solution procedures result in the computational preservation of total energy for undamped systems. The present methodology is also extended to model the dynamics of deployment/retrieval of the flexible members. A moving spatial grid corresponding to the configuration of a deployed rigid beam is employed as a reference for the dynamic variables. A transient integration scheme which accurately accounts for the deforming spatial grid is derived from a space-time finite element discretization of a Hamiltonian variational statement. The computational results of this general deforming finite element beam formulation are compared to reported results for a planar inverse-spaghetti problem.

Downer, Janice Diane↗
NASA NTRS1986

Correlation of data on strain accumulation adjacent to the San Andreas Fault with available models

Theoretical and numerical studies of deformation on strike slip faults were performed and the results applied to geodetic observations performed in the vicinity of the San Andreas Fault in California. The initial efforts were devoted to an extensive series of finite element calculations of the deformation associated with cyclic displacements on a strike-slip fault. Measurements of strain accumulation adjacent to the San Andreas Fault indicate that the zone of strain accumulation extends only a few tens of kilometers away from the fault. There is a concern about the tendency to make geodetic observations along the line to the source. This technique has serious problems for strike slip faults since the vector velocity is also along the fault. Use of a series of stations lying perpendicular to the fault whose positions are measured relative to a reference station are suggested to correct the problem. The complexity of faulting adjacent to the San Andreas Fault indicated that the homogeneous elastic and viscoelastic approach to deformation had serious limitations. These limitation led to the proposal of an approach that assumes a fault is composed of a distribution of asperities and barriers on all scales. Thus, an earthquake on a fault is treated as a failure of a fractal tree. Work continued on the development of a fractal based model for deformation in the western United States. In order to better understand the distribution of seismicity on the San Andreas Fault system a fractal analog was developed. The fractal concept also provides a means of testing whether clustering in time or space is a scale-invariant process.

Turcotte, Donald L.↗
NASA NTRS1997

Static Aeroelastic Analysis of Transonic Wind Tunnel Models Using Finite Element Methods

A computational method for accurately predicting the static aeroelastic deformations of typical transonic transport wind tunnel models is described. The method utilizes a finite element method (FEM) for predicting the deformations. Extensive calibration/validation of this method was carried out using a novel wind-off wind tunnel model static loading experiment and wind-on optical wing twist measurements obtained during a recent wind tunnel test in the National Transonic Facility (NTF) at NASA LaRC. Further validations were carried out using a Navier-Stokes computational fluid dynamics (CFD) flow solver to calculate wing pressure distributions about several aeroelastically deformed wings and comparing these predictions with NTF experimental data. Results from this aeroelastic deformation method are in good overall agreement with experimentally measured values. Including the predicted deformations significantly improves the correlation between CFD predicted and experimentally measured wing & pressures.

Hooker, John R.↗
NASA NTRSDate not supplied

Multiscale Mesh Adaptation for Transonic Aeroelastic Flutter Problems

This work applies multiscale mesh adaptation with refine to reduce spatial discretization error of aeroelastic computational fluid dynamics (CFD) simulations. Benchmark flutter models, such as the pitch and plunge NACA64A-010 airfoil and the benchmark supercritical wing, are studied with both a linearized frequency-domain solver and time-marching CFD coupled to a modal structural solver in FUN3D. The undeformed NASA Common Research Model (CRM), an aeroelastic jig shape variant of the CRM, is also studied with the linearized frequency-domain approach. For these cases, the adaptation process converges to comparable flutter predictions to hand-generated meshes but with smaller node counts. However the additional disciplines of the linearized frequency-domain analysis, the mesh deformation, and the unsteady finite-volume solver create robustness challenges that need to be addressed before it can be applied as a fully automated process for complex transonic aeroelastic problems. In particular, negative volumes are observed to be an issue for FUN3D’s linear elasticity mesh deformation solver when moving the adapted meshes.

Aeroelasticity↗
NASA NTRS1991

The width of fault zones in a brittle-viscous lithosphere: Strike-slip faults

A fault zone in an ideal brittle material overlying a very weak substrate could, in principle, consist of a single slip surface. Real fault zones have a finite width consisting of a number of nearly parallel slip surfaces on which deformation is distributed. The hypothesis that the finite width of fault zones reflects stresses due to quasistatic flow in the ductile substrate of a brittle surface layer is explored. Because of the simplicity of theory and observations, strike-slip faults are examined first, but the analysis can be extended to normal and thrust faulting.

Parmentier, E. M.↗
NASA NTRS2002

Warpage of Large Curved Composite Panels due to Manufacturing Anomalies

This paper discusses the influences of a misaligned layer, a resin-rich slightly thicker layer, and a small thermal gradient on the thermally-induced deformations of large curved composite panels during cooldown from their cure temperature. The deformations represent warpage of the panels due to anomalies that occur during layup, consolidation, and cure. Two-dimensional finite element analyses are used The deformations are categorized as to their impact on circumferential and twist warpage metrics. The results are intended to highlight the sensitivity of manufactured panel shape to the various unwanted effects that can occur during manufacturing.

Starnes, James H., Jr.↗