近似平衡多项式加速度动力显式算法

利用已知初始时刻的信息，建立一种可以取到任意阶高精度的多项式加速度单步隐式算法。在该隐式方法中，待采解方程纽系数矩阵中质量阵的系数远远大于阻尼阵和剐度阵的系数，略去非对角阻尼阵和非对角刚度阵对方程组的影响，得到一种近似平衡多项式加速度动力显式计算方法。此方法的精度主要由加速度多项式插值的项数、步长、质量阵的每件数、质量刚度比（质量阵和刚度阵的范数之比）决定。在此基础上给出了这种算法的通式，进行了精度分析，结果表明：如果时间步长h足够短，n次加速度近似平衡动力显式算法的精度可以达到O（hn＋1）。算例采用5次加速度近似平衡显式算法，计算结果的精确性证明了本算法的可行性。     

传递矩阵－多项式法的改进及其应用

本文提出一种计算转子－滑动轴承－挤压油膜阻尼器转子系统的阻尼固有频率及稳定性的方法，即将传递矩奉多项式法改进后，应用于这类具有分叉结点的系统，直接得到系统的特征多项式，然后用Ｂａｉｒｓｔｏｗ－Ｎｅｗｔｏｎ僻因子法求出全部特征根。本方法具有占内存少、计算速度快、数值稳定等特点，文中对一转子模型进行了计算，其结果与实测基本相符。     

修正荷载的Wilson-θ法的稳定性与精度分析

Wilson-θ法求得的位移、速度与加速度不满足t时刻的动力平衡方程。提出修正荷载的Wilson-θ法:增加一个荷载δF(t),使得t时刻的动力平衡方程得以满足;将-θ′δF(t)作为荷载加入t+Δt时刻的计算中,当θ′=0时,修正荷载的Wilson-θ法退化为Wilson-θ法。对应于不同的θ′值,在无条件稳定的前提下,θ的取值范围也不同。定义了逐步积分法中的计算误差。计算结果表明:计算误差与θ值成正相关,当θ′=0.6,无条件稳定的θ为最小值1.24,因而θ′=0.6,θ=1.24时,计算误差最小,建议在计算中采用。保持Wilson-θ法无条件稳定、几乎不增加计算量的条件下,修正荷载的Wilson-θ法可以提高计算精度。     

基于拉格朗日插值的无网格直接配点法和稳定配点法

由于无网格法中大多数近似函数均为有理式,不具有Kronecker delta性质,因此难以精确地施加本质边界条件.边界误差较大容易导致整个求解域求解结果精度低,甚至引起数值不稳定现象.文章在无网格直接配点法和稳定配点法中引入拉格朗日插值函数作为形函数,构建了拉格朗日插值配点法(LICM)和拉格朗日插值稳定配点法(SLICM).由于拉格朗日插值具有Kronecker delta性质,可以像有限元法一样简单而精确地施加本质边界条件,提高这两种方法的数值求解精度.稳定配点法基于子域对强形式方程进行积分,可以满足高阶积分约束,即可以保证形函数在积分形式下也满足高阶一致性条件,实现精确积分.同时,进行子域积分还可以减少离散矩阵的条件数,从而提高算法的稳定性.进一步提高拉格朗日插值稳定配点法的精度和稳定性.通过数值算例验证这两种方法的精度、收敛性和稳定性,结果表明基于拉格朗日插值的配点法的精度优于基于重构核近似的配点法,拉格朗日插值稳定配点法的精度和稳定性均优于拉格朗日插值配点法.     

参数多项式方法的扩展——参数分布形状特征式法

在深入分析参数多项式方法数学原理的基础上，提出了参数分布形状特征式法。根据流动问题的物理本质及大量试验数据，总结出参数分布的形状特征式，将其代入流动基本方程，就可以快速而准确地求解流场。它是参数多项式方法的一种扩展，具有重要的应用前景。     

A Polynomial Approach to the Piecewise Hyperbolic Method

In this paper, a local (third-order accurate) shock capturing method for hyperbolic conservation laws is presented. The method has been made with the same idea as the PHM method, but with a simpler reconstruction. A comparison with the classic high order methods is discussed.     

脱层屈曲的高阶摄动解

