Search NASA⌕ Search

SEARCH · Search NASA

Results for “rate-dependent material”

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 37 records · Page 2

Elastoviscoplastic snap-through behavior of shallow arches subjected to thermomechanical loads

The problem of snap-through buckling of clamped shallow arches under thermomechanical loads is investigated. The analysis is based on nonlinear kinematic relations and nonlinear rate-dependent unified constitutive equations. A finite element approach is employed to predict the, in general, inelastic buckling behavior. The construction material is alloy B1900 + Hf, which is commonly utilized in high-temperature environments. The effect of several parameters is assessed. These parameters include the rise parameter and temperature. Comparison between elastic and elastoviscoplastic responses is also presented.

Simitses, George J.↗

Three-dimensional deformation analysis of two-phase dislocation substructures

Three-dimensional deformation analysis of two-phase dislocation substructures was carried out, extending the Qian and Fan (1991) approach to 3D stress-strain fields by using the Budiansky and Wu (1962) criterion for strain compatibility between the 'hard' and 'soft' regions. The result is a rate-dependent viscoplastic theory, named the dislocation substructure viscoplasticity (DSV), which incorporates a self-consistent effect of dislocation substructure on material response. An algorithm developed for numerical implementation of the DSV theory is presented.

Freed, Alan D.↗

Thermo-Poro-Mechanical Modeling of RTV Intumescence

Room temperature vulcanizing (RTV) silicone is a high-temperature adhesive used as a gap-filler between heatshield tiles in numerous entry missions. Its propensity to intumesce, or swell upon exposure to heat, is a well-known effect that needs to be carefully quantified during design. At tile interfaces of charring ablators, intumescence, combined with differential recession, could cause the gap filler to protrude past the ablator outer mold line, forming a “fence”. Fencing can in turn cause transition to turbulence of the flow wetting the heat shield, leading to augmented surface heating. Recent experiments conducted at the Plasmatron X facility, the high enthalpy wind tunnel of the Center for Hypersonics and Entry Systems Studies, have shown prominent fencing of RTV gap fillers in PICA, under both nitrogen and air plasmas. Similar observations are well known in the arcjet literature. Further experiments under controlled environment, performed using in situ X-ray micro-computed tomography (micro-CT) at the Advanced Light Source (ALS), have shown heating rate-dependent swelling and shrinkage of RTV during pyrolysis. To simulate RTV intumescence, a novel model was introduced in the Porous Materials Analysis Toolbox based on OpenFOAM, PATO, to account for pore-pressure buildup within both closed- and open-pores. The governing equation for the thermo-poro-mechanical response were developed, assuming linear elasticity for the charring silicone. A new multi-pyrolysis model that tracks non-monotonic advancement of material properties with pyrolysis was proposed. This model addresses the limitations of state-of-the-art ablator models to capture the different stages of thermal degradation and coupled thermomechanics. Swelling of RTV was simulated using the new thermo-poro-mechanical model and compared against in situ micro-CT data. Results showed good agreement in intumescence height and temperature profiles at all heating rates, indicating that the key factor contributing to RTV swelling is the internal pressure build-up within closed- and open-pores. As RTV is cured into a soft (rubbery) compound with low-porosity and permeability, initial temperature increase and pyrolysis gas production cause a significant increase of internal pressure, causing a pronounced volume growth. As thermal degradation progresses, rigidization of the silicone occurs due to char hardening which counteract volume shrinkage after gas pressure relief. Overall, our model shows that accounting for changes in properties such porosity, permeability and key thermomechanical coefficients is crucial for capturing the RTV volume change during ablation and enable a predictive capability for heatshield tile interface response. A plan for future calibration of thermomechanical properties evolution during degradation is discussed, as a key next step to close the new model.

RTV↗

Thermo-Poro-Mechanical Modeling of RTV Intumescence

