弹性力学及有限元
弹性力学基础及有限单元法

第一章1、弹性力学的任务是什么弹性力学的任务是分析各种结构物或其构件在弹性阶段的应力和位移,校核它们是否具有所需的强度和刚度,并寻求或改进它们的计算方法。
2、弹性力学的基本假设是什么?为什么要采用这些假设?(1) 假设物体是连续的——物体内部由连续介质组成,物体中没有空隙,因此物体中的应力、应变、位移等量是连续的•可以用坐标的连续函数表示。
实际上,所有的物体均由分子构成,但分子的大小及分子间的距离与物体的尺寸相比是很微小的,故可以不考虑物体内的分个构造。
根据这个假设所得的结果与实验结果是符合的。
(2) 假设物体是匀质的和各向同性的一一物体内部各点与各方向上的介质相同,因此,物体各部分的物理性质是相同的。
这样,物体的弹性常数(弹性模量、泊松比)不随位置坐标和方向而变化。
钢材由微小结晶体组成,晶体本身是各向异性的、但由于晶体很微小而排列又不规则,按其材料的平均性质,可以认为钢材是各向同性的。
木材不是各向同性的。
(3) 假设物体是完全弹性的一一物体在外加因家(裁荷、温度变化等)的作用下发生变形,在外加固素去除后,物体完全恢复其原来形状而没有任何剩余变形。
同时还假定材料服从胡克定律,即应力与形变成正比。
(4) 假设物体的变形是很小的——在载荷或温度变化等的作用下,物体变形而产生的位移,与物体的尺寸相比,是很微小的。
在研究物体受力后的平衡状态时,可以不考虑物体尺寸的改变。
在研究物体的应变时,可以赂去应变的乘积,因此,在微小形变的情况下弹性理论中的微分方程将是线性的。
(5) 假设物体内无初应力一一认为物体是处于自然状态,即在载荷或温度变化等作用之前,物体内部没合应力。
也就是说,出弹性理论所求得的应力仅仅是由于载荷或温度变化等所产生的。
物体中初应力的性质及数值与物体形成的历史有关。
若物体中有韧应力存在,则由弹性理论所求得的应力加上初应力才是物体中的实际应力。
上面基本假设中•假设(4)是属于几何假设,其他假设是属于物理假设。
弹性力学与有限元的关系

弹性力学和有限元关系:
弹性力学也称弹性理论,主要研究弹性体在外力作用或温度变化等外界因素下所产生的应力、应变和位移,从而解决结构或机械设计中所提出的强度和刚度问题。
在研究对象上,弹性力学同材料力学和结构力学之间有一定的分工。
材料力学基本上只研究杆状构件;结构力学主要是在材料力学的基础上研究杆状构件所组成的结构,即所谓杆件系统;而弹性力学研究包括杆状构件在内的各种形状的弹性体。
在有限单元法中,选择节点位移作为基本未知量时称为位移法;选择节点力作为基本未知量时称为力法;取一部分节点力和一部分节点位移作为基本未知量时称为混合法。
位移法易于实现计算自动化,所以,在有限单元法中位移法应用范围最广。
当采用位移法时,物体或结构物离散化之后,就可把单元总的一些物理量如位移,应变和应力等由节点位移来表示。
这时可以对单元中位移的分布采用一些能逼近原函数的近似函数予以描述。
通常,有限元法我们就将位移表示为坐标变量的简单函数。
这种函数称为位移模式或位移函数。
弹性力学边值问题及有限元法(PPT)

