企业官网建设流程全解析

简介：UDEC中文指导说明（详尽版）是一份面向岩土工程初学者的离散元数值模拟教程，围绕通用离散元代码（UDEC）的基本原理与操作展开。内容涵盖软件安装启动、块体划分、材料模型定义、边界条件设置、初始平衡迭代等关键环节，并给出建模思路与模拟过程的系统讲解，适合需要快速上手UDEC、建立完整学习框架的读者。资源包内共1个doc文档，大小1.88MB，属于单文件电子文档，便于下载后直接阅读或打印。文档章节结构完整，包含引言、总论、与其他方法比较、一般特性、应用领域、概念与术语、命令语法及UDEC应用基础等模块，既能帮助初学者理解离散元理论，也能作为日常建模操作时的参考手册。目前已有2468人学习下载。通过这份资料，读者可以掌握从模型初始化到结果分析的全流程，为后续开展边坡稳定性、采矿巷道、隧道工程等复杂岩土数值模拟打下基础。

1. 拿到UDEC第一周，我建议你先别急着建模

做岩石力学和边坡工程的人，大概率绕不开UDEC这个软件。它全称是Universal Distinct Element Code，通用离散元程序，属于Itasca系列里主打“块体离散元”的二维分析工具。和FLAC、3DEC师出同门，但在处理节理岩体、块体滑动、结构面控制型破坏这些问题上，UDEC有它不可替代的位置。

为什么这么说？因为连续介质软件（比如FLAC）把岩体当成一个整体，节理、裂隙、层面这些弱面只能靠“弱化参数”去间接模拟，本质上是“等效”，不是“真实”。而UDEC的思路完全不同：它把岩体切成一个个独立的块体，块体之间通过接触来传递力和位移，节理就是实实在在的接触面。这样当边坡失稳时，块体可以真正地滑出去、翻倒、脱离，这是连续介质方法很难做到的。

我第一次用UDEC是在做某水电边坡的节理岩体稳定性分析，那会儿UEDC教程少得可怜，英文手册几百页，啃得头昏脑涨。后来陆陆续续用了四五年，踩坑无数，才慢慢摸清楚这个软件的正确打开方式。这篇东西不是照着手册念，纯粹是我个人把这些年用UDEC的经验、踩过的坑、反复试出来的参数习惯整理出来，给刚开始接触UEDC的同行一个参考。

需要注意的一点：UDEC在网上能搜到的信息相对分散，中文的系统性资料少，很多同学卡在“界面怎么用”“命令怎么写”“为什么不收敛”这些地方。下面这些内容，我就按照一个数值模拟项目的正常流程顺序来写：从拿到一个实际问题，到建模、设参数、求解、后处理，再到遇到问题怎么排查。中间穿插大量实测经验和具体命令，尽量让这篇文章能当作一份详细的中文操作参考书来用。

2. 界面与项目组织：建模之前先把这些摸清楚

2.1 UDEC的两种操作方式

UDEC有两种操作方式：一种是图形界面（UGUI），通过菜单点选完成建模和计算；另一种是命令流模式，在控制台逐条输入或批量执行命令。年轻人第一次接触时往往习惯用界面，毕竟所见即所得。但说实话，我的建议是：界面操作用来熟悉功能、快速查看模型状态，真正的工程分析应该用命令流。

原因有三条。第一，UDEC里面的很多高级控制参数、特殊设置藏得很深，菜单里根本找不到，必须用命令才能访问。第二，工程模型不是一次就能跑通的，调参数、改边界、加支护都是反复迭代的过程，界面操作每次都要重复点选，效率太低。第三，命令流本身就是一份完整的数据档案，三个月后你翻回来看，还能清楚知道当时模型是怎么建的、参数怎么设的，界面操作的记录很难保留。

2.2 启动后第一件事：设置工作目录

这个细节太容易被忽略了。UDEC计算会生成一堆中间文件，比如保存文件、历史记录文件、日志文件。如果你的工作目录一团乱，不同模型的文件混在一起，用不了多久就会搞混，甚至出现打开保存文件后模型内容和预期完全不同的情况。

我的习惯是：每一个分析项目建一个独立的文件夹，然后在文件夹内按“原始数据”“模型文件”“计算结果”“报告”分设子目录。每次启动UDEC时，通过控制台的cd命令切换到这个项目目录，或者直接在该目录下打开UDEC。千万别用默认目录跑完全部项目，后面整理数据时你会哭的。

2.3 单位制的混乱问题

