有限元-第1章

有限元-第1章
有限元-第1章

有限元理论基础

有限元理论基础

有限元理论基础 2.1 数值模拟技术 2.1.1数值模拟技术简介 在工程技术领域中许多力学问题和场问题,实质上就是在一定的边界条件下求解一些微分方程。对于少数简单问题,人们可以通过建立它们的微分方程与边界约束求出该问题的解析解。但是对于比较复杂的数学方程问题以及不规则的边界条件通过激吻戏法往往难以求解,而需要借助各种数值模拟方法活的相应的工程数值解,这就是所谓的数值模拟技术。 在实际工程领域中,用数值模拟技术可以对复杂的工程结构进行受力和响应分析,这样可以在设计或者加工前预知实体结构工作状态下的大概情况。 目前在工程实际应用中,常用的数值求解方法有:有限单元法、有限差分法、边界元等但从实用性和使用范围来说,有限单元法则是随着计算机技术的发展而被广泛应用的一种行之有效的数值计算方法。 2.2.2 有限元法 有限元法是一种基于能量原理的数值计算

方法,是解决工程实际问题的一种有效的数值计 算工具。它是里茨法的另一种表示形式,它可应用里茨法分析的所有弹性理论。 限元法是处理连续的结构体离散或有限个单元集合,也就是将连续的求解域离散为一定数量的单元集合体。且每个单元都具有一定的节点,相邻单元通过节点相互连续,同时使用等效节点力代替作用于单元上的力和选定场函数的节点值作为基本未知量。并在每一单元中假设一个近似插值函数以表示单元中场函数的分布规律:进而利用力学中的某些变分原理去建立用以求解节点未知量的有限元法方程,从而将一个连续域中的无限自由度问题化为离散域中的有限自由度问题。求解后,可利用解出的节点值和设定的插值函数确定整个单元集体上的场函数。有限元求解问题中的单元分析:t t t a k F= 式中::t F单元节点作用力。 t K:单元刚度矩阵。 t a:单元节点位移。 通过单元分析确定单元刚度矩阵,建立单元节点作用力和单元为伊关系。有限元求解问题时建立 的结构整体平衡方程:P KU=

有限元分析理论基础

有限元分析概念 有限元法:把求解区域看作由许多小的在节点处相互连接的单元(子域)所构成,其模型给出基本方程的分片(子域)近似解,由于单元(子域)可以被分割成各种形状和大小不同的尺寸,所以它能很好地适应复杂的几何形状、复杂的材料特性和复杂的边界条件 有限元模型:它是真实系统理想化的数学抽象。由一些简单形状的单元组成,单元之间通过节点连接,并承受一定载荷。 有限元分析:是利用数学近似的方法对真实物理系统(几何和载荷工况)进行模拟。并利用简单而又相互作用的元素,即单元,就可以用有限数量的未知量去逼近无限未知量的真实系统。 线弹性有限元是以理想弹性体为研究对象的,所考虑的变形建立在小变形假设的基础上。在这类问题中,材料的应力与应变呈线性关系,满足广义胡克定律;应力与应变也是线性关系,线弹性问题可归结为求解线性方程问题,所以只需要较少的计算时间。如果采用高效的代数方程组求解方法,也有助于降低有限元分析的时间。 线弹性有限元一般包括线弹性静力学分析与线弹性动力学分析两方面。 非线性问题与线弹性问题的区别: 1)非线性问题的方程是非线性的,一般需要迭代求解; 2)非线性问题不能采用叠加原理; 3)非线性问题不总有一致解,有时甚至没有解。 有限元求解非线性问题可分为以下三类:

