Modeling and waveform optimization of stick–slip micro-drives using the method of dimensionality reduction

In this paper, the dynamics of piezo-actuated stick–slip micro-drives are studied experimentally and theoretically. First, the stick–slip-based force-generating test stand is introduced, and experimental results are presented. Then, a numerical model is formulated which explicitly includes the dynamics of normal and tangential properties of the contact areas in the frictional driving elements of the drive. The contact forces are simulated using the method of dimensionality reduction. We show that the experimentally observed behavior can be described without using any fitting parameters or assuming any generalized laws of friction if the explicit contact mechanics of the frictional contacts is taken into account. Furthermore, an even simpler model of the drive is developed to get a qualitative understanding of the system. It is employed to gain a new actuation method, which reduces the vibrations of the drive’s runner and therefore enhances its performance.     

Understanding root cause of stick–slip vibrations in deep drilling with drag bits

In search for the root cause of stick–slip, a mode of torsional vibrations of a drilling assembly, a linear stability analysis of coupled axial–torsional vibrations has been carried out. It has been shown that in a rotary drilling system with axial and torsional degree of freedom two distinct modes of self-excited vibrations are present: axial and torsional. These axial (torsional) modes of vibrations are due to resonance between the cutting forces acting at the bit and the axial (torsional) natural modes of drillstring vibrations. It has been demonstrated that although axial and torsional modes of vibrations do affect each other the underlying mechanisms driving these modes of vibrations are completely different. In particular, the only driving mechanism of the axial vibrations is the regenerative effect, while there are two distinct mechanisms that drive the torsional vibrations: (i) the cutting action of the bit, and (ii) the wearflat/rock interaction. Moreover, in the case of the torsional vibrations the regenerative effect plays only a secondary role. The results of the present study indicate that the axial compliance can play a stabilizing role. In particular, the stabilizing role of the axial compliance increases as the ratio of the torsional to the axial natural frequency of the drillstring vibrations decreases.     

An adaptive finite element method for viscoplastic flows in a square pipe with stick–slip at the wall

This paper presents the numerical solution of non-linear yield stress phenomena by using a new mixed anisotropic auto-adaptive finite element method. The Poiseuille flow of a Bingham fluid with slip yield boundary condition at the wall is considered. Despite its practical interest, for instance for pipeline flows of yield-stress fluids such as concrete and cements, this problem has not been addressed yet to our knowledge. The case of a pipe with a square section has been investigated in detail. The computations cover the full range of the two main dimensionless numbers and exhibit complex flow patterns: all the different flow regimes are completely identified.     

Transition from stick to slip in Hertzian contact with “Griffith” friction: The Cattaneo–Mindlin problem revisited

Classically, the transition from stick to slip is modelled with Amonton–Coulomb law, leading to the Cattaneo–Mindlin problem, which is amenable to quite general solutions using the idea of superposing normal contact pressure distributions – in particular superposing the full sliding component of shear with a corrective distribution in the stick region. However, faults model in geophysics and recent high-speed measurements of the real contact area and the strain fields in dry (nominally flat) rough interfaces at macroscopic but laboratory scale, all suggest that the transition from ‘static’ to ‘dynamic’ friction can be described, rather than by Coulomb law, by classical fracture mechanics singular solutions of shear cracks. Here, we introduce an ‘adhesive’ model for friction in a Hertzian spherical contact, maintaining the Hertzian solution for the normal pressures, but where the inception of slip is given by a Griffith condition. In the slip region, the standard Coulomb law continues to hold. This leads to a very simple solution for the Cattaneo–Mindlin problem, in which the “corrective” solution in the stick area is in fact similar to the mode II equivalent of a JKR singular solution for adhesive contact. The model departs from the standard Cattaneo–Mindlin solution, showing an increased size of the stick zone relative to the contact area, and a sudden transition to slip when the stick region reaches a critical size (the equivalent of the pull-off contact size of the JKR solution). The apparent static friction coefficient before sliding can be much higher than the sliding friction coefficient and, for a given friction fracture “energy”, the process results in size and normal load dependence of the apparent static friction coefficient. Some qualitative agreement with Fineberg's group experiments for friction exists, namely the stick–slip boundary quasi-static prediction may correspond to the arrest of their slip “precursors”, and the rapid collapse to global sliding when the precursors arrest front has reached about half the interface may correspond to the reach of the “critical” size for the stick zone.     

On the Hamiltonian Hopf Bifurcations in the 3D Hénon–Heiles Family

