Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “torque”

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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Error field measurements with rotating RMP fields for DIII-D H-mode

3D magnetic sensors are employed to identify the amplitude and toroidal phase of error fields (EF) by analyzing the torque balance for magnetic islands entrained by rotating resonant magnetic perturbations (RMPs) in DIII-D H-mode plasmas. This technique of torque balance allows for efficient error field identification, offering a valuable tool for scenario-specific and optimized error field compensation (EFC) and requires only magnetic diagnostics. The torque balance used in this work includes the contribution from electromagnetic torque due to error fields, wall response, RMP fields, and viscous. Results show that viscous torque plays a crucial role, particularly during locked modes and H-mode plasmas, ensuring accurate data fits with lower residuals. The torque balance technique reveals that the L- and H-mode plasmas have distinct EF configurations, and consistent with a model-based EF analysis including MHD response in IPEC and the SURFMN EF simulation. This technique shows great robustness in measuring the intrinsic EF amplitude regardless of its amplitude or toroidal phase. Repeated discharges with EFC disparities exhibit consistent results of intrinsic error field within a reasonable range near the “standard” error field compensation. Additionally, the use of a rotating n = 1 resonant magnetic perturbation offers the advantage of reducing disruption risks by entraining saturated magnetic islands. These findings are instrumental for optimizing EF correction in fusion devices, thereby enhancing tearing mode suppression and overall plasma stability.

3D magnetic sensors↗
DOE OSTI2020

Current-driven transverse domain wall oscillations in perpendicular spin-valve structures

Spin-transfer-driven oscillations of a transverse domain wall confined to a perpendicular spin-valve structure are investigated using a one-dimensional model. The stack consists of a polarizer, nonmagnetic spacer, soft free layer, and pinned magnetic layer. It is found that the domain-wall oscillation frequency is a nonmonotonic, highly asymmetrical function of applied electric current, showing a strong dependence on the current direction and the relative strengths of the interfacial and bulk spin-transfer torques. Micromagnetic analysis reveals that the surprising and atypical oscillator response is due to an interplay between the interfacial spin-transfer torque, the bulk spin-transfer torque, the exchange torque, and the damping torque. The underlying physical and material responses are examined, including the important role of the domain-wall twist. The competitions between the involved torques under different operating conditions suggest that the oscillator could serve as a model system to investigate magnetic and spintronic phenomena at the nanoscale. Finally, the observed current-dependent twisting of the free-layer magnetization about the axis of precession may further be found interesting for investigations of the interaction between spin-polarized current and chiral spin structures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2000

Postural control model interpretation of stabilogram diffusion analysis

Collins and De Luca [Collins JJ. De Luca CJ (1993) Exp Brain Res 95: 308-318] introduced a new method known as stabilogram diffusion analysis that provides a quantitative statistical measure of the apparently random variations of center-of-pressure (COP) trajectories recorded during quiet upright stance in humans. This analysis generates a stabilogram diffusion function (SDF) that summarizes the mean square COP displacement as a function of the time interval between COP comparisons. SDFs have a characteristic two-part form that suggests the presence of two different control regimes: a short-term open-loop control behavior and a longer-term closed-loop behavior. This paper demonstrates that a very simple closed-loop control model of upright stance can generate realistic SDFs. The model consists of an inverted pendulum body with torque applied at the ankle joint. This torque includes a random disturbance torque and a control torque. The control torque is a function of the deviation (error signal) between the desired upright body position and the actual body position, and is generated in proportion to the error signal, the derivative of the error signal, and the integral of the error signal [i.e. a proportional, integral and derivative (PID) neural controller]. The control torque is applied with a time delay representing conduction, processing, and muscle activation delays. Variations in the PID parameters and the time delay generate variations in SDFs that mimic real experimental SDFs. This model analysis allows one to interpret experimentally observed changes in SDFs in terms of variations in neural controller and time delay parameters rather than in terms of open-loop versus closed-loop behavior.

NASA Program Biomedical Research and Countermeasur↗
DOE OSTI2025

Thermal and electric multidomain dynamic model for integration of power grid distribution with behind-the-meter devices

