Dynamics of gas bubbles in viscoelastic fluids. II. Nonlinear viscoelasticity

The nonlinear oscillations of a spherical, acoustically forced gas bubble in nonlinear viscoelastic media are examined. The constitutive equation [Upper-Convective Maxwell (UCM)] used for the fluid is suitable for study of large-amplitude excursions of the bubble, in contrast to the previous work of the authors which focused on the smaller amplitude oscillations within a linear viscoelastic fluid [J. S. Allen and R. A. Roy, J. Acoust. Soc. Am. 107, 3167-3178 (2000)]. Assumptions concerning the trace of the stress tensor are addressed in light of the incorporation of viscoelastic constitutive equations into bubble dynamics equations. The numerical method used to solve the governing system of equations (one integrodifferential equation and two partial differential equations) is outlined. An energy balance relation is used to monitor the accuracy of the calculations and the formulation is compared with the previously developed linear viscoelastic model. Results are found to agree in the limit of small deformations; however, significant divergence for larger radial oscillations is noted. Furthermore, the inherent limitations of the linear viscoelastic approach are explored in light of the more complete nonlinear formulation. The relevance and importance of this approach to biomedical ultrasound applications are highlighted. Preliminary results indicate that tissue viscoelasticity may be an important consideration for the risk assessment of potential cavitation bioeffects.     

Nonlinear waves in a fluid-filled thin viscoelastic tube

In the present paper the propagation property of nonlinear waves in a thin viscoelastic tube filled with incom-pressible inviscid fluid is studied.The tube is considered to be made of an incompressible isotropic viscoelastic material described by Kelvin-Voigt model.Using the mass conservation and the momentum theorem of the fluid and radial dynamic equilibrium of an element of the tube wall,a set of nonlinear partial differential equations governing the prop-agation of nonlinear pressure wave in the solid-liquid coupled system is obtained.In the long-wave approximation the nonlinear far-field equations can be derived employing the reductive perturbation technique (RPT).Selecting the expo-nent α of the perturbation parameter in Gardner-Morikawa transformation according to the order of viscous coefficient η,three kinds of evolution equations with soliton solution,i.e.Korteweg-de Vries (KdV)-Burgers,KdV and Burgers equations are deduced.By means of the method of traveling-wave solution and numerical calculation,the propagation properties of solitary waves corresponding with these evolution equations are analysed in detail.Finally,as a example of practical application,the propagation of pressure pulses in large blood vessels is discussed.     

Jerky flow model as a basis for research in deformation instabilities

A phenomenological jerky flow model was developed in which macroscale plastic strain rates are defined by dislocation kinetics. The model takes into account destructive processes governed by shear and bulk defect accumulation. At the heart of the model lie equations of solid mechanics and relaxation-type constitutive equations. A loaded elastoplastic solid is treated as a nonlinear dynamic system whose evolution, according to synergetic laws, is much contributed by negative and positive feedbacks expressed, respectively, through constitutive equations of the first group (relaxation equations) and constitutive equations of the second group (kinetic equations for deformation defect and damage accumulation rates). The negative feedback stabilizes deformation by relaxation, bringing the process to some local dynamic equilibrium. The positive feedback destabilizes deformation, driving the system to a critical state. Numerical experiment was performed in 2D and 3D statements. Statistical analysis of stress fluctuations about the average trend shows that the jerky flow model of an elastoplastic medium demonstrates evolution characteristic of nonlinear dynamic systems: through states of dynamic chaos and self-organized criticality to a global catastrophe.     

聚合物充填过程中裹气现象的数值模拟

针对聚合物充填过程中的裹气现象,采用一种有限元（FEM）-间断有限元（DG）耦合算法对其进行数值模拟。对于自由运动界面,采用水平集（Level Set）方法进行捕捉;用XPP（eXtended Pom-Pom）本构模型来描述黏弹性流体的流变行为。采用有限元-间断有限元耦合算法求解统一的流场方程,并采用隐式间断有限元求解XPP本构方程、Level Set及其重新初始化方程。数值结果与文献中的实验结果及模拟结果吻合较好,验证了数值方法的稳定性及准确性。分析带有非规则嵌件型腔内,注射速度及浇口尺寸对裹气现象的影响,裹气容易出现在较高注射速度及较小浇口的情形。     

On spectral relaxation method approach for steady von Kármán flow of a Reiner-Rivlin fluid with Joule heating, viscous dissipation and suction/injection

