有限元基本方程推导
- 格式:docx
- 大小:16.20 KB
- 文档页数:3
第二章有限单元法的基本原理作为一种比较成熟的数值计算方法,有限元的数学基础是变分原理。
经过半个过世纪的发展,它的数学基础已经比较完善。
从数学角度分析,有限元法是以变分原理和剖分插值为基础的数值计算方法。
它广泛的应用于解算各种类型的偏微分方程,特别对椭圆型方程,因为椭圆型方程的边值问题等价于适当的变分问题,即能量积分的级值问题。
通过变分,导出相应的泛涵,再把作用域从几何上剖分为足够小的单元,这样就能够用简单的图形去拟合复杂的边界,用简单的初等函数去模拟单元的性质。
在解算中先对每个单元进行分析,后在通过连接单元的节点对作用域的整体进行分析,就是对泛涵求极值,从而把一个复杂的偏微分方程求解问题,变成解线形代数方程组的问题。
尽管这样会出现大量的未知数,由于采用了矩阵分析的方法,总体上很有规律,适合编制程序用计算机完成。
通常的数学考虑包括这些: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 是体积力向量。
第二章 有限元分析基本理论有限元法的基本思路是将一个连续求解区域分割成有限个不重叠且按一定方式相互连接在一起的子域(单元),利用在每一个单元内假设的近似函数来分片地表示全求解域上待求的未知场函数。
单元内的场函数通常由未知场函数或其导数在单元各个节点的数值和其插值函数来近似表示。
这样,未知场函数或其导数在各个节点上的数值即成为未知量(自由度)。
根据单元在边界处相互之间的连续性,将各单元的关系式集合成方程组,求出这些未知量,并通过插值函数计算出各个单元内场函数的近似值,从而得到全求解域上的近似解。
有限元将一个连续的无限自由度问题变成离散的有限自由度问题进行求解。
如果将区域划分成很细的网格,也即单元的尺寸变得越来越小,或随着单元自由度的增加及插值函数精度的提高,解的近似程度将不断被改进。
如果单元是满足收敛要求的,近似解最后可收敛于精确解。
2.1 有限元分析的基本概念和计算步骤首先以求解连续梁为例,引出结构有限元分析的一些基本概念和计算步骤。
如图2-1,连续梁承受集中力矩作用。
将结构离散为三个节点,两个单元。
结构中的节点编号为1、2、32.1.1单元分析在有限元分析过程中,第一步是进行结构离散,并对离散单元进行分析,分析的目的是得到单元节点的力与位移的关系。
单元分析的方法有直接法和能量法,本节采用直接法。
从连续梁中取出一个典型单元e ,左边为节点i ,右边为节点j 。
将节点选择在支承点处,单元两端只产生转角位移e i θ、ej θ,顺时针转动为正。
独立的单元杆端内力为弯矩i m 、j m ,顺时针为正。
记:{}e j i eu ⎭⎬⎫⎩⎨⎧=θθ为单元e 的节点位移向量;{}ej i em m f ⎭⎬⎫⎩⎨⎧=为单元e 的杆端力向量。
根据结构力学位移法可得如下平衡方程:⎪⎭⎪⎬⎫+=+=e j e e i e e j ej e e i e e i k k m k k m θθθθ22211211 (2-1)式中:ee e e ee i k k i k k 2412212211====,lEIi e =,EI 、l 分别为单元e 的抗弯刚度和长度。
有限元单元方程推导过程1.引言有限元分析是一种数值计算方法,用于求解结构力学、流体动力学等领域的物理问题。
在有限元分析中,有限元单元是构成整个有限元模型的基本单元,通过推导有限元单元的方程,可以实现对结构或系统的精确分析和计算。
本文将从有限元方法的基本原理出发,详细介绍有限元单元方程的推导过程。
2.有限元方法基本原理有限元方法是将连续的物理问题离散化,转化为有限个代表性元素的集合,通过对每个元素施加适当的边界条件和力学方程,最终得到整个系统的解。
有限元方法通过有限元单元之间的相互作用,从而模拟整个系统的行为。
3.有限元单元的概念有限元单元是有限元模型中最小的离散单元,它是对实际的结构或系统进行离散化的结果。
不同的物理问题和结构,可以采用不同类型的有限元单元进行离散化,如梁单元、壳单元、板单元等。
4.有限元单元方程的一般形式有限元单元方程的一般形式可以表示为:\[K_{e}U_{e}=F_{e}\]其中\(K_{e}\)为有限元单元的刚度矩阵,\(U_{e}\)为有限元单元的位移矢量,\(F_{e}\)为有限元单元的荷载矢量。
5.有限元单元方程推导的基本步骤有限元单元方程的推导主要包括以下几个基本步骤:5.1 单元刚度矩阵的推导首先需要根据有限元单元的几何形状和材料性质,推导出单元刚度矩阵。
单元刚度矩阵可以通过对单元内部的应变能量或者应力-应变关系进行积分得到。
5.2 单元位移矢量的表示在推导单元方程过程中,需要选择合适的位移矢量表示方式,可以采用基函数展开的方法,将位移矢量表示为一组未知系数乘以基函数的线性组合形式。
5.3 单元荷载矢量的求解单元荷载矢量是由外部施加的荷载和边界条件共同决定的,在推导单元方程的过程中需要将这些荷载转化为局部坐标系下的形式,并利用位移矢量的表示方式,将荷载矢量表达为位移矢量和未知系数的线性组合。
5.4 单元方程的组装需要将单元刚度矩阵、位移矢量和荷载矢量组装成完整的单元方程,可以通过坐标变换或者有限元单元之间的关系对单元方程进行组装。