As renewable energy sources like solar and wind power become more integrated into the grid, coordinated control of behind-the-meter devices is crucial for enhancing grid flexibility and reliability and for meeting cost targets, with standardized models being developed to support this transition. The increasing flexibility and uncertainty of integrated renewable energy grids, along with interactions between various subsystems, make traditional steady-state modeling insufficient to capture transient and dynamic behaviors. Current models (e.g., composite load and battery equivalent models) focus on thermodynamic or electrical characteristics but overlook critical electromechanical interactions. This limits the ability to share performance information for grid services and hampers fast dynamic simulations. In addition, motor stalling is usually triggered by a fault event and attributed to the characteristics of the mechanical torque of the motor, resulting in absorption of a large amount of reactive power during the stalling period. Further, this significant withdrawal of reactive power will deteriorate the dynamic voltage stability of power grids and cause delayed voltage recovery. Therefore, an in-depth modeling of the thermodynamics or mechanical torque is essential to study the impacts of the realistic torque characteristics of those behind-the-meter devices on power system voltage stability. This study developed a dynamic multidomain model for building HVAC systems, such as air-source heat pumps, to simulate their thermal and electrical responses to grid transients. The model can accurately predict power metrics with a mean absolute percentage error of 10%, by validating against with power system computer-aided design performance data. Case studies demonstrate the model capability of capturing the transient response to sudden voltage changes, rapid load fluctuations, and system shutdowns respectively. During a sudden voltage drop (30% for 0.1s), a fully loaded heat pump’s motor speed dropped, continued declining, and shut down after 3.6s, with severe power oscillations and a torque spike. A partially loaded unit experienced temporary oscillations but stabilized. Under higher building loads, compressor speed increased from 64% to 100%, with power and torque rising before stabilizing. In safety-triggered shutdowns, power decreased after minor fluctuations, and torque briefly spiked before dropping to zero.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Advanced Permanent Magnet Generator Topologies Using Multimaterial Shape Optimization and 3D Printing: Preprint

A vast majority of utility-scale wind turbine generators in the United States are dependent on foreign- sourced rare-earth permanent magnets that are vulnerable to supply chain uncertainties. Many small wind original equipment manufacturers are motivated to pursue continuous improvements to the generator design to lower the material and production costs and improve performance by lowering cogging torque and increasing the efficiency. Traditional design and manufacturing offer limited opportunities. In this work, we demonstrate advanced design approaches for a 15-kW baseline wind turbine generator by making use of recent progress in three-dimensional (3D) printing of polymer- bonded magnets and, electrical and structural steel. We explore three methods of magnet parametrization using Bezier curves resulting in symmetric, asymmetric and multimaterial magnet designs. We employ a multiphysics approach combining parametric computer-aided design modeling, finite- element analysis and targeted sampling to identify novel designs with more opportunities for reducing rare-earth material, improving efficiency and minimizing cogging torque. The results show that asymmetric-pole design and multimaterial-pole designs offer a greater opportunity to minimize rare-earth magnet materials by up to 35% with similar performance as the baseline generator, suggesting newer opportunities with design freedom beyond traditional limits of symmetry and as allowed by 3D printing.

Bezier curves↗
DOE OSTI2022

Important powertrain dynamics for developing models for control of connected and automated electrified vehicles

Connected and Automated Vehicles (CAV) technology presents significant opportunities for energy saving in the transportation sector. CAV technology forecasts vehicle and powertrain power needs under various terrain, ambient, and traffic conditions. Integration of the CAV technology in Hybrid Electric Vehicles (HEVs) provides the opportunity for optimal vehicle operation. Indeed, Hybrid Electric Vehicle powertrains present high degrees of flexibility and possibility for choosing optimum powertrain modes based on the predicted traction power needs. In modeling complex CAV powertrain dynamics, the modeler needs to consider short-time scale powertrain dynamics, such as engine transients, and hysteresis of mode-switching for a multi-mode HEV. Therefore, the powertrain dynamics essential for developing powertrain controllers for a class of connected HEVs is presented. To this end, control-oriented powertrain dynamic models for a test vehicle consisting of full electric, hybrid, and conventional engine operating modes are developed. The resulting powertrain model can forecast vehicle traction torque and energy consumption for the specified prediction horizon of the test vehicle. The model considers different operating modes and associated energy penalty terms for mode switching. Thus, the vehicle controller can determine the optimum powertrain mode, torque, and speed for forecasted vehicle operation via utilizing connectivity data. The powertrain model is validated against the experimental data and shows prediction error of less than 5% for predicting vehicle energy consumption. The model is used to create energy penalty maps that can be used for CAV control, for example fuel penalty map for engine torque changes (10–40 Nm) at each engine speed. The results of model-based optimization show optimum switching delays ranging from 0.4 to 1.4 s to avoid hysteresis in mode switching.

