关于有限元分析法及其应用举例 摘要：本文主要介绍有限元分析法，作为现代设计理论与方法的一种，已经在 众多领域普遍使用。介绍了它的起源和国内外发展现状。阐述了有限元法的基 本思想和设计方法。并从实际出发，例举了有限元法的一个简单应用———啤 酒瓶的应力分析和优化,表明了利用有限元分析法的众多优点。随着计算机的 发展，基于有限元分析方法的软件开发越来越多。本文也在其软件开发方面进 行阐述，并简单介绍了一下主流软件的发展情况和使用范围。并就这一领域的 未来发展趋势进行阐述。 关键词：有限元分析法软件啤酒瓶 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年代初兴起的应用数学、现代力学及计算机科学相互渗透、综合利用的边缘科学。有限元法最初应用在工程科学技术中，用于模拟并且解决工程力学、热学、电磁学等物理问题。对于过去用解析方法无法求解的问题和边界条件及结构形状都不规则的复杂问题，有限元法则是一种有效的分析方法。有限元法的基本思想是先将研究对象的连续求解区域离散为一组有限个且按一定方式相互联结在一起的单元组合体。由于单元能按不同的联结方式进行组合，且单元本身又可以有不同形状，因此可以模拟成不同几何形状的求解小区域；
有限元分析及应用作业报告
目录 有限元分析及应用作业报告....................................... I 目录 ........................................................ II 试题1 . (1) 一、问题描述 (1) 二、几何建模与分析 (2) 三、第1问的有限元建模及计算结果 (2) 四、第2问的有限元建模及计算结果 (7) 五、第3问的有限元建模及计算结果 (13) 六、总结和建议 (16) 试题5 (17) 一、问题的描述 (17) 二、几何建模与分析 (18) 三、有限元建模及计算结果分析 (18) 四、总结和建议 (26) 试题6 (27) 一、问题的描述 (27) 二、几何建模与分析 (27) 三、有限元建模及计算结果分析 (27) 五、总结和建议 (35)
试题1 一、问题描述 图示无限长刚性地基上的三角形大坝，受齐顶的水压力作用，试用三节点常应变单元和六节点三角形单元对坝体进行有限元分析，并对以下几种计算方案进行比较： 1)分别采用相同单元数目的三节点常应变单元和六节点三角形单元计算； 2)分别采用不同数量的三节点常应变单元计算； 3)当选常应变三角单元时，分别采用不同划分方案计算。 图1-1模型示意图及划分方案
二、几何建模与分析 图1-2力学模型 由于大坝长度>>横截面尺寸，且横截面沿长度方向保持不变，因此可将大坝看作无限长的实体模型，满足平面应变问题的几何条件；对截面进行受力分析，作用于大坝上的载荷平行于横截面且沿纵向方向均匀分布，两端面不受力，满足平面应变问题的载荷条件。因此该问题属于平面应变问题，大坝所受的载荷为面载荷，分布情况及方向如图1-2所示，建立几何模型，进行求解。 假设大坝的材料为钢，则其材料参数：弹性模量E=2.1e11，泊松比σ=0.3 三、第1问的有限元建模 本题将分别采用相同单元数目的三节点常应变单元和六节点三角形单元计算。1）设置计算类型：两者因几何条件和载荷条件均满足平面应变问题，故均取Preferences为Structural 2）选择单元类型：三节点常应变单元选择的类型是PLANE42（Quad 4node42），该单元属于是四节点单元类型，在网格划分时可以对节点数目控制使其蜕化为三节点单元；六节点三角形单元选择的类型是PLANE183（Quad 8node183），该单元属于是八节点单元类型，在网格划分时可以对节点数目控制使其蜕化为六节点单元。因研究的问题为平面应变问题，故对Element behavior（K3）设置为plane strain。 3）定义材料参数：按以上假设大坝材料为钢，设定：ANSYS Main Menu: Preprocessor →Material Props →Material Models →Structural →Linear →Elastic →Isotropic →input EX:2.1e11, PRXY:0.3 → OK 4）生成几何模型： a. 生成特征点：ANSYS Main Menu: Preprocessor →Modeling →Create →Keypoints→In Active CS→依次输入三个点的坐标：
有限元分析概念 有限元法：把求解区域看作由许多小的在节点处相互连接的单元（子域）所构成，其模型给出基本方程的分片（子域）近似解，由于单元（子域）可以被分割成各种形状和大小不同的尺寸，所以它能很好地适应复杂的几何形状、复杂的材料特性和复杂的边界条件 有限元模型：它是真实系统理想化的数学抽象。由一些简单形状的单元组成，单元之间通过节点连接，并承受一定载荷。 有限元分析：是利用数学近似的方法对真实物理系统（几何和载荷工况）进行模拟。并利用简单而又相互作用的元素，即单元，就可以用有限数量的未知量去逼近无限未知量的真实系统。 线弹性有限元是以理想弹性体为研究对象的，所考虑的变形建立在小变形假设的基础上。在这类问题中，材料的应力与应变呈线性关系，满足广义胡克定律；应力与应变也是线性关系，线弹性问题可归结为求解线性方程问题，所以只需要较少的计算时间。如果采用高效的代数方程组求解方法，也有助于降低有限元分析的时间。 线弹性有限元一般包括线弹性静力学分析与线弹性动力学分析两方面。 非线性问题与线弹性问题的区别： 1）非线性问题的方程是非线性的，一般需要迭代求解； 2）非线性问题不能采用叠加原理； 3）非线性问题不总有一致解，有时甚至没有解。 有限元求解非线性问题可分为以下三类：
1）材料非线性问题 材料的应力和应变是非线性的，但应力与应变却很微小，此时应变与位移呈线性关系，这类问题属于材料的非线性问题。由于从理论上还不能提供能普遍接受的本构关系，所以，一般材料的应力与应变之间的非线性关系要基于试验数据，有时非线性材料特性可用数学模型进行模拟，尽管这些模型总有他们的局限性。在工程实际中较为重要的材料非线性问题有：非线性弹性（包括分段线弹性）、弹塑性、粘塑性及蠕变等。 2）几何非线性问题 几何非线性问题是由于位移之间存在非线性关系引起的。 当物体的位移较大时，应变与位移的关系是非线性关系。研究这类问题一般都是假定材料的应力和应变呈线性关系。它包括大位移大应变及大位移小应变问题。如结构的弹性屈曲问题属于大位移小应变问题，橡胶部件形成过程为大应变问题。 3）非线性边界问题 在加工、密封、撞击等问题中，接触和摩擦的作用不可忽视，接触边界属于高度非线性边界。 平时遇到的一些接触问题，如齿轮传动、冲压成型、轧制成型、橡胶减振器、紧配合装配等，当一个结构与另一个结构或外部边界相接触时通常要考虑非线性边界条件。 实际的非线性可能同时出现上述两种或三种非线性问题。
第一章 绪论 有限元发展过程： 有限元法在西方起源于收音机和导弹的结构设计，发表这方面文章最早而且最有影响的是西德教授，于1954—1955年间分阶段在《Aircraft Engineering 》上发表上许多有关这方面的论文，并在此基础上写成了《能量原理与结构分析》，此书内容提供了有限元法的理论基础。美国的、 、 和等人于1956年发表了了篇题为《复杂结构的刚度和挠度分析》一文，此文提出了计算复杂结构刚度影响系数的方法，并说明了如何利用计算机进行分析。美国于1960年在一篇介绍平面应力分析的论文中，首先提出了有限元的名字。1965年英国及其合作者解决了将有限元法应用于所有场的问题，使有限元法的应用更加广泛。 有限元法的基本思路： 有限元法的基本思路和基本原理以结构力学中的位移法为基础，把复杂的结构或连续体看成为有限个单元的组合，各单元彼此在节点处连续而组成整体，把连续体分成有限个单元和节点，称之为离散化，先对单元进行特性分析，然后根据各单元在节点处的平衡协调条件建立方程，综合后作整体分析。 这样一分一合，先离散再综合的过程，就把复杂结构或连续体的计算问题转化为简单单元的分析与综合问题。 有限元分析中可采取三种方法： 位移法——取节点位移作为基本未知数 力 法——取节点力作为基本未知数 混合法—— 有限元法分析过程： 1、结构离散化（单元划分） 2、选择位移模式 为了能用节点位移表示单元体的位移、应变和应力，在分析连续体时，必须对单元中位移的分布做出一定的假定，也就是假定位移是坐标的某种简单函数，这种函数称为位移模式或位移函数（形函数）。 {}[]{}e u N δ= (1) 3、分析单元的力学特性 （1）利用几何方程：由位移表达式导出用点位移表示单元应变的关系式 {}[]{} e εδ=B {}ε为单元内任一点的应变列阵 (2) 非线性有限元 线性有限元 几何非线性 材料非线性 有限元
《有限元分析及应用》课程标准 课程代码：汽车学分：3 建议课时数：64 英文名称： 适用专业：计算机辅助设计与分析 先修课程：《计算机辅助设计》 课程团队负责人及成员：陈良萍、刘宏强、王云、赵静、李蕾、黄艺、史俊玲、 毛新 1.课程定位和设计思路 1.1课程定位 本课程是为计算机辅助设计与分析专业本科生开设的一门专业核心课程，重点介绍有限元法的基本原理和方法、一些成熟的有限元软件功能和简单的分析步骤，同时结合工程实际，为他们进一步学习或实际应用及参加科研工作开辟道路。其任务是通过先修课程中所学知识的综合运用和新知识的获取，使学生初步掌握现代设计中的一种重要方法，开阔视野，提高能力，以适应科学技术发展的要求。 1.2设计思路 在教学中，首先通过力学中的矩阵位移法思想的对比教学，引出连续介质力学有限单元法的学习重点在于单元的插值函数如何构造。这因为，虽说矩阵位移法是对杆系结构而言的，但其结构的离散化和组建整体刚度方程的思想完全可以借鉴到连续介质力学，它们的不同点只是在单元刚度矩阵的建立；而不同单元类型的单元刚度矩阵的建立，又取决于对应单元插值函数的构造。这样处理，不但使学生抓住了本课程的教学重点，而且对有限单元法的整体思想有了宏观上掌握；起到主动学习而非被动接受的作用。在单元构造的教学中，理论学习的重点在于常规单元的介绍；通过常规单元介绍插值函数的完备性与收敛性等。接之，介绍高次单元、等参单元等教学内容。在理论教学中，强调数学论证的严谨性和工程应用的适应性。
结合工程实例教学，拓宽学生数值分析方面的应用能力在课内对不同的单元类 型进行介绍时，及时抓住不同单元在应用中的对比教学与其适用性，并结合工程实例介绍单元类型的合理选取和单元网格的合理划分等。为学生在实际问题的数值分析中如何选定单元和剖分单元奠定了一定的基础和经验。 2.工作任务和课程目标 2.1工作任务 由于采用有限单元法的分析计算软件大多已商业化，而熟悉应用这些中的常规软件也应是本门课程的主要教学内容。在课内学生学会使用软件建立分析模型的基本步骤，其中包括分析模型抽象、几何模型绘制、单元网格划分、材料定义、边界条件定义、方程求解方法等。因课内教学时数的不足，学生应利用课余时间学习，以提高对实际问题的数值分析能力。 2.2课程目标 从教学思想和方法上对原课程进行改革，使学生从较高层次上理解有限元方法的实质，掌握有限元分析的工具，并具备初步处理工程问题的能力;使该课程成为具有较宽口径和较大覆盖面的、面向计算机辅助设计方面的专业基础课；注意课程体 系的整体优化，强调课程的深度、广度与应用。 3.教学方针落实情况
有限元法基本原理与应用 班级机械2081 姓名方志平 指导老师钟相强 摘要：有限元法的基础是变分原理和加权余量法，其基本求解思想是把计算域划分为有限个互不重叠的单元，在每个单元内，选择一些合适的节点作为求解函数的插值点，将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式，借助于变分原理或加权余量法，将微分方程离散求解。采用不同的权函数和插值函数形式，便构成不同的有限元方法。 关键词：有限元法；变分原理；加权余量法；函数。 Abstract：Finite element method is based on the variational principle and the weighted residual method, the basic idea is to solve the computational domain is divided into a finite number of non-overlapping units, each unit, select some appropriate function for solving the interpolation node points as , the differential variables rewritten or its derivative by the variable value of the selected node interpolation functions consisting of linear expressions, by means of variational principle or weighted residual method, the discrete differential equations to solve. Different forms of weight functions and interpolation functions, it constitutes a different finite element method. Keywords：Finite element method; variational principle; weighted residual method; function。 引言 有限元方法最早应用于结构力学，后来随着计算机的发展慢慢用于流体力学的数值模拟。在有限元方法中，把计算域离散剖分为有限个互不重叠且相互连接的单元，在每个单元内选择基函数，用单元基函数的线形组合来逼近单元中的真解，整个计算域上总体的基函数可以看为由每个单元基函数组成的，则整个计算域内的解可以看作是由所有单元上的近似解构成。在河道数值模拟中，常见的有限元计算方法是由变分法和加权余量法发展而来的里兹法和伽辽金法、最小二乘法等。根据所采用的权函数和插值函数的不同，有限元方法也分为多种计算格式。从权函数的选择来说，有配置法、矩量法、最小二乘法和伽辽金法，从计算单元网格的形状来划分，有三角形网格、四边形网格和多边形网格，从插值函数的精度来划分，又分为线性插值函数和高次插值函数等。不同的组合同样构成不同的有限元计算格式。对于权函数，伽辽金(Galerkin)法是将权函数取为逼近函数中的基函数；最小二乘法是令权函数等于余量本身，而内积的极小值则为对代求系数的平方误差最小；在配置法中，先在计
有限元分析及应用大作业 作业要求： 1）个人按上机指南步骤至少选择习题中3个习题独立完成，并将计算结果上交； 也可根据自己科研工作给出计算实例。 2）以小组为单位完成有限元分析计算； 3）以小组为单位编写计算分析报告； 4）计算分析报告应包括以下部分： A、问题描述及数学建模； B、有限元建模（单元选择、结点布置及规模、网格划分方案、载荷及边界 条件处理、求解控制） C、计算结果及结果分析（位移分析、应力分析、正确性分析评判） D、多方案计算比较（结点规模增减对精度的影响分析、单元改变对精度的 影响分析、不同网格划分方案对结果的影响分析等） 题一：图示无限长刚性地基上的三角形大坝，受齐顶的水压力作用，试用三节点常应变单元和六节点三角形单元对坝体进行有限元分析，并对以下几种计算方案进行比较： 1）分别采用相同单元数目的三节点常应变单元和六节点三角形单元计算；（注意ANSYS中用四边形单元退化为三节点三角形单元） 2）分别采用不同数量的三节点常应变单元计算； 3）当选常应变三角单元时，分别采用不同划分方案计算。 解：1.建模： 由于大坝长度>>横截面尺寸，且横截面沿长度方向保持不变，因此可将大坝看作无限长的实体模型，满足平面应变问题的几何条件；对截面进行受力分析，作
用于大坝上的载荷平行于横截面且沿纵向方向均匀分布，两端面不受力，满足平面应变问题的载荷条件。因此该问题属于平面应变问题，大坝所受的载荷为面载荷，分布情况P=98000-9800*Y；建立几何模型，进行求解；假设大坝的材料为钢，则其材料参数：弹性模量E=2.1e11，泊松比σ=0.3； 2：有限元建模过程： 2.1 进入ANSYS : 程序→ANSYS APDL 15.0 2.2设置计算类型: ANSYS Main Menu: Preferences →select Structural →OK 2.3选择单元类型: ANSYS Main Menu: Preprocessor →Element Type→Add/Edit/Delete →Add →select Solid Quad 4node 182(三节点常应变单元选择Solid Quad 4node 182，六节点三角形单元选择Solid Quad 8node 183)→OK (back to Element Types window) →Option →select K3: Plane Strain →OK→Close (the Element Type window) 2.4定义材料参数： ANSYS Main Menu: Preprocessor →Material Props →Material Models →Structural →Linear →Elastic →Isotropic →input EX:2.1e11, PRXY:0.3 →OK 2.5生成几何模型： 生成特征点： ANSYS Main Menu: Preprocessor →Modeling →Create →Keypoints→In Active CS →依次输入四个点的坐标：input:1(0,0),2(10,0),3(1,5),4(0.45,5) →OK 生成坝体截面： ANSYS Main Menu: Preprocessor →Modeling →Create →Areas →Arbitrary →Through KPS →依次连接四个特征点，1(0,0),2(6,0),3(0,10) →OK 2.6 网格划分： ANSYS Main Menu: Preprocessor →Meshing →Mesh Tool→(Size Controls) lines: Set →依次拾取两条直角边:OK→input NDIV: 15 →Apply→依次拾取斜边:OK →input NDIV: 20 →OK →(back to the mesh tool window)Mesh:Areas, Shape: tri, Mapped →Mesh →Pick All (in Picking Menu) →Close( the Mesh Tool window) 2.7 模型施加约束： 给底边施加x和y方向的约束： ANSYS Main Menu: Solution →Define Loads →Apply →Structural →Displacement →On lines →pick the lines →OK →select Lab2:UX, UY →OK 给竖直边施加y方向的分布载荷： ANSYS 命令菜单栏: Parameters →Functions →Define/Edit →1) 在下方的下拉列表框内选择x ,作为设置的变量；2) 在Result窗口中出现{X},写入所施加的载荷函数： 98000-9800*{Y}；3) File>Save(文件扩展名：func) →返回：Parameters →Functions →Read from file：将需要的.func文件打开，参数名取meng，它表示随之将施加的载荷→OK →ANSYS Main Menu: Solution →Define Loads →Apply →Structural →Pressure →On Lines →拾取竖直边；OK →在下拉列表框中，选择：Existing table →OK →选择需要的载荷为meng参数名→OK 2.8 分析计算： ANSYS Main Menu: Solution →Solve →Current LS →OK(to close the solve Current Load
第五章有限元素方法
§5.1有限元素方法的基本思想 有限元素法是一套求解微分方程的系统化数值计算方法。它比传统解法具有理论完整可靠，物理意义直观明确，适应性强，形式单纯、规范，解题效能强等优点。 从数学上来说, 有限元素方法是基于变分原理。它不象差分法那样直接去解偏微分方程, 而是求解一个泛函取极小值的变分问题。有限元素法是在变分原理的基础上吸收差分格式的思想发展起来的。 采用有限元素法还能使物理特性基本上被保持, 计算精度和收敛性进一步得到保证。 有限元素法优点： - 降低实验所需成本 - 減少試验对象的变异困难 - 方便参数控制 - 可获得实验无法获得的信息
有限元素法基本概念： 元素(element)，节点(node)，连結元素 有限元素法的基本思想: ?实际的物理問題很难利用单一的微分方程式描述，更无法順利求其解析解. ?有限元素法是将复杂的几何外型結构的物体切割成许多简单的几何形状称之为元素. ?元素与与元素间以“节点”相连. ?由于元素是简单的几何形状，故可以順利地写出元素的物理方程式,並求得节点上的物理量. ?采用內插法求得元素內任意点的物理量.
§5.2二维场的有限元素方法 1． 场域划分的约定 三角形元素。三角形元素越小，场域的分割就越细，计算的精度就会越高。因而在实际应用中是按精度的要求来决定场域内各处三角形元素的大小。 一般规定每个三角形元素的三个边的边长尽量地接近，尽量避免三角形元素具有大的钝角，一般最长的一条边不得大于最短边的三倍。 在分割场域时要求各三角形元素之间只能以顶点相交，即两相邻的三角形元素有两个公共的顶点及一条等长的公共边。不能把一个三角形的顶点取在另一个三角形的边上。 划分时还应当注意要尽量地使由相邻边界节点之间的线段所近似构成的曲线足够光滑。 如果在场域D内有不同的介质，则需要将介质的交面线选为分割线。
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 有限元法的思想 有限元法是求解连续介质力学问题的数值方法，更一般意义是一种分析结构问题和连续场数学物理问题的数值方法。 有限元法的基本思想是离散化和分片插值。 即把连续的几何机构离散成有限个单元，并在每一个单元中设定有限个节点，从而将连续体看作仅在节点处相连接的一组单元的集合体，同时选定场函数的节点值作为基本未知量并在每一单元中假设一个近似插值函数以表示单元中场函数的分布规律，再建立用于求解节点未知量的有限元方程组，从而将一个连续域中的无限自由度问题转化为离散域中的有限自由度问题。 求解得到节点值后就可以通过设定的插值函数确定单元上以至个集合体上的场函数。对每个单元，选取适当的插值函数，使得该函数在子域内部、在子域分界面上以及子域与外界面上都满足一定的条件。单元组合体在已知外载荷作用下处于平衡状态时，列出一系列以节点、位移为未知量的线性方程组，利用计算机解出节点位移后，再用弹性力学的有关公式，计算出各单元的应力、应变，当各单元小到一定程度，那么它就代表连续体各处的真实情况。
第9章受内外压筒体的有限元建模与应力变形分析(Project 2) 计算分析模型如图9-1 所示, 习题文件名: cylinder。 X (a) σO=100N/mm2 σI =200N/mm2 γ =7.85g/cm3 μ =0.3 E =210000N/mm2 (b) 图9-1 计算分析模型 9.1进入ANSYS 程序→ANSYSED 6.1ed →Interactive →change the working directory into yours→input Initial jobname: cylinder→Run 9.2 设置计算类型 ANSYS Main Menu: Preferences…→select Structural →OK 9.3 选择单元类型 ANSYS Main Menu: Preprocessor → Element Type →Add/Edit/Delete… → Add… →select Solid Quad 4node 42 →Apply →select Solid Brick 8node 45 → OK → Close (the Element
Types window) 9.4定义材料参数 ANSYS Main Menu: Preprocessor →Material Props →Materials Models →Structural→Lineal →Elastic→Isotropic…→input EX:2.1e5, PRXY:0.3→ OK 关闭材料定义窗口 9.5构造筒体模型 ?生成模型截平面 ANSYS Main Menu: Preprocessor →Modeling→Create →Keypoints →In Active CS… →按次序输入横截平面的十个特征点和旋转对称轴上两点坐标（十个特征点：(300,0,0), (480,0,0), (480,100,0), (400,100,0), (400,700,0), (480,700,0), (480,800,0), (300,800,0), (300,650,0), (300,150,0)，对称轴上两点：(0,0,0), (0,800,0)）（每次输入完毕，用Apply结束，0可以不输入） →Cancel (back to Create window) →-Areas- Arbitrary → Through KPs →依次连接截面边线上的十个特征点(注意在选完第10点后结束，不要再选第1点)→ OK ?对平面进行网格划分 ANSYS Main Menu: Preprocessor →Meshing→Mesh Tool →(Size Controls) Globl: Set →input SIZE (element edge length): 50 →OK (back to MeshTool window)→Mesh → Pick All (in Picking Menu) → Close( the MeshTool window) ?用旋转法生成筒体模型 ANSYS Main Menu: Preprocessor →Modeling→Operate →Extrude→Elem Ext Opts→select TYPE:SOLID 45→Element sizing options for extrusion No. Elem divs: 1→OK (back to Extrude window)→Areas →About Axis →Pick All(in Picking Menu)→OK→Pick the two keypoints (11,12) of the Symmetrical Axis → OK→input ARC: 90; NSEG: 3→ OK 9.6 模型加位移约束 ANSYS Main Menu: Solution→Define Loads →Apply→Structural→Displacement ?两截面分别加Z, X方向的约束 ANSYS Utility Menu: Select → Entities…→Nodes → By Location →select X coordinates →input 0→ OK (back to Displacement window)→On Nodes → Pick All(in Picking Menu) → select Lab2:UX →OK →ANSYS Utility Menu: Select → Everything ANSYS Utility Menu: Select → Entities…→ Nodes → By Location →select Z coordinates →input 0→ OK (back to Displacement window)→On Nodes →Pick All(in Picking Menu) → select Lab2:UZ →OK →ANSYS Utility Menu: Select →Everything ?底面加Y方向的约束 ANSYS Utility Menu: Select → Entities… → Nodes → By Location →select Y coordinates →input 0→ OK (back to Displacement window)→On Nodes →Pick All(in Picking Menu) →
有限元原理及工程应用 ——大作业 学院：机械工程学院 班级：硕4002班 小组成员：李追3114001089 陈草3114001080 2015.5.19
作业题目： 利用有限元方法对简支梁问题进行求解，梁的横截面为矩形，其约束情况如图1所示。 已知梁的几何尺寸和物理参数如下： （1）几何尺寸：长度40cm L =，截面尺寸2cm 02cm .b t ?=?； （2）物理参数：弹性模量70E =GPa ，泊松比0.3ν=，密度3 =2700kg/m ρ。 图1．梁及其横截面示意图 要求： (1) 至少划分五个节点（四个单元）； (2) 给出单元节点信息； (3) 给出单元刚度矩阵和质量矩阵； (4) 给出总刚度矩阵和总质量矩阵； (5) 求出梁各界固有频率及振型（五阶）； (6) 将所得结果与理论值进行对比，验证方法的可行性。
解：由有限元知识，根据Rayleigh-Ritz 法，解有限元分为四步：建立离散化、单元分析、形 成总体方程、解方程，具体步骤如下： （1）建立离散化 这里我们将矩形截面简支梁等分四等分，即分为六节点的五个杆单元，如图2所示： 每个单元尺寸40 cm=8cm 55 L l = =，这里只考虑杆在竖直平面的弯曲，每个节点只有y 方向位移和绕z 轴的旋转自由度。 （2）单元分析 构造一组Lagrange 插值基函数，在本节点值为1，其他节点值为0。从Rayleigh-Ritz 法可以看到，插值函数要p 次可微，最高阶导数出现在应变能表达式中；同样，我们可以这一原则适用于基函数的选择以及形状函数，否则我们将无法正确计算应变能当我们使用有限元逼近方法。 梁的弯曲问题，应变能计算公式： ()2 220,12L z v x t U EI dx x ???= ???? ? （1-1） 其中，E 为弹性模量，I z 为截面惯性矩。从公式可知，位移函数必须连续，并且二阶导数平方可积。 如图3，是一维杆单元模型，每个节点两个自由度，该单元含有四个自由度，即（,,,i zi j zj v v θθ）。本题中我们采用三次多项式插值函数： ()231234u x x x x αααα=+++ （1-2）
有限元方法理论及其应用
1 课程论文：弹性力学有限元位移法原理（30分） 撰写一篇论文，对有限元位移法的原理作一般性概括和论述。要求论文论及但不 限于下列内容：1）弹性力学有限元位移法的基本思想和数学、力学基础；2）有限元法求解的原理和过程，推导计算列式；对基本概念和矩阵符号进行解释和讨论；3）等参单元的概念、原理和应用。 1.1 对一维杆单元有限元形式的理解 将一维杆单元分成三段加以推导，并应用驻值条件0p D ?∏=?，我们得到节点的平衡 方程[K]{D}{R}=，即： 12 2341100112106012112600118u u AE cL u L u -?? ???? ?? ????--??????= ??????--??????????-???? ?? 我对此提出了几点疑问： 1) 为什么边界条件u 1=0，就要划去刚度矩阵[K]中对应的行列再解方程？ 2) 为什么刚度矩阵[K]会奇异？ 3) 为什么平衡方程本身是矛盾的，而加上边界条件u 1=0之后就能解出一个唯一的近似解？ 4) 为什么刚度矩阵[K]是对称的？ 下面我谈谈自己的理解:节点平衡方程是在u 1不定的前提下，假设单元内位移都是线性变化推导出来的，由此u 1相当于一个不确定的定值约束，再加上中间两个节点的连续性要求，系统实际上只有三个独立的自由度（广义坐标）。 对于第一个问题，其实刚度矩阵[K]中的元素不是一成不变的，相反它是伴随边界条件动态变化的。当u 1=0时由刚度矩阵的推导过程可以知道，刚度矩阵的第一行和第一列都会变为0，所以此时第一行和第一列对于求解方程是没有作用的。 对于第二个问题，由于系统自由度（广义坐标）只有三个，而我们的方程却列出了四个，显然这四个方程不可能线性无关，所以刚度矩阵奇异。
《有限元分析及应用》配书盘 曾攀 (清华大学机械工程系) 说明 该配书盘针对《有限元分析及应用》一书中有关有限元分析的自主程序开发、与ANSYS平台的衔接、基于ANSYS的有限元建模、基于MARC的有限元建模的章节，提供相应的电子材料及文档，以便在进行实际编程和应用国际著名商业软件进行建模和分析时参考。电子文档材料包括三大部分：(1)有限元分析源程序(f,c,ANSYS衔接)；(2) 四类问题有限元分析的操作指南(ANSYS,MARC)；(3) ANSYS一般性帮助文件。具体的文件目录和清单如下。 在目录/有限元分析源程序(f,c,ANSYS衔接)/中有以下内容 (1) 使用说明文件 自主程序开发使用说明(fortran,C,ANSYS平台衔接).pdf (2 ) 在子目录/fortran源程序及与ANSYS衔接(FEM2D)/中有以下文件 源程序文件： FEM2D.FOR 程序需读入的数据文件： BASIC.IN(模型的基本信息文件，需手工生成) NODE_ANSYS.IN (节点信息文件，可由ANSYS前处理导出，或手工生成) ELEMENT_ANSYS.IN(单元信息文件，可由ANSYS前处理导出，或手工生成)程序输出的数据文件： DATA.OUT (一般结果文件) FOR_POST.DAT(专供ANSYS进行后处理的结果数据文件) 与ANSYS后处理衔接的接口程序： USER_POST.LOG(在ANSYS中进行后处理的命令流文件) (3 ) 在子目录/c源程序及与ANSYS衔接(JIEKOU)/中有以下文件 源程序文件： JIEKOU.CPP 程序需读入的数据文件： NODE_ANSYS.IN(从ANSYS前处理导出的节点信息文件) ELEMENT_ANSYS.IN(从ANSYS前处理导出的单元信息文件) INPUT.DAT(包含除网格划分信息之外的所有前处理信息) 程序输出的数据文件：
有限元法理论及应用大作业 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）不能，位移函数应该包括常数项和一次项。
3.5 ANSYS软件加载、求解、后处理技术 3.5.1 ANSYS 3.5.1 ANSYS 荷载概述荷载概述 在这一节中将讨论: 有限元分析软件及应用 8 有限元分析软件及应用 8 A. 载荷分类 3.5 ANSYS 软件加载、求解、后处理技术 3.5 ANSYS 软件加载、求解、后处理技术 B. 加载 C. 节点坐标系 D. 校验载荷 孙瑛 孙瑛 E. 删除载荷 哈哈尔尔滨滨工工业业大学空大学空间结间结构研构研究中心究中心 2010秋 2010秋 SSRC SSRC 1/ 76 S Space pace S Stru truc ctu ture re R Res esear earc ch h C Center enter, H , HI IT, T, CH CHIN INA A
理技术 A. 载荷分类 B. 加载 A. 载荷分类 B. 加载 ANSYS中的载荷可分为: 可在实体模型或 FEA 模型节点和单元上加载自由度DOF - 定义节点的自由度( DOF )值结构分析_ 沿单元边界均布的压力 沿线均布的压力 位移集中载荷 - 点载荷结构分析_力面载荷 - 作用在表面的分布载荷结构分析_压力 在关键点处 在节点处约 约束体积载荷 - 作用在体积或场域内热分析_ 体积膨胀、内生 束 成热、电磁分析_ magnetic current density等实体模型 FEA 模型惯性载荷 - 结构质量或惯性引起的载荷重力、角速度等 在关键点加集中力在节点加集中力 SSR SSRC C SSR SSRC C 2/ 76 3/ 76 S Space pace S Stru truc ctu ture re R Res esear earc ch h C Center enter, H , HI IT, T, CH CHIN INA A S Space pace S Stru truc ctu ture re R Res esear earc ch h C Center enter, H , HI IT, T, CH CHIN INA A