非线性弹性杆中纵波传播过程的数值模拟

基于Hamilton空间体系的多辛理论研究了包含材料非线性效应和几何弥散效应的非线性弹性杆中纵波传播问题。导出了其Bridges意义下的多辛形式及其多种守恒律,并构造了等价于Box多辛格式的隐式多辛格式对不考虑材料非线性效应和几何弥散效应、只考虑材料非线性效应、只考虑几何弥散效应、同时考虑材料非线性效应和几何弥散效应四种情况下不同截面参数的圆杆中的纵波传播过程进行数值模拟,模拟结果不仅全面地反映了非线性效应和几何弥散效应对纵波传播的影响,而且也反映了多辛方法的两大优点:精确的保持多种守恒律和良好的长时间数值行为。     

辛体系下含对称破缺因素动力学系统的近似守恒律

自冯康先生创立Hamilton系统辛几何算法以来,诸如辛结构和能量守恒等守恒律逐渐成为动力学系统数值分析方法有效性的检验标准之一。然而,诸如阻尼耗散、外部激励与控制和变参数等对称破缺因素是实际力学系统本质特征,影响着系统的对称性与守恒量。因此,本文在辛体系下讨论含有对称破缺因素的动力学系统的近似守恒律。针对有限维随机激励Hamilton系统,讨论其辛结构;针对无限维非保守动力学系统、无限维变参数动力学系统、Hamilton函数时空依赖的无限维动力学系统和无限维随机激励动力学系统,重点讨论了对称破缺因素对系统局部动量耗散的影响。上述结果为含有对称破缺因素的动力学系统的辛分析方法奠定数学基础。     

Analyzing dynamic response of non-homogeneous string fixed at both ends

To survey the local conservation properties of complex dynamic problems, a structure-preserving numerical method, named as generalized multi-symplectic method, is proposed to analyze the dynamic response of a non-homogeneous string fixed at both ends. Firstly, based on the multi-symplectic idea, a generalized multi-symplectic form derived from the vibration equation of a non-homogeneous string is presented. Secondly, several conservation laws are deduced from the generalized multi-symplectic form to illustrate the local properties of the system. Thirdly, a centered box difference scheme satisfying the discrete local momentum conservation law exactly, named as a generalized multi-symplectic scheme, is constructed to analyze the dynamic response of the non-homogeneous string fixed at both ends. Finally, numerical experiments on the generalized multi-symplectic scheme are reported. The results illustrate the high accuracy, the good local conservation properties as well as the excellent long-time numerical behavior of the generalized multi-symplectic scheme well.     

On a “deterministic” explanation of the stochastic resonance phenomenon

The cantilevered carbon nanotube is a traditional model in the design of some precise nano-oscillators. The sensitivity that is denoted by the quality factor defined as the ratio of the energy stored to the energy dissipated by losses in the oscillator is an important index for the nano-oscillators. In this paper, the small mass impacting on the single-walled carbon nanotube is simplified as an impact load and the governing equation of the transverse oscillation for the nanotube is established firstly. Based on the structure-preserving idea, the generalized multi-symplectic formulations of the governing equation are constructed. The oscillation of the damping nanotube under the transverse impact load with different amplitudes is simulated after the small numerical dissipation and the good convergence of the scheme constructed is verified. From the numerical results, the effects of the induced tension and the energy dissipation are investigated. More importantly, the high precision and the feasibility of the numerical approach proposed in this paper for reproducing the energy dissipation in the carbon nanotube oscillator are illustrated, which gives a new way to investigate the properties of the carbon nanotube oscillators.     

Geometrical nonlinear waves in finite deformation elastic rods

IntroductionSomenewphenomenaofnonlinearwavesinthesolidmediumsuchasshockwave ,solitarywaveetc.arepaidmoreattentiontoincreasinglybyresearchersbecausetheytakeonalotofimportantproperties.ItistheoreticallyanalyzedinRefs.[1 -6]thattheformationmechanismsofshockwaveandsolitarywaveintheelasticthinrodsaswellastheirpropagationproperties.TheexistenceofsolitarywaveintheelasticmediumsuchasarodandaplatehasbeenverifiedinRef.[7]byexperiments.Shockwaveandsolitarywavearesteadilypropagatingtraveling_wavesgenerat…     

变截面弹性直杆纵振动分析的小波——DQ法

在经典微分求积(DQ)法基础上, 根据多分辨分析理论, 以尺度函数为基础构造插 值基函数, 形成了新的微分方程边值问题的求解方法------小波--DQ法, 并应用该方法分析了变截面弹性直杆的纵振动问题, 给出了其频率方程, 计算出了不同参数下固 支--固支, 自由--自由楔形直杆和锥形直杆的固有频率, 数值结果表明该方法是一个简单易行高精度的方法, 该方法可以推广应用于其他力学问题的数值分析.     

Longitudinal waves in one dimensional non-uniform waveguides

