Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nonlinear viscoelasticity”

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 91 records · Page 5

Elastoviscoplastic buckling behavior of simply supported columns

The elastoviscoplastic buckling behavior of simply supported columns is investigated. The analysis is based on an incremental formulation that can deal with nonlinear geometrical and constitutive relations. Both viscoelastic and viscoplastic models are discussed, and Bodner-Partom's constitutive equations are used in the viscoplastic case. A numerical example is given to demonstrate the buckling behavior of columns under various temperature and loading levels.

Song, Y.↗

Viscoelastic/damage modeling of filament-wound spherical pressure vessels

A model of the viscoelastic/damage response of a filament-wound spherical vessel used for long-term pressure containment is developed. The matrix material of the composite system is assumed to be linearly viscoelastic. Internal accumulated damage based upon a quadratic relationship between transverse modulus and maximum circumferential strain is postulated. The resulting nonlinear problem is solved by an iterative routine. The elastic-viscoelastic correspondence is employed to produce, in the Laplace domain, the associated elastic solution for the maximum circumferential strain which is inverted by the method of collocation to yield the time-dependent solution. Results obtained with the model are compared to experimental observations.

Hackett, Robert M.↗

A creep-rupture model of filament-wound spherical pressure vessels

The creep-rupture model is that of a quasi-isotropic filament-wound spherical pressure vessel, subjected to internal pressurization. The matrix material of the composite system is assumed to be linearly viscoelastic. Internal damage resulting from the relaxation of the matrix and the corresponding increase in microcracks is represented by a functional relationship between circumferential strain and transverse modulus. The numerical solution to this nonlinear problem is an iterative technique, whereby the elastic-viscoelastic correspondence principle is employed. In the Laplace domain, the associated elastic solution is obtained and this solution is inverted by the multidata method to yield the time-dependent solution.

Dozier, Jan D.↗

A kinematically driven anisotropic viscoelastic constitutive model applied to tires

Aircraft tires are composite structures manufactured with viscoelastic materials such as carbon black filled rubber and nylon cords. When loaded they experience large deflections and moderately large strains. Detailed structural models of tires require the use of either nonlinear shell or nonlinear three dimensional solid finite elements. Computational predictions of the dynamic response of tires must consider the composite viscoelastic material behavior in a realistic fashion. We describe a modification to a nonlinear anisotropic shell finite element so it can be used to model viscoelastic stresses during general deformations. The model is developed by introducing internal variables of the type used to model elastic strain energy. The internal variables are strains, curvatures, and transverse shear angles which are in a one-to-one correspondence with the generalized coordinates used to model the elastic strain energy for nonlinear response. A difference-relaxation equation is used to relate changes in the observable strain field to changes in the internal strain field. The internal stress state is introduced into the equilibrium equations by converting it to nodal loads associated with the element's displacement degrees of freedom. In this form the tangent matrix in the Newton-Raphson solution algorithm is not modified from its form for the nonlinear statics problem. Only the gradient vector is modified and the modification is not computationally costly. The existing finite element model for the Space Shuttle nose gear tire is used to provide examples of the algorithm. In the first example, the tire's rim is displaced at a constant rate up to a fixed value. In the second example, the tire's rim is enforced to follow a saw tooth load and unload curve to generate hysteresis loops.

Johnson, Arthur R.↗

Linear and nonlinear mechanical properties of a series of epoxy resins

The linear viscoelastic properties have been measured for a series of bisphenol-A-based epoxy resins cured with the diamine DDS. The linear viscoelastic master curves were constructed via time-temperature superposition of frequency dependent G-prime and G-double-prime isotherms. The G-double-prime master curves exhibited two sub-Tg transitions. Superposition of isotherms in the glass-to-rubber transition (i.e., alpha) and the beta transition at -60 C was achieved by simple horizontal shifts in the log frequency axis; however, in the region between alpha and beta, superposition could not be effected by simple horizontal shifts along the log frequency axis. The different temperature dependency of the alpha and beta relaxation mechanisms causes a complex response of G-double-prime in the so called alpha-prime region. A novel numerical procedure has been developed to extract the complete relaxation spectra and its temperature dependence from the G-prime and G-double-prime isothermal data in the alpha-prime region.

Curliss, D. B.↗

Egg yolk as a model for gelation: From rheometry to flow physics

