Stability of large flexible damped spacecraft modeled as elastic continua

Vibrational stability of a large flexible, structurally damped spacecraft subject to large rigid body rotations is analysed modelling the system as an elastic continuum. Using solution of rigid body attitude motion under torque free conditions and modal analysis, the vibrational equations are reduced to ordinary differential equations with time-varying coefficients. Stability analysis is carried out using Floquet theory and Sonin-Polya theorem. The cases of spinning and non-spinning spacecraft idealized as a flexible beam plate undergoing simple structural vibration are analysed in detail. The critical damping required for stabilization is shown to be a function of the spacecraft's inertia ratio and the level of disturbance.

Statistical analysis of spectral data: A methodology for designing an intelligent monitoring system for the diabetic foot

Early detection of (pre-)signs of ulceration on a diabetic foot is valuable for clinical practice. Hyperspectral imaging is a promising technique for detection and classification of such (pre-)signs. However, the number of the spectral bands should be limited to avoid overfitting, which is critical for pixel classification with hyperspectral image data. The goal was to design a detector/classifier based on spectral imaging (SI) with a small number of optical bandpass filters. The performance and stability of the design were also investigated. The selection of the bandpass filters boils down to a feature selection problem. A dataset was built, containing reflectance spectra of 227 skin spots from 64 patients, measured with a spectrometer. Each skin spot was annotated manually by clinicians as "healthy" or a specific (pre-)sign of ulceration. Statistical analysis on the data set showed the number of required filters is between 3 and 7, depending on additional constraints on the filter set. The stability analysis revealed that shot noise was the most critical factor affecting the classification performance. It indicated that this impact could be avoided in future SI systems with a camera sensor whose saturation level is higher than 106, or by postimage processing.

Stability and transition in the flow of polymer solutions

The linear stability analysis of a plane Couette flow of viscoelastic fluid have been studied with the emphasis on two dimensional disturbances with wave number k similar to Re-1/2, where Re is Reynolds number based on maximum velocity and channel width. We employ three models to represent the dilute polymer solution: the classical Oldroyd-B model, the Oldroyd-B model with artificial diffusivity and the non-homogeneous polymer model. The result of the linear stability analysis is found to be sensitive to the polymer model used. While the plane Couette flow is found to be stable to infinitesimal disturbances for the first two models, the last one exhibits a linear instability.

Analysis of solar radiation pressure induced coupled liberations of a gravity stabilized axisymmetric satellite

This paper presents an analysis of solar radiation pressure induced coupled librations of gravity stabilized cylindrical spacecraft with a special reference to geostationary communication satellites. The Lagrangian approach is used to obtain the corresponding equations of motion. The solar induced torques are assumed to be free of librational angles and are represented by their Fourier expansion. The response and periodic solutions are obtained through linear and nonlinear analyses, using the method of harmonic balance in the latter case. The stability conditions are obtained using Routh-Hurwitz criteria. To establish the ranges of validity the analytic response is compared with the numerical solution. Finally, values of the system parameters are suggested to make the satellite behave as desired. Among these is a possible approach to subdue the solar induced roll resonance. It is felt that the approximate analysis presented here should significantly reduce the computational efforts involved in the design and stability analysis of the systems.

Steady-State Stability of Current Mode Active-Clamp ZVS DC-DC Converters

Active-clamp dc-dc converters are pulsewidth-modulated converters having two switches featuring zero-voltage switching at frequencies beyond 100 kHz. Generalized equivalent circuits valid for steady-state and dynamic performance have been proposed for the family of active-clamp converters. The active-clamp converter is analyzed for its dynamic behavior under current control in this paper. The steady-state stability analysis is presented. On account of the lossless damping inherent in the active-clamp converters, it appears that the stability region in the current-controlled active-clamp converters get extended for duty ratios, a little greater than 0.5, unlike in conventional hard-switched converters. The conventional graphical approach fails to assess the stability of current-controlled active-clamp converters due to the coupling between the filter inductor current and resonant inductor current. An analysis that takes into account the presence of the resonant elements is presented to establish the condition for stability. This method correctly predicts the stability of the current-controlled active-clamp converters. A simple expression for the maximum duty cycle for subharmonic free operation is obtained. The results are verified experimentally.

