Engineering Papers⌕ Search

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.

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

At least 541 records · Page 30

An approach for relating the results of quantitative nondestructive evaluation to intrinsic properties of high-performance materials

One of the most difficult problems the manufacturing community has faced during recent years has been to accurately assess the physical state of anisotropic high-performance materials by nondestructive means. In order to advance the design of ultrasonic nondestructive testing systems, a more fundamental understanding of how ultrasonic waves travel and interact within the anisotropic material is needed. The relationship between the ultrasonic and engineering parameters needs to be explored to understand their mutual dependence. One common denominator is provided by the elastic constants. The preparation of specific graphite/epoxy samples to be used in the experimental investigation of the anisotropic properties (through the measurement of the elastic stiffness constants) is discussed. Accurate measurements of these constants will depend upon knowledge of refraction effects as well as the direction of group velocity propagation. The continuing effort for the development of improved visualization techniques for physical parameters is discussed. Group velocity images are presented and discussed. In order to fully understand the relationship between the ultrasonic and the common engineering parameters, the physical interpretation of the linear elastic coefficients (the quantities that relate applied stresses to resulting strains) are discussed. This discussion builds a more intuitional understanding of how the ultrasonic parameters are related to the traditional engineering parameters.

Miller, James G.↗

Automated Finite Element Analysis of Elastically-Tailored Plates

A procedure for analyzing and designing elastically tailored composite laminates using the STAGS finite element solver has been presented. The methodology used to produce the elastic tailoring, namely computer-controlled steering of unidirectionally reinforced composite material tows, has been reduced to a handful of design parameters along with a selection of construction methods. The generality of the tow-steered ply definition provides the user a wide variety of options for laminate design, which can be automatically incorporated with any finite element model that is composed of STAGS shell elements. Furthermore, the variable stiffness parameterization is formulated so that manufacturability can be assessed during the design process, plus new ideas using tow steering concepts can be easily integrated within the general framework of the elastic tailoring definitions. Details for the necessary implementation of the tow-steering definitions within the STAGS hierarchy is provided, and the format of the ply definitions is discussed in detail to provide easy access to the elastic tailoring choices. Integration of the automated STAGS solver with laminate design software has been demonstrated, so that the large design space generated by the tow-steering options can be traversed effectively. Several design problems are presented which confirm the usefulness of the design tool as well as further establish the potential of tow-steered plies for laminate design.

Jegley, Dawn C.↗

Aeroelastic Modeling of X-56A Stiff-Wing Configuration Flight Test Data

Aeroelastic stability and control derivatives for the X-56A Multi-Utility Technology Testbed (MUTT), in the stiff-wing configuration, were estimated from flight test data using the output-error method. Practical aspects of the analysis are discussed. The orthogonal phase-optimized multisine inputs provided excellent data information for aeroelastic modeling. Consistent parameter estimates were determined using output error in both the frequency and time domains. The frequency domain analysis converged faster and was less sensitive to starting values for the model parameters, which was useful for determining the aeroelastic model structure and obtaining starting values for the time domain analysis. Including a modal description of the structure from a finite element model reduced the complexity of the estimation problem and improved the modeling results. Effects of reducing the model order on the short period stability and control derivatives were investigated.

Grauer, Jared A.↗

A multigrid solver for the semiconductor equations

We present a multigrid solver for the exponential fitting method. The solver is applied to the current continuity equations of semiconductor device simulation in two dimensions. The exponential fitting method is based on a mixed finite element discretization using the lowest-order Raviart-Thomas triangular element. This discretization method yields a good approximation of front layers and guarantees current conservation. The corresponding stiffness matrix is an M-matrix. 'Standard' multigrid solvers, however, cannot be applied to the resulting system, as this is dominated by an unsymmetric part, which is due to the presence of strong convection in part of the domain. To overcome this difficulty, we explore the connection between Raviart-Thomas mixed methods and the nonconforming Crouzeix-Raviart finite element discretization. In this way we can construct nonstandard prolongation and restriction operators using easily computable weighted L(exp 2)-projections based on suitable quadrature rules and the upwind effects of the discretization. The resulting multigrid algorithm shows very good results, even for real-world problems and for locally refined grids.

Bachmann, Bernhard↗

