数学总导航 · 本课程目录 · ← 上一章 · 下一章 →
两端固定的晾衣绳,挂上重量后弯成什么形状?#
将绳子拉紧,两端固定,均匀挂上轻小负载。用 u(t) 表示向下位移,一个归一化的小挠度模型是 −u′′=1、u(0)=u(1)=0。
解是 u(t)=t(1−t)/2,中点下垂最多,为 1/8。这个数依赖我们选择的单位与归一化张力;它只是用来展示计算机制,不是任何实绳的直接工程预测。
如果重量分布粗糙甚至分段变化,直接要求每一点都存在二阶导数会很别扭。弱形式把要求转成整体平衡:让每一种小的试探位移都满足相应的积分关系。接下来会逐步说明这个转写过程,以及为什么它保证唯一解。
从绳子的方程走到积分平衡,中间不能省#
先假设 u,v 足够光滑,且试探位移 v 在两端为零。把 −u′′=f 乘 v 再积分:−∫u′′v=∫fv。分部积分后,左边为 −[u′v]01+∫u′v′。
因为 v 端点为零,边界项消失,留下 ∫u′v′=∫fv。现在两边只要求一阶弱导数与函数平方可积,就可以把这条关系延伸为弱解定义。
这不是说所有经典方程都可无条件倒推回来;从弱解得到更强的经典正则性还需额外条件。这里先保证弱意义下问题完整且稳定。
把能量增量完整展开#
能量 E(v)=21∫∣v′∣2−∫fv。在候选 u 周围改变 εh,有
E(u+εh)−E(u)=ε(∫u′h′−∫fh)+2ε2∫∣h′∣2.若 u 为最小点,一阶系数对所有 h 必须为零,便得到弱形式。反过来,一旦一阶系数全零,余下二次项非负,说明 u 真是全局最小点。固定边界又排除非零常数方向,使最小点唯一。
用一个形状模板就能算出答案#
先只允许 uN(t)=ct(1−t)。代入能量,∫01(1−2t)2dt=1/3,∫01t(1−t)dt=1/6;当 f=1,
E(c)=6c2−6c=61(c−1/2)2−241.因此 c=1/2 最优。这个形状模板恰好包含精确解,所以一次就找到答案。一般载荷下,一个模板不够,就增加更多基函数,得到 Galerkin 线性系统。
强制性到底防住了什么#
强制性要求能量二次部分至少按 α∥u∥2 增长,防止解逃到很远但代价不增加。连续性则控制相互作用不会突然放大。
Lax–Milgram 把这两种控制结合起来,给存在、唯一和稳定性。Céa 引理继续告诉你:数值解与真解的差,不会比所选有限维空间能做到的最佳近似差太多。理论并不替你选择网格,却能明确指出改善近似应往哪里努力。
带你做:载荷加倍,解为什么也加倍#
对弱问题 a(u,v)=F(v),如果载荷变为 2F,线性性使 2u 满足新方程;唯一性保证它就是新解。因此在线性小挠度模型中,重量加倍,下垂也加倍。
若载荷发生小扰动 ΔF,两解之差满足 a(Δu,v)=ΔF(v)。取 v=Δu,由强制性得
α∥Δu∥2≤∥ΔF∥∥Δu∥,∥Δu∥≤α∥ΔF∥.这份估计直接回答“称重误差会带来多少形状误差”,而不是只保证一个抽象解存在。
真实绳子在大变形时可能不再线性,材料也可能改变张力。模型中的加倍规律依赖线性前提;数学推导清楚,才能知道实际应用时要核验哪一种近似是否仍合理。
回头整理:定义、工具与练习#
不直接追问曲线的二阶导数#
考虑一根固定两端的弦:
−u′′=f在 (0,1),u(0)=u(1)=0.如果 f 不连续,要求处处存在经典二阶导数可能太严格。把方程乘测试函数后积分,转移一次导数,得到
∫01u′v′dt=∫01fvdt对所有 v∈H01(0,1).现在只需一阶弱导数平方可积。这叫弱形式,边界条件通过函数空间编码。它不是随意放宽,而是保留积分意义下的全部平衡关系。
Lax–Milgram 定理#
以下使用实 Hilbert 空间版本。设 a:H×H→R 为双线性形式,满足
∣a(u,v)∣≤M∥u∥∥v∥,a(v,v)≥α∥v∥2,α>0.后者称强制性。对每个有界线性泛函 F,存在唯一 u∈H 使 a(u,v)=F(v) 对全部 v 成立,并有 ∥u∥≤∥F∥/α。
证明路线:Riesz 表示把 a 写为 ⟨Au,v⟩。强制性给 ∥Au∥≥α∥u∥,故 A 单射且像闭。若 y 垂直于像,取 v=y 得 a(y,y)=0,强制性迫使 y=0,故像还稠密,从而满射。范数界由把 v=u 代回得到。
复数形式应改用半双线性并要求 Rea(v,v)≥α∥v∥2,本章不混用两种约定。
核验 Poisson 问题的每个条件#
取 H=H01(0,1),用 ∥u∥H=∥u′∥2,由 Poincaré 不等式知它与原 H1 范数等价,故完备。
令 a(u,v)=∫u′v′,则有界常数和强制常数都可取 1。若 f∈L2,F(v)=∫fv 满足
∣F(v)∣≤∥f∥2∥v∥2≤∥f∥2∥v′∥2.因此弱解存在唯一,而且 ∥u′∥2≤∥f∥2。这同时给了输入扰动的稳定性估计。
当 f=1 时,显式解为 u(t)=t(1−t)/2。把它代入弱形式可验证;对一般 L2 输入,即使没有漂亮公式,定理仍给结论。
能量最小化#
当 a 对称时,同一个 u 最小化
E(v)=21a(v,v)−F(v).展开可得 E(u+h)−E(u)=a(h,h)/2≥α∥h∥2/2。所以驻点不仅局部最优,而且唯一全局最优。若 a 不对称,Lax–Milgram 仍可成立,但不能用这个能量泛函直接产生原方程。
更一般的直接法:证明能量有下界并强制;取极小化序列;用自反性提取弱收敛子列;保证可行集弱闭、能量弱下半连续;最后用严格凸性得唯一。每一步有独立条件。
Galerkin:回到有限个未知数#
选有限维子空间 HN=span{ϕ1,…,ϕN},设 uN=∑cjϕj,要求 a(uN,ϕi)=F(ϕi)。得到刚度矩阵 Aij=a(ϕj,ϕi)。
Galerkin 正交性为 a(u−uN,vN)=0。Céa 引理给出
∥u−uN∥H≤αMvN∈HNinf∥u−vN∥H.证明把误差的强制性界与正交性结合,再用双线性形式的有界性即可。它告诉我们:离散解的误差与所选空间的最佳近似误差同阶,泛函分析真正连接了存在性与数值算法。
练习与解答#
- 对 f=1、H1=span{t(1−t)},求 Galerkin 系数。
- 若删去零边界条件,只用导数范数会怎样?
解答:1. ∫(1−2t)2dt=1/3,∫t(1−t)dt=1/6,所以系数为 1/2,恰好恢复精确解。2. 非零常数的导数范数为零,不再是范数,强制性与唯一性都需要重新处理。
数学总导航 · 本课程目录 · ← 上一章 · 下一章 →