Nonlinear vibrations and stability of elastic and viscoelastic systems under random stationary loads

The paper deals with numerical analysis of nonlinear vibrations of viscoelastic systems under a stochastic action in the form of a Gaussian stationary process with rational spectral density. The analysis is based on numerical simulation of the original stationary process, numerical solution of the differential equations describing the motion of the system, and computation of the maximum Lyapunov exponent if the stability of this motion is studied. An example of a plate subjected to a random stationary load applied in its plane is used to consider specific issues concerning the application of the proposed method and the peculiarities of the behavior of geometrically nonlinear elastic and viscoelastic stochastic systems. Special attention is paid to the interaction of a deterministic periodic action and a stochastic action from the viewpoint of stability of the system motion. It is shown that in some cases imposing a “colored” noise may stabilize an unstable system subjected to a periodic load.     

Incompressible Fluids in Thin Domains with Navier Friction Boundary Conditions (I)

This study is motivated by problems arising in oceanic dynamics. Our focus is the Navier–Stokes equations in a three-dimensional domain Ω_ɛ, whose thickness is of order O(ɛ) as ɛ → 0, having non-trivial topography. The velocity field is subject to the Navier friction boundary conditions on the bottom and top boundaries of Ω_ɛ, and to the periodicity condition on its sides. Assume that the friction coefficients are of order O(ɛ^3/4) as ɛ → 0. It is shown that if the initial data, respectively, the body force, belongs to a large set of H¹(Ω_ɛ), respectively, L²(Ω_ɛ), then the strong solution of the Navier–Stokes equations exists for all time. Our proofs rely on the study of the dependence of the Stokes operator on ɛ, and the non-linear estimate in which the contributions of the boundary integrals are non-trivial.     

A method for calculating the lyapunov exponent spectrum of a periodically excited non-autonomous dynamical system

The relation between the Lyapunov exponent spectrum of a periodically excited non-autonomous dynamical system and the Lyapunov exponent spectrum of the corresponding autonomous system is given and the validity of the relation is verified theoretically and computationally. A direct method for calculating the Lyapunov exponent spectrum of non-autonomous dynamical systems is suggested in this paper, which makes it more convenient to calculate the Lyapunov exponent spectrum of the dynamical system periodically excited. Following the definition of the Lyapunov dimensionD _L ^(A) of the autonomous system, the definition of the Lyapunov dimensionD _L of the non-autonomous dynamical system is also given, and the difference between them is the integer 1, namely,D _L ^(A) −D_L=1. For a quasi-periodically excited dynamical system, similar conclusions are formed. Project supported by the National Natural Science Foundation of China (No. 19772027), the Science Foundation of Shanghai Municipal Commission of Education (99A01) and the Science Foundation of Shanghai Municipal Commission of Science and Technology (No. 98JC14032).     

The moment Lyapunov exponent of a Timoshenko beam under bounded noise excitation

The moment Lyapunov exponents and the Lyapunov exponent of a two-dimensional system under bounded noise excitation are studied in this paper. The method regular perturbation is applied to obtain the small noise expansion of the pth moment Lyapunov exponent and the Lyapunov exponent. The results are applied to the study of the almost-sure and moment stability of the stationary solutions of the elastic beam subjected to the stochastic axial load. The boundaries of the almost-sure and moment stability of the elastic beam as the function of the damping coefficient and characteristics of the stochastic force are obtained.     

BV Solutions of the Semidiscrete Upwind Scheme

  We consider the semidiscrete upwind scheme

Asymptotic computation of the inhomogeneity energy in a body in an external stress field

We consider a problem on an ellipsoidal inhomogeneity in an infinitely extended homogeneous isotropic elastic medium. The inhomogeneity differs from the ambient body in the elastic moduli (Poisson’s ratio ν and shear modulus μ) and in that it has intrinsic strains. We use the equivalent inclusion method to write out expressions for the Helmholtz and Gibbs free energy of the inhomogeneity as quadratic forms in the intrinsic strains and strains at infinity. The general expressions for the coefficients of these quadratic forms are written out as three rank four tensors characterizing the contribution to the energy by the plastic strain (ɛ _p ²), by the strain at infinity (ɛ ₀²), and (only for the Gibbs energy) by the cross term ɛ ⁰ ɛ _p .     

Stability of Degenerate Stationary Waves for Viscous Gases

