用切比雪夫多项式研究短梁的弯曲计算

考虑剪切效应,利用切比雪夫多项式构造严格满足表面切应力边界条件的轴向位移表达式,建立了短梁弯曲问题的新理论．利用奇异函数把作用在短梁上的复杂外载荷表示为分布载荷,推导出了短梁弯曲时的截面正应力公式及挠曲线表达式．把采用切比雪夫多项式推导出短梁的弯曲计算公式计算结果与弹性理论计算结果进行比较,可知该方法的计算精度较高．研究结果表明：在复杂外载荷作用下,当长高比小于等于6时,剪切变形对梁的弯曲挠度影响较大,而当长高比小于3时,剪切变形对梁的弯曲应力影响较大;因此建议采用切比雪夫多项式方法给出的挠度表达式、弯曲应力进行计算,因为切比雪夫多项式方法不但给出了复杂外载荷作用下梁截面挠度、弯曲应力的计算通式,而且该方法具有计算过程简便、精度高的优点．     

带功能梯度材料的压电底层中周期裂纹对SH波的散射

本文研究了压电材料底层中周期裂纹对SH波的散射,通过渗透边界条件和界面上连续边界条件,将问题转化为一组带Hilbter核的奇异积分方程。利用利用切比雪夫多项式逼近方法求解Hilbter核的奇异积分方程,给出了标准动应力强度因子和电位移强度因子的表达式。最后通过数值算例说明了几何参数、物性参数,入射波频率和振幅等对强度因子的影响.     

岩土力学参数概率分布的切比雪夫多项式推断

提出了较大样本岩土力学参数概率分布的切比雪夫多项式逼近法。基于数值逼近原理,直接根据试验样本矩,运用切比雪夫多项式推断岩土力学参数的概率密度函数,并用精度较高的K—S检验法,从理论上证明所求密度函数的正确性和实用性。该方法直接根据试验样本信息和统计方法推断,而不是事先假定成经典的理论概率分布,因此数学和物理意义更加充分。通过对各种经典分布曲线（正态分布、指数分布、对数正态分布等）数值检验,结果表明所得到的逼近表达式有很好的拟合性能。根据样本数据得出的某岩石抗压强度概率密度函数,与实际统计所得分布频率非常接近,可以满足岩土工程可靠性分析的要求。     

基于动态响应参数的点阵桁架夹芯板脱焊损伤检测数值方法研究

点阵桁架夹芯板作为一种特殊的超轻多孔材料,具有广阔应用前景。针对其加工和服役过程中经常发生的脱焊损伤,本文提出一种基于结构振动特性进行金字塔型点阵桁架夹芯板无损评价的方法。首先,基于有限元理论对夹芯板在不同损伤状态(损伤个数、位置)下的振动特性进行数值模拟;然后基于计算所得振动特性参数,运用切比雪夫多项式逼近法计算相应均布载荷曲面的曲率,进而通过此参数与原始未损伤模型参数作差进行损伤识别;最后运用曲面拟合算法解决该检测方法对原始未损伤模型参数的依赖性问题。数值计算结果表明:四周简支金字塔型点阵桁架夹芯板中存在损伤时,损伤前后结构损伤位置处的均布载荷曲面的曲率差值比较明显,基于此差分方法可精确识别各种损伤状况下夹芯板脱焊损伤的位置和数量;利用曲面拟合算法可以消除损伤结构参数中因损伤造成的尖点,从而不需以原始未损伤模型参数为依据即可精确识别损伤。     

微分求积法在计算功能梯度Timoshenko梁临界荷载中的应用研究

提出了一种改进型等比数列布点方式研究轴向功能梯度Timoshenko变截面梁的屈曲临界荷载。首先基于Timoshenko梁理论,建立了求解功能梯度材料Timoshenko变截面梁屈曲临界荷载的变系数常微分方程,由微分求积法理论将其转化为标准型的广义代数特征值问题,再采用QR法求解该代数特征方程组,可一次性地计算出轴向功能梯度Timoshenko变截面梁的屈曲临界荷载。数值计算结果表明:当梁上区间单元划分段数N取28时,采用改进型等比数列布点方式和切比雪夫多项式根布点方式时,由微分求积法(DQM)获得的屈曲临界荷载数值解计算精度等价且与实际值完全吻合,证明了本文方法的可行性和计算精度;当N取8和12时,采用切比雪夫多项式根作为布点方式的计算值与实际值误差较大甚至失真,而采用等比数列变步长布点方式时,公比q为控制算法精度的控制参数,通过调整公比q可获得精确值,相对于切比雪夫多项式根作为布点方式,这一优势十分明显。     

