武汉理工大学有限元理论基础精品课件

合集下载

有限元课件第2章-单元分析精选全文完整版

有限元课件第2章-单元分析精选全文完整版

a1
1 2A
ui uj
xi xj
yi yj
um xm ym
1
a2
1 2A
1
ui uj
yi yj
1 um ym
(2-14)
1
a3
1 2A
1
xi xj
ui uj
1 xm um
式中, A为三角形单元的面积,有
1 A 11
2
xi xj
yi yj
(2-15)
1 xm ym
y
m(7)
i(2)
j(1)
x
特别指出:为使求得面积的值为正值,本单元节点号
(2-16)
1 2A
[(ai
bi
x
ci
y)
i
(a
j
b
j
x
c
j
y)
j
(am
bm
x
cm
y)
m
]

Ni
1 2A
(ai
bi x ci y)
(i, j, m)
(2-18)
位移模式(2-16)可以简写为
u Niui N ju j Nmum Ni i N j j Nm m (2-19)
式(2-19)中的Ni、Nj、Nm是坐标的函数,反应 了单元的位移形态,称为单元位移函数的形函数。数
Ni 0
0 Ni
Nj 0
0 Nj
Nm 0
0 Nm
u v
j j
缩写为
um
vm
{ f } [N ]{ } (2-20)
[N]为形函数矩阵,进一步写成分块形式:
[N ] [[ Ni ] [N j ] [Nm ]]

有限元经典PPT第4章

有限元经典PPT第4章

Pii Kiiui
Ki1u1 Ki2u2 Kiiui K u i,i1 i1
ui
n
Kiiui Kiiui
Kiju j
4.1.2 平面应力问题有限元的基本思想和瑞雷-里兹法
v3 f3y
3
u3
f3x
f1y v1 u1
1 f1x
v2 f2y u2
2 f2x
给定一个三角形单元和作用在角点上 的六个力,要求得六个角点的位移。 或者是要求三角形角点发生指定的位 移,在三角形三个角点如何加力?
很显然,问题的精确解很困难。采用 瑞雷-里兹法求近似式解
e号单元的三个节点I,j,k的力对应的 力的平衡方程是第2i-1,2i;2j-1,2j;2k1,2k个平衡方程
e号单元的三个节点I,j,k的位移是第 2i-1,2i;2j-1,2j;2k-1,2k个未知数
弹性模量:E 横截面积:A
1
1 L
2
2L
3
局部系单元刚度阵:
k
1
EA L
1 -1
-1
1
2 集成总刚:
0 1
解得:
ux uy
L EA
3.8284L
EA
i
j
第一类位移条件:
Ki1u1 Ki2u2 Kiiui Ki1ui1
ui 0
令: Kij 0 i j
m
vi 0
Kii 1
um 0
Pi 0
ui 0
第二类位移条件:um um
大数
充大数法: Kii Kii
第一步:求转换矩阵
k2
EA 1 2L -1
-1
1
P
cos 0
T sin

武汉理工大学ansys讲义ppt课件

武汉理工大学ansys讲义ppt课件

