Estimation of the dynamical properties of polyurethane foam through use of Prony series

A system identification procedure is formulated for estimation of parameters associated with a dynamic model of a single-degree-of-freedom foam-mass system. The foam is modelled as a linear viscoelastic material, whose constitutive law is expressed by an exponential hereditary relaxation kernel. The identification procedure is based on modelling the free response of the system as a Prony series (sum of exponentials terms) and fitting this Prony series to the data. This estimated response model is then utilized to estimate the parameters in the system model based on an explicit solution of the model. The procedure is analyzed for its reliability under different sources of error and uncertainties, such as the presence of weak components and experimental noise, and some modifications are evaluated to improve the robustness of the procedure. Finally, the procedure is applied to experimental data to obtain relevant stiffness, viscous and viscoelastic parameters associated with the system. Variations in values of these parameters as a function of static compression are also investigated.     

Rayleigh-Lamb wave propagation on a fractional order viscoelastic plate

A previous study of the authors published in this journal focused on mechanical wave motion in a viscoelastic material representative of biological tissue [Meral et al., J. Acoust. Soc. Am. 126, 3278-3285 (2009)]. Compression, shear and surface wave motion in and on a viscoelastic halfspace excited by surface and sub-surface sources were considered. It was shown that a fractional order Voigt model, where the rate-dependent damping component that is dependent on the first derivative of time is replaced with a component that is dependent on a fractional derivative of time, resulted in closer agreement with experiment as compared with conventional (integer order) models, such as those of Voigt and Zener. In the present study, this analysis is extended to another configuration and wave type: out-of-plane response of a viscoelastic plate to harmonic anti-symmetric Lamb wave excitation. Theoretical solutions are compared with experimental measurements for a polymeric tissue mimicking phantom material. As in the previous configurations the fractional order modeling assumption improves the match between theory and experiment over a wider frequency range. Experimental complexities in the present study and the reliability of the different approaches for quantifying the shear viscoelastic properties of the material are discussed.     

基于分数阶流变模型的铁基块体非晶合金黏弹性行为研究

采用分数阶黏弹单元替代经典模型中的黏壶, 结合非晶合金在外加载荷作用下的微观结构演化, 建立了以分数阶微积分表示的非晶合金黏弹性本构模型. 并根据Hertz弹性理论及分数阶黏弹性本构模型, 推导了块体非晶合金在纳米压痕球形压头下的位移与载荷及时间关系式. 基于推导的解析式, 对铁基块体非晶合金在表观弹性区的纳米压痕位移与载荷及时间曲线进行了非线性拟合分析. 相较于整数阶模型, 分数阶模型不仅具有较高的拟合精度, 其拟合参数能敏锐地反应加载速率对块体非晶合金黏弹性行为的影响, 且参数的变化规律与载荷作用下非晶合金微观结构演化呈现出较强的相关性.     

Analysis of two colliding fractionally damped spherical shells in modelling blunt human head impacts

The collision of two elastic or viscoelastic spherical shells is investigated as a model for the dynamic response of a human head impacted by another head or by some spherical object. Determination of the impact force that is actually being transmitted to bone will require the model for the shock interaction of the impactor and human head. This model is indended to be used in simulating crash scenarios in frontal impacts, and provide an effective tool to estimate the severity of effect on the human head and to estimate brain injury risks. The model developed here suggests that after the moment of impact quasi-longitudinal and quasi-transverse shock waves are generated, which then propagate along the spherical shells. The solution behind the wave fronts is constructed with the help of the theory of discontinuities. It is assumed that the viscoelastic features of the shells are exhibited only in the contact domain, while the remaining parts retain their elastic properties. In this case, the contact spot is assumed to be a plane disk with constant radius, and the viscoelastic features of the shells are described by the fractional derivative standard linear solid model. In the case under consideration, the governing differential equations are solved analytically by the Laplace transform technique. It is shown that the fractional parameter of the fractional derivative model plays very important role, since its variation allows one to take into account the age-related changes in the mechanical properties of bone.     

THE DYNAMICS OF A VIBROMACHINE WITH PARAMETRIC EXCITATION

The dynamics of a vibration machine with piecewise linear elastic ties under parametric harmonic excitation is investigated. Different designs of elastic elements with periodically time-varying elasticity are described. Specific non-linear features of parametric oscillations in the system under study are revealed (the invariance of parametric vibration regime to possible disturbance of phase co-ordinates, conditions of limitedness of amplitude of parametric vibrations, spectral features of non-linear parametric regimes, etc.). By the utilization of these non-linear effects, a procedure for the design of the main parameters of a parametric vibromachine is proposed.     