基于Chebyshev展开的区间穿孔板超材料分析

目前对于声学超材料的传输特性分析和优化大多是基于确定的数值和确定的模型,然而在实际工程和结构设计中存在大量材料自身特性和几何物理参数的不确定性.如果忽略这些不确定变量对声学超材料传输特性分析和优化过程的影响,得到的结果可能不正确.针对这一现状,拟将切比雪夫区间模型引入多层穿孔板超材料,提出多层穿孔板超材料声学透射率的区间切比雪夫展开-蒙特卡洛模拟法(interval Chebyshev expansionMonte Carlo simulation method,ICE-MCSM).该方法采用截断切比雪夫多项式近似拟合多层穿孔板超材料的声学透射率响应曲线,构造声学透射率响应曲线的切比雪夫代理模型;然后采用蒙特卡洛模拟法(Monte Carlo simulation method,MCSM)随机生成一定数量的不确定区间变量的样本数据点,并将生成的不确定区间变量样本数据点代入切比雪夫代理模型,预测单个不确定区间变量和多个不确定区间变量条件下的多层穿孔板超材料声学透射率区间的上界和下界.数值分析结果表明,ICE-MCSM预测的声学透射率变化区间的上界和下界与直接蒙特卡洛法(direct Monte Carlo simulation method,DMCSM)预测的声学透射率上界和下界的结果非常接近.与DMCSM相比,ICE-MCSM具有更高的计算效率.因此,ICE-MCSM可有效且高效地分析不确定区间变量条件下多层穿孔板超材料声学透射率传输特性,具有良好的工程应用前景.     

求非线性动力系统周期解的切比雪夫多项式法

周期运动是一种在客观世界中普遍存在的运动形式,它与混沌运动之间存在十分密切的关系,因而具有很重要的研究价值。利用切比雪夫多项式的若干良好性质,对自治非线性动力系统进行分析,将状态矢量在主周期上展开为切比雪夫多项式的形式,从而将原问题转变为非线性代数方程组的求解问题,得出一种可以方便、迅速地获得周期轨道近似多项式表达式的方法。此方法不依赖于小参数假设,可以用于分析强非线性问题,而且对参数激励系统同样有效。在计算机条件允许时,对高维系统也能迅速、精确地得到其周期轨道的近似多项式表达式。以三维Rossler系统和五维非线性磁浮转子系统周期轨道的计算为例,通过与四阶Runge-Kutta数值积分结果比较,说明此方法的精确、高效性。     

切比雪夫谱有限条

文[1]用有限条法分析薄板弯曲时,将薄板分成若干条带,条的端部以振动梁函数等特殊函数为基函数,对不同的边界条件,取不同基函数,给计算和程序编制带来麻烦。本文将以切比雪夫多项式作为有限条端部的基函数,並通过建立边界约束方程来统一处理各种边界条件,边界约束方程可简化有限条的刚度矩阵,计算结果与理论解较吻合。     

一种新的求解非线性系统周期解方法

本文利用切比雪夫多项式的若干良好性质,对非自治强非线性动力系统进行分析。将状态矢量在主周期上先展开谐波级数的形式,再将各谐波展开为切比雪夫多项式的形式,从而将求周期解的问题转变为非线性代数方程组的求解问题,得出一种可以方便、迅速地获得近似周期解的解析方法。此方法不依赖于小参数假设,可以用于分析强非线性问题和高维问题,而且对参数激励系统同样有效。以Duffing系统周期解的计算为例,通过与标准谐波平衡方法和四阶Runge-Kutta数值积分结果比较,说明此方法的有效性。     

断裂过程区应力应变场的数学分析

本文将条状分离区看作是位错连续分布区域,利用复变函数理论及切比雪夫多项式级数展开提出了一种统一的数学分析的方法,用以对均匀材料裂纹前方分离过程区的非线性效应进行有效的,精确的解析分析。避免以前各种模型一般都需要借助于有限元法或其他离散方法才能得到数值解答的难点。本文利用这一新方法,对断裂过程区是内聚力模型,条状屈服区模型的裂纹问题,受约束金属薄膜的裂纹问题,进行了分析计算,并与已有的结果进行了比较,得到了满意的结果。     

配点型区间有限元法

