SEARCH · Engineering Papers

Results for “stiff problem”

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 595 records · Page 33

All sources NASA NTRS DOE OSTI

NASA NTRS1986

Inclusion of elastically connected members in exact buckling and frequency calculations

A standard stiffness matrix procedure which permits any combination of rigid, elastic, pinned or sliding connections of the degrees of freedom at the ends of a member to the nodes of its parent structure is described, in order to show how easily it can be extended to allow an existing algorithm to be used to ensure that no eigenvalues of the parent structure can be missed even when 'exact' member theory is used. The eigenvalues are the natural frequencies of undamped free vibration analyses or the critical load factors of buckling problems. The method preserves the exactness of the member theory and an efficient method for computer application is indicated. The theory also permits any combination of rigid, elastic, pinned or sliding connections between the freedoms of a substructure and those of its parent structure.

Williams, F. W.↗

DOE OSTI2022

Hydromechanical impact of basement rock on injection-induced seismicity in Illinois Basin

Abstract The common explanation of observed injection-induced microseismicity is based on the measured stress state at the injection interval and the assumption that it remains the same in the vicinity. We argue here that representing the stress state in different geologic formations over the injection site with the single Mohr’s circle is insufficient due to local stratigraphic features and contrast in compressibilities of the involved formations. The role of hydromechanical coupling in the microseismic response is also crucial for the proper assessment of the problem. Thoroughly monitored Illinois Basin Decatur Project revealed the majority of CO 2 injection-associated microseismic events being originated in the crystalline basement. Even though basement faults can serve as the conduits for fluid flow—the predicted pressure increase seems to be insufficient to trigger seismicity. To address this issue, accurate laboratory measurements of rock properties from the involved formations are conducted. The pre-injection stress state and its evolution are evaluated with the hydromechanically coupled numerical model. It appears that the presence of an offset in a stiff competent layer affects the stress state in its vicinity. Therefore, both the pre-injection stress state and its evolution during the fluid injection should be addressed during the induced seismicity assessment.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1987

The effects of pin elasticity, clearance, and friction on the stresses in a pin-loaded orthotropic plate

The effects of pin elasticity, clearance, and friction on the stresses in a pin loaded orthotropic plate are studied. The effects are studied by posing the problem as a planar contact elasticity problem, the pin and the plate being two elastic bodies which interact through contact. Coulomb friction is assumed, the pin loads the plate in one of its principal material directions, and the plate is infinite in extent. A collocation scheme and interaction, in conjunction with a complex variable series solution, are used to obtain numerical results. The contact region between the plate and pin is unknown and must be solved for as part of the solution. The same is true of the region of friction induced no slip. Two pin stiffnesses, two clearance levels, two friction levels and two laminates, a (0/+ or - 45/90)s and a (02/+ or - 45)s, are studied. The effects of pin elasticity, clearance, and friction on the load capacity of the plate are assessed by comparing the load capacity of the plate with the capacity when the pin is rigid, perfectly fitting, and frictionless.

Hyer, M. W.↗

NASA NTRS1972

Thermochemical nonequilibrium in atomic hydrogen at elevated temperatures

A numerical study of the nonequilibrium flow of atomic hydrogen in a cascade arc was performed to obtain insight into the physics of the hydrogen cascade arc. A rigorous mathematical model of the flow problem was formulated, incorporating the important nonequilibrium transport phenomena and atomic processes which occur in atomic hydrogen. Realistic boundary conditions, including consideration of the wall electrostatic sheath phenomenon, were included in the model. The governing equations of the asymptotic region of the cascade arc were obtained by writing conservation of mass and energy equations for the electron subgas, an energy conservation equation for heavy particles and an equation of state. Finite-difference operators for variable grid spacing were applied to the governing equations and the resulting system of strongly coupled, stiff equations were solved numerically by the Newton-Raphson method.

Scott, R. K.↗

NASA NTRS1987

Effects of transverse shearing flexibility on the postbuckling of plates loaded by inplane shear

