第2章有限元分析过程概要精选PPT.ppt

《第2章有限元分析过程概要精选PPT.ppt》由会员分享，可在线阅读，更多相关《第2章有限元分析过程概要精选PPT.ppt（60页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第第2章有限元分析过程概要章有限元分析过程概要1第1页，此课件共60页哦有限元分析过程的概要有限元分析过程的概要 本章先通过一个简单的实例，采用直接的推导方法，逐步展示本章先通过一个简单的实例，采用直接的推导方法，逐步展示有限元分析的基本流程，从中可以了解有限元方法的思路形成有限元分析的基本流程，从中可以了解有限元方法的思路形成过程，以及如何由具体的求解步骤归纳出一种通用的标准求解过程，以及如何由具体的求解步骤归纳出一种通用的标准求解方法。方法。第二章第二章2第2页，此课件共60页哦2.1 有限元分析的目的和概念有限元分析的目的和概念 o任何具有一定使用功能的构件任何具有一定使用功能的构件(称

2、为变形体称为变形体(deformed body)都是由满足要求的材料所制造的，在设计阶段，就需都是由满足要求的材料所制造的，在设计阶段，就需要对该构件在可能的外力作用下的内部状态进行分析，以便核要对该构件在可能的外力作用下的内部状态进行分析，以便核对所使用材料是否安全可靠，以避免造成重大安全事故。对所使用材料是否安全可靠，以避免造成重大安全事故。o描述可承力构件的力学信息一般有三类：描述可承力构件的力学信息一般有三类：(1)构件中因承载在任意位置上所引起的移动构件中因承载在任意位置上所引起的移动(称为称为位移位移(displacement)；(2)构件中因承载在任意位置上所引起的变形状态构件中

3、因承载在任意位置上所引起的变形状态(称称为为应变应变(strain)；(3)构件中因承载在任意位置上所引起的受力状态构件中因承载在任意位置上所引起的受力状态(称为称为应应力力(stress)；3第3页，此课件共60页哦2.1 有限元分析的目的和概念有限元分析的目的和概念o若该构件为简单形状，且外力分布也比较单一，如：杆、梁、若该构件为简单形状，且外力分布也比较单一，如：杆、梁、柱、板就可以采用材料力学的方法，一般都可以给出解析公柱、板就可以采用材料力学的方法，一般都可以给出解析公式，应用比较方便；式，应用比较方便；o但对于几何形状较为复杂的构件却很难得到准确的结果，甚但对于几何形状较为复杂的构

4、件却很难得到准确的结果，甚至根本得不到结果。至根本得不到结果。o有限元分析的目的有限元分析的目的 就是针对具有任意复杂几何形状变形体，完整获取在复杂就是针对具有任意复杂几何形状变形体，完整获取在复杂外力作用下它内部的准确力学信息，即求取该变形体的三类力外力作用下它内部的准确力学信息，即求取该变形体的三类力学信息学信息(位移、应变、应力位移、应变、应力)。4第4页，此课件共60页哦2.1 有限元分析的目的和概念有限元分析的目的和概念o大型液压机机架的设计过程大型液压机机架的设计过程5第5页，此课件共60页哦2.1 有限元分析的目的和概念有限元分析的目的和概念o有限元分析模型与得到的变形状况有限元

5、分析模型与得到的变形状况 6第6页，此课件共60页哦2.1 有限元分析的目的和概念有限元分析的目的和概念o为什么采用有限元方法就可以针对具有任意复杂几何形状的结为什么采用有限元方法就可以针对具有任意复杂几何形状的结构进行分析，并能够得到准确的结果呢？构进行分析，并能够得到准确的结果呢？o这时因为有限元方法是基于这时因为有限元方法是基于“离散逼近离散逼近(discretized approximation)”的基本策略，可以采用较多数量的简的基本策略，可以采用较多数量的简单函数的组合来单函数的组合来“近似近似”代替非常复杂的原函数。代替非常复杂的原函数。7第7页，此课件共60页哦2.1 有限元分

6、析的目的和概念有限元分析的目的和概念o一个复杂的函数，可以通过一系列的基底函数一个复杂的函数，可以通过一系列的基底函数(base function)的组合来的组合来“近似近似”，也就是函数逼近，其中有两，也就是函数逼近，其中有两种典型的方法：种典型的方法：(1)基于全域的展开基于全域的展开(如采用傅立叶级数展开如采用傅立叶级数展开)，(2)基于子域基于子域(sub-domain)的分段函数的分段函数(pieces function)组合组合(如采用分段线性函数的连接如采用分段线性函数的连接)。8第8页，此课件共60页哦【典型例题典型例题】2.1(1)p一个一维函数的两种展开方式的比较一个一维函

7、数的两种展开方式的比较 p 设有一个一维函数，分析它的展开与逼近形式设有一个一维函数，分析它的展开与逼近形式p解答：首先考虑基于全域的展开形式，如采用傅立叶级数解答：首先考虑基于全域的展开形式，如采用傅立叶级数(Fourier series)展开，则有展开，则有(2-1)9第9页，此课件共60页哦【典型例题典型例题】2.1(1)o第二种是基于子域第二种是基于子域 上的分段展开形式，若采用线上的分段展开形式，若采用线性函数性函数,有有o这两种函数展开的方式如图这两种函数展开的方式如图2-2所示。所示。(2-2)10第10页，此课件共60页哦【典型例题典型例题】2.1(1)基于全域的函数展开与逼近

8、基于全域的函数展开与逼近 11第11页，此课件共60页哦【典型例题典型例题】2.1(1)基于子域的函数展开与逼近基于子域的函数展开与逼近 12第12页，此课件共60页哦【典型例题典型例题】2.1(1)o对于第一种的函数逼近方式，就是力学分析中的经典瑞利对于第一种的函数逼近方式，就是力学分析中的经典瑞利-里兹方法里兹方法(Rayleigh-Ritz principle)的思想，而针对的思想，而针对第二种的函数逼近方式，就是现代力学分析中的有限元方法的思第二种的函数逼近方式，就是现代力学分析中的有限元方法的思想，其中的分段就是想，其中的分段就是“单元单元”的概念。的概念。o基于分段的函数描述具有非

9、常明显的优势：基于分段的函数描述具有非常明显的优势：(1)可以将原函可以将原函数的复杂性数的复杂性“化繁为简化繁为简”，使得描述和求解成为可能，使得描述和求解成为可能，(2)所采用的简单函数可以人工选取，因此，可取最简单的线所采用的简单函数可以人工选取，因此，可取最简单的线性函数，或取从低阶到高阶的多项式函数，性函数，或取从低阶到高阶的多项式函数，(3)可以将原可以将原始的微分求解变为线性代数方程。始的微分求解变为线性代数方程。o但分段的做法可能会带来的问题有：但分段的做法可能会带来的问题有：(1)因采用了因采用了“化繁为化繁为简简”，所采用简单函数的描述的能力和效率都较低，所采用简单函数的描

10、述的能力和效率都较低，(2)由于简单函数的描述能力较低，必然使用数量众多的分段来进行由于简单函数的描述能力较低，必然使用数量众多的分段来进行弥补，因此带来较多的工作量。弥补，因此带来较多的工作量。13第13页，此课件共60页哦【典型例题典型例题】2.1(1)o综合分段函数描述的优势和问题，只要采用综合分段函数描述的优势和问题，只要采用功能完善的软件以及能够进行高速处理的计功能完善的软件以及能够进行高速处理的计算机，就可以完全发挥算机，就可以完全发挥“化繁为简化繁为简”策略的策略的优势，有限元分析的概念就在于此。优势，有限元分析的概念就在于此。14第14页，此课件共60页哦2.2 一维阶梯杆结构

11、问题的求解 o【典型例题】2.2(1)15第15页，此课件共60页哦【典型例题】2.2(1)o1D阶梯杆结构问题的材料力学求解阶梯杆结构问题的材料力学求解 16第16页，此课件共60页哦【典型例题】2.2(1)o由杆件的平衡关系可知由杆件的平衡关系可知17第17页，此课件共60页哦【典型例题】2.2(1)o由于材料是弹性的，由虎克定律由于材料是弹性的，由虎克定律(Hooke law)有有 18第18页，此课件共60页哦【典型例题】2.2(1)19第19页，此课件共60页哦【典型例题】2.2(1)20第20页，此课件共60页哦【典型例题】2.2(1)o讨论：以上完全按照材料力学的方法，将对象进行

12、分解来获讨论：以上完全按照材料力学的方法，将对象进行分解来获得问题的解答，它所求解的基本力学变量是力（或应力），得问题的解答，它所求解的基本力学变量是力（或应力），由于以上问题非常简单，而且是静定问题，所以可以直接求由于以上问题非常简单，而且是静定问题，所以可以直接求出，但对于静不定问题，则需要变形协调方程出，但对于静不定问题，则需要变形协调方程(compatibility equation)，才能求解出应力变量，在构，才能求解出应力变量，在构建问题的变形协调方程时，则需要一定的技巧；建问题的变形协调方程时，则需要一定的技巧；o若采用位移作为首先求解的基本变量，则可以使问题的求解变若采用位移作

13、为首先求解的基本变量，则可以使问题的求解变得更规范一些。得更规范一些。o下面就基于下面就基于A、B、C三个点的位移来进行以上问题的求解。三个点的位移来进行以上问题的求解。21第21页，此课件共60页哦【典型例题】2.2(2)o1D阶梯杆结构的节点位移求解及平衡关系阶梯杆结构的节点位移求解及平衡关系 o分别画出每个节点的分离受力图分别画出每个节点的分离受力图22第22页，此课件共60页哦【典型例题】2.2(2)杆杆杆杆23第23页，此课件共60页哦【典型例题】2.2(2)节点节点A节点节点B节点节点C24第24页，此课件共60页哦【典型例题】2.2(2)o将节点将节点A、B、C的平衡关系写成一个

14、方程组的平衡关系写成一个方程组25第25页，此课件共60页哦【典型例题】2.2(2)o改写成矩阵形式改写成矩阵形式26第26页，此课件共60页哦【典型例题】2.2(2)o将材料弹性模量和结构尺寸代入将材料弹性模量和结构尺寸代入(2-23)方程中方程中27第27页，此课件共60页哦【典型例题】2.2(2)o这样得到的结果与这样得到的结果与【典型例题典型例题】2.2(1)所得到的结果完全一所得到的结果完全一致。致。28第28页，此课件共60页哦讨论：讨论：o还可以将式还可以将式(2-23)写成写成o外力矩阵外力矩阵 29第29页，此课件共60页哦讨论：讨论：o内力矩阵内力矩阵30第30页，此课件共

15、60页哦讨论：讨论：o由方程由方程(2-23)可知，这是一个基于节点可知，这是一个基于节点A、B、C描述的全结描述的全结构的平衡方程，该方程的特点为：构的平衡方程，该方程的特点为：o(a)基本的力学参量为节点位移和节点力。基本的力学参量为节点位移和节点力。o o(b)直接给出全结构的平衡方程，而不是象直接给出全结构的平衡方程，而不是象【典型例题典型例题】2.2(1)那样，需要针对每一个杆件去进行递推。那样，需要针对每一个杆件去进行递推。o o(c)在获得节点位移变量后，其它力学参量（如应变和应力）在获得节点位移变量后，其它力学参量（如应变和应力），都可以分别求出。，都可以分别求出。31第31页

16、，此课件共60页哦【典型例题】2.2(3)o1D阶梯杆结构基于位移求解的通用形式阶梯杆结构基于位移求解的通用形式 o为了将方程为了将方程(2-23)写成更规范、更通用的形式，用来求解写成更规范、更通用的形式，用来求解【典型例题典型例题】2.2(1)所示结构的更一般的受力状况，下面在所示结构的更一般的受力状况，下面在式式(2-23)的基础上，直接推导出通用平衡方程。的基础上，直接推导出通用平衡方程。32第32页，此课件共60页哦【典型例题】2.2(3)o将式(2-23)写成 33第33页，此课件共60页哦【典型例题】2.2(3)o再将其分解为两个杆件之和，即写成再将其分解为两个杆件之和，即写成(

17、2-31)34第34页，此课件共60页哦【典型例题】2.2(3)o上式左端的第上式左端的第1项实质为项实质为 o杆件杆件中的左节点的内力和右节点的内力中的左节点的内力和右节点的内力(2-32)35第35页，此课件共60页哦【典型例题】2.2(3)o上式左端的第上式左端的第2项实质为项实质为o杆件杆件中的左节点的内力和右节点的内力中的左节点的内力和右节点的内力(2-33)36第36页，此课件共60页哦【典型例题】2.2(3)o可以看出：方程可以看出：方程(2-31)的左端就是杆件的左端就是杆件的内力表达和杆件的内力表达和杆件的内力表达之和，这样就将原来的基于节点的平衡关系，变为通的内力表达之和，

18、这样就将原来的基于节点的平衡关系，变为通过每一个杆件的平衡关系来进行叠加。过每一个杆件的平衡关系来进行叠加。o这里就自然引入单元的概念，即将原整体结构进行这里就自然引入单元的概念，即将原整体结构进行“分段分段”，以，以划分出较小的划分出较小的“构件构件”，每一个，每一个“构件构件”上具有节点，还可上具有节点，还可以基于节点位移写出该以基于节点位移写出该“构件构件”的内力表达关系，这样的的内力表达关系，这样的“构构件件”就叫做单元就叫做单元(element)，它意味着在几何形状上、节点，它意味着在几何形状上、节点描述上都有一定普遍性和标准性描述上都有一定普遍性和标准性o只要根据实际情况将单元表达

19、式中的参数（如材料常数、几何只要根据实际情况将单元表达式中的参数（如材料常数、几何参数）作相应的代换，它就可以广泛应用于这一类构件参数）作相应的代换，它就可以广泛应用于这一类构件(单元单元)的描述。的描述。37第37页，此课件共60页哦【典型例题】2.2(3)o从式从式(2-32)和式和式(2-33)可以看出，虽然它们分别用来描述杆件可以看出，虽然它们分别用来描述杆件和杆件和杆件的，但它们的表达形式完全相同，因此本质上是一样。的，但它们的表达形式完全相同，因此本质上是一样。o实际上，它们都是杆单元实际上，它们都是杆单元(bar element)。o可以将杆单元表达为如图可以将杆单元表达为如图2

20、-7所示的标准形式。所示的标准形式。38第38页，此课件共60页哦【典型例题】2.2(3)o将单元节点位移写成将单元节点位移写成 o将单元节点外力写成将单元节点外力写成 39第39页，此课件共60页哦【典型例题】2.2(3)o该单元节点内力为该单元节点内力为 o它将与单元的节点外力相平衡它将与单元的节点外力相平衡40第40页，此课件共60页哦【典型例题】2.2(3)o进一步表达为进一步表达为单元内力与外力的单元内力与外力的平衡方程平衡方程(equilibrium equation)，单元的刚度矩阵单元的刚度矩阵(stiffness matrix)刚度矩阵中的刚度刚度矩阵中的刚度系数系数(sti

21、ffness coefficient)41第41页，此课件共60页哦2.3有限元分析的基本流程有限元分析的基本流程 o1D三连杆结构的有限元分析过程三连杆结构的有限元分析过程 42第42页，此课件共60页哦2.3有限元分析的基本流程有限元分析的基本流程o(1)节点编号和单元划分节点编号和单元划分 o(2)计算各单元的单元刚度方程计算各单元的单元刚度方程 43第43页，此课件共60页哦2.3有限元分析的基本流程有限元分析的基本流程o单元单元的刚度方程为的刚度方程为 44第44页，此课件共60页哦2.3有限元分析的基本流程有限元分析的基本流程o单元单元的刚度方程为的刚度方程为 45第45页，此课件

22、共60页哦2.3有限元分析的基本流程有限元分析的基本流程o单元单元的刚度方程为的刚度方程为 46第46页，此课件共60页哦2.3有限元分析的基本流程有限元分析的基本流程o(3)组装各单元刚度方程组装各单元刚度方程 o由于整体结构是由各个单元按一定连接关系组合而成的，因此，由于整体结构是由各个单元按一定连接关系组合而成的，因此，需要按照节点的对应位置将以上方程（需要按照节点的对应位置将以上方程（2-40）、（）、（2-41）、）、（2-42）进行组装，以形成一个整体刚度方程，即）进行组装，以形成一个整体刚度方程，即 47第47页，此课件共60页哦2.3有限元分析的基本流程有限元分析的基本流程48

23、第48页，此课件共60页哦2.3有限元分析的基本流程有限元分析的基本流程49第49页，此课件共60页哦2.3有限元分析的基本流程有限元分析的基本流程o(4)处理边界条件并求解 50第50页，此课件共60页哦2.3有限元分析的基本流程有限元分析的基本流程o求解上述方程，有求解上述方程，有 o(5)求支反力求支反力 由方程组（由方程组（2-46）的最后一行方程）的最后一行方程51第51页，此课件共60页哦2.3有限元分析的基本流程有限元分析的基本流程o(6)求各个单元的其它力学量（应变、应力）52第52页，此课件共60页哦2.3有限元分析的基本流程有限元分析的基本流程53第53页，此课件共60页哦

24、2.3有限元分析的基本流程有限元分析的基本流程o有限元分析的基本流程图示有限元分析的基本流程图示 54第54页，此课件共60页哦2.4有限元分析的特点有限元分析的特点 o有限元分析的最大特点就是标准化和规范化，这种特点使得大规有限元分析的最大特点就是标准化和规范化，这种特点使得大规模分析和计算成为可能，当采用了现代化的计算机以及所编制的模分析和计算成为可能，当采用了现代化的计算机以及所编制的软件作为实现平台时，则复杂工程问题的大规模分析就变为了现软件作为实现平台时，则复杂工程问题的大规模分析就变为了现实。实。o实现有限元分析标准化和规范化的载体就是单元，这就需实现有限元分析标准化和规范化的载体

25、就是单元，这就需要我们构建起各种各样的具有代表性的单元，一旦有了这要我们构建起各种各样的具有代表性的单元，一旦有了这些单元，就好像建筑施工中有了一些标准的预制构件些单元，就好像建筑施工中有了一些标准的预制构件(如如梁、楼板等梁、楼板等)，可以按设计要求搭建出各种各样的复杂结，可以按设计要求搭建出各种各样的复杂结构，如图构，如图2-11所示。所示。55第55页，此课件共60页哦2.4有限元分析的特点有限元分析的特点56第56页，此课件共60页哦2.4有限元分析的特点有限元分析的特点o有限元分析中常用的一些典型单元有限元分析中常用的一些典型单元57第57页，此课件共60页哦2.4有限元分析的特点有

26、限元分析的特点o有限元方法的主要任务就是对常用的各种单元（包括有限元方法的主要任务就是对常用的各种单元（包括1D、2D、3D问题的单元）构造出相应的单元刚度矩阵问题的单元）构造出相应的单元刚度矩阵o当然，如果还采用如当然，如果还采用如【典型例题典型例题】2.2(2)所示的直接法来进行所示的直接法来进行构造，会非常烦琐构造，会非常烦琐o而采用能量原理（如：虚功原理或最小势能原理）来建立而采用能量原理（如：虚功原理或最小势能原理）来建立相应的平衡关系则比较简单，这种方法可以针对任何类型相应的平衡关系则比较简单，这种方法可以针对任何类型的单元进行构建，以得到相应的刚度矩阵。的单元进行构建，以得到相应的刚度矩阵。58第58页，此课件共60页哦2.5 本章要点本章要点n函数的全域逼近与子域分段逼近的概念函数的全域逼近与子域分段逼近的概念 n将对象离散为一些标准的单元将对象离散为一些标准的单元 n将单元的节点位移作为基本力学变量将单元的节点位移作为基本力学变量 n建立单元的节点的平衡关系建立单元的节点的平衡关系 n单元刚度矩阵以及直接求取的方法单元刚度矩阵以及直接求取的方法 n有限元分析的标准流程有限元分析的标准流程(离散化、单元描述、整体组装、离散化、单元描述、整体组装、问题求解问题求解)59第59页，此课件共60页哦结结 束束60第60页，此课件共60页哦