An axially symmetric perturbed isotropic harmonic oscillator undergoes several bifurcations as the parameter adjusting the relative strength of the two terms in the cubic potential is varied. We show that three of these bifurcations are Hamiltonian Hopf bifurcations. To this end we analyse an appropriately chosen normal form. It turns out that the linear behaviour is not that of a typical Hamiltonian Hopf bifurcation as the eigen-values completely vanish at the bifurcation. However, the nonlinear structure is that of a Hamiltonian Hopf bifurcation. The result is obtained by formulating geometric criteria involving the normalized Hamiltonian and the reduced phase space.     

Distribution of stable islands within chaotic areas in the non-hyperbolic and hyperbolic regimes in the Hénon–Heiles system

Nonlinear Dynamics - We provide rigorous computer-assisted proofs of the existence of different dynamical objects, like stable families of periodic orbits, bifurcations and stable invariant tori...     

Application of a modified Lindstedt–Poincaré method in coupled TDOF systems with quadratic nonlinearity and a constant external excitation

An analytical approach is developed for the nonlinear oscillation of a conservative, two-degree-of-freedom (TDOF) mass-spring system with serial combined linear–nonlinear stiffness excited by a constant external force. The main idea of the proposed approach lies in two categories, the first one is the transformation of two nonlinear differential equations of a two-mass system using suitable intermediate variables into a single nonlinear differential equation. Another is the treatment a quadratic nonlinear oscillator (QNO) by the modified Lindstedt–Poincaré (L-P) method presented recently by the authors. The first-order and second-order analytical approximations for the modified L-P method are established for the QNOs with satisfactory results. After solving the nonlinear differential equation, the displacements of two-mass system can be obtained directly from the governing linear second-order differential equation. Unlike the common perturbation method, the modified L-P method is valid for weak as well as strong nonlinear oscillation systems. On the other hand, the new approach yields simple approximate analytical expressions valid for small as well as large amplitudes of oscillation. In short, this new approach yields extended scope of applicability, simplicity, flexibility in application, and avoidance of complicated numerical integration as compared to the previous approaches such as the perturbation and classical harmonic balance methods. Two examples of nonlinear TDOF mass-spring systems excited by a constant external force are selected and the approximate solutions are verified with the exact solutions derived from the Jacobi elliptic function and also the numerical fourth-order Runge–Kutta solutions.     

A new study of chaotic behavior and the existence of Feigenbaum’s constants in fractional-degree Yin–Yang Hénon maps

hing dynamics in a square tank are numerically investigated when the tank is subjected to horizontal, narrowband random ground excitation. The natural frequencies of the two predominant sloshing modes are identical and therefore 1:1 internal resonance may occur. Galerkin’s method is applied to derive the modal equations of motion for nonlinear sloshing including higher modes. The Monte Carlo simulation is used to calculate response statistics such as mean square values and probability density functions (PDFs). The two predominant modes exhibit complex phenomena including “autoparametric interaction” because they are nonlinearly coupled with each other. The mean square responses of these two modes and the liquid elevation are found to differ significantly from those of the corresponding linear model, depending on the characteristics of the random ground excitation such as bandwidth, center frequency and excitation direction. It is found that the direction of the excitation is a significant factor in predicting the mean square responses. The frequency response curves for the same system subjected to equivalent harmonic excitation are also calculated and compared with the mean square responses to further explain the phenomena. Changing the liquid level causes the peak of the mean square response to shift. Furthermore, the risk of the liquid overspill from the tank is discussed by showing the three-dimensional distribution charts of the mean square responses of liquid elevations.     

Influence of friction and asymmetric freeplay on the limit cycle oscillation in aeroelastic system: An extended Hénon’s technique to temporal integration

Nonlinear effects such as friction and freeplay on the control surfaces can affect aeroelastic dynamics during flight. In particular, these nonlinearities can induce limit cycle oscillations (LCO), changing the system stability, and because of this it is essential to employ computational methods to predict this type of motion during the aircraft development cycle. In this context, the present article presents a matrix notation for describing the Hénon’s method used to reduce errors when considering piecewise linear nonlinearities in the numerical integration process. In addition, a new coordinate system is used to write the aeroelastic system of equations. The proposal defines a displacement vector with generalized and physical variables to simplify the computational implementation of the Hénon’s technique. Additionally, the article discusses the influence of asymmetric freeplay and friction on the LCO of an airfoil with control surface. The results show that the extended Hénon’s technique provides more accurate LCO predictions, that friction can change the frequency and amplitude of these motions, and the asymmetry of freeplay is important to determine the LCO behavior.     

A method for the estimation of the film thickness and plate tilt angle in thin film misaligned plate–plate rheometry

