程序员视角的计算力学：一套好懂好用的认知框架-尧图网络科技

程序员视角的计算力学：一套好懂好用的认知框架

- - - 做过多年工业仿真与三维图形软件开发，从极限平衡、关键块体、有限元前处理，到三维点云、网格处理、智能爆破软件CAD设计、Gis数字地球，踩过不少坑之后，慢慢攒出了这套跨工具的统一认知框架。

先唠3个底层共识，想通了后面都顺

所有力学现象，根子都是分子间的作用力和能量平衡，一层层放大，就成了我们看到的工程现象。但不用啥都从分子开始算，不同尺度有不同的玩法，够用就行。
世界本质上全是“散的”：岩石是矿物颗粒堆的，水是水分子凑的，空气是气体分子飘的。我们常说的“连续介质”，其实是个偷懒的宏观近似——当研究的东西，比组成它的粒子大得多的时候，忽略那些小缝隙，计算效率能翻好几个数量级，误差还在工程能接受的范围里。
站在写代码的角度，整个物理世界都能用「几何+能量+拓扑」搭出来：
几何搭出各种单元的样子，单元之间无非就是吸收能量、释放能量；能量一变，单元的XYZ位置就跟着动；大量单元一起位移，宏观上就成了变形、裂隙、断裂、滑移、旋转。
固液气三态全适用，最后落到计算机里，无非就是粒子、结构化网格、非结构化网格这几种表达形式。
注：这是仿真视角的简化，牺牲了极致微观的准确性，但工程上特别好使。

三个核心支柱：几何 · 能量 · 拓扑

几何是载体：不管是原子、颗粒、网格节点还是岩块，都能抽象成带XYZ坐标的几何单元，只是粒度大小不一样。
能量是唯一驱动力：单元之间的相互作用，本质就是吸能、放能。能量差变成力，推着单元位移、旋转，最后都往能量最低的状态跑，就稳定了。
补充：实际还有摩擦、阻尼这些能量耗散，变成热量跑掉了，工程上一般做简化处理。
拓扑是连接规则：单元之间粘不粘、连不连，全看拓扑。
- 拓扑没断 → 只是变形（弹性、塑性）
- 拓扑断了 → 就是破坏（断裂、破碎、分开了）

从微观到宏观，四个尺度各玩各的

不用强行从最底层往上推，干啥活就用对应尺度的方法。

分子/原子级（埃~纳米）
研究化学键、分子作用力，用分子动力学、第一性原理。纯搞机理的，工程开发基本用不上，知道有这么回事就行。
细观级（微米~毫米）
研究矿物颗粒、微裂纹、胶结物，用颗粒DEM、细观损伤力学。相当于微观和宏观之间的桥，解释宏观参数是怎么来的。
岩块级（厘米~米）
就是实验室里做实验的那块完整石头/铁块，看成均质的，用FEM、FDM算。主要是给工程计算提供基础参数。
工程岩体级（米~公里）
真正干活的尺度，有节理、断层、破碎带这些结构面，两种玩法：
- 连续等效法：把节理“抹平”，当成整体用FEM/FDM算，省事
- 块体离散元：顺着结构面切成一块块，算它们怎么滑、怎么撞

两大体系：连续 vs 非连续

先讲判断标准，遇到问题先套这个

看尺度：研究对象远大于组成单元 → 用连续；尺寸差不多 → 必须用非连续
看目标：想看整体应力、流速分布 → 用连续；想看破碎、滑移、碰撞过程 → 用非连续

补充：两者不是死的，过程中能互相转，比如完整岩石裂了，就从连续变成非连续了。

一、连续介质体系

就是把东西当成一整块无缝的，用“场”来描述，适合没大规模破碎分离的情况。

1. 固态分支（岩土/结构核心）

底层理论：固体力学三大方程（平衡、几何、本构），本构就是材料独有的应力应变规律，分弹性、弹塑性、流变、损伤这些。
主流方法：
- FEM有限元：最通用、精度最高，固体领域绝对主流
- FDM有限差分：显式求解，大变形、动力问题天然稳
- VEM虚拟单元：FEM升级版，前沿科研为主
开源工具：CalculiX、Code_Aster、OpenSees（岩土专属，本构特别全）
常用场景：算应力、边坡稳定、大坝受力、洞室围岩分析

