1有限元分析概述

有限元分析非变形体（刚体）变形体1.1概述各力学学科的分支的关系（对象、变量、方程、求解途径）变形体科学研究的目的：定量获取所研究对象的所有信息确定对象定义参量或变量获得定量关系推广到该类问题的任意情形任意变形体力学分析的基本变量（信息量的表达方式）及方程有限元方法的思路及发展1.2有限元法的进展与应用有限元法不仅能应用于结构分析，还能解决归结为场问题的工程问题，从二十世纪六十年代中期以来，有限元法得到了巨大的发展，为工程设计和优化提供了有力的工具。1.2.1算法与有限元软件从二十世纪60年代中期以来，进行了大量的理论研究，不但拓展了有限元法的应用领域，还开发了许多通用或专用的有限元分析软件。理论研究的一个重要领域是计算方法的研究，主要有：大型线性方程组的解法，非线性问题的解法，动力问题计算方法。目前应用较多的通用有限元软件如下表所列：软件名称简介MSC/Nastran著名结构分析程序，最初由NASA研制MSC/Dytran动力学分析程序MSC/Marc非线性分析软件ANSYS通用结构分析软件ADINA非线性分析软件ABAQUS非线性分析软件1.2.2应用实例有限元法已经成功地应用在以下一些领域：固体力学，包括强度、稳定性、震动和瞬态问题的分析；传热学；电磁场；流体力学。结构与焊缝布置焊接残余应力分析（用Sysweld完成）焊接过程的温度分布与轴向残余应力应用案例一飞机液压泵振动强度数值模拟•项目简述某型歼击机液压泵传动轴频繁发生断裂事故，严重影响飞行安全。而该液压泵及传动轴的静强度和地面试车测试均没有任何问题，由于条件限制无法在飞机上进行动态测试。因此希望通过仿真计算解决问题•项目挑战–由37个零部件构成的装配体，模型复杂–系统刚度受转动和装配等因素的影响大–载荷谱非常复杂，计算点多，计算量大•解决方案–合理简化模型(如用只承压杆代替轴承)–通过改变接触刚度模拟装配刚度变化–传动轴疲劳寿命计算液压泵第4、5、11阶模态液压泵几何模型和有限元网格模型整体几何模型内部转动部分几何模型整体网格模型内部转动部分网格模型多段瞬态载荷谱采用子模型技术提高分析精度改进前后的疲劳寿命分布•项目描述–VDO电子助力转向系统–需要提高力矩传感器的灵敏度，以便更精确地执行驾驶员的转向要求•项目挑战–初始设计的扭转变形钢片几乎没有信号输出，无法实现扭矩传感•解决方案–通过结构分析发现原始设计的缺陷–第一次改进设计，效果很好，但由于结构尺寸过大，基本不实用–经过30多次方案改进，最后获得了一个非常满意的设计（传感器电路仿真也在ANSYS里一起完成）应用案例二电子助力转向系统传感器优化设计力矩传感器电子电机电子助力转向系统原始扭转变形钢片设计此处固定一个陶瓷厚膜电阻来“测量”扭转引起的应变，电阻值的变化就反映了扭矩的变化原始设计的缺陷是：在电阻贴片处的两个主应变大小几乎一样但方向相反，导致电阻的变化相互抵消，因此几乎没有信号输出！第一次改进设计的应变分布状态非常良好（基本上只有第一主应变，其它主应变很小），扭转引起的电阻变化很大，传感效果好。但结构宽度太大，无法集成在转向系统中，实用性差第一次改进设计最终设计专利产品VDO公司官方杂志封面•项目描述–核反应堆压力容器必须满足密封设计的变形限制要求，防止密封失效–结构的密封性能与结构形式、结构尺寸、材料性能及载荷等因素有着直接的关系，需要了解各因素对密封性能的影响程度及影响趋势，便于针对性地改进密封性能•项目挑战–需解决螺栓预紧、螺栓连接结构接触非线性分析、三维耦合热弹塑性接触密封分析等技术•解决方案–三维模型简化–首先进行螺栓预紧，然后施加内压、温度瞬态载荷以及压紧弹簧力、密封环反弹力等，同时考虑接触及接触传热，进行热弹塑性分析–应用概率设计功能，进行密封性能敏感性分析，确定对密封性能影响相对显著的因素应用案例三压力容器密封性能敏感性分析几何和网格模型螺栓应力对各输入参数的敏感性法兰内环槽变形对各输入参数的敏感性法兰外环槽变形对各输入参数的敏感性计算结论：设计中如螺栓应力计算结果超出规范限制，则可使用弹性模量更低的材料或减小螺栓预紧力以降低螺栓应力；如因法兰密封环槽分离量过大，已不能满足密封要求时，在不能降低压力的情况下可适当增大斜面倾角的角度•项目描述–用于F-15飞机的弹射座椅改进设计–需要计算在弹射和前向碰撞两种最大载荷状态下的座椅可靠性•项目挑战–100多个零部件，模型极其复杂–载荷施加非常困难•解决方案–在Workbench环境下使用Mechanical软件，利用其双向参数链接功能输入CAD模型，并自动创建零部件的装配接触–利用Workbench高级网格处理能力–利用Workbench先进的加载功能（如空间质量点、远程等效力等）–与CAD协同进行结构改进和优化设计应用案例四模块化弹射座椅系统级设计和分析座椅载荷示意图载荷和边界条件施加改进前后的结构应力分布座椅弹射状态结构强度计算座椅载荷示意图载荷和边界条件施加结构应力分布前向碰撞状态F-15弹射座椅从Pro/E传入Workbench的几何模型和Workbench划分的网格•项目描述–商用飞机机身大量板结构都承受剪切和挤压载荷，在极限载荷下的非线性屈曲使蒙皮和加强筋都发生塑性失效–用仿真分析代替部分实验指导设计•项目挑战–结构复杂，板壳连接形式多种多样–板壳厚度变化和中心偏移处理困难–材料类型多（金属板、多层板、纤维增强金属复合板等）•解决方案–应用APDL进行快速参数化建模–ANSYS先进板壳单元及其建模技术–利用接触功能模拟各种连接形式–采用复合材料、Hill各向异性材料等应用案例五空客飞机机身板壳结构极限载荷分析商用飞机上的板壳结构及其剪切/挤压屈曲失效实验机身板壳结构模型（铆接、点焊、粘连结构）变厚度板壳及其中心偏移连接结构细节处理不同载荷下的屈曲失效形态分析结果（右）与实验测试（左）的对比•项目描述–S波段多普勒气象雷达–天线支撑结构刚强度校核及改进•项目挑战–风载荷的准确获得及施加–结构连接形式复杂•解决方案–采用ANSYSCFX进行流体分析，获得风载荷压力分布数据–利用ANSYS流固耦合计算功能自动将风压值传入结构分析中–综合运用梁、壳、表面效应等单元建立结构模型，并采用自由度耦合模拟方向转动和俯仰转动轴承连接–进行结构静力、模态和谐响应分析应用案例六气象雷达天线流固耦合优化设计气象雷达结构ANSYSCFX流体分析模型及天线面上流体压力分布流体网格自适应流体分析的风压载荷自动传到结构模型上结构分析网格模型静力分析结果（位移和应力分布）模态分析和谐响应分析结果有限元分析的基本方法1）建立实际工程问题的计算模型利用几何、载荷的对称性简化模型建立等效模型2）选择适当的分析工具侧重考虑以下几个方面：物理场耦合问题大变形网格重划分3）前处理（Preprocessing）建立几何模型（GeometricModeling，自下而上，或基本单元组合）有限单元划分（Meshing）与网格控制有限元分析的基本方法4）求解（Solution）给定约束（Constraint）和载荷（Load）求解方法选择计算参数设定5）后处理（Postprocessing）后处理的目的在于分析计算模型是否合理，提出结论。用可视化方法（等值线、等值面、色块图）分析计算结果，包括位移、应力、应变、温度等；最大最小值分析；特殊部位分析。有限元法的基本思路有限元法的基本思路可以归结为：将连续系统分割成有限个分区或单元，对每个单元提出一个近似解，再将所有单元按标准方法组合成一个与原有系统近似的系统。1.3有限元法的计算步骤有限元法的计算步骤归纳为以下三个基本步骤：网格划分，单元分析，整体分析。网格划分有限元法的基础是用有限个单元体的集合来代替原有的连续体。因此首先要对弹性体进行必要的简化，再将弹性体划分为有限个单元组成的离散体。单元之间通过单元节点相连接。由单元、结点、结点连线构成的集合称为网格。通常把三维实体划分成4面体或6面体单元的网格，平面问题划分成三角形或四边形单元的网格。单元分析对于弹性力学问题，单元分析，就是建立各个单元的节点位移和节点力之间的关系式。由于将单元的节点位移作为基本变量，进行单元分析首先要为单元内部的位移确定一个近似表达式，然后计算单元的应变、应力，再建立单元中节点力与节点位移的关系式。整体分析对由各个单元组成的整体进行分析，建立节点外载荷与结点位移的关系，以解出结点位移，这个过程为整体分析。241102mA242101mAmLL1.021MPaEE521102NF1003用有限元法求图示受拉阶梯杆的位移和应力。已知杆截面面积，，各段杆长；材料弹性模量，作用于杆端的拉力。1.单元划分根据材料力学的平面假设，等截面受拉杆的同一截面的不同点上可认为具有相同的位移和应力，即位移只与截面的轴向坐标有关，所以可将阶梯杆看作由两个”一维单元”组成，同一单元内截面积及材料特性不变，并用线段表示一维单元。最简单的情况是，每一单元有两个结点，它们分别位于单元两端。相邻两单元由公共结点联接。因此，受拉杆简化为由两个一维单元和三个结点构成的有限元模型。如图b所示。图中①和②是单元号，1、2、3是结点号。取结点位移作为基本末知量，应力由求得的结点位移算出。2.确定单元插值函数（形函数）有限元法将整个求解域离散为一系列仅靠公共结点联接的单元，而每一单元本身却视为光滑的连续体。单元内任一点的场变量（如本算例中的位移）可由本单元的结点值根据场变量在单元中的假定分布规律（插值函数）插值求得。本例中，每个单元有两个结点，采用线性插值方式是适宜的。如图c是一典型单元图。两结点号分别为i和j。设单元中坐标为x处的场变量为xe，结点场变量值分别记为eiej和。根据线性插值关系得：eejeijiejeiEiejejijieiijjiijeiejeiexxLxxLxxxxxxxxxxxxxxx（1）形函数矩阵的分量数目应与单元自由度数相等。对于两自由度线性插值单元由式（1）可知形函数矩阵的两分量为eijejiLxxxLxxx(2)3.单元方程（单元结点位移与结点力的关系）由等截面杆变形与拉力的关系（虎克定律）得到：ejeiejeeeeiejeieeeFLEAFLEA（3）eiFejF式中，和分别为作用于单元e的结点i和结点j的结点力。式（3）写成矩阵形式为或简记为ejeiejeieeeFFLEA1111（4）eeeFK（5）4.单元组集------建立总体方程组ij为获得总体方程组，必须先将单元方程按照局部自由度（和）和总体自由度（Ф1Ф2Ф3）的对应关系进行扩展。即单元1的扩展方程为00000110112211321111FFLEA（6）式中，各项上角码表示单元序号；下角码表示自由度总体序号。单元2的扩展方程为23223212220110110000FFLEA（7）由于相邻单元公共结点上的基本场变量（位移）想同，所以可将扩展后的各单元方程组相加，将式（6）和（7）相加得32132122222222222211111111111