实心、开口及空心多边形截面的认证Saint-Venant扭转刚度：一种基于闭式解析公式的工程截面分析精确算术原始–对偶方法

《International Journal of Engineering Science》：Certified Saint-Venant torsional stiffness of solid, open and hollow polygonal sections: An exact-arithmetic primal–dual method based on closed-form analytical formulae for engineering section analysis

　　

Saint-Venant扭转刚度是棱柱形结构构件扭转设计的基础，然而对于实践中遇到的大多数截面——冷弯、焊接及组合型材——不存在闭式解，设计依赖于特定族的手册公式（薄壁截面的薄条求和公式、闭口截面的Bredt公式），其精度未被量化。研究人员提出一种原始–对偶流程，将扭转常数Jt（刚度为G Jt）封闭在一个双侧认证区间内：经典的Ritz下界由基于平衡多项式应力场构造的上界补充，且试函数空间的选择使得每个能量积分均化为有理三角形上的有理矩。对于任何具有有理顶点坐标的多边形——即在CAD平面上绘制的任何截面——区间的两个端点及其宽度均为精确有理数，无浮点舍入误差；在有理三角剖分下，同一构造覆盖非凸开口剖面（C型、L型、T型截面）及多重连通空心截面。工程师获得的是一个有保障且可精确复现的参考值：在区间包围不存在闭式的Jt处，其宽度度量残余不确定性，既可验证有限元与边界元扭转求解器，亦可界定设计公式的误差。该方法将单位正方形封闭在围绕Jt=0.140577的宽度约10-4的有理区间内；表明在宽厚比b/t=1时薄条公式高估紧凑矩形达+137%，参考槽钢高估+0.6%至+1.3%，而随着壁厚增加Bredt公式低估矩形管-1.5%至-9.2%——认证区间精确锁定了符号相反的误差。作为副产品，一个精确判据遴选出等边三角形是唯一容许初等Prandtl解（Prandtl solution）的三角形。

研究背景方面，棱柱杆件的均匀扭转问题归结为截面上的Prandtl问题，Saint-Venant扭转常数_Jt作为纯几何属性决定了扭矩与扭率的关系_Mt=_GΘ_Jt，其物理扭转刚度为_{G Jt}。目前仅有圆、椭圆、等边三角形及少数特殊曲线边界截面存在初等闭式解，矩形截面仅存缓慢收敛的Fourier级数解；工程中广泛采用的薄条求和公式与Bredt公式均为渐近设计公式，其精度未量化且无误差界。现有数值方法（有限元、边界元）给出的是数值估计而非带保证的区间，浮点运算引入舍入误差且不可精确复现。由于冷弯、焊接及组合型材在CAD中以有理顶点坐标的多边形表示，亟需一种能提供精确、可认证且可复现的_Jt参考值的理论方法，用以验证数值代码并评估手册公式误差，这构成了本研究开展的动因。研究人员开展了基于精确算术的原始–对偶变分框架研究，对有理顶点多边形截面构造_Jt的双侧认证区间，得出该区间端点为精确有理数、宽度度量能量误差的结论，明确了常见设计公式的系统偏差方向与量级，并证明等边三角形是唯一具有初等Prandtl解的三角形，该研究发表于《International Journal of Engineering Science》，为工程截面扭转分析提供了可认证的精确基准。

作者采用的关键技术方法包括：基于Prandtl薄膜比拟（Prandtl membrane analogy）的变分原理，建立原始（运动许可）与对偶（静力平衡）两类试验函数空间；原始场取为边界函数_B0=∏_i?_i与自由多项式_qp之积_vp=_{B0 qp}以恒等满足零边界条件，对偶场取为_σp=(-x,-y)+?^⊥ψ_p以保持散度-2的平衡条件；利用有理多边形的边线仿射形式（affine form）?_i为正系数的性质，将能量积分化为有理三角形上的多项式矩（polynomial moments）并在精确算术下计算Ritz下界与Prager–Synge型上界；扩展至非凸与空心截面时采用有理三角剖分（rational triangulations）与分片多项式空间，结合Ehrhart拟多项式理论处理有理多胞形矩的精确计算。

研究结果部分保留原文小标题并说明如下：

1. Introduction

研究人员回顾了Saint-Venant扭转问题的Prandtl表述（式1），定义了归一化Prandtl应力函数φ≡?（normalised Prandtl stress function）及其与量纲应力函数Φ的关系Φ=_GΘ?，阐明截面域Ω、边界?Ω、梯度?、Laplacian Δ与旋转梯度?^⊥ψ=(?_yψ,-?_xψ)的物理与几何含义，指出非均匀扭转的翘曲常数_Cw不在本研究范围；介绍了多边形边线的仿射形式?_i（正对内点）及其乘积_B0作为边界函数的构造，说明原始场_vp与对偶场_σp的选取如何在构造上满足运动与静力容许条件，从而将约束变分问题转化为关于_qp与ψ_p系数的无约束问题，且所有积分在有理系数下保持精确算术。

2. State of the art

研究人员梳理了Prandtl设置、经典解及变分原理给出的Ritz下界历史，指出Pólya–Szeg?、Prager–Synge、Trefftz的上界框架已知，但本研究贡献在于将其与有理多胞形面积矩（area moments）结合并在精确算术下实现双侧认证；对比了有限元位移法给出从上收敛的上界、平衡应力法给出下界的现有实践，阐明本研究在任意多项式阶与网格下无需渐近论证即保证不等式成立、端点为精确有理数且无舍入扰动、可逐位复现而不依赖重跑代码的三个差异；指出能量宽度（energy width）具Prager–Synge型恒等式，可作为任意数值解的a posteriori误差估计量；归纳工程收益为求解器验证、手册公式误差控制与安全设计取值三方面，并概述论文从初等解判据、凸截面认证、算例、非凸与空心扩展的组织结构。

3. When does a section have an elementary solution: an exact criterion

研究人员针对三角形截面取边界函数_B0=λ₁λ₂λ₃（λ_i为重心坐标，barycentric coordinates），推导其Laplacian Δ_B0=2[(g₁·g₂)λ₃+(g₁·g₃)λ₂+(g₂·g₃)λ₁]，其中g_i=?λ_i为常向量且|g_i|=1/h_i=a_i/(2A)；通过分析梯度V=?[Δ(_B0)]为零的条件，导出|g₁|=|g₂|=|g₃|等价于边长a₁=a₂=a₃，即仅等边三角形使Δ_B0为常数，从而存在初等Prandtl解?=c_B0；对于非等边三角形无此类初等边界函数，但认证方法仍可提供保证区间，同时指出圆与椭圆对应二次_B0的类似结论。

讨论与结论部分总结：研究人员强调所提原始–对偶精确算术方法不为给出新的_Jt闭式公式，而是将经典变分界限与有理几何矩结合，首次在精确有理数端点下认证_Jt区间及其能量宽度，适用于CAD导出的任意有理顶点多边形（凸、非凸开口如C/L/T型、多重连通空心）。结论表明：单位正方形的_Jt被认证在宽度~10^-4的有理区间内（中心0.140577）；薄条公式在方截面宽厚比1时高估137%，在参考槽钢高估0.6%–1.3%，Bredt公式随壁厚增加低估矩形管1.5%–9.2%，认证区间精确标定这些异号的系统误差；等边三角形被精确判据遴选为唯一具初等Prandtl解的三角形；该认证区间为有限元与边界元扭转求解器提供可复现的精确基准，并为工程手册公式的误差量化提供理论依据。