第六讲 有限元

合集下载

有限元方法的基本原理

有限元方法的基本原理

有限元方法的基本原理
有限元方法是一种数值分析方法,用于求解复杂结构的力学问题。

其基本原理如下:
1. 将结构离散化:首先将结构分割成许多小的单元(有限元),每个单元可视作一个简单的结构部件。

这样可以将原始连续结构的复杂问题简化为每个小单元的简单问题。

2. 定义弯曲关系:对每个单元建立力学模型,包括定义材料的弹性模量、泊松比、截面积等力学性质参数。

3. 建立单元的位移方程:利用有限元方法,采用适当的形函数,建立每个单元的位移方程,一般为不定位移分析。

4. 组装全局方程:将所有单元的位移方程组装成整个结构的全局方程。

5. 求解方程组:通过数值方法(如高斯消元法、迭代法等),求解结构的位移和应力等力学量。

6. 分析结果:根据结构的位移和应力等力学量,可对结构的强度、刚度、振动等进行分析和评价。

有限元方法的基本原理是将复杂结构的力学问题通过离散化处理,化为易于计算的小单元问题,再通过数值方法求解整个结构的力学行为。

6-第六讲-杆件有限元分析

6-第六讲-杆件有限元分析

3 0 1 2 k ji 2 2 k jj 3
补充-整体分析
补充-整体分析
整体刚度矩阵的建立 (3) 编码法
杆件有限元分析:案例
四 杆 桁 架 结 构
杆件有限元分析:案例
(1) 结构的离散化与编号
四 杆 桁 架 结 构
对该结构进行自然离散,节点编号和单元编号如上图所示,有关节点 和单元的信息见表1至表3。
杆件有限元分析:案例
(1) 结构的离散化与编号
四 杆 桁 架 结 构
对该结构进行自然离散,节点编号和单元编号如上图所示,有关节点 和单元的信息见表1至表3。
杆件有限元分析:案例
(2) 各个单元的矩阵描述
四 杆 桁 架 结 构
由于所分析的结构包括有斜杆,所以必须在总体坐标下对节点 位移进行表达,所推导的单元刚度矩阵也要进行变换,各单元 经坐标变换后的刚度矩阵如下。
基本概念 回顾:杆件有限元分析
(3) 单元应变场的表达 由弹性力学中的几何方程,有1D问题的应变
1D 杆 单 元
其中
叫做几何矩阵(strain-displacement matrix)。
基本概念 回顾:杆件有限元分析
(4) 单元应力场的表达 由弹性力学中的物理方程,有1D问题的应力
1D 杆 其中 单 元 叫做应力矩阵(stress-displacement matrix)
2D 杆 单 元
回顾:杆件有限元分析
(2)整体坐标系下的单元刚度方程
2D 杆 单 元
(3)整体坐标系下的单元应力
补充-整体分析
整体刚度矩阵的建立 (1) 位移转换法
补充-整体分析
整体刚度矩阵的建立 (1) 位移转换法
补充-整体分析

第六讲热温度场有限元

第六讲热温度场有限元

8-2 平面稳态温度场的有限元法
上式第一部分为内部单元的温度刚阵:
J e k [(bi2 ci2 )Ti (bib j ci c j )Tj (bibk ci ck )Tk ] Ti 4 A
对于内部单元的温度刚阵,i,j,k三点轮换,记为矩阵形式:
J e Ti bi2 ci2 J e k T j 4 A J e Tk
2T 2T 2 0 2 x y T k (Ta T ) n 在内 在1上
据变分原理,此问题等价于求泛函J[T(x,y)]的极值函数, 参考相关教材,可得上述热传导作为欧拉方程的相应泛 函:
k T 2 T 2 1 2 J [T ( x, y)] [( ) ( ) ]dxdy ( T TaT )ds 2 2 x y 1
考虑物体不包含内热源的情况,各向同性材料中的温度场满 足Laplace方程
2T 2T 2T 2 2 2 0 x y z
8-1 温度场问题的基本方程
在分析稳态热传导问题时,不需要考虑物体的初始温度分布 对最后的稳定温度场的影响,因此不必考虑温度场的初始条 件,而只需考虑换热边界条件。计算稳态温度场实际上是求 解偏微分方程的边值问题。温度场是标量场,将物体离散成 有限单元后,每个单元结点上只有一个温度未知数,比弹性 力学问题要简单。进行温度场计算时有限单元的形函数与弹 性力学问题计算时的完全一致,单元内部的温度分布用单元 的形函数,由单元结点上的温度来确定。由于实际工程问题 中的换热边界条件比较复杂,在许多场合下也很难进行测量, 如何定义正确的换热边界条件是温度场计算的一个难点。
bi b j ci c j b2 c 2 j j