Room temperature vulcanizing (RTV) silicone is a high-temperature adhesive used as a gap-filler between heatshield tiles in numerous entry missions. Its propensity to intumesce, or swell upon exposure to heat, is a well-known effect that needs to be carefully quantified during design. At tile interfaces of charring ablators, intumescence, combined with differential recession, could cause the gap filler to protrude past the ablator outer mold line, forming a “fence”. Fencing can in turn cause transition to turbulence of the flow wetting the heat shield, leading to augmented surface heating. Recent experiments conducted at the Plasmatron X facility, the high enthalpy wind tunnel of the Center for Hypersonics and Entry Systems Studies, have shown prominent fencing of RTV gap fillers in PICA, under both nitrogen and air plasmas. Similar observations are well known in the arcjet literature. Further experiments under controlled environment, performed using in situ X-ray micro-computed tomography (micro-CT) at the Advanced Light Source (ALS), have shown heating rate-dependent swelling and shrinkage of RTV during pyrolysis. To simulate RTV intumescence, a novel model was introduced in the Porous Materials Analysis Toolbox based on OpenFOAM, PATO, to account for pore-pressure buildup within both closed- and open-pores. The governing equation for the thermo-poro-mechanical response were developed, assuming linear elasticity for the charring silicone. A new multi-pyrolysis model that tracks non-monotonic advancement of material properties with pyrolysis was proposed. This model addresses the limitations of state-of-the-art ablator models to capture the different stages of thermal degradation and coupled thermomechanics. Swelling of RTV was simulated using the new thermo-poro-mechanical model and compared against in situ micro-CT data. Results showed good agreement in intumescence height and temperature profiles at all heating rates, indicating that the key factor contributing to RTV swelling is the internal pressure build-up within closed- and open-pores. As RTV is cured into a soft (rubbery) compound with low-porosity and permeability, initial temperature increase and pyrolysis gas production cause a significant increase of internal pressure, causing a pronounced volume growth. As thermal degradation progresses, rigidization of the silicone occurs due to char hardening which counteract volume shrinkage after gas pressure relief. Overall, our model shows that accounting for changes in properties such porosity, permeability and key thermomechanical coefficients is crucial for capturing the RTV volume change during ablation and enable a predictive capability for heatshield tile interface response. A plan for future calibration of thermomechanical properties evolution during degradation is discussed, as a key next step to close the new model.

silicone intumescence↗

Impact ejecta on the moon

The partitioning of energy and the distribution of the resultant ejecta on the moon is numerically modeled using a Eulerian finite difference grid. The impact of an iron meteoroid at 15 km/sec on a gabbroic anorthosite lunar crust is examined. The high speed impact induced flow is described over the entire hydrodynamic regime from a time where the peak pressures are 6 Mbar until the stresses everywhere in the flow are linearly elastic, and less than 5 kbar. Shock-induced polymorphic phase changes, (plagioclase and pyroxene to hollandite and perovskite), and the subsequent reversion to low pressure phases are demonstrated to enhance shock wave attenuation. A rate-dependent equation of state is used for describing the hysteretic effect of the phase change. Ballistic equations for a spherical planet are then applied to material with net velocity away from the moon.

Okeefe, J. D.↗

Impact Simulation of Composite Panels for Aerospace Applications

The aerospace research community has been working towards the development of a comprehensive framework for the design of composite structures at the component, subassembly, assembly, and system levels. Such a framework requires accurate, efficient and easy-to-use mathematical and computational models that are in sync with physical reality. The research work in this paper is designed to help an aerospace structural engineer move towards this goal by simulating dynamic impact to composite panels. First, the existing orthotropic visco-elastic visco-plastic material model in LS-DYNA commercial finite element program, MAT_COMPOSITE_TABULATED_PLASTICITY_DAMAGE (MAT_213) is enhanced with the development and use of rate-dependent damage data. Using state-of-the-art experimental procedures, impact validation test data are created. Second, the versatility of the constitutive model is examined by its ability to simulate a low velocity (27 ft/s, 8 m/s) impact test, and two high velocity impact tests (385–417 ft/s, 117–127 m/s) that straddle the ballistic limit of the projectile and the composite panels. Third, these impact events are modeled using three different finite element types – under-integrated and fully-integrated solid elements, and fully-integrated thin shell elements. Insight into accuracy versus numerical efficiency is gained by comparing the performance of solid versus thin shell elements. The numerical results show that accurate predictions of these impact events can be obtained with reasonable effort and care.

Simulations↗

Thermoviscoplastic analysis of engine cowl leading edge subjected to oscillating shock-shock interaction

A finite element thermoviscoplastic analysis method, which employs a unified constitutive model proposed by Bodner and Partom, is used to predict rate-dependent nonlinear structural behavior. The method is evaluated by predicting stress-strain behavior of a uniaxially loaded bar of nickel-based superalloy (B1900 + Hf) material. The method is used to predict the time-dependent thermoviscoplastic response of a B1900 + Hf leading edge subjected to oscillating shock-shock interaction loading. Viscoplastic analysis shows that the leading edge experiences significant plastic straining. The plastic region increases with cyclic loading in the high heat flux area.

