A fractional differential constitutive model for dynamic stress intensity factors of an anti-plane crack in viscoelastic materials

Fractional differential constitutive relationships are introduced to depict the history of dynamic stress inten- sity factors （DSIFs） for a semi-infinite crack in infinite viscoelastic material subjected to anti-plane shear impact load. The basic equations which govern the anti-plane deformation behavior are converted to a fractional wave-like equation. By utilizing Laplace and Fourier integral transforms, the fractional wave-like equation is cast into an ordinary differential equation （ODE）. The unknown function in the solution of ODE is obtained by applying Fourier transform directly to the boundary conditions of fractional wave-like equation in Laplace domain instead of solving dual integral equations. Analytical solutions of DSIFs in Laplace domain are derived by Wiener-Hopf technique and the numerical solutions of DSIFs in time domain are obtained by Talbot algorithm. The effects of four parameters α, β, b1, b2 of the fractional dif- ferential constitutive model on DSIFs are discussed. The numerical results show that the present fractional differential constitutive model can well describe the behavior of DSIFs of anti-plane fracture in viscoelastic materials, and the model is also compatible with solutions of DSIFs of anti-plane fracture in elastic materials.     

Solitary wave solutions of the nonlinear generalized Pochhammer–Chree and regularized long wave equations

In this paper, we implement some fast and high accuracy numerical algorithms to obtain the solitary wave solutions of generalized Pochhammer?CChree (PC) and regularized long wave (RLW) equations. We employ the discrete Fourier transform to discretize the original partial differential equations (PDEs) in space and obtain a system of ordinary differential equations (ODEs) in Fourier space which will be solved with fourth order time-stepping methods. The proposed methods are fast and accurate due to the use of the fast Fourier transform in combination with explicit fourth-order time stepping methods. For RLW equation we investigate the propagation of a single solitary and interaction of two and three solitary waves. Moreover, three invariants of motion (mass, energy, and momentum) are evaluated to determine the conservation properties of the problem, and the numerical schemes lead to accurate results. The numerical results are compared with analytical solutions and with those of other recently published methods to confirm the accuracy and efficiency of the presented schemes.     

Higher order dispersive effects in regularized Boussinesq equation

In this paper, we consider the higher order Boussinesq (HBq) equation which models the bi-directional propagation of longitudinal waves in various continuous media. The equation contains the higher order effects of frequency dispersion. The present study is devoted to the numerical investigation of the HBq equation. For this aim a numerical scheme combining the Fourier pseudo-spectral method in space and a Runge–Kutta method in time is constructed. The convergence of semi-discrete scheme is proved in an appropriate Sobolev space. To investigate the higher order dispersive effects and nonlinear effects on the solutions of HBq equation, propagation of single solitary wave, head-on collision of solitary waves and blow-up solutions are considered.     

层状饱和土Biot固结问题状态空间法

针对饱和多孔介质空间非轴对Biot固结问题，引入状态变量，构造了两组相比独立的状态变量方程，利用Fourier级数和Laplace-Hankel变换，将状态变量方程转换为两组一阶常微分方程组，提出了均质饱和多孔介质空间非轴对称Biot固结问题的传递矩阵，得到以状态变量和传递矩阵乘积的形式表示的均质饱和多孔介质空间非轴对称Biot固结问题的解，利用层间完全接触的条件，可得到N层饱和多孔介质空间非轴对称Biot固结问题的一般解析表达式，文中考虑几种不同的边界条件，分析了两个算例，数值结果表明该方法具有较高的计算精度和良好的计算稳定性。     

结构粘弹性分析的线法

本文将以现代常微分方程求解器为支撑软件的线法推广应用于含时间变量的粘弹性力学问题。基基本思想是先用拦普拉斯积分变换把问题的控制偏微分方程组和定解条件进行交换,得到拉氏变换空间中的定常问题,然后用传统的线法借助于常微分方程求解器得到该变换空间中的数值解,再应用拉氏数值逆变换求得原时空域中的数值解。文中还给出方法的基本原理和实施过程。算例表明,采用拉普拉斯积分变换于线法分析结构粘弹性力学等含时间变量的     

