Multi-rotor with suspended load: System dynamics and control toolbox

There is an increasing demand for Unmanned Aerial Systems (UAS) to carry suspended loads as this can provide significant benefits to several applications in agriculture, law enforcement and construction. The load impact on the underlying system dynamics should not be neglected as significant feedback forces may be induced on the vehicle during certain flight manoeuvres. The constant variation in operating point induced by the slung load also causes conventional controllers to demand increased control effort. Much research has focused on standard multi-rotor position and attitude control with and without a slung load. However, predictive control schemes, such as Nonlinear Model Predictive Control (NMPC), have not yet been fully explored. To this end, we present a novel controller for safe and precise operation of multi-rotors with heavy slung load in three dimensions. The paper describes a System Dynamics and Control Simulation Toolbox for use with MATLAB/SIMULINK which includes a detailed simulation of the multi-rotor and slung load as well as a predictive controller to manage the nonlinear dynamics whilst accounting for system constraints. It is demonstrated that the controller simultaneously tracks specified waypoints and actively damps large slung load oscillations. A linear-quadratic regulator (LQR) is derived and control performance is compared. Results show the improved performance of the predictive controller for a larger flight envelope, including aggressive manoeuvres and large slung load displacements. The computational cost remains relatively small, amenable to practical implementations.

Improving PLDA speaker verification using WMFD and linear-weighted approaches in limited microphone data conditions

This paper proposes the addition of a weighted median Fisher discriminator (WMFD) projection prior to length-normalised Gaussian probabilistic linear discriminant analysis (GPLDA) modelling in order to compensate the additional session variation. In limited microphone data conditions, a linear-weighted approach is introduced to increase the influence of microphone speech dataset. The linear-weighted WMFD-projected GPLDA system shows improvements in EER and DCF values over the pooled LDA- and WMFD-projected GPLDA systems in inter-view-interview condition as WMFD projection extracts more speaker discriminant information with limited number of sessions/ speaker data, and linear-weighted GPLDA approach estimates reliable model parameters with limited microphone data.

A hybrid of modified PSO and local search on a multi-robot search system

Particle swarm optimization (PSO), a new population based algorithm, has recently been used on multi-robot systems. Although this algorithm is applied to solve many optimization problems as well as multi-robot systems, it has some drawbacks when it is applied on multi-robot search systems to find a target in a search space containing big static obstacles. One of these defects is premature convergence. This means that one of the properties of basic PSO is that when particles are spread in a search space, as time increases they tend to converge in a small area. This shortcoming is also evident on a multi-robot search system, particularly when there are big static obstacles in the search space that prevent the robots from finding the target easily; therefore, as time increases, based on this property they converge to a small area that may not contain the target and become entrapped in that area.Another shortcoming is that basic PSO cannot guarantee the global convergence of the algorithm. In other words, initially particles explore different areas, but in some cases they are not good at exploiting promising areas, which will increase the search time.This study proposes a method based on the particle swarm optimization (PSO) technique on a multi-robot system to find a target in a search space containing big static obstacles. This method is not only able to overcome the premature convergence problem but also establishes an efficient balance between exploration and exploitation and guarantees global convergence, reducing the search time by combining with a local search method, such as A-star.To validate the effectiveness and usefulness of algorithms,a simulation environment has been developed for conducting simulation-based experiments in different scenarios and for reporting experimental results. These experimental results have demonstrated that the proposed method is able to overcome the premature convergence problem and guarantee global convergence.

New graph-based algorithms to efficiently solve large scale open pit mining optimisation problems

In the mining optimisation literature, most researchers focused on two strategic-level and tactical-level open-pit mine optimisation problems, which are respectively termed ultimate pit limit (UPIT) or constrained pit limit (CPIT). However, many researchers indicate that the substantial numbers of variables and constraints in real-world instances (e.g., with 50-1000 thousand blocks) make the CPIT’s mixed integer programming (MIP) model intractable for use. Thus, it becomes a considerable challenge to solve the large scale CPIT instances without relying on exact MIP optimiser as well as the complicated MIP relaxation/decomposition methods. To take this challenge, two new graph-based algorithms based on network flow graph and conjunctive graph theory are developed by taking advantage of problem properties. The performance of our proposed algorithms is validated by testing recent large scale benchmark UPIT and CPIT instances’ datasets of MineLib in 2013. In comparison to best known results from MineLib, it is shown that the proposed algorithms outperform other CPIT solution approaches existing in the literature. The proposed graph-based algorithms leads to a more competent mine scheduling optimisation expert system because the third-party MIP optimiser is no longer indispensable and random neighbourhood search is not necessary.

