2程序设计与软件应用.ppt

第2章 有限元程序设计方法 2.1 程序基本框图 1、输入基本数据（结构描述）： （1）控制数据：如结点总数、单元总数、约束条件总数等； （2）结点数据：如结点编号、结点坐标、约束条件等； （3）单元数据：如单元编号、单元结点序号、单元的材料特性、几何特性等； （4）载荷数据：包括集中载荷、分布载荷等。 2、单元分析 （1）各单元的bi,ci(i,j,m) ， 面积A； （2）应变矩阵[B]，应力矩阵[S]； （3）单元刚度矩阵[k]； （4）单元等价载荷列向量[F]。 总刚存贮 全矩阵存贮法：不利于节省计算机的存贮空间，很少采用。K[i,j] 对称三角存贮法：存贮上三角或下三角元素。 半带宽存贮法 ：存贮上三角形（或下三角形）半带宽以内的元素 。 一维压缩存贮法 ：半带宽存贮中仍包含了许多零元素。存贮每一行的第一个非零元素到主对角线元素。 例：计算下图半带宽。 4、引入约束条件 手算时采用去行列法 ，而计算机编程时采用乘大数法。即：指定结点位移对应的主对角元素乘上一个大数 ，同时将{P}中对应元素换为结点位移指定值与扩大了的主对角线元素的乘积 。 5、线性方程组求解 求解方法常用：GAUSS消元法，QR分解法等。其程序在一些专著中列出（例如见：徐士良编。FORTRAN常用算法程序集。清华大学出版社，1992）。在此不作详细介绍，其方法参阅[数值分析]有关书籍。 6、单元应力 节点位移求单元应力。首先整体节点位移变换成单元节点位移，然后再用物理方程求单元应力。 1、 Adina 麻省理工学院K.J. Bathe教授 ，用于结构、流体、流固耦合、热、热结构耦合分析，材料模型有金属、土壤与岩石、塑料、橡胶、织物、木材、陶瓷和混凝土等等。带裂纹传播的断裂力学 、ABS 防抱死系统分析、燃油泵及减震器分析、压电分析土壤固结、热机械分析 。 2、 ABAQUS 美国ABAQUS公司1978 推出，其特点是非线性功能强大，非线性求解效率和能力得到世界公认，单元有八类433种。特色：退火成型过程、海洋工程结构、水下冲击分析 、疲劳分析 3、ANSYS 美国ANSYS开发，集结构、流体、电场、磁场、声场分析及耦合分析于一体 。单元百余种。非线性功能较差，集成别家产品。 4、MSC/ NASTRAN 1966年美国国家航空航天局 （NASA）招标开发 ，它在航空航天领域有着崇高的地位。经过40余年的发展，用户从最初的航空航天领域，逐步发展到国防、汽车、造船、机械制造、兵器、铁道、电子、石化、能源材料工程、科研教育等各个领域，成为用户群最多、应用最为广泛的有限元分析软件。 5、SAP2000/ETABS 土木结构工程设计分析（含中国规范) 6、Midas CIVIL/Gen Civil桥梁为主, Midas/Gen 建筑 7、PKPM 结构分析与设计 2.3、ANSYS算例 网格划分 2.3 提高计算精度的方法 （1） 计算结果的整理 计算结果包括位移和应力两个方面。在位移方面，一般无须进行整理工作。应力结果则需要整理。通常认为计算出的应力是三角形单元形心处的应力。而相邻单元之间的应力存在突变，甚至正、负符号都不相同。为了由计算结果推算出结构内某一点的接接实际的应力，必须通过某种平均计算。通常可采用两单元平均法或绕结点平均法。 平均法整理单元应力 两单元平均法：把两个相邻单元中的常应力加以平均，用来表示公共边界中点处的应力。 绕结点平均法：把环绕某一结点的各单元常应力加以平均，用以表示该结点的应力。在内结点效果较好，而在边界结点可能很差，一般改为应由内结点的应力外推计算出来。 （2）网格的细分 通过网格的细分，使每个单元的面积缩小，那么尽管每个单元是应变、常应力单元，仍可较好地反映结构中的应力变化，使得到的解答收敛于问题的精确解。 （3）网格合理布局 根据应力梯度使网格的布局合理化。即在梯度大的区域网格密些，梯度小的区域应稀些。密、稀网格之间应逐步过渡。 （4）改用高阶单元 受集中力的悬臂梁，采用128个三结点三角形常应变单元，以及3个八结点四边形高阶单元结果。由图可见，采用高阶元的计算精度比常应变元高得多。 * 开始 输入基本数据 计算单元刚度矩阵 形成总体刚度矩阵 形成结点荷载向量 引入约束条件 求解方程组，输出结点位移 计算单元应力，输出结果 结束 开始 输入基本数据 计算单元刚度矩阵 形成总体刚度矩阵 形成结点荷载向量 引入约束条件 求解方程组，输出结点位移 计算单元应力，输出结果 结束 3、系统分析 （1）整体刚度矩阵[K]的组装； （2）整体载荷列阵{P}的形成； [