(3)在边界上部分给定面力,部分给定位移——混合边界条件
基本解法
弹性力学边值问题——基本方程+边界条件
给定作用在物体全部边界或内部的外界作用(包括温度影响、 外力等),求解物体内由此产生的应力场和位移场。
具体地说,对物体内每一点,当它处在弹性阶段,其应力分 量、应变分量、位移分量等15个未知函数要满足平衡方程、几 何力程、本构方程这15个泛定方程,在边界上并要满足给定的 全部边界条件。
通过与原问题基本方程及边界条件等效的变分原理,建立求 解的代数方程组,求解有限个节点上的场变量值
用有限个节点场变量值插值得到全求解域任意位置的场变量
单元内近似函程形式必须一样 单元内近似函数一般取Lagrange多项式
单元位移函数
对三角形单元,假定单元内的位移分量是坐标的线性函数
x
x
xy
y
xz
z
Fbx
0
yx
x
y
y
yz
z
Fby
0
zx
x
zy
y
z
z
Fbz
0
平衡方程的意义
受力而平衡的弹性体内 各应力之间(及其与体 力之间)的相互制约关 系
几何方程
x
u x
y
v y
z
w z
xy
u y
v x
xy
v z
w y
xy
w x
u z
应变与位移之间的关系, 以及应变之间的关系
物理方程
也叫本构方程
应力应变之间的关系
x
E(1 ) (1 )(1 2)
( x
1
y
弹性力学及有限元法:第2章 弹性力学中的若干典型问题及基本解法的讨论

2
2.1 平面问题
平面问题是工程实际中最常遇到的问题,许多工程实际问题 都可以简化为平面问题来进行求解。平面问题一般可以分为两类, 一类是平面应力问题,另一类是平面应变问题。
平面应力问题
平面应变问题
3
2.1.1 平面应力问题
平面应力问题的特征:
(1)所研究的对象在z方向上的尺寸很小(即呈平板状);
x y
2 xy
I3 0
因此,求解平面应力状态下主应力的方程为
3 I1 2 I2 0 解出的平面应力状态下的主应力具体为式
1, 2
x
y
2
[(
x
y
2
)2
1
2 xy
]
2
3 0
7
(2.6) (2.7) (2.8)
8
2.1.2 平面应变问题
平面应变问题的特征:
(1)如图2-2所示,当物体z方向上的 y
对于轴对称问题,采用圆柱坐标r、、z比采用直角坐标x、y、
z方便得多。这是因为,当以弹性体的对称轴为z轴时(如图2-3所 示),则所有的应力分量、应变分量和位移分量都将只是r和 z的
函数,而与无关(即不随变化)。
z
(a)
o
x
13
2.2 空间轴对称问题
C
dz
PB
z
A
r
dr
d
o
d
r
r
z
z
r
z r z
dz
C
Z
z
z z
dz
r
z
r
r
z
dr
dz
r rz
R
(同济大学)第1讲_弹性力学及有限元方法概述

有限元分析
的一般规律物体在空间的位置随时间的改变
对象内容
任务
对象内容
任务
概述
ANSYS 静力分析z起重机械有限元应用
整机模态分析
车辆安全性
工件淬火3.06 min 时的温度、组织分布(NSHT3D)
同济大学
同济大学
金属反挤压成型:温度分布和变化铸造成型:温度变化和气泡
速度
压力导流管分析
超音速飞行压力分布汽车气动分析
高速导弹气动
同济大学
两根热膨胀系数不同的棒焊接在一起,加热后的变形情况
子结构方法分析大型结构的早期应用法
梁单元
建模时充分利用重复性。
9第2章弹性力学平面问题及空间问题有限元

假定的位移函数是多项式,它是连续函数,可以肯定,在单元内部位移函数是单值连续的。由于单 元的位移函数 u 、 v 都是坐标 x 、 y 的线性函数,在单元边界上位移也是线性变化的,两个相邻单元在 公共节点上具有相同的节点位移,因而相邻单元在公共边界上位移连续,即协调条件得到满足。 由上面分析可以看出,三角形常应变单元的位移模式可以保证计算结果的收敛。
px
py
px
py ]
T
(2-1-7b)
(2 )若在 jm 边上受线性分布的水平方向的面力,它在 j 点的集度为 q ,在 m 点的集度为零 (如图 2-5) 。可预计由该面力求得的等效节点载荷只有 R xj 、
R xm ,其余节点载荷分量必为零。
将 jm 边上的分布面力写成 s 的函数,为
s { p} [ (1 ) q 0]T l 在 jm 边上的形函数也需用变量 s 表示,根据形函数的含义,
Ve
[k ii ] [k ij ] [ k im ] [k ji ] [k ij ] [k jm ] [k mi ] [ k mj ] [k mm ]
式中, t 为单元的厚度,当单元划分得足够小时,可以认为每个单元的厚度 t 为常值。子阵为
(2-1-5)
[k rs ] [ Br ]T [ D][B s ]tA
101
二、 单元刚度矩阵 1、单元几何矩阵 [ B ] 有了单元的位移模式,利用平面问题的几何方程求得应变分量
0 x x u e e 0 { } [ L][ N ]{} [B ]{} y y v xy y x
弹性力学与有限元完整版

