Search NASA⌕ Search

SEARCH · Search NASA

Results for “constitutive relations”

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

Effective constitutive relations for the microstructure of periodic frames

Effective micro-structure properties of two-dimensional large repetitive frame-like structures are derived using straightforward combinations of energy methods of linear transformations. Once the actual structure is identified symmetry considerations are used in order to identify its independent property constants. The actual values of these constants are constructed according to a building block format which is carried out in the three consecutive steps: (1) all basic planar lattices are identified, (2) effective continuum properties are derived for each of these planar basic grids using energy methods, and (3) orthogonal transformations are finally used to determine the contribution of each basic set to the overall effective continuum properties of the structure.

Nayfeh, A. H.↗

Modeling of nonlinear constitutive relations of woven ceramic composites

The linear and nonlinear tensile behavior of SiC/SiC composites is examined. This analysis is based on the observed stress-strain behavior and damage evolution from a series of loading and unloading tensile tests. Stress distribution of the composite with transverse cracks in the transverse yarns is examined, and the nonlinearity of the composite stress-strain relation due to the presence of transverse cracks is predicted by application of a maximum failure strain criterion for the transverse yarn.

Kuo, Wen-Shyong↗

Modeling of Wall-Bounded Complex Flows and Free Shear Flows

Various wall-bounded flows with complex geometries and free shear flows have been studied with a newly developed realizable Reynolds stress algebraic equation model. The model development is based on the invariant theory in continuum mechanics. This theory enables us to formulate a general constitutive relation for the Reynolds stresses. Pope was the first to introduce this kind of constitutive relation to turbulence modeling. In our study, realizability is imposed on the truncated constitutive relation to determine the coefficients so that, unlike the standard k-E eddy viscosity model, the present model will not produce negative normal stresses in any situations of rapid distortion. The calculations based on the present model have shown an encouraging success in modeling complex turbulent flows.

Shih, Tsan-Hsing↗

Variational Asymptotic Homogenization of Finitely Deformed Viscoelastic-Viscoplastic Composites

The objective of this paper is to develop a constitutive model for finitely deformed viscoelastic-viscoplastic materials and a micromechanics approach to homogenizing composites consisting of such materials. The development of the constitutive model involves establishing a thermodynamic framework based on finite strain theory, developing a viscoelasticity and a viscoplasticity model based on the thermodynamic framework, developing a radial return algorithm based on a classic framework, and deriving a closed-from incremental constitutive relation in the spatial configuration. The development of the micromechanics approach involves pulling-back the above constitutive relation to the material configuration, formulating a variational statement with the resulting constitutive relation, discretizing variational statement in a finite-dimensional space, and solving the discretized variational statement using an Euler–Newton predictor–corrector method. The constitutive model is calibrated via monotonic uniaxial compression tests on a polymer, and the calibrated model is validated by comparing its predictions with the cyclic test data. It is shown capable of characterizing viscoelasticity, viscoplasticity, and complex loading paths. The micromechanics approach’s capabilities are demonstrated through homogenizing a unidirectional fiber-reinforced composite, subjected to uniaxial, biaxial, and shear loading, at different strain rates. It is demonstrated to be capable of handling rate dependence and complex loading paths. The present framework can be further improved by implementing more sophisticated viscoelasticity and viscoplasticity models in future work.

Finite element analysis↗

Nucleation and growth of rolling contact failure of 440C bearing steel

