应力圆（莫尔圆）的绘制与解读
在固体力学中，**应力圆（莫尔圆）**是分析二维应力状态、求取主应力、最大剪应力以及应力方向的强有力工具。通过几何作图，可直观地揭示应力分量之间的关系，帮助工程师快速判断材料的受力情况。本文将系统阐述莫尔圆的绘制步骤、数学原理以及实际解读方法，并给出典型数值示例。
平面应力状态：在 $x$–$y$ 平面内，任一点的应力张量可表示为
[
\boldsymbol{\sigma}=
\begin{bmatrix}
\sigma_x & \tau_{xy}\
\tau_{xy} & \sigma_y
\end{bmatrix}
]其中 $\sigma_x,\sigma_y$ 为法向应力，$\tau_{xy}$ 为剪切应力（满足 $\tau_{xy}=\tau_{yx}$）。
应力变换公式：当坐标系逆时针旋转角度 $\theta$ 时，新的应力分量为
[
\begin{aligned}
\sigma_{x'} &= \frac{\sigma_x+\sigma_y}{2}
+\frac{\sigma_x-\sigma_y}{2}\cos2\theta
+\tau_{xy}\sin2\theta,\[4pt]
\sigma_{y'} &= \frac{\sigma_x+\sigma_y}{2}
-\frac{\sigma_x-\sigma_y}{2}\cos2\theta
-\tau_{xy}\sin2\theta,\[4pt]
\tau_{x'y'} &= -\frac{\sigma_x-\sigma_y}{2}\sin2\theta
+\tau_{xy}\cos2\theta .
\end{aligned}
]莫尔圆的几何意义：在 $(\sigma,\tau)$ 坐标系中，以 $\displaystyle \sigma_{\text{avg}}=\frac{\sigma_x+\sigma_y}{2}$ 为圆心，半径
[
R=\sqrt{\left(\frac{\sigma_x-\sigma_y}{2}\right)^2+\tau_{xy}^2}
]绘制的圆即为莫尔圆。圆上任一点对应于某一坐标系旋转角度 $\theta$ 的应力状态。
2. 绘制莫尔圆的步骤
确定已知应力分量
- 给定 $\sigma_x,\sigma_y,\tau_{xy}$（注意符号约定，拉伸为正，剪切正方向通常取逆时针为正）。
计算圆心与半径
- 圆心：$\displaystyle C=\sigma_{\text{avg}}=\frac{\sigma_x+\sigma_y}{2}$。
- 半径：$\displaystyle R=\sqrt{\left(\frac{\sigma_x-\sigma_y}{2}\right)^2+\tau_{xy}^2}$。
在 $\sigma$–$\tau$ 坐标系中作图
- 横坐标为法向应力 $\sigma$，纵坐标为剪切应力 $\tau$。
- 标出点 $A(\sigma_x,\tau_{xy})$ 与点 $B(\sigma_y,-\tau_{xy})$，两点必在圆上且互为直径端点。
绘制圆
- 以 $C$ 为圆心，以 $R$ 为半径画完整圆。
标注主应力与最大剪应力
- 主应力：圆上与 $\tau$ 轴交点，即 $\sigma_{1,2}=C\pm R$（对应 $\tau=0$）。
- 最大剪应力：圆上与 $\sigma$ 轴交点，即 $\tau_{\max}=R$（对应 $\sigma=C$）。
求取对应的方向角
- 由几何关系 $\tan 2\theta = \dfrac{2\tau_{xy}}{\sigma_x-\sigma_y}$，可得到使 $\tau_{x'y'}$ 为零（即主应力方向）的坐标系旋转角 $\theta_p$。
- 最大剪应力方向相差 $45^\circ$，即 $\theta_{\tau\max}= \theta_p + 45^\circ$。
3. 示例计算
已知某平面应力状态为
[
\sigma_x = 80\ \text{MPa},\qquad
\sigma_y = 20\ \text{MPa},\qquad
\tau_{xy}=30\ \text{MPa}.
]
3.1 圆心与半径
[
\begin{aligned}
C &= \frac{80+20}{2}=50\ \text{MPa},\
R &= \sqrt{\left(\frac{80-20}{2}\right)^2+30^2}
=\sqrt{30^2+30^2}=42.43\ \text{MPa}.
\end{aligned}
]
3.2 主应力
[
\sigma_1 = C+R = 92.43\ \text{MPa},\qquad
\sigma_2 = C-R = 7.57\ \text{MPa}.
]
3.3 最大剪应力
[
\tau_{\max}=R=42.43\ \text{MPa}.
]
3.4 主应力方向
[
\tan 2\theta_p = \frac{2\tau_{xy}}{\sigma_x-\sigma_y}
= \frac{60}{60}=1
;\Rightarrow;
2\theta_p = 45^\circ ;\Rightarrow; \theta_p = 22.5^\circ .
]
因此，主应力方向相对于原 $x$ 轴逆时针为 $22.5^\circ$，最大剪应力方向为 $22.5^\circ+45^\circ=67.5^\circ$。
3.5 作图要点
- 在 $\sigma$–$\tau$ 平面上标出点 $A(80,30)$ 与点 $B(20,-30)$。
- 以 $C=50$ 为横坐标中心，半径 $42.43$ 画圆。
- 圆与 $\tau=0$ 交点即为 $\sigma_1,\sigma_2$，与 $\sigma=50$ 交点即为 $\pm\tau_{\max}$。
4. 莫尔圆的工程意义
- 失效判定：在屈服准则（如 von Mises、Tresca）中，只需读取圆上相应的应力值即可快速判断是否超过材料屈服强度。
- 应力分布分析：通过比较不同截面或不同加载工况的莫尔圆，可直观评估应力集中程度。
- 设计优化：若主应力差距过大，可通过改变结构形状或加载方式，使莫尔圆半径减小，从而降低最大剪应力，提高结构安全性。
5. 常见误区与注意事项
- 符号统一：绘图时必须保持拉伸为正、剪切正方向（逆时针）一致，否则圆的相对位置会出现镜像错误。
- 平面应力假设：莫尔圆仅适用于平面应力或平面应变问题，三维应力需分别绘制三个互相垂直的莫尔圆。
- 数值精度：在手工作图时，半径的计算误差会直接影响主应力和最大剪应力的读数，建议使用计算器或软件辅助。
6. 小结
莫尔圆通过简单的几何构造，将二维应力变换的代数关系转化为直观的图形表达。掌握圆心、半径的求法以及主应力、最大剪应力的读取技巧，可在材料强度分析、结构设计和失效预测中发挥重要作用。通过本文的步骤演示与数值示例，读者应能够独立完成莫尔圆的绘制、解读，并将其应用于实际工程问题的快速评估。