1)材料非线性问题 材料的应力和应变是非线性的,但应力与应变却很微小,此时应变与位移呈线性关系,这类问题属于材料的非线性问题。由于从理论上还不能提供能普遍接受的本构关系,所以,一般材料的应力与应变之间的非线性关系要基于试验数据,有时非线性材料特性可用数学模型进行模拟,尽管这些模型总有他们的局限性。在工程实际中较为重要的材料非线性问题有:非线性弹性(包括分段线弹性)、弹塑性、粘塑性及蠕变等。 2)几何非线性问题 几何非线性问题是由于位移之间存在非线性关系引起的。 当物体的位移较大时,应变与位移的关系是非线性关系。研究这类问题一般都是假定材料的应力和应变呈线性关系。它包括大位移大应变及大位移小应变问题。如结构的弹性屈曲问题属于大位移小应变问题,橡胶部件形成过程为大应变问题。 3)非线性边界问题 在加工、密封、撞击等问题中,接触和摩擦的作用不可忽视,接触边界属于高度非线性边界。 平时遇到的一些接触问题,如齿轮传动、冲压成型、轧制成型、橡胶减振器、紧配合装配等,当一个结构与另一个结构或外部边界相接触时通常要考虑非线性边界条件。 实际的非线性可能同时出现上述两种或三种非线性问题。

有限元基础

有限元方法基础吉林大学汽车工程学院

引论 一、结构分析的目的 1、克服传统设计方法的不足。 2、为改进设计、优化设计提供指导。 以往的设计大都是基于经验的,尽管基于经验的设计方法在以往的设计中取得了巨大的成功,并且今后的设计仍然将以经验设计为主。但也必须承认基于经验的设计也存在许多不足。首先经验的积累需要时间,有时也不一定都可靠。 此外,通常基于经验的设计能解决行不行的问题,很难解决优不优的问题,尽管为了解决结构优化的问题,也有人提出了一些新的设计理念,如加法设计原则。解决上述问题的根本手段就是采用合理的、科学的方法对上述性能进行分析。

二、数值分析与实验分析的比较 分析方法可分为理论计算和实验测试两大类。 1、基于实验的分析方法 指通过的实验测试获取需要的性能参数的方法。这种方法获取不同的性能参数需要采用不同的测试方法、仪器设备和辅助实验装置。如:强度实验,可以采用电阻应变片及应变仪、光弹涂膜或云纹栅、应变涂料等; 扭转与弯曲刚度实验则需要专门的实验台等等。 实验方法的最大优点是工程上普遍认为,实验结果真实可靠,通常被当作产品最终定型的权威性依据。

实验方法的优点: 是工程上普遍认为,实验结果真实可靠, 通常被当作产品最终定型的权威性依据。 实验方法的缺点: 1)实验一定要在样品或样机试制之后才能进行,无法与设计并行。 2)成本高、周期长,通常只适合批量生产的产品。 3)可以获得的数据量有限,无法对设计提供更多的指导,更无法进行结构优化。 4)受实验手段的限制,有些参数无法测准。

2、基于理论计算的分析方法 指通过理论分析或数值计算获取所需的性能参数的分析方法 优点: 1)经济、快捷,成本低、周期短。(与实验相比)2)一次分析可以获得大量的数据。 3)可以与设计同步进行。 4)可以配合优化算法,对设计进行优化。 缺点:分析结果受模型质量、算法理论及边界条件等诸多因素影响。

有限元分析理论基础

有限元分析概念 有限元法:把求解区域瞧作由许多小的在节点处相互连接的单元(子域)所构成,其模型给出基本方程的分片(子域)近似解,由于单元(子域)可以被分割成各种形状与大小不同的尺寸,所以它能很好地适应复杂的几何形状、复杂的材料特性与复杂的边界条件 有限元模型:它就是真实系统理想化的数学抽象。由一些简单形状的单元组成,单元之间通过节点连接,并承受一定载荷。 有限元分析:就是利用数学近似的方法对真实物理系统(几何与载荷工况)进行模拟。并利用简单而又相互作用的元素,即单元,就可以用有限数量的未知量去逼近无限未知量的真实系统。 线弹性有限元就是以理想弹性体为研究对象的,所考虑的变形建立在小变形假设的基础上。在这类问题中,材料的应力与应变呈线性关系,满足广义胡克定律;应力与应变也就是线性关系,线弹性问题可归结为求解线性方程问题,所以只需要较少的计算时间。如果采用高效的代数方程组求解方法,也有助于降低有限元分析的时间。 线弹性有限元一般包括线弹性静力学分析与线弹性动力学分析两方面。 非线性问题与线弹性问题的区别: 1)非线性问题的方程就是非线性的,一般需要迭代求解; 2)非线性问题不能采用叠加原理; 3)非线性问题不总有一致解,有时甚至没有解。 有限元求解非线性问题可分为以下三类:

1)材料非线性问题 材料的应力与应变就是非线性的,但应力与应变却很微小,此时应变与位移呈线性关系,这类问题属于材料的非线性问题。由于从理论上还不能提供能普遍接受的本构关系,所以,一般材料的应力与应变之间的非线性关系要基于试验数据,有时非线性材料特性可用数学模型进行模拟,尽管这些模型总有她们的局限性。在工程实际中较为重要的材料非线性问题有:非线性弹性(包括分段线弹性)、弹塑性、粘塑性及蠕变等。 2)几何非线性问题 几何非线性问题就是由于位移之间存在非线性关系引起的。 当物体的位移较大时,应变与位移的关系就是非线性关系。研究这类问题一般都就是假定材料的应力与应变呈线性关系。它包括大位移大应变及大位移小应变问题。如结构的弹性屈曲问题属于大位移小应变问题,橡胶部件形成过程为大应变问题。 3)非线性边界问题 在加工、密封、撞击等问题中,接触与摩擦的作用不可忽视,接触边界属于高度非线性边界。 平时遇到的一些接触问题,如齿轮传动、冲压成型、轧制成型、橡胶减振器、紧配合装配等,当一个结构与另一个结构或外部边界相接触时通常要考虑非线性边界条件。 实际的非线性可能同时出现上述两种或三种非线性问题。 有限元理论基础

有限元法的理论基础

有限元法的理论基础 有限元法是一种离散化的数值计算方法,对于结构分析而言,它的理论基础是能量原理。能量原理表明,在外力作用下,弹性体的变形、应力和外力之间的关系受能量原理的支配,能量原理与微分方程和定解条件是等价的。下面介绍有限元法中经常使用的虚位移原理和最小势能原理。 1.虚位移原理 虚位移原理又称虚功原理,可以叙述如下:如果物体在发生虚位移之前所受的力系是平衡的(物体内部满足平衡微分方程,物体边界上满足力学边界条件),那么在发生虚位移时,外力在虚位移上所做的虚功等于虚应变能(物体内部应力在虚应变上所做的虚功)。反之,如果物体所受的力系在虚位移(及虚应变)上所做的虚功相等,则它们一定是平衡的。可以看出,虚位移原理等价于平衡微分方程与力学边界条件。所以虚位移原理表述了力系平衡的必要而充分的条件。 虚位移原理不仅可以应用于弹性性力学问题,还可以应用于非线性弹性以及弹塑性等非线性问题。 2.最小势能原理 最小势能原理可以叙述为:弹性体受到外力作用时,在所有满足位移边界条件和变形协调条件的可以位移中,真实位移使系统的总势能取驻值,且为最小值。根据最小势能原理,要求弹性体在外力作用下的位移,可以满足几何方程和位移边界条件且使物体总势能取最小值的条件去寻求答案。最小势能原理仅适用于弹性力学问题。 2.2有限元法求解问题的基本步骤 弹性力学中的有限元法是一种数值计算方法,对于不同物理性质和数学模型的问题,有限元法的基本步骤是相同的,只是具体方式推导和运算求解不同,有限元求解问题的基本步骤如下。 2.2.1问题的分类 求解问题的第一步就是对它进行识别分析,它包含的更深层次的物理问题是什么?比如是静力学还是动力学,是否包含非线性,是否需要迭代求解,要从分析中得等到什么结果等。对这些问题的回答会加深对问题的认识与理解,直接影响到以后的建模与求解方法的选取等。 2.2.2建模 在进行有限元离散化和数值求解之值,我们为分析问题设计计算模型,这一步包括决定哪种特征是所要讨论的重点问题,以便忽略不必要的细节,并决定采用哪种理论或数学公式描述结果的行为。因此,我们可以忽略几何不规则性,把一些载荷看做是集中载荷,并把某些支撑看做是固定的。材料可以理想化为线弹性和各向同性的。根据问题的维数、载荷以及理论化的边界条件,我们能够决定采用梁理论、板弯曲理论、平面弹性理论或者一些其他分析理论描述结构性能。在求解中运用分析理论简化问题,建立问题的模型。 2.2.3连续体离散化 连续体离散化,习惯上称为有限元网络划分,即将连续体划分为有限个具有规则形状的单元的集合,两相邻单元之间只通过若干点相互连接,每个连接点称为节点。单元节点的设置、性质、数目等应视问题的性质、描述变形的需要和计算精度而定,如二维连续体的单元可为三角形、四边形,三维连续体的单元可以是四面体、长方体和六面体等。为合理有效地表示连续体,需要适当选择单元的类型、数目、大小和排列方式。 离散化的模型与原来模型区别在于,单元之间只通过节点相互连接、相互作用,而无其他连接。因此这种连接要满足变形协调条件。离散化是将一个无限多自由度的连续体转化为一个有限多自由度的离散体过程,因此必然引起误差。主要有两类:建模误差和离散化误差。

