厚壁圆筒或管道中的应力计算
注塑机筒疲劳强度计算的设计准则
2023年 第49卷·71·作者简介：袁卫明（1971-），男，本科，正高级工程师，副总工程师，现从事塑料注射成型机研发设计。
收稿日期：2023-07-200 引言注塑机机筒是注射机构中的重要零部件，在工作中其要承载注射高压的冲击，当前注塑机的注射压力已从传统的170 MPa 发展到270 MPa 以上。
面对机筒在高压和超高压中出现的失效现象，沿用传统的注塑机筒强度理论[1]，不能圆满解释机筒失效的实际现象。
本文从厚壁圆筒的弹塑性力学理论[2]分析研究注塑机筒的工作特性，阐述以往用弹性力学角度分析研究机筒强度的局限性，提出了符合实际的注塑机筒疲劳强度的设计准则，并用实例加以论证。
1 厚壁圆筒1.1 厚壁圆筒的应力分析根据厚壁圆筒体[3]的应力变形特点，我们假设将厚壁圆筒看成是由许多个薄壁圆筒相互连在一起所组成，如图1所示，当厚壁圆筒内径承受内压力后，其组成的各层薄壁圆筒由里至外逐步受力，其变形受到里层薄壁圆筒的约束和受到外层薄壁圆筒的限制，因此各个单元薄壁圆筒体都会受到内外侧变形的约束和限制所引起的均布压力作用，从里往外各层薄壁圆筒体的变形被受到的约束和限制是不同的，环向应力沿壁厚方向分布是不均匀的，这是厚壁圆筒形变和应力的一个基本特点。
厚壁圆筒应力、应变的另一个特点是：由于厚壁圆筒是由多个薄壁圆筒组成，在多层材料变形的相互约束和限制下，沿径向方向产生了径向应力，沿壁厚方向径向应力分布是不均匀的。
厚壁圆筒和薄壁圆筒注塑机筒疲劳强度计算的设计准则袁卫明，成明祥(德清申达机器制造有限公司，浙江 湖州 313205)摘要：传统注塑机筒强度设计理论未能合理解释回答在实际中产生的一些失效现象问题，对比厚壁圆筒的力学分析，确认判断注塑机筒沿用以往的设计理论具有局限性和适用范围。
通过引用分析目前在厚壁圆筒中较常用的弹塑性强度理论设计观点，结合实例，提出了符合实际的注塑机筒强度理论的设计准则。
平面问题作业--厚壁圆筒应力分析
图 厚壁圆筒问题
问题描述及要求
如图所示为一厚壁圆筒，其内半径r 1=50 mm ，外半径r 2=100 mm ，作用在内孔上的压力p=10 MPa ，无轴向压力，轴向长度很大可视为无穷。
材料参数：2e11（弹性模量），泊松比：0.3；计算厚壁圆筒的径向应力σr 和切向应力σt 沿半径r 方向的分布。
根据材料力学的知识，σr 、σt 沿r 方向的分布的解析解为
⎪
⎪⎭
⎫
⎝⎛--=2222
1
2221r 1r r r r p
r σ ⎪⎪⎭
⎫ ⎝⎛+-=222212221t 1r r r r p r σ
提示：该问题符合平面应变问题的条件，故可以简化平面应变问题进行分析。
另外，根据对称性，可取圆筒的四分之一并施加垂直于对称面的约束进行分析。
利用路径操作。
（1）
步骤：
1、定义单元类型
Ok
options
2、定义材料属性
3、创建模型
4、划分单元
Size controls--lines--set
apply
拾取圆弧边输入20
mesh
5、施加约束
apply
拾取左边线
6、施加载荷
7、求解
8、显示单元
Plot--elements
9、定义路径
顺次拾取下边线结点
Plot paths Map onto path
10、作路线图
11、结果。
带裂纹厚壁圆筒应力强度因子的几种计算方法
第 ２ ６卷 第 ４期
２０ ０ ２年 ８月
南
京 理
工
大
学
学
报
、
Ｖｏ ． ６ １２ Ｎｏ ４ ．
Ａｕ ｇ．２ ２ ００
ｖ ｒ ｉｙ ｏｆＳｃ ｅ ｅ ａｎ Ｔｅ ｈ ｏｙ Ｊ ｕ ｎ ｌｏ ｎ ｉ Ｕｎｉ ｅ ｓ ｔ ｉｎｃ ｄ ｃ ｎｏｌ ｇ ｏ ｒ ａ ｆＮａ ｊｎ
４ ３１
式中， ： 志 【 ３—４
｝ 平 应 本 文 以后 的公 式 都 是 按 平 面应 变 公 式 而 得 到 的 。 面力
，
平 面 应 变
图 １ 带 裂 纹厚 壁 筒 模 型 图
Ｆｉ １ Ｍ ｏ ｅ ｆ ｔ ｉ ｋ ｗａ Ｉ ｇ． ｄ Ｉ ｈｃ ｌ ｏ ｃ ｌ ｄ ｒ ｗｉｈ ｃ ａ ｋ ｙ ｉ ｅ ｔ ｒ ｃ ｎ
收稿 日期 ：０ １ ３—３ ２ ０ —０ ０
国防 科技 预研 行业 基 金资 助项 目 陈爱 军 男 ３ ０岁 博 士生
（ １ ）
维普资讯
总第 １ ５期 ２
陈爱 军
徐
诚
胡 小秋
带 裂纹 厚 壁 圆筒应 力 强度 因子 的几种 计算 方 法
关键词
分类号
应 力 强 度 因子 ， 限元 分 析 ， 壁 壳体 有 厚
Ｏ ４ ． ，Ｔ ０ ３ ６ １ Ｊ３ １
承 受 内压 的厚 壁 筒 在 工 程 中广 泛 使 用 ， ： 业 管 道 、 管 炮 管 等 等 。 壁 筒 在 工 作 时 ， 如 工 枪 厚
由 于 高 压 的循 环 作 用 ， 内壁 常 产 生 裂 纹 ， 续 工 作 ， 由于 裂 纹 的扩 展 而 导 致 破 坏 ， 了工 作 继 会 为
厚壁圆筒应力分析
温度变化引起的弹性热应力
热应力
构件之间热变形 的相互约束
构件热变形受到 外界约束
构件内部温度 分布不均匀
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。
压力容器应力分析_厚壁圆筒弹性应力分析
工程上一般将设计压力在10≤p≤100MPa之间的压力容器称为高压容器，而将100MPa压力以上的称为超高压容器。
高压容器不但压力高，而且同时伴有高温，例如合成氨就是在15～32MPa压力和500℃高温下进行合成反应。
一般来说，高压和超高压容器的径比K > 1.2，称此类容器为“厚壁容器”。
本章讨论的对象，是厚壁圆筒型容器。
承受压力载荷或者温差载荷的厚壁圆筒容器，其上任意点的应力，是三向应力状态。
即存在经向应力（又称轴向应力）、周向应力和径向应力。
针对厚壁筒的应力求解，将在平衡方程、几何方程、物理方程三个方面进行分析。
2.2.1 弹性应力－压力载荷引起的弹性应力（1）轴向(经向)应力ϭz222200002200002220()1i z i i i i i i i z i iP P FP P p R p R F R R p R p R p p KR K R R K R σππππσ−=−=⋅−⋅=−−−⋅===−−径比(2) 周向应力ϭ和径向应力ϭrθ三对截面：一对圆柱面，相距dr一对纵截面，相差dθ一对横截面，长度为1Ϭz作用在横截面上Ϭr作用在圆柱面上Ϭθ作用在纵截面上平衡方程（沿径向列平衡方程）()()112sin 102r r r d d r dr d rd dr θθσσθσθσ++⋅−⋅−⋅=sin 22d d θθ≈略去高阶无穷小，并使得到平衡方程r r d r drθσσσ−=几何方程()r w dw wdwdr drε+−==径向应变周向应变()r w d rd wrd r θθθεθ+−==上述表达式是Lame 在1833年推得的，又称为Lame 公式。
当仅有内压时，p o =0，有()222222211111112i o i o r z i z r p R K r p R K r p K θθσσσσσσ⎧⎛⎞=⋅−⎪⎜⎟−⎝⎠⎪⎪⎛⎞⎪=⋅+⎜⎟⎪−⎝⎠⎨⎪⎛⎞=⋅⎪⎜⎟−⎝⎠⎪⎪=+⎪⎩246810010********σθ R i / σθ R oK可见，当K 越大时，应力的分布就越不均匀。
厚壁圆筒应力分析
厚壁圆筒应力分析1、概述K>1.2的壳体成为厚壁圆筒。
厚壁容器承压的应力特点有（此处不考虑热应力）：一、不能忽略径向应力，应做三向应力分析；二、厚壁容器的应力在厚度方向不是均匀分布，而是应力梯度。
所以，在求解的时候需要联立几何方程、物理方程、平衡方程才能确定厚壁各点的应力大小。
2、解析解一、内压为i p ，外压为0p 的厚壁圆筒，需要求出径向应力r σ、周向应力θσ和轴向应力z σ，其中轴向应力z σ不随半径r 变化。
（1）几何方程如图所示，取内半径r ，增量为dr 的一段区域两条弧边的径向位移为ω和ωωd +，其应变的表达式为：r rd rd d r drd dr d r ωθθθωεωωωωεθ=-+==-+=))(（周向应力：径向应力：（1）θσ对r 求导，得：()θθσσωωωωωσ-=⎪⎭⎫⎝⎛-=-='⎪⎭⎫ ⎝⎛=r rr 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 rrr σσσθ=- （10） 联立（7）（10）得0322=+drd dr d r r r σσ （11）对（11）进行求得r σ，在代入（10）得22rBA rB A r +=-=θσσ （12） 其中A 、B 是两个积分常数，要求A ，B 需要两个方程，根据内外壁边界条件0,,p R r p R r r i r i -==-==σσ （13）将（13）代入（12）得：22020202202002)(ii i i i i RR R R p p B R R R p R p A --=--=（14）最后剩下z σ未求出，最后在轴向用平衡方程，内力等于外力。
厚壁圆筒应力分析
厚壁圆筒应力分析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.只受内压作用：（1）在厚壁圆筒中，筒体处于三向应力状态，其中θσ为拉应力，r σ为压应力，且沿壁厚非均匀分布；而z σ介于θσ和r σ之间，即2r z θσσσ+=，且沿壁厚均匀分布。
（2）在筒体内壁面处，θσ、r σ的绝对值比外壁面处为大，其中θσ具有最大值，且恒大于内压力i p ，其危险点将首先在内壁面上产生。
（3）θσ沿壁厚分布随径比K 值的增加趋向更不均匀，不均匀度为内、外壁周向应力之比，即2()1()2io r R r R K θθσσ==+=。
显然，不均匀度随2K 成比例，可见K 值愈大，应力分布愈不均匀。
当内壁材料开始屈服时，外壁材料远小于屈服限，因此筒体材料的强度不能得到充分的利用。
由此可知，用增加筒体壁厚（即增加K 值）的方法来降低厚壁圆筒的内壁应力，只在一定范围内有效，而内压力接近或超过材料的许用应力时，增加厚度是完全无效的。
为了提高筒壁材料的利用率，有效的办法是改变应力沿壁厚分布的不均匀性，使其趋于均化。
2.往往采用组合圆筒或单层厚壁圆筒自增强处理技术，以提高筒体的弹性承载能力。
3.温差应力：厚壁圆筒的厚壁可能从内表面或外表面被加热，由于筒壁较厚，并有一定的热阻，在筒体的内、外壁之间存在温度差，温度较高部分因受热而引起膨胀变形，同时受到温度较低部分的约束，从而使前者受压缩，而后者受拉伸，出现了温差应力或称热应力。
（1）厚壁圆筒中，温差应力与温度差t ∆成正比，而与温度本身的绝对值无关，因此在圆筒内壁或外壁进行保温以减小内、外壁的温度差，可以降低厚壁圆筒的温差应力。
（2）温差应力的分布规律为三向应力沿壁厚均为非均匀分布，其中，轴向应力是环（周）向应力与径向应力之和，即t t t z r θσσσ=+ ；在内、外壁面处，径向应力为零，轴向应力和环（周）向应力分别相等，且最大应力发生在外壁面处。
（3）温差应力是由于各部分变形相互约束而产生的，因此应力达到屈服极限而屈服时，温差应力不但不会继续增加，而且在很大程度上会得到缓和，这就是温差应力的自限性，它属于二次应力。
厚壁圆筒问题ansys求解操作
厚壁圆筒问题（平面应变问题）
[本例提示] 介绍了利用 ANSYS 进行静力学分析的方法、步骤和过程，并对将空间问题简
化为平面问题的条件、方法进行了简单的介绍。
概述
平面问题 所谓平面问题指的是弹性力学的平面应力问题和平面应变问题。 当结构为均匀薄板，作用在板上的所有面力和体力的方向均平行于板面，而且不沿厚度 方向发生变化时，可以近似认为只有平行于板面的三个应力分量 σx、σy 、τxy 不为零，所以这 种问题就被称为平面应力问题。 设有无限长的柱状体，在柱状体上作用的面力和体力的方向与横截面平行，而且不沿长 度长度发生变化。此时，可以近似认为只有平行于横截面的三个应变分量 ε x、ε y 、γ xy 不为零， 所以这种问题就被称为平面应变问题。 对称性 当结构具有对称面而载荷也对称于该对称面时，可以利用该对称性，取结构的一半进行 分析，并且约束掉对称面上垂直方向的位移，从而减少了计算工作量。
2
厚壁圆筒问题 ansys 求解
图 7-4
单元类型库对话框
图 7-5
单元选项对话框
图 7-6
材料模型对话框
钮，然后关闭图 7-6 所示的对话框。 （注意：从解析公式中可以看出，径向应力 σr 和切向应力
3
厚壁圆筒问题 ansys 求解
图 7-7
材料特性对话框
图 7-8
创建面对话框
σt 与弹性模量、泊松比无关，但是，这两个参数在有限元分析中却是必须的。 ） 7.3.4 创建实体模型 拾取菜单 Main Menu→Preprocessor →Modeling→Create→ Areas →Circle→ By Dimensions 。 弹出的图 7-8 所示的对话框，在“RAD1” 、 “RAD2” 、 “T HETA2”文本框中分别输入 0.1、0.05 和 90，单击“Ok” 按钮。 7.3.5 划分单元 拾取菜单 Main Menu→Preprocessor →Meshing→MeshTool。弹出的图 7-9 所示的对话框，
第三章第四节2--厚壁圆筒-应力
bardcdr dr r
微单元体
r dr
b
c
a
d 2
dr
d
r
r
d 2
厚壁圆筒
图3-17 厚壁圆筒微元体受力情况
在圆筒体半径为r处，以相距dr的二环向截面及夹角 d
的二径向截面截取任一微元体，其微元体在轴向的长度为1。
由于轴向应力对径向应力的平衡没有影响，所以图中未标出
轴向应力。
根据半径r方向力的平衡条件，有：
d 1-32KK 22-p1
相应对载荷的限制为： 或
p
K2 -1
2K 2
pmax
K2 -1
2K2
当 K 时 p m ， a0 x .5 ，其含义是，
对厚壁圆筒，其壁厚的无限增加只能换来允许承受载 荷的有限增加。即用增加壁厚来增大承载能力是有限 和有条件的。在应力低的筒体外壁处增大壁厚，对筒 体提高承载能力作用不大，甚至造成浪费或其他问题。
比较厚壁圆筒应力计算公式与薄壁圆筒壳应力计算公式，对 了解圆筒壳应力计算公式的精确度和适用范围是十分有益的。 以环向应力为例，圆筒壳环向薄膜应力为：
p R p 2((0 0 R R R R ii))2K( -1 1 Kp )
式中，R为圆筒壳平均半径。 若以厚壁圆筒应力公式进行计算，其最大环向应力为：
多层板厚壁筒体及绕带筒体的采用，可以有效 地避开单层厚壁筒体的上述局限性。
（二）、根据弹性失效准则，厚壁圆筒的承压能力是根据内壁的
强度条件决定的
承内压厚壁圆筒的应力最大部位是在内壁壁面处，根据工
程上常用的弹性失效准则，应力最大部位的应力强度达到极限值
时，结构即失去了承载能力。因而，按第三强度理论建立的内壁
第二章_厚壁圆筒的弹塑性应力分析1
c3
2
E
Ri2 pi Ro2
Ro2 po Ri2
c1
1
E
c3
1
E
Ri2 pi Ro2
Ro2 po Ri2
c2
1
E c4
1
E
Ri2 Ro2 ( pi po ) Ro2 Ri2
（c）
12/2/2019
厚壁圆筒的弹塑性应力分析
d d rd d rr 1 r(r) (2-10)
12/2/2019
厚壁圆筒的弹塑性应力分析
Page - 17
二、厚壁圆筒的应力和位移解
本节采用位移法求解在均匀内、外压作用下的厚壁
圆筒。将几何方程式代入物理方程式，得出用位移 分量表示的物理方程
r
E
du (
1 dr
（2-20）
12/2/2019
厚壁圆筒的弹塑性应力分析
Page - 28
（3）两端封闭同时受轴向刚性约束的筒体（高压管 道或厚壁圆筒无限长）
轴向变形受到约束，
z 0
z 2C32Ri2R po2 i R Ro i2 2po
C1
(12E)(1)C3
(12)(1)
12/2/2019
厚壁圆筒的弹塑性应力分析
Page - 32
z 2c3 E Z
c3
E 1
Z c1 1 2
c4
1
E
c2
厚壁圆筒的弹塑性应力分析
（2-14）
厚壁圆筒的弹塑性分析
厚壁圆筒的弹塑性分析弹塑性分析是一种结构分析方法，适用于材料在一定强度范围内既具有弹性行为又具有塑性行为的情况。
厚壁圆筒是一种常见的结构，广泛应用于工程中，如汽车零部件、压力容器等。
本文将介绍厚壁圆筒的弹塑性分析方法，并结合一个具体的例子进行说明。
厚壁圆筒的弹性分析是指在圆筒内外受到压力作用时圆筒的变形和应力分布的计算。
在弹性阶段，材料的应力-应变关系是线性的，可以通过胡克定律描述。
在塑性阶段，材料的应力-应变关系是非线性的，需要采用本构关系来描述。
首先，我们来介绍圆筒的几何参数。
厚壁圆筒可以由内外半径分别为R1和R2的圆柱体围成，圆柱体的高度为h。
此外，圆筒的材料有一个屈服强度σy，用于描述材料的塑性行为。
对于厚壁圆筒，弹性阶段的计算相对简单。
在内外压力P的作用下，圆筒的应变可以通过应力与材料的弹性模量E之间的关系得到。
圆筒的轴向应变εr可以通过胡克定律得到：εr=σr/E其中，σr是圆筒轴向应力，E是材料的弹性模量。
圆筒的周向应变、轴向切变应变可以根据几何关系得到。
在弹性阶段，应力满足柯西－格林弹性方程：σr=λ(εr+εθ)+2μεrσθ=λ(εr+εθ)+2μεθτrz = μ(εr - εθ)其中，λ和μ是材料的拉梅常数，可以通过杨氏模量E和泊松比ν计算得到。
当圆筒的应力达到屈服强度σy时，就进入了塑性阶段。
在塑性阶段，应力与应变之间的关系通过本构关系来描述。
常用的本构关系包括线性硬化本构关系、塑性截面变形本构关系等。
本文以线性硬化本构关系为例进行说明。
线性硬化本构关系假设材料的塑性应变是线性增加的。
圆筒中心的塑性应力σp和塑性应变εp可以通过以下方程计算：σp=σyεp=(σr-σy)/E*H其中，E*是圆筒在弹性阶段的等效弹性模量，H是圆筒的等效刚度。
对于给定的压力P，可以通过迭代法来确定圆筒的应力和应变分布。
首先假设圆筒是在弹性阶段，在初始状态下计算应力和应变分布。
然后，通过本构关系计算塑性应力和塑性应变分布。
厚壁圆筒__弹塑性力学知识
arr
r rb
内表面产生压缩的切向残余应力，当再次加载时，产生 的切向应力被抵消一部分，可提高圆筒的弹性极限压力。
五、位移分量（平面应变状态）
(1 )a 2 p b 2 1. 弹性阶段： u (1 2 )r 2 2 E (b a ) r
e s ij s ij s ij r
sr
r 2s s b 2
2b
2
2 1 r
ρ rb
sq
r 2s s b 2
2b 2
2 1 r
r r s r s s ln p* a r ρ sq s s 1 l n p * a a a2 p * b2 a2 p* b2 1 2 1 2 s re 2 sq e 2 2 2 b a r b a r r b2 p a2 1 2 a r r s s ln 2 2 a b a r r sr s sr 2 a 2 p b2 2b 2 b 2 a 2 1 r 2 r r b
r: 弹塑性分界面的半径。
ss ds r dr r
s r s s lnr C
C p s s lna
塑性区：a r r
屈服条件： sq – sr= ss
r s r s s ln p a r sq s s 1 ln p a
塑性区的应力分量是静定的。
边界条件：
sr sr
r a p1 r b
p2
一、弹性分析
2.3_厚壁圆筒应力分析
轴向应力
pi R p0 R z A 2 R Ri
2 i 2 0 2 0
称Lamè （拉美）公式
（2-34）
14
2.3 厚壁圆筒应力分析
2.3.1 弹性应力（续）
过程设备设计
当仅有内压或外压作用时，拉美公式可以简化，此时，厚壁圆筒 应力值和应力分布分别如表2－1和图2－17 表2-1 厚壁圆筒的筒壁应力值
to
厚壁圆筒各处的热应力见表2-2， Et Pt 表中 21 厚壁圆筒中热应力分布如图2-20所示。
28
2.3 厚壁圆筒应力分析
2.3.1 弹性应力（续） 表2-2 厚壁圆筒中的热应力
过程设备设计
热应 力
rt
r 任意半径 处
圆筒内壁K r K 处 0
圆筒外壁 K r 1 处 0
＝A
（2-25）
6
2.3 厚壁圆筒应力分析
2.3.1 弹性应力（续）
过程设备设计
2. 周向应力与径向应力 由于应力分布的不均匀性，进行应力分析时，必须从微元体着 手，分析其应力和变形及它们之间的相互关系。 a. 微元体
b. 平衡方程
c. 几何方程 ：微元体位移与应变之间的关系。（用位移法求解） d. 物理方程：弹性范围内，微元体的应变与应力的关系 e. 平衡、几何和物理方程综合—求解应力的微分方程 （求解微分方程，积分，边界条件定常数）
z
1 pi 2 K 1
K2 po 2 K 1
15
2.3 厚壁圆筒应力分析
2.3.1 弹性应力（续）
过程设备设计
z
K2 1 m ax pi K2 1
m in pi
压力容器应力分析-厚壁圆筒应力分析
•2.3 厚壁圆筒应力分析
•③除 外，其它应力沿壁厚的不均匀程度与径比K值有关。 • 以 为例，外壁与内壁处的 • 周向应力 之比为： • K值愈大不均匀程度愈严重， • 当内壁材料开始出现屈服时， 外壁材料则没有达到屈服， • 因此筒体材料强度不能得到充分的利用。
•2.3 厚壁圆筒应力分析
•二、温度变化引起的弹性热应力
•一、压力载荷引起的弹性应力 • 1、轴向（经向）应力
•对两端封闭的圆筒，横截面在变形后仍保持平面。所 以，假设轴向应力沿壁厚方向均匀分布，得：
•＝ A •（2-25）
•2.3 厚壁圆筒应力分析
• 2、周向应力与径向应力 •由于应力分布的不均匀性，进行应力分析时，必须从微元体 着手，分析其应力和变形及它们之间的相互关系。
•厚壁圆筒中热应力及其分布的规律为： •① 热应力大小与内外壁温差成正比
• 取决于壁厚，径比K值愈大 值也愈大，表2-2中的
•
值也愈大。
•②热应力沿壁厚方向是变化的
•2.3 厚壁圆筒应力分析
•3、内压与温差同时作用引起的弹性应力
•（2-39 ）
•具体计算公式见表2-3，分布情况见图2-21。
•2.3 厚壁圆筒应力分析
•a. 微元体
•b. 平衡方 程 •c. 几何方程 (位移－应变）
•d. 物理方程（应变－应力）
•e. 平衡、几何和物理方程综合—求解应力的微分方程 • （求解微分方程，积分，边界条件定常数）
•2.3 厚壁圆筒应力分析
•a. 微元体 •如图2-15(c)、(d)所示，由圆柱面mn、m1n1和纵截面mm1、nn1 组成，微元在轴线方向的长度为1单位。 •b. 平衡方程
•（2-26）
•2.3 厚壁圆筒应力分析
厚壁圆筒模型在均布应力作用下位移计算
中国建筑标准知识中国建筑标准知识，那可是相当有趣又实用的东西呢。
咱先说说这建筑标准为啥重要。
你看，盖房子就跟做衣服似的。
做衣服你得有个尺寸标准吧，不能袖子长到拖地，裤腿短得像七分裤还不对称。
盖房子也是一样，要是没有标准，这房子可能歪歪扭扭的，就像个喝醉酒的大汉，站都站不稳。
房间的大小、门窗的位置，都得按照标准来，这样住着才舒服、安全。
比如说，卧室要是太小了，就像把人塞在一个小盒子里，那多难受啊；客厅太大了，又感觉空落落的，像个大仓库。
建筑标准就是那个给房子“量体裁衣”的尺子。
咱再聊聊建筑材料的标准。
建筑材料就好比是房子的肌肉和骨骼。
木材得是干燥的、结实的，要是用了那种湿哒哒、软绵绵的木材，就像用棉花糖来搭房子，风一吹就散架了。
钢材呢，得符合一定的强度标准，不然就像用软面条来做支撑，那房子还不得塌成一堆废墟啊。
水泥也有它的标准，标号合适的水泥才能让房子牢固得像一座小山，要是标号不对，就像用泥巴糊墙，根本经不住风吹雨打。
房屋的采光标准也很有讲究。
房子里要是没有足够的采光，就像住在山洞里一样，黑乎乎的，人的心情都会变得压抑。
采光好的房子，白天阳光能洒进来，那感觉就像给房子注入了活力，房间里亮堂堂的，人在里面住着也觉得神清气爽。
窗户的大小、朝向，这些都是根据建筑标准来确定的，可不是随便开个洞就行的。
这就好比人的眼睛，窗户的位置和大小决定了房子能接收到多少光线，就像眼睛的大小和睁开的程度决定了人能看到多少东西一样。
还有建筑的防火标准呢。
这可是关乎生命财产安全的大事。
如果房子没有按照防火标准来建造，一旦着火，那可就像干柴堆里丢了个火把，一下子就烧得精光。
防火门得能真正起到防火的作用，不能是那种徒有其表的“纸老虎”。
防火墙也得足够坚固，能够阻挡火势的蔓延。
这就像给房子穿上了一层防火的铠甲，让住在里面的人有足够的时间逃生。
建筑的抗震标准更是不能忽视。
咱们国家有些地方地震比较频繁，房子要是不抗震，那就像纸糊的灯笼，一摇就碎。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
厚壁圆筒或管道中的应力计算
（1）概述
当厚壁管或圆柱体受到内部和外部压力时，会在壁中产生环箍和纵向应力。
（2）轴向方向应力
σa = (p i r i2 - p o r o2 )/(r o2 - r i2) (1)
σa=轴向应力（MPa，psi）
pi=管道或圆柱体中的内部压力（MPa，psi）
p o=管道或圆柱体中的外部压力（MPa，psi）
r i=管道或圆柱体的内径（mm，in）
r o=管子或圆柱体的外半径（mm，in）
（3）周向应力-环向应力
圆周方向上的应力——环向应力——在管或圆筒壁上的一点上可以表示为：
σc = [(p i r i2 - p o r o2) / (r o2 - r i2)] - [r i2 r o2 (p o - p i) / (r2 (r o2 r i2))] (2)其中：
σc=周向应力（MPa，psi）
r=管道或圆筒壁中点的半径（mm，in）（r i<r<r o）
r=r i时的最大应力（管道或圆柱体内部）
(4)合成应力
气缸壁中单个点的组合应力不能通过使用矢量加法的单个矢量来描述。
相反，可以使用描述两个物理向量之间的线性连接的应力张量（矩阵）。
径向应力
管壁或圆筒壁中某一点处的径向应力可以表示为：
σr= [(p i r i2 - p o r o2) / (r o2 - r i2)] + [r i2 r o2 (p o - p i) / (r2 (r o2 - r i2))] (3) r=r o时的最大应力（管道或圆柱体外部）
(5)示例-厚壁圆筒中的应力
在内径为200mm（半径为100mm）、外径为400mm（半径为200mm）的圆柱体中，相对于外部压力存在100MPa的压力。
轴向应力可计算为:
σa=（（（100 MPa）（100 mm）2-（0 MPa）（200 mm）2）/（（200 mm =33.3 MPa
内壁（100 mm）的周向应力（环向应力）可计算为:
σc=[（（100 MPa）（100 mm）2-（0 MPa）（200 mm）2）/（200 mm
=167 MPa
内壁（100 mm）的径向应力可计算为:
σr=[（（100 MPa）（100 mm）2-（0 MPa）（200 mm）2）/（200 mm
=-100MPa。