Z面 X面
•②应力符号意义
•正应力: 由法线方向确定
x、 y、 z
•剪应力: xy
作用面
作用方向
•符号规定:
正面上与坐标轴正向一致,为正;
负面上与坐标轴负向一致,为正。
正面 负面
Z面
X面
•③剪应力互等定理
xy yx
相等
yz zy
xz zx
4. 完全弹性假设
应力和应变之间存在一一对应关系,与时间及变形历史无关。满 足胡克定理。
5. 小变形假设
在弹性体的平衡等问题讨论时,不考虑因变形所引起的几何尺寸 变化,使用物体变形前的几何尺寸来替代变形后的尺寸。采用这 一假设,在基本方程中,略去位移、应变和应力分量的高阶小量 ,使基本方程成为线性的偏微分方程组。
大小和方向不同。
体力分量:将体力沿三个坐标轴xyz 分解,用X
、Y、Z表示,称为体力分量。
符号规定:与坐标轴方向一致为正,反之为负
。 应该注意的是:在弹性力学中,体力是指单位
体积的力 。
体力的因次:[力]/[长度]^3
表示:F={X Y Z}
② 面力
与体力相似,在物体表面上任意一点P 所受面力的大小 和方向,在P点区域取微小面积元素△S ,
压力,物体之间的接触力等。
集中力——作用物体一点上的力。(在弹性力学中一
般所受体力的大小和方向,在P点区域取
一微小体积元素△V, 设△V 的体力合力为△F,则
△V 的平均体力为
当△V 趋近于0, 则为P点的体力
体力是矢量:一般情况下,物体每个点体力的
第一篇 弹性力学
第一章 弹性力学基本方程
1.1 绪论 1.2 弹性力学的基本假定 1.3 几个基本概念 1.4 弹性力学基本方程
第2章_弹性力学基础及有限元法的基本原理1

W U
当外力的形式是多样的时,外力的虚功等于:
W f Pc f Pv dV f Ps dS
T T T v s
• 1.4 平面问题定义
严格地讲,任何结构都是空间的。对于某些特殊情 况,空间问题可以转化为平面问题。
(1)平面应力问题 满足条件: 1)几何条件 厚度尺寸远远小于截面尺寸; 2)载荷条件 载荷平行于板平面且沿厚度方向均匀 分布,而板平面不受任何外力作用。
1)位移函数 分片插值→ 假设一种函数来表示单元位移分布 一般选取多项式(简单而且易求导)
可用于离散的单元: • 三角形单元; • 矩形单元; • 不规则四边形单元。 DOF 节点的自由度:节点所具有的位移分量的数量。 一个单元所有节点的自由度总和称为单元自由度。 (1)单元参数只能通过节点传递到相邻单元 (2)单元和节点必须统一编号
2.2 单元分析(位移、应力、应变) 任务:形成单元刚度矩阵,建立单元特性方程 因此必须建立坐标系,如下图:
1D问题的弹性模量
E杨氏弹性模量
泊松比是指材料在单向受拉或受压时,横向正应变与轴向 正应变的绝对值的比值,也叫横向变形系数,它是反映材 料横向变形的弹性常数。 若在弹性范围内加载,横向应变εx与纵向应变εy之间存 在下列关系: εx=- νεy 式中ν为材料的一个弹性常数,称为泊松比。泊松比是 量纲为一的量。 可以这样记忆:空气的泊松比为0,45#钢0.3,水的泊松 比为0.5,中间的可以推出。
• 未知数 应力 6个+应变 6个+位移 3个=15个 • 方程个数 平衡方程 3个+几何方程6个+物理方程6个=15个 原则上可以根据15个方程求出15个未知物理量 但实际求解时先求出一部分再通过方程求解剩下的。 目前有限元法主要采用的是位移法,以三个位移 分量为基本未知量。位移-应变-应力,应力和外力平衡
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
(2) 面力
Q —— 面力分布集度(矢量) F lim S 0 S
F Xi Yj Z k
—— 作用于物体表面单位面积上的外力
z
Q
X Y Z —— 面力矢量在坐标轴上投影
单位: 1N/m2 =1Pa (帕)
Z
k i
x O j
X
S Y
y
1MN/m2 = 106Pa = 1MPa (兆帕)
yx dzdx zx zx dz dxdy zx dxdy Xdxdydz 0 z
由其余两个平衡方程
可以得出与之相似的两个方 程。化简,除以dxdydz,得
z
F
F
y
0
和
0
空间问题的平衡微分方程 (纳维叶方程)
如物体处于运动状态,根据达朗伯(d’Alembert)原理,在体 力项中引入惯性力:
第2章
§2-1 §2-2 §2-3 §2-4 §2-5
应力分析与平衡方程
一点的应力状态、应力张量 主应力与主剪应力、应力张量不变量 八面体应力、应力强度 应力球张量和应力偏张量 平衡方程——应力和外力的关系
§2.1 一点的应力状态、应力张量
基本概念: 外力、应力、形变、位移。 1. 外力
体力、面力
程
一.几何方程
PA=dx PB=dy C C’ P P’ A A’ B B’ PC=dz
研究在oxy平面内投影的 变形,
O
x
一点的变形 线段的伸长或缩短; 线段间的相对转动; 考察P点邻域内线段的变形:
P
u
P
B
u u dx x
dx A
v
dy
A
v v dx x
PA dx
变形前
y’
y
e1’
x’ x
数学上,对坐标变换时服从一定坐标变换式的9个数所定 义的量叫二阶张量,应力张量通常表示为:
4. 斜面上的应力
作斜面abc垂直于x’轴,该斜面上的应力矢量为P。P在 坐标系下的三个分量为Px, Py 和Pz ,则
z
由斜面应力(Cauchy)公式
z’
Pz
P
y’ Py
Px
y
x’ x
xy
xy
yx
y
x
x
在用应力莫尔圆时必须用此规定求解问题
3. 应力张量
考察物体内任一点o,设oxyz为旧坐标系,其 单位矢量为e1、e2、e3,相应的应力分量为
z z’
设ox’y’z’为新坐标,其单位矢 量为e1’、e2’ 、e3’ 。相应的 应力分量为
e3’ e1
e3 e2’ e2
剪应力互等定理
O x
xz xy y yx y yz x zx zy z
y
yx
zx
zy
yz
应力符号的意义:
第2个下标 y 表示τ的方向. 应力正负号的规定: 正应力—— 拉为正,压为负。 剪应力—— 坐标正面上,与坐标正向一致时为正; 坐标负面上,与坐标正向相反时为正。
由外力引起的在 P点的某一面上内力分布集度 应力的法向分量 应力的切向分量
P ΔA
ΔQ
n
(法线)
应力分量 单位:
—— 正应力 —— 剪应力
与面力相同
MPa (兆帕)
应力关于坐标连续分布的
( x, y, z ) ( x, y , z )
(2) 一点的应力状态
通过一点P 的各个面上应力状况的集合 —— 称为一点的应力状态 x面的应力: x , xy , xz y面的应力: z面的应力:
同样方法研究另外两平面yoz和zox上投影线元的变形可 得到类似的方程。综合起来,得弹性力学几何方程。也 称柯西(Cauchy)方程
几何方程
说明:
(1) 几何方程反映任一点的位移(3个分量)与该点应变(6个分量)
间的关系,是弹性力学的基本方程之一。
(2) 当 位移分量u、v 、 w已知,则6个应变分量可完全确定;
八面体剪应力为:
八面体剪应力对于塑性理论具有重要意义,为了使用方便, 将它乘以 ,并称之为应力强度,用符号 来表示, 即
§2.4 应力球张量和应力偏张量
描述一点应力状态的9个应力分量构成一个对称应力张量
其中
称为应力张量的分量。
引入平均应力 则
应力张量可分解为两个张量之和
简写为 式中
运动微分方程
第三章
应变分析与几何方程
第二节有关力学基本概念描述已知: *在载荷作用下,物体的形状和位置要发生变化, *力学中用应变来度量一点形状的改变;用位移来 度量一点位置的改变. 如已知物体中每一点的位移,则受载物体的位 置和形状均可确定.即位移与应变之间存在一定 的关系.
描述位移与应变之间关系的方程称为几何方
称
为应力偏量,
为应力球形张量,
为单位张量。
球形张量是代表各向均匀拉伸或压缩的应力状态。
球形张量应力(静水应力)作用下,物体只产生各向 相同的线应变而无剪应变。对应物体的体积改变,而形状 不变。
应力偏量代表各面正应力中偏离静水应力的量,是正应力 之和为零的应力状态。该应力状态下,物体的体积不改变 而形状改变。 静水压力实验研究表明,在均匀受力情况下,即使应力达到 很大值,材料也不产生塑性变形。 故:应力球形张量不产生材料的塑性变形; 应力偏量是产生塑性变形的真正原因。
的平衡方程 M ee ' 0
整理,并略去微量后,得 yz zy 同样可以得出 xy yx , zx xz
剪应力互等定理
列出x轴方向的力的平衡方程
F
x
0
x x dx dydz x dydz yx yx dy dzdx x y
现在主应力空间里,考察通过物体内任一点M这 的一个微分面,该微分面的外法向n与三个应力 主轴呈等倾斜。这样的微分面共有8个,它们可组 成一个包含点M在内的无限小的正八面体,如图 所示。这些微分面上的应力,就称为八面体应力。
由于这些斜面的法线的方向余弦的绝对值都相等:
同时有: 于是得:
带入正应力的计算公式,可得八面体正应力为:
xy
第1个下标 x 表示τ所在面的法线方向;
z
z
O x
xz xy y yx y yz x zx zy z
y
yx
zx
zy yz
与材力中剪应力τ正负号规定的区别: 规定使得单元体顺时转的剪应力τ 为正,反之为负。
yx y
x
y
xy yx
y , yx , yz
z , zx , zy
x xy xz 用矩阵表示: yx y yz zx zy z
其中,只有6个量独立。
z
z
xy yx yz zy zx xz
PB dy
变形后
y
P
P A
ห้องสมุดไป่ตู้
u v u u dx x
v v dy y
u
A
v v dx x
B
u u dy y B v v dy y
u dy y
B
注:这里略去了二阶以上高阶无穷小量。
PA的正应变: u u dx u u x x x dx PB的正应变: v v dy v y v y dy y P点的剪应变:
3. 形变 (1) 一点形变的度量
形变 —— 物体的形状改变 (1)线段长度的改变 ——用线(正)应变ε度量 (2)两线段间夹角的改变。 ——用剪应变γ度量 (剪应变——两垂直线段夹角(直角)的改变量)
三个方向的线应变:
三个平面内的剪应变:
x , y , z xy , yz , zx
上述方程为 的齐次线性方程组, 且常数项都为 零。因为: ,故 不能同时为零, 所以方程组的系数行列式应为零,即
将行列式展开,得到求解主应力 应力张量 的特征方程。
的三次方程,称为
式中
设特征方程的三个根为
展开后有
,则
(特征方程)
比较上两式,有
对一个给定的应力状态,其主应力的大小和方向是确定的, 不随坐标系的变换而变化。故 是不随坐标系的变换 而变化的量,称为应力张量不变量。 分别称为应力张量的第一、第二、第三不变量。
由此可见,过某点的任意斜面上的应力分量,都可以用过该 点的平行于坐标面的微分面上的9个应力分量来表示。写成 矩阵的形式,即:
斜面上的总应力为: 斜面上的正应力为:
§2.2 主应力与应力张量不变量
设斜截面外法线方向为 为
,它的方向余弦
应力矢量P在坐标轴上的投影为:
将上式展开
当斜面法线方向满足上述方程时,该斜面上只有正应力,没有剪应 力,称该平面为主平面;主平面上的正应力称为主应力;主应力方 向(即主平面法线方向)称为主方向。
z C
z
A
O
应变的正负: 线应变: 伸长时为正,缩短时为负;
剪应变: 以直角变小时为正,变大时为负; x
x P
y
B y
§2.5 平衡(运动)微分方程
在物体内的任意一点P,割取一
个微小的平行六面体,棱边的长度
分别为PA=dx,PB=dy,PC=dz。 首先,以连接六面体前后两面
中心的直线 ee ' 为矩轴,列出力矩
反之,已知6个应变分量,不能确定位移分量。 (∵积分需要确定积分常数,由边界条件决定。)