有限元分析及应用作业报告
目录 有限元分析及应用作业报告....................................... 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→依次输入三个点的坐标：
关于有限元分析法及其应用举例 摘要：本文主要介绍有限元分析法，作为现代设计理论与方法的一种，已经在 众多领域普遍使用。介绍了它的起源和国内外发展现状。阐述了有限元法的基 本思想和设计方法。并从实际出发，例举了有限元法的一个简单应用———啤 酒瓶的应力分析和优化,表明了利用有限元分析法的众多优点。随着计算机的 发展，基于有限元分析方法的软件开发越来越多。本文也在其软件开发方面进 行阐述，并简单介绍了一下主流软件的发展情况和使用范围。并就这一领域的 未来发展趋势进行阐述。 关键词：有限元分析法软件啤酒瓶 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年代初兴起的应用数学、现代力学及计算机科学相互渗透、综合利用的边缘科学。有限元法最初应用在工程科学技术中，用于模拟并且解决工程力学、热学、电磁学等物理问题。对于过去用解析方法无法求解的问题和边界条件及结构形状都不规则的复杂问题，有限元法则是一种有效的分析方法。有限元法的基本思想是先将研究对象的连续求解区域离散为一组有限个且按一定方式相互联结在一起的单元组合体。由于单元能按不同的联结方式进行组合，且单元本身又可以有不同形状，因此可以模拟成不同几何形状的求解小区域；
有限元分析及其应用-2010 思考题： 1、有限元法的基本思想是什么？有限元法的基本步骤有那些？其中“离散”的含义是什 么？是如何将无限自由度问题转化为有限自由度问题的？ 答：基本思想：几何离散和分片插值。 基本步骤：结构离散、单元分析和整体分析。 离散的含义：用假想的线或面将连续物体分割成由有限个单元组成的集合，且单元之间仅在节点处连接，单元之间的作用仅由节点传递。当单元趋近无限小，节点无限多，则这种离散结构将趋近于实际的连续结构。 2、有限元法与经典的差分法、里兹法有何区别？ 区别：差分法：均匀离散求解域，差分代替微分，要求规则边界，几何形状复杂精度较低； 里兹法：根据描述问题的微分方程和相应的定解构造等价的泛函表达式，求得近似解； 有限元：基于变分法，采用分片近似进而逼近总体的求解微分方程的数值计算方法。 3、一根单位长度重量为q的悬挂直杆，上端固定，下端受垂直向下的外力P，试 1）建立其受拉伸的微分方程及边界条件； 2）构造其泛函形式； 3）基于有限元基本思想和泛函求极值构造其有限元的计算格式（即最小势能原理）。4、以简单实例为对象，分别按虚功原理和变分原理导出有限元法的基本格式（单元刚度矩 阵）。 5、什么是节点力和节点载荷？两者有何区别？ 答：节点力：单元与单元之间通过节点相互作用 节点载荷：作用于节点上的外载 6、单元刚度矩阵和整体刚度矩阵各有何特点？其中每个矩阵元素的物理意义是什么（按自 由度和节点解释）？ 答：单元刚度矩阵：对称性、奇异性、主对角线恒为正 整体刚度矩阵：对称性、奇异性、主对角线恒为正、稀疏性、带状性。 Kij，表示j节点产生单位位移、其他节点位移为零时作用i节点的力，节点力等于节点位移与单元刚度元素乘积之和。 7、单元的形函数具有什么特点？有哪些性质？ 答：形函数的特点：Ni为x，y的坐标函数，与位移函数有相同的阶次。 形函数Ni在i节点的值为1，而在其他节点上的值为0； 单元内任一点的形函数之和恒等于1； 形函数的值在0~1间变化。 8、描述弹性体的基本变量是什么？基本方程有哪些组成？ 答：基本变量：外力、应力、应变、位移 基本方程：平衡方程、几何方程、物理方程、几何条件 9、何谓应力、应变、位移的概念？应力与强度是什么关系？ 答：应力：lim△Q/△A=S △A→0 应变：物体形状的改变 位移：弹性体内质点位置的变化 10、问题的微分方程提法、等效积分提法和泛函变分提法之间有何关系？何谓“强形 式”？何谓“弱形式”，两者有何区别？建立弱形式的关键步骤是什么？
The finite element analysis Finite element method, the solving area is regarded as made up of many small in the node connected unit (a domain), the model gives the fundamental equation of sharding (sub-domain) approximation solution, due to the unit (a domain) can be divided into various shapes and sizes of different size, so it can well adapt to the complex geometry, complex material properties and complicated boundary conditions Finite element model: is it real system idealized mathematical abstractions. Is composed of some simple shapes of unit, unit connection through the node, and under a certain load. Finite element analysis: is the use of mathematical approximation method for real physical systems (geometry and loading conditions were simulated. And by using simple and interacting elements, namely unit, can use a limited number of unknown variables to approaching infinite unknown quantity of the real system. Linear elastic finite element method is a ideal elastic body as the research object, considering the deformation based on small deformation assumption of. In this kind of problem, the stress and strain of the material is linear relationship, meet the generalized hooke's law; Stress and strain is linear, linear elastic problem boils down to solving linear equations, so only need less computation time. If the efficient method of solving algebraic equations can also help reduce the duration of finite element analysis. Linear elastic finite element generally includes linear elastic statics analysis and linear elastic dynamics analysis from two aspects. The difference between the nonlinear problem and linear elastic problems: 1) nonlinear equation is nonlinear, and iteratively solving of general; 2) the nonlinear problem can't use superposition principle; 3) nonlinear problem is not there is always solution, sometimes even no solution. Finite element to solve the nonlinear problem can be divided into the following three categories: 1) material nonlinear problems of stress and strain is nonlinear, but the stress and strain is very small, a linear relationship between strain and displacement at this time, this kind of problem belongs to the material nonlinear problems. Due to theoretically also cannot provide the constitutive relation can be accepted, so, general nonlinear relations between stress and strain of the material based on the test data, sometimes, to simulate the nonlinear material properties available mathematical model though these models always have their limitations. More important material nonlinear problems in engineering practice are: nonlinear elastic (including piecewise linear elastic, elastic-plastic and viscoplastic, creep, etc. 2) geometric nonlinear geometric nonlinear problems are caused due to the nonlinear relationship between displacement. When the object the displacement is larger, the strain and displacement relationship is nonlinear relationship. Research on this kind of problem Is assumes that the material of stress and strain is linear relationship. It consists
有限元法基本原理与应用 班级机械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．简述有限单元法的应用范围 答：①工程地质现象机制的研究；②工程区岩体应力边界条件或区域构造力的反馈；③工程岩土体位移场和应力场的模拟；④岩土体稳定性模拟 2．简述有限元单元法的基本原理 答：有限元单元法是随着电子计算机的发展而迅速发展起来的一种现代计算方法。它是50年代首先在连续体力学领域----飞机结构静，动态特性分析中应用的一种由此奥的数分析方法，随后很快广泛的应用于求解热传导。电磁场、流体力学等连续性问题。有限元分析计算的思路和做法可归纳如下： ①物体离散化 将整个工程结构离散为由各个单元组成的计算模型，这一步称作单元剖分。离散散后单元与单元之间利用单元的节点相互连接起来；单元节点的设置、性质、树木等应是问题的性质，描述变形形态的需要和计算进度而定（一般情况但愿划分月息则描述变形情况月精确，及月接近实际变形，但计算两越大）。所以有限元中分析的结构已不是原有的物体或结构物，而是同新材料的由众多单元以一定方式连接成的离散物体。这样，用有限元分析计算所获得的结果只是近似的。如果划分单元数目非常多而又合理，则所获得的结果就与实际情况相符合。 ②单元特性分析 A．选择位移模式 在有限单元法中，选择节点位移为基本未知量称为位移法；选择节点力作为基本未知量时称为力法；取一部分节点力和一部分节点位移作为基本未知量时称为混合法。位移法易于实现计算机自动化，所以，在有限单元法中位移法应用范围最广。当采用位移法时，物体或结构离散化之后，就可把单元总的一些物理量如位移，应变和应力等由节点位移来表示。这时可以对单元中位移的分布采用一些能逼近原原函数的近似函数予以描述。通常，有限元法我们就将位移作为坐标变量的简单函数。这种函数称为位移模式或位移函数，如y=a其中a 是待定系数，y是与坐标有关的某种函数。 B．分析但愿的力学性质 根据单元的材料性质、形状、尺寸、节点数目、位置及其含义等，找出单元节点力和节点位移的关系式，折中单元分析中的关键一部。此时需要应用弹性力学中的几何方程和物理方程来来建立力和位移的方程式，从而导出单元刚度矩阵，这是有限元法的基本步骤之一。C．计算等效节点力 物体离散化后，假定力是通过节点从一个单元传递到另一个单元。但是，对于实际的连续题，力是从单元的公共边传递到另一个单元中去的。因而，这种作用在单元辩解上的表面力、体积力和集中力都需要等效的移动节点上去，也就是用等效的节点力来代替所有作用在单元上的力。 ③单元组集 利用结构力的平衡条件和边界条件把各个单元按原来的结构重新连接起来，形成整体的有限元方程 ④求解未知节点位移 解有限元方程式得出位移。这里，可以根据方程的具体特点来选择合适的计算方法。通过上述分析，可以看出，有限单元法的基本思想是“一分一合”，分是为了进行单元分
《有限元分析及应用》课程标准 课程代码：汽车学分：3 建议课时数：64 英文名称： 适用专业：计算机辅助设计与分析 先修课程：《计算机辅助设计》 课程团队负责人及成员：陈良萍、刘宏强、王云、赵静、李蕾、黄艺、史俊玲、 毛新 1.课程定位和设计思路 1.1课程定位 本课程是为计算机辅助设计与分析专业本科生开设的一门专业核心课程，重点介绍有限元法的基本原理和方法、一些成熟的有限元软件功能和简单的分析步骤，同时结合工程实际，为他们进一步学习或实际应用及参加科研工作开辟道路。其任务是通过先修课程中所学知识的综合运用和新知识的获取，使学生初步掌握现代设计中的一种重要方法，开阔视野，提高能力，以适应科学技术发展的要求。 1.2设计思路 在教学中，首先通过力学中的矩阵位移法思想的对比教学，引出连续介质力学有限单元法的学习重点在于单元的插值函数如何构造。这因为，虽说矩阵位移法是对杆系结构而言的，但其结构的离散化和组建整体刚度方程的思想完全可以借鉴到连续介质力学，它们的不同点只是在单元刚度矩阵的建立；而不同单元类型的单元刚度矩阵的建立，又取决于对应单元插值函数的构造。这样处理，不但使学生抓住了本课程的教学重点，而且对有限单元法的整体思想有了宏观上掌握；起到主动学习而非被动接受的作用。在单元构造的教学中，理论学习的重点在于常规单元的介绍；通过常规单元介绍插值函数的完备性与收敛性等。接之，介绍高次单元、等参单元等教学内容。在理论教学中，强调数学论证的严谨性和工程应用的适应性。
结合工程实例教学，拓宽学生数值分析方面的应用能力在课内对不同的单元类 型进行介绍时，及时抓住不同单元在应用中的对比教学与其适用性，并结合工程实例介绍单元类型的合理选取和单元网格的合理划分等。为学生在实际问题的数值分析中如何选定单元和剖分单元奠定了一定的基础和经验。 2.工作任务和课程目标 2.1工作任务 由于采用有限单元法的分析计算软件大多已商业化，而熟悉应用这些中的常规软件也应是本门课程的主要教学内容。在课内学生学会使用软件建立分析模型的基本步骤，其中包括分析模型抽象、几何模型绘制、单元网格划分、材料定义、边界条件定义、方程求解方法等。因课内教学时数的不足，学生应利用课余时间学习，以提高对实际问题的数值分析能力。 2.2课程目标 从教学思想和方法上对原课程进行改革，使学生从较高层次上理解有限元方法的实质，掌握有限元分析的工具，并具备初步处理工程问题的能力;使该课程成为具有较宽口径和较大覆盖面的、面向计算机辅助设计方面的专业基础课；注意课程体 系的整体优化，强调课程的深度、广度与应用。 3.教学方针落实情况
有限元分析及应用大作业 作业要求： 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
机械与汽车学院 曹国强 主要内容： 1、有限元法的基本思想。 2、结构力学模型的简化和结构离散化。 3、有限元法的实施过程。 一、有限元法的基本思想 有限元法是随着计算机的发展而发展起来的一种有效的数值方法。其基本思想是：将连续的结构分割成数目有限的小单元体（称为单元），这些小单元体彼此之间只在数目有限的指定点（称为节点）上相互连接。用这些小单元体组成的集合体来代替原来的连续结构。再把每个小单元体上实际作用的外载荷按弹性力学中的虚功等效原理分配到单元的节点上，构成等效节点力，并按结构实际约束情况决定受约束节点的约束。这一过程称为结构的离散化。其次，对每个小单元体选择一个简单的函数来近似地表示其位移分量的分布规律，并按弹性力学中的变分原理建立起单元节点力和节点位移之间的关系（单元刚度方程），最后，把全部单元的节点力和节点位移之间的关系组集起来，就得到了一组以结构节点位移为未知量的代数方程组（总体刚度方程），同时考虑结构的约束情况，消去那些结构节点位移为零的方程，再由最后的代数方程组就可求得结构上有限个离散节点的各位移分量。求得了结构上各节点的位移分量之后，即可按单元的几何方程和物理方程求得各单元的应变和应力分量。 有限元法的实质就是把具有无限个自由度的连续体，理想化为有限个自由度的单元的集合体，使问题简化为适合于数值解法的结构型问题。 经典解法（解析法）与有限元法的区别 解析法 { } 建立一个描述连续体性质的偏微分方程组 有限元解法 连续体 数目增加到∞ 大小趋于0 微元 有限元 离散化 （单元分析）集合 总体分析 求得近似解
二、结构力学模型的简化和结构离散化 （一）结构力学模型的简化 用有限元法研究实际工程结构问题时，首先要从工程实际问题中抽象出力学模型，即要对实际问题的边界条件、约束条件和外载荷进行简化，这种简化应尽可能地反映实际情况，不至于使简化后的解答与实际差别过大，但也不要带来计算上的过分复杂，在力学模型的简化过程中，必须判断实际结构的问题类型，是二维问题还是三维问题。如果是平面问题，是平面应力问题，还是平面应变问题。同时还要搞清楚结构是否对称，外载荷大小和作用位置，结构的几何尺寸和力学参数（弹性模量E、波松比μ等）。 （二）结构的离散化 将已经简化好的结构力学模型划分成只在一些节点连续的有限个单元，把每个单元看成是一个连续的小单元体，各单元之间只在一些点上互相联结，这些点称作节点，每个单元体称为一个单元。用只在节点处连接的单元的集合体代替原来的连续结构，把外载荷按虚功等效原理移置到有关受载的节点上，构成节点载荷，把连续结构进行这样分割的过程称为结构的离散化。现举例说明。 设一平面薄板，中间有一个园孔，其左端固定，右端受面力载荷q，试对其进行有限元分割和力学模型简化。
第1讲抛物问题有限元方法 1、椭圆问题有限元方法 考虑椭圆问题边值问题： （1） 问题（1）的变分形式：求使满足 （2） 的性质，广义解的正则性结果。 区域的剖分，矩形剖分，三角剖分，剖分规则，正则剖分条件，拟一致剖分条件。 剖分区域上分片次多项式构成的有限元空间。 的逼近性质，逆性质： 这里，为的插值逼近。 问题（2）的有限元近似：求使满足 （3） （3）的解唯一存在，且满足。 （3）的解所满足的矩阵方程（离散方程组）形式： （4） 刚度矩阵的由单元刚度矩阵组装而成。 模误差分析：由（2）－（3）可得 （5） 由（5）可首先得到 则得到 （6） －模误差分析 设满足 用与此方程做内积，由（5）式和插值逼近性质得到 再利用模误差估计结果，得到 （7） 最优阶误差估计和超收敛估计概念。 当与时间相关时（为抛物问题准备），由（5）式得 （8） 利用（7），类似分析可得 （9） 2、抛物问题半离散有限元方法 考虑抛物型方程初边值问题：
（10） （10）的变分形式：求使满足 （11） （11）的半离散有限元近似：求使满足 （12） 令，代入（12），依次取可导出常微分方程组： （13） 其中为质量矩阵，K为刚度矩阵。。 求解常微分方程组（13），得到代回的表达式，即得半离散有限元解。 定理1.问题（12）的解唯一存在且满足稳定性估计： （14） 证明：在（12）中取得到 整理为（注意是正定的） 对此式积分，证毕。 误差分析。引进解的椭圆投影逼近：满足 （15） 根据椭圆问题的有限元结果可知 （16） 分解误差： 的估计由（16）式给出，只须估计。 由（11），（12）和（15）知，满足 取，类似稳定性论证可得 （17） 可取为的投影，插值逼近等。 由（17）式，三角不等式和（16），得到 （18） 3、抛物问题全离散有限元近似 剖分时间区间：。 引进差分算子： 规定，当为连续函数时，，则有 由此得到 （19） （20） 定义问题（11）的全离散向后Euler有限元近似：求，使满足 （21） 将代入（21）可导出全离散方程组 （22）
第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）