Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nonlinear programming”

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,009 records · Page 56

Global/local methods for probabilistic structural analysis

A probabilistic global/local method is proposed to reduce the computational requirements of probabilistic structural analysis. A coarser global model is used for most of the computations with a local more refined model used only at key probabilistic conditions. The global model is used to establish the cumulative distribution function (cdf) and the Most Probable Point (MPP). The local model then uses the predicted MPP to adjust the cdf value. The global/local method is used within the advanced mean value probabilistic algorithm. The local model can be more refined with respect to the g1obal model in terms of finer mesh, smaller time step, tighter tolerances, etc. and can be used with linear or nonlinear models. The basis for this approach is described in terms of the correlation between the global and local models which can be estimated from the global and local MPPs. A numerical example is presented using the NESSUS probabilistic structural analysis program with the finite element method used for the structural modeling. The results clearly indicate a significant computer savings with minimal loss in accuracy.

Millwater, H. R.↗

A History of Full-Scale Aircraft and Rotorcraft Crash Testing and Simulation at NASA Langley Research Center

This paper summarizes 2-1/2 decades of full-scale aircraft and rotorcraft crash testing performed at the Impact Dynamics Research Facility (IDRF) located at NASA Langley Research Center in Hampton, Virginia. The IDRF is a 240-ft.-high steel gantry that was built originally as a lunar landing simulator facility in the early 1960's. It was converted into a full-scale crash test facility for light aircraft and rotorcraft in the early 1970 s. Since the first full-scale crash test was preformed in February 1974, the IDRF has been used to conduct: 41 full-scale crash tests of General Aviation (GA) aircraft including landmark studies to establish baseline crash performance data for metallic and composite GA aircraft; 11 full-scale crash tests of helicopters including crash qualification tests of the Bell and Sikorsky Advanced Composite Airframe Program (ACAP) prototypes; 48 Wire Strike Protection System (WSPS) qualification tests of Army helicopters; 3 vertical drop tests of Boeing 707 transport aircraft fuselage sections; and, 60+ crash tests of the F-111 crew escape module. For some of these tests, nonlinear transient dynamic codes were utilized to simulate the impact response of the airframe. These simulations were performed to evaluate the capabilities of the analytical tools, as well as to validate the models through test-analysis correlation. In September 2003, NASA Langley closed the IDRF facility and plans are underway to demolish it in 2007. Consequently, it is important to document the contributions made to improve the crashworthiness of light aircraft and rotorcraft achieved through full-scale crash testing and simulation at the IDRF.

Jackson, Karen E.↗

International Space Station 2A Array Modal Analysis

On December 9th 2009, the International Space Station (ISS) 2A solar array mast experienced prolonged longeron shadowing during a Soyuz undocking. Analytical reconstruction of induced thermal and dynamic structural loads showed an exceedance of the mast buckling limit. Possible structural damage to the solar array mast could have occurred during this event. A Low fidelity video survey of the 2A mast showed no obvious damage of the mast longerons or battens. The decision was made to conduct an on-orbit dynamic test of the 2A array on December 18th, 2009. The test included thruster pluming on the array while photogrammetry data was recorded. The test was similar to other Dedicated Thruster Firings (DTFs) that were performed to measure structural frequency and damping of a solar array. Results of the DTF indicated lower frequency mast modes than model predictions, thus leading to speculation of mast damage. A detailed nonlinear analysis was performed on the 2A array model to assess possible solutions to modal differences. The setup of the parametric nonlinear trade study included the use of a detailed array model and the reduced mass and stiffness matrices of the entire ISS being applied to the array interface. The study revealed that the array attachment structure is nonlinear and thus was the source of error in the model prediction of mast modes. In addition, a detailed study was performed to determine mast mode sensitivity to mast longeron damage. This sensitivity study was performed to assess if the ISS program has sufficient instrumentation for mast damage detection.

Laible, Michael↗

Electric Field Conjugation in the Presence of Model Uncertainty

