轮对系统的Hopf分岔研究

针对轮对系统中的非线性动力学问题, 本文基于Hopf分岔代数判据得到考虑陀螺效应的轮对系统Hopf分岔点解析表达式, 即轮对系统蛇形失稳的线性临界速度解析表达式. 基于分岔理论得到轮对系统的第一、第二Lyapunov系数表达式, 并结合打靶法分别得到不同纵向刚度下, 考虑陀螺效应与不考虑陀螺效应的轮对系统分岔图. 通过对比有无陀螺效应的轮对系统分岔图发现, 在同一纵向刚度下, 考虑陀螺效应的轮对系统线性临界速度和非线性临界速度均大于不考虑陀螺效应的轮对系统, 即陀螺效应可以提高轮对系统的运动稳定性. 基于Bautin分岔理论, 以纵向刚度和纵向速度作为参数, 分别得到考虑陀螺效应和不考虑陀螺效应的轮对系统, 从亚临界Hopf分岔到超临界Hopf分岔, 再从超临界Hopf分岔到亚临界Hopf分岔的迁移机理拓扑图. 通过对比有、无陀螺效应的轮对系统Bautin分岔拓扑图发现, 陀螺效应将改变轮对系统的退化Hopf分岔点, 但对于轮对系统Bautin分岔拓扑图的影响不大.      

Some basic questions on mechanosensing in cell–substrate interaction

Cells constantly probe their surrounding microenvironment by pushing and pulling on the extracellular matrix (ECM). While it is widely accepted that cell induced traction forces at the cell–matrix interface play essential roles in cell signaling, cell migration and tissue morphogenesis, a number of puzzling questions remain with respect to mechanosensing in cell–substrate interactions. Here we show that these open questions can be addressed by modeling the cell–substrate system as a pre-strained elastic disk attached to an elastic substrate via molecular bonds at the interface. Based on this model, we establish analytical and numerical solutions for the displacement and stress fields in both cell and substrate, as well as traction forces at the cell–substrate interface. We show that the cell traction generally increases with distance away from the cell center and that the traction-distance relationship changes from linear on soft substrates to exponential on stiff substrates. These results indicate that cell adhesion and migration behaviors can be regulated by cell shape and substrate stiffness. Our analysis also reveals that the cell traction increases linearly with substrate stiffness on soft substrates but then levels off to a constant value on stiff substrates. This biphasic behavior in the dependence of cell traction on substrate stiffness immediately sheds light on an existing debate on whether cells sense mechanical force or deformation when interacting with their surroundings. Finally, it is shown that the cell induced deformation field decays exponentially with distance away from the cell. The characteristic length of this decay is comparable to the cell size and provides a quantitative measure of how far cells feel into the ECM.     

A semi-analytical approach to Biot instability in a growing layer: Strain gradient correction,weakly non-linear analysis and imperfection sensitivity

Many experimental works have recently investigated the dynamics of crease formation during the swelling of long soft slabs attached to a rigid substrate. Mechanically, the spatially constrained growth provokes a residual strain distribution inside the material, and therefore the problem is equivalent to the uniaxial compression of an elastic layer.The aim of this work is to propose a semi-analytical approach to study the non-linear buckling behaviour of a growing soft layer. We consider the presence of a microstructural length, which describes the effect of a simple strain gradient correction in the growing hyperelastic layer, considered as a neo-Hookean material. By introducing a non-linear stream function for enforcing exactly the incompressibility constraint, we develop a variational formulation for performing a stability analysis of the basic homogeneous solution. At the linear order, we derive the corresponding dispersion relation, proving that even a small strain gradient effect allows the system to select a critical dimensionless wavenumber while giving a small correction to the Biot instability threshold. A weakly non-linear analysis is then performed by applying a multiple-scale expansion to the neutrally stable mode. By applying the global conservation of the mechanical energy, we derive the Ginzburg–Landau equation for the critical single mode, identifying a pitchfork bifurcation. Since the bifurcation is found to be subcritical for a small ratio between the microstructural length and the layer׳s thickness, we finally perform a sensitivity analysis to study the effect of the initial presence of a sinusoidal imperfection on the free surface of the layer. In this case, the incremental solution for the stream function is written as a Fourier series, so that the surface imperfection can have a cubic resonance with the linear modes. The solutions indicate the presence of a turning point close to the critical threshold for the perfect system. We also find that the inclusion of higher modes has a steepening effect on the surface profile, indicating the incipient formation of an elastic singularity, possibly a crease.     

