两点边值问题的一种精细求解方法

将求解域均匀离散,由状态参量在相邻结点间的精细积分关系式,确定一组代数方程;并将其写成矩阵形式,代入边界条件后,代数方程组的系数矩阵可化为块三对角形式。针对这一特性,给出了一种高效的递推消元算法。由于没有离散误差,该方法具有较高的精度,不仅适用于任意边界的常规两点边值问题,还适用于奇异摄动边值问题。数值算例充分证明了本文方法的精度和效率。     

求解非线性动力系统周期解的改进打靶法

针对有周期解的动力系统边值问题可以转化为初值问题这一特点,改进了周期解的打靶 法数值求解. 在计算边界条件代数方程关于待定初值参数导数的过程中利用前一次 Runge-Kutta方法计算得到的节点函数值并通过再次利用Runge-Kutta方法获得了该导数值. 用此方法求解了Duffing方程及非线性转子---轴承系统的周期解,用Floquet理论判断了 周期解的稳定性,与普通打靶法作了比较,验证了方法的有效性.     

滑动最小二乘插值函数配点法

给出了利用滑动最小二乘法构造加权残值法中试函数的方法，对试函数中的基函数以及权函数的选取提出了建议；该试函数适用于任何定解问题，采用配点法求出试函数中的系数，进而可得到定解问题的近似解，利用该试函数对简支板的挠曲，悬臂梁的弯曲，以及中心具有小圆孔的大板的均匀拉伸等三个例子进行了数值计算，并与理论结果进行对比，同时还检验了该法的精度对结点数，配点数，以及结点影响半径的依赖情况，结果表明，该试函数适用于多种边值问题，且精度高，该法简化了选择试函数的过程，尤其适用于工程中的各种数值计算。     

含参数常微分方程边值问题的改进COLSYS算法

对常微分方程求解器ＣＯＬＳＹＳ的算法进行了分析、研究和改进,使其能够对参数方程和常微分方程的联立体系直接求解,一方面,免除了目前不得不采用的平凡常微分方程技巧所带来的大量不必要的计算,提高了效率;另一方面,由于减小了问题的规模,使条件数降低,从而改进了数值稳定性和精度。     

分段常微分方程边值问题的改进COLSYS算法

对分段定义的常微分方程边值问题的 COLSYS算法进行了改进 ,提高了其计算效率及数值稳定性 ,并且可以节省内存 ,这些优势随着分段数的增加而迅速加大。     

一种新的求解对流扩散边值问题的无网格方法

基于配点法和楔形基函数,提出了一种新的求解对流扩散边值问题的无网格方法。通过一维和二维的问题验证了该数值方法的可行性;并根据数值算例和分析,可以看到该数值方法能达到满意的收敛效果。该数值方法的隐格式形式能够有效地消除对流占优问题的数值振荡现象,是一种真正的无网格方法。     

求解Helmholtz方程基于核重构思想的最小二乘配点法

基于核重构思想构造近似函数，将配点法和最小二乘原理相结合对微分方程进行离散， 建立了Helmholtz方程的最小二乘配点格式，并分别研究了Helmholtz方程的波传播问题和 边界层问题. 通过数值算例可以发现，给出的数值计算结果非常接近于精确解，计算精度明显高于SPH 法的数值结果，且随着节点数目的增加，其精确度越来越高，具有良好的收敛性.     

光力学问题中非线性参数的一种求解方法

文中提出用非线性最小二乘法,从全场光力学条纹图中决定非线性参数。     

改进的移动最小二乘法

近年来发展的无网格方法大多采用移动员小二乘法来构造试函数，而应用移动最小二乘法形成的方程组有时会是病态的甚至奇异的，从而限制了它的发展和应用。本文采用带权正交函数作为基函数对移动最小二乘法做了改进，避免出现病态方程组，且在计算过程中不需要进行短阵求逆运算，提高了计算速度。之后，借鉴牛顿法、平衡法和摄动法对由移动最小二乘法得到的非线性代数方程组提出了新的求解方法。     

焦散线法求解裂纹应力强度因子的两种改进方法

为提高焦散线法求解应力强度因子的精度，提出了两种基于多点非线性最小二乘法的改进方法。文中推导了控制方程，编制了程序，阐述了实验技术。作为应用实例，对三种典型裂纹试件作了测试与处理。结果表明，用本文方法处理的结果比文［１］经典法的结果有很大改进。作者还对图像处理用于焦散线法应力强度因子的提取作了初步探讨［５     

改进虚拟边界算法在超声速流动问题求解中的应用

