王勖成《有限单元法》学习总结

合集下载

计算力学 有限单元法 清华大学 王勖成

计算力学 有限单元法 清华大学 王勖成

0.5 本课程的内பைடு நூலகம்与要求
1. 计算力学(1)包括: 第1章 有限单元法的理论基础 ——加权余量法和变分原理 第2章 弹性力学问题有限单元法的 一般原理和表达格式
第3章 线性代数方程组的解法 第4章 单元和插值函数的构造 第5章 等参单元和数值积分 第6章 有限单元法应用中的若干实际考虑
2. 学习的形式 1)课堂讲授(48学时+8学时辅导) 2)作业习题(概念和算法的练习) 3)自选论文(程序) 程序实践,要求写出程序报告。
增压风洞的第1阶模态
f = 10.36 Hz
0.3.2 对各种复杂材料细观结构的 分析
Voronoi单元的特点
充分考虑夹杂的微 观结构特点
夹杂随机分布 夹杂大小正态分布 夹杂形状任意
能计算宏观结构变 形 便于分析微观结构 变化对宏观性能的 影响
单元脱层过程的模拟
脱层断裂准则 网格重划分技术
复杂复合材料的损伤模拟
等著名学者著教材
计算力学(1)
—— 有限元方法与数值分析
授课教师:牛莉莎
Chap.0 绪论
计算力学的定义
以现代力学、应用数学为基础,以计算机及 其技术为工具,以求解现代工程和科技中的 力学问题为目标。研究离散化理论和求解方 程的一门应用基础性学科。它伴随计算机的 出现而兴起,现仍在快速发展的学科。
课程设置目标及主要内容
课程设置的目标
使学生掌握扎实的理论基础, 掌握有限单元法的原理、方法和计算机执行; 另一方面培养学生很强的实践动手能力, 包括建模、求解、使用和开发程序的能力。
0.1. 有限元方法 Finite Element Method
0.1.1 有限元方法是一种 关于微分方程的数值解的方法; 微分方程的数值解的方法

有限单元法课后习题全部解答_王勖成

有限单元法课后习题全部解答_王勖成

d 2w dx2
δ
dw dx