Viscoelastic analysis of multi-body systems using the finite element method

A study of the effect of viscoelastic material damping on the dynamic response of multibody systems, consisting of interconnected rigid, elastic and viscoelastic components, is presented. The motion of each elastic or viscoelastic body is identified by using three sets of modes: rigid body, reference and normal modes. Rigid body modes describe translation and large angular rotation of a body reference. Reference modes are the result of imposing the body-axis conditions. Normal modes define the deformation of the body relative to the body reference. Constraints between different components are formulated by using a set of non-linear algebraic equations that can be introduced to the dynamic formulation by using a Lagrange multiplier technique or can be utilized to eliminate dependent co-ordinates by partitioning the constraint Jacobian matrix. In developing the system equations of motion of the viscoelastic component, an assumption of a linear viscoelastic model is made. A Kelvin-Voigt model is employed, wherein the stress is assumed to be proportional to the strain and its time derivative. The formulation yields a constant damping matrix and the damping forces depend only on the local deformation; thus, no additional coupling between the reference and elastic co-ordinates appears in the formulation when considering the viscoelastic effects. It is demonstrated, by a numerical example, that the viscoelastic material damping can have a significant effect on the dynamic response of multibody systems.     

MHD Maxwell flow modeled by fractional derivatives with chemical reaction and thermal radiation

Heat transfer in a time-dependent flow of incompressible viscoelastic Maxwell fluid induced by a stretching surface has been investigated under the effects of heat radiation and chemical reaction. The magnetic field is applied perpendicular to the direction of flow. Velocity, temperature, and concentration are functions of z and t for the modeled boundary-layer flow problem. To have a hereditary effect, the time-fractional Caputo derivative is incorporated. The pressure gradient is assumed to be zero. The governing equations are non-linear, coupled and Boussinesq approximation is assumed for the formulation of the momentum equation. To solve the derived model numerically, the spatial variables are discretized by employing the finite element method and the Caputo-time derivatives are approximated using finite difference approximations. It reveals that the fractional derivative strengthens the flow field. We also observe that the magnetic field and relaxation time suppress the velocity. The lower Reynolds number enhances the viscosity and thus motion weakens slowly. The velocity initially decreases with increasing unsteadiness parameter δ. Temperature is an increasing function of heat radiation parameter but a decreasing one for the volumetric heat absorption parameter. The increasing value of the chemical reaction parameter decreases concentration. The Prandtl and Schmidt numbers adversely affect the temperature and concentration profiles respectively. The fractional parameter changes completely the velocity profiles. The Maxwell fluids modeled by the fractional differential equations flow faster than the ordinary fluid at small values of the time t but become slower for large values of the time t.     

Size-dependent geometrically nonlinear free vibration analysis of fractional viscoelastic nanobeams based on the nonlocal elasticity theory

In recent decades, mathematical modeling and engineering applications of fractional-order calculus have been extensively utilized to provide efficient simulation tools in the field of solid mechanics. In this paper, a nonlinear fractional nonlocal Euler–Bernoulli beam model is established using the concept of fractional derivative and nonlocal elasticity theory to investigate the size-dependent geometrically nonlinear free vibration of fractional viscoelastic nanobeams. The non-classical fractional integro-differential Euler–Bernoulli beam model contains the nonlocal parameter, viscoelasticity coefficient and order of the fractional derivative to interpret the size effect, viscoelastic material and fractional behavior in the nanoscale fractional viscoelastic structures, respectively. In the solution procedure, the Galerkin method is employed to reduce the fractional integro-partial differential governing equation to a fractional ordinary differential equation in the time domain. Afterwards, the predictor–corrector method is used to solve the nonlinear fractional time-dependent equation. Finally, the influences of nonlocal parameter, order of fractional derivative and viscoelasticity coefficient on the nonlinear time response of fractional viscoelastic nanobeams are discussed in detail. Moreover, comparisons are made between the time responses of linear and nonlinear models.     

Using wavelet multi-resolution nature to accelerate the identification of fractional order system

Because of the fractional order derivatives, the identification of the fractional order system(FOS) is more complex than that of an integral order system(IOS). In order to avoid high time consumption in the system identification, the leastsquares method is used to find other parameters by fixing the fractional derivative order. Hereafter, the optimal parameters of a system will be found by varying the derivative order in an interval. In addition, the operational matrix of the fractional order integration combined with the multi-resolution nature of a wavelet is used to accelerate the FOS identification, which is achieved by discarding wavelet coefficients of high-frequency components of input and output signals. In the end, the identifications of some known fractional order systems and an elastic torsion system are used to verify the proposed method.     