Dynamic analysis of mathematical model of ethanol fermentation with gas stripping

In this paper, a mathematical model for ethanol fermentation with gas stripping is investigated. Firstly, the model with continuous substrate input is taken. We study the existence and local stability of two equilibrium points. According to Poincare–Bendixson Theorem, the sufficient condition for the globally asymptotical stability of positive equilibrium point is obtained, which implies that we can get stable ethanol product. Secondly, we study the model with impulsive substrate input and obtain the sufficient condition for the local stability of cell-free periodic solution by using the Floquet’s theory of impulsive differential equation and small-amplitude perturbation skills. In a certain limiting case, it is shown that a nontrivial periodic solution emerges via a supercritical (subcritical) bifurcation. Finally, our results are confirmed by means of numerical simulation.     

Slow Sinusoidal Modulation Through Bifurcations: The Effect of Additive Noise

We analyse the response of two oscillators with quadratic coupling that exhibit an internal resonance. With sinusoidal excitation, as the excitation amplitude increases, a bifurcation in the response occurs. The response of one oscillator changes from linear variation with excitation amplitude to a constant saturated value. The other mode changes from zero to large amplitude, the change sometimes being quite rapid as the excitation amplitude is very slowly increased. We consider slow sinusoidal variation of the excitation amplitude through this bifurcation. Noise must now be included in the model, as even very low-level amplitude noise can affect critically the value at which the jump occurs. Amplitude modulation can extend the range over which the zero response of one oscillator occurs, causing an effective stabilisation of that form of the response; noise on the other hand is destabilizing. We analyse these competing effects using matched asymptotic expansions. A nested set of three expansions is needed to describe the rapid jump; the innermost expansion describes how noise triggers the rapid jump. Excellent comparisons are obtained with numerical simulations. The analytic results can be used to find ranges and frequencies of the modulation and noise levels that control the system so that the zero-amplitude solution is maintained effectively at zero even into parameter ranges where the autonomous zero-amplitude solution is locally unstable.     

Supercritical as well as subcritical Hopf bifurcation in nonlinear flutter systems

The Hopf bifurcations of an airfoil flutter system with a cubic nonlinearity are investigated, with the flow speed as the bifurcation parameter. The center manifold theory and complex normal form method are Used to obtain the bifurcation equation. Interestingly, for a certain linear pitching stiffness the Hopf bifurcation is both supercritical and subcritical. It is found, mathematically, this is caused by the fact that one coefficient in the bifurcation equation does not contain the first power of the bifurcation parameter. The solutions of the bifurcation equation are validated by the equivalent linearization method and incremental harmonic balance method.     

Quasi-periodicity and chaos in a differentially heated cavity

Convective flows of a small Prandtl number fluid contained in a two-dimensional vertical cavity subject to a lateral thermal gradient are studied numerically. The chosen geometry and the values of the material parameters are relevant to semiconductor crystal growth experiments in the horizontal configuration of the Bridgman method. For increasing Rayleigh numbers we find a transition from a steady flow to periodic solutions through a supercritical Hopf bifurcation that maintains the centro-symmetry of the basic circulation. For a Rayleigh number of about ten times that of the Hopf bifurcation, the periodic solution loses stability in a subcritical Neimark–Sacker bifurcation, which gives rise to a branch of quasiperiodic states. In this branch, several intervals of frequency locking have been identified. Inside the resonance horns the stable limit cycles lose and gain stability via some typical scenarios in the bifurcation of periodic solutions. After a complicated bifurcation diagram of the stable limit cycle of the 1:10 resonance horn, a soft transition to chaos is obtained. PACS 44.25.+f, 47.20.Ky, 47.52.+j     

Post-bifurcation analysis of a thin-walled hyperelastic tube under inflation

