动力系统双变量重构的行为

利用重构变换的牙可比行列式，研究了常微方程系统双变量重构的行为，发现重构的拓朴性质依赖于重构维数和延迟时间，及依赖于双变量本身的稳定性质，由于这种稳定性的影响，使得双变量时间序列的重构不一定优于单变量的重构。利用条件Ｌｙａｐｕｎｏｖ指数分析了这种影响。数值结果还表明，双变量重构时嵌入的充分条件与Ｔａｋｅｎｓ指数分析了这种影响。数值结果还表明，双变量重构时嵌入的ａｋｅｎｓ定理一致。文中以Ｒｏｓｓｌｅ     

一维格子结构动力分析的传递矩阵法

本文将传递矩阵法推广应用于分析一维格子结构的波传播和动力响应特性。一个格子结构的元件可分为主元件和次元件,传递矩阵沿主元件形成并考虑次元件的作用。文中通过例子说明形成一个周期单元传递矩阵的方法,指出利用传递矩阵计算无限或半无限长格子结构波传播的传播常数及有限长格子结构固有频率和频率响应函数的方法。作为数例,文中计算了一维平面格子结构的传播常数和频率响应函数。     

高维周期解支的数值延续算法

本文提出了描述动力系统周期解的一个定解问题，给出其的充要条件，并基于样条函数在Ｇａｕｓｓ点的配置法和伪弧长算法，发展了适合高维、复杂构形、具有极值分侧周期解支的数延续算法。计算了若干非线性常微和偏微分方程组的周期解支，以说明本文数值方法在解题规模，精度上的有效性。     

欧拉动弯曲问题的混沌控制

用周期激振力法、参数共振微扰法和延时反馈控制方法对欧拉动弯曲问题的混沌行为进行了控制，通过计算受控系统的Ｌｙａｐｕｎｏｖ指数来确定控制参数，得到受控后系统的高周期、低周期和准周期等稳定周期振动结果，。     

实现有限元并行分布计算的一种新策略

在几种典型的计算网络上，给出了实现有限元并行分布计算的一种全新策略。它对子结构的划分方式没有任何限制，使结构划分方式对通讯不产生任何影响，并利用所谓的Σ通讯完成有关迭代计算。这种策略广泛适用于多项式加速法的并行分布迭代计算，使有限元并行分布计算的算法及程序与具体的计算网络有很好的分离性，同时也很大程度地保留了已有串行有限元算法及程序的优点。以预处理的共轭斜量法为例，在ＩｎｍｏｓＴ８００Ｔｒａｎｓｐｕｔｅｒ系统上实现了有限元并行分布计算。通过数值算例，验证了本文方法的可行性与有效性。     

基于奇异谱分析的BDS卫星钟差周期项提取

在建立星座自主时间基准时，必须扣除卫星钟差的周期波动，以避免将其引入系统时间。为准确地扣除周期波动，提出一种基于奇异谱分析（SSA）的BDS卫星钟差周期项提取方法。首先，经SSA分解获得钟差信号的多个重构成分；然后，引入重标极差分析方法计算各重构成分的Hurst指数，根据Hurst指数辨识钟差的低频信号主导分量和高频信号随机分量，从而实现周期项的重建提取。对不同类型的BDS卫星钟进行了分析和研究，结果表明：所提方法能够准确提取钟差周期项，所提取周期项的频谱特征比钟差多项式拟合残差的频谱更为清晰；利用多项式拟合方法扣除周期项后，万秒稳和日稳平均分别提升4.1%和2.4%，而利用所提方法扣除周期项后，万秒稳和日稳平均分别提升20.9%和34.1%，提高了卫星钟频率稳定度，为星座自主时间基准的建立提供了基础。     

用格子Boltzmann模型模拟Lotka-Volterra系统

用格子Boltzmann方法模拟了Lotka-Voltcrra系统的动力学行为．通过使用多尺度技术,我们得到了扩散反应系统。以及零扩散的反应系统。从而得到常微分方程系统，这个理论结果表明该格子Boltzmann模型包含了Lotka-Volterra系统的非线性行为，数值结果表明附加的人为扰动可以改变某些初始分布。     

直角和贴体坐标系下两种格子Boltzmann方法的比较研究

