离散正交法(DOM)离散点的正交分析与程序设计.pdfVIP

第23卷第6期 计算力学学报 V01-23，No．6 Chinese of Mechanics December2006 2006年12月 JournalComputational 白象忠“”， 田振国1’2 (1．燕山大学建筑工程与力学学院，秦皇岛066004； 2．中国科学院力学研究所国家非线性连续介质力学重点实验室(LNM)，北京100080) 摘要：对能够求解一系列线性常微分方程组边值问题的数值计算方法——离散正交法(DOM)进行了离散点 的正交分析，给出了计算机实现数值计算的程序设计原理与计算流程图，指出了该方法能够克服传统计算方法 由于所求函数的快速增长所引起的边界效应和局部效应的缺点，给出了得到稳定计算过程的可能性。为了推广 应用，文中介绍了离散正交法的基本原理、实现方法和求解过程，讨论了采用离散正交法来求解非线性问题的处 理方法。并且以承受均布载荷的环形板为例，将采用离散正交法的计算结果与经典解作了对比。 关键词：离散正交法；边值问题；常微分方程组；数值方法 中图分类号：0441．3 文献标识码：A 理方案，并且以承受均布载荷作用的环形板为例， l 引 言 将采用离散正交法计算的结果与经典解作了对比。 自离散正交法[1]问世以来，由于计算技术逐步 应用离散正交法不仅可以求解线性问题，亦可 发展，该法至今尚未在解决板壳非线性边值问题中 以求解非线性问题。对于同时含有时间变量的动 得到广泛应用，在数值计算领域中，也未能充分发 力学问题，按照时间变量或者位置变量建立差分格 挥作用。众所周知，当求解板壳非线性问题时，通 式，仍旧可以应用离散正交法进行求解[4]。 常采用迭代过程，迭代过程中每一步的线性边值问 2离散正交法的基本原理 题都是可解的。因此，准线性化的方法是迭代过程 的基础，即计算过程中的每一步都必须是利用上一 离散正交法是对一系列线性常微分方程组的 步的结果作为下一步的基础[2’3]。这样，作为求解 边值问题进行求解的数值方法。对于一些工程问 一系列线性常微分方程组边值问题的数值计算方 题，例如，圆盘上作用轴对称载荷的弯曲问题，圆柱 壳的轴对称变形问题，都可以用常微分方程描述： 法——离散正交法(DOM)，可以成为用来解决板 jp 壳非线性边值问题的一种有效的计算手段。 iLFilt=A(z)f+q(z)(口≤X≤6) (1) U．L 为了推广应用该法的计算机实施，使计算者能 边界条件为 够正确掌握正交点的确定和计算结果的选取，顺利 D1f(n)一d1 (2) 地根据计算目的编制出计算程序来，本文对该法进 D2f(6)=d2 (3) 行了离散点的正交化分析；并在此基础上，给出了 程序设计原理和计算流程图；指出了采用该方法， qz…q。)T为竹维纵向矢量微分方程的非齐次项， 能够克服在传统计算方法中由于待求函数的快速 A(z)为扎阶方阵，