Search NASA⌕ Search

SEARCH · Search NASA

Results for “Problem structuring methods”

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 1,099 records · Page 61

High temperature composite analyzer (HITCAN) user's manual, version 1.0

This manual describes 'how-to-use' the computer code, HITCAN (HIgh Temperature Composite ANalyzer). HITCAN is a general purpose computer program for predicting nonlinear global structural and local stress-strain response of arbitrarily oriented, multilayered high temperature metal matrix composite structures. This code combines composite mechanics and laminate theory with an internal data base for material properties of the constituents (matrix, fiber and interphase). The thermo-mechanical properties of the constituents are considered to be nonlinearly dependent on several parameters including temperature, stress and stress rate. The computation procedure for the analysis of the composite structures uses the finite element method. HITCAN is written in FORTRAN 77 computer language and at present has been configured and executed on the NASA Lewis Research Center CRAY XMP and YMP computers. This manual describes HlTCAN's capabilities and limitations followed by input/execution/output descriptions and example problems. The input is described in detail including (1) geometry modeling, (2) types of finite elements, (3) types of analysis, (4) material data, (5) types of loading, (6) boundary conditions, (7) output control, (8) program options, and (9) data bank.

Lackney, J. J.↗

Automatic Determination of the Conic Coronal Mass Ejection Model Parameters

Characterization of the three-dimensional structure of solar transients using incomplete plane of sky data is a difficult problem whose solutions have potential for societal benefit in terms of space weather applications. In this paper transients are characterized in three dimensions by means of conic coronal mass ejection (CME) approximation. A novel method for the automatic determination of cone model parameters from observed halo CMEs is introduced. The method uses both standard image processing techniques to extract the CME mass from white-light coronagraph images and a novel inversion routine providing the final cone parameters. A bootstrap technique is used to provide model parameter distributions. When combined with heliospheric modeling, the cone model parameter distributions will provide direct means for ensemble predictions of transient propagation in the heliosphere. An initial validation of the automatic method is carried by comparison to manually determined cone model parameters. It is shown using 14 halo CME events that there is reasonable agreement, especially between the heliocentric locations of the cones derived with the two methods. It is argued that both the heliocentric locations and the opening half-angles of the automatically determined cones may be more realistic than those obtained from the manual analysis

Pulkkinen, A.↗

Applications of the theory of optimal control of distributed-parameter systems to structural optimization

An extension of classical methods of optimal control theory for systems described by ordinary differential equations to distributed-parameter systems described by partial differential equations is presented. An application is given involving the minimum-mass design of a simply-supported shear plate with a fixed fundamental frequency of vibration. An optimal plate thickness distribution in analytical form is found. The case of a minimum-mass design of an elastic sandwich plate whose fundamental frequency of free vibration is fixed. Under the most general conditions, the optimization problem reduces to the solution of two simultaneous partial differential equations involving the optimal thickness distribution and the modal displacement. One equation is the uniform energy distribution expression which was found by Ashley and McIntosh for the optimal design of one-dimensional structures with frequency constraints, and by Prager and Taylor for various design criteria in one and two dimensions. The second equation requires dynamic equilibrium at the preassigned vibration frequency.

Armand, J. P.↗

Numerical simulation of shock/turbulent boundary layer interaction

Most flows of aerodynamic interest are compressible and turbulent. However, our present knowledge on the structures and mechanisms of turbulence is mostly based on incompressible flows. In the present work, compressibility effects in turbulent, high-speed, boundary layer flows are systematically investigated using the Direct Numerical Simulation (DNS) approach. Three-dimensional, time-dependent, fully nonlinear, compressible Navier-Stokes equations were numerically integrated by high-order finite-difference methods; no modeling for turbulence is used during the solution because the available resolution is sufficient to capture the relevant scales. The boundary layer problem deals with fully-turbulent compressible flows over flat geometries. Apart from its practical relevance to technological flows, turbulent compressible boundary layer flow is the simplest experimentally realizable turbulent compressible flow. Still, measuring difficulties prohibit a detailed experimental description of the flow, especially in the near-wall region. DNS studies provide a viable means to probe the physics of compressible turbulence in this region. The focus of this work is to explore the paths of energy transfer through which compressible turbulence is sustained. The structural similarities and differences between the incompressible and compressible turbulence are also investigated. The energy flow patterns or energy cascades are found to be directly related to the evolution of vortical structures which are generated in the near-wall region. Near-wall structures, and mechanisms which are not readily accessible through physical experiments are analyzed and their critical role on the evolution and the behavior of the flow is documented extensively.

