薄膜比拟与开口截面扭转：混塔水平缝张开抗扭承载力推导

1. 从一道没有规范条文的验算题说起

上个月校核一个混合塔筒接缝的极限承载力工况，卡在了一个很尴尬的位置：水平缝在极端风况下局部张开，同时塔顶还有不小的扭矩。翻遍手头规范，都是按"全截面受压、全截面传递剪力"来验算，张开状态没有现成公式。最后只能从理论开始重新捋，捋完发现，薄膜比拟、开口截面扭转、混塔水平缝张开抗扭这三件事，看起来分别是弹性力学、钢结构、混凝土塔筒构造三门课里的内容，实际上是同一条逻辑链上的三个节点。

如果你做风电塔架、高耸结构或者薄壁构件设计，这篇文章就是把这条逻辑链完整展开：先讲清楚薄膜比拟为什么能统一理解扭转问题；再说开口截面的抗扭为什么天生就弱，弱在哪个物理量上；最后落到混塔水平缝张开的抗扭承载力推导，给出一个工程上可以直接套用的公式和算例。研究生和刚入行的结构工程师，也能从里面对"从第一性原理出发做推导"有一个具体感知。

所有公式我都给到能照着算一遍的程度，凡是涉及参数取值的地方，都会说明工程上惯用的范围，以及我个人在项目里的取值习惯。如果只想要结果做验算，可以直接跳到第5节，但我还是建议把前几节的物理图像扫一眼，否则公式里每一项代表什么、边界条件在哪，套错了自己都发现不了。

2. 薄膜比拟的核心图像：剪应力就是薄膜山坡的坡度

2.1 扭转方程与薄膜方程的对应关系

圣维南扭转问题里，截面上的剪应力可以用应力函数来描述，记作φ，它满足平面范围内的泊松方程：

∇²φ = ∂²φ/∂x² + ∂²φ/∂y² = -2Gθ

G是剪切模量，θ是单位长度扭转角，边界条件是截面外边界上φ等于零。

这个方程在任意截面形状下都不好解解析解。普朗特尔很聪明，他发现了一个物理类比：取一块薄膜，周边固定在和截面外边界形状一致的边框上，然后施加微小均布压力p。在小挠度范围内，薄膜挠度z满足：

∇²z = ∂²z/∂x² + ∂²z/∂y² = -p/T

T是薄膜张力。两个方程数学形式完全一样，所以可以直接做变量替换。令φ = C·z，代入扭转方程得到C·(-p/T) = -2Gθ，从而C = 2GθT/p。也就是说，扭转应力函数在截面内的分布，等于薄膜在均布压力下的挠曲面形状。

这是整个薄膜比拟的核心：不用解复杂的扭转偏微分方程，只要能想象出薄膜挠曲面长什么样，扭转问题的数学结构就已经被完整复刻了。

2.2 剪应力、扭矩与薄膜坡度、体积的关系

应力函数和剪应力分量的关系是：

τ_zx = ∂φ/∂y，τ_zy = -∂φ/∂x

用薄膜比拟替换之后，剪应力就正比于薄膜的坡度。换句话说，某一点的剪应力大小正比于该处薄膜的坡度，剪应力方向沿着薄膜等高线的法线方向。坡度陡的地方剪应力大，坡度缓的地方剪应力小，薄膜最高点附近坡度为零，对应截面扭心或剪应力为零的位置。

扭矩和应力函数的关系则是一个积分式：

M_t = 2∫∫φ dA

代入薄膜比拟，扭矩正比于薄膜覆盖的体积。你可以想象一座"薄膜山"，山的体积决定总扭矩，山的局部坡度决定局部剪应力。这就是整个类比最直观的两条结论，后面的开口截面抗扭强弱、混塔水平缝张开的退化逻辑，全都可以用这两句话解释。

2.3 为什么这个"老工具"今天仍不过时

有人会问，现在有限元随便一跑就有数值解，为什么还要研究一百多年前的薄膜比拟？我的观点是：数值解给你一个数字，但给不了你判断这个数字是否合理的直觉。薄膜比拟提供的是抗扭问题的几何直觉，概念设计阶段估算截面抗扭能力、校核粗算结果、判断有限元输出的剪应力分布是否可疑，它都很管用。

尤其对于"剪应力流是否连续"这类问题，薄膜比拟能在十几秒内给出定性判断。闭口截面和开口截面，同样外轮廓尺寸，扭转刚度差几百倍，有限元跑出来你很惊讶，但薄膜比拟一看薄膜形状就明白了。这个优势是有限元代替不了的。