Finite element analysis and experimental study on dynamic properties of a composite beam with viscoelastic damping

This paper investigates the frequency dependent viscoelastic dynamics of a multifunctional composite structure from finite element analysis and experimental validation. The frequency-dependent behavior of the stiffness and damping of a viscoelastic material directly affects the system's modal frequencies and damping, and results in complex vibration modes and differences in the relative phase of vibration. A second order three parameter Golla–Hughes–McTavish (GHM) method and a second order three fields Anelastic Displacement Fields (ADF) approach are used to implement the viscoelastic material model, enabling the straightforward development of time domain and frequency domain finite elements, and describing the frequency dependent viscoelastic behavior. Considering the parameter identification a strategy to estimate the fractional order of the time derivative and the relaxation time is outlined. Agreement between the curve fits using both the GHM and ADF and experiment is within 0.001 percent error. Continuing efforts are addressing the material modulus comparison of the GHM and the ADF model. There may be a theoretical difference between viscoelastic degrees of freedom at nodes and elements, but their numerical results are very close to each other in the specific frequency range of interest. With identified model parameters, numerical simulation is carried out to predict the damping behavior in its first two vibration modes. The experimental testing on the layered composite beam validates the numerical predication. Experimental results also show that elastic modulus measured from dynamic response yields more accurate results than static measurement, such as tensile testing, especially for elastomers.     

Stochastic response of van der Pol oscillator with two kinds of fractional derivatives under Gaussian white noise excitation

This paper aims to investigate the stochastic response of the van der Pol(VDP) oscillator with two kinds of fractional derivatives under Gaussian white noise excitation.First,the fractional VDP oscillator is replaced by an equivalent VDP oscillator without fractional derivative terms by using the generalized harmonic balance technique.Then,the stochastic averaging method is applied to the equivalent VDP oscillator to obtain the analytical solution.Finally,the analytical solutions are validated by numerical results from the Monte Carlo simulation of the original fractional VDP oscillator.The numerical results not only demonstrate the accuracy of the proposed approach but also show that the fractional order,the fractional coefficient and the intensity of Gaussian white noise play important roles in the responses of the fractional VDP oscillator.An interesting phenomenon we found is that the effects of the fractional order of two kinds of fractional derivative items on the fractional stochastic systems are totally contrary.     

NON-LINEAR DYNAMIC ANALYSIS OF THE TWO-DIMENSIONAL SIMPLIFIED MODEL OF AN ELASTIC CABLE

The non-linear behavior of an elastic cable subjected to a harmonic excitation is investigated in this paper. Using Garlerkin's method and method of multiple scales, the discrete dynamical equations and a set of first order non-linear differential equations are obtained. A numerical simulation is used to obtain the steady state response and stable solutions. Finally the coupled dynamic features between the out-planar pendulation and the in-planar vibration of an elastic cable are analyzed.     

Noether symmetry for fractional Hamiltonian system

Fractional Hamiltonian systems within combined Riemann-Liouville fractional order derivative and combined Caputo fractional order derivative are established. Then Noether quasi-symmetry and conserved quantity for the fractional Hamiltonian systems are presented. Thirdly, perturbation to Noether quasi-symmetry and adiabatic invariant are studied. Several special cases are discussed in each section. And finally, two applications, i.e., the fractional Lotka biochemical oscillator model and the fractional isotropic harmonic oscillator model, are discussed to illustrate the results and methods.     

Analysis of the flow of non-Newtonian visco-elastic fluids in fractal reservoir with the fractional derivative

In both the oil reservoir engineering and seepage flow mechanics, heavy oil with relaxation property shows non-Newtonian rheological characteristics. The relationship between shear rate g& and shear stress t is nonlinear. Because of the relaxation phenomena of heavy oil flow in porous media, the equation of motion can be written as[1] 2,rrvpqkppqtrrtll秏骣+=-+琪抖桫 (1) where lv and lp are velocity relaxation and pressure retardation times. For most porous media, the above motion equation (1)…     

Sparse identification method of extracting hybrid energy harvesting system from observed data