UDEC默认采用国际单位制（SI）：长度用米，质量用千克，时间用秒，应力用帕斯卡。这套单位系统本身不复杂，但很多实际地质资料的单位并不是SI的，比如抗拉强度用兆帕、容重用千牛每立方米、变形模量用吉帕，这就容易埋下单位混乱的隐患。

我见过一个比较典型的错误：把室内试验得到的凝聚力（单位为千帕）直接输进去，但弹性模量用的是吉帕，结果算出来的破坏模式完全不对，还找不出问题在哪。所以建模前一定要先统一单位，我自己的做法是先在纸上把下面这些参数按SI单位列清楚：

• 密度：千克/立方米
• 重力加速度：9.8米/平方秒
• 弹性模量：帕斯卡
• 抗拉/抗压强度：帕斯卡
• 凝聚力：帕斯卡
• 内摩擦角：度
• 节理法向刚度/切向刚度：帕斯卡/米

顺手说一句，KPa换算成Pa是乘1000，GPa换算成Pa是乘10的9次方，这个不用记，但要保留换算过程在计算书里，方便核查。单位统一之后再开始建模，可以省掉后面大量排查的时间和精力。

3. 几何建模与块体切割的完整思路

3.1 做模型前先回答三个问题

很多初学者打开UDEC就急着画边界、切块体，结果模型跑起来要么不收敛，要么结果没法解释。我的经验是，动手操作之前先回答三个问题：这模型是二维的哪个剖面？模拟范围多大才合理？控制性结构面有哪些、产状如何？

第一个问题决定你在哪个平面上建模。UDEC本质是二维程序，它模拟的是某个剖面的力学行为，而不是三维空间的完整响应。比如边坡分析通常取主剖面，隧洞围岩稳定通常取垂直于洞轴线的断面。剖面选得不好，后面的结果会系统性地偏离现场实际。

第二个问题涉及边界范围。模型太小，边界效应会严重影响开挖或加载区的响应；模型太大，计算时间不划算。一般经验是：重点关注区域的边界到模型边界的距离，至少要有3到5倍的关键尺寸。比如隧洞直径10米，模型横向至少留30到50米；边坡高度50米，坡脚到模型底部也至少往下延伸150米以上。

第三个问题最容易被忽略，但也最致命。UDEC的核心价值就在于模拟结构面的控制作用，如果节理的产状、间距、延展性搞错了，后面算得再精细也是白算。建议把现场测绘得到的节理分组数据、玫瑰花图、极点投影图在工作前整理好，模型里的每一组节理都要有对应的地质依据。

3.2 几何建模流程我习惯用什么命令

UDEC建模常用的几何命令主要有：边界生成（polygon block或block edge）、块体切割（crack或jset）、删除/合并块体等。下面是一条最基础但很完整的建模命令序列：

; 设置模型边界，生成一个矩形区域 block polygon -50,-60 50,-60 50,60 -50,60 ; 生成随机节理组，方向60度，间距5米，倾角85度 ; 再生成一组，方向-30度，间距8米，倾角95度 jset angle 60 spacing 5 dip 85 jset angle -30 spacing 8 dip 95 ; 对块体进行切割 crack ; 查看生成结果 plot block

这里要注意，jset命令里面的angle参数是节理组的方向角，dip是节理在剖面内的倾角修正参数，很多新手容易搞混。默认情况下jset生成的节理是一组平行但位置随机的裂隙，spacing控制间距，angle控制走向。如果要模拟真实的岩体裂隙网络，还需要配合generate命令里的统计分布选项，比如节理间距服从负指数分布、迹长服从对数正态分布等。这些在jset里可以用stat选项来设置，想做得更接近真实的可以深入研究。

；上面这个例子是最简单的均质节理切割，真实工程中一般要分区域设节理组，比如洞室围岩的松动区、原岩区用不同参数，做法是先按区域切块，再单独加载节理组。

3.3 块体生成的坑

crack命令执行后，UDEC会根据已有的几何边界把模型切成一个个独立的块体。表面上看起来很简单，但实际操作中经常遇到几个问题。

最常见的问题是crack之后出现大量重叠或非常扁平的块体。原因往往是节理间距和角度设置不合理，导致节理之间夹角太小，切出来的块体长宽比严重失衡。这种畸形块体在后续计算中特别容易造成接触探测错误、收敛困难。

