下载文档原格式
材料力学有限元法知识点总结材料力学是一门研究物质内部结构、性质和变形行为的学科，而有限元法则是一种在工程和科学领域中广泛应用的数值计算方法。
有限元法可以将一个复杂的实体划分为无数小的单元，通过对这些小单元进行分析和计算，最终得到整个实体的力学性质和行为。
本文将对材料力学有限元法的一些核心概念和知识点进行总结。
1. 有限元法基础概念有限元法基于将实际连续的物体离散为有限数量的单元，通过计算每个单元的受力、变形等性质，再通过组合这些单元的结果来近似整个物体的行为。
它包含以下几个基础概念：1.1 单元（Element）：有限元法中的基本组成单元，可以是一维的线段、二维的三角形或四边形，或三维的四面体、六面体等。
1.2 节点（Node）：单元的角点或边上的点，用于定义单元之间的连接关系和边界条件。
1.3 自由度（Degree of Freedom）：每个节点与力学性质相关的物理量，如位移、应力等。
根据问题的不同，在每个节点上可能有一个或多个自由度。
1.4 单元刚度矩阵（Element Stiffness Matrix）：描述单元内部受力和变形关系的矩阵，在有限元法中通过组合所有单元的刚度矩阵来得到整个系统的刚度矩阵。
1.5 全局刚度矩阵（Global Stiffness Matrix）：由所有单元刚度矩阵组合而成的整个系统的刚度矩阵，用于计算节点的位移和应力。
2. 有限元法的数学原理有限元法的数学原理主要基于以下两个方面：2.1 变分原理（Variational Principle）：有限元法的数学基础是根据变分原理推导实现的。
它通过对结构的势能进行变分并进行最小化，得到满足结构力学行为和边界条件的位移和应力场。
2.2 加权残差法（Weighted Residuals Method）：有限元法通过将变分原理中的势能函数展开为一系列基函数的线性组合，并使用权重函数对残差进行加权求和的方式进行近似。
这样可以将求解连续问题转化为离散问题，进而进行数值计算。
《有限元基础及应用》课程教学大纲一、课程基本信息二、课程目标（一）总体目标：有限元法是求解复杂工程问题进行数值模拟非常有效的方法，是现代数字化科技的一种重要基础性原理。
将它应用于科学研究中，可以成为探究物质客观规律的先进手段；将它应用于工程技术中，可成为工程设计和分析的可靠工具。
有限元法已经成为机械工程、车辆工程、航空航天工程、土木建筑等专业的必修课或选修课，有限元商用软件也是广大工程技术人员从事产品开发、设计、分析，以及生产服务的重要工具。
通过本课程的学习使同学们掌握有限元分析方法的基础知识和原理；掌握大型有限元分析软件（ANSYS）的使用；有限元方法的实际应用：能够针对具有复杂几何形状的变形体完整获取复杂外力作用下它内部准确力学信息，在准确进行力学分析的基础上，可以对所研究对象进行强度、刚度等方面的判断，以便对研究结构进行静态、动态的强度和刚度分析、参数设计以及结构优化设计。
内容由浅入深，通俗易懂，结合实践应用分析，培养学生理论联系实际和解决实际问题的能力。
（二）课程目标：课程目标1：掌握有限元方法的基本原理，分析过程和步骤，形函数的构造方法，以及针对不同维度、不同结构准确选择合适的单元的技巧；课程目标2：掌握有限元分析方法，具有对不同工程问题建立相应力学模型再选取适合的有限元模型离散，最后得到高精度低成本的数值模拟结果；课程目标3：利用有限元原理和应用软件（ANSYS），能够针对车辆结构中具有复杂几何形状的零部件完整获取复杂外力作用下其内部的准确力学信息（位移、应力和应变），并能根据强度、刚度、稳定性及疲劳等进行分析判断结构的安全性，具有分析和解决工程实际问题的能力；课程目标4：掌握大型商用有限元软件（ANSYS）对车辆结构部件的静力学、动力学和多物理场耦合问题进行数值模拟和分析。
能够了解不同单元的适用范围以及有限元方法数值模拟的局限性。
（三）课程目标与毕业要求、课程内容的对应关系本课程支撑专业培养计划中毕业要求1、2、3、5。
有限元理论基础有限元方法的基础是变分原理和加权余量法，其基本求解思想是把计算域划分为有限个互不重叠的单元，在每个单元内，选择一些合适的节点作为求解函数的插值点，将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式，借助于变分原理或加权余量法，将微分方程离散求解。
釆用不同的权函数和插值函数形式，便构成不同的有限元方法。
4.加权余量法：是指采用使余量的加权函数为零求得微分方程近似解的方法称为加权余量法。
(Weighted residual method WRM)是一种直接从所需求解的微分方程及边界条件出发，寻求边值问题近似解的数学方法。
加权余量法是求解微分方程近似解的一种有效的方法。
设问题的控制微分方程为：在V域内厶(")-八0 (5.1.1)在S 边界上〃(“)-& = 0 (5.1.2)式中：L、B——分别为微分方程和边界条件中的微分算子；f、g ——为与未知函数u无关的己知函数域值；u——为问题待求的未知函数当弄!J用力u权余•肚法求近丁以解首先在求耳军域上理立一个T式閑数H 一般兵升如下形式:仁土CN=NC(5.1.3)T M式中：c{----------- 彳寺定系数. 也可称为广义坐标；N：--- 取白完备函冬攵*S线.性无关的基函孕攵°由于〃一般只圮彳守求函缨攵U的近1以耳岂因u匕将式(5 1.3) 代入式(5 1 1)牙口式(5 1.2)后将诃•不誉斯兄，昔迅:| R] = L(flb— f在V域内\R B =B(^~g在S 边界上("14)城然 & 、尽反映了r式函竽攵与实解之问的偏差. 它丁门分另U称做内召卩牙口边界余覺。
若在域\'内引入内部权函数硏，在边界S上引入边界权函数W B 则可理立11个消除余甘的条件.一般可农示为：L兀W B1R B dS = 0 (/ = L2.L ,〃) (51-5)• V • S不同的权函数幵；和jr R反映了不同的消除余•眩的准则。
有限单元法的数学基础1、引言有限元方法归根结底是一种数值计算方法，它有严格的数学证明作为其近似的客观性和合理性的保证。
力学问题最终归结为一组微分方程的边值问题或者初值问题抑或是混合问题。
比如弹性静力学最终归结为L-N 方程的微分提法。
在很难或者根本不可能得到所得方程的理论解的情况下，究竟用什么样的方法才能得到方程的近似解（这种近似解已经能够满足实际工程的需要），在这种情况下，二十世纪五六十年代由结构力学家进而由数学家提出和证明了这种思想方法的合理性。
有限元方法产生于力学计算，但是，它本质上并不是力学的专利。
世间万物的变化过程很多都可以通过微分方程特别是偏微分方程来描述，也就是说，微分方程是很多现象和过程的数学结构，而大多数的微分方程是不能得到理论解的，这时候就可以使用有限元方法来求其近似解，因为有限元方法是求解微分方程（组）的数值计算方法。
它适用于力学的微分方程，也同样适用于其它领域的相应的微分方程的数值求解。
2、有限元方法数学根源对于一个给定的微分方程定解问题，为了求其近似解，我们可以使用Ritz 方法和Galerkin 方法。
下面分别阐述这两种方法，然后讨论有限元方法和他们的关系。
（1） Ritz 法Ritz 法源于最小势能原理，设H 是可分的Hilbert 空间，在H 中取有限维空间Sn ，它是由N 个线性无关向量12,,,N φφφ 张成，即：121,,(,,)NN n n i i N N i S C C C C R ωωφ=⎧⎫≡=∀∈⎨⎬⎩⎭∑用N S 代替H ，在N S 上求泛函J(w)的极值,即求N U ∈N S ，使得()N J U =min ()N N S N J ωω∈实际上寻求N U 只需通过解一个线性方程组1()(,)()02J D F ωωωω=-≥D--------双线性形式 F--------线性泛函1NN i i i C ωφ==∑111,111()(,)()21(,)()2N N NN i i i i i i i i i NN i j i j i ii j i J D C C F C D C C F C ωφφφφφφ====== =-∑∑∑∑∑－因此，()N J ω是一个以12,,,N C C C 为未知数（自变量）的二次多项式12(,,,)N j C C C ，如果二次项的系数矩阵,1,2,,[(,)]i j i j N D φφ= 是正定的，那么12(,,,)N j j C C C = 在N+1维空间是一个开口向上的椭球抛物面，它有且只有一个极（最）小值点，所谓在N S 上求()N J ω的极值，就是确定00012,,,N C C C ，使得：00012(,,,)N j C C C ＝1000,,12min (,,,)N C C R N j C C C ∈极值条件：ijC ∂∂|00012,,,N C C C ＝0 （1,,i N = ） 得:01()()ni ji i i D CF φφφ==∑ (1,,i N = )即：00012[,,,]T N C C C C = 适合方程组：KC=F11[(),,()]T F F F φφ=112111222212(,)(,)(,)(,)(,)(,)(,),(,),,(,)N N N N N N D D D D D D K D D D φφφφφφφφφφφφφφφφφφ⎛⎫ ⎪ ⎪= ⎪ ⎪⎝⎭，，，，，， 。
§5-5数值积分１、问题的提出在上一节中对等参元进行单元分析时要进行下列积分： (i) 单元刚度矩阵(ii)体积力的等效结点力(iii)边界力的等效结点力(iv)温升载荷的等效结点力式（5-4-5）~（5-4-8）分别归结为计算以下两种形式的积分对于上述积分仅在单元的形状十分规则的情况下才能得到解析的结果（精确值），一般情况只能用数值积分方法（主要是高斯求积法）求近似值。
虽然数值积分是“被迫“采用的，但后来发现：有选择地控制积分点的个数和位置，可以方便地实现我们的某些特殊意图。
这样一来，数值积分就成为有限元分析的一个重要组成部分，以至本来可以精确积分的三角形单元也常常采用数值积分。
２、数值积分的基本概念任何积分工作取决于三个要素：给定的积分区间，给定的被积函数，具体的积分方法。
下面以一维情况为例介绍数值积分的基本概念 (i) 梯形法函数()x f 在区间(a,b)的积分可以表达为 ()()ini ibax f W dx x f I ∑⎰=≈=1⎰⎰⎰---111111),()(dxdxy x f dx x f 、 [][][][][][][]ηξd d J t B E B tdxdyB E B k T Te det 1111⎰⎰⎰⎰--=={}[][]ηξσd d J t f f N td f f N r y x T y x T eV det 1111⎭⎬⎫⎩⎨⎧=⎭⎬⎫⎩⎨⎧=⎰⎰⎰⎰--{}[]{}ηξσγd Jd t B T det 01111T ⎰⎰--={}[]()()⎪⎪⎭⎫ ⎝⎛=⎭⎬⎫⎩⎨⎧=⎰⎰⎰--dy y f dx x f tds q p N r T 1111,ΓΓ(5-4-5) (5-4-8) (5-4-7) (5-4-6)i W ：权系数； i x ：积分样点；()i x f ：积分样点的函数值。
梯形法的求积公式为其中，1--=n ab h ，而a b W ni i -=∑=1(ii) 当被积函数为n-1次多项式P n-1(x )时，则由n 个样点及其样点值（x i , P n-1(x i ),i=1,n ）可以精确重构这个多项式，从而可以得到精确解。
有限元法有限元法是一套求解微分方程的系统化数值计算方法，它比传统解法具有理论完整可靠，物理意义直观明确，解题效能强等优点，特别是由于这种方法适应性强，形式单纯、规范，所以近年来在电子计算机的配合下，已推广应用到很多工程技术部门和某些科学领域。
本章是从应用的角度来介绍有限元法的基本知识，首先通过典型的位移法阐述有限元法的一般原理与解算过程，然后叙述了剖分单元的技巧，最后介绍与有限元法有关的弹性力学问题。
常用符号规定如下(括号内为力学术语或释例)：Ω，表示区域及其边界。
表示区域Ω的单元及其边界。
表示单元的第i个顶点，简记作节点i。
表示系数(刚度)矩阵。
()表示单元的系数(刚度)矩阵。
(x，y，z)表示总体的直角坐标。
(）表示单元的局部坐标。
(，，），（，，，）等表示单元的自然坐标。
(x，y ，)表示节点i的直角坐标。
(u，v，w）表示一组待定函数（分别为沿x，y，z方向的位移分量），其列矢量表示为u。
1(u，v，w)表示（u，v，w）在单元上的插值函数，其列矢量表示为u。
(u，v，w）表示节点i的函数（位移）值。
{u，v，w}表示节点i的一组参数值，即函数直到某阶导数在节点i上的值按一定次序排成的列矢量{u}。
例如{u}= {u，v，w}=（u，u，u，u，v，v，v，v，w，w，w，w）式中τ表示转置。
{u，v，w}表示{u，v，w}按单元的节点序号排成的列矢量，表示为{u}。
等表示单元的型函数。
{R}表示n次多项式中含变量x，y，z各项按一定次序排成的列矢量，并以表示其中第k个分量。
例如二元二次多项式{}表示n 次多项式中各项相应的系数排成的列矢量，并以表示其中第k个分量。
例如对于{}，{}={f，g，h}表示与节点参数相对应的一组已知函数及其导数在节点i上的值排成的列矢量。
2{f，g，h}表示{f，g，h}按单元的节点序号排成的列矢量。
，或放在定义或公式之后表示其中函数u，v，w，变量x，y，z或下标i，j，k作循环替换后，该定义或公式仍然成立。
有限元方法求解微分方程，特别是椭圆型边值问题的一种离散化方法，其基础是变分原理和剖分逼近。
有限元方法是传统的里茨－加廖金方法的发展，并融会了差分法的优点，处理上统一，适应能力强，已广泛应用于科学与工程中庞大复杂的计算问题。
作为有限元方法出发点的变分原理，是表达物理基本定律的一种普遍形式。
其表述可概括如下：给出一个依赖物理状态v的变量J(v)（v是函数,J(v)在数学上称为泛函），同时给出J(v)的容许函数集V,即一切可能的物理状态,则真实的状态是V中使J(v)达到极小值的函数。
剖分逼近是有限元离散化的手段，把问题的整体（即求解域）剖分为有限个基本块，称为"单元"，然后通过单元上的插值逼近，得到一个结构简单的函数集，称为"有限元空间",它一般是容许函数集V的子集或有某种联系。
有限元方法就是在这个有限元空间中寻找J(v)的极小解作为近似解。
典型问题为具体说明有限元方法，讨论二维有界域Ω上的椭圆型方程, (1)变系数β表示介质不均匀。
物理学中许多平衡态或定常态问题都可归结为这个典型方程。
与方程(1)相配的有如下三类边界条件：第一类：；第二类：；第三类：。
这里的φ、g及α均为定义在边界дΩ上的已知函数，表示外法向导数，第二类边界条件是第三类当α=0时的特例。
为说明有限元方法能统一处理复杂的情况，假定讨论的问题是混合边值，并且介质有间断,即дΩ分成Г0和Г1两部分，分别有边界条件, (2),(3)β（x,y)有间断线，把Ω分为Ω-,Ω+两部分，在间断线上微分方程(1)无定义，而代之以接触条件, (4)及表示间断线上分别指向Ω+及Ω-的法向导数。
变分原理与微分方程(1)及附加条件(2)、(3)、(4)的边值问题相对应的是物理学中的极小能量原理。
构造"能量积分"并取J(v)的容许函数集V为一切满足边界条件(2)且一阶偏导数平方可积的函数,则使J(v)达到极小值的u,即,(6)也必满足方程(1)及(2)、(3)、(4)。
有限元基础知识归纳有限元知识点归纳1.、有限元解的特点、原因？答：有限元解一般偏小，即位移解下限性原因：单元原是连续体的一部分，具有无限多个自由度。
在假定了单元的位移函数后，自由度限制为只有以节点位移表示的有限自由度，即位移函数对单元的变形进行了约束和限制，使单元的刚度较实际连续体加强了，因此，连续体的整体刚度随之增加，离散后的刚度较实际的刚度K为大，因此求得的位移近似解总体上将小于精确解。
2、形函数收敛准则（写出某种单元的形函数，并讨论收敛性）P49(1)在节点i处Ni=1,其它节点Ni=0;(2)在单元之间，必须使由其定义的未知量连续；(3)应包含完全一次多项式；(4)应满足∑Ni=1以上条件是使单元满足收敛条件所必须得。
可以推证，由满足以上条件的形函数所建单元是完备协调的单元，所以一定是收敛的。
4、等参元的概念、特点、用时注意什么？（王勖成P131）答：等参元—为了将局部坐标中几何形状规则的单元转换成总体（笛卡尔）坐标中的几何形状扭曲的单元，以满足对一般形状求解域进行离散化的需要，必须建立一个坐标变换。
即：为建立上述的变换，最方便的方法是将上式表示成插值函数的形式，即：其中m是用以进行坐标变换的单元节点数，xi,yi,zi是这些结点在总体（笛卡尔）坐标内的坐标值，Ni’称为形状函数，实际上它也是局部坐标表示的插值函数。
称前者为母单元，后者为子单元。
还可以看到坐标变换关系式和函数插值表示式：在形式上是相同的。
如果坐标变换和函数插值采用相同的结点，并且采用相同的插值函数，即m=n，Ni’=Ni,则称这种变换为等参变换。
5、单元离散？P42答：离散化既是将连续体用假想的线或面分割成有限个部分，各部分之间用有限个点相连。
每个部分称为一个单元，连接点称为结点。
对于平面问题，最简单、最常用的离散方式是将其分解成有限个三角形单元，单元之间在三角形顶点上相连。
这种单元称为常应变三角形单元。
常用的单元离散有三节点三角形单元、六节点三角形单元、四节点四边形单元、八节点四边形单元以及等参元。
（3001）《固体力学》专业综合考试内容:一、结构力学1. 静定与静不定杆系结构内力求解-力法；2. 杆系结构的位移解法（直接刚度法）；3.薄壁工程梁理论；4. 结构力学中的能量原理及其应用。
包括应变能及余应变能、虚功原理、最小位能原理、最小余能原理、卡氏定理、单位载荷定理、叠加原理与位移互等定理。
二、弹性力学1. 平面问题的基本理论。
包括平面问题的平衡微分方程、几何方程、物理方程、及应力函数、逆解法与半逆解法等；2. 平面问题的直角坐标解答。
包括多项式解答、矩形梁的纯弯曲、简支梁受均布荷载等；3. 平面问题的极坐标解答。
包括极坐标中的平衡微分方程、几何方程、物理方程、轴对称应力和相应的位移等；4. 空间问题的基本理论。
包括空间问题的平衡微分方程、主应力与应力主向、几何方程、物理方程、轴对称问题及球对称问题的基本方程等；5. 弹性力学中的变分解法。
弹性体的位移变分方程及最小势能原理，基于最小势能原理的近似解法，应力变分方程及最小余能原理，基于最小余能原理的近似解法。
三、有限元方法及应用1．有限元法中的变分原理。
包括最小位能原理、最小余能原理、广义变分原理、修正变分原理等。
2. 有限元方法的一般性原理和表达格式； 3. 单元和插值函数的构造；4. 等参元和数值积分。
四、结构振动理论1．单自由度系统自由振动。
包括能量法，无阻尼自由振动，有阻尼衰减振动。
2.单自由度系统强迫振动。
包括简谐激励下的响应，强迫振动的复指数解法，频响函数与频域分析。
3.多自由度系统的振动。
包括拉格朗日方程，实模态分析，复模态分析，假设模态法。
参考书目：1. 结构力学基础，黄其青，王生楠编，西安：西北工业大学出版社2. 飞行器结构力学，王生楠主编，西安：西北工业大学出版社3. 弹性力学，吴家龙编著，北京：高等教育出版社4. 有限元法中的变分原理，王生楠编，西安：西北工业大学出版社5. 有限单元法王勖成编著，北京：清华大学出版社，20036. 振动理论及应用方同，薛璞著，西安：西北工业大学出版社。
有限单元法知识点总结1. 有限元法概述有限单元法（Finite Element Method ，简称FEM）是一种数值分析方法，适用于求解工程结构、热传导、流体力学等领域中的强耦合、非线性、三维等问题，是一种求解偏微分方程的数值方法。
有限元法将连续的物理问题抽象为由有限数量的简单几何单元（例如三角形、四边形、四面体、六面体等）组成的离散模型，通过对单元进行适当的数学处理，得到整体问题的近似解。
有限元法广泛应用于工程、材料、地球科学等领域。
2. 有限元法基本原理有限元法的基本原理包括离散化、加权残差法和形函数法。
离散化是将连续问题离散化为由有限数量的简单单元组成的问题，建立有限元模型。
加权残差法是选取适当的残差形式，并通过对残差进行加权平均，得到弱形式。
形函数法是利用一组适当的形函数来表示单元内部的位移场，通过形函数的线性组合来逼近整体位移场。
3. 有限元法的步骤有限元法的求解步骤包括建立有限元模型、建立刚度矩阵和载荷向量、施加边界条件、求解代数方程组和后处理结果。
建立有限元模型是将连续问题离散化为由简单单元组成的问题，并确定单元的连接关系。
建立刚度矩阵和载荷向量是通过单元的应变能量和内力作用，得到整体刚度矩阵和载荷向量。
施加边界条件是通过给定位移或力的边界条件，限制未知自由度的取值范围。
求解代数方程组是将有限元模型的刚度方程和载荷方程组成一个大型代数方程组，通过数值方法求解。
后处理结果是对数值结果进行处理和分析，得到工程应用的有用信息。
4. 有限元法的元素类型有限元法的元素类型包括结构单元、板壳单元、梁单元、壳单元、体单元等。
结构单元包括一维梁单元、二维三角形、四边形单元、三维四面体、六面体单元。
板壳单元包括各种压力单元、弹性单元、混合单元等。
梁单元包括梁单元、横梁单元、大变形梁单元等。
壳单元包括薄壳单元、厚壳单元、折叠单元等。
体单元包括六面体单元、锥体单元、八面体单元等。
5. 有限元法的数学基础有限元法的数学基础包括变分法、能量方法、有限元插值等。
有限元方法基础教程第三版答案第二单元1. 诉述有限元法的定义答:有限元法是近似求解一般连续场问题的数值方法2. 有限元法的基本思想是什么答:首先,将表示结构的连续离散为若干个子域,单元之间通过其边界上的节点连接成组合体。
其次,用每个单元内所假设的近似函数分片地表示求解域内待求的未知厂变量。
3. 有限元法的分类和基本步骤有哪些答:分类:位移法、力法、混合法;步骤:结构的离散化,单元分析,单元集成,引入约束条件,求解线性方程组,得出节点位移。
4. 有限元法有哪些优缺点答:优点:有限元法可以模拟各种几何形状复杂的结构,得出其近似解;通过计算机程序,可以广泛地应用于各种场合;可以从其他CAD软件中导入建好的模型;数学处理比较方便,对复杂形状的结构也能适用;有限元法和优化设计方法相结合,以便发挥各自的优点。
缺点:有限元计算,尤其是复杂问题的分析计算,所耗费的计算时间、内存和磁盘空间等计算资源是相当惊人的。
对无限求解域问题没有较好的处理办法。
尽管现有的有限元软件多数使用了网络自适应技术,但在具体应用时,采用什么类型的单元、多大的网络密度等都要完全依赖适用者的经验。
5. 梁单元和平面钢架结构单元的自由度由什么确定答:由每个节点位移分量的总和确定6. 简述单元刚度矩阵的性质和矩阵元素的物理意义答:单元刚度矩阵是描述单元节点力和节点位移之间关系的矩阵单元刚度矩阵中元素aml的物理意义为单元第L个节点位移分量等于1,其他节点位移分量等于0时,对应的第m个节点力分量。
7. 有限元法基本方程中的每一项的意义是什么P14 答:Q——整个结构的节点载荷列阵(外载荷、约束力);整个结构的节点位移列阵;结构的整体刚度矩阵,又称总刚度矩阵。
8. 位移边界条件和载荷边界条件的意义是什么答:由于刚度矩阵的线性相关性不能得到解,引入边界条件,使整体刚度矩阵求的唯一解。
9. 简述整体刚度矩阵的性质和特点P14 答:对称性;奇异性;稀疏性;对角线上的元素恒为正。
