结构有限元设计相关问题探索.pptVIP

结构有限元程序设计相关问题;有限元结构分析程序理论及模型化;问题的求解模型描述;有限元法的求解思路;方程系数 称为刚度系数，如写成矩阵形式，就是单元刚度阵，记作 ;对于平面内梁单元：有两个出口节点，每一个节点3个自由度，2个线变位，1个角变位。出口位移向量为 ;说明;;2 弹性理论的基本公式 ;在三维坐标系下，应变-位移关系为 ;将应力-应变、应变-位移关系带入平衡方程，不计初始应变，可以得到 ;力学中两个重要要原理;弹性力学的变分原理;弹性体总势能的具体表达式 ;3 假定位移场的有限元公式推导;设单元内任一点位移由下式给出 ;对于有 个单元的结构，结构的总势能为 ;结构的总势能成为 ;整理结构总势能的表达式得到 ;由最小总势能原理，处于平衡位置的位移 满足下面的方程组 ;4 结构模型化及程序模块划分;模型化的两方面内容;工程实际对结构分析程序功能要求;多目标程序的核心问题;有限元程序都的数据管理方式;有限元程序都的数据管理方式;IDT(IBSW+1)：文件的第1个字，其中存文件在内存中占有的长度LEN，以2字节为单位。 IDT(IBSW+2)：文件的第2个字，记为IPT。内存中前面IPT个字不是立即可用的资源；而后面LEN－IPT个字是立即可用的资源。 IDT(IBSW+3)：文件的第3个字，记为IV。表示文件尚有IV个字的资源可供使用。 IDT(IBSW+4)：文件的第4个字，记为ISTAT。冻结标记，当为10000时表示冻结。此时只可使用立即可用的资源。所有记录原位置不动。 IDT(IBSW+5)：文件的第5个字，为本文件的记录数NREC。 IDT(IBSW+6)：文件的第6个字，记为JFILBS。文件的库基址。由此可以找出本文件在库BANK中导向器存区的基址。 IDT(IBSW+7)：文件的第7个字，本文件导向器的长度，记为LDIREC。 说明：IBSW是本文件在内存IDT数组中存区的基址。所谓基址是指向存区之前的地址数。 ;有限元分析数据库（模型/载荷）;5 程序的单元分析模型;结构体系方面;有限元离散后，未知数一般是节点位移。由于结构构造方面的要求，往往有些单元刚度很大而其它一些单元刚度很小，相比之下这些刚度大的单元不能发生变形。因此这些刚度大的单元出口位移就不是相互独立的了（具有相关性）。结构的自由度与独立位移就成为一个重要问题。 复杂结构结构常常需要多种单元组合进行模拟，由于不同单元出口位移性态不同，如何处理不同单元之间（如梁单元和杆单元）的连接也是一个重要问题。有时不同结构部件相互连接并不是完全固接，其计算模型的选取也是需要处理的问题。 ;6.1 结构节点位移规格数;由于有不同的结构形式以及不同的边界条件或约束条件，节点任一个位移都可能有以下5种规格，分别以0，1，2，3，4的代码表示;上述节点位移规格数都是在结构全局坐标系中描写的，这个坐标系称作绝对坐标系。结构节点坐标是由绝对坐标系来定位的。 位移规格数的引入一方面可以处理非齐次节点位移，另一方面，它不仅能协助把指定为零及几何可动位移排除在总位移向量之外，还可以处理指定位移和相关位移，并也将其排除在总位移向量之外。 这对扩大解题范围、减少独立未知数，避免方程病态，从而提高求解精度有很大作用。;6.3主－从关系及其应用;新的解决方案;7 结构总刚度阵;可以将此划分为三个部分;总刚度阵的组织、未知数的编排 ;3.2 ;由于总刚度阵是对称正定阵，可以选择修正的平方根方法，或LDLT方法求解。首先将总刚度阵分解为 的形式，其中 是对角阵，而 则是对角元素等于1.0且上半矩阵全为零元素的下半三角阵。;可采用“一维变带宽”的方案进行总刚度阵的组织。即如图3－24中的矩阵，在对角线左面有部分区域有不为零的元素，但在左侧的边界线之左还有一区全为零元素。由于是对称矩阵。因此只要存储左半部分的元素就可以了。 一维变带宽矩阵的实际存储方式是将各行首尾相连，成为一个一维数组。该一维数组当然很长，常用的方法是用对角元地址数组IDIAG[1:ND]辨别每一行的分界线，如图3－25所示，这样可以将带外的零元素排除在外。 但这个一维数组的总刚度阵在内存一般是放不下的，因此需要放在BANK库中。则计算机整型字的限制32767可能不够对角元地址之用。在这里不采用IDIAG数组，而是存放每一行的半带宽（不包括对角元在内）。 ;;总未知数编排对于半带宽有重大影响;总刚度阵的装配 ;对于任一个单元，出口位移与结构总位移的关系可由下式给出 ;总刚度阵的求解;8 荷载模式与荷载工况;有限元模型介绍;工况与载荷介绍;载荷类型;工况类型;载荷 工况;这些荷载来源不一，但是在结构实际工作中却一般是同时作用的。对于结构的线性计算理论来说，完全可以对于各个单种因素的荷载进行求解；