下载文档原格式
1、角点必须按顺时针方向排列；2、Crack 命令用于产生块体中单一直线特征的裂缝。
裂缝由端点坐标（x1，y1）和（x2，y2）所确定。
3、Jset 命令则是自动节理组生成器。
根据所给定的特征参数（即倾角、迹长、岩桥长度、间距和空间位置）产生一组裂缝。
4、round d---d是圆角距离，建议在block命令前指定圆角长度。
5、DELETE 命令，能从模型中删除一个块体。
例如，为了删除槽口块体，delete range 4.5，5.5 8，10。
6、GEN命令激活三角形网格有限单元自动生成器。
命令GEN edge v 将作用于任意形状的块体。
其v值定义三角形单元的最大边长，即v值越小，块体中的单元越小。
应当注意的是：具有高的边长比值的块体并不能产生单元，其极限的比重近似为1:10。
7、采用命令GEN quad v，指定模型为塑性材料模型的单元。
该类型的单元提供了对于塑性问题的精确解。
然而，GEN quad 命令可能对某些形状的块体不起作用。
在此情况下，应当采用GEN edge8、Change 命令改变块体为指定的变形块体。
Cons=0意味着模型块体材料被移出或开挖。
Cons=1 改变块体为各向同性弹性特性；而Cons=3则改变块体为摩尔－库仑模型，考虑塑性特性。
缺省值为所有变形体则自动改变为Cons=1。
P219、cha nge jcons=2，所以不连续结构面的缺省模型是Jcons=2。
10、可用以下命令检查材料号Plot block mat12、INSITU命令用来初始化应力。
采用该命令，可以赋值初始应力。
13、hist xvel 5, 5 hist ydisp 0, 11 第一个是记录位移坐标（x=5，y=5）附近结点x方向的速度，而第二个是记录接近坐标（x=0，y=11）位置处y方向的位移。
14、set grav 0.0 , -9.81第一个是x方向的加速度，第二个值为y方向的加速度为9.81m/sec2（向下作用）。
通用离散元用户指导（U D E C 3.1）2004.9目录1 引言 (1)1.1 总论 (1)1.2 与其他方法的比较 (2)1.3 一般特性 (2)1.4 应用领域 (3)2 开始启动 (4)2.1 安装和启动程序 (4)2.1.7 内存赋值 (4)2.1.9 运行UDEC (5)2.1.10 安装测试程序 (5)2.2 简单演示－通用命令的应用 (5)2.3 概念与术语 (6)2.4 UDEC模型：初始块体的划分 (8)2.5 命令语法 (9)2.6 UDEC应用基础 (10)2.6.1 块体划分 (10)2.6.2 指定材料模型 (16)2.6.2.1 块体模型 (16)2.6.2.2 节理模型 (17)2.6.3 施加边界条件和初始条件 (19)2.6.4 迭代为初始平衡 (21)2.6.5 进行改变和分析 (24)2.6.6 保存或恢复计算状态 (25)2.6.7 简单分析的总结 (25)2.8 系统单位 (26)3 用UDEC求解问题 (27)3.1 一般性研究 (27)3.1.1 第1步：定义分析模型的对象 (28)3.1.2 第2步：产生物理系统的概念图形 (28)3.1.3 第3步：建造和运行简单的理想模型 (28)3.1.4 第4步：综合特定问题的数据 (29)3.1.5 第5步：准备一系列详细的运行模型 (29)3.1.6 第6步：进行模型计算 (29)3.1.7 第7步：提供结果和解释 (30)3.2 产生模型 (30)3.2.1 确定UDEC模型合适的计算范围 (30)3.2.2 产生节理 (32)3.2.2.1 统计节理组生成器 (32)3.2.2.2 VORONOI多边形生成器 (34)3.2.2.3 例子 (34)3.2.3 产生内部边界形状 (35)3.3 变形块体和刚体的选择 (38)3.4 边界条件 (42)3.4.1 应力边界 (42)3.4.1.1 施加应力梯度 (43)3.4.1.2 改变边界应力 (44)3.4.1.3 打印和绘图 (44)3.4.1.4 提示和建议 (45)3.4.2 位移边界 (46)3.4.3 真实边界－选择合理类型 (46)3.4.4 人工边界 (46)3.4.4.1 对称轴 (46)3.4.4.2 截取边界 (46)3.4.4.3 边界元边界 (49)3.5 初始条件 (50)3.5.1 在均匀介质中的均匀应力：无重力 (50)3.5.2 无节理介质中具有梯度变化的应力：均匀材料 (51)3.5.3 无节理介质中具有梯度变化的应力：非均匀材料 (51)3.5.4 具有非均匀单元的密实模型 (52)3.5.5 随模型变化的初始应力 (53)3.5.6 节理化介质的应力 (54)3.5.7 绘制应力等值线图 (55)3.6 加载与施工模拟 (57)3.7 选择本构模型 (62)3.7.1 变形块体材料模型 (63)3.7.2 节理材料模型 (64)3.7.3 合理模型的选择 (65)3.8 材料性质 (71)3.8.1 岩块性质 (71)3.8.1.1 质量密度 (71)3.8.1.2 基本变形性质 (71)3.8.1.3 基本强度性质 (72)3.8.1.4 峰后效应 (73)3.8.1.5 现场性质参数的外延 (77)3.8.2 节理性质 (80)3.9 提示和建议 (81)3.9.1 节理几何形状的选择 (81)3.9.2 设计模型 (81)3.9.3 检查模型运行时间 (82)3.9.4 对允许时间的影响 (82)3.9.5 单元密度的考虑 (83)3.9.6 检查模型响应 (83)3.9.7 检查块体接触 (83)3.9.8 应用体积模量和剪切模量 (83)3.9.9 选择阻尼 (84)3.9.10 给块体和节理模型指定模型和赋值 (84)3.9.11 避免圆角误差 (85)3.9.12 接触嵌入 (85)3.9.13 非联结块体 (86)3.9.14 初始化变量 (86)3.9.15 确定坍塌荷载 (86)3.9.16 确定安全系数 (86)3.10 解释 (88)3.10.1 不平衡力 (88)3.10.2 块体/网格结点的速度 (88)3.10.3 块体破坏的塑性指标 (89)3.11 模拟方法 (89)3.11.1 有限数据系统模拟 (89)3.11.2 混沌系统的模拟 (90)3.11.3 局部化、物理的不稳定性和应力路径 (91)1 引言1.1 总论通用离散元程序（UDEC，Universal Distinct Element Code）是一个处理不连续介质的二维离散元程序。
(完整版)UDEC实例详解编辑整理：尊敬的读者朋友们：这里是精品文档编辑中心，本文档内容是由我和我的同事精心编辑整理后发布的，发布之前我们对文中内容进行仔细校对，但是难免会有疏漏的地方，但是任然希望（(完整版)UDEC实例详解）的内容能够给您的工作和学习带来便利。
同时也真诚的希望收到您的建议和反馈，这将是我们进步的源泉，前进的动力。
本文可编辑可修改，如果觉得对您有帮助请收藏以便随时查阅，最后祝您生活愉快业绩进步，以下为(完整版)UDEC实例详解的全部内容。
（完整版)UDEC实例详解编辑整理：张嬗雒老师尊敬的读者朋友们：这里是精品文档编辑中心，本文档内容是由我和我的同事精心编辑整理后发布到文库，发布之前我们对文中内容进行仔细校对，但是难免会有疏漏的地方,但是我们任然希望（完整版)UDEC实例详解这篇文档能够给您的工作和学习带来便利。
同时我们也真诚的希望收到您的建议和反馈到下面的留言区,这将是我们进步的源泉,前进的动力。
本文可编辑可修改,如果觉得对您有帮助请下载收藏以便随时查阅，最后祝您生活愉快业绩进步,以下为〈（完整版）UDEC实例详解> 这篇文档的全部内容。
实例1 荷载作用下边坡稳定问题朔准线黄河大桥，边坡高105m，层理间距取3m，岩层视倾角6°；主节理间距取2m，视倾角85°，计算模型如下:DEM计算模型相关计算参数取值如下：桥基长15m，宽12m，荷载180000kN；灰岩弹性模量E=4×104MPa，泊松比μ=0.23，岩体密度γ=2.7×103kg/m3。
newround 0.05set delc offblock 0 0 0 50 50 50 62。
6 115。
2 68.9 123。
9 85。
5 132。
5 ＆105.5 140.1 125.4 156.4 180 156。
4 180 0jset 85,0 300，0 0 0 3，0 (50 50）jset 6 0 300 0 0 0 4 0 （50 50)gen auto 5change jmat=1 range angle 84 86change jmat=2 range angle 5 7prop mat=1 den=2700 b=1。
科技写作学习心得（共5篇）第一篇：科技写作学习心得《科技论文写作》学习心得经过两周时间的学习，《科技论文写作》这门课已经结束，在此谈一下个人学习心得。
《科技论文写作》主要讲授科技论文的特点、作用，科技论文写作过程中的选题、资料检索、表达方式及语言的运用、科技创新思维以及学术论文的写作与发表，学位论文的写作与答辩等内容。
所以对我们单证研究生而言是极为实用的，也是极为重要的。
学习了这门课程，我感觉自己在以下三个方法的水平得到了提高：1.知识方面：通过本课程的学习，掌握科技论文写作的基本原则及常用方法；掌握将试验研究结果加工整理为符合规范要求的科技论文的方法；掌握论文摘要、正文和参考文献的写作技巧，等等。
2.能力方面：通过课程教学，提高进行科技论文写作的能力；具备将试验研究结果加工整理为符合规范要求的科技论文的能力,以及具备科技论文编辑的基本能力等。
3.素质方面：通过课程教学，培养了严谨求实、理论联系实际的科学态度；培养了独立地获取知识和综合运用知识的能力；培养了分析问题和解决问题的能力；提高了综合素质。
通过学习这门课程，我学到了好多东西。
首先，我明白了毕业论文的特点与基本要求以及毕业论文的写作宗旨。
然后，知道怎么去选择毕业论文题目，以及明白了毕业论文的选题原则和毕业论文的选题方法。
其次，知道在毕业论文写作的前期需要做一些准备。
例如，收集、筛选和整理资料；选择论文的类型；拟定论文的写作大纲。
再其次，毕业论文的项目构成及说明，它又包括毕业论文的项目构成，标题、目录、摘要、关键词和注释与参考文献。
最后，就是正文的写作，主要内容是引言、主体、结论和论文的自查。
科技论文区别于其他文体的特点：科技论文是创新性科学技术研究工作成果的科学论述，是某些理论性、实验性或观测性新知识的科学记录，是某些已知原理应用于实际中取得新进展、新成果的科学总结。
它是反映文稿的科学依据和著者尊重他人研究成果而向读者提供文中引用有关资料的出处，或为了节约篇幅和叙述方便，提供在论文中提及而没有展开的有关内容的详尽文本。
计算机学习总结2022计算机学习总结范文7篇2022计算机学习总结范文1在计算机的发展史上，有三件事是我觉得最有转折性意义的，第一就是冯诺依曼进行的计算机逻辑设计，他运用其非凡的分析、综合能力和深厚的数学基础编写出了EDVAC方案，对EDVAC的总体结构和逻辑设计的优化起到了关键的作用。
EDVAC方案当中有两个重大的改进：一是用二进制代替十进制，便于电子元件表示数据，简化了计算器的运算速度；二是提出了“存储程序”的概念，程序和数据都存放在存储器中，实现了基于程序的计算机自动执行，实现了程序执行中的“条件转移”。
这两种改变基本上奠定了现代计算机的基本模型，时至今日，现代的计算机依然没有突破冯诺依曼所设计的基本结构。
第二件重大的事是微处理器的出现，微处理器是将运算器和控制器集成在一起的大规模/超大规模集成电路芯片。
微处理器的出现不仅提高了计算机的运算速度，而且实现了计算机的微型化，因此计算机得以进入广大人民群众的生活里，成为他们生活、娱乐、工作的一个必不可少的工具，从而拓宽了计算机的使用领域和使用价值。
但在我看来，计算机微型化的影响不仅仅只是使人们的生活更加便利，而是让更多有思想的人接触到了这一改变未来的伟大发明，当他们接触到这项发明后被吸引，他们的兴趣被激发，他们愿意把自己的时间、智慧投入这项发明里，因此，计算机热潮开始席卷而来，每天都会有无数的创意、更新出现。
这种热潮、风气大大地推动了计算机的完善，最有代表性的就是Windows的诞生。
第三件就是互联网的出现，其中影响最大的毫无疑问就是因特网。
计算机网络是自主计算机的互联集合，主要包括个人区域网、局域网、城域网、广域网、互联网等形式。
因特网得到广泛的应用，为人们的工作、学习提供了网络新闻、信息检索、电子邮件、博客、文件传输、电子商务等多种服务，渐渐地改变了人们的日常生活，甚至影响着科技的发展。
各个不同地区之间的文化得以相互传播。
2022计算机学习总结范文2通过近三个多月的学习，我收获了很多意外的惊喜。
关于DCS操作学习总结（通用3篇）DCS操作学习总结1 首先感谢蔡总和郁主任给我提供了这么好的一个学习机会去上海参加横河CS3000 DCS系统组态的培训。
使我能在更好的环境中更加系统的学习我公司所使用到的DCS系统。
这是我毕业大约5年时间第一次脱产学习，我非常珍惜这次学习机会，在十多天的学习生活中，我始终保持较高的学习热情，争取能更好的为公司服务。
2017年3月4日，我受公司安排前往上海横河电机中国培训中心学习CS3000 DCS系统组态。
在这十多天的学习过程中，通过课程学习，了解了DCS系统的硬件构成、软件安装，以及项目创建、常规反馈、顺序控制和人机界面定义等内容。
DCS是分布式控制系统的英文缩写（Distributed Control System），在国内自控行业又称之为集散控制系统。
是相对于集中式控制系统而言的一种新型计算机控制系统，它是在集中式控制系统的基础上发展、演变而来的。
它是一个由过程控制级和过程监控级组成的以通信网络为纽带的多级计算机系统，综合了计算机，通信、显示和控制等4种技术，其基本思想是分散控制、集中操作、分级管理、配置灵活以及组态方便。
首先，DCS的骨架—系统网络，它是DCS的基础和核心。
由于网络对于DCS整个系统的实时性、可靠性和扩充性，起着决定性的作用，因此各厂家都在这方面进行了精心的设计。
对于DCS的系统网络来说，它必须满足实时性的要求，即在确定的时间限度内完成信息的传送。
这里所说的“确定”的时间限度，是指在无论何种情况下，信息传送都能在这个时间限度内完成，而这个时间限度则是根据被控制过程的实时性要求确定的。
因此，衡量系统网络性能的指标并不是网络的速率，即通常所说的每秒比特数（bps），而是系统网络的实时性，即能在多长的时间内确保所需信息的传输完成。
系统网络还必须非常可靠，无论在任何情况下，网络通信都不能中断，因此多数厂家的DCS均采用双总线、环形或双重星形的网络拓扑结构。
《软件单元测试》读书札记目录一、内容概览 (2)二、软件单元测试概述 (3)1. 软件单元测试定义与重要性 (4)2. 软件单元测试发展历程 (5)3. 软件单元测试基本原则 (6)三、软件单元测试基础概念 (7)1. 单元测试的级别 (9)2. 测试方法 (10)3. 测试工具与框架 (12)四、软件单元测试流程 (14)1. 测试计划制定 (14)2. 测试用例设计 (15)3. 测试环境搭建 (17)4. 测试执行与监控 (18)5. 测试报告编写 (19)五、软件单元测试的实战技巧 (20)1. 测试数据的准备与处理 (22)2. 单元测试中的常见问题及解决方案 (23)3. 单元测试与集成测试的协同工作 (23)4. 持续集成与自动化测试实践 (25)六、软件单元测试的最新趋势与挑战 (26)1. 人工智能与软件测试的融合 (28)2. 云计算环境下的软件测试挑战 (30)3. 大数据时代的软件测试技术创新 (31)七、软件单元测试的未来发展展望 (32)1. 测试技术的不断更新换代 (34)2. 自动化测试工具的进一步完善 (35)3. 测试人员的技能提升与职业发展路径 (37)八、结语及心得体会 (38)1. 本书阅读总结及感悟 (39)2. 对软件单元测试的进一步思考与实践计划 (41)一、内容概览引言：简要介绍了软件单元测试的重要性，以及为什么需要进行单元测试。
同时阐述了单元测试的基本概念，为后续内容做了铺垫。
单元测试基础：详细介绍了单元测试的基本概念、原则、流程和测试方法。
包括单元测试的定义、目的、测试范围、测试策略等，为读者打下了坚实的基础。
测试技术与工具：详细讲解了常用的软件单元测试技术和工具，如测试框架、断言库、模拟和伪造技术等。
这些技术和工具在单元测试中发挥着重要作用，有助于提高测试效率和质量。
实战案例：通过多个实际项目的案例分析，展示了如何应用单元测试的流程和技巧。
docker实训总结万能版500字心得【最新版4篇】目录（篇1）一、Docker 简介二、Docker 实训的目标与收获三、Docker 的安装与配置四、Docker 的实战应用五、Docker 的未来发展趋势与个人展望正文（篇1）一、Docker 简介Docker 是一种开源的容器技术，它能让开发者打包他们的应用以及依赖包到一个可移植的容器中，然后发布到任何流行的 Linux 机器或Windows 机器上，也可以实现虚拟化。
容器是完全使用沙箱机制，相互之间不会有任何接口。
二、Docker 实训的目标与收获通过 Docker 实训，我旨在掌握 Docker 的基本概念、安装与配置方法，以及在实际项目中的应用。
实训过程中，我学会了如何创建、运行和管理 Docker 容器，熟悉了 Dockerfile 和 Docker-compose 等工具的使用，同时也了解了一些 Docker 的高级特性。
通过这次实训，我对Docker 技术有了更加深入的理解和实践经验，对未来的工作有很大帮助。
三、Docker 的安装与配置在开始 Docker 实训之前，首先需要安装 Docker。
根据操作系统的不同，安装方法也有所区别。
在 Linux 系统中，可以通过运行以下命令来安装 Docker：```sudo apt-get updatesudo apt-get install docker.io```在 Windows 系统中，需要从 Docker 官网下载安装包并运行安装程序。
安装完成后，可以通过运行`docker version`命令来检查 Docker 的版本信息。
接下来是 Docker 的配置。
可以通过编辑`/etc/docker/daemon.json`文件来设置 Docker 的默认参数，例如：```{"log-level": "warn","storage-driver": "overlay2"}```四、Docker 的实战应用在实际项目中，Docker 的应用非常广泛。
UDEC 3.0 版本3．利用UDEC解决的问题本章就UDEC在处理岩石力学工程问题中的使用给予说明。
在3.1部分（节）中，就实施地质力学分析可行步骤作了简述。
接着从3.2节到3.10节就有关在任何模型创建和解决方法中具体（特定）方面的检测中作了介绍。
内容包括：●模型的产生（3.2节）；●刚性块或可变形块分析的选择（3.3节）；●边界与初始条件（3.4节和3.8节）；●负载连续模型（3.6节）；●块体和节理构成的模型选择和材料特性（3.7节和3.8节）；●改善模型效率的方法（3.9节）；●结果说明（3.10节）。
最后，地质力学领域建模思想将在3.11节中受到检验；在这领域中的初学建模者可能希望首先参考这节。
在地质力学中建模方法与其他的工程领域中有显著差异，例如结构工程学。
在进行任何地质力学分析时应注意这一点很重要。
3.1概述（总的手段）地质工程建模过程涉及到特定的考虑和不同于上述提到的构造物体的设计理念。
对于在岩土内或上面的建筑和巷道的分析和设计必须用小到微小具体位置的数据，并注意变形和强度特性可能的相当变化。
在一个岩石或土地位置取得完全的现场数据是不可能的。
例如，有关应力、岩性和不连续性的信息最多仅能部分了解。
由于有关对于设计必须的输入数据是有限的，那么地质力学数学模型应主要用来理解影响整个系统行为的主导机理过程。
一但系统行为被理解，那么就可以具体（特定）为设计过程开发出简单的预测方案。
该方法适应于地质力技术工程，在这里总是缺乏好的数据，但在其他方面的应用中如果有充足数据，在设计上直接运用UDEC也许是可以的。
当程序提供适当的数据，利用UDEC图3.1建模位置范围UDEC可以用在一种完全可预测的方式中（如图3.1右侧部分）或作为一个“数值实验室”来验证观点（如左侧部分）。
决定使用类型是场地情况（和预算）而不是程序。
如果有充足的高质量数据可找到，UDEC可给出良好的预测。
既然多数UDEC应用适合于仅有较小数据可找到的情况，本节讨论这种推荐方法用来处理数学模型，好象是一个实验室实验。
1.把图形保存下来，能在AUTOCAD(图片)中打开、编辑plot block cable red supp ye stru bl;显示块体、锚杆（红色）、支架（黄色）、梁（蓝色）的图形set plot dxf(jpg) 256；设置图形为256色set out c:\ss.dxf(jpg) (set out 1.dxf 这个文件保存在当前带数值模拟的文件夹里)copy c:\ss.dxf(jpg) ( copy 1.dxf 这个文件保存在当前带数值模拟的文件夹里)；把这图形以ss.dxf文件保存在C磁盘下2.把数据导出set log onprint pline 1 ydis;把pline 1的y方向位移的数据导出到UDEC目录下udec.txt文件中set log off3.plot block stress ;显示块体的应力plot block dis ;显示块体的位移(有x和y方向)plot block pl ;显示块体的塑性区(plastic)save xx.save ;保存计算结果res xx.save ;调用4.液压支架的命令supp xc yc wid l seg n mat j;(xc,yc)是指中心点坐标,l指支架的宽度,n指分段数,mat j指支柱材料性质为j prop mat j sup_kn -1 ;sup_kn指支架的刚度,-1与表1相对应del range x1 x2 y1 y2 ;挖掘范围(x1, y1)(x1 ,y2)(x2, y2)(x2, y1)的块体,由支架支撑table 1 0 42.e6 0.05 5.0e6 0.1 6.0e6;表1表示的是液压支架的(P-DS)特性曲线5 巷道施工中断面加梯子梁的命令(见图1)stru gen xc yc np 100 fa a thetra b mat=16 thick=0.2 ;100 指分100段,a b 指角度prop mat=16 st_ymod=13.5e9 st_prat=0.14 st_den=7800prop mat=16 st_yield=6e7 st_yresid=6e7 st_ycomp=2.5e7prop mat=16 if_kn=1.35e9 if_ks=1.35e9 if_tens=0 if_fric=18 if_coh=0;interface-界面例子讲解:建模(以米为单位,;后为解释部分)round=0.1 ;方块的圆角块半径为0.1米set ovtol=1.0 ;块体与块体之间相互嵌入量最大值为1米bl 0,0 0,26 50,26 50,0 ;在(x1, y1)(x1 ,y2)(x2, y2)(x2, y1)生成块体范围crack 0,25 50,25 ;在(x1, y1)(x2 ,y2)两点间画直线jregion id 1 0,0 0,10 50,10 50,0 delete ;删除此区域的块体jset 90,0 2.5,0 2.5,0 3,0 0,0 range jregion 1 ;jset A,0 a,0 b,0 c,0 x0,y0 range jregion 1jset 90,0 2.5,0 2.5,0 3,0 1.5,2.5 range jregion 1 可画成列的线段(见图2)jset 0,0 26,0 0,0 2.5,0 0,0 range jregion 1pa ;当程度运行至此时暂停,可以看看你所建的部分模型,用continue继续运行下面部分save t.save ;建模保存在t.save中岩层赋属性及原岩力平衡计算res t.save ;调用已建好的模型gen quad 12 ;定义块体最大变形,若没有此语句,刚所有块体均为刚性块体zone model mo range 0,180 0,45.2 ;在范围(x1, y1)(x1 ,y2)(x2, y2)(x2, y1)的块体符合库仑准则change jcons=2 range 0,180 0,45.2 ;节理面间接触-coulmb滑移;下面主要是讲岩层赋属性change （范围） mat=1 range reg 0,0 0,3 180,3 180,0 ;岩层1的范围 change mat=2 range reg 0,3 0,4.5 180,4.5 180,3change mat=3 range reg 0,4.5 0,10.2 180,10.2 180,4.5 change mat=4 range reg 0,10.2 0,11.2 180,11.2 180 10.2 prop （赋值的意思） mat=1 dens=2500 ;岩层1的密度prop mat=2 dens=2500 prop mat=3 dens=1300zone （块体） k=2e10, g=1.1e10, fric=30,coh=2e6,ten=4e6 range mat=1 ;岩层1的块体力学参数(见表1) zone k=8.7e9, g=4.2e9, fric=25,coh=1.5e6,ten=1.5e6 range mat=2 zone k=13.05e9, g=6.3e9, fric=31,coh=4e6,ten=2e6 range mat=3change jmat(节理)=1 range reg 0,0 0,3 180,3 180,0 ;在这范围1内的块体间相互接触的接触面力学参数 change jmat=2 range reg 0,3 0,4.5 180,4.5 180,3change jmat=3 range reg 0,4.5 0,10.2 180,10.2 180,4.5prop jmat=1 jkn=7e9, jks=2e9, jcoh=0.1e6 , jfric=0,jten=0.1e6 ;在这范围1内的块体间相互接触的接触面力学参数 prop jmat=2 jkn=6e9, jks=2e9, jcoh=1e6 , jfric=20,jten=1e6 prop jmat=3 jkn=6e9, jks=2e9, jcoh=1e6 , jfric=20,jten=1e6;工作面埋深550米set gravity 0,-10 ;地下岩层主要受重力,还有构造应力(水平应力=垂直应力*侧压系数)bound stress 0,0,-1.375e7 range 0 180 45.1 45.3 ; boundary stress sxx0 ,sxy0, syy0 range x1 x2 y1 y2;(550—煤层到上边界的距离)*2500*(-10)= -1.375e7 insitu str -1.5144e7 0 -1.262e7 szz=-1.5144e7 ygrad 3e4 0 2.5e4 zgrad 0 3e4;initi asituate stress 初始设置应力 ,侧压系数为1.2;syy=(550+煤层到下边界的距离)*2500*(-10)= -1.262e7（在下边界上施加的力）-1.262e7*1.2=-1.5144e7 ,y 方向的梯度为=1*2500*10=2.5e4, x 、z 方向的梯度为=2.5e4*1.2=3.0e4bound xvel=0 range -0.1 0.1 0 45.2 ;固定左边界(xvel-也就是x 方向的速度),见图3 bound xvel=0 range 179.9 180.1 0 45.2 ;固定右边界bound yvel=0 range 0 180 -0.1 0.1 ;固定下边界（因为上边界有上覆岩层故不需固定了） solve ;计算save fyuanyan1.save注: insitu str sxx sxy syy szz=-1.6638e8 ygrad sxxy sxyy syyy zgrad szzx szzy (xgrad sxxx sxyx syyy,上面没有写,表示在x 方向没有变化 )(σx=sxx ,σy=syy, σz=szz,τxy=sxy ,τxz=sxz, τyz=syz 弹性力学;grad 表示梯度,即在此方向的变化量) (xgrad sxxx sxyx syyy ygrad sxxy sxyy syyy zgrad szzx szzy(sxx=sxx0+sxxx.x+sxxy .y,syy=syy0+syyx.x+syyy.y ,szz=szz0+szzx.x+szzy.y) (μ—侧压系数,sxx=syy . μ,szz=syy . μ,sxxy=syyy . μ,szzy=syyy. μ)(τ=c+σ.tg φ,c 值,反映岩石剪切时的粘结阻力,故称岩石的内聚力(或粘结力) ,φ,值反映岩石剪切时摩擦阻力, 的大小,故称岩石的内摩擦角,tg φ,相当于摩擦系数f 。
可得用c 、φ两个数在应力圆中判断某种应力状态下的岩石是否产生破坏，通常，岩石愈坚硬，c 、φ两个值超大，反之亦然。
1G=103M=109()21E G υ=+，()312EK υ=-】锚杆的效果显不著。
round=0.01set ovtol=3bl 0,0 0,45.2 180,45.2 180,0;基本底3m(1.5*1)crack 0,3 180,3crack 57,3 57,0crack 73.8,3 73.8,0jregion id 11 0,0 0,3 57,3 57,0 deletejset 0,0 57,0 0,0 1,0 0,0 range jregion 11jset 90,0 1,0 1,0 1.5,0 0,0 range jregion 11jset 90,0 1,0 1,0 1.5,0 0.75,1 range jregion 11jregion id 12 57,0 57,3 73.8,3 73.8,0 deletejset 0,0 16.8,0 0,0 0.25,0 57,0 range jregion 12jset 90,0 0.25,0 0.25,0 0.5,0 57,0 range jregion 12jset 90,0 0.25,0 0.25,0 0.5,0 57.25,0.25 range jregion 12 jregion id 13 73.8,0 73.8,3 180,3 180,0 deletejset 0,0 106.2,0 0,0 1,0 73.8,0 range jregion 13jset 90,0 1,0 1,0 1.5,0 73.8,0 range jregion 13jset 90,0 1,0 1,0 1.5,0 74.55,1 range jregion 13;直接底1.5m(1*0.5/0.5*0.25)crack 0,4.5 180,4.5crack 57,3 57,4.5crack 73.8,3 73.8,4.5;前57米(1*0.5)jregion id 21 0,3 0,4.5 57,4.5 57,3 deletejset 0,0 57,0 0,0 0.5,0 0,3 range jregion 21jset 90,0 0.5,0 0.5,0 1,0 0.5,3 range jregion 21jset 90,0 0.5,0 0.5,0 1,0 0,3.5 range jregion 21;57到73.8米(0.5*0.25)jregion id 22 57,3 57,4.5 73.8,4.5 73.8,3 deletejset 0,0 16.8,0 0,0 0.25,0 57,3 range jregion 22jset 90,0 0.25,0 0.25,0 0.5,0 57,3 range jregion 22jset 90,0 0.25,0 0.25,0 0.5,0 57.25,3.25 range jregion 22 ;后73.8到180米(1*0.5)jregion id 23 73.8,3 73.8,4.5 180,4.5 180,3 deletejset 0,0 106.2,0 0,0 0.5,0 73.8,3 range jregion 23jset 90,0 0.5,0 0.5,0 1,0 73.8,3 range jregion 23jset 90,0 0.5,0 0.5,0 1,0 74.3,3.5 range jregion 23;5.7米的煤层，采3.5米crack 0,10.2 180,10.2crack 0,9.5 57,9.5crack 73.8,9.5 180,9.5crack 57,4.5 57,10.2crack 73.8,4.5 73.8,10.2crack 63,4.5 63,7.75crack 67.8,4.5 67.8,10.2crack 57,7.75 67.8,7.75crack 67.8,8 180,8;前57米(1*0.5)jregion id 31 0,4.5 0,9.5 57,9.5 57,4.5 deletejset 0,0 57,0 0,0 0.5,0 0,4.5 range jregion 31jset 90,0 0.5,0 0.5,0 1,0 0,4.5 range jregion 31jset 90,0 0.5,0 0.5,0 1,0 0.5,5 range jregion 31jregion id 30 0,9.5 0,10.2 57,10.2 57,9.5 deletejset 90,0 0.7,0 0.7,0 1,0 0,9.5 range jregion 30;57到73.8米(0.5*0.25)jregion id 321 57,4.5 57,7.75 63,7.75 63,4.5 deletejset 0,0 6,0 0,0 0.25,0 57,4.5 range jregion 321jset 90,0 0.25,0 0.25,0 0.5,0 57,4.5 range jregion 321 jset 90,0 0.25,0 0.25,0 0.5,0 57.25,4.75 range jregion 321 jregion id 322 63,4.5 63,7.75 67.8,7.75 67.8,4.5 delete jset 0,0 3.2,0 0,0 0.25,0 63,4.5 range jregion 322jset 90,0 0.25,0 0.25,0 0.5,0 63,4.5 range jregion 322 jset 90,0 0.25,0 0.25,0 0.5,0 63.25,4.75 range jregion 322 jregion id 323 57,7.75 57,10.2 67.8,10.2 67.8,7.75 delete jset 0,0 10.8,0 0,0 0.25,0 57,7.75 range jregion 323jset 90,0 0.25,0 0.25,0 0.5,0 57.25,7.75 range jregion 323 jset 90,0 0.25,0 0.25,0 0.5,0 57,8 range jregion 323 jregion id 324 67.8,4.5 67.8,8 73.8,8 73.8,4.5 deletejset 0,0 6,0 0,0 0.25,0 67.8,4.5 range jregion 324jset 90,0 0.25,0 0.25,0 0.5,0 68,4.5 range jregion 324jset 90,0 0.25,0 0.25,0 0.5,0 68.25,4.75 range jregion 324 jregion id 325 67.8,8 67.8,10.2 73.8,10.2 73.8,8 delete jset 0,0 6,0 0,0 0.25,0 67.8,4.5 range jregion 325jset 90,0 0.25,0 0.25,0 0.5,0 68,8 range jregion 325jset 90,0 0.25,0 0.25,0 0.5,0 68.25,8.25 range jregion 325 ;后73.8到180米(1*0.5)jregion id 331 73.8,4.5 73.8,9.5 180,9.5 180,4.5 delete jset 0,0 106.2,0 0,0 0.5,0 73.8,4.5 range jregion 331jset 90,0 0.5,0 0.5,0 1,0 74.3,4.5 range jregion 331jset 90,0 0.5,0 0.5,0 1,0 73.8,5 range jregion 331jregion id 300 73.8,9.5 73.8,10.2 180,10.2 180,9.5 delete jset 90,0 0.7,0 0.7,0 1,0 74.3,9.5 range jregion 300;直接顶1m(1.5*1)crack 0,11.2 180,11.2crack 57,10.2 57,11.2crack 73.8,10.2 73.8,11.2;前57米(1*0.5)jregion id 41 0,10.2 0,11.2 57,11.2 57,10.2 delete;jset 0,0 57,0 0,0 0.5,0 0,10.2 range jregion 41jset 90,0 1,0 1,0 1.5,0 0.75,10.2 range jregion 41;jset 90,0 0.5,0 0.5,0 1,0 0,10.7 range jregion 41;57到73.8米(0.5*0.25)jregion id 42 57,10.2 57,11.2 73.8,11.2 73.8,10.2 delete jset 0,0 16.8,0 0,0 0.25,0 57,10.2 range jregion 42jset 90,0 0.25,0 0.25,0 0.5,0 57.25,10.2 range jregion 42 jset 90,0 0.25,0 0.25,0 0.5,0 57,10.45 range jregion 42;后73.8到180米(1.5*1)jregion id 43 73.8,10.2 73.8,11.2 180,11.2 180,10.2 delete ;jset 0,0 106.2,0 0,0 0.5,0 73.8,10.2 range jregion 43;jset 90,0 0.5,0 0.5,0 1,0 73.8,10.2 range jregion 43jset 90,0 1,0 1,0 1.5,0 73.8,10.2 range jregion 43;基本顶10m(12*5)crack 0,21.2 180,21.2jregion id 5 0,11.2 0,21.2 180,21.2 180,11.2 deletejset 0,0 180,0 0,0 5,0 0,11.2 range jregion 5jset 90,0 5,0 5,0 12,0 0,11.2 range jregion 5jset 90,0 5,0 5,0 12,0 6,16.2 range jregion 5;上覆岩层24m(6*4)crack 0,45.2 180,45.2jregion id 6 0,21.2 0,45.2 180,45.2 180,21.2 deletejset 0,0 180,0 0,0 4,0 0,21.2 range jregion 6jset 90,0 4,0 4,0 5,0 3,21.2 range jregion 6jset 90,0 4,0 4,0 5,0 0,25.2 range jregion 6save inir.saveres inir.savegen quad 1.6 1.2 range 0,180 0,11.2gen quad 13 5.5 range 0,180 11.2,45.2zone model mo range 0,180 0,45.2 ;( change cons=3 range 0,180 0,45.2)change jcons=2 range 0,180 0,45.2;设定不同岩层、煤层change mat=1 range reg 0,0 0,3 180,3 180,0change mat=2 range reg 0,3 0,4.5 180,4.5 180,3change mat=3 range reg 0,4.5 0,10.2 180,10.2 180,4.5change mat=4 range reg 0,10.2 0,11.2 180,11.2 180 10.2change mat=5 range reg 0,11.2 0,21.2 180,21.2 180,11.2change mat=6 range reg 0,21.2 0,45.2 180,45.2 180,21.2change jmat=1 range reg 0,0 0,3 180,3 180,0change jmat=2 range reg 0,3 0,4.5 180,4.5 180,3change jmat=3 range reg 0,4.5 0,10.2 180,10.2 180,4.5change jmat=4 range reg 0,10.2 0,11.2 180,11.2 180 10.2change jmat=5 range reg 0,11.2 0,21.2 180,21.2 180,11.2change jmat=6 range reg 0,21.2 0,45.2 180,45.2 180,21.2prop mat=1 dens=2500prop mat=2 dens=2500prop mat=3 dens=1390prop mat=4 dens=2500prop mat=5 dens=2700prop mat=6 dens=2500zone k=2e10, g=1.1e10, fric=40, coh=8.5e6, ten=5.6e6 range mat=1 (prop mat=1 k=2e10, g=1.1e10, fric=40, coh=8.5e6, ten=5.6e6 )zone k=17e9, g=8.9e9, fric=34, coh=6e6, ten=2.4e6 range mat=2zone k=13.05e9, g=6.3e9, fric=28, coh=4e6, ten=2.1e6 range mat=3 zone k=2e10, g=1.1e10, fric=33, coh=6.2e6, ten=2.3e6 range mat=4 zone k=2.4e10, g=1.3e10, fric=40, coh=9.5e6, ten=5.4e6 range mat=5 zone k=2e10, g=1.1e10, fric=35, coh=8.4e6, ten=4e6 range mat=6prop jmat=1 jkn=14e9, jks=8.7e9, jcoh=6.5e6 , jfric=35,jten=0prop jmat=2 jkn=9.5e9, jks=6.4e9, jcoh=4.8e6 , jfric=28,jten=0prop jmat=3 jkn=6e9, jks=4.7e9, jcoh=2.8e6 , jfric=20,jten=0prop jmat=4 jkn=8.9e9, jks=5.9e9, jcoh=4.9e6, jfric=27, jten=0prop jmat=5 jkn=16e9, jks=9.5e9, jcoh=7.8e6 , jfric=36,jten=0prop jmat=6 jkn=14e9, jks=8.7e9, jcoh=6.5e6 , jfric=29,jten=0;工作面埋深550米set gravity 0,-10bound stress 0,0,-1.375e7 range 0 180 45.1 45.3insitu str -1.722e7 0 -1.435e7 szz=-1.722e7 range 0 180 21.2 45.2insitu ygrad 3e4 0 2.5e4 zgrad 0 3e4 range 0 180 21.2 45.2insitu str -1.7544e7 0 -1.462e7 szz=-1.7544e7 range 0 180 11.2 21.2insitu ygrad 3.24e4 0 2.7e4 zgrad 0 3.24e4 range 0 180 11.2 21.2insitu str -1.7574e7 0 -1.4645e7 szz=-1.7574e7 range 0 180 10.2 11.2insitu ygrad 3e4 0 2.5e4 zgrad 0 3e4 range 0 180 10.2 11.2insitu str -1.7669076e7 0 -1.472423e7 szz=-1.7669076e7 range 0 180 4.5 10.2 insitu ygrad 1.668e4 0 1.39e4 zgrad 0 1.668e4 range 0 180 4.5 10.2insitu str -1.7804076e7 0 -1.483673e7 szz=-1.7804076e7 range 0 180 0 4.5 insitu str 3e4 0 2.5e4 zgrad 0 3e4 range 0 180 0 4.5bound xvel=0 range -0.1 0.1 0 45.2bound xvel=0 range 179.9 180.1 0 45.2bound yvel=0 range 0 180 -0.1 0.1solvesave yuanyan.savecal kh.txtres yuanyan.savereset dispreset vel;开掘巷道并支护巷道del range reg 63,4.5 63,7.75 67.8,7.75 67.8,4.5;加顶梯子梁(顶梁4.6米)stru gen xc 65.4 yc 4.95 np 100 fa 0 theta 180 mat=16 thick 0.2 connect prop mat=16 st_ymod=13.5e9 st_prat=0.14 st_den=7800prop mat=16 st_yield=60e6 st_yresid=60e6 st_ycomp=35e6prop mat=16 if_kn=1.35e9 if_ks=1.35e9 if_tens=0 if_fric=18 if_coh=0;两帮梯子梁(2.75米);stru gen xc 65.4 yc 6.15 np 100 fa 150.64 theta 59.62 mat=17 thick 0.2;prop mat=17 st_ymod=21e9 st_prat=0.15 st_den=7800;prop mat=17 st_yield=40e6 st_yresid=40e6 st_ycomp=40e6;prop mat=17 if_kn=1e9 if_ks=1e9 if_tens=0 if_fric=50 if_coh=0;stru gen xc 65.4 yc 6.15 np 100 fa 329.75 theta 59.62 mat=18 thick 0.2;prop mat=18 st_ymod=21e9 st_prat=0.15 st_den=7800;prop mat=18 st_yield=40e6 st_yresid=40e6 st_ycomp=40e6;prop mat=18 if_kn=1e9 if_ks=1e9 if_tens=0 if_fric=50 if_coh=0;顶锚杆(L2.4米,6根)cable (63.2,7.75) (62.379,10.005) 5 10 314e-6 11 4e4 connectcable (64.08,7.75) (64.08,10.15) 5 10 314e-6 11 4e4 connectcable (64.96,7.75) (64.96,10.15) 5 10 314e-6 11 4e4 connectcable (65.84,7.75) (65.84,10.15) 5 10 314e-6 11 4e4 connectcable (66.72,7.75) (66.72,10.15) 5 10 314e-6 11 4e4 connectcable (67.6,7.75) (68.421,10.005) 5 10 314e-6 11 4e4 connect;左帮锚杆(L2米,4根)cable (63,4.95) (61,4.95) 5 10 314e-6 11 4e4 connectcable (63,5.8) (61,5.8) 5 10 314e-6 11 4e4 connectcable (63,6.65) (61,6.65) 5 10 314e-6 11 4e4 connectcable (63,7.5) (61.121,8.184) 5 10 314e-6 11 4e4 connect;右帮锚杆(L2米,4根)cable (67.8,4.95) (69.8,4.95) 5 10 314e-6 11 4e4 connectcable (67.8,5.8) (69.8,5.8) 5 10 314e-6 11 4e4 connectcable (67.8,6.65) (69.8,6.65) 5 10 314e-6 11 4e4 connectcable (67.8,7.5) (69.679,8.184) 5 10 314e-6 11 4e4 connect;锚杆参数20*2000prop m=10 cb_dens 7500 cb_ycomp 430e6 cb_yield 260e3 cb_ymod 1.2e11 ;锚固体的参数prop m=11 cb_kbond=6.3e9 cb_sbond=6e6cable (64.5,7.75) (64.5,14.95) 5 12 314e-6 13 4e4cable (65.8,7.75) (65.8,14.95) 5 12 314e-6 13 4e4;锚索参数20*7200prop m=12 cb_dens 7500 cb_ycomp 700e6 cb_yield 500e3 cb_ymod 1.2e11 ;锚固体的参数prop m=13 cb_kbond=6.3e9 cb_sbond=6e6chang mat=7 range region 60,4.5 60,10.2 71,10.2 71,4.5prop mat=7 dens=1390zone k=14.05e9, g=6.9e9,fric=37, coh=4.4e6, ten=2.1e6 range mat=7solvesave kh.savecall kgx.txtres kh.savereset dispreset vel;采3.5米煤层del range reg 67.8,4.5 67.8,8 167,8 167,4.5stru del range 67.7,68 4.5,7.75cable del range 67.8,70 4.5 7.8def aaxx=68.58supp xx 4.6 wid 1.5 seg 3 mat 1 endcommandxx=xx+1.55end_loopendaasupp 65.4 4.8 angle 90 wid 4.2 seg 3 mat 2 prop m 1 sup_kn -1 sup_tmax=1000prop m 2 sup_kn -2table 1 0 4.2e6 0.05 5.0e6 0.1 6.0e6table 2 0 2.5e6 0.2 3.2e6set pline 60,11.2 72,11.2 12set pline 63,7.75 67.8,7.75 12set pline 60,11.2 72,11.2 12set pline 63,4.5 63,7.75 8set pline 67.8,4.5 67.8,7.75 8;顶板下沉量hist ydis 63.3,7.75hist ydis 67.5,7.75hist ydis 65.4,7.75;底板鼓起量hist ydis 63.3,4.5hist ydis 67.5,4.5hist ydis 65.4,4.5;左帮水平移近量hist xdis 63,5hist xdis 63,6.1hist xdis 63,7.2;右帮水平移近量hist xdis 67.8,5hist xdis 67.8,6.1hist xdis 67.8,7.2step 60000save kgx.savecall kgs.txtres kgx.save;放顶煤2.2米def aaaxx=73.23supp xx 4.6 delendcommandxx=xx+1.55end_loopendaaadel range reg 157,4.5 157,8 180,8 180,4.5del range reg 72.5,8 72.5,10.2 180,10.2 180,8set pline 60,11.2 72,11.2 12set pline 63,7.75 67.8,7.75 12set pline 60,11.2 72,11.2 12set pline 63,4.5 63,7.75 8set pline 67.8,4.5 67.8,7.75 8;顶板下沉量hist ydis 63.3,7.75hist ydis 67.5,7.75hist ydis 65.4,7.75;底板鼓起量hist ydis 63.3,4.5hist ydis 67.5,4.5hist ydis 65.4,4.5;左帮水平移近量hist xdis 63,5hist xdis 63,6.1hist xdis 63,7.2;右帮水平移近量hist xdis 67.8,5hist xdis 67.8,6.1hist xdis 67.8,7.2set sup_delsolvesave kgs.savecall kgss.txt用的较多的是:1.changeCH ANGE keyword <keyword><ra nge. . . >Block, joint and cable element material characteristics are prescribed and changed with the CHANGE command. All blocks with centroids lying within the optional range (see Section 1.1.3) have block material characteristics changed. Likewise, all joints with contact coordinates lying within the optional range, or cables with nodes lying within the optional range, have material characteristics changed. If no range is specified, all blocks and joints will have characteristics changed according to the keywords given below. Extended zone models anduser-defined models are assigned using the extended ZONE command. User-defined joint models are defined using the JMODEL command. Local storage joint models are assigned using the JOINT command.The following keywords are used to change characteristics.1. Block Characteristicscons nConstitutive number n is assigned to designated deformable blocks(see Table 1.2).m at nMaterial property number n is assigned to designated rigid or deformableblocks. (All blocks initially default to mat = 1. The maximum value for n is 50.)Table 1.2 Constitutive models for deformable blockscons Model Description0 null material (The null model is used to model excavated material.The stresses within the null block are automatically set to zero.)1 linearly elastic, isotropic (default)3 elastic/plastic, Mohr-Coulomb failure (This model should be used withcaution since accurate solution to plasticity problems requires thatthe triangular zoning have a gridpoint at the centroid of each block.GENERATE quad zoning should be used whenever possible to improveplasticity analyses. However, no significant errors have been noted inproblems for which the above criterion has not been met.)6 elastic/plastic, Drucker-Prager failure (The same caution discussed abovefor the Mohr-Coulomb model (cons = 3) also applies for the Drucker-Prager model.)CH ANGENOTES:1. Constitutive models may also be assigned to regions of zones within blocks withthe ZONE command.2. All block constitutive models are described in Section 2 in Theory and Background.2. Joint Characteristicsjc ons nConstitutive number n is assigned to designated contacts (see Table1.3).jm at nMaterial property number n is assigned to designated contacts. (All contacts initially default to jmat = 1. The maximum value for n is 50.)Table 1.3 Joint constitutive modelsjcons Model Description1 point contact elastic/plastic with Coulomb slip failure (units are [force/displacement] for contact stiffnesses, and [force] for cohesion and tension)2 joint area* contact elastic/plastic with Coulomb slip failure (units are [stress] for cohesion and tension, and[stress/displacement] for joint stiffnesses) (default)3 continuously yielding joint model (see Section 3 in Theory and Background for a detailed explanation)5 same as jcons = 2, except that the internal fracture flag is set for each joint segment when joint shear or tensile strength is exceeded. If the fracture flag is set, residual values for friction, cohesion and tension are used in all subsequent calculations.7 optional Barton-Bandis (BB) joint model. See Section 3 in Special Features for details.* NOTE: The minimum joint area is limited to twice the rounding length, so thatit is not necessary to specify point contact properties if both point contacts andarea contacts occur between blocks.3. Cable Characteristicscab le matg <mats>Material property number matg is assigned to designated cable nodes.Material property number mats is assigned to designated cable elements.Midpoint of element must lie within the range to be changed.(All cable elements and nodes default to mat=1. The maximum valuefor matg or mats is 50.)4. Domain Characteristicsdm at nMaterial property number n is assigned to designated domains. (Alldomains initially default to dmat=1. The maximum value for n is 50.)CAB LE x1 y1 x2 y2 npoint mats matg <preten><keyword>delete <range>Execution of this command creates reinforcing elements that explicitly model theshear behavior of a grout annulus. These elements are generated between endpoints(x1,y1) and (x2,y2).NOTE: For excavation problems, if point (x1,y1) is inside the excavation periphery,the first nodal point will be on the excavation periphery. The point (x2,y2) shouldalways be located in the rock mass. This reinforcing logic can only be used withdeformable blocks, and the CABLE command must be invoked after the GENERATEcommand. Reinforcement properties and grout properties are specified and storedusing the PROPERTY command. The following parameters are also required.matg material property number for grout. (The PROPERTY commandshould be used to specify cb kbond and cb sbond for grout.)NOTE: The units for cb kbond and cb sbond are [force / cable length/ displacement] and [force / cable length], respectively.mats material property number for reinforcing (i.e., steel). (The PROPERTYcommand should be used to specify cb ymod, cb yield, cb ycomp,cb fstrain, and cb density for the material.)npoint number of lumped mass nodal points (npoint ≥2)Pre-tensioning of the reinforcing can be specified by providing a value (in unitsof force) for the optional parameter preten. Always use the PRINT property cablecommand to check property assignment.The following keywords may be used.connect The optional connect keyword will shift the cable node end closest to a structural (beam) element node to coincide with the beam node.The structural element nodes must be created before the cable nodes.The connection is not allowed to fail. The common, or shared, nodewill not appear in the list of cable nodes—it does appear in the list of structural element nodes. The cable element list will show the commonnode ID number at one end of the cable element and the normalcable node at the other. Beam and cable nodes are entirely separateentities. This is because being that the beam nodes are in contactwith the surface of a block, whereas the cable nodes are connected tothe model in the interior of a zone. All common nodes are thereforestructural element nodes, rather than cable nodes. As a consequence,it is possible for the cable element information display to indicateidentical node IDs at either end when the connect keyword is used.UDEC Version 4.0delete The optional delete keyword will remove all cable elements from the specified range. If no range is given, all cable elements will bedeleted.extend The extend keyword causes the end of the current cable to connect to the end node of an existing cable element.The cable reinforcing model is described in Section 1 in Special Features.The material properties for segments of an existing cable element can be changed with the CHANGE cable command.Plot ca ble <keyword><n>plots location of cable elements. For most of the keywords (shown below),the identification number n for the cable (defined by the CABLEcommand) may be used to plot only the cable elements associatedwith the number. The corresponding value of n for the cable maybe found by using the command PLOT cable number. The following keywords may be used.af ail cable element axial failureax ial <n> axial forcePL OT ca ble ele mentele ment cable element numbersfa il shows failure modes in color.gf ail cable node grout failureno de cable node numbersnu mber the identification number of the cable group used, for example, for line plotssd isp nodal point displacement vectorssh ear <n> relative shear force between cable node and host material st rain <n> axial strain in cable elementssvel nodal point velocity vectorsxd isp <n> x-displacement of cable nodesxvel <n> x-velocity of cable nodesyd isp <n> y-displacement of cable nodesyvel <n> y-velocity of cable nodesNOTE:(1) The line plot switch can be used to plot a cable variable as a line plot (see plot category 4).(2) The sense of the cable variable plots can be reversed by giving the switch yrev after the cable variable keyword.Plot stru ct <keyword><n>plots location of structural (beam) elements and associated variables. (See the STRUCT command and Section 1 in Special Features.) For certain keywords (shown below), the identification number, n, for the beam (defined by the STRUCT command) may be used to plot only the structural elements associated with that number. The corresponding value of n for the beam may be found by using the command PLOT struct number. The following keywords may be used.af ail structural element axial failureav el <n> angular velocity of elementax ial <n> axial forceif ail interface failurein ormal <n> normal force at structural interfaceint erface structural interface locationis hear <n> shear force at structural interfacemo ment <n> momentnu mber ID number of the element group used, for example, forline plotssd isp nodal point displacement vectorssh ear <n> shear forcesv el velocity vectors of structural nodesth ick structural element is plotted with thickness rather thana single linexd isp <n> x-displacement of structural nodesxv el <n> x-velocity of structural nodesyd isp <n> y-displacement of structural nodesyv el <n> y-velocity of structural nodesNOTE:(1) The line plot switch can be used to plot a structural element variableas a line plot (see plot category 4).(2) The sense of structural element variable plots can be reversed bygiving the switch yrev after the structural element keyword.Plot su pport structural support elementsPrint ca ble information on cable-reinforcing elements. The output is divided into cable-element information and cable-node information. (SeeSection 1 in Special Features for further descriptions.) The columnheadings for cable elements are:(1) cable element address(2) cable element segment identification number(3) cable node identification number — first node(4) cable node identification number — second node(5) material number for cable element(6) cross-sectional area of cable element(7) axial force in cable element (tension is negative)(8) length of cable element(9) axial strain in cable element (tension is negative)(10) failure indicator for element:0 elastic1 at yield in tension2 elastic, tension yield in past3 at yield in compression4 elastic, compression yield in past5 failure in extensional strainThe column headings for cable nodes are:(1) cable node address(2) cable node identification number(3) x-coordinate of cable node(4) y-coordinate of cable node(5) x-velocity of cable node(6) y-velocity of cable node(7) x-displacement of cable node(8) y-displacement of cable node(9) material number for cable node(10) x-force applied to node(11) y-force applied to node(12) shear force in grout(13) node fixed in x if = 1, else 0(14) node fixed in y if = 1, else 0(15) failure indicator for node:0 elastic1 yield2 elastic, yield in past(16) address of zone in which node is currently locatedPlot stru ct keywordinformation on structural elements. The following keywords are available.el ement <keyword>information on structural element segments. The following optional keywords are available to print information separately.d isp displacements. The column headings are:(1) structural element address(2) x-coordinate of midpoint(3) y-coordinate of midpoint(4) x-displacement of element(5) y-displacement of elementf orce forces. The column headings are:(1) structural element ID number(2) axial force (compression positive)。
