轴向变速运动粘弹性弦线横向振动的复模态Galerkin方法

在考虑初始张力和轴向速度简谐涨落的情况下,利用含预应力三维变形体的运动方程,建立了轴向变速运动弦线横向振动的非线性控制方程,材料的粘弹性行为由Kelvin模型描述．利用匀速运动线性弦线的模态函数构造了变速运动非线性弦线复模态Galerkin方法的基底函数,并借助构造出来的基底函数研究了复模态Galerkin方法在轴向变速运动粘弹性弦线非线性振动分析中的应用．数值结果表明,复模态Galerkin方法相比实模态Galerkin方法对变系数陀螺系统有较高的收敛速度．     

轴向运动粘弹性板的横向振动特性

研究了轴向运动粘弹性矩形薄板的动力特性和稳定性问题．从二维粘弹性微分型本构关系出发,建立了轴向运动粘弹性板的运动微分方程．采用微分求积法,对四边简支、一对边简支一对边固支两种边界条件下粘弹性板的无量纲复频率进行了数值计算．分析了薄板的长宽比、无量纲运动速度及材料的无量纲延滞时间对其横向振动及稳定性的影响．     

非线性粘性柱的稳定性和混沌运动

研究了受轴向周期力作用的各向同性简支柱的动力学稳定性。假定粘弹性材料满足Ｌｅａｄｅｒｍａｎ非线性本构关系。导出运动方程为非线性微分－积分方程,并利用Ｇａｌｅｒｋｉｎ方法简化为非线性微分－积分方程。应用平均法进行了稳定性分析,并用数值结果进行验证。数值结果还表明系统可能存在混沌运动。     

小垂度粘弹性索非线性响应及振动主动控制

研究具有初始应力的小垂度粘弹性索的非线性动态响应及振动主动控制。在假定索材料的本构关系为一般微分本构类型的基础上,建立小垂度粘弹性索的运动微分方程;应用Galerkin方法将其转化为可用Runge-Kutta数值积分方法求解的一系列三阶非线性常微分方程。在仅考虑面内的横向振动及忽略非线性的情况下得到了连续状态空间中的状态方程,将状态方程离散为差分方程形式,并用矩阵指数来逐步近似状态转移矩阵;基于二次性能指标的最小化得到了最优的控制力与状态向量。最后通过数值仿真研究说明了粘性参数对索动态响应的影响。     

非线性粘弹性柱的稳定性和混沌运动

研究了受轴向周期力作用的各向同性简支柱的动力学稳定性。假定粘弹性材料满足Lea-derman非线性本构关系。导出运动方程为非线性偏微分-积分方程,并利用Galerkin方法简化为非线性微分-积分方程。应用平均法进行了稳定性分析,并用数值结果进行验证。数值结果还表明系统可能存在混沌运动。     

轴向变速运动弦线的非线性振动的稳态响应及其稳定性

研究具有几何非线性的轴向运动弦线的稳态横向振动及其稳定性．轴向运动速度为常平均速度与小简谐涨落的叠加．应用Hamilton原理导出了描述弦线横向振动的非线性偏微分方程．直接应用于多尺度方法求解该方程．建立了避免出现长期项的可解性条件．得到了近倍频共振时非平凡稳态响应及其存在条件．给出数值例子说明了平均轴向速度、轴向速度涨落的幅值和频率的影响．应用Liapunov线性化稳定性理论,导出倍频参数共振时平凡解和非平凡解的不稳定条件．给出数值算例说明相关参数对不稳定条件的影响．     

具有分数导数型本构关系的粘弹性柱的动力稳定性

研究简支的受轴向周期激励的粘弹性柱动力稳定性,柱的材料满足分数导数型本构关系.建立了描述粘弹性柱动力学行为的弱奇异性Volterra积分-偏微分方程,利用Galerkin方法将其化归为弱奇异性Volterra积分-常微分方程.利用平均化方法的思想给出了粘弹性柱运动稳定状态的存在性条件.给出一种新的计算方法,克服了存储整个响应历史数据的困难,并给出了数值算例,计算结果与解析方法的结论比较吻合.     

粘弹性运动带动力响应分析

基于Kelvin粘弹性材料本构模型及带运动方程,建立了运动带非线性动力学分析模型.基于该模型和Lie群分析方法推导了匀速运动及简谐运动带线性问题的解析解;基于该非线性模型的数值仿真讨论了运动带材料参数、带稳态运动速度、扰动速度对系统动态响应的影响.结果表明:1)当带匀速运动时,无论系统是线性还是非线性,运动带横向振动"频率"都随着带运动稳态速度增加而减小.2)随着材料粘性增加,系统耗散能力逐渐增强,动态响应逐渐减小.3)当带运动速度简谐波动时,系统动态响应随扰动速度增大而增大.扰动频率对带横向振动影响较大.     

服从OldroydB型微分模型的粘弹性流体问题的数值解法

本文对服从OldroydB型微分模型的粘弹性流体问题给出了一种数值逼近算法．该算法对压力方程采用标准混合有限元方法,对速度方程采用并行非重叠区域分解方法和特征线法．这种并行算法在子区域上用Galerkin方法,通过积分平均方法显式地给出内边界的数值流．在本文最后还给出了该算法的最优L^2。一误差估计．     

分数积分的一种数值计算方法及其应用

