新闻资讯

暖气片钢板为什么只鼓不裂：格里菲斯公式√(2Eγ/πa)里的π来自半椭圆，钢的γ被塑性放大了上千倍

2026年11月9日0

一、打压时暖气片先鼓，玻璃却一碰就碎

做过打压试验的人都有一个印象：1.5 毫米厚的钢板水道被压到 0.8 兆帕，板面会明显地鼓出来一块，像被人用手指从里面顶了一下，但怎么压都不裂；而同样厚度的一块玻璃，磕一下就顺着裂纹走到底。两者的弹性模量只差三倍（200 对 70 吉帕），抗拉强度差得也不离谱，行为却是"鼓"和"碎"两个世界。

这个差别不在强度里，在一个几乎没人注意的希腊字母里：γ，表面能。1920 年格里菲斯正是靠它，把"材料为什么远达不到理论强度"这件事算清楚。今天把他那条 σ = √(2Eγ/πa) 一次拆完：分子里的 2 和分母里的 π 各是哪来的，以及为什么这个公式直接套到钢板上会算出荒唐的答案。

二、能量账本：一个涨得像面积，一个涨得像长度

取一块板，中间一条贯穿裂纹，总长 2a，板厚取单位 1。远处拉应力 σ。

第一笔账是释放的弹性能。裂纹张开等于把材料里原本被拉紧的那部分松开，可以证明这部分能量是 Ue = πσ²a²/E——注意它是 a 的平方，因为裂纹张开释放的能量正比于"受影响区域的面积"。

第二笔账是新增的表面能。裂纹每往前推一点，就多出两块新表面，每块面积正比于长度。整条裂纹上下两个面、两个裂尖，总共是 Us = 4aγ，它是 a 的一次方。

把两笔账相减，总能量 U = 4aγ − πσ²a²/E。一个线性增长、一个平方增长，平方项起先很小、后来居上，中间必有一个最高点。对 a 求导令其为零：dU/da = 4γ − 2πσ²a/E = 0，解出

σf = √(2Eγ / πa)

平方项和线性项相减取极值，这是所有"1/2 次方"公式的共同出身：能量按体积算、门槛按面积算，两者对 a 的幂次差一阶，极值处自然剩一个 a 的一次方在分母，开方之后就变成 √a。这个结构在暖气立管压屈那篇里也出现过，只不过那边是 π²EI/L² 与热膨胀位移量相除，开出来的是 L 的一次方。

三、π 不是圆整来的，是裂纹张开成半椭圆的面积

很多人以为分母那个 π 是"圆管、圆孔"留下的习惯写法，其实它跟管道毫无关系。

把裂纹两个面张开的位移写出来：弹性力学给的结果是 w(x) = (2σ/E)·√(a² − x²)，x 从 −a 到 a。这是一条椭圆的轮廓：中间张得最开，到两个尖端收成零。现在把应力乘位移沿整条裂纹积分，就是"把这些面重新合上需要做多少功"，也就是释放出来的能量：

∫−aa √(a² − x²) dx = πa²/2

这个积分值就是半径 a 的半个圆的面积。π 是从这个积分里出来的，它描述的是"裂纹张开后是一个什么形状"，不是圆整系数。同样一个 π，也出现在 Inglis 1913 年解椭圆孔的那条公式里：孔边最大应力 σtip = σ(1 + 2√(a/ρ))，ρ 是尖端曲率半径。孔越扁（ρ 越小）尖端应力越高，ρ 趋于零时趋于无穷。

格里菲斯的聪明之处在于：他绕开了这个无穷大。既然尖端应力算不出来，那就不算应力，改算能量——无穷大的应力场里，能量其实是有限的，而且可以用一条简单的抛物线关系表达。这是断裂力学真正的起点。

四、分子那个 2 是哪来的：4 除以 2

把 2Eγ/πa 里的常数一个个对回去，会发现那个 2 既不是"两个 E"，也不是"两个 γ"：

4 来自表面能那一项：上下两个面 × 两个裂尖，每个裂尖推进 da 都要新开两段长 da 的面，所以系数是 4。

2 来自求导：a² 求导掉下来一个 2。4 ÷ 2 = 2，于是写成 2Eγ/πa。

换句话说，只要裂纹换成"单边裂纹"（只有一个裂尖），4 就变成 2，分子也就变成 Eγ——这就是为什么不同文献里这个公式有时写 2Eγ/πa、有时写 Eγ/πa，差别全在"数了几个裂尖"，而不是物理上有什么不同。

五、把钢的数带进去：算出来 16 兆帕，荒谬

