第8章 特征值问题的变分原理8.1 Sturm-Liouville微分方程与特征值在求解微分方程、结构的稳定性或者求结构的固有频率时，我们经常会遇到下面的微分算子 (8.1.1)其中都是已知的函数，，那么方程 (8.1.2)称为Sturm-Liouville方程，其中权函数，当且仅当在的一个零测度集上等号成立当给定了齐次边界条件后，只有某一些特定的才能使得该方程有非零解, 使得该方程有非零解的称为特征值，相应的解称为特征函数常见的边界条件为(1) 两端固定∶2) 两端自由∶3) 一端固定、另一端自由∶或， 或我们可以在复函数空间中定义一个内积运算为 容易证明，，即是对称(自伴)算子如果记 ()为Sturm-Liouville方程的特征值是相应的特征函数也就是说 (8.1.3)那么，对于特征值和特征函数，我们可以得到以下一些性质：1. 所有特征值是实的 若是一组特征值和特征函数，即 (8.1.4)则也是一组特征值和特征函数，即 (8.1.5)将(8.1.4) 乘、(8.1.5) 乘相减并积分可得 2. 特征函数正交性 由算子的对称性 另一方面，由于是算子的特征向量，所以有 因此` 当时，要求上式成立，只有 当时，若是的两个线性无关的特征向量，选择 代替，满足正交性要求。
对于有多个线性无关特征向量的重特征值问题，也可类似处理(Schmit正交化)这样，我们总是可以选择合适的特征函数，使得 (8.1.6)也就是说可以把特征函数单位正交化否则，我们只要把得到的特征函数作下面的变换就可以3. 特征函数的富里叶展开对于任意一个连续函数，均可以用Sturm-Liouville算子的特征函数进行富里叶展开 (8.1.7)其中 这里，严格的证明我们不去讨论8.2 Sturm-Liouville特征值问题的Rayleigh原理根据上面Sturm-Liouville方程的算子及内积定义，对于任意的函数,定义 (8.2.1)我们称该泛函为算子的Rayleigh商定理8.1 上述定义的Rayleigh商与算子的特征值有 (8.2.2)这里表示泛函取驻值，此时得到的为对应特征值的特征函数。
当时，所有特征值，此时取最小特征值,则(8.2.2) 变成 (8.2.3)证明: 对任何一个函数，按Sturm-Liouville方程的特征函数进行富里叶展开其中那么由于 因此有 从而即 当时，式(8.2.3) 可由性质(2)得到 定理8.1称为Sturm-Liouville特征值问题的Rayleigh原理8.3 特征值问题的Rayleigh-Ritz法 根据Rayleigh原理，Ritz提出了求解Sturm-Liouville微分方程特征值的近似计算方法:首先把特征值问题转化为变分问题(8.2.2)，然后再用数值方法来求解该变分问题令 其中是待定的常数，是选定的一系列基函数，它们满足指定的边界条件在实际应用中最好从一组完备的函数系中来选取基试算函数，如幂函数，三角函数等将的表达式代入的定义中，可以得到 其中 而且矩阵和是对称的要使得上面的取到最小值，那么必定要求满足 从而得到 也就是 这是一个(广义)代数特征值问题，可以通过迭代方法，SVD方法或者其他数值方法来求解。
例8.1 解：其特征值为，特征函数为Rayleigh商为如果取近似函数为，那么，代入的表达式中得到 它比真实的稍大如果取近似函数为，代入的表达式中得到，它和真实的几乎相等8.4 Sturm-Liouville四阶微分方程的特征值问题 Sturm-Liouville四阶微分算子为 (8.4.1)特征方程为 (8.4.2)这里边界条件为每端各取下列两个边条件(1) 或者；(2) 或者与该方程特征值问题等价的变分问题为 8.5 结构的稳定性结构的平衡状态可以分为三类: 稳定平衡、不稳定平衡、随遇平衡在工程中经常会遇到结构失稳问题：(1) 细长杆受压，当压力从零开始增加时，杆件保持为直线，当压力到达一定值的时候，杆件被压弯，产生较大的变形2) 板条或者工字梁在最大抗弯刚度平面内弯曲，当载荷到达一定值时，会发生侧向弯曲与扭转3) 圆柱壳的失稳定义8.1（稳定性）设结构处于某一平衡状态，受到任一微小扰动后而稍微离开原平衡位置当扰动消失后，如果结构能回到原来位置，则称此平衡状态为稳定平衡状态；如果结构可能继续偏离，不能回到原来位置，则称此平衡状态为不稳定平衡状态。
介于稳定平衡和不稳定平衡之间的过渡平衡，称为临界平衡状态，简称临界状态注：这里我们用临界平衡替代前述的随遇平衡这是两个不同的概念，因为临界平衡可能是随遇平衡，也可能是稳定平衡或不稳定平衡，它在工程上更有用定理8.2 设为系统的总势能，是结构的位移函数，则其平衡点必定满足，并且（1）对于任意，，则必定是稳定的平衡点；（2）至少存在一个非零的，使得，则必定是不稳定平衡点；（3）对于任意，，并且至少有一个非零的使得不等式中的等号成立，则必定是系统的临界平衡点定理8.3 平衡点稳定的充要条件是使得总势能取严格极小值例8.2 图8.1（a）为一刚性压杆（不变形），承受中心压力为，底端为铰支座，顶端有弹簧系数为的水平弹簧支承 （a） （b）图8.1例8.2图解：当为竖直时，系统能平衡，这是原始的平衡形式现在考虑倾斜位置是否还存在新的平衡状态为此，写出平衡条件（水平方向） 即 这个方程有两个解 这正好是原始平衡状态（Ⅰ），另一解为 这是新的平衡路径（Ⅱ）将这些解画在图（b）上，显然分支点为 分支点将原始平衡路径Ⅰ分成两段：前段上的点属于稳定平衡，而后段属于不稳定平衡。
而在新的平衡路径上，当载荷减少时倾角反而增大，所以也是属于不稳定平衡对于这类具有不稳定分支的完善体系，进行稳定性验证时要特别小心，一般应考虑初始缺陷（初曲率、偏心）的影响现在按小挠度理论分析所谓小挠度理论，是将平衡方程中位移分量按小量线性化处理由于 代入方程（a）得 故两个平衡态为 ，前者是原始平衡状态，后者是新的平衡状态图8.2 小挠度理论结果显然，按小挠度理论计算出的分支点与大挠度完全一样；但对分支点以后的情形，小挠度给出的随遇平衡是一种假象例8.3 考虑图8.3（a）单自由度非完善体系，刚性杆有初倾角，其余同上例 (a) (b) (c)图8.3例8.3图解：在图（b）中，平衡条件（水平方向）为 （a）所以 （b）对于不同的，曲线画在图（c）上曲线上有极值点，为此令，解得 代入式（b）得到相应的极值载荷为 （c）显然，初倾角越大，临界载荷就越小现在用小挠度理论设，从而 代入式（b）得 （d）对于所有的，上述的曲线均是以作为渐近线的，即很显然，式（d）中的曲线无极值。
稳定性问题的出现是由于整个结构总势能的正定性破坏能造成的，所以只要讨论正定性破坏的条件，就可以得到所需的临界载荷1）是平衡点要满足的条件，表明的正定性被破坏所以，或有非零解存在可以得到临界载荷的值2）一般我们考虑的是小应变、大变形问题换言之，由于应变量是小量，所以应变能可以按原来的公式计算；但外力势能，则要按变形后的位置计算例8.4 用能量法计算前面例子中的临界载荷解：点弹簧的弹性势能为 当转过角度时（图8.2（b）），外力沿方向移动了： 所以外力的势能为 总势能为 由得 从其非零解条件可得 例8.5图8.4 两端简支的压杆首先来考虑两端是简支的压杆选取坐标系如图所示，则其弯矩为 用表示的平衡方程（适用于小挠度） 引入记号 则方程的通解为（当时） 由的条件可得，再将的条件代入上式可得 为使方程有非零解（第二平衡路径），必须有 因此 取最小的作为临界载荷选，即 这就是我们要求的两端简支梁的临界载荷8.6 求压杆临界载荷的变分方法考虑一根直梁，在轴向力作用下的失稳问题由于没有横向载荷，挠度满足方程 （8.6.1）梁的边界条件为每端各取下列两个条件∶(1) 或者 ；(2) 或者 .显然是该方程的解。
对于一些特殊的，该方程有非零解，我们称的这些值为压杆的临界载荷，相应的为特征函数对应的变分形式为 （8.6.2）8.7 临界载栽荷的Rayleigh-Ritz法以两端简支的压杆稳定性为例由于作用，其在微元上所作用的弯矩为 当挠度由变到时，这一弯矩所做的功为 从而对整根梁来说，力所做的功为 因为, 式中处的项自动消失这样外力势能为 而总势能为 设挠度曲线可以近似写成 这里满足边界条件，为已知的插值函数系，为待定参数，代入总势能表达式得 式中 从可以得到一组要满足的代数方程 。