Unity中实现船只与海浪物理交互：从Gerstner波到浮力模拟

1. 项目概述：从静态模型到动态交互的沉浸感飞跃

在Unity里做海面场景，很多朋友可能都经历过：费劲找了一个漂亮的水体Shader，调了半天参数，海面波光粼粼，看着挺像那么回事。然后你兴冲冲地把做好的船模拖进去，结果发现船像个铁疙瘩一样杵在海面上，纹丝不动，或者只是僵硬地上下平移，瞬间就出戏了。这种“水上漂”的割裂感，是很多海洋场景项目从“好看”迈向“可信”的关键瓶颈。

“随波起伏的船只”这个需求，核心要解决的就是物理可信度与视觉沉浸感的问题。它不仅仅是让船体简单地跟随一个正弦波上下运动，而是要模拟船体作为一个刚体，与复杂、动态的海面进行实时、双向的物理交互。船头应该劈开浪花，船身会根据波浪的走向产生横摇和纵摇，船尾会留下扰动的水痕。这背后涉及到海浪的数学模拟、刚体物理计算、以及两者之间高效的数据耦合。

我最近在为一个海洋探索类的独立游戏项目搭建这套系统，从最基础的波浪公式开始，到整合物理引擎，再到优化性能，踩了不少坑，也总结出一套相对稳定、效果不错的方案。这个教程将分享我的实现路径，并提供一套可复用的核心脚本与Shader资源思路。无论你是想做一个简单的航海demo，还是开发复杂的海战游戏，这套关于船只与海浪交互的底层逻辑都能为你提供一个坚实的起点。

2. 核心思路拆解：分离、采样与力的艺术

实现船只随波起伏，最忌讳的就是把船和海浪的代码逻辑搅在一起。一个清晰、解耦的架构是后期调试和效果迭代的基础。我的核心思路可以概括为三步：海浪高度场生成 -> 船体采样与姿态计算 -> 物理力施加。

2.1 海浪模拟：不止于视觉的高度场

首先，我们必须建立一个海浪高度场。这个“场”是一个数学函数，对于世界坐标系中的任意一个水平位置(x, z)，它都能返回该点当前时刻的海平面高度y。这是所有计算的基础。

常见的模拟方法有几种：

• Gerstner波：这是游戏和图形学中最经典的海浪模型。它的优势在于不仅能计算高度，还能通过波函数对水平坐标进行偏移，从而模拟出波峰尖锐、波谷宽阔的真实波形，并且能自然地计算出每个点的法线（用于光照）和切线（用于流动贴图）。
• FFT（快速傅里叶变换）海洋：通过叠加大量不同频率、方向的正弦波，并进行逆FFT来生成高度场。效果非常逼真，能模拟广阔的、统计学的海面，但计算量较大，通常用于高端PC或主机游戏。
• 贴图位移：使用一张预先计算好的、循环播放的高度图来驱动顶点位移。性能极好，但动态性和变化性较弱，适合移动平台或风格化水面。

对于大多数需要实时交互的中小型项目，叠加多个Gerstner波是一个在效果和性能上取得绝佳平衡的方案。你可以通过调整每个波的振幅、波长、速度和方向，来组合出从平静湖面到惊涛骇浪的不同海况。

关键点：我们模拟海浪的首要目的，不是为了渲染（那是Shader的工作），而是为了在CPU端或Compute Shader中生成一个可供查询的高度场数据。船只物理系统需要频繁地向这个“场”查询其所在位置的海面高度和法线。

2.2 船体交互：多采样点与浮力模拟

让一个长方体（船体）简单地采样其中心点下方的海浪高度并移动，效果会很假。真实的船只会因为形状和重量分布，与海面产生多个接触点。

我的做法是，在船体上预设多个采样点（Sample Points）。这些点通常沿着船体的龙骨线（从船头到船尾的中线）和船舷分布。在每一帧物理更新时：

1. 获取每个采样点在世界空间中的水平坐标(x, z)。
2. 向海浪高度场查询该坐标处的当前海面高度y_water和海面法线。
3. 比较采样点当前的世界高度y_point与y_water。如果y_point低于y_water，说明该点浸入水中，需要受到浮力；反之则没有。