Engineering↗
DOE OSTI2022

Optimizing the fuel efficiency of an opposed piston engine for electric power generation

This paper investigates the optimal crankshaft motion for an opposed piston (OP) engine in a novel hybrid architecture to maximize fuel efficiency. The OP engine was selected for this work due to its inherent thermodynamic benefits and the balanced nature of the engine which can achieve downsizing through reducing the number of cylinders rather than the individual cylinder volume. The typical geartrain required on an OP engine was exchanged for two electric motors, reducing friction loss and decoupling the crankshafts. Using the motors to control the crankshaft motion profiles, this architecture introduces capabilities to dynamically vary compression ratio, combustion volume, and scavenging dynamics. To leverage these opportunities, an optimization scheme was developed utilizing nonlinear optimization of a 0-D model to compute the crankshaft motion profile that maximizes the work generated by the system. This optimization was then iteratively coupled with a high fidelity model which supplies the cylinder flow boundary conditions. This iterative approach reduces the model complexity used in the optimal control problem (OCP) while capturing the gas exchange dynamics critical to the 2-stroke cycle of the OP engine. By using the rate of change of motor torque as the input to the OCP, the torque fluctuation in a single cycle can be limited to ensure tracking feasibility. The results show crankshaft velocity slows during the compression stroke and conversely accelerates during the expansion stroke, reducing the peak motor torque required for control and thus reducing the motor losses. The extended residence time at top dead center, however, leads to an increase in heat transfer, illustrating the trade-off between the work extraction efficiency and the indicated engine efficiency.

Engineering↗
NASA NTRS1978

Steering laws for double-gimbal control-moment gyros

Derivation of CMG steering laws was considered. The type of CMG considered was a double gimbal one which has a constant magnitude momentum vector. The equations for the generation of torque of an inertially fixed CMG are presented. To realize a particular commanded torque, the gimbal rates of the CMG must be specified. These gimbal rates are shown to be functions of gimbal angle orientation and the commanded torque. The conventional cross-product steering law was derived by minimizing a function which was the square of the magnitude of the error between the commanded torque and the realizable torque of the gyro. This steering law was severely dependent upon the inner gimbal angle. When the inner gimbal reaches 90 degrees, the cross product law requires infinite outer gimbal rates to realize finite torque values.

Winder, S. W.↗
NASA NTRS1982

A digital optical torquemeter for high rotational speed applications

A digital optical torquemeter system designed for applications at high rotational speeds was fabricated and tested for zero stability at speeds up to 20,000 rpm. Data obtained in a spin rig and with simulated inputs demonstrate that the system is capable of measuring torque bar twist to within 0.03 degrees at speeds of 30,000 rpm. The optical system uses fiber optic bundles to transmit light to the torque bar and to silicon avalanche detectors. The system is microcomputer based and provides measurements of average torque and torque as a function of angular shaft position. The torquemeter requires no bearings or other contact between the rotating torque bar and the nonrotating optics, and tolerates movement of the torque bar as large as 1 mm relative to the optics.

Lesco, D. J.↗
NASA NTRS1988

GOES dynamic propagation of attitude

The spacecraft in the next series of Geostationary Operational Environmental Satellites (GOES-Next) are Earth pointing and have 5-year mission lifetimes. Because gyros can be depended on only for a few years of continuous use, they will be turned off during routine operations. This means attitude must, at times, be determined without benefit of gyros and, often, using only Earth sensor data. To minimize the interruption caused by dumping angular momentum, these spacecraft have been designed to reduce the environmental torque acting on them and incorporate an adjustable solar trim tab for fine adjustment. A new support requirement for GOES-Next is that of setting the solar trim tab. Optimizing its setting requires an estimate of the unbalanced torque on the spacecraft. These two requirements, determining attitude without gyros and estimating the external torque, are addressed by replacing or supplementing the gyro propagation with a dynamic one, that is, one that integrates the rigid body equations of motion. By processing quarter-orbit or longer batches, this approach takes advantage of roll-yaw coupling to observe attitude completely without Sun sensor data. Telemetered momentum wheel speeds are used as observations of the unbalanced external torques. GOES-Next provides a unique opportunity to study dynamic attitude propagation. The geosynchronous altitude and adjustable trim tab minimize the external torque and its uncertainty, making long-term dynamic propagation feasible. This paper presents the equations for dynamic propagation, an analysis of the environmental torques, and an estimate of the accuracies obtainable with the proposed method.

