有限元分析及应用第四章
有限元方法及应用
![有限元方法及应用](https://img.taocdn.com/s3/m/171296bf1a37f111f1855b81.png)
第四章有限元法及应用FEM(The Finite Element Method) and Application主要内容有限元法的基本概念 有限元法的分类有限元法原理有限元法的求解步骤 常用有限元软件简介第一节概述第一节概述一、为什么需要有限元?)随着生产的发展,不断要求设计高质量、高水平的大型、复杂和精密的机械和工程结构。
在实践中人们也逐渐认识到要达到正确的、高水平的设计,就必须预先通过有效的计算手段确切的了解即将诞生的机械和工程结构在未来工作时的应力、应变及位移等情况,从多种可能的方案中去选择合乎要求的方案。
)但是,传统的一些方法往往难于完成对工程实际问题的有效分析。
弹性力学的经典理论对于几何上复杂、不规则边界、有裂缝或厚度突变以及几何非线性、材料非线性等问题往往解决很困难;优化设计、可靠性设计等也难或根本无法解决。
蓄水后大坝的位移与应变情况、地震时大坝的位移与应变情况等三峡大坝的受力情况温度场分布航天飞机飞行中的受热分析导弹、飞机飞行的流体动力学分析流场分布磁场分布分析卫星、飞船在轨运行时磁场的影响传统方法在处理这类问题时,往往要对一个实际的物理系统作出多种假设,比如形状假设、连续性假设、物体的各项同性假设,然后通过经典理论方法得出问题的解析解,这种解析解从形式上看,可以得出关于实际问题的连续解,比如用方程描述三峡大坝某一点的位移和应变,但这样的解析解往往和实际情况有比较大的偏差。
这对于精度要求不高的领域是可以的,但对于有些领域,就不能满足实际的需要了。
同时,实际中常常要遇到一些几何上复杂、不规则边界、有裂缝或厚度突变以及几何非线性、材料非线性的物理系统,对这些系统经典理论解决起来相当困难,有时甚至无法解决,也就是无法求得解析解。
因此,寻求离散数值分析法就成了必由之路。
常用的数值分析法有两种:差分法和有限元法。
差分法是在传统方法的基础上,将传统方法建立的微分方程中的微分dx、dy、dz变成差分Δx,Δy,Δz,从而把微分方程变成代数方程,用一步步迭代的方法,逐步求出物理系统中各个离散点的物理量,用差分离散解代替连续解。
有限元分析及应用课件
![有限元分析及应用课件](https://img.taocdn.com/s3/m/caeb313230b765ce0508763231126edb6e1a7617.png)
设置材料属性、单元类型等参数。
求解过程
刚度矩阵组装
根据每个小单元的刚度,组装成全局的刚度矩阵。
载荷向量构建
根据每个节点的外载荷,构建全局的载荷向量。
求解线性方程组
使用求解器(如雅可比法、高斯消元法等)求解线性方程组,得到节点的位移。
后处理
01
结果输出
将计算结果以图形、表格等形式输 出,便于观察和分析。
有限元分析广泛应用于工程领域,如结构力学、流体动力学、电磁场等领域,用于预测和优化结构的 性能。
有限元分析的基本原理
离散化
将连续的求解域离散化为有限 个小的单元,每个单元具有特
定的形状和属性。
数学建模
根据物理问题的性质,建立每 个单元的数学模型,包括节点 力和位移的关系、能量平衡等。
求解方程
通过建立和求解线性或非线性 方程组,得到每个节点的位移 和应力分布。
PART 05
有限元分析的工程应用实 例
桥梁结构分析
总结词
桥梁结构分析是有限元分析的重要应用之一,通过模拟桥梁在不同载荷下的响应,评估 其安全性和稳定性。
详细描述
桥梁结构分析主要关注桥梁在不同载荷(如车辆、风、地震等)下的应力、应变和位移 分布。通过有限元模型,工程师可以预测桥梁在不同工况下的行为,从而优化设计或进
刚性问题
刚性问题是有限元分析中的一种 特殊问题,主要表现在模型中某 些部分刚度过大,导致分析结果 失真
刚性问题通常出现在大变形或冲 击等动态分析中,由于模型中某 些部分刚度过高,导致变形量被 忽略或被放大。这可能导致分析 结果与实际情况严重不符。
解决方案:为避免刚性问题,可 以采用多种方法进行优化,如采 用更合适的材料模型、调整模型 中的参数设置、采用更精细的网 格等。同时,可以采用多种方法 对分析结果进行验证和校核,以 确保其准确性。
有限元原理与应用
![有限元原理与应用](https://img.taocdn.com/s3/m/a11c77d60740be1e650e9af7.png)
第二节 平面刚架有限元法
二、单元分析
第二节 平面刚架有限元法
二、单元分析
第二节 平面刚架有限元法
二、单元分析
第二节 平面刚架有限元法
二、单元分析
第二节 平面刚架有限元法
二、单元分析
第二节 平面刚架有限元法
三、坐标变换
第二节 平面刚架有限元法
三、坐标变换
第二节 平面刚架有限元法
三、坐标变换
四 载荷移置
第二节 平面问题有限元法
四 载荷移置
第二节 平面问题有限元法
五 约束处理
第二节 平面问题有限元法
五 约束处理
第二节 平面问题有限元法
五 约束处理
第二节 平面问题有限元法
五 约束处理
第二节 平面问题有限元法
六 求解线方程组
七 计算其它物理量
第二节 平面问题有限元法
八 计算结果处理
第二节 轴对称问题有限元法
二、单元分析
第二节 轴对称问题有限元法
第二节 轴对称问题有限元法
第二节 轴对称问题有限元法
第二节 轴对称问题有限元法
第二节 轴对称问题有限元法
三、单元刚度矩阵
第二节 轴对称问题有限元法
三、单元刚度矩阵
第二节 轴对称问题有限元法
三、单元刚度矩阵
第二节 轴对称问题有限元法
第二节 平面问题有限元法
3 总刚矩阵的特点
第二节 平面问题有限元法
3 总刚矩阵的特点
第二节 平面问题有限元法
四 载荷移置
第二节 平面问题有限元法
四 载荷移置
第二节 平面问题有限元法
四 载荷移置
第二节 平面问题有限元法
四 载荷移置
第二节 平面问题有限元法
第四章有限元单元法
![第四章有限元单元法](https://img.taocdn.com/s3/m/39656e5232687e21af45b307e87101f69e31fb11.png)
¾ ADINA
★ FLUENT
¾ ANSYS
★ SAP
(3) 有限单元法的未来
应用需求:技术革新、设计理论、 制造方法
基础产业:汽车、船舶、冶金、飞机 高新产业:航天、微电系统、纳米器
件:
(3) 有限单元法的未来 待发展的方面
1) 新的材料本构模型和单元型式
2) 结构在复杂环境条件下的全寿命过程响应分析
人类认识自然的得力助手
力学或工程领域求解问题的两大法宝:解 析法和离散法
4.1.1 有限单元法的基本概念
解析法:
它从研究连续体中无限小的微分体入手,得出 描述连续体性质的微分方程。然后根据边界条 件、初始条件可解得一个通解。这个解可给出连 续体内任一点上所求参数的值。 核心是微分方程。 微分方程的建立过程是近似的,而微分方程的 求解过程是精确的。
E
• 长度分别为l1、l2 • 桁架的铰链处受到外力 X1、Y1、X2、Y2、X3、Y3
• 在1点和3点固定铰支 求解内力
铰接桁架
求解过程:
(1). 将结构划分成典型单元的集合——离散化(节点力、节点位移)
--求解过程:
(2).分析每个单元上节点力和节点位移之间的关系 ――单元特性分析
轧机牌坊(三维实体问题,弹性)
稳静态结构问题实例
出钢机部件分析(三维壳体问题,弹性)
稳静态结构问题实例
万向接轴叉头(三维实体问题,弹性)
稳静态结构问题实例
轧钢机刚度(三维实体问题,多体接触)
稳静态结构问题实例
沧州铁狮子(三维不规则实体,弹性)
稳静态结构问题实例
轧制过程仿真(三维实体,弹塑性、接触)
.500E+09
10000
ANSYS有限元分析入门与应用指南
![ANSYS有限元分析入门与应用指南](https://img.taocdn.com/s3/m/4d67afb2f9c75fbfc77da26925c52cc58bd69089.png)
ANSYS有限元分析入门与应用指南第一章:ANSYS有限元分析概述ANSYS是一种常用于工程领域的有限元分析软件,主要用于对各种结构进行力学分析、流体动力学分析、热传导分析等。
本章将对ANSYS的基本原理、工作流程和应用领域进行介绍。
1.1 ANSYS的基本原理ANSYS基于有限元方法,将实际结构或系统离散为有限数量的单元,通过对单元进行各种物理特性的分析,最终得到整个结构的行为。
有限元方法是一种数值分析方法,可以有效解决传统方法难以处理的复杂问题。
1.2 ANSYS的工作流程ANSYS的工作流程包括几个关键步骤:前处理、求解和后处理。
前处理阶段主要负责模型的建立和单元网格的划分,求解阶段进行物理场的计算和求解,后处理阶段对结果进行可视化和分析。
1.3 ANSYS的应用领域ANSYS可应用于各个工程领域,如固体力学、流体力学、热传导、电磁场等。
在航空航天、汽车工程、建筑结构、电子设备等领域都有广泛的应用。
第二章:ANSYS建模与前处理在使用ANSYS进行有限元分析之前,需要对模型进行建模和前处理工作。
本章将介绍ANSYS建模的基本方法和前处理的必要步骤。
2.1 模型建立ANSYS提供了多种建模方法,包括几何建模、CAD导入、脚本编程等。
用户可以根据需要选择合适的建模方法,对模型进行几何设定。
2.2 材料定义和属性设置在进行有限元分析之前,需要为材料定义材料性质和属性。
ANSYS提供了多种材料模型,用户可以根据具体需求进行选择和设置。
2.3 网格划分网格划分是有限元分析中非常重要的一步,它决定了模型的离散精度和计算效果。
ANSYS提供了多种单元类型和划分算法,用户可以根据需要进行合理的网格划分。
第三章:ANSYS求解与后处理在进行前处理完成后,就可以进行有限元分析的求解和后处理了。
本章将介绍ANSYS的求解方法和后处理功能。
3.1 求解方法ANSYS提供了多种求解方法,如直接法、迭代法等。
根据模型的复杂程度和求解要求,用户可以选择合适的方法进行求解。
第四讲-机械结构的有限元分析
![第四讲-机械结构的有限元分析](https://img.taocdn.com/s3/m/91c38863caaedd3383c4d343.png)
5.线性静力分析 线性静力分析注意问题 线性静力分析
NORTHEASTERN UNIVERSITY
单元类型必须指定为线性结构单元类型; 材料属性是线性的,可以是各向同性或正交各向 异性、常量或与温度相关的量; 必须定义弹性模量和泊松比; 对于象重力一样的惯性载荷,必须要定义能计算 出质量的参数,如密度等; 对热载荷,必须要定义热膨胀系数; 对应力、应变感兴趣的区域,网格划分比对位移 感兴趣的区域要密。
10 材料非线性
材料的应力一应变关系是非线性的,影响材料 这种非线性关系的因素可能有:加载历史(如弹 一塑性响应)、环境状况(如温度)、加载的时间 总量(如蠕变)等。 ( ) ANSYS的非线性分析能力主要包括:弹塑性分 析、超弹性分析、蠕变分析等。
弹性力学及有限元
11状态非线性
NORTHEASTERN UNIVERSITY
弹性力学及有限元
3. 结构静力分析定义
NORTHEASTERN UNIVERSITY
计算在固定不变的载荷作用下结构的位移、应力、 应变和力,它不考虑变化的惯性力和阻尼的影响。 计算那些固定不变的惯性载荷对结构的影响,以及 那些可以近似为等价静力作用的随时间变化载荷。 静力分析既可以是线性的也可以是非线性的。
Ansys中的常用杆元
NORTHEASTERN UNIVERSITY
三维空间杆元Link8: : 三维空间杆元
弹性力学及有限元
3.平面梁元
NORTHEASTERN UNIVERSITY
平面梁元:为二维空间单元,每个节点有三个自 由度:一个轴向自由度,一个横向自由度(挠度)和 一个旋转自由度(转角),适用于受弯曲、拉压组合 作用的刚架结构。 这类单元适用于平面刚架结构,即在刚架结构中, 所有构件横截面的主惯性轴之一与刚架所受载荷 在同一平面内,如机床的主轴、导轨等 。
有限元分析及应用
![有限元分析及应用](https://img.taocdn.com/s3/m/36709102866fb84ae45c8d1f.png)
方程
变量
求解 方法
积分法
积分法
解析法、半 解析法、半 解析法 解析法
4
2015-06-13
基本假设
材料力学
1. 连续性假设 2. 均匀性假设 3 各向同性假设 3. 4. 小变形条件
1. 2 任意变形体力学分析的基本变量及方程 研究对象:任意形状的变形体 弹性力学
1. 假定物体是连续的 2. 假定物体是完全弹性的 3. 假定物体是均匀的 4. 假定物体是各向同性的 5. 假定位移和形变是微小 的 满足前4条为理想弹性体
8
2015-06-13
7
2015-06-13
1.3 有限元方法的思路及发展过程 思路: 以计算机为工具,分析任意变形体以获得所有 力学信息,并使得该方法能够普及、简单、高效、 方便,一般人员可以使用。 发展过程:
如: 屋顶结构
杆
支承
真实结构 (物理模型) 铰接
实现过程
连续体
离散体
力学模型
两个重要工具: (1)在理论推导方面采用矩阵方法 (2)在计算中采用了计算机
x yx bx 0 x y xy y b y 0 y x
并且所有力学变量都是x, y的函数,不随 的函数 不随z变化,即有 变化 即有
但
zx
0,
z
E
zy
0
( x v y v u x y
4) 边界条件
25 2015-06-13
l x m yx p x l xy m y p y
u u 和 v v
26
2015-06-13
有限单元法原理及应用简明教程ppt课件
![有限单元法原理及应用简明教程ppt课件](https://img.taocdn.com/s3/m/172af1a480c758f5f61fb7360b4c2e3f5627251f.png)
(a) 瞬变结构
(b) 分离体分析
(c) 平衡状态分析
图2-32 瞬变结构
24
第二章 结构几何构造分析
(2) 两刚片规则 两刚片用三根既不完全平行也不交于同一点的链杆 相联,所得结构是几何不变结构。
(a) 铰与链杆连接两刚片 (b) 三链杆连接两刚片 图2-33 两刚片连接规则
25
第二章 结构几何构造分析
章
生刚体位移时,称之为几何不变结构或几何稳定结构,
节
反之则称为几何可变结构或几何不稳定结构。几何可
目 录
变结构不能承受和传递载荷。对结构进行几何构造分
析也是能够对工程结构作有限单元法分析的必要条件。
11
第二章 结构几何构造分析
(a) 结构本身可变 (b) 缺少必要的约束条件 (c) 约束汇交于一点 图2-1 几何可变结构
节
何不变结构上,由增加二元体而发展的结构,是一个
目
几何不变结构。铰接三角形是最简单的几何不变结构。
录
图2-31 铰接三角形
23
第二章 结构几何构造分析
结构的特征是:当它受载荷作用时会产生微小的 位移, 但位移一旦发生后, 即转变成一几何不变结 构,但结构的内力可能为无限大值或不定值,这样的 结构称为瞬变结构。显然,瞬变结构在工程结构设计 中应尽量避免。
(5) 约束处理,求解系统方程
(6) 其它参数计算
4
第一章 概述
图1-2 工程问题有限单元法分析流程
5
第一章 概述
1.3 工程实例
返 回 章 节 目 录
(a) 铲运机举升工况测试
(b) 铲运机工作装置插入工况有限元分析
图1-3 WJD-1.5型电动铲运机
有限元法_精品文档
![有限元法_精品文档](https://img.taocdn.com/s3/m/7a97ff9f51e2524de518964bcf84b9d528ea2c9f.png)
12
一、有限元法的基本概念
1.什么是有限元法
我们实际要处理的对象都是连续体,在传统设 计思维和方法中,是通过一些理想化的假定后,建 立一组偏微分方程及其相应的边界条件,从而求出 在连续体上任一点上未知量的值。
25
4)具有灵活性和适用性,适应性强(它可以把形状 不同、性质不同的单元组集起来求解,故特别适 用于求解由不同构件组合的结构,应用范围极为 广泛。它不仅能成功地处理如应力分析中的非均 匀材料、各向异性材料、非线性应力应变以及复 杂的边界条件等问题,且随着其理论基础和方法 的逐步完善,还能成功地用来求解如热传导、流 体力学及电磁场领域的许多问题)
21
对于一个具体的工程结构,单元的划分越小, 求解的结果就越精确,同时,其计算工作量也就越 大。
梯子的有限元模型不到100个方程; 在ANSYS分析中,一个小的有限元模型可能有几 千个未知量,涉及到的单元刚度系数几百万个。 单元划分的精细程度,取决于工程实际对计算 结果精确性的要求。
22
4)有限元分析 有限元分析就是利用数学近似的方法对真实
5)在具体推导运算过程中,广泛采用了矩阵方法。
26
2.有限元法的作用 1)减少模型试验的数量(计算机模拟允许对大量
的假设情况进行快速而有效的试验); 2)模拟不适合在原型上试验的设计(例如:器官
移植、人造膝盖); 3)节省费用,降低设计与制造、开发的成本; 4)节省时间,缩短产品开发时间和周期; 5)创造出高可靠性、高品质的产品。
因为点是无限多的,存在无限自由度的问题, 很难直接求解这种偏微分方程用来解决实际工程问 题,因此需要采用近似方法来处理。
杆梁结构有限元分析(第四章)汇总.
![杆梁结构有限元分析(第四章)汇总.](https://img.taocdn.com/s3/m/ffbca306a300a6c30c229fea.png)
可以验证:它满足位移边界条件。这是一个待定函数,也称为试函 数,所谓该函数是待定的,就是因为它中间有一个待定系数,这就需要 通过一个原理来确认它,下面由虚功原理来进行确认。基于式(4-13)的试 函数,则它的应变、虚位移以及虚应变为
( x) c u ( x) c.x ( x) c
x 0
E
du dx
d x 0 dx
x
x
力边界条件BC(p)
( x ) | x l
F px A
4.2 杆件有限元分析的标准化标准与算例
从求解思路来说,可以有两类方法来对该问题进行求解,即: 直接求解法:可以由3个方程来直接求解3个变量;
间接法(试函数):选取变量(位移)作为基本的待求变量,将其它变量都
E c c去δc后,有解
1D问题的最小势能原理求解 先介绍最小势能原理的基本表达式。设有满足位移边界条件BC(u)的许 可位移场,计算该系统的势能为
(u ) U W
其中U为应变能,W为外力功,对于如图4-2所示的算例,有
U (u ( x)) (u ( x))d
面讨论该问题的力学描述与求解。
图4-2 一端固定的拉杆
4.2 杆件有限元分析的标准化标准与算例
基本变量: 由于该问题是沿x方向的一维问题,因此只有沿x方向的基本变量, 即定义沿x方向移动为位移: 定义:沿x方向移动为位移:
沿x方向的相对伸长(或缩短)量为应变:
沿x方向的单位横截面上的受力为应力: 基本方程:
1 2 x x
W P u( x l )
4.2 杆件有限元分析的标准化标准与算例
对于包含有待定系数的试函数而言,真实的位移函数应使得该系统的 势能取极小值,即
《有限元分析及应用》PPT课件
![《有限元分析及应用》PPT课件](https://img.taocdn.com/s3/m/ea89dd4dcd1755270722192e453610661ed95aef.png)
有限元法的孕育过程及诞生和发展
牛顿(Newton)
莱布尼茨(Leibniz G. W.)
6
大约在300年前,牛顿和莱布尼茨发明了积 分法,证明了该运算具有整体对局部的可加 性。虽然,积分运算与有限元技术对定义域 的划分是不同的,前者进行无限划分而后者 进行有限划分,但积分运算为实现有限元技 术准备好了一个理论基础。
u y
dy
vB
v
v y
dy
66
在小变形的前提下,∠A’P’A1很小,可以认 为,线段PA位移后的绝对伸长,可以用线段两 端点沿x轴的位移之差来表示,即:。
PA PA
uA
uP
u
u x
dx u
u x
dx
从而线段PA的正应变
x为:。 x
PA PA PA
u dx x
dx
u x
v
dy
同理线段PB的正应变
y
dy
zy
1 2
zy
z
dz
0
略去微量项,得 yz zy
MY 0 zx xz
MZ 0
xy yx
剪切力互等定律
53
二维问题:平衡微分方程
x yx X 0
x y xy y Y 0 x y
剪切力互等定律
xy yx
54
应力边界条件
四面微分体Mabc
55
效的力系所代替,只能产生局部应力的改变,而在离
这一面积稍远处,其影响可以忽略不计。
60
61
62
均匀分布载荷作用 下的平板,应力分 布是均匀的。
材料力学中的拉伸 应力计算公式就是 圣维南原理应用的 结论。
63
一对集中力F/2作 用点区域仍然有比 较大的应力梯度变 化,但是比等效力
《弹性力学与有限元》第4章梁的有限元分析
![《弹性力学与有限元》第4章梁的有限元分析](https://img.taocdn.com/s3/m/f8721f0b844769eae009edac.png)
j
l
ix
(4-2)
利用材料力学梁的理论很容易求得节点力与位移的关系:
Fi
=
du dx
|x=0 =
− EAα 2
=
EA ui
− l
uj
(4-3)
Fj
=
du dx
|x=l =
EAα 2
=
u EA j
− l
u i
{} {} 将结果写成矩阵形式为:
P
e
=
⎡⎣K
⎤e ⎦
a
e ,其中
{ } { } P
e
=
⎧⎪ ⎨
( ) v = v 0 = α
i
1
( ) θ = θ 0 = − dv | = −α
i
dx x =0
2
() v = v l = α + α l + α l 2 + α l 3
j
1
2
3
4
() θ = θ l = − dv | = −α − 2α l − 3α l 2
j
dx x =l
2
3
4
由(4-6)式可解得以下常数:
4. 2 杆和杆系
当结构长度尺寸比两个截面方向的尺寸大很多的时候,这类结构称为杆件。在杆件结构 中,垂直与长度方向的截面称为横截面,横截面中心的连线称为轴线。如果杆的轴线示直线, 则称为直杆;如果轴线是曲线,则称为曲杆;如果各个截面的尺寸和形状不变,称为等截面 杆,反之称为变截面杆。
杆件的结构可以范围桁杆和梁两类。和其它结构采用铰链连接的杆称为桁杆,桁杆的连 接处可以自由转动,因此这类结构只能承受拉压作用,内部应力为拉应力。影响应力的几何 因素主要是截面面积,与截面形状无关;和其它结构采用固定连接的杆称为梁,梁的连接处 不能自由转动,因此梁不仅能承受拉压,而且能承受完全和扭转作用。这类杆件的内部应力 的状况比较复杂,应力大小和分布不仅与截面有关,而且与截面的形状和方位有很大关系, 建立有限元模型时,这两类杆件可以用相应的杆单元和梁单元离散。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
则称ϕ1、ϕ2Lϕ n 线性相关;
(ii) 若 c1ϕ1 + c2ϕ 2 + L + cnϕ n ≡ 0
仅当
c1
才成立,则称
ϕ=1c、2
=L= ϕ2Lϕ
cn
n
≡0
线性无关。
(2) 线性空间的维数
若线性空间E满足
(i)任意 n+1 个元素一定线性相关。
(ii)存在着 n 个线性无关的元素。
则称线性空间E的维数为 n。
a ⋅ b = a ⋅ b ⋅ cosα ≤ a ⋅ b
上式为 Euclid 空间的三角不等式,此式仅是 Schwarz 不等式的一个特例。 5、收敛性与完备性 (1)收敛性
∀ 点列{xn } ∈E(赋范线性空间),若存在
lim xn − x0 = 0
n →∞
则,x0 称为点列{xn }的强极限,读作:{xn }强收敛于 x0 ,注意模的定义不同收敛的涵
c1ϕ1 + c2ϕ 2
c1ϕ1′ + c2ϕ 2′
第 1 页 共 17 页
有限元分析与应用
霍战鹏
也在(a, b)上连续。所有函数本身及一阶导数都在(a, b)上连续的函数组成一种线性空
间,记作 C1[a, b]。 例4 Rn n 维欧氏空间是线性空间,R2(二维平面), R3(三维空间)是 n 维欧氏空
形的项点为结点,以结点处的函数值对单元内的位移场进行分片线性插值。根据第 3-4 节的
分析可知,对于这样定义的函数 u(x,y)在Ω上连续,且积分
y
∫∫ ∫∫ ∫∫ Ω
u 2dxdy
、
Ω
∂u ∂x
2 dxdy
、
Ω
∂u ∂y
2 dxdy
存在。
在单元边界和结点处,通常意义下 ∂u ∂x 、∂u ∂y 不存在。
义导数的空间就是完备的线性空间,由于在通过结点时 u′ 不连续,有限跳跃量, u′′ 在结点
处为δ—函数。δ—函数本身可积,但平方不可积。我们可以得到这样的结论:对于分段线
性插值的函数 u(x)而言:u∈H1(0, L),但 u∉H2(0, L)。 (2)设Ω为一平面二维区域,将Ω分成若干三角形的片(单元)(图 4-3)。取各三角
例1 在平面R2内,向量 x(x1, x2)可以有下列三种模的定义
x= 2
x12
+
x
2 2
x 1
=
x1
+
x2
{ } x ∞
= max
x1 ,
x2
例2 设 x, y ∈E 则,x-y 的模可以表示这两个元素的“接近程度”,若在 R2 空间中的
两个元素,x (1,1), y (2, 4) 可以有如下模的定义
(可加性) (u, u)=0。
定义了内积的线性空间称为内积空间。
例 ∀ u(x), v(x) ∈C2[a, b] 至少存在以下四种形式的内积
(u,
v)
=
b
∫a
u
⋅
vdx
(u,
v)
=
b
∫a
[u
⋅
v
+
u′
⋅
v′]dx
(u,
v)
=
∫b a
[u
⋅
v
+
u′
⋅
v′
+
u′′
⋅
v′′]dx
(u,
v)
=
b
∫a
ρ
(x
)⋅
b
u = ∫ ρu ⋅ udx a
而 u、v 正交则意味着
b
(u, v) = ∫ ρu ⋅ vdx = 0
a
即通常理解的 u、v 以ρ为权正交。
4、Schwarz 不等式
设 u、v 是内积空间的两个元素,t 为任一实数,则 tu-v 也是内积空间的一个元素,显
然,它自身的内积
(tu − v,tu − v) ≥ 0
许多数学物理问题的解存在于某一类函数空间中。换句话来说为了得到有意义的 解,必需明确解的存在空间。即对组成解的函数的类型作一限定。值得注意的是,如果 限制的过于严格会将一些解排除在外,限制过宽可能导致解无意义。有限元解的存在空 间为索伯列夫空间,是 Hilbert 空间的子空间。
§4-3 索伯列夫空间 HK
( ) ≤ 2 c12ϕ12 (x) + c2 2ϕ 2 2 ( x) 2
所以 c1ϕ1 ( x) + c2ϕ 2 (x) 也是(a, b)上平方可积的函数。所有(a, b)上平方可积的函数
组成一个线性空间,记作 L2 (a, b) 。 例3 C1[a,b]
若 ϕ1 (x) 、ϕ1′(x) 、ϕ 2 (x)、ϕ 2′ (x) 在[a, b]上连续,则
故定义在[a, b]上的所有连续函数组成一个线性空间。记作 C[a, b]。
例2 L2(a,b)
若ϕ1 ( x) 、ϕ 2 ( x)是(a, b)上平方可积的函数,即
存在,则
b
b
∫ ∫ ϕ1 (x)2 dx , ϕ 2 (x)2 dx
a
a
(c1ϕ1 ( x) + c2ϕ 2 (x))2
= c12ϕ12 ( x) + c2 2ϕ 2 2 (x) + 2c1c2ϕ1 ( x)ϕ 2 ( x)
则模
b
∫ u = u ⋅ udx = u L2 a
而 u、v 正交则意味着 b ∫ u ⋅ vdx = 0 a
例如在[0, π]上当 m≠n 时 sinmx 与 sinnx 正交。
第 4 页 共 17 页
有限元分析与应用
霍战鹏
例2 若ρ>0 且积分 存在,也可定义 则模
b
∫ ρu ⋅ vdx
a
b
(u, v) = ∫ ρu ⋅ vdx a
(iii) 对任意 x、y∈E 有 ‖x+y‖≤‖x‖+‖y‖ (此式又称三角不等式)。
只要满足这些要求,均可作为一种模的定义。‖x‖描述了元素 x 的“大小”,‖x-y‖
第 2 页 共 17 页
有限元分析与应用
霍战鹏
则描述了两个元素 x、y 之间的“距离”。设真实解为 u,有限元解为 uh,如果当网格无限细 分时有‖u-uh‖→0,则说明有限元解收敛于真实解。模的定义不同,收敛的意义也不同。
y (2,4)
x−y = 2
(x1 − )y1 2 + (x2 − ) y2 2 =
10
x− y 1
=
x1 − y1
+
x2
− y2
=4
{ } x − y ∞ = max x1 − y1 , x2 − y2 =3
可见模的定义不同,其意义不同, 在实数域内,
10 3
x1 (1,1)
模 x 与绝对值 x 是等价的。
(iii) (α + β )x = αx + βx (iv) α (x + y) = αx + αy
若K为实数域则E称为实线性空间,K为复数域则E称为复线性空间。
例1 C[a,b]
若 ϕ1 ( x)、ϕ 2 ( x) 是[a, b]上的连续函数,则 c1ϕ1 + c2ϕ 2 也是[a, b]上的连续函数。
u − v = (u − v,u − v)
3、正交性
内积空间与一般线性空间的不同之处是可以用内积来定义两个元素之间的正交关系,函
数之间的“正交”。
若( u、v)=0 则称 u、v 正交。显然,模及正交性涵义取决于内积的定义。
例1 若积分 存在,则可定义
b
∫ u ⋅ vdx
a
b
(u, v) = ∫ u ⋅ vdx a
u
⋅
vdx
其中 ρ(x)是[a, b]上的给定函数。
2、内积模 对于内积空间,可以直接利用内积来定义元素的模
u = (u,u)
模的定义并不是唯一的,如果确有必要,我们仍然可以再定义其他形式的模(例如一致 模)。但对于内积空间用得最多的当然是由内积定义的模。
在内积空间E中, ∀u,v ∈ E ,u 与 v 之间的距离可用内积模表示为
例1 若 ϕ1、ϕ2Lϕ n 线性无关,则所有形式为
c1ϕ1 + c2ϕ 2 + L + cnϕ n
的试探函数组成 n 维线性空间。而所有形式为
n
∑ u = uiϕi i =1 n
∑ v = viϕ i i =1
的位移场则组成 2n 维线性空间。
例2 由于可以找出任意多个线性无关的连续函数(例如:1、x、x 2 L x n L ),
图4-1
(2)两种常用的模
下面是最常见、也是具有代表性的两种模
① 一致模 若 u∈C[a, b],则 u 必在[a, b]上取到最大值和最小值,故可按如下方式定
义一致模‖u‖∞。
u ∞ = max u a ≤ x≤b
② L2 模
若 u∈L2(a, b) 则
b
∫ L2 = u 2dx a
存在,可以按如下方式定义 L2 模
间的特例。
例5 Pn(x) [a,b] 定义在[a,b]上的 n 次多项式 Pn(x) [a,b]⊂C[a,b] 2、线性空间的维数
构成线性空间。
(1)线性相关与线性无关
设 ϕ1、ϕ2Lϕ n 为线性空间E的 n 个元素 (i)若存在不全为零的常数 c 1、 c 2… c n 使得
c1ϕ1 + c2ϕ 2 + L + cnϕ n ≡ 0
v 的内积,记作(u、v),且满足
(i) ( u、v)=(v、u)
(对称性)
第 3 页 共 17 页
有限元分析与应用