提出了一种改进的虚拟单元浸没边界法, 并与一种高阶格式的有限差分算法相结合, 运用于求解超声速复杂几何绕流问题.该算法的核心思想在于在固体边界的内部和外部分别施加满足边界关系的作用点, 使得几何边界离散更加细化, 起到了壁面附近网格局部加密的作用.采用源空间内流体点作为反距离插值算法的重构点, 有效避免了插值点数目过少而与作用点相重合情况.通过对二维激波反射现象 (马赫数为 2.81) 和三维超声速球体绕流问题 (马赫数为 1.2) 的数值模拟, 与实验结果对比表明, 本文改进算法相对一般的虚拟边界法来说能显著提高数值精度, 减小计算误差.计算结果揭示了球体绕流中剪切层、压缩波系和尾迹的相互作用导致自由剪切层失稳的机理.剪切层厚度和湍流雷诺脉动经历了线性增长、大幅度震荡和小幅度波动三个阶段, 导致剪切层表面褶皱因子变化呈指数规律增长.其湍流结构表现出明显的各向异性, 具体在流向雷诺正应力在湍流脉动中占主导地位, 激波的压缩作用对不同方向雷诺正应力的影响存在空间迟滞效应.      

A MIXED ELECTRIC BOUNDARY VALUE PROBLEM FOR AN ANTI－PLANE PIEZOELECTRIC CRACK

The analytical continuation method is adopted to solve a mixed electric boundary value problem for a piezoelectric medium under anti-plane deformation. The crack face is partly conductive and partly impermeable. The results show that the stress intensity factor is identical with the mode III stress intensity factor independent of the conducting length. But the electric field and the electric displacement are dependent on the electric boundary conditions on the crack faces and are singular not only at the crack tips but also at the junctures between the impermeable part and conducting portions.     

一种改进的位移反分析方法

为了发展可用于各种岩土工程及工程地质问题的一般的位移反分析系统方法, 本文在边界控制方法的基础上进行了改进。其目的是对于其各种复杂的系统问题确定最重要的力学参数。     

两种材料契合弹性力学问题的虚边界元-配点法

提出两种材料契合弹性力学问题的虚边界元－配点法.对于材料不同的区域分别采用各自的基本解,这样就避免了一般边界元采用Ｈｅｔｅｎｙｉ’ｓ基本解的局限性和麻烦．编制相应程序,通过实例将计算结果同理论解进行了比较,表明该方法是非常有效的．     

AN IMPROVED HYBRID BOUNDARY NODE METHOD IN TWO-DIMENSIONAL SOLIDS

The hybrid boundary node method （HBNM） is a promising method for solving boundary value problems with the hybrid displacement variational formulation and shape functions from the moving least squares（MLS） approximation. The main idea is to reduce the dimensionality of the former and keep the meshless advantage of the latter. Following its application in solving potential problems, it is further developed and numerically implemented for 2D solids in this paper. The rigid movement method is employed to solve the hyper-singular integrations. Numerical examples for some 2D solids have been given to show the characteristics. The computation results obtained by the present method are in excellent agreement with the analytical solution. The parameters that influence the performance of this method are studied through numerical examples.     

改进的奇异边界法模拟三维位势问题

奇异边界法是与基本解法相对应的一种边界型无网格数值离散方法. 该方法提出了源点强度因子的概念, 克服了传统基本解方法中最复杂最头疼的虚拟边界问题.基于边界元法中处理奇异积分的数值处理技术, 导出了源点强度因子的解析表达式, 提出了改进的无网格奇异边界法, 并进一步将该方法应用于三维位势问题. 该方法消除了传统方法中样本点的选取, 在不增加计算量的前提下, 极大地提高了奇异边界法的计算精度与稳定性.      

NUMERICAL METHOD FOR SOLVING LINEAR BOUNDARY VALUE PROBLEMS BY THE CHEBYSHEVτ_METHOD

A new τ_method is presented for the two dimensional linear boundary value problems. Theoretical and numerical analyses are presented. These results indicate that our method works nicely and efficiently.     

Timoshenko梁谐响应求解新方法探索

在空间域上采用只与结点有关的无网格方法离散,在时间域上采用精细积分方法求 解. 无网格离散过程中,利用伽辽金积分等效弱形式代替微分形式的控制方程,并 用修正变分原理满足位移边界条件,采用移动最小二乘法求解离散的形函数,把形 函数代入等效积分弱形式得到离散的二阶方程;精细积分过程中非齐次项采 用Romberg积分. 同时给出了两种不同边界条件的谐响 应求解的两个数值算例,得到了精确的数值结果.     

瞬态声振耦合问题的边界元分析及实验研究

本文探求任意外形结构在外力作用下的瞬态声场特性,研究声辐射机理,用时域边界元技术讨论外力、结构振动响应和声辐射三者之间的关系,并用实验结果加以验证,为瞬态声振耦合研究提供一种新的数值分析方法。