遇到这种情况，我的处理办法是：先执行crack remove命令把节理移除，调整节理参数后再重新切割。如果只是局部出现了个别畸形块体，可以手动用block delete把不合理的块体删掉，或者用block join把相近的块体合并。这个“建模-检查-修正”的循环我基本每次都要走好几遍，是正常的，别指望一步到位。

另外还有一个细节：crack之后的裂缝是“有长度无宽度”的，也就是说初始状态下节理是不开度的。如果需要模拟有初始开度的节理（比如张开裂隙），需要对特定节理赋予初始位移或开度值，命令层面可以用change配合相关节理参数来实现，但这个只在特殊问题中用到，常规分析可以不管。

4. 材料参数与边界条件：怎么设才靠谱

4.1 块体本构模型怎么选

UDEC中块体材料默认是刚性的，也就是块体本身不变形，所有变形都发生在节理上。这个假设在处理块体强度远高于节理强度、破坏以结构面控制为主的岩体时是合理的，比如硬岩边坡、块状岩体隧洞。

但如果你的问题是软岩、破碎岩体，块体的自身变形不可忽略，那就必须把块体设置为变形体（deformable），这时需要给块体赋予弹性模量、泊松比、密度等参数，并选择相应的本构模型。UDEC支持的本构模型包括各向同性弹性模型（elastic isotropic）、摩尔-库仑塑性模型（mohr）、应变软化模型（strain-softening）等，与FLAC3D用法类似。

我的选型经验是：优先用刚性块体做参数敏感性分析、机理研究阶段；进入工程方案比选和定量评价阶段，再把块体改成变形体。因为刚性块体计算速度明显更快，参数调试周期短。但要注意，刚性块体无法考虑块体自身的破碎，如果现场岩体极破碎，模型结果会偏“乐观”。

4.2 节理本构与参数，决定结果合理性

节理本构模型的选择和参数取值，是整个UDEC模拟中最关键、也是最难的一步。UDEC中最常用的是库仑滑动模型（Coulomb slip model），它把节理力学行为简化为法向刚度、切向刚度、凝聚力和摩擦角的组合。这个模型虽然简单，但工程上多数问题用它是够用的，关键是把参数取准。

这里先贴一段典型的节理参数赋值命令：

prop mat=1 jkn=5e9 jks=2e9 jfric=35 jcoh=0.1 jten=0.0

其中jkn是节理法向刚度（Pa/m），jks是切向刚度，jfric是摩擦角（度），jcoh是节理凝聚力（Pa），jten是节理抗拉强度（Pa）。

关于刚度参数，很多初学者会照搬手册上的经验值，但我的习惯是：刚度不宜取太大。节理刚度越大，计算时步越小，计算耗时越长，而且容易出现接触振动。一个朴素的判断准则是：节理法向刚度应当远大于周围岩体的等效刚度，但又不至于大到引发数值问题。如果你做的是层状岩体边坡，节理间距1米，岩体弹性模量5GPa，那jkn取5e9附近是合适的，因为它相当于把单位长度上的单条节理刚度与实际岩柱刚度匹配起来。如果jkn取到5e12，那结果也不会更准，计算时间反而暴涨，完全没必要。

接触黏聚力和抗拉强度对结果影响很大。实际节理大多是不能承受拉应力的，所以jten一般取0。但节理面有粗糙度，存在一定“表观黏聚力”，这可以从室内直剪试验的峰值强度包线反推。没有试验数据时，保守取jcoh=0也是可以的，算出来的结果偏安全但趋势合理。

4.3 边界条件，应力场比你想的重要

UDEC建模后默认边界是固定的（即给定位移为零）。但真实工程中，模型边界处于连续岩体中，其应力状态主要由地应力场控制。所以除非是边坡浅层分析、边界远离影响区等特殊情况，否则我强烈建议做地应力初始化，再开挖或加载。

地应力初始化的常用做法是：先关掉节理的本构行为（让它变成刚性），只对块体材料施加自重应力，并侧向施加一个kk系数控制的水平应力（即侧压力系数），然后求解到平衡。求解平衡后再把节理激活，开始正式的施工阶段模拟。

这里给一个简化的命令示例：

; 预设侧压力系数为1.2 insitu stress -8e6,0,-5e6 szz=-8e6 ; 设置边界条件：固定左右边界水平位移，底部竖直位移固定 boundary xvel=0 range x=-50, -50 boundary xvel=0 range x=50, 50 boundary yvel=0 range y=-60, -60 ; 施加重力 set gravity 0,-9.8 ; 弹性求解至初始平衡 solve elastic