In this study we use the spectral relaxation method (SRM) for the solution of the steady von Kármán flow of a Reiner-Rivlin fluid with Joule heating and viscous dissipation. The spectral relaxation method is a new Chebyshev spectral collocation based iteration method that is developed from the Gauss-Seidel idea of decoupling systems of equations. In this work, we investigate the applicability of the method in solving strongly nonlinear boundary value problems of von Kármán flow type. The SRM results are validated against previous results present in the literature and with those obtained using the bvp4c, a MATLAB inbuilt routine for solving boundary value problems. The study highlights the accuracy and efficiency of the proposed SRM method in solving highly nonlinear boundary layer type equations.     

Dynamic characteristics of traveling waves for a rotating laminated circular plate with viscoelastic core layer

The dynamic governing equations and the corresponding boundary conditions for a rotating thin laminated circular plate with a viscoelastic core layer are derived in this paper based on the Hamilton principle. The analysis on dynamic features of the forward and Backward Traveling Waves for the rotating laminated plate is performed by means of Galerkin's method. The frequency-dependent complex modulus model for describing the constitutive behavior of the viscoelastic core layer is employed. The dynamic characteristics of frequencies and dampings of traveling waves for the rotating plate are obtained numerically. The effects of geometrical and material parameters on the critical speed of the rotating laminated plate with viscoelastic core are discussed in detail.     

Nonlinear dynamics of higher-dimensional system for an axially accelerating viscoelastic beam with in-plane and out-of-plane vibrations

In this paper, the bifurcations and chaotic motions of higher-dimensional nonlinear systems are investigated for the nonplanar nonlinear vibrations of an axially accelerating moving viscoelastic beam. The Kelvin viscoelastic model is chosen to describe the viscoelastic property of the beam material. Firstly, the nonlinear governing equations of nonplanar motion for an axially accelerating moving viscoelastic beam are established by using the generalized Hamilton’s principle for the first time. Then, based on the Galerkin’s discretization, the governing equations of nonplanar motion are simplified to a six-degree-of-freedom nonlinear system and a three-degree-of-freedom nonlinear system with parametric excitation, respectively. At last, numerical simulations, including the Poincare map, phase portrait and Lyapunov exponents are used to analyze the complex nonlinear dynamic behaviors of the axially accelerating moving viscoelastic beam. The bifurcation diagrams for the in-plane and out-of-plane displacements via the mean axial velocity, the amplitude of velocity fluctuation and the frequency of velocity fluctuation are respectively presented when other parameters are fixed. The Lyapunov exponents are calculated to identify the existence of the chaotic motions. From the numerical results, it is indicated that the periodic, quasi-periodic and chaotic motions occur for the nonplanar nonlinear vibrations of the axially accelerating moving viscoelastic beam. Observing the in-plane nonlinear vibrations of the axially accelerating moving viscoelastic beam from the numerical results, it is found that the nonlinear responses of the six-degree-of-freedom nonlinear system are much different from that of the three-degree-of-freedom nonlinear system when all parameters are same.     

Dynamics of gas bubbles in viscoelastic fluids. I. Linear viscoelasticity

The nonlinear oscillations of spherical gas bubbles in linear viscoelastic fluids are studied. A novel approach is implemented to derive a governing system of nonlinear ordinary differential equations. The linear Maxwell and Jeffreys models are chosen as the fluid constitutive equations. An advantage of this new formulation is that, when compared with previous approaches, it facilitates perturbation methods and numerical investigations. Analytical solutions are obtained using a multiple scale perturbation method and compared with the Newtonian results for various Deborah numbers. Numerical analysis of the full equations supports the perturbation analysis, and further reveals significant differences between the viscoelastic and Newtonian cases. Differences in the oscillation phase and harmonic structure characterize some of the viscoelastic effects. Subharmonic excitations at particular fluid parameters lead to a discrete group modulation of the radial excursions; this appears to be a unique, previously undiscovered phenomenon. Implications for medical ultrasound applications are discussed in light of these current findings.     

Multi-pulse chaotic dynamics in non-planar motion of parametrically excited viscoelastic moving belt

