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 541 records · Page 30

Stress and deformation modeling of multiple rotary combustion engine trochoid housings

This paper documents the development of the capability to produce finite element models of alternate trochoid housing configurations. The effort needed to produce these models is greatly reduced by the use of a newly developed specialized finite element preprocessor which is described. The results of static stress comparisons conducted on a Mazda trochoid housing are presented. Planned future development of this modeling capability to operational situations is also presented.

Lychuk, W. M.↗

Passive control of buckling deformation via Anderson Localization Phenomenon

Buckling problems of two types of multi-span elastic plates with transverse stiffeners are considered using a method based on the finite difference calculus. The discreteness of the stiffeners is accounted for. It is found that the torsional rigidity of the stiffener plays an important role in the buckling mode pattern. When the torsional rigidity is properly adjusted, the stiffener can act as an isolator of deformation for the structure at buckling so that the deflection is only limited to a small area.

Elishakoff, Isaac↗

The EST Model for Predicting Progressive Damage and Failure of Open Hole Bending Specimens

Progressive damage and failure in open hole composite laminate coupons subjected to flexural loading is modeled using Enhanced Schapery Theory (EST). Previous studies have demonstrated that EST can accurately predict the strength of open hole coupons under remote tensile and compressive loading states. This homogenized modeling approach uses single composite shell elements to represent the entire laminate in the thickness direction and significantly reduces computational cost. Therefore, when delaminations are not of concern or are active in the post-peak regime, the version of EST presented here is a good engineering tool for predicting deformation response. Standard coupon level tests provides all the input data needed for the model and they are interpreted in conjunction with finite element (FE) based simulations. Open hole bending test results of three different IM7/8552 carbon fiber composite layups agree well with EST predictions. The model is able to accurately capture the curvature change and deformation localization in the specimen at and during the post catastrophic load drop event.

damage↗

Constraints on crustal rheology and age of deformation from models of gravitational spreading in Ishtar Terra, Venus

Gravitational spreading is expected to lead to rapid relaxation of high relief due to the high surface temperature and associated weak crust on Venus. In this study, we use new Magellan radar and altimetry data to determine the extent of gravitational relaxation in Ishtar Terra, which contains the highest relief on Venus as well as areas of extremely high topographic slope. Within Ishtar Terra the only mountain belts found on Venus, Akna, Danu, Freyja, and Maxwell Montes, nearly encircle the smooth, high (3-4 km) plateau of Lakshmi Planum. Finite-element models of this process give expected timescales for relaxation of relief and failure at the surface. From these modeling results we attempt to constrain the strength of the crust and timescales of deformation in Ishtar Terra. Below we discuss observational evidence for gravitational spreading in Ishtar Terra, results from the finite-element modeling, independent age constraints, and implications for the rheology and timing of deformation.

Smrekar, Suzanne E.↗

Finite-Element Methods for Real-Time Simulation of Surgery

Two finite-element methods have been developed for mathematical modeling of the time-dependent behaviors of deformable objects and, more specifically, the mechanical responses of soft tissues and organs in contact with surgical tools. These methods may afford the computational efficiency needed to satisfy the requirement to obtain computational results in real time for simulating surgical procedures as described in Simulation System for Training in Laparoscopic Surgery (NPO-21192) on page 31 in this issue of NASA Tech Briefs. Simulation of the behavior of soft tissue in real time is a challenging problem because of the complexity of soft-tissue mechanics. The responses of soft tissues are characterized by nonlinearities and by spatial inhomogeneities and rate and time dependences of material properties. Finite-element methods seem promising for integrating these characteristics of tissues into computational models of organs, but they demand much central-processing-unit (CPU) time and memory, and the demand increases with the number of nodes and degrees of freedom in a given finite-element model. Hence, as finite-element models become more realistic, it becomes more difficult to compute solutions in real time. In both of the present methods, one uses approximate mathematical models trading some accuracy for computational efficiency and thereby increasing the feasibility of attaining real-time up36 NASA Tech Briefs, October 2003 date rates. The first of these methods is based on modal analysis. In this method, one reduces the number of differential equations by selecting only the most significant vibration modes of an object (typically, a suitable number of the lowest-frequency modes) for computing deformations of the object in response to applied forces.

Basdogan, Cagatay↗

Multi-scale characterization and simulation of impact welding between immiscible Mg/steel alloys

