Time-Frequency Analysis of Nonlinear Systems: The Skeleton Linear Model and the Skeleton Curves

The joint time-frequency analysis method is adopted to study the nonlinear behavior varying with the instantaneous response for a class of S.D.O.F nonlinear system. A time-frequency masking operator, together with the conception of effective time-frequency region of the asymptotic signal are defined here. Based on these mathematical foundations, a so-called skeleton linear model (SLM) is constructed which has similar nonlinear characteristics with the nonlinear system. Two skeleton curves are deduced which can indicate the stiffness and damping in the nonlinear system. The relationship between the SLM and the nonlinear system, both parameters and solutions, is clarified. Based on this work a new identification technique of nonlinear systems using the nonstationary vibration data will be proposed through time-frequency filtering technique and wavelet transform in the following paper.

Analysis of a Rubber Cone Tensioned by a Concentrated Force

Using the constitutive equation of a rubber-like materials given by Gao (1997), this paper investigates the problem of a cone under tension of a concentrated force at its apex. Under consideration is the axial-symmetry case and the large strain is taken into account. The stress strain fields near the apex are obtained by both asymptotic analysis and finite element calculation. The two results are consistent well. When the cone angle is 180 degrees, the solution becomes that of non-linear Boussinesq's problem for tension case.

Nonlinear Faraday Waves in a Parametrically Excited Circular Cylindrical Container

In the cylindrical coordinate system, a singular perturbation theory of multiple-scale asymptotic expansions was developed to study single standing water wave mode by solving potential equations of water waves in a rigid circular cylinder, which is subjec

Mesofracture mechanics: a necessary link

Classical fracture mechanics is based on the premise that small scale features could be averaged to give a larger scale property such that the assumption of material homogeneity would hold. Involvement of the material microstructure, however, necessitates different characteristic lengths for describing different geometric features. Macroscopic parameters could not be freely exchanged with those at the microscopic scale level. Such a practice could cause misinterpretation of test data. Ambiguities arising from the lack of a more precise range of limitations for the definitions of physical parameters are discussed in connection with material length scales. Physical events overlooked between the macroscopic and microscopic scale could be the link that is needed to bridge the gap. The classical models for the creation of free surface for a liquid and solid are oversimplified. They consider only the translational motion of individual atoms. Movements of groups or clusters of molecules deserve attention. Multiscale cracking behavior also requires the distinction of material damage involving at least two different scales in a single simulation. In this connection, special attention should be given to the use of asymptotic solution in contrast to the full field solution when applying fracture criteria. The former may leave out detail features that would have otherwise been included by the latter. Illustrations are provided for predicting the crack initiation sites of piezoceramics. No definite conclusions can be drawn from the atomistic simulation models such as those used in molecular dynamics until the non-equilibrium boundary conditions can be better understood. The specification of strain rates and temperatures should be synchronized as the specimen size is reduced to microns. Many of the results obtained at the atomic scale should be first identified with those at the mesoscale before they are assumed to be connected with macroscopic observations. Hopefully, "mesofracture mechanics" could serve as the link to bring macrofracture mechanics closer to microfracture mechanics.

The vector fields admitting one-parameter spatial symmetry group and their reduction

For a n-dimensional vector fields preserving some n-form, the following conclusion is reached by the method of Lie group. That is, if it admits an one-parameter, n-form preserving symmetry group, a transformation independent of the vector field is constructed explicitly, which can reduce not only dimesion of the vector field by one, but also make the reduced vector field preserve the corresponding ( n - 1)-form. In partic ular, while n = 3, an important result can be directly got which is given by Me,ie and Wiggins in 1994.

Preliminary report on the energy balance for nonlinear oscillations

In this paper, a reliable technique for calculating angular frequencies of nonlinear oscillators is developed. The new algorithm offers a promising approach by constructing a Hamiltonian for the nonlinear oscillator. Some illustrative examples are given. (C) 2002 Published by Elsevier Science Ltd.

Attractor crisis and bursting in a fluid flow with two no-slip directions

Characteristic burtsing behavior is observed in a driven, two-dimensional viscous flow, confined to a square domain and subject to no-slip boundaries. Passing a critical parameter value, an existing chaotic attractor undergoes a crisis, after which the flow initially enters a transient bursting regime. Bursting is caused by ejections from and return to a limited subdomain of the phase space, whereas the precrisis chaotic set forms the asymptotic attractor of the flow. For increasing values of the control parameter the length of the bursting regime increases progressively. Passing another critical parameter value, a second crisis leads to the appearance of a secondary type of bursting, of very large dynamical range. Within the bursting regime the flow then switches in irregular intervals from the primary to the secondary type of bursting. Peak enstrophy levels for both types of bursting are associated to the collapse of a primary vortex into a quadrupolar state.

A new method for evaluation of stress intensities for interface cracks

A new method is presented for calculating the values of K-I and K-II in the elasticity solution at the tip of an interface crack. The method is based on an evaluation of the J-integral by the virtual crack extension method. Expressions for calculating K-I and K-II by using the displacements and the stiffness derivative of the finite element solution and asymptotic crack tip displacements are derived. The method is shown to produce very accurate solutions even with coarse element mesh.

Higher-order analysis of crack-tip fields in power-law hardening materials