A Comparative Study of Liapunov Stability Analyses of Flexible Damped Satellites

A literal Liapunov stability analysis of a spacecraft with flexible appendages often requires a division of the associated dynamic potential into as many dependent parts as the number of appendages. First part of this paper exposes the stringency in the stability criteria introduced by such a division and shows it to be removable by a “reunion policy.” The policy enjoins the analyst to piece together the sets of criteria for each part. Employing reunion the paper then compares four methods of the Liapunov stability analysis of hybrid dynamical systems illustrated by an inertially coupled, damped, gravity stabilized, elastic spacecraft with four gravity booms having tip masses and a damper rod, all skewed to the orbital plane. The four methods are the method of test density function, assumed modes, and two and one-integral coordinates. Superiority of one-integral coordinate approach is established here. The design plots demonstrate how elastic effects delimit the satellite boom length.

Direct stability evaluation of power systems with detailed generator models using structure-preserving energy functions

Energy-based direct methods for transient stability analysis are potentially useful both as offline tools for planning purposes as well as for online security assessment. In this paper, a novel structure-preserving energy function (SPEF) is developed using the philosophy of structure-preserving model for the system and detailed generator model including flux decay, transient saliency, automatic voltage regulator (AVR), exciter and damper winding. A simpler and yet general expression for the SPEF is also derived which can simplify the computation of the energy function. The system equations and the energy function are derived using the centre-of-inertia (COI) formulation and the system loads are modelled as arbitrary functions of the respective bus voltages. Application of the proposed SPEF to transient stability evaluation of power systems is illustrated with numerical examples.

Stability of the viscous flow of a fluid through a flexible tube

The stability of Hagen-Poiseuille flow of a Newtonian fluid of viscosity eta in a tube of radius R surrounded by a viscoelastic medium of elasticity G and viscosity eta(s) occupying the annulus R < r < HR is determined using a linear stability analysis. The inertia of the fluid and the medium are neglected, and the mass and momentum conservation equations for the fluid and wall are linear. The only coupling between the mean flow and fluctuations enters via an additional term in the boundary condition for the tangential velocity at the interface, due to the discontinuity in the strain rate in the mean flow at the surface. This additional term is responsible for destabilizing the surface when the mean velocity increases beyond a transition value, and the physical mechanism driving the instability is the transfer of energy from the mean flow to the fluctuations due to the work done by the mean flow at the interface. The transition velocity Gamma(t) for the presence of surface instabilities depends on the wavenumber k and three dimensionless parameters: the ratio of the solid and fluid viscosities eta(r) = (eta(s)/eta), the capillary number Lambda = (T/GR) and the ratio of radii H, where T is the surface tension of the interface. For eta(r) = 0 and Lambda = 0, the transition velocity Gamma(t) diverges in the limits k much less than 1 and k much greater than 1, and has a minimum for finite k. The qualitative behaviour of the transition velocity is the same for Lambda > 0 and eta(r) = 0, though there is an increase in Gamma(t) in the limit k much greater than 1. When the viscosity of the surface is non-zero (eta(r) > 0), however, there is a qualitative change in the Gamma(t) vs. k curves. For eta(r) < 1, the transition velocity Gamma(t) is finite only when k is greater than a minimum value k(min), while perturbations with wavenumber k < k(min) are stable even for Gamma--> infinity. For eta(r) > 1, Gamma(t) is finite only for k(min) < k < k(max), while perturbations with wavenumber k < k(min) or k > k(max) are stable in the limit Gamma--> infinity. As H decreases or eta(r) increases, the difference k(max)- k(min) decreases. At minimum value H = H-min, which is a function of eta(r), the difference k(max)-k(min) = 0, and for H < H-min, perturbations of all wavenumbers are stable even in the limit Gamma--> infinity. The calculations indicate that H-min shows a strong divergence proportional to exp (0.0832 eta(r)(2)) for eta(r) much greater than 1.