This paper is concerned with the asymptotic stability of degenerate stationary waves for viscous gases in the half space. We discuss the following two cases: (1) viscous conservation laws and (2) damped wave equations with nonlinear convection. In each case, we prove that the solution converges to the corresponding degenerate stationary wave at the rate t ^−α/4 as t → ∞, provided that the initial perturbation is in the weighted space L²_a=L²(\mathbb R₊; (1+x)^a dx){L^2_\alpha=L^2({\mathbb R}_+;\,(1+x)^\alpha dx)} . This convergence rate t ^−α/4 is weaker than the one for the non-degenerate case and requires the restriction α < α_*(q), where α_*(q) is the critical value depending only on the degeneracy exponent q. Such a restriction is reasonable because the corresponding linearized operator for viscous conservation laws cannot be dissipative in L²_a{L^2_\alpha} for α > α^*(q) with another critical value α^*(q). Our stability analysis is based on the space–time weighted energy method in which the spatial weight is chosen as a function of the degenerate stationary wave.     

Crack linkup: An experimental analysis

TheT _ɛ ^* integral was used to assess stable crack growth and crack linkup in 0.8 mm thick 2024-T3 aluminum tension specimens with multiple site damage (MSD) under monotonic and cyclic loads. TheT _ɛ ^* values were obtained directly from the recorded moiré fringes on the fracture specimens with and without MSD. TheT _ɛ ^* resistance curves of these fracture specimens of different geometries were in excellent agreement with each other. The results suggest thatT _ɛ ^* is a material parameter which can be used to characterize crack growth and linkup in the absence of large overloading.T _ɛ ^* based crack growth and net-section-yield based crack linkup criteria for MSD specimens are proposed. The crack tip opening angle (CTOA) criterion can also be used to correlate crack growth larger than 2 mm.     

Rigorous Results in Steady Finger Selection in Viscous Fingering

This paper concerns the existence of a steadily translating finger solution in a Hele-Shaw cell for small but non-zero surface tension (ɛ²). Though there are numerous numerical and formal asymptotic results for this problem, we know of no mathematically rigorous results that address the selection problem. We rigorously conclude that for relative finger width λ in the range , with small, analytic symmetric finger solutions exist in the asymptotic limit of surface tension if and only if the Stokes constant for a relatively simple nonlinear differential equation is zero. This Stokes constant S depends on the parameter and earlier calculations by a number of authors have shown it to be zero for a discrete set of values of a. The methodology consists of proving the existence and uniqueness of analytic solutions for a weak half-strip problem for any λ in a compact subset of (0, 1). The weak problem is shown to be equivalent to the original finger problem in the function space considered, provided we invoke a symmetry condition. Next, we consider the behavior of the solution in a neighborhood of an appropriate complex turning point for the restricted case , for some . This turning point accounts for exponentially small terms in ɛ, as ɛ→0⁺ that generally violate the symmetry condition. We prove that the symmetry condition is satisfied for small ɛ when the parameter a is constrained appropriately. (Accepted July 4, 2002 Published online January 15, 2003) Communicated by F. OTTO     

Homogenization of a Boundary-Value Problem with Varying type of Boundary Conditions in a Thick Two-Level Junction

We consider a mixed boundary-value problem for the Poisson equation in a plane two-level junction Ω _ɛ that is the union of a domain Ω₀ and a large number 2N of thin rods with thickness of order ɛ = (N ⁻¹). Depending on their lengths, the thin rods are divided into two levels. In addition, the rods from each level are ɛ-periodically alternated. Inhomogeneous Neumann boundary conditions are given on the vertical sides of the thin rods of the first level, and homogeneous Dirichlet boundary conditions are given on the vertical sides of the rods of the second level. We investigate the asymptotic behavior of a solution of this problem as ɛ → 0 and prove a convergence theorem and the convergence of the energy integral. __________ Translated from Neliniini Kolyvannya, Vol. 8, No. 2, pp. 241–257, April–June, 2005.     

Gelation behavior of polysaccharide-based interpenetrating polymer network (IPN) hydrogels

We report the preparation and rheological characterization of interpenetrating polymer network (IPN) hydrogels made from alginate and hydrophobically modified ethyl hydroxyl ethyl cellulose (HMEHEC). To our knowledge, there have been no studies of the gelation behavior of IPNs. We found that the rheology of these systems can be easily tuned, with the elastic modulus of the IPN strongly dependent on the relative ratio of HMEHEC to alginate. The sol–gel transition of these systems was found to satisfy the Winter–Chambon criterion for gelation at various crosslinker densities. From the power law relationship of the dynamic moduli (G ^′ ~G ^″ ~ω ⁿ), the exponent n appears to be dependent on both the crosslinker density and relative amount of two polymers. The value of n was found to be ~0.5 for all samples for stoichiometric amounts of crosslinker. The effect of molecular weight of HMEHEC on the gel point and viscoelastic exponent has also been reported. Alginate seems to dominate the kinetics of the process but the effect of high molecular weight HMEHEC on the gel point, especially at lower proportion was also evident.     