We consider the problem of bulging, or necking, of an infinite thin-walled hyperelastic tube that is inflated by an internal pressure, with the axial stretch at infinity maintained at unity. We present a simple procedure that can be used to derive the bifurcation condition and to determine the near-critical behaviour analytically. It is shown that there is a bifurcation with zero mode number and that the associated axial variation of near-critical bifurcated configurations is governed by a first-order differential equation that admits a locally bulging or necking solution. This result suggests that the corresponding bifurcation pressure can be identified with the so-called initiation pressure which featured in recent experimental studies. This is supported by good agreement between our theoretical predictions and one set of experimental data. It is also shown that the Gent material model can support both bulging and necking solutions whereas the Varga and Ogden material models can only support bulging solutions. Relevance of the present method to the study of non-linear wave propagation in a fluid-filled distensible tube is also discussed.     

Limit cycle oscillation suppression of 2-DOF airfoil using nonlinear energy sink

This paper presents a novel mechanical attachment, i.e., nonlinear energy sink （NES）, for suppressing the limit cycle oscillation （LCO） of an airfoil. The dynamic responses of a two-degree-of-freedom （2-DOF） airfoil coupled with an NES are studied with the harmonic balance method. Different structure parameters of the NES, i.e., mass ratio between the NES and airfoil, NES offset, NES damping, and nonlinear stiffness in the NES, are chosen for studying the effect of the LCO suppression on an aeroelastic system with a supercritical Hopf bifurcation or subcritical Hopf bifurcation, respectively. The results show that the structural parameters of the NES have different influence on the supercritical Hopf bifurcation system and the subcritical Hopf bifurcation system.     

轮对非线性动力学模型稳定性和分岔特性研究

为了探究轮对系统的横向失稳问题,考虑了陀螺效应和一系悬挂阻尼的影响作用,建立非线性轮轨接触关系的轮对动力学模型,研究轮对系统的蛇行稳定性、Hopf分岔特性及迁移转化机理.通过稳定性判据获得了轮对系统失稳临界速度.采用中心流形定理和规范型方法对轮对动力学模型进行化简,得到与轮对系统分岔特性相同的一维复变量方程,理论推导求得轮对系统的第一Lyapunov系数的表达式,根据其符号即可判断轮对系统的Hopf分岔类型.讨论了不同参数对轮对系统Hopf分岔临界速度的影响,探究了轮对系统的超临界、亚临界Hopf分岔域在二维参数空间的分布规律.利用数值模拟得到轮对系统的3种典型Hopf分岔图,验证了轮对系统超临界、亚临界Hopf分岔域分布规律的正确性.结果表明,轮对系统的临界速度随着等效锥度的增大而减小,随着一系悬挂的纵向刚度和纵向阻尼的增大而增大,随着纵向蠕滑系数的增大呈先增大后减小.系统参数变化会引起轮对系统Hopf分岔类型发生改变,即亚临界与超临界Hopf分岔相互迁移转化.轮对系统Hopf分岔域在二维参数空间的分布规律对于轮对系统参数匹配和优化设计具有一定的指导意义.     

齿轮传动系统共存吸引子的不连续分岔

大量的多吸引子共存是引起齿轮传动系统具有丰富动力学行为的一个重要因素.多吸引子共存时,运动工况的变化以及不可避免的扰动都可能导致齿轮传动系统在不同运动行为之间跳跃变换,对整个机器产生不良的影响.目前,一些隐藏的吸引子没有被发现,共存吸引子的分岔演化规律没有被完全揭示.考虑单自由度直齿圆柱齿轮传动系统,构建由局部映射复合的Poincaré映射,给出Jacobi矩阵特征值计算的半解析法.应用数值仿真、延拓打靶法和Floquet特征乘子求解共存吸引子的稳定性与分岔,应用胞映射法计算共存吸引子的吸引域,讨论啮合频率、阻尼比和时变激励幅值对系统动力学的影响,揭示齿轮传动系统倍周期型擦边分岔、亚临界倍周期分岔诱导的鞍结分岔和边界激变等不连续分岔行为.倍周期分岔诱导的鞍结分岔引起相邻周期吸引子相互转迁的跳跃与迟滞,使倍周期分岔呈现亚临界特性.鞍结分岔是共存周期吸引子出现或消失的主要原因.边界激变引起混沌吸引子及其吸引域突然消失,对应周期吸引子的分岔终止.     