在分析Taylor展开``点逼近'区间有限元法不足的基础上, 提出了基于Chebyshev第一类正交多项式全局逼近目标函数的配点型区间有限元法. 该方法不需要计算目标函数对不确定性变量的灵敏度, 不要求不确定性变量的变化范围为小区间, 并适合求解目标函数为不确定变量非线性函数的情形. 目标函数正交展开式的系数采用Gauss-Chebyshev求积公式得到,故需要在不确定性变量所在区间内配置高斯积分点. 计算目标函数在高斯点的取值是该方法的主要工作量, 当不确定性变量数为m, 并选用高斯十点法进行积分时, 需要对系统进行12$m$次分析. 算例表明, 在其他区间有限元法失效的情况下, 配点型区间有限元法依然能够得到几乎精确的区间界限.      

Free vibration analysis of elastic rods using global collocation

This paper investigates the performance of a novel global collocation method for the eigenvalue analysis of freely vibrated elastic structures when either basis or shape functions are used to approximate the displacement field. Although the methodology is generally applicable, numerical results are presented only for rods in which one-dimensional basis functions in the form of a power series, as well as equivalent Lagrange, Bernstein or Chebyshev polynomials are used. The new feature of the proposed methodology is that it can deal with any type of boundary conditions; therefore, the cases of two Dirichlet as well as one Dirichlet and one Neumann condition were successfully treated. The basic finding of this work is that all these polynomials lead to results identical to those obtained by the power series expansion; therefore, the solution depends on the position of the collocation points only.     

High-order theory for arched structures and its application for the study of the electrostatically actuated MEMS devices

A high-order theory for arched rods and beams based on expansion of the two-dimensional (2D) equations of elasticity into Legendre’s polynomials series has been developed. The 2D equations of elasticity have been expanded into Legendre’s polynomials series in terms of a thickness coordinate. Thereby, all equations of elasticity including Hooke’s law have been transformed to corresponding equations for coefficients of Legendre’s polynomials expansion. Then system of differential equations in term of displacements and boundary conditions for the coefficients of Legendre’s polynomials expansion coefficients has been obtained. Cases of the first and second approximations have been considered in details. For obtained boundary-value problems, a finite element method has been used and numerical calculations have been done with COMSOL Multiphysics and MATLAB. Developed theory has been applied for study pull-in instability and stress–strain state of the electrostatically actuated micro-electro-mechanical Systems.     

弹性连接旋转柔性梁动力学分析

采用Chebyshev 谱方法对考虑根部连接弹性的平面内旋转柔性梁动力学特性进行研究. 基于Gauss-Lobatto 节点与Chebyshev 多项式方法对柔性梁变形场进行离散,通过投影矩阵法施加固定及弹性连接边界条件. 利用Chebyshev 谱方法获得了系统固有频率和模态振型数值解,通过与有限元方法及加权残余法的比较,验证了方法的有效性. 分析了弹性连接刚度、角速度比率、系统径长比及梁的长细比等参数对系统固有频率及模态振型的影响. 研究发现:由于系统弯曲模态、拉伸模态的频率随各参数的变化规律不一致,将出现频率转向与振型转换现象;随着弹性连接刚度、角速度比率及系统径长比的增大,低阶弯曲模态频率增大并超过高阶拉伸模态频率,随着梁的长细比的增大,低阶拉伸模态频率增大并超过高阶弯曲模态频率.      

Finite element analysis of incompressible laminar boundary layer flows

A numerical procedure was developed to solve the two-dimensional and axisymmetric incompressible laminar boundary layer equations using the semi-discrete Galerkin finite element method. Linear Lagrangian, quadratic Lagrangian, and cubic Hermite interpolating polynomials were used for the finite element discretization; the first-order, the second-order backward difference approximation, and the Crank-Nicolson method were used for the system of non-linear ordinary differential equations; the Picard iteration and the Newton-Raphson technique were used to solve the resulting non-linear algebraic system of equations. Conservation of mass is treated as a constraint condition in the procedure; hence, it is integrated numerically along the solution line while marching along the time-like co-ordinate. Among the numerical schemes tested, the Picard iteration technique used with the quadratic Lagrangian polynomials and the second-order backward difference approximation case turned out to be the most efficient to achieve the same accuracy. The advantages of the method developed lie in its coarse grid accuracy, global computational efficiency, and wide applicability to most situations that may arise in incompressible laminar boundary layer flows.     

Three-dimensional vibrations of annular thick plates with linearly varying thickness