A 'two-body' elasto-plastic finite element model of 2-dimensional rolling and rolling-plus-sliding was developed to treat the effect of surface irregularities. The model consists of a smooth cylinder in contact with a semi-infinite half-space that is either smooth or fitted with one of 0.4 microns deep or 7 microns deep groove, or a 0.4 microns high ridge-like asperity. The model incorporates elastic-linear-kinematic hardening-plastic (ELKP) and non-linear-kinematic hardening-plastic (NLKP) material constitutive relations appropriate for hardened bearing steel and the 440C grade. The calculated contact pressure distribution is Hertzian for smooth body contact, and it displays intense, stationary, pressure spikes superposed on the Hertzian pressure for contact with the grooved and ridged surface. The results obtained for the 0.4 microns deep groove compare well with those reported by Elsharkawy and Hamrock for an EHD lubricated contact. The effect of translating the counterface on the half space as opposed to indenting the half space with the counter face with no translation is studied. The stress and strain values near the surface are found to be similar for the two cases, whereas they are significantly different in the subsurface. It is seen that when tiny shoulders are introduced at the edge of the groove in the finite element model, the incremental plasticity and residual stresses are significantly higher in the vicinity of the right shoulder (rolling direction is from left to right) than at the left shoulder. This may explain the experimental observation that the spall nucleation occurs at the exit end of the artificially planted indents. Pure rolling calculations are compared with rolling + sliding calculations. For a coefficient of friction, mu = 0.1, the effect of friction is found to be small. Efforts were made to identify the material constitutive relations which best describe the deformation characteristics of the bearing steels in the initial few cycles. Elastic-linear-kinematic hardening-plastic (ELKP) material constitutive relations produce less net plastic deformation in the initial stages for a given stress, than seen in experiments. A new set of constitutive relations: non-linear-kinematic hardening-plastic (NLKP) was used. This material model produces more plasticity than the ELKP model and shows promise for treating the net distortions in the early stages. Techniques for performing experimental measurements that can be compared with the finite element calculations were devised. The measurements are being performed on 9mm-diameter, 440C steel cylindrical rolling elements in contact with 12.5 mm-diameter, 52100 steel balls in a 3-ball-rod fatigue test machine operating at 3600 RPM. Artificial, 7 microns deep, indents were inserted on the running track of the cylindrical rolling elements and profilometer measurements of these indents made, before and after the rolling. These preliminary measurements show that the indents are substantially deformed plastically in the process of rolling. The deformations of the groove calculated with the finite element model are comparable to those measured experimentally.

Gupta, V.↗

Interpretation and modeling of the averaged equations for a fluid-solid flow

A self-consistent derivation of the conservation laws is given for flows of a fluid-solid mixture. A unified analytical framework for obtaining constitutive relations is provided. This analysis uses a control volume/control surface approach that is widely used in fluid mechanics. All terms in the governing equations and the constitutive relations are written in terms of the mass-weighted averages except solid concentration. It is believed that the mass-weighted average is the natural bridge between micromechanics and constitutive relations. The derived momentum equations contain terms that differ from all existing models except that of Prosperetti and Jones (1984). However, their assumptions are not needed here. Special attention is given to the solid phase pressure. The physical basis of the previously assumed form for this pressure (Givler 1987) becomes clear. A number of related phenomena are also discussed. These include the anti-diffusion and anisotropic normal stresses. The energy equations are also different from existing models.

Shen, Hayley H.↗

Finite element analysis of elasto-plastic soils. Report no. 4: Finite element analysis of elasto-plastic frictional materials for application to lunar earth sciences

The behavior of finite element models employing different constitutive relations to describe the stress-strain behavior of soils is investigated. Three models, which assume small strain theory is applicable, include a nondilatant, a dilatant and a strain hardening constitutive relation. Two models are formulated using large strain theory and include a hyperbolic and a Tresca elastic perfectly plastic constitutive relation. These finite element models are used to analyze retaining walls and footings. Methods of improving the finite element solutions are investigated. For nonlinear problems better solutions can be obtained by using smaller load increment sizes and more iterations per load increment than by increasing the number of elements. Suitable methods of treating tension stresses and stresses which exceed the yield criteria are discussed.

Marr, W. A., Jr.↗

Strength, Fracture Toughness, Fatigue, and Standardization Issues of Free-standing Thermal Barrier Coatings

Strength, fracture toughness and fatigue behavior of free-standing thick thermal barrier coatings of plasma-sprayed ZrO2-8wt % Y2O3 were determined at ambient and elevated temperatures in an attempt to establish a database for design. Strength, in conjunction with deformation (stress-strain behavior), was evaluated in tension (uniaxial and trans-thickness), compression, and uniaxial and biaxial flexure; fracture toughness was determined in various load conditions including mode I, mode II, and mixed modes I and II; fatigue or slow crack growth behavior was estimated in cyclic tension and dynamic flexure loading. Effect of sintering was quantified through approaches using strength, fracture toughness, and modulus (constitutive relations) measurements. Standardization issues on test methodology also was presented with a special regard to material's unique constitutive relations.

Choi, Sung R.↗

Deformation and Tensile Cyclic Fatigue of Plasma-Sprayed ZrO2-8wt% Y2O3 Thermal Barrier Coatings

Deformation (constitutive relations) of free-standing, thick thermal barrier coatings of sprayed ZrO2-8Wt% Y2O3 was determined at ambient temperature in both pure tension and pure compression using cylindrical bar test specimens. The material exhibited both significant nonlinearity and hysteresis in its load-strain curves, The load-strain relations in four-point uniaxial flexure were determined from tension and compression sides and were compared with individual pure tension and compression constitutive data. Effect of sintering on deformation behavior was significant, resulting in a dramatic change in constitutive relation. Cyclic fatigue testing of the coating material in tension-tension at room temperature showed an insignificant susceptibility to fatigue, similar to the slow crack growth behavior of the material in flexure in 800 C air.

