Stress Analysis of a Square Elastic Body Under Biaxial Loading Using Airy Stress Functions

本研究通过推导满足双调和方程及边界条件的艾里应力函数闭式解，对单轴及双轴压缩载荷下正方形弹性体的应力分布进行了分析性研究，其结果与实验光弹性数据高度吻合。

要点

想象一下，你有一块完美的正方形橡胶。现在，想象有人用力按压它的顶端和底端，或者同时向四个边施压。那块橡胶内部会发生什么？压力是均匀地扩散开来，还是会在某些奇怪的地方被挤压？

这篇论文就像是一个数学水晶球，让我们无需切开橡胶块，也无需使用昂贵的计算机模拟，就能精确观察到其内部正在发生的一切。作者是来自东京农业大学的研究人员，他们利用了一个经典的数学工具——**艾里应力函数（Airy Stress Function）**来解开这个谜题。

以下是他们研究工作的通俗易懂版解析：

问题所在：正方形很棘手

长期以来，科学家们一直知道如何计算圆形物体（比如被挤压的硬币）内部的应力。这就像是在用一个圆形的框架解谜，数学逻辑流动得非常顺畅。但当形状变成正方形时，数学变得非常复杂。由于角点和直边的存在，很难找到一个完美的、精确的公式。通常情况下，工程师必须依赖计算机程序（如有限元分析）来提供近似的答案。

这篇论文说：“让我们为正方形寻找一个精确的答案。”

方法：一种“应力配方”

为了解决这个问题，作者使用了一种特殊的数学配方（艾里应力函数）。你可以把这个配方想象成一把万能钥匙，它能自动平衡材料内部所有的作用力，使它们不会四散崩裂。

1. 分解过程： 他们将边缘承受的复杂压力分解成一系列简单的波（就像池塘里的涟漪）。
2. 无穷级数： 他们写出了一个公式，通过累加成千上万个这样微小的波来构建出总体的应力图景。
3. 调节旋钮： 他们必须不断调整每个波的“音量”（数学系数），直到边缘处的压力与预期的目标（无论是强力挤压还是平滑挤压）完全吻合。

研究结果：他们的发现

1. “简单模式”测试：
首先，他们用一个简单案例测试了他们的数学模型：向四周均匀施压。正如预期的那样，内部的应力是完全均匀的。这证明了他们的“配方”运行正确。

2. “挤压”测试（单轴加载）：
接下来，他们模拟了仅对顶端和底端施压的情况（类似于巴西坚果实验）。

• 令人惊讶的发现： 在圆盘中，中间部分的张力（拉伸力）是完全笔直且均匀的。但在正方形中，作者发现顶端和底端附近的应力并不是平坦的。因为正方形拥有角点和直边，材料抵抗挤压的方式各不相同，从而在受力位置产生了一个“凹陷”或局部的应力变化。
• 证据： 他们将自己的数学推导与现实世界中受压塑料的照片（称为光弹性实验）以及计算机模拟进行了对比。他们的数学“水晶球”与现实照片几乎完美契合。

3. “双重挤压”（双轴加载）：
最后，他们观察了同时对顶/底边和左/右边施压时会发生什么。

• 他们发现，内部的应力变成了两种挤压力的复杂混合体。根据你在内部观察的位置不同，最强应力与最弱应力之间的“差异”也会随之改变。这就像混合两种不同的颜料；结果取决于你在混合物中的具体位置。

为什么这很重要（根据论文所述）

作者并不是在声称这项研究明天就能治愈疾病或建造新桥梁。相反，他们是在提供一个金标准参考。

• 基准： 正如需要一把尺子来检查卷尺是否准确一样，这种精确的数学解是用来检查计算机模拟是否工作正确的必要工具。
• 洞察力： 它揭示了正方形材料行为中隐藏的细节，而这些细节是圆形物体数学模型所无法捕捉到的。它表明，物体的形状（正方形对比圆形）实际上会改变在你指尖下方的应力流动方式。

简而言之，这篇论文为我们提供了一张精确的、关于正方形材料内部不可见力量的地图，证明了即使是在如此简单的形状中，物理学也可以表现得如此复杂且独特。