有限元理论基础

有限元理论基础 2.1 数值模拟技术 2.1.1数值模拟技术简介 在工程技术领域中许多力学问题和场问题,实质上就是在一定的边界条件下求解一些微分方程。对于少数简单问题,人们可以通过建立它们的微分方程与边界约束求出该问题的解析解。但是对于比较复杂的数学方程问题以及不规则的边界条件通过激吻戏法往往难以求解,而需要借助各种数值模拟方法活的相应的工程数值解,这就是所谓的数值模拟技术。 在实际工程领域中,用数值模拟技术可以对复杂的工程结构进行受力和响应分析,这样可以在设计或者加工前预知实体结构工作状态下的大概情况。 目前在工程实际应用中,常用的数值求解方法有:有限单元法、有限差分法、边界元等但从实用性和使用范围来说,有限单元法则是随着计算机技术的发展而被广泛应用的一种行之有效的数值计算方法。 2.2.2 有限元法 有限元法是一种基于能量原理的数值计算方法,是解决工程实际问题的一种有效的数值计算工具。它是里茨法的另一种表示形式,它可应用里茨法分析的所有弹性理论。 限元法是处理连续的结构体离散或有限个单元集合,也就是将连续的求解域离散为一定数量的单元集合体。且每个单元都具有一定的节点,相邻单元通过节点相互连续,同时使用等效节点力代替作用于单元上的力和选定场函数的节点值作为基本未知量。并在每一单元中假设一个近似插值函数以表示单元中场函数的分布规律:进而利用力学中的某些变分原理去建立用以求解节点未知量的有限元法方程,从而将一个连续域中的无限自由度问题化为离散域中的有限自由度问题。求解后,可利用解出的节点值和设定的插值函数确定整个单元集体上的场函数。 有限元求解问题中的单元分析:t t t a k F = 式中::t F 单元节点作用力。 t K :单元刚度矩阵。 t a :单元节点位移。 通过单元分析确定单元刚度矩阵,建立单元节点作用力和单元为伊关系。有限元求解问题时建立的结构整体平衡方程:P KU = 式中:P —结构整体等效点力载荷 K —结构总体刚度矩阵 U —结构节点位移阵列 单元内力的计算:t DBa =σ 式中:D —弹性矩阵 P —应变矩阵 整个结构的有限元分析就是一句上述方程而进行的具体的有限元求解过程如图