Hybrid energy harvesters under external excitation have complex dynamical behavior and the superiority of promoting energy harvesting efficiency. Sometimes, it is difficult to model the governing equations of the hybrid energy harvesting system precisely, especially under external excitation. Accompanied with machine learning, data-driven methods play an important role in discovering the governing equations from massive datasets. Recently, there are many studies of data-driven models done in aspect of ordinary differential equations and stochastic differential equations (SDEs). However, few studies discover the governing equations for the hybrid energy harvesting system under harmonic excitation and Gaussian white noise (GWN). Thus, in this paper, a data-driven approach, with least square and sparse constraint, is devised to discover the governing equations of the systems from observed data. Firstly, the algorithm processing and pseudo code are given. Then, the effectiveness and accuracy of the method are verified by taking two examples with harmonic excitation and GWN, respectively. For harmonic excitation, all coefficients of the system can be simultaneously learned. For GWN, we approximate the drift term and diffusion term by using the Kramers-Moyal formulas, and separately learn the coefficients of the drift term and diffusion term. Cross-validation (CV) and mean-square error (MSE) are utilized to obtain the optimal number of iterations. Finally, the comparisons between true values and learned values are depicted to demonstrate that the approach is well utilized to obtain the governing equations for the hybrid energy harvester under harmonic excitation and GWN.     

Non-linear parameter estimation with Volterra series using the method of recursive iteration through harmonic probing

Volterra series provides a platform for non-linear response representation and definition of higher order frequency response functions (FRFs). It has been extensively used in non-parametric system identification through measurement of first and higher order FRFs. A parametric system identification approach has been adopted in the present study. The series response structure is explored for parameter estimation of polynomial form non-linearity. First and higher order frequency response functions are extracted from the measured response harmonic amplitudes through recursive iteration. Relationships between higher order FRFs and first order FRF are then employed to estimate the non-linear parameters. Excitation levels are selected for minimum series approximation error and the number of terms in the series is controlled according to convergence requirement. The problem of low signal strength of higher harmonics is investigated and a measurability criterion is proposed for selection of excitation level and range of excitation frequency. The procedure is illustrated through numerical simulation for a Duffing oscillator. Robustness of the estimation procedure in the presence of measurement noise is also investigated.     

The equation of state of a microinhomogeneous medium and the frequency dependence of its elastic nonlinearity

In the framework of a rheological model, a nonlinear dynamic equation of state of a microinhomogeneous medium containing nonlinear viscoelastic inclusions is derived. The frequency dependences of the effective nonlinear parameters are determined for the difference frequency and second harmonic generation processes in the case of a quadratic elastic nonlinearity. It is shown that the frequency dependence of the nonlinear elasticity of the medium is governed by the linear relaxation response of the inclusions at the primary excitation frequency, as well as by the relaxation of the inclusions at the nonlinear generation frequencies.     

More on finite element modeling of damped composite systems

Finite element procedures are developed and verified for layered beams and rings having either continuously or discontinuously constrained viscoelastic damping layers. The two configurations considered are (1) a three-layered sandwich beam or ring (closed curved beam) consisting of two thin elastic layers with a viscoelastic core in between, and (2) a damped composite made of a thin-walled elastic structure having a finite number of mass segments or elastic segments adhered to it by a viscoelastic material. Viscoelastic material dependence on frequency and temperature is accounted for. Numerical predictions of transverse driving point impedances agree very well with available experimental data.     

Forced response of neutrally buoyant inflated viscoelastic cantilevers to ocean waves

The dynamics of the neutrally buoyant inflated viscoelastic cantilevers constituting a submarine detection system is investigated. Thin shell theory is used to account for the stresses arising due to the internal pressure. A significant feature of the analysis is the use of the reduced shell equation which is similar in form to that for a vibrating beam with rotary effects. The forcing function in the form of surface wave excitation consists of a fundamental frequency and its second harmonic. Both the effects of apparent inertia and viscous drag are accounted for. The highly complicated non-linear, coupled equations are analyzed numerically. Use of the reduced form of the shell equations appears to avoid the problems of numerical instability and convergence reported by several investigators. The amount of information generated is rather enormous; however, for conciseness, only a few of the typical data, sufficient to establish trends, are presented. The results suggest that for the case of simple harmonicexcitation, the non-linear hydrodynamic drag introduces no superharmonic components into the response. The analysis provides valuable information concerning the system parameters leading to critical response and hence should prove useful in the design of inflatable structural members.     

Stochastic stability of a fractional viscoelastic column under bounded noise excitation

The stability of a viscoelastic column under the excitation of stochastic axial compressive load is investigated in this paper. The material of the column is modeled using a fractional Kelvin–Voigt constitutive relation, which leads to that the equation of motion is governed by a stochastic fractional equation with parametric excitation. The excitation is modeled as a bounded noise, which is a realistic model of stochastic fluctuation in engineering applications. The method of stochastic averaging is used to approximate the responses of the original dynamical system by a new set of averaged variables which are diffusive Markov vector. An eigenvalue problem is formulated from the averaged equations, from which the moment Lyapunov exponent is determined for the column system with small damping and weak excitation. The effects of various parameters on the stochastic stability and significant parametric resonance are discussed and confirmed by simulation results.