Approximate stationary solution and stochastic stability for a class of differential equations with parametric colored noise

This paper aims to study a class of differential equations with parametric Gaussian colored noise. We present the general framework to get the solvability conditions of the approximate stationary probability density function, which is determined by the Fokker-Planck-Kolmogorov (FPK) equations. These equations are derived using the stochastic averaging method and the operator theory with the perturbation technique. An illustrative example is proposed to demonstrate the procedure of our proposed method. The analytical expression of approximate stationary probability density function is obtained. Numerical simulation is carried out to verify the analytical results and excellent agreement can be easily found. The FPK equation for the probability density function of order ε ⁰ is used to examine the almost-sure stability for the amplitude process. Finally, the stability in probability of the amplitude process is investigated by Lin and Cai’s method.     

Periodic solutions of dynamical systems

Summary In this paper we look for T-periodic solutions of dynamical systems. Particularly we consider the system whereU ɛC ¹(ℝ ⁿ x x ℝ, ℝ),U(x, t + T)=U(x,t) ∀ x ℝ ⁿ , ∀t ɛ ℝ T>0. We assume that the problem is asymptotically linear with a bounded nonlinearity. Under a resonance assumption, we find a multiplicity of T-periodic solutions for T large enough.

---

Sommario In questo lavoro si cercano soluzioni periodiche di periodo T assegnato di sistemi dinamici. In particolare si considera un sistema di n equazioni differenziali del secondo ordine del tipo doveU ɛC ¹(ℝ ⁿ x x ℝ, ℝ),U(x, t + T)=U(x,t) ∀ x ℝ ⁿ , ∀t ɛ ℝ T>0. Nel caso in cui il problema sia asintoticamente lineare, con termine nonlineare limitato e in condizioni di risonanza, troviamo che esiste tale che per il sistema ha una molteplicità di soluzioni.

---

---

Presented at the VII A.I.M.E.T.A. and supported by M.P.I. (40% and 60%).     

On near-wall hot-wire measurements

A specially constructed hot-wire probe was used to obtain very near-wall velocity measurements in both a fully developed turbulent channel flow and flat plate boundary layer flow. The near-wall hot-wire probe, having been calibrated in a specially constructed laminar flow calibration rig, was used to measure the mean streamwise velocity profile, distributions of streamwise and spanwise intensities of turbulence and turbulence kinetic energy k in the viscous sublayer and beyond; these distributions compare very favorably with available DNS results obtained for channel flow. While low Reynolds number effects were clearly evident for the channel flow, these effects are much less distinct for the boundary layer flow. By assuming the dissipating range of eddy sizes to be statistically isotropic and the validity of Taylor's hypothesis, the dissipation rate ɛ _iso in the very near-wall viscous sublayer region and beyond was determined for both the channel and boundary layer flows. It was found that if the convective velocity U _c in Taylor's hypothesis was assumed to be equal to the mean velocity Uˉ at the point of measurement, the value of (ɛ⁺ _iso)₁ thus obtained agrees well with that of (ɛ ⁺)_DNS for y ⁺ ≥ 80 for channel flow; this suggests the validity of assuming U _c=Uˉ and local isotropy for large values of y ⁺. However, if U _c was assumed to be 10.6u _τ , the value of (ɛ⁺ _iso)₂ thus obtained was found to compare reasonably well with the distribution of (ɛ⁺ _iso)_DNS for y ⁺≤ 15. Received: 31 May 1999/Accepted: 20 December 1999     

Random impacts of a complex damped system

We investigated the random impacts of a complex damped system. Firstly the interested deterministic complex damped system was revisited and the unstable periodic attractors could be found by means of Poincaré map, time evolution and phase plot since the top Lyapunov exponent could not be applied to decide the unstable states of the proposed system. Secondly the stochastic complex damped system was examined and random impacts would be discovered, namely, the initial deterministic system will be stabilized using the stochastic force properly. The top Lyapunov exponent versus the noise intensity will be observed and one can find the change of dynamical behaviors from instability to stability. Also we implemented Poincaré map analysis, time history and phase plot to confirm the obtained results of top Lyapunov exponent, and we can find excellent agreement between these results. Therefore random noise can be applied to control the dynamical behaviors.     

Minimizing Sequences Selected via Singular Perturbations, and their Pattern Formation

