Identification of stiffness, dissipation and input parameters of randomly excited non-linear systems: Capability of restricted potential models (RPM)

A dynamic identification technique in the time domain for time invariant systems under random external forces is presented. This technique is based on the use of the class of restricted potential models (RPM), which are characterized by a non-linear stiffness and a special form of damping, that is a product of the input power spectral density (PSD) matrix and the velocity gradient of a non-linear function of the total mechanical energy. By applying stochastic differential calculus and by specific analytical manipulations, some algebraic equations, depending on the response statistics and on the mechanic parameters that characterize RPM, are obtained. These equations can be used for the dynamic identification of the above mechanic parameters once the response statistics of the system to be identified are evaluated. The proposed technique allows one to identify single-degree-of-freedom or multi-degrees-of-freedom systems in the case of unmeasurable input. Further, the probabilistic characteristics of the external forces can be completely estimated in terms of PSD matrix.     

The free energy of mixing for n-variant martensitic phase transformations using quasi-convex analysis

The construction of effective models for materials that undergo martensitic phase transformations requires usable and accurate functional representations for the free energy density. The general representation of this energy is known to be highly non-convex; it even lacks the property of quasi-convexity. A quasi-convex relaxation, however, does permit one to make certain estimates and powerful conclusions regarding phase transformation. The general expression for the relaxed free energy is however not known in the n-variant case. Analytic solutions are known only for up to 3 variants, whereas cases of practical interests involve 7-13 variants. In this study we examine the n-variant case utilizing relaxation theory and produce a seemingly obvious but very powerful observation regarding a lower bound to the quasi-convex relaxation that makes practical evolutionary computations possible. We also examine in detail the 4-variant case where we explicitly show the relation between three different forms of the free energy of mixing: upper bound by lamination, the Reuß lower bound, and a lower estimate of the -measure bound. A discussion of the bounds and their utility is provided; sample computations are presented for illustrative purposes.     

Bifurcation and universal unfolding for a rotating shaft with unsymmetrical stiffness

The 1/2 subharmonic resonance bifurcation and universal unfolding are studied for a rotating shaft with unsymmetrical stiffness. The bifurcation behavior of the response amplitude with respect to the detuning parameter was studied for this class of problems by Xiao et al. Obviously, it is highly important to research the bifurcation behavior of the response amplitude with respect to the unsymmetry of stiffness for this problem. Here, by means of the singularity theory, the bifurcation and universal unfolding of amplitude with respect to the unsymmetrical stiffness parameter are discussed. The results indicate that it is a high codimensional bifurcation problem with codimension 5, and the universal unfolding is given. From the mechanical background, we study four forms of two parameter unfoldings contained in the universal unfolding. The transition sets in the parameter plane and the bifurcation diagrams are plotted. The results obtained in this paper show rich bifurcation phenomena and provide some guidance for the analysis and design of dynamic buckling experiments of this class of system, especially, for the choice of system parameters. The project supported by the National Natural Science Foundation of China (19990510), the National Key Basic Research Special Foundation (G1998020316) and Liuhui Center for Applied Mathematics, Nankai University and Tianjin University     

A phenomenological model of electrorheological fluids

The fundamental assumption of the paper is that the extra stress tensor of an electrorheological fluid is an isotropic tensor valued function of the rate of strain tensor D and the vector n (which characterizes the orientation and length N of the fibers formed by application of an electric field). The resulting constitutive equation for is supplemented by the solution of the previously studied time evolution equation for n. Plastic behavior for the shear and normal stresses is predicted. Anticipating that the action of increasing shear rate is i) to orient the fibers more and more in the direction of flow and ii) simultaneously to break up the fibers leads to the conclusion that for the same behavior is encountered as without an electric field. Using realistically possible approximation formulas for the dependence of and N on leads to the Bingham behavior for and power law behavior for large shear rates.

Computing Green's function of elasticity in a half-plane with impedance boundary condition

This Note presents an effective and accurate method for numerical calculation of the Green's function G associated with the time harmonic elasticity system in a half-plane, where an impedance boundary condition is considered. The need to compute this function arises when studying wave propagation in underground mining and seismological engineering. To theoretically obtain this Green's function, we have drawn our inspiration from the paper by Durán et al. (2005), where the Green's function for the Helmholtz equation has been computed. The method consists in applying a partial Fourier transform, which allows an explicit calculation of the so-called spectral Green's function. In order to compute its inverse Fourier transform, we separate as a sum of two terms. The first is associated with the whole plane, whereas the second takes into account the half-plane and the boundary conditions. The first term corresponds to the Green's function of the well known time-harmonic elasticity system in (cf. J. Dompierre, Thesis). The second term is separated as a sum of three terms, where two of them contain singularities in the spectral variable (pseudo-poles and poles) and the other is regular and decreasing at infinity. The inverse Fourier transform of the singular terms are analytically computed, whereas the regular one is numerically obtained via an FFT algorithm. We present a numerical result. Moreover, we show that, under some conditions, a fourth additional slowness appears and which could produce a new surface wave. To cite this article: M. Durán et al., C. R. Mecanique 334 (2006).     

Stability of a shallow arch with one end moving at constant speed

