有限元法
精讲solidworks有限元分析步骤
2013-08-29 17:31 by:有限元来源:广州有道有限元 1. 软件形式： ㈠. SolidWorks的内置形式： ◆COSMOSXpress——只有对一些具有简单载荷和支撑类型的零件的静态分析。 ㈡. SolidWorks的插件形式： ◆COSMOSWorks Designer——对零件或装配体的静态分析。 ◆COSMOSWorks Professional——对零件或装配体的静态、热传导、扭曲、频率、掉落测试、优化、疲劳分析。 ◆COSMOSWorks Advanced Professional——在COSMOSWorks Professional的所有功能上增加了非线性和高级动力学分析。 ㈢. 单独发行形式： ◆COSMOS DesignSTAR——功能与COSMOSWorks Advanced Professional相同。 2. 使用FEA的一般步骤： FEA=Finite Element Analysis——是一种工程数值分析工具，但不是唯一的数值分析工具！其它的数值分析工具还有：有限差分法、边界元法、有限体积法… ①建立数学模型——有时，需要修改CAD几何模型以满足网格划分的需要， （即从CAD几何体→FEA几何体），共有下列三法： ▲特征消隐：指合并和消除在分析中认为不重要的几何特征，如外圆角、圆边、标志等。▲理想化：理想化是更具有积极意义的工作，如将一个薄壁模型用一个平面来代理（注：如果选中了“使用中面的壳网格”做为“网格类型”，COSMOSWorks会自动地创建曲面几何体）。 ▲清除：因为用于划分网格的几何模型必须满足比实体模型更高的要求。如模型中的细长面、多重实体、移动实体及其它质量问题会造成网格划分的困难甚至无法划分网格—这时我们可以使用CAD质量检查工具（即SW菜单: Tools→Check…）来检验问题所在，另外含有非常短的边或面、小的特征也必须清除掉（小特征是指其特征尺寸相对于整个模型尺寸非常小！但如果分析的目的是找出圆角附近的应力分布，那么此时非常小的内部圆角应该被保留）。 ②建立有限元模型——即FEA的预处理部分，包括五个步骤： ▲选择网格种类及定义分析类型（共有静态、热传导、频率…等八种类别）——这时将产生一个FEA算例，左侧浏览器中之算例名称之后的括号里是配置名称； ▲添加材料属性: 材料属性通常从材料库中选择，它不并考虑缺陷和表面条件等因素，与几何模型相比，它有更多的不确定性。
有限元理论方法
关于有限元分析法及其应用举例 摘要：本文主要介绍有限元分析法，作为现代设计理论与方法的一种，已经在 众多领域普遍使用。介绍了它的起源和国内外发展现状。阐述了有限元法的基 本思想和设计方法。并从实际出发，例举了有限元法的一个简单应用———啤 酒瓶的应力分析和优化,表明了利用有限元分析法的众多优点。随着计算机的 发展，基于有限元分析方法的软件开发越来越多。本文也在其软件开发方面进 行阐述，并简单介绍了一下主流软件的发展情况和使用范围。并就这一领域的 未来发展趋势进行阐述。 关键词：有限元分析法软件啤酒瓶 Abstract：This thesis mainly introduces the finite element analysis, as a modern design theory and methods used widely in in most respects. And this paper introduces its origins and development in world. It also expounds the basic thinking and approach of FEM..Proceed from the actual situation,this text holds the a simple application of finite-element method———the analysis and optimized of an beer bottle and indicate the the numerous benefits of finite element analysis .As computers mature and based on the finite element analysis of the software development is growing. This article introduces its application in the software development aspects as well, and briefly states the development and scope of the mainstream software. And it’s also prospect future development tendency in this area . Key: Finite Element Analysis Software Beer bottle 0 绪论 有限元法(Finite Element Method，FEM)，是计算力学中的一种重要的方法，它是20世纪50年代末60年代初兴起的应用数学、现代力学及计算机科学相互渗透、综合利用的边缘科学。有限元法最初应用在工程科学技术中，用于模拟并且解决工程力学、热学、电磁学等物理问题。对于过去用解析方法无法求解的问题和边界条件及结构形状都不规则的复杂问题，有限元法则是一种有效的分析方法。有限元法的基本思想是先将研究对象的连续求解区域离散为一组有限个且按一定方式相互联结在一起的单元组合体。由于单元能按不同的联结方式进行组合，且单元本身又可以有不同形状，因此可以模拟成不同几何形状的求解小区域；
有限元分析理论基础
有限元分析概念 有限元法：把求解区域看作由许多小的在节点处相互连接的单元（子域）所构成，其模型给出基本方程的分片（子域）近似解，由于单元（子域）可以被分割成各种形状和大小不同的尺寸，所以它能很好地适应复杂的几何形状、复杂的材料特性和复杂的边界条件 有限元模型：它是真实系统理想化的数学抽象。由一些简单形状的单元组成，单元之间通过节点连接，并承受一定载荷。 有限元分析：是利用数学近似的方法对真实物理系统（几何和载荷工况）进行模拟。并利用简单而又相互作用的元素，即单元，就可以用有限数量的未知量去逼近无限未知量的真实系统。 线弹性有限元是以理想弹性体为研究对象的，所考虑的变形建立在小变形假设的基础上。在这类问题中，材料的应力与应变呈线性关系，满足广义胡克定律；应力与应变也是线性关系，线弹性问题可归结为求解线性方程问题，所以只需要较少的计算时间。如果采用高效的代数方程组求解方法，也有助于降低有限元分析的时间。 线弹性有限元一般包括线弹性静力学分析与线弹性动力学分析两方面。 非线性问题与线弹性问题的区别： 1）非线性问题的方程是非线性的，一般需要迭代求解； 2）非线性问题不能采用叠加原理； 3）非线性问题不总有一致解，有时甚至没有解。 有限元求解非线性问题可分为以下三类：
1）材料非线性问题 材料的应力和应变是非线性的，但应力与应变却很微小，此时应变与位移呈线性关系，这类问题属于材料的非线性问题。由于从理论上还不能提供能普遍接受的本构关系，所以，一般材料的应力与应变之间的非线性关系要基于试验数据，有时非线性材料特性可用数学模型进行模拟，尽管这些模型总有他们的局限性。在工程实际中较为重要的材料非线性问题有：非线性弹性（包括分段线弹性）、弹塑性、粘塑性及蠕变等。 2）几何非线性问题 几何非线性问题是由于位移之间存在非线性关系引起的。 当物体的位移较大时，应变与位移的关系是非线性关系。研究这类问题一般都是假定材料的应力和应变呈线性关系。它包括大位移大应变及大位移小应变问题。如结构的弹性屈曲问题属于大位移小应变问题，橡胶部件形成过程为大应变问题。 3）非线性边界问题 在加工、密封、撞击等问题中，接触和摩擦的作用不可忽视，接触边界属于高度非线性边界。 平时遇到的一些接触问题，如齿轮传动、冲压成型、轧制成型、橡胶减振器、紧配合装配等，当一个结构与另一个结构或外部边界相接触时通常要考虑非线性边界条件。 实际的非线性可能同时出现上述两种或三种非线性问题。
有限元法的基本思想及计算 步骤
有限元法的基本思想及计算步骤 有限元法是把要分析的连续体假想地分割成有限个单元所组成的组合体，简称离散化。这些单元仅在顶角处相互联接，称这些联接点为结点。离散化的组合体与真实弹性体的区别在于：组合体中单元与单元之间的联接除了结点之外再无任何关联。但是这种联接要满足变形协调条件，即不能出现裂缝，也不允许发生重叠。显然，单元之间只能通过结点来传递内力。通过结点来传递的内力称为结点力，作用在结点上的荷载称为结点荷载。当连续体受到外力作用发生变形时，组成它的各个单元也将发生变形，因而各个结点要产生不同程度的位移，这种位移称为结点位移。在有限元中，常以结点位移作为基本未知量。并对每个单元根据分块近似的思想，假设一个简单的函数近似地表示单元内位移的分布规律，再利用力学理论中的变分原理或其他方法，建立结点力与位移之间的力学特性关系，得到一组以结点位移为未知量的代数方程，从而求解结点的位移分量。然后利用插值函数确定单元集合体上的场函数。显然，如果单元满足问题的收敛性要求，那么随着缩小单元的尺寸，增加求解区域内单元的数目，解的近似程度将不断改进，近似解最终将收敛于精确解。 用有限元法求解问题的计算步骤比较繁多，其中最主要的计算步骤为： 1）连续体离散化。首先，应根据连续体的形状选择最能完满地描述连续体形状的单元。常见的单元有：杆单元，梁单元，三角形单元，矩形单元，四边形单元，曲边四边形单元，四面体单元，六面体单元以及曲面六面体单元等等。其次，进行单元划分，单元划分完毕后，要将全部单元和结点按一定顺序编号，每个单元所受的荷载均按静力等效原理移植到结点上，并在位移受约束的结点上根据实际情况设置约束条件。 2）单元分析。所谓单元分析，就是建立各个单元的结点位移和结点力之间的关系式。现以三角形单元为例说明单元分析的过程。如图1所示，三角形有三个结点i，j，m。在平面问题中每个结点有两个位移分量u，v和两个结点力分量F x，F y。三个结点共六个结点位移分量可用列
有限元分析理论基础
有限元分析概念 有限元法:把求解区域瞧作由许多小的在节点处相互连接的单元(子域)所构成,其模型给出基本方程的分片(子域)近似解,由于单元(子域)可以被分割成各种形状与大小不同的尺寸,所以它能很好地适应复杂的几何形状、复杂的材料特性与复杂的边界条件 有限元模型:它就是真实系统理想化的数学抽象。由一些简单形状的单元组成,单元之间通过节点连接,并承受一定载荷。 有限元分析:就是利用数学近似的方法对真实物理系统(几何与载荷工况)进行模拟。并利用简单而又相互作用的元素,即单元,就可以用有限数量的未知量去逼近无限未知量的真实系统。 线弹性有限元就是以理想弹性体为研究对象的,所考虑的变形建立在小变形假设的基础上。在这类问题中,材料的应力与应变呈线性关系,满足广义胡克定律;应力与应变也就是线性关系,线弹性问题可归结为求解线性方程问题,所以只需要较少的计算时间。如果采用高效的代数方程组求解方法,也有助于降低有限元分析的时间。 线弹性有限元一般包括线弹性静力学分析与线弹性动力学分析两方面。 非线性问题与线弹性问题的区别: 1)非线性问题的方程就是非线性的,一般需要迭代求解; 2)非线性问题不能采用叠加原理; 3)非线性问题不总有一致解,有时甚至没有解。 有限元求解非线性问题可分为以下三类:
1)材料非线性问题 材料的应力与应变就是非线性的,但应力与应变却很微小,此时应变与位移呈线性关系,这类问题属于材料的非线性问题。由于从理论上还不能提供能普遍接受的本构关系,所以,一般材料的应力与应变之间的非线性关系要基于试验数据,有时非线性材料特性可用数学模型进行模拟,尽管这些模型总有她们的局限性。在工程实际中较为重要的材料非线性问题有:非线性弹性(包括分段线弹性)、弹塑性、粘塑性及蠕变等。 2)几何非线性问题 几何非线性问题就是由于位移之间存在非线性关系引起的。 当物体的位移较大时,应变与位移的关系就是非线性关系。研究这类问题一般都就是假定材料的应力与应变呈线性关系。它包括大位移大应变及大位移小应变问题。如结构的弹性屈曲问题属于大位移小应变问题,橡胶部件形成过程为大应变问题。 3)非线性边界问题 在加工、密封、撞击等问题中,接触与摩擦的作用不可忽视,接触边界属于高度非线性边界。 平时遇到的一些接触问题,如齿轮传动、冲压成型、轧制成型、橡胶减振器、紧配合装配等,当一个结构与另一个结构或外部边界相接触时通常要考虑非线性边界条件。 实际的非线性可能同时出现上述两种或三种非线性问题。 有限元理论基础
solidworks进行有限元分析的一般步骤
1.软件形式： ㈠. SolidWorks的内置形式： ◆COSMOSXpress——只有对一些具有简单载荷和支撑类型的零件的静态分析。 ㈡. SolidWorks的插件形式： ◆COSMOSWorks Designer——对零件或装配体的静态分析。 ◆COSMOSWorks Professional——对零件或装配体的静态、热传导、扭曲、频率、掉落测试、优化、疲劳分析。 ◆COSMOSWorks Advanced Professional——在COSMOSWorks Professional的所有功能上增加了非线性和高级动力学分析。 ㈢. 单独发行形式： ◆COSMOS DesignSTAR——功能与COSMOSWorks Advanced Professional相同。 2.使用FEA的一般步骤： FEA=Finite Element Analysis——是一种工程数值分析工具，但不是唯一的数值分析工具！其它的数值分析工具还有：有限差分法、边界元法、有限体积法… ①建立数学模型——有时，需要修改CAD几何模型以满足网格划分的需要， （即从CAD几何体→FEA几何体），共有下列三法： ▲特征消隐：指合并和消除在分析中认为不重要的几何特征，如外圆角、圆边、标志等。▲理想化：理想化是更具有积极意义的工作，如将一个薄壁模型用一个平面来代理（注：如果选中了“使用中面的壳网格”做为“网格类型”，COSMOSWorks会自动地创建曲面几何体）。▲清除：因为用于划分网格的几何模型必须满足比实体模型更高的要求。如模型中的细长面、多重实体、移动实体及其它质量问题会造成网格划分的困难甚至无法划分网格—这时我们可以使用CAD质量检查工具（即SW菜单: Tools→Check…）来检验问题所在，另外含有非常短的边或面、小的特征也必须清除掉（小特征是指其特征尺寸相对于整个模型尺寸非常小！但如果分析的目的是找出圆角附近的应力分布，那么此时非常小的内部圆角应该被保留）。 ②建立有限元模型——即FEA的预处理部分，包括五个步骤： ▲选择网格种类及定义分析类型（共有静态、热传导、频率…等八种类别）——这时将产生一个FEA算例，左侧浏览器中之算例名称之后的括号里是配置名称； ▲添加材料属性: 材料属性通常从材料库中选择，它不并考虑缺陷和表面条件等因素，与几何模型相比，它有更多的不确定性。 ◇右键单击“实体文件夹”并选择“应用材料到所有”——所有零部件将被赋予相同的材料属性。 ◇右键单击“实体文件夹”下的某个具体零件文件夹并选择“应用材料到所有实体”——某个零件的所有实体(多实体)将被赋予指定的材料属性。 ◇右键单击“实体文件夹”下具体零件的某个“Body”并选择“应用材料到实体”——只有
有限元软件介绍和比较
有限元软件介绍和比较 一、msc/patran+nastran, ansys, abaqus 三者的比较 俺最喜欢的是msc/patran+nastran，因为当年国内飞机公司最先引进的就是nastran，其菜单式的操作，比用手写有限元程序，爽多了！！特别是建立飞机这类巨大型结构，可以说，只有patran的建模最强！！（有人在仿真说abaqus能建整个飞机模型，哈哈，吹牛不上税，就凭其目前功能，要花一百年！！） 另外，msc财大气粗，其教程是手把手式，航空上最常用的有限元分析，都有现成的例题，step by step，傻瓜都会很快地入门！！由于其广泛应用于航空航天/汽车工业，所以，至今为止，如果要学CAE软件，俺认为应首选msc/patran+nastran。 与patran+nastran相比，ansys的界面就低了一些，操作也没有patran舒服。不过，差别不是很大。ansys据俺的体会，唯一的强项就是多场耦合。其他的功能， msc/patran+nastran都有。不过，ansys的apdl语言比较高级，是其最大优势，或者说，msc 应向这一方向发展！！不过，apdl最开始学也很费事，得一条一条查，一条一条记，这个过程没有两三个月下不来。由此，ansys的清爽度比msc差一些。 abaqus，如果自己用手编写过有限元程序的，入门应该不难。其命令格式，跟自己用手编程序一个套路。abaqus的强项是其分析功能很全面，特别是非线性部分，基本上都包含了。abaqus最大的缺点是上手慢，其教程太差，除了几本手册，基本上等于没有教程。要学abaqus，其时间要比msc, ansys长多了！！现在看，学abaqus实在没什么省时间的方法（比如它的 training lecture，一本250$，买来一看，气晕俺，还没手册说得详细！！），所以唯一的笨方法就是要看手册啦。(如果说msc是windows点鼠标时代的水平，abaqus就是敲dos命令的原始时代。不过，如果愣要用非线性分析，而nastran/ansys都没用，也只能用abaqus了。估计几年后，其CAE应能发展patran的水平，其教程应有step by step的水平。否则，为了一个非线性，多花数倍的时间，实在不爽！！或者说，花一辈子时间，才会用其中一部分功能，真可谓生也有涯，学也无涯，以有涯学无涯，不如不学算了！！ 二、MSC.PATRAN和ANSYS比较 MSC.PATRAN最早由美国宇航局(NASA)倡导开发的, 是工业领域最著名的并行框架式有限元前后处理及分析系统，其开放式、多功能的体系结构可将工程设计、工程分析、结果评估、用户化身和交互图形界面集于一身，构成一个完整 CAE集成环境。 ANSYS软件是融结构、流体、电场、磁场、声场分析于一体的大型通用有限元分析软件。由世界上最大的有限元分析软件公司之一的美国ANSYS开发，它能与多数CAD软件接口，实现数据的共享和交换，如Pro/Engineer, NASTRAN, Algor, I－DEAS, AutoCAD等，是现代产品设计中的高级CAD工具之一。 在建立复杂模型上ANSYS不如PATRAN，但PATRAN很繁琐。ANSYS比较适合于教学和科研，但ANSYS的求解效率确实不如NASTRAN。所以NASTRAN比较适合于工程。比较如下： 1、PATRAN界面层次分明，建模思路清晰；ANSYS界面菜单重叠、繁杂、互相覆盖，建模思路交替杂乱，条理不清。 2、PATRAN在一个界面内完成所有的同类模型（Geo. Fem BC. Mat. Prop.等各自为一类）操作。而ANSYS要重复打开和关闭多个相互重叠覆盖的界面，才能完成一个特征的创建和参数的输入等操作，非常烦琐。 3、PATRAN将计算任务提交给NASTRAN在后台运算后，在前台PATRAN仍然可以进行各种建模操作。而ANSYS提交了计算任务后，就不能再使用其前后处理功能。ANSYS的使用效率就大大地降低。
ANSYS 有限元分析基本流程
第一章实体建模 第一节基本知识 建模在ANSYS系统中包括广义与狭义两层含义，广义模型包括实体模型和在载荷与边界条件下的有限元模型，狭义则仅仅指建立的实体模型与有限元模型。建模的最终目的是获得正确的有限元网格模型，保证网格具有合理的单元形状，单元大小密度分布合理，以便施加边界条件和载荷，保证变形后仍具有合理的单元形状，场量分布描述清晰等。 一、实体造型简介 1．建立实体模型的两种途径 ①利用ANSYS自带的实体建模功能创建实体建模： ②利用ANSYS与其他软件接口导入其他二维或三维软件所建立的实体模型。 2．实体建模的三种方式 (1)自底向上的实体建模 由建立最低图元对象的点到最高图元对象的体，即先定义实体各顶点的关键点，再通过关键点连成线，然后由线组合成面，最后由面组合成体。 (2)自顶向下的实体建模 直接建立最高图元对象，其对应的较低图元面、线和关键点同时被创建。 (3)混合法自底向上和自顶向下的实体建模 可根据个人习惯采用混合法建模，但应该考虑要获得什么样的有限元模型，即在网格划分时采用自由网格划分或映射网格划分。自由网格划分时，实体模型的建立比较1e单，只要所有的面或体能接合成一体就可以：映射网格划分时，平面结构一定要四边形或三边形的面相接而成。 二、ANSYS的坐标系 ANSYS为用户提供了以下几种坐标系，每种都有其特定的用途。 ①全局坐标系与局部坐标系：用于定位几何对象(如节点、关键点等)的空间位置。 ②显示坐标系：定义了列出或显示几何对象的系统。 ③节点坐标系：定义每个节点的自由度方向和节点结果数据的方向。 ④单元坐标系：确定材料特性主轴和单元结果数据的方向。 1．全局坐标系 全局坐标系和局部坐标系是用来定位几何体。在默认状态下，建模操作时使用的坐标系是全局坐标系即笛卡尔坐标系。总体坐标系是一个绝对的参考系。ANSYS提供了4种全局坐标系：笛卡尔坐标系、柱坐标系、球坐标系、Y-柱坐标系。4种全局坐标系有相同的原点，且遵循右手定则，它们的坐标系识别号分别为：0是笛卡尔坐标系(cartesian)，1是柱坐标系 (Cyliadrical)，2是球坐标系(Spherical)，5是Y-柱坐标系(Y-aylindrical)，如图2-1所示。
有限元方法理论及其应用
1 课程论文：弹性力学有限元位移法原理（30分） 撰写一篇论文，对有限元位移法的原理作一般性概括和论述。要求论文论及但不限于下列内容：1）弹性力学有限元位移法的基本思想和数学、力学基础；2）有限元法求解的原理和过程，推导计算列式；对基本概念和矩阵符号进行解释和讨论；3）等参单元的概念、原理和应用。 1.1 对一维杆单元有限元形式的理解 我对此提出了几点疑问： 1)为什么边界条件u1=0，就要划去刚度矩阵[K]中对应的行列再解方程？ 2)为什么刚度矩阵[K]会奇异？ 3)为什么平衡方程本身是矛盾的，而加上边界条件u1=0之后就能解出一 个唯一的近似解？ 4)为什么刚度矩阵[K]是对称的？ 下面我谈谈自己的理解:节点平衡方程是在u1不定的前提下，假设单元内位移都是线性变化推导出来的，由此u1相当于一个不确定的定值约束，再加上中间两个节点的连续性要求，系统实际上只有三个独立的自由度（广义坐标）。 对于第一个问题，其实刚度矩阵[K]中的元素不是一成不变的，相反它是伴随边界条件动态变化的。当u1=0时由刚度矩阵的推导过程可以知道，刚度矩阵的第一行和第一列都会变为0，所以此时第一行和第一列对于求解方程是没有作用的。 对于第二个问题，由于系统自由度（广义坐标）只有三个，而我们的方程却列出
了四个，显然
这四个方程不可能线性无关，所以刚度矩阵奇异。 对于第三个问题，首先我们应该明确方程区别于等式，虽然左右两边都是用“=”连接，但是方程只在特殊条件下取得定解。由于平衡方程是在没有约束的条件下推导出来的，显然它不可能满足等式要求。宏观上看，系统在没有外部约束，而又施加有外力，显然系统会产生加速度而绝不会平衡。所以平衡方程本身是矛盾的。而加上边界条件之后，不但满足了平衡的前提，还改变了矩阵的结构和性质，所以有解。但是，由于我们提前假设了位移线性变化，相当于人为对单元施加了额外约束，让位移按照我们假设的规律变化，所以得到的解是过刚的近似解。但对于方程本身而言是精确解。 对于第四个问题，其力学的作用机理类似于作用力与反作用力，由于刚度矩阵不表征方向，所以其大小是相等的。 1.2 有限元法的思想 有限元法是求解连续介质力学问题的数值方法，更一般意义是一种分析结构问题和连续场数学物理问题的数值方法。 有限元法的基本思想是离散化和分片插值。 即把连续的几何机构离散成有限个单元，并在每一个单元中设定有限个节点，从而将连续体看作仅在节点处相连接的一组单元的集合体，同时选定场函数的节点值作为基本未知量并在每一单元中假设一个近似插值函数以表示单元中场函数的分布规律，再建立用于求解节点未知量的有限元方程组，从而将一个连续域中的无限自由度问题转化为离散域中的有限自由度问题。 求解得到节点值后就可以通过设定的插值函数确定单元上以至个集合体上的场函数。对每个单元，选取适当的插值函数，使得该函数在子域内部、在子域分界面上以及子域与外界面上都满足一定的条件。单元组合体在已知外载荷作用下处于平衡状态时，列出一系列以节点、位移为未知量的线性方程组，利用计算机解出节点位移后，再用弹性力学的有关公式，计算出各单元的应力、应变，当各单元小到一定程度，那么它就代表连续体各处的真实情况。
板结构有限元分析实例详解
板结构有限元分析实例详解1：带孔平板结构静力分析本节介绍带孔平板结构静力分析问题，同时介绍布尔操作的基本用法。 8.3.1 问题描述与分析 有孔的矩形平板，左侧边缘固定，长400mm，宽200 mm，厚度为10 mm，圆孔在板的正中心，半径为40 mm，左侧全约束，右侧边缘均布应力1MPa，如图8.7所示。求板的变形、位移及应力变化情况。（材料的材料属性为：弹性模量为300000 MPa，剪切模量为0.31。） 图8.7 带孔的矩形平板 由于小孔处边缘不规则，本文采用PLANE82高阶平面单元进行分析。 8.3.2 求解过程 8.3.2.1 定义工作目录及文件名 启动ANSYS Mechanical APDL Product Launcher窗口，如图8.8所示。在License下 拉选框中选择ANSYS Multiphysics产品，在Working Directory输入栏中输入工作目 录：C:\ANSYS12.0 Structural Finite Elements Analysis and Practice\Chapter 8\8-1，在Job Name一栏中输入工作文件名：Chapter8-1。以上参数设置完毕后，单 击Run按钮运行ANSYS。
图8.8 ANSYS设置窗口菜单 可以先在目标文件位置建立工作目录，然后单击Browse按钮选择工作目录；也 可以通过单击Browse按钮选择工作文件名。 8.3.2.2 定义单元类型和材料属性 选择Main Menu>Preferences命令，出现Preferences for GUI Filtering对话框， 如图8.9所示，在Individual discipline(s) to show in the GUI中勾选Structural，过滤掉ANSYS GUI菜单中与结构分析无关的选项，单击OK按钮关闭该对话框。 图8.9 Preferences for GUI Filtering对话框
有限元法的概述
有限元法的概述 有限元方法(Finite Element Method)是力学，数学物理学，计算方法，计算机技术等多种学科综合发展和结合的产物。在人类研究自然界的三大科学研究方法(理论分析，科学试验，科学计算)中，对于大多数新型领域，由于科学理论和科学实践的局限性，科学计算成为一种最重要的研究手段。在大多数工程研究领域，有限元方法是进行科学计算的重要方法之一;利用有限元方法几乎可以对任意复杂的工程结构进行分析，获取结构的各种机械性能信息，对工程结构进行评判，对工程事故进行分析。有限元法在设计过程中有极为关键的作用。 人们对各种力学问题进行分析求解，其方法归结起来可以分为解析法(Analytical Method)和数值法(Numeric Method).如果给定一个问题，通过一定的推导可以用具体的表达式来获得问题的解答，这样的求解方法就称为解析法。但是由于实际结构物的复杂性，除了少数极其简单的问题外，绝大多数科学研究和工程计算问题用解析法求解式极其困难的。因此，数值法求解便成为了一种不可替代的广泛应用的方法，并取得了不断的发展，如有限元法，有限差分法，边界元方法等都是属于数值求解方法。其中有限元法式 20 世纪中期伴随着计算机技术的发展而迅速发展起来的一种数值分析方法，它的数学逻辑严谨，物理概念清晰，应用非常广泛，能活灵活现处理和求解各种复杂的问题。有限元方法采用矩阵式来表达基本公式，便于计算机编程，这些优点赋予了它强大的生命力。 有限元方法的实质是将复杂的连续体划分成为有限多个简单的单元体，化无限自由度问题为优先自由度问题，将连续场函数的(偏)微分方程的求解问题转化为有限个参数的代数方程组的求解问题。用有限元方法分析工程结构的问题时，将一个理想体离散化后，如何保证其数值的收敛性和稳定性是有限元理论讨论的主要内容之一，而
有限元d 分析与介绍
有限元分析（FEA，Finite Element Analysis）利用数学近似的方法对真实物理系统（几何和载荷工况）进行模拟。还利用简单而又相互作用的元素，即单元，就可以用有限数量的未知量去逼近无限未知量的真实系统。 有限元分析是用较简单的问题代替复杂问题后再求解。它将求解域看成是由许多称为有限元的小的互连子域组成，对每一单元假定一个合适的(较简单的）近似解，然后推导求解这个域总的满足条件(如结构的平衡条件），从而得到问题的解。这个解不是准确解，而是近似解，因为实际问题被较简单的问题所代替。由于大多数实际问题难以得到准确解，而有限元不仅计算精度高，而且能适应各种复杂形状，因而成为行之有效的工程分析手段。 有限元是那些集合在一起能够表示实际连续域的离散单元。有限元的概念早在几个世纪前就已产生并得到了应用，例如用多边形（有限个直线单元）逼近圆来求得圆的周长，但作为一种方法而被提出，则是最近的事。有限元法最初被称为矩阵近似方法，应用于航空器的结构强度计算，并由于其方便性、实用性和有效性而引起从事力学研究的科学家的浓厚兴趣。经过短短数十年的努力，随着计算机技术的快速发展和普及，有限元方法迅速从结构工程强度分析计算扩展到几乎所有的科学技术领域，成为一种丰富多彩、应用广泛并且实用高效的数值分析方法。 有限元方法与其他求解边值问题近似方法的根本区别在于它的近似性仅限于相对小的子域中。20世纪60年代初首次提出结构力学计算有限元概念的克拉夫（Clough）教授形象地将其描绘为：“有限元法=Rayleigh Ritz法+分片函数”，即有限元法是Rayleigh Ritz法的一种局部化情况。不同于求解（往往是困难的）满足整个定义域边界条件的允许函数的Rayleigh Ritz法，有限元法将函数定义在简单几何形状（如二维问题中的三角形或任意四边形）的单元域上（分片函数），且不考虑整个定义域的复杂边界条件，这是有限元法优于其他近似方法的原因之一。
1有限元法简介
1有限元法简介 1.1有限单法的形成 在工程技术领域内，经常会遇到两类典型的问题。其中的第一类问题，可以归结为有限个已知单元体的组合。例如，材料力学中的连续梁、建筑结构框架和桁架结构。我们把这类问题，称为离散系统。如图1-1所示平面桁架结构，是由6个承受轴向力的“杆单元”组成。尽管离散系统是可解的，但是求解图1-2所示这类复杂的离散系统，要依靠计算机技术。 图1-1 平面桁架系统
图1-2 大型编钟“中华和钟”的振动分析及优化设计（曾攀教授） 第二类问题，通常可以建立它们应遵循的基本方程，即微分方程和相应的边界条件。例如弹性力学问题，热传导问题，电磁场问题等。由于建立基本方程所研究的对象通常是无限小的单元，这类问题称为连续系统。 图1-3 V6引擎的局部 下面是热传导问题的控制方程与换热边界条件： t T c Q z T z y T y x T x ??=+??? ??????+??? ? ??????+??? ??????ρλλλ （1- 1） 初始温度场也可以是不均匀的，但各点温度值是已知的： () 00 x,y,z T T t == （1- 2） 通常的热边界有三种，第三类边界条件如下形式： ()f T-T h n T λ=??- （1- 3） 尽管我们已经建立了连续系统的基本方程，由于边界条件的限制，通常只能得到少数简单问题的精确解答。对于许多实际的工程问题，还无法给出精确的解答，例如，图1-3所示V6引擎在工作中的温度分布。这为解决这个困难，工程师们和数学家们提出了许多近似方法。 在寻找连续系统求解方法的过程中，工程师和数学家从两个不同的路线得到了相同的结果，即有限元法。有限元法的形成可以回顾到二十世纪50年代，来源于固体力学中矩阵结构法的发展和工程师对结构相似性的直觉判断。从固体力学的角度来看，桁架结构等标准离散系统与人为地分割成有限个分区后的连续系统在结构上存在相似性。 1956年M..J.Turner, R.W.Clough, H.C.Martin, L.J.Topp 在纽约举行的航空学会年会上介
有限元法理论及应用参考答案分析
有限元法理论及应用大作业 1、试简要阐述有限元理论分析的基本步骤主要有哪些？ 答：有限元分析的主要步骤主要有： （1）结构的离散化，即单元的划分； （2）单元分析，包括选择位移模式、根据几何方程建立应变与位移的关系、根据虚功原理建立节点力与节点位移的关系，最后得到单元刚度方程； （3）等效节点载荷计算； （4）整体分析，建立整体刚度方程； （5）引入约束，求解整体平衡方程。 2、有限元网格划分的基本原则是什么？指出图示网格划分中不合理的地方。 题2图 答：一般选用三角形或四边形单元，在满足一定精度情况，尽可能少一些单元。 有限元划分网格的基本原则: 1.拓扑正确性原则。即单元间是靠单元顶点、或单元边、或单元面连接 2.几何保持原则。即网络划分后，单元的集合为原结构近似 3.特性一致原则。即材料相同，厚度相同 4.单元形状优良原则。单元边、角相差尽可能小 5.密度可控原则。即在保证一定精度的前提下，网格尽可能的稀疏一些。（a）(b)中节点没有有效的连接，且（b）中单元边差相差很大。 （c）中没有考虑对称性，单元边差很大。 3、分别指出图示平面结构划分为什么单元？有多少个节点？多少个自由度？
题3图 答：（a ）划分为杆单元， 8个节点，12个自由度。 （b ）划分为平面梁单元，8个节点，15个自由度。 （c ）平面四节点四边形单元，8个节点，13个自由度。 （d ）平面三角形单元，29个节点，38个自由度。 4、什么是等参数单元？。 答：如果坐标变换和位移插值采用相同的节点，并且单元的形状变换函数与位移插值的形函数一样，则称这种变换为等参变换，这样的单元称为等参单元。 5、在平面三节点三角形单元中，能否选取如下的位移模式，为什么？ (1). ?????++=++=2 65432 21),(),(y x y x v y x y x u αααααα (2). ?????++=++=2 65242 3221),(),(y xy x y x v y xy x y x u αααααα 答：（1）不能，因为位移函数要满足几何各向同性，即单元的位移分布不应与人为选取的 坐标方位有关，即位移函数中的坐标x,y 应该是能够互换的。所以位移多项式应按巴斯卡三角形来选择。 （2）不能，位移函数应该包括常数项和一次项。
有限元法的理论基础
有限元法的理论基础 有限元法是一种离散化的数值计算方法，对于结构分析而言，它的理论基础是能量原理。能量原理表明，在外力作用下，弹性体的变形、应力和外力之间的关系受能量原理的支配，能量原理与微分方程和定解条件是等价的。下面介绍有限元法中经常使用的虚位移原理和最小势能原理。 1.虚位移原理 虚位移原理又称虚功原理，可以叙述如下：如果物体在发生虚位移之前所受的力系是平衡的(物体内部满足平衡微分方程，物体边界上满足力学边界条件），那么在发生虚位移时，外力在虚位移上所做的虚功等于虚应变能（物体内部应力在虚应变上所做的虚功）。反之，如果物体所受的力系在虚位移（及虚应变）上所做的虚功相等，则它们一定是平衡的。可以看出，虚位移原理等价于平衡微分方程与力学边界条件。所以虚位移原理表述了力系平衡的必要而充分的条件。 虚位移原理不仅可以应用于弹性性力学问题，还可以应用于非线性弹性以及弹塑性等非线性问题。 2.最小势能原理 最小势能原理可以叙述为：弹性体受到外力作用时，在所有满足位移边界条件和变形协调条件的可以位移中，真实位移使系统的总势能取驻值，且为最小值。根据最小势能原理，要求弹性体在外力作用下的位移，可以满足几何方程和位移边界条件且使物体总势能取最小值的条件去寻求答案。最小势能原理仅适用于弹性力学问题。 2.2有限元法求解问题的基本步骤 弹性力学中的有限元法是一种数值计算方法，对于不同物理性质和数学模型的问题，有限元法的基本步骤是相同的，只是具体方式推导和运算求解不同，有限元求解问题的基本步骤如下。 2.2.1问题的分类 求解问题的第一步就是对它进行识别分析，它包含的更深层次的物理问题是什么？比如是静力学还是动力学，是否包含非线性，是否需要迭代求解，要从分析中得等到什么结果等。对这些问题的回答会加深对问题的认识与理解，直接影响到以后的建模与求解方法的选取等。 2.2.2建模 在进行有限元离散化和数值求解之值，我们为分析问题设计计算模型，这一步包括决定哪种特征是所要讨论的重点问题，以便忽略不必要的细节，并决定采用哪种理论或数学公式描述结果的行为。因此，我们可以忽略几何不规则性，把一些载荷看做是集中载荷，并把某些支撑看做是固定的。材料可以理想化为线弹性和各向同性的。根据问题的维数、载荷以及理论化的边界条件，我们能够决定采用梁理论、板弯曲理论、平面弹性理论或者一些其他分析理论描述结构性能。在求解中运用分析理论简化问题，建立问题的模型。 2.2.3连续体离散化 连续体离散化，习惯上称为有限元网络划分，即将连续体划分为有限个具有规则形状的单元的集合，两相邻单元之间只通过若干点相互连接，每个连接点称为节点。单元节点的设置、性质、数目等应视问题的性质、描述变形的需要和计算精度而定，如二维连续体的单元可为三角形、四边形，三维连续体的单元可以是四面体、长方体和六面体等。为合理有效地表示连续体，需要适当选择单元的类型、数目、大小和排列方式。 离散化的模型与原来模型区别在于，单元之间只通过节点相互连接、相互作用，而无其他连接。因此这种连接要满足变形协调条件。离散化是将一个无限多自由度的连续体转化为一个有限多自由度的离散体过程，因此必然引起误差。主要有两类：建模误差和离散化误差。
机械结构有限元分析
机械结构有限元分析 有限元分析软件ANSYS在机械设计中的应用 摘要：在机械设计中运用ANSYS软件进行有限元分析是今后机械设计发展的必然趋势，将有限元方法引入到机械设计课程教学中，让学生参与如何用有限元法来求解一些典型零件的应力，并将有限元结果与教材上的理论结果进行对照。这种新的教学方法可以大大提高学生的学习兴趣，增强学生对专业知识的理解和掌握，同时还可以培养学生的动手能力。在机械设计课程教学中具有很强的实用价值。 关键词：机械设计有限元 Ansys 前言：机械设计课程是一门专业基础课，其中很多教学内容都涉及到如何求取零件的应力问题，比如齿轮、v带、螺栓等零件。在传统的教学过程中，都是根据零件的具体受力情况按材料力学中相应的计算公式来求解。比如，在求解齿轮的接触应力时，是把齿轮啮合转化为两圆柱体的接触，再用公式求解。这些公式本身就比较复杂，还要引入各种修正参数，因此我们在学习这些内容时普遍反映公式难记，学习起来枯燥乏味，而且很吃力。 近年来有限元法在结构分析中应用越来越广泛，因此如果能将这种方法运用到机械设计课程中，求解一些典型零件的应力应变，并将分析结果和教材上的理论结果进行对比，那么无论是对于提高学生学习的热情和积极性，增强对重点、难点知识的理解程度，还是加强学生的计算机水平都是一件非常有益的事情。 由于直齿圆柱齿轮的接触强度计算是机械设计课程中的一个重要内容，齿轮强度的计算也是课程中工作量最繁琐的部分。下面就以渐开线直齿圆柱齿轮的齿根弯曲疲劳强度的计算为例，探讨在机械设计课程中用ANSYS软件进行计算机辅助教学的步骤和方法，简述如何将有限元方法应用到这门课程的教学中。 1．传统的直齿圆柱齿轮齿根弯曲疲劳强度的计算 传统方法把轮齿看作宽度为b的矩形截面的悬臂梁。因此齿根处为危险剖面，它可用30。切线法确定。如图l所示。 作与轮齿对称中心线成30。角并与齿根过渡曲线相切的切线，通过两切点作平行与齿轮轴线 的剖面，即齿根危险剖面。理论上载荷应由同时啮合的多对齿分担，但为简化计算，通常假设全部载荷作用于齿顶来进行分析，另用重合度系数E对齿根弯曲应力予以修正。 由材料力学弯曲应力计算方法求得齿根最大弯曲应力为：
(完整版)有限元法的基本原理
第二章有限元法的基本原理 有限元法吸取了有限差分法中的离散处理内核，又继承了变分计算中选择试探函数并对区域积分的合理方法。有限元法的理论基础是加权余量法和变分原理，因此这里首先介绍加权余量法和变分原理。 2.1等效积分形式与加权余量法 加权余量法的原理是基于微分方程等效积分的提法，同时它也是求解线性和非线性微分方程近似解的一种有效方法。在有限元分析中，加权余量法可以被用于建立有限元方程，但加权余量法本身又是一种独立的数值求解方法。 2.1.1 微分方程的等效积分形式 工程或物理学中的许多问题，通常是以未知场函数应满足的微分方程和边界条件的形式提出来的，可以一般地表示为未知函数u 应满足微分方程组 12()()()0A A A ?? ?== ? ??? u u u M （在Ω内） （2-1） 域Ω可以是体积域、面积域等，如图2-1所示。同时未知函数u 还应满足边界条件 12()()()0B B B ?? ?== ? ??? u u u M （在Γ内） （2-2） 要求解的未知函数u 可以是标量场（例如压力或温度），也可以是几个变量组成的向量场（例如位移、应变、应力等）。A ，B 是表示对于独立变量（例如空间坐标、时间坐标等）的微分算子。微分方程数目应和未知场函数的数目相对应，因此，上述微分方程可以是单个的方程，也可以是一组方程。所以在以上两式中采用了矩阵形式。 以二维稳态的热传导方程为例，其控制方程和定解条件如下： ()()()0A k k q x x y y φφφ????=++=???? （在Ω内） （2-3）
0()0q B k q n φφφφφ?-=Γ?=??-=Γ???（在上）（在上） （2-4） 这里φ表示温度（在渗流问题中对应压力）；k 是流度或热传导系数（在渗流问题中对应流度/K μ）；φ和q 是边界上温度和热流的给定值（在渗流问题中分别对应边界上的压力和边界上的流速）；n 是有关边界Γ的外法线方向；q 是源密度（在渗流问题中对应井的产量）。 在上述问题中，若k 和q 只是空间位置的函数时，问题是线性的。若k 和q 是φ及其导数的函数时，问题则是非线性的。 由于微分方程组（2-1）在域Ω中每一点都必须为零，因此就有 1122()(()())0u d v A u v A u d ΩΩ Ω≡++Ω≡? ?T V A L （2-5） 其中 12v V v ?? ?= ? ??? M （2-6） 其中V 是函数向量，它是一组和微分方程个数相等的任意函数。 式（2-5）是与微分方程组（2-1）完全等效的积分形式。我们可以说，若积分方程对于任意的V 都能成立，则微分方程（2-1）必然在域内任一点都得到满足。同理，假如边界条件（2-2）亦同时在边界上每一点都得到满足，对于一组任意函数，下式应当成立 1122 ()(()())0u d v B u v B u d ΓΓΓ≡++Γ≡??VB L 因此积分形式 ()()0u d u d ΓΓ Ω+Γ=??T T V A V B 对于所有的V 和V 都成立是等效于满足微分方程（2-1）和边界条件（2-2）。我们把（2-7）式称为微分方程的等效积分形式。 2.1.2等效积分的“弱”形式 在一般情况下，对（2-7）式进行分部积分得到另一种形式： ()()()()0T T v d v d ΩΓ Ω+Γ=??C D u E F u （2-8） 其中C ，D ，E ，F 是微分算子，它们中所包含的导数的阶数较（2-7）式的低，这样对函数u 只需要求较低阶的连续性就可以了。在（2-8）式中降低连续性要求是以提高V 和V 的连续性要求为代价的，由于原来对V 和V （在（2-7）式中）并无连续性要求，但是适当提高对其连续性的要求并不困难，因为它们是可以选择的已知函数。这种降低对函数u 连续性要求的作法在近似计算中，尤其是在有限单元法中是十分重要的。（2-8）式称为微分方程