d 3w dx3
δ
w
L 0
= 0
∫ 1.5 如有一问题的泛函= 为 Π(w)
L
EI
0 2
d 2w dx2
2
+
kw2 2
+ qwdx ,其中 E,
I,
k 是常数,q
是给定函数,w 是未知函数,试导出原问题的微分方程和边界条件.
∫ = δΠ(w)
L 0
x3 L2
)
+
a2 (x2

x3 L
)
+
x3 L3
(1)
x3 上式中的最后一项 L3 前面没有待定系数,这是由于使用了在 x=L 处φ=1 的强制边界条件。
从物理意义上说,相当于给定边界条件的解为齐次方程的通解加一个特解的缘故。将(1)
式代入教材(1.2.26)式,得到残量:
R(
x)
=
a1 (−6
x L2
+
Q
δφ
dΩ
+
Γ−Γq
k
∂φ ∂n
δφ
dΩ

Γq
αφ

q

k
∂φ ∂n
δφ d
Γ
欧拉方程: k
∂2φ ∂x2
+
k
∂2φ ∂y 2
+
Q
=0
Γφ
自然边界: αφ

q

k
∂φ ∂n
=0
Γ

Γq
强制边界:
k
∂φ ∂n
=0
习题 1.8: 板弯曲问题的平衡方程为:

有限单元法

有限单元法

有限单元法人们常说:“教学有法,教无定法。

”的确,要提高语文课堂教学的质量,要提高学生的语文素养,教师不能一味地把知识灌输给学生,而应该为学生营造轻松、自主、开放的课堂氛围,从而提高学生的学习兴趣。

如何将课堂活动落到实处?语文老师们苦思冥想,找出了许多种教学方法,但这些教学方法都存在一个共同的问题:一节课下来,学生的知识似乎没有增加多少,他们好像只懂得了听讲,对知识点不求甚解,效果可见不佳。

那么怎样才能让学生在有限的时间内既扎实基础又培养能力呢?有限单元法可以助你一臂之力。

这就是有限单元法。

在上《夏天里的成长》这篇课文时,我把全班分成了三组,每一组负责查阅《大自然的语言》《夏天里的成长》和《童年的水墨画》三篇课文。

每个小组安排一名组员负责摘抄三篇课文中具有代表性的段落,并把它们进行分类整理,写出自己的感受。

这一环节引导学生在课外对课文进行深入地了解,发挥了课本学习的延伸作用。

《夏天里的成长》一课中,安排了三次关于“蝉鸣”的交流讨论,我告诉学生“不同的季节会听到不同的蝉声,我们所熟悉的蝉声就来自这个春天……”“请大家拿出各自的工具书,通过字典或百度来了解一下‘蝉’这个字的含义。

”“‘鸣’的古意是什么?”通过交流与探讨,同学们纷纷表示会收集“鸣”的资料,丰富自己的知识。

整个过程轻松愉快,活跃了课堂气氛,培养了学生读书的好习惯。

除了这些,我还用了有限单元法设计了“一石激起千层浪”这一环节,精心创设教学情境,使学生置身于具体的情境之中,受到熏陶,得到启迪。

在交流讨论时,有同学提出“有的蝉是好几年才叫一次的,一辈子就叫一回,也有的蝉在一年中的不同时候都叫……那么蝉为什么叫的次数不同呢?”面对这样的问题,我们没有急于给出答案,而是鼓励学生继续查阅资料,多思考,相信他们肯定会带着这个问题走进下一课。

这一环节的设计巧妙利用了网络资源,拓宽了学生的视野,开阔了学生的思路,学生仿佛一下子解开了心中的疑惑,收获良多。

有限元基础学习心得

有限元基础学习心得

有限元基础学习心得一、问题:1、在开始安装软件时无法正常安装。

2、一些输入符号上的错误,如2.1e11,习惯上输入成了2.1ell,说明对物理意义并不是很清楚。

3、只是按照步骤一步一步往下走,不应该单纯只追求结果,应该要弄懂每一步都是什么意思。

但是现在做完之后根本不知道错在哪一步。

4、老师在课堂上讲过的坝体的载荷分布问题,应该是水深处压力,F应该修改为10000(0.45-X),这样计算的结果会合理一些。

5、英文界面的问题。

6、在操作时要细心,不能丢三落四,尽量独自完成练习,但是可以与同学做学习心得上的交流。

7、操作时不记得要经常保存。

8、对于有限元基本思想的理解不深(为什么要划分网格,ANSYS不是有限元分析的唯一软件)。

9、在生成几何模型时提前划分网格的一处有哪些,局部坐标系的用处有哪些。

二、建议1、希望老师可以推荐几本好的教材,学习起来比较得心应手。

2、希望可以多安排一些上机练习,练习量比较少,进步不大。

(这样理论学习上应该会有很大提高。

)3、上机时指导更加详细一些,一些问题还是有一些难度的。

4、讲课的速度开始时有些快,示范操作时速度慢一些,有一些同学可能会跟不上。

5、上课时多讲解一些操作方面的知识(特别是网格划分和结果显示,以及选择合适的单元类型的方法),增加一些对实际问题的分析和解决实例。

6、希望老师可以将软件及课程中出现的重要单词罗列出来,具体操作步骤的意义可以挑典型例题加以讲解,适当做一些总结。

7、希望老师可以在重要章节可以多重复几遍,加深印象。

8、建议老师安排同学们分组进行一些没有操作步骤提示的问题。

9、上机作业可能会存在抄袭现象。

10、对于用矩阵表达的一些公式的意义多加以讲解。

11、希望可以增加一些弹性力学的讲解。

12、希望老师能在作业每个操作步骤里添加一些解释性的说明。

13、希望可以多讲解一些船舶建模的基本方法以及它与桥梁建模之间的区别。

三、经验\感受:1、建议同学们在遇到问题时最好能记下来,积累经验,避免犯同样的错误。

有限元读书报告

有限元读书报告

有限元读书报告1.有限元的基本理论在当前科学和工程技术的发展和研究中,有限元分析方法是应用最广泛的数值方法,它最早由clough在20世纪60年代提出了“有限单元法”的概念,研究人员们以此为基础不断的探索与创新,经过40年的发展从有限元法的基本概念演化出了一种新的数值分析方法。

有限元分析法把连续体的全求解域看成是由许多个子域组成,对全求解域进行离散,再对各个子域单元上分片假定一个合适的近似解,最后推导全求解域的满足条件建立方程,解出方程即可。

在确定工程和物理问题的数学模型后,使用有限元方法计算模型。

基本思想可以概括如下:1.把连续体的全求解域看成是由许多个子域组成的,并对其进行离散,一个连续体是通过各个单元边界上的节点互连组合成的。

2.在每个元素上假设近似函数,然后使用这些近似函数来表示解域中的未知字段变量。

每个单元中假定的近似函数通常由单元每个节点处的未知场函数值及其相应的插值函数表示。

我们知道,在这些节点上,场函数的值是相同的,所以它们可以作为数值解中的基本未知数。

然后就可以求出原始的场函数无穷多自由度的求解问题转化为场函数节点值的有限自由度的求解问题。

3.基于原问题的数学模型,采用等效加权法或变分原理建立有限元解方程,并用数值方法求解。

从上面所述的有限元法的基本思路中可以得到其具有以下四个特性:1.适应性,表现在其适用于复杂几何模型中;2.可应用性,表现于其在各种物理问题中的使用;3.可靠性,表现为其建立于严格的理论基础上;4.高效性,表现为其特别适合计算机的编程和执行。

由于上述四个特点,有限元法已成为应用最广泛的数值方法。

2.有限元的发展趋势纵观当今国际上cae软件的发展情况,可以看出有限元分析方法的一些发展趋势:2.1与CAD软件的无缝集成当今有限元分析软件的一个发展趋势是与通用cad软件的集成使用,即在用cad软件完成部件和零件的造型设计后,能直接将模型传送到cae软件中进行有限元网格划分并进行分析计算,如果分析的结果不满足设计要求则重新进行设计和分析直到令人满意,这大大提高了设计水平和效率。

(完整word版)有限单元法读书报告

(完整word版)有限单元法读书报告

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

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

关键词:有限单元法;插值函数;网格划分;实例分析1 有限单元法概述1。

1 有限单元法的简介有限单元法[1]是应用局部的近似解来建立整个定义域的解的一种方法。

先把注意力集中在单个单元上,进行上述所谓的单元分析.基本前提是每一单元要尽可能小,以致其边界值在整个边界上的变化也是小的。

这样,边界条件就能取某一在结点间插值的光滑函数来近似,在单元内也容易建立简单的近似解。

因此,比起经典的近似法,有限元法具有明显的优越性。

比如经典的Ritz法,要求选取一个函数来近似描述整个求解区域中的位移,并同时满足边界条件,这是相当困难的.而有限元法采用分块近似,只需对一个单元选择一个近似位移函数,且不必考虑位移边界条件,只须考虑单元之间位移的连续性即可.对于具有复杂几何形状或材料、荷载有突变的实际结构,不仅处理简单,而且合理适宜.1.2 有限单元法的基本方法简介有限单元法,是一种有效解决数学问题的解题方法.在有限元方法中,把计算域离散剖分为有限个互不重叠且相互连接的单元,在每个单元内选择基函数,用单元基函数的线形组合来逼近单元中的真解,整个计算域上总体的基函数可以看为由每个单元基函数组成的,则整个计算域内的解可以看作是由所有单元上的近似解构成。

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

根据所采用的权函数和插值函数的不同,有限元方法也分为多种计算格式.从权函数的选择来说,有配置法、矩量法、最小二乘法和伽辽金法,从计算单元网格的形状来划分,有三角形网格、四边形网格和多边形网格,从插值函数的精度来划分,又分为线性插值函数和高次插值函数等.不同的组合同样构成不同的有限元计算格式。

第三章 有限单元法

第三章 有限单元法

第3章有限单元法在工程技术领域内,工程师常常运用数学和力学的知识将实际问题抽象成它们应遵循的基本方程(常微分方程或偏微分方程)和相应的边界条件。

对于大多数的工程技术问题,由于物体的几何形状和载荷作用方式是很复杂的,除了方程性质比较简单且几何边界相当规则的少数问题之外,试图按经典的弹性力学和塑性力学方法获得解析解是十分困难的,甚至是不可能的。

为了克服这种困难,有两条解决途径:一是引入简化假设,将方程和边界条件简化为能够处理的问题,从而得到它在简化状态下的解答。

这种方法只在有限的情况下可行,因为过多的简化将可能导致不正确的甚至错误的答案。

另一条解决途径就是数值解法,如有限差分法、边界元法、有限单元法和离散元法等。

对于非线性问题,有限单元法更为有效,且已经出现了许多通用程序。

有限单元法的主要优点是:①建立于严格理论基础上的可靠性。

因为用于建立有限元方程的变分原理或加权余量法在数学上已被证明是微分方程和边界条件的等效积分形式。

只要原问题的数学模型是正确的,同时用来求解有限元方程的算法是稳定、可靠的,如果单元满足收敛准则,则近似解最后收敛于原数学模型的精确解;②适应性强,应用范围广,不仅能成功地分析具有复杂边界条件、非线性、非均质材料、动力学等难题,而且还可以推广到解答数学方程中的其它边值问题,如热传导、电磁场、流体力学等问题;③适合计算机实现的高效性。

由于有限元分析的各个步骤可以表达成规范化的矩阵形式,最后导致求解方程可以统一为标准的矩阵代数问题,特别适合计算机的编程和执行。

已经出现了许多大型结构分析通用程序,如:NASTRAN、ASKA、ADINA、ANSYS、ABAQUS等,可以直接应用。

这些优点使有限单元法得到了广泛的应用和发展。

3.1有限单元法分析的基本步骤在工程或物理问题的数学模型(基本变量、基本方程、求解域和边界条件等)确定以后,有限单元法作为对其进行分析的数值计算方法的基本步骤如下:(1) 离散化一个复杂的弹性体可以看成是由无限个质点组成的连续体,它具有无限个自由度。

有限元基础知识归纳

有限元基础知识归纳

有限元知识点归纳1.、有限元解的特点、原因?答:有限元解一般偏小,即位移解下限性原因:单元原是连续体的一部分,具有无限多个自由度。

在假定了单元的位移函数后,自由度限制为只有以节点位移表示的有限自由度,即位移函数对单元的变形进行了约束和限制,使单元的刚度较实际连续体加强了,因此,连续体的整体刚度随之增加,离散后的刚度较实际的刚度K为大,因此求得的位移近似解总体上将小于精确解。

2、形函数收敛准则(写出某种单元的形函数,并讨论收敛性)P49(1)在节点i处N i=1,其它节点N i=0;(2)在单元之间,必须使由其定义的未知量连续;(3)应包含完全一次多项式;(4)应满足∑Ni=1以上条件是使单元满足收敛条件所必须得。

可以推证,由满足以上条件的形函数所建单元是完备协调的单元,所以一定是收敛的。

4、等参元的概念、特点、用时注意什么?(王勖成P131)答:等参元—为了将局部坐标中几何形状规则的单元转换成总体(笛卡尔)坐标中的几何形状扭曲的单元,以满足对一般形状求解域进行离散化的需要,必须建立一个坐标变换。

即:为建立上述的变换,最方便的方法是将上式表示成插值函数的形式,即:其中m是用以进行坐标变换的单元节点数,xi,yi,zi是这些结点在总体(笛卡尔)坐标内的坐标值,Ni’称为形状函数,实际上它也是局部坐标表示的插值函数。

称前者为母单元,后者为子单元。

还可以看到坐标变换关系式和函数插值表示式:在形式上是相同的。

如果坐标变换和函数插值采用相同的结点,并且采用相同的插值函数,即m=n,Ni’=Ni,则称这种变换为等参变换。

5、单元离散?P42答:离散化既是将连续体用假想的线或面分割成有限个部分,各部分之间用有限个点相连。

每个部分称为一个单元,连接点称为结点。

对于平面问题,最简单、最常用的离散方式是将其分解成有限个三角形单元,单元之间在三角形顶点上相连。

这种单元称为常应变三角形单元。

常用的单元离散有三节点三角形单元、六节点三角形单元、四节点四边形单元、八节点四边形单元以及等参元。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

在所有区域内连续可导的并在边界上满足给定位 移条件式的可能位移中,真实位移使系统的总位能取 驻值。即在所有可能位移中,真实位移使系统总位能 取最小值,因此称为最小位能原理。
二、有限元法理论基础-加权余量法和变分原理
2.4 弹性力学的基本方程和变分原理 2.4.4 线弹性力学的变分原理 最小余能原理: 系统的总位能是弹性体余能和外力余能之和:
经分部积分后的等效积分“弱”形式为:
上式第一项代表虚应力在应变上所作的功,第二项代表虚位移约束反 力在给定位移上所作的虚功,称为虚应力原理。
二、有限元法理论基础-加权余量法和变分原理
2.4 弹性力学的基本方程和变分原理 2.4.4 线弹性力学的变分原理 最小位能原理: 系统的总位能是弹性体变形位能和外力位能之和:
几何方程的矩阵形式为:
二、有限元法理论基础-加权余量法和变分原理
2.4 弹性力学的基本方程和变分原理 2.4.1 弹性力学基本方程的矩阵形式: 物理方程——应力与应变关系: 弹性力学中应力-应变之间的转换关系也称弹性关 系。对各向同性的线弹性材料,应力通过应变的表达 式可用矩阵形式表示为: 其中,D为弹性矩阵,它完全取决于弹性体材料的 弹性模量和泊松比。
一、绪论
1.4 有限元法的未来:
① 为了真实地模拟新材料和新结构的行为,需要发 展新的材料本构模型和单元形式; ② 为了分析和模拟各种类型和形式的结构在复杂载 荷工况和环境作用下的全寿命过程响应,需要发展新 的数值分析方案; ③ 有限元软件和CAD/CAM/CAE等软件系统共同集成有限元法理论基础-加权余量法和变分原理
2.4 弹性力学的基本方程和变分原理 2.4.1 弹性力学基本方程的矩阵形式: 几何边界条件: 弹性体在边界上单位面积的内力等于:
边界条件矩阵形式为: 把边界力学方程记为一般形式,则有:
二、有限元法理论基础-加权余量法和变分原理
2.4 弹性力学的基本方程和变分原理 2.4.1 弹性力学基本方程的矩阵形式: 弹性体应变能和余能: 单位体积的应变能(应变能密度)为:
二、有限元法理论基础-加权余量法和变分原理
2.4 弹性力学的基本方程和变分原理 2.4.1 弹性力学基本方程的矩阵形式: 力的边界条件: 弹性体在边界上单位面积的内力等于在边界上已 知弹性体单位面积上作用的面积力,即: 设边界外法线的方向余弦为 体的内力为: ,则边界上弹性
边界条件矩阵形式为:
王勖成《有限单元法》
(学习总结)
汇报人:XXX 时 间:XXX
2017/4/12 1
内容提纲
一、绪论 二、有限元法的理论基础-加权余量法和变分原理 三、弹性力学问题有限元方法的一般原理和表达式 四、单元和插值函数的构造 五、等参元与数值积分 六、有限元法运用中的若干实际考虑
七、线性代数方程组的解法
二、有限元法理论基础-加权余量法和变分原理
2.4 弹性力学的基本方程和变分原理 2.4.1 弹性力学基本方程的矩阵形式: 平衡方程: 弹性体V域内任一点沿坐标轴x,y,z方向的平衡方程 为:
平衡方程矩阵形式:
二、有限元法理论基础-加权余量法和变分原理
2.4 弹性力学的基本方程和变分原理 2.4.1 弹性力学基本方程的矩阵形式: 几何方程——应变与位移关系: 在微小位移和微小变形情况下,略去位移导数的 高次幂,则应变向量和位移向量间的关系为:
八、有限元分析计算机程序
一、绪论
1.1 有限元法要点: ① 将一个表示结构或连续体的求解域离散为若干个 子域(单元),并通过他们边界上的节点相互联接成 为组合体; ② 用每个单元内在所假设的近似函数来分片地表示 全求解域内待求的未知场变量。而每个单元内的近似 函数由未知场函数(或及其导数)在单元各结点上的 数值和与其对应的插值函数来表达(此表达式为矩阵 形式); ③ 通过和原问题数学模型(基本方程、边界条件) 等效的变分原理或加权余量法,建立求解基本未知量 (场函数结点值)的代数方程组或者场微分方程组。
应变能是个正定函数,只有当弹性体内所有的点 都没有应变时,应变能才为零。 单位体积的余能(余能密度)为: 余能也是正定函数,在线弹性力学中弹性体的应 变能等于余能。
二、有限元法理论基础-加权余量法和变分原理
2.4 弹性力学的基本方程和变分原理 2.4.2 弹性力学基本方程的张量形式: 平衡方程: 张量形式的平衡方程为:
一、绪论
1.2 有限元法特性:
① 对于复杂几何构型的适应性(单元在空间可以是一维、二 维或三维的,而每一种单元可以有不同形状); ② 对各种物理问题的可应用性(用单元内近似函数分片地表示 全求解域的未知场函数,并为限制场函数所满足的方程形式, 也为限制各个单元所对应方程必须是相同的形式); ③ 建立于严格理论基础上的可靠性(用于建立有限元方程的变 分原理或加权余量法在数学上以证明是微分方程和边界条件的 等效积分形式); ④ 适合计算机实现的高校性(有限元分析的各个步骤可以表达 成规范的矩阵形式,最后导致求解方程可以统一为标准的矩阵 代数问题,特别适合计算机编程和执行)。
其扩展形式为:
二、有限元法理论基础-加权余量法和变分原理
2.4 弹性力学的基本方程和变分原理 2.4.2 弹性力学基本方程的张量形式: 几何方程: 张量形式的几何方程为:
其扩展形式为:
二、有限元法理论基础-加权余量法和变分原理
2.4 弹性力学的基本方程和变分原理 2.4.2 弹性力学基本方程的张量形式: 物理方程: 张量形式的物理方程为:
在所有在弹性体内满足平衡方程,在边界上满足 力的边界条件的可能盈利中,真实应力使系统的总余 能取驻值,真实位移使系统总位能取最小值类同步骤, 证明在所有应力中,真实应力使系统总余能取最小值, 因此称为最小余能原理。
二、有限元法理论基础-加权余量法和变分原理
2.4 弹性力学的基本方程和变分原理 2.4.4 线弹性力学的变分原理 弹性力学变分原理的能量上、下界: 根据能量平衡,应变能应等于外力功,因此得到 弹性系统的总位能与总余能之和为零:
二、有限元法理论基础-加权余量法和变分原理
2.3 变分原理和里兹方法 2.3.2 里兹方法: 设未知函数的近似解由一族带有待定参数的试探 函数来近似表示,即: 泛函的变分为零相当于将泛函对所包含的待定参 数进行全微分,并令所得方程等于零,即: 由于 是任意的,满足上式时必然有 都等于零。这是与待定系数a的个数相等的方程组, 用以求近似解的经典方法叫做里兹法。
二、有限元法理论基础-加权余量法和变分原理
2.3 变分原理和里兹方法 2.3.1线性、自伴随微分方程变分原理的建立: 泛函的极值性: 对于2m阶微分方程,含0~m-1阶导数的边界条件 称为强制边界条件,近似函数应事先满足。含m~2m1 阶导数的边界条件成为自然边界条件,近似函数不 必事先满足。 设近似场函数 ,则 其中, 是真正的泛函, 是等效积分伽辽金提法 的弱形式,应有:
二、有限元法理论基础-加权余量法和变分原理
2.2 微分方程的等效积分形式和加权余量法 2.2.3 基于等效积分形式的近似方法—加权余量法: 采用使余量的加权积分为零来求微分方程近似解 的方法称为加权余量法。 根据对权函数W的不同选择可得到不同的加权余 量计算方法,常用的方法有: ① 配点法: ② 子域法:在n个子域内W=I,在子域意外W=0。即 强迫余量在n个子域的积分为零。 ③ 最小二乘法: 使 最 小,即 。 ④ 力矩法:
里兹法的实质是从一族假定解中寻求满足泛函变分的最好 解,显然近似解的精度与试探函数的选择有关。
二、有限元法理论基础-加权余量法和变分原理
2.3 变分原理和里兹方法 2.3.2 里兹方法: 当 n 趋于无穷时,近似解收敛于微分方程精确解 的条件如下: ① 试探函数应取自完备函数系列。满足此要求的试 探函数称为是完备的; ② 试探函数应满足 连续性要求,即表示泛函的场 函数最高的微分阶数是m时,试探函数0~m-1阶导 数应是连续的,以保证泛函中的积分存在,满足 此要求的试探函数称为是协调的。
二、有限元法理论基础-加权余量法和变分原理
1.3 变分原理和里兹方法 2.3.1线性、自伴随微分方程变分原理的建立: 泛函的构造: 原问题的微分方程和边界条件的等效积分的伽辽 金法等效于它的变分原理,即原问题的微分方程和边 界条件等效于泛函的变分等于零,亦即泛函取驻值。 反之,如果泛函取驻值则等效于满足问题的微分方程 和边界条件。而泛函可以通过原问题的等效积分的伽 辽金提法得到的,并称这样的到的变分原理为自然变 分原理。
其扩展形式为:
二、有限元法理论基础-加权余量法和变分原理
2.4 弹性力学的基本方程和变分原理 2.4.2 弹性力学基本方程的张量形式: 力的边界条件: 张量形式的力的边界条件为:
其扩展形式为:
二、有限元法理论基础-加权余量法和变分原理
2.4 弹性力学的基本方程和变分原理 2.4.3 平衡方程和几何方程的等效积分“弱”形式 虚位移原理: 平衡方程和力边界条件为: 其等效积分为:
二、有限元法理论基础-加权余量法和变分原理
2.3 变分原理和里兹方法 2.3.1线性、自伴随微分方程变分原理的建立: 线性、自伴随微分算子: 若微分方程: L为微分算子,若 ,则为线性。 L(u)与任意函数的内积: 若 ,则算子为自伴随的。 泛函的构造: 原问题微分方程和边界条件: 与上式等效的伽辽金法: 若: ,其中: 为 原问题的泛函。
二、有限元法理论基础-加权余量法和变分原理
2.1 微分方程的等效积分形式和加权余量法 2.1.3 基于等效积分形式的近似方法—加权余量法: 假设未知函数u可以采用近似函数表示,近似函 数是一族有待定参数的已知函数,一般形式为:
通常n取有限项的近似解不能精确满足微分方程 式和边界条件,故产生残差R,即: 把等效积分形式写成余量形式:
一、绪论
1.3 有限元法的发展和现状:
相关文档
最新文档