Search NASA⌕ Search

SEARCH · Search NASA

Results for “viscoelastic modeling”

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 235 records · Page 13

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↗

The Extensional Rheology of Non-Newtonian Materials

The evolution of the transient extensional stresses in dilute and semi-dilute viscoelastic polymer solutions are measured with a filament stretching rheometer of a design similar to that first introduced by Sridhar, et al. The solutions are polystyrene-based (PS) Boger fluids that are stretched at constant strain rates ranging from 0.6 less than or equal to epsilon(0) less than or equal to 4s(exp -1) and to Hencky strains of epsilon greater than 4. The test fluids all strain harden and Trouton ratios exceeding 1000 are obtained at high strains. The experimental data strain hardens at lower strain levels than predicted by bead-spring FENE models. In addition to measuring the transient tensile stress growth, we also monitor the decay of the tensile viscoelastic stress difference in the fluid column following cessation of uniaxial elongation as a function of the total imposed Hencky strain and the strain rate. The extensional stresses initially decay very rapidly upon cessation of uniaxial elongation followed by a slower viscoelastic relaxation, and deviate significantly from FENE relaxation predictions. The relaxation at long times t is greater than or equal to 5 s, is compromised by gravitational draining leading to non-uniform filament profiles. For the most elastic fluids, partial decohension of the fluid filament from the endplates of the rheometer is observed in tests conducted at high strain rates. This elastic instability is initiated near the rigid endplate fixtures of the device and it results in the progressive breakup of the fluid column into individual threads or 'fibrils' with a regular azimuthal spacing. These fibrils elongate and bifurcate as the fluid sample is elongated further. Flow visualization experiments using a modified stretching device show that the instability develops as a consequence of an axisymmetry-breaking meniscus instability in the nonhomogeneous region of highly deformed fluid near the rigid endplate.

Spiegelberg, Stephen H.↗

Reversible and Irreversible Time-Dependent Behavior of GRCop-84

A series of mechanical tests were conducted on a high-conductivity copper alloy, GRCop-84, in order to understand the time dependent response of this material. Tensile, creep, and stress relaxation tests were performed over a wide range of temperatures, strain rates, and stress levels to excite various amounts of time-dependent behavior. At low applied stresses the deformation behavior was found to be fully reversible. Above a certain stress, termed the viscoelastic threshold, irreversible deformation was observed. At these higher stresses the deformation was observed to be viscoplastic. Both reversible and irreversible regions contained time dependent deformation. These experimental data are documented to enable characterization of constitutive models to aid in design of high temperature components.

viscoelasticity↗

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.↗

Vertical movements following a dip-slip earthquake

The results of finite element (FEM) computations based on elastic dislocation theory are presented. It is shown that the ratio of fault depth to the lithospheric thickness is a significant factor in the sign of the rebound. Subsidence occurs if the fault penetrates less than approximately 0.68 of the way through the lithosphere. An uplift occurs if the fault cuts deeper. The FEM employed consisted of a 546 node, 500 element grid and a quasi-plane strain algorithm. The situation modeled consisted of a 30 deg dip fault in an initially planar elastic lithosphere overlying a Maxwell viscoelastic asthenosphere. It is suggested that postseismic rebound may be measurable over a several year interval, thus providing data on lithosphere thicknesses.

Melosh, H. J.↗

Bend stress relaxation and tensile primary creep of a polycrystalline alpha-SiC fiber

Understanding the thermomechanical behavior (creep and stress relaxation) of ceramic fibers is of both practical and basic interest. On the practical level, ceramic fibers are the reinforcement for ceramic matrix composites which are being developed for use in high temperature applications. It is important to understand and model the total creep of fibers at low strain levels where creep is predominantly in the primary stage. In addition, there are many applications where the component will only be subjected to thermal strains. Therefore, the stress relaxation of composite consituents in such circumstances will be an important factor in composite design and performance. The objective of this paper is to compare and analyze bend stress relaxation and tensile creep data for alpha-SiC fibers produced by the Carborundum Co. (Niagara Falls, NY). This fiber is of current technical interest and is similar in composition to bulk alpha-SiC which has been studied under compressive creep conditions. The temperature, time, and stress dependences will be discussed for the stress relaxation and creep results. In addition, some creep and relaxation recovery experiments were performed in order to understand the complete viscoelastic behavior, i.e. both recoverable and nonrecoverable creep components of these materials. The data will be presented in order to model the deformation behavior and compare relaxation and/or creep behavior for relatively low deformation strain conditions of practical concern. Where applicable, the tensile creep results will be compared to bend stress relaxation data.

Hee Man, Yun↗

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.↗

Modeling and Simulation of Vehicle Dynamics on the Surface of Phobos

In this paper we analyze the dynamics of a spacecraft in proximity of Phobos by developing the equations of motion of a test mass in the Phobos rotating frame using a model based on circularly-restricted three body problem, and by analyzing the dynamics of a ATHLETE hopper vehicle interacting with the soil under different soil-interaction conditions. The main conclusion of the numerical studies is that the system response is dominated by the stiffness and damping parameters of the leg springs, with the soil characteristics having a much smaller effect. The system simulations identify ranges of parameters for which the vehicle emerges stably (relying only on the passive viscoelastic damper at each leg) or unstably (needing active attitude control) from the hop. The implication is that further experimental and possibly computational modeling work, as well as site characterization (from precursor missions) will be necessary to obtain validated performance models.

soil-structure interaction↗

Wave transmission characteristics and anisotropy of canine carotid arteries

A method was developed to generate and record three types of small amplitude waves (pressure, torsion and axial) in the exposed carotid artery of anesthetized dogs. The pressure waves were studied with the aid of miniature pressure transducers; electro-optical tracking units monitored the axial and circumferential surface displacements. Results from 6 dogs are presented in the form of the phase velocities and attenuation of three types of waves. The data demonstrate incompatibility with an isotropic elastic model for the mechanical behavior of the artery. The measured damping appears to be primarily due to the viscoelastic properties of the vessel wall material.

Moritz, W. E.↗

Numerical analysis of some problems related to the mechanics of pneumatic tires: Finite deformation/rolling contact of a viscoelastic cylinder and finite deformation of cord-reinforced rubber composites

The formulation and numerical analysis of several problems related to the behavior of pneumatic tires are considered. These problems include the general rolling contact problem of a rubber-like viscoelastic cylinder undergoing finite deformations and the finite deformation of cord-reinforced rubber composites. New finite element models are developed for these problems. Numerical results obtained for several representative cases are presented.

Oden, J. T.↗

Relation of ultrasonic energy loss factors and constituent properties in unidirectional composites

A model is developed relating composite constituents properties with ultrasonic energy loss factors for longitudinal waves propagating in the principal directions of a unidirectional graphite/epoxy fiber composite. All the constituents are assumed to behave as linear viscoelastic materials with energy dissipation properties defined by loss factors. It is found that by introducing a new constituent called the interface material, the composite and constituent properties can be brought into consistency with simple series and parallel models. An expression relating the composite loss factors to the loss factors of the constituents is derived and its coefficients are evaluated.

Williams, J. H., Jr.↗

Exact finite element method analysis of viscoelastic tapered structures to transient loads

A general method is presented for determining the dynamic torsional/axial response of linear structures composed of either tapered bars or shafts to transient excitations. The method consists of formulating and solving the dynamic problem in the Laplace transform domain by the finite element method and obtaining the response by a numerical inversion of the transformed solution. The derivation of the torsional and axial stiffness matrices is based on the exact solution of the transformed governing equation of motion, and it consequently leads to the exact solution of the problem. The solution permits treatment of the most practical cases of linear tapered bars and shafts, and employs modeling of structures with only one element per member which reduces the number of degrees of freedom involved. The effects of external viscous or internal viscoelastic damping are also taken into account.

Spyrakos, Constantine Chris↗

Ultrasonic wave propagation loss factor in composite in terms of constituent properties

A model is developed to relate the energy loss factors of ultrasonic longitudinal waves propagating in the principal directions of a unidirectional graphite fiber composite to the composite constituent properties. All the constituents are assumed to behave as linear viscoelastic materials with energy dissipation properties defined by loss factors. It is found that by introducing a new constituent called the interface material, the composite and constituent properties can be brought into consistency with simple series and parallel models. An expression relating the composite loss factors to the loss factors of the constituents is derived and its coefficients are evaluated.

Williams, J. H., Jr.↗

The viscoelastic behavior of notched glassy polymers

In the bulk, glassy polymers exhibit a nonlinear viscoelastic response during deformation. Stress or strain induced damage (i.e. crazing, microshear banding) results in the production of nonrecoverable work and observed nonlinearity. Stress or strain dependent shift factors have been used to mathematically model the mechanical behavior of these polymers. Glassy polymers that have been notched, may exhibit very different load displacement response compared to the same material under bulk deformation. If a sharp notch is introduced into the body then loaded, the load displacement trace may appear to be single-valued in the absence of viscoelasticity and crack growth. This suggests the volume of damaged material is small compared to the overall dimensions of the specimen. The ability to produce a single-valued load-load-line displacement trace through the use of the Correspondence Principle may prove to be useful for fracture of viscoelastic materials.

Crook, R. A.↗

Mechanical property of solithane 113

Models to study failure under stress of Solithane 113 and set of equations for stress, strain, and energy responses of linear viscoelastic systems for wide spectrum of load inputs

ENERGY EQUATION↗

Crustal deformation in great California earthquake cycles

Periodic crustal deformation associated with repeated strike slip earthquakes is computed for the following model: A depth L (less than or similiar to H) extending downward from the Earth's surface at a transform boundary between uniform elastic lithospheric plates of thickness H is locked between earthquakes. It slips an amount consistent with remote plate velocity V sub pl after each lapse of earthquake cycle time T sub cy. Lower portions of the fault zone at the boundary slip continuously so as to maintain constant resistive shear stress. The plates are coupled at their base to a Maxwellian viscoelastic asthenosphere through which steady deep seated mantle motions, compatible with plate velocity, are transmitted to the surface plates. The coupling is described approximately through a generalized Elsasser model. It is argued that the model gives a more realistic physical description of tectonic loading, including the time dependence of deep slip and crustal stress build up throughout the earthquake cycle, than do simpler kinematic models in which loading is represented as imposed uniform dislocation slip on the fault below the locked zone.

Li, Victor C.↗

The effects of physical aging at elevated temperatures on the viscoelastic creep on IM7/K3B

Physical aging at elevated temperature of the advanced composite IM7/K3B was investigated through the use of creep compliance tests. Testing consisted of short term isothermal, creep/recovery with the creep segments performed at constant load. The matrix dominated transverse tensile and in-plane shear behavior were measured at temperatures ranging from 200 to 230 C. Through the use of time based shifting procedures, the aging shift factors, shift rates and momentary master curve parameters were found at each temperature. These material parameters were used as input to a predictive methodology, which was based upon effective time theory and linear viscoelasticity combined with classical lamination theory. Long term creep compliance test data was compared to predictions to verify the method. The model was then used to predict the long term creep behavior for several general laminates.

Gates, Thomas S.↗