3. 开口截面的扭转承载力：薄膜体积一眼看穿的"弱抗扭"

3.1 开口薄壁截面的自由扭转模型

工字钢、槽钢、角钢这类开口薄壁截面，在自由扭转计算里可以看成若干个狭长矩形条的组合。每个窄条宽度为b、厚度为t，且b远大于t，则单条承担的扭矩为：

M_ti = (1/3)Gθ b t³

整个截面的扭转刚度是各矩形条之和：

I_t = Σ (1/3)b_i t_i³

这里有一个经常被忽略的前提："自由扭转"意味着截面可以自由翘曲。实际构件两端如果受约束，会产生约束扭转，截面里出现附加正应力，抗扭能力会比自由扭转公式算出来的高。但做承载力下限估计时，按自由扭转验算是偏保守的，所以开口截面的抗扭校核通常都从这组公式开始。

单个窄条内最大剪应力出现在长边的中点：

τ_max = Gθ t

把总扭矩M_t = GθI_t代进去，整个截面的最大剪应力等于：

τ_max = M_t·t_max / Σ[(1/3)b_i t_i³]

t_max取构成截面中的最厚矩形条的厚度，它决定整个截面的抗扭上限。这套公式在很多钢结构设计手册里出现，但不少手册只给了系数没讲清来源。它就是从薄膜比拟里面"坡度正比剪应力"推出来的。

3.2 数值算例：同尺寸开口截面与闭口截面的差距

为了有直观感受，用一个尺寸接近实际工程的热轧工字钢和闭口箱形截面做对比。工字钢取截面高度500mm、翼缘宽200mm、翼缘厚16mm、腹板厚10mm，忽略圆角。按公式算：

施加10 kN·m的扭矩，最大剪应力：

τ_max = 10×10⁶×16 / 7.02×10⁵ ≈ 227.8 MPa

扭转角：

θ = M_t/(GI_t) ≈ 10×10⁶ / (79000×7.02×10⁵) ≈ 0.18 rad/m ≈ 10.3°/m

这个结果非常惊人。10 kN·m对一根500mm高的构件来说是个不大的扭矩，但开口截面居然产生超过200MPa的剪应力和每米10度的扭转角，材料基本接近屈服了。

再按同样外包尺寸做一个壁厚10mm的闭口箱形截面，用Bredt公式估算：

I_t,closed = 4A_m² / ∮(ds/t)

按中心线围成面积和周长计算，I_t可以到2.86×10⁸ mm⁴左右，大约是同一个工字钢开口截面的400倍。同样的10 kN·m扭矩下，壁内剪应力只有约5 MPa，扭转角小到可以忽略。

3.3 薄膜比拟对此的解释：闭口是气球，开口是水沟

为什么差这么多？薄膜比拟一句话就能解释。闭口截面的薄膜像顶着均布压强的气球，中间整体鼓起来，薄膜覆盖的体积非常大，对应扭转刚度大；剪应力流沿壁厚闭合流动，应力分布相对均匀。开口截面则像在一条细长水沟里铺薄膜，虽然局部坡度可能很陡，但总体体积很小，刚度自然小，局部剪应力反而大。

这个图像对理解风电混塔水平缝的退化机制极其关键：闭口截面的抗扭依赖剪应力流在截面环向路径上"绕一圈"。这个环一旦在某处被打断，闭口就退化成开口，抗扭性能断崖式下跌。混塔水平缝张开后的抗扭问题，本质上就是这个问题。

4. 风电混塔水平缝：从闭口环流到开口化退化的物理图景

4.1 混塔水平缝的构造与传力路径

风电混合塔筒下部是预制混凝土节段，上部是钢塔段。混凝土节段之间通过水平接缝连接，接缝处通常布置竖向预应力筋，把上下节段压紧。接缝面有的做成粗糙界面，有的设剪力键，有的直接平接，靠摩擦传递水平剪力。

正常设计状态下，预应力提供的法向压应力使得接缝面处于全截面受压状态，水平剪力、包括扭矩引起的环向剪力，全部由界面摩擦力承担。扭矩作用下的传力机制，本质上和闭口薄壁截面的剪应力流相似：混凝土筒壁的环向剪应力沿圆周"绕一圈"形成闭合环路。这个环流路径的完整性，是整个抗扭体系成立的前提。水平缝一旦在某处张开，环流路径被切断，抗扭机制立刻从闭口环流退化为局部摩擦加钢筋销栓，抗扭承载能力会明显下降。