本文用摄动法求解了三维复合材料中有纤维搭桥的圆形脱层的屈曲问题。搭桥纤维被简化为连续分布于脱层面上的线性弹簧。用Taylor级数展开法确定了每一级摄动方程用Bessel函数表达的解析解．给出了脱层屈曲临界载荷、挠度曲线与脱层半径、搭桥因子等的封闭关系式。根据这一结果并考虑脱层只能向一个方向挠曲的约束条件，可以进一步探讨脱层屈曲的特征。     

A Class of High Order Nonlocal Operators

We study a class of nonlocal operators that may be seen as high order generalizations of the well known nonlocal diffusion operators. We present properties of the associated nonlocal functionals and nonlocal function spaces including nonlocal versions of Sobolev inequalities such as the nonlocal Poincaré and nonlocal Gagliardo–Nirenberg inequalities. Nonlocal characterizations of high order Sobolev spaces in the spirit of Bourgain–Brezis–Mironescu are provided. Applications of nonlocal calculus of variations to the well-posedness of linear nonlocal models of elastic beams and plates are also considered.     

Nonlinear-dispersive Shallow-water Models of High Order Approximation

The class of nonlinear-dispersive (NLD) models of hydrodynamics was considered, the comparison of a series of hydrodynamic models was performed, and the selection of the most satisfactory model for carrying out economical and effective finite-difference schemes adequately modeling a wide class of natural wave processes was made.     

弹塑性扭转问题具多项式基的径向点插值无网格法

对于弹塑性扭转问题描述的椭圆变分不等式,采用具多项式基的径向点插值法无网格方法与Uzawa方法耦舍,得到了带松弛因子的离散迭代算法,并给出了数值算例,分析了参数对结果的影响.通过与有限元法比较,表明该方法是求解弹塑性扭转问题的有效的方法之一.     

高阶多维半离散中心迎风格式及其应用

提出了求解多维对流-扩散方程的四阶半离散中心迎风格式。该格式以中心加权基本无振荡(CWENO)重构为基础,同时考虑到在Riemann扇内波传播的局部速度,从而更加准确地估计出了局部Riemann扇的宽度,最终既回避了网格的交错,又降低了格式的数值粘性,建立了介于迎风格式和中心格式之间的半离散中心迎风格式。本文还将该四阶半离散中心迎风格式与涡度-流函数方法相结合,有效地求解了二维不可压Euler方程组和Navier-Stokes方程组。     

一类高精度TVD差分格式及其应用

构造了一维非线性双曲型守恒律的一个新的高精度、高分辨率的守恒型TvD差分格式。其构造思想是：首先，将计算区间划分为若干个互不相交的小区间，再根据精度要求等分小区间，通过各细小区间上的单元平均状态变量，重构各细小区间交界面上的状态变量，并加以校正；其次，利用近似Riemann解计算细小区间交界面上的数值通量，并结合高阶Runge—Kutta TVD方法进行时间离散，得到了高精度的全离散方法。证明了该格式的TVD特性。该格式适合于使用分量形式计算而无须进行局部特征分解。通过计算几个典型的问题，验证了格式具有高精度、高分辨率且计算简单的优点。     

基于时序变形预测的数字图像相关加速方法

利用数字图像相关求解连续变形物体的位移场和应变场时,会遇到处理速度非常慢的问题,原因是相关算法计算量大且忽略了物体在时间轴上的变形规律。本文提出了一种基于变形预测的数字图像相关方法,该方法利用物体在时间轴上的变形规律,通过已经得到的变形值来预测后面时刻的变形初值。首先分析了物体在时间轴上的变形规律,然后结合实际应用对其进行修正,得到物体下一步变形的初值,最后通过NR方法(Newton-Raphson Method)得到物体的位移场和应变场。通过计算机模拟和金属试件拉伸实验并结合GPU(Graphic Processing Unit)编程验证了该方法的有效性。计算机模拟和金属试件拉伸实验的计算结果表明,该方法能够在计算精度保持不变的情况下,使计算速度提高4~7倍。     

转台数字锁相系统最高极限速率的确定方法

数字锁相速率系统实质是一个非线性系统。确定其最高速率极限是设计者十分关心的问 题。本文采用求锁相环非线性系统极限环频率的方法确定了系统的速率极限，并给出了计算公 式。用该方法对某转台速率系统进行了计算，结果证明了此方法的有效性。