有限元的基本原理

有限元的基本原理

有限元的基本原理
有限元方法是一种数值计算方法,常用于求解工程问题中的连续介质力学问题。

其基本原理是将复杂的连续介质分割成有限数量的简单几何形状的子域,称为有限元,然后利用数学方法和计算机技术对每个有限元进行离散化处理。

基于有限元原理,我们可以得到以下步骤:
1. 离散化:将连续的物理问题离散化为有限个由节点和单元组成的网格,在每个单元上选择适当的方程形式。

2. 建立本构方程:根据材料的力学性质,建立适当的本构关系表达式,将其转化为数学方程。

3. 单元形函数:在每个有限元上选择适当的单元形函数,将物理问题转换为离散问题。

4. 求解:对离散化后的方程进行求解,得到节点的未知位移。

5. 后处理:根据得到的位移信息,计算相应的应力和应变,以及其他感兴趣的物理量。

有限元方法的精度和收敛性与网格的划分有关,更精细的网格可以得到更准确的结果,但也会增加计算量。

因此,有限元方法是一个权衡计算效率和精度的方法。

有限元方法广泛应用于结构力学、流体力学、电磁场等领域的
建模和仿真中,可以有效地分析和解决各种工程问题。

其应用范围涉及机械、航空航天、汽车、建筑、电子等多个工程领域,为工程设计和优化提供了有力的工具。

第六讲 空间问题的有限元方法

第六讲 空间问题的有限元方法
( 2)
和 u, 、

表示 i 和 两个 节 点沿局
部坐 标系 三个 坐标 轴 的节 点力 ,对 应 于 局部 坐标 系三 个坐 标 轴 的节 点 力矩 为
M j、 M j、 M io x y z

,、

, m3 3

同样 ,f , 和 两个 节点 相对局 部坐标 系 的线位
维普资讯
技 术 与讲 座 Tc il s n e naL s hc o e
第六讲 空间问题 的有限元方法
Th ii e e Fnt Elme tMe h d f r3 dme so s e n to o i n in
陈 乐 生
标系 ,原点 为节 点 i 连线 为 ,Y 与横截 面 , 一 , 两个惯 性主 轴平 行 。 了确 定 面 平面 的位置 ,必 为 须在 而 平 面上给 定一 个参考 点 k ,当然 k 点不 能 在 f上 。于是 ,我们 可 以得到局 部坐标 系与整 体 , 坐标 系之 间的变换 关系 如下 :
( 州大 学机 械 工程及 自动化 学院 ) 福
前 面 几讲 我们 介 绍 了平面 问题 的 有 限元分 析
— —


ห้องสมุดไป่ตู้
— —
方 法 及其 原 理 ,关 于空 间 问题 的有 限元分 析这 里 只介绍 工程上 常见 的空 间粱单元 和轴对 称 问题 。
, l



mI

称 为轴对 称 问题 ,这在 工程上 是一类 常见 的 问题 , 或 者说 许 多 问题 工程 实 际 问题可 以简 化 为轴 对 称 问题加 以解决 。 ‘

有限元分析第六章

有限元分析第六章

第六章 非第二章至第五章的讨论以最小势能原理为基础,要求在单元内假设的位移场(试探函数)满足协调条件(在不同的单元内可以假设不同的的位移场)。

满足协调条件的单元,它们的收敛性等问题已在第四章中做了研究。

等参数单元就是目前处理二阶问题时应用最广的一种协调单元。

此外,还有一些单元,它们不满足协调条件,但仍可以收敛到真实解,这类单元称为非协调单元,可以看成是对等参数单元的一种改进,目的在于:在计算量增加不多的情况下,使单元的实际精度有所改善。

对于四阶问题(例如板、壳),协调条件要求单元之间位移和位移的一阶导数(转角)连续。

