有限元分析第九章.docVIP

应用程序 §9-1概述 有限元方法的成功和广泛应用推动了应用程序的进步和发展。随着第三代电子计算机的出现，从２０世纪６０年代未期到７０年代未期，一些大型和特大型通用程序，例如ADINA、ANSYS、ASAS、ASKA、BERSAFE、BOSOR、ELAS、MARC、MARCS、NASTRAN、NEPSAP、PASEC、SAP、STRUDL等相继问世，并在使用中不断扩充新的版本。这些程序占据了相当在的一批用户市场，利用这些程序除可以分析一般的线性问题外，还可以处理非线性连续介质、塑性、断裂力学、温升应力和蠕变、船舶结构、冲击波传播、焊接问题、粘弹性、结构优化、热传导、桥梁系统、复合材料、气动弹性、结构稳定性、海洋结构等问题[12]。这种程序中有的（ASKAT 和NASTRAN）有30~40万条语句，编制和维护这些程序要化费大量的人力。 在我国由于计算机条件的限制，大型分析程序的研制工作开展得比较迟，70年代未期以后有了较快的发展，研制了如：JIGFEX（结构分析系统）、HAJKS（航空结构分析程序系统）SPSWZ用户程序软件包系统、SDL程序系统、DDDU（多单元、多约束、多工况结构优化设计程序系统）。此外，还从国外引入了SAP－５、SAP－６、SAP－７、ADINA、ASKA等程序。从而使我国的分析水平大大前进了一步。 近十年来，微型机性能提高很快，微型机价格低、易于普及，用于微型机的结构分析软件的研制工作，在七十年代中期已在国外引起一些人的注意。在国内，JIGFEX程序系统和部分程序已经移植到微机上；1982年从美国引入的SAP－81程序也扩展成SAP－8４ SAP－5程序的一些情况。SAP－5程序1974~1976年在美国完成，1979年引入我国，是一个用于线性静力、动力分析的结构分析通用程序。这个程序在引入我国时，就国内来说是先进的，赢得了相当大的一批用户市场，国内很多计算站都装有这个程序。此外，由于这个程序的文本是公开的，便于一些研究生把它解体、修改后用于自己论文课题的研究。这里简单介绍 SAP程序的主要目的有两个：其一使读者从中了解一个通用分析程序大体要具备那些功能；其二是从中体会各方面的知识（本课程所讲内容及其它课程中的讲解内容）如何综合一个整体。此外，如果有条件，把它作为教学程序安排一次上机实习，对建立有限元方法的整体印象也是有利的。 §9-2 SAP－5程序的主要功能 静力分析 计算结构在静载荷下各节点的位移和各单元的内力或应力。载荷可以是体积力、边界力。节点集中力和温升载荷。材料常数（E、ν）可以是温度的函数。总体刚度矩阵等带宽方式存贮，且利用对称性，只存贮上三角部分的有关元素，以分块的形式记在磁盘文件上（图９－１）， 在组装总刚阵时每次在内存中组装两块，组装好后记入磁盘，再组装另外两块。解方程时先将两块读入内存进行消元，消好一块后记入另一个磁盘文件，再读入一块继续进行消元。当作用于结构的载荷分为几种工况时，将各载荷向量（即每种工况下的方程右端项）拼成一个矩阵，同时消元求解。 动力分析 动力分析是解以下形式的广义特征值问题： 求使（9-2-1）有非零解的ω及相应的非零解{φ},，ω称为固有频率、{φ}为相应的振型，它们组成一个特征对。特征对的总个数等于系统自由度个数，但工程中感兴趣的往往是最低的若干阶固有频率ω1、ω2、……、ωNF，相应的振型为{φ1}、{φ2}、……、{φNF}，有时相邻的若干阶固有频率相等（称为重根）或十分接近（称为密集根），这两种情况都可能给特征值问题的求解造成困难。 SAP－5程序要求总刚阵[K]正定，提供了两种解（9-2-1）的方法：行列式搜索法和子空间迭代法。其中的行列式搜索法求解速度较快，可以六阶以下（包括六阶）的重根问题；但要求内存空间能容纳总刚度矩阵[K]和总质量阵[M]，因而限制了算题的规模；此外，当要求的特征对个数较多（例如数十阶时），行列式搜索法会遇到数值上的困难。子空间迭代法是一种适应能力很强的方法，它的总刚阵分块记在磁盘文件上，在求解过程中只把两块调入内存空间，因而算题的规模可以比较大；缺点是求解速度比较慢，且不排除“漏掉”某一阶特征对的可能性。 一般情况下，由程序根据算题规模自行确定算法，但在有些情况下（例如算题规模不大，但要求特征对的个数却比较多），也可以人为地选择子空间迭代法。 时间历程响应 即解以下形式的方程组 其中[K]为刚度矩阵、[C]为阻尼矩阵，[M]为质量矩阵、{R(t)}为任意形式的外干扰力。外激励除了以干扰力的方式出现外，也允许以基础的运动（平动或转动）的方式出现。 SAP－5程序为（9-2-2）提供了两种解法：振型迭加法和直接积分法。其中振型迭加法适用于外激励频谱比较单纯（例如简谐干扰）的情况，绝大部分计算量将化费在求解固有频率和特征值）上