元法概念意义与应用

合集下载

有限元法

有限元法

李中秋20111323 热能一班第一章有限元法简介有限元法是求取复杂微分方程近似解的一种非常有效的工具,是现代数字化科技的一种重要基础性原理。

将它用于在科学研究中,可成为探究物质客观规律的先进手段。

将它应用于工程技术中,可成为工程设计和分析的可靠工具。

1.1 有限元法发展简史早在1870年,英国科学家Rayleigh就采用假想的“试函数”来求解复杂的微分方程,1909年Ritz将其发展成为完善的数值近似方法,为现代有限元方法打下坚实基础。

20世纪40年代,由于航空事业的飞速发展,设计师需要对飞机结构进行精确的设计和计算,便逐渐在工程中产生了的矩阵力学分析方法;1943年,Courant 发表了第一篇使用三角形区域的多项式函数来求解扭转问题的论文;1956年波音公司的Turner,Clough,Martin和Topp在分析飞机结构时系统研究了离散杆、梁、三角形的单元刚度表达式;1960年Clough在处理平面弹性问题,第一次提出并使用“有限元方法”(finite element met hod)的名称;1955年德国的Argyris出版了第一本关于结构分析中的能量原理和矩阵方法的书,为后续的有限元研究奠定了重要的基础,1967年Zienkiewicz和Cheung出版了第一本有关有限元分析的专著;1970年以后,有限元方法开始应用于处理非线性和大变形问题;我国的一些学者也在有限元领域做出了重要的贡献,如胡海昌于1954提出了广义变分原理[8],钱伟长最先研究了拉格朗日乘子法与广义变分原理之间关系,钱令希在20世纪五十年代就研究了力学分析的余能原理,冯康在20世纪六十年代就独立地、并先于西方奠定了有限元分析收敛性的理论基础。

1.2基本概念1.2.1 有限单元数值计算的思路是将复杂问题简单化,求近似解。

即将复杂的结构分解成若干相对简单的构件或部件,分别分析,然后求解。

而且这种近似解可以收敛于问题的精确解。

二元函数求极限的积分换元法综述

二元函数求极限的积分换元法综述

二元函数求极限的积分换元法综述在高等数学中,求二元函数的极限是一个非常重要的概念。

对于一些复杂的函数,直接求解其极限可能会比较困难。

而积分换元法是一种常用的有效方法,可以简化二元函数极限的求解过程。

本文将对积分换元法在求解二元函数的极限中的应用进行综述。

一、积分换元法简介积分换元法是一种常用的积分求解技巧,它通过引入新的变量替代原变量,从而将原积分转化为更加简单的形式。

在二元函数求极限中,我们可以借鉴积分换元法的思想,将原二元函数转化为与之等价的更容易求解的函数形式。

二、二元函数求极限的积分换元法步骤1. 确定变量替换对于给定的二元函数,我们首先需要确定合适的变量替换。

通常情况下,我们选择将二元函数中的一个自变量表示为另一个自变量的函数形式。

2. 进行变量替换根据确定的变量替换,我们将原二元函数中的自变量进行对应的替换。

这样可以将原二元函数转化为只含有一个变量的函数。

3. 求解极限通过变量替换,我们得到了一个只含有一个变量的函数。

接下来,我们可以使用常规的一元函数求极限的方法,对这个函数进行求解。

4. 还原变量在求解极限后,我们需要将之前引入的新变量还原为原二元函数的自变量。

这样可以得到最终的极限结果。

三、实例分析以求解二元函数 f(x,y) = sin(x^2 + y^2) / (x^2 + y^2) 在点 (0,0) 处的极限为例,综合使用积分换元法进行求解。

1. 确定变量替换我们可以将 x^2 + y^2 表示为 r^2,其中 r 表示点 (x, y) 到原点的距离。

2. 进行变量替换根据变量替换 r^2 = x^2 + y^2,我们将原二元函数中的自变量进行替换。