在第七章中将会看到,实现上述协调条件不是件容易的事,而且为此要增加相当大的计算量,因而人们在自编程序中常常对非协调单元感兴趣。

本章只讨论二阶问题,主要包括:非协调元的构造和分析方法,非协调元的理论基础(显然不能再利用最小势能原理),收敛判别方法。

这些结论对四阶问题同样适用。

从关于非协调元的讨论中,读者可以看到,有限元方法有了坚实的数学基础以后,在构造方法时思路可以开阔很多。

§6-1Wilson 非协调元Wilson 非协调元可以看成是由等参数单元演变来的单元,现以二维情况为例。

1、母体单元 形函数 母体单元ê:边长为2的正方形 自然坐标:ξ、η取四个角点为节点,在单元内的序号为1~4。

形函数2、实际单元 e可看成母体单元ê经变换F 得到利用上面定义的形函数,坐标的变换可写成其中(x i , y i )为实际单元中节点的坐标。

至此,还看不出Wilson 非协调单元与上一章介绍的等参数单元之间的差别。

3、单元内假设位移场图6-1图6-2) (4~1)1)(1(41),(=++=i N i i i ηηξξηξe eF →ˆ: ∑∑====4141i i i i iiy N y x Nx )1()1(),()1()1(),(242341222141ηαξαηξηαξαηξ-+-+=-+-+=∑∑==i i ii i iv Nv u Nu (6-1-1)同四节点等参元相比,单元内假定的位移场多了四项:它们有如下特性:(1) 不影响节点处的位移值,故称αl 为非节点自由度或单元的“内自由度”。

有限元分析与应用 第6讲、等参单元

有限元分析与应用 第6讲、等参单元
i
我们可以看到,位移插值函数公式(3)和 坐标变换公式(4)具有完全相同的形式,它们 用同样数目的对应节点值作为参数,并有完 全相同的形状函数 N (ζ ,η ), 作为这些节点 值前面的系数,我们称具有这种特点的单元 为等参数单元
i
等参变换步骤: 等参变换步骤
1找变换 x = x(ξ ,η ), y = (ξ ,η ) ,使x0y面上的任意四边形变成在 上的边长为2的正方形.
1 4 1 4 1 4 1 4
(1 − ξ )(1 − η ) (1 + ξ )(1 − η ) (1 + ξ )(1 + η ) (1 − ξ )(1 + η )
利用节点处得(ξ,η)坐标,上式可以写成统一得形式:
1 Ni (ξ ,η ) = (1 + ξiξ )(1 + ηiη ) 4
其中(ξi,ηi)为
ξoη 面
2在 ξoη 面上构造多项式插值函数 N k (ξ ,η ) 满足µ = ∑ N k (ξ ,η )µ k
3再变回xoy即: µ = ∑ N k (ξ ( x, y ) η ( x, y ))µ k = ∑ N k ( x, y )µ k 由于在 ξoη 面交界两测 u是连续的,xoy 面上也同样连续,但现在 N k (x, y )已经不 再是x,y的多项式了.
等参数单元平面问题变换的有限元格式
前面讲的建立有限元计算格式的推导过程中,前几步的主要 目的是求出以节点位移表示的单元位移插值函数,或求出单元形 状函数,后几步的主要目的是求出单元刚度矩阵,然后是用已知节 点位移计算应力。对于等参数单元,上面得到了四节点四边形等 参数单元的形状函数,下面主要讨论单元刚度矩阵的形成。 单元应变—单元位移—节点位移之间的关系. 由平面问题几何方程和位移插值公式(3)有:

第六讲热传导过程有限元分析

第六讲热传导过程有限元分析

V
(c
u t
u
kx
u x
u
x
ky
u y
u )dV
y
Q udV
V
q q0 ud
瞬态热传导有限元分析 ➢工程背景
一个长方形截面的冷空气通道,几何模型如下图所示。假设在垂直于纸面的方向上,通道内的初 始温度为0℃。通道的导热系数为0.044W/m·℃,比热和密度的乘积为1J/(m3·℃),内壁维持在0℃, 外壁与流体发生对流交换,且与周围环境间的热换系数为10 W/m·℃,环境温度为30℃,求3s后通道 壁面中的温度和热流密度。
a场体单元材料参数图
a场边界单元材料参数
b场单元材料参数
➢前处理
点击工具栏中“前处理”按钮进入GID。 注:进入GID后要进行ELAB1.0的数据转化data→problemtype→ELAB 几何建模: 首先建立一个小的矩形面,利用gid中copy命令中的拉伸功能建立如下图所示的几何模型,详细步骤 可以参考《有限元分析基础与应用》相关章节。
温度场u分布云图
热流场x方向分布云图
热流场y方向分布云图
➢有限元语言描述文件
为生成该问题有限元计算的所有程序源代码,针对之前的ELAB有限元分析得到的微分方程弱 形式,ELAB软件提供简洁的有限元语言描述文件,包括微分方程描述文件、多物理场描述文件以 及求解命令流控制文件。
针对该问题的有限元描述文件包括heatxy.fde(温度场fde文件), hfxy.fde(热流场fde文件), heat.mdi, heat.gcn ✓微分方程描述文件heatxy.fde(温度场fde文件)
START a l1: BFT SOLVC a SOLVSTR b a gidres(coor0); if (stop==0) goto l1;
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

二. 弹性力学研究问题的基本方法
1. 微分原理——取微小体积元dxdydz进行分析。需用到( 偏)微分方程, 2. 三大类方程 1)力的平衡方程 2)几何变形方程 3)材料的物理方程 3. 边界条件的考虑 1)位移方面 2)外力方面 三. 2维的情况 1.力的平衡方程 1).在物体内(x,y)点处取一微元dxdy, 厚 度为单位长度 注意: 四个侧面上力的表示
v v y (v dy v ) / dy y y
x
u
u u u dy u u y y y ad边转角: tan v v y v dy dy v 1 y 1 y y u
x
y
u x v y Biblioteka y v u x y
3.物理方程 ——广义Hooke定律
1 ( x y ) E 1 y ( y x ) E xy 2(1 ) xy xy G E
x
E ( x y ) 2 1 E y ( y x ) 2 1 E xy G xy xy 2(1 )
② y方向——合力为0
( y y y dy)dx ( y )dx ( xy xy x dx)dy xy dy Ydxdy 0
y
y y
Y——y方向单位体积所受的体积力 上式化简为:
y y