This communication describes a general procedure for the estimation of the true gap and tilt angle in commercially available torsional flow plate–plate rheometers by simply measuring the torque and normal force acting on the plates when shearing a Newtonian fluid.     

Nonlinear investigation of anti-convection and Rayleigh–Benard convection in systems with heat release at the interface

Anti-convection and Rayleigh–Benard convection generated by the joint action of external heating and heat sources (sinks) on the interface in layers with finite thicknesses are studied. Numerical simulations of the finite-amplitude convective regimes have been mage for the real two-liquid system (silicone oil 10 cs – ethylenglycol), convenient for the performance of experiments. The nonlinear boundary value problem was solved by means of the finite-difference method. Anti-convective structures in fluid systems subject to anti-convective instability only in the presence of heat sources (sinks) on the interface, have been obtained. This new type of the anti-convective motion appears in the case where one layer is strongly heated from above, while the temperature gradient in another layer is very weak.     

The incentive effect of venture capital in bilateral partnership systems with the bias mono-stable Cobb–Douglas utility

Nonlinear Dynamics - Partnerships, between multiple sides that share goals and strive for mutual benefit, are ubiquitous both between and within the enterprises, and competition and cooperation are...     

$$\pmb {H_\infty }$$ H ∞ fusion estimation of time-delayed nonlinear systems with energy constraints: the finite-horizon case

Nonlinear Dynamics - The fusion estimation issue of sensor networks is investigated for nonlinear time-varying systems with energy constraints, time delays as well as packet loss. For the addressed...     

On the use of the model proposed by Leonov for the explanation of a secondary plateau of the loss modulus in heterogeneous polymer–filler systems with agglomerates

The purpose of the presented work was to test the capability of the model proposed by Leonov (J Rheol 34:1039–1068, 1990) for the prediction of secondary plateaus on the storage and loss moduli during small-amplitude oscillatory shear flow experiments on filled or heterogeneous polymer melts. Though the occurrence of a plateau on the storage modulus can be well explained in the frame of a filler network, a plateau on the loss modulus can hardly be described with the classical models. In the Leonov model, the continuum of dissipative processes is attributed to the rupture of flocs of particles. Experiments with polyolefins filled with magnesium hydroxide show that there is a clear connection between the amount of agglomerates and the occurrence of a plateau on the loss modulus. However, the value of the critical strain for floc rupture that can be calculated from the experiment shows that the processes responsible for the low-frequency dissipation are rather changes of configuration within the agglomerates than floc rupture. These processes are not described by the Leonov model, and the predicted strain dependence of the plateau is not observed experimentally.     

Karhunen–Loève Galerkin method with decimated sampling technique for the simulation of complex fluids defined in the phase space

In complex fluids, solute molecules with structural length scales much larger than atomic are dispersed in solvents of simple fluids such as water. The rheological properties of complex fluids are determined by dynamics of solute molecules which can be modeled by the Fokker–Planck equation defined in a six-dimensional phase space. In the present investigation, we devise a method of efficient simulation of complex fluid flows employing the Karhunen–Loève Galerkin (KLG) method. Adopting the decimated sampling of solvent flow fields, a reduced-order model for the Fokker–Planck equation is obtained, which can be employed for the the simulation of complex fluids with a decent computer time. As a specific example, we consider a flow of dilute polymeric liquids over a cylinder, whose constitutive equation is the FENE (finitely extensible nonlinear elastic) model. It is found that the KLG method with the decimated sampling technique yields accurate results at a computational cost less than a hundredth of that for the numerical simulation using the Fokker–Planck equation. The KLG method supplemented by the decimated sampling technique is an efficient method of coarse-graining for equations of complex fluids defined in the phase space.     

On the Application of the Levenberg–Marquardt Method in Conjunction with an Explicit Runge–Kutta and an Implicit Rosenbrock Method to Assess Burning Velocities from Confined Deflagrations

The potential of the Levenberg–Marquardt method combined with an explicit Runge–Kutta method for non-stiff systems, and, an implicit Rosenbrock method for stiff systems to investigate burning velocities using explosion bombs was explored. The implementation of this combination of methods was verified on three benchmark test problems, and, by the application of two integral balance models to laminar hydrogen-air and methane-air explosions. The methodology described here was subsequently applied to quantify the coefficients of a turbulent burning velocity correlation for a methane-air explosion in the decaying flow field of the standard 20-litre explosion sphere. The outcome of this research indicates that the usefulness of the 20-litre sphere can be extended beyond the measurement of practical explosion parameters. When combined with the methodology in this paper, turbulent burning velocity correlations can be assessed in different parts of the Borghi-diagram.