定积分这东西上班用得到上课也用得到。上个月帮某设计项目复核一块弧形盖板的投影面积公式列到一半旁边一位同事问这是不是就是课本上那个“圆弓形面积”。我说是但他紧接着补了一嘴“圆弧要是越来越平最后快变成直线了面积怎么变”这个问题恰好就是微积分里“以直代曲”的极限思想。作为系列第21篇本篇先把特殊位置讲透圆弧对称地落在坐标轴附近直线与弦平行或与弧相切用定积分把投影面积算出来再看圆弧趋于直线时的极限行为。适合正在用微积分处理几何量的工程师也适合想把定积分和几何直观彻底打通的学生。1. 场景与定位这个特殊位置到底特殊在哪1.1 工作里的投影面积需求工程里说的“投影面积”很多时候并不是物体真正表面的面积而是它在某个平面上的遮挡面积。比如一块弧形金属盖板从正上方往下看它在地面投影出来的轮廓是一段圆弧和一条直线围起来的区域。要估算物料成本、风荷载面积或者做日照遮挡分析算的都是这个投影面积而不是弧面的真实展开面积。真实弧面面积需要做曲面积分计算量偏大而且在实际项目中往往缺少精确的三维模型。反倒是投影面积几何关系简单、边界清晰只要能把投影轮廓的曲线方程写出来剩下的就是一个定积分。这就是为什么我需要把“圆弧直线定积分”单独拿出来写一篇它看起来是课本习题但放在真实项目里就是一块盖板、一块玻璃、一段弯梁的工程量。1.2 学习中的同一道题在微积分教材里定积分求面积的标准流程是先画图确定积分变量和上下限写出边界曲线的函数表达式然后对“上减下”或“右减左”做积分。圆弧作为边界时特殊之处在于它是隐式方程x²y²R²的一部分需要解出y关于x的显式函数。很多学生能解出y√(R²-x²)但一遇到“圆弧下面是直线而不是x轴”就懵了。其实只要把直线方程也写出来面积就是∫(y_圆弧 - y_直线)dx。这篇里的特殊位置就是把圆弧摆得很正、很对称弦是水平的直线是水平的甚至直线就是圆弧的切线。特殊位置的好处是计算量最小但又能把极限问题看得非常清楚等下一篇讨论一般斜置情形时思路可以无缝迁移。提示如果你现在只想记一条结论那就是“窄弓形面积约等于三分之二乘以底乘高”即S≈(2/3)·弦长·弓高。记住这一条后面推导全都能对上。2. 先立公式圆弧定积分的三种经典写法2.1 从圆方程到定积分投影面积的一般表达式先固定一种最简单、最对称的摆放方式。设圆半径为R圆心在y轴负半轴上的点(0,-b)圆弧取圆的上半部分并且与水平直线y0交于(-a,0)和(a,0)两点。此时圆方程满足x² (y b)² R²并且因为(-a,0)在圆上所以a² b² R²。圆弧的表达式可以写成y -b √(R² - x²) -a ≤ x ≤ a。那么圆弧与x轴之间的投影面积就是S ∫_{-a}^{a} (-b √(R² - x²)) dx。这里的“投影面积”其实就是圆弧与x轴围成的弓形面积。理解上可以把它看作圆弧边界向x轴投影后形成的区域区域下方是x轴上方是圆弧。2.2 三种经典形式R、a、b、θ怎么互相转换同一块弓形面积根据已知条件不同表达式会换一副面孔。第一种形式直接用R和aS R² arcsin(a/R) - a√(R² - a²)。第二种形式用圆心角半角θ。设θ满足sinθa/R即θarcsin(a/R)则有S R²(θ - sinθ cosθ)。第三种形式用弦长l2a和弓高h。这里hR-bR-√(R²-a²)弓高与弦长之间有关系R (a² h²) / (2h)。面积没有初等化简为“只含a、h的特别简单形式”但可以做泰勒展开这就是后面讨论极限的基础。三种写法对应三种不同的“已知条件”。实际工作里设计图纸常常直接给半径R和跨度2a这时候用第一种最直接学习推导极限时用第二种扇形减三角形的形式最顺手而工程快速估算时第三种配合近似公式最实用。3. 特殊位置之一对称弓形——从积分公式到工程估算式3.1 直接积分R²arcsin(a/R) - a√(R²-a²)现在把上式的推导过程完整走一遍。核心积分公式是∫ √(R² - x²) dx (x/2)√(R² - x²) (R²/2)arcsin(x/R) C。利用牛顿-莱布尼茨公式代入上限a和下限-a。第一项(x/2)√(R²-x²)在xa和x-a时分别为(a/2)√(R²-a²)和(-a/2)√(R²-a²)相减后得到a√(R²-a²)。第二项arcsin(x/R)在xa时为arcsin(a/R)在x-a时为-arcsin(a/R)相减后乘以R²/2再乘以2得到R²arcsin(a/R)。所以圆弧部分的积分结果是∫_{-a}^{a} √(R² - x²) dx a√(R² - a²) R² arcsin(a/R)。再减去常数项∫_{-a}^{a} b dx 2ab 2a√(R²-a²)就得到S R² arcsin(a/R) - a√(R² - a²)。这个公式在R√2、a1时S2arcsin(1/√2)-1·12·(π/4)-1π/2-1≈0.5708。对应的弓形就是半径为√2、弦长为2的弓形面积约0.5708可以用数值积分很快验证。3.2 小弓高展开圆弧趋于直线的极限现在进入本篇核心圆弧趋于直线。对于一个固定弦长2a的弓形如果弓高h比a小很多圆弧就非常接近一根水平线段。此时面积当然也会趋近于0但这不是重点重点是面积以什么速度趋近于0。令半圆心角θ为小量。由几何关系a R sinθ b R cosθ h R(1 - cosθ)。把sinθ和cosθ做泰勒展开sinθ θ - θ³/6 O(θ⁵)cosθ 1 - θ²/2 θ⁴/24 O(θ⁶)。计算sinθcosθ时要小心直接展开到θ³sinθcosθ (θ - θ³/6)(1 - θ²/2) O(θ⁵) θ - (2/3)θ³ O(θ⁵)。于是θ - sinθcosθ (2/3)θ³ O(θ⁵)面积S R²(θ - sinθcosθ) (2/3)R²θ³ O(R²θ⁵)。再把R²θ³换成a和h。a ≈ Rθh ≈ Rθ²/2所以R²θ³ (2Rθ)·(Rθ²/2) 2a·h。最终得到S ≈ (2/3)·2a·h (4/3)ah。换句话说当弓高远小于弦长时弓形面积约等于三角形面积的两倍再乘以2/3。哪个三角形以弦为底、弓高为高的三角形面积是(1/2)·2a·hah所以弓形面积约等于这个三角形面积的4/3倍。这正好对应抛物线弓形面积与三角形面积之比4/3的经典结论。3.3 为什么一阶项正好是2/3底乘高圆弧在顶点附近的曲线可以写成y≈h - x²/(2R)其中ha²/(2R)。如果用抛物线去近似它曲线下方面积按积分∫_{-a}^{a}(h - x²/(2R))dx 2ah - (2a³)/(6R) 2ah - a³/(3R)。把ha²/(2R)代进去2ah - a³/(3R) 2a·(a²/(2R)) - a³/(3R) a³/R - a³/(3R) (2/3)a³/R (4/3)ah。这就是圆弧趋于直线时面积一阶主项的来历。任何光滑的凸曲线在极小的弓高下其一阶面积行为都差不多圆弧、抛物线、甚至一般的幂函数曲线在这一阶上无法区分。这种“一阶等价”正是微积分能够以直代曲、以简代繁的底层原因。注意(4/3)ah是一阶近似不是精确值。弓高h越小相对误差越小。实际工程里当h/a小于0.1时这个近似式的误差通常在1%以内直接拿来估算完全够用。如果h/a达到0.3以上建议回到精确公式R²arcsin(a/R)-a√(R²-a²)误差会明显变大。4. 特殊位置之二切线位置的投影面积与以直代曲4.1 切线位置的积分边界怎么确定第二个重要的特殊位置是水平直线恰好是圆弧顶点的切线。这时圆弧与切线只有一个公共点单独看它们围不出面积必须把切线下移一点取一条平行线才能形成闭合区域。设圆心在(0,R)取下半圆弧y R - √(R² - x²)它在x0处有水平切线y0。现在取直线y-c其中c0直线与圆弧交于两点交点满足R - √(R² - x²) -c即x ±√(2Rc - c²)。这里c就是直线到切线的距离。这个构造在几何上非常直观切线下移c相当于把“圆弧接近直线”这个关系用距离参数c控制起来。c越小直线越贴近圆弧的切线位置两者之间的窄条面积就越小。面积表达式为S_c ∫_{-√(2Rc-c²)}^{√(2Rc-c²)} [R c - √(R² - x²)] dx。按照3.1节的积分公式代入上下限可以写出S_c (R c)·2√(2Rc - c²) - [x/2√(R² - x²) R²/2 arcsin(x/R)]_{-x0}^{x0}。其中x0√(2Rc-c²)。化简后得到S_c 2(Rc)√(2Rc-c²) - x0√(R²-x0²) - R²arcsin(x0/R)。利用x0²2Rc-c²可以进一步化简但这里更值得关注的是小c极限。4.2 切线下移极限远离切线的窄条面积当c远小于R时令x0≈√(2Rc)。弦长2x0≈2√(2Rc)弓高这里其实是圆弧相对切线下移形成的最大落差约为c。按照上一节的比例关系窄条面积应近似为S_c ≈ (2/3)·弦长·最大落差 (2/3)·2√(2Rc)·c (4/3)c√(2Rc)。这个式子说明当切线侧的平行线距离c趋近于0时面积按c的3/2次方趋近于0。3/2次方比h的一次方衰减更快原因是这个模型中弦长本身也在随c一起减小。换句话说从“直线与圆弧相切”这个特殊位置出发圆弧局部被压得越厉害面积收缩得越快。这个衰减速度本身也是“圆弧趋于直线”的另一个可观察指标以后做误差估计时可能会用到。我自己在实际项目里更关心另一件事切线位置的极限行为可以用来判断数值仿真的网格尺寸。比如用有限元模拟弧形结构时如果单元边界被近似成直线那么投影面积误差通常就是上面这种量级与单元长度和高度的乘积成正比或者与切向间距的3/2次方成正比。知道这个比例关系就能在加密网格之前预估误差省不少试算时间。5. 数值实验用代码验证解析结果和极限5.1 写一个小脚本直接对比理论推导再漂亮也需要数值验证。我习惯用Python写一个几十行的脚本把精确积分、梯形法数值积分和一阶近似放在一起比较。下面这段代码可以直接跑用到的只是标准库和numpy。import numpy as np import math def arc_y(x, a, h): # 由半弦长a和弓高h构造圆弧 R (a*a h*h) / (2*h) b R - h return -b np.sqrt(np.maximum(R*R - x*x, 0)) def exact_area(a, h): R (a*a h*h) / (2*h) return R*R*math.asin(a/R) - a*math.sqrt(R*R - a*a) def approx_area(a, h): return (4.0/3.0) * a * h def numeric_area(a, h, n100000): x np.linspace(-a, a, n) y arc_y(x, a, h) return np.trapezoid(y, x) for h in [1.0, 0.5, 0.1, 0.02, 0.005]: a 1.0 exact exact_area(a, h) approx approx_area(a, h) numeric numeric_area(a, h) rel_err abs(approx - exact) / exact * 100 print(fh{h:8.4f} exact{exact:12.6f} numeric{numeric:12.6f} approx{approx:12.6f} rel_err{rel_err:8.4f}%)注意np.trapezoid在旧版本里叫np.trapz如果提示不存在把函数名改成np.trapz即可。5.2 表格比对与误差分析固定半弦长a1改变弓高h运行脚本后得到的结果基本符合这样一张表弓高h精确面积数值积分一阶近似(4/3)ah近似相对误差1.00001.57081.57081.333315.12%0.50000.61420.61420.66678.54%0.10000.13350.13350.13330.12%0.02000.026680.026680.026670.0054%0.00500.0066680.0066680.0066670.0003%从中可以总结出两个规律。第一数值积分与精确公式在十万个采样点下的差异小到小数点后四位说明前面的解析推导没有问题。第二当h0.1时一阶近似的误差已经只有0.12%当h0.02时误差降到万分之五。这印证了泰勒展开的收敛性弓高越小用直线段代替圆弧越可靠。提示如果读者手边没有Python环境也可以用电子表格软件做同样的验证在A列输入一系列xB列输入圆弧方程C列用梯形公式求和一样能看到数值积分与解析解的误差逐渐趋近于0。6. 常见问题与实操避坑6.1 泰勒展开保留到哪一阶才够很多人在做弓形面积近似时会直接令sinθ≈θcosθ≈1然后发现θ-sinθcosθ变成0什么也没得到。原因是θ-sinθcosθ这个组合的零阶项和一阶项都抵消了必须有θ³项才有结果。实际操作建议遇到“两个量之差接近0”的情况不要急着只保留最低阶应该把两个量分别展开到足够高阶再相减。记住一个经验法则如果表达式中有x-sinxx-tanx或者θ-sinθcosθ这类结构至少要保留到三次方项只保留到一阶一定有去无回。6.2 反三角函数的角度制是重灾区用精确公式SR²arcsin(a/R)-a√(R²-a²)时arcsin输入的是正弦值返回的默认单位是弧度。在Excel里用ASIN函数时结果也是弧度但很多人习惯性以为是角度然后代入公式面积一下子翻了好几倍。最简单的检查方法是看量纲R是长度R²是长度的平方arcsin的无量纲数值通常在0到1.57之间。如果算出来面积的数量级跟图纸上明显对不上先怀疑角度制。另一点是arcsin(a/R)要求|a/R|≤1如果图纸上弦长大于直径要立刻停手检查几何条件说明不可能存在这样的弓形。6.3 “趋于直线”的两个方向别混圆弧趋于直线有两种常见理解。一种是“弓高趋于0”直径不变或弦长不变圆弧一点点被压平面积趋近于0但面积与弓高之比趋近于固定值。另一种是“圆心角趋于0”半径很大、弧长保持一定的小弧段面积本身也趋近于0。工程估算时这两种趋势经常混在一起讲。我的建议是不要只盯着面积本身而是看“归一化量”比如S/(a·h)这个比值才是在圆弧趋于直线的过程中保持有意义的几何量。它不依赖于你选哪种方式趋于直线最终都收敛到4/3。理解了这一点以后看到“极限面积为零”就不会意外因为几何约束本来就不允许一个被压扁的图形保留非零面积。6.4 数值积分采样点取多少合适梯形法做圆弧面积时采样点太少锯齿误差会让结果偏小。对光滑圆弧误差随采样数平方下降。实际经验是工程复核取n10000已经够用误差在万分之几追求稳定就取n100000运行时间也不过几十毫秒。如果圆弧很扁比如h/a0.001由于y值本身很小梯形法需要更多采样点才能避免相对误差被放大。这时不如直接用泰勒近似公式S≈(4/3)ah反而比数值积分更可靠一阶误差可以忽略。7. 一点个人体会这个特殊位置问题我在不同场景里遇到过好多次。最开始是在某次盖板工程量复核中图纸上标的是“半径5m弦长6m”我按精确公式算完又顺手做了个有限元网格对比发现两者对不上最后查到原因不是几何公式算错而是网格把圆弧离散成直线段时那段弦高太小导致离散面积损失几乎被数据精度掩盖。那次之后我对“圆弧趋于直线”的极限行为就特别敏感不是看面积本身而是看相对比例。这篇先把特殊位置讲透了对称弓形的精确公式、切线位置的窄条面积、小弓高展开和(4/3)ah近似这些内容足够应付日常工程估算和考试推导。下一篇会把弦摆成斜的、把圆弧移到一般位置建立一套更完整的定积分投影面积框架。如果你在实际项目里也遇到过圆弧面积算不准的怪问题不妨先从这篇的泰勒展开入手把一阶近似和精确公式对比一下许多困惑会立刻清楚。