This paper presents buckling and postbuckling results for plates loaded by inplane shear. The buckling results have been plotted to show the effects of thickness on the stress coefficient for aluminum plates. Results are given for various length-to-width ratios. Postbuckling results for thin plates with transverse shearing flexibility are compared to results from classical theory. The problems considered are the postbuckling response of plates in shear made of aluminum and of a + or - 45 deg graphite-epoxy laminate. Thus the materials are isotropic and orthotropic, respectively. The plates are considered to be long with side edges simply supported, with various inplane edge conditions, and the plates are subject to a constant shearing displacement along the side edges. Characteristic curves presenting the average shear stress resultant as a function of the applied displacement are given. These curves indicate that change in inplane edge conditions influence plate postbuckling stiffness and that transverse shearing is important in some cases.

Stein, Manuel↗

NASA NTRS1997

Dynamic Reconstruction and Multivariable Control for Force-Actuated, Thin Facesheet Adaptive Optics

The Multiple Mirror Telescope (MMT) under development at the University of Arizona takes a new approach in adaptive optics placing a large (0.65 m) force-actuated, thin facesheet deformable mirror at the secondary of an astronomical telescope, thus reducing the effects of emissivity which are important in IR astronomy. However, The large size of the mirror and low stiffness actuators used drive the natural frequencies of the mirror down into the bandwidth of the atmospheric distortion. Conventional adaptive optics takes a quasi-static approach to controlling the, deformable mirror. However, flexibility within the control bandwidth calls for a new approach to adaptive optics. Dynamic influence functions are used to characterize the influence of each actuator on the surface of the deformable mirror. A linearized model of atmospheric distortion is combined with dynamic influence functions to produce a dynamic reconstructor. This dynamic reconstructor is recognized as an optimal control problem. Solving the optimal control problem for a system with hundreds of actuators and sensors is formidable. Exploiting the circularly symmetric geometry of the mirror, and a suitable model of atmospheric distortion, the control problem is divided into a number of smaller decoupled control problems using circulant matrix theory. A hierarchic control scheme which seeks to emulate the quasi-static control approach that is generally used in adaptive optics is compared to the proposed dynamic reconstruction technique. Although dynamic reconstruction requires somewhat more computational power to implement, it achieves better performance with less power usage, and is less sensitive than the hierarchic technique.

Grocott, Simon C. O.↗

DOE OSTI2023

Analysis of Neural Networks as Random Dynamical Systems

In this report we present our findings and outcomes of the NNRDS (analysis of Neural Networks as Random Dynamical Systems) project. The work is largely motivated by the analogy of a large class of neural networks (NNs) with a discretized ordinary differential equation (ODE) schemes. Namely, residual NNs, or ResNets, can be viewed as a discretization of neural ODEs (NODEs) where the NN depth plays the role of the time evolution. We employ several legacy tools from ODE theory, such as stiffness, nonlocality, autonomicity, to enable regularization of ResNets thus improving their generalization capabilities. Furthermore, armed with NN analysis tools borrowed from the ODE theory, we are able to efficiently augment NN predictions with uncertainty overcoming wellknown dimensionality challenges and adding a degree of trust towards NN predictions. Finally, we have developed a Python library QUiNN (Quantification of Uncertainties in Neural Networks) that incorporates improved-architecture ResNets, besides classical feed-forward NNs, and contains wrappers to PyTorch NN models enabling several major classes of uncertainty quantification methods for NNs. Besides synthetic problems, we demonstrate the methods on datasets from climate modeling and materials science.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1976

On the problem of stress singularities in bonded orthotropic materials

The problem of stress singularities at the leading edge of a crack lying in the neighborhood of a bimaterial interface in bonded orthotropic materials is considered. The main objective is to study the effect of material orthotropy on the singular behavior of the stress state when the crack touches or intersects the interface. The results indicate that, due to the large number of material constants involved, in orthotropic materials, the power of stress singularity as well as the stress intensity factor can be considerably different than that found in the isotropic materials with the same stiffness ratio perpendicular to the crack.

Erdogan, F.↗

NASA NTRS2007