Biringen, Sedat↗

Assessment of impact damage of composite rocket motor cases

This contract reviewed the available literature on mechanisms of low velocity impact damage in filament wound rocket motor cases, MDE methods to quantify damage, critical coupon level test methods, manufacturing and material process variables and empirical and analytical modeling off impact damage. The critical design properties for rocket motor cases are biaxial hoop and axial tensile strength. Low velocity impact damage is insidious because it can create serious nonvisible damage at very low impact velocities. In thick rocket motor cases the prevalent low velocity impact damage is fiber fracture and matrix cracking adjacent to the front face. In contrast, low velocity loading of thin wall cylinders induces flexure, depending on span length and the flexure induces delamination and tensile cracking on the back face wall opposed to impact occurs due to flexural stresses imposed by impact loading. Important NDE methods for rocket motor cases are non-contacting methods that allow inspection from one side. Among these are vibrothermography, and pulse-echo methods based on acoustic-ultrasonic methods. High resolution techniques such as x-ray computed tomography appear to have merit for accurate geometrical characterization of local damage to support development of analytical models of micromechanics. The challenge of coupon level testing is to reproduce the biaxial stress state that the full scale article experiences, and to determine how to scale the composite structure to model full sized behavior. Biaxial tensile testing has been performed by uniaxially tensile loading internally pressurized cylinders. This is experimentally difficult due to gripping problems and pressure containment. Much prior work focused on uniaxial tensile testing of model filament wound cylinders. Interpretation of the results of some studies is complicated by the fact that the fabrication process did not duplicate full scale manufacturing. It is difficult to scale results from testing subscale cylinders since there are significant differences in out time of the resins relative to full scale cylinder fabrication, differences in hoop fiber tensioning and unsatisfactory coupon configurations. It appears that development of a new test method for subscale cylinders is merited. Damage tolerance may be improved by material optimization that uses fiber treatments and matrix modifications to control the fiber matrix interface bonding. It is difficult to develop process optimization in subscale cylinders without also modeling the longer out times resins experience in full scale testing. A major breakthrough in characterizing the effect of impact damage on residual strength, and understanding how to scale results of subscale evaluations, will be a sound micromechanical model that described progressive failure of the composite. Such models will utilize a three dimensional stress analysis due to the complex nature of low velocity impact stresses in thick composites. When these models are coupled with non-contact NDE methods that geometrically characterize the damage and acoustic methods that characterize the effective local elastic properties, accurate assessment of residual strength from impact damage may be possible. Directions for further development are suggested.

Paris, Henry G.↗

Analysis methods and preliminary design study

The shuttle orbiter body flap was studied to define loads and temperatures and to assess various structural concepts. Laminated-structure analysis capability was required. An appraisal of the available structural analysis programs, NASTRAN and SPAR, indicates that neither had all the capabilities required. The thermal analysis program, MITAS solved the problem once a model was generated but model generation and verification was laborious and transfer of temperatures to the structural program usually required interpolation. Therefore it was decided to incorporate both capabilities in the SPAR finite element program.

Blackburn, C. L.↗

The Pack Method for Compressive Tests of Thin Specimens of Materials Used in Thin-Wall Structures

The strength of modern lightweight thin-wall structures is generally limited by the strength of the compression members. An adequate design of these members requires a knowledge of the compressive stress-strain graph of the thin-wall material. The "pack" method was developed at the National Bureau of Standards with the support of the National Advisory Committee for Aeronautics to make possible a determination of compressive stress-strain graphs for such material. In the pack test an odd number of specimens are assembled into a relatively stable pack, like a "pack of cards." Additional lateral stability is obtained from lateral supports between the external sheet faces of the pack and outside reactions. The tests seems adequate for many problems in structural research.

