Asymptotic expansions of Bessel, Anger and Weber transformations

In this paper we explore higher order numerical quadrature for the integrations of systems containing Bessel, Anger and Weber functions. The method is constructed by finding the approximate solution of the differential equation and truncating the asymptotic series. Numerical examples based on theoretical results are presented to demonstrate the efficiency of the method.     

Szegő–Lobatto quadrature rules

Gauss-type quadrature rules with one or two prescribed nodes are well known and are commonly referred to as Gauss–Radau and Gauss–Lobatto quadrature rules, respectively. Efficient algorithms are available for their computation. Szeg? quadrature rules are analogs of Gauss quadrature rules for the integration of periodic functions; they integrate exactly trigonometric polynomials of as high degree as possible. Szeg? quadrature rules have a free parameter, which can be used to prescribe one node. This paper discusses an analog of Gauss–Lobatto rules, i.e., Szeg? quadrature rules with two prescribed nodes. We refer to these rules as Szeg?–Lobatto rules. Their properties as well as numerical methods for their computation are discussed.     

Performance and numerical behavior of the second‐order scheme of precise time‐step integration for transient dynamic analysis

Spurious high‐frequency responses resulting from spatial discretization in time‐step algorithms for structural dynamic analysis have long been an issue of concern in the framework of traditional finite difference methods. Such algorithms should be not only numerically dissipative in a controllable manner, but also unconditionally stable so that the time‐step size can be governed solely by the accuracy requirement. In this article, the issue is considered in the framework of the second‐order scheme of the precise integration method (PIM). Taking the Newmark‐β method as a reference, the performance and numerical behavior of the second‐order PIM for elasto‐dynamic impact‐response problems are studied in detail. In this analysis, the differential quadrature method is used for spatial discretization. The effects of spatial discretization, numerical damping, and time step on solution accuracy are explored by analyzing longitudinal vibrations of a shock‐excited rod with rectangular, half‐triangular, and Heaviside step impact. Both the analysis and numerical tests show that under the framework of the PIM, the spatial discretization used here can provide a reasonable number of model types for any given error tolerance. In the analysis of dynamic response, an appropriate spatial discretization scheme for a given structure is usually required in order to obtain an accurate and meaningful numerical solution, especially for describing the fine details of traction responses with sharp changes. Under the framework of the PIM, the numerical damping that is often required in traditional integration schemes is found to be unnecessary, and there is no restriction on the size of time steps, because the PIM can usually produce results with machine‐like precision and is an unconditionally stable explicit method. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

黏弹性阻尼作用下轴向运动Timoshenko梁振动特性的研究

黏弹性阻尼一直是轴向运动系统的研究热点之一.以往研究轴向运动系统大都没有考虑黏弹性阻尼的影响.但在工程实际中, 存在黏弹性阻尼的轴向运动体系更为普遍.本文研究了黏弹性阻尼作用下轴向运动Timoshenko梁的振动特性.首先, 采用广义Hamilton原理给出了轴向运动黏弹性Timoshenko梁的动力学方程组和相应的简支边界条件.其次, 应用直接多尺度法得到了轴速和相关参数的对应关系, 给出了前两阶固有频率和衰减系数在黏弹性作用下的近似解析解.最后, 采用微分求积法分析了在有无黏弹性作用下前两阶固有频率和衰减系数随轴速的变化; 给出了前两阶固有频率和衰减系数在黏弹性作用下的近似数值解, 验证了近似解析解的有效性.结果表明: 随着轴速的增大, 梁的固有频率逐渐减小.梁的固有频率和衰减系数随着黏弹性系数的增大而逐渐减小, 其中衰减系数与黏弹性系数成正比关系, 黏弹性系数对第一阶衰减系数和固有频率的影响很小, 对第二阶衰减系数和固有频率的影响较大.      

稀薄流非线性模型方程离散速度坐标法有限差分解