这样可以得到性质更简单的新的函数 f(r) = sin(r^2) / r^2。

3. 求解极限通过变量替换,我们将二元函数的极限转化为一元函数的极限。

对新函数 f(r) 使用一元函数求极限的方法,我们得到lim(r→0) f(r) = 1。

偏微分方程的有限元法

偏微分方程的有限元法
第3页/共106页
第五章 偏微分方程的有限元法
有限元法特点有限元法的物理意义直观明确,理论完整可靠。 因为变分原理描述了支配物理现象的物理学中的最小作用原理(如力学中的最小势能原理)。 优异的解题能力。有限元法对边界几何形状复杂以及媒质物理性质变异等复杂物理问题求解上,有突出优点: ① 不受几何形状和媒质分布的复杂程度限制。 ②不必单独处理第二、三类边界条件。 ③ 离散点配置比较随意,通过控制有限单元剖分密度和单元插值函数的选取,可以充分保证所需的数值计算精度。
有限元法于上世纪50年代首先在力学领域-----飞机结构的静、动态特性分析中得到应用,随后很快广泛的应用于求解热传导、电磁场、流体力学等连续性问题。有限元法主要用于求解拉普拉斯方程和泊松方程所描述的各类物理场中。
第1页/共106页
第五章 偏微分方程的有限元法
有限元法---变分原理
第4页/共106页
5.1 泛函与变分原理
数学上,通常自变量与因变量间的关系称为函数,而泛函则是函数集合的函数,也就是函数的函数,即自变量为函数,而不是变量。
5.1.1 泛函的定义 泛函通常是指一种定义域为函数,而值域为实数的“函数”。 设C是函数的集合,B是实数集合。如果对C中的任一元素y(x),在B中都有一个元素J与之对应,则称J为y(x)的泛函,记为J[y(x)]。
5.1.3 泛函的变分
定义最简泛函
F(x,y,y’)称为泛函的“核函数”
泛函的变分
最简泛函: 核函数只包含自变量 x、未知函数y(x)以及导数y’(x)
第9页/共106页
5.1 泛函与变分原理
利用二元函数的泰勒展开
第10页/共106页
5.1 泛函与变分原理
其中
分别称为泛函的一阶变分和二阶变分。

元法的讲解

元法的讲解

第四章求解导热问题的有限单元法第4.1节概述第4.2节泛函变分原理第4.3节有限单元法第4.1节概述粗略地讲:有限元法是获得微分方程近似解的一种方法,是一种适合计算机来求解的数值计算方法。

(元素特性方程和总体合成方程的建立可以采用直接法,变分法,加权余数法和能量平衡法等四种方法之一,所以粗略地说有限元法是获得微分方程近似解的一种方法也有道理)比较严格的定义:有限单元法是求解泛函变分问题的一种近似方法。

那么这两种说法有什么联系,或者说是共同之处呢?变分和微分是对未知函数的不同描述,同一连续介质问题往往都可以找到微分和变分的等价表达方式。

变分和微分几乎是同时发展起来的两个数学分支,其目的是相同的,都是求解未知函数,但是方法上有很大差别。

在已知边界条件的情况下,求微分方程的精确解析虽然已有完整的理论,但是真正能解出的只有极少数的几种简单情况,因为在很多情况下,微分方程并不存在初等函数解析解。

(对于各种各样的映射,初等函数的表达能力实在太有限了,初等函数包括:冥函数、指数函数、对数、三角函数,以及它们的四则运算等。

)由于寻求微分方程的初等函数解析解有困难,所以我们在前一章讲述了微分方程的近似解法,即差分法。

泛函变分原理虽然也可以用解析法(即积分)求得未知函数,但是因为有很多被积函数根本无法找到初等原函数,也就不能积分,尤其是对于二维和三维问题,解析法更加困难。

所以我们也要寻求泛函变分的近似解法。