Stiffness and Fatigue Life Estimator for Polymer Composite Laminates Using Machine Learning

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, Python-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 established. Results show that both neural net type algorithms can provide an excellent estimate of initial laminate stiffness as well as fatigue life of eight-ply symmetric polymer matrix composite laminates (PMCs). RNNs are better able to capture the shape of the fatigue curve of a laminate. The resulting tool and GUI can be very useful for system level studies to obtain an estimate of desired properties and life of PMC composite laminates. Further, the associated surrogate models can 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 micromechanics-based multiscale analyses a viable industrial tool for large scale structural problems.

multiscale analysis↗

Stiffness and Fatigue Life Estimator for Polymer Composite Laminates Using Machine Learning

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, Python-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 established. Results show that both neural net type algorithms can provide an excellent estimate of initial laminate stiffness as well as fatigue life of eight-ply symmetric polymer matrix composite laminates (PMCs). RNNs are better able to capture the shape of the fatigue curve of a laminate. The resulting tool and GUI can be very useful for system level studies to obtain an estimate of desired properties and life of PMC composite laminates. Further, the associated surrogate models can 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 micromechanics-based multiscale analyses a viable industrial tool for large scale structural problems.

multiscale analysis↗

Influence of thermal distortion on the anomalous behavior of a gravity gradient satellite

Consideration of the problem of determining the influence of solar heating on the low-frequency dynamic stability of a gravity-gradient satellite. The scope of the effort includes the development of a computer program that uses a quasi-static approach for describing a satellite's unsteady orbital motion. The essential feature of the quasi-static approach is that the satellite's long, slender booms reach their thermal equilibrium position solely as a result of thermal bending without consideration of the dynamic effects of boom vibration. A nonlinear analytical model and a corresponding computer program have been developed to study the effect of absorptivity, sun angle, and hinge stiffness on the stability of a satellite's motion.

Goldman, R. L.↗

Some Divergence Characteristics of Low-Aspect-Ratio Wings at Transonic and Supersonic Speeds

The problem of chordwise, or camber, divergence at transonic and supersonic speeds is treated with primary emphasis on slender delta wings having a cantilever support at the trailing edge. Experimental and analytical results are presented for four wing models having apex half-angles of 5 deg, 10 deg, 15 deg, and 20 deg. A Mach number range from 0.8 to 7.3 is covered. The analytical results include calculations based on small-aspect-ratio theory, lifting-surface theory, and strip theory. A closed-form solution of the equilibrium equation is given, which is based on low-aspect-ratio theory but which applies only to certain stiffness distributions. Also presented is an iterative procedure for use with other aerodynamic theories and with arbitrary stiffness distribution.

Woolston, Donald S.↗

Mass-additive modal test method for verification of constrained structural models

A method for deriving constrained or fixed-base modes and frequencies from free-free modes of a structure with mass-loaded boundaries is developed. Problems associated with design and development of test fixtures can be avoided with such an approach. The analytical methodology presented is used to assess applicability of the mass-additive method for three types of structures and to determine the accuracy of derived constrained modes and frequencies. Results show that mass loading of the boundaries enables local interface modes to be measured within a desired frequency bandwidth, thus allowing constrained modes to be derived with considerably fewer free-free modes than for unloaded boundaries. Good convergence was obtained for a simple beam and a truss-like Shuttle payload, both of which had well-spaced modes and stiff interface support structures. Slow convergence was obtained for a space station module prototype, a shell-like structure having high modal density.

Admire, John R.↗

The Dynamics of Solid Propellant Rocket Motors

From the point of view of dynamics, a solid propellant rocket motor is a unique structure in that it is composed of a substantial mass of propellant material case-bonded to a relatively massless, thin-walled cylinder. The mechanical properties of the propellant are such that it contributes little to the stiffness of the composite structure but it does contribute to the dynamic characteristics of the structure because of its mass. Furthermore, due to the viscoelastic character of the propellant, it can be expected to provide considerable damping to the system. At this point in the development of the art of design of solid propellant rocket motors there are no clear-cut or well-founded methods to evaluate quantitatively the contributions of the propellant to the dynamic responses of the composite structure. It is clear that unless such methods are devised, it will be difficult to arrive at accurate missile designs. In this paper a survey of the dynamic problems of solid propellant rocket motors is presented starting from the simplest model thereof and proceeding, step-by-step, to the consideration of more sophisticated and realistic models. The consideration is restricted to infinitesimal deformations of propellant grains with linear mechanical properties. Substantial progress has been achieved towards the solution of many important dynamical problems. We shall attempt to summarize the pertinent developments and indicate along which lines we feel future study should proceed. A few illustrative solutions are included in areas wherein progress has been substantial.