Aitchison, C S↗

Studies in tilt rotor VTOL aircraft aeroelasticity, volume 2

Two methods for natural mode vibration analysis are discussed. The first consists of a direct approach based on a finite element representation of the complete structure as an entity. The mass and stiffness matrices for the complete structure are assembled by properly combining the mass and stiffness matrices of the individual elements into which the structure has been divided. The second approach is that of component mode synthesis. This method is based on the concept of synthesizing the natural modes of the complete structure from modes of conveniently difined substructures, or components, into which the structure has been partitioned. In this way the expedient of reducing the system degrees of freedom, and thus the size of the eigenvalue problem, can be introduced by partial modal synthesis.

Kvaternik, R. G.↗

Computational fracture mechanics

Numerical procedures for accurate determination of elastic stress intensity factors for the general two-dimensional crack problem are reviewed. The elastic perfectly plastic state of crack tip deformation is studied by a finite element procedure. Elastic-plastic fields in the immediate vicinity of a crack tip are determined numerically by finite element procedures based on asymptotic studies of crack tip singularities in plastic materials. The small-scale yielding problem is modeled, and expressions for crack tip opening displacement, shear singularity amplitude, and plastic zone extent are derived. A finite element solution to the large-scale yielding of a circumferentially cracked round tension bar is obtained. The three-dimensional aspects of flawed structures and numerical methods of treating them are studied. Ductile fracture mechanisms, in particular crack tip fracture on the microscale, are discussed.

Rice, J. R.↗

Use of asymptotic methods in vibration analysis

The derivation of dynamic differential equations, suitable for studying the vibrations of rotating, curved, slender structures was examined, and the Hamiltonian procedure was advocated for this purpose. Various reductions of the full system are displayed, which govern the vibrating troposkien when various order of magnitude restrictions are placed on important parameters. Possible advantages of the WKB asymptotic method for solving these classes of problems are discussed. A special case of this method is used illustratively to calculate eigenvalues and eigenfunctions for a flat turbine blade with small flexural stiffness.

Ashley, H.↗

Application of a systematic finite-element model modification technique to dynamic analysis of structures

A systematic finite-element model modification technique has been applied to two small problems and a model of the main wing box of a research drone aircraft. The procedure determines the sensitivity of the eigenvalues and eigenvector components to specific structural changes, calculates the required changes and modifies the finite-element model. Good results were obtained where large stiffness modifications were required to satisfy large eigenvalue changes. Sensitivity matrix conditioning problems required the development of techniques to insure existence of a solution and accelerate its convergence. A method is proposed to assist the analyst in selecting stiffness parameters for modification.

Robinson, J. C.↗

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

N-body models

The experimental discovery of an overstability in the central regions of galaxies is briefly discussed, and numerical methods for integrating orbits are briefly addressed. The overstability manifests itself as a growing amplitude in the orbit of a galaxy's nucleus about its mass centroid. This finding may complicate studies of the topological properties of orbits and studies of the bifurcation structure of orbits. A sample problem is used to illustrate the importance of a Liouville theorem in N-body calculations.

Miller, R. H.↗

Robustness of solutions to a benchmark control problem

The stochastic robustness of solutions to a benchmark control design problem has been analyzed. The analysis quantifies the controllers' stability and performance robustness with structured uncertainties in up to six system parameters. The analysis provides insights about system response that are not readily derived from other robustness criteria, providing a common ground for judging controllers produced by alternative methods.

Stengel, Robert F.↗

Stochastic search in structural optimization - Genetic algorithms and simulated annealing

An account is given of illustrative applications of genetic algorithms and simulated annealing methods in structural optimization. The advantages of such stochastic search methods over traditional mathematical programming strategies are emphasized; it is noted that these methods offer a significantly higher probability of locating the global optimum in a multimodal design space. Both genetic-search and simulated annealing can be effectively used in problems with a mix of continuous, discrete, and integer design variables.