有限元的基础理论

§1有限元的基础理论 §1-1 概述 有限元法是一种数值计算的近似方法。早在40年代初期就已有人提出,但当时由于没有计算工具而搁置,一直到50年代中期,高速数字电子计算机的出现和发展为有限元法的应用提供了重要的物质条件,才使有限元法得以迅速发展。 有限元法在西方起源于飞机和导弹的结构设计,发表这方面文章最早而且最有影响的是西德的J.H.Argyris教授,于1954–1955年间,他在《Aircraft engineering》上发表了许多有关这方面的论文,并在此基础上写成了《能量原理与结构分析》,此书成为有限元法的理论基础。美国的M.T.Turner,R.W.Clough,H.C.Martin和L.J.Topp等人于1956年发表了一篇题为《复杂结构的刚度和挠度分析》一文,此文提出了计算复杂结构刚度影响系数的方法,说明了如何利用计算机进行分析。美国教授R.W.Clough于1960年在一篇介绍平面应力分析的论文中,首次提出了有限元法的名字。1965年英国的O.C.Zienliewice教授及其合作者解决了将有限元应用于所有场的问题,使有限元法的应用范围更加广泛。 有限元法的优点很多,其中最突出的优点是应用范围广。发展至今,不仅能解决静态的、平面的、最简单的杆系结构,而且还可以解决空间问题、板壳问题、结构的稳定性问题、动力学问题、弹塑性问题和粘弹性问题、疲劳和脆性断裂问题以及结构的优化设计问题。而且不论物体的结构形式和边界条件如何复杂,也不论材料的性质和外载荷的情况如何,原则上都能应用。 §1-2 有限元的基础理论 有限元法的基本思路和基本原则以结构力学中的位移法为基础,把复杂的结构或连续体看成有限个单元的组合,各单元彼此在节点处连接而组成整体。把连续体分成有限个单元和节点,称为离散化。先对单元进行特性分析,然后根据各节点处的平衡和协调条件建立方程,综合后作整体分析。这样一分一合,先离散再综合的过程,就是把复杂结构或连续体的计算问题转化为简单单元的分析与综合的问题。因此,一般的有限元解法包括三个主要步骤:离散化、单元分析、整体分析。 §1-2-1 离散化 一个复杂的弹性体可以看作由无限个质点组成的连续体。为了进行解算,可以将此弹性体简化为有限个单元组成的集合体,这些单元只在有限个节点上铰接,因此,这集合体只具有有限个自由度,这就为解算提供了可能。有无限个质点的连续体转化为有限个单元的集合体,就称为离散化。 §1-2-2 单元分析 单元分析首先要进行单元划分。在工程结构中,一般采用四种类型的基本单元,即标量单元、线单元(杆、梁单元)、面单元和体单元。四种基本单元的若干例子及各单元节点自由度(节点位移)表示在图(1-1)中。而单元划分一般注意下面几点: 一、从有限元本身来看,单元划分的越细,节点布置得越多,计算的结果越精确。但计算时间和计算费用的增加。所以在划分单元时对应兼顾这两个方面。 二、在边界比较曲折,应力比较集中,应力变化较大的地方,单元应划分的细点,而在应力变化平缓处单元划分的大些。单元由小到大应逐渐过渡。 三、对于三角形单元,三条边长应尽量接近,不应出现钝角,以免计算出现较大的偏差。对于矩形单元,长度和宽度也不应相差过大。 四、任意一个三角形单元的角点必须同时也是相邻单元边上的角点,而不能是相邻单元边上的内点。划分其他单元时也应遵循此原则。

有限元基础理论复习资料