Stability of non-parabolic flow in a flexible tube

Flows with velocity profiles very different from the parabolic velocity profile can occur in the entrance region of a tube as well as in tubes with converging/diverging cross-sections. In this paper, asymptotic and numerical studies are undertaken to analyse the temporal stability of such 'non-parabolic' flows in a flexible tube in the limit of high Reynolds numbers. Two specific cases are considered: (i) developing flow in a flexible tube; (ii) flow in a slightly converging flexible tube. Though the mean velocity profile contains both axial and radial components, the flow is assumed to be locally parallel in the stability analysis. The fluid is Newtonian and incompressible, while the flexible wall is modelled as a viscoelastic solid. A high Reynolds number asymptotic analysis shows that the non-parabolic velocity profiles can become unstable in the inviscid limit. This inviscid instability is qualitatively different from that observed in previous studies on the stability of parabolic flow in a flexible tube, and from the instability of developing flow in a rigid tube. The results of the asymptotic analysis are extended numerically to the moderate Reynolds number regime. The numerical results reveal that the developing flow could be unstable at much lower Reynolds numbers than the parabolic flow, and hence this instability can be important in destabilizing the fluid flow through flexible tubes at moderate and high Reynolds number. For flow in a slightly converging tube, even small deviations from the parabolic profile are found to be sufficient for the present instability mechanism to be operative. The dominant non-parallel effects are incorporated using an asymptotic analysis, and this indicates that non-parallel effects do not significantly affect the neutral stability curves. The viscosity of the wall medium is found to have a stabilizing effect on this instability.

The effect of boundaries on the plane Couette flow of granular materials: a bifurcation analysis

The tendency of granular materials in rapid shear flow to form non-uniform structures is well documented in the literature. Through a linear stability analysis of the solution of continuum equations for rapid shear flow of a uniform granular material, performed by Savage (1992) and others subsequently, it has been shown that an infinite plane shearing motion may be unstable in the Lyapunov sense, provided the mean volume fraction of particles is above a critical value. This instability leads to the formation of alternating layers of high and low particle concentrations oriented parallel to the plane of shear. Computer simulations, on the other hand, reveal that non-uniform structures are possible even when the mean volume fraction of particles is small. In the present study, we have examined the structure of fully developed layered solutions, by making use of numerical continuation techniques and bifurcation theory. It is shown that the continuum equations do predict the existence of layered solutions of high amplitude even when the uniform state is linearly stable. An analysis of the effect of bounding walls on the bifurcation structure reveals that the nature of the wall boundary conditions plays a pivotal role in selecting that branch of non-uniform solutions which emerges as the primary branch. This demonstrates unequivocally that the results on the stability of bounded shear how of granular materials presented previously by Wang et al. (1996) are, in general, based on erroneous base states.

Stability of fluid flow past a membrane

The stability of fluid flow past a membrane of infinitesimal thickness is analysed in the limit of zero Reynolds number using linear and weakly nonlinear analyses. The system consists of two Newtonian fluids of thickness R* and H R*, separated by an infinitesimally thick membrane, which is flat in the unperturbed state. The dynamics of the membrane is described by its normal displacement from the flat state, as well as a surface displacement field which provides the displacement of material points from their steady-state positions due to the tangential stress exerted by the fluid flow. The surface stress in the membrane (force per unit length) contains an elastic component proportional to the strain along the surface of the membrane, and a viscous component proportional to the strain rate. The linear analysis reveals that the fluctuations become unstable in the long-wave (alpha --> 0) limit when the non-dimensional strain rate in the fluid exceeds a critical value Lambda(t), and this critical value increases proportional to alpha(2) in this limit. Here, alpha is the dimensionless wavenumber of the perturbations scaled by the inverse of the fluid thickness R*(-1), and the dimensionless strain rate is given by Lambda(t) = ((gamma) over dot* R*eta*/Gamma*), where eta* is the fluid viscosity, Gamma* is the tension of the membrane and (gamma) over dot* is the strain rate in the fluid. The weakly nonlinear stability analysis shows that perturbations are supercritically stable in the alpha --> 0 limit.