Quantification of stability in an agility drill using linear and nonlinear measures of variability

This study implemented linear and nonlinear methods of measuring variability to determine differences in stability of two groups of skilled (n = 10) and unskilled (n = 10) participants performing 3m forward/backward shuttle agility drill. We also determined whether stability measures differed between the forward and backward segments of the drill. Finally, we sought to investigate whether local dynamic stability, measured using largest finite-time Lyapunov exponents, changed from distal to proximal lower extremity segments. Three-dimensional coordinates of five lower extremity markers data were recorded. Results revealed that the Lyapunov exponents were lower (P < 0.05) for skilled participants at all joint markers indicative of higher levels of local dynamic stability. Additionally, stability of motion did not differ between forward and backward segments of the drill (P > 0.05), signifying that almost the same control strategy was used in forward and backward directions by all participants, regardless of skill level. Furthermore, local dynamic stability increased from distal to proximal joints (P < 0.05) indicating that stability of proximal segments are prioritized by the neuromuscular control system. Finally, skilled participants displayed greater foot placement standard deviation values (P < 0.05), indicative of adaptation to task constraints. The results of this study provide new methods for sport scientists, coaches to characterize stability in agility drill performance.

Study of robust H(infinity) filtering application in loosely coupled INS/GPS system

Since a celebrate linear minimum mean square (MMS) Kalman filter in integration GPS/INS system cannot guarantee the robustness performance, a H(infinity) filtering with respect to polytopic uncertainty is designed. The purpose of this paper is to give an illustration of this application and a contrast with traditional Kalman filter. A game theory H(infinity) filter is first reviewed; next we utilize linear matrix inequalities (LMI) approach to design the robust H(infinity) filter. For the special INS/GPS model, unstable model case is considered. We give an explanation for Kalman filter divergence under uncertain dynamic system and simultaneously investigate the relationship between H(infinity) filter and Kalman filter. A loosely coupled INS/GPS simulation system is given here to verify this application. Result shows that the robust H(infinity) filter has a better performance when system suffers uncertainty; also it is more robust compared to the conventional Kalman filter.

Finite Element Solution Of Surface-Tension Driven Flows In Laser Surface-Melting

Lasers are very efficient in heating localized regions and hence they find a wide application in surface treatment processes. The surface of a material can be selectively modified to give superior wear and corrosion resistance. In laser surface-melting and welding problems, the high temperature gradient prevailing in the free surface induces a surface-tension gradient which is the dominant driving force for convection (known as thermo-capillary or Marangoni convection). It has been reported that the surface-tension driven convection plays a dominant role in determining the melt pool shape. In most of the earlier works on laser-melting and related problems, the finite difference method (FDM) has been used to solve the Navier Stokes equations [1]. Since the Reynolds number is quite high in these cases, upwinding has been used. Though upwinding gives physically realistic solutions even on a coarse grid, the results are inaccurate. McLay and Carey have solved the thermo-capillary flow in welding problems by an implicit finite element method [2]. They used the conventional Galerkin finite element method (FEM) which requires that the pressure be interpolated by one order lower than velocity (mixed interpolation). This restricts the choice of elements to certain higher order elements which need numerical integration for evaluation of element matrices. The implicit algorithm yields a system of nonlinear, unsymmetric equations which are not positive definite. Computations would be possible only with large mainframe computers.Sluzalec [3] has modeled the pulsed laser-melting problem by an explicit method (FEM). He has used the six-node triangular element with mixed interpolation. Since he has considered the buoyancy induced flow only, the velocity values are small. In the present work, an equal order explicit FEM is used to compute the thermo-capillary flow in the laser surface-melting problem. As this method permits equal order interpolation, there is no restriction in the choice of elements. Even linear elements such as the three-node triangular elements can be used. As the governing equations are solved in a sequential manner, the computer memory requirement is less. The finite element formulation is discussed in this paper along with typical numerical results.

Spin-state equilibria in the LCoO3 system from 59Co NMR