4.2 水平缝张开的本质原因：弯矩控制，扭矩叠加

需要澄清一个容易混淆的点：水平缝张开并不是扭矩直接拉开的。接缝面的法向应力主要由轴力和弯矩决定。正常工况下，预应力产生的压应力大于使用荷载产生的拉应力，接缝保持受压闭合。到了极限工况，比如极端风况或地震作用，倾覆弯矩产生的拉应力超过预应力压应力，受拉侧接缝就张开了。

所以更准确的说法是：弯矩导致水平缝出现"张开区"，而在这个已经张开的几何状态下，截面还要同时承担扭矩。设计关心的问题是，水平缝张开之后，抗扭承载力还剩多少。工程上正是这种"弯矩先劈开、扭矩同时作用"的复合受力让设计变得棘手。

4.3 张开区与闭合区的划分：一套可操作的几何参数

定义张开区为接缝面中法向应力为正（拉应力）的扇区，设其总圆心角为2α，闭合区为其余扇区，总圆心角为2π-2α。α可以从截面应力分布直接求出。

环形截面在轴力N和弯矩M_y共同作用下的法向应力为：

σ(z) = -N/A + M_y·z/I_y

z是沿弯曲方向到截面中心的投影距离，受拉侧取正。I_y是截面对弯曲轴的惯性矩，环形截面取：

I_y = π(R⁴-r⁴)/4

令σ=0可以得到中性轴位置z0，从而确定张开角：

cosα = z0/R，其中 z0 = (N/A)·(I_y/M_y)

当z0≤-R时全截面受压，α=0；当z0≥R时全截面受拉，α=π。这个关系把宏观内力（轴力、弯矩）和张开的几何状态联系起来了，工程上很实用：拿到一个极限工况的内力组，就能算出水平缝张开到了什么程度。

这里要特意提醒一句：α不是独立设计参数，它由特定工况下的N、M决定，不能单独取最大弯矩对应的α、再叠加另一个工况的最大扭矩。必须用同一个荷载组合的内力组。

5. 水平缝张开状态的抗扭承载力推导：摩擦项与销栓项

5.1 基本假定

推导之前先把假设写清楚，后面用公式才有底：

1. 接缝面为圆环截面，外半径R、内半径r，壁厚沿环向均匀。
2. 水平缝因弯矩作用已存在对称张开区，总圆心角为2α。
3. 闭合区内混凝土界面处于受压状态，摩擦剪应力按 τ_f = μσ_n 计算；σ_n取预应力提供的平均法向应力。
4. 张开区内混凝土界面脱离接触，摩擦贡献为零；但贯穿接缝的预应力筋仍有销栓作用，每根销栓抗剪承载力为V_s。
5. 不计混凝土本身抗拉强度，不计钢筋的斜向受拉分量，结果偏保守。

在这些假定下，水平缝张开状态的总抗扭承载力由两项组成：

M_u = M_f + M_s

M_f是闭合区的摩擦抗扭，M_s是张开区预应力筋的销栓抗扭。

5.2 第一项：闭合区摩擦抗扭的推导

在闭合区任取一个微元面dA，其半径为ρ，摩擦剪应力为μσ_n，沿环向。微元对截面形心的抗扭贡献为：

dM_f = μσ_n · ρ · dA

极坐标下dA = ρdφdρ，所以：

M_f = μσ_n ∫_r^R∫_闭合区 ρ² dφ dρ

角度积分范围是闭合区的总圆心角2π-2α，因此：

M_f = μσ_n · [(R³-r³)/3] · (2π-2α)

为了对照，先列出完全未张开时的抗扭承载力：

M_0 = (2π/3)μσ_n(R³-r³)

于是闭合区摩擦项可以写成：

M_f = M_0 · (1-α/π)

这是核心式子之一。它说明摩擦抗扭能力大致随张开角线性下降。张开区总圆心角60°，摩擦抗扭剩5/6；张开区总圆心角180°，摩擦抗扭只剩一半。

M_0这个完好状态值可以改写为更工程化的形式：

M_0 = μN_pre·R_eq，R_eq = (2/3)(R³-r³)/(R²-r²)

