尽管某些国家或地区的人（例如印度）上厕所时不用厕纸（卫生纸），但对大多数体面的现代人来说，厕纸是不可或缺的生活用品。而人们每天使用的这些厕纸中，超过 75% 的都是由原生树木制造的（图 1）。
为了满足人们对厕纸的需求，全球每天需要生产 2.8 亿卷卫生纸，砍伐 8.2 万棵树。巨大的厕纸需求使自然资源和环境背负上了沉重的负担。
在 2018 年的一项调查中显示，我国人均生活用纸消费量处于全球较低水平，仅为 4.4 公斤（图 2）。而自我标榜环保的美国、德国、英国等西方发达国家，无一例外，都是厕纸消费大国。美国人均年厕纸消费量最高，达到了 12.7 公斤。其次是德国和英国，也都超过了 10 公斤。
相较之下，古代人就比较环保。在古代，纸还不那么地普及，只能用来书写的时候，人们用来擦屁股的物品可谓五花八门（图 3）。富人用羊毛、织物、草纸等柔软物品；穷人则用树叶、竹片、干草、小石块、农作物秸秆、果皮等，附近有水源的则用水洗。
而我们现在普遍使用的卷状厕纸直到 1890 年才出现。现代卷状厕纸通常中间有个空心纸芯（卷筒），上面一圈圈地缠绕着卫生纸。卫生纸宽度为 10 厘米左右，总长度一般在十几到几百米不等。具体长度取决于卷状厕纸的内、外半径和纸的厚度。相信很多人都好奇过一卷厕纸到底有多长？本文将建立多个模型来回答这个问题。
问题
如图 5 所示，已知卷状厕纸外半径为 ，内半径为 ，一层纸厚为 ，厕纸卷的层数 = ( - )/，要求卷状厕纸的总长度 的表达式。
模型
等面积模型
一种最简单、也最容易想到的方法是利用厕纸卷展开前后侧面积相等来计算。
如图 6 左所示，未展开的厕纸卷侧面可近似为一个外半径为 ，内半径为 的圆环，其面积为如图 6 右所示，如果把厕纸卷展开，则侧面圆环变成长度为 ，厚度为 的矩形。因此，侧面积还可以表示为由以上两式得到的面积相等，可得一卷厕纸的总长度为这就算完了？模型和结果也太简单了吧，显然满足不了我们“弃简从繁”的需求！接下来我们尝试一些更“高级”、更有趣的模型。
同心圆模型
除了上面连小学生都会的等面积模型，还可以通过加和每一圈纸的周长来计算厕纸的总长度。为了简化问题，我们将厕纸卷近似为一圈圈的同心圆（图 7）。
由于相邻层圆环半径相差 ，第 1、2、、 层半径分别为：显然，厕纸的总长度可以表示为所有层圆环的长度之和，即：注意，上式给出的结果与等面积模型的完全一样。这是因为同心圆模型与等面积模型没有实质差别（图 8）。
在同心圆模型中，如果每个圆环的周长乘以纸张厚度 就转化为了面积。因此，周长的加和本质上就是面积的加和。
螺旋线模型
上一个模型中，我们假设厕纸卷一圈圈半径的增大是离散的。实际上，卷纸半径随着纸的缠绕角度逐渐增大更符合直观。
数学上将旋转半径随旋转角度匀速地增大而产生的轨迹称为阿基米德螺旋线，自然界中也有很多类似螺旋线的结构（图 9）。因此可将卷纸假设为阿基米德螺旋线。
阿基米德螺旋线在极坐标下可表示为 ，其中 是缠绕角度， 是纸张在该角度的半径。纸卷缠绕的总角度为 ，半径由 变为 ，即 = ， = 。因此 的形式如下将整个螺旋线分成 段弧，第 段弧 的长度为 ，, , 。由于纸张的厚度相对于其长度来讲非常小，弧 对应的角度 较小，曲率半径变化 。当 （即 0），螺旋线的总长度可近似为因此，则整个卷纸的长度可以由以下积分近似我们惊奇的发现，上式给出的结果与前两个模型的完全一样。但实际上，上式给出的并不是螺旋线的精确长度，螺旋线的精确长度可由光滑曲线的弧长公式求得：注意，上式中的 正是我们之前为了近似而忽略掉的曲率半径变化 。上式较为复杂，但考虑到 ，上式可近似为将以上结果与等面积模型（或同心圆模型）结果对比，不难发现，螺旋线模型的结果仅比等面积模型多出一项 ，而多出的这一项显然远小于 2 。因此，螺旋线模型的数值结果应与前两个模型没有明显差别，我们将在稍后的实例计算中验证这一推论。
更懂卷模型
同心圆模型假设一圈圈厕纸的半径是离散的，而螺旋线模型假设厕纸半径是连续均匀变化的，这两个模型都不是厕纸卷的完美近似。真正的厕纸卷并不是以同心圆或螺旋线的方式缠绕在纸芯上的。
如图 12 左所示，将厚度为 的纸的一端固定于圆柱（纸芯）上，然后旋转圆柱，纸逐渐缠绕在圆柱上。圆柱旋转的角度接近完整一周（还差 ）时，纸张的缠绕半径都是相同的，这形成了圆环的一部分。在旋转角度 2 - 角后，纸张与圆柱相切并从该切点直接延伸到接触点。然后以接触点为中心，纸张转过 角。至此，圆柱完成了一圈的旋转，也形成了圆柱上的第一圈纸。
仔细观察，发现第一圈纸可分为三部分：圆环的一部分（绿色），矩形（红色）和扇形（蓝色）。这三部分是由角 所分割开来的， 是直角三角形的一个角，其大小为其对边长度即切点到接触点的距离，可由勾股定理求得：当圆柱继续转动，纸张一层一层的绕在圆柱上，最终形成图 12 右的结果。不难发现，每一层纸都由一段部分圆环（绿色）、一段矩形（红色）和一段两蓝线所夹区域（蓝色）组成。两蓝线的夹角与 相等。最终形成的卷纸也可以分为三部分：圆环的一部分（绿色），矩形（红色）和扇形（蓝色）。绿色区域是部分圆环，其面积为红色区域是矩形，其面积为蓝色区域是扇形，半径为 - - ，其面积为因此，由等面积可以求得卷纸的总长度为考虑到 = 和 = ，上式可化为将以上结果与等面积模型（或同心圆模型）结果对比，不难发现，上式的结果仅比等面积模型多出右侧后两项，而多出的这两项显然远小于 。因此，由上式计算得到的数值结果应与等面积模型没有明显差别，我们将接下来的实例计算中验证这一推论。
结果
为了验证模型，我们以某品牌传统型卷状厕纸为例，计算其长度并与其包装上的标注长度比较。为了应用模型计算该型号厕纸长度，需要知道该型号厕纸卷的相关参数。经测量，该型号未使用全新厕纸的参数（图 13）如下：外径和内径的均值分别为 110 mm 和 50 mm，因此，外半径和内半径分别为 = 110/2 = 55 mm 和 = 50/2 = 25 mm；叠在一起的 10 层厕纸总厚度约为 3.2 mm，因此一层厕纸的厚度 = 3.2/10 = 0.32 mm，层数为 = (55-25)/0.32 94。
将这些参数分别代入同心圆（等面积）、螺旋线和更懂卷模型，得到各模型算得的厕纸长度分别为 23562 mm，23562 mm 和 23556 mm。几种模型给出的结果并没有什么明显差别，都约为 23.6 m。而实际该厕纸包装上标注的长度为 23 m（图 13）。模型计算出的厕纸长度与包装上的标注长度之间相差仅 2.6%。
结论
针对一卷厕纸有多长这个问题，本文分别建立了等面积、同心圆、螺旋线和更懂卷模型。这四个模型依次由易到难（适合不同层次的读者），但都仅需要知道厕纸卷的内外半径和一层纸的厚度就可以计算出厕纸的总长度。通过对各模型的厕纸长度表达式对比发现，四个模型给出的结果相同或相近。为了验证模型，我们将这四个模型应用到某型号的厕纸上。四个模型给出的结果几乎没有差别，都约为 23.6 m，这与该型号厕纸包装上的标注长度（23 m）仅相差 2.6%。
通过比较四个模型，我们发现，尽管螺旋线和更懂卷模型比等面积和同心圆要复杂得多，但复杂的模型并没有得到明显更好的结果。因此，对于仅关心厕纸长度的人来讲，等面积模型就已经是一个足够简单和足够精确的模型。而对于关心厕纸到底是如何卷在纸芯上的人，更懂卷模型是一个很好的演示。
本文针对的问题是二维的，所建立的模型也都是二维的。现实生活中还存在类似的三维问题，如线轴上缠绕的线（图 14）。但实际上，本文的模型也可以稍作扩展来解决这些三维问题。例如对于线轴问题，可将线轴上每一层线想象成一层纸，而每一层线又都是三维螺旋。关于线轴上的线有多长，这里就不再展开讨论，有兴趣的读者可参考本文模型自行建模计算。
附录
点击文章左下角“阅读原文”，在原推文的附录处可以获取本文 PDF 版和附件。
参考资料
[1]
标准起草小组. 《绿色产品评价纸和纸制品》国家标准编制说明, 2017:
[2]
大模头. 为什么我鼓励大家不擦屁股? 2021:
[3]
Martin Armstrong. The u.s. leads the world in toilet paper consumption, 2018:
[4]
Peter R Johnston. How long is my toilet roll?–a simple exercise in mathematical modelling. International Journal of Mathematical Education in Science and Technology, 44(6):938–950, 2013.
[5]UT CALCULUS. Arc length of polar curves:
[6]Jan. Calculating the length of the paper on a toilet paper roll, 2016.:
[7]Tork. Tork conventional toilet roll advanced –2-ply, 2020:
转载内容仅代表作者观点
不代表中科院物理所立场
来源：数学模型
编辑：牧鱼
