用Unity和C#手把手教你实现一个简单的社会力模型（Social Force Model）模拟器

用Unity和C#实现社会力模型的人群模拟系统

在游戏开发和虚拟仿真领域，模拟真实人群行为一直是个令人着迷的挑战。想象一下，当你在开发一个大型开放世界游戏时，如何让NPC在集市中自然地流动？或者在建筑设计中，如何预测紧急情况下的人群疏散路径？这正是社会力模型(Social Force Model)大显身手的地方。

社会力模型由物理学家Dirk Helbing在2000年提出，它将行人间的互动简化为物理力的作用。不同于传统路径规划算法，这个模型能自发产生复杂的群体行为模式——从有序流动到恐慌逃散。本文将带你从零开始，在Unity中实现一个完整的社会力模拟系统，包括核心算法、可视化调试和参数调优技巧。

1. 项目准备与环境搭建

首先创建一个新的3D Unity项目（2021 LTS或更新版本）。我们将使用C#脚本来实现所有逻辑，不需要额外插件。在Hierarchy中创建以下基本结构：

SimulationManager (空对象) ├── Agents (空对象，用于存放所有行人实例) └── Obstacles (空对象，用于存放障碍物)

安装必要的包：

• Unity.Mathematics：高性能数学运算
• Burst：提升计算效率
• （可选）UnityEditor.UI：用于调试界面

创建三个核心C#脚本：

1. SocialForceAgent.cs- 行人个体逻辑
2. SocialForceSimulation.cs- 全局管理
3. Obstacle.cs- 障碍物定义

提示：启用Burst编译可以显著提升大规模人群模拟的性能，但调试时建议先关闭

2. 社会力模型的核心实现

2.1 定义Agent基础属性

打开SocialForceAgent.cs，我们先定义行人的基本物理属性：

[System.Serializable] public class AgentParams { public float mass = 80f; // 行人质量(kg) public float radius = 0.3f; // 碰撞半径(m) public float desiredSpeed = 1.2f; // 期望速度(m/s) public float relaxationTime = 0.5f; // 速度调整时间常数 } public class SocialForceAgent : MonoBehaviour { [Header("Movement Parameters")] public AgentParams parameters; public Vector2 currentVelocity; public Vector2 desiredDirection; private Vector2 _position; private Rigidbody2D _rb; }

2.2 实现三种核心作用力

社会力模型的核心是三种力的叠加计算：

1. 自驱力：使行人保持期望速度和方向

Vector2 CalculateSelfDrivingForce() { Vector2 desiredVelocity = desiredDirection * parameters.desiredSpeed; return (desiredVelocity - currentVelocity) / parameters.relaxationTime; }

2. 行人排斥力：避免与其他行人碰撞

Vector2 CalculateAgentRepulsion(SocialForceAgent other) { Vector2 direction = _position - other._position; float distance = direction.magnitude; float minDistance = parameters.radius + other.parameters.radius; if (distance > minDistance) return Vector2.zero; float magnitude = Mathf.Exp((minDistance - distance) / 0.1f); return direction.normalized * magnitude; }

3. 障碍物排斥力：避开墙壁和障碍物

Vector2 CalculateObstacleForce(Obstacle obstacle) { Vector2 closestPoint = obstacle.GetClosestPoint(_position); Vector2 direction = _position - closestPoint; float distance = direction.magnitude - parameters.radius; if (distance > 1f) return Vector2.zero; float magnitude = Mathf.Exp(-distance / 0.2f); return direction.normalized * magnitude; }

2.3 整合力计算与运动更新

在FixedUpdate中整合所有力的作用：

void FixedUpdate() { Vector2 totalForce = Vector2.zero; // 自驱力 totalForce += CalculateSelfDrivingForce() * parameters.mass; // 行人排斥力 foreach (var other in SimulationManager.Instance.agents) { if (other != this) totalForce += CalculateAgentRepulsion(other); } // 障碍物排斥力 foreach (var obstacle in SimulationManager.Instance.obstacles) { totalForce += CalculateObstacleForce(obstacle); } // 更新速度 currentVelocity += totalForce * Time.fixedDeltaTime; // 限制最大速度 if (currentVelocity.magnitude > parameters.desiredSpeed * 1.5f) { currentVelocity = currentVelocity.normalized * parameters.desiredSpeed * 1.5f; } // 更新位置 _rb.velocity = currentVelocity; }

3. 可视化与调试系统

3.1 实时可视化力场

为了直观调试，我们可以添加力场可视化：

void OnDrawGizmosSelected() { // 绘制期望方向 Gizmos.color = Color.green; Gizmos.DrawRay(transform.position, desiredDirection * 0.5f); // 绘制当前速度 Gizmos.color = Color.blue; Gizmos.DrawRay(transform.position, currentVelocity * 0.5f); // 绘制排斥力范围 Gizmos.color = new Color(1,0,0,0.1f); Gizmos.DrawSphere(transform.position, parameters.radius * 3f); }

3.2 创建调试UI

添加简单的UI来实时调整参数：

[Header("Debug Controls")] public Slider massSlider; public Slider speedSlider; void Start() { massSlider.onValueChanged.AddListener(v => parameters.mass = v); speedSlider.onValueChanged.AddListener(v => parameters.desiredSpeed = v); }

4. 高级优化技巧

4.1 空间分区优化

当行人数量超过100时，需要优化碰撞检测：

// 在SimulationManager中添加 private SpatialGrid _grid; void Update() { _grid.Clear(); foreach (var agent in agents) { _grid.Add(agent); } } // 在Agent中修改排斥力计算 foreach (var other in _grid.GetNeighbors(_position, 2f)) { if (other != this) totalForce += CalculateAgentRepulsion(other); }

4.2 Burst加速计算

使用Unity的Burst编译器提升性能：

[BurstCompile] public struct ForceCalculationJob : IJobParallelFor { [ReadOnly] public NativeArray<Vector2> positions; [ReadOnly] public NativeArray<float> radii; public NativeArray<Vector2> velocities; public void Execute(int index) { // 并行计算逻辑 } }

4.3 典型参数配置参考

下表展示了不同场景下的推荐参数：

5. 实战案例：地铁站人流模拟

让我们构建一个完整的地铁站场景：

1. 场景布置：

   • 创建站台（长20m×宽5m）
   • 添加4根立柱作为障碍物
   • 设置两个入口和两个出口

2. 行人设置：

   • 生成100个行人
   • 设置70%的人目标为出口A，30%为出口B
   • 随机分布初始速度（0.8-1.3m/s）

3. 特殊行为处理：

// 在Agent中添加群体行为 void UpdateDesiredDirection() { if (Random.value < 0.01f) { // 1%几率改变目标 desiredDirection = (newExitPosition - _position).normalized; } // 跟随前方行人 RaycastHit2D hit = Physics2D.Raycast(_position, desiredDirection, 3f); if (hit.collider != null && hit.collider.CompareTag("Agent")) { desiredDirection = Vector2.Lerp(desiredDirection, hit.collider.GetComponent<SocialForceAgent>().currentVelocity.normalized, 0.1f); } }

4. 性能优化统计：

在实现过程中，有几个关键发现值得注意：首先，排斥力的指数衰减系数对模拟真实性影响极大——过强会导致行人像弹球一样反弹，过弱则会发生不自然的穿透。经过多次测试，0.1-0.3的范围最适合大多数场景。其次，在转角处添加轻微的向心力可以显著改善行人转弯的流畅度。