Lubrication theory for micropolar fluids and its application to a journal bearing with finite length

In this paper, the field equation of micropolar fluid with general lubrication theory assumptions is simplified into two systems of coupled ordinary differential equation. The analytical solutions of velocity and microrotation velocity are obtained. Micropolar fluid lubrication Reynolds equation is deduced. By means of numerical method, the characteristics of a finitely long journal bearing under various dynamic parameters, geometrical parameters and micropolar parameters are shown in curve form. These characteristics are pressure distribution, load capacity, coefficient of flow flux and coefficient of friction. Practical value of micropolar effects is shown, so micropolar fluid theory further closes to engineering application.     

双周期加筋微穿孔板的振动响应及声透射的研究

本文主要研究了水下无穷大双周期加筋微穿孔薄板,在平面声波斜入射下的振动响应和声透射,并提出了一种半解析半数值的计算方法。利用微穿孔板的声阻抗以及薄板表面的振速边界条件,建立了加筋穿孔薄板的振动方程,并根据傅立叶变换及空间波数法将振动位移表达为波数分量的迭加形式。采用数值计算的方法对波数分量进行求解并通过傅里叶逆变换,最终得到了双周期加筋穿孔薄板的振动响应及透射系数。通过与Takahashi穿孔板声压结果的对比,证明了本方法的正确性。在算例中,分析了加强筋及穿孔率对薄板结构的振动和声透射的影响。     

Analysis of bending waves in beam on viscoelastic random foundation using wavelet technique

This paper deals with the mathematical model of dynamic behaviour of the beam resting on viscoelastic random foundation. It is considered by assuming the modulus of subgrade reaction to be a homogeneous random function of space variable. The problem is governed by the fourth-order differential equation with random parameters. The main results of this article are the approximate analytical solutions for the displacement field, variance and dynamic-stiffness coefficient. It has been made a comparison of numerical results obtained by using two different methods: Adomian’s decomposition and Bourret’s approximation. The special method of finding inverse Laplace transform based on the wavelet theory is adopted and used in numerical examples. For making numerical calculations and plots the programs in MATHEMATICA have been prepared.     

Nonlinear vibration of a two-mass system with nonlinear stiffnesses

An analytical approach is developed for nonlinear free vibration of a conservative, two-degree-of-freedom mass–spring system having linear and nonlinear stiffnesses. The main contribution of the proposed approach is twofold. First, it introduces the transformation of two nonlinear differential equations of a two-mass system using suitable intermediate variables into a single nonlinear differential equation and, more significantly, the treatment a nonlinear differential system by linearization coupled with Newton’s method and harmonic balance method. New and accurate higher-order analytical approximate solutions for the nonlinear system are established. After solving the nonlinear differential equation, the displacement of two-mass system can be obtained directly from the governing linear second-order differential equation. Unlike the common perturbation method, this higher-order Newton–harmonic balance (NHB) method is valid for weak as well as strong nonlinear oscillation systems. On the other hand, the new approach yields simple approximate analytical expressions valid for small as well as large amplitudes of oscillation unlike the classical harmonic balance method which results in complicated algebraic equations requiring further numerical analysis. In short, this new approach yields extended scope of applicability, simplicity, flexibility in application, and avoidance of complicated numerical integration as compared to the previous approaches such as the perturbation and the classical harmonic balance methods. Two examples of nonlinear two-degree-of-freedom mass–spring system are analyzed and verified with published result, exact solutions and numerical integration data.     

基于Fourier级数的时变周期系数Riccati微分方程精细积分

结合Fourier级数展开方法,本文提出了基于精细积分的时变周期系数Riccati微分方程求解高效算法.首先,利用Fourier级数展开方法将周期系统表示成三角级数形式,在一个积分步内使用精细积分方法得到对应Hamilton系统状态转移矩阵的表达式.然后,通过Riccati变换的方法,得到含有状态转移矩阵的时变周期系数Riccati微分方程解的递推格式.本文方法充分利用了方程本身的周期性特点,文中的数值算例表明算法具有计算效率高、结果可靠等优势.     