Markley, F. Landis↗
NASA NTRS1995

High performance stepper motors for space mechanisms

Hybrid stepper motors are very well adapted to high performance space mechanisms. They are very simple to operate and are often used for accurate positioning and for smooth rotations. In order to fulfill these requirements, the motor torque, its harmonic content, and the magnetic parasitic torque have to be properly designed. Only finite element computations can provide enough accuracy to determine the toothed structures' magnetic permeance, whose derivative function leads to the torque. It is then possible to design motors with a maximum torque capability or with the most reduced torque harmonic content (less than 3 percent of fundamental). These later motors are dedicated to applications where a microstep or a synchronous mode is selected for minimal dynamic disturbances. In every case, the capability to convert electrical power into torque is much higher than on DC brushless motors.

Sega, Patrick↗
NASA NTRS1995

Dynamic analysis of astronaut motions in microgravity: Applications for Extravehicular Activity (EVA)

Simulations of astronaut motions during extravehicular activity (EVA) tasks were performed using computational multibody dynamics methods. The application of computational dynamic simulation to EVA was prompted by the realization that physical microgravity simulators have inherent limitations: viscosity in neutral buoyancy tanks; friction in air bearing floors; short duration for parabolic aircraft; and inertia and friction in suspension mechanisms. These limitations can mask critical dynamic effects that later cause problems during actual EVA's performed in space. Methods of formulating dynamic equations of motion for multibody systems are discussed with emphasis on Kane's method, which forms the basis of the simulations presented herein. Formulation of the equations of motion for a two degree of freedom arm is presented as an explicit example. The four basic steps in creating the computational simulations were: system description, in which the geometry, mass properties, and interconnection of system bodies are input to the computer; equation formulation based on the system description; inverse kinematics, in which the angles, velocities, and accelerations of joints are calculated for prescribed motion of the endpoint (hand) of the arm; and inverse dynamics, in which joint torques are calculated for a prescribed motion. A graphical animation and data plotting program, EVADS (EVA Dynamics Simulation), was developed and used to analyze the results of the simulations that were performed on a Silicon Graphics Indigo2 computer. EVA tasks involving manipulation of the Spartan 204 free flying astronomy payload, as performed during Space Shuttle mission STS-63 (February 1995), served as the subject for two dynamic simulations. An EVA crewmember was modeled as a seven segment system with an eighth segment representing the massive payload attached to the hand. For both simulations, the initial configuration of the lower body (trunk, upper leg, and lower leg) was a neutral microgravity posture. In the first simulation, the payload was manipulated around a circular trajectory of 0.15 m radius in 10 seconds. It was found that the wrist joint theoretically exceeded its ulnal deviation limit by as much as 49. 8 deg and was required to exert torques as high as 26 N-m to accomplish the task, well in excess of the wrist physiological limit of 12 N-m. The largest torque in the first simulation, 52 N-m, occurred in the ankle joint. To avoid these problems, the second simulation placed the arm in a more comfortable initial position and the radius and speed of the circular trajectory were reduced by half. As a result, the joint angles and torques were reduced to values well within their physiological limits. In particular, the maximum wrist torque for the second simulation was only 3 N-m and the maximum ankle torque was only 6 N-m.

Newman, Dava J.↗
NASA NTRS1996

Solar and Heliospheric Observatory (SOHO) Flight Dynamics Simulations Using MATLAB (R)

This paper describes a study to verify onboard attitude control laws in the coarse Sun-pointing (CSP) mode by simulation and to develop procedures for operational support for the Solar and Heliospheric Observatory (SOHO) mission. SOHO was launched on December 2, 1995, and the predictions of the simulation were verified with the flight data. This study used a commercial off the shelf product MATLAB(tm) to do the following: Develop procedures for computing the parasitic torques for orbital maneuvers; Simulate onboard attitude control of roll, pitch, and yaw during orbital maneuvers; Develop procedures for predicting firing time for both on- and off-modulated thrusters during orbital maneuvers; Investigate the use of feed forward or pre-bias torques to reduce the attitude handoff during orbit maneuvers - in particular, determine how to use the flight data to improve the feed forward torque estimates for use on future maneuvers. The study verified the stability of the attitude control during orbital maneuvers and the proposed use of feed forward torques to compensate for the attitude handoff. Comparison of the simulations with flight data showed: Parasitic torques provided a good estimate of the on- and off-modulation for attitude control; The feed forward torque compensation scheme worked well to reduce attitude handoff during the orbital maneuvers. The work has been extended to prototype calibration of thrusters from observed firing time and observed reaction wheel speed changes.

