关于经典Appell-Hamel例——————分析力学札记之九

 不可积分微分约束系统动力学有Appell-Chetaev模型和Vacco模型,对经典Appell-Hamel例两个模型有同样解.     

单点子域积分与多点子域积分

基于单点子域精细积分的思想,针对抛物线型热传导方程初边值问题,提出了多点子域积分的概念,推出了一种多点子域积分的FTCS格式。该格式为显格式,并证明其为无条件稳定。数值算例表明,多点子域积分的FTCS格式具有比单点子域积分的FTCS格式收敛速度快的特点。     

泊松方程中域积分化为边界积分的方法

本文讨论了二维和三维泊松方程中域积分化为边界积分的方法。对于形如x~ig_x(y,z)、y~ig_x(x,z)和z~ig_z(x,y)的荷载给出了域积分转化为边界积分的正确公式。而对于复杂荷载,利用泰勒展开将域积分近似地转化为边界积分并给出了误差估计。计算结果表明利用本文方法可大大节省计算时间。因此,本文方法是一种十分有效的方法。     

层流边界层的有限差分计算

前言在解决边界层的问题时有两种方法,一种称为精确方法,对边界层基本微分方程(Prandtl方程)进行积分求解,另一种则称为近似方法,是从动量积分方程(Karman 方程)或能量积分方程来求解.在过去,用分析方法来解基本微分方程是颇困难而复杂的.Blasius 引入了速度剖面相似的假设使偏微分方程简化为常微分方...     

一类指数矩阵函数及其应用

研究了一阶常微分方程组特解的精细积分方法. 针对非齐次项为多项式、指数函数以及二者的乘积的情况,在Duhamel积分形式特解的基础上,引入了一类指数矩阵函数. 通过该类函数的线性组合即可表达出非齐次方程的特解. 建立了该类指数矩阵函数的一种高效递推算法,并在此基础上实现了特解的精细积分. 由于特解的积分过程能充分利用通解精细积分过程的中间量,因此两个精细积分过程能有机地结合起来,形成了一种高效、统一的广义精细积分法. 对上述递推算法做了进一步优化,并给出了通用的计算公式.算例结果证明了该方法的有效性.      

一种基于直接计算高阶奇异积分的断裂力学双边界积分方程分析法

相对于有限元法,边界单元法在求解断裂问题上有着独特的优势,现有的边界单元法中主要有子区域法和双边界积分方程法.采用一种改进的双边界积分方程法求解二维、三维断裂问题的应力强度因子,对非裂纹边界采用传统的位移边界积分方程,只需对裂纹面中的一面采用面力边界积分方程,并以裂纹间断位移为未知量直接用于计算应力强度因子.采用一种高阶奇异积分的直接法计算面力边界积分方程中的超强奇异积分;对于裂纹尖端单元,提供了三种不同形式的间断位移插值函数,采用两点公式计算应力强度因子.给出了多个具体的算例,与现存的精确解或参考解对比,可得到高精度的计算结果.      

薄体位势问题边界元法中的解析积分算法

薄体结构的数值分析是边界元法的难点问题之一。该文导出了一种完全解析积分算法，用这种算法计算了薄体平面位势问题边界元法中出现的几乎弱奇异、强奇异和超奇异积分。当边界离散为一系列线性单元，边界积分方程离散计算的积分可归纳为三种形式。对薄体问题，源点与积分单元距离通常相距很近，这些积分产生显著几乎奇异性，直接采用常规高斯积分不能有效计算。为此该文导出了这些几乎奇异积分的全解析计算公式。按源点与单元的距离是否为零，公式分两种情况。新算法采用全解析积分公式处理几乎奇异积分，首先精确计算出薄体问题边界未知位势和法向位势梯度，然后再进一步计算了域内点的物理参量。算例表明该文算法可处理狭长比为1．E-08的薄体问题，显示了边界元法分析薄体问题具有独特的优势。     

基于平均值定理和点积分方案的自然单元法及其程序实现

自然单元法是一种基于自然邻接点插值求解偏微分方程的无网格数值方法.它使用Voronoi图或Delaunay三角形作为背景积分网格,使用几何测度构造插值点形函数并形成刚度矩阵.平均值定理定义在未知函数定义域内任何球心(或圆心)的值等于球面(或圆周)上值的平均或加权平均,对于未知函数所满足的平衡方程是充分必要的.因此用平均值定理和点积分方案将求解域内平均应变值由散度定理转化为区域周界上的环路积分,改进传统的积分格式.算例表明,这一积分方案能进一步精简计算量和提高计算效率,是一种自适应的数值计算方法.     

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

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

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

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