有限元知识点归纳 1.、有限元解的特点、原因？ 答：有限元解一般偏小，即位移解下限性 原因：单元原是连续体的一部分，具有无限多个自由度。在假定了单元的位移函数后，自由度限制为只有以节点位移表示的有限自由度，即位移函数对单元的变形进行了约束和限制，使单元的刚度较实际连续体加强了，因此，连续体的整体刚度随之增加，离散后的刚度较实际的刚度K为大，因此求得的位移近似解总体上将小于精确解。 2、形函数收敛准则（写出某种单元的形函数，并讨论收敛性）P49 (1)在节点i处N i=1,其它节点N i=0; (2)在单元之间，必须使由其定义的未知量连续； (3)应包含完全一次多项式； (4)应满足∑Ni=1 以上条件是使单元满足收敛条件所必须得。可以推证，由满足以上条件的形函数所建单元是完备协调的单元，所以一定是收敛的。 4、等参元的概念、特点、用时注意什么？（王勖成P131） 答：等参元—为了将局部坐标中几何形状规则的单元转换成总体（笛卡尔）坐标中的几何形状扭曲的单元，以满足对一般形状求解域进行离散化的需要，必须建立一个坐标变换。即： 为建立上述的变换，最方便的方法是将上式表示成插值函数的形式，即： 其中m是用以进行坐标变换的单元节点数，xi,yi,zi是这些结点在总体（笛卡尔）坐标内的坐标值，Ni’称为形状函数，实际上它也是局部坐标表示的插值函数。称前者为母单元，后者为子单元。 还可以看到坐标变换关系式和函数插值表示式：在形式上是相同的。如果坐标变换和函数插值采用相同的结点，并且采用相同的插值函数，即m=n，Ni’=Ni,则称这种变换为等参变换。 5、单元离散？P42 答：离散化既是将连续体用假想的线或面分割成有限个部分，各部分之间用有限个点相连。每个部分称为一个单元，连接点称为结点。对于平面问题，最简单、最常用的离散方式是将其分解成有限个三角形单元，单元之间在三角形顶点上相连。这种单元称为常应变三角形单元。常用的单元离散有三节点三角形单元、六节点三角形单元、四节点四边形单元、八节点四边形单元以及等参元。 6、数值积分，阶次选择的基本要求？ 答：通常是选用高斯积分 积分阶次的选择—采用数值积分代替精确积分时，积分阶数的选取应适当，因为它直接影响计算精度，计算工作量。选择时主要从两方面考虑。一是要保证积分的精度，不损失收敛性；二是要避免引起结构总刚度矩阵的奇异性，导致计算的失败。
有限元分析概念 有限元法：把求解区域看作由许多小的在节点处相互连接的单元（子域）所构成，其模型给出基本方程的分片（子域）近似解，由于单元（子域）可以被分割成各种形状和大小不同的尺寸，所以它能很好地适应复杂的几何形状、复杂的材料特性和复杂的边界条件 有限元模型：它是真实系统理想化的数学抽象。由一些简单形状的单元组成，单元之间通过节点连接，并承受一定载荷。 有限元分析：是利用数学近似的方法对真实物理系统（几何和载荷工况）进行模拟。并利用简单而又相互作用的元素，即单元，就可以用有限数量的未知量去逼近无限未知量的真实系统。 线弹性有限元是以理想弹性体为研究对象的，所考虑的变形建立在小变形假设的基础上。在这类问题中，材料的应力与应变呈线性关系，满足广义胡克定律；应力与应变也是线性关系，线弹性问题可归结为求解线性方程问题，所以只需要较少的计算时间。如果采用高效的代数方程组求解方法，也有助于降低有限元分析的时间。 线弹性有限元一般包括线弹性静力学分析与线弹性动力学分析两方面。 非线性问题与线弹性问题的区别： 1）非线性问题的方程是非线性的，一般需要迭代求解； 2）非线性问题不能采用叠加原理； 3）非线性问题不总有一致解，有时甚至没有解。 有限元求解非线性问题可分为以下三类：
1）材料非线性问题 材料的应力和应变是非线性的，但应力与应变却很微小，此时应变与位移呈线性关系，这类问题属于材料的非线性问题。由于从理论上还不能提供能普遍接受的本构关系，所以，一般材料的应力与应变之间的非线性关系要基于试验数据，有时非线性材料特性可用数学模型进行模拟，尽管这些模型总有他们的局限性。在工程实际中较为重要的材料非线性问题有：非线性弹性（包括分段线弹性）、弹塑性、粘塑性及蠕变等。 2）几何非线性问题 几何非线性问题是由于位移之间存在非线性关系引起的。 当物体的位移较大时，应变与位移的关系是非线性关系。研究这类问题一般都是假定材料的应力和应变呈线性关系。它包括大位移大应变及大位移小应变问题。如结构的弹性屈曲问题属于大位移小应变问题，橡胶部件形成过程为大应变问题。 3）非线性边界问题 在加工、密封、撞击等问题中，接触和摩擦的作用不可忽视，接触边界属于高度非线性边界。 平时遇到的一些接触问题，如齿轮传动、冲压成型、轧制成型、橡胶减振器、紧配合装配等，当一个结构与另一个结构或外部边界相接触时通常要考虑非线性边界条件。 实际的非线性可能同时出现上述两种或三种非线性问题。
有限元知识点汇总 第一章 1、何为有限元法？其基本思想是什么？ 》有限元法是一种基于变分法而发展起来的求解微分方程的数值计算方法。 》基本思想：化整为零，化零为整 2、为什么说有限元法是近似的方法，体现在哪里？ 》有限元法的基本思想是几何离散和分片插值； 》用离散单元的组合来逼近原始结构，体现了几何上的近似；用近似函数逼近未知量在单元内的真实解，体现了数学上的近似；利用与问题的等效的变分原理建立有限元基本方程，又体现了明确的物理背景。 3、单元、节点的概念？ 》单元：把参数单元划分成网格，这些网格就称为单元。 》节点：网格间相互连接的点称为节点。 4、有限元法分析过程可归纳为几个步骤？ 》3大步骤；——结构离散化；——单元分析；——整体分析。 5、有限元方法分几种？本课程讲授的是哪一种？ 》有限元方法分3种；——位移法、力法、混合法。 》本课程讲授的：位移法 6、弹性力学的基本变量是什么？何为几何方程、物理方程及虚功方程？弹性矩阵的特点？》弹性力学的基本变量是——{外力、应力、应变、位移} 》几何方程——{描述弹性体应变分量与位移分量之间关系的方程} 》物理方程——{描述应力分量与应变分量之间的关系} 》虚功方程——{描述内力和外力的关系的方程} 》弹性矩阵特点——{ } 7、何为平面应力问题和平面应变问题？ 》平面应力问题——{满足（1）几何条件——所研究的是一根很薄的等厚度薄板，即一个方向上的几何尺寸远远小于其余两个面上的几何尺寸;（2）载荷条件——作用于薄板上的载荷平行于板平面且沿厚度方向均匀分布，而在两板面上无外力作用} 》平面应变问题——{满足（1）几何条件——所研究的是长柱体，即长度方向的尺寸远远大于横截面的尺寸，且横截面沿长度方向不变；（2）载荷条件——作用于长柱体结构上的载荷平行于横截面且沿纵向方向均匀分布，两端面不受力} 第二章 7、形函数的特点？ 》1形函数Ni再节点i处等于1，在其他节点上的值等于0，对于Nj、Nm也有同样的性质。》2在单元内任一点的各形函数之和等于1，即Ni+Nj+Nm=1 8、单元刚度矩阵的性质？ 》1 K^e中每个元素都有明确的物理意义，每个元素都是一个刚度系数，他是单位节点位移分量所引起的节点力分量 》2 k^e是对称矩阵，具有对称性。 》3 K^e的每一行或每一列元素之和为零，是奇异矩阵
有限元方法基础吉林大学汽车工程学院
引论 一、结构分析的目的 1、克服传统设计方法的不足。 2、为改进设计、优化设计提供指导。 以往的设计大都是基于经验的，尽管基于经验的设计方法在以往的设计中取得了巨大的成功，并且今后的设计仍然将以经验设计为主。但也必须承认基于经验的设计也存在许多不足。首先经验的积累需要时间，有时也不一定都可靠。 此外，通常基于经验的设计能解决行不行的问题，很难解决优不优的问题，尽管为了解决结构优化的问题，也有人提出了一些新的设计理念，如加法设计原则。解决上述问题的根本手段就是采用合理的、科学的方法对上述性能进行分析。
二、数值分析与实验分析的比较 分析方法可分为理论计算和实验测试两大类。 1、基于实验的分析方法 指通过的实验测试获取需要的性能参数的方法。这种方法获取不同的性能参数需要采用不同的测试方法、仪器设备和辅助实验装置。如：强度实验，可以采用电阻应变片及应变仪、光弹涂膜或云纹栅、应变涂料等； 扭转与弯曲刚度实验则需要专门的实验台等等。 实验方法的最大优点是工程上普遍认为，实验结果真实可靠，通常被当作产品最终定型的权威性依据。
实验方法的优点： 是工程上普遍认为，实验结果真实可靠， 通常被当作产品最终定型的权威性依据。 实验方法的缺点： 1）实验一定要在样品或样机试制之后才能进行，无法与设计并行。 2）成本高、周期长，通常只适合批量生产的产品。 3）可以获得的数据量有限，无法对设计提供更多的指导，更无法进行结构优化。 4）受实验手段的限制，有些参数无法测准。
2、基于理论计算的分析方法 指通过理论分析或数值计算获取所需的性能参数的分析方法 优点： 1）经济、快捷，成本低、周期短。（与实验相比）2）一次分析可以获得大量的数据。 3）可以与设计同步进行。 4）可以配合优化算法，对设计进行优化。 缺点：分析结果受模型质量、算法理论及边界条件等诸多因素影响。
【MA TLAB 算例】3.3.7(2) 三梁平面框架结构的有限元分析 (Beam2D2Node) 如图3-19所示的框架结构，其顶端受均布力作用，结构中各个 截面的参数都为：113.010Pa E =?，746.510I m -=?，426.810A m -=?。试基 于MA TLAB 平台求解该结构的节点位移以及支反力。 图3-19 框架结构受一均布力作用 解答：对该问题进行有限元分析的过程如下。 （1） 结构的离散化与编号 将该结构离散为3个单元，节点位移及单元编号如图3-20所示， 有关节点和单元的信息见表3-5。 （a ） 节点位移及单元编号
（b）等效在节点上的外力 图3-20 单元划分、节点位移及节点上的外载 （2）各个单元的描述 首先在MA TLAB环境下，输入弹性模量E、横截面积A、惯性矩I和长度L，然后针对单元1，单元2和单元3，分别二次调用函数Beam2D2Node_ElementStiffness，就可以得到单元的刚度矩阵k1(6×6)和k2(6×6)，且单元2和单元3的刚度矩阵相同。 >> E=3E11; >> I=6.5E-7; >> A=6.8E-4; >> L1=1.44; >> L2=0.96; >> k1=Beam2D2Node_Stiffness(E,I,A,L1); >> k2=Beam2D2Node_Stiffness(E,I,A,L2); （3）建立整体刚度方程 将单元2和单元3的刚度矩阵转换成整体坐标下的形式。由于该结构共有4个节点，则总共的自由度数为12，因此，结构总的刚度矩阵为KK(12×12)，对KK清零，然后两次调用函数Beam2D2Node_Assemble进行刚度矩阵的组装。 >> T=[0,1,0,0,0,0;-1,0,0,0,0,0;0,0,1,0,0,0;0,0,0,0,1,0;0,0,0,-1,0,0;0,0,0,0,0,1] ; >> k3=T'*k2*T; >> KK=zeros(12,12); >> KK=Beam2D2Node_Assemble(KK,k1,1,2);
有限元法基本原理与应用 班级机械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)法是将权函数取为逼近函数中的基函数；最小二乘法是令权函数等于余量本身，而内积的极小值则为对代求系数的平方误差最小；在配置法中，先在计
一、命令流 举例： 有一长为 100mm 的矩形截面梁，截面为 10X1mm ，与一规格为 20mmX7mmX10mm 的实体连接， 约束实体的端面， 在梁端施加大小为 3N 的 y 方向的压力， 梁与实体都为一材 料，弹性模量为 30Gpa ，泊松比为 0.3 。本例主要讲解梁与实体连接处如何利用耦合及约束 方程进行处理。 命令流如下： FINI /CLE LSEL,S,LOC,X,21,130 ! 选择梁线 LATT,1,2,2 ! 指定梁的单元属性 LESIZE,ALL,,,10 !指定梁上的单元份数 LMESH,ALL !划分梁单元 VSEL,ALL !选择所有实体 VATT,1,1,1 ! 设置实体的单元属性 ESIZE,1 !指定实体单元尺寸 MSHAPE,0,2D ! 设置实体单元为 2D MSHKEY,1 !设置为映射网格划分方法 VMESH,ALL ! 划分实体单元 ALLS !全选 FINI !退出前处理 /FILNAME,BEAM_AND_SOLID_ELEMENTS_CONNECTION ! 定义工作文件名 /TITLE,COUPLE_AND_CONSTRAINT_EQUATION ! 定义工作名 /PREP7 ET,1,SOLID95 ET,2,BEAM4 MP,EX,1,3E4 MP,PRXY,1,0.3 R,1 R,2,10.0,10/12.0,1000/12.0,10.0,1.0 BLC4,,,20,7,10 WPOFFS,0,3.5 WPROTA,0,90 VSBW,ALL WPOFFS,0,5 WPROTA,0,90 VSBW,ALL WPCSYS,-1 K,100,20,3.5,5 K,101,120,3.5,5 L,100,101 !进入前处理 !定义实体单元类型为 SOLID95 ! 定义梁单元类型为 BEAM4 !定义材料的弹性模量 !定义泊松比 !定义实体单元实常数 !定义梁单元实常数 !创建矩形块为实体模型 !将工作平面向 Y 方向移动 3.5 !将工作平面绕 X 轴旋转 !将实体沿工作平面剖开 !将工作平面向 Y 方向移动 !将工作平面绕 X 轴旋转 !将实体沿工作平面剖开 90 度 5 90 度 !将工作平面设为与总体笛卡儿坐标一致 !创建关键点 !创建关键点 !连接关键点生成梁的线实体
Matlab有限元分析20140226 为了用Matlab进行有限元分析，首先要学会Matlab基本操作，还要学会使用Matlab进行有限元分析的基本操作。 1. 复习：上节课分析了弹簧系统 x 推导了系统刚度矩阵
2. Matlab有限元分析的基本操作 （1）单元划分（选择何种单元，分成多少个单元，标号）（2）构造单元刚度矩阵（列出…） （3）组装系统刚度矩阵（集成整体刚度矩阵） （4）引入边界条件（消除冗余方程） （5）解方程 （6）后处理（扩展计算）
3. Matlab有限元分析实战【实例1】
分析： 步骤一：单元划分
>>k1=SpringElementStiffness(100)
a) 分析SpringAssemble库函数 function y = SpringAssemble(K,k,i,j) % This function assembles the element stiffness % matrix k of the spring with nodes i and j into the % global stiffness matrix K. % function returns the global stiffness matrix K % after the element stiffness matrix k is assembled. K(i,i) = K(i,i) + k(1,1); K(i,j) = K(i,j) + k(1,2); K(j,i) = K(j,i) + k(2,1); K(j,j) = K(j,j) + k(2,2); y = K; b) K是多大矩阵？ 今天的系统刚度矩阵是什么？ 因为 11 22 1212 k k k k k k k k - ?? ?? - ????--+ ?? 所以 1000100 0200200 100200300 - ?? ?? - ????-- ??？
【ANSYS 算例】4.7.1(3) 基于3节点三角形单元的矩形薄板分析 如图4-20所示为一矩形薄平板，在右端部受集中力100 000N F =作用，材料常数为：弹性模量7110Pa E =?、泊松比1/3μ=，板的厚度为0.1m t =，在ANSYS 平台上，按平面应力问题完成相应的力学分析。 (a) 问题描述 (a) 有限元分析模型 图4–20 右端部受集中力作用的平面问题(高深梁) 解答 在ANSYS 平台上，完成的分析如下。 1. 基于图形界面的交互式操作(step by step) (1) 进入ANSYS(设定工作目录和工作文件) 程序 → ANSYS Interactive →Working directory (设置工作目录) →Initial jobname (设置工作文件名): 2D3Node →Run → OK (2) 设置计算类型 ANSYS Main Menu : Preferences… → Structural → OK (3) 选择单元类型 ANSYS Main Menu : Preprocessor →Element Type →Add/Edit/Delete… →Add… →Solid ：Quad 4node 42 →OK (返回到Element Types 窗口) → Options… →K3: Plane Strs w/thk(带厚度的平面应力问题) →OK →Close (4) 定义材料参数 ANSYS Main Menu : Preprocessor →Material Props →Material Models →Structural →Linear →Elastic → Isotropic: EX:1.0e7 (弹性模量)，PRXY: 0.33333333 (泊松比) → OK → 鼠标点击该窗口右上角的“ ”来关闭该窗口 (5) 定义实常数以确定平面问题的厚度 ANSYS Main Menu: Preprocessor →Real Constant s… →Add/Edit/Delete →Add →Type 1→ OK →Real Constant Set No: 1 (第1号实常数), THK: 0.1 (平面问题的厚度) →OK →Close (6) 生成单元模型 生成4个节点 ANSYS Main Menu: Preprocessor →Modeling → Create → Nodes → On Working Plane →输入节点1的x,y,z 坐标(2,1,0)，回车→输入节点2的x,y,z 坐标(2,0,0)，回车→输入节点3的x,y,z 坐标(0,1,0)，回车→输入节点4的x,y,z 坐标(0,0,0)，回车→OK 定义单元属性 ANSYS Main Menu: Preprocessor →Modeling → Create → Elements → Elem Attributes →Element type number:1 →Material number:1→Real constant set number:1 →OK 生成单元 ANSYS Main Menu: Preprocessor →Modeling → Create → Elements → User Numbered → Thru Nodes →Number to assign to element:1→Pick nodes:2,3,4→OK →Number to assign to element:2→Pick nodes:3,2,1→OK (7) 模型施加约束和外载 左边两个节点施加X,Y 方向的位移约束 ANSYS Main Menu: Solution → Define Loads → Apply →Structural → Displacement → On
1、有限元这种数值计算方法起源于20世纪50年代中期航空工程中飞机结构的矩阵分析。 2、有限单元法的基本思想:在力学模型上将一个原来连续的物体离散成为有限个具有一定 大小的单元,这些单元仅在有限个节点上相连接,并在节点上引进等效力以代替实际作用于单元上的外力。 3、节点:网格间相互连接的点。 4、边界:网格与网格的交界线。 5、有限元的优点:①理论基础简明,物理概念清晰,且可在不同的水平上建立起对该法的 理解②具有灵活性和适用性,应用范围极为广泛③该法在具体推导运算中,广泛采用了矩阵方法。 6、有限单元法分类(从选择基本未知量的角度:位移法(以节点位移为基本未知量,通用 性广、力法(以节点力、混合法(一部分以节点位移,另一部分以节点力 7、有限元法分析计算的基本步骤:①结构的离散化②单元分析(选择位移模式,建立单元 刚度方程,计算等效节点力③整体分析④求解方程,得出节点位移⑤由节点位移计算单元的应变与应力。 8、单元划分:将某个机械结构划分为由各种单元组成的计算模型。 9、有限元法基本近似性------几何近似。
10、弹性力学的任务:分析弹性体在受外力作用并处于平衡状态下产生的应力、应变和位移状态及其相互关系等。 11、弹性力学假设所研究的物体是连续的、完全弹性的、均匀的、各向同性的、微小变形的和无初应力的 12、外力:体力(分布在物体体积内的力---重力、惯性力、电磁力面力(分布在物体表面上的力---流体压力、接触力、风力 13、应力:物体受外力作用,或由于温度有所改变,其内部发生的内力。σ={ σx σy σz τx τy τz } = [σx σy σz τx τy τz ]T 14、应变:物体受到外力作用时,其形状发生改变时的形变。---长度和角度。 ε={ εx εy εz γx γy γz } = [εx εy εz γx γy γz ]T 15、位移:弹性体在载荷作用下,不仅会发生形变,还将产生位移,即弹性体位置 的移动。 δ={u v w }=[u v w ]T 16:、变形协调条件:设想在变形前,把弹性体分为许多微小立方单元体。变形后,每个单元体都产生任意变形而变成一些六面体。可能发生这样的情况,这些六面体
1. 诉述有限元法的定义 答：有限元法是近似求解一般连续场问题的数值方法 2. 有限元法的基本思想是什么 答：首先，将表示结构的连续离散为若干个子域，单元之间通过其边界上的节点连接成组合体。其次，用每个单元内所假设的近似函数分片地表示求解域内待求的未知厂变量。 3. 有限元法的分类和基本步骤有哪些 答：分类：位移法、力法、混合法；步骤：结构的离散化，单元分析，单元集成，引入约束条件，求解线性方程组，得出节点位移。 4. 有限元法有哪些优缺点 答：优点：有限元法可以模拟各种几何形状复杂的结构，得出其近似解；通过计算机程序，可以广泛地应用于各种场合；可以从其他CAD软件中导入建好的模型；数学处理比较方便，对复杂形状的结构也能适用；有限元法和优化设计方法相结合，以便发挥各自的优点。 缺点：有限元计算，尤其是复杂问题的分析计算，所耗费的计算时间、内存和磁盘空间等计算资源是相当惊人的。对无限求解域问题没有较好的处理办法。尽管现有的有限元软件多数使用了网络自适应技术，但在具体应用时，采用什么类型的单元、多大的网络密度等都要完全依赖适用者的经验。 5. ?梁单元和平面钢架结构单元的自由度由什么确定 答：每个节点上有几个节点位移分量，就称每个节点有几个自由度 6. ?简述单元刚度矩阵的性质和矩阵元素的物理意义 答：单元刚度矩阵是描述单元节点力和节点位移之间关系的矩阵 单元刚度矩阵中元素aml的物理意义为单元第L个节点位移分量等于1，其他节点位移分量等于0时，对应的第m个节点力分量。 7. 有限元法基本方程中的每一项的意义是什么 答：整个结构的节点载荷列阵（外载荷、约束力），整个结构的节点位移列阵，结构的整体刚度矩阵，又称总刚度矩阵。 8. 位移边界条件和载荷边界条件的意义是什么 答：由于刚度矩阵的线性相关性不能得到解，从而引入边界条件。 9. ?简述整体刚度矩阵的性质和特点 答：对称性；奇异性；稀疏性；对角线上的元素恒为正。 11. 简述整体坐标的概念 答：单元刚度矩阵的坐标变换式把平面刚架的所有单元在局部坐标系X’Y’Z’下的单元刚度矩阵变换到一个统一的坐标系xOy下，这个统一的坐标系xOy称为整体坐标系。 13. 简述平面钢架问题有限元法的基本过程 答：力学模型的确定，结构的离散化，计算载荷的等效节点力，计算各单元的刚度矩阵，组集整体刚度矩阵，施加边界约束条件，求解降价的有限元基本方程，求解单元应力，计算结果的输出。 14. 弹性力学的基本假设是什么。 答：连续性假定，弹性假定，均匀性和各向同性假定，小变形假定，无初应力假定。 15.弹性力学和材料力学相比，其研究方法和对象有什么不同。 答：研究对象：材料力学主要研究杆件，如柱体、梁和轴，在拉压、剪切、弯曲和扭转等作用下的应力、形变和位移。弹性力学研究各种形状的弹性体，除杆件外，还研究平面体、空间体，板和壳等。因此，弹性力学的研究对象要广泛得多。研究方法：弹性力学和材料力学
第一章结构分析及有限元分析基础知识 注：摘自《NX知识工程应用技术——CAD/CAE篇》 洪如瑾编译 清华大学出版社 [目标] 本章将简述结构分析及有限元分析的基础知识，为学习与应用结构分析做好准备，包括： ※ 结构与结构分析定义 ※ 结构的线性静态分析 ※ 材料行为与故障 ※ 有限元分析的基本概念 ※ 有限元模型 1.1结构分析基础知识 1.1.1结构基本概念 1.结构定义 结构可以定义为一个正承受作用的载荷处于平衡中的系统。平衡条件意味着结构是不移动的。一个自由的支架不是一个结构，它未被连接到任一物体上并无载荷作用与它。仅当它附着到外部世界，并且有作用力、压力或力矩时，支架成为一个结构。 例如横跨江面的大桥就是一个普通的结构，一个支架通过它的支撑连接到地面上，桥的重量是在结构上的一种载荷（力）。当汽车通过桥时，附加的力作用于桥的不同位置。 一个好的结构必须满足以下标准： （1） 当预期的载荷作用时，结构必须不出现故障。这个似乎是显而易见的，并意味着结构必须是“强度足够的”。故障意味着结构破裂、分离、弯曲，以及支撑作用载荷失败。 注意：考虑到意外的载荷，通常在设计中提供安全余量。余量常常利用安全因素来描述。例如，如果在结构上期待载荷是10 000磅，规定安全因素是2.0，则结构将设计成能经受住20 000磅载荷。 （2） 当载荷作用时，结构必须不产生过分变形。这意味着结构必须“刚度足够”。 变形可接受的极限（弯曲度、挠度、拉伸等）取决于特定情况。例如，在通常住宅中的地板由足够的吊带支撑，以防止当人在地板岸上行走时有“柔软”的感觉。 （3） 在它的服务生命周期，结构的行为应不会恶化。这意味着结构必须“足够耐用”，必须考虑环境影响和“磨损与破裂”。如果一座桥假定维持50年，则桥的设计必须提供整个50年寿命的结构完整性与充分的安全余量。2.结构分析 结构分析是用于决定一个结构是否将正确完成任务的工程分析过程。结构将在某些方式中进行模拟和求解描述它的行为的数学方程。分析可以人工方法或用计算机方法来完成。 结构分析的结果（答案）用于评估性能，摘要如下： （1）“强度足够吗？”：应力必须是在一可接受的范围内。 （2）“刚度足够吗？”：位移必须是在一可接受的范围内。 （3）“耐用度足够？”：对一个长的疲劳周期应力必须足够低。
有限元分析基础教程Fundamentals of Finite Element Analysis (ANSYS算例) 曾攀 清华大学 2008-12
有限元分析基础教程曾攀 有限元分析基础教程 Fundamentals of Finite Element Analysis 曾攀 (清华大学) 内容简介 全教程包括两大部分，共分9章；第一部分为有限元分析基本原理，包括第1章至第5章，内容有：绪论、有限元分析过程的概要、杆梁结构分析的有限元方法、连续体结构分析的有限元方法、有限元分析中的若干问题讨论；第二部分为有限元分析的典型应用领域，包括第6章至第9章，内容有：静力结构的有限元分析、结构振动的有限元分析、传热过程的有限元分析、弹塑性材料的有限元分析。本书以基本变量、基本方程、求解原理、单元构建、典型例题、MATLAB程序及算例、ANSYS算例等一系列规范性方式来描述有限元分析的力学原理、程序编制以及实例应用；给出的典型实例都详细提供有完整的数学推演过程以及ANSYS实现过程。本教程的基本理论阐述简明扼要，重点突出，实例丰富，教程中的二部分内容相互衔接，也可独立使用，适合于具有大学高年级学生程度的人员作为培训教材，也适合于不同程度的读者进行自学；对于希望在MATLAB程序以及ANSYS平台进行建模分析的读者，本教程更值得参考。 本基础教程的读者对象：机械、力学、土木、水利、航空航天等专业的工程技术人员、科研工作者。
目录 [[[[[[\\\\\\ 【ANSYS算例】3.3.7(3) 三梁平面框架结构的有限元分析 1 【ANSYS算例】4.3.2(4) 三角形单元与矩形单元的精细网格的计算比较 3 【ANSYS算例】5.3(8) 平面问题斜支座的处理 6 【ANSYS算例】6.2(2) 受均匀载荷方形板的有限元分析9 【ANSYS算例】6.4.2(1) 8万吨模锻液压机主牌坊的分析(GUI) 15 【ANSYS算例】6.4.2(2) 8万吨模锻液压机主牌坊的参数化建模与分析(命令流) 17 【ANSYS算例】7.2(1) 汽车悬挂系统的振动模态分析(GUI) 20 【ANSYS算例】7.2(2) 汽车悬挂系统的振动模态分析(命令流) 23 【ANSYS算例】7.3(1) 带有张拉的绳索的振动模态分析(GUI) 24 【ANSYS算例】7.3(2) 带有张拉的绳索的振动模态分析(命令流) 27 【ANSYS算例】7.4(1) 机翼模型的振动模态分析(GUI) 28 【ANSYS算例】7.4(2) 机翼模型的振动模态分析(命令流) 30 【ANSYS算例】8.2(1) 2D矩形板的稳态热对流的自适应分析(GUI) 31 【ANSYS算例】8.2(2) 2D矩形板的稳态热对流的自适应分析(命令流) 33 【ANSYS算例】8.3(1) 金属材料凝固过程的瞬态传热分析(GUI) 34 【ANSYS算例】8.3(2) 金属材料凝固过程的瞬态传热分析(命令流) 38 【ANSYS算例】8.4(1) 升温条件下杆件支撑结构的热应力分析(GUI) 39 【ANSYS算例】8.4(2) 升温条件下杆件支撑结构的热应力分析(命令流) 42 【ANSYS算例】9.2(2) 三杆结构塑性卸载后的残余应力计算(命令流) 45 【ANSYS算例】9.3(1) 悬臂梁在循环加载作用下的弹塑性计算(GUI) 46 【ANSYS算例】9.3(2) 悬臂梁在循环加载作用下的弹塑性计算(命令流) 49 附录 B ANSYS软件的基本操作52 B.1 基于图形界面(GUI)的交互式操作(step by step) 53 B.2 log命令流文件的调入操作(可由GUI环境下生成log文件) 56 B.3 完全的直接命令输入方式操作56 B.4 APDL参数化编程的初步操作57
有限元法基础试题（A ） 一、填空题（5×2分） 1.1单元刚度矩阵e T k B DBd Ω = Ω? 中，矩阵B 为__________，矩阵D 为___________。 1.2边界条件通常有两类。通常发生在位置完全固定不能转动的情况为_______边界，具体指定有限的非零值位移的情况，如支撑的下沉，称为_______边界。 1.3内部微元体上外力总虚功： ()(),,,,e x x xy y bx xy x y y by d W F u F v dxdy δστδτσδ??=+++++??+(),,,,x x y y xy y x u v u u dxdy σδσδτδδ??+++??的表达式中，第一项为____________________的虚功，第二项为____________________的虚功。 1.4弹簧单元的位移函数1N +2N =_________。 1.5 ij k 数学表达式：令j d =_____，k d =_____，k j ≠，则力i ij F k =。 二、判断题（5×2分） 2.1位移函数的假设合理与否将直接影响到有限元分析的计算精度、效率和可靠性。（ ） 2.2变形体虚功原理适用于一切结构（一维杆系、二维板、三位块体）、适用于任何力学行为的材料（线性和非线性），是变形体力学的普遍原理。 （ ） 2.3变形体虚功原理要求力系平衡，要求虚位移协调，是在“平衡、协调”前提下功的恒等关系。 （ ） 2.4常应变三角单元中变形矩阵是x 或y 的函数。 （ ） 2.5 对称单元中变形矩阵是x 或y 的函数。 （ ） 三、简答题（26分） 3.1列举有限元法的优点。（8分） 3.2写出有限单元法的分析过程。（8分） 3.3列出3种普通的有限元单元类型。（6分） 3.4简要阐述变形体虚位移原理。（4分） 四、计算题（54分） 4.1对于下图所示的弹簧组合，单元①的弹簧常数为10000N/m ，单元②的弹簧常数为20000N/m ，单元③的弹簧常数为10000N/m ，确定各节点位移、反力以及单元②的单元力。（10分） 4.2对于如图所示的杆组装，弹性模量E 为10GPa ，杆单元长L 均为2m ，横截面面积A 均为2×10-4m 2，弹簧常数为2000kN/m ，所受荷载如图。采用直接刚度法确定节点位移、作用力和单元②的应力。（10分）
有限元法理论及应用大作业 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.1 有限元方法的发展历史、现状和将来 一，历史 有限元法的起源应追溯到上世纪40年代（20世纪40年代）。1943年R.Courant从数学的角度提出了有限元法的基本观点。50年代中期在对飞机结构的分析中，诞生了结构分析的矩阵方法。1960年R.W.Clough在分析弹性力学平面问题时引入了“Finite Element Method”这一术语，从而标志着有限元法的思想在力学分析中的广泛推广。 60、70年代计算机技术的发展，极大地促进了有限元法的发展。具体表现在： 1）由弹性力学的平面问题扩展到空间、板壳问题。 2）由静力平衡问题——稳定性和动力学分析问题。 3）由弹性问题——弹塑性、粘弹性等问题。 二，现状 现在有限元分析法的应用领域已经由开始时的固体力学，扩展到流体力学、传热学和电磁力学等多个传统的领域。已经形成了一种非常成熟的数值分析计算方法。大型的商业化有限元分析软件也是层出不穷，如： SAP系列的代表SAP2000（Structure Analysis Program） 美国安世软件公司的ANSYS大型综合有限元分析软件 美国航天航空局的NASTRAN系列软件 除此以外，还有MASTER、ALGO、ABIQUES、ADINA、COSMOS等。 三，将来 有限元的发展方向最终将和CAD的发展相结合。运用“四个化”可以概括其今后的发展趋势。那就是：可视化、集成化、自动化和网络化。 §1.2 有限元法的特点 机械零构件的受力分析方法总体说来分为解析法和数值法两大类。如大家学过的材料力学、结构力学等就是经典的解析力学分析方法。在这些解析力学方法中，弹性力学的分析方法在数学理论上是最为严谨的一种分析方法。 其解题思路是：从静力、几何和物理三个方面综合考虑，建立描述弹性体的平衡、应力、应变和位移三者之间的微分方程，然后考虑边界条件，从而求出微分方程的解析解。其最大的有点就是，严密精确。缺点就是微分方程的求解困难，很多情况下，无法求解。 数值方法是一种近似的计算方法。具体又分为“有限差分法”和“有限元法”。 “有限差分法”是将得到的微分方程离散成近似的差分方程。通过对一系列离散的差分
有限元分析及其应用-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、问题的微分方程提法、等效积分提法和泛函变分提法之间有何关系？何谓“强形式”？何谓“弱形式”，两者有何区别？建立弱形式的关键步骤是什么？答：强弱的区分在于是否完全满足物理模型的条件。所谓强形式，是指由于物理模型的复杂性，各种边界条件的限制，使得对于所提出的微分方程，对所需要求得的解的要求太强。也