1-1 工程和科学中典型问题
所研究问题的数学建模
力学模型
P
(平面应力问题)
P
有限元模型
第一章 有限元方法简介
1-2 有限元方法定义
有限元方法(The Finite Element Method, FEM)是计算机问世 以后迅速发展起来的一种分析方法。众所周知,每一种自然现 象的背后都有相应的物理规律,对物理规律的描述可以借助相 关的定理或定律表现为各种形式的方程(代数、微分、或积 分)。这些方程通常称为控制方程(Governing equation)。针 对实际的工程问题推导这些方程并不十分困难,然而,要获得 问题的解析的数学解却很困难。人们多采用数值方法给出近似 的满足工程精度要求的解答。有限元方法就是一种应用十分广 泛的数值分析方法。
的美国Ansoft公司。ANSYS公司成为全球最大的仿真软件公司。 目前,ANSYS整个产品线包括 结构分析(ANSYS Mechanical)
系列, 流体动力学(ANSYS CFD(FLUENT/CFX))系列,电子设计 (ANSYS ANSOFT)系列以及ANSYS Workbench和 EKM等。产品 广泛应用于航空、航天、电子、车辆、船舶、交通、通信、建筑、
节点
单元
2.有限元法物理上的理解
☆将连续体分割(离散)为有 限个、且按一定方式相互联结 在一起的小单元的组合体(单 元之间在节点处铰接)。用该 离散结构(单元组合体)近似 代替原来的连续体。
☆如果合理地求出各小单元的力学特性,就可以求出单元组合 体(离散结构)的力学特性,从而在给定的载荷和约束条件 下求出各节点的位移,求出各单元的应力。
间,又做不出什么东西,没有成就感,容易产生心理疲劳,缺
乏耐心。“苦中作乐”应是学ANSYS所必须保持的一种良好

有限元的核心思想和基本概念ppt课件(共17张PPT)

有限元的核心思想和基本概念ppt课件(共17张PPT)
内力:在外力作用下,物体内部不同部分 应力=内力/横截面面积
➢ 内力:在外力作用下,物体内部不同部分之间的相互作用力。
MSC-NASTRAN软件在航空航天领域有着很高的位置,目前世界上规模最大的有限元分析系统。
之间的相互作用力。物体横截面上的合力。 4、由单一构造场求解开展到耦合场问题的求解
构造力学:研讨有许多杆件组成的杆系的内力,位移。 许多商业化有限元分析软件都开发了和著名的CAD软件〔例如Pro/ENGINEER、Unigraphics、SolidEdge、SolidWorks、IDEAS、
第14页,共17页。
➢ 4、由单一构造场求解开展到耦合场问题的 求解
➢ 如今用于求解构造线性问题的有限元方法 和软件曾经比较成熟,开展方向是构造非 线性、流体动力学和耦合场问题的求解。 例如当流体在弯管中流动时,流体压力会 使弯管产生变形,而管的变形又反过来影 响到流体的流动……这就需求对构造场和 流场的有限元分析结果交叉迭代求解,即 所谓\"流固耦合\"的问题。
3、由求解线性问题开展到求解非线性问题 弹塑性阶段:去除外力物体不能恢复到外力作用前的外形。
用户自定义流场边境条件、用户自定义构 2、更为强大的网格处置才干 (技术难题,关键步骤)
给用户一个开放的环境,允许用户根据本人的实践情况对软件进展扩展,包括用户自定义单元特性、用户自定义资料、用户自定义流场
造断裂判据和裂纹扩展规律等等。 边境条件、用户自定义构造断裂判据和裂纹扩展规律等等。
杆件:长度远远大于横截面高度的构件。
内力,位移。
第4页,共17页。
➢ 应力:物体横截面上单位面积上的内力。 ➢ 应力=内力/横截面面积 ➢ 应变:单位长度上的位移。 ➢ 应变=位移/构件长度 ➢ 弹性阶段:去除外力物体还能恢复到外力

有限元基础教学课件PPT

有限元基础教学课件PPT

ε E T u (几何线性)
为梯度矢
ε u 一一对应,多连通域中未必一一对应. 在单连通域中:
31
§0.2 应力分析
取P点处一微平行六面体与xyz平行, 决定P点应力状态的6个分量记为
ζ x y z yz zx xy
f f x fy fz
T
T
ε E u,
T
u : u u : P E ν ζ

p


