物理模拟：从一根毛线到一块数字布料

可视化演示

展示质点弹簧布料网格受风力变形、Verlet积分约束求解与软体四面体碰撞变形，切换模式观察物理模拟。

布料基石：质点与弹簧网络

织毛衣的第一步是起针。奶奶在针上排出一排针脚，每个针脚只需要记住两件事：我此刻在哪儿，谁在拉着我。布料模拟的第一步也是“起针”——把布面离散成一张网格，每个网格交点叫质点，相邻质点之间用一根弹簧拴住。质点负责记录位置、速度和“轻重”（质量），弹簧负责记住两个质点之间的初始距离，也就是“原长”。

这一步的灵感来自生活：整块布的行为，等于所有针脚互相拉扯的合谋。网格越密，布越细腻，但计算量也越大。最常用的拓扑是四边形网格：每个内部质点与上下左右四个邻居相连。不过四边形抗剪切能力差，一拉就斜——所以工程上会再加对角弹簧，像毛衣的加固线一样，防止布面“扯歪”。

质点弹簧模型（Mass-Spring）是整个布料模拟的基石，它的能量由每根弹簧的形变决定。弹性势能公式如下：

就像起针决定了毛衣的版型，质点和弹簧构成了布料的“骨”。针脚密，布细腻；针脚疏，布粗糙。接下来，我们要看看这些毛线到底是怎么用力的。

弹簧动力学：胡克定律与阻尼

毛线拉长后会回弹，但不会永远弹跳下去——因为空气和线自身的摩擦在“刹车”。弹簧的这一体两面，正是模拟的日常：回弹交给胡克定律，刹车交给阻尼。如果只有回弹没有阻尼，布会像装了永动机一样荡个没完。

胡克定律说：弹簧的恢复力与形变量成正比，方向与拉伸方向相反。k 是弹性系数，形变量越大，回拉力越猛。而阻尼项与速度成正比：质点动得越快，阻力越大，方向与运动方向相反。把两项叠在一起：

有了力，质点就会动。但“动”是连续的，计算机只能离散地推进——这正是下一节“时间推进器”要解决的事。

数值积分：显式欧拉与 Verlet

织毛衣不能一步到位，得一行一行往上织。物理模拟也一样：把时间切成一小段一小段，每段叫一个时间步 Δt。每走一步，就根据当前状态推算下一帧。这个过程在数学上叫数值积分。

最直观的推进方法是显式欧拉：用当前速度直接推位置，用当前力直接推速度。但它有个著名的毛病——步长稍大就“爆炸”：位置越过平衡点太多，下一帧回弹更猛，最后整块布崩成雪花。布料模拟界更爱用 Verlet 积分：干脆不显式存速度，而是用前两帧的位置“外推”当前帧：

一针一针往前推，时间步就是“每一针的宽度”：针脚太宽，布就变形；太细，织得慢。但光有推进还不够——布不能穿过身体，这就轮到“版型规则”登场了。

碰撞与约束：投影的艺术

毛衣穿在身上，不能穿到身体里面去。道理简单，可数字世界里的“身体”和“布”都没有实体，谁都不拦着谁。于是物理引擎想出一个粗暴又聪明的办法：约束——规定某些关系必须满足，不满足就直接“掰”回去。

最典型的约束是距离约束：一根弹簧两端质点的间距要尽量保持原长。我们不调整力，而是直接计算误差，然后沿两点连线方向把位置修正到合法值。这种“直接改位置”的框架叫 PBD（Position Based Dynamics），是游戏引擎里最流行的布料方案。它的核心公式是约束修正量：

这就像织毛衣时随时对照版型修一修：哪里鼓了按回去，哪里紧了松一松——不纠结于力，直接改位置。碰撞检测也同理：质点穿透地面，就把它投影回地面上方。有了这套规则，布才真正“穿得住”。

软体物理：有限元与压力体

如果往毛衣里塞一团棉花，按压会凹、松手会弹——这就是软体。布是“面”，软体是“体”。软体不仅要像布那样抗拉伸，还要对抗一个更重要的东西：体积变化。捏一个果冻，它内部会被挤压，压强升高，然后把你的手指顶回去。

模拟软体有两条主流路线。有限元法（FEM） 把物体切成许多小四面体，逐个分析每个单元的形变与应力，精确但计算量惊人；压力体模型 则把物体内部当成一个气球腔，用气压抵抗挤压。游戏里常用后者：软体表面每个小面元都有面积 A，内部压强 P 沿外法线方向推出去：

就像棉花填充让毛衣立体，软体物理让虚拟物体有了“体积感”。一压一弹，都是压强在说话。从果冻怪物到软体机器人，这套模型撑起了半个“有弹性”的世界。

性能优化：空间哈希与 GPU 并行

一群织工同时织一件毛衣，得先分好片区：甲织左袖，乙织右袖，最后拼接。可要是所有人的手都能伸到同一个格子里，就打架了。布料模拟里最大的开销是碰撞检测：布料的每一点都要确认自己没穿进别的东西。朴素做法是每对质点都测一遍，复杂度 O(n²)，一万个质点就是五千万次运算，直接卡死。

空间哈希 把空间切成均匀的正方格子，每个质点只跟自己所在格子和周围 9 个格子里的邻居比对。查邻居的平均成本从 O(n) 降到 O(1)，总复杂度变成 O(n)。格子的坐标经过一个哈希函数，映射到一维数组：