We study the Cahn-Hilliard energy E _ɛ(u) over the unit square under the constraint of a constant mass m with (ɛ > 0) and without ɛ= 0) interfacial energy. Minimizers of E ₀(u) have no preferred pattern and we select patterns via sequences of conditionally critical points of E _ɛ(u) converging to minimizers as ɛ tends to zero. Those critical points are not minimizers if the singular limit has no minimal interface. We obtain them by a global bifurcation analysis of the Euler-Lagrange equations for E _ɛ(u) where the mass m is the bifurcation parameter. We make use of the symmetry of the unit square, and the elliptic maximum principle, in turn, implies that the location of maxima and minima is fixed for all solutions on global branches. This property is used to guarantee the existence of a singular limit and to verify the Weierstrass-Erdmann corner condition which proves its minimizing property. Accepted January 21, 2000?Published online November 24, 2000     

Averaging the acoustics equations for a viscoelastic material with channels filled with a viscous compressible fluid

A mathematical model describing small oscillations of a combined medium consisting of an ɛ-periodic porous viscoelastic material and a viscous compressible fluid filling the pores is considered. For this model an effective averaged model is developed and the convergence of the solutions of prelimiting problems, as ɛ → 0, to the solution of the averaged problem in the L ² space norm is proved.     

Identification of Heterogeneous Constitutive Parameters in a Welded Specimen: Uniform Stress and Virtual Fields Methods for Material Property Estimation

Local strain data obtained throughout the entire weld region encompassing both the weld nugget and heat affected zones (HAZs) are processed using two methodologies, uniform stress and virtual fields, to estimate specific heterogeneous material properties throughout the weld zone. Results indicate that (a) the heterogeneous stress–strain behavior obtained by using a relatively simple virtual fields model offers a theoretically sound approach for modeling stress–strain behavior in heterogeneous materials, (b) the local stress–strain results obtained using both a uniform stress assumption and a simplified uniaxial virtual fields model are in good agreement for strains ɛ _xx < 0.025, (c) the weld nugget region has a higher hardening coefficient, higher initial yield stress and a higher hardening exponent, consistent with the fact that the steel weld is overmatched and (d) for ɛ _xx > 0.025, strain localization occurs in the HAZ region of the specimen, resulting in necking and structural effects that complicate the extraction of local stress strain behavior using either of the relatively simple models.

Lyapunov exponent and stochastic stability of quasi-non-integrable Hamiltonian systems

An n degree-of-freedom (DOF) non-integrable Hamiltonian system subject to light damping and weak stochastic excitation is called quasi-non-integrable Hamiltonian system. In the present paper, the stochastic averaging of quasi-non-integrable Hamiltonian systems is briefly reviewed. A new norm in terms of the square root of Hamiltonian is introduced in the definitions of stochastic stability and Lyapunov exponent and the formulas for the Lyapunov exponent are derived from the averaged Itô equations of the Hamiltonian and of the square root of Hamiltonian. It is inferred that the Lyapunov exponent so obtained is the first approximation of the largest Lyapunov exponent of the original quasi-non-integrable Hamiltonian systems and the necessary and sufficient condition for the asymptotic stability with probability one of the trivial solution of the original systems can be obtained approximately by letting the Lyapunov exponent to be negative. This inference is confirmed by comparing the stability conditions obtained from negative Lyapunov exponent and by examining the sample behaviors of averaged Hamiltonian or the square root of averaged Hamiltonian at trivial boundary for two examples. It is also verified by the largest Lyapunov exponent obtained using small noise expansion for the second example.     

Reversible Bifurcation of Homoclinic Solutions in Presence of an Essential Spectrum

We consider bifurcations of a class of infinite dimensional reversible dynamical systems which possess a family of symmetric equilibria near the origin. We also assume that the linearized operator at the origin L_ɛ has an essential spectrum filling the entire real line, in addition to the simple eigenvalue at 0. Moreover, for parameter values ɛ < 0 there is a pair of imaginary eigenvalues which meet in 0 for ɛ = 0, and which disappear for ɛ > 0. The above situation occurs for example when one looks for travelling waves in a system of superposed perfect fluid layers, one being infinitely deep. We give quite general assumptions which apply in such physical examples, under which one obtains a family of bifurcating solutions homoclinic to every equilibrium near the origin. These homoclinics are symmetric and decay algebraically at infinity, being approximated at main order by the Benjamin–Ono homoclinic. For the water wave example, this corresponds to a family of solitary waves, such that at infinity the upper layer slides with a uniform velocity, over the bottom layer (at rest).