极坐标系下Fourier-Legendre谱元方法与有限差分法数值扩散的比较

提出一种Fourier-Legendre谱元方法用于求解极坐标系下的Navier-Stokes方程,其中极点所在单元的径向采用Gauss-Radau积分点,避免了r=0处的1/r坐标奇异性。时间离散采用时间分裂法,引入数值同位素模型跟踪同位素的输运过程验证数值模拟的精度,分别利用谱元法和有限差分法的迎风差分格式求解匀速和加速坩埚旋转流动中的同位素方程。计算结果表明,有限差分法中的一阶迎风差分格式存在严重的数值假扩散,二阶迎风差分格式的数值结果较精确,增加节点可以有效地缓解数值扩散。然而,谱元法具有以较少节点得到高精度解的优势。     

比例边界坐标插值方法在谱元法中的应用——无穷域Euler方程的数值模拟

将比例边界坐标插值方法引入谱元法, 构成比例边界谱单元, 对无穷域Euler方程进行数值模拟.阐述了比例边界谱单元的基本使用方法以及基于比例边界谱元的Runge-Kutta间断Galerkin方法求解Euler方程的过程;计算了无穷域圆柱和NACA0012翼型绕流问题, 并与已有结果进行了比较, 显示了计算结果的正确性.用基于比例边界谱元的间断Galerkin方法求解无穷域Euler方程时, 最多只需将求解域划分为2个子域, 避免了一般谱方法将求解域划分为9个或者27个子域的麻烦. 比例边界谱单元为无穷域Euler方程的直接求解提供了一个可供参考的方法.     

非正交曲线坐标的叶轮机械三元跨声速流方程和特征理论

本文使用张量运算,导出非正交曲线坐标下三元弱守恒型方程、时间相关流面迭代方程及特征理论.它们对于求解叶轮机械三元带激波跨声速流场,进行周向、径向畸变的数值分析及颤振不稳定气动力计算等都是有价值的. 1.非正交曲线坐标下的三元弱守恒型方程 叶轮机械非正交曲线坐标三元流方程     

高阶边界元奇异积分的一种通用高效计算方法

针对三维边界元法中曲面单元上的(弱、强、超)奇异积分提出了一种通用高效的计算方法。经极坐标变换,将奇异积分转化为常规积分;采用数值方法计算Cauchy主值积分和Hadamard有限项积分系数;引入保角变换和反曲变换消除因单元畸形或因积分点靠近单元边界而引起的周向积分奇异性。该方法可以统一处理(弱、强、超)奇异积分,并且只需要知道核函数的奇异阶数和少数几个点上的被积函数值,不依赖于积分和函数的具体选取;所需的积分点少,精度高,并且受单元畸形程度影响较小,稳定性好。采用该方法计算了声学和弹性力学中的典型奇异积分,并结合二阶Nystrm方法求解了弹性力学的边界积分方程,验证了方法的高精度和高效性。本文数值积分程序可向作者索取。     

三维问题边界元法中几乎奇异积分的正则化算法

当源点靠近边界单元时,边界积分方程通常存在几乎奇异积分的计算难题.基于三角形单元,将源点到单元的距离与单元特征长度比值定义为接近度,用于度量边界单元中积分奇异性的程度.将单元上的面积分在局部的极坐标系ρθ下表示,利用一些初等函数的积分公式,获得对变量ρ作单层积分的解析表达式.几乎强奇异和超奇异面积分被转化为沿单元围道上一系列线积分,而Gauss数值积分能够有效计算这些线积分.应用该算法分析三维弹性薄壁结构获得了成功.     

轴向槽动压滑动轴承非线性油膜力解析模型

本文中提出了一种求解流体润滑轴向槽径向滑动轴承非线性油膜力的解析模型.采用油膜气穴边界条件,基于Sturm-Liouville理论,求解了非线性油膜的压力分布.为了便于求解油膜动压润滑的Reynolds方程,将油膜压力函数分解为特解和通解相加的形式,润滑油膜的破裂位置通过连续性条件确定.运用分离变量法,将特解的压力分布分解为周向分离函数和轴向分离函数相加的形式,周向分离函数运用Sommerfeld变换求解.将通解的压力分布分解为周向分离函数和轴向分离函数相乘的形式.采用变量代换,将周向分离函数方程转化为Sturm-Liouville型方程,根据边界条件求得本征值和本征函数系,进而得到通解的周向压力分布;通过求解微分方程,得出轴向分离函数为含本征值的双曲正切函数.在油膜完备区域,对油膜压力分布的解析表达式进行积分,从而求得有限宽轴向槽径向滑动轴承非线性油膜力.计算结果表明:本文中提出的方法和有限差分法的结果吻合得较好,验证了本文中所提出解析模型的正确有效性.     

非均质问题中的无网格边界单元法