Pandey, Ajay K.↗

"Ultra"-Fast Fracture Strength of Advanced Structural Ceramic Materials Studied at Elevated Temperatures

The accurate determination of inert strength is important in reliable life prediction of structural ceramic components. At ambient temperature, the inert strength of a brittle material is typically regarded as free of the effects of slow crack growth due to stress corrosion. Therefore, the inert strength can be determined either by eliminating active species, especially moisture, with an appropriate inert medium, or by using a very high test rate. However, at elevated temperatures, the concept or definition of the inert strength of brittle ceramic materials is not clear, since temperature itself is a degrading environment, resulting in strength degradation through slow crack growth and/or creep. Since the mechanism to control strength is rate-dependent viscous flow, the only conceivable way to determine the inert strength at elevated temperatures is to utilize a very fast test rate that either minimizes the time for or eliminates slow crack growth. Few experimental studies have measured the elevated-temperature, inert (or "ultra"-fast fracture) strength of advanced ceramics. At the NASA Lewis Research Center, an experimental study was initiated to better understand the "ultra"-fast fracture strength behavior of advanced ceramics at elevated temperatures. Fourteen advanced ceramics - one alumina, eleven silicon nitrides, and two silicon carbides - have been tested using constant stress-rate (dynamic fatigue) testing in flexure with a series of stress rates including the "ultra"-fast stress rate of 33 000 MPa/sec with digitally controlled test frames. The results for these 14 advanced ceramics indicate that, notwithstanding possible changes in flaw populations as well as flaw configurations because of elevated temperatures, the strength at 33 000 MPa/sec approached the room-temperature strength or reached a higher value than that determined at the conventional test rate of 30 MPa/sec. On the basis of the experimental data, it can be stated that the elevated-temperature, inert strength of an advanced ceramic material can be defined as the strength where no slow crack growth takes place at the temperature.

Choi, Sung R.↗

Light thermal structures and materials for high speed flight

Over the last twenty years, unified viscoplastic constitutive models have evolved to meet this need. These constitutive models provide a means for representing a material's response from the elastic through the plastic range including strain-rate dependent plastic flow, creep, and stress relaxation. Rate-dependent plasticity effects are known to be important at elevated temperatures. The purpose of this paper is to describe computational and experimental research programs underway at the Light Thermal Structures Center focused on the investigation of the response of structures and materials to local heating. In the first part of the paper, finite element thermoviscoplastic analysis is highlighted. In the second part of the paper, the thermal-structures experimental program is outlined.

Thornton, Earl A.↗

Explicit 3D continuum fracture modeling with smooth particle hydrodynamics