浮力的计算遵循阿基米德原理的简化版：浮力大小与浸入体积成正比，方向垂直于海平面（即沿着海面法线方向）。我们可以用采样点浸入的深度(y_water - y_point)来近似估算该点贡献的浮力大小。最后，将所有采样点产生的浮力，汇总计算出一个总的作用力和扭矩（旋转力），施加到船体的Rigidbody上。

2.3 物理集成：Rigidbody与力的艺术

这里是实现“真实感”的关键，也是新手最容易出错的地方。绝对不能直接修改船体Transform的position和rotation！必须通过Unity的物理引擎（PhysX）来驱动。

我们需要为船只添加一个Rigidbody组件，并合理设置质量（Mass）、阻力（Drag/Angular Drag）和碰撞体（Collider）。在FixedUpdate中（这是物理更新的周期）：

1. 遍历所有采样点，计算当前帧的总浮力矢量totalForce和总扭矩totalTorque。
2. 使用Rigidbody.AddForceAtPosition(force, position)方法。这个方法允许你在世界空间的特定位置（采样点位置）施加一个力，物理引擎会自动计算出这个力对刚体产生的线性运动（移动）和角运动（旋转）影响。这是实现逼真横摇纵摇的核心API。
3. 适当增加Rigidbody的线性阻力和角阻力，可以模拟水对船体的阻尼效应，让运动不会过于“弹跳”，更显沉稳。

通过这种方式，船体的运动完全由物理引擎驱动，它会与Unity环境中的其他物理对象（比如其他船只、漂浮物）自然交互，碰撞检测也完全正常。

3. 核心模块实现详解

下面，我将分步骤拆解关键代码模块。我会提供核心代码片段并解释其原理，你可以根据自己项目的具体需求进行调整。

3.1 Gerstner波高度场的C#实现

我们首先创建一个WaveManager单例类，用于管理和计算全局的海浪高度场。它存储了多个Gerstner波的参数。