In this paper, we present an exact higher-order asymptotic analysis on the near-crack-tip fields in elastic-plastic materials under plane strain, Mode I. A four- or five-term asymptotic series of the solutions is derived. It is found that when 1.6 < n less-than-or-equal-to 2.8 (here, n is the hardening exponent), the elastic effect enters the third-order stress field; but when 2.8< n less-than-or-equal-to 3.7 this effect turns to enter the fourth-order field, with the fifth-order field independent. Moreover, if n>3.7, the elasticity only affects the fields whose order is higher than 4. In this case, the fourth-order field remains independent. Our investigation also shows that as long as n is larger than 1.6, the third-order field is always not independent, whose amplitude coefficient K3 depends either on K1 or on both K1 and K2 (K1 and K2 arc the amplitude coefficients of the first- and second-order fields, respectively). Firmly, good agreement is found between our results and O'Dowd and Shih's numerical ones[8] by comparison.

Higher-order analysis of crack tip fields in elastic power-law hardening materials

A HIGHER-ORDER asymptotic analysis of a stationary crack in an elastic power-law hardening material has been carried out for plane strain, Mode 1. The extent to which elasticity affects the near-tip fields is determined by the strain hardening exponent n. Five terms in the asymptotic series for the stresses have been derived for n = 3. However, only three amplitudes can be independently prescribed. These are K1, K2 and K5 corresponding to amplitudes of the first-, second- and fifth-order terms. Four terms in the asymptotic series have been obtained for n = 5, 7 and 10; in these cases, the independent amplitudes are K1, K2 and K4. It is found that appropriate choices of K2 and K4 can reproduce near-tip fields representative of a broad range of crack tip constraints in moderate and low hardening materials. Indeed, fields characterized by distinctly different stress triaxiality levels (established by finite element analysis) have been matched by the asymptotic series. The zone of dominance of the asymptotic series extends over distances of about 10 crack openings ahead of the crack tip encompassing length scales that are microstructurally significant. Furthermore, the higher-order terms collectively describe a spatially uniform hydrostatic stress field (of adjustable magnitude) ahead of the crack. Our results lend support to a suggestion that J and a measure of near-tip stress triaxiality can describe the full range of near-tip states.

Singular behaviour near the tip of a sharp V-notch in a power law hardening material

This paper presents an asymptotic analysis of the near-tip stress and strain fields of a sharp V-notch in a power law hardening material. First, the asymptotic solutions of the HRR type are obtained for the plane stress problem under symmetric loading. It is found that the angular distribution function of the radial stress sigma(r) presents rapid variation with the polar angle if the notch angle beta is smaller than a critical notch angle; otherwise, there is no such phenomena. Secondly, the asymptotic solutions are developed for antisymmetric loading in the cases of plane strain and plane stress. The accurate calculation results and the detailed comparisons are given as well. All results show that the singular exponent s is changeable for various combinations of loading condition and plane problem.

Inhomogeneous behavior of antiplane shear crack in softening material with local unloading

Examined in this work is the anti-plane stress and strain near a crack in a material that softens beyond the elastic peak and unloads on a linear path through the initial state. The discontinuity in the constitutive relation is carried into the analysis such that one portion of the local solution is elliptic in character and the other hyperbolic. Material elements in one region may cross over to another as the loading is increased. Local unloading can thus prevail. Presented are the inhomogeneous character of the asymptotic stress and strain in the elliptic and hyperbolic region, in addition to the region in which the material elements had experienced unloading. No one single stress or strain coefficient would be adequate for describing crack instability.

A finite-element analysis of viscoelastically damped sandwich plates

A finite element analysis associated with an asymptotic solution method for the harmonic flexural vibration of viscoelastically damped unsymmetrical sandwich plates is given. The element formulation is based on generalization of the discrete Kirchhoff theory (DKT) element formulation. The results obtained with the first order approximation of the asymptotic solution presented here are the same as those obtained by means of the modal strain energy (MSE) method. By taking more terms of the asymptotic solution, with successive calculations and use of the Padé approximants method, accuracy can be improved. The finite element computation has been verified by comparison with an analytical exact solution for rectangular plates with simply supported edges. Results for the same plates with clamped edges are also presented.

Analysis of flexural vibration of viscoelastically damped sandwich plates

The simplified governing equations and corresponding boundary conditions of flexural vibration of viscoelastically damped unsymmetrical sandwich plates are given. The asymptotic solution of the equations is then discussed. If only the first terms of the asymptotic solution of all variables are taken as an approximate solution, the result is identical with that obtained from the Modal Strain Energy (MSE) Method. As more terms of the asymptotic solution are taken, the successive calculations show improved accuracy. With the natural frequencies and the modal loss factors of a damped sandwich plate known, one can calculate the response of the plate to various loads providing a reliable basis for engineering design.

Two-time scales inner solutions and motion of a geostrophic vortex

Based on the authors' previous work, in this paper the systematical analyses on the motion and the inner solutions of a geostrophic vortex have been presented by means of thematched asymptotic expansion method with multiple time scales (S/gh001/2 and α S/gh001/2) and space scales. It has been shown that the leading inner solutions to the core structure in two-time scales analyses are identified with the results in normal one-time scale analyses. The time averages of the first-order solutions on short time variable τ are the same as the first-order solutions obtained in one normal time scale analyses. The geostrophic vortex induces an oscillatory motion in addition to moving with the background flow. The period, amplitude andthe deviation from the mean trajectory depend on the core structure and the initial conditions. The velocity of the motion of vortex center varies periodically and the time average of the velocity on short time variable τ is equal to the value of the local mean velocity.