厚壁圆筒应力分析
06_压力容器应力分析_厚壁圆筒弹塑性应力分析
2.3 厚壁圆筒应力分析
2.3.2 弹塑性应力
因为“弹性筒”内壁面同时也是“塑性筒”的外 壁面，所以在交界面上( r=Rc )，也满足 Mises 条件
r R
c
r r R
c
2 s 3
联立上述三式，得到弹、塑性区界面压力pc的另一表达 式如下
pc
s R R
2.3.4 提高屈服承载能力的措施
（2）高压厚壁筒提高屈服承载能力的措施
2.3 厚壁圆筒应力分析
2.3.4 提高屈服承载能力的措施
下图为经过自增强处理后，单层厚壁筒中的应力 分布情况。自增强法最早出现于20世纪初，首先应用 于炮筒的制造。目前已经应用于石油化工中的高压及 超高压容器、超高压管道、超高压压缩机气缸等。
残余应力的计算是依据“卸载定理”的，参见教 材。该部分须掌握残余应力的分布图。
2.3 厚壁圆筒应力分析
2.3.2 弹塑性应力
2.3 厚壁圆筒应力分析
（1）爆破过程 OA：弹性变形 AB：进入屈服 BC：屈服并强化 CD：爆破 pc:塑性垮塌压力， 工程上称为爆破 压力。
2.3.3 屈服压力和爆破压力
2.3 厚壁圆筒应力分析
（2）理想弹塑性材料
2.3.2 弹塑性应力
对于理想弹塑性材料，忽略材料的硬化阶段，同 时认为材料的屈服极限为常数。
2.3 厚壁圆筒应力分析
（3）塑性失效准则
2.3.2 弹塑性应力
筒体为理想弹塑性材料，当屈服区扩展至外壁 面，使筒体整体屈服，此时承受的内压力为筒体承 受的最高极限载荷。 （4）屈服条件 当材料从弹性阶段进入理想塑性阶段时，应满 足一定的条件，以此来判定材料是否进入屈服阶段， 此条件称为“屈服条件”（屈服失效判据）。 常用的屈服条件有：Tresca屈服条件和Mises 屈服条件。
厚壁圆筒应力分析
温度变化引起的弹性热应力
热应力
构件之间热变形 的相互约束
构件热变形受到 外界约束
构件内部温度 分布不均匀
2.3 厚壁圆筒应力分析
2、厚壁圆筒的热应力
◆厚壁圆筒中的热应力由平衡方程、几何方程和物理方程， 结合边界条件求解。
◆当厚壁圆筒处于对称于中心轴且沿轴向不变的温度场时， 稳态传热状态下，三向热应力的表达式为：
应变
径向应变、轴向应变和周向应变
分析方法
8个未知数，只有2个平衡方程，属静不定问 题，需平衡、几何、物理等方程联立求解。
2.3 厚壁圆筒应力分析
2.3.1 弹性应力
p0
研究在内压、 外压作用下， 厚壁圆筒中的 应力。
图2-15 厚壁圆筒中的应力
2.3 厚壁圆筒应力分析
2.3.1 弹性应力
•以轴线为z轴建立圆柱坐标。 •求解远离两端处筒壁中的三向应力。
2.3 厚壁圆筒应力分析
c. 几何方程（续） 径向应变 周向应变 变形协调方程
（2-27） （2-28）
2.3 厚壁圆筒应力分析
d. 物理方程
（2-29）
2.3 厚壁圆筒应力分析
e. 平衡、几何和物理方程综合—求解应力的微分方程
（2－33）
2.3 厚壁圆筒应力分析
边界条件为：当
时，
；
当
时，
。
由此得积分常数A和B为：
2.3 厚壁圆筒应力分析
周向应力 径向应力
轴向应力
称Lamè（拉美）公式
（2-34）
2.3 厚壁圆筒应力分析
表2-1 厚壁圆筒的筒壁应力值
2.3 厚壁圆筒应力分析
(a)仅受内压
(b)仅受外压
图2-17 厚壁圆筒中各应力分量分布
厚壁圆筒应力分析剖析
厚壁圆筒应力分析剖析一、应力分析方法1.在应力分析中，通常采用静力学的方法，根据力学定律对厚壁圆筒进行应力分析。
2.厚壁圆筒的应力分析可以分为轴向应力、周向应力和切向应力三个方向上的应力分析。
二、应力计算公式1.轴向应力：σa=(P·r)/t其中，σa表示轴向应力，P表示圆筒受到的内外压力，r表示圆筒内径，t表示圆筒壁厚。
2.周向应力：σc=(P·r)/(2t)其中，σc表示周向应力。
3. 切向应力：τ = (P · ri) / t其中，τ 表示切向应力，ri 表示圆筒中心点到任意一点的径向距离。
三、实例分析假设有一个内径为 10cm，外径为 15cm，壁厚为 2cm 的厚壁圆筒，内外压力分别为 5MPa 和 10MPa。
现对该厚壁圆筒进行应力分析。
1.轴向应力：根据公式σa = (P · r) / t，代入 P = 5MPa，r = 7.5cm，t =2cm，计算得σa = (5×7.5) / 2 = 18.75MPa。
同理，代入 P = 10MPa，r = 7.5cm，t = 2cm，计算得σa =(10×7.5) / 2 = 37.5MP a。
2.周向应力：根据公式σc = (P · r) / (2t)，代入 P = 5MPa，r = 7.5cm，t= 2cm，计算得σc = (5×7.5) / (2×2) = 9.375MPa。
同理，代入 P = 10MPa，r = 7.5cm，t = 2cm，计算得σc =(10×7.5) / (2×2) = 18.75MPa。
3.切向应力：根据公式τ = (P · ri) / t，代入 P = 5MPa，ri = 7.5cm，t =2cm，计算得τ = (5×7.5) / 2 = 18.75MPa。
同理，代入 P = 10MPa，ri = 7.5cm，t = 2cm，计算得τ =(10×7.5) / 2 = 37.5MPa。
压力容器厚壁圆筒的弹塑性应力分析
未来发展方向和前景展望
THANK YOU
汇报人：XX
有限元法的优缺点及其在 工程实践中的应用案例
厚壁圆筒的弹塑性应力分析中的材料模型
理想弹塑性模型：假设材料在受力过程中遵循胡克定律，忽略材料的应变率效应 和温度效应。
弹塑性有限元法：将厚壁圆筒离散化为有限个单元，每个单元的应力应变关系通 过弹塑性本构方程描述。
增量理论：基于增量形式的本构方程，考虑了前一次加载时残留在材料中的应力 场对当前加载的影响。
厚壁圆筒的弹塑性应力 分析的未来发展
PART 01 添加章节标题
PART 02
压力容器厚壁圆 筒的弹塑性应力
分析概述
压力容器厚壁圆筒的结构特点
厚壁圆筒由金属材料制成，具有高强度和耐腐蚀性能。 厚壁圆筒的结构设计应满足压力容器的工艺要求和使用条件。 厚壁圆筒的厚度通常较大，以承受内压和其他附加载荷。 厚壁圆筒的制造过程中需要进行焊接、热处理、无损检测等质量控制措施。
PART 06
厚壁圆筒的弹塑 性应力分析的未
来发展
新型材料对厚壁圆筒弹塑性应力分析的影响
新型材料的出现将改变厚壁圆筒的弹塑性应力分析的边界条件和载荷条件。 新型材料的力学性能对厚壁圆筒的弹塑性应力分析的精度和可靠性提出了更高的要求。 新型材料的加工制造技术将促进厚壁圆筒的弹塑性应力分析方法的改进和发展。 未来将有更多的新型材料应用于厚壁圆筒的制造，需要进一步研究这些材料对弹塑性应力分析的影响。
提高压力容器的安裂而引起的安全事故 为压力容器的设计、制造和使用提供科学依据
PART 03
厚壁圆筒的弹塑 性应力分析方法
有限元法在厚壁圆筒弹塑性应力分析中的应用
有限元法的定义和原理
厚壁圆筒的弹塑性应力分 析的数学模型
05_压力容器应力分析_厚壁圆筒弹性应力分析
2.3 厚壁圆筒应力分析
三对截面取出微元体：
2.3.1 弹性应力
② 周向应力σϴ 径向应力σr
一对圆柱面，相距dr，σr作用于该面上。 一对纵截面，相差dθ，σϴ作用于该面上。 一对横截面，长度为1， σz作用于该面上。
根据轴对称性， σϴ和σr仅与r有关。
2.3 厚壁圆筒应力分析
Hale Waihona Puke 2.3.1 弹性应力周向应变
(mn m' n' )
(r w)d rd w rd r
上述二式为(2-27)式。
2.3 厚壁圆筒应力分析
对周向应变求导，有
2.3.1 弹性应力
dw r w d 1 dw w dr 2 dr r r dr r 1 r r
上式又称广义虎克定律。
(2 29)
2.3 厚壁圆筒应力分析
2.3.1 弹性应力
平衡、几何和物理方程综合－求解应力的微分方程 由物理方程(2-29)式，可得
(1 ) r ( r ) E d d r d 1 dr E dr dr
2.3 厚壁圆筒应力分析
2.3 厚壁圆筒应力分析
合成氨反应器 结构示意图
2.3 厚壁圆筒应力分析
合成氨高压反应器制造安装
2.3 厚壁圆筒应力分析
(1) 薄壁容器力学分析模型
(2) 厚壁容器力学分析模型
针对厚壁筒的应力求解，将在平衡方程、几何方 程、物理方程三个方面进行分析。
2.3 厚壁圆筒应力分析
由上面方程组可导出下列“二阶齐次变系数微分方程”
d r 3 d r 0 2 dr r dr
第三节-厚壁圆筒应力分析
厚壁圆筒应力分析3.3.1弹性应力 3.3.2弹塑性应力3.3.3屈服压力和爆破压力33.4提高屈服承载能力的措施3.3.1弹性应力 3.3.2弹塑性应力一、弹塑性应力描述弹塑性疗壁圆筒的儿何与载荷参数：尺，/>; RJ;陽P () 本小节的U 的：求弹性区和塑性区里的应力假设：a.理想弹塑性材料b.圆筒体只取远离边缘区第三节 厚壁圆筒应力分析内压t 塑性区t2-22处于弹塑性状态的厚壁圆筒图2-23理想弹•塑性材料的应力■应变关系1、塑性区应力平衡方程：刃-旦drMises屈服失效判据:CF e-丐=—=丁2联立积分，得<T r=-^trJnr+Ar = &：6=-Pi内壁边界条件，求出A后带回上式得将r = R e: cr r= -p c代入(2-42)得2 ! R<p(=--a s ln-+Pl结论：① b = pjbj②q, cr^=/(lnr) rt,③cr：=-(b「+ b&) H const (区别:弹区cr. =-© + b&) =const )2 2弹性区内壁处于屈服状态:（刃）Y一（6）“ =眉$Kc=Ro/Rc(2-46)(2-26) (2-40) (2-41)将(2-42)带入(2-40)得(2-42 )(2-43)(2-44 )(2-45 ) 山表2J拉美公式得出:与2-45联立导出弹性区与塑性区交界面的pi与Rc的关系Pi =由（2-34）式（以代代替门）得若按屈雷斯卡（H.Tresca）屈服失效判据，也可导岀类似的上述各表达式。
各种应力表达式列于表2-4中结论：② 6 a d=f(r) rT->(r z. T,与「无关二、残余应力肖厚壁圆筒进入弹塑性状态后卸除内爪力pi —残余应力思考：残余应力是如何产生的卸载定理：卸载时应力改变量Ab = b-b和应变的改变量△£ = £-£之间存在着弹性关系= 图2・24。
厚壁圆筒应力分析
厚壁圆筒应力分析1、概述K>1.2的壳体成为厚壁圆筒。
厚壁容器承压的应力特点有（此处不考虑热应力）：一、不能忽略径向应力，应做三向应力分析；二、厚壁容器的应力在厚度方向不是均匀分布，而是应力梯度。
所以，在求解的时候需要联立几何方程、物理方程、平衡方程才能确定厚壁各点的应力大小。
2、解析解一、内压为i p ，外压为0p 的厚壁圆筒，需要求出径向应力r σ、周向应力θσ和轴向应力z σ，其中轴向应力z σ不随半径r 变化。
（1）几何方程如图所示，取内半径r ，增量为dr 的一段区域两条弧边的径向位移为ω和ωωd +，其应变的表达式为：r rd rd d r dr d dr d r ωθθθωεωωωωεθ=-+==-+=))(（周向应力：径向应力：（1） θσ对r 求导，得：()θθσσωωωωωσ-=⎪⎭⎫ ⎝⎛-=-='⎪⎭⎫ ⎝⎛=r r r dr d r r r dr d r dr d 112 （2） （2）物理方程根据胡克定理表示为:[]z Eσσμσεθθ+-=r (1 （3） 两式相减，消去z σ得：[]θθσσμεε-+=r E）（1-r []z r E σσμσεθ+-=(1r （4） 将（4）代入（2）得：[]）z r Edr d σσμσεθθ+-=(1 （5） 对（3）的θε求导得，z σ看做常数：⎪⎭⎫ ⎝⎛-=dr r d dr d E dr d σμσεθθ1 （6） 联立（5）、（6）得：[]θθθσσμσμσ-)1-r rdr d dr d +=（ （7） （3）平衡方程如图所示，沿径向和垂直径向建立坐标系，把θσ向x 轴和y 轴分解，得：⎪⎭⎫ ⎝⎛=-+2sin 2θθd p p p r dr r （8） 其中()θσσd dr r d p r r dr r ++=+)（ （9）θσrd p r r =由于θd 很小，22sin θθd d ≈⎪⎭⎫⎝⎛，略去二阶微量r r d d σ，得 drd r r r σσσθ=- （10） 联立（7）（10）得0322=+drd dr d r r r σσ （11）对（11）进行求得r σ，在代入（10）得22r B A r BA r +=-=θσσ（12） 其中A 、B 是两个积分常数，要求A ，B 需要两个方程，根据内外壁边界条件00,,p R r p R r r ir i -==-==σσ（13）将（13）代入（12）得：22020202202002)(i ii i i i R R RR p p B R R R p R p A --=--=（14）最后剩下z σ未求出，最后在轴向用平衡方程，内力等于外力。
厚壁圆筒应力分析剖析
厚壁圆筒应力分析剖析厚壁圆筒是一种常见的结构，广泛应用于各个领域，比如压力容器、热交换器等。
在使用厚壁圆筒的过程中，必须进行应力分析，以确保结构的安全性和可靠性。
首先，研究厚壁圆筒的应力分析需要考虑以下几个方面。
1.圆筒的几何形状：厚壁圆筒是由外径、厚度和长度组成的。
这些几何参数会影响圆筒内部的应力分布情况。
2.材料特性：圆筒的材料特性直接影响其应力分布。
研究厚壁圆筒时，通常会考虑材料的弹性模量和泊松比等参数。
3.加载条件：圆筒的应力分布受外部载荷的影响。
载荷的形式可以是压力、温度、重力等。
加载条件的确定对于应力分析至关重要。
接下来，我们将详细介绍厚壁圆筒的应力分析方法。
1.内外压力分析：考虑厚壁圆筒内外的压力差异。
当内外压力相等时，圆筒应力较小。
当内压大于外压时，圆筒将会受到较大的应力。
2.纵向应力分析：厚壁圆筒在纵向方向上承受的应力主要为轴向拉应力。
如果存在压力差，则拉应力沿厚度逐渐增加。
3.周向应力分析：在周向上，厚壁圆筒受到的应力主要为周向拉应力。
当圆筒内外压力不平衡时，周向应力将会增加。
4.切应力分析：切应力是圆筒内部的剪切应力分量。
在圆筒壁厚度的不同位置，切应力的大小也会有所不同。
5.应力分布图：为了更好地理解厚壁圆筒的应力分布情况，可以绘制应力分布图。
这样可以直观地了解不同部位的应力分布情况，以便进行结构优化。
总结一下，厚壁圆筒的应力分析对于确保结构安全性至关重要。
通过分析内外压力、纵向应力、周向应力和切应力，可以更好地理解圆筒的应力分布情况。
通过应力分布图，可以更直观地了解圆筒不同部位的应力情况，从而进行优化设计。
在实际工程中，应力分析的结果可以用来指导材料的选择、结构的设计以及使用中的安全操作。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
厚壁圆筒应力分析
1、概述
K>1.2的壳体成为厚壁圆筒。
厚壁容器承压的应力特点有（此处不考虑热应力）：一、不能忽略径向应力，应做三向应力分析；二、厚壁容器的应力在厚度方向不是均匀分布，而是应力梯度。
所以，在求解的时候需要联立几何方程、物理方程、平衡方程才能确定厚壁各点的应力大小。
2、解析解
一、内压为i p ，外压为0p 的厚壁圆筒，需要求出径向应力r σ、周向应力θσ和轴向应力z σ，其中轴向应力z σ不随半径r 变化。
（1）几何方程
如图所示，取内半径r ，增量为dr 的一段区域两条弧边的径向位移为ω和ωωd +，其应变的表达式为：
r rd rd d r dr
d dr d r ω
θθθωεω
ωωωεθ=
-+==
-+=
))(（周向应力：径向应力：（1）
θσ对r 求导，得：
()θθσσωωωω
ωσ-=⎪⎭⎫
⎝⎛-=-='⎪⎭⎫ ⎝⎛=r r
r dr d r r r dr d r dr d 112 （2） （2）物理方程 根据胡克定理表示为
[]z E
σσμσεθθ+-=
r (1
（3） 两式相减，消去z σ得：
[]θθσσμεε-+=
r E ）
（1-r []z r E
σσμσεθ+-=(1r
（4） 将（4）代入（2）得：
[]）
z r E
dr d σσμσεθθ+-=(1
（5） 对（3）的θε求导得，z σ看做常数：
⎪⎭
⎫
⎝⎛-=dr r d dr d E dr d σμσεθθ1 （6） 联立（5）、（6）得：
[]θθθσσμσμσ-)
1-r r
dr d dr d +=
（ （7） （3）平衡方程
如图所示，沿径向和垂直径向建立坐标 系，把θσ向x 轴和y 轴分解，得：
⎪⎭
⎫ ⎝⎛=-+2sin 2θθd p p p r dr r （8）
其中
()θσσd dr r d p r r dr r ++=+)（ （9）
θσrd p r r =
由于θd 很小，2
2sin θ
θd d ≈⎪⎭⎫
⎝⎛，略去二阶微量r r d d σ，得 dr
d r
r
r σσσθ=- （10） 联立（7）（10）得
0322=+dr
d dr d r r r σσ （11）
对（11）进行求得r σ，在代入（10）得
2
2
r
B
A r
B A r +=-
=θσσ （12） 其中A 、B 是两个积分常数，要求A ，B 需要两个方程，根据内外壁边界条件
0,,p R r p R r r i r i -==-==σσ （13）
将（13）代入（12）得：
22020
202202
002)(i
i i i i i R
R R R p p B R R R p R p A --=
--=
（14）
最后剩下z σ未求出，最后在轴向用平衡方程，内力等于外力。
)S 220z z i R R F -==（内力πσσ （15） 0202i R p R p F i ππ-=外力 （16）
联立（15）（16）得：
A R R R p R p R R p R p R i
i i i i i z =--=--=2
2020022200202)(πππσ （17） 所以得到内外压作用下的拉美公式：
周向应力：2
2202
020********
)(r R R R R p p R R R p R p i i i i i i --+--=θσ （18） 径向应力：222
0202022020021
)(r R R R R p p R R R p R p i i i i i i r -----=σ 轴向应力：2
2020
02i
i i z R R R p R p --=σ
利用拉美公式，可以得到如下解析解：
设厚壁圆筒受内压Mpa p i 4=，受外压Mpa p o 1=，圆筒内径
mm i 120R =，外径mm o 160R =。
则得到如下解析解
3、数值解
有限元模型及其假设
并且得到总的应力分布：
由于筒体为对称结构，所以用轴对称模型。
同样设 设厚壁圆筒受内压Mpa p i 4=，受外压Mpa p o 1=，
圆筒内径mm i 120R =，壁厚t 为40mm 高度H 为200mm ，弹性模量E 为Mpa 1025⨯，泊松比为10Mpa 。
模型如右图所示： 得到，厚壁圆筒只在受内压和外压作用时的三向应力，在ANSYS 中得到如下图所示：
实体模型及网格划分图
和拉美公式的到的解析解进行比较，下面是利用Excel整理得到的厚壁圆筒三向应力有限元数值解和拉美公式解析解比较图，如下图所示：
从图中看出，ANSYS的数值解和拉美公式解析解几乎一样，进一步验证了拉美公式的正确性。