采用格子Boltzmann方法(LBM)和改进的插值格子Boltzmann方法(GILBM)研究了45°斜方腔的顶盖驱动流和Roach通道内的流动特性,并与基准解进行了对比。结果表明,对于45°斜方腔的顶盖驱动流,当雷诺数较小时,两种方法的计算结果与基准解吻合较好;但当雷诺数较大时,采用LBM的计算结果准确性降低,而基于GILBM方法得到的结果准确度升高,且计算稳定性好。对于Roach通道内的流体流动而言,两种方法的计算精度和复杂边界的复杂程度与雷诺数大小有关。根据流场边界形状的复杂程度,网格划分与计算精确度的不同要求,两种方法各有利弊。     

用概率型模型测量格子气的粘性系数

在Ｂｏｌｔｚｍａｎｎ假设下，提出一种概率型格子气模型，并通过对平面Ｐｏｓｅｕｉｌｌｅ流的模拟，用数值方式“测量”了格子气体的粘性系数。数值结果同理论上得到的解析粘性系数进行了定量比较，得到了非常符合的结果。     

实现有限元并行分布计算的一种新策略

在几种典型的计算网络上，给出了实现有限元并行分布计算的一种全新策略。它对子结构的划分方式没有任何限制，使结构划分方式对通讯不产生任何影响，并利用所谓的∑通讯完成有关迭代计算。这种策略广泛适用于多项式加速法的并行分布迭代计算，使有限元并行分布计算的算法及程序与具体的计算网络有很好的分离性，同时也很大程度地保留了已有串行有限元算法及程序的优点。以预处理的共轭斜量法为例，在Ｉｎｍｏｓ Ｔ８００ Ｔｒａｎ     

基于核重构思想的配点型无网格方法的研究--一维问题

无网格方法按其离散原理可分为Galerkin型、配点型等。其中Galerkin型无网格方法的实施需要背景网格，不属于真正的无网格法；配点型无网格方法的实施不需要背景网格，是真正的无网格法。本文首先介绍了重构核点法的基本原理，然后基于核重构思想，与配点法相结合，以一维问题为例，研究了配点型无网格方法，对该方法构造过程中的近似函数及其导数的计算、修正函数的计算及方法的实现等问题进行了探讨。并结合若干典型算例，检验了其计算精度与收敛姓。     

格子Boltzmann方法模拟多孔介质惯性流的边界条件改进

格子Boltzmann方法可以有效地模拟水动力学问题,边界处理方法的选择对于可靠的模拟计算至关重要.本文基于多松弛时间格子Boltzmann模型开展了不同边界条件下,周期对称性结构和不规则结构中流体流动模拟,阐述了不同边界条件的精度和适用范围. 此外,引入一种混合式边界处理方法来模拟多孔介质惯性流, 结果表明:对于周期性对称结构流动模拟,体力格式边界条件和压力边界处理方法是等效的,两者都能精确地捕捉流体流动特点; 而对于非周期性不规则结构,两种边界处理方法并不等价,体力格式边界条件只适用于周期性结构;由于广义化周期性边界条件忽略了垂直主流方向上流体与固体格点的碰撞作用,同样不适合处理不规则模型;体力-压力混合式边界格式能够用来模拟周期性或非周期性结构流体流动,在模拟多孔介质流体惯性流时,比压力边界条件有更大的应用优势,可以获得更大的雷诺数且能保证计算的准确性.      

格子Boltzmann方法模拟多孔介质惯性流的边界条件改进

格子Boltzmann方法可以有效地模拟水动力学问题,边界处理方法的选择对于可靠的模拟计算至关重要.本文基于多松弛时间格子Boltzmann模型开展了不同边界条件下,周期对称性结构和不规则结构中流体流动模拟,阐述了不同边界条件的精度和适用范围.此外,引入一种混合式边界处理方法来模拟多孔介质惯性流,结果表明:对于周期性对称结构流动模拟,体力格式边界条件和压力边界处理方法是等效的,两者都能精确地捕捉流体流动特点;而对于非周期性不规则结构,两种边界处理方法并不等价,体力格式边界条件只适用于周期性结构;由于广义化周期性边界条件忽略了垂直主流方向上流体与固体格点的碰撞作用,同样不适合处理不规则模型;体力–压力混合式边界格式能够用来模拟周期性或非周期性结构流体流动,在模拟多孔介质流体惯性流时,比压力边界条件有更大的应用优势,可以获得更大的雷诺数且能保证计算的准确性.     

爆炸下球壳变形空间周期分布的理论计算方法

