SEARCH · Search NASA

Results for “linear 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 199 records · Page 11

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Nonlinear Architectures Can Alter the Dynamics of Polymer–Nanoparticle Composites

Polymer nanocomposites exhibit remarkable physical properties that are attractive for many applications. These systems have been so far investigated by using linear polymer chains; the role of polymer matrix architecture on local dynamics, bulk rheology and nanoparticle motion remains unexplored. Here, using quasielastic neutron scattering, bulk rheology and x-ray photon correlation spectroscopy, we investigated nanocomposites with spherical silica nanoparticles well dispersed in poly (ethylene oxide) matrices having different architectures (specifically linear, stars and hyperbranched). The results reveal that the mechanical reinforcement of the nanocomposites with the nonlinear polymers can be altered by orders of magnitude with respect to the conventional nanocomposite with the linear polymer. Additionally, the polymer compactness and inter-penetrability are found to play crucial roles in determining their bulk rheology. At microscopic level, average segmental dynamics is remarkably slowed down by the attractive NPs in the matrices of high degree of branching whereas no significant effect is observed in linear polymer matrix at the same NP loading. In addition, the nanoscale dynamics of particles in the compact nonlinear matrices exhibit strong decoupling from the bulk viscoelasticity, allowing their fast relaxation even at ≈ 30% by volume. These results provide an experimental evidence that macromolecular architecture is a powerful new tool for tuning the bulk rheological properties as well as the nanoscale dynamics of PNCs without the need for changing polymer molecular weight, nanoparticle size, shape, loading or dispersion state.

36 MATERIALS SCIENCE↗

NASA NTRS1993

Numerical simulation of the non-Newtonian mixing layer

This work is a continuing effort to advance our understanding of the effects of polymer additives on the structures of the mixing layer. In anticipation of full nonlinear simulations of the non-Newtonian mixing layer, we examined in a first stage the linear stability of the non-Newtonian mixing layer. The results of this study show that, for a fluid described by the Oldroyd-B model, viscoelasticity reduces the instability of the inviscid mixing layer in a special limit where the ratio (We/Re) is of order 1 where We is the Weissenberg number, a measure of the elasticity of the flow, and Re is the Reynolds number. In the present study, we pursue this project with numerical simulations of the non-Newtonian mixing layer. Our primary objective is to determine the effects of viscoelasticity on the roll-up structure. We also examine the origin of the numerical instabilities usually encountered in the simulations of non-Newtonian fluids.

Azaiez, Jalel↗

NASA NTRS2013

Granular Media-Based Tunable Passive Vibration Suppressor

and vibration suppression device is composed of statically compressed chains of spherical particles. The device superimposes a combination of dissipative damping and dispersive effects. The dissipative damping resulting from the elastic wave attenuation properties of the bulk material selected for the granular media is independent of particle geometry and periodicity, and can be accordingly designed based on the dissipative (or viscoelastic) properties of the material. For instance, a viscoelastic polymer might be selected where broadband damping is desired. In contrast, the dispersive effects result from the periodic arrangement and geometry of particles composing a linear granular chain. A uniform (monatomic) chain of statically compressed spherical particles will have a low-pass filter effect, with a cutoff frequency tunable as a function of particle mass, elastic modulus, Poisson fs ratio, radius, and static compression. Elastic waves with frequency content above this cutoff frequency will exhibit an exponential decay in amplitude as a function of propagation distance. System design targeting a specific application is conducted using a combination of theoretical, computational, and experimental techniques to appropriately select the particle radii, material (and thus elastic modulus and Poisson fs ratio), and static compression to satisfy estimated requirements derived for shock and/or vibration protection needs under particular operational conditions. The selection of a chain of polymer spheres with an elastic modulus .3 provided the appropriate dispersive filtering effect for that exercise; however, different operational scenarios may require the use of other polymers, metals, ceramics, or a combination thereof, configured as an array of spherical particles. The device is a linear array of spherical particles compressed in a container with a mechanism for attachment to the shock and/or vibration source, and a mechanism for attachment to the article requiring isolation (Figure 1). This configuration is referred to as a single-axis vibration suppressor. This invention also includes further designs for the integration of the single-axis vibration suppressor into a six-degree-of-freedom hexapod "Stewart"mounting configuration (Figure 2). By integrating each singleaxis vibration suppressor into a hexapod formation, a payload will be protected in all six degrees of freedom from shock and/or vibration. Additionally, to further enable the application of this device to multiple operational scenarios, particularly in the case of high loads, the vibration suppressor devices can be used in parallel in any array configuration.