介绍了一种不需要内部网格计算非均匀介质问题的边界元算法.该算法是建立在一种能将任何区域积分转换成边界积分的径向积分转换法基础上,首先用对应各向同性问题的基本解来建立以正规化位移表示的非均质问题的积分方程,然后用径向积分转换法将出现在积分方程中的区域积分转换成边界积分,从而形成不需要使用内部网格来计算区域积分的纯边界元算法.与其它无网格法相比,此方法需要很少的内部点,有些问题甚至不需要内部点都能得到满意的结果,因此,可以计算大型的三维非均匀介质工程问题.由于此方法继承了边界元和无网格算法的优点,因而具有广阔的发展前景.     

三维变系数热传导问题边界元分析中几乎奇异积分计算

在边界积分的数值计算过程中,当源点离积分单元很近时,边界积分就会具有几乎奇异性,此时不能直接用高斯数值积分公式计算几乎奇异积分。本文以三维非均质热传导问题为例,介绍了一种计算几乎奇异边界积分的新方法。首先,采用Newton-Raphson迭代算法确定积分单元上离源点最近的点;然后,将积分单元上任意一点的坐标在最近点处展开成泰勒级数,并计算源点到积分单元任意点的距离;最后,将距离函数代入几乎奇异边界积分中,并运用指数变换方法导出积分单元上几乎奇异积分的计算公式。文中给出了两个非均质热传导问题的算例来验证所述方法的正确性、有效性和稳定性。     

折线型裂纹对SH波的动力响应

利用Fourier积分变换方法,得出了无限平面中用裂纹位错密度函数表示的单裂纹散射场．根据无穷积分的性质,把单裂纹的散射场分解为奇异部分和有界部分．利用单裂纹的散射场建立了折线裂纹在SH波作用下的Cauchy型奇异积分方程．根据折线裂纹散射场和所得的积分方程讨论了裂纹在折点处的奇性应力及折点处的奇性应力指数．利用所得的奇性应力定义了折点处的应力强度因子．对所得Cauchy型奇积分方程的数值求解,可得裂纹端点和折点处的动应力强度因子。     

初应力对单向复合材料圆柱板周向波的影响

基于“增量变形力学”理论,研究了径向和轴向均匀初应力作用下单向复合材料圆柱板中周向波的传播特性,应用Legendre多项式方法求解了耦合波动方程。讨论了单向复合材料纤维方向分别为周向和轴向时,初应力对圆柱板中的周向类Lamb波和SH波的影响。数值分析结果表明初应力对周向类Lamb波和周向SH波的影响是非常不同的;轴向初应力对频散曲线、位移和应力分布的影响与径向初应力的影响也不相同。     

Acoustic scattering by a sphere in the time domain

A sound pulse is scattered by a sphere leading to an initial–boundary value problem for the wave equation. A method for solving this problem is developed using integral representations involving Legendre polynomials in a similarity variable and Volterra integral equations. The method is compared and contrasted with the classical method, which uses Laplace transforms in time combined with separation of variables in spherical polar coordinates.     

A numerical investigation of a spectral‐type nodal collocation discontinuous Galerkin approximation of the Euler and Navier–Stokes equations

In this paper, we focus on the applicability of spectral‐type collocation discontinuous Galerkin methods to the steady state numerical solution of the inviscid and viscous Navier–Stokes equations on meshes consisting of curved quadrilateral elements. The solution is approximated with piecewise Lagrange polynomials based on both Legendre–Gauss and Legendre–Gauss–Lobatto interpolation nodes. For the sake of computational efficiency, the interpolation nodes can be used also as quadrature points. In this case, however, the effect of the nonlinearities in the equations and/or curved elements leads to aliasing and/or commutation errors that may result in inaccurate or unstable computations. By a thorough numerical testing on a set of well known test cases available in the literature, it is here shown that the two sets of nodes behave very differently, with a clear advantage of the Legendre–Gauss nodes, which always displayed an accurate and robust behaviour in all the test cases considered.Copyright © 2012 John Wiley & Sons, Ltd.     

Multidomain pseudospectral methods for nonlinear convection-diffusion equations

Multidomain pseudospectral approximations to nonlinear convection-diffusion equations are considered. The schemes are formulated with the Legendre-Galerkin method, but the nonlinear term is collocated at the Legendre/Chebyshev-Gauss-Lobatto points inside each subinterval. Appropriate base functions are introduced so that the matrix of the system is sparse, and the method can be implemented efficiently and in parallel. The stability and the optimal rate of convergence of the methods are proved. Numerical results are given for both the single domain and the multidomain methods to make a comparison.     

三维边界元法高阶元几乎奇异积分半解析法

分析了三维边界元法高阶曲面单元几何特征,定义接近度来表征源点与积分单元的接近程度.利用源点在积分单元上的垂足点建立局部极坐标系,构造与几乎奇异积分核函数具有相同奇异性的近似函数.从奇异积分核函数中扣除其近似函数,分离出积分核中主导的奇异函数部分,将奇异积分分解为规则核函数和奇异核函数两项积分.规则核函数积分应用常规Gauss数值积分计算,奇异核函数积分在局部极坐标系ρθ下分离积分变量ρ和θ,对ρ积分建立解析计算列式,对θ积分应用常规Gauss数值积分计算,从而对三维位势问题高阶边界单元几乎强奇异和几乎超奇异积分建立一种新的半解析算法.给出了若干温度场算例,采用边界元法高阶单元几乎奇异积分半解析法计算了近边界内点位势和位势梯度,并与线性单元正则化算法计算结果对比,结果证明提出的半解析法计算几乎奇异面积分和薄壁结构更加高效.      