Hajela, Prabhat↗

On domains of convergence in optimization problems

Numerical optimization algorithms require the knowledge of an initial set of design variables. Starting from an initial design x(sup 0), improved solutions are obtained by updating the design iteratively in a way prescribed by the particular algorithm used. If the algorithm is successful, convergence is achieved to a local optimal solution. Let A denote the iterative procedure that characterizes a typical optimization algorithm, applied to the problem: Find x belonging to R(sup n) that maximizes f(x) subject to x belonging to Omega contained in R(sup n). We are interested in problems with several local maxima (x(sub j))(sup *), j=1, ..., m, in the feasible design space Omega. In general, convergence of the algorithm A to a specific solution (x(sub j))(sup *) is determined by the choice of initial design x(sup 0). The domain of convergence D(sub j) of A associated with a local maximum (x(sub j))(sup *) is a subset of initial designs x(sup 0) in Omega such that the sequence (x(sup k)), k=0,1,2,... defined by x(sup k+1) = A(x(sup k)), k=0,1,... converges to (x(sub j))(sup *). The set D(sub j) is also called the basin of attraction of (x(sub j))(sup *). Cayley first proposed the problem of finding the basin of attraction for Newton's method in 1897. It has been shown that the basin of attraction for Newton's method exhibits chaotic behavior in problems with polynomial objective. This implies that there may be regions in the feasible design space where arbitrarily close starting points will converge to different local optimal solutions. Furthermore, the boundaries of the domains of convergence may have a very complex, even fractal structure. In this paper we show that even simple structural optimization problems solved using standard gradient based (first order) algorithms exhibit similar features.

Diaz, Alejandro R.↗

A Neural Network Aero Design System for Advanced Turbo-Engines

An inverse design method calculates the blade shape that produces a prescribed input pressure distribution. By controlling this input pressure distribution the aerodynamic design objectives can easily be met. Because of the intrinsic relationship between pressure distribution and airfoil physical properties, a Neural Network can be trained to choose the optimal pressure distribution that would meet a set of physical requirements. Neural network systems have been attempted in the context of direct design methods. From properties ascribed to a set of blades the neural network is trained to infer the properties of an 'interpolated' blade shape. The problem is that, especially in transonic regimes where we deal with intrinsically non linear and ill posed problems, small perturbations of the blade shape can produce very large variations of the flow parameters. It is very unlikely that, under these circumstances, a neural network will be able to find the proper solution. The unique situation in the present method is that the neural network can be trained to extract the required input pressure distribution from a database of pressure distributions while the inverse method will still compute the exact blade shape that corresponds to this 'interpolated' input pressure distribution. In other words, the interpolation process is transferred to a smoother problem, namely, finding what pressure distribution would produce the required flow conditions and, once this is done, the inverse method will compute the exact solution for this problem. The use of neural network is, in this context, highly related to the use of proper optimization techniques. The optimization is used essentially as an automation procedure to force the input pressure distributions to achieve the required aero and structural design parameters. A multilayered feed forward network with back-propagation is used to train the system for pattern association and classification.

Sanz, Jose M.↗

Parallelization of the Flow Field Dependent Variation Scheme for Solving the Triple Shock/Boundary Layer Interaction Problem

A parallelized version of the Flowfield Dependent Variation (FDV) Method is developed to analyze a problem of current research interest, the flowfield resulting from a triple shock/boundary layer interaction. Such flowfields are often encountered in the inlets of high speed air-breathing vehicles including NASA's Hyper-X. In order to resolve the complex shock structure and to provide adequate resolution for boundary layer computations of the convective heat transfer from surfaces inside the inlet, models containing over 500,000 nodes are needed. Efficient parallelization of the computation is essential to obtaining the results in a timely manner. Results from different parallelization schemes, based upon multi-threading and message passing, as implemented on multiple processor supercomputers and on distributed workstations are compared.

Schunk, Greg↗