Dillon, Robert P.↗

DOE OSTI2020

Viscoelastic Material Characterization and Modeling of Photovoltaic Module Packaging Materials for Direct Finite-Element Method Input

Numerical tools, such as the finite-element method, are increasingly used to design and evaluate the photovoltaic (PV) modules, providing for the reduction of development time and improved performance and reliability. However, high-fidelity material models are necessary to accurately model the complex structural behavior of the involved packaging materials. A common simplification used in recent years is to model the polymer materials (i.e., encapsulant and backsheet) as linear elastic, which will lead to inaccurate results. Therefore, in this work, we present a thorough characterization of the time- and temperature-dependent mechanical response of predominant PV module encapsulant and backsheet materials. Based on this material characterization, we developed and experimentally validated generalized Maxwell models to describe each material's viscoelastic response. In addition, we included measurements of the coefficient of thermal expansion and presented all material models in such a fashion for direct input into commercial finite-element method modeling software.

FEM↗

NASA NTRS1994

Modeling of periodic great earthquakes on the San Andreas fault: Effects of nonlinear crustal rheology

We analyze the cycle of great earthquakes along the San Andreas fault with a finite element numerical model of deformation in a crust with a nonlinear viscoelastic rheology. The viscous component of deformation has an effective viscosity that depends exponentially on the inverse absolute temperature and nonlinearity on the shear stress; the elastic deformation is linear. Crustal thickness and temperature are constrained by seismic and heat flow data for California. The models are for anti plane strain in a 25-km-thick crustal layer having a very long, vertical strike-slip fault; the crustal block extends 250 km to either side of the fault. During the earthquake cycle that lasts 160 years, a constant plate velocity v(sub p)/2 = 17.5 mm yr is applied to the base of the crust and to the vertical end of the crustal block 250 km away from the fault. The upper half of the fault is locked during the interseismic period, while its lower half slips at the constant plate velocity. The locked part of the fault is moved abruptly 2.8 m every 160 years to simulate great earthquakes. The results are sensitive to crustal rheology. Models with quartzite-like rheology display profound transient stages in the velocity, displacement, and stress fields. The predicted transient zone extends about 3-4 times the crustal thickness on each side of the fault, significantly wider than the zone of deformation in elastic models. Models with diabase-like rheology behave similarly to elastic models and exhibit no transient stages. The model predictions are compared with geodetic observations of fault-parallel velocities in northern and central California and local rates of shear strain along the San Andreas fault. The observations are best fit by models which are 10-100 times less viscous than a quartzite-like rheology. Since the lower crust in California is composed of intermediate to mafic rocks, the present result suggests that the in situ viscosity of the crustal rock is orders of magnitude less the rock viscosity determined in the laboratory.

Reches, Ze'ev↗

NASA NTRS1996

A viscoelastic higher-order beam finite element

A viscoelastic internal variable constitutive theory is applied to a higher-order elastic beam theory and finite element formulation. The behavior of the viscous material in the beam is approximately modeled as a Maxwell solid. The finite element formulation requires additional sets of nodal variables for each relaxation time constant needed by the Maxwell solid. Recent developments in modeling viscoelastic material behavior with strain variables that are conjugate to the elastic strain measures are combined with advances in modeling through-the-thickness stresses and strains in thick beams. The result is a viscous thick-beam finite element that possesses superior characteristics for transient analysis since its nodal viscous forces are not linearly dependent an the nodal velocities, which is the case when damping matrices are used. Instead, the nodal viscous forces are directly dependent on the material's relaxation spectrum and the history of the nodal variables through a differential form of the constitutive law for a Maxwell solid. The thick beam quasistatic analysis is explored herein as a first step towards developing more complex viscoelastic models for thick plates and shells, and for dynamic analyses. The internal variable constitutive theory is derived directly from the Boltzmann superposition theorem. The mechanical strains and the conjugate internal strains are shown to be related through a system of first-order, ordinary differential equations. The total time-dependent stress is the superposition of its elastic and viscous components. Equations of motion for the solid are derived from the virtual work principle using the total time-dependent stress. Numerical examples for the problems of relaxation, creep, and cyclic creep are carried out for a beam made from an orthotropic Maxwell solid.