An analytical study of time-delayed control of friction-induced vibrations in a system with a dynamic friction model

We investigate the control of friction-induced vibrations in a system with a dynamic friction model which accounts for hysteresis in the friction characteristics. Linear time-delayed position feedback applied in a direction normal to the contacting surfaces has been employed for the purpose. Analysis shows that the uncontrolled system loses stability via. a subcritical Hopf bifurcation making it prone to large amplitude vibrations near the stability boundary. Our results show that the controller achieves the dual objective of quenching the vibrations as well as changing the nature of the bifurcation from subcritical to supercritical. Consequently, the controlled system is globally stable in the linearly stable region and yields small amplitude vibrations if the stability boundary is crossed due to changes in operating conditions or system parameters. Criticality curve separating regions on the stability surface corresponding to subcritical and supercritical bifurcations is obtained analytically using the method of multiple scales (MMS). We have also identified a set of control parameters for which the system is stable for lower and higher relative velocities but vibrates for the intermediate ones. However, the bifurcation is always supercritical for these parameters resulting in low amplitude vibrations only.     

BIFURCATIONS OF AIRFOIL IN INCOMPRESSIBLE FLOW

Bifurcations of an airfoil with nonlinear pitching stiffness in incompressible flow are investigated. The pitching spring is regarded as a spring with cubic stiffness. The motion equations of the airfoil are written as the four dimensional one order differential equations. Taking air speed and the linear part of pitching stiffness as the parameters, the analytic solutions of the critical boundaries of pitchfork bifurcations and Hopf bifurcations are obtained in 2 dimensional parameter plane. The stabilities of the equilibrium points and the limit cycles in different regions of 2 dimensional parameter plane are analyzed. By means of harmonic balance method, the approximate critical boundaries of 2-multiple semi-stable limit cycle bifurcations are obtained, and the bifurcation points of supercritical or subcritical Hopf bifurcation are found. Some numerical simulation results are given.     

Qualitative study of cavitated bifurcation for a class of incompressible generalized neo-Hookean spheres

IntroductionIn 1 958,GentandLindleyobservedthephenomenonofsuddenvoidnucleationinsolidsexperimentallyintensioningahomogenousclose_grainedvulcanizedrubbercylinderforthefirsttime.ButthemathematicalmodelonvoidnucleationandgrowthhasnotbeendescribedasabifurcationproblembasedonthetheoryofnonlinearelasticmechanicsbyBall[1]until1 982 .Inrecentyears,manyinvestigationshavebeenmadeonthisaspect.Theproblemofcavitatedbifurcationforincompressibleisotropichyperelasticmaterialswithpower_lawtypehasbeeninvestig…     

The Effect of Elastic Surface Coating on the Finite Deformation and Bifurcation of a Pressurized Circular Annulus

A plane-strain theory of an elastic solid coated with a thin elastic film on part or all of its boundary was developed recently by Steigmann and Ogden (1997a). In this paper the theory is applied to the (plane-strain) problem of a thick-walled circular cylindrical tube which is subject to both internal and external pressure and which has an elastic coating on one or both of its circular cylindrical boundaries. The effect of the coating on the symmetrical response of the annular cross-section of the tube is determined first. It is noted, in particular, that while the pressure may exhibit a maximum followed by a minimum during inflation for an uncoated tube it may be a monotonic increasing function of the radius for a coated tube with coating elastic modulus sufficiently large. Next, the possibility of bifurcation from a symmetrical configuration is examined and again the influence of the coating is analysed. The effect of a coating on the outer boundary is compared with that on the inner boundary. Specifically, during compression, coating on the outer boundary delays bifurcation compared with the uncoated case. On the other hand, when the coating is on the inner boundary, bifurcation is either delayed or advanced relative to the uncoated situation depending on the values of the bending stiffness and tube thickness parameters. Generally, bifurcation is delayed by an increase in the magnitude of the bending stiffness of the coating at fixed values of the other parameters. This revised version was published online in August 2006 with corrections to the Cover Date.     