The free vibration of annular thick plates with linearly varying thickness along the radial direction is studied, based on the linear, small strain, three-dimensional (3-D) elasticity theory. Various boundary conditions, symmetrically and asymmetrically linear variations of upper and lower surfaces are considered in the analysis. The well-known Ritz method is used to derive the eigen-value equation. The trigonometric functions in the circumferential direction, the Chebyshev polynomials in the thickness direction, and the Chebyshev polynomials multiplied by the boundary functions in the radial direction are chosen as the trial functions. The present analysis includes full vibration modes, e.g., flexural, thickness-shear, extensive, and torsional. The first eight frequency parameters accurate to at least four significant figures for five vibration categories are obtained. Comparisons of present results for plates having symmetrically linearly varying thickness are made with others based on 2-D classical thin plate theory, 2-D moderate thickness plate theory, and 3-D elasticity theory. The first 35 natural frequencies for plates with asymmetrically linearly varying thickness are compared to the finite element solutions; excellent agreement has been achieved. The asymmetry effect of upper and lower surface variations on the frequency parameters of annular plates is discussed in detail. The first four modes of axisymmetric vibration for completely free circular plates with symmetrically and asymmetrically linearly varying thickness are plotted. The present results for 3-D vibration of annular plates with linearly varying thickness can be taken as benchmark data for validating results from various plate theories and numerical methods.     

Wave propagation analysis in anisotropic and inhomogeneous uncracked and cracked structures using pseudospectral finite element method

The work deals with the development of an effective numerical tool in the form of pseudospectral method for wave propagation analysis in anisotropic and inhomogeneous structures. Chebyshev polynomials are used as basis functions and Chebyshev–Gauss–Lobatto points are used as grid points. The formulation is implemented in the same way as conventional finite element method. The element is tested successfully on a variety of problems involving isotropic, orthotropic and functionally graded material (FGM) structures. The formulation is validated by performing static, free vibration and wave propagation analysis. The accuracy of the element in predicting stresses is compared with conventional finite elements. Free vibration analysis is carried out on composite and FGM beams and the computational resources saved in each case are presented. Wave propagation analysis is carried out using the element on anisotropic and inhomogeneous beams and layer structures. Wave propagation in thin double bounded media over long propagating distances is studied. Finally, a study on scattering of waves due to embedded horizontal and vertical cracks is carried out, where the effectiveness of modulated pulse in detecting small cracks in composites and FGMs has been demonstrated.     

A Fourier–Chebyshev spectral collocation method for simulating flow past spheres and spheroids

An accurate Fourier–Chebyshev spectral collocation method has been developed for simulating flow past prolate spheroids. The incompressible Navier–Stokes equations are transformed to the prolate spheroidal co‐ordinate system and discretized on an orthogonal body fitted mesh. The infinite flow domain is truncated to a finite extent and a Chebyshev discretization is used in the wall‐normal direction. The azimuthal direction is periodic and a conventional Fourier expansion is used in this direction. The other wall‐tangential direction requires special treatment and a restricted Fourier expansion that satisfies the parity conditions across the poles is used. Issues including spatial and temporal discretization, efficient inversion of the pressure Poisson equation, outflow boundary condition and stability restriction at the pole are discussed. The solver has been validated primarily by simulating steady and unsteady flow past a sphere at various Reynolds numbers and comparing key quantities with corresponding data from experiments and other numerical simulations. Copyright © 1999 John Wiley & Sons, Ltd.     

A trigonometric series expansion method for the Orr-Sommerfeld equation

A trigonometric series expansion method and two similar modified methods for the Orr-Sommerfeld equation are presented. These methods use the trigonometric series expansion with an auxiliary function added to the highest order derivative of the unknown function and generate the lower order derivatives through successive integrations. The proposed methods are easy to implement because of the simplicity of the chosen basis functions. By solving the plane Poiseuille flow(PPF), plane Couette flow(PCF), and Blasius boundary layer flow with several homogeneous boundary conditions,it is shown that these methods yield results with the same accuracy as that given by the conventional Chebyshev collocation method but with better robustness, and that obtained by the finite difference method but with fewer modal number.     

Stability of unstably stratified shear flow between parallel plates *

The linear stability of unstably stratified shear flows between two horizontal parallel plates has been investigated. The eigenvalue problem was solved numerically by making use of the expansion method in Chebyshev polynomials, and critical Rayleigh numbers were obtained accurately in the Reynolds number range of [0.01,100]. It was found that the critical Rayleigh number for two-dimensional disturbances increases with an increase of the Reynolds number. The result strongly supports previous stability analyses except for the analysis by Makino and Ishikawa (1985) in which a decrease of the critical Rayleigh number was obtained. For some cases, a discontinuity in the critical wavenumber occurs, due to the development of two extrema in the neutral stability boundary.