Johnson, Arthur R.↗

NASA NTRS1987

Comparison of measured temperatures, thermal stresses and creep residues with predictions on a built-up titanium structure

Temperature, thermal stresses, and residual creep stresses were studied by comparing laboratory values measured on a built-up titanium structure with values calculated from finite-element models. Several such models were used to examine the relationship between computational thermal stresses and thermal stresses measured on a built-up structure. Element suitability, element density, and computational temperature discrepancies were studied to determine their impact on measured and calculated thermal stress. The optimum number of elements is established from a balance between element density and suitable safety margins, such that the answer is acceptably safe yet is economical from a computational viewpoint. It is noted that situations exist where relatively small excursions of calculated temperatures from measured values result in far more than proportional increases in thermal stress values. Measured residual stresses due to creep significantly exceeded the values computed by the piecewise linear elastic strain analogy approach. The most important element in the computation is the correct definition of the creep law. Computational methodology advances in predicting residual stresses due to creep require significantly more viscoelastic material characterization.

Jenkins, Jerald M.↗

DOE OSTI2024

Chain dynamics of alternating polymers P(C n EG 4 )

Alternating or sequence defined polymers attract the attention of an increasing number of researchers recently. Due to their different blocks, they are very customizable and material properties can be tuned. In this publication, we present dynamical studies with focus on polymer dynamics, investigated by rheology and fast field cycling (FFC) NMR. The molecular weight dependencies of the relaxation time and zero-shear viscosity could hint to entanglement effects; however, the spectral shape of the viscoelastic data resembles the polymer dynamics of unentangled melts. Furthermore, taking both techniques into account, i.e. rheology and FFC NMR leads to the conclusion that the relaxation spectra of the alternating polymer are similar to those of linear polymers.

36 MATERIALS SCIENCE↗

NASA NTRS1993

The analytical representation of viscoelastic material properties using optimization techniques

This report presents a technique to model viscoelastic material properties with a function of the form of the Prony series. Generally, the method employed to determine the function constants requires assuming values for the exponential constants of the function and then resolving the remaining constants through linear least-squares techniques. The technique presented here allows all the constants to be analytically determined through optimization techniques. This technique is employed in a computer program named PRONY and makes use of commercially available optimization tool developed by VMA Engineering, Inc. The PRONY program was utilized to compare the technique against previously determined models for solid rocket motor TP-H1148 propellant and V747-75 Viton fluoroelastomer. In both cases, the optimization technique generated functions that modeled the test data with at least an order of magnitude better correlation. This technique has demonstrated the capability to use small or large data sets and to use data sets that have uniformly or nonuniformly spaced data pairs. The reduction of experimental data to accurate mathematical models is a vital part of most scientific and engineering research. This technique of regression through optimization can be applied to other mathematical models that are difficult to fit to experimental data through traditional regression techniques.

Hill, S. A.↗

NASA NTRS1985

Transient polar motions and the nature of the asthenosphere for short time scales

A uniformly valid mathematical formalism is developed to study the secular motions of the rotational axis of a layered viscoelastic earth due to seismic excitation. The changes required for implementing the formulation within the framework of the faulting problem. The rationale of adopting the chosen nrheological model, which contains a low-viscosity zone beneath the lithosphere and is based on linear Maxell constitutive relationship, is discussed. The impact of this low-viscosity channel on thhe two families of relaxation time, governing both isostatic readjustment and rotational processes, is considered. It is found that the polar motions depend sensitively on the viscosity structure of the asthenosphere and not at all on the underlying mantle. A gloal low-velocity zone with short-term asthenospheric viscosities less than about 5 x 10 to the 18th Pa-s and widths greater than 50 km is ruled out.

Boschi, E.↗

DOE OSTI2023

Bead–Spring Simulation of Ionomer Melts—Studying the Effects of Chain-Length and Associating Group Fraction on Equilibrium Structure and Extensional Flow Behavior