NASA Langley's Approach to the Sandia's Structural Dynamics Challenge Problem

The objective of this challenge is to develop a data-based probabilistic model of uncertainty to predict the behavior of subsystems (payloads) by themselves and while coupled to a primary (target) system. Although this type of analysis is routinely performed and representative of issues faced in real-world system design and integration, there are still several key technical challenges that must be addressed when analyzing uncertain interconnected systems. For example, one key technical challenge is related to the fact that there is limited data on target configurations. Moreover, it is typical to have multiple data sets from experiments conducted at the subsystem level, but often samples sizes are not sufficient to compute high confidence statistics. In this challenge problem additional constraints are placed as ground rules for the participants. One such rule is that mathematical models of the subsystem are limited to linear approximations of the nonlinear physics of the problem at hand. Also, participants are constrained to use these models and the multiple data sets to make predictions about the target system response under completely different input conditions. Our approach involved initially the screening of several different methods. Three of the ones considered are presented herein. The first one is based on the transformation of the modal data to an orthogonal space where the mean and covariance of the data are matched by the model. The other two approaches worked solutions in physical space where the uncertain parameter set is made of masses, stiffnesses and damping coefficients; one matches confidence intervals of low order moments of the statistics via optimization while the second one uses a Kernel density estimation approach. The paper will touch on all the approaches, lessons learned, validation 1 metrics and their comparison, data quantity restriction, and assumptions/limitations of each approach. Keywords: Probabilistic modeling, model validation, uncertainty quantification, kernel density

Horta, Lucas G.↗

NASA NTRS1989

Structural optimization of rotor blades with straight and swept tips subject to aeroelastic constraints

The main conclusions obtained in the present study are summarized. Their application to the structural optimization of a helicopter blade should be limited by the assumptions used in obtaining the numerical results presented here. The optimum design procedure described here is very efficient, and can produce improved designs with a very limited number of precise analyses. The method of constructing the approximate problem is such that previously conducted aeroelastic analyses can be reused in a new optimization problem. For example, if an optimization study is preceded by a parametric study in which the effect of various combinations of blade design parameters is examined, all the aeroelastic analyses performed for the parametric study can be reutilized in the optimization study. This is not possible when the approximate problem is built from Taylor series expansions. The results of the optimization are quite sensitive to the aeroelastic stability margins required of the blade. In the optimization of case 2, changing the aeroelastic stability constraints from simply requiring that the blade be stable in hover, to requiring that the stability margins be maintained during the course of the optimization, reduced the gains in n/rev vibration levels by more than 50 percent. The introduction of tip sweep can reduce the n/rev vertical hub shears beyond the level that can be obtained by just modifying the mass and stiffness distributions of the blade.

Friedmann, Peretz P.↗

DOE OSTI2021

Dissipative magnetohydrodynamics for nonresistive relativistic plasmas: An implicit second-order flux-conservative formulation with stiff relaxation

Based on a 14-moment closure for nonresistive (general-) relativistic viscous plasmas, here we describe a new numerical scheme that is able to handle all first-order dissipative effects (heat conduction, bulk and shear viscosities), as well the anisotropies induced by the presence of magnetic fields. The latter is parametrized in terms of a thermal gyrofrequency or, equivalently, a thermal Larmor radius and allows to correctly capture the thermal Hall effect. By solving an extended Israel-Stewart-like system for the dissipative quantities that enforces algebraic constraints via stiff-relaxation, we are able to cast all first-order dissipative terms in flux-divergence form. This allows us to apply traditional high-resolution shock capturing methods to the equations, making the system suitable for the numerical study of highly turbulent flows. We present several numerical tests to assess the robustness of our numerical scheme in flat spacetime. The 14-moment closure can seamlessly interpolate between the highly collisional limit found in neutron star mergers, and the highly anisotropic limit of relativistic Braginskii magnetohydrodynamics appropriate for weakly collisional plasmas in black-hole accretion problems. We believe that this new formulation and numerical scheme will be useful for a broad class of relativistic magnetized flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1975