This paper investigates the multi-pulse global bifurcations and chaotic dynamics for the nonlinear, non-planar oscillations of the parametrically excited viscoelastic moving belt using an extended Melnikov method in the resonant case. Using the Kelvin-type viscoelastic constitutive law and Hamilton's principle, the equations of motion are derived for the viscoelastic moving belt with the external damping and parametric excitation. Applying the method of multiple scales and Galerkin's approach to the partial differential governing equation, the four-dimensional averaged equation is obtained for the case of 1:1 internal resonance and primary parametric resonance. From the averaged equations obtained, the theory of normal form is used to derive the explicit expressions of normal form with a double zero and a pair of pure imaginary eigenvalues. Based on the explicit expressions of normal form, the extended Melnikov method is used for the first time to investigate the Shilnikov-type multi-pulse homoclinic bifurcations and chaotic dynamics. The paper demonstrates how to employ the extended Melnikov method to analyze the Shilnikov-type multi-pulse homoclinic bifurcations and chaotic dynamics of high-dimensional nonlinear systems in engineering applications. Numerical simulations show that for the nonlinear non-planar oscillations of the viscoelastic moving belt, the Shilnikov-type multi-pulse chaotic motions can occur. Overall, both theoretical and numerical studies suggest that the chaos for the Smale horseshoe sense in motion exists.     

Large-amplitude nonlinear vibrations of a Mooney–Rivlin rectangular membrane

An analysis of the linear and nonlinear vibration response and stability of a pre-stretched hyperelastic rectangular membrane under harmonic lateral pressure and finite initial deformations is presented in this paper. Geometric nonlinearity due to finite deformations and material nonlinearity associated with the hyperelastic constitutive law are taken into account. The membrane is assumed to be made of an isotropic, homogeneous, and incompressible Mooney–Rivlin material. The results for a neo-Hookean material are obtained as a particular case and a comparison of these two constitutive models is carried out. First, the exact solution of the membrane under a biaxial stretch is obtained, being this initial stress state responsible for the membrane stiffness. The equations of motion of the pre-stretched membrane are then derived. From the linearized equations, the natural frequencies and mode shapes of the membrane are analytically obtained for both materials. The natural modes are then used to approximate the nonlinear deformation field using the Galerkin method. A detailed parametric analysis shows the strong influence of the stretching ratios and material parameters on the linear and nonlinear oscillations of the membrane. Frequency–amplitude relations, resonance curves, and bifurcation diagrams, are used to illustrate the nonlinear dynamics of the membrane. The present results are compared favorably with the results evaluated for the same membrane using a nonlinear finite element formulation.     

Bifurcation analysis and amplitude equations for viscoelastic convective fluids

Summary The possible bifurcations of a convective instability in viscoelastic fluid are studied. The viscoelastic behaviour is modelized by means of the Oldroyd type fluid whose parameters can be adjusted to suit a large class of polymeric fluids. We analyse in some detail bifurcations of codimension one (stationary or oscillatory convection) and codimension two for such kind of fluids. By a weak nonlinear analysis, the coefficients of the amplitude equations corresponding to the different bifurcations are also determined. It has been found that the nature of the convective solution depends crucially on both the viscoelastic parameters and the constitutive equation used to describe the fluid.     

Dynamics of two-dimensional composites of elastic and viscoelastic layers

Models of frequency response, acoustic transmission, and transient wave propagation are presented for a two-dimensional composite of elastic and viscoelastic layers, simply supported at the two boundaries. The three models adopt transfer matrices to relate state variables over the two faces of a layer. In the frequency domain, a viscoelastic constitutive law is derived by nonlinear fitting a Padé series to measured data of complex shear modulus. For an elastic material, the eigenproblem admits positive real eigenvalues and their negatives. For a viscoelastic material, it admits positive complex eigenvalues and their negative conjugates. The imaginary part of the eigenvalue acts as a velocity-dependent viscous damper. Modal analysis solving transient response utilizes the complex eigenquantities and the static-dynamic superposition method.     

考虑频散效应的一维非线性地震波数值模拟

非线性理论是解决地学问题的重要手段, 充分认识非线性波动特征有助于解释实际观测资料中的一些特殊地震现象. 本文基于Hokstad改造的非线性本构方程, 利用交错网格有限差分法实现了固体介质中一维非线性地震波数值模拟; 采用通量校正传输方法消除非线性数值模拟中波形振荡, 提高模拟精度. 通过与解析解的对比验证了模拟结果的正确性. 研究结果显示了非线性系数对地震波的传播有重要影响, 并且当取适当的非线性和频散系数时, 地震波表现出孤立波的传播特性. 最后分析了不同的非线性地震波在固体介质中的传播演化特征.     

