Search NASA⌕ Search

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.

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

At least 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

Numerical solution methods for viscoelastic orthotropic materials

Numerical solution methods for viscoelastic orthotropic materials, specifically fiber reinforced composite materials, are examined. The methods include classical lamination theory using time increments, direction solution of the Volterra Integral, Zienkiewicz's linear Prony series method, and a new method called Nonlinear Differential Equation Method (NDEM) which uses a nonlinear Prony series. The criteria used for comparison of the various methods include the stability of the solution technique, time step size stability, computer solution time length, and computer memory storage. The Volterra Integral allowed the implementation of higher order solution techniques but had difficulties solving singular and weakly singular compliance function. The Zienkiewicz solution technique, which requires the viscoelastic response to be modeled by a Prony series, works well for linear viscoelastic isotropic materials and small time steps. The new method, NDEM, uses a modified Prony series which allows nonlinear stress effects to be included and can be used with orthotropic nonlinear viscoelastic materials. The NDEM technique is shown to be accurate and stable for both linear and nonlinear conditions with minimal computer time.

Gramoll, K. C.↗
NASA NTRS2017

Non-Linear Analysis of the NASA Super Pressure Balloons: Whole Flight Simulations

Tensys have a long-established background in the shape generation and load analysis of architectural stressed membrane structures. Founded upon their inTENS finite element analysis suite, these activities have broadened to encompass 'lighter than air' structures such as aerostats, airships, hybrid air vehicles and stratospheric balloons. Tensys have acted as consultants to the NASA Super Pressure Balloon (SPB) Program since 2004. Previous papers have focussed upon the application of inTENS to the overall structural and stability analysis of pumpkin type balloons as used in the SPB Program. Particular emphasis has been placed upon the ability to study both stress and stability at all stages of a flight. As the program has developed, increasing modelling fidelity has been introduced as design refinement has moved emphasis from the overall shape to the performance of individual details. Examples include the introduction of a contact capability into inTENS to better represent the separate cap layer of film. Localised investigations have considered the consequence of local geometric anomalies introduced during fabrication and the effects of debonding of the PBO tendon within its sleeve. Analysis to date has used a material model for the polyethylene shell film developed by Dr Rand of Winzen Engineering, supported by a program of fine resolution material tests at the Balloon Research and Development Laboratory (BRDL) at the GSFC Wallops Island Facility. Based upon work first reported by Schapery, this model has been incorporated into so-called 'snapshot' analyses within inTENS. For a given elapsed time point and temperature distribution, non-linear material properties are iteratively updated for the current stress state in each individual element until changes in those material properties are insignificant. This process has to be iterative as the film stresses and material properties are interdependent. Attempts to derive a time-stepping incremental viscoelastic capability using the Schapery Rand model encountered problems when dealing with the varying temperatures associated with balloon deployment, pressurisation and diurnal behaviour. An alternative approach has been developed by Pellegrino et al at Caltech, again with support from the NASA Balloon Program Office. This is a large strain non-linear viscoelastic model that includes film out-of-plane mechanical and thermal effects. Working from a unified base, this model can be utilised for time-stepping analyses in both a modulus or compliance mode, or a combination of both. This new capability is intended to enhance the current NASA balloon design process. This paper presents the implementation of the Caltech model into a specialist finite element suite, inTENS and its application to whole flight simulations.

Wakefield, David↗
DOE OSTI2023

Anticipating gelation and vitrification with medium amplitude parallel superposition (MAPS) rheology and artificial neural networks

Abstract Anticipating qualitative changes in the rheological response of complex fluids (e.g., a gelation or vitrification transition) is an important capability for processing operations that utilize such materials in real-world environments. One class of complex fluids that exhibits distinct rheological states are soft glassy materials such as colloidal gels and clay dispersions, which can be well characterized by the soft glassy rheology (SGR) model. We first solve the model equations for the time-dependent, weakly nonlinear response of the SGR model. With this analytical solution, we show that the weak nonlinearities measured via medium amplitude parallel superposition (MAPS) rheology can be used to anticipate the rheological aging transitions in the linear response of soft glassy materials. This is a rheological version of a technique called structural health monitoring used widely in civil and aerospace engineering. We design and train artificial neural networks (ANNs) that are capable of quickly inferring the parameters of the SGR model from the results of sequential MAPS experiments. The combination of these data-rich experiments and machine learning tools to provide a surrogate for computationally expensive viscoelastic constitutive equations allows for rapid experimental characterization of the rheological state of soft glassy materials. We apply this technique to an aging dispersion of Laponite ® clay particles approaching the gel point and demonstrate that a trained ANN can provide real-time detection of transitions in the nonlinear response well in advance of incipient changes in the linear viscoelastic response of the system.