粘弹性固体的精细积分有限元算法

粘弹性固体本构方程的数学表达式分为微分型和积分型两种，其数值求解主要是时域上离散计算。文中从微分型表达式出发导出其状态空间方程的数学表达式，通过严格推导论证了它与微、积分型表达式的等价性；引入状态空间方程，从而利用精细积分格式来求解粘弹性固体本构方程；给出了粘弹性固体本构方程的精细积分有限元算法，为求解粘弹性固体本构方程的数值解提供了一个新的途径，具有计算简便，求解精度高等优点。     

Application of a modified Lindstedt–Poincaré method in coupled TDOF systems with quadratic nonlinearity and a constant external excitation

An analytical approach is developed for the nonlinear oscillation of a conservative, two-degree-of-freedom (TDOF) mass-spring system with serial combined linear–nonlinear stiffness excited by a constant external force. The main idea of the proposed approach lies in two categories, the first one is the transformation of two nonlinear differential equations of a two-mass system using suitable intermediate variables into a single nonlinear differential equation. Another is the treatment a quadratic nonlinear oscillator (QNO) by the modified Lindstedt–Poincaré (L-P) method presented recently by the authors. The first-order and second-order analytical approximations for the modified L-P method are established for the QNOs with satisfactory results. After solving the nonlinear differential equation, the displacements of two-mass system can be obtained directly from the governing linear second-order differential equation. Unlike the common perturbation method, the modified L-P method is valid for weak as well as strong nonlinear oscillation systems. On the other hand, the new approach yields simple approximate analytical expressions valid for small as well as large amplitudes of oscillation. In short, this new approach yields extended scope of applicability, simplicity, flexibility in application, and avoidance of complicated numerical integration as compared to the previous approaches such as the perturbation and classical harmonic balance methods. Two examples of nonlinear TDOF mass-spring systems excited by a constant external force are selected and the approximate solutions are verified with the exact solutions derived from the Jacobi elliptic function and also the numerical fourth-order Runge–Kutta solutions.     

Normal forms for NFDEs with parameters and application to the lossless transmission line

A method for the computation of normal forms for neutral functional differential equations (NFDEs) with parameters is developed by considering an extension of phase space, based on the method of computing normal forms for FDEs with parameters previously introduced by Faria. The Hopf bifurcation of the differential difference equation is considered as an example of a circuit involving a lossless transmission line. The direction and stability of the bifurcating periodic solutions are also determined. Finally, numerical simulations are carried out to support the analytic results. This research is supported by the NNSF of China.     

A Numerical Method for Wave Propagation in Viscoelastic Stratified Porous Media

This paper presents a numerical method, a transmission matrix method, for the wave propagation in viscoelastic stratified saturated porous media. The wave propagation in saturated media, based on Biot theory, is a coupled problem. In this stratified three-dimensional model we do the Laplace transform for the time variable and the Fourier transform for the horizontal space coordinate. The original problem is transformed into ordinary differential equations with six independent unknown variables, which are only the function of the coordinate of depth. Thus, we get a transmission matrix of the wave problem for each layer. In the process of solution we use numerical method to calculate the eigenvalues and the eigenvectors of the transmission matrices. In the first step of the solution process we can obtain the wave field in the transformed space. The fast Fourier transform (FFT) method is used to do the inverse Laplace and the inverse Fourier transforms to get the solution in the time space. The detailed formulae are derived and some numerical examples are given.     

强非线性动力系统的频率增量法

提出一类强非线性动力系统的暧时频率增量法,将描述动力系统的二阶常微分方程,化为以相位为自变量、瞬廛频率为未知函数的积分方程;用谐波平衡原理,将求解瞬时频率的积分问题,归结为求解以频率增量的Fourier系数为独立变量的线性代数方程组;给出了若干例子。     

An analytical symplectic approach to the vibration analysis of orthotropic graphene sheets