泛函变分的近似解法包括里兹法和有限元法(里兹法是有限元法的前身),这两种方法的原理完全相同,即:构造一个近似的初等函数,用近似的初等函数去逼近未知函数。

因为任何未知函数都可以找到它的近似初等函数(如:包含待定系数的多项式或三角函数),所以从根本上克服了解析法(无法找到初等原函数)的局限性—牺牲极小的理论计算精度,却换回了对大量复杂二维和三维工程问题的适用性。

微分方程的近似解法:差分法泛函变分的近似解法:里兹法,有限元法第4.2节泛函变分原理一、泛函的概念(借助讲解)二、变分的概念借助普通函数微分的概念,用类比法讲解三、泛函的极值条件借助普通函数的极值条件,用类比法讲解四、里兹法(补充内容,但是很重要)泛函变分的近似解法一、泛函的概念通过教材§泛函的概念:函数的函数泛函与普通函数的区别就在于:函数的自变量是数;而泛函的自变量则是函数,泛函的定义域由具有一定条件的一组函数组成。

有限元法和应用总结课件

有限元法和应用总结课件

线弹性有限元
线弹性有限元是以理想弹性体为研究对象旳, 所考虑旳变形建立在小变形假设旳基础上。在 此类问题中,材料旳应力与应变呈线性关系, 满足广义胡克定律;应力与应变也是线性关系, 线弹性问题可归结为求解线性方程问题,所以 只需要较少旳计算时间。假如采用高效旳代数 方程组求解措施,也有利于降低有限元分析旳 时间。
平面单元划分原则
• 1.单元形状:常用单元形状有三角形单元、矩形单元和等 参数单元。他们旳特点是单元旳节点数越多,其计算精 度越高,三角形单元与等参数单元可适应任意边界。
• 2.划分原则: • 1)划分单元旳个数,视计算机要求旳精度和计算机容量
而定,单元分得越多,块越小其精度越高,但需要旳计 算机容量越大,所以,须根据实际情况而定。 • 2)划分单元旳大小,可根据部位不同有所不同,在位 移或应力变化大旳部位取得单元要小;在位移或应力变 化小旳部位取得单元要大,在边界比较平滑旳部位,单 元可大。
移,另一部分基本未知量为节点力。
*8.有限元法分析过程(续)
• 有限元位移法计算过程旳系统性、规律性强,尤 其合适于编程求解。一般除板壳问题旳有限元应 用一定量旳混正当外,其他全部采用有限元位移 法。所以,一般不做尤其申明,有限元法指旳是 有限元位移法。
• 有限元分析旳后处理主要涉及对计算成果旳加工 处理、编辑组织和图形表达三个方面。它能够把 有限元分析得到旳数据,进一步转换为设计人员 直接需要旳信息,如应力分布状态、构造变形状 态等,而且绘成直观旳图形,从而帮助设计人员 迅速旳评价和校核设计方案。
• 虚位移原理是平衡方程和力旳边界条件旳等效积 分旳“弱”形式;
• 虚应力原理是几何方程和位移边界条件旳等效积 分“弱”形式。
3.虚功原理(续)

有限元方法的发展及应用

有限元方法的发展及应用

有限元⽅法的发展及应⽤有限元⽅法的发展及应⽤摘要:有限元法是⼀种⾼效能、常⽤的计算⽅法。

有限元法在早期是以变分原理为基础发展起来的,所以它⼴泛地应⽤于以拉普拉斯⽅程和泊松⽅程所描述的各类物理场中。

⾃从1969年以来,某些学者在流体⼒学中应⽤加权余数法中的迦辽⾦法或最⼩⼆乘法等同样获得了有限元⽅程,因⽽有限元法可应⽤于以任何微分⽅程所描述的各类物理场中,⽽不再要求这类物理场和泛函的极值问题有所联系。

基本思想:由解给定的泊松⽅程化为求解泛函的极值问题。

1有限元法介绍1.1有限元法定义有限元法(FEA,Finite Element Analysis)的基本概念是⽤较简单的问题代替复杂问题后再求解。