Ionomers are associative polymers with diverse applications ranging from selective membranes and high-performance adhesives to abrasion- and chemical-resistant coatings, insulation layers, vacuum packaging, and foamed sheets. Within equilibrium melt, the ionic or associating groups are known to form thermally reversible, associative clusters whose presence can significantly affect the system’s mechanical, viscoelastic, and transport properties. It is, thus, of great interest to understand how to control such clusters’ size distribution, shape, and stability through the designed choice of polymer architecture and the ionic groups’ fraction, arrangement, and interaction strength. In this work, we represent linear associating polymers using a Kremer–Grest type bead–spring model and perform large-scale MD simulations to explore the effect of polymer chain-length (l) and fraction (f s ) of randomly placed associating groups on the size distribution and stability of formed clusters. We consider different chain-lengths (below and above entanglement), varying fractions of associating groups (represented by ‘sticky’ beads) between 5 and 20%, and a fixed sticky–sticky nonbond interaction strength of four times that between regular non-associating beads. For all melts containing associating groups the equilibrium structure factor S(q) displays a signature ionomer peak at low wave vector q whose intensity increases with increasing f s and l. The average cluster size Nc increases with f s . However, the effect of chain-length on Nc appears to be pronounced only at higher values of f s . Under extensional flows, the computed stress (and viscosity) is higher at higher f s and l regardless of strain rate. Beyond a critical strain rate, we observe fragmentation of the associative clusters, which has interesting effects on the stress/viscous response.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Range and strength of mechanical interactions of force dipoles in elastic fiber networks

Mechanical forces generated by myosin II molecular motors drive diverse cellular processes, most notably shape change, division and locomotion. These forces may be transmitted over long range through the cytoskeletal medium – a disordered, viscoelastic network of biopolymers. The resulting cell size scale force chains can in principle mediate mechanical interactions between distant actomyosin units, leading to self-organized structural order in the cell cytoskeleton. Inspired by such force transmission through elastic structures in the cytoskeleton, we consider a percolated fiber lattice network, where fibers are represented as linear elastic elements that can both bend and stretch, and the contractile activity of myosin motors is represented by force dipoles. Then, by using a variety of metrics, we show how two such contractile force dipoles interact with each other through their mutual mechanical deformations of the elastic fiber network. As a prelude to two-dipole interactions, we quantify how forces propagate through the network from a single anisotropic force dipole by analyzing clusters of nodes connected by highly strained bonds, as well as through the decay rate of strain energy with distance from a force dipole. We show that predominant fiber bending screens out force propagation, resulting in reduced and strongly network configuration-dependent dipole interactions. On the other hand, stretching-dominated networks support longer-ranged inter-dipole interactions that recapitulate the predictions of linear elasticity theory. By characterizing the differences between tensile and compressive force propagation in the fiber network, we show how inter-dipole interaction depends on the dipoles’ mutual separation and orientation. The resulting elastic interaction energy may mediate a force between multiple distant dipoles, leading to their self-organization into ordered configurations. In conclusion, this provides a potential pathway for active mechanical force-driven structural order in elastic biopolymer networks.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS1987

Rate-dependent mode I interlaminar crack growth mechanisms in graphite/epoxy and graphite/PEEK

In this paper the mode I fracture behavior of graphite/epoxy and graphite/PEEK composites is examined over four decades of crosshead rates (0.25-250 mm/min). Straight-sided double-cantilever-beam specimens consisting of unidirectional laminates were tested at room temperature. For graphite/epoxy the load-deflection response was linear to fracture, and stable slow crack growth initiating at the highest load level was observed for all rates tested. In contrast, mode I crack growth in the graphite/PEEK material was often unstable and showed stick-slip behavior. Subcritical crack growth occurring prior to the onset of fracture was observed at intermediate displacement rates. A mechanism for the fracture behavior of the graphite/PEEK material (based on viscoelastic, plastic, and microcrack coalescence in the process zone) is proposed and related to the observed rate-dependent phenomena.

Gillespie, J. W., Jr.↗

DOE OSTI2020

Periodic training of creeping solids

Significance It is well appreciated that many disordered materials deform their shape irreversibly (plastically) under an external load (e.g., memory foam). Here, we show that this plasticity can be exploited to train materials to develop novel elastic responses by straining them periodically. By applying different periodic strains to a common viscoelastic material, we are able to design a number of different responses. These include a maximally negative Poisson’s ratio, bistable behavior, and nonlocal bond-specific responses. In contrast to computer-aided design, we rely on plasticity to self-organize the system in response to local stresses. This approach shows promise to achieve an unprecedented control over behavior at large strains well beyond the linear-response regime.