R_eq是等效摩擦半径，对薄壁环形截面约等于内外半径的平均值。于是手算逻辑变成：总预压力乘以等效半径得到抗扭力臂，再乘摩擦系数，就是完好状态的抗扭能力。

5.3 第二项：张开区预应力筋的销栓抗扭

张开区混凝土面已经脱开，但竖向贯穿的预应力筋还在。接缝平面内出现相对滑移时，预应力筋相当于被横向剪切，这就是钢筋混凝土里的销栓作用。单根钢筋的销栓抗剪承载力近似为：

V_s = α_s·f_y·A_s

A_s是单根钢筋面积，f_y是屈服强度，α_s是销栓折减系数。销栓作用受周围混凝土局部承压控制，α_s不能取得太高。工程上一般在0.10到0.20之间，我习惯取0.15；有可靠灌浆套筒约束时可以适当上调，但没有试验依据前不建议超过0.25。

预应力筋沿圆周均匀布置时，总数N_s中落入张开区的根数为：

n_open = N_s·(α/π)

所有张开区钢筋的销栓力近似作用在预应力筋的布置半径ρ_s处，沿环向，于是：

M_s = V_s·ρ_s·n_open = V_s·ρ_s·N_s·(α/π)

张开角越大，进入张开区的钢筋越多，销栓项贡献越大。这和摩擦项随α增大而降低正好是此消彼长的关系。

5.4 总承载力公式与两种压应力模型

把两项合并，得到水平缝张开状态下的抗扭承载力估算公式：

M_u = M_0(1-α/π) + V_s·ρ_s·N_s·(α/π)

这个式子结构很简洁：M_0是完好状态抗扭能力，第二项是剩余机制的补充。

但这里有一个隐藏的坑。前面假定σ_n取原始平均预应力、不做重分布，得到摩擦项随张开角线性下降。实际上，当一部分截面张开后，预压力会重新分配到残余的闭合区面积上，闭合区平均压应力会升高。如果按压应力完全重分布计算：

σ_n' = N_pre/[π(R²-r²)(1-α/π)] = σ_n/(1-α/π)

M_f' = μσ_n'·[(R³-r³)/3]·(2π-2α) = M_0

也就是说，如果压应力在闭合区均匀重分布并且摩擦系数保持不变，摩擦抗扭项不会下降，始终等于M_0。初看反直觉，理解起来也简单：承压面积变小了，但单位面积摩擦力变大了，摩擦总力等于μ乘以总预压力，总预压力没变，摩擦抗扭总量自然不变。

这里的关键是"压应力均匀重分布"这个假设是否成立。实际环形截面的应力重分布并不均匀，越接近张开边界应力集中越明显；而且闭合区面积缩到很小时，局部压应力会逼近混凝土承压极限，摩擦滑移关系也不再是简单的库仑模型。所以当α接近π/2以上时，重分布模型已经明显失真。

工程上建议这样处理：

下限：不重分布模型，M_u = M_0(1-α/π)+M_s，用于设计校核。 上限：完全重分布模型，M_u = M_0+M_s(α/π)，用于理解机理、评估上限。 真实承载力落在这两条线之间。方案阶段用下限校核，不足时再做精细化有限元验证。

这个"两条线包络"的处理方式，比给一个看起来精确的单一安全系数更符合工程实际。

5.5 完整算例：一座5MW级混塔接缝

给一个能直接复算的完整算例。某混塔水平缝：

• 外半径R = 3.0m，内半径r = 2.7m，壁厚300mm
• 有效预压力N_pre = 20000kN
• 接缝面摩擦系数μ = 0.6
• 预应力筋总数N_s = 64根，单根有效面积A_s = 150mm²，f_y = 1670MPa，布置半径ρ_s = 2.82m
• 销栓折减系数α_s = 0.15
• 极限工况由弯矩确定的张开角α = π/6，即总张开圆心角60°

完好状态：

A = π(R²-r²) = π×1.71 = 5.373m²

σ_n = 20000/5.373 ≈ 3.72MPa

M_0 = 0.6×3.72×10⁶×(7.317/3)×2π ≈ 34.2 MN·m

用下限模型，张开60°时的摩擦项：

M_f = 34.2×(1-1/6) = 28.5 MN·m

单筋销栓承载力：

V_s = 0.15×1670×150 = 37575N ≈ 37.6kN

销栓项：

M_s = 37.6×10³×2.82×64×(1/6) ≈ 1.13 MN·m

总抗扭承载力下限：

