在固体力学中，**应力圆（莫尔圆）**是分析二维应力状态、求取主应力、最大剪应力以及应力方向的强有力工具。通过几何作图，可直观地揭示应力分量之间的关系，帮助工程师快速判断材料的受力情况。本文将系统阐述莫尔圆的绘制步骤、数学原理以及实际解读方法，并给出典型数值示例。
平面应力状态：在 $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. 小结
莫尔圆通过简单的几何构造，将二维应力变换的代数关系转化为直观的图形表达。掌握圆心、半径的求法以及主应力、最大剪应力的读取技巧，可在材料强度分析、结构设计和失效预测中发挥重要作用。通过本文的步骤演示与数值示例，读者应能够独立完成莫尔圆的绘制、解读，并将其应用于实际工程问题的快速评估。