xy x
yx yx
x
1 [ x ( y z )] E 1 y [ y ( z x )] E 1 z [ z ( x y )] E
xy
xy
yz
zx
2(1 ) yz G E 2(1 ) zx zx
0 u x v x
1 x E y 2 xy 1 0
物理 方程
1 ( x y ) E 1 y ( y x ) E 2(1 ) xy xy xy
1
1 0
x 0 y xy 1 2 2(1 ) 0
平面应变问题弹性矩阵D
§4.2 平面问题的有限元模型
平面问题有限元分析的特点
1. 连续体的离散化——划分单元与节点 ● 在杆件、刚架结构结构中,一根杆、梁可作为上个单元,单元之 间的端点就是节点。这属于——自然分割(自然划分)。 ●对于连续体,则需要人为将其分割成单元,单元之间以有限个节 点连接。
y y
1.平面应力问题 — z 0, z 0 1). 研究对象——等厚薄板。 2).外载 ——沿xoy平面,在z方向无变化。 3).应力 ——在板上下面
t
t 2
x
z
z |
z
t 2
0,
zx |
z
0,
zy |
z
t 2
0
由于板很薄,故认为在整个板内,有:
平衡 方程
xy yx
x
几何 方程
y
u x v y
v u x y
xy
x
x x y 0 xy y
E ( x y ) 2 1 E y ( y x ) 1 2 E xy G xy xy 2(1 )
b-c上的 正应力
a-d上的 正应力
经过dx后正 应力增量
2).平衡关系 ① x方向——合力为0
yx x ( x dx)dy ( x )dy ( yx dy)dx yxdx Xdxdy 0 x y
注:X——x方向单位体积所受的体积力
上式化简为:
x yx X 0 x y
G E
yz
G

