有限元法数学基础
- 格式:ppt
- 大小:2.95 MB
- 文档页数:33
材料力学有限元法知识点总结材料力学是一门研究物质内部结构、性质和变形行为的学科,而有限元法则是一种在工程和科学领域中广泛应用的数值计算方法。
有限元法可以将一个复杂的实体划分为无数小的单元,通过对这些小单元进行分析和计算,最终得到整个实体的力学性质和行为。
本文将对材料力学有限元法的一些核心概念和知识点进行总结。
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):有限元法通过将变分原理中的势能函数展开为一系列基函数的线性组合,并使用权重函数对残差进行加权求和的方式进行近似。
这样可以将求解连续问题转化为离散问题,进而进行数值计算。
有限单元法的数学基础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 φφφφφφφφφφφφφφφφφφ⎛⎫ ⎪ ⎪= ⎪ ⎪⎝⎭,,,,,, 。
有限元基础知识
嘿,朋友们!今天咱要来聊聊有限元基础知识啊,这可真是个超有意思的东西!
你们有没有玩过拼图游戏呀?有限元就有点像把一个复杂的东西,比如一个机器零件啦,拆分成好多好多小的部分,就像拼图的小块块一样。
比如说,你想想看一辆汽车,它那么复杂,要是直接去研究它可太难了。
但通过有限元,咱就可以把它分成一个个小区域,分别去分析、去理解,这不就简单多了嘛!
有限元就像是给我们一个探索复杂世界的秘密武器!它让那些看似遥不可及、搞不懂的东西变得清晰起来。
你知道吗?工程师们经常用这个方法来解决各种各样的问题呢!比如设计更牢固的桥梁,或者让飞机飞得更安全、更稳定。
就好比有一座摇摇欲坠的老桥,工程师们就可以用有限元方法,一点一点地分析每个部分,找出问题所在,然后想办法加固它,让它重新变得坚固可靠。
这多了不起啊!
那有限元具体是咋工作的呢?简单来说,就是先划分网格,这就像是给那个复杂的东西画格子。
然后再对这些小格子进行计算和分析。
就好像你在做数学题一样,一步步算出答案。
“哎呀,这听起来好难啊!”你可能会这么说。
但别害怕呀!一开始可能觉得有点难理解,但只要你深入进去,就会发现它的奇妙之处。
而且现在有好多软件可以帮我们进行有限元分析呢,超方便的!
总之,有限元基础知识是个非常有用、非常有趣的东西!它就像一把钥匙,能帮我们打开复杂工程世界的大门,让我们更好地去理解和创造。
大家赶紧去探索一下吧,相信你们一定会爱上它的!。
第二章有限单元法的基本原理作为一种比较成熟的数值计算方法,有限元的数学基础是变分原理。
经过半个过世纪的发展,它的数学基础已经比较完善。
从数学角度分析,有限元法是以变分原理和剖分插值为基础的数值计算方法。
它广泛的应用于解算各种类型的偏微分方程,特别对椭圆型方程,因为椭圆型方程的边值问题等价于适当的变分问题,即能量积分的级值问题。
通过变分,导出相应的泛涵,再把作用域从几何上剖分为足够小的单元,这样就能够用简单的图形去拟合复杂的边界,用简单的初等函数去模拟单元的性质。
在解算中先对每个单元进行分析,后在通过连接单元的节点对作用域的整体进行分析,就是对泛涵求极值,从而把一个复杂的偏微分方程求解问题,变成解线形代数方程组的问题。
尽管这样会出现大量的未知数,由于采用了矩阵分析的方法,总体上很有规律,适合编制程序用计算机完成。
通常的数学考虑包括这些:1)从古典变分方法原理去定义微分方程边值问题的广义解以及在古典变分方法的框架对有限元进行理论分析。
2)保证偏微分方程边值问题的提法正确,即要求解存在、唯一和稳定,即保证数值解法是可靠的。
3)有限元中重要的一点是采用了分块多项式插值函数,因此,有限元的误差估计转化为插值逼近的误差估计问题。
4)有限元的收敛性和误差估计。
由于本文是应用有限元的理论解决大地测量中的问题,因此,这里将不讨论上叙问题,而是从固体力学的基本方程出发,通过虚功原理建立起离散化的有限元方程。
另外,还以八节点六面体单元为例,简要叙述了实际中最常用的等参单元的概念及其数值变化的一些公式。
§2.1 弹性力学基本方程有限元法中经常要用到弹性力学的基本方程,这里写出这些方程的矩阵表达式。
2-1-1、平衡方程对任意一点的受力情况分析,沿坐标轴方向x, y ,z分解得到平衡方程0*00000000=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎣⎡∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂z y xxz yz xy z y x F F F z yzz x y z y x τττσσσ 记为: 0=+F A σ其中A 是微分算子,F 是体积力向量。
一:有限元的基本思想有限元法的基本思想是将连续的求解区域离散为一组有限个、且按一定方式相互联结在一起的单元的组合体。
由于单元能按不同的联结方式进行组合,且单元本身又可以有不同形状,因此可以模型化几何形状复杂的求解域。
通常有限元法都遵循以下基本步骤:物体的离散化:离散化是有限元法的基础,这就是依据结构的实际情况,选择合适的单元形状、类型、数目、大小以及排列方式,将拟分析的物体假想地分成有限个分区或分块的集合体。
假设这些单元在处于它们边界上的若干个离散节点处相互连接,这些节点的位移将是该问题的基本未知参数。
挑选形函数或插值函数:选择一组函数,通常是多项式,最简单的情况是位移的线性函数。
这些函数应当满足一定条件,该条件就是平衡方程,它通常是通过变分原理得到的,可由每个“有限单元”的节点位移唯一地确定该单元中的位移状态。
确定单元的性质:确定单元性质就是对单元的力学性质进行描述。
确定了单元位移后,可以很方便地利用几何方程和物理方程求得单元的应变和应力。
一般用单元的刚度矩阵来描述单元的性质,确定单元节点力与位移的关系。
组成物体的整体方程组:组成物体的整体方程组就是由已知的单元刚度矩阵和单元等效节点载荷列阵集成表示整个物体性质的结构刚度矩阵和结构载荷列阵,从而建立起整个结构己知量-------总节点载荷与整个物体未知量-------总节点位移的关系。
解有限元方程和辅助计算:引入强制边界条件,解方程得到节点位移。
一般整体方程组往往数目庞大,可能是几十个、几百个,以至于成千上万个。
对于这些方程组需要一定的计算数学方法解出其未知量。
然后,根据实际问题进行必要的辅助计算。
完整的有限元的求解过程如下图所示:二:有限元的数学方法从更广泛的观点看,有限元法的数学基础是变分原理。
根据变分原理发展而来的有限元法,在求解微分方程方面得到了广泛的应用。
变分原理是表达物理基础定律的一种普遍形式,其表达课概括如下:给出一个依赖物理状态v 的变量()J v ,同时给出()J v 的容许函数集v ,即一切可能的物理状态,则真是的状态是v 中使()J v 达到极小值的函数。