物体表面 u , 取未知函数 u ,经代换
: E DE u f 0 : u : u u
T
Px, y, z
: P E ν DET u (位移表示的应力边界条件)
14
应用领域:机械工程
(a) 铲运机举升工况测试
(b) 铲运机插入工况有限元分析
WJD-1.5型电动铲运机
15
液压挖掘机
(a) KOMATSU液压挖掘机
(b) 某液压挖掘机动臂有限元分析
16
驾驶室受侧向力 应力云图
接触问题结构件 应力云图
17
液压管路速度场分布云图
磨片热应力云图
支架自由振动云图
称为弹性矩阵
34
ζ Dε 或 ε D 1ζ
1 1 1 D E 0 0 0 1 0 0 0 1 0 0 0 0 0 0 21 0 0 0 0 0 0 21 0 0 21 0 0 0 0
i 1
RB
m
(Gu g ) 0
i 1
m
为了消除残差,通常引进内部权函数 WI 和边界权函 数WB ,将它们分别与 RI 和 RB 相乘,列出消除内部残 值方程式及消除边界方程式分别如下: RIWI dv 0 V C j ( j 1,2,, n) m S RBWB ds 0

有限元基础-上课件

有限元基础-上课件

总结词
有限元方法在电-磁场分析中能够模拟电磁 场分布和相互作用,为电磁装置设计提供精 确的预测。
详细描述
有限元方法在电-磁场分析中,能够考虑电 场强度、磁场强度、电流等参数,以及电磁 场与物质的相互作用。这为电磁装置设计提 供了精确的预测,如变压器、电动机、发电 机等的设计,以确保其性能和稳定性。
06
04
有限元方法的基本步骤
选取单元体与划分网格
选取单元体
选择适合问题特性的单元体,通常选 择容易解析和计算的几何形状,如三 角形、矩形等。
划分网格
将问题域分解成由单元体组成的网格 ,每个单元体之间通过节点相连。
建立单元体的刚度矩阵与质量矩阵
建立刚度矩阵
根据单元体的力学特性和边界条件,建立单元体的刚度矩阵,反映了单元体抵 抗变形的能力。
热传导分析
总结词
有限元方法在热传导分析中能够模拟热 量的传递和分布,为热工设计和优化提 供依据。
VS
详细描述
有限元方法在热传导分析中,能够考虑热 量的产生、传递和分布,以及材料热物理 性质的影响。这为热工设计和优化提供了 依据,如电子设备、机械零件、建筑保温 等的设计,以实现高效、稳定的热管理。
电-磁场分析
弹性力学本构方程
本构方程的数学表述
01
描述了材料的应力应变关系。
线弹性本构
02
材料在受力后会发生形变,但这种形变是可逆的,与应力大小
成正比。
非线性本构
03
材料在受力后发生的形变与应力大小不成正比,呈现出非线性
关系。
弹性力学边界条件与初始条件
边界条件
物体在边界上受到的力或位移约 束。
初始条件
物体在初始时刻的位移和速度状 态。

有限元法基础ppt课件

有限元法基础ppt课件