除了空间划分，GPU 并行 是另一大杀器：质点和弹簧彼此独立，天然适合并行——每个线程管几个质点，一帧里跑好几轮迭代。再加上细节层次：远处用粗网格、近处用细网格。十万质点的披风，就这样在游戏里 60 帧飘了起来。

引擎实战：稳定条件与参数调优

奶奶会按季节选毛线：夏天细棉线透气，冬天粗毛线保暖。引擎里的布料也有“线号”——一堆参数等着你调。Unity 的 Cloth、UE 的 Chaos Cloth 都暴露了刚度（stiffness）、阻尼（damping）、迭代次数、质量等旋钮。调之前，必须先懂一条“生死线”：显式积分的稳定条件。

如果时间步长太大，弹簧力会把质点推过头，回弹再推更远，最终爆炸。临界步长由质量和刚度决定：

有了这条公式，调参顺序就清晰了：先定刚度，满足手感；再算临界步长，不够就加密子步；最后微调阻尼，去掉余振。工程上还有个常用技巧：多次小步迭代——把一步拆成几个子步，既绕过稳定性限制，又保持效率。毛线选不好毛衣就织不好，参数配不好布料就变成爆炸粒子特效。

常见伪影：抖动、拉伸与锁定

织得太松，毛衣领口越穿越大；织得太紧，缩水硬邦邦。布料模拟也有两个对应毛病：拉伸（超弹）和抖动（jitter）。拉伸来自约束太软或迭代不足：弹簧没被“掰”到位，布越拉越长；抖动来自步长太大或误差来回震荡：画面像筛糠。还有一个更隐蔽的问题——锁定：弯曲约束太强，布料硬得像纸板，怎么都不肯自然打褶。

约束投影的迭代收敛速度能解释很多现象：

知道毛病出在“松紧”，就能对症下药：加迭代、减步长、微调阻尼。若还嫌慢，用过松弛（over-relaxation）加速收敛，或者换 XPBD 消除过冲。一块平静妥帖的数字布料，就织出来了。

AI 应用场景

Python 演示

import mathimport random  class Particle:    def __init__(self, x, y, pinned=False):        self.x = x        self.y = y        self.ox = x        self.oy = y        self.pinned = pinned     def step(self, dt, gravity):        if self.pinned:            return        vx = (self.x - self.ox) * 0.98        vy = (self.y - self.oy) * 0.98        self.ox, self.oy = self.x, self.y        self.x += vx        self.y += vy + gravity * dt * dt  class Spring:    def __init__(self, a, b, rest):        self.a = a        self.b = b        self.rest = rest     def solve(self, particles):        p1 = particles[self.a]        p2 = particles[self.b]        dx = p2.x - p1.x        dy = p2.y - p1.y        dist = math.hypot(dx, dy)        if dist < 1e-9:            return        diff = (dist - self.rest) / dist * 0.5        if not p1.pinned:            p1.x += dx * diff            p1.y += dy * diff        if not p2.pinned:            p2.x -= dx * diff            p2.y -= dy * diff  class Cloth:    def __init__(self, cols=12, rows=9, gap=12.0):        self.particles = []        self.springs = []        for r in range(rows):            for c in range(cols):                x = c * gap                y = r * gap                jitter = random.uniform(-0.3, 0.3)                is_pinned = (r == 0)                self.particles.append(Particle(x + jitter, y, pinned=is_pinned))        for r in range(rows):            for c in range(cols):                i = r * cols + c                if c < cols - 1:                    self._add_spring(i, i + 1)                if r < rows - 1:                    self._add_spring(i, i + cols)     def _add_spring(self, i, j):        a = self.particles[i]        b = self.particles[j]        rest = math.hypot(b.x - a.x, b.y - a.y)        self.springs.append(Spring(i, j, rest))     def simulate(self, dt=0.016, gravity=500.0, iterations=4):        for p in self.particles:            p.step(dt, gravity)        for _ in range(iterations):            for s in self.springs:                s.solve(self.particles)        for p in self.particles:            if p.y > 300:                p.y = 300                p.oy = 300  def main():    cloth = Cloth(12, 9, 12.0)    for frame in range(60):        cloth.simulate()        if frame % 10 == 0:            xs = [p.x for p in cloth.particles]            ys = [p.y for p in cloth.particles]            width = max(xs) - min(xs)            print("frame=%3d width=%6.2f  bottom_y=%6.2f" % (frame, width, max(ys)))  if __name__ == "__main__":    main()

1. is_pinned = (r == 0)第一行质点固定，相当于把毛衣领口挂在衣架上，其余部分才能自由下垂。
2. vx = (self.x - self.ox) * 0.98用“当前位置减上一帧位置”还原速度，0.98 是隐式阻尼，避免永动。
3. self.y += vy + gravity * dt * dtVerlet 积分的位移更新，重力在这一小步里贡献一个二次位移。
4. diff = (dist - self.rest) / dist * 0.5距离误差归一化，两个端点各摊一半修正量。
5. for _ in range(iterations):迭代次数越多，弹簧越“听话”，布料越不容易拉伸。
6. if p.y > 300:地面碰撞的暴力投影：不许穿到地板下面，同时重置上一帧位置防止抖振。

常见误区

面试题