技术摘要

技术摘要：基于 Airy 应力函数研究双轴载荷下正方形弹性体的应力分析

问题陈述
本研究旨在通过解析方法，确定在单轴和双轴条件下，承受集中载荷和分布压缩载荷的正方形、均匀、各向同性且线弹性体的应力分布。虽然经典弹性理论为圆形和椭圆形区域（例如直径加载的圆盘）提供了闭式解，但对于多边形几何形状（尤其是正方形）的解析处理仍然较少，这是由于缺乏径向对称性以及应用实际边界条件的复杂性所致。尽管有限元法（FEM）常用于此类问题，但其本质上是近似解。本研究旨在填补这一空白，为正方形域提供精确或半解析解，这类几何形状在实际工程场景中经常遇到，例如间接抗拉强度评估（如巴西试验）和多轴材料测试。

方法论
本分析基于二维线性弹性理论，利用 Airy 应力函数 ($\phi$) 方法。该方法将控制平衡方程简化为单一的双调和方程 ($\nabla^4\phi = 0$)，在不存在体力的情况下，该方程本质上满足平衡条件。

1. 公式化： 问题定义在一个边长为 $2L$、中心位于原点的正方形区域内。外部法向应力 $\sigma_1(y)$ 和 $\sigma_2(x)$ 分别施加于侧向和垂直边缘，所有边界上的剪应力均为零（无摩擦或无牵引条件）。
2. 级数展开： 为了处理任意应力分布，将边界条件展开为余弦傅里叶级数。Airy 应力函数被构建为一个包含多项式项以及双曲与三角函数（$\cosh$, $\sinh$, $\cos$）无穷级数的可分离解。
3. 系数确定： 通过强制执行边界条件来确定级数展开中的未知系数。这一过程产生了一个耦合的线性方程组。
4. 数值稳定性： 为了解决高阶项（$n, m \gg 1$）系数发散的问题，作者引入了重标度变量和修正系数。这种变换确保了数值稳定性，使得无穷系统可以通过截断至有限阶数 ($N_{max}$) 并通过矩阵求逆进行求解。

主要贡献

• 广义解析框架： 本文推导出了在任意对称双轴载荷下正方形弹性体应力场的闭式半解析解，这类几何形状在经典的坐标分离系统中难以处理。
• 收敛性分析： 研究展示了级数解的收敛行为。研究表明，截断阶数 $N_{max} = 128$ 足以提供足够的数值精度，且 $N_{max}=128$ 与 $256$的结果重合，相对误差为$10^{-4}$。
• 验证： 该方法通过已知解（如均匀双轴拉伸，恢复空间均匀应力场）进行了验证，并针对集中载荷情况，与已知的光弹性条纹图样及有限元分析（FEM）进行了对比。

结果

• 单轴压缩： 在集中单轴压缩下，所得应力场与先前观察到的光弹性条纹图样高度一致。主应力差 ($\Delta\sigma$) 等值线与实验干涉条纹对齐。
• 几何效应： 一个显著的发现是沿加载线 ($x=0$) 的拉应力分量 ($\sigma_{xx}$) 的行为。与圆形域中拉应力保持恒定的情况不同，正方形几何结构由于周围材料提供的足够阻力，在载荷施加点附近诱发了局部压缩。这导致靠近边缘处出现非均匀的 $\sigma_{xx}$ 分布，这是圆形域解中未捕捉到的特征。
• 双轴压缩： 对于双轴载荷，该解代表了两个正交压缩场景的叠加。分析揭示了主应力差的空间变化，这种变化取决于域内的特定位置以及施加载荷的比率。

意义
作者将这项工作定位为理论验证、基准测试和教学用途的有力工具。通过提供精确或半解析表达式，本研究为验证计算模型（如 FEM）以及解释涉及正方形或矩形试样的机械测试装置中的实验数据提供了参考标准。该框架被视为未来扩展到矩形、各向异性或分层域，以及动态和三维弹性问题的垫脚石，同时并未声称在分析的载荷配置范围之外具有直接的工业应用价值。