Spin-state equilibria in the whole set of LCoO3 (where L stands for a rare-earth metal or Y) have been investigated with the use of 59Co NMR as a probe for the polycrystalline samples (except Ce) in the temperature interval 110-550 K and frequency range 3- 11.6 MHz. Besides confirming the coexistence of the high-spin—low-spin state in this temperature range, a quadrupolar interaction of ∼0.1 -0.5 MHz has been detected for the first time from 59Co NMR. The NMR line shape is found to depend strongly on the relative magnitude of the magnetic and quadrupolar interactions present. Analysis of the powder pattern reveals two basically different types of transferred hyperfine interaction between the lighter and heavier members of the rare-earth series. The first three members of the lighter rare-earth metals La, Pr (rhombohedral), and Nd (tetragonal), exhibit second-order quadrupolar interaction with a zero-asymmetry parameter at lower temperatures. Above a critical temperature TS (dependent on the size of the rare-earth ion), the quadrupolar interaction becomes temperature dependent and eventually gives rise to a first-order interaction thus indicating a possible second-order phase change. Sm and Eu (orthorhombic) exhibit also a second-order quadrupolar interaction with a nonzero asymmetry parameter ((η∼0.47)) at 300 K, while the orthorhombic second-half members (Dy,..., Lu and Y) exhibit first-order quadrupolar interaction at all temperatures. Normal paramagnetic behavior, i.e., a linear variation of Kiso with T-1, has been observed in the heavier rare-earth cobaltites (Er,..., Lu and Y), whereas an anomalous variation has been observed in (La,..., Nd)CoO3. Thus, Kiso increases with increasing temperature in PrCoO3 and NdCoO3. These observations corroborate the model of the spin-state equilibria in LCoO3 originally proposed by Raccah and Goodenough. A high-spin—low-spin ratio, r=1, can be stabilized in the perovskite structure by a cooperative displacement of the oxygen atoms from the high-spin towards the low-spin cation. Where this ordering into high- and low-spin sublattices occurs at r=1, one can anticipate equivalent displacement of all near-neighbor oxygen atoms towards a low-spin cobalt ion. Thus the heavier LCoO3 exhibits a small temperature-independent first-order quadrupolar interaction. Where r<1, the high- and low-spin states are disordered, giving rise to a temperature-dependent second-order quadrupolar interaction with an anomalous Kiso for the lighter LCoO3.

On the direct parallel solution of systems of linear equations: New algorithms and systolic structures

The paper presents two new algorithms for the direct parallel solution of systems of linear equations. The algorithms employ a novel recursive doubling technique to obtain solutions to an nth-order system in n steps with no more than 2n(n −1) processors. Comparing their performance with the Gaussian elimination algorithm (GE), we show that they are almost 100% faster than the latter. This speedup is achieved by dispensing with all the computation involved in the back-substitution phase of GE. It is also shown that the new algorithms exhibit error characteristics which are superior to GE. An n(n + 1) systolic array structure is proposed for the implementation of the new algorithms. We show that complete solutions can be obtained, through these single-phase solution methods, in 5n−log2n−4 computational steps, without the need for intermediate I/O operations.

Asymptotic behavior of a hierarchical system of learning automata

Learning automata arranged in a two-level hierarchy are considered. The automata operate in a stationary random environment and update their action probabilities according to the linear-reward- -penalty algorithm at each level. Unlike some hierarchical systems previously proposed, no information transfer exists from one level to another, and yet the hierarchy possesses good convergence properties. Using weak-convergence concepts it is shown that for large time and small values of parameters in the algorithm, the evolution of the optimal path probability can be represented by a diffusion whose parameters can be computed explicitly.

Analysis of an internally mountable accelerometer balance system for use with non-isotropic models in shock tunnels

Knowledge of drag force is an important design parameter in aerodynamics. Measurement of aerodynamic forces at hypersonic speed is a challenge and usually ground test facilities like shock tunnels are used to carry out such tests. Accelerometer based force balances are commonly employed for measuring aerodynamic drag around bodies in hypersonic shock tunnels. In this study, we present an analysis of the effect of model material on the performance of an accelerometer balance used for measurement of drag in impulse facilities. From the experimental studies performed on models constructed out of Bakelite HYLEM and Aluminum, it is clear that the rigid body assumption does not hold good during the short testing duration available in shock tunnels. This is notwithstanding the fact that the rubber bush used for supporting the model allows unconstrained motion of the model during the short testing time available in the shock tunnel. The vibrations induced in the model on impact loading in the shock tunnel are damped out in metallic model, resulting in a smooth acceleration signal, while the signal become noisy and non-linear when we use non-isotropic materials like Bakelite HYLEM. This also implies that careful analysis and proper data reduction methodologies are necessary for measuring aerodynamic drag for non-metallic models in shock tunnels. The results from the drag measurements carried out using a 60 degrees half angle blunt cone is given in the present analysis.