Wave approach is used to analyze the longitudinal wave motion in one dimensional non-uniform waveguides. With assumptions of constant wave velocity and no wave conversion, there exist four types of non-uniform rods and corresponding traveling wave solutions are investigated. The obtained results indicate that the kinetic energy is preserved as a constant and the wave amplitude is inversely proportional to square root of the cross-sectional area of the rod. Under certain condition, there exists a cut-off frequency for the rod with variation in geometric or material properties, below which waves do not propagate along the non-uniform rod. For the rod with arbitrary variable cross-section, the conclusions are similar if the wave frequency is high enough. And a series solution of the wave motion is presented.     

Multi-symplectic Runge-Kutta methods for Landau-Ginzburg-Higgs equation

Nonlinear wave equations have been extensively investigated in the last sev- eral decades. The Landau-Ginzburg-Higgs equation, a typical nonlinear wave equation, is studied in this paper based on the multi-symplectic theory in the Hamilton space. The multi-symplectic Runge-Kutta method is reviewed, and a semi-implicit scheme with certain discrete conservation laws is constructed to solve the first-order partial differential equations （PDEs） derived from the Landau-Ginzburg-Higgs equation. The numerical re- sults for the soliton solution of the Landau-Ginzburg-Higgs equation are reported, showing that the multi-symplectic Runge-Kutta method is an efficient algorithm with excellent long-time numerical behaviors.     

Wave-based defect detection and interwire friction modeling for overhead transmission lines

In this study, the feasibility of continuous, online monitoring of power lines using ultrasonic waves is considered. Local and global wave-based approaches for wire break detection in overhead transmission lines are presented. Both methods use a sending/receiving transducer to generate an ultrasonic, longitudinal, elastic wave in the cable. Defects in the cable cause a portion of the incident ultrasonic wave to be reflected back to the transducer, which when received, can be used to identify the presence of the defect. Although the transducers can only be attached to the surface of the cable, subsurface wires can also be interrogated since elastic energy spreads to these wires through friction contact. This study also explores how the elastic energy of a propagating wave becomes distributed among contacting rods via friction contact. This work focuses specifically on a two-rod system in which the wave energy from an excited “active” rod is transmitted to a neighboring “passive” rod through friction contact. An energy-based model is used to approximate the time average elastic wave power in the two rods as a function of propagation distance. Power predictions from the energy-based model compare well with experimental measurements and finite element simulations.     

Structure-preserving properties of three differential schemes for oscillator system

A numerical method for the Hamiltonian system is required to preserve some structure-preserving properties. The current structure-preserving method satisfies the requirements that a symplectic method can preserve the symplectic structure of a finite dimension Hamiltonian system, and a multi-symplectic method can preserve the multi-symplectic structure of an infinite dimension Hamiltonian system. In this paper,the structure-preserving properties of three differential schemes for an oscillator system are investigated in detail. Both the theoretical results and the numerical results show that the results obtained by the standard forward Euler scheme lost all the three geometric properties of the oscillator system, i.e., periodicity, boundedness, and total energy,the symplectic scheme can preserve the first two geometric properties of the oscillator system, and the St¨ormer-Verlet scheme can preserve the three geometric properties of the oscillator system well. In addition, the relative errors for the Hamiltonian function of the symplectic scheme increase with the increase in the step length, suggesting that the symplectic scheme possesses good structure-preserving properties only if the step length is small enough.     

Modeling and simulation of martensitic phase transitions with a triple point

A framework for modeling complex global energy landscapes in a piecewise manner is presented. Specifically, a class of strain-dependent energy functions is derived for the triple point of Zirconia (ZrO₂), where tetragonal, orthorhombic (orthoI) and monoclinic phases are stable. A simple two-dimensional framework is presented to deal with this symmetry breaking. An explicit energy is then fitted to the available elastic moduli of Zirconia in this two-dimensional setting. First, we use the orbit space method to deal with symmetry constraints in an easy way. Second, we introduce a modular (piecewise) approach to reproduce or model elastic moduli, energy barriers and other characteristics independently of each other in a sequence of local steps. This allows for more general results than the classical Landau theory (understood in the sense that the energy is a polynomial of invariant polynomials). The class of functions considered here is strictly larger. Finite-element simulations for the energy constructed here demonstrate the pattern formation in Zirconia at the triple point.     

Reflection and transmission of elastic waves at an elastic/porous solid saturated by two immiscible fluids

Wave propagation in a porous elastic medium saturated by two immiscible fluids is investigated. It is shown that there exist three dilatational waves and one transverse wave propagating with different velocities. It is found that the velocities of all the three longitudinal waves are influenced by the capillary pressure, while the velocity of transverse wave does not at all. The problem of reflection and refraction phenomena due to longitudinal and transverse wave incident obliquely at a plane interface between uniform elastic solid half-space and porous elastic half-space saturated by two immiscible fluids has been analyzed. The amplitude ratios of various reflected and refracted waves are found to be continuous functions of the angle of incidence. Expression of energy ratios of various reflected and refracted waves are derived in closed form. The amplitude ratios and energy ratios have been computed numerically for a particular model and the results obtained are depicted graphically. It is verified that during transmission there is no dissipation of energy at the interface. Some particular cases have also been reduced from the present formulation.     