取 E = 200 吉帕，铁的纯净表面能 γs ≈ 2 焦耳/平方米，裂纹深度 a = 1 毫米：

σf = √(2 × 2×1011 × 2 / (π × 10−3)) = √(2.546×1014) ≈ 15.96 兆帕

15.96 兆帕意味着什么？意思是"1 毫米深的裂纹，只要 16 兆帕就能把它撕开"。可 1.5 毫米的钢板，光屈服门槛就是 235 兆帕，工作应力还有 10 兆帕上下。如果这个数是真的，凡是带 1 毫米划伤的暖气片，出厂打压就该当场炸开。现实里它连渗都不渗。

结论只有一个：γ 用错了。钢的 γ 不是 2。

六、γ 的第二层：塑性功把它顶高了几个数量级

奥罗万在 1948 年补上了这一层。裂纹往前推的时候，尖端那一小块材料先屈服、被拉长、最后才断开——断开它要的不只是"造两个新面"的 2 焦耳，还要把这一小团金属整体塑性变形一遍的功。所以真正该代入的是

γeff = γs + γp，而 γp ≫ γs

对金属，γp 的教科书量级是 103 焦耳/平方米，也就是表面能的上千倍。同一个公式，只把 γ 从 2 换成 1000，γ 涨了 500 倍，σ 涨 √500 = 22.4 倍，从 15.96 跳到 357 兆帕。这一下就合理了：357 兆帕高于屈服 235 兆帕，说明钢板不是"脆断"，而是先大面积屈服、裂尖钝化，然后才谈断裂。

还能再往上走一层。用本站算过的延性断裂韧度反推：γeff = JIC / 2。Q235 母材的 JIC 在 220 ~ 398 千焦/平方米（见地暖打压气代水那篇），代进去 γeff = 1.1×105 ~ 2.0×105 焦耳/平方米，是纯净表面能的五万到十万倍。此时

σf = √(2 × 2×1011 × 1.1×105 / (π × 10−3)) ≈ 3742 兆帕

3742 兆帕已经远超钢的抗拉强度 400 兆帕——这个数本身就是在宣布判据失效。线弹性断裂力学的全部推导都建立在"裂尖周围是弹性的"这个前提上，而 3742 兆帕的应力根本不可能存在，材料在 235 兆帕就开始屈服了。前提被自己否掉，结论自然不能再用。

γ 取值	物理含义	1 毫米裂纹的 σf	是否可信
2 J/m²	纯净表面能（玻璃、理想脆性）	15.96 MPa	玻璃可信，钢不可信
1000 J/m²	教科书塑性功量级	357 MPa	量级正确
1.1×10⁵ J/m²	由 Q235 的 JIC 反推	3742 MPa	超抗拉，判据失效

七、天花板：平面应变门槛要 453 毫米厚

既然线弹性判据在钢上会失效，那到底什么条件下它还能用？ASTM E399 给了一条硬门槛：试样的厚度必须满足

B ≥ 2.5 × (KIC / σy)²

代入 Q235 的 KIC ≈ 100 兆帕·√米、σy = 235 兆帕：B ≥ 2.5 × (100/235)² = 0.453 米 = 453 毫米。

这个数字很扎人。它说的是：要想测出"真正的、与厚度无关的"平面应变断裂韧性 KIC，试样至少要有 45 厘米厚。暖气片的板厚是 1.5 毫米，比门槛薄了 300 倍。所以任何一份写着"Q235 的 KIC = 100"的资料，其试样都厚得不像话；把它直接拿来判 1.5 毫米的板，等于用另一个世界的参数。薄板只能用平面应力下的 Kc（随厚度变化）、J 积分或 CTOD，这也是为什么钢板硬度报告只写 HV 一样——薄件能做的检测，和厚件能做的检测，从来不是一套。

还有一条能救回线弹性的路：小范围屈服。裂纹尖端塑性区尺寸 rp = K²/(2πσy²)。取实际工况 K ≈ 0.95 兆帕·√米，算得 rp ≈ 2.6 微米，比 1.5 毫米板厚小 577 倍。只要塑性区小到这么一点，外面的应力场仍然是弹性的，线弹性公式就还成立——这就是为什么低应力下的暖气片，还能用 K 去判断会不会裂。

八、两条式子其实是同一条

上面用的是能量式。工程上更常用的是应力强度式：K = Y·σ·√(πa)，Y 是形状因子，中心裂纹取 1，边缘裂纹取 1.12（准确值 1.1215，来自自由表面对裂纹张开的修正级数首项，不是某个解析闭式常数）。