Choi, Sung R.↗

Fracture, failure, and fragmentation

A general approach is required for describing matter of behavior when the failure is likely to involve growth and coalescence of a large number of fractures. Failures of this kind appear frequently in rapid dynamic processes, particularly in the formation of spall fragments. An approach to formulating constitutive relations that accounts for the opening, shear and growth of an ensemble of cracks is discussed. The approach accounts for plastic flow accompanying fragmentation. The resulting constitutive relations were incorporated into a Lagrangian computer program. A theoretical approach to coalescence is described. The simplest formulation uses a linear Liouville equation, with crack growth limited by the mean free path of cracks, assumed constant. This approach allows for an anisotropic distribution of cracks. An alternative approach in which the decrease of the mean free path with increasing crack size is accounted for, but the crack distribution is assumed isotropic is described. A reduction of the governing Liouville equation to an ordinary differential equation of third order is possible, and the result can be used to determine how mean free path decreases with increasing crack size.

Dienes, J. K.↗

Analysis of metal-matrix composite structures. I - Micromechanics constitutive theory. II - Laminate analyses

The micromechanical constitutive theory is used to examine the nonlinear behavior of continuous-fiber-reinforced metal-matrix composite structures. Effective lamina constitutive relations based on the Abouli micromechanics theory are presented. The inelastic matrix behavior is modeled by the unified viscoplasticity theory of Bodner and Partom. The laminate constitutive relations are incorporated into a first-order deformation plate theory. The resulting boundary value problem is solved by utilizing the finite element method. Attention is also given to computational aspects of the numerical solution, including the temporal integration of the inelastic strains and the spatial integration of bending moments. Numerical results the nonlinear response of metal matrix composites subjected to extensional and bending loads are presented.

Arenburg, R. T.↗

Crater modification by gravity - A mechanical analysis of slumping

An analysis of the stability of a crater from a mechanical point of view indicates that the observed slumping of lunar craters requires a perfectly plastic constitutive relation for the lunar surface rock. The angle of internal friction of this material must be less than a few degrees. The evidence for a perfectly plastic constitutive relation is described, and a simplified model of a crater in a perfectly plastic medium is used to investigate the nature of its collapse. Crater stability is found to depend principally upon a dimensionless parameter equal to the product of the density, acceleration of gravity, and depth divided by the yield strength. Criteria for stability, 'slope failures', and 'floor failures' are reported.

Melosh, H. J.↗

Seasonal and Interannual Variations of Ice Sheet Surface Elevation at the Summit of Greenland: Observed and Modeled

Observed seasonal and interannual variations in the surface elevation over the summit of the Greenland ice sheet are modeled using a new temperature-dependent formulation of firn-densification and observed accumulation variations. The observed elevation variations are derived from ERS (European Remote Sensing)-1 and ERS-2 radar altimeter data for the period between April 1992 and April 1999. A multivariate linear/sine function is fitted to an elevation time series constructed from elevation differences measured by radar altimetry at orbital crossovers. The amplitude of the seasonal elevation cycle is 0.25 m peak-to-peak, with a maximum in winter and a minimum in summer. Inter-annually, the elevation decreases to a minimum in 1995, followed by an increase to 1999, with an overall average increase of 4.2 cm a(exp -1) for 1992 to 1999. Our densification formulation uses an initial field-density profile, the AWS (automatic weather station) surface temperature record, and a temperature-dependent constitutive relation for the densification that is based on laboratory measurements of crystal growth rates. The rate constant and the activation energy commonly used in the Arrhenius-type constitutive relation for firn densification are also temperature dependent, giving a stronger temperature and seasonal amplitudes about 10 times greater than previous densification formulations. Summer temperatures are most important, because of the strong non-linear dependence on temperature. Much of firn densification and consequent surface lowering occurs within about three months of the summer season, followed by a surface build-up from snow accumulation until spring. Modeled interannual changes of the surface elevation, using the AWS measurements of surface temperature and accumulation and results of atmospheric modeling of precipitation variations, are in good agreement with the altimeter observations. In the model, the surface elevation decreases about 20 cm over the seven years due to more compaction driven by increasing summer temperatures. The minimum elevation in 1995 is driven mainly by a temporary accumulation decrease and secondarily by warmer temperatures. However, the overall elevation increase over the seven years is dominated by the accumulation increase in the later years.