The Wide-Field Infrared Survey Telescope (WFIRST) is a 2.4-meter-diameter space telescope NASA program. The payload will include a coronagraph instrument (CGI). The CGI designs under development use deformable mirrors (DM) to create a point spread function (PSF) with a dark region around the obscured star object. Electric field conjugation (EFC) is an iterative nonlinear optimization procedure that uses measurements of the electric field at the image to determine the DM displacements to modify the PSF to create the region of high contrast in the image. EFC requires a numerical model of the coronagraph to calculate the Jacobian of the system, which is used, along with regularization, to solve for the DM displacements for each iteration of the nonlinear optimization. Ideally, the coronagraph is aligned and calibrated, and the calibration data are used in the numerical model for calculating the Jacobian. However, calibration and alignment measurements always contain uncertainty resulting in calibration error. Therefore, the Jacobian calculated from the numerical model is not an exact representation of the physical coronagraph, and the resulting DM solution for an EFC iteration does not have the exact impact on the electric field of the coronagraph as predicted by the EFC. The result is slow convergence, and, as will be shown, drives the use of more restrictive regularization. Using Monte Carlo trials, we investigate the effect of calibration error on EFC convergence and regulatization. Comparison to results from the High Contrast Imaging Testbed Hybrid Lyot Coronagraph are also presented.

Marx, David↗

Importance of helicopter dynamics to the mathematical model of the helicopter

A mathematical model of the helicopter requires appropriate representation of the constituent elements of rotor dynamics. General-purpose programs that model a variety of configurations for a broad range of operating conditions result in varying and incompatible levels of sophistication. Analysis of specific dynamic problems facilitates the identification of configuration parameters which determine system behavior. For the present analysis, the nonlinear equations of a torsionally rigid hingeless rotor are linearized about an equilibrium condition to determine flap-lag stability characteristics in hover. A collocation method was used to obtain the coupled natural frequencies and modes. These modes allow exact treatment of the effect of elastic coupling which more than compensates for the destabilizing inertial coupling. The sensitivity of damping to the number of modes was found to be small, and reasonable accuracy was obtained the first flapwise and edgewise coupled modes. The range of destabilizing precone was found to be small.

White, W. F., Jr.↗

Model verification of large structural systems

A computer program for the application of parameter identification on the structural dynamic models of space shuttle and other large models with hundreds of degrees of freedom is described. Finite element, dynamic, analytic, and modal models are used to represent the structural system. The interface with math models is such that output from any structural analysis program applied to any structural configuration can be used directly. Processed data from either sine-sweep tests or resonant dwell tests are directly usable. The program uses measured modal data to condition the prior analystic model so as to improve the frequency match between model and test. A Bayesian estimator generates an improved analytical model and a linear estimator is used in an iterative fashion on highly nonlinear equations. Mass and stiffness scaling parameters are generated for an improved finite element model, and the optimum set of parameters is obtained in one step.

Lee, L. T.↗

A comparison of linear acoustic theory with experimental noise data for a small-scale hovering rotor

Linear acoustic calculations based on full aerodynamic data as input are presented and compared with measured cases reported by Boxwell et al. (1978). The full aerodynamic data are obtained using three programs giving radial loading, chordwise loading, and chordwise position of transition. It is shown that in the theoretical results the most significant noise source mechanism is due to blade thickness. Thus the conclusions of Boxwell et al. as to the importance of nonlinearities around the blades are upheld. These conclusions concern the width, shape and the level of the acoustic pressure calculated from linear acoustic theory. Some of the approximations involved in the application of acoustic analogy using quadrupole sources are discussed. It is necessary that the near- and far-field problems of rotating blades be treated together as shown for the case of an oscillating sphere.

Farassat, F.↗

The determination of atmospheric trace gases using the correlation interferometer technique

The theoretical feasibility inherent in the correlation interferometry technique is demonstrated. High signal to noise and increased information content in a given partial interferogram are attractive features of the technique. The separation of interfering gas species is overcome by a precise measurement program for spectral line parameters of methane and other species. There are optical, mechanical, and electronic difficulties with the interferometry hardware design. Flight data from Cessna 402 over southeastern Virginia during April 1978 yielded amplitude fluctuations that are explainable only as internal nonlinearities. Misalignment of the AGC/interferometer optical axis also contributed to a large error source in the flight data. It is concluded that no useful data on gas burdens was obtained from these flights.

Kindle, E. C.↗

Real-time hybrid computer simulation of a small turboshaft engine and control system

The development of an analytical model of a small turboshaft engine designed for helicopter propulsion systems is described. The model equations were implemented on a hybrid computer system to provide a real time nonlinear simulation of the engine performance over a wide operating range. The real time hybrid simulation of the engine was used to evaluate a microprocessor based digital control module. This digital control module was developed as part of an advanced rotorcraft control program. After tests with the hybrid engine simulation the digital control module was used to control a real engine in an experimental program. A hybrid simulation of the engine's electrical hydromechanical control system was developed. This allowed to vary the fuel flow and torque load inputs to the hybrid engine simulation for simulating transient operation. A steady-state data and the experimental tests are compared. Analytical model equations, analog computer diagrams, and a digital computer flow chart are included.

