Two-dimensional finite-difference equations for shallow spherical shells.
Linear static deformation and buckling of shallow spherical shells under asymmetric load by solving two dimensional finite difference equations
Liepins, A. A.↗
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.
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↗
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.↗
Modular incorporation of deformable spine and 3D neck musculature into a simplified human body finite element model
Explore the source record for details and available documents.
Mitesh Lalwala↗
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.↗
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↗
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.↗
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.↗
Analysis of delamination related fracture processes in composites
This is a final report that summarizes the results achieved under this grant. The first major accomplishment is the development of the sublaminate modeling approach and shear deformation theory. The sublaminate approach allows the flexibility of considering one ply or groups of plies as a single laminated unit with effective properties. This approach is valid when the characteristic length of the response is small compared to the sublaminate thickness. The sublaminate approach was validated comparing its predictions with a finite element solution. A shear deformation theory represents an optimum compromise between accuracy and computational effort in delamination analysis of laminated composites. This conclusion was reached by applying several theories with increasing level of complexity to the prediction of interlaminar stresses and strain energy release rate in a double cracked-lap-shear configuration.
Armanios, Erian A.↗
Postseismic deformation and stress diffusion due to viscoelasticity and comments on the modified Elsasser model
Finite element techniques have been used to investigate the far-field deformation and stress changes due to asthenospheric viscoelastic relaxation following a dip-slip earthquake. The diffusion of extensional stress toward the subduction zone following a thrust earthquake on land is qualitatively consistent with the modified Elsasser model as proposed by Rydelek and Sacks (1988, 1990) to explain the coupled occurrence of land and sea earthquakes near Japan. However, the magnitude of the diffusing tensional signal is significantly smaller. The nominal model consists of a partially faulted elastic lithosphere overlying a viscoelastic substrate. Other models consider thin channel flow, rupturing of the entire elastic lithosphere, and changes in the depth of faulting. While some of these changes have significant impact on the magnitude and spatial features of the stress and deformation field, the far-field stress remains small. Numerical experiments demonstrate that the assumptions of the modified Elsasser model accentuate uniaxial deformation. When these assumptions are replicated in the finite element calculations, reasonable agreement between the models is achieved.
Cohen, Steven C.↗
Supercurrent-induced topological phase transitions
We show that finite current in superconductors can induce topological phase transitions, as a result of the deformation of the quasiparticle spectrum by a finite center-of-mass (COM) momentum of the Cooper pairs. To show the wide applicability of this mechanism, we examine the topological properties of three prototypical systems: the Kitaev chain, s-wave superconductors, and d-wave superconductors. Here we introduce a finite COM momentum as an external field corresponding to the supercurrent and show that all the models exhibit current-induced topological phase transitions. We also discuss the possibility of observing the phase transitions in experiments and the relation to the other finite COM momentum pairing states.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
On the suppression of zero energy deformation modes
Based on the Hellinger-Reissner principle and the deformation energy due to assumed stresses and displacements, the problem of the kinematic deformation modes in assumed stress hybrid/mixed finite elements has been examined. Basic schemes are developed for the choice of assumed stress terms that will suppress all kinematic deformation modes. Quadrilateral membrane and axisymmetric elements, and three-dimensional hexahedral elements, are used to illustrate the suggested procedure.
Pian, T. H. H.↗
Design, Manufacture, and Testing of an Open-Source Benchmark Composite Hydrokinetic Turbine Blade: Preprint
In a trend toward clean energy alternatives, recent years have seen great strides in the marine energy space. Consequently, there is a pressing need for the design, development, and validation of novel energy harvesting technologies such as hydrokinetic devices, which capture kinetic energy from waves, tides, and currents. However, these devices span numerous concepts and designs that often lack solid benchmark research that can be freely referenced throughout their development. This work focuses on the design process of an open-source composite hydrokinetic turbine blade for a three-bladed marine turbine rotor assembly with a diameter of 2.5 m. The proposed blade consists of two structural composite skins that are bonded with an adhesive and filled with a foam core. This study also explores and contrasts the efficiency and resolution of low-fidelity rapid design methodologies and comprehensive high-fidelity approaches in the context of blade design, modeling, and analysis efforts, a key objective in this research. Blade hydrodynamic loads were modeled and applied to finite-element blade models to study deformations and potential failure. Ongoing and upcoming efforts will result in blade manufacture and structural testing at the National Renewable Energy Laboratory. In future work, multiple blades will be deployed at the Living Bridge site at the University of New Hampshire and will be compared to rigid aluminum blades of the same geometry, developed by Sandia National Laboratories. Ultimately, this research will lay foundational groundwork for researchers and manufacturers, establishing a baseline composite blade design that will serve as a benchmark in the development of future hydrokinetic turbine blades.
blade design↗
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.↗
Manufacturing Distortions of Curved Composite Panels
This papa briefly discusses the influences of through-thickness thermal expansion, a misaligned ply, and a resin-rich slightly thicker ply on the deformations of a curved composite laminate during cool down from tbc cure temperature. Both two-dimensional and three-dimensional level finite-element analyses are used. The deformations are categorized as to radial and tangential deformations and twist, and for each of the three influences, these deformations are quantified. An additional outcome of the study is an indication of the level of analysis needed to study each of these three influences.
Starnes, James H., Jr.↗
Modeling deformation and failure in AlSi-polyester abradable sealcoating material using microstructure-based finite element simulation
A plasma-sprayed aluminum-silicon (AlSi)/polyester coating is applied in modern gas turbine engines as an abradable sealcoating to maintain tight clearances between the rotating blades and the static casing. While running the engine, the rotating blades “rub” with the abradable coating, which results in extreme strain rate (up to 10 6 ) dynamics and a high-temperature environment. Due to the difficulty of collecting direct measurements, predictive computational models are important for analyzing the deformation and failure of the abradable material, to help meet the design target of avoiding damage to the blade tip and maintaining high fuel efficiency. In this research, a microstructure-based finite element (FE) computational model was developed to capture the complex mechanical behavior of the AlSi/polyester microstructure. The model is based on a virtual representative-volume-element (RVE) of a metal-polymer microstructure, reconstructed from x-ray computed tomography. It models the plastic deformation of, and damage to, each AlSi and polyester constituents, as well as the failure at their interface. The model was calibrated and validated with uniaxial tension and compression experiments, conducted at two temperatures (298 K and 533 K) at an applied strain rate of 10 3 - 10 4 s -1 . The material exhibited strongly asymmetric tension-compression behavior and a sensitivity to temperature, which was well captured by the model. The model was further applied to investigate changes in mechanical behavior due to variations in constituents’ volume fractions, which provides guidance to the microstructural design of AlSi/polyester abradable materials. The model is expected to facilitate the development of improved abradable materials by bypassing the conventional trial-and-error approach and extensive testing requirements.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Modeling of deformation of battery cells using thick shell element formulation
Here, we describe a new approach for modeling nonlinear deformation and stress distribution of battery cells using a new thick shell finite element formulation with a through-thickness calculation of stresses and strains that satisfy equilibrium conditions. Battery cells are transversely layered materials that contain numerous thin layers in a repeating sequence. The layers are made of materials with significantly different mechanical responses. Explicit discretization of the layers is computationally impractical except for very small domains, while homogenized material models cannot account for stress and strain variations and partition through the thickness. The new formulation allows for calculation of the transverse and interlayer stresses based on the material properties of the individual cell layers. The application to problems where the transverse stresses have strong influence is also possible.
25 ENERGY STORAGE↗