我的实践：我自己做岩土三维可视化，核心用 VTK 做后处理渲染，配合 osgEarth 做大场景GIS地形展示；前后处理的几何建模、网格剖分，本质都是在给连续介质的FEM/FDM计算搭几何载体。

2. 液态分支（流体/渗流）

底层理论：普通流体用N-S方程，岩土渗流用达西定律简化。
主流方法：
- FVM有限体积法：CFD领域绝对主流，守恒性好
- LBM格子玻尔兹曼：介观粒子法，GPU跑特别快，复杂孔隙天生适配
  我的实践:我自己写了个小工具用这个模拟爆炸,爆破孔按时许进行连续爆破后,能量之间的互相影响.
- FEM有限元：算慢流、渗流还行，强对流场景不稳
开源工具：OpenFOAM（CFD标杆）、FluidX3D、OpenGeoSys（地质多场专用）
常用场景：坝体渗流、降雨入渗、地下水运移

3. 气态分支

可压缩流体，用FVM算，工具OpenFOAM、SU2。岩土里基本用不上，顶多算爆炸冲击波、巷道通风。

二、非连续介质体系

拆成一个个独立单元，各自按牛顿定律运动，靠接触力互相影响，专门对付破碎、滑移、翻滚的场景。

1. 颗粒类非连续

一个个小颗粒，用DEM离散单元法，靠接触刚度、摩擦、可断裂粘结这些规则运行。
开源工具：Yade（岩土DEM标杆）、LIGGGHTS（工业颗粒流）
常用场景：岩石破碎机理、泥石流、颗粒流

2. 块体类非连续（岩土最有特色的方向）

顺着结构面切成大岩块，用块体DEM或者DDA（块体自身还能变形）。
落地方案：
- 开源求解器：Yade（刚性块体模块）
- 商业标杆：3DEC（闭源，行业标杆，用来做结果对标）
我的实践：我自己的 XPlote 系列软件里，危岩体滑落模拟就是按这个思路做的：先用几何切割算法沿结构面拆块，再用 Bullet 物理引擎算动力学和碰撞，最后用 osgEarth 做大场景实时渲染。
实测下来，工程级场景用刚性块体假设完全够用，性价比比 DDA 高太多，不需要硬上可变形块体。
常用场景：岩质边坡失稳、危岩体滑落分裂、地下洞室掉块

复杂问题：跨体系耦合

遇到多物理场、全过程的问题，就混着来：

流固耦合（FSI）：流体用FVM + 固体用FEM/FDM，比如浪拍大坝、渗流驱动岩体变形
连续-离散耦合（FEM-DEM）：先算连续体损伤开裂，裂了之后转成离散块体算运动，比如完整岩体从裂到塌的全过程
流体-离散耦合（CFD-DEM）：流体带着颗粒/块体跑，比如水下滑坡、泥石流远程运动

遇到任何力学问题，五步定位法

定尺度：工程问题直接上宏观方法，搞机理再往下挖
判介质：选连续还是非连续体系
分物态：固、液、气对应不同的控制方程
选方法：按精度、效率、变形大小，挑数值方法
找工具：对应分支选开源求解器，或者自己搭引擎

放到代码里，本质是什么

分子尺度：几何对应分子的空间结构，拓扑对应化学键
颗粒尺度：一个个点/球就是粒子，没有固定连接，只有碰到了才相互作用，天生适合破碎、流动
变形与断裂的几何本质：
- 弹性/塑性变形：就是节点挪位置，网格的连接关系一点没变，所以FEM算起来特别快
- 断裂/破碎/分裂：连接关系断了，一个整体变成好几个独立对象，对应代码里就是三维网格拓扑重构，生成新的几何实体
我的实践：这也是我做三维切割、网格重建最深的感受：所谓的岩体分裂，代码里本质就是网格拓扑关系的断开，生成新的独立几何实体；很多仿真的难点，说到底不是力学难，是几何拓扑的处理难。