两式并不独立。能量释放率 G = πσ²a/E，而 K² = E·G，把 G 代进去立刻得到 K = σ√(πa)——它们是一件事的两种写法。令 K 等于门槛 Kc = √(2Eγ)，再令它等于 σ√(πa)：

σ√(πa) = √(2Eγ) ⇒ σ = √(2Eγ/πa)

闭合。所以打压漏量那篇用 K = 1.12σ√(πa) 反算临界裂纹长度，和本篇用能量式反算，得到的是同一根门槛，只是单位一个是"兆帕·√米"、一个是"焦耳/平方米"。

九、倍率同异：缺陷长度比材料韧性更贵

把三个变量分别动一下，看 σf 怎么变：

裂纹深度 a 翻 4 倍 ⇒ σ 掉一半（1/√4 = 0.5）。

γ 涨 500 倍（2 → 1000）⇒ σ 涨 22.4 倍（√500）。

材料从玻璃换成钢，E 涨 2.86 倍 ⇒ σ 只涨 1.69 倍（√2.86）。

三者完全不同步。这说明一件事：在格里菲斯公式里，最贵的不是材料有多强，而是有没有那一条 1 毫米的缝。a 在分母、在根号里，动一点点就成倍变化；而想让 σ 涨 22 倍，得把塑性功整个抬上去 500 倍。工程上"提高材料韧性"和"把缺陷做小"哪个更划算，答案在这个式子里是明摆着的。

十、落回暖气片：临界裂纹要 1.1 米

取实际工况：工作压力 0.8 兆帕，水道当量直径 40 毫米，板厚 1.5 毫米，环向应力 σ = pD/(2t) = 0.8 × 40 / (2 × 1.5) = 10.67 兆帕（这组参数与1.5 毫米壁厚那篇一致）。反解临界裂纹：

ac = 2Eγeff / (πσ²) = 2 × 2×1011 × 1000 / (π × 10.67²×1012) = 1.12 米

1.12 米，比整片暖气片的对角线还长。换成 K 形式也一样：一条 2 毫米深的焊趾裂纹，K = 1.12 × 10.67 × √(π × 0.002) = 0.95 兆帕·√米，离 100 兆帕·√米的门槛差 105 倍。

这就是那条最实用的结论：0.8 兆帕的压力，不可能把暖气片压裂，只能把它压渗。鼓起来是弹性变形，裂需要一条长得离谱的缺陷，而渗只需要一条 0.35 微米宽的贯通通道。两者相差的倍数，大到让人误以为"鼓了就是快炸了"。

反过来看为什么铸铁暖气片怕敲：灰铸铁几乎没有塑性，γeff 只有几十焦耳/平方米。同样 0.5 毫米缺陷，取 γeff = 50、E = 120 吉帕，σf = √(2 × 1.2×1011 × 50 / (π × 5×10−4)) ≈ 87 兆帕，低于它 200 兆帕的抗拉强度——先断后屈服。钢板同样 0.5 毫米缺陷算出来是 505 兆帕，高于抗拉 400 兆帕，先屈服后断。同一个公式、同一个几何，只换 γ，行为就翻转了。铸铁件表面那些圆角处的应力集中之所以危险，就是因为它同时压低了缺口的临界尺寸。

十一、现场能用的四条

一、裂纹长度是平方级敏感。a 翻 4 倍，承载力掉一半。发现 1 毫米的焊趾裂纹，比发现一道划痕重要得多，别用"反正不漏"来安慰自己。

二、钢的韧性来自塑性，冷天会掉。温度降到脆性转变区，γp 掉一个数量级，σf 就掉 √10 = 3.16 倍。这就是为什么冬天打压、冬天搬运的破损率明显更高，而不是夏天。

三、焊趾和热影响区是 γ 最低的地方。焊缝的断裂韧度只有母材的四分之一到五分之一，而形状因子又因为"边缘裂纹"从 1 涨到 1.12。强度标号管的是屈服和拉断，管不到这个位置。

四、薄板别套 KIC 表。1.5 毫米比平面应变门槛薄 300 倍，看到"Q235 KIC=100"就直接套用，算出来的安全系数是假的。薄件要判，得用 Kc、J 或 CTOD。

十二、本篇不算

本篇只算准静态、常温、单一裂纹的能量门槛。不涉及：疲劳裂纹在交变应力下怎么一步步长（Paris 公式）、腐蚀介质里 γ 怎么被电化学过程持续压低、以及管道内压波动带来的应力腐蚀开裂。这三件事各有一套常数，都不能用这条式子代。