Sensitivity of Hingeless Rotor Blade Flap-lag Stability in Hover to Analytical Modelling Assumptions

Prediction of flap-lag stability using a single bending mode for each degree-of-freedom is examined in the case in which the bending modes are assumed to be the same in the flap and lag directions and are independent of pitch angle and stiffness distribution. It is shown that this model gives results analogous to those obtained by Ormiston employing a rigid blade model with the blade and hub stiffness represented by springs in the limiting cases of the elastic coupling parameter R = 0 and 1. For intermediate values of R the results are shown to be quite different. The mode shape assumptions are shown to result in what is referred to as the parallel spring model in contrast to Ormiston's model which is referred to as a series spring model. The similarities and differences between these two models are developed in some details. The differences between these two models are examined for various typical rotor blade characteristics. Other aspects of the sensitivity of this problem are also considered.

Curtiss, H. C., Jr.↗

NASA NTRS1986

Analysis of Delamination Growth from Matrix Cracks in Laminates Subjected to Bending Loads

A major source of delamination damage in laminated composite materials is from low-velocity impact. In thin composite laminates under point loads, matrix cracks develop first in the plies, and delaminations then grow from these cracks at the ply interfaces. The purpose of this study was to quantify the combined effects of bending and transverse shear loads on delamination initiation from matrix cracks. Graphite-epoxy laminates with 90 deg. plies on the outside were used to provide a two-dimensional simulation of the damage due to low-velocity impact. Three plate bending problems were considered: a 4-point bending, 3-point bending, and an end-clamped center-loaded plate. Under bending, a matrix crack will form on the tension side of the laminate, through the outer 90 deg. plies and parallel to the fibers. Delaminations will then grow in the interface between the cracked 90 deg. ply and the next adjacent ply. Laminate plate theory was used to derive simple equations relating the total strain energy release rate, G, associated with the delamination growth from a 90 deg. ply crack to the applied bending load and laminate stiffness properties. Three different lay-ups were tested and results compared. Test results verified that the delamination always formed at the interface between the cracked 90 deg. ply and the next adjacent ply. Calculated values for total G sub c from the analysis showed good agreement for all configurations. The analysis was able to predict the delamination onset load for the cases considered. The result indicated that the opening mode component (Mode I) for delamination growth from a matrix crack may be much larger than the component due to interlaminar shear (Mode II).

Murri, G. B.↗

NASA NTRS1985

Parametric study of stress-intensity factors in bonded composite stringer panels

Stress-intensity factors are determined for an infinite cracked orthotropic sheet adhesively bonded to an orthotropic stringer. Since the stringer is modeled as a semi-infinite sheet, the solution is most appropriate for a crack tip located near a stringer edge. Both adherends are treated as homogeneous, orthotropic media which are representative of many fiber-reinforced composite materials. The complex variable theory of elasticity was used to obtain a set of integral equations describing the problem. The integral equations are replaced by an equivalent sheet of algebraic equations, which are solved to obtain the shear stress distribution in the adhesive layer. From these adhesive stresses, the stress-intensity factors are found. A parametric study is conducted to determine the sensitivity of the system to material properties and specimen configuration. Unless the crack tip is very close to or under the stringer, the stress-intensity factor is approximately that of the unstiffened sheet. However, as the crack propagates beneath the stringer, the stress-intensity factor decreases significantly. Increasing the stringer stiffness or the adhesive stiffness also decreases the stress-intensity factor.

Bigelow, C. A.↗

NASA NTRS1988

Analysis of delamination growth from matrix cracks in laminates subjected to bending loads