The Modified Ghost Fluid Method Applied to Fluid-Elastic Structure Interaction

In this work, the modified ghost fluid method is developed to deal with 2D compressible fluid interacting with elastic solid in an Euler-Lagrange coupled system. In applying the modified Ghost Fluid Method to treat the fluid-elastic solid coupling, the Navier equations for elastic solid are cast into a system similar to the Euler equations but in Lagrangian coordinates. Furthermore, to take into account the influence of material deformation and nonlinear wave interaction at the interface, an Euler-Lagrange Riemann problem is constructed and solved approximately along the normal direction of the interface to predict the interfacial status and then define the ghost fluid and ghost solid states. Numerical tests are presented to verify the resultant method.     

Higher-order stochastic averaging to study stability of a fractional viscoelastic column

Stochastic stability of a fractional viscoelastic column axially loaded by a wideband random force is investigated by using the method of higher-order stochastic averaging. By modelling the wideband random excitation as Gaussian white noise and real noise and assuming the viscoelastic material to follow the fractional Kelvin–Voigt constitutive relation, the motion of the column is governed by a fractional stochastic differential equation, which is justifiably and uniformly approximated by an averaged system of Itô stochastic differential equations. Analytical expressions are obtained for the moment Lyapunov exponent and the Lyapunov exponent of the fractional system with small damping and weak random fluctuation. The effects of various parameters on the stochastic stability of the system are discussed.     

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.     

Theoretical and experimental study of the nonlinear resonance vibration of cementitious materials with an application to damage characterization

This paper presents a theoretical and experimental study of the nonlinear flexural vibration of a cement-based material with distributed microcracks caused by an important deterioration mechanism, alkali-silica reaction (ASR). The general equation of motion is derived for the flexural vibration of a slender beam with the nonlinear hysteretic constitutive relationship for consolidated materials, and then an approximate formula for excitation-dependent resonance frequency is obtained. A downward shift of the resonance frequency is related to the nonlinearity parameters defined in the constitutive relationship. Vibration experiments are conducted on standard mortar bar samples undergoing progressive ASR damage. The absolute nonlinearity parameters are determined from these experimental results using the theoretical solution in order to investigate their dependence on the damage state of the material. With the progress of the ASR damage, the absolute value of the hysteresis nonlinearity parameter increases by as much as six times from the intact (undamaged) state in the sample with highly reactive aggregate; this is in contrast to a change of about 16% in the linear resonance frequency. It is demonstrated that the combined theoretical and experimental approach developed in this research can be used to quantitatively characterize ASR damage in mortar samples and other cement-based materials.     

Nonlinear rheological models for structured interfaces

The GENERIC formalism is a formulation of nonequilibrium thermodynamics ideally suited to develop nonlinear constitutive equations for the stress-deformation behavior of complex interfaces. Here we develop a GENERIC model for multiphase systems with interfaces displaying nonlinear viscoelastic stress-deformation behavior. The link of this behavior to the microstructure of the interface is described by including a scalar and a tensorial structural variable in the set of independent surface variables. We derive an expression for the surface stress tensor in terms of these structural variables, and a set of general nonlinear time evolution equations for these variables, coupling them to the deformation field. We use these general equations to develop a number of specific models, valid for application near equilibrium, or valid for application far beyond equilibrium.     

离散牛顿正则化方法及应用

把离散Newton法和解不适定问题的正则化方法结合起来,给出了离散Newton正则化方法的迭代格式,并给出了这种迭代格式的收敛性分析的结果。最后考虑了这种方法在微分方程反问题上的应用,数值计算结果表明了这种方法的有效性。     

Dynamic stability in parametric resonance of axially accelerating viscoelastic Timoshenko beams

This paper investigates dynamic stability of an axially accelerating viscoelastic beam undergoing parametric resonance. The effects of shear deformation and rotary inertia are taken into account by the Timoshenko thick beam theory. The beam material obeys the Kelvin model in which the material time derivative is used. The axial speed is characterized as a simple harmonic variation about the constant mean speed. The governing partial-differential equations are derived from Newton's second law, Euler's angular momentum principle, and the constitutive relation. The method of multiple scales is applied to the equations to establish the solvability conditions in summation and principal parametric resonances. The sufficient and necessary condition of the stability is derived from the Routh-Hurvitz criterion. Some numerical examples are presented to demonstrate the effects of related parameters on the stability boundaries.