应用插值小波解第二类Fredholm积分方程

引入了一种解第二类Fredholm积分方程的新的数值算法,该数值方法利用插值小波变换将积分方程转化成线性方程组并求解,经过变换后得到的线性方程组的矩阵是一个稀疏的带状矩阵.数值算例表明,与传统算法比较该方法计算量小,并且具有较高的精度.     

三维数值流形方法及其积分区域的确定算法

该文采用长方体单元作为三维数值流形方法的数学覆盖 ,在此覆盖下物理覆盖和物体的交集是刚度矩阵的积分区域 ,采用的是Hammer积分法 .为了求解刚度矩阵 ,针对Hammer积分方法 ,该文提出了一种确定积分区域的算法和算法代码     

第一类卷积型Volterra积分方程的快速配置边值方法

针对第一类卷积型Volterra积分方程的数值解,研究其快速算法;基于特殊的多步配置方法,利用未计算的近似值,构造了高阶数值格式;通过格式,将原积分方程离散为线性方程组,其中系数矩阵可分解为Toeplitz矩阵和稀疏矩阵;利用快速Fourier变换计算该线性方程组,运算量为O(NlongN);数值例子验证了方法的高效性。     

一类具弱奇性核Volterra积分方程的配置法求解

研究一类具弱奇性核Volterra积分方程的配置法求解.  利用压缩映射定理证明了该类方程解的存在唯一性, 构造了求解这类方程的配置算法, 并对算法进行误差分析, 数值实验结果验证了理论的正确性. 该数值方法可应用于更一般的非线性Volterra积分方程.     

非线性矩阵方程数值解的广义哈密顿算法

采用具有近二阶收敛速度的算法计算一类非线性矩阵方程的数值解.根据矩阵方程的解的特征,提出一个基于正定矩阵流形几何结构的广义哈密顿算法.进而比较广义哈密顿算法与经典的多步定常迭代方法的计算行为.最后通过数值模拟表明广义哈密顿算法具有更快的收敛速度.      

三维弹性力学问题中有限元方程的预处理方法

针对三维弹性问题中有限元方程的数值求解，建立了一类简单且实用的代数多重网格预处理共轭梯度法(AMG-CG法)，详细描述了相应代数多重网格方法的粗化技术及网格转移算子的构造．由于该预处理方法能有效地降低刚度矩阵的条件数，使刚度矩阵的谱分布更集中，从而大大提高了计算效率．数值结果表明，AMGCG法对求解三维弹性问题有限元方程是十分有效和健壮的。     

非线性分数阶Volterra积分-微分方程的SCW数值方法

推导第二类Chebyshev小波(SCW)分数阶算子矩阵,利用SCW算子矩阵方法求解了一类非线性分数阶Volterra积分-微分方程.此方法将分数阶积分-微分方程转化成非线性代数方程组求解,可以简化分数阶方程的求解,所得到的数值结果表明该方法是有效和精确的.     

二维Laplace方程Neumann问题直接边界积分方程的Galerkin解法

对任意形状区域的二维Laplace方程△u（x）=0的Neumann问题,用Green公式和基本解-1/2ln｜x-y｜推导得出与之等价的直接边界识分方程,采用直接边界积分方程的Galerkin解法来解该第二类Fredholm积分方程,在进行边界离散化处理时采用常单元。为了提高数值计算的误差精度,在形成线性代数方程组的刚度矩阵元素时,对二重积分的内层积分采用精确积分表达式,外层积分使用Gauss数值积分,数值实验表明该方法的有效性和实用性。     

基于精细时程积分的结构动力响应降维分析

利用指数矩阵精细算法及状态方程直接积分法,讨论了求解动力响应问题的时程积分方式。通过选择代数精度高的Cotes积分,得出了计算精度非常高的动力响应结果。采用减缩主从自由度的精细时程积分算法对动力方程进行降维积分,通过保留指定的主自由度,删除其余的自由度来减小质量阵、阻尼阵和刚度阵的维数,既降低了指数矩阵的维数又保持了必要的计算精度,使指数矩阵分解所需时间大为降低。数值算例表明所给方法在保障求解精度的前提下具有很高的求解效率。     

Laplace方程的Galerkin边界元解法

Galerkin方法是基于变分原理基础上的一种把微分方程或积分方程转化为等价的变分方程。通过离散变分方程求原方程数值解的数值计算方法。把Laplace方程的边值问题转化为边界积分方程后，通过与边界积分方程等价的变分形式，采用线性单元，利用Galerkin边界元方法求解。在计算单元刚度矩阵时，对二重积分的第一重使用精确积分，第二重使用数值积分，从而有效克服了奇异积分的计算，数值算例验证了Galerkin方法误差的理论结果。     