2(1 ) xy E
给定一定的边界条,由上述方程, 从理论上讲, 就能求出各应力, 应变, 位移分量. 由于方程中存在偏微分,故求解较为因难,只有少数一些非 常简单的受力和边界下能得到解析解,工程中一般采用近似或数值计 算方法——有限差分法或有限元法。
五. 弹性力学平面问题 工程中有许多构件与受力状态可以简化为二维情况处理, 这就是弹性力学 的平面问题. 平面问题有两种情况:
x ( x , y ) 2 x ( x , y ) 2 x ( x dx, y ) x ( x , y ) dx dx 2 x 2! x
略去高阶项
( x , y ) x ( x dx , y ) x ( x , y ) x dx x
G E
平面应力问题与平面应变问题的平衡方程和几何方程均相同,物理方 程形式相同,系数要作变换。 平面应 力问题
E E 1 2


1

平面应 变问题
平面应变问题物理方程也可表示为:
1 x E (1 ) y xy (1 )(1 2 ) 1 0
x
四. 3维的情况
位移: 应力:
应变:
u, v , w
x , y , z , xy , yz , zx
x , y , z , xy , yz , zx
x yx zx X 0 x y z y xy zy Y 0 y x z z yz xz Z 0 z y x xy yx , yz zy , xz zx
x
u u dx x方向位移—— x v v dx y方向位移—— x
x方向应变=ab在X方向变形量/ab长度=ba在X方向位移差/ ab长度
故:
x (u
u u dx u) / dx x x
v
ab边转角:
v v v dx v v x tan x x u u u dx dx u 1 x 1 x x x
§4.1 弹性力学的基础知识
一. 弹性力学的研究内容及基本假设 1. 材料力学的研究内容——规则形体(杆,轴,梁等)应力与变形规律。 2. 弹性力学的研究——不规则形体应力与变形规律。 3. 弹性力学对研究对象的基本假设——5点 1). 连续性——物质内部无空穴,可用连续函数描述。 2).均匀性——物体内各位置的物质具有相同的性质。 3).各向同性——在同一位置的各方向具有相同的性质。 4).线弹性——变形与外力呈线性关系,外力去除后,物体恢复原状。 5).小变形——物体变形量远小于其尺寸,分析时可略去高阶小量(二 阶以上)
xy dxdy
得:
xy 2x
dx dy yxdxdy
2
yx 2y
dxdy2 0
xy yx
剪应力互等定理
平衡方程
xy yx
2.几何变形方程
x yx X 0 x y y yx Y 0 y x
d
1.力的平衡方程
剪应力下标的含义:
xy
剪应力与y轴平行
剪应力所在平面垂直于x轴
2.几何变形方程
x
y
z
v y w z
u x
xy

yz
xz
v u x y w v y z u w z x
3.物理方程 ——广义Hooke定律
1 1 1 2 x [ x ( y z )] { x [ y ( x y )]} ( x y) E E E 1 2 1 1 1 y [ y ( x z )] { y [ x ( x y )]} ( y x) E E E 1 2(1 ) xy xy xy
y
y y
yx yx
x
d
c
xy xy
x x
①. 面上的应力表示成两个分量:正应 力和剪应力
②同一则面上应力均匀公布 ③侧面位置变化dx(dy)后,应力有增量
xy
a ( x, y )
b
y yx
x
●关于应力增量的计算
由Talyor公式展开:
2). a-d边 d点移至d’: y方向位移—— x方向位移——
v dy y u u dy y v v v dy y
y
u u dy y
c
d
a d
( x, y )
dx
b
y方向应变=ad在y方向变形量/ad长度 v dy
c
b
=da 在Y方向位移差位移差/ ad长度
a
故:
Y 0
x
d
c xy xy
x x
③ 力矩平衡——对中点的合力矩为0
xy
xy a
( x, y )
b
y yx
x
yx dx dx dy dy ( xy dx)dy xy dy ( yx dy)dx yxdx 0 x 2 2 y 2 2
1). 研究对象——Z方向尺寸较大,故Z向的位移受限。认为Z向无应 变,但有应力。
2).外载 ——沿xoy平面,在z方向无变化。 3).应力 在垂直于Z方向,取一载面,
w 0,
y
zx 0, zy 0
但 z 0,
zx 0, zy 0,
4).基本量 位移:u, v 应力: x , y , xy
z 0,
4).基本量 位移:
zx 0,
zy 0
zx 0,
zy 0, 但 z 0
应变: x , y , xy
相关文档
最新文档