A nonlocal continuum orthotropic plate model is proposed to study the vibration behavior of single-layer graphene sheets (SLGSs) using an analytical symplectic approach.A Hamiltonian system is established by introduc-ing a total unknown vector consisting of the displacement amplitude,rotation angle,shear force,and bending moment. The high-order governing differential equation of the vibra-tion of SLGSs is transformed into a set of ordinary differential equations in symplectic space.Exact solutions for free vibra-tion are obtianed by the method of separation of variables without any trial shape functions and can be expanded in series of symplectic eigenfunctions. Analytical frequency equations are derived for all six possible boundary con-ditions. Vibration modes are expressed in terms of the symplectic eigenfunctions.In the numerical examples,com-parison is presented to verify the accuracy of the proposed method. Comprehensive numerical examples for graphene sheets with Levy-type boundary conditions are given.A para-metric study of the natural frequency is also included.     

Non-axisymmetrical vibration of elastic circular plate on layered transversely isotropic saturated ground

The non-axisymmetrical vibration of elastic circular plate resting on a layered transversely isotropic saturated ground was studied.First,the 3-d dynamic equations in cylindrical coordinate for transversely isotropic saturated soils were transformed into a group of governing differential equations with 1-order by the technique of Fourier ex- panding with respect to azimuth,and the state equation is established by Hankel integral transform method,furthermore the transfer matrixes within layered media are derived based on the solutions of the state equation.Secondly,by the transfer matrixes,the general solutions of dynamic response for layered transversely isotropic saturated ground excited by an arbitrary harmonic force were established under the boundary conditions, drainage conditions on the surface of.ground as well as the contact conditions.Thirdly, the problem was led to a pair of dual integral equations describing the mixed boundary- value problem which can be reduced to the Fredholm integral equations of the second kind solved by numerical procedure easily.At the end of this paper,a numerical result concerning vertical and radical displacements both the surface of saturated ground and plate is evaluated.     

Approximation of generalized bounded solutions of evolution equations with unbounded operator

We substantiate a parametrization method for a differential equation in a Banach space with an unbounded operator coefficient. We propose an algorithm for finding bounded generalized solutions with an arbitrary order of accuracy.     

Problem of periodic solutions for neutral functional differential equation

ＩｎｔｒｏｄｕｃｔｉｏｎＬｅｔＣ(ｋ- 1)2π =ｈ(ｔ) ｜ｈ :Ｒ →Ｒｉｓ (ｋ -1 )＿ｔｈｏｒｄｅｒｃｏｎｔｉｎｕｏｕｓｄｉｆｆｅｒｅｎｔｉａｂｌｅａｎｄｈ(ｔ+ 2π) ≡ｈ(ｔ) ,　　Ｃ2π =ｈ(ｔ) ｜ｈ :Ｒ →Ｒｉｓｃｏｎｔｉｎｕｏｕｓａｎｄｈ(ｔ+ 2π) ≡ｈ(ｔ) ,　　‖ｈ(ｔ)‖ =ｓｕｐｔ∈ [0 ,2π] ｜ｈ(ｔ) ｜ ,　　‖ｈ(ｔ)‖Ｐｋ- 1 =ｍａｘ‖ｈ(ｔ)‖ ,‖ｈ′(ｔ)‖ ,… ,‖ｈ(ｋ- 1) (ｔ)‖ ,　　ｘ(ｍ) (ｔ+ ·) (θ) =ｘ(ｍ) (ｔ+θ)　　θ∈Ｒ (ｍ =0 ,1 ,2 ,… ,ｋ-1 ) .Ｃｌｅａｒｌｙ ,ｘ(ｍ) (ｔ + ·) ∈Ｃ2π, …     

Non-linear oscillations of a third-order differential equation

An investigation is conducted into the behavior of the solutions of a third-order non-linear differential equation which is characterized by a non-linearity depending solely upon the Euclidean norm of the associated phase space. The non-linearity represents a central restoring force, which has important applications in modern control theory. For small non-linearities, the existence of a limit cycle is established by a fixed point technique, the approach to the limit cycle is approximated by averaging methods, and the periodic solution is harmonically represented by perturbation. Computer solutions of the differential equation are provided in order to reinforce the analysis. Some related differential equations are discussed including one in which the periodic solution is explicitly prescribed.