Headrick, R. D.↗
NASA NTRS1995

Traction Drives for Zero Stick-Slip Robots, and Reaction Free, Momentum Balanced Systems

Two differential (dual input, single output) drives (a roller-gear and a pure roller), and a momentum balanced (single input, dual output) drive (pure roller ) were designed, fabricated, and tested. The differential drives are each rated at 295 rad/sec (2800 rpm) input speed, 450 N-m (4,000 in-lbf) output torque. The momentum balanced drive is rated at 302 rad/sec (2880 rpm) input speed, and dual output torques of 434N-m (3840 in-lbf). The Dual Input Differential Roller-Gear Drive (DC-700) has a planetary roller-gear system with a reduction ratio (one input driving the output with the second input fixed) of 29.23: 1. The Dual Input Differential Roller Drive (DC-500) has a planetary roller system with a reduction ratio of approximately 24:1. Each of the differential drives features dual roller-gear or roller arrangements consisting of a sun, four first row planets, four second row planets, and a ring. The Momentum Balanced (Grounded Ring) Drive (DC-400) has a planetary roller system with a reduction ratio of 24:1 with both outputs counterrotating at equal speed. Its single roller cluster consists of a sun, five first and five second row planets, a roller cage or spider and a ring. Outputs are taken from both the roller cage and the ring which counterrotate. Test results reported for all three drives include angular and torque ripple (linearity and cogging), viscous and Coulomb friction, and forward and reverse power efficiency. Of the two differential drives, the Differential Roller Drive had better linearity and less cogging than did the Differential Roller-Gear Drive, but it had higher friction and lower efficiency (particularly at low power throughput levels). Use of full preloading rather than a variable preload system in the Differential Roller Drive assessed a heavy penalty in part load efficiency. Maximum measured efficiency (ratio of power out to power in) was 95% for the Differential Roller-Gear Drive and 86% for the Differential Roller Drive. The Momentum Balanced (Grounded Ring) Drive performed as expected kinematically. Reduction r-atios to the two counterrotating outputs (design nominal=24:1) were measured to be 23.98:1 and 24.12:1 at zero load.. At 25ONm (2200 in-lbf) output torque the ratio changed 2% due to roller creep. This drive was the smoothest of all three as determined from linearity and cogging tests, and maximum measured efficiency (ratio of power out to power in) was 95%. The disadvantages of full preloading as comvared to variable preload were apparent in this drive as in the Differential Roller Drive. Efficiencies at part load were low, but improved dramatically with increases in torque. These were consistent with friction measurements which indicated losses primarily from Coulomb friction. The initial preload level setting was low so roller slip was encountered at higher torques during testing.

Anderson, William J.↗
NASA NTRS1998

Magellan Aerodynamic Characteristics During the Termination Experiment Including Thruster Plume-Free Stream Interaction

Results are presented on the aerodynamic characteristics of the Magellan spacecraft during the October 1994 Termination Experiment, including the effects of the thruster engine exhaust plumes upon the molecular free stream around the spacecraft and upon the aerodynamics coefficients. As Magellan passed through the Venusian atmosphere, the solar arrays were turned in opposite directions relative to the free stream creating a torque on the spacecraft. The spacecraft control system was programmed to counter the effects of this torque with attitude control engines to maintain an inertially fixed attitude. The orientation and reaction engine telemetry returned from Magellan are used to create a model of the aerodynamic torques. Geometric models of the Magellan spacecraft are analyzed with the aid of both free molecular and Direct Simulation Monte Carlo codes. The simulated aerodynamic torques determined are compared to the measured torques. The Direct Simulation Monte Carlo method is also used to model the attitude engine exhaust plumes, the free stream disturbance caused by these plumes, and the resulting torques acting on the spacecraft compared to no-exhaust plume cases. The effect of the exhaust plumes was found to be sufficiently large that thrust reversal is possible.