A major source of delamination damage in laminated composite materials is from low-velocity impact. In thin composite laminates under point loads, matrix cracks develop first in the plies, and delaminations then grow from these cracks at the ply interfaces. The purpose of this study was to quantify the combined effects of bending and transverse shear loads on delamination initiation from matrix cracks. Graphite-epoxy laminates with 90 deg plies on the outside were used to provide a two-dimensional simulation of the damage due to low-velocity impact. Three plate bending problems were considered: a 4-point bending, 3-point bending, and an end-clamped center-loaded plate. Under bending, a matrix crack will form on the tension side of the laminate, through the outer 90 deg plies and parallel to the fibers. Delaminations will then grow in the interface between the cracked 90 deg ply and the next adjacent ply. Laminate plate theory was used to derive simple equations relating the total strain energy release rate, G, associated with the delamination growth from a 90 deg ply crack to the applied bending load and laminate stiffness properties. Three different lay-ups were tested and results compared. Test results verified that the delamination always formed at the interface between the cracked 90 deg ply and the next adjacent ply. Calculated values for total G sub c from the analysis showed good agreement for all configurations. The analysis was able to predict the delamination onset load for the cases considered. The result indicated that the opening mode component (Mode I) for delamination growth from a matrix crack may be much larger than the component due to interlaminar shear (Mode II).

Murri, Gretchen Bostaph↗

NASA NTRSDate not supplied

Aeroelastic Analysis of Highly Flexible Wings with Linearized Frequency-Domain Aerodynamics

Flutter analysis of configurations with geometric structural nonlinearities typically is done with time-domain analysis. The results from this process are computationally expensive and can yield cumbersome results that may be difficult to manage/interpret. Compared to time-domain methods, frequency-domain flutter analysis can provide additional insight into the characteristics of a flutter stability problem. By linearizing the aeroelastic problem about the nonlinear equilibrium state, this work applies frequency-domain aeroelastic analysis to the Pazy wing, the subject of the Large Deformation Working Group in the Aeroelastic Prediction Workshop. Generalized aerodynamic forces (GAFs) are computed with both a doublet-lattice method and a computational fluid dynamics solver at a range of reduced frequencies as well as a range of dynamic pressures to account for the dependence of the mode shapes on the nonlinear equilibrium state. These GAFs are used in a p-k flutter solver, which is modified to handle the dependence of the stiffness matrix and GAFs on the dynamic pressure.

Bret K Stanford↗

NASA NTRS1992

Optimal matrix approximants in structural identification

Problems of model correlation and system identification are central in the design, analysis, and control of large space structures. Of the numerous methods that have been proposed, many are based on finding minimal adjustments to a model matrix sufficient to introduce some desirable quality into that matrix. In this work, several of these methods are reviewed, placed in a modern framework, and linked to other previously known ideas in computational linear algebra and optimization. This new framework provides a point of departure for a number of new methods which are introduced here. Significant among these is a method for stiffness matrix adjustment which preserves the sparsity pattern of an original matrix, requires comparatively modest computational resources, and allows robust handling of noisy modal data. Numerical examples are included to illustrate the methods presented herein.

Beattie, C. A.↗

NASA NTRS2023

Prediction of Stiffness and Fatigue Lives of Polymer Matrix Composite Laminates Using Artificial Neural Networks

Machine learning (ML) models are increasingly being used in many engineering fields due to the advancements in ML algorithms and availability of high-speed computing power. One of the most popular ML class of models is artificial neural networks (ANN). ML is increasingly being used in the design and analysis of composite materials and structures, specifically in the constitutive modeling of composite materials with the focus on greatly accelerating multiscale analyses of composite materials and structures through development of surrogate models. Towards that end, both Python and MATLAB-based neural nets have been developed to predict initial stiffness and fatigue life of an eight-ply symmetric polymer matrix composite laminate. Two types of neural networks, a Multilayer Perceptron (MLP) and a Recurrent Neural Network (RNN), have been developed for both platforms. Results show that the both neural net types can provide an excellent estimate of initial stiffness as well as fatigue life of eight-ply symmetric polymer matrix composite laminate. RNNs are better able to capture the shape of the fatigue curve of a laminate. This tool can be very useful for system level studies to obtain an estimate of desired properties and life of PMC composite laminates. The associated surrogate models could also be used in composite multiscale analyses to replace the actual physics-based calculations at lower scales and thereby significantly increase the computational efficiency of such analyses and thus make multiscale analyses a viable industrial tool for large scale structural problems.

Composite↗