有限单元法
一、数值模拟方法概述 二、有限单元法简介 三、有限单元法分析步骤 四、利用有限元软件进行工程分析
一、数值模拟方法概述
工程技术领域中的许多力学问题和场问题,如固 体力学中的位移场、应力场分析、电磁学中的电磁 分析、振动特性分析、热力学中的温度场分析,流 体力学中的流场分析等,都可以归结为在给定边界 条件下求解其控制方程的问题。
结构矩阵分析方法认为:整体结构可以看作是由有限 个力学小单元相互连接而组成的集合体,每个单元的 力学特征可以看作建筑物的砖瓦,装配在一起就能提 供整体结构的力学特性。
结构矩阵分析方法分析的结构本身都明显地由杆件组 成,杆件的特征可通过经典的位移法分析建立。
虽然矩阵位移法整个分析方法和步骤都与有限单元法 相似,也是用矩阵来表达、用计算机来求解,但是它 与目前广泛应用的有限单元法是有本质区别的。
❖ 国际上早在20世纪50年代末、60年代初就投入大量的人力和 物力开发具有强大功能的有限元分析程序。其中最为著名的是 由美国国家宇航局(NASA)在1965年委托美国计算科学公司 和贝尔航空系统公司开发的NASTRAN有限元分析系统。该系 统发展至今已有几十个版本,是目前世界上规模最大、功能最 强的有限元分析系统。
有限元法
既可以分析杆系结构,又分析非杆系的连续 体结构。
三、有限单元法简介
有限单元法的常用术语:
有限元模型 是真实系统理想化的数学抽象。
定义
真实系统
有限元模型
自由度(DOFs- degree of freedoms)
自由度(DOFs) 用于描述一个物理场的响应特性。
UY ROTY
ROTZ UZ
UX ROTX
目前在工程技术领域内常用的数值模拟方法有: 1、有限单元法FEM( Finite Element Method) 2、边界元法BEM(Boundary Element Method ) 3、有限差分法FDM( Finite Difference Method 4、离散单元法DEM(Discrete Element Method) 其中有限单元法是最具实用性和应用最广泛的。

《有限元基础》课件

《有限元基础》课件
广泛适用性
有限元方法可以应用于各种物理问题和工程领域 ,如结构力学、流体力学、热传导、电磁场等。
高效性
有限元方法采用分块逼近的方式,将整体问题分 解为多个子问题,从而大大降低了问题的规模和 复杂度,提高了计算效率。
精度可控制
通过选择足够小的离散元尺寸和足够多的元数目 ,可以控制求解的精度,使得结果更加精确可靠 。
有限元方法对初值和边界条件 的选取比较敏感,不同的初值 和边界条件可能导致截然不同 的结果。
高阶偏微分方程的离散化 困难
对于一些高阶偏微分方程,有 限元方法的离散化过程可能会 变得相当复杂和困难。
有限元方法的发展趋势
并行化和高性能计算
随着计算机技术的发展,有限元方法的计算效率和精度得到了极大的提高。未来,随着并行化和高性能计算技术的进 一步发展,有限元方法的计算效率将会得到进一步提升。
02
有限元的数学基础
线性代数基础知识
向量与矩阵
介绍向量的基本概念、向量的运算、矩阵的表示和基 本运算。
线性方程组
阐述线性方程组的基本概念、解法以及在有限元分析 中的应用。
特征值与特征向量
介绍特征值和特征向量的概念、计算方法以及在有限 元分析中的应用。
变分法基础知识
变分法的基本概念
阐述变分法的基本思想、定义和定理,以及在 有限元分析中的作用。
弱收敛与弱*收敛
03
介绍弱收敛和弱*收敛的概念、性质以及在有限元分析中的应用

03
有限元方法的基本步骤
问题的离散化
总结词
将连续的问题离散化,将连续体划分为有限个小的单元,每个单元称为有限元 。
详细描述
在有限元方法中,首先需要对实际问题进行离散化,即将连续的问题划分为有 限个小的单元,每个单元称为有限元。离散化的目的是将连续的物理量近似为 离散的数值,以便进行数值计算。
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
弹性力学也是研究弹性体在外力作用下的平衡和运动, 以及由此产生的应力和变形。
2、研究的对象:有相同也有区别。
材料力学基本上只研究杆、梁、柱、轴等杆状构件,即 长度远大于宽度和厚度的构件。弹性力学虽然也研究杆 状构件,但还研究材料力学无法研究的板与壳及其它实 体结构,即两个尺寸远大于第三个尺寸,或三个尺寸相 当的构件。
果更精确,因而应用的范围更广泛。
但是,弹性力学也有其固有的弱点。由于研究对象 的变形状态较复杂,处理的方法又较严谨,因而解 算问题时,往往需要冗长的数学运算。为了简化计 算,便于数学处理,它保留了材料力学中关于材料 性质的假定。
10
弹性力学中关于材料性质的假定
(1) 物体是连续的,
即物体整个体积内部被组成这种物体的介质填满,不留 任何空隙。这样,物体内的一些物理量,如应力、应变、 位移等等才可以用座标的连续函数来表示。
20
应 变
体素的变形可分为两类:一是长度的变化,二是角度的变化。
任一线素的长度的变化与原有长度的比值称为线应变(或
称正应变),用符号 来表示。沿坐标轴的线应变,则加上 相应的角码,分别用 x、y、z 来表示。当线素伸长时,其线
应变为正。反之,线素缩短时,其线应变为负。这与正应力 的正负号规定相对应。
剪应力
加上两个角码,前一个角码表明作用面垂直于哪 一个坐标轴,后一个角码表明作用方向沿着哪一
个坐标轴。例如,剪应力 xy是作用在垂直于X轴
的面上而沿着y轴方向作用的。
15
2-2 应力的概念
应力的正负
如果某一个面上的外法线是沿着坐标轴的正方向, 这个面上的应力就以沿坐标轴正方向为正,沿坐标 轴负方向为负。
线素AB的正应变为:
x
=
(u
u dx) u x
dx
=
u x
同理,AD的正应变为:
(v v dy) v
y =
y dy
= v 2y2
求剪应变 xy ,也就是线素AB与AD之间的直角的改变
X向线素AB的转角 a , Y向线素AD的转角 b
y
u u dy yvv ydyD"b D'
D
C
A' u
沿坐标轴正向为正,负向为负。
体力,是分布于物体体积内的外力,如重力、磁力、
惯性力等。
单位体积内的体力亦可分解为三个成分,用记号Px、Py、 Pz表示,
沿坐标轴正向为正,负向为负。
弹性体受外力以后,其内部将产生应力。
13
2-2 应力的概念
为了研究物体内某点P的应力,考虑一个弹性体内微小的平 行六面体PABC的受力情况,该微元体称为体素。
3
2-1 材料力学与弹性力学
有限单元法
— 本课程中所指的是有限单元法在弹性 力学问题中的应用。因此要用到弹性力学的 某些基本概念和基本方程。本章将简单介绍 这些概念和方程,作为弹性力学有限单元法 的预备知识。
4
2-1 材料力学与弹性力学
弹性力学 — 区别与联系 — 材料力学
1、研究的内容:基本上没有什么区别。
任意两个原来彼此正交的线素,在变形后其夹角的变化 值称为角应变或剪应变,用符号 来表示。两坐标轴之间 的角应变,则加上相应的角码,分别用 xy、 yz、 zx 来表示。 规定当夹角变小时为正,变大时为负,与剪应力的正负号规 定相对应
(正的 xy 引起正的 xy ,等等)。
21
应变分量与位移分量的关系
11
(4) 物体是各向同性的,
也就是说物体内每一点各个不同方向的物理性质和机械 性质都是相同的。
(5) 物体的变形是微小的,
亦即当物体受力以后,整个物体所有各点的位移都远小 于物体的原有尺寸,因而应变和转角都远小于1,
这样,在考虑物体变形以后的平衡状态时,可以用变形 前的尺寸来代替变形后的尺寸,而不致有显著的误差;
2 yzdxdz
dy 2
2 zydxdy
dz 2
=
0
简化得
yz = zy
剪应力互等 xy = yx, yz = zy, zx = xz (2-1)
17
平衡微分方程
当物体在外力作用下保 持静止时,称物体处于 平衡状态。
弹性体中的应力不是任 意的,必须满足静力平 衡条件。
在单元体处于三维应力 作用下,根据微元体所 受合力为零的条件,可 以导出直角坐标系中的 三维平衡方程式:
并且,在考虑物体的变形时,应变和转角的平方项或乘 积项都可以略去不计,这就使得弹性力学中的微分方程 都成为线性方程。
12
2-2 应力的概念
作用于弹性体的外力(或称荷载)可能有两种:
表面力,是分布于物体表面的力,如静水压力,一 物体与另一物体之间的接触压力等。
单位面积上的表面力通常分解为平行于座标轴的三个成 分,用记号qx,qy,qz 来表示,
ABCD---A'B'C'D'
求线素AB、AD的正应变 x、 y
y
u ?u dy ?y
v ?v dy ?y
D" b D'
D
C
A' u
C'
B'
a
v ?v dx
?x
dy
v
A
B
B"
u ?u dx
dx
?x
0
x
图 1-5
,用位移分量来表示:
A点在X方向的位移分 量为u;
B点在X方向的位移:
u u = u u dx x
一般说来,弹性体内各点的应力状态都不相同,因此, 描述弹性体内应力状态的上述六个应力分量并不是常量,而 是坐标x、y、z的函数。
六个应力分量的总体,可以用一个列矩阵 s 来表示:
s x
s
y
s
=
sxzy
=
s
x
sy
sz
xy
yz
zx T
(2-2)
yz
19
zx
2-3 位移及应变、几何方程、刚体 位移
弹性体在受外力以后,还将发生变形。物体的变形 状态,一般有两种方式来描述:
1、给出各点的位移;2、给出各体素的变形。
弹性体内任一点的位移,用此位移在x、y、z
三个坐标轴上的投影u、v、w来表示。以沿坐标轴
正方向为正,沿坐标轴负方向为负。这三个投影称
为位移分量。一般情况下,弹性体受力以后,各点
的位移并不是定值,而是坐标的函数。
图 2-4
Z
Y X
PA=dx,PB=dy,PC=dz 每一个面上的应力分
解为一个正应力和两 个剪应力,分别与三 个坐标轴平行
s 正应力
剪应力
14
2-2 应力的概念
s 正应力
s 为了表明这个正应力的作用面和作用方向,加 上一个角码,例如,正应力 x是作用在垂直 于x轴的面上同时也沿着X轴方向作用的。
完全确定该点的应变分量,它们就称为该点的应变分量。
六个应变分量的总体,可以用一个列矩阵 来表示:
x
0
0
x
y
0
=
xzy
yz
zx
=
0
y 0
y 0
x
z
0
z
u
v
=
B
0
w
y
(2-3-2)
z
0
x
26
刚体位移
由几何方程(2-3)可见,当弹性体的位移分量完全确定时, 应变分量是完全确定的。反过来,当应变分量完全确定 时,位移分量却不完全确定;这是因为,具有确定形状 的物体,可能发生不同的刚体位移。为了说明这一点,
b
=
v x
u y
0
x
图 2-5
24
以上是考察了体素在XOY一个平面内的变形情
况,
x
=
u x
y
=
v y
xy
=
a
b
=
v x
u y
同样方法来考察体素在XOZ和YOZ平面内的变 形情况,可得:
z
=
w z

yz
=
v z
w y

zx
=
w x
u z
联立得到几何方程,表明应变分量与位移分量之
间的关系。
x
=
s x
x
xy
y
xz
z
Px
=
0
s y
y
yx
x
yz
z
Py
=
0
s z
z
zy
y
zx
x
Pz
=
0
xy = yx xz = zx yz = zy
18
应力分量
可以证明:如果 s x、s y、s z、 xy、 yz、 zx 这六个 量在P点是已知的,就可以求得经过该点的任何面上的正应 力和剪应力,因此,这六个量可以完全确定该点的应力状态, 它们就称为在该点的应力分量。
u x

y
=
v y

z
=
w z
xy
=
u y
v ,
x
yz
=
v z
w y

zx
=
w x
u
z
(2-3-1)
25
应变分量矩阵
可以证明,如果弹性体内任一点,已知这三个垂直方
向的正应变及其相应的三个剪应变,则该点任意方向
的正应变和任意二垂直线间的剪应变均可求出,当然
也可求出它的最大和最小正应变。因此,这六个量可以
5
弹性力学 — 区别与联系 — 材料力学
3、研究的方法:有较大的区别。
相关文档
最新文档