SEARCH · Search NASA
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.
Theoretical paper presents solutions of equations that describe polymeric aging. Solutions apply under such loading conditions as constant strain (stress relaxation), constant strain rate, and stress relaxation with sudden change of cross-link density. Theoretical models form framework for predicting chemomechanical aging behavior of elastomers and polymers.
Explore the source record for details and available documents.
A new mathematical approach for quantifying the mechanical properties of elastomeric materials under biaxial loads is presented. Specific equations relating principal strains and principal true stresses for a homogeneous, isotropic, and elastic material are proposed that resemble the conventional Hooke's law of linear elasticity. The predicted stresses are compared to those from three different sets of experimental data and to stresses from three different theories based on Rivlin's work-of-deformation approach. The proposed approach is considered to be of benefit to design engineers involved in a broad range of rubber products.
Abstract not provided.
Applying smart materials in aeropropulsion systems may improve the performance of aircraft engines through a variety of vibration, noise, and shape-control applications. To facilitate the experimental characterization of these smart structures, researchers have been focusing on developing analytical models to account for the coupled mechanical, electrical, and thermal response of these materials. One focus of current research efforts has been directed toward incorporating a comprehensive thermal analysis modeling capability. Typically, temperature affects the behavior of smart materials by three distinct mechanisms: Induction of thermal strains because of coefficient of thermal expansion mismatch 1. Pyroelectric effects on the piezoelectric elements; 2. Temperature-dependent changes in material properties; and 3. Previous analytical models only investigated the first two thermal effects mechanisms. However, since the material properties of piezoelectric materials generally vary greatly with temperature (see the graph), incorporating temperature-dependent material properties will significantly affect the structural deflections, sensory voltages, and stresses. Thus, the current analytical model captures thermal effects arising from all three mechanisms through thermopiezoelectric constitutive equations. These constitutive equations were incorporated into a layerwise laminate theory with the inherent capability to model both the active and sensory response of smart structures in thermal environments. Corresponding finite element equations were formulated and implemented for both the beam and plate elements to provide a comprehensive thermal effects modeling capability.
A constitutive equation was developed to predict the inelastic thermomechanical response of single crystal turbine blades. These equations are essential for developing accurate finite element models of hot section components and contribute significantly to the understanding and prediction of crack initiation and propagation. The method used was limited to unified state variable constitutive equations. Two approaches to developing an anisotropic constitutive equation were reviewed. One approach was to apply the Stouffer-Bodner representation for deformation induced anisotropy to materials with an initial anisotropy such as single crystals. The second approach was to determine the global inelastic strain rate from the contribution of the slip in each of the possible crystallographic slip systems. A three dimensional finite element is being developed with a variable constitutive equation link that can be used for constitutive equation development and to predict the response of an experiment using the actual specimen geometry and loading conditions.
In this work, we investigate the applicability and relative merit of discrete and continuous versions of physics-informed neural network (PINN) methods for learning unknown governing equations or constitutive relations in a nonlinear dynamical system. In the case of unknown dynamics, entire right-hand-side (RHS) equations of the ordinary differential equations are unknown. In the case of unknown constitutive relations, however, the RHS equations are known up to the specification of constitutive relations (that may depend on the state of the system). We use a deep neural network to model unknown governing equations or constitutive relations. The discrete PINN approach combines classical multistep discretization methods for dynamical systems with neural-network-based machine learning methods. On the other hand, the continuous versions utilize deep neural networks to minimize the residual function for the continuous governing equations. We use the case of a fedbatch bioreactor system to study the effectiveness of these approaches and discuss conditions for their applicability. Our results indicate that the accuracy of the trained neural network models is much higher for the cases where we only have to learn a constitutive relation instead of all dynamics. This finding corroborates the well-known fact from scientific computing that building as much structural information as is available into an algorithm can enhance its efficiency and/or accuracy.
The applicability of a thermodynamic constitutive theory of deformation to the prediction of primary creep and creep strain relaxation behavior in metals is examined. Constitutive equations derived from the theory are subjected to a parametric analysis in order to determine the influence of several parameters on the curve forms generated by the equations. A computer program is developed which enables the solution of a generalized constitutive equation using experimental data as input. Several metals were tested to form a data base of primary creep and relaxation behavior. The extent to which these materials conformed to the constitutive equation showed wide variability, with the alloy Ti-6Al-4V exhibiting the most consistent results. Accordingly, most of the analysis is concentrated upon data from that alloy, although creep and relaxation data from all the materials tested are presented. Experimental methods are outlined as well as some variations in methods of analysis. Various theoretical and practical implications of the work are discussed.
Potential gas turbine applications will expose polymer matrix composites to very high strain rate loading conditions, requiring an ability to understand and predict the material behavior under extreme conditions. Specifically, analytical methods designed for these applications must have the capability of properly capturing the strain rate sensitivities and nonlinearities that are present in the material response. The Ramaswamy-Stouffer constitutive equations, originally developed to analyze the viscoplastic deformation of metals, have been modified to simulate the nonlinear deformation response of ductile, crystalline polymers. The constitutive model is characterized and correlated for two representative ductile polymers. Fiberite 977-2 and PEEK, and the computed results correlate well with experimental values. The polymer constitutive equations are implemented in a mechanics of materials based composite micromechanics model to predict the nonlinear, rate dependent deformation response of a composite ply. Uniform stress and uniform strain assumptions are applied to compute the effective stresses of a composite unit cell from the applied strains. The micromechanics equations are successfully verified for two polymer matrix composites. IM7/977-2 and AS4/PEEK. The ultimate strength of a composite ply is predicted with the Hashin failure criteria that were implemented in the composite micromechanics model. The failure stresses of the two composite material systems are accurately predicted for a variety of fiber orientations and strain rates. The composite deformation model is implemented in LS-DYNA, a commercially available transient dynamic explicit finite element code. The matrix constitutive equations are converted into an incremental form, and the model is implemented into LS-DYNA through the use of a user defined material subroutine. The deformation response of a bulk polymer and a polymer matrix composite are predicted by finite element analyses. The results compare reasonably well to experimental values, with some discrepancies. The discrepancies are at least partially caused by the method used to integrate the rate equations in the polymer constitutive model.
The material strength and strain rate effects associated with the hypervelocity impact problem were considered. A yield criterion involving the second and third invariants of the stress deviator and a strain rate sensitive constitutive equation were developed. The part of total deformation which represents change in shape is attributable to the stress deviator. Constitutive equation is a means for analytically describing the mechanical response of a continuum under study. The accuracy of the yield criterion was verified utilizing the published two and three dimensional experimental data. The constants associated with the constitutive equation were determined from one dimensional quasistatic and dynamic experiments. Hypervelocity impact experiments were conducted on semi-infinite targets of 1100 aluminum, 6061 aluminum alloy, mild steel, and commercially pure lead using spherically shaped and normally incident pyrex projectiles.
The limiting shear stress shear rheological model was applied to property measurements pursuant to the use of the constitutive equation and the application of the constitutive equation to elastrohydrodynamic (EHD) traction. Experimental techniques were developed to subject materials to isothermal compression which is similar to the history the materials were subjected to in EHD contacts. In addition, an apparatus was developed for measuring the shear stress-strain behavior of solid lubricating materials. Four commercially available materials were examined under pressure. They exhibit elastic and limiting shear stress behavior similar to that of liquid lubricants. The application of the limiting shear stress model to traction predictions was extended employing the primary materials properties measured in the laboratory. The shear rheological model was also applied to a Grubin-like EHD inlet analysis for predicting film thicknesses when employing the limiting shear stress model material behavior.
Fraction-order viscoelastic (FOV) material models have been proposed and studied in 1D since the 1930's, and were extended into three dimensions in the 1970's under the assumption of infinitesimal straining. It was not until 1997 that Drozdov introduced the first finite-strain FOV constitutive equations. In our presentation, we shall continue in this tradition by extending the standard, FOV, fluid and solid, material models introduced in 1971 by Caputo and Mainardi into 3D constitutive formula applicable for finite-strain analyses. To achieve this, we generalize both the convected and co-rotational derivatives of tensor fields to fractional order. This is accomplished by defining them first as body tensor fields and then mapping them into space as objective Cartesian tensor fields. Constitutive equations are constructed using both variants for fractional rate, and their responses are contrasted in simple shear. After five years of research and development, we now possess a basic suite of numerical tools necessary to study finite-strain FOV constitutive equations and their iterative refinement into a mature collection of material models. Numerical methods still need to be developed for efficiently solving fraction al-order integrals, derivatives, and differential equations in a finite element setting where such constitutive formulae would need to be solved at each Gauss point in each element of a finite model, which can number into the millions in today's analysis.
A previously developed analytical formulation has been modified in order to more accurately account for the effects of hydrostatic stresses on the nonlinear, strain rate dependent deformation of polymer matrix composites. State variable constitutive equations originally developed for metals have been modified in order to model the nonlinear, strain rate dependent deformation of polymeric materials. To account for the effects of hydrostatic stresses, which are significant in polymers, the classical J2 plasticity theory definitions of effective stress and effective inelastic strain, along with the equations used to compute the components of the inelastic strain rate tensor, are appropriately modified. To verify the revised formulation, the shear and tensile deformation of two representative polymers are computed across a wide range of strain rates. Results computed using the developed constitutive equations correlate well with experimental data. The polymer constitutive equations are implemented within a strength of materials based micromechanics method to predict the nonlinear, strain rate dependent deformation of polymer matrix composites. The composite mechanics are verified by analyzing the deformation of a representative polymer matrix composite for several fiber orientation angles across a variety of strain rates. The computed values compare well to experimentally obtained results.
Piezoelectric material inherently possesses coupling between electrostatics and structural dynamics. Utilizing linear piezoelectric theory results in an intrinsically coupled pair of piezoelectric constitutive equations. One equation describes the direct piezoelectric effect where strains produce an electric field and the other describes the converse effect where an applied electrical field produces strain. The purpose of this study is to compare finite element analysis and experiments of a thin plate with bonded piezoelectric material. Since an isotropic plate in combination with a thin piezoelectric layer constitutes a special case of a laminated composite, the classical laminated plate theory is used in the formulation to accommodated generic laminated composite panels with multiple bonded and embedded piezoelectric layers. Additionally, the von Karman large deflection plate theory is incorporated. The formulation results in laminate constitutive equations that are amiable to the inclusion of the piezoelectric constitutive equations yielding in a fully electro-mechanically coupled composite laminate. Using the finite element formulation, the governing differential equations of motion of a composite laminate with embedded piezoelectric layers are derived. The finite element model not only considers structural degrees of freedom (d.o.f.) but an additional electrical d.o.f. for each piezoelectric layer. Comparison between experiment and numerical prediction is performed by first treating the piezoelectric as a sensor and then again treating it as an actuator. To assess the piezoelectric layer as a sensor, various uniformly distributed pressure loads were simulated in the analysis and the corresponding generated voltages were calculated using both linear and nonlinear finite element analyses. Experiments were carried out by applying the same uniformly distributed loads and measuring the resulting generated voltages and corresponding maximum plate deflections. It is found that a highly nonlinear relationship exists between maximum deflection and voltage versus pressure loading. In order to assess comparisons of predicted and measured piezoelectric actuation, sinusoidal excitation voltages are simulated/applied and maximum deflections are calculated/measured. The maximum deflection as a function of time was determined using the linear finite elements analysis. Good correlation between prediction and measurement was achieved in all cases.
A number of titanium matrix composite (TMC) systems are currently being investigated for high-temperature air frame and propulsion system applications. As a result, numerous computational methodologies for predicting both deformation and life for this class of materials are under development. An integral part of these methodologies is an accurate and computationally efficient constitutive model for the metallic matrix constituent. Furthermore, because these systems are designed to operate at elevated temperatures, the required constitutive models must account for both time-dependent and time-independent deformations. To accomplish this, the NASA Lewis Research Center is employing a recently developed, complete, potential-based framework. This framework, which utilizes internal state variables, was put forth for the derivation of reversible and irreversible constitutive equations. The framework, and consequently the resulting constitutive model, is termed complete because the existence of the total (integrated) form of the Gibbs complementary free energy and complementary dissipation potentials are assumed a priori. The specific forms selected here for both the Gibbs and complementary dissipation potentials result in a fully associative, multiaxial, nonisothermal, unified viscoplastic model with nonlinear kinematic hardening. This model constitutes one of many models in the Generalized Viscoplasticity with Potential Structure (GVIPS) class of inelastic constitutive equations.
The combination of thermal and mechanical loading expected in practice means that constitutive equations of metal matrix composites must be developed which deal with time-independent and time-dependent irreversible deformation. Also, the internal state of composites is extremely complicated which underlines the need to formulate macroscopic constitutive equations with a limited number of state variables which represent the internal state at the micro level. One available method for calculating the macro properties of composites in terms of the distribution and properties of the constituent materials is the method of homogenization whose formulation is based on the periodicity of the substructure of the composite. A homogenization procedure was developed which lends itself to the use of the finite element procedure. The efficiency of these procedures, to determine the macroscopic properties of a composite system from its constituent properties, was demonstrated utilizing an aluminum plate perforated by directionally oriented slits. The selection of this problem is based on the fact that, extensive experimental results exist, the macroscopic response is highly anisotropic, and that the slits provide very high stress gradients which severely test the effectiveness of the computational procedures. Furthermore, both elastic and plastic properties were investigated so that the application to practical systems with inelastic deformation should be able to proceed without difficulty. The effectiveness of the procedures was rigorously checked against experimental results and with the predictions of approximate calculations. Using the computational results it is illustrated how macroscopic constitutive equations can be expressed in forms of the elastic and limit load behavior.
The results presented here are part of an ongoing research program, to develop strain rate dependent deformation and failure models for the analysis of polymer matrix composites subject to high strain rate impact loads. A micromechanics approach is employed in this work, in which state variable constitutive equations originally developed for metals have been modified to model the deformation of the polymer matrix, and a strength of materials based micromechanics method is used to predict the effective response of the composite. In the analysis of the inelastic deformation of the polymer matrix, the definitions of the effective stress and effective inelastic strain have been modified in order to account for the effect of hydrostatic stresses, which are significant in polymers. Two representative polymers, a toughened epoxy and a brittle epoxy, are characterized through the use of data from tensile and shear tests across a variety of strain rates. Results computed by using the developed constitutive equations correlate well with data generated via experiments. The procedure used to incorporate the constitutive equations within a micromechanics method is presented, and sample calculations of the deformation response of a composite for various fiber orientations and strain rates are discussed.
The formulation of rheological constitutive equations—models that relate internal stresses and deformations in complex fluids—is a critical step in the engineering of systems involving soft materials. While data-driven models provide accessible alternatives to expensive first-principles models and less accurate empirical models in many engineering disciplines, the development of similar models for complex fluids has lagged. The diversity of techniques for characterizing non-Newtonian fluid dynamics creates a challenge for classical machine learning approaches, which require uniformly structured training data. Consequently, early machine-learning based constitutive equations have not been portable between different deformation protocols or mechanical observables. Here, we present a data-driven framework that resolves such issues, allowing rheologists to construct learnable models that incorporate essential physical information, while remaining agnostic to details regarding particular experimental protocols or flow kinematics. These scientific machine learning models incorporate a universal approximator within a materially objective tensorial constitutive framework. By construction, these models respect physical constraints, such as frame-invariance and tensor symmetry, required by continuum mechanics. We demonstrate that this framework facilitates the rapid discovery of accurate constitutive equations from limited data and that the learned models may be used to describe more kinematically complex flows. This inherent flexibility admits the application of these “digital fluid twins” to a range of material systems and engineering problems. We illustrate this flexibility by deploying a trained model within a multidimensional computational fluid dynamics simulation—a task that is not achievable using any previously developed data-driven rheological equation of state.