Impact phenomena shaped our solar system. As usual for most solar system processes, the scales are far different than we can address directly in the laboratory. Impact velocities are often much higher than we can achieve, sizes are often vastly larger, and most impacts take place in an environment where the only gravitational force is the mutual pull of the impactors. The Smooth Particle Hydrodynamics (SPH) technique has been applied in the past to the simulations of giant impacts. In these simulations, the colliding objects were so massive (at least a sizeable fraction of the Earth's mass) that material strength was negligible compared to gravity. This assumption can no longer be made when the bodies are much smaller. To this end, we have developed a 3D SPH code that includes a strength model to which we have added a von Mises yielding relation for stresses beyond the Hugoniot Elastic Limit. At the lower stresses associated with brittle failure, we use a rate-dependent strength based on the nucleation of incipient flaws whose number density is given by a Weibull distribution. Following Grady and Kipp and Melosh et al., we introduce a state variable D ('damage'), 0 less than D less than 1, which expresses the local reduction in strength due to crack growth under tensile loading. Unfortunately for the hydrodynamics, Grady and Kipp's model predicts which fragments are the most probable ones and not the ones that are really formed. This means, for example, that if a given laboratory experiment is modeled, the fragment distribution obtained from the Grady-Kipp theory would be equivalent to a ensemble average over many realizations of the experiment. On the other hand, the hydrodynamics itself is explicit and evolves not an ensemble average but very specific fragments. Hence, there is a clear incompatibility with the deterministic nature of the hydrodynamics equations and the statistical approach of the Grady-Kipp dynamical fracture model. We remedy these shortcomings by making the incipient flaw distribution explicit, i.e., particles carry activation strains which are distributed at random with a probability of occurrence given by the Weibull distribution. If the local principal axis strain exceeds this limit, damage starts to grow. By growing explicit cracks together with statistical cracks (damage) at the sub-particle scale, we ensure that material strength and fragmentation is independent of model resolution. We tested our scheme by simulating laboratory impact experiments on basalt spheres.

Benz, W.↗

A Computationally-Efficient, Multi-Mechanism Based Framework for the Comprehensive Modeling of the Evolutionary Behavior of Shape Memory Alloys

The report summarizes the accomplishments made during the 4-year duration of the project. Here, the major emphasis is placed on the different tasks performed by the two research teams; i.e., the modeling activities by the University of Akron (UA) team and the experimental and neutron diffraction studies conducted by the University of Central Florida (UCF) team, during this 4-year period. Further technical details are given in the upcoming sections by UA and UCF for each of the milestones/years (together with the corresponding figures and captions).The project majorly involved the development, validation, and application of a general theoretical model that is capable of capturing the nonlinear hysteretic responses, including pseudoelasticity, shape memory effect, rate-dependency, multi-axiality, asymmetry in tension versus compression response of shape memory alloys. Among the targeted goals for the SMA model was its ability to account for the evolutionary character response (including transient and long term behavior under sustained cycles) for both conventional and high temperature (HT) SMAs, as well as being able to simulate some of the devices which exploit these unique material systems. This required extensive (uniaxial and multi-axial) experiments needed to guide us in calibrating and characterizing the model. Moreover, since the model is formulated on the theoretical notion of internal state variables (ISVs), neutron diffraction experiments were needed to establish the linkage between the micromechanical changes and these ISVs. In addition, the design of the model should allow easy implementation in large scale finite element application to study the behavior of devices making use of these SMA materials under different loading controls. Summary of the activities, progress/achievements made during this period is given below in details for the University of Akron and the University (Section 2.0) of Central Florida (Section 3.0).

Alloys↗

Advances and trends in structures and dynamics; Proceedings of the Symposium, Washington, DC, October 22-25, 1984

Among the topics discussed are developments in structural engineering hardware and software, computation for fracture mechanics, trends in numerical analysis and parallel algorithms, mechanics of materials, advances in finite element methods, composite materials and structures, determinations of random motion and dynamic response, optimization theory, automotive tire modeling methods and contact problems, the damping and control of aircraft structures, and advanced structural applications. Specific topics covered include structural design expert systems, the evaluation of finite element system architectures, systolic arrays for finite element analyses, nonlinear finite element computations, hierarchical boundary elements, adaptive substructuring techniques in elastoplastic finite element analyses, automatic tracking of crack propagation, a theory of rate-dependent plasticity, the torsional stability of nonlinear eccentric structures, a computation method for fluid-structure interaction, the seismic analysis of three-dimensional soil-structure interaction, a stress analysis for a composite sandwich panel, toughness criterion identification for unidirectional composite laminates, the modeling of submerged cable dynamics, and damping synthesis for flexible spacecraft structures.

Noor, A. K.↗

A New Local Debonding Model with Application to the Transverse Tensile and Creep Behavior of Continuously Reinforced Titanium Composites

A new, widely applicable model for local interfacial debonding in composite materials is presented. Unlike its direct predecessors, the new model allows debonding to progress via unloading of interfacial stresses even as global loading of the composite continues. Previous debonding models employed for analysis of titanium matrix composites are surpassed by the accuracy, simplicity, and efficiency demonstrated by the new model. The new model was designed to operate seamlessly within NASA Glenn's Micromechanics Analysis Code with Generalized Method of Cells (MAC/GMC), which was employed to simulate the time- and rate-dependent (viscoplastic) transverse tensile and creep behavior of SiC/Ti composites. MAC/GMC's ability to simulate the transverse behavior of titanium matrix composites has been significantly improved by the new debonding model. Further, results indicate the need for a more accurate constitutive representation of the titanium matrix behavior in order to enable predictions of the composite transverse response, without resorting to recalibration of the debonding model parameters.

Bednarcyk, Brett A.↗

Analysis of thermomechanical fatigue of unidirectional titanium metal matrix composites

Thermomechanical fatigue (TMF) data was generated for a Ti-15V-3Cr-3Al-3Sn (Ti-15-3) material reinforced with SCS-6 silicon carbide fibers for both in-phase and out-of-phase thermomechanical cycling. Significant differences in failure mechanisms and fatigue life were noted for in-phase and out-of-phase testing. The purpose of the research is to apply a micromechanical model to the analysis of the data. The analysis predicts the stresses in the fiber and the matrix during the thermal and mechanical cycling by calculating both the thermal and mechanical stresses and their rate-dependent behavior. The rate-dependent behavior of the matrix was characterized and was used to calculate the constituent stresses in the composite. The predicted 0 degree fiber stress range was used to explain the composite failure. It was found that for a given condition, temperature, loading frequency, and time at temperature, the 0 degree fiber stress range may control the fatigue life of the unidirectional composite.

Mirdamadi, M.↗

An experimental investigation on the tensile moduli and strengths of graphite/epoxy laminates

The results of a series of tensile tests on some graphite/epoxy laminates, at rates varying from 0.002 to 2 in./min are examined. The loads were applied at various angles to the fiber directions in each case. The rate-dependent behavior of the stress-strain response is assessed. Evidence is presented to indicate that failure first occurs on inner plies, and that, in some cases, moduli increase with increasing stress (or strain) level. Lamination theory is used to predict the moduli, and comparisons with experiment are given. This theory is also used in conjunction with three failure theories to predict ultimate strengths (with varying degrees of success). Further, two approaches to ply unloading after first-ply failure are used and discussed. One is a standard method found in the literature while the other is a proposed 'strength-of-materials' type of technique which is computationally much simpler.

Yeow, Y. T.↗

Slow growth of cracks in a rate sensitive Tresca solid.

Extension of the classical theory of fracture to viscoelastic and elastic-plastic materials in which the plasticity effects are confined to a narrow band encompassing the crack front. It is suggested that the Griffith-Irwin criterion of fracture, which requires that the energy release rate computed for a given boundary value problem equal the critical threshold, ought to be replaced by a differential equation governing the slow growth of a crack prior to the onset of rapid propagation. A new term which enters the equation of motion in the dissipative media is proportional to the energy lost within the end sections of the crack, and thus it reflects the extent of inelastic behavior of a solid. A concept of 'apparent' surface energy is introduced to account for the geometry-dependent and the rate-dependent phenomena which influence toughness of an inelastic solid. Three hypotheses regarding the condition for fracture in the subcritical range of load are compared.

Wnuk, M. P.↗

Mechanical properties and morphology of poly(etheretherKetone)

Mechanical properties and morphology of poly(etheretherketone) (PEEK) were studied for samples having different thermal histories. Isothermal and rate-dependent crystallization were studied to ascertain the relationship between crystallinity/morphology and processing condition. Degree of crystallinity and microstructure were controlled by cooling the melt at different rates, ranging from quenching to slowly cooling, and by annealing amorphous material above the glass transition temperature Tg. It is found that degree of crystallinity was not as important as processing history in determining the room temperature mechanical properties. Samples with the same degree of crystallinity had very different tensile properties, depending on rate of cooling from the melt. All samples yielded by shear band formation and necked down. Quenched films had the largest breaking strains, drawing to 270 percent. Slowly cooled films exhibited ductile failure at relatively low strains. Best combined mechanical properties were obtained from semicrystalline films cooled at intermediate rates from the melt.

Cebe, Peggy↗

A crystallographic model for nickel base single crystal alloys

The purpose of this research is to develop a tool for the mechanical analysis of nickel-base single-crystal superalloys, specifically Rene N4, used in gas turbine engine components. This objective is achieved by developing a rate-dependent anisotropic constitutive model and implementing it in a nonlinear three-dimensional finite-element code. The constitutive model is developed from metallurgical concepts utilizing a crystallographic approach. An extension of Schmid's law is combined with the Bodner-Partom equations to model the inelastic tension/compression asymmetry and orientation-dependence in octahedral slip. Schmid's law is used to approximate the inelastic response of the material in cube slip. The constitutive equations model the tensile behavior, creep response and strain-rate sensitivity of the single-crystal superalloys. Methods for deriving the material constants from standard tests are also discussed. The model is implemented in a finite-element code, and the computed and experimental results are compared for several orientations and loading conditions.

Dame, L. T.↗