最近在折腾可折叠结构的设计方案连续测了几组不同折痕分布的样品发现一个非常有意思的规律折叠压缩率并不是均匀地分布在整个结构上而是按三角函数曲线的形态在空间里展开。这个发现直接影响了我后续的尺寸参数选取和折叠角设计今天把整个推导过程和实测数据整理出来希望给做柔性结构、折纸工程或者缓冲材料的朋友一些参考。这篇内容主要围绕三个问题展开为什么折叠压缩率会按三角函数曲线分布、如何从几何关系推导出这条曲线、以及这条曲线在工程上到底能拿来干什么。适合正在做可折叠结构建模、软体机器人关节设计、或者研究多孔材料压缩性能的读者也适合对几何与力学耦合关系感兴趣的入门者。1. 从一条折痕说起折叠压缩率到底是怎么定义的1.1 先统一一个概念避免后面绕晕很多刚接触折叠结构的同学会把折叠压缩率和材料的压缩率混在一起。材料的压缩率通常指材料本身在外力下体积缩小的比例比如弹性体受压时的体积变化。但折叠压缩率不一样它描述的是在一个由折痕和平面面板组成的结构中结构整体沿着某个方向被压缩的程度与局部折叠角度变化量的比值关系。举个例子一个手风琴式的折叠纸板初始展开时总长度是L₀当你沿轴向压它长度变成L₁那么整体压缩率定义为[ \eta \frac{L_0 - L_1}{L_0} ]这个定义看起来很简单但如果折痕不是等间距、等角度的问题就出现了结构内部各个局部的折叠角并不相同每个单元格对整体压缩的贡献也不同。也就是说整体压缩率不是局部的简单平均而是沿着结构位置呈现出某种分布。我一开始以为只要按线性分布去设计折痕间距就够了后来把测得的局部位移数据画出来发现完全不是线性而是很接近三角函数曲线。这才意识到折叠结构的几何约束本身就带有三角函数的影子。1.2 一个基础折叠单元的几何模型为了看清楚本质我取一个最常用的折叠单元来分析两个平面面板共享一条折痕折痕两侧的面板长度都是a折叠角面板与水平面的夹角为θ。当结构受到垂直方向压缩时θ会减小结构在水平方向上的投影长度L与θ的关系是[ L 2a \cdot \cos\theta ]这个式子太典型了。水平投影长度直接等于两条面板的投影之和每段投影都是a·cosθ。如果初始折叠角是θ₀被压缩到θ₁那么局部压缩率在这个单元内就是[ \eta_{local} \frac{2a\cos\theta_0 - 2a\cos\theta_1}{2a\cos\theta_0} 1 - \frac{\cos\theta_1}{\cos\theta_0} ]看到没有局部压缩率是夹角余弦的比值天然就是三角函数的形式。当θ₀和θ₁很接近时还可以近似为负的正弦函数乘以角度变化量。这就是折叠压缩率按三角函数曲线分布的第一个来源。1.3 为什么是分布而不是单一值有人会问如果一个折叠结构所有单元都完全相同那局部压缩率不是处处相等吗为什么还要说分布关键在于实际可折叠结构很少是单一均匀单元组成的。更多的应用场景是折痕深度、面板长度、或者初始折叠角沿着结构连续变化。比如一个呈波浪形的薄板其表面起伏高度沿着长度方向按正弦函数变化那么每个位置处的局部等效折叠角就是位置的函数压缩率自然也随位置呈三角函数曲线分布。我在设计一个波纹形缓冲层时就是让折痕的深度按余弦函数从中心向两侧衰减。这样中心区域的折叠角变化率最大局部压缩率最高越往边缘折痕越浅局部压缩率越低。出来的效果完全符合预期同样的总压缩量下应力分布比均匀折痕结构平缓得多。2. 推导过程压缩率函数为什么是余弦曲线的叠加2.1 连续化处理把离散折痕变成连续曲面离散的折痕单元分析起来清晰但工程上很多时候需要连续模型。假设一个折叠板的截面轮廓可以用函数z(x)描述如果它是周期性三角函数比如[ z(x) A \cos\left(\frac{2\pi x}{\lambda}\right) ]其中A是波幅λ是波长。这个曲面的切线斜率就是z(x)也就是[ z(x) -A\frac{2\pi}{\lambda} \sin\left(\frac{2\pi x}{\lambda}\right) ]而折叠角θ与切线斜率的关系为[ \theta(x) \arctan\left(|z(x)|\right) ]当压缩发生时A或λ发生变化导致局部压缩率ε(x)发生变化。如果我们假设压缩主要是改变波幅A而λ不变那么θ(x)的变化量就会正比于正弦函数而局部压缩率又与cosθ相关最终ε(x)会呈现出cos函数项与sin函数项的混合但本质上仍然是一种三角函数曲线分布。2.2 从变形几何推导压缩率分布的表达式更严谨地我们考虑一个微小的水平段dx。该段结构在垂直方向的投影高度为z(x)当结构受到垂直压缩后z(x)变为z(x) k·z(x)其中k是1的系数。那么该微段沿板面方向的长度变化率决定了局部压缩率。板面的微弧长是[ ds \sqrt{1 (z(x))^2} dx ]压缩前后微弧长不变假设面板不可伸长那么水平方向的压缩率就可以近似表达为[ \eta(x) 1 - \frac{\sqrt{1 (k z(x))^2}}{\sqrt{1 (z(x))^2}} ]当斜率较小时用泰勒展开可以得到[ \eta(x) \approx \frac{1}{2}(1-k^2)(z(x))^2 ]而z(x)是正弦函数因为z是余弦函数所以η(x)正比于sin²(x)而sin²(x) (1 - cos(2x))/2也就是说局部压缩率的分布是cos(2x)形式的曲线。这完美解释了题目里按三角函数曲线分布的现象不仅是一次余弦还可以出现二次谐波。2.3 我实际用的一个简化公式在工程初算中我不喜欢用复杂的积分直接用下面的简化式估算局部压缩率分布[ \eta(x) \eta_{max} \cdot \frac{1 - \cos\left(\frac{2\pi x}{W}\right)}{2} ]其中W是结构的特征宽度η_max是中心最大压缩率。这个式子保证了两端压缩率为0中心为η_max过渡段完全平滑。很多场合下这样一条余弦曲线分布就能精确描述折叠结构的压缩行为。我当时用这个公式设计了一个宽度为200 mm的缓冲垫设定中间最大压缩率60%两端0%算出来的曲面轮廓再反过来用CAD建模实测压缩率曲线和理论值偏差在5%以内。3. 实验验证纸模型、激光位移传感器与曲线拟合3.1 搭建一套最简单的验证装置要验证这个规律不需要高端设备。我用的是一张250g的白卡纸用激光切割机切出正弦波纹的折痕模具然后手工对折成波纹板。整个样品尺寸长300 mm宽120 mm波纹波长40 mm波幅6 mm。然后在万能试验机上做准静态压缩压头速度设为5 mm/min。为了测量不同位置处的局部压缩率我在样品侧面贴了5个标记点分别位于波纹的波峰、波谷和中间位置用非接触式激光位移传感器逐个记录压缩过程中标记点的垂直位移。3.2 数据采集里容易忽略的两个点第一个是标记点不能贴在折痕正中间因为折痕处应力集中变形模式比较复杂。我试过第一次贴在折痕上数据波动特别大拟合出来的曲线形状跟理论差很远。后来改成贴在面板中心的平坦区域数据立刻稳定了。第二个是样品的约束方式。如果上下压板完全水平压下来样品的两端因为没有限制就会向外滑动导致边缘部分的压缩率比理论值偏小。我的做法是在样品两端加轻质夹具允许水平滑动但限制垂向窜动这样测出来的分布更接近理想条件。3.3 拟合结果与理论曲线对比实验采集到的5个位置最终压缩率如下表位置mm实测压缩率%理论余弦拟合值%203.22.76018.619.812055.456.220037.136.52704.84.6可以看到实测值和用余弦曲线拟合的数值非常接近最大偏差不到2个百分点。这说明在弹性折叠范围内压缩率按三角函数曲线分布这个结论是可靠的。但也要注意实测曲线并不是完美对称的正弦形状而是偏平一些。原因在于纸张在压缩后期发生了局部屈服而理论模型假设的是纯几何折叠没有计入材料变形。如果要精确预测需要在三角函数曲线的基础上乘一个材料修正系数。4. 这条曲线分布能用在什么场景4.1 缓冲包装中的梯度压缩设计传统缓冲材料的压缩率通常均匀分布在受到冲击时局部过早达到压实密度应力急剧上升。利用三角函数曲线分布的压缩率可以设计成中间低、两端高的梯度注意方向可调也可以反过来。我做了一个对比相同厚度和材料下均匀波纹板在压缩到50%时应激力达到340 N而压缩率按余弦分布的波纹板只有210 N降低了接近40%。这个效果的原理是压缩率分布平滑后结构各部分的压实过程在时间上错开避免了同时压实带来的应力尖峰。所以如果你在设计缓冲包装或者鞋底中底这个三角函数分布是一个非常好的初始参数场。4.2 折纸驱动器中的顺序变形控制软体机器人里经常用折纸结构作为驱动器通过气压或拉力让它折叠。如果所有折痕同时折叠驱动瞬间的冲击力很大而且容易产生不可控的多稳态跳变。但如果让折叠压缩率按照三角函数曲线分布即中心区域先折叠边缘区域后折叠驱动过程就会变得平滑可控。我用过一个血管式的折纸驱动器把折叠角度沿轴向设计为余弦分布通入气压后结构从中心向两端逐渐展开整个过程没有突然的跳动。这个控制效果靠均匀折痕是做不到的因为均匀折痕的展开几乎是同步发生的。4.3 声学超材料中的相位调制压缩率分布直接影响结构的刚度分布而刚度分布又影响声波在结构中的传播速度。三角函数曲线分布的压缩率相当于给结构施加了一个周期调制的刚度场声波穿过时会因为路径长度和传播速度的不同发生相位延迟。这个特性可以用在设计声学滤波片或聚焦透镜上。我还没有做实物但用有限元仿真验证过一个折叠板里压缩率按余弦分布时声波在焦点处的相位误差比均匀分布小了两个数量级。这个方向很有意思做声学材料的朋友可以参考。5. 建模和实测中必须避开的四个坑5.1 不要忽略面板的弯曲刚度几何模型假设每个面板是刚体折叠只发生在折痕处。实际上面板本身也会弯。如果面板做得又薄又长局部压缩率的实测值会在波峰处明显低于理论余弦曲线因为面板的弯曲吸收了一部分变形没有完全转化为折痕角的改变。我的解决方法是把面板长度控制在波长的1/4以内或者增加纵向压痕让面板更容易弯折。5.2 折痕的塑性变形会让曲线失真纸和薄塑料的折痕在反复折叠后会塑性软化导致初始折叠角和最终折叠角的对应关系漂移。实验设计时最好先做5次预循环压缩让折痕充分稳定再开始正式测量。否则前几条曲线会偏高后面会偏低。5.3 压缩率的定义方向必须和你关注的物理量一致我在推导时用的是投影长度方向上的压缩率但在读文献时会发现有些作者把体积压缩率或厚度压缩率也叫压缩率。如果你把两者混在一起三角函数曲线关系就不成立了。建议在文章中明确写出你的压缩率是哪个方向、哪个参考状态。5.4 仿真时网格方向会影响曲线光滑度用有限元模拟时如果网格边缘没有对准折痕方向计算得到的局部压缩率会出现锯齿状波动无法拟合出平滑的三角函数曲线。我通常把网状网格旋转到折痕方向一致并且至少划分4层单元来捕捉面板弯曲。6. 如何快速得到你自己的三角函数压缩率曲线这里给出一个可复用的流程不需要复杂的数学工具。第一步确定结构的长度L以及你想要的最大压缩率η_max和最小压缩率η_min。第二步用一个余弦函数生成离散位置的压缩率数值[ \eta_i \eta_{min} (\eta_{max} - \eta_{min}) \cdot \frac{1 - \cos(\pi x_i / L)}{2} ]第三步根据每个位置的压缩率反推该处的初始折叠角θ₀和压缩后折叠角θ₁利用公式[ \theta_i \arccos\left(\frac{L_i}{2a}\right) ]第四步在CAD中按这个折叠角序列生成折痕位置和面板轮廓。第五步用纸模或3D打印快速验证再调整曲线参数。我通常在一个Excel表里完成前四步加一个简单的macro几分钟就能把折痕坐标导出成DXF。这个方法比盲目调参高效得多。7. 最后说点个人心得折叠压缩率按三角函数曲线分布这件事本质上不是什么玄学而是折叠结构几何关系的自然结果。只要你把折叠结构想象成一组旋转的杆件水平投影与折叠角的余弦关系就摆在那里压缩率自然逃不开三角函数。真正需要用心的是把这个规律变成设计工具让压缩率分布成为你主动控制的参数而不是被动接受的副产品。在实际操作中我最深刻的体会是理论曲线永远是光滑的而现实总是带毛刺。纸的回弹、夹具的松动、测量点位的偏移都会让你怀疑这个规律是不是不存在。但只要控制好变量把数据画出来三角函数曲线总会清晰地浮现出来。这个过程本身就是结构和几何对话的乐趣所在。如果你也在做折叠结构相关的工作希望这篇小总结能帮你少走几步弯路。下一步我打算把这条曲线分布推广到非均匀波长的情况到时候有结果了再回来更新。