Structure-preserving approach for infinite dimensional nonconservative system

The current structure-preserving theory, including the symplectic method and the multisymplectic method, pays most attention on the conservative properties of the continuous systems because that the conservative properties of the conservative systems can be formulated in the mathematical form. But, the nonconservative characteristics are the nature of the systems existing in engineering. In this letter, the structure-preserving approach for the infinite dimensional nonconservative systems is proposed based on the generalized multi-symplectic method to broaden the application fields of the current structure-preserving idea. In the numerical examples,two nonconservative factors, including the strong excitation on the string and the impact on the cantilever, are considered respectively. The vibrations of the string and the cantilever are investigated by the structure-preserving approach and the good long-time numerical behaviors as well as the high numerical precision of which are illustrated by the numerical results presented.     

Waves in Nonlocal Elastic Solid with Voids

In this paper, the governing relations and equations are derived for nonlocal elastic solid with voids. The propagation of time harmonic plane waves is investigated in an infinite nonlocal elastic solid material with voids. It has been found that three basic waves consisting of two sets of coupled longitudinal waves and one independent transverse wave may travel with distinct speeds. The sets of coupled waves are found to be dispersive, attenuating and influenced by the presence of voids and nonlocality parameters in the medium. The transverse wave is dispersive but non-attenuating, influenced by the nonlocality and independent of void parameters. Furthermore, the transverse wave is found to face critical frequency, while the coupled waves may face critical frequencies conditionally. Beyond each critical frequency, the respective wave is no more a propagating wave. Reflection phenomenon of an incident coupled longitudinal waves from stress-free boundary surface of a nonlocal elastic solid half-space with voids has also been studied. Using appropriate boundary conditions, the formulae for various reflection coefficients and their respective energy ratios are presented. For a particular model, the effects of non-locality and dissipation parameter (\(\tau \)) have been depicted on phase speeds and attenuation coefficients of propagating waves. The effect of nonlocality on reflection coefficients has also been observed and shown graphically.     

TRIAL FUNCTION METHOD FOR THE LONGITUDINAL VIBRATION OF RODS WITH VARIABLE CROSS-SECTIONS

给出了一种试探函数法,并研究了变截面杆的纵振动问题. 先给出振动控制方程的特殊函数形式的试探解,然后要求此解满足控制方程,反过来确定了控制方程各种可能的系数函数（即截面变化函数）并得到了控制方程的精确解. 作为例子,给出了一种变截面杆在3 种边界条件下的频率方程,计算出了固有频率. 研究表明,试探函数法简单、直接,适合于研究变截面杆的纵振动问题. 对于杆扭转振动、薄膜振动以及管中波传播等问题,该方法同样有推广应用价值.     

Variation in the length of perfectly elastic rods under torsion

On the basis of elastic constitutive relations that reflect geometrically nonlinear second-order effects, we refine the theory of torsion of rectilinear rods of an arbitrary transverse cross-section. In particular, we obtain a universal formula, independent of the material properties, that determines the longitudinal strain arising as the rod undergoes free torsion. According to this formula, the length of a rod made of an isotropic perfectly elastic material can, in contrast to the traditional concepts, either increase or decrease as the rod undergoes torsion. Moreover, the variation in the length depends only on the geometry of the transverse cross-section.     

On thermoelastic dielectrics

Constitutive equations for a linear thermoelastic dielectric are derived from the energy balance equation assuming dependence of the stored energy function on the strain tensor, the polarization vector, the polarization gradient tensor and entropy. A method is indicated for constructing a hierarchy of constitutive equations for materials with arbitrary symmetry by introducing various thermodynamic potentials. Maxwell's relations are constructed for the thermodynamic potential W^L. The entropy inequality is used to obtain stability conditions for an elastic dielectric in equilibrium under prescribed boundary constraints. Frequencies are explicitly determined for a plane wave propagating along the x₁-axis in an infinite centro-symmetric isotropic thermoelastic dielectric.     

A linear theory of straight elastic rods

This paper is based on the work of Green & Laws who have given a general thermodynamical theory of rods which is valid for any material. Here, starting with the general non-linear theory of elastic rods, we derive a linear theory allowing for thermal effects. The resulting free energy as a quadratic function of kinematic variables is restricted by certain symmetry conditions. The basic equations then separate into four groups, two for flexure, one for torsion and one for extension of the rod with temperature effects occurring only in the latter group. Wave propagation along an infinite rod is considered. There are two wave speeds for each type of flexure, two for torsion and three for isothermal extension and all wave speeds depend on the wave length.