应力成分输入
中心 C: — MPa
半径 R: — MPa
θ_p: —°
应力单元 ↔ 莫尔应力圆 实时联动
σx 面
σy 面
当前面 (σθ,τθ)
主平面
最大剪切面
理论与主要公式
$$\sigma_{1,2} = \frac{\sigma_x+\sigma_y}{2} \pm \sqrt{\left(\frac{\sigma_x-\sigma_y}{2}\right)^2 + \tau_{xy}^2}$$
主应力:\(\sigma_1 \geq \sigma_2\);莫尔圆中心 \(C=(\sigma_x+\sigma_y)/2\),半径为根号项。
$$\tau_{max} = \sqrt{\left(\frac{\sigma_x-\sigma_y}{2}\right)^2 + \tau_{xy}^2}$$
最大剪应力出现在与主应力面相差 45° 的平面上。
$$\tan 2\theta_p = \frac{2\tau_{xy}}{\sigma_x - \sigma_y}$$
主应力方向角 \(\theta_p\) 为 \(x\) 轴到主应力方向的角度。
什么是莫尔应力圆
🙋
老师,我经常在材料力学里看到“应力状态”这个词,到底是什么意思啊?
🎓
简单来说,就像你用手指去戳一块橡皮泥,橡皮泥内部某个点在不同方向上感受到的“挤压”或“拉扯”的力,这个整体情况就是应力状态。在实际工程中,比如一个桥梁的钢梁,它内部每一点都同时承受着来自不同方向的应力。你可以试着在模拟器里拖动σx、σy和τxy这三个滑块,它们分别代表x方向的压应力、y方向的压应力和剪切应力,改变它们就能看到这个点“感受”到的力在如何变化。
🙋
诶,真的吗?那“主应力”又是什么?听起来好复杂。
🎓
别怕,其实很简单。想象一下,你手里拿着一张有皱纹的纸,总可以找到两个互相垂直的特殊方向,沿着这两个方向,纸只被纯粹地拉伸或压缩,没有任何“撕扯”(剪切)。这两个方向上的正应力就是主应力。在模拟器里,你输入任意一组应力,它都会自动帮你算出这两个主应力σ₁和σ₂,并且用圆上的两个端点直观地标出来。你试试把剪切应力τxy设为零,会发现主应力方向就和坐标轴重合了。
🙋
哦!那旁边那个最大剪应力τmax又是干嘛的?这个很重要吗?
🎓
非常重要!对于很多材料,尤其是金属,破坏往往不是从拉断开始,而是从“剪坏”开始的。比如一根螺栓被拧得太紧时,内部最大的剪切应力如果超过了材料的承受极限,它就可能被剪断。τmax就是材料内部可能出现的最大剪切应力值。在模拟器的圆上,它就是圆的最高点。你试着把σx和σy设成一样,然后只改变τxy,会发现τmax就等于τxy,而且这时候材料最容易发生剪切破坏。
物理模型与关键公式
莫尔应力圆的核心是平面应力状态的坐标变换。给定一点在x-y坐标系下的应力分量(σx, σy, τxy),可以求出任意倾斜截面上的正应力σ和剪应力τ。所有这些可能的(σ, τ)点构成了一个圆——莫尔圆。
$$ \sigma_{1,2}= \frac{\sigma_x+\sigma_y}{2}\pm \sqrt{\left(\frac{\sigma_x-\sigma_y}{2}\right)^2 + \tau_{xy}^2}$$
其中,σ₁和σ₂是两个主应力(最大和最小正应力)。公式中的根号项就是莫尔圆的半径R。圆心位于横坐标(正应力轴)上的 $(\frac{\sigma_x+\sigma_y}{2}, 0)$ 处。
最大剪应力τmax在数值上等于莫尔圆的半径R,它作用在与主平面成45度角的截面上。
$$ \tau_{max}= R = \sqrt{\left(\frac{\sigma_x-\sigma_y}{2}\right)^2 + \tau_{xy}^2}$$
主应力的方向角θp(从x轴逆时针转到σ₁方向)由下式确定:$ \tan 2\theta_p = \frac{2\tau_{xy}}{\sigma_x - \sigma_y} $。这个角度在模拟器中也会实时计算并显示。
主应力与最大剪应力公式
给定平面应力状态($\sigma_x, \sigma_y, \tau_{xy}$)时,主应力(剪应力为零的面上的正应力)与最大剪应力为:
$\sigma_{1,2} = \dfrac{\sigma_x+\sigma_y}{2} \pm R, \qquad R=\sqrt{\left(\dfrac{\sigma_x-\sigma_y}{2}\right)^2 + \tau_{xy}^2}$
$\tau_{max} = R = \dfrac{\sigma_1-\sigma_2}{2}, \qquad \tan 2\theta_p = \dfrac{2\tau_{xy}}{\sigma_x-\sigma_y}$
$\sigma_1$ 为最大主应力,$\sigma_2$ 为最小主应力,$\theta_p$ 为主方向。最大剪应力面与主应力面成45°,其大小等于主应力差的一半(莫尔圆半径 $R$)。本模拟器可由输入应力即时计算 $\sigma_1, \sigma_2, \tau_{max}, \theta_p$。
应力变换与莫尔圆作图
单元体旋转角 $\theta$ 后,面上的应力由应力变换方程给出。
$\sigma_\theta = \dfrac{\sigma_x+\sigma_y}{2} + \dfrac{\sigma_x-\sigma_y}{2}\cos 2\theta + \tau_{xy}\sin 2\theta$
$\tau_\theta = -\dfrac{\sigma_x-\sigma_y}{2}\sin 2\theta + \tau_{xy}\cos 2\theta$
$(\sigma_\theta, \tau_\theta)$ 的轨迹即为 $\sigma$–$\tau$ 平面上的莫尔应力圆。作图步骤:
- 圆心:$C=\left(\dfrac{\sigma_x+\sigma_y}{2},\ 0\right)$ 在 $\sigma$ 轴上。
- 代表点:绘出 $X(\sigma_x,\ \tau_{xy})$ 与 $Y(\sigma_y,\ -\tau_{xy})$。
- 圆:以 $XY$ 为直径作圆(半径 $R$)。
- 主应力:圆与 $\sigma$ 轴交点为 $\sigma_1, \sigma_2$,圆顶为 $\tau_{max}$。
- 2θ 法则:物理面转 $\theta$,圆上点沿同向转 $2\theta$。
在屈服准则中的应用
由莫尔圆得到的主应力与最大剪应力可直接用于屈服/破坏准则。
| 准则 | 屈服条件(平面应力) | 适用 |
| 最大剪应力(Tresca) | $\tau_{max}=\dfrac{\sigma_1-\sigma_2}{2}\le \dfrac{\sigma_Y}{2}$ | 延性金属,偏保守 |
| 畸变能(von Mises) | $\sqrt{\sigma_1^2-\sigma_1\sigma_2+\sigma_2^2}\le \sigma_Y$ | 延性金属,与试验吻合好 |
| 最大主应力(Rankine) | $\sigma_1 \le \sigma_Y$ | 脆性材料(铸铁、混凝土) |
延性材料因剪切屈服,宜用 Tresca 或 von Mises;脆性材料因最大拉应力破坏,宜用最大主应力准则。莫尔圆可推广到三维($\sigma_1,\sigma_2,\sigma_3$),最大剪应力为最大圆半径 $(\sigma_1-\sigma_3)/2$。
现实世界中的应用
土木与结构工程:在分析房屋的梁、柱或桥梁的关键连接部位(如焊接点、螺栓连接处)时,工程师使用莫尔圆来确定该点最危险的应力(主应力)和最容易发生剪切滑移的方向(最大剪应力面),从而进行安全校核。
机械设计与失效分析:设计齿轮的齿根、轴的键槽或压力容器的开孔附近时,这些地方应力状态复杂。通过莫尔圆分析,可以判断零件是更可能发生脆性断裂(由最大拉应力控制)还是塑性屈服(由最大剪应力控制),进而选择合适的材料和设计形状。
地质与岩土工程:分析地壳中某一点的应力状态,用于预测断层的滑动方向(与最大剪应力方向相关)或评估隧道、边坡的稳定性。莫尔圆是理解岩石剪切破坏准则(如库仑准则)的基础工具。
材料测试与科学研究:在实验室进行双轴拉伸或压缩试验时,莫尔圆用于将测得的载荷转换到材料内部晶粒的滑移面上,帮助研究者理解材料的微观变形机制和宏观力学性能之间的联系。
常见误解与注意事项
首先,需要明确“莫尔应力圆以二维(平面应力)为前提”。例如,它非常适合用于板或薄壁结构的分析,但要完整表达实体(三维)零件内部复杂的应力状态,则需要三个莫尔圆。可以认为本工具显示的是其中一个圆(由最大/最小主应力确定的圆)。
其次,剪切应力符号(方向)的设置错误非常常见。$\tau_{xy}$ 的正负不会改变圆的中心和半径,但圆会绘制在上方或下方,从而导致主应力作用面的角度(2θ)产生180°偏差。例如,可以比较 $\sigma_x=100$, $\sigma_y=0$, $\tau_{xy}=30$ 和 $\tau_{xy}=-30$ 的情况。主应力值相同,但圆会上下翻转。在实际应用中,如果坐标系定义未严格统一,可能会导致在错误的方向进行加固。
另外,“冯·米塞斯应力无法直接从莫尔圆读取”也是一个关键点。本工具会计算并显示两者,但冯·米塞斯应力是基于剪切应变能的等效应力,是包含三维影响的概念。在二维平面应力状态下,可通过 $\sigma_{vM} = \sqrt{\sigma_1^2 - \sigma_1\sigma_2 + \sigma_2^2}$ 由主应力计算得出。莫尔圆是求解主应力的手段,无法将冯·米塞斯应力本身绘制在圆上。
具体计算示例
以钢结构焊接节点为例:输入σx=120 MPa(水平拉力)、σy=30 MPa(竖向压力)、τxy=45 MPa(剪切力)。莫尔圆的圆心为[(120+30)/2, 0]=(75, 0),半径R=√[(120-30)²/4+45²]=√[2025+2025]≈63.6 MPa。主应力计算:σ₁=75+63.6=138.6 MPa,σ₂=75-63.6=11.4 MPa。最大剪切应力τmax=(138.6-11.4)/2=63.6 MPa,主应力方向角θp=0.5×arctan[2×45/(120-30)]=0.5×arctan(1)≈22.5°。
计算示例:焊接接头的主应力评价
根据焊接接头的实测应力(σ_x = 120 MPa、σ_y = 40 MPa、τ_xy = 60 MPa)求主应力:
- 主应力:σ_1,2 = (σ_x+σ_y)/2 ± √[((σ_x−σ_y)/2)² + τ_xy²]
- σ_1 = 80 + √(1600+3600) = 80 + 72.1 ≈ 152 MPa
- σ_2 = 80 − 72.1 ≈ 7.9 MPa
- 最大剪切应力:τ_max = (σ_1−σ_2)/2 ≈ 72 MPa
von Mises 应力 = √(σ_1²−σ_1σ_2+σ_2²) ≈ 148 MPa 确认安全系数。