SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Technique computes stresses resulting from axisymmetric transient thermal loading in circular solid-propellant grain section with circular ports. Propellant is assumed to be linear, thermal rheologically simple, viscoelastic material; material properties are represented by exponential series in time.
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.
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.
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.
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.
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.
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.
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.
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.
Viscoelastic theories evaluated for ability to correlate linear dynamic data with nonlinear viscosity and normal stress data
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.
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
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.
Lagrangian turbulence in continua is introduced as a failure of Euclidian metric of a material system of coordinates, and it is shown that the Ricci tensor of such a system of coordinates may be used as a measure of L-turbulence. The conditions which are sufficient, though not necessary, for instability of an initially infinitesimal disturbance of this tensor are determined for fluid, elastic, and viscoelastic bodies. The results are obtained from the linearized equations, and therefore the growth of unstable solutions may be bounded.
Anticipating the implementation of advanced SiC/SiC composites into turbine section components of future aero-propulsion engines, the primary objective of this on-going study is to develop physics-based analytical and finite-element modeling tools to predict the effects of constituent creep on SiC/SiC component service life. A second objective is to understand how to possibly manipulate constituent materials and processes in order to minimize these effects. Focusing on SiC/SiC components experiencing through-thickness stress gradients (e.g., airfoil leading edge), prior NASA creep modeling studies showed that detrimental residual stress effects can develop globally within the component walls which can increase the risk of matrix cracking. These studies assumed that the SiC/SiC composites behaved as isotropic viscoelastic continuum materials with creep behavior that was linear and symmetric with stress and that the creep parameters could be obtained from creep data as experimentally measured in-plane in the fiber direction of advanced thin-walled 2D SiC/SiC panels. The present study expands on those prior efforts by including constituent behavior with non-linear stress dependencies in order to predict such key creep-related SiC/SiC properties as time-dependent matrix stress, constituent creep and content effects on composite creep rates and rupture times, and stresses on fiber and matrix during and after creep.
A linearized stability analysis has been applied to a fluid flowing in a gravity field between horizontal planes in Couette flow under conditions such that the temperature of the bottom plane exceeds that of the top. It is shown that, under conditions likely to be encountered with polymer solutions, oscillatory instabilities will not be controlling. Criteria are offered for ascertaining when an analysis based upon a second-order fluid model may be expected to yield physically meaningful results. It is also shown that, for the fluid model considered, critical conditions for stability are not changed when disturbances which vary in the flow direction are substituted for those which are a function of the coordinate transverse to the flow.
Well-posed models and computational algorithms are developed and analyzed for control of a class of partial differential equations that describe the motions of thermo-viscoelastic structures. An abstract (state space) framework and a general well-posedness result are presented that can be applied to a large class of thermo-elastic and thermo-viscoelastic models. This state space framework is used in the development of a computational scheme to be used in the solution of a linear quadratic regulator (LQR) control problem. A detailed convergence proof is provided for the viscoelastic model and several numerical results are presented to illustrate the theory and to analyze problems for which the theory is incomplete.
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.