36 MATERIALS SCIENCE↗

NASA NTRS1988

Optimal discrete-time LQR problems for parabolic systems with unbounded input: Approximation and convergence

An abstract approximation and convergence theory for the closed-loop solution of discrete-time linear-quadratic regulator problems for parabolic systems with unbounded input is developed. Under relatively mild stabilizability and detectability assumptions, functional analytic, operator techniques are used to demonstrate the norm convergence of Galerkin-based approximations to the optimal feedback control gains. The application of the general theory to a class of abstract boundary control systems is considered. Two examples, one involving the Neumann boundary control of a one-dimensional heat equation, and the other, the vibration control of a cantilevered viscoelastic beam via shear input at the free end, are discussed.

Rosen, I. G.↗

NASA NTRS1996

Expendable Launch Vehicle Studies

Analytical support studies of expendable launch vehicles concentrates on the stability of the dynamics during launch especially during or near the region of maximum dynamic pressure. The in-plane linearized dynamic equations of a generic launch vehicle with multiple flexible bending and fuel sloshing modes are developed. The design of a robust LQR controller based on the reduced order system is accomplished using the parameter perturbation technique. The ELV modeling and analysis team has spent the past three years working on the theoretical development and application of sensitivity analysis to solve eigenvalue problems associated with structural dynamics. Specific application areas include stochastic vibrations, viscously damped vibrations and nonlinear dynamics. In stochastic linear vibrations, sensitivity analysis methods have been developed which determined the expectation and variance of the response for given probability density functions of various stochastic parameters. Examples worked include beam vibrations with up to six stochastic parameters; and the eigensensitivity results differ little from those obtained, with much effect, from Monte Carlo techniques. For viscously damped vibrations, eigensensitivity analysis has given excellent results for homogeneous beams, modeled by either quadratic or quartic eigenvalue equations. Specific applications include Kevin and Maxwell-type viscoelastic beams.

Bainum, P. M.↗

NASA NTRS1985

Ultrasonic characterization of changes in viscoelastic properties of epoxy during cure

A technique for using the longitudinal velocity (LV) of an ultrasonic wave to monitor the extent of cross linking (CL) during the cure of thermosetting resins is described. The method was developed by monitoring the rate of change in LV during the cure of a bisphenol-A epoxy resin with an amine adduct. The experiment included variations in the temperature and stoichiometry in order to express the rate of change in terms of the reaction kinetics. The pulse-echo method was used with a single transducer operating at 20 MHz. Numerical models were defined to account for the acoustic response of a single layer, the attenuation and the reflection coefficient. A linear relationship was observed between the inverse of the temperature and the log of the rate of change in the velocity, supporting the theory that the velocity could be used to monitor the extent of the cross-linking reaction. An activation energy of 11.9 kcal/mole was calculated for the mixture being investigated.

Winfree, W. P.↗

DOE OSTI2026

Nonlinear viscoelastic response of silicone additively manufactured direct ink write (DIW) foams under repetitive compression

To investigate dynamic fatigue behavior of foam in military protective applications, such as helmets, additively manufactured (AM) foams were compressively strained into the plateau region using a reduced design of experiments. A simple power law was found to govern the decline in dynamic stiffness (complex modulus) as the foams underwent the purchase order requirement of 10,000 cycles of small deformation in the plateau region. This rate of decline was newly found to correlate with the degree of nonlinearity in the material’s deformation, quantified using total harmonic distortion. Materials with low nonlinearity exhibited relatively stable stiffness across cycles, while those with high nonlinearity experienced greater losses. The observed nonlinearity depended on both applied stress and strain rate. A strong linear correlation (R 2 = 0.78) was identified between second-order nonlinearity and the time-dependent stiffness response. Two lattice structures were examined: face-centered tetragonal (FCT) and simple cubic (SC). The SC material exhibited higher total harmonic distortion (5%) and lower stiffness retention than the FCT (2%). In conclusion, these results suggest that for cyclic compression applications in a wide variety of industries such as packaging, personal protective equipment, or aerospace, selecting materials with lower stress and greater structural uniformity can enhance the stability of dynamic performance.

42 ENGINEERING↗