命令中insitu stress的格式是stress xx, yy, xy，如果你是三维习惯的话需要特别注意它只有三个分量，别想着填zz，因为UDEC不关心面外应力。但程序允许填szz值，它对块体内部的应力计算有一定影响，但不会参与面内的位移求解。

很多人的模型结果怪异，比如开挖后围岩完全不破坏、或者一开始就大片屈服，原因往往是初始应力场没有正确还原。初始应力搞不准，后面一切分析都是空中楼阁。

5. 求解与收敛：那些折磨人的参数和消息

5.1 时间步长与阻尼：理解UDEC的“迭代时钟”

UDEC采用显式时步积分法，这意味着计算过程是逐步推进的。每一步的时间步长由所有单元和接触中最小的特征尺寸决定，块体越细、刚度越大，时间步长就越小，整体计算步数就越多。

显式法的好处是程序内存占用低、能处理大变形和非连续问题；代价是“迭代次数天生巨大”。一个中等规模的模型，几十万个时步是很常见的事。如果每个时步都很小，再叠加模型块体数量大，那计算一天一夜也不稀奇。

阻尼的作用是消耗系统动能，让求解过程最终收敛到静力平衡状态。UDEC提供了多种阻尼类型，其中局部阻尼（local damping）最常用，因为它的优势是对不同频率的振动都能有效衰减，且不影响最终静力解。solve命令会自动使用默认的局部阻尼系数（0.8），一般来说不需要动它。但如果系统持续振荡、不收敛，可以尝试增大阻尼值到0.9左右，或者改用瑞利阻尼配合粘性边界。

5.2 什么才算“收敛”

初学者最常见的疑问是：solve命令运行半天，控制台不停刷消息，到底什么时候才收敛？

UDEC默认的收敛判据是：最大不平衡力与平均不平衡力的比值降至某一阈值以下（一般默认1e-5），或者是最大节点速度低于某一值。你可以在命令里显式控制：

solve ratio 1e-5

意思是当最大不平衡力与代表性内力的比值小于1e-5时停止计算。

另一个实用的做法是设置最大时步数上限：

cycle 20000

防止程序无限跑下去。计算完成后，一定要查看不平衡力曲线和最大不平衡力值，如果曲线衰减得很好而且最终值很小，说明模型已经达到了较好的静力平衡。如果曲线发散或一直振荡，那就要返回检查模型。

此外我常用的一个判断技巧是查看关键点的位移历史记录。位移历史曲线如果在计算后期已经趋于水平直线，表明该点不再运动，模型基本稳定；如果还在持续线性增长，说明滑面已经形成、边坡可能正在发生渐进破坏，这时候模型不“收敛”恰恰是真实物理行为的反映——它本来就不该停住。

5.3 计算不收敛怎么办

如果模型怎么算都不收敛，教你三个排查方向。

第一，检查有没有“飞块”。UDEC迭代时如果某个块体速度异常增大，往往是因为它的周围约束不够或者接触参数太软。出现这种情况时，先用plot vel看速度分布，找出速度异常的块体，再检查该处的接触和边界条件。

第二，检查接触穿透。当接触刚度不足而外力过大时，块体可能互相“嵌入”，导致求解发散。可以增加接触刚度，或者检查块体之间是否存在不合理的重叠几何（上面提到的畸形块体）。

第三，检查阻尼设置。局部阻尼的值如果比默认值小很多，系统可能会长期振荡不收敛；比默认值大很多，又可能过度压制真实的不稳定。别轻易动阻尼，除非你清楚自己在干什么。

另外还有一点容易被忽略：建模时如果存在完全孤立的块体（俗称“悬浮块体”），它们在受力时容易发生绕定轴旋转、数值发疯，最好在建完模型后用plot block配合颜色区分检查一下孤立块体，必要时手动删除或给予约束。

6. 后处理与结果解读：图要能讲出工程故事

后处理阶段看似技术含量不高，但很多人的分析都栽在这上面：模型能算，图也能出，但结论说不清楚。以下三个操作是我每次必做的。

6.1 云图：位移与塑性区

用下面命令可以快速查看位移和塑性区分布：

plot block disp plot block plastic

plot block disp看的是位移场，能直观反映变形范围、滑动方向、影响区域。一个合格的位移云图要么有明显的位移集中带（对应潜在滑面位置），要么有清晰的变形分区（比如开挖卸载回弹区、应力集中区）。如果云图处处均匀、看不出层次变化，多半是地应力或者开挖卸荷设置有问题，回去查。