DYNAMIC RESPONSE↗

Graph theory approach to the eigenvalue problem of large space structures

Graph theory is used to obtain numerical solutions to eigenvalue problems of large space structures (LSS) characterized by a state vector of large dimensions. The LSS are considered as large, flexible systems requiring both orientation and surface shape control. Graphic interpretation of the determinant of a matrix is employed to reduce a higher dimensional matrix into combinations of smaller dimensional sub-matrices. The reduction is implemented by means of a Boolean equivalent of the original matrices formulated to obtain smaller dimensional equivalents of the original numerical matrix. Computation time becomes less and more accurate solutions are possible. An example is provided in the form of a free-free square plate. Linearized system equations and numerical values of a stiffness matrix are presented, featuring a state vector with 16 components.

Reddy, A. S. S. R.↗

Toward Ultralight High Strength Structural Materials via Collapsed Carbon Nanotube Bonding

The growing commercial availability of carbon nanotube (CNT) macro-assemblies such as sheet and yarn is making their use in structural composite components increasingly feasible. However, the mechanical properties of these materials continue to trail those of state-of-the-art carbon fiber composites due to relatively weak inter-tube load transfer. Forming covalent links between adjacent CNTs promises to mitigate this problem, but it has proven difficult in practice to introduce them chemically within densified and aligned CNT materials due to their low permeability. To avoid this limitation, this work explores the combination of pulsed electrical current, temperature, and pressure to introduce inter-CNT bonds. Reactive molecular dynamics simulations identify the most probable locations, configurations, and conditions for inter-nanotube bonds to form. This process is shown to introduce covalent linkages within the CNT material that manifest as improved macroscale mechanical properties. The magnitude of this effect increases with increasing levels of prealignment of the CNT material, promising a new synthesis pathway to ultralight structural materials with specific strengths and stiffnesses exceeding 1 and 100 GPa/(g/cu.cm), respectively.

Jensen, Benjamin D.↗

Parametric studies of stitching effectiveness for preventing substructure disbond

A methodology is desired that will allow a designer to select appropriate amounts of through-thickness reinforcement needed to meet design requirements. The goal is to use a relatively simple analysis to minimize the amount of testing that needs to be performed, and to make test results from simple configurations applicable to more general structures. Using this methodology one should be able to optimize the selection of stitching materials, the weight of the yarn, and the stitching density. The analysis approach is to treat substructure disbond as a crack propagation problem. In this approach, the stitches have little influence until a delamination begins to grow. Once the delamination reaches, or extends beyond a stitch, the stitch serves to reduce the strain-energy-release-rate (G) at the crack tip for a given applied load. The reduced G can then be compared to the unstitched materials toughness to predict the load required to further extend the crack. The current model treats the stitch as a simple spring which responds to displacements in the vertical (through-thickness) direction. In concept, this approach is similar to that proposed by other authors. Test results indicate that the model should be refined to include the shearing stiffness of the stitch. The strain-energy-release-rate calculations are performed using a code which uses interconnected higher-order plates to model built-up composite cross-sections. When plates are stacked vertically, the interfacial tractions between the plates can be computed. The plate differential equations are solved in closed-form. The code, called SUBLAM, was developed as part of this section in one dimension. Because of this limitation, rows of stitches are treated as a two-dimensional sheet. The spring stiffness of a row of stitches can be estimated from the stitch material, weight, and density. As a practical and conservative approach, we can assume that the stitch is bonded until a crack passes the stitch location. After the crack passes, it is fully bonded. A series of tests were performed to exercise this methodology and incorporated an attached flange such that the sudden change in thickness initiated a delamination. The analysis was used to estimate the materials' critical G from that of the unstitched specimens. With this data, a prediction was made for the load required to delaminate the stitched specimens. Using the methodology, design charts have been created for simplified geometries. These charts give stitch force, along with G(sub 1) and G(sub 2) as as function of the stitch spring stiffness. Using the charts, it should be possible to determine the stitch spring stiffness and strength required to reduce the G to a desired level. From these parameters, the actual stitching material, weight, and density can be computed.

Flanagan, Gerry↗

Solution of the subsynchronous whirl problem in the high-pressure hydrogen turbomachinery of the Space Shuttle Main Engine