The effect of boundaries on the plane Couette flow of granular materials: a bifurcation analysis

The tendency of granular materials in rapid shear ow to form non-uniform structures is well documented in the literature. Through a linear stability analysis of the solution of continuum equations for rapid shear flow of a uniform granular material, performed by Savage (1992) and others subsequently, it has been shown that an infinite plane shearing motion may be unstable in the Lyapunov sense, provided the mean volume fraction of particles is above a critical value. This instability leads to the formation of alternating layers of high and low particle concentrations oriented parallel to the plane of shear. Computer simulations, on the other hand, reveal that non-uniform structures are possible even when the mean volume fraction of particles is small. In the present study, we have examined the structure of fully developed layered solutions, by making use of numerical continuation techniques and bifurcation theory. It is shown that the continuum equations do predict the existence of layered solutions of high amplitude even when the uniform state is linearly stable. An analysis of the effect of bounding walls on the bifurcation structure reveals that the nature of the wall boundary conditions plays a pivotal role in selecting that branch of non-uniform solutions which emerges as the primary branch. This demonstrates unequivocally that the results on the stability of bounded shear flow of granular materials presented previously by Wang et al. (1996) are, in general, based on erroneous base states.

Steady-State Stability of Current-Mode Active-Clamp ZVS DC-DC Converters

Active-clamp dc-dc converters are pulsewidth-modulated converters having two switches featuring zero-voltage switching at frequencies beyond 100 kHz. Generalized equivalent circuits valid for steady-state and dynamic performance have been proposed for the family of active-clamp converters. The active-clamp converter is analyzed for its dynamic behavior under current control in this paper. The steady-state stability analysis is presented. On account of the lossless damping inherent in the active-clamp converters, it appears that the stability region in the current-controlled active-clamp converters get extended for duty ratios, a little greater than 0.5 unlike in conventional hard-switched converters. The conventional graphical approach fails to assess the stability of current-controlled active-clamp converters, due to the coupling between the filter inductor current and resonant inductor current. An analysis that takes into account the presence of the resonant elements is presented to establish the condition for stability. This method correctly predicts the stability of the current-controlled active-clamp converters. A simple expression for the maximum duty cycle for subharmonic-free operation is obtained. The results are verified experimentally.

Effect of strength - Matric suction variations on stability of a natural slope

Stability analysis of residual soil slopes are now increasingly being performed with the incorporation of the matric suction component of strength. The matric suction (u(a)-u(w)) component of shear strength is known as apparent cohesion. The relation between matric suction and apparent cohesion (c(app)) may be linear or non-linear. The impact of type of apparent strength versus matric suction relationship on the stability of an unsaturated residual soil slope is examined in this study. Results of the study showed that the factor of safety values were unaffected by the nature of the strength versus matric suction relationship for the residual soil slope examined. This was so as contribution from the effective stress- strength component to the factor of safety predominated over the contribution made by the apparent strength component.

Torque Angle Loop Analysis of Synchronous Machines

This paper reports the dynamic stability analysis of a single machine infinite bus system through torque angle loop analysis and forms an extension of the work on Block diagrams and torque angle loop analysis of synchronous machines reported by I. Nagy [3]. It aims to incorporate in the machine model, the damper windings (one on each axis) and to compare the dynamic behaviour of the system with and without damper windings. The effect of using different stabilizing signals (viz. active power and speed deviations) on the dynamic performance is analysed and the significant effect of damper windings on the dynamic behaviour of the system is highlighted.