提出了一种只需要存储部分历史数据的分数积分的数值计算方法,并给出了误差估计。这种方法可对包含分数积分和分数导数的积分-微分方程进行较长时间的数值计算,克服了存储全部历史数据的困难,并能对计算误差进行控制。作为应用,给出了具有分数导数型本构关系的粘弹性Timoshenko梁的动力学行为研究的控制方程,利用分离变量法讨论梁在简谐激励作用下的动力响应,然后用新提出的数值方法对控制方程进行数值计算,数值计算结果和理论结果进行了比较,它们比较吻合。     

Nonlinear evolution of the travelling waves at the surface of a thin viscoelastic falling film

The problem of hydrodynamic instability of a thin condensate viscoelastic liquid film flowing down on the outer surface of an axially moving vertical cylinder is investigated. In order to improve the accuracy of numerical results, the viscoelastic and heat transfer parameters have been included into the governing equations. Also, the analytical solutions are obtained by utilizing the long-wave perturbation method. The influence of some physical parameters is discussed in both linear and nonlinear steps of the problem. It has been revealed that the stability of the film flow is weakened when the radius of cylinder and the temperature difference are reduced. Moreover, it is found that the increment of down-moving motion of the cylinder can enhance the flow stability. Further, the thin film flow can be destabilized by the viscoelastic property. The results show that both supercritical stability and subcritical instability can take place within the film flow system given appropriate conditions. Moreover, the absence of Reynolds number leads to an obvious difference in the behavior of some physical parameters.     

An upwind center difference parallel method and numerical analysis for the displacement problem with moving boundary

A nonlinear coupled mathematical system of two‐phase seepage flow displacement is discussed in this paper including an elliptic equation for the pressure and a convection‐dominated diffusion equation for the saturation. In fact, the boundary of an underground region where the fluid flows through is nonstationary. So a moving boundary should be considered. The saturation equation is convection‐dominated, therefore the method of upwind finite difference is introduced for the accurate computation. The upwind approximation could eliminate numerical oscillation and strong stability is shown. Since the computational work of saturation is larger than the pressure, the authors apply a parallel method, decomposing the whole domain into several nonoverlapping subdomains, to simplify the computation. A domain decomposition method coupled with upwind differences is presented for the saturation. The pressure equation is discretized by a five‐point center finite difference method. By using a transformation and defining new inner products and norms, error estimates in l² norm is discussed. Finally, two experimental tests are given to illustrate the efficiency and accuracy of the parallel algorithm.     

A Mildly Non-linear System of Parabolic Differential Equations Arising from Mass Transfer with Chemical Reaction in Electrochemistry

A system of three connected parabolic equations is studied.The first equation is the straight forward diffusion equationin one space dimension and the solution can be written down.The remaining two cannot be solved analytically but it is interestingto observe that a solution does exist for their difference.By considering the problem as a moving boundary value probleman approximate solution is obtained by a finite difference technique.An analysis of stability is performed and numerical resultsfor a specific chemical reaction are presented.     

The Conservation and Convergence of Two Finite Difference Schemes for KdV Equations with Initial and Boundary Value Conditions

Korteweg-de Vries equation is a nonlinear evolutionary partial differential equation that is of third order in space. For the approximation to this equation with the initial and boundary value conditions using the finite difference method, the difficulty is how to construct matched finite difference schemes at all the inner grid points. In this paper, two finite difference schemes are constructed for the problem. The accuracy is second-order in time and first-order in space. The first scheme is a two-level nonlinear implicit finite difference scheme and the second one is a three-level linearized finite difference scheme. The Browder fixed point theorem is used to prove the existence of the nonlinear implicit finite difference scheme. The conservation, boundedness, stability, convergence of these schemes are discussed and analyzed by the energy method together with other techniques. The two-level nonlinear finite difference scheme is proved to be unconditionally convergent and the three-level linearized one is proved to be conditionally convergent. Some numerical examples illustrate the efficiency of the proposed finite difference schemes.     

Non-linearly parametric resonances of an axially moving viscoelastic sandwich beam with time-dependent velocity

Non-linearly parametric resonances of an axially moving viscoelastic sandwich beam are investigated in this paper. The beam is moving with a time-dependent velocity, namely a harmonically varied velocity about the mean velocity. The partial differential equation is discretized into nonlinear ordinary differential equations via the method of Galerkin truncation and then the steady-state response is obtained using the method of multiple scales, an approximate analytical method. The tuning equations are obtained by eliminating secular terms and the amplitude of the vibration is derived from the tuning equations expressed in polar form, and two bifurcation points are obtained as well. Additionally, the stability conditions of trivial and nontrivial solutions are analyzed using the Routh–Hurwitz criterion. Eventually, the effects of various parameters such as the thickness of core layer, mean velocity, initial tension, and the amplitude of axially moving velocity on amplitude–frequency response curves and unstable regions are investigated.     

轴向运动三参数黏弹性梁弱受迫振动的渐近分析

研究了轴向运动三参数黏弹性梁的弱受迫振动．建立了轴向运动三参数黏弹性梁受迫振动的控制方程．使用多尺度法渐近分析了运动梁的稳态响应,导出了解稳定性边界方程、稳态振幅的表达式以及稳态响应非零解的存在条件．依据Routh-Hurwitz定律决定了非线性稳态响应非零解的稳定性．