Responses of Duffing oscillator with fractional damping and random phase

This paper aims to investigate dynamic responses of stochastic Duffing oscillator with fractional-order damping term, where random excitation is modeled as a harmonic function with random phase. Combining with Lindstedt–Poincaré (L–P) method and the multiple-scale approach, we propose a new technique to theoretically derive the second-order approximate solution of the stochastic fractional Duffing oscillator. Later, the frequency–amplitude response equation in deterministic case and the first- and second-order steady-state moments for the steady state in stochastic case are presented analytically. We also carry out numerical simulations to verify the effectiveness of the proposed method with good agreement. Stochastic jump and bifurcation can be found in the figures of random responses, and then we apply Monte Carlo simulations directly to obtain the probability density functions and time response diagrams to find the stochastic jump and bifurcation. The results intuitively show that the intensity of the noise can lead to stochastic jump and bifurcation.     

Asymmetric bifurcations of thick-walled circular cylindrical elastic tubes under axial loading and external pressure

In this paper, we consider bifurcation from a circular cylindrical deformed configuration of a thick-walled circular cylindrical tube of incompressible isotropic elastic material subject to combined axial loading and external pressure. In particular, we examine both axisymmetric and asymmetric modes of bifurcation. The analysis is based on the three-dimensional incremental equilibrium equations, which are derived and then solved numerically for a specific material model using the Adams–Moulton method. We assess the effects of wall thickness and the ratio of length to (external) radius on the bifurcation behaviour.     

Instability conditions due to structural nonlinearities in regenerative chatter

Chatter is an instability condition in machining processes characterized by nonlinear behavior, such as the presence of limit cycles, jump phenomenon, subcritical Hopf and period doubling bifurcations. Although the use of nonlinear techniques has provided a better understanding of chatter, neither a unifying model nor an exact solution has yet been developed due to the intricacy of the problem. This work proposes a weakly nonlinear model with square and cubic terms in both structural stiffness and regenerative terms, to represent self-excited vibrations in machining. An approximate solution is derived by using the method of multiple scales. In addition, a qualitative analysis of the effect of the nonlinear parameters on the stability of the system is performed. The structural cubic term gives a better representation of the nonlinear behavior, whereas the square term represents a distant attractor in the stability chart. Instability due to subcritical Hopf bifurcations is established in terms of the eigenvalues of the model in normal form. An important contribution of this analysis is the representation of hysteresis in terms of new lobes within the conventional stability limits, useful in restoring stability. This analysis leads to a further understanding of the nonlinear behavior of regenerative chatter.     

Nonlinear chatter with large amplitude in a cylindrical plunge grinding process

The phenomenon that the stable smooth grinding process coexists with chatter vibrations with large amplitudes in a cylindrical plunge grinding process is investigated in this paper. In the analyzed dynamic model, the workpiece and the grinding wheel involved in the grinding process are regarded as a slender hinged-hinged Euler?CBernoulli beam and a damped spring mass system, respectively, and the contact force between the two is treated as the main factor that affects the dynamic behaviors of the process. Called regenerative force, the contact force represents the interaction with regenerative effects between the workpiece and the wheel. To clarify the relation between the force and the dynamical behaviors in the grinding process, all the effects of the system parameters related to the interaction, such as the grinding stiffness, the rotation speeds of the workpiece and the wheel, on the dynamic motions of the process are studied. To this end, the eigenvalues analysis is firstly carried out to find the chatter-free-region, in which the smooth grinding process is stable and the chatter vibration may be absent. And then the nonlinear chatter vibrations when the values of concerned parameter leave the chatter-free region are predicted numerically. It is interesting that both the supercritical and subcritical Hopf bifurcations are found on the same boundary of the chatter-free region. As we know, there must be a zone in the chatter-free region where the stable smooth grinding process coexists with the chatter vibration when the subcritical one arises and the switching point between the supercritical and the subcritical ones is a Bautin bifurcation point mathematically. Thus, the Bautin bifurcation analysis is performed to scan the subregion in which the smooth grinding process is not unconditional stable anymore.