早期研究提出了对振动叠加应变增长现象的解剖式分析方法，进而发现爆炸加载下带扰动源球壳上的弯曲波和壳体变形呈空间周期分布的规律。参考Timoshenko梁的弯曲理论，基于平截面假定和壳体发生较小的弯曲变形的假设，推导出球壳上弯曲波波速和波长的关系，计算得到最短弯曲波和与膜振动频率相近的弯曲波的波速，还结合早期研究提出的壳体变形分布周期与弯曲波波速的关系，计算得到了壳体变形空间分布的周期。结果表明:（1）理论计算结果与数值仿真结果基本吻合，其中弯曲波波速的计算结果与数值仿真结果相差在15%以内，壳体变形空间分布周期的计算结果与数值仿真结果相差在12%以内；（2）弯曲波波长越短，波速越快，当波长无限短时，波速趋于极限值，约为声速的0.574倍。本计算方法为解剖式分析方法提供了一定的理论依据。     

平面Poisenille流动的二维Hopf分叉

本文用直接数值模拟的方法计算了二维Ｐｏｉｓｅｕｉｌｌｅ流动中扰动波的演化问题。得到了二维平衡态，在一定的波数下，Ｒｅ３９５０时，这种平衡态，将变得不稳定，模拟发现出现第二周期解，即二次分叉。     

平面Poisenille流动的二维Hopf分叉

本文用直接数值模拟的方法计算了二维Ｐｏｉｓｅｕｉｌｌｅ流动中扰动波的演化问题。得到了二维平衡态，在一定的波数下，Ｒｅ３９５０时，这种平衡态，将变得不稳定，模拟发现出现第二周期解，即二次分叉。     

含孔正交各向异性复合材料板的动应力集中

运用一种改进的非结构化四边形格子法,对含孔正交各向异性板条受面内冲击拉伸时弹性应力波的传播过程和孔边的动应力集中进行了研究．非结构化格子法采用与有限元类似的网格剖分方法,并基于围绕每个节点的积分平衡方程,并自然满足复杂边界的自由边界条件．计算中不需存储刚度矩阵,因而计算速度快、效率高、节省内存,在解决应力波传播问题中具有显著的优越性．通过对多种工况进行数值模拟,分析了材料的各向异性性质、纤维方向、孔径比、加载脉冲周期等参数对孔边动应力的影响,得到了一些规律性的结果．并与现有实验结果进行对比,验证了该方法的有效性．     

三维非均匀介质中弹性波传播的数值模拟

提出了一种三维非均匀介质中弹性波传播数值模拟的方法,文中称为三维格子法。该算法是二维格子法（一种二维非均匀介质中P－SV波传播的数值模拟算法）向三维非均匀介质情况的推广。在空间离散上该文方法与有限元方法类似,容许根据连续体的形状和介质分界面任意剖面网格,且自然满足自由表面边界条件。不同于常规有限差分法在各个节点上满足动力学微分方程,该算法通过满足各节点周围格子的整体平衡（积分平衡方程）来对问题进行求解,三维格子法所需的计算机内存及计算耗时与同阶精度的规则网格有限差分法相当。算例表明,该文提出的三维格子法具有较高的精度且可很好地模拟三维复杂形状地表对弹性波的反射和绕射。     

sine—Gordon方程的动力学行为讨论

本文以ｓｉｎｅ－Ｃｏｒｄｏｎ方程为全了把无穷维动力系统理论和非一动力学方法结合起来讨论无穷维动力学系统的动力学行为的一种方法，首先证明了有限维广义渐惯性流形的存在性。然后给出了该流形上常微分方程组，且讨论了解的不变子空间结构。最后用摄动和数值计算讨论了动力学行为，所得结果与（６）所得的一致。这种方法也许为时军蝗分析提供了一条有效途径     

挤压油膜阻尼器—滑动轴承—刚性转子系统的稳定性及分岔行为

对挤压油膜阻尼器－滑动轴承－转子系统的稳定性及分岔行为进行了研究，由于该动力系统为一强非线性系统，具有复杂的非线性现象。本文采用Ｆｌｏｑｕｅｔ理论对其周期解的稳定性进行了计算分析：随着系统参数的变化，该系统将出现稳态周期解、准周期分岔、倍周期分岔。文中也对系统平衡点的稳定性进行了分析，讨论了Ｈｏｐｆ分岔行为。