有限元基础理论复习 第一章:有限元法及ANSYS概述 1.CAE的概念是什么?(P1) CAE即计算机辅助工程,指工程设计中的分析计算与仿真。 2.有限单元法的基本思想是什么?(P2) 有限单元法的基本思想是将物体(即连续的求解域)离散成有限个且按一定方式相互联结在一起的单元的组合,来模拟或逼近原来的物体,从而将一个连续的无限自由度问题简化为离散的有限自由度问题求解的一种数值分析法。 3.单元、节点概念的定义是什么?(P2) 网格划分中每一个小的块体称为单元。 确定单元形状、单元之间相互联结的点称为节点。 4.节点力与节点载荷的区别是什么?(P2) 单元上节点处的结构内力为节点力,外力(有集中力、分析力等)为节点载荷。故一个是内力,一个是外力。 第二章:有限元法基础理论 1.平面应力问题与平面应变问题的区别是什么?(P25) 恒有δz=0,τzx =τxz=0, τzy=τyz=0,不为0的应力分量为δx,δy,τxy,这种问题称为平面应力问题。 恒有w=0,εz=γyz=γzx=0,不为0的应力分量为εx,εy,γxy,这种问题就称为平面应变问题。 2.轴对称问题有什么特征?它和平面应力问题的主要区别是什么?(P34) 轴对称应力问题的特征是如果弹性体的几何形状、约束条件及载荷都对称于某一轴,则所有的位移、应变及应力也对称于此轴。 与平面应力问题不同的是:单元体为圆环体,单元之间由结圆铰接,节点力为结圆上的均布力,单元边界为回转面。 3.什么是等参数单元?(P40) 等参数变换即坐标变换和单元内德场函数采用相同数目的节点参数及相同的插值函数,等参数变换的单元称之为等参数单元。 4.介绍虚位移原理和最小势能原理?(P44) 虚位移原理:如果在虚位移发生之前,物体处于平衡状态,那么在虚位移发生时,外力所做的虚功等于物体的虚应变能。 最小势能原理:在所有满足边界条件的协调(连续)位移中,那些满足平衡条件的位移使物体势能取驻值,即δПp=δU-δV=0,对于线性弹性体,势能取最小值。 5.计算题 第三章:ANSYS建模 1.什么是工作平面?(P97) 光标在屏幕上是一个点,在空间实际上代表一条直线。为了用光标拾取一个点,必须要有一个假想的平面与该直线相交,这样才能唯一地确定空间中的一个点,该平面即工作平面。 2.ANSYS内有哪几种坐标系,适用哪些场合?(P 95) 1.全局和局部坐标系统:在空间坐标定位几何项(如节点和关键点) 2.显示坐标系统:决定列出和显示几何项的坐标系。 3.节点坐标系统:定义每个节点自由度的方向及节点计算结果的定位。 4.单元坐标系统:定位材料特性及单元计算结果数据。

第一章:有限元基本概念

第一章 有限元方法的基本原理与概念 有限元方法是求解数学物理问题的一种强有力的工具,它在变分法理论的基础上吸收了有限差分方法的思想而发展起来的。随着计算求解技术的发展成为研究学研究和工程设计等部门中的一种有效的手段。为了使读者更好的了解这一方法,本章主要介绍数学物理问题的一般提法,以及有限元方法所能解决问题的一些背景及其应用范围。本章还通过对有限差分法、变分法和加权残数法的介绍和比较,让读者进一步了解有限元方法的基本原理及其优越性。 §1.1 偏微分方程与定解问题 自然界中的许多物理问题在某些情况下可以通过适当的微分方程或偏微分方程描述,通常把这类方程称为数学物理方程。这些方程在给定的特定边界条件和初始条件中可求得解答。数学物理方程的求解方法基本可以分为两大类型:一类称之为解析法,例如分离变量法、特殊函数法、积分变换法等,利用这种方法可得到问题的精确解;另一类为数值法,例如有限差分法、数值积分方法、边界元法、有限元法等,利用这种方法可得到问题的近似解。不论使用哪种方法都需要对数学物理方程及其定解条件(初值和边值)应有基本了解,故本节对此方面的基本概念作一简要介绍。 §1.1.1偏微分方程的分类 数学物理问题一般由偏微分方程控制。方程中偏导数的最高阶次称之为方程的阶;未知函数多于一个的偏微分方程组成的方程组称为偏微分方程组。如果方程中包含的其它函数或其偏导数的项都是未知函数的一次式,则称之为线性偏微分方程;否则为非线性偏微分方程。在非线性方程中,未知函数的最高阶偏导数对于未知函数是线性的称之为拟线性偏微分方程,如果方程中存在有与未知函数无关项称之为非齐次的,否则为齐次方程。 在数学物理方程中最为典型的一维波动方程 222 2 2(,)u u a f x t x t ??+=?? (1.1.1) 为二阶非齐次线性方程。

