SEARCH · Engineering Papers
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.
Computer program synthesizes a passive network by minimizing the difference in desired and actual frequency response. The program solves for the critical points of the error function /weighted least squares fit between calculated and desired frequency response/ by the multivariable Newton-Raphson method with components constrained to an admissible region.
Computer program uses Bernoullis formula and Newton-Raphson method to provide steady state fluid flow analysis of line pressure drop in a system with six outlets for each of two main storage tanks. Program flexibility arises in the ease with which changes in the fluid line geometry can be made.
Computer program for transfer orbit trajectory optimization based on Lagrange multipliers, Newton-Raphson method, and power series
A Newton-Raphson method of iteration was used in evaluating the radial and axial projection of the distance between the ball center and the outer raceway groove curvature center (V and W). Fatigue life evaluations were made. The similar analysis of a conventional bearing can be directly obtained from the arched bearing analysis by simply letting the amount of arching be zero (g = 0) and not considering equations related to the unloaded half of the outer race. The analysis was applied to a 150-mm angular contact ball bearing. Results for life, contact loads, and angles are shown for a conventional bearing (g = 0) and two arched bearings (g = 0.127 mm (0.005 in.), and 0.254 mm (0.010 in.)). The results indicate that an arched bearing is highly desirable for high speed applications. In particular, for a DN value of 3 million (20,000 rpm) and an applied axial load of 4448 N (1000 lb), an arched bearing shows an improvement in life of 306 percent over that of a conventional bearing. At 4.2 million DN (28,000 rpm), the corresponding improvement is 340 percent. It was also found for low speeds, the arched bearing does not offer the advantages that it does for high speed applications.
Thrust-load analysis of a 150-mm angular contact ball bearing, taking into account centrifugal forces but omitting gyroscopics, elastohydrodynamics, and thermal effects. A Newton-Raphson method of iteration is used to evaluate the radial and axial projection of the distance between the ball center and the outer raceway groove curvature center. Fatigue life of the bearing is evaluated. Results for life, contact loads, and angles are given for a conventional bearing and two arched bearings.
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.
The longitudinal aerodynamic derivatives of the F-8 aircraft with supercritical wing were obtained from flight data by a parameter-extraction algorithm at Mach numbers of 0.8, 0.9, and 0.98. A set of derivatives were obtained from which calculated aircraft responses were correlated almost identically with actual flight responses. In general, the trends of the extracted derivatives obtained by the algorithm agreed with those obtained by a Newton-Raphson method and with preliminary data from the Langley 8-foot transonic pressure tunnel. The wind-tunnel damping derivatives were, however, substantially higher than the converged damping derivatives possibly because of Reynolds number differences between flight and model tests.
A modified Newton-Raphson method has been used to determine the parameters in equations describing the vortex flow to obtain a best match with flight measurements of the flow behind the C-5A airplane. The flight measurements were made using a specially instrumented T-33 airplane which passed as closely as possible to the centers of the trailing vortices at several distances behind the C-5A airplane. The flight measurements were transformed to flow velocity relative to an inertial frame of reference. The assumed form of the flow consisted of the superposition of two counterrotating, finite core vortices. The positions of the vortex centers, their total circulation, the effective eddy viscosity and measurement bias were the parameters adjusted. The assumed form of vortex flow fit well the measured velocities for the numerous sets of data, both flaps up and down for the C-5A airplane. The resulting values of total circulation, however, were about two-thirds that expected of a wing with an elliptical loading. A partial explanation of the less than expected circulation is a dip in the spanwise lift distribution at the airplane's center line. The distance between the trailing vortices at the smallest times encountered is somewhat less than that expected for an elliptical wing loading.
A theoretical investigation was conducted to determine the importance of the lateral acceleration (beta) derivatives in the extraction of lateral-directional stability derivatives for swept wing airplanes at high angles of attack. Representative values of lateral acceleration derivatives in yaw and roll (Cn beta and Cl beta) were used in a computer program to generate representative flight motions at several angles of attack and altitudes. The computer-generated motions were then subjected to a parameter identification process based on a modified Newton-Raphson method. Two identification techniques were evaluated, one which included the beta derivatives and one which neglected them. The results of the study indicate that omission of the beta derivatives from mathematical models used in the derivative-extraction techniques can produce erroneous values for the lateral-directional stability derivatives particularly at high angles of attack, where the beta derivatives are large. The largest errors occur in the dynamic derivatives, but large errors may also occur in the static derivatives for cases in which the beta derivatives have large effects on the flight motions of the airplane. In addition, the resulting identified mathematical models provide poor motion prediction as well as erroneous predictions of dynamic modal characteristics. These results strongly indicate that the effects of beta derivatives should be considered in any attempt to extract lateral-directional aerodynamic parameters at high angles of attack.
An efficient algorithm is presented for the solution of Kepler's equation f(E)=E-M-e sin E=0, where e is the eccentricity, M the mean anomaly and E the eccentric anomaly. This algorithm is based on simple initial approximations that are cubics in M, and an iterative scheme that is a slight generalization of the Newton-Raphson method. Extensive testing of this algorithm has been performed on the UNIVAC 1108 computer. Solutions for 20,000 pairs of values of e and M show that for single precision, 42.0% of the cases require one iteration, 57.8% two and 0.2% three. For double precision one additional iteration is required.
A system identification technique is developed for classes of distributed systems using finite element approximations. Vibrating systems represented by partial differential equations have physical parameters associated with mass, stiffness, and damping distributions which need to be known in order to properly control and design mathematical models of the system. In order to identify these parameters a weighted least-squares algorithm and modified Newton-Raphson method is used for the identification process. The theory and technique is demonstrated by estimating the system parameters of a vibrating cantilever beam made up of several different structural properties.
An inviscid discrete vortex model, with newly derived expressions for the tangential velocity imposed at the separation points, is used to investigate the symmetric and asymmetric vortex separation on cones and tangent ogives. The circumferential locations of separation are taken from experimental data. Based on a slender body theory, the resulting simultaneous nonlinear algebraic equations in a cross-flow plane are solved with Broyden's modified Newton-Raphson method. Total force coefficients are obtained through momentum principle with new expressions for nonconical flow. It is shown through the method of function deflation that multiple solutions exist at large enough angles of attack, even with symmetric separation points. These additional solutions are asymmetric in vortex separation and produce side force coefficients which agree well with data for cones and tangent ogives.
Presented are the research results from the research grant entitled: Active Control of Robot Manipulators, sponsored by the Goddard Space Flight Center (NASA) under grant number NAG-780. This report considers a class of robot manipulators based on the closed-kinematic chain mechanism (CKCM). This type of robot manipulators mainly consists of two platforms, one is stationary and the other moving, and they are coupled together through a number of in-parallel actuators. Using spatial geometry and homogeneous transformation, a closed-form solution is derived for the inverse kinematic problem of the six-degree-of-freedom manipulator, built to study robotic assembly in space. Iterative Newton Raphson method is employed to solve the forward kinematic problem. Finally, the equations of motion of the above manipulators are obtained by employing the Lagrangian method. Study of the manipulator dynamics is performed using computer simulation whose results show that the robot actuating forces are strongly dependent on the mass and centroid locations of the robot links.
The forced steady response of a single degree of freedom system involving a large nonlinearity, represented by unsymmetric piecewise-linear stiffness, is determined by a harmonic balance Newton-Raphson method with the application of the fast Fourier transformation (FFT) algorithm. All possible subharmonic, harmonic, and superharmonic responses are sought. The responses of two systems involving subharmonic and superharmonic dominant motions are calculated by the newly developed method and compared with previously published findings. The results obtained by using this method reveal the details of the response of the systems more efficiently than previous methods.
This report presents results from the research grant entitled Active Control of Robot Manipulators, funded by the Goddard Space Flight Center, under Grant NAG5-780, for the period July 1, 1988 to January 1, 1989. An analysis is presented of a 6 degree-of-freedom robot end-effector built to study telerobotic assembly of NASA hardware in space. Since the end-effector is required to perform high precision motion in a limited workspace, closed-kinematic mechanisms are chosen for its design. A closed-form solution is obtained for the inverse kinematic problem and an iterative procedure employing Newton-Raphson method is proposed to solve the forward kinematic problem. A study of the end-effector workspace results in a general procedure for the workspace determination based on link constraints. Computer simulation results are presented.
Based on the Newton-Raphson method, this paper presents recursive algorithms that are rapidly convergent and more stable for modeling the equivalent continuous-time (discrete-time) model from the available discrete-time (continuous-time) model for a fixed sampling period. The newly developed recursive algorithms relax the constraints imposed upon the existing model conversion algorithms, and, thus, enhance the applications of microprocessors and associated microelectronics to digital control systems. A practical example is presented to demonstrate the effectiveness of the proposed procedures.
An optical processor to solve partial differential equations for computational fluid dynamics applications is considered. This application is new and original for optical processors. The algorithms that are used are optical realizations of the Newton-Raphson method for nonlinear equations and a new optical LU direct decomposition and Gauss-Seidel iterative solution to the resultant linear algebraic equations. These algorithms are used to solve Burger's equation (a specific form of the momentum equation in fluid dynamics). The nonlinear equations provide 1-D velocity data at each time step. Simulation results of optical processing with these algorithms on computational fluid dynamics data is included.
This paper deals with a class of robot manipulators built based on the kinematic chain mechanism (CKCM). This class of CKCM manipulators consists of a fixed and a moving platform coupled together via a number of in-parallel actuators. A closed-form solution is derived for the inverse kinematic problem of a six-degre-of-freedom CKCM manipulator designed to study robotic applications in space. Iterative Newton-Raphson method is employed to solve the forward kinematic problem. Dynamics of the above manipulator is derived using the Lagrangian approach. Computer simulation of the dynamical equations shows that the actuating forces are strongly dependent on the mass and centroid of the robot links.