《有限元分析第九章15718.pdf》由会员分享，可在线阅读，更多相关《有限元分析第九章15718.pdf（6页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、-.z.第九章 应用程序 9-1 概述 有限元方法的成功和广泛应用推动了应用程序的进步和开展。随着第三代电子计算机的出现，从世纪年代未期到年代未期，一些大型和特大型通用程序，例如 ADINA、ANSYS、ASAS、ASKA、BERSAFE、BOSOR、ELAS、MARC、MARCS、NASTRAN、NEPSAP、PASEC、SAP、STRUDL 等相继问世，并在使用中不断扩大新的版本。这些程序占据了相当在的一批用户市场，利用这些程序除可以分析一般的线性问题外，还可以处理非线性连续介质、塑性、断裂力学、温升应力和蠕变、船舶构造、冲击波传播、焊接问题、粘弹性、构造优化、热传导、桥梁系统、复合材料、
2、气动弹性、构造稳定性、海洋构造等问题12。这种程序中有的ASKAT 和 NASTRAN有 3040 万条语句，编制和维护这些程序要化费大量的人力。在我国由于计算机条件的限制，大型分析程序的研制工作开展得比较迟，70 年代未期以后有了较快的开展，研制了如：JIGFE*构造分析系统、HAJKS航空构造分析程序系统SPSWZ 用户程序软件包系统、SDL 程序系统、DDDU多单元、多约束、多工况构造优化设计程序系统。此外，还从国外引入了 SAP、SAP、SAP、ADINA、ASKA等程序。从而使我国的分析水平大大前进了一步。近十年来，微型机性能提高很快，微型机价格低、易于普及，用于微型机的构造分析软件
3、的研制工作，在七十年代中期已在国外引起一些人的注意。在国内，JIGFE*程序系统和局部程序已经移植到微机上；1982 年从美国引入的 SAP81 程序也扩展成 SAP8版本，下面着重介绍 SAP5 程序的一些情况。SAP5 程序 19741976 年在美国完成，1979 年引入我国，是一个用于线性静力、动力分析的构造分析通用程序。这个程序在引入我国时，就国内来说是先进的，赢得了相当大的一批用户市场，国内很多计算站都装有这个程序。此外，由于这个程序的文本是公开的，便于一些研究生把它解体、修改后用于自己论文课题的研究。这里简单介绍 SAP 程序的主要目的有两个：其一使读者从中了解一个通用分析程序大
4、体要具备那些功能；其二是从中体会各方面的知识本课程所讲内容及其它课程中的讲解内容如何综合一个整体。此外，如果有条件，把它作为教学程序安排一次上机实习，对建立有限元方法的整体印象也是有利的。9-2 SAP5 程序的主要功能、静力分析 计算构造在静载荷下各节点的位移和各单元的内力或应力。载荷可以是体积力、边界力。节点集中力和温升载荷。材料常数E、可以是温度的函数。总体刚度矩阵等带宽方式存贮，且利用对称性，只存贮上三角局部的有关元素，以分块的形式记在磁盘文件上图，在组装总刚阵时每次在内存中组装两块，组装好后记入磁盘，再组装另外两块。解方程时先将两块读入内存进展消元，消好一块后记入另一个磁盘文件，再读
5、入一块继续进展消元。当作用于构造的载荷分为几种工况时，将各载荷向量即每种工况下的方程右端项拼成一个矩阵，同时消元求解。、动力分析 动力分析是解以下形式的广义特征值问题：求使9-2-1有非零解的 及相应的非零解,，称为固有频率、为相应的振型，它对 称 图 MK29-2-1-.z.们组成一个特征对。特征对的总个数等于系统自由度个数，但工程中感兴趣的往往是最低的假设干阶固有频率 1、2、NF，相应的振型为1、2、NF，有时相邻的假设干阶固有频率相等称为重根或十分接近称为密集根，这两种情况都可能给特征值问题的求解造成困难。SAP5 程序要求总刚阵K正定，提供了两种解9-2-1的方法：行列式搜索法和子空
6、间迭代法。其中的行列式搜索法求解速度较快，可以六阶以下包括六阶的重根问题；但要求内存空间能容纳总刚度矩阵K和总质量阵M，因而限制了算题的规模；此外，当要求的特征对个数较多例如数十阶时，行列式搜索法会遇到数值上的困难。子空间迭代法是一种适应能力很强的方法，它的总刚阵分块记在磁盘文件上，在求解过程中只把两块调入内存空间，因而算题的规模可以比较大；缺点是求解速度比较慢，且不排除“漏掉*一阶特征对的可能性。一般情况下，由程序根据算题规模自行确定算法，但在有些情况下 例如算题规模不大，但要求特征对的个数却比较多，也可以人为地选择子空间迭代法。、时间历程响应 即解以下形式的方程组 其中K为刚度矩阵、C为阻
7、尼矩阵，M为质量矩阵、R(t)为任意形式的外干扰力。外鼓励除了以干扰力的方式出现外，也允许以根底的运动平动或转动的方式出现。SAP5 程序为9-2-2提供了两种解法：振型迭加法和直接积分法。其中振型迭加法适用于外鼓励频谱比较单纯例如简谐干扰的情况，绝大局部计算量将化费在求解固有频率和特征值上。直接积分法则适用于外鼓励频谱比较复杂例如，短时间冲击力的情况，地震响应可作为这一类，到底选用哪种解法由用户自己决定。SAP5 程序只处理初始条件初位移和初速度为零的情况，而地震就属于这种情况，要能处理任意值的初始条件困难并不大，但初始条件的输入数据量有时会相当可观。分析过程中将求得各时间步上的节点位移和单
8、元应力或内力，但这样做数据输出量将很大。为了突出重点。程序允许用户指定需要输出的节点和单元以及输出方式。、响应谱分析 本功能为分析根底做平动情况下构造的响应节点位移，单元内力或应力。与历程响应的不同之处在于；根底的运动不是以时间历程、而是以位移谱或加速度谱的形式给定。而输出的位移和应力则按标准加以处理，作为各节点位移和各单元应力的“最大响应输出。如果输入地震谱曲线则可分析构造的地震响应；假设输入路面谱则可分析车辆的响应。分析响应的方法为振型迭加法。、频率响应 本功能为分析根底做简谐平动时构造的稳态响应。根底的频率可以在一定范围内扫描。节点位移和单元应力的“最大值将被找出，故得到的即为各节点位移
9、和单元应力的“幅频曲线。本功能中分析响应的方法也为振型迭加法。以上是 SAP5 程序的根本分析功能。下面再介绍几种常用到的辅助功能。、数据检查 用户使用 SAP5 程序算题需要输入算题的原始数据，数据的任何过失都可以导致上机失败，浪费了珍贵的机时，而输入数据的过失几乎是难以防止的。针对这种情况，SAP5程序设置了数据检查功能，它将以尽可能少的计算量不形成单刚、不组装总刚、不求解方程完成下述几项任务：检查输入数据有无格式上的错误，并将输入数据按其应有的意义加标题打印输出，供用户进一步检查。tRUMUCUK(9-2-2)-.z.检查算题的规模是否为所用机器的内存所容许。并将判断结论告知用户 求得总
10、刚阵的阶数和半带宽，供用户估算正式计算所需的机时，费用和需要的外存空间。、带宽优化 本功能是将构造的节点重新编号，以选出一种带宽较小的编号方案，这样可以大大减少数据的存贮空间和解方程的计算量。输出节点位移时仍按原来用户的编号输出。、构造绘图 当计算机配有绘图仪时，可以绘制以下输出图形：检查态只绘制构造网格图，用户可根据这种图形检查节点坐标或单元拓朴数据有无错误。进展静力分析时，绘制一份原网格图和各式况的变形图。进展动力分析时，绘制一份原网格图和所求各阶振型的振型图。上述所有图形可以是整个构造的投影图，也可以是局部构造的放大投影图。程序将自动调整比例系数并告之用户。以上是 SAP5 程序的主要功
11、以能。也是一般通用分析程序需要具备的根本功能。9-3 SAP5 程序的单元库 SAP5 程序的单元库包括了以下 11 种单元、空间桁架单元 这是一种两个节点的等截面直杆单元，它的力学模型同第节中的杆单元一样，区别仅在于单元可以在三维空间中取任意的方位，每个节点有三个自由度。单元的节点自由度为 单元刚度矩阵为 。、空间梁单元 这是一种两个节点的直梁单元，比中的梁单元复杂。单元可以在三维空间中取任意方位。每个节点有六个自由度，单元节点自由度为 单元刚度矩阵为 。在单元分析中综合考虑了拉压、扭转、弯曲、剪切变形。、空间膜单元 这是一个四节点四边形等参数单元，单元可以在三维空间中取任意方位，但四个节点
12、必须在同一平面内，单元处于平面应力状态。单元可以退化为常应变三角单元，也可以增加非协调项成为二维 Wilson 非协调单元，第六章中已经证明，当单元为矩形时可以通过分片检验，在这种情况下非协调项有助于改善单元内精度。单元材料可以是正交各向异性材料。图。、二维有限元 这是 平面内的四节点等参数单元。可以退化为常应变三角元。也可以加非协调项成为 Wilason 非协调元，本单元可以作为平面应变单元、或平面应力单元或轴对称单元使用图 66y,v*,u z,w i j 图*,u,*z,w,z i j 图 y,v,y 1212y,v*,u i l 图5 k j y,v*,u z,w i l 图4 k j
13、-.z.、八节点三维 Wilson 非协调单元 这种单元是在八节点三维等参数单元位移场中参加了九个非协调项，即 当单元为长方体时可以通过分片检验，非协调项有助于改善单元内的精度。本单元有 24 个外自由度 9 个内自由度，单元分析中先形成 的单元刚度矩阵，再凝聚掉 9个内自由度。以 的单元刚度矩阵进入总刚度矩阵图、板壳单元 本单元为LCCT9板单元或4 LCCT9个拼成的四边形单元，也可以用 LCCT9 平板单元来逼近壳体，本单元可以在空间取任何方位，用常应变三角元描述平面应力局部。用 LCCT9 单元描述弯曲局部。每个节点有六个自由度，但绕中面法线的转动自由度必须加以约束。单元的材料可以是正
14、交各向异性材料图、边界单元 这种单元实际是一根拉压或扭转弹簧。可以用来描述弹性支承。当弹簧刚度取很大数值时，还可以实现指定位移或斜支承约束条件。、八二十一节点三维等参数单元 本单元节点个数可在八二十一范围内选择，当节点个数取到二十时，单元可以有较好的精度和适应能力，单元的面也可以是曲面，组成单元的材料可以是正交各向异性材料图。、伪单元 这是一个由 用户自行定义的只 包括两个节点的单元。单元刚度矩阵的 元素 由 用 户 输入。利用这个单元可以 处理仅涉及两个节 点的约束关系，例 如光 滑 接 触 等 图。、读入单刚单元 8 7 6 5 4 3 2 1 图 33332424y,v,y*,u,*z,
15、w,z i j 图 j y,v,y*,u,*z,w,z i l j k o y,v*,u z,w 图8 图-.z.这也是一个由用户自行定义的单元，但节点个数不限，涉及到的自由度不超过，单元刚度矩阵的元素由用户输入，利用这个单元可以处理涉及到三个或三个以上节点的位移约束关系。、管单元 本单元可以是直管单元或圆弧形弯管单元。每个单元两个节点，十二个自由,管材料为各向同性材料。当使用这种单元时，规定：总体坐标的 轴必须铅垂向上。图 SAP5 程序原文本序号为的单元暂缺。引入我国后。根据实际需要有的单位寻单元库进展了扩大，目前，有的把四 九节点等参数单元作为序号为的单元，有的把八 十四节点混合刚度三维
16、元作为序号为的单元，尚未来得及加以统一。9-使 用SAP5 程序的方法 使用SAP5 程序的一般要经过以下几个主要步骤：、构造的离散化 将待分析的构造离散化为单元。一个构造可能要用到SAP5 程序的几种单元。每一种单元可编为一组。在各组内将单元从开场顺序编号。、节点编号和坐标 将所有节点从开场顺序加以编号，建立总体直角坐标系。总体坐标系轴的选择一般没有限制。但有些情况下例如使用管单元、三维 Wilson 非协调单元、轴对称单元时对总体坐标系的 y 轴或 z 轴有所限制。计算出各节点在总体坐标系中的坐标。、建立数据文件 将节点的坐标和约束信息、单元的几何、材料、载荷和拓朴信息。构造载荷信息及各种
17、控制数据按使用手册的规定准备好输入数据，并以批处理或会话方式建立数据文件。、数据的检查和修改 一般利用数据检查状态数据加以反复检查。改正其中的错误，并利用总体刚度矩阵的阶数和带宽对正式计算态所需要的 CPU 时间和外存容量进展估计。、正式计算 当确认原始数据无误后。即可转为正式计算状态。得到分析结果。参考文献 1 西阿莱，“有限元素法的数值分析，蒋尔雄等译，*科学技术出版1978。2 斯特朗，费克斯，“有限元法分析，*俊芝等译，科学1983。3 姜礼尚，庞之恒，“有限元方法及其理论根底。人民教育1979。4 陈传淼，“有限元方法及其提高精度的分析，*科学技术(1982)5 欧阳，马文华等，“弹
18、性 塑性 有限元,*科学技术(1983)6 RD Cook,“有限元分析的概念和应用，程耿东译，科学技术 7 SSRAD,“THE FINITE ELEMENT METHOD IN ENGINEERING(1980)8 RGlowinski,EY Rodin,OC Zienkiewicz,“Energy Methods in Finite Element Analysis,(1978)9 JEAKIN,“Application and Implementation of Finite Element Methods(1982)*,u,*z,w,z i j 图 y,v,y-.z.10 I霍拉德，
19、贝尔，“有限单元法 在应力分析中的应用，凌复华译，国防工业(1978)。11 董平，J罗赛托斯，“有限单元法根本方法与实施，*圣坤等译，国防工业,(1979)12 GPBazeley,YKCheung,BMIrons,CCZienkiewicz,“Triangular elements in plate bending-conforming and nonconforming solutions,Wright-patterson I(1965)13 HGSchaeffer,“A review of the international Symposium on structural mechanics Software,puters and Structures Vol.8 No.5 1978 14 J.N.Reddy:AN INTRODUCTION TO THE FINITE ELEMENT METHOD.McGraw-Hill Book pany.1984.15 王勖成，邵敏：“有限单元法根本原理和数值分析方法清华大学，1997 16 谭建国：使用 Ansys6.0 进展有限元分析大学，2002.5.17*鸿庆，任侠：构造有限元分析中国铁道，2000