它是起源于20世纪50年代末60年代初兴起的应⽤数学、现代⼒学及计算机科学相互渗透、综合利⽤的边缘科学。

有限元法的基本思想是将求解域看成是由许多称为有限元的⼩的互连⼦域组成,对每⼀单元假定⼀个合适的(较简单的)近似解,然后推导求解这个域总的满⾜条件(如结构的平衡条件),从⽽得到问题的解。

这个解不是准确解,⽽是近似解,因为实际问题被较简单的问题所代替。

由于⼤多数实际问题难以得到准确解,⽽有限元不仅计算精度⾼,⽽且能适应各种复杂形状,因⽽成为⾏之有效的⼯程分析⼿段。

有限元法最初应⽤在⼯程科学技术中,⽤于模拟并且解决⼯程⼒学、热学、电磁学等物理问题。

1.2有限元法优缺点有限元⽅法是⽬前解决科学和⼯程问题最有效的数值⽅法,与其它数值⽅法相⽐,它具有适⽤于任意⼏何形状和边界条件、材料和⼏何⾮线性问题、容易编程、成熟的⼤型商⽤软件较多等优点。

(1)概念浅显,容易掌握,可以在不同理论层⾯上建⽴起对有限元法的理解,既可以通过⾮常直观的物理解释来理解,也可以建⽴基于严格的数学理论分析。

(2)有很强的适⽤性,应⽤范围极其⼴泛。

它不仅能成功地处理线性弹性⼒学问题、费均质材料、各向异性材料、⾮线性应⽴-应变关系、⼤变形问题、动⼒学问题已及复杂⾮线性边界条件等问题,⽽且随着其基本理论和⽅法的逐步完善和改进,能成功地⽤来求解如热传导、流体⼒学、电磁场等领域的各类线性、⾮线性问题。

换元法后新方程和原方程

换元法后新方程和原方程

换元法后新方程和原方程【原创实用版】目录1.换元法的概念和作用2.换元法后新方程的生成3.新方程与原方程的关系4.换元法在实际问题中的应用正文一、换元法的概念和作用换元法是数学中一种常用的方法,它通过将一个复杂数学问题转换为另一个相对简单的问题,使得问题得以简化和解决。

换元法可以使问题中的变量和关系更加清晰,从而更容易找到问题的解决方案。

这种方法在微积分、代数和几何等领域都有广泛的应用。

二、换元法后新方程的生成在使用换元法时,我们通常会引入一个新的变量,用以表示原问题中的某个变量或关系的变化。

例如,如果我们有一个关于 x 的方程,我们可以通过引入一个新的变量 y 来表示 x 的平方,从而将原方程转化为一个关于 y 的新方程。

这个新方程通常会更容易求解,因为它表达了原问题中的变量和关系之间的关系。

三、新方程与原方程的关系新方程与原方程之间的关系是紧密相连的。

新方程的解通常可以转化为原方程的解,而原方程的解也可以通过新方程得到。

换元法实际上是将原方程转化为一个新方程,从而将问题转化为一个更容易解决的形式。

通过这种方式,我们可以更容易地找到问题的解决方案。

四、换元法在实际问题中的应用换元法在实际问题中有广泛的应用。

例如,在解决物理问题时,我们通常会使用换元法来简化问题。

通过引入新的变量,我们可以将复杂的物理问题转化为更容易解决的数学问题。

同样,在解决经济问题时,我们也会使用换元法来简化问题。

通过引入新的变量,我们可以更好地理解经济问题的本质,并找到问题的解决方案。

总的来说,换元法是一种非常有用的数学方法,它可以帮助我们简化复杂的问题,并更容易地找到问题的解决方案。

不定积分第二种换元法

