发布时间 : 星期二 文章UDEC3.1中文说明更新完毕开始阅读7818c6114431b90d6c85c706
1 引 言
1.1 总 论
通用离散元程序（UDEC，Universal Distinct Element Code）是一个处理不连续介质的二维离散元程序。UDEC用于模拟非连续介质（如岩体中的节理裂隙等）承受静载或动载作用下的响应。非连续介质是通过离散的块体集合体加以表示。不连续面处理为块体间的边界面，允许块体沿不连续面发生较大位移和转动。块体可以是刚体或变形体。变形块体被划分成有限个单元网格，且每一单元根据给定的“应力－应变”准则，表现为线性或非线性特性。不连续面发生法向和切向的相对运动也由线性或非线性“力－位移”的关系控制。在UDEC中，为完整块体和不连续面开发了几种材料特性模型，用来模拟不连续地质界面可能显现的典型特性。UDEC是基于“拉格朗日”算法很好地模拟块体系统的变形和大位移。
UDEC包含了功能强大的程序语言FISH函数。借助于FISH函数，用户可以编写自己的功能函数，扩展UDEC的应用功能。FISH函数为简化分析，适应特殊要求的UDEC的用户，提供了一个强有力的工具。
UDEC采用的离散单元法理论由Cundall（1971）首次提出，至今已经过了20多年的发展。在1985年，Cundall博士和Itasca公司在IBM系列兼容微机上开发了UDEC工程计算应用程序。该软件为建立数以千块模型的高速计算而设计。基于浮点运算速度的优势和低成本的内置RAM，用UDEC程序可大大地提高了计算大规模问题的能力。例如，在具有4MB RAM的微机上，UDEC能够求解2500个刚体（或1000个具有8个自由度变形体）的模型。该模型的求解速度大约为每分钟200次。在RAM确定的情况下，其计算速度是与模型的块体数量成线性关系。
对于典型的模型，约1500个刚体（或500个变形体）或更少，在UDEC中采用的显式解法，大约需要2000～4000计算步可以获得问题的解。例如，一个500个变形块体的模型，计算4000步大约需要6 min。因此，典型的工程问题用UDEC计算仅需几十分钟或几个小时。
UDEC是一个命令驱动（而不是菜单驱动）的计算程序。尽管菜单驱动程序易于初次学习，但在UDEC中所提供的命令驱动结构具有如下优点：
1、输入的“语言”是基于可识别的文字命令，使你易于识别每一个命令的作用（例如BOUNDARY命令，是指施加模型的边界条件）。
2、工程模拟通常是按照系列施工顺序构成 ―― 即，构造原岩应力，施加作用的荷载、开挖隧道、安装支护等。一系列（从文件或键盘上）输入命令完全对应于实际的施工顺序。
3、根据文本编辑器，很容易对UDEC数据文件进行编辑和修改。几个数据文件能
－1－
相互连接，进行多个问题求解，这对于进行参数的灵敏度分析是十分有用的。
4、命令驱动结构允许用户开发前后处理程序，控制UDEC必要的输入/ 输出。用户可以为一系列UDEC的模拟，编写节理模拟函数，产生特定的节理结构。可采用FISH程序语言，并插入到输入的文件中，使计算很容易实现。
1.2 与其他方法的比较
对于UDEC程序，一个共同的问题是，UDEC是一个有限元程序还是离散元程序？他们的主要区别是什么？UDEC程序与其他程序有何关系？为回答上述问题，下面将给予解释。
许多有限元、边界单元和拉格朗日有限差分程序都具有“界面单元”或“节理单元”，使程序能够模拟问题中的不连续面，扩大程序的应用范围。然而，他们的公式在一个或多个方面通常受到限制：首先，当考虑很多相互切割的节理就可能打乱系统的逻辑关系；其次，不可能自动识别新的接触面进行自动考虑；第三，计算公式可能有小位移和无转动条件限制，所以通常适用连续介质的程序。
术语“离散单元法”（Discrete element method）意味着： （a）允许离散块体发生有限的位移和转动，包括完全脱离； （b）在计算过程中，自动识别新的接触面。
在不连续介质中，如果没有第一个属性，程序不可能产生某些重要的机理。如果没有第二个特性，程序将限制在事先已知的相互作用的有限块体数。离散元法（Distinct element method）是由Cundall和Strack（1979）采用变形接触和显式、时间域的初始运动方程（而不是变换，块体方程）提出的特殊的离散单元法程序。
离散单元法的计算机程序主要有以下四类：
1、Distinct Element Programs － 该类程序采用显式时间步直接进行运动方程的求解。块体可以是刚体或变形体（通过细分成单元）；接触面是可变形的。UDEC就属此类。
2、Modal Methods － 该类方法类似于刚体离散单元法，但对于变形体采用模型叠加技术。
3、Discontinuous Deformation Analysis － 接触是刚体，块体可以是刚体或变形体。通过迭代算法可以获得非嵌入条件；块体变形性基于应变模型的叠加。
4、Momentum-Exchange Methods－接触面和块体都是刚体：块体接触面在瞬时碰撞的过程中惯性矩发生交换，可以表征滑动和摩擦特性。
1.3 一般特性
UDEC主要用于岩石边坡的渐进破坏研究及评价岩体的节理、裂隙、断层、层面对地下工程和岩石基础的影响。UDEC对研究不连续特征的潜在破坏模型是十分理想的工具。
－2－
当地质结构特征明显且易于明确描述的情况适宜使用该程序进行分析。UDEC开发了人工或自动节理生成器，用以模拟产生岩体中一组或多组不连续面。在模型中，可以产生变化范围较大的节理模式。屏幕绘图工具允许用户随时观看节理模型。在最后确定所选择的节理模型前，能容易进行调整与修改。
也可以获得不同的节理材料特性。基本模型是指定节理弹性刚度、摩擦角、粘结力、张拉强度和剪胀特性的库仑滑动准则。对该模型的改进包括随着位移的发展而粘结力和张拉强度的降低弱化。在此还可获得一个比较复杂的模拟连续屈服的节理模型，用以模拟弱化为累积塑性剪切位移函数的连续变化特性。作为一个选择模型，还可获得Barton－Bandis节理模型。节理模型和性质参数也可分别赋给单一节理或节理组。应当注意，即使地质图上所显示的节理为直线段，节理的几何粗糙度也可以通过节理材料模型加以表征。
UDEC的块体可以是刚体或变形体。对于变形块体，开发了包括用于开挖模拟的空模型（null）、应变硬化/软化的剪切屈服破坏模型以及非线性不可逆的剪切破坏和压缩模型。因此，块体能被用来模拟回填、土体介质以及完整岩石。
UDEC的基本公式假设为二维平面应变模型。此条件涉及断面保持为定值，并在平行于该断面的平面上作用荷载的无限长结构。所以，非连续面也被假设为平面特性。另外，UDEC提供了一个平面应力问题的选择。对于平面应变分析，如果在垂直于平面方向的应力?zz，为最大或最小主应力，在垂直于平面方向， 块体可能出现塑性屈服，
UDEC的显式求解算法允许进行动态或静态分析。对于动态计算，用户指定的速度或应力波可作为外部的边界条件或者内部激励直接输入到模型中。一个简单的动态波型库也可以获取。UDEC为动力分析设计了自由边界条件。
在静态分析中，包括了应力（力）和固定位移（速度为零）两种边界条件。边界条件在不同的位置可以是不同的。同时，在UDEC中还可以获得边界元边界，用于模拟无限弹性边界。也可以获得半平面解用来描述自由面效应。
UDEC还能够模拟通过模型中的孔隙和不连续面的流体流动。在此认为块体是不可渗透的。岩体的渗透率取决于节理的力学变形。也能够进行力学－流体全耦合分析。反过来，节理水压也将影响力学特性。流体被处理为平行板的粘性流。
程序中的结构单元可用于模拟岩体加固和工程表面支护。加固包括端部锚固、全长锚索和锚杆。表面支护模拟诸如喷射混凝土、混凝土衬砌和其他形式的隧道支护。 UDEC包含一个强有力的程序语言，FISH，能够使用户定义新的变量和函数。FISH是一个编辑器。通过UDEC数据文件进入程序被翻译并储存在内存中。
1.4 应用领域
UDEC最初是为节理岩石边坡的稳定性分析开发的。对于块体不连续公式和运动方
－3－
程（包括惯性项）采用显式时间步求解方法，便于块状岩体边坡的渐进破坏分析和大变形运动研究。
UDEC常用于采矿工程，已经进行了深部地下采矿洞室的静态与动态分析。洞室围岩破坏诱发的断裂、滑移是用UDEC分析研究的实例之一。通过在模型的边界施加动应力或速度波研究爆破影响。地震诱发的断层滑移也通过采用连续屈服节理模型进行了研究。结构单元已经用于模拟全长岩锚和喷射混凝土的各种岩体加固系统。
UDEC还应用于地下结构和深部高辐射废料的储存研究领域。通过应用热模型，UDEC已经应用于模拟与核废料相关的热荷载效应。
UDEC在作为一个计算设计工具，仍受到一定的限制。然而，程序较适用于研究节理效应的潜在破坏机理。节理岩体特性是一个“有限数据系统”－即，在很大程度上内部结构和应力状态是未知和不可知的。因此，建立一个完备的节理模型是不可能的。而且，UDEC是一个二维程序，除了特殊情况外，不可能表征具有三维结构的节理模型。不过，应用UDEC程序，可以从现象学的角度研究节理岩体地下工程开挖响应。该方法可加深岩石力学设计中各种不同现象的相互影响的理解。采用这种方法，工程师能够通过识别地下工程可能产生不可接受的变形或加载导致的破坏机理，从而揭示工程所潜在的诸多问题。
值得注意的是，UDEC程序对于模拟颗粒流动或动态分析火山喷发是不适宜的。对于该类研究，可以采用PFC2D程序。
2 开始启动
2.1 安装和启动程序
本节为首次使用UDEC的用户提供指导。如果你熟悉该程序仅仅是偶尔使用，你会发现本节尤其是2.6节对于改变你原有印象是有帮助的。UDEC程序共有65个主命令，有接近400个关键词。 2.1.7 内存赋值
UDEC自动调节内存大小达到8MB。可以通过下列命令查询、改变内存值： Uedc m Uedc 14 Print mem
如果更多的内存可以获得，其内存能够通过应用环绕磁盘文件获得额外内存。表2.2给出了最大块体数与所需内存的关系。
－4－