M_u = 28.5 + 1.13 = 29.6 MN·m

继续看张开角更大的情形。总张开圆心角180°，即α=π/2时：

M_f = 34.2×0.5 = 17.1 MN·m

M_s = 37.6×10³×2.82×64×0.5 ≈ 3.39 MN·m

M_u ≈ 20.5 MN·m，比完好状态下降约40%。即使极端到全截面张开，摩擦项归零，销栓项也只有6.78 MN·m，大约是完好状态的20%。

这个算例给出了一个重要工程结论：水平缝张开后的抗扭承载力主要靠闭合区摩擦，预应力筋的销栓只能作为安全裕度，不能作为主要抗扭机制。设计时如果濒临张开状态，首先要控制张开角，其次再谈加强预应力筋。

6. 工程应用边界、参数取值与完善方向

6.1 这个公式能用在哪，不能用在哪

这个简化公式适合用于：混塔方案阶段的截面抗扭校核、水平缝张开状态下承载力的快速估算、有限元结果合理性判断。对预制混凝土塔筒节段之间的干接缝尤其适用，因为界面行为接近摩擦模型。

不适用范围也要说清楚：带复杂剪力键、接缝面开槽灌浆、接缝材料有明显粘结强度的节点，不能直接用。这些节点的抗剪机制除了摩擦还有机械咬合和粘结，摩擦模型会低估实际承载力，但作为下限问题不大。反过来，如果接缝面有长期疲劳磨光、或者做了降低摩擦的处理，摩擦系数要下调，否则公式可能高估，必须按实测数据修正。

另外一个限制是公式只给承载力，不给变形。如果关心张开状态下的扭转刚度，需要建立滑移与剪应力的本构关系，把完整的M-θ曲线求出来，这超出了简化模型的范畴。实际工程中，如果抗扭刚度影响整机动力学特性，建议直接做带接触的实体有限元分析。

6.2 参数取值的经验与坑

摩擦系数μ是整个公式里最敏感的参数。混凝土与混凝土之间的干接缝，粗糙度良好时取0.6到0.7没有问题，但要区分界面状态：粗糙界面，凹凸深度大，取高值；平滑界面甚至高精度水磨面，摩擦系数要降到0.4到0.5。长期风致微振会让界面磨光，动摩擦系数比静摩擦低，疲劳验算里要用动摩擦值。如果水平缝存在灌浆或环氧粘结层，界面失效模式会从摩擦滑移变成粘结剪切破坏，摩擦模型本身就不适用了。

销栓折减系数α_s取0.15时已经包含了局部混凝土压碎的影响。预应力筋如果在接缝附近有波纹管或灌浆缺陷，销栓承载力还会进一步下降。对于后张预应力体系，钢绞线在波纹管内的初始偏心也会影响销栓行为。这些不确定因素都要在安全系数里消化掉。

6.3 和有限元、试验的衔接思路

做精细有限元分析时，建议在接缝面设置接触对，法向采用硬接触并允许张开，切向采用带摩擦系数的库仑摩擦模型，预应力筋用工况初始应力。后处理重点看两个量，一是水平缝张开角的实际值与α线性公式算出来的差异，二是张开区预应力筋销栓剪力与V_s估计量是否一致。这两点直接校验简化公式的误差来源。

试验方面，如果有条件做缩尺接缝扭剪试验，建议把张开角做成控制变量，测量不同α下的M_u。试验结果大概率落在我前面说的下限和上限包络之间，α越大，越靠近下限。

6.4 几点个人体会

遇到规范条文覆盖不了的问题，回到基础力学原理去重新推一遍，这种方法在工程上其实很高效。薄膜比拟这种"老工具"的价值恰恰在于它把抽象方程变成几何图像，画出来、想明白，很多问题就没有那么难了。

推导承载力公式时，明确写出假设比堆公式重要得多，因为后面所有讨论都围绕假设展开。工程上准备两个模型也同样重要：一个保守下限用于设计，一个理想上限用于理解机制，两者之间的差距就是风险所在，也是精细化分析的工作空间。

我在实际项目里用这套公式做过初步判断，结果方向都是准的。但我不会建议在一个风电塔筒的关键节点上只依赖手算公式下最终结论，重要节点最终都应该有精细化分析或试验支撑。反过来，如果简化结果和精细结果差得很远，先不要急着质疑有限元，回头检查α取值和摩擦系数是否统一，大多数时候问题出在模型参数传递上，而不是推导本身。