不定积分第二种换元法
通过简单的积分问题,如 $int x^2 dx$,展示如何使用第二种换元法进行求解。解析过程中强调变量代 换和积分区间的调整。
复杂实例解析
总结词
复杂实例展示了方法的实际应用
详细描述
选取具有挑战性的不定积分问题,如 $int frac{e^x}{x} dx$,逐步展示如何通过第二种 换元法化简积分,并最终得出答案。
扩展微积分的应用范围
掌握第二种换元法后,学生可以在更广泛的 领域应用微积分知识,解决实际问题。
在其他数学领域的应用
在实变函数中的应用
实变函数是研究实数范围上的函数的数学分 支,第二种换元法在实变函数中也有广泛的 应用。
在复变函数中的应用
复变函数是研究复数范围内函数的数学分支, 其中许多问题可以通过第二种换元法得到解 决。
在第二种换元法中,首先需要选择一个适当的换元函数,通常是为了简化被积函数的形式。然后确定新变量的范 围,将原不定积分中的自变量替换为新变量。接着将被积函数转化为新变量的函数,最后根据新变量的范围计算 不定积分的结果。
04
第二种换元法实例解析
简单实例解析
总结词
简单实例有助于理解基本概念和方法
详细描述
THANKS
感谢观看
03
第二种换元法原理
第二种换元法的定义
总结词
不定积分的第二种换元法是通过引入新的变量来简化不定积分的过程。
详细描述
不定积分的第二种换元法是一种基于变量替换的方法,通过选择适当的换元函 数,将原不定积分转化为更易于计算的形式。
第二种换元法的适用范围
总结词
第二种换元法适用于被积函数难以直接积分的情况,尤其是含有根号或三角函数 的不定积分。
意义
不定积分第二种换元法的意义在于,它提供了一种有 效的工具来解决一些难以处理的不定积分问题。在实 际应用中,许多物理、工程和科学问题都需要解决不 定积分,而第二种换元法可以帮助我们更准确地计算 这些不定积分,从而为解决实际问题提供更可靠的数 学支持。此外,不定积分第二种换元法也是数学理论 体系的重要组成部分,它推动了数学的发展和进步。
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

元法概念意义与应用 TYYGROUP system office room 【TYYUA16H-TYY-TYYYUA8Q8-有限元法概论、意义与应用班级: 2013信息姓名:张正学号指导老师:曾伟梁摘要:有限元法的基础是变分原理和加权余量法,其基本求解思想是把计算域划分为有限个互不重叠的单元,在每个单元内,选择一些合适的节点作为求解函数的插值点,将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式,借助于变分原理或加权余量法,将微分方程离散求解。

采用不同的权函数和插值函数形式,便构成不同的有限元方法。

关键词:有限元法;变分原理;加权余量法;函数。

Abstract:Finite element method is based on the variational principle and the weighted residual method, the basic idea is to solve the computational domain is divided into a finite number of non-overlapping units, each unit, select some appropriate function for solving the interpolation node points as , the differential variables rewritten or its derivative by the variable value of the selected node interpolation functions consisting of linear expressions, by means of variational principle or weighted residual method, the discrete differential equations to solve. Different forms of weight functions and interpolation functions, it constitutes a different finite element method.Keywords:Finite element method; variational principle; weighted residual method; function。

引言随着现代科学技术的发展,人们正在不断建造更为快速的交通工具、更大规模的建筑物、更大跨度的桥梁、更大功率的发电机组和更为精密的机械设备。

这一切都要求工程师在设计阶段就能精确地预测出产品和工程的技术性能,需要对结构的静、动力强度以及温度场、流场、电磁场和渗流等技术参数进行分析计算。

例如分析计算高层建筑和大跨度桥梁在地震时所受到的影响,看看是否会发生破坏性事故;分析计算核反应堆的温度场,确定传热和冷却系统是否合理;分析涡轮机叶片内的流体动力学参数,以提高其运转效率。

这些都可归结为求解物理问题的控制偏微分方程式往往是不可能的。

近年来在计算机技术和数值分析方法支持下发展起来的有限元分析(FEA,Finite Element Analysis)方法则为解决这些复杂的工程分析计算问题提供了有效的途径。