Vaporizing foil actuator spot welding method is used in this paper to join magnesium alloy AZ31 and uncoated high-strength steel DP590, which are typically considered as un-weldable due to their high physical property disparities, low mutual solubility, and the lack of any intermetallic phases. Characterization results from scanning electron microscopy (SEM) and high-resolution transmission electron microscopy (HRTEM) of the weld interface indicate that the impact creates an Mg nanocrystalline interlayer with abundant Fe particles. The interlayer exhibits intact bonding with both DP590 and AZ31 substrates. To investigate the fundamental bond formation mechanisms at the interface, a finite element (FE)-based process simulation is first performed to calculate the local temperature and deformation at the interface under the given macroscopic experimental condition. Finally, taking the FE results at the interface as inputs, molecular dynamics (MD) simulations are conducted to study the interlayer formation at the Mg/Fe interface during the impact and cooling. The results found a high velocity shearing-induced mechanical mixing mechanism that mixes Mg/Fe atoms at the interface and creates the interlayer, leading to the metallurgical bond between Mg/steel alloys.

36 MATERIALS SCIENCE↗

Analytical modeling of structure-soil systems for lunar bases

The study of the behavior of granular materials in a reduced gravity environment and under low effective stresses became a subject of great interest in the mid 1960's when NASA's Surveyor missions to the Moon began the first extraterrestrial investigation and it was found that Lunar soils exhibited properties quite unlike those on Earth. This subject gained interest during the years of the Apollo missions and more recently due to NASA's plans for future exploration and colonization of Moon and Mars. It has since been clear that a good understanding of the mechanical properties of granular materials under reduced gravity and at low effective stress levels is of paramount importance for the design and construction of surface and buried structures on these bodies. In order to achieve such an understanding it is desirable to develop a set of constitutive equations that describes the response of such materials as they are subjected to tractions and displacements. This presentation examines issues associated with conducting experiments on highly nonlinear granular materials under high and low effective stresses. The friction and dilatancy properties which affect the behavior of granular soils with low cohesion values are assessed. In order to simulate the highly nonlinear strength and stress-strain behavior of soils at low as well as high effective stresses, a versatile isotropic, pressure sensitive, third stress invariant dependent, cone-cap elasto-plastic constitutive model was proposed. The integration of the constitutive relations is performed via a fully implicit Backward Euler technique known as the Closest Point Projection Method. The model was implemented into a finite element code in order to study nonlinear boundary value problems associated with homogeneous as well as nonhomogeneous deformations at low as well as high effective stresses. The effect of gravity (self-weight) on the stress-strain-strength response of these materials is evaluated. The calibration of the model is performed via three techniques: (1) physical identification, (2) optimized calibration at the constitutive level, and (3) optimized calibration at the finite element level (Inverse Identification). Activities are summarized in graphic and outline form.

Macari-Pasqualino, Jose Emir↗

GPS Monitoring of Surface Change During and Following the Fortuitous Occurrence of the M(sub w) = 7.3 Landers Earthquake in our Network

Accomplishments: (1) Continues GPS monitoring of surface change during and following the fortuitous occurrence of the M(sub w) = 7.3 Landers earthquake in our network, in order to characterize earthquake dynamics and accelerated activity of related faults as far as 100's of kilometers along strike. (2) Integrates the geodetic constraints into consistent kinematic descriptions of the deformation field that can in turn be used to characterize the processes that drive geodynamics, including seismic cycle dynamics. In 1991, we installed and occupied a high precision GPS geodetic network to measure transform-related deformation that is partitioned from the Pacific - North America plate boundary northeastward through the Mojave Desert, via the Eastern California shear zone to the Walker Lane. The onset of the M(sub w) = 7.3 June 28, 1992, Landers, California, earthquake sequence within this network poses unique opportunities for continued monitoring of regional surface deformation related to the culmination of a major seismic cycle, characterization of the dynamic behavior of continental lithosphere during the seismic sequence, and post-seismic transient deformation. During the last year, we have reprocessed all three previous epochs for which JPL fiducial free point positioning products available and are queued for the remaining needed products, completed two field campaigns monitoring approx. 20 sites (October 1995 and September 1996), begun modeling by development of a finite element mesh based on network station locations, and developed manuscripts dealing with both the Landers-related transient deformation at the latitude of Lone Pine and the velocity field of the whole experiment. We are currently deploying a 1997 observation campaign (June 1997). We use GPS geodetic studies to characterize deformation in the Mojave Desert region and related structural domains to the north, and geophysical modeling of lithospheric behavior. The modeling is constrained by our existing and continued GPS measurements, which will provide much needed data on far-field strain accumulation across the region and on the deformational response of continental lithosphere during and following a large earthquake, forming the basis for kinematic and dynamic modeling of secular and seismic-cycle deformation. GPS geodesy affords both regional coverage and high precision that uniquely bear on these problems.

Miller, M. Meghan↗

Finite element flow-thermal-structural analysis of aerodynamically heated leading edges