In this paper we consider a shallow arch with rise parameter h, free of lateral loading, but subject to prescribed end motion e with constant speed c. Attention is focused on finding out whether dynamic snap-through will occur. Quasi-static analysis is first performed to identify all equilibrium configurations and their stability properties when e and h are specified. If the arch is stretched quasi-statically, it will be straightened up and no snap-through will occur. However, when the speed c is not negligible it is possible for the arch to snap to the other side dynamically. Careful analysis shows that the only possible situation when dynamic snap-through may occur is and . In this case, to prevent dynamic snap-through to occur the end speed c must not exceed a critical speed, which is a function of e and h. The minimum critical stretching speed is found to be 25.9 for all possible combinations of e and h.     

Exact static and dynamic critical loads of a sinusoidal arch under a point force at the midpoint

This paper derives the exact static and dynamic critical loads for a pinned sinusoidal arch under a concentrated force at the midpoint. For quasi-static loading, the exact critical load can be derived analytically when the rise parameter h is greater than 4.81. In the case when the concentrated force is applied suddenly, exact dynamic critical load can be formulated when the equilibrium configurations exist.     

Effective flow of a viscous liquid through a helical pipe

We study the flow of a viscous fluid through a pipe with helical shape parameterized with , where the small parameter stands for the distance between two coils of the helix. The pipe has small cross-section of size . Using the asymptotic analysis of the microscopic flow described by the Navier–Stokes system, with respect to the small parameter that tends to zero, we find the effective fluid flow described by an explicit formula of the Poisseuile type including a small distorsion due to the particular geometry of the pipe. To cite this article: E. Marušić-Paloka, I. Pažanin, C. R. Mecanique 332 (2004).

Résumé

On considère un écoulement dans un tube de section circulaire et de forme hélicoïdale paramétré par , où est la distance entre deux tours de la spirale. Le rayon de la section du tube est lui aussi supposé égal à . A partir de l'écoulement microscopique décrit par le système de Navier–Stokes et en utilisant l'analyse asymptotique par rapport à ce petit paramètre on obtient l'écoulemment effectif décrit par une formule explicite de type Poiseuille associée à une petite déviation due à la géometrie du tube. Pour citer cet article : E. Marušić-Paloka, I. Pažanin, C. R. Mecanique 332 (2004).     

Ši’lnikov Chaos in the Generalized Lorenz Canonical Form of Dynamical Systems

This paper studies the generalized Lorenz canonical form of dynamical systems introduced by elikovský and Chen [International Journal of Bifurcation and Chaos 12(8), 2002, 1789]. It proves the existence of a heteroclinic orbit of the canonical form and the convergence of the corresponding series expansion. The ilnikov criterion along with some technical conditions guarantee that the canonical form has Smale horseshoes and horseshoe chaos. As a consequence, it also proves that both the classical Lorenz system and the Chen system have ilnikov chaos. When the system is changed into another ordinary differential equation through a nonsingular one-parameter linear transformation, the exact range of existence of ilnikov chaos with respect to the parameter can be specified. Numerical simulation verifies the theoretical results and analysis.     

A perturbation-incremental method for strongly non-linear non-autonomous oscillators

A perturbation-incremental method is extended for the analysis of strongly non-linear non-autonomous oscillators of the form , where g(x) and are arbitrary non-linear functions of their arguments, and ε can take arbitrary values. Limit cycles of the oscillators can be calculated to any desired degree of accuracy and their stabilities are determined by the Floquet theory. Branch switching at period-doubling bifurcation along a frequency-response curve is made simple by the present method. Subsequent continuation of an emanating branch is also discussed.     

Collision operator including large-angle scattering for moderately coupled plasma

A new approach has been developed to treat the large-angle as well as the small-angle binary collisions in high temperature and high density plasmas when the test particle distribution function f_α is even function about the test velocity and the relations of the mass and the velocity between the test particles and the field particles are satisfied with m_αm_β (such as electron–ion collision or Lorentz-gas model) and . With the approach, the Boltzmann collision operator is derived to be suitable for the plasma considered as weakly coupled (Coulomb logarithm ) and moderately coupled , i.e., for the electron–ion coupling constant Γ_ei<0.1. The modified collision operator has a direct and practical connection to the Rosenbluth potentials, the new reduced electron–ion collision operator differs from the original Fokker–Planck operator for Coulomb collisions by terms of order . Moreover, some calculations of relaxation rate and transport properties are given for new reduced electron–ion collision operator that shows corrections.     

Homoclinic Twist Bifurcation in a System of Two Coupled Oscillators

We analyse the dynamics of two identical Josephson junctions coupled through a purely capacitive load in the neighborhood of a degenerate symmetric homoclinic orbit. A bifurcation function is obtained applying Lin's version of the Lyapunov–Schmidt reduction. We locate in parameter space the region of existence of n-periodic orbits, and we prove the existence of n-homoclinic orbits and bounded nonperiodic orbits. A singular limit of the bifurcation function yields a one-dimensional mapping which is analyzed. Numerical computations of nonsymmetric homoclinic orbits have been performed, and we show the relevance of these computations by comparing the results with the analysis.