Hart, C. E.↗

Deterministic Multiaxial Creep and Creep Rupture Enhancements for CARES/Creep Integrated Design Code

High temperature and long duration applications of monolithic ceramics can place their failure mode in the creep rupture regime. A previous model advanced by the authors described a methodology by which the creep rupture life of a loaded component can be predicted. That model was based on the life fraction damage accumulation rule in association with the modified Monkman-Grant creep rupture criterion. However, that model did not take into account the deteriorating state of the material due to creep damage (e.g., cavitation) as time elapsed. In addition, the material creep parameters used in that life prediction methodology, were based on uniaxial creep curves displaying primary and secondary creep behavior, with no tertiary regime. The objective of this paper is to present a creep life prediction methodology based on a modified form of the Kachanov-Rabotnov continuum damage mechanics (CDM) theory. In this theory, the uniaxial creep rate is described in terms of sum, temperature, time, and the current state of material damage. This scalar damage state parameter is basically an abstract measure of the current state of material damage due to creep deformation. The damage rate is assumed to vary with stress, temperature, time, and the current state of damage itself. Multiaxial creep and creep rupture formulations of the CDM approach are presented in this paper. Parameter estimation methodologies based on nonlinear regression analysis are also described for both, isothermal constant stress states and anisothermal variable stress conditions This creep life prediction methodology was preliminarily added to the integrated design code CARES/Creep (Ceramics Analysis and Reliability Evaluation of Structures/Creep), which is a postprocessor program to commercially available finite element analysis (FEA) packages. Two examples, showing comparisons between experimental and predicted creep lives of ceramic specimens, are used to demonstrate the viability of Ns methodology and the CARES/Creep program.

Jadaan, Osama M.↗

Deterministic and Probabilistic Creep and Creep Rupture Enhancement to CARES/Creep: Multiaxial Creep Life Prediction of Ceramic Structures Using Continuum Damage Mechanics and the Finite Element Method

High temperature and long duration applications of monolithic ceramics can place their failure mode in the creep rupture regime. A previous model advanced by the authors described a methodology by which the creep rupture life of a loaded component can be predicted. That model was based on the life fraction damage accumulation rule in association with the modified Monkman-Grant creep ripture criterion However, that model did not take into account the deteriorating state of the material due to creep damage (e.g., cavitation) as time elapsed. In addition, the material creep parameters used in that life prediction methodology, were based on uniaxial creep curves displaying primary and secondary creep behavior, with no tertiary regime. The objective of this paper is to present a creep life prediction methodology based on a modified form of the Kachanov-Rabotnov continuum damage mechanics (CDM) theory. In this theory, the uniaxial creep rate is described in terms of stress, temperature, time, and the current state of material damage. This scalar damage state parameter is basically an abstract measure of the current state of material damage due to creep deformation. The damage rate is assumed to vary with stress, temperature, time, and the current state of damage itself. Multiaxial creep and creep rupture formulations of the CDM approach are presented in this paper. Parameter estimation methodologies based on nonlinear regression analysis are also described for both, isothermal constant stress states and anisothermal variable stress conditions This creep life prediction methodology was preliminarily added to the integrated design code CARES/Creep (Ceramics Analysis and Reliability Evaluation of Structures/Creep), which is a postprocessor program to commercially available finite element analysis (FEA) packages. Two examples, showing comparisons between experimental and predicted creep lives of ceramic specimens, are used to demonstrate the viability of this methodology and the CARES/Creep program.

Jadaan, Osama M.↗

Computational Simulation of the High Strain Rate Tensile Response of Polymer Matrix Composites

A research program is underway to develop strain rate dependent deformation and failure models for the analysis of polymer matrix composites subject to high strain rate impact loads. Under these types of loading conditions, the material response can be highly strain rate dependent and nonlinear. State variable constitutive equations based on a viscoplasticity approach have been developed to model the deformation of the polymer matrix. The constitutive equations are then combined with a mechanics of materials based micromechanics model which utilizes fiber substructuring to predict the effective mechanical and thermal response of the composite. To verify the analytical model, tensile stress-strain curves are predicted for a representative composite over strain rates ranging from around 1 x 10(exp -5)/sec to approximately 400/sec. The analytical predictions compare favorably to experimentally obtained values both qualitatively and quantitatively. Effective elastic and thermal constants are predicted for another composite, and compared to finite element results.