有限元法是一种高效能、常用的计算方法.有限元法在早期是以变分原理为基础发展起来的,所以它广泛地应用于以拉普拉斯方程和泊松方程所描述的各类物理场中(这类场与泛函的极值问题有着紧密的联系)。

自从1969年以来,某些学者在流体力学中应用加权余数法中的迦辽金法(Galerkin)或最小二乘法等同样获得了有限元方程,因而有限元法可应用于以任何微分方程所描述的各类物理场中,而不再要求这类物理场和泛函的极值问题有所联系.有限元方法最早应用于结构力学,后来随着计算机的发展慢慢用于流体力学的数值模拟。

在有限元方法中,把计算域离散剖分为有限个互不重叠且相互连接的单元,在每个单元内选择基函数,用单元基函数的线形组合来逼近单元中的真解,整个计算域上总体的基函数可以看为由每个单元基函数组成的,则整个计算域内的解可以看作是由所有单元上的近似解构成。

在河道数值模拟中,常见的有限元计算方法是由变分法和加权余量法发展而来的里兹法和伽辽金法、最小二乘法等。

根据所采用的权函数和插值函数的不同,有限元方法也分为多种计算格式。

从权函数的选择来说,有配置法、矩量法、最小二乘法和伽辽金法,从计算单元网格的形状来划分,有三角形网格、四边形网格和多边形网格,从插值函数的精度来划分,又分为线性插值函数和高次插值函数等。

不同的组合同样构成不同的有限元计算格式。

对于权函数,伽辽金(Galerkin)法是将权函数取为逼近函数中的基函数;最小二乘法是令权函数等于余量本身,而内积的极小值则为对代求系数的平方误差最小;在配置法中,先在计算域内选取N个配置点。

令近似解在选定的N个配置点上严格满足微分方程,即在配置点上令方程余量为0。

插值函数一般由不同次幂的多项式组成,但也有采用三角函数或指数函数组成的乘积表示,但最常用的多项式插值函数。

有限元插值函数分为两大类,一类只要求插值多项式本身在插值点取已知值,称为拉格朗日(Lagrange)多项式插值;另一种不仅要求插值多项式本身,还要求它的导数值在插值点取已知值,称为哈密特(Hermite)多项式插值。

单元坐标有笛卡尔直角坐标系和无因次自然坐标,有对称和不对称等。

常采用的无因次坐标是一种局部坐标系,它的定义取决于单元的几何形状,一维看作长度比,二维看作面积比,三维看作体积比。

在二维有限元中,三角形单元应用的最早,近来四边形等参元的应用也越来越广。

对于二维三角形和四边形电源单元,常采用的插值函数为有Lagrange插值直角坐标系中的线性插值函数及二阶或更高阶插值函数、面积坐标系中的线性插值函数、二阶或更高阶插值函数等。

一、有限元法的孕育过程及诞生和发展大约在300年前,牛顿和莱布尼茨发明了积分法,证明了该运算具有整体对局部的可加性。

虽然,积分运算与有限元技术对定义域的划分是不同的,前者进行无限划分而后者进行有限划分,但积分运算为实现有限元技术准备好了一个理论基础。

在牛顿之后约一百年,着名数学家高斯提出了加权余值法及线性代数方程组的解法。

这两项成果的前者被用来将微分方程改写为积分表达式,后者被用来求解有限元法所得出的代数方程组。

在18世纪,另一位数学家拉格郎日提出泛函分析。

泛函分析是将偏微分方程改写为积分表达式的另一途经。

在19世纪末及20世纪初,数学家瑞雷和里兹首先提出可对全定义域运用展开函数来表达其上的未知函数。

1915年,数学家伽辽金提出了选择展开函数中形函数的伽辽金法,该方法被广泛地用于有限元。

1943年,数学家库朗德第一次提出了可在定义域内分片地使用展开函数来表达其上的未知函数。