有限元法的基础理论

一、里兹法与迦辽金法(摘自电磁场有限元方法 金建铭) 1. 里兹法 里兹法是一种变分方法,其中边值问题用变分表达式(也称泛函)表示,泛函的极小值对应于给定边界条件下的控制微分方程。通过求泛函相对于其变量的极小值可得到近似解。 2. 伽辽金法 伽辽金法属于残数加权方法类型,它通过对微分方程的残数求加权的方法得到方程的解。 若u 是方程的近似解,将u 代入方程可得到非零的残数: r Lu f =- u 的最佳近似应能使残数r 在Ω内所有点上有最小值。残数加权方法要求: 0i i R rd ωΩ =Ω=? 这里i R 表示残数的加权积分,i ω是所选的加权函数。 在伽辽金法中,加权函数与近似解展开中所用的函数相同。通常,这样可得到最精确的 解。 二、有限元方法 里兹法和伽辽金法中,在整个解域内找出能表示或至少近似表示问题真实解的试探函数是非常重要的。然而对于许多问题,这个步骤是十分困难的,对二维和三维问题尤其如此。为此,我们可将整个区域划分成小子域,并应用定义在每个子域上的试探函数。因为子域是小区域,因而在每一子域内函数的变化不大,所以定义在子域上的试探函数通常比较简单。这正是有限元法的基本思想。应用里兹法的过程通常称为里兹有限元法或变分有限元法,而应用伽辽金方法的过程通常称为伽辽金有限元方法。 有限元法与经典里兹法和伽辽金法的不同之处是在试探函数的公式上。在经典里兹法和伽辽金法中,试探函数由定义在全域上的一组基函数组成。这种组合必须能够(至少近似)表示真实解,也必须满足适当的边界条件。在有限元法中,试探函数是由定义在组成全域的子域上的一组基函数构成。因为子域很小,所以定义在子域上的基函数能够十分简单。 三、关于形函数(摘自有限元法在电磁计算中的应用 张榴晨) 对于一个待求的微分方程,用一组线性独立的尝试函数i ψ和待定系数i C 来表示方程的近似解,并用加权余数法(迦辽金法)来求解这些待定系数。求解待定系数的代数方程组为: 1 []1,2,,n i j i j i d C q d j n ψψψΩ Ω =??Ω=Ω =∑? ? 这里j ψ为所选择的加权函数,应用迦辽金法时,所选取的加权函数即为尝试函数。 有限元中应用的尝试函数代表了单元上近似解的一种插值关系,它决定了近似解在单元上的形状。因此尝试函数在有限元法中又称为形函数。对于一维有限元来说,形函数为一个直线段;对一维高阶有限元来说,形函数为一个曲线段;对二维一阶有限元来说,形函数为一个平面;对二维高阶有限元来说,形函数为一个曲面;三维有限元来说,形函数为多维平面或曲面。选择形函数时可以使一个任意元上的函数只与该元所对应的节点势函数值有关,而与其它各点的值无关。 1. 一维有限元

相关文档
最新文档