p自适应分析求解快速迭代算法的研究

提出一种结构系统ｐ自适应有限元分析快速求解算法．该方法以预处理共轭梯度法为基础，采用了逐单元矩阵处理技术和基于单元的分块预处理矩阵技术；提出以高斯消元法求得ｐ＝１或２时问题的解作为高阶问题初始解的求解策略，避免了显式组装系统总体刚度矩阵，具有内存需求少、收敛速度快的优点．     

薄板稳定性分析中的一个精化不协调元

基于直接刚度法建立的十二参数任意四边形薄板精化不协调单元,采用三次和线性插值函数的某种组合构造了其单元几何刚度阵,并着重对几何刚度阵进行精化,给出了精化几何刚度阵显式。数值结果表明,其收敛速度快、精度高、易于实施。     

有限元单刚矩阵计算的神经网络法

固体力学有限元计算问题可以采用神经网络进行实时求解，神经网络可用改进的Ｈｏｐｆｉｅｌｄ网络。求解之前必须确定信号网络的联接权矩阵（相应结构刚度矩阵）。分析了单元刚度矩阵的组成，提出用ＢＰ网络来实现单刚矩阵的实时计算方法。分别用传统有限元法和ＢＰ网络法求解了一三角形单元的单刚矩阵元素，结果表明ＢＰ网络法计算具有实时性，精度满足工程要求。     

粘弹性单元的刚度矩阵

:本文分析了粘弹性材料的应力应变关系 ,导出了粘弹性单元的刚度矩阵 ,该刚度矩阵便于用直接刚度法装配结构刚度矩阵。文中还给出了改写后的动力方程 ,这种动力方程非常适合在设置粘弹性单元的结构中进行时程分析。     

梁式结构非线性有限元分析及其在微型机上的实现

本文讨论梁式结构几何非性静力分析的有限元方法,推导了空间梁大位移切线刚度矩阵,研究了在微型机上实施的方法。为了提高计算效率,采取了分别管理单元线性刚度阵及非线性刚度阵、一次形成线性总刚度阵而多次修改总刚度阵、按内存资源自动分块法修改总刚度阵的策略。节省了计算时间及内存资源,从而在微型机上实现了结构非线性分析。在IBM-PC-XT、机上编制了结构非线性静力分析有限元程序NONDJ-W。给出了对梁、平面柔性刚架、空间柔性刚架的计算实例。有解析近似解的算例表明用有限元的计算结果具有一定的精确度。     

覆盖位移函数对数值流形方法刚度矩阵的影响

基于数值流形方法的覆盖位移函数理论,从覆盖位移函数对刚度矩阵形成关系的角度出发分析了其对刚度矩阵求解的影响。通过分析可知,刚度矩阵中矩阵元素数值递增速度与覆盖位移函数中绝对坐标x、y相关。并提出一种改进的覆盖位移函数来改善在采用高阶覆盖位移函数情况下刚度矩阵的局部刚度过大的问题,以利于矩阵求解计算。最后通过算例验证了该方法是有效的。     

基于虚拟质量的飞机颤振快速分析方法

针对研制初期,飞机结构设计参数存在严重不确定性,颤振设计工作难以快速有效推进的问题,开展了方案设计阶段飞机颤振快速分析方法研究.采用虚拟质量法,通过在原结构一定位置施加虚拟质量,建立可覆盖局部结构参数变化的统一模态振型,然后将此虚拟质量从质量矩阵中移除,以防止影响原结构的动力学特性,并采用虚拟质量模态振型对不同参数下的...     

变截面箱形梁畸变分析的刚度法

利用箱形梁的畸变微分方程和弹性地基梁挠曲微分方程之间的比拟关系 ,推导出变截面箱形梁畸变分析的单元刚度矩阵 ,该单元刚度矩阵与一般杆系结构的单元刚度矩阵完全相似 ,使得箱形梁桥畸变分析可按一般杆系结构的刚度法进行 .     

关于复模态分析法的推广

把复模态分析法推广到刚度矩阵可以奇导且质量矩阵、阻尼矩阵和刚度矩阵均不必为对称的情形。利用常微分方程组理论在较一般条件下求出了线性有阻尼多自由度振动系统对任意外激励的精确响应。在较特殊的条件下，与复模态分析法的结果是一致。     

含有负刚度夹杂复合材料稳定性问题

Lyapunov间接方法可用于判断含有负刚度夹杂复合材料的稳定性.本文从动力学方程推导Lyapunov间接方法判据.在不考虑阻尼的情况下,可以从质量矩阵和刚度矩阵相对特征值的符号来判断含有负刚度夹杂复合材料的稳定性.本文通过算例论证该方法的可行性.