Zwally, H. Jay↗

Constitutive Modeling of Superalloy Single Crystals and Directionally Solidified Materials

A unified viscoplastic constitutive relation based on crystallographic slip theory was developed for the deformation analysis of nickel base face centered cubic superalloy single crystals at elevated temperature. The single crystal theory is embedded in a self consistent method to derive a constitutive relation for a directionally solidified material comprised of a polycrystalline aggregate of columnar cylindrical grains. One of the crystallographic axes of the cylindrical crystals points in the columnar direction while the remaining crystallographic axes are oriented at random in the basal plane perpendicular to the columnar direction. These constitutive formulations are coded in FORTRAN for use in nonlinear finite element and boundary element programs.

Walker, K. P.↗

Reynolds-Averaged Navier-Stokes Computations of the NASA Juncture Flow Model Using FUN3D and OVERFLOW

Two Reynolds-averaged Navier-Stokes codes, FUN3D and OVERFLOW, are used to assess the capability of Spalart-Allmaras-based turbulence models to predict the flow over the NASA Juncture Flow model. Both free-air and in-tunnel simulations are performed. While the tunnel walls have some influence, it is found to be relatively minor in the juncture region of interest. Results from the two codes are found to be consistent with each other in attached flow regions, but results in the area of separation still show grid and code sensitivity, even on grids as large as 400 million unknowns. Nonetheless, it is possible to draw conclusions regarding the model capabilities. Without a quadratic constitutive relation, the linear model predicts a separation size that is too large. The inclusion of a nonlinear quadratic constitutive relation improves results significantly, but separation still occurs too far upstream. Prediction of turbulent normal stress differences play a key role in this flow. Although the nonlinear model makes better predictions in this regard, they could still be improved, particularly in the streamwise normal component near the wall.

C L Rumsey↗

Particle Distribution in Suspension Shear Flow

An earlier proposed constitutive relation for normal stresses originated by random particle fluctuations is used to describe a joint effect of thermal and shear-induced fluctuations on concentrational distributions in suspension flow. Averaged products of components of the fluctuation velocity are evaluated on a basis of the rational mechanics approach combined with a simple kinematic consideration. The equation of momentum conservation of the dispersed phase of a suspension closed with this constitutive relation is applied to unidirectional shear flow in the gravity field and to rotational Couette flow. Coupling of the thermal and shear-induced fluctuations results in that the ability of shear flow to suspend particles has a minimum at a certain particle size, all other things being equal. The developed model provides also for a reasonable explanation of particle distributions observed in Couette flow. The approach based on the consideration of momentum balance for the dispersed phase is proved to lead to an effective equation of convective diffusion of the suspended particles. Coefficients of mutual diffusion due to both thermal and shear-induced fluctuations are drastically different from corresponding self-diffusivities as regards both their scaling and their concentrational dependence.

Buyevich, Yu A.↗

A simple element for aeroelastic analysis of undamaged and damaged wings

The effects of material damage on the aeroelastic response of an anisotropic composite wing are investigated analytically. The wing is modeled as a Timoshenko beam of arbitrary cross section, applying a 24-DOF one-dimensional beam element and Cartesian coordinates. The formulations for the element stiffness matrix, the strain-displacement relations, the box-beam constitutive relations, the contributions of stringers and web, and the warping correction factor are outlined; the application of strip theory in deriving the aerodynamic model is briefly explained; and numerical results are presented in extensive graphs. It is shown that asymmetry and damage can induce bending-stretching coupling phenomena and thereby affect the stability of the aircraft either positively or negatively. The need to take such effects into account in the design of advanced composite wings is stressed.

Kapania, Rakesh K.↗

Thermo-elasto-viscoplastic analysis of problems in extension and shear

The problems of extension and shear behavior of structural elements made of carbon steel and subjected to large thermomechanical loads are investigated. The analysis is based on nonlinear geometric and constitutive relations, and is expressed in a rate form. The material constitutive equations are capable of reproducing all nonisothermal, elasto-viscoplastic characteristics. The results of the test problems show that: (1) the formulation can accommodate very large strains and rotations; (2) the model incorporates the simplification associated with rate-insensitive elastic response without losing the ability to model a rate-temperature dependent yield strength and plasticity; and (3) the formulation does not display oscillatory behavior in the stresses for the simple shear problem.

Riff, R.↗