从一般非线性Bo ltzm ann方程出发,发展并实现了一套适于大范围K nudsen数稀薄流问题数值模拟的统一算法。采用BGK模型和Shakov模型近似碰撞项,进而引入两个二速度无量纲简化分布函数,通过关于分子速度第三分量取矩积分,将三速度单一模型方程变换为二速度微分方程组。基于G auss-H erm ite积分公式和正交多项式G auss积分公式,借助离散速度坐标法消除简化模型方程对分子速度空间的连续依赖性,从相空间到物理空间得到一组带源项双曲守恒离散方程,并给出其显式和隐式二阶迎风TVD有限差分解。以二维圆柱A r气体超声速绕流算例,验证了数值算法的有效性,比较分析了漫反射和镜面反射两种气体分子壁面反射模型的计算结果。     

一类带弱奇异核非线性偏积分微分方程的全离散有限元

1引言我们将研究下面一类带弱奇异核非线性偏积分微分方程的数值解:u_t-▽·(a(u)▽u)-integral from n=0 to tβ(t-s)△u(s)ds=f(u),x∈Ω,t∈(?),(1.1) u(·,t)=0,x∈(?)Ω,t∈J,(1.2) u(·,0)=v(x),x∈Ω,(1.3)其中Ω为平面上的凸角域,J=(0,T],α和f为R上的光滑函数,满足0相似文献   

数字正交检波技术在主动探测中的应用研究

提出基于光学系统的“猫眼”效应,应用单频调制的连续激光,对扫描式红外热成像目标进行有效探测的新方法。分析了光学系统的后向反射特性,针对回波信号特征,提出利用回波信号的相位信息,采用直接采样的数字正交检波技术实现对低信噪比回波信号的可靠检测,并在外场试验中应用该方案对实际目标回波信号获得了满意的检测结果。     

变截面欧拉梁自由振动分析的重采样微分求积法

采用重采样微分求积法求解了变截面欧拉梁的自由振动问题。推导了变截面梁的控制方程离散格式,采用重采样矩阵方法对边界条件进行处理,给出了变截面梁自由振动算法。采用本文方法对不同类型截面形式和不同边界条件的变截面梁进行自由振动分析,并和其他解法进行比较。计算结果表明,本文方法可以适用于不同变截面类型和不同边界条件,计算精度与解析解吻合良好,具有良好的收敛性能。在同等精度条件下网格点数少于现有计算方法。重采样转换矩阵边界处理方法相比于传统边界处理方法具有更快的收敛性能。     

基于微分求积有限元法的双层微板系统振动特性研究

基于修正的偶应力理论和两变量精化的剪切变形理论,建立了由Winkler-Pasternak连续弹性夹层连接的双层微板系统的自由振动模型,着重推导了系统异步振动的运动微分方程和势能泛函。融合Gauss-Lobatto求积准则和微分求积准则构造了具有C¹连续性的微分求积有限元。通过与已有文献进行对比,验证了数值方法的有效性。详细讨论了各种因素对系统同步和异步振动特性的影响。结果表明,系统的自由振动特性对材料尺度参数、长宽比、长厚比以及边界条件呈现出依赖性;弹性夹层刚度仅对系统异步振动产生作用;随着模态阶次的增大,材料尺度参数和弹性夹层刚度对异步振动频率和模态的影响变得显著。     

Unified two-phase nonlocal formulation for vibration of functionally graded beams resting on nonlocal viscoelastic Winkler-Pasternak foundation

A nonlocal study of the vibration responses of functionally graded (FG) beams supported by a viscoelastic Winkler-Pasternak foundation is presented. The damping responses of both the Winkler and Pasternak layers of the foundation are considered in the formulation, which were not considered in most literature on this subject, and the bending deformation of the beams and the elastic and damping responses of the foundation as nonlocal by uniting the equivalently differential formulation of well-posed strain-driven (ε-D) and stress-driven (σ-D) two-phase local/nonlocal integral models with constitutive constraints are comprehensively considered, which can address both the stiffness softening and toughing effects due to scale reduction. The generalized differential quadrature method (GDQM) is used to solve the complex eigenvalue problem. After verifying the solution procedure, a series of benchmark results for the vibration frequency of different bounded FG beams supported by the foundation are obtained. Subsequently, the effects of the nonlocality of the foundation on the undamped/damping vibration frequency of the beams are examined.