plot block plastic看的是塑性区分布。塑性区代表当前应力状态已经达到屈服条件的地方。在边坡分析中，如果塑性区从坡脚贯通到坡顶，意味着已经形成连续滑动面，安全系数就很低了。注意区分拉破坏、剪破坏和压破坏在图例中的颜色，不要张冠李戴。

6.2 安全系数怎么估计

UDEC本身没有内置强度折减法命令，但可以通过多次折减强度参数、观察计算收敛状态来获得边坡安全系数。

具体做法是：对节理参数jcoh和jfric同时除以折减系数F（比如先F=1.0跑一次，再F=1.2跑一次，逐步增大），每次折减后求解，观察边坡是否产生大范围位移或是否无法收敛。当折减到某个值，模型出现持续的、加速的位移，即认为达到了临界状态，这个F就近似是安全系数。

这种方法虽然原始，但在UDEC框架下非常实用，尤其适合节理控制型边坡，比连续方法的等效强度折减更逼近真实破坏机制。要注意的是每次折减后都要重新求解至平衡，并记录坡顶和坡脚几个监测点的位移增量，以位移突变作为判断临界值的主要依据。

6.3 历史记录是最后的“照妖镜”

我几乎在每次模拟中都会设置几个关键点的历史记录，命令如下：

hist id=1 ydisp 20,30 hist id=2 ydisp 20,35 plot hist 1 2

hist可以记录节点位移、速度、应力、不平衡力等随时间步的变化曲线。面板上多设几个测点后，你就能看到不同位置的响应时序关系和耦合关系，这对判断破坏演化过程很有价值。

举个小例子：某次我做含软弱夹层边坡的稳定性分析时，坡脚测点的位移曲线早早就出现了加速趋势，而坡顶测点还在缓慢移动。如果不看历史曲线，只凭最终云图很可能误判为“整体滑动”，但历史曲线清楚地显示这是一场从坡脚渐进发展至坡顶的破坏过程。两种解读对工程决策的影响完全不同。所以，历史记录不是可有可无的附加功能，它是理解计算过程、解读结果机制的关键工具。

7. 避坑与实战技巧汇总

7.1 建模精度与计算效率的平衡

UDEC的块体切割越细、节理越多，模型就越接近真实岩体结构，但计算时间也会急剧膨胀。如何平衡，我的经验是：在模型的关键区域（开挖影响区、潜在滑移带）加密节理，在远场区域适当增大节理间距，用稀疏网格包裹模型边界，既能保证核心区域的精度，又能有效控制总体计算量。

7.2 存档习惯：保存文件别只存一个

算了几小时突然程序崩溃、保存文件损坏，是每个数值模拟工作者最崩溃的时刻。UDEC支持save命令保存正在运行的模型状态，建议在每次求解阶段切换前都保存一次，并采用带序号的文件名，比如model_step1.sav、model_step2.sav。千万别一直覆盖同一个文件，否则返工都没法返。

7.3 常见报错的快速排查

围绕“不回收敛”“无法计算”“结果异常”三类问题，我整理了一个极简易的排查表：

7.4 从“能跑通”到“可信赖”，再走三步

最后说一句掏心窝的话：模型能算通、能出图，这只是第一步。要让计算结果真正可信、可用于工程决策，我建议每个项目至少再做三件事：

第一步，参数敏感性分析。将最不确定的参数（通常是节理刚度、摩擦角）上下浮动一定范围跑几组，看结果变化有多大。如果结果的敏感性过高，说明当前参数获取精度不够，需要加强试验或反演。

第二步，对比监测数据。如果项目有变形监测值，把模拟位移和实测位移放一起对比。差一个数量级基本可以断定参数或边界条件有问题，就要回头调整。数值模型的价值不在于复现，而在于预测；但预测的“底座”需要有实测数据来校验。

第三步，写清假设。每一步参数取值背后依据的试验报告、经验公式、工程类比都要在计算说明书中记录。数值分析报告如果通篇只有图表没有参数依据，是站不住脚的。这块工作不性感，但却是让模型从“学术玩具”变成“工程语言”的关键一步。

我用UDEC这些年，最大的体会是：这软件并不难学，难的是建立工程直觉，理解它算出来的东西到底在说什么。多跑模型、多对比实测、多怀疑自己的参数，慢慢就会形成一套属于自己的判断体系。上面的这些命令和经验，都是一次次算崩、熬夜调参之后攒出来的，希望能让你少走点弯路。

本文还有配套的精品资源，点击获取