Evaluation of a remote drafting system for regulating sheep access to supplement

Remote drafting technology now available for sheep makes possible targeted supplementation of individuals within a grazing flock. This system was evaluated by using 68 Merino wethers grazing dry-season, native Mitchell grass pasture (predominantly Astrebla spp.) as a group and receiving access to lupin grain through a remote drafter 0, 1, 2, 4 or 7 days/week for 8 weeks. The sole paddock watering point was separately fenced and access was via a one-way flow gate. Sheep exited the watering point through a remote drafter operated by solar power and were drafted by radio frequency identification (RFID) tag, according to treatment, either back into the paddock or into a common supplement yard where lupins were provided ad libitum in a self-feeder. Sheep were drafted into the supplement yard on only their first time through the drafter during the prescribed 24-h period and exited the supplement yard via one-way flow gates in their own time. The remote drafter operated with a high accuracy, with only 2.1% incorrect drafts recorded during the experimental period out of a total of 7027 sheep passes through the remote drafter. The actual number of accesses to supplement for each treatment group, in order, were generally less than that intended, i.e. 0.02, 0.69, 1.98, 3.35 and 6.04 days/week. Deviations from the intended number of accesses to supplement were mainly due to sheep not coming through to water on their allocated day of treatment access, although some instances were due to incorrect drafts. There was a non-linear response in growth rate to increased frequency of access to lupins with the growth rate response plateauing at similar to 3 actual accesses per week, corresponding to a growth rate of 72.5 g/head. day. This experiment has demonstrated the application of the remote drafting supplementation system for the first time under grazing conditions and with the drafter operated completely from solar power. The experiment demonstrates a growth response to increasing frequency of access to supplement and provides a starting point with which to begin to develop feeding strategies to achieve sheep weight-change targets.

A conditionally linearized Monte Carlo filter in non-linear structural dynamics

State and parameter estimations of non-linear dynamical systems, based on incomplete and noisy measurements, are considered using Monte Carlo simulations. Given the measurements. the proposed method obtains the marginalized posterior distribution of an appropriately chosen (ideally small) subset of the state vector using a particle filter. Samples (particles) of the marginalized states are then used to construct a family of conditionally linearized system of equations and thus obtain the posterior distribution of the states using a bank of Kalman filters. Discrete process equations for the marginalized states are derived through truncated Ito-Taylor expansions. Increased analyticity and reduced dispersion of weights computed over a smaller sample space of marginalized states are the key features of the filter that help achieve smaller sample variance of the estimates. Numerical illustrations are provided for state/parameter estimations of a Duffing oscillator and a 3-DOF non-linear oscillator. Performance of the filter in parameter estimation is also assessed using measurements obtained through experiments on simple models in the laboratory. Despite an added computational cost, the results verify that the proposed filter generally produces estimates with lower sample variance over the standard sequential importance sampling (SIS) filter.

Radiation induced depinning of a charge density wave system

The frequency-dependent response of a pinned charge density wave is considered in terms of forced vibration of an oscillator held in an anharmonic well. It is shown that the effective pinning-frequency can be reduced by applying a d.c. field. If a strong a.c. field, superposed on a d.c. field is applied on such a system “jumps” can be observed in the frequency dependent response of the system. The conditions at which these “jumps” occur are investigated with reference to NbSe3. The possibility of observing such phenomena in other systems like superionic conductors, non-linear dielectrics like ferroelectrics is pointed out. The characteristics are expressed in terms of some “scaled variables” — in terms of which the characteristics show a universal behaviour.

Approximate analysis of coupled non-linear non-conservative systems subjected to step-function excitation

The paper deals with the approximate analysis of non-linear non-conservative systems oftwo degrees of freedom subjected to step-function excitation. The method of averaging of Krylov and Bogoliubov is used to arrive at the approximate equations for amplitude and phase. An example of a spring-mass-damper system is presented to illustrate the method and a comparison with numerical results brings out the validity of the approach.