第80夜 | 压杆稳定7：一端固支一端铰支0.7背后的真相
地点：力学研讨室。人物：
小生：翻看着前几天的推导草稿，眼神中带着求知欲。 老师傅：在黑板上画下了一根一端固支、一端铰支的压杆。 
1. 难缠的铰支端：多出来的未知力
小生：师傅，前两夜咱们用微积分严丝合缝地证明了“一端固支一端自由（ ）”和“两端固支（ ）”。今天终于要面对最后一个硬骨头了——“一端固支一端铰支”压杆！那个工程中常用的 到底是怎么算出来的？我已经迫不及待想列平衡方程了。
老师傅：很好，保持这种探索欲。我们还是以底端固支点为原点建立向上的 轴。 这次顶端变成了铰支端，小生，你分析一下，顶端的约束情况和之前有什么不同？
小生：铰支端允许转动，所以没有顶端约束弯矩。但是，它不允许水平移动！
老师傅：没错。那么当压杆受压发生微弯时，为了阻止顶端发生水平位移，铰支座会产生什么样的反力？
小生：必然会产生一个水平方向的约束反力！我们假设这个水平 反力为 。
老师傅：好，现在你在任意截面 处切开，取上半部分研究，试着列出合力矩为 0 的平衡方程。规定全局的正方向为逆时针为正，顺时针为负。
小生：我来看看……截面上假设有正向内弯矩 （使压杆凸向右侧），在取上半部分分离体时，这个内弯矩作用在切口底部，方向是顺时针的，所以记为 。 上半部分的外力有两个：向下的压力 和向左的水平 反力 （假设它向左）。 力 的作用线偏离当前截面形心的距离是 （因为顶端挠度为0），它相对截面形心产生顺时针的力矩，记为 。 力 距离当前截面的竖直高度是 ，它向左作用，相对截面形心产生逆时针的力矩，记为 。 根据合力矩等于 0：
移项得到弯矩表达式：
2. 非齐次方程的三重变身
老师傅：完美。接下来，把它代入挠曲线近似微分方程 ，并化简成标准形式。
小生：代入进去：
把含有 的项移到左边，并令 ：
师傅，等号右边不再是单纯的常数了，而是一个关于 的一次函数。这是一个非齐次方程。齐次解我知道，依然是 。但是特解怎么找？
老师傅：既然等号右边是一个关于 的一次函数，你可以大胆猜测，特解应该也是一个什么样的函数？
小生：应该也是个一次函数！因为一次函数的二阶导数也是0呀！我观察一下，如果令特解刚好等于右边那项除以 ，也就是 ，求它的二阶导数刚好是 0。
老师傅：你把它代入方程左边验证一下。
小生：代入左边…… ，完美相等！这就是我们要找的特解！ 所以，非齐次方程的通解就是齐次解加上特解：
3. 边界条件的代入与线性方程组的推导
老师傅：通解写得很利索。现在，面对通解里的三个未知数 、 以及未知力比例 ，我们需要代入边界条件，去确定它们的值。你先把导数 准备好。
小生：求导得到： 。 这根杆子有三个边界条件： 底端固支： 时， ；且 。 顶端铰支： 时， 。
我把它们一个个代进去： 第一个： 第二个： 第三个：
老师傅：好，现在你拥有了一个包含三个未知量（ ）的三元一次齐次线性方程组。我们把它写成矩阵形式。
小生：提取系数，矩阵就是：
师傅，这在线性代数里是一个标准的齐次线性方程组。对于这种方程组，通常有两种情况。 第一种情况，是 , , 。
老师傅：这种情况在数学上叫“平凡解”。如果套用在我们这根压杆上，就意味着压杆的挠度 处处为 0，也就是压杆根本没有发生变形，还是笔直的。这显然不是我们想要寻找的失稳状态。我们需要压杆真的发生弯曲失稳，也就是变形不全为 0。这就要求这个方程组必须存在“非零解”（或者叫非平凡解）。而在代数中，一个齐次方程组要有非零解的充要条件，需要满足什么条件？
小生：我想想，线性代数里面讲过，这种情况需要系数矩阵的行列式必须等于 0！
老师傅：解释得非常透彻。平凡解代表着安全的直线平衡状态，而非平凡解的出现，就是压杆失稳的数学宣告。这就是决定压杆能否失稳的关键条件。展开这个行列式，看看你能得到什么？
小生：我按第一行展开： 化简后就是： 两边同时除以 ，得到一个极致简洁的超越方程：
4. 超越方程的解与 0.7 的真相
小生：师傅，我卡住了。 这种既有三角函数又有一次函数的超越方程，代数上根本求不出精确解啊！
老师傅：代数求不出，但可以通过数值分析或图解方法。为了看得更清楚，我们可以先把它抽象成一个标准的数学方程：求解 。
小生：我明白了，把方程两边看作两个独立的函数： 和 。在坐标系里画出正切曲线和一条过原点、倾角为 45 度的直线。它们的交点，就是方程的根！
老师傅：没错。排除掉不代表失稳的原点 ，通过图形和数值计算，它们的第一象限最小非零交点发生在哪里？
小生：我查一下……最小非零解大约是 （弧度）。
老师傅：很好。现在，我们把 代回来，也就是：
小生：得到了这个数值，剩下的就水到渠成了！ 把 代进去：
两边平方，把 移过去：
为了看清楚 ，我们要把它凑成 的形式！ 既然 ，那 大约等于 。而 刚好是 ！ 所以：
老师傅：大功告成！你看到了吗，那个在工程手册里被查阅了无数遍的近似值 ，它的根源其实是数学方程 的一个无理数解 。 从合力矩的物理洞察，到非齐次微分方程的建立，再到矩阵行列式的化简，最后通过数值求解揭开真相。小生，你用微积分的铁锤，把压杆稳定性的所有基石都亲手夯实了。