Cestero, Francisco J.↗
NASA NTRS2000

Formation of Large Regular Satellites of Giant Planets in an Extended Gaseous Nebula: Satellite Migration And Survival - 2

Using an optically thick inner disk and an extended, optically thin outer disk as described in Mosqueira and Estrada, we compute the torque as a function of position in the subnebula, and show that although the torque exerted on the satellite is generally negative, which leads to inward migration as expected, there are regions of the disk where the torque is positive. For our model these regions of positive torque correspond roughly to the locations of Callisto and Iapetus. Though the outer location of zero torque depends on the (unknown) size of the transition region between the inner and outer disks, the result that Saturn's is found much farther out (at approximately 3r(sub c, sup S) where r(sub c, sup S) is Saturn's centrifugal radius) than Jupiter's (at approximately 2r(sub c, sup J), where r(sub c, sup J) is Jupiter's centrifugal radius) is mostly due to Saturn's less massive outer disk, and larger Hill radius. For a satellite to survive in the disk the timescale of satellite migration must be longer than the timescale for gas dissipation. For large satellites (approximately 1000 km) migration is dominated by the gas torque. We consider the possibility that the feedback reaction of the gas disk caused by the redistribution of gas surface density around satellites with masses larger than the inertial mass causes a large drop in the drift velocity of such objects, thus improving the likelihood that they will be left stranded following gas dissipation. We adapt the inviscid inertial mass criterion to include gas drag, and m-dependent non-local deposition of angular momentum.

Mosqueira, I.↗
NASA NTRS2005

Integrated Power and Attitude Control System Demonstrated With Flywheels G2 and D1

On September 14, 2004, NASA Glenn Research Center's Flywheel Development Team experimentally demonstrated a full-power, high-speed, two-flywheel system, simultaneously regulating a power bus and providing a commanded output torque. Operation- and power-mode transitions were demonstrated up to 2000 W in charge and 1100 W in discharge, while the output torque was simultaneously regulated between plus or minus 0.8 N-m. The G2 and D1 flywheels--magnetically levitated carbon-fiber wheels with permanent magnet motors--were used for the experiment. The units were mounted on an air bearing table in Glenn's High Energy Flywheel Facility. The operational speed range for these tests was between 20,000 and 60,000 rpm. The bus voltage was regulated at 125 V during charge and discharge, and charge-discharge and discharge-charge transitions were demonstrated by changing the amount of power that the power supply provided between 300 and 0 W. In a satellite system, this would be the equivalent of changing the amount of energy that the solar array provides to the spacecraft. In addition to regulating the bus voltage, we simultaneously controlled the net torque produced by the two flywheel modules. Both modules were mounted on an air table that was restrained by a load cell. The load cell measured the force on the table, and the torque produced by the two flywheels on the table could be calculated from that measurement. This method was used to measure the torque produced by the modules, yielding net torques from -0.8 to 0.8 N-m. This was the first Glenn demonstration of the Integrated Power and Attitude Control System (IPACS) at high power levels and speeds.

Jansen, Ralph H.↗
NASA NTRS2007

Analytical Approach Validation for the Spin-Stabilized Satellite Attitude

An analytical approach for spin-stabilized spacecraft attitude prediction is presented for the influence of the residual magnetic torques and the satellite in an elliptical orbit. Assuming a quadripole model for the Earth s magnetic field, an analytical averaging method is applied to obtain the mean residual torque in every orbital period. The orbit mean anomaly is used to compute the average components of residual torque in the spacecraft body frame reference system. The theory is developed for time variations in the orbital elements, giving rise to many curvature integrals. It is observed that the residual magnetic torque does not have component along the spin axis. The inclusion of this torque on the rotational motion differential equations of a spin stabilized spacecraft yields conditions to derive an analytical solution. The solution shows that the residual torque does not affect the spin velocity magnitude, contributing only for the precession and the drift of the spin axis of the spacecraft. The theory developed has been applied to the Brazilian s spin stabilized satellites, which are quite appropriated for verification and comparison of the theory with the data generated and processed by the Satellite Control Center of Brazil National Research Institute. The results show the period that the analytical solution can be used to the attitude propagation, within the dispersion range of the attitude determination system performance of Satellite Control Center of Brazil National Research Institute.

Zanardi, Maria Cecilia F. P. S.↗