Egg yolks are an excellent model for studying sol-gel transitions, particularly the power law viscoelasticity that defines the critical point of gelation. However, prior studies lack comprehensive datasets and fail to visualize flow behavior linked to temperature and time-dependent linear and nonlinear rheology. Here, we present a detailed dataset characterizing egg yolk viscoelasticity across temperature, time, and forcing amplitude using oscillatory shear, step strain, step stress, and constant high strain rate. Novel protorheology visualizations link rheological properties with observable flow behavior. Our findings highlight the nuanced determination of the critical gel point, emphasizing observation timescale dependencies. We compare methods to identify critical temperatures for gelation, including power law viscoelasticity, moduli crossover, diverging zero-shear viscosity, and emerging equilibrium elastic modulus, while visualizing flow consequences near these transitions. Egg yolk is an accessible, realistic, and nontoxic material relevant to the physicist and the chef alike, making it ideal for understanding the rheology of critical gels. By integrating protorheology photos and videos with rigorous rheometric data, we deepen the understanding of critical gels, with broader impacts for studying other materials with sol-gel transitions.

Marsh, Maxwell C. [Department of Mechanical Scienc↗

Interlaminar R-Curve Effects in Thermoplastic Composites: Experiments and Model Correlations

The interlaminar fracture behavior of thermoplastic composites under Mode-I and Mode-II static loading is investigated in this research. The material system under consideration is a carbon fiber reinforced low melt semi-crystalline resin TC1225 LMPAEK T700G (T700/LMPAEK)†. Characterization experiments were conducted to measure the load-displacement response, peak load, and propagation of delamination from a pre-implanted Kapton® insert. The Mode-I response was characterized using a double cantilever beam (DCB) test according to the ASTM D5528-22 standard. It was observed that the propagation exhibited significant resistance to delamination growth, i.e., an R-Curve effect. The test data indicated that the steady-state fracture toughness, G SS I , was approximately twice the mean initiation fracture toughness, G Ic . The Mode-II response of an end notch flexure (ENF) specimen tested according to the ASTM D7905-19 standard exhibited an unexpected stable behavior upon attaining peak load. Prior to attaining peak load, the Mode-II response exhibited some non-linearity, which may be attributed to nonlinear bulk material response, e.g., hyper-elasticity, viscoelasticity, plasticity, etc., and potentially friction at the load/support rollers and R-Curve effects. Finite element models based on the Cohesive Zone Method (CZM) and the Virtual Crack Closure Technique (VCCT) were used to predict the DCB and ENF specimen responses. DCB analysis results obtained using G SS I predicted propagation consistent with experimental data. However, results obtained with the as characterized G Ic were inaccurate in predicting delamination growth onset and peak load. ENF analysis results obtained with the as-characterized Mode-II fracture toughness, G IIc , were inaccurate in predicting the peak load and propagation. Hence, R-Curve effects were included in both DCB and ENF simulations, which enabled improved predictions. The experimentally determined Mode-I R-Curve was used for the DCB simulations, whereas an R-Curve was derived for Mode-II assuming that the change in compliance of the specimen was associated with crack growth only.

Prabhakar M Rao↗

Interlaminar R-Curve Effects in Thermoplastic Composites: Experiments and Model Correlations

The interlaminar fracture behavior of thermoplastic composites under Mode-I and Mode-II static loading is investigated in this research. The material system under consideration is a carbon fiber reinforced low melt semi-crystalline resin TC1225 LMPAEK T700G (T700/LMPAEK)†. Characterization experiments were conducted to measure the load-displacement response, peak load, and propagation of delamination from a pre-implanted Kapton® insert. The Mode-I response was characterized using a double cantilever beam (DCB) test according to the ASTM D5528-22 standard. It was observed that the propagation exhibited significant resistance to delamination growth, i.e., an R-Curve effect. The test data indicated that the steady-state fracture toughness, G SS I , was approximately twice the mean initiation fracture toughness, G Ic . The Mode-II response of an end notch flexure (ENF) specimen tested according to the ASTM D7905-19 standard exhibited an unexpected stable behavior upon attaining peak load. Prior to attaining peak load, the Mode-II response exhibited some non-linearity, which may be attributed to nonlinear bulk material response, e.g., hyper-elasticity, viscoelasticity, plasticity, etc., and potentially friction at the load/support rollers and R-Curve effects. Finite element models based on the Cohesive Zone Method (CZM) and the Virtual Crack Closure Technique (VCCT) were used to predict the DCB and ENF specimen responses. DCB analysis results obtained using G SS I predicted propagation consistent with experimental data. However, results obtained with the as characterized G Ic were inaccurate in predicting delamination growth onset and peak load. ENF analysis results obtained with the as-characterized Mode-II fracture toughness, G IIc , were inaccurate in predicting the peak load and propagation. Hence, R-Curve effects were included in both DCB and ENF simulations, which enabled improved predictions. The experimentally determined Mode-I R-Curve was used for the DCB simulations, whereas an R-Curve was derived for Mode-II assuming that the change in compliance of the specimen was associated with crack growth only.

Prabhakar Rao↗

Thermomechanical response of Gr/Pi composites

The effects of temperature changes upon the stresses and strains in composite laminates having carbon fibers in a polyimide matrix were evaluated. Composites having laminae in which there were nonlinear stress-strain relations for stresses transverse to the fibers and for axial shear stresses were treated. Material properties were considered to be temperature dependent. Separately, the effects of laminae viscoelastic response were also treated. The results suggest that, for this material, nonlinearities due either to stress or time dependent effects do not appear to be of major practical importance for conventional high temperature composite structures.

Rosen, B. W.↗

Thermoviscoelastic characterization and prediction of Kevlar/epoxy composite laminates

The thermoviscoelastic characterization of Kevlar 49/Fiberite 7714A epoxy composite lamina and the development of a numerical procedure to predict the viscoelastic response of any general laminate constructed from the same material were studied. The four orthotropic material properties, S sub 11, S sub 12, S sub 22, and S sub 66, were characterized by 20 minute static creep tests on unidirectional (0) sub 8, (10) sub 8, and (90) sub 16 lamina specimens. The Time-Temperature Superposition-Principle (TTSP) was used successfully to accelerate the characterization process. A nonlinear constitutive model was developed to describe the stress dependent viscoelastic response for each of the material properties. A numerical procedure to predict long term laminate properties from lamina properties (obtained experimentally) was developed. Numerical instabilities and time constraints associated with viscoelastic numerical techniques were discussed and solved. The numerical procedure was incorporated into a user friendly microcomputer program called Viscoelastic Composite Analysis Program (VCAP), which is available for IBM PC type computers. The program was designed for ease of use. The final phase involved testing actual laminates constructed from the characterized material, Kevlar/epoxy, at various temperatures and load level for 4 to 5 weeks. These results were compared with the VCAP program predictions to verify the testing procedure and to check the numerical procedure used in the program. The actual tests and predictions agreed for all test cases which included 1, 2, 3, and 4 fiber direction laminates.

Gramoll, K. C.↗

Physics augmented machine learning discovery of composition-dependent constitutive laws for 3D printed digital materials

Multi-material 3D printing, particularly through polymer jetting, enables the fabrication of digital materials by mixing distinct photopolymers at the micron scale within a single build to create a composite with tunable mechanical properties. Here, this work presents an integrated experimental and computational investigation into the composition-dependent mechanical behavior of 3D printed digital materials. We experimentally characterize five formulations, combining soft and rigid UV-cured polymers under uniaxial tension and torsion across three strain and twist rates. The results reveal nonlinear and rate-dependent responses that strongly depend on composition. To model this behavior, we develop a physics-augmented neural network (PANN) that combines a partially input convex neural network (pICNN) for learning the composition-dependent hyperelastic strain energy function with a quasi-linear viscoelastic (QLV) formulation for time-dependent response. The pICNN ensures convexity with respect to strain invariants while allowing non-convex dependence on composition. To enhance interpretability, we apply $L_0$ sparsification. For the time-dependent response, we introduce a multilayer perceptron (MLP) to predict viscoelastic relaxation parameters from composition. The proposed model accurately captures the nonlinear, rate-dependent behavior of 3D printed digital materials in both uniaxial tension and torsion, achieving high predictive accuracy for interpolated material compositions. This approach provides a scalable framework for automated, composition-aware constitutive model discovery for multi-material 3D printing.

Constitutive modeling↗

Stress-controlled medium-amplitude oscillatory shear (MAOStress) of PVA–Borax

We report the first-ever complete measurement of MAOStress material functions, which reveal that stress can be more fundamental than strain or strain rate for understanding linearity limits as a function of Deborah number. The material used is a canonical viscoelastic liquid with a single dominant relaxation time: polyvinyl alcohol (PVA) polymer solution cross-linked with tetrahydroborate (Borax) solution. We outline experimental limit lines and their dependence on geometry and test conditions. These MAOStress measurements enable us to observe the frequency dependence of the weakly nonlinear deviation as a function of stress amplitude. The observed features of MAOStress material functions are distinctly simpler than MAOStrain, where the frequency dependence is much more dramatic. The strain-stiffening transient network model was used to derive a model-informed normalization of the nonlinear material functions that accounts for their scaling with linear material properties. Moreover, we compare the frequency dependence of the critical stress, strain, and strain-rate for the linearity limit, which are rigorously computed from the MAOStress and MAOStrain material functions. While critical strain and strain-rate change by orders of magnitude throughout the Deborah number range, critical stress changes by a factor of about 2, showing that stress is a more fundamental measure of nonlinearity strength. This work extends the experimental accessibility of the weakly nonlinear regime to stress-controlled instruments and deformations, which reveal material physics beyond linear viscoelasticity but at conditions that are accessible to theory and detailed simulation.

Mechanics↗

Deformations and strains in a thick adherend lap joint

Displacement fields in a thick adherend lap joint were measured by high-sensitivity moire interferometry. Contour maps of in-plane U and V displacements were obtained across adhesive and adherend surfaces. Loads ranged from a modest load to a near-failure load. Quantitative results are given for displacements and strains in the adhesive and along the adhesive/adherend boundary lines. The results show nearly constant shear strain in the adhesive, nonlinear strains as a function of load or average shear stress, and viscoelastic or time-dependent response. Longitudinal normal strains in the adhesive are nearly two orders of magnitude less than the shear strains. With its subwavelength displacement resolution and high spatial resolution, moire interferometry is especially well suited for deformation studies of adhesive joints.

Post, D.↗

Calculating Dynamic Shear Moduli Of Polymers

Improved nonlinear mathematical model developed to fit experimental data on relaxation of stresses in viscoelastic materials. Developed specifically to study shear moduli of rubbery solid rocket propellants, also used to characterize polymers in general, other viscoelastic materials, and composites of viscoelastic materials. Facilitates and improves accuracy of analysis and numerical simulation of mechanical behavior of structural components (for example, tires) made of such materials.

Moore, Carleton J.↗

Quasi-static extension of a tensile crack contained in a viscoelastic-plastic solid

Final stretch criterion of failure is applied to the problem of quasi-static extension of a crack embedded in an elastic-plastic or viscoelastic-plastic matrix. The slow growth under subcritical conditions in a rate-sensitive Tresca solid is shown to be a superposition of creep rupture and McClintock's ductile growth. This type of growth occurs at subcritical magnitude of the imposed K-factor and can be accounted for only through a recognition of inelastic properties of solids. In the subcritical range there is no unique value for K sub c independent of geometrical configuration and flaw size. Not only the produced states of stress and strain are dependent on the loading path, but also the material resistance to fracture turns out to be a function of the history of loading that precedes catastrophic failure. A nonlinear integro-differential equation of motion is derived for a crack progressing through a viscoelastic medium with some limited ability to plastic flow. Examples of numerical integration are given incorporating both monotonic and cyclic loading programs.

Wnuk, M. P.↗

Strains and tilts on crustal blocks

Geodetic work done over the past century is examined to investigate block motion in areas of intense tectonic deformation, with special attention to geodetic releveling data obtained for Japan. Problems in interpreting strainmeter and tiltmeter to study the behavior of crustal blocks and block boundaries are discussed. Block dimensions of 5 to 50 km seem to occur frequently in regions of intense tectonic activity. Block boundaries are often considerably weaker than contiguous crustal blocks, resulting in a concentration of strains at boundaries and decrease of strains within blocks. Block to block variations in tilt phase and magnitude are observed. There is some evidence that block boundaries may exhibit strain-dependent elastic properties or respond viscoelastically, possibly accounting for the slow transmission of tectonic deformation reported in Japan. If nonlinear behavior is characteristic of regions fragmented by crustal blocks, it will generally not be possible to apply a site correction factor based on the observed distortion of tidal or seismic strains in interpreting secular strains.

Bilham, R. G.↗

Deformations and strains in adhesive joints by moire interferometry

Displacement fields in a thick adherend lap joint and a cracked lap shear specimen were measured by high sensitivity moire interferometry. Contour maps of in-plane U and V displacements were obtained across adhesive and adherent surfaces. Loading sequences ranged from modest loads to near-failure loads. Quantitative results are given for displacements and certain strains in the adhesive and along the adhesive/adherend boundary lines. The results show nonlinear displacements and strains as a function of loads or stresses and they show viscoelastic or time-dependent response. Moire interferometry is an excellent method for experimental studies of adhesive joint performance. Subwavelength displacement resolution of a few micro-inches, and spatial resolution corresponding to 1600 fringes/inch (64 fringes/mm), were obtained in these studies. The whole-field contour maps offer insights not available from local measurements made by high sensitivity gages.

Post, D.↗