强Guass-seidel法求解有限元模型修正问题

有限元模型修正是一类特殊的二次反特征值问题.我们将有限元模型修正看成二次规划问题来解决,并采用非线性Gauss-Seidel方法来求解其相应的Lagrange对偶函数.最后,给出的数值实验说明方法的有效性.

L-矩阵的预条件方法及其比较定理

给出了解线性方程组Ax=b的预条件Guass-Seidel法,讨论了对于不可约的L-矩阵应用这种方法的收敛性并得到了比较定理.此外,给出了收敛最快时的系数取值.通过数值例子说明该文提出的预条件Guass-Seidel法是有效的.

α-块对角占优矩阵与两类迭代法的收敛性

探讨两类α-块对角占优矩阵,利用矩阵非奇异性和迭代矩阵的谱半径估计,证明当系数矩阵为这两类块对角占优矩阵时,线性方程组Ax=b的块Jacobi迭代法和块Guass Seidel迭代法均收敛.

基于PVM的线性方程组的一种网上并行迭代算法

针对基于PVM的桌面PC机联网而成的网络并行计算环境中，处理机的运算速度较快，而处理机间的通信相对较慢的实际情况，提出了求解线性方程组的一种分组Guass-Seidel并行迭代算法，该算法将线性方程组的增广矩阵按行分块储存在各处理机，每台处理机分别对各自的块采用Guass-Seidel迭代法进行迭代计算，其处理机间的通信较少，实现容易。并用1～24台桌面PC机联成的局域网，在PVM 3．4 on Windows2000，VC 6．0并行计算平台上编程对该算法进行了数值试验，试验结果表明，该算法较传统的Jacobi并行迭代算法和传统的Guass—Seidel并行迭代算法更优越。

虚拟手术中带摩擦的碰撞力计算及其物理仿真研究

虚拟手术训练中,碰撞力的计算必不可少,常见的碰撞力计算中,经常会忽略碰撞体之间的摩擦。提出一种带有摩擦的碰撞力计算方法。首先,建立计算机图形学物理仿真中最为常用的粒子系统,总体的常微分方程通过隐式欧拉方法求解,通过碰撞检测找出潜在的碰撞点,建立碰撞约束和摩擦约束数学模型,再通过类高斯-赛德尔算法求解出摩擦接触力。然后分别进行10次单接触面和双接触面的接触碰撞实验,分析其相应时间和稳定后起伏平均值,再进行100次多接触点碰撞实验,分析其平均运算时间。实验表明,在单接触面和双接触面接触的情况下,模型碰撞的平均响应时间为0.02 s,接触力稳定后的平均起伏范围为±0.02个单位。在多接触点接触的实验中,11个接触点接触时,平均每个接触点的运算时间为1.9 ms,满足物理仿真的实时性要求。