这实际上就是有限元的做法。

所以,到这时为止,实现有限元技术的第二个理论基础也已确立。

“有限单元”,这样的名词。

此后,这样的叫法被大家接受,有限元技术从此正式诞生。

1990年10月美国波音公司开始在计算机上对新型客机B-777进行“无纸设计”,仅用了三年半时间,于1994年4月第一架B-777就试飞成功,这是制造技术史上划时代的成就,其中在结构设计和评判中就大量采用有限元分析这一手段。

在有限元分析的发展初期,由于其基本思想和原理的“简单”和“朴素”,以至于许多学术权威都对其学术价值有所鄙视,国际着名刊物Journal of Applied Mechanics 许多年来都拒绝刊登有关于有限元分析的文章。

然而现在,有限元分析已经成为数值计算的主流,不但国际上存在如ANSYS等数种通用有限元分析软件,而且涉及到有限元分析的杂志也有几十种之多。

有限元方法与其他求解边值问题近似方法的根本区别在于它的近似性仅限于相对小的子域中。

20世纪60年代初首次提出结构力学计算有限元概念的克拉夫(Clough)教授形象地将其描绘为:“有限元法=Rayleigh Ritz法+分片函数”,即有限元法是Rayleigh Ritz法的一种局部化情况。

不同于求解(往往是困难的)满足整个定义域边界条件的允许函数的Rayleigh Ritz法,有限元法将函数定义在简单几何形状(如二维问题中的三角形或任意四边形)的单元域上(分片函数),且不考虑整个定义域的复杂边界条件,这是有限元法优于其他近似方法的原因之一。

有限元方法(FEM)的基础是变分原理和加权余量法,其基本求解思想是把计算域划分为有限个互不重叠的单元,在每个单元内,选择一些合适的节点作为求解函数的插值点,将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式,借助于变分原理或加权余量法,将微分方程离散求解。

采用不同的权函数和插值函数形式,便构成不同的有限元方法。

有限元方法最早应用于结构力学,后来随着计算机的发展慢慢用于流体力学的数值模拟。

在有限元方法中,把计算域离散剖分为有限个互不重叠且相互连接的单元,在每个单元内选择基函数,用单元基函数的线形组合来逼近单元中的真解,整个计算域上总体的基函数可以看为由每个单元基函数组成的,则整个计算域内的解可以看作是由所有单元上的近似解构成。

在河道数值模拟中,常见的有限元计算方法是由变分法和加权余量法发展而来的里兹法和伽辽金法、最小二乘法等。

根据所采用的权函数和插值函数的不同,有限元方法也分为多种计算格式。

从权函数的选择来说,有配置法、矩量法、最小二乘法和伽辽金法,从计算单元网格的形状来划分,有三角形网格、四边形网格和多边形网格,从插值函数的精度来划分,又分为线性插值函数和高次插值函数等。

不同的组合同样构成不同的有限元计算格式。

对于权函数,伽辽金(Galerkin)法是将权函数取为逼近函数中的基函数;最小二乘法是令权函数等于余量本身,而内积的极小值则为对代求系数的平方误差最小;在配置法中,先在计算域内选取N个配置点。

令近似解在选定的N个配置点上严格满足微分方程,即在配置点上令方程余量为0。

插值函数一般由不同次幂的多项式组成,但也有采用三角函数或指数函数组成的乘积表示,但最常用的多项式插值函数。

有限元插值函数分为两大类,一类只要求插值多项式本身在插值点取已知值,称为拉格朗日(Lagrange)多项式插值;另一种不仅要求插值多项式本身,还要求它的导数值在插值点取已知值,称为哈密特(Hermite)多项式插值。

单元坐标有笛卡尔直角坐标系和无因次自然坐标,有对称和不对称等。

常采用的无因次坐标是一种局部坐标系,它的定义取决于单元的几何形状,一维看作长度比,二维看作面积比,三维看作体积比。

相关文档
最新文档