using System.Collections.Generic; using UnityEngine; [System.Serializable] public class GerstnerWave { public float amplitude; // 振幅 (A): 波峰到静水面的高度 public float wavelength; // 波长 (L): 两个波峰之间的距离 public float speed; // 速度 (S): 波前进的快慢 public Vector2 direction; // 方向 (D): 标准化后的(x,z)方向向量 public float steepness; // 陡峭度 (Q): 控制波峰的尖锐程度，通常 <= 1/(amplitude * waveCount) private float _frequency; private float _phaseConstant; public float Frequency => _frequency; public float PhaseConstant => _phaseConstant; public void UpdateDerivedParams() { // 根据物理公式推导：频率 w = 2π / L, 相位常数 φ = S * 2π / L // 深水波色散关系简化：ω² = gk，这里我们用速度来控制相位，更直观 _frequency = 2 * Mathf.PI / wavelength; _phaseConstant = speed * _frequency; } } public class WaveManager : MonoBehaviour { public static WaveManager Instance; public List<GerstnerWave> waves = new List<GerstnerWave>(); public float timeScale = 1.0f; // 全局时间缩放，用于控制海浪动画速度 private float _currentTime = 0f; void Awake() { if (Instance != null && Instance != this) { Destroy(this.gameObject); } else { Instance = this; } foreach (var wave in waves) { wave.UpdateDerivedParams(); } } void Update() { _currentTime += Time.deltaTime * timeScale; } // 核心方法：获取世界坐标(x,z)处，海浪的高度(y)和法线(normal) public void GetWaveInfo(Vector3 worldPos, out float height, out Vector3 normal) { height = 0f; Vector3 tangent = new Vector3(1, 0, 0); // 初始切线 Vector3 binormal = new Vector3(0, 0, 1); // 初始副切线 Vector2 posXZ = new Vector2(worldPos.x, worldPos.z); foreach (var wave in waves) { float freq = wave.Frequency; float phase = wave.PhaseConstant; Vector2 dir = wave.direction.normalized; // 计算点积和相位 float dot = Vector2.Dot(dir, posXZ); float theta = freq * dot + phase * _currentTime; // 计算正弦余弦值，避免重复计算 float cosTheta = Mathf.Cos(theta); float sinTheta = Mathf.Sin(theta); // Gerstner波公式核心部分 float QA = wave.steepness * wave.amplitude; height += wave.amplitude * sinTheta; // 水平偏移（Gerstner波特性） float horizontalOffset = QA * cosTheta; Vector2 offset = dir * horizontalOffset; // 对切线和副切线的贡献（用于计算法线） float WA = wave.amplitude * freq; float factor = WA * cosTheta; tangent += new Vector3(-dir.x * dir.x * factor, dir.x * factor, -dir.x * dir.y * factor); binormal += new Vector3(-dir.x * dir.y * factor, dir.y * factor, -dir.y * dir.y * factor); } // 通过切线和副切线的叉积计算法线 normal = Vector3.Cross(binormal, tangent).normalized; // 注意：这里计算出的法线是未归一化的，但叉积后再归一化结果正确 } }

这个类提供了GetWaveInfo方法，它是船只脚本查询海面信息的唯一接口。在Shader中，我们也可以用完全相同的公式来计算顶点偏移和法线，确保视觉和逻辑一致。

3.2 船只浮力计算器脚本

接下来是附着在船体上的核心脚本BoatBuoyancy。

using System.Collections.Generic; using UnityEngine; public class BoatBuoyancy : MonoBehaviour { public Rigidbody boatRigidbody; public float buoyancyStrength = 2.0f; // 浮力系数，用于微调 public float waterDrag = 1.0f; // 水中线性阻力 public float waterAngularDrag = 1.5f; // 水中旋转阻力 public List<Transform> samplePoints = new List<Transform>(); // 采样点空物体 // 调试用，在Scene视图中绘制采样点和力 public bool debugDraw = true; public float debugForceScale = 0.1f; private float _originalDrag; private float _originalAngularDrag; void Start() { if (boatRigidbody == null) boatRigidbody = GetComponent<Rigidbody>(); if (boatRigidbody == null) { Debug.LogError("BoatBuoyancy requires a Rigidbody on the same GameObject or assigned."); enabled = false; return; } _originalDrag = boatRigidbody.drag; _originalAngularDrag = boatRigidbody.angularDrag; // 如果没有手动指定采样点，尝试在子物体中寻找 if (samplePoints.Count == 0) { var points = GetComponentsInChildren<Transform>(); // 这里可以添加一个命名规则，例如所有名为“BuoyancyPoint”的子物体 // 简单起见，这里只是警告 Debug.LogWarning("No sample points assigned for buoyancy. Please assign transforms in the inspector."); } } void FixedUpdate() { if (WaveManager.Instance == null) return; bool isUnderwater = false; float totalSubmergedVolume = 0f; // 估算总浸入体积，可用于高级效果 foreach (var point in samplePoints) { Vector3 worldPos = point.position; WaveManager.Instance.GetWaveInfo(worldPos, out float waterHeight, out Vector3 waterNormal); float pointHeight = worldPos.y; float depth = waterHeight - pointHeight; // 如果采样点在水面以下 if (depth > 0) { isUnderwater = true; // 简化浮力计算：力的大小与浸入深度成正比，方向为水面法线方向 // 更精确的模型可以考虑采样点代表的体积 float buoyancyForceMagnitude = buoyancyStrength * depth * Physics.gravity.magnitude; Vector3 buoyancyForce = waterNormal * buoyancyForceMagnitude; // 在采样点位置施加浮力 boatRigidbody.AddForceAtPosition(buoyancyForce, worldPos); // 调试绘制 if (debugDraw) { Debug.DrawLine(worldPos, worldPos + buoyancyForce * debugForceScale, Color.cyan); } totalSubmergedVolume += depth; // 粗略估算 } } // 根据是否在水下，调整阻力模拟水的阻尼效果 if (isUnderwater) { boatRigidbody.drag = waterDrag; boatRigidbody.angularDrag = waterAngularDrag; } else { // 完全离开水面时恢复原阻力（如空气阻力） boatRigidbody.drag = _originalDrag; boatRigidbody.angularDrag = _originalAngularDrag; } } // 在编辑器Scene视图中可视化采样点 void OnDrawGizmosSelected() { if (!debugDraw) return; Gizmos.color = Color.green; foreach (var point in samplePoints) { if (point != null) Gizmos.DrawSphere(point.position, 0.1f); } } }

脚本使用要点：

1. 采样点布置：在船体模型下创建多个空物体（GameObject），作为采样点。通常沿龙骨线（船底中线）从船头到船尾布置3-5个点，在船体两侧再各布置1-2个点。点的数量和位置直接影响船的稳定性和摇晃行为，需要反复调试。
2. 参数调试：buoyancyStrength是关键。值太小，船会沉；值太大，船会像软木塞一样剧烈弹跳。waterDrag和waterAngularDrag用于模拟水的粘滞阻力，让船的运动有“重量感”和“阻尼感”，不会无限摇晃。
3. AddForceAtPosition：这是魔法发生的地方。在浸入水中的采样点位置施加一个向上的力，物理引擎会自动计算出这个力如何使船体旋转（产生横摇、纵摇）和移动。

3.3 与渲染同步的海浪Shader（概念）

为了让船体与海浪的视觉效果完美融合，我们需要一个与WaveManager使用相同参数的海面Shader。这里给出一个Unlit Shader的核心片段，展示如何在顶点着色器中应用Gerstner波。

// 这是一个简化的Shader结构，实际项目可能需要包含光照、镜面反射等 Shader "Custom/OceanGerstner" { Properties { _MainTex ("Texture", 2D) = "white" {} _Color ("Color", Color) = (0.1, 0.3, 0.6, 1) // 波浪参数可以通过脚本传递，这里也可以定义一些覆盖参数 _WaveAmplitude ("Wave Amplitude", Float) = 0.5 _WaveLength ("Wave Length", Float) = 10 // ... 其他属性 } SubShader { Tags { "RenderType"="Opaque" } LOD 100 Pass { CGPROGRAM #pragma vertex vert #pragma fragment frag #include "UnityCG.cginc" // 假设波浪参数通过全局数组从C#脚本传递进来 // 在实际项目中，你可能需要使用StructuredBuffer或全局数组 uniform float4 _WaveParams[10]; // [振幅， 波长， 速度， 陡峭度] uniform float2 _WaveDirections[10]; uniform float _GlobalTime; struct appdata { float4 vertex : POSITION; float2 uv : TEXCOORD0; }; struct v2f { float2 uv : TEXCOORD0; float4 vertex : SV_POSITION; float3 worldNormal : TEXCOORD1; }; // Gerstner波计算函数（与C#端保持一致） float3 GerstnerWave(float3 worldPos, float amplitude, float wavelength, float speed, float steepness, float2 direction) { float freq = 2 * UNITY_PI / wavelength; float phase = speed * freq; float dotProd = dot(direction, worldPos.xz); float theta = freq * dotProd + phase * _GlobalTime; float cosT = cos(theta); float sinT = sin(theta); float QA = steepness * amplitude; float height = amplitude * sinT; float horizontal = QA * cosT; float3 offset; offset.x = direction.x * horizontal; offset.y = height; offset.z = direction.y * horizontal; return offset; } v2f vert (appdata v) { v2f o; // 将模型顶点转换到世界空间 float3 worldPos = mul(unity_ObjectToWorld, v.vertex).xyz; float3 waveOffset = float3(0,0,0); float3 tangent = float3(1,0,0); float3 binormal = float3(0,0,1); // 叠加多个波（示例叠加两个） for(int i = 0; i < 2; i++) { float amp = _WaveParams[i].x; float len = _WaveParams[i].y; float spd = _WaveParams[i].z; float stp = _WaveParams[i].w; float2 dir = _WaveDirections[i]; // 调用Gerstner函数，这里简化了法线计算 waveOffset += GerstnerWave(worldPos, amp, len, spd, stp, dir); } // 应用波浪偏移 worldPos += waveOffset; o.vertex = mul(UNITY_MATRIX_VP, float4(worldPos, 1.0)); o.uv = v.uv; // 法线计算（简化版，实际需要更精确的切/副切计算） // 可以使用屏幕空间导数或传递计算好的法线 o.worldNormal = normalize(float3(0,1,0)); // placeholder return o; } fixed4 frag (v2f i) : SV_Target { // 简单颜色输出 fixed4 col = _Color; return col; } ENDCG } } }

重要提示：在实际项目中，为了保持CPU（物理）和GPU（渲染）计算的一致性，最佳实践是将波浪参数（如振幅、方向）定义在一个ScriptableObject或静态配置类中，然后同时提供给WaveManager（C#）和Shader（通过MaterialPropertyBlock或全局Shader变量）。这样可以确保船“感受”到的波浪和眼睛“看到”的波浪是完全同步的，避免出现船在波谷而视觉在波峰的错位现象。

4. 性能优化与高级技巧

当你的海面上有多艘船，或者波浪非常复杂时，性能可能成为瓶颈。以下是一些优化思路和进阶技巧。

4.1 性能优化策略

1. 采样点优化：不是所有船都需要几十个采样点。对于小型、简单的船只，4-6个关键点可能就足够了。可以通过脚本动态减少远处或不重要船只的采样点更新频率（如每2-3帧更新一次）。
2. 波浪查询批量化：WaveManager.GetWaveInfo在每帧会被每个采样点调用多次。如果波的数量很多，计算量不小。可以考虑：
   • 使用Compute Shader：将所有的采样点坐标打包成一个Buffer，在Compute Shader中并行计算所有点的高度和法线，结果读回CPU或直接用于GPU。这是处理大量交互对象（如数百个漂浮物）的最高效方法。
   • 预计算或简化：对于远处或背景中的海面，可以使用较低复杂度的波浪公式，甚至用一张动态的高度图纹理来采样，而不是实时计算Gerstner波。
3. LOD（细节层次）系统：为船只的浮力计算设置LOD。近距离高精度（多采样点，每帧更新），远距离低精度（少采样点，降低更新频率，甚至用简单的正弦波动画代替物理模拟）。
4. 对象池与缓存：对于频繁生成和销毁的漂浮物（如炮弹溅落的水花、掉落的货物），使用对象池管理其浮力脚本，避免频繁的GetComponent和内存分配。

4.2 增强真实感的进阶技巧

1. 船首尾浪与航迹：可以通过在船头、船尾下方设置额外的粒子发射器或网格生成器来模拟。根据船体与水面的相对速度，控制粒子的发射速率和大小。更高级的做法是使用RenderTexture绘制航迹贴图到水面上。
2. 动态阻力：目前的waterDrag是常数。更真实的模拟是，阻力应与船体浸入水中的横截面积和速度的平方成正比。你可以根据采样点估算的浸没横截面，动态调整Rigidbody.drag。
3. 推进力与转向：不要直接给Rigidbody施加扭矩来转向。更好的方法是在船尾两侧各设置一个“推力点”，通过独立控制这两个点向后的力的大小差来实现转向（类似坦克的差速转向），这样更符合真实船舶的螺旋桨或舵效原理。
4. 海浪对船的阻力：目前船只会被动地被海浪推着走。可以增加一个计算，当船体侧面（船舷采样点）与波浪运动方向垂直时，受到一个侧向的阻力，模拟海浪拍击船身的效果。

5. 常见问题与调试实录

在开发过程中，我遇到了不少典型问题，这里记录下排查思路和解决方法。

5.1 船只表现异常问题排查表

5.2 调试与可视化技巧

1. 善用Debug.Draw：如代码所示，在FixedUpdate中使用Debug.DrawLine绘制每个采样点上的浮力矢量（用青色线表示）。这能让你直观地看到哪些点在水下，受力方向如何，是调试浮力逻辑最强大的工具。
2. 场景视图调试：在Scene视图运行时，你可以选择船体，查看其Rigidbody组件的速度、角速度矢量，帮助你理解船体的运动状态。
3. 自定义Editor工具：可以为BoatBuoyancy脚本编写一个简单的Editor脚本，在Inspector中增加一个按钮，一键在船体周围按规则生成采样点空物体，提升工作流效率。
4. 参数动画曲线调试：将关键参数如buoyancyStrength声明为public，在游戏运行时直接在Inspector中调整并观察效果，Unity会记录下修改过程，非常方便。

这套系统搭建完毕后，你的船只就不再是场景中一个孤立的模型，而是真正融入动态海洋环境的一部分。你可以通过调整波浪参数来创造从风平浪静到暴风骤雨的各种海况，观察船只的不同反应。这种底层物理交互带来的真实感，是任何简单的动画都无法比拟的。