Goldberg, Robert K.↗

Indirect synthesis of multidegree-of-freedom transient systems

The indirect synthesis method is developed and shown to be capable of leading a near-optimal design of multidegree-of-freedom and multidesign-element transient nonlinear dynamical systems. The basis of the approach is to select the open design parameters such that the response of the portion of the system being designed approximates the limiting performances solution. The limiting performance problem can be formulated as one of linear programming by replacing all portions of the system subject to transient disturbances by control forces and supposing that the remaining portions are linear as are the overall kinematic constraints. One then selects the design parameters that respond most closely to the limiting performance solution, which can be achieved by unconstrained curve-fitting techniques.

Chen, Y. H.↗

Computing Stress, Stability, And Vibration Of Shells

BOSOR4 computer program developed as comprehensive program for analysis of stress, stability, and vibration of complex, branched shells of revolution made of elastic materials. Used to analyze prismatic shells and panels. Performs large-deflection axisymmetric stress analysis, small-deflection nonsymmetric stress analysis, modal vibration analysis with axisymmetric nonlinear prestress included, and buckling analysis with axisymmetric or nonsymmetric prestress. One main advantage, provisions for such realistic engineering details as eccentric load paths, internal supports, arbitrary branching conditions, and "library" of wall constructions. Based on finite-difference energy method and offers very rapid convergence with increasing numbers of mesh points. Written in FORTRAN 77.

Bushnell, David↗

Transfer Functions Via Laplace- And Fourier-Borel Transforms

Approach to solution of nonlinear ordinary differential equations involves transfer functions based on recently-introduced Laplace-Borel and Fourier-Borel transforms. Main theorem gives transform of response of nonlinear system as Cauchy product of transfer function and transform of input function of system, together with memory effects. Used to determine responses of electrical circuits containing variable inductances or resistances. Also possibility of doing all noncommutative algebra on computers in such symbolic programming languages as Macsyma, Reduce, PL1, or Lisp. Process of solution organized and possibly simplified by algebraic manipulations reducing integrals in solutions to known or tabulated forms.

Can, Sumer↗

A Disturbance Rejection Approach to Actuator and Sensor Placement

For various reasons as discussed for instance in, the selection of actuator and sensor positions is still ad hoc. This is especially true for flexible structures where many candidate configurations can exist. This study is an attempt to make the selection process more methodical. One approach to actuator and sensor placement is to optimize a closed loop performance metric directly by selecting the actuators, sensors, and controller gains simultaneously. This direct approach makes sense if the desired closed loop performance is well defined. Since the individual actuator and sensor contributions to the closed loop performance metric is complex, the solution strategy usually employs non linear programming with many design and numerical iterations. A second approach is to select actuators and/or sensors based on open loop properties so that closed loop performance is indirectly optimized. Since the individual sensor and actuator contributions to the open loop metric is simple, nonlinear optimization is usually not needed. This approach will suggest efficient actuator and sensor configurations for any type of control law. The method suggested in this study falls into the latter class of approaches.

Lim, K. B.↗

Variable-Speed Simulation of a Dual-Clutch Gearbox Tiltrotor Driveline

This investigation explores the variable-speed operation and shift response of a prototypical two-speed dual-clutch transmission tiltrotor driveline in forward flight. Here, a Comprehensive Variable-Speed Rotorcraft Propulsion System Modeling (CVSRPM) tool developed under a NASA funded NRA program is utilized to simulate the drive system dynamics. In this study, a sequential shifting control strategy is analyzed under a steady forward cruise condition. This investigation attempts to build upon previous variable-speed rotorcraft propulsion studies by 1) including a fully nonlinear transient gas-turbine engine model, 2) including clutch stick-slip friction effects, 3) including shaft flexibility, 4) incorporating a basic flight dynamics model to account for interactions with the flight control system. Through exploring the interactions between the various subsystems, this analysis provides important insights into the continuing development of variable-speed rotorcraft propulsion systems.

DeSmidt, Hans↗

Application of symbolic/numeric matrix solution techniques to the NASTRAN program

The matrix solving algorithm of any finite element algorithm is extremely important since solution of the matrix equations requires a large amount of elapse time due to null calculations and excessive input/output operations. An alternate method of solving the matrix equations is presented. A symbolic processing step followed by numeric solution yields the solution very rapidly and is especially useful for nonlinear problems.

Buturla, E. M.↗