Reducing three-dimensional elasticity problems to two-dimensional problems by approximating stresses and displacements by legendre polynomials

Shell equations are constructed in orthogonal curvilinear coordinates using approximations of stresses and displacements by Legendre polynomials. The order of the constructed system of differential equations is independent of whether stresses and displacements or their combination are specified on the shell surfaces, which provides the correct formulation of the surface conditions in terms of both displacements and stresses. This allows the system of differential equations of laminated shells to be constructed using matching conditions for displacements and stresses on the contact surfaces. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 48, No. 3, pp. 179–190, May–June, 2007.     

Integration of B-spline geometry and ANCF finite element analysis

The goal of this investigation is to introduce a new computer procedure for the integration of B-spline geometry and the absolute nodal coordinate formulation (ANCF) finite element analysis. The procedure is based on developing a linear transformation that can be used to transform systematically the B-spline representation to an ANCF finite element mesh preserving the same geometry and the same degree of continuity. Such a linear transformation that relates the B-spline control points and the finite element position and gradient coordinates will facilitate the integration of computer aided design and analysis (ICADA). While ANCF finite elements automatically ensure the continuity of the position and gradient vectors at the nodal points, the B-spline representation allows for imposing a higher degree of continuity by decreasing the knot multiplicity. As shown in this investigation, a higher degree of continuity can be systematically achieved using ANCF finite elements by imposing linear algebraic constraint equations that can be used to eliminate nodal variables. The analysis presented in this study shows that continuity of the curvature vector and its derivative which corresponds in the cubic B-spline representation to zero knot multiplicity can be systematically achieved using ANCF finite elements. In this special case, as the knot multiplicity reduces to zero, the recurrence B-spline formula causes two segments to automatically blend together forming one cubic segment defined on a larger domain. Similarly in this special case, the algebraic constraint equations required for the C ³ continuity convert two ANCF cubic finite elements to one finite element, demonstrating the strong relationship between the B-spline representation and the ANCF finite element representation. For the same order of interpolation, higher degree of continuity at the finite element interface can lead to a coarser mesh and to a lower dimensional model. Using the B-spline/ANCF finite element transformation developed in this paper, the equations of motion of a finite element mesh that represents exactly the B-spline geometry can be developed. Because of the linearity of the transformation developed in this investigation, all the ANCF finite element desirable features are preserved; including the constant mass matrix that can be used to develop an optimum sparse matrix structure of the nonlinear multibody system dynamic equations.     

Spline sectorial elements

In this paper, the quadratic and cubic splines local interpolation on a sectorial element in polar coordinates is discussed and a class of spline sectorial elements for analyses of plane and thin, problems are presented. A reasonable treatment of the assumed displacement fields for elements with nodes at the origin, (r=0) has been made so that the elements can not only characterize the geometrical properties at the origin but also remove the singularity of strains and stresses there. Some numerical examples are given to show the efficiency of the proposed elements.     

A NEW SEMI-ANALYTIC ALGORITHM OF NEARLY SINGULAR INTEGRALS IN HIGH ORDER BOUNDARY ELEMENT ANALYSIS OF 3D POTENTIAL

分析了三维边界元法高阶曲面单元几何特征,定义接近度来表征源点与积分单元的接近程度.利用源点在积分单元上的垂足点建立局部极坐标系,构造与几乎奇异积分核函数具有相同奇异性的近似函数.从奇异积分核函数中扣除其近似函数,分离出积分核中主导的奇异函数部分,将奇异积分分解为规则核函数和奇异核函数两项积分.规则核函数积分应用常规Gauss数值积分计算,奇异核函数积分在局部极坐标系ρθ下分离积分变量ρ和θ,对ρ积分建立解析计算列式,对θ积分应用常规Gauss数值积分计算,从而对三维位势问题高阶边界单元几乎强奇异和几乎超奇异积分建立一种新的半解析算法.给出了若干温度场算例,采用边界元法高阶单元几乎奇异积分半解析法计算了近边界内点位势和位势梯度,并与线性单元正则化算法计算结果对比,结果证明提出的半解析法计算几乎奇异面积分和薄壁结构更加高效.     

Legendre积分法在随机有限元法中的应用

将Legendre积分法应用于随机结构的有限元分析,针对非线性问题,建立基于Legendre积分法的随机有限元列式。选择不同的Legendre积分点数目进行算例分析,并用Monte—Carlo法的计算进行对比研究,考察该方法的有效性。计算结果表明本文提出的Legendre积分随机有限元有很高的计算效率,在精度上,较少的积分点在一阶矩、二阶矩计算上即有较高的精度,在积分点数较多时,三阶矩、四阶矩也有较高的精度。