Subsynchronous whirl of the high-pressure fuel turbopump limited operation of the Space Shuttle Main Engine for some months in early 1976. The means by which this problem was successfully attacked is of particular interest to the rotor-dynamics community, not only because this machinery was the highest power-to-weight ratio known (77,000 hp...760 pounds), but because of the multiple forcing functions involved and the means, both analytical and experimental, which were utilized in separating variables, pointing toward successful solutions, and evaluating results. The general means of identifying fundamental characteristics, analyzing data, and conducting computer investigations are delineated. The results of analysis and testing are discussed. Since whirl inception occurs at a shaft speed greater than twice the first system critical, this was increased by stiffening the shaft and bearing supports; the system damping and system stiffness was additionally increased by proper interstage seal design to the extent that the instability threshold is now beyond the operating range.

Ek, M. C.↗

A magnetic damper for first mode vibration reduction in multimass flexible rotors

Many rotating machines such as compressors, turbines and pumps have long thin shafts with resulting vibration problems, and would benefit from additional damping near the center of the shaft. Magnetic dampers have the potential to be employed in these machines because they can operate in the working fluid environment unlike conventional bearings. An experimental test rig is described which was set up with a long thin shaft and several masses to represent a flexible shaft machine. An active magnetic damper was placed in three locations: near the midspan, near one end disk, and close to the bearing. With typical control parameter settings, the midspan location reduced the first mode vibration 82 percent, the disk location reduced it 75 percent and the bearing location attained a 74 percent reduction. Magnetic damper stiffness and damping values used to obtain these reductions were only a few percent of the bearing stiffness and damping values. A theoretical model of both the rotor and the damper was developed and compared to the measured results. The agreement was good.

Kasarda, M. E. F.↗

On the application of ENO scheme with subcell resolution to conservation laws with stiff source terms

Two approaches are used to extend the essentially non-oscillatory (ENO) schemes to treat conservation laws with stiff source terms. One approach is the application of the Strang time-splitting method. Here the basic ENO scheme and the Harten modification using subcell resolution (SR), ENO/SR scheme, are extended this way. The other approach is a direct method and a modification of the ENO/SR. Here the technique of ENO reconstruction with subcell resolution is used to locate the discontinuity within a cell and the time evolution is then accomplished by solving the differential equation along characteristics locally and advancing in the characteristic direction. This scheme is denoted ENO/SRCD (subcell resolution - characteristic direction). All the schemes are tested on the equation of LeVeque and Yee (NASA-TM-100075, 1988) modeling reacting flow problems. Numerical results show that these schemes handle this intriguing model problem very well, especially with ENO/SRCD which produces perfect resolution at the discontinuity.

Chang, Shih-Hung↗

Tools for Designing and Analyzing Structures

Structural Design and Analysis Toolset is a collection of approximately 26 Microsoft Excel spreadsheet programs, each of which performs calculations within a different subdiscipline of structural design and analysis. These programs present input and output data in user-friendly, menu-driven formats. Although these programs cannot solve complex cases like those treated by larger finite element codes, these programs do yield quick solutions to numerous common problems more rapidly than the finite element codes, thereby making it possible to quickly perform multiple preliminary analyses - e.g., to establish approximate limits prior to detailed analyses by the larger finite element codes. These programs perform different types of calculations, as follows: 1. determination of geometric properties for a variety of standard structural components; 2. analysis of static, vibrational, and thermal- gradient loads and deflections in certain structures (mostly beams and, in the case of thermal-gradients, mirrors); 3. kinetic energies of fans; 4. detailed analysis of stress and buckling in beams, plates, columns, and a variety of shell structures; and 5. temperature dependent properties of materials, including figures of merit that characterize strength, stiffness, and deformation response to thermal gradients

Luz, Paul L.↗

Rotor-to-stator rub vibration in centrifugal compressor

One example of excessive vibration encountered during loading of a centrifugal compressor train (H type compressor with HP casing) is discussed. An investigation was made of the effects of the dynamic load on the bearing stiffness and the rotor-bearing system critical speed. The high vibration occurred at a "threshold load," but the machine didn't run smoothly due to rubs even when it had passed through the threshold load. The acquisition and discussion of the data taken in the field as well as a description of the case history which utilizes background information to identify the malfunction conditions is presented. The analysis shows that the failures, including full reverse precession rub and exact one half subharmonic vibration, were caused by the oversize bearings and displacement of the rotor center due to foundation deformation and misalignment between gear shafts, etc. The corrective actions taken to alleviate excessive vibration and the problems which remain to be solved are also presented.

Gao, J. J.↗