Hypersonic vehicles operate in a hostile aerothermal environment which has a significant impact on their aerothermostructural performance. Significant coupling occurs between the aerodynamic flow field, structural heat transfer, and structural response creating a multidisciplinary interaction. A long term goal of the Aerothermal Loads Branch at the NASA Langley Research Center is to develop a computational capability for integrated fluid, thermal and structural analysis of aerodynamically heated structures. The integrated analysis capability includes the coupling between the fluid and the structure which occurs primarily through the thermal response of the structure, because: (1) the surface temperature affects the external flow by changing the amount of energy absorbed by the structure, and (2) the temperature gradients in the structure result in structural deformations which alter the flow field and attendant surface pressures and heating rates. In the integrated analysis, a finite element method is used to solve: (1) the Navier-Stokes equations for the flow solution, (2) the energy equation of the structure for the temperature response, and (3) the equilibrium equations of the structure for the structural deformation and stresses.

Dechaumphai, Pramote↗

Simplified and refined finite element approaches for determining stresses and internal forces in geometrically nonlinear structural analysis

Two methods for determining stresses and internal forces in geometrically nonlinear structural analysis are presented. The simplified approach uses the mid-deformed structural position to evaluate strains when rigid body rotation is present. The important feature of this approach is that it can easily be used with a general-purpose finite-element computer program. The refined approach uses element intrinsic or corotational coordinates and a geometric transformation to determine element strains from joint displacements. Results are presented which demonstrate the capabilities of these potentially useful approaches for geometrically nonlinear structural analysis.

Robinson, J. C.↗

Powder spreading, densification, and part deformation in binder jetting additive manufacturing

Binder jetting additive manufacturing (AM) can print complex structures in economical and scalable manner. Binder jetting AM comprises of deposition and weak binding of particles, known as green part, at room temperature and subsequent binder removal and sintering densification at high temperatures. However, during the densification (i.e., sintering), the part significantly deforms due to volume shrinkage. The deformation during sintering is difficult to predict, which prevents the widespread application of this technology. In this research, powder spreading simulation using discrete element method (DEM) was performed first to capture local variations in powder bed configuration. Second, finite element method (FEM) with a phenomenological continuum constitutive model was used to predict part shrinkage during the sintering process. DEM simulation showed variations in packing density, particle segregation, formation of uneven powder bed surface, and shift in particle size distribution (PSD). The sintering simulation modeled part deformation with a reasonable accuracy of 3% for solid-state sintering and intermediate liquid phase sintering. A demonstration case with non-uniform initial packing density showed that inhomogeneous green part density and PSD should be accounted for prediction of part deformation in binder jetting AM.

36 MATERIALS SCIENCE↗

A multiscale model to understand the interface chemistry, contacts, and dynamics during lithium stripping

A reversible Li-metal electrode, paired with a solid electrolyte, is critical for attaining higher energy density and safer batteries beyond the current lithium-ion cells. A stable stripping process may be even harder to attain as the stripping process will remove Li-atoms from the surface, and naturally reduce surface contact area, if not self-corrected by other mechanisms, such as diffusion and plastic deformation under an applied external stack pressure. Here, we capture these mechanisms occurring at multiple length- and time- scales, i.e., interface interactions, vacancy hopping, and plastic deformation, by integrating density functional theory (DFT) simulations, kinetic Monte Carlo (KMC), and continuum finite element method (FEM). By assuming the self-affine nature of multiscale contacts, we predict the steady-state contact area as a function of stripping current density, interface wettability, and stack pressure. We further estimate the exponential increase of overpotential due to contact area loss to maintain the same stripping current density. We demonstrate that a lithiophilic interface requires less stack pressure to reach the same steady-state contact area fraction than a lithiophobic interface. A “tolerable steady-state” contact area loss for maintaining stable stripping is estimated at 20 %, corresponding to a 10 % increase in overpotential. To constrain contact loss within the tolerance, the required stack pressure is 0.1, 0.5, and 2 times the yield strength of lithium metal for three distinct interfaces, lithiophilic Li/lithium oxide(Li2O), Li/lithium lanthanum zirconium oxide(LLZO), and lithiophoblic Li/lithium fluoride(LiF), respectively. The modeling results agree with experiments on the impact of the stack pressure quantitatively, while the discrepancy in stripping rate sensitivity is attributed to the simplifying interface interaction in our simulations. Overall, this multiscale simulation framework demonstrates the importance of electrochemical-mechanical coupling in understanding the dynamics of the Li/SE interface during stripping.

Feng, Min↗

Thin liquid film as an optical nonlinear-nonlocal medium and memory element in integrated optofluidic reservoir computer

Understanding light–matter interaction lies at the core of our ability to harness physical effects and to translate them into new capabilities realized in modern integrated photonics platforms. Here, we present the design and characterization of optofluidic components in an integrated photonics platform and computationally predict a series of physical effects that rely on thermocapillary-driven interaction between waveguide modes and topography changes of optically thin liquid dielectric film. Our results indicate that this coupling introduces substantial self-induced phase change and transmittance change in a single channel waveguide, transmittance through the Bragg grating waveguide, and nonlocal interaction between adjacent waveguides. We then employ the self-induced effects together with the inherent built-in finite relaxation time of the liquid film, to demonstrate that the light-driven deformation can serve as a reservoir computer capable of performing digital and analog tasks, where the gas–liquid interface operates both as a nonlinear actuator and as an optical memory element.

42 ENGINEERING↗

Comprehensive Material Characterization and Simultaneous Model Calibration for Improved Computational Simulation Credibility

Computational simulation is increasingly relied upon for high-consequence engineering decisions, and a foundational element to solid mechanics simulations is a credible material model. Our ultimate vision is to interlace material characterization and model calibration in a real-time feedback loop, where the current model calibration results will drive the experiment to load regimes that add the most useful information to reduce parameter uncertainty. The current work investigated one key step to this Interlaced Characterization and Calibration (ICC) paradigm, using a finite load-path tree to incorporate history/path dependency of nonlinear material models into a network of surrogate models that replace computationally-expensive finite-element analyses. Our reference simulation was an elastoplastic material point subject to biaxial deformation with a Hill anisotropic yield criterion. Training data was generated using either a space-filling or adaptive sampling method, and surrogates were built using either Gaussian process or polynomial chaos expansion methods. Surrogate error was evaluated to be on the order of 10 ⁻5 and 10 ⁻3 percent for the space-filling and adaptive sampling training data, respectively. Direct Bayesian inference was performed with the surrogate network and with the reference material point simulator, and results agreed to within 3 significant figures for the mean parameter values, with a reduction in computational cost over 5 orders of magnitude. These results bought down risk regarding the surrogate network and facilitated a successful FY22-24 full LDRD proposal to research and develop the complete ICC paradigm.

36 MATERIALS SCIENCE↗

Nonapplicability of linear finite element programs to the stress analysis of tires

A static finite element stress analysis of an inflated radial car tire was carried out. The deformed shape of the sidewall presents outward bulging. The analysis of a homogeneous isotropic toroidal shell shows that the problem is common to all solids of this type. The study suggests that the geometric stiffness due to the inflation pressure has to be taken into account. Also, the resulting large displacements make it necessary for the geometry to be updated at each load step.

Durand, M.↗

Exploiting symmetries for efficient post-buckling analysis of composite plates

Two aspects of the post-buckling analysis of composite plates are considered in this paper. The first pertains to identifying the different types of symmetry exhibited by the pre- and post-buckling responses of rectangular and skewed composite plates. A procedure is presented for exploiting these symmetries in finite-element and finite-difference analyses. The procedure can be used with existing programs having multipoint constraint capability. The second aspect pertains to the post-buckling response of biaxially loaded composite plates. Numerical results from finite-element analysis are presented which show the range of parameters for which shear deformation and/or anisotropy are important in buckling and post-buckling analysis.

Noor, A. K.↗

Surface deformation and geoid anomalies over single and double-layered convective systems

Using a primitive variable formulation of the finite-element method, the differences in the surface observables, such as topography and geoid, produced by single- and double-layered thermal convection, were compared. Both constant and depth-dependent viscosities have been considered. For the same Rayleigh number, larger surface perturbations are produced by single-cell convection. For the same Nusselt number, the magnitudes of the surface observables are greater for double-layered convection. For the same surface heat-flux, surface topographies have similar magnitudes, but the relative amplitudes of geoid anomalies depend greatly on the style of viscosity stratification. This difference in the geoid between the two systems increases with greater surface heat-flow, regardless of viscosity structure.

Koch, M.↗

Space ten-meter telescope (STMT) - Structural and thermal feasibility study of the primary mirror

The structural and thermal behavior of a ten-meter primary mirror for a space optical/near-IR telescope in geosynchronous orbit is studied. The glass-type lightweighted mirror is monolithic, of the double arch type, and is supported at only three points. The computer programs SSPTA (thermal), NASTRAN (finite element), and ACCOS V (optical) are used in sequence to determine the temperature, deformation, and optical performance of the mirror. A mirror temperature of 130 K or less appears to be obtainable by purely passive means. With a fused silica or standard Zerodur blank, thermally-induced deformation is unacceptable and cannot be fully corrected by an active secondary mirror over the desired field. Either active thermal control or a blank of lower thermal expansion coefficient would be required.

Bely, Pierre Y.↗