Lennon, Kyle R. (ORCID:0000000212515461)↗
NASA NTRS1979

Mechanics of fracture - Fundamentals and some recent developments

An overview is presented of the fundamental aspects of and recent developments in fracture mechanics. Reference is made to linear elastic fracture mechanics including the state of stresses and displacements in the vicinity of cracks, effects of crack geometry and orientation on stress intensity factors, energy balance of Griffith, Irwin's stress intensity concept, and linear elastic fracture mechanics testing for fracture toughness. Other aspects of this paper include the non-linear behavior of materials and their influence on fracture mechanics parameters, consideration of viscoelasticity and plasticity, non-linear fracture toughness parameters as C.O.D., R-curve and J-integral, and a non-linear energy method, proposed by Liebowitz. Finite element methods applied to fracture mechanics problems are indicated. Also, consideration has been given to slow crack growth, dynamic effects on K(IC), Sih's criterion for fracture, Lee and Liebowitz's criterion relating crack growth with plastic energy, and applications of fracture mechanics to aircraft design. Suggestions are offered for future research efforts to be undertaken in fracture mechanics.

Liebowitz, H.↗
NASA NTRS1975

Rate and time dependent failure of structural adhesives

Studies on two adhesives (Metlbond 1113 and 1113-2) identified as having important applications in the bonding of composite materials are presented. A testing program to ascertain stress-strain, strain-rate, time, yield, and/or failure behavior of these materials in bulk form using uniaxial tensile constant strain-rate, creep, and relaxation tests is described. The stress-strain behavior of each material is shown to be significantly rate dependent. A rate dependent stress whitening (crazing) phenomenon occurs prior to either yield or fracture. A region of linear elasticity, a region of viscoelasticity, and the onset of yielding are identified in the stress-strain behavior. The linear elastic limit and the yield point are shown to be rate dependent and agree well with an empirical equation proposed by Ludwik. A creep to failure phenomenon is shown to exist and is correlated with a delayed yield equation proposed by Crochet. Analytical predictions based on a modified Bingham model are shown to agree well with experimental stress-strain strain-rate data. Analytical predictions based on a modified Ramberg-Osgood equation are also shown for comparison purposes.

Brinson, H. F.↗
NASA NTRS1995

Expendable launch vehicle studies

Analytical support studies of expendable launch vehicles concentrate on the stability of the dynamics during launch especially during or near the region of maximum dynamic pressure. The in-plane dynamic equations of a generic launch vehicle with multiple flexible bending and fuel sloshing modes are developed and linearized. The information from LeRC about the grids, masses, and modes is incorporated into the model. The eigenvalues of the plant are analyzed for several modeling factors: utilizing diagonal mass matrix, uniform beam assumption, inclusion of aerodynamics, and the interaction between the aerodynamics and the flexible bending motion. Preliminary PID, LQR, and LQG control designs with sensor and actuator dynamics for this system and simulations are also conducted. The initial analysis for comparison of PD (proportional-derivative) and full state feedback LQR Linear quadratic regulator) shows that the split weighted LQR controller has better performance than that of the PD. In order to meet both the performance and robustness requirements, the H(sub infinity) robust controller for the expendable launch vehicle is developed. The simulation indicates that both the performance and robustness of the H(sub infinity) controller are better than that for the PID and LQG controllers. The modelling and analysis support studies team has continued development of methodology, using eigensensitivity analysis, to solve three classes of discrete eigenvalue equations. In the first class, the matrix elements are non-linear functions of the eigenvector. All non-linear periodic motion can be cast in this form. Here the eigenvector is comprised of the coefficients of complete basis functions spanning the response space and the eigenvalue is the frequency. The second class of eigenvalue problems studied is the quadratic eigenvalue problem. Solutions for linear viscously damped structures or viscoelastic structures can be reduced to this form. Particular attention is paid to Maxwell and Kelvin models. The third class of problems consists of linear eigenvalue problems in which the elements of the mass and stiffness matrices are stochastic. dynamic structural response for which the parameters are given by probabilistic distribution functions, rather than deterministic values, can be cast in this form. Solutions for several problems in each class will be presented.

Bainum, Peter M.↗
NASA NTRS1984

Secular rotational motions and the mechanical structure of a dynamical viscoelastic earth

A survey is presented of analytical methods for computing the linear responses of the rotational axis of a layered viscoelastic earth to surface loading. Theoretical research in this area is first summarized, and the differences between the mechanical boundary conditions to be applied at the interface separating the upper and lower mantles for an adiabatically and chemically stratified mantle are discussed. Some examples of polar wander and secular variation of the spin rate from glacial excitation are presented for various types of chemical and viscosity stratifications. The effects of an artificial density jump at the base of the lithosphere in models are examined, and certain issues concerning the fluid tidal Love number for different types of density stratification are addressed. The meaning of effective plate thickness over geological time scales for rotational dynamics is discussed.

Yuen, D. A.↗
DOE OSTI2023

Peri-Net-Pro: the neural processes with quantified uncertainty for crack patterns

Abstract This paper develops a deep learning tool based on neural processes (NPs) called the Peri-Net-Pro, to predict the crack patterns in a moving disk and classifies them according to the classification modes with quantified uncertainties. In particular, image classification and regression studies are conducted by means of convolutional neural networks (CNNs) and NPs. First, the amount and quality of the data are enhanced by using peridynamics to theoretically compensate for the problems of the finite element method (FEM) in generating crack pattern images. Second, case studies are conducted with the prototype microelastic brittle (PMB), linear peridynamic solid (LPS), and viscoelastic solid (VES) models obtained by using the peridynamic theory. The case studies are performed to classify the images by using CNNs and determine the suitability of the PMB, LBS, and VES models. Finally, a regression analysis is performed on the crack pattern images with NPs to predict the crack patterns. The regression analysis results confirm that the variance decreases when the number of epochs increases by using the NPs. The training results gradually improve, and the variance ranges decrease to less than 0.035. The main finding of this study is that the NPs enable accurate predictions, even with missing or insufficient training data. The results demonstrate that if the context points are set to the 10th, 100th, 300th, and 784th, the training information is deliberately omitted for the context points of the 10th, 100th, and 300th, and the predictions are different when the context points are significantly lower. However, the comparison of the results of the 100th and 784th context points shows that the predicted results are similar because of the Gaussian processes in the NPs. Therefore, if the NPs are employed for training, the missing information of the training data can be supplemented to predict the results.

Mathematics↗
NASA NTRS1972

Governing equations for vibrating constrained-layer damping sandwich plates and beams.

A simple differential equation is derived to describe constrained-layer damping in nonsymmetric sandwich plates and beams composed of isotropic and homogeneous layers. The natural boundary conditions related to this equation are determined and some typical numerical results obtained by this equation are given. The equation is valid within the linear theories of elasticity and viscoelasticity in the absence of any constraints on thicknesses, positions, symmetries, and densities of the layers.

Yan, M.-J.↗
NASA NTRS1985

Bounding solutions of geometrically nonlinear viscoelastic problems

Integral transform techniques, such as the Laplace transform, provide simple and direct methods for solving viscoelastic problems formulated within a context of linear material response and using linear measures for deformation. Application of the transform operator reduces the governing linear integro-differential equations to a set of algebraic relations between the transforms of the unknown functions, the viscoelastic operators, and the initial and boundary conditions. Inversion either directly or through the use of the appropriate convolution theorem, provides the time domain response once the unknown functions have been expressed in terms of sums, products or ratios of known transforms. When exact inversion is not possible approximate techniques may provide accurate results. The overall problem becomes substantially more complex when nonlinear effects must be included. Situations where a linear material constitutive law can still be productively employed but where the magnitude of the resulting time dependent deformations warrants the use of a nonlinear kinematic analysis are considered. The governing equations will be nonlinear integro-differential equations for this class of problems. Thus traditional as well as approximate techniques, such as cited above, cannot be employed since the transform of a nonlinear function is not explicitly expressible.

Stubstad, J. M.↗
NASA NTRS1986

Bounding solutions of geometrically nonlinear viscoelastic problems

Integral transform techniques, such as the Laplace transform, provide simple and direct methods for solving viscoelastic problems formulated within a context of linear material response and using linear measures for deformation. Application of the transform operator reduces the governing linear integro-differential equations to a set of algebraic relations between the transforms of the unknown functions, the viscoelastic operators, and the initial and boundary conditions. Inversion either directly or through the use of the appropriate convolution theorem, provides the time domain response once the unknown functions have been expressed in terms of sums, products or ratios of known transforms. When exact inversion is not possible approximate techniques may provide accurate results. The overall problem becomes substantially more complex when nonlinear effects must be included. Situations where a linear material constitutive law can still be productively employed but where the magnitude of the resulting time dependent deformations warrants the use of a nonlinear kinematic analysis are considered. The governing equations will be nonlinear integro-differential equations for this class of problems. Thus traditional as well as approximate techniques, such as cited above, cannot be employed since the transform of a nonlinear function is not explicitly expressible.

Stubstad, J. M.↗
NASA NTRS1982

Postseismic deformation due to subcrustal viscoelastic relaxation following dip-slip earthquakes

The deformation of the Earth following a dip-slip earthquake is calculated using a three layer rheological model and finite element techniques. The three layers are an elastic upper lithosphere, a standard linear solid lower lithosphere, and a Maxwell viscoelastic asthenosphere-a model previously analyzed in the strike-clip case (Cohen, 1981, 1982). Attention is focused on the magnitude of the postseismic subsidence and the width of the subsidence zone that can develop due to the viscoelastic response to coseismic reverse slip. Detailed analysis for a fault extending from the surface to 15 km with a 45 deg dip reveals that postseismic subsidence is sensitive to the depth to the asthenosphere but is only weakly dependent on lower lithosphere depth. The greatest subsidence occurs when the elastic lithosphere is about 30 km thick and the asthenosphere lies just below this layer (asthenosphere depth = 2 times the fault depth). The extremum in the subsidence pattern occurs at about 5 km from the surface trace of the fault and lies over the slip plane. In a typical case after a time t = 30 tau (tau = Maxwell time) following the earthquake the subsidence at this point is 60% of the coseismic uplift. Unlike the horizontal deformation following a strike slip earthquake, significant vertical deformation due to asthenosphere flow persists for many times tau and the magnitude of the vertical deformation is not necessarily enhanced by having a partially relaxing lower lithosphere.

Cohen, S. C.↗
NASA NTRS1984

Postseismic deformation due to subcrustal viscoelastic relaxation following dip-slip earthquakes

The deformation of the earth following a dip-slip earthquake is calculated using a three layer rheological model and finite element techniques. The three layers are an elastic upper lithosphere, a standard linear solid lower lithosphere, and a Maxwell viscoelastic asthenosphere - a model previously analyzed in the strike-clip case (Cohen, 1981, 1982). Attention is focused on the magnitude of the postseismic subsidence and the width of the subsidence zone that can develop due to the viscoelastic response to coseismic reverse slip. Detailed analysis for a fault extending from the surface to 15 km with a 45 deg dip reveals that postseismic subsidence is sensitive to the depth to the asthenosphere but is only weakly dependent on lower lithosphere depth. The greatest subsidence occurs when the elastic lithosphere is about 30 km thick and the asthenosphere lies just below this layer (asthenosphere depth = 2 times the fault depth). The extremum in the subsidence pattern occurs at about 5 km from the surface trace of the fault and lies over the slip plane. In a typical case after a time t = 30 tau (tau = Maxwell time) following the earthquake, the subsidence at this point is 60 percent of the coseismic uplift. Unlike the horizontal deformation following a strike slip earthquake, significant vertical deformation due to asthenosphere flow persists for many times tau and the magnitude of the vertical deformation is not necessarily enhanced by having a partially relaxing lower lithosphere. Previously announced in STAR as N83-13683

Cohen, S. C.↗
NASA NTRS2005

Biaxial Stress Limit for ULDB Film

The current ULDB design applies stress to the shell film biaxially to control creep in the latitudinal direction. The recent change in design paradigm, from a uniaxial to biaxial stress state basis, arose from a new perspective that the biaxial loading can control strain in both principal surface dimensions as discussed below. The current ULDB project path was thus enabled by a more thorough understanding of the nonlinear viscoelastic properties of the shell film material, linear low-density polyethylene (LLDPE). Although a very similar material is also used in NASA zero-pressure (ZPB) and long-duration balloons (LDB), the different stress state requires a new approach to shell material qualification.

Sterling, W. J.↗
DOE OSTI2021

Learning viscoelasticity models from indirect data using deep neural networks

In this study, we propose a novel approach to model viscoelasticity materials, where rate-dependent and non-linear constitutive relationships are approximated with deep neural networks. We assume that inputs and outputs of the neural networks are not directly observable, and therefore common training techniques with input–output pairs for the neural networks are inapplicable. To that end, we develop a novel computational approach to both calibrate parametric and learn neural-network-based constitutive relations of viscoelasticity materials from indirect displacement data in the context of multiple-physics systems. We show that limited displacement data holds sufficient information to quantify the viscoelasticity behavior. We formulate the inverse computation – modeling viscoelasticity properties from observed displacement data – as a PDE-constrained optimization problem and minimize the error functional using a gradient-based optimization method. The gradients are computed by a combination of automatic differentiation and implicit function differentiation rules. The effectiveness of our method is demonstrated through numerous benchmark problems in geomechanics and porous media transport.

97 MATHEMATICS AND COMPUTING↗