UDEC60常用命令集锦
- 格式:doc
- 大小:583.50 KB
- 文档页数:9
UDEC60常用命令集锦rest falll.savdelete range -2,2 -2,2reset disp reset hist hist un bal hist ydis 0,2 step 2000plot hold block stress disp ;位移归0edgeExample 3.2 EJfta rtrundit^ latglh on crack getteraiionro D .2;ro D . 1bi 0.0 o f io i0f in lo.o cr 0,3.09.7f10]f the rounding length ia reduced (o O. I. then the crack will be located where sfyecified.I lie SET edgeThisconimana, the u^er 匚日n set a small rnunding length for solution accuracy bur avoid blocks with small edg亡lengths and, consequently, adversel} high aspect ratios (see below). For exajnpk. ifcommand allows the user to Ueiine a minimum blxKvd*亡kngch nianually. WithLhe SET edge 0.4 and ROUND 0.1 arc $pe(;ifk也then block 、tnal吐山an G4 VrilL ilol be cruaicMl. and the rounding length for blocks ill be OJ. These coinnunds be gi\ en beLor^tlie BLOCK coinmaiid.Joints jregionjset angle 135 trace 1 gap 2 spacing 1.732 origin 0.259,0.966 range jreg 1 ;jregion 2 jointsjset angle jset angle 0 trace 10 trace 1gap 2 spacing 1?732 origin Q a0 range jreg 2gap 2 spacing 1.732 origin 1?5.0.866 range jreg 2jset angle 60 trace 1 gap 2 spacing 1?732 origin 1.0 range jreg 2jset angle 60 trace 1 gap 2 spacing 1?732 origin. 1.1.732 range jreg 2Joints that are joined are assigned high values for normal and shear stiffness. It is tempting to give VCQ high values for stiffnesses to prevent movement along a fictitious joint. 1 lowcvcr, the timestep calculation in UDEC is based upon stiffnesses; the response (and solution convergence) will be very slow if very high siifTncsses are spexified. Tlic lowest MifTiiess consisieni with small joint deformaiion is used when joints are joined. The rule-of-thumb is that joint stifTnesses t I n and k s. should be set to a factor times the equivalent stiffness of the stiffest neighboring zone. The equation to calculate fictitious joint normal stiffness related to the equivalent stiffness (expressed in stress-per-distance units) of a zone in the normal direction is of the form k n = f actor x max(3.1) Figure 3.14 Zone dimension used ui stiffness calculationSystems of Units2.S Systems of L'nilsUDEC accepts any con listen f of engineering units, Examples of consis:enr sets of units forh阳it parameters are shown in Table 2.5. The user should be very careful when converting from one stem of uniK to unoiher. An excellent reference on the subject or units and conversion between the Imperiikl and SI s\ stems can be iimnd in the Juunuti。
UDEC实例翻译与命令解析翻译:珠穆朗玛1 地震诱发地层坍塌Seismic-Induced Groundfall1.1 问题描述本例展示使用UDEC模拟分析地震诱发地层坍塌的一类的问题,模型见图1.1,该模型基于加拿大安大略省萨德伯里市鹰桥公司弗雷则矿34-1-554切割断面的一个剖面图的结构和尺寸. 用二维平面应变模型代表垂直于超采轴向方向的平面效应,超采面高5m,宽10m.假定两个连续节理交叉平面分析:一个角度为45度,另一个为-9度,两者节理间距均为5m,为了演示的目的,一个近似垂直的“虚拟节理”也被添加到块体内开挖面顶部以增强不稳定性。
围岩参数来自试验室平均测试数值,假定岩石块体参数如下:假定块体仅具有弹性行为,节理假定符合库伦滑动准则,选择典型的教课书数值作为节理参数,如下:初始应力状态按各向同性估计为24Mpa(假定垂直荷载由覆盖深度大约800m的岩层产生)。
1.2 UDEC分析UDEC模拟顺序分三个阶段,首先,模型在初始应力状态下进行无超采固结.其次,进行开挖并且模型循环至平衡状态.本阶段超采面周围的应力分布见图 1.2.超采正上方和下方的块体滑动后稳定.在第三阶段.估计了两个不同的峰值速度的地震事件.对所有地震模拟,在问题域的外周边界引入粘滞边界用以消除波的反射.从而模拟有限的岩体,地震事件用施加到模型顶部y方向的正弦应力波表现.应力波被叠加到已存在的初始地应力上.在第一个模拟中,施加1.25Mpa的峰值应力,应当注意的是,由于粘滞边界条件实际是在模型顶部, 施加的有效影响应力应该是1.25 MPa/2, or 0.625 MPa.0.02秒后的开挖面拱顶的应力分布见图1.3,两点的位移被监测,1点位于开挖面的左角,点2位于拱顶块体的右角, 图1.4的位移时间曲线显示两点本质上是弹性反应.本例关心的问题是在模型顶部施加的速度和计算速度的对比,下面的公式可以用以估计施加的波速.使用这个方程,施加的最大波速大概是0.04m/sec,图1.5显示的峰值波速小于0.06m/sec. 估计的波速和监测波速的不同在于使用的围岩模量.而是没有考虑节理变形的相等变形模量.在第二个案例中,施加应力波峰值12.5 Mpa(有效应力6.25Mpa).0.02秒后的开挖拱顶应力分布见图1.6.该图显示出拱顶岩体不受力,表面该块体已经松散并正在下落.对于关心的问题,后来三个时间的几何体和应力分布见图1.8至图1.10.在问题的顶部预测的波速(从上面的方程)是0.4m/sec.从模型中计算的波速见图1.11,再次,由于使用的是原岩弹性模量而不是岩体的变形模量导致预测和监测的波速之间的差异.1.3节包含了该模型的数据列表,该列表包含了一个FISH函数(show)被用来创建坍塌的动画文件,每隔0.02秒俘获一个显示的图片.通过改变FISH参数time_int可以改变动画帧的间隔.视图的总数也可以通过改变snap_shot的数值进行改变.为了显示80帧的显示图片而创建的该电影文件需要大概13MB的硬盘空间.1.3 数据文件列表Example 1.1 SEISMIC.DATtitleSEISMIC INDUCED ROOF COLLAPSE 地震诱发拱顶坍塌;round 0.01; define original boundary of modeled region 定义模型区域的原始边界block -25,-20 -25,20 25,20 25,-20; generate joint pattern over entire original region 在整个原始区域生成节理形态jregion id 1 -25,-25 -25,25 25,25 25,-25jset 45,0 200,0 0,0 5.0,0 (0,0) range jreg 1jset -9,0 200,0 0,0 5.0,0 (0,0) range jreg 1; put in joints needed for the later excavation 为了后面开挖而设置的节理crack -5.01,-2.51 5.01,-2.51crack -5.01, 2.51 5.01, 2.51crack -5,-2.5 -5,2.5crack 5,-2.5 5,2.5crack 2.25,2.5 1.93,5.0; generate fdef zones and assign joint properties (mat=1 & jmat=1;default) 生成单元和设置节理参数generate edge 9.0 range -30,30 -30,30prop mat=1 d=0.00300 k=39060 g=31780prop jmat=1 jkn=20000 jks=20000prop jmat=1 jf=30.0; apply boundary conditions and initial conditions to 在地应力下施加边界条件和初始条件; consolidate model under field stressesbound stress=-24.0, 0.0, -24.0 ygrad=-.03 0 -.03insitu stress=-24.0, 0.0, -24.0 ygrad=-.03 0 -.03bound yvel 0.0 range -26,26 -21,-19grav 0.0 -10.0; track the x-displacement, and y-displacement over time 追踪位移hist solvehist xdis=0,7 ydis=0,7 type 1solve rat 1e-5; save consolidated statesave seismic1.sav; make excavationdelete range -5,5 -2.5,2.5solve rat 1e-5; save excavated statesave seismic2.sav;rest seismic2.sav; apply seismic load from top (peak velocity=0.04 m/sec);; set up nonreflecting boundarybound mat=1bound xvisc range -26 -23 -21 21bound xvisc range 23 26 -21 21bound xvisc yvisc range -26 26 -21 -19bound xvisc yvisc range -26 26 19 21; apply sinusoidal stress wavebound stress 0 0 -1.25 yhist=cos(100.0,0.0195) range -26 26 19 21;reset time hist disp rothist ydis (-4.48,2.57)hist ydis (0,2.57) yvel (0,2.57) yvel (4,2.57) yvel(-4.48,2.57)hist yvel (0,20) yvel (25,10) yvel (25,-10) yvel (0,-20)hist yvel (-25,-10) yvel (-25,10)hist sxx (25,10) sxx (25,-10) sxx (-25,-10) sxx (-25,10)hist syy (0,20);damp 0.1 1.0 mass; 0.02 sec.cyc time 0.02save seismic3.sav;rest seismic2.sav; apply seismic load from top (peak velocity=0.4 m/sec); set up nonreflecting boundarybound mat=1bound xvisc range -26 -23 -21 21bound xvisc range 23 26 -21 21bound xvisc yvisc range -26 26 -21 -19bound xvisc yvisc range -26 26 19 21; apply sinusoidal stress wavebound stress 0 0 -12.5 yhist=cos(100.0,0.0195) range -26 26 19 21reset time hist disphist ydis (-4.48,2.57)hist ydis (0,2.57) yvel (0,2.57) yvel (4,2.57) yvel(-4.48,2.57)hist yvel (0,20) yvel (25,10) yvel (25,-10) yvel (0,-20)hist yvel (-25,-10) yvel (-25,10)hist sxx (25,10) sxx (25,-10) sxx (-25,-10) sxx (-25,10)hist syy (0,20);damp 0.1 1.0 masssave seismov.sav;; 0.02 sec. —————————————————————————————————————cyc time 0.02save seismic4.sav; 0.25 sec.cyc time 0.23save seismic5.sav; 0.50 sec.cyc time 0.25save seismic6.sav; 0.75 seccyc time 0.25save seismic7.sav;rest seismov.sav; make a movie of the groundfall;wind -12 12 -12 12set ovtol 0.05plot block vel max 2.0 blue stress max 50movie onmovie file = seismic.dcxmovie step 1000step 400003 隧道支护荷载Tunnel Support Loading3.1 问题陈述本例模拟展示了UDEC在检查衬砌隧道方面的应用,着重强调了荷载在混凝土衬砌中的发展,本例也解释了模拟连续建造操作中独立阶段的模拟程序.隧道系统的理想几何体见图3.1.系统包含在海床下大约70m(中线)深度,中线间距12m 的两个隧道, 初始水位在隧道中线上方110m处.服务隧道直径5.24m,衬砌厚度37cm.主隧道直径8.22m,衬砌厚度46cm.服务隧道先于主隧道开挖和衬砌.随后设置主隧道衬砌,水位上升增加到100m.—————————————————————————————————————施工顺序是:(1)开挖服务隧道excavation of the service tunnel;(2)衬砌服务隧道lining of the service tunnel;(3)开挖主隧道excavation of the main tunnel;(4)衬砌主隧道lining of the main tunnel; and(5)升高水位raising of the water level.分析的目的是评价每个施工阶段服务隧道和主隧道支护状况.本例的材料参数见下:岩体——开挖隧道的围岩参数为:弹性模量 elastic modulus 0.89 GPa泊松比Poisson’s ratio 0.35单轴抗压强度uniaxial compressive strength 3.5 MPa粘聚力 cohesion 1 MPa密度 density 1340 kg/m3混凝土衬砌——弹性模量为 24 GPa ,泊松比为0.19. 假定衬砌为线弹性材料。
服务器指令大全服务器指令大全本文档为服务器指令大全,包含了一系列常用的服务器指令及其详细使用方法。
以下是各个章节的细化内容。
第1章:服务器管理指令1.1 登录服务器指令登录服务器指令用于远程登录服务器,具体指令包括:- ssh:用于以SSH协议登录远程服务器。
- telnet:用于以Telnet协议登录远程服务器。
1.2 用户管理指令用户管理指令用于管理服务器上的用户,包括:- useradd:用于添加新用户。
- userdel:用于删除用户。
- passwd:用于设置用户密码。
1.3 进程管理指令进程管理指令用于管理服务器上的进程,常用指令有:- ps:用于查看当前运行的进程列表。
- kill:用于终止指定进程。
1.4 系统信息指令系统信息指令用于查看服务器的基本信息,包括:- uname:用于查看系统内核和硬件信息。
- df:用于查看磁盘空间使用情况。
- free:用于查看内存使用情况。
第2章:文件管理指令2.1 文件操作指令文件操作指令用于对文件进行操作,包括:- ls:用于列出目录内容。
- cd:用于切换当前工作目录。
- cp:用于复制文件或目录。
- mv:用于移动文件或目录。
- rm:用于删除文件或目录。
2.2 文件权限指令文件权限指令用于设置文件的权限和所有权,包括:- chmod:用于修改文件或目录的权限。
- chown:用于修改文件或目录的所有者。
- chgrp:用于修改文件或目录的所属组。
2.3 文件搜索指令文件搜索指令用于在服务器中搜索文件,包括:- find:用于按照指定条件搜索文件。
- grep:用于在文件中搜索指定内容。
第3章:网络管理指令3.1 网络配置指令网络配置指令用于配置服务器的网络设置,包括:- ifconfig:用于查看和配置网络接口信息。
- route:用于查看和配置路由表信息。
3.2 网络状态指令网络状态指令用于查看服务器的网络连接状态,包括:- netstat:用于显示网络状态信息。
目录实验一以太网交换机基本配置 (1)实验二以太网端口配置实验 (7)实验三利用TFTP管理交换机配置 (13)实验四虚拟局域网VLAN (16)实验五生成树配置 (25)实验六802.1x和AAA配置 (38)实验七路由器基本配置 (445)实验八PPP配置 (51)实验九FR配置 (56)实验十静态路由协议配置 (64)实验十一RIP协议配置 (68)实验十二OSPF协议配置 (74)实验十三访问控制列表配置 (88)实验十四地址转换配置 (95)实验十五DHCP配置 (101)实验十六升级路由器或交换机的操作系统 (116)实验一以太网交换机基本配置一、交换机常用命令配置模式1 业务描述(1)Quidway系列产品的系统命令采用分级保护方式,命令被划分为参观级、监控级、配置级、管理级4个级别,简介如下:✧参观级:网络诊断工具命令(ping、tracert)、从本设备出发访问外部设备的命令(包括:Telnet客户端、RLogin)等,该级别命令不允许进行配置文件保存的操作。
✧监控级:用于系统维护、业务故障诊断等,包括display、debugging命令,该级别命令不允许进行配置文件保存的操作。
✧配置级:业务配置命令,包括路由、各个网络层次的命令,这些用于向用户提供直接网络服务。
✧管理级:关系到系统基本运行,系统支撑模块的命令,这些命令对业务提供支撑作用,包括文件系统、FTP、TFTP、XModem下载、配置文件切换命令、电源控制命令、备板控制命令、用户管理命令、命令级别设置命令、系统内部参数设置命令等。
(2)命令视图:系统将命令行接口划分为若干个命令视图,系统的所有命令都注册在某个(或某些)命令视图下,只有在相应的视图下才能执行该视图下的命令:各命令视图的功能特性、进入各视图的命令等的细则:◆命令视图功能特性列表2 配置参考命令(1)命令行在线帮助在任一命令视图下,键入“?”获取该命令视图下所有的命令及其简单描述。
;newtitle moni ;定义名称moni块体都有“圆角”,其目的在于避免块体悬挂在有棱角的节点上。
由于块体悬挂引起应力集中。
然而,圆角值存在与模型有关的上限值。
对于变形块体,最大圆角长度应当不超过块体平均棱长的1%。
圆角长度可以如下命令加以改变:round dround=0.025;d 是圆角距离(缺省值是d=0.5)。
模型中的所有圆角长度都是相同的。
For example, if the commands SET edge 0.4and ROUND 0.1 are specified, then block edge lengths smaller than 0.4 will not be created, and the rounding length for blocks will be 0.1. These commands must be given before the BLOCK command.set ovtol=1;块体与块体之间相互嵌入量最大值为1米bl 0,0 0,100 300,100 300,0 ;定义范围,四点坐标,顺时针方向cr 0,15 300,15;Crack 命令用于产生块体中单一直线特征的裂缝,裂缝由端点坐标(x1,y1)和(x2,y2)所确定。
jreg id 1 0,0 0,15 300,15 300,0 delete ;定义节理,命名1,定义范围,四点坐标,delete为常规语言jset 90,0 5,0 5,0 5,0 0,0 rang jreg 1;Jset 命令则是自动节理组生成器。
根据所给定的特征参数(即倾角、迹长、岩桥长度、间距和空间位置)产生一组裂缝。
定义节理角度90°;节理的长度;节理的距离,即纵向间隔;横向距离,即横向隔5一个;起始点坐标为0,0jset 90,0 5,0 5,0 5,0 2.5,5 rang jreg 1 ;定义节理的另外一项gen quad 10 6 range 0 300 0 100 ;定义块体最大变形,若没有此语句,刚所有块体均为刚性块体;automatic generation of diagonally opposed triangular zones to improve plastic flow calculation.(对角三角形区域的自动生成改善塑性流动计算) Parameters xw and yw are zone widths in the x-and y-directions zone model mo range 0 300 0 100 ;定义摩尔库伦模型的范围,X的范围,Y的范围change jcons=5 range 0 300 0 100 ;定义节理的某个属性change mat=1 range reg 0,0 0,15 300,15 300,0 ;定义物质1的范围,即赋予这个范围为物质1,只赋名,没有定义属性change mat=2 range reg 0,15 0,20 300,20 300,15change jmat=1 range reg 0,0 0,15 300,15 300,0 ;定义节理1的范围set jmatdf 2 ;定义节理属性,使符合摩尔库伦准则save ch.sav ;保存文件cal prop-25.txt ;调出文件res ch.sav ;调出前面保存的文件接着计算;参数设置prop mat=1 dens=2720 ;物质的密度zone k=12.12e9 g=10.26e9 fric=33 coh=5.77e6 ten=3.6e6 range mat=1;定义物质1的体积模量、剪切模量、内摩擦角、粘聚力、抗拉强度prop jmat=1 jkn=11e8 jks=12e7 jfric=12 jcoh=0 jten=0;定义节理的属性法向刚度、切向刚度、内摩擦角、粘聚力、抗拉强度;边界条件boun stress 0 0 0 range 0,300 99.9 100.1 ;定义边界条件,三向应力为0,上边界,范围,纵向波动范围boun stress 0 0 -8.6095397e6 range ;定义边界应力条件,竖直方向的应力为,负号表示方向向下insitu str -7.5e6 0 -15e6 szz -7.5e6 ygrad 1.1018e4 0 2.203607e4 ;定义应力属性,sxx方向即水平应力,sxy方向即剪切应力,syy方向即竖直应力,szz方向的应力,ygrad应力梯度set gravity 0,-10 ;重力加速度bound xvel=0 range -0.1 0.1 0 100 ;X位移边界,位移波动范围,范围bound xvel=0 range 299.9 300.1 0 100 ;X位移边界,同上bound yvel=0 range -0.1 300.1 -0.1,0.1 ;Y位移边界hist solve_ratio type 1 ;设置不平衡力的精度,普适solve ;计算save 25-1.sav ;保存cal 25-2.txt ;调出文件res 25-1.sav ;调出前面保存的文件del 80,110 15,20;永久删除这个范围,即开挖,不可再充填,X的范围,Y的范围model null range ;变性删除这个范围,即之后可以充填物质,X的范围,Y的范围model m range ;充填这个范围,solve ;计算save 25-301.sav ;再保存,周而复始set pline 25,20 300,20 10 ;定义测线,起始点坐标,测线分段, 前两组数字表示坐标最后一个数表示布置了多少测点reset histset log on ;定义hist,保存在一个位置,在程序中输入hist1即可调出hist xdisp 10,20 ;测该点的X方向位移hist sxx 10,20 ;测该点X方向的应力hist syy 10,20 ;测该点Y方向的应力reset dispreset hist ;位移清零,历史记录清零his unbalsolve ratio 1e-5 ;记录不平衡力,计算直到精度达到1e-5ch cons=3 range 77.5,80 0,3ch mat=7 range 77.5,80 0,3 ;应变软化模型,这个是做充填用的cable (73.4,2.91) (71.15,4.99) 10 12 380e-6 11 ;锚杆,起始点坐标,后面四个数值可固定prop mat 12 cb_dens 7800 cb_ycomp 430e6 cb_yield 160e3 cb_ymod1.3e11prop mat 11 cb_kbond 6.3e9 cb_sbond 6e5;定义锚杆的属性,密度、抗拉强度、屈服强度density for block materical 1 is zero cannot cycle----块体密度未设range 范围;density 密度;volume 体积;regedit 调出注册表;pl bound pl 显示塑性区;pl bound dis 显示位移矢量图continue pa以后继续算number 可以不关udec,继续使用udecres xx.save (先调用文件)set pline x1,y1 x2,y2 n (n--观测线分的段数)set log on ;固定的set log 文件名print+pline+n+ syy sxx or xdis (观测线名称任意数字)set log off1) pl bl hist 1,2,... 观测点位移2) pl bl nu(mber) 显示块体标号3) pl bl cab red stru red cab--锚杆stru--梁4) pl bl sxx/syy/szz 显示应力图5) pl bl dis 位移矢量图6) set pl dxf 256 把UDEC 图引入CAD 换行copy 文件名自定.dxf块体力学参数k--体积模量;g--剪切模量;fric--摩擦角;coh--粘结力(内聚力);ten--抗拉强度;d--体积力接触面力学参数jkn--法向刚度;jks--切向强度;jc--粘结力;jf--摩擦角;jt--抗拉强度1) D--质量密度F--摩擦角B--体积模量; C--内聚力(粘结力) S--剪切模量2) 应力--正号代表张力,负号代表压力应变--正的应变代表伸长,负的应变代表压缩重力--正号的重力物质往下拉,负号的重力将物质往上提;内存赋值udec mudec 14print memhist yvel (20,20) type 1; type 是在屏幕上以指定的间隔显示其值step 一般典型问题要2000~4000次循环。
;newtitle moni ;定义名称moni块体都有“圆角”,其目的在于避免块体悬挂在有棱角的节点上。
由于块体悬挂引起应力集中。
然而,圆角值存在与模型有关的上限值。
对于变形块体,最大圆角长度应当不超过块体平均棱长的1%。
圆角长度可以如下命令加以改变:round dround=0.025;d 是圆角距离(缺省值是d=0.5)。
模型中的所有圆角长度都是相同的。
For example, if the commands SET edge 0.4and ROUND 0.1 are specified, then block edge lengths smaller than 0.4 will not be created, and the rounding length for blocks will be 0.1. These commands must be given before the BLOCK command.set ovtol=1;块体与块体之间相互嵌入量最大值为1米bl 0,0 0,100 300,100 300,0 ;定义范围,四点坐标,顺时针方向cr 0,15 300,15;Crack 命令用于产生块体中单一直线特征的裂缝,裂缝由端点坐标(x1,y1)和(x2,y2)所确定。
jreg id 1 0,0 0,15 300,15 300,0 delete ;定义节理,命名1,定义范围,四点坐标,delete为常规语言jset 90,0 5,0 5,0 5,0 0,0 rang jreg 1;Jset 命令则是自动节理组生成器。
根据所给定的特征参数(即倾角、迹长、岩桥长度、间距和空间位置)产生一组裂缝。
定义节理角度90°;节理的长度;节理的距离,即纵向间隔;横向距离,即横向隔5一个;起始点坐标为0,0jset 90,0 5,0 5,0 5,0 2.5,5 rang jreg 1 ;定义节理的另外一项gen quad 10 6 range 0 300 0 100 ;定义块体最大变形,若没有此语句,刚所有块体均为刚性块体;automatic generation of diagonally opposed triangular zones to improve plastic flow calculation.(对角三角形区域的自动生成改善塑性流动计算) Parameters xw and yw are zone widths in the x-and y-directions zone model mo range 0 300 0 100 ;定义摩尔库伦模型的范围,X的范围,Y的范围change jcons=5 range 0 300 0 100 ;定义节理的某个属性change mat=1 range reg 0,0 0,15 300,15 300,0 ;定义物质1的范围,即赋予这个范围为物质1,只赋名,没有定义属性change mat=2 range reg 0,15 0,20 300,20 300,15change jmat=1 range reg 0,0 0,15 300,15 300,0 ;定义节理1的范围set jmatdf 2 ;定义节理属性,使符合摩尔库伦准则save ch.sav ;保存文件cal prop-25.txt ;调出文件res ch.sav ;调出前面保存的文件接着计算;参数设置prop mat=1 dens=2720 ;物质的密度zone k=12.12e9 g=10.26e9 fric=33 coh=5.77e6 ten=3.6e6 range mat=1;定义物质1的体积模量、剪切模量、内摩擦角、粘聚力、抗拉强度prop jmat=1 jkn=11e8 jks=12e7 jfric=12 jcoh=0 jten=0;定义节理的属性法向刚度、切向刚度、内摩擦角、粘聚力、抗拉强度;边界条件boun stress 0 0 0 range 0,300 99.9 100.1 ;定义边界条件,三向应力为0,上边界,范围,纵向波动范围boun stress 0 0 -8.6095397e6 range ;定义边界应力条件,竖直方向的应力为,负号表示方向向下insitu str -7.5e6 0 -15e6 szz -7.5e6 ygrad 1.1018e4 0 2.203607e4 ;定义应力属性,sxx方向即水平应力,sxy方向即剪切应力,syy方向即竖直应力,szz方向的应力,ygrad应力梯度set gravity 0,-10 ;重力加速度bound xvel=0 range -0.1 0.1 0 100 ;X位移边界,位移波动范围,范围bound xvel=0 range 299.9 300.1 0 100 ;X位移边界,同上bound yvel=0 range -0.1 300.1 -0.1,0.1 ;Y位移边界hist solve_ratio type 1 ;设置不平衡力的精度,普适solve ;计算save 25-1.sav ;保存cal 25-2.txt ;调出文件res 25-1.sav ;调出前面保存的文件del 80,110 15,20;永久删除这个范围,即开挖,不可再充填,X的范围,Y的范围model null range ;变性删除这个范围,即之后可以充填物质,X的范围,Y的范围model m range ;充填这个范围,solve ;计算save 25-301.sav ;再保存,周而复始set pline 25,20 300,20 10 ;定义测线,起始点坐标,测线分段, 前两组数字表示坐标最后一个数表示布置了多少测点reset histset log on ;定义hist,保存在一个位置,在程序中输入hist1即可调出hist xdisp 10,20 ;测该点的X方向位移hist sxx 10,20 ;测该点X方向的应力hist syy 10,20 ;测该点Y方向的应力reset dispreset hist ;位移清零,历史记录清零his unbalsolve ratio 1e-5 ;记录不平衡力,计算直到精度达到1e-5ch cons=3 range 77.5,80 0,3ch mat=7 range 77.5,80 0,3 ;应变软化模型,这个是做充填用的cable (73.4,2.91) (71.15,4.99) 10 12 380e-6 11 ;锚杆,起始点坐标,后面四个数值可固定prop mat 12 cb_dens 7800 cb_ycomp 430e6 cb_yield 160e3 cb_ymod1.3e11prop mat 11 cb_kbond 6.3e9 cb_sbond 6e5;定义锚杆的属性,密度、抗拉强度、屈服强度density for block materical 1 is zero cannot cycle----块体密度未设range 范围;density 密度;volume 体积;regedit 调出注册表;pl bound pl 显示塑性区;pl bound dis 显示位移矢量图continue pa以后继续算number 可以不关udec,继续使用udecres xx.save (先调用文件)set pline x1,y1 x2,y2 n (n--观测线分的段数)set log on ;固定的set log 文件名print+pline+n+ syy sxx or xdis (观测线名称任意数字)set log off1) pl bl hist 1,2,... 观测点位移2) pl bl nu(mber) 显示块体标号3) pl bl cab red stru red cab--锚杆stru--梁4) pl bl sxx/syy/szz 显示应力图5) pl bl dis 位移矢量图6) set pl dxf 256 把UDEC 图引入CAD 换行copy 文件名自定.dxf块体力学参数k--体积模量;g--剪切模量;fric--摩擦角;coh--粘结力(内聚力);ten--抗拉强度;d--体积力接触面力学参数jkn--法向刚度;jks--切向强度;jc--粘结力;jf--摩擦角;jt--抗拉强度1) D--质量密度F--摩擦角B--体积模量; C--内聚力(粘结力) S--剪切模量2) 应力--正号代表张力,负号代表压力应变--正的应变代表伸长,负的应变代表压缩重力--正号的重力物质往下拉,负号的重力将物质往上提;内存赋值udec mudec 14print memhist yvel (20,20) type 1; type 是在屏幕上以指定的间隔显示其值step 一般典型问题要2000~4000次循环。
《Udec 命令总结》精华版1. 安装、打开、保存(1) 安装:①执行Udec 3.1→将Crack文件中的内容替换;②复制Udec.exe,粘贴为快捷方式→属性,目标,加入空格256【开始内存8M,将内存改为256】,可以复制快捷方式至桌面或硬盘。
(2) 打开:①直接在udec> 命令行输入;②写好程序,udec>命令输入call,然后将*.txt文件拖入命令行,执行(3) 保存:输入save d:\kaicai.sav, 调用命令rest,将kaicai.sav 拖入命令行【!!!文件名最好不用汉字,有时候不识别】注释:如果保存为save d:\111\kaicai.sav, 注意其中111文件必须提前建好,否则无法保存或者保存错误(4) 操作:Ctrl+Z选中图像可以放大,Ctrl+Z 双击复原,屏幕中会出现十字叉,按住鼠标左键不放,移动光标直到你满意的窗口为止;pause暂停,此时可以察看任何信息;continue继续调用下面程序段。
ESC可以随时进行停止,但不能继续;英文分号; 表示注释不运行命令。
2. 基本命令2.1 基本设置Udec>n【new刷新窗口,从新调用一个程序,修改后的*.txt文件必须输入n,重新运行文本文件】Udec>title 【或heading代表标题,后面紧跟标题的名称。
如:hang dao mo ni 或济宁三号井围岩变形破坏规律研究】Udec> round d 【“圆角”命令,Udec中所有的块体都有圆角,目的是为防止块体悬挂在有棱角的节点上,由于块体悬挂将产生应力集中。
d指块体与块体之间的圆角半径,默认值是0.5,其值要求小于模型中最小块体的最短那条边长的二分之一,最大圆角长度不能超过块体平均棱长的1%。
在block 命令前指定圆角长度。
如:round 0.05】圆角图1 10×10块体圆角Udec>set ovtol 0.5 【此命令是指层与层之间的嵌入厚度,当提示为“overlap too large”时就需要修改此值更大一些,可以显示设计的块体,plot overlap!嵌入太大的原因可能为块体强度太小】Udec>set log onprintset log off 【命令用于导出数据,可以将数据导出至硬盘f:\】命令:set log f:\yuanyan.logprint pline 1 syy 【记录测线1应力】print pline 2 ydisset log off;设置观测线set pline 1 x1,y1 x2,y2 n (n--观测线分的段数)set pline 25,20 300,20 10 ;定义测线,起始点坐标,测线分段2.2 图形划分;块体命令Udec> block x1,y1 x2,y2 x3,y3 x4,y4【建立模型框架,四个坐标角点必须按“顺时针”排列,也可以为五个坐标点。
服务器指令大全服务器指令大全:一、常用服务器指令:1.连接到远程服务器:ssh [用户名][服务器IP地质]2.查看当前路径:pwd3.切换目录:cd [目录路径]4.查看文件列表:ls5.查看文件内容:cat [文件名]6.创建新文件:touch [文件名]7.删除文件:rm [文件名]8.复制文件:cp [源文件] [目标文件] 9.移动文件:mv [源文件] [目标文件] 10.创建新目录:mkdir [目录名]11.删除目录:rmdir [目录名]12.查看进程列表:ps -ef13.终止进程:kill [进程ID]14.查看服务器负载:top15.重启服务器:reboot二、网络相关指令:1.查看网络接口信息:ifconfig2.查看网络连接状态:netstat3.查看路由表:route4.查看DNS配置:cat /etc/resolv.conf 5.测试网络连通性:ping [目标IP地质]三、安全相关指令:1.查看当前登录用户:who2.查看系统登录日志:cat /var/log/auth.log3.修改用户密码:passwd [用户名]4.查看用户权限:ls -l [文件/目录]5.修改文件权限:chmod [权限设置] [文件/目录] 6.查看防火墙状态:ufw status7.开启防火墙:ufw enable8.关闭防火墙:ufw disable四、数据库相关指令:1.登录数据库:mysql -u [用户名] -p2.创建数据库:CREATE DATABASE [数据库名];3.删除数据库:DROP DATABASE [数据库名];4.查看数据库列表:SHOW DATABASES;5.创建表格:CREATE TABLE [表名] ([列名1] [数据类型1], [列名2] [数据类型2],);6.插入数据:INSERT INTO [表名] ([列名1], [列名2],) VALUES ([值1], [值2],);7.查询数据:SELECT FROM [表名];8.更新数据:UPDATE [表名] SET [列名1]=[新值1], [列名2]=[新值2],WHERE [条件];9.删除数据:DELETE FROM [表名] WHERE [条件];附件:无附件。
rest falll.savdelete range -2,2 -2,2reset disp reset hist hist un bal hist ydis 0,2 step 2000plot hold block stress disp ;位移归0edgeExample 3.2 EJfta rtrundit^ latglh on crack getteraiionro D .2;ro D . 1bi 0.0 o f io i0f in lo.o cr 0,3.09.7f10]f the rounding length ia reduced (o O. I. then the crack will be located where sfyecified.I lie SET edgeThisconimana, the u^er 匚日n set a small rnunding length for solution accuracy bur avoid blocks with small edg亡lengths and, consequently, adversel} high aspect ratios (see below). For exajnpk. ifcommand allows the user to Ueiine a minimum blxKvd*亡kngch nianually. WithLhe SET edge 0.4 and ROUND 0.1 arc $pe(;ifk也then block 、tnal吐山an G4 VrilL ilol be cruaicMl. and the rounding length for blocks ill be OJ. These coinnunds be gi\ en beLor^tlie BLOCK coinmaiid.Joints jregionjset angle 135 trace 1 gap 2 spacing 1.732 origin 0.259,0.966 range jreg 1 ;jregion 2 jointsjset angle jset angle 0 trace 10 trace 1gap 2 spacing 1・732 origin Q a0 range jreg 2gap 2 spacing 1.732 origin 1・5.0.866 range jreg 2jset angle 60 trace 1 gap 2 spacing 1・732 origin 1.0 range jreg 2jset angle 60 trace 1 gap 2 spacing 1・732 origin. 1.1.732 range jreg 2Joints that are joined are assigned high values for normal and shear stiffness. It is tempting to give VCQ high values for stiffnesses to prevent movement along a fictitious joint. 1 lowcvcr, the timestep calculation in UDEC is based upon stiffnesses; the response (and solution convergence) will be very slow if very high siifTncsses are spexified. Tlic lowest MifTiiess consisieni with small joint deformaiion is used when joints are joined. The rule-of-thumb is that joint stifTnesses t I n and k s. should be set to a factor times the equivalent stiffness of the stiffest neighboring zone. The equation to calculate fictitious joint normal stiffness related to the equivalent stiffness (expressed in stress-per-distance units) of a zone in the normal direction is of the formk n = f actor x max(3.1) Figure 3.14 Zone dimension used ui stiffness calculationSystems of Units2.S Systems of L'nilsUDEC accepts any con listen f of engineering units, Examples of consis:enr sets of units forh阳it parameters are shown in Table 2.5. The user should be very careful when converting from one stem of uniK to unoiher. An excellent reference on the subject or units and conversion between the Imperiikl and SI s\ stems can be iimnd in the Juunuti </J Petndeuni lec /iimlti^y (DecCEiiber 1977). No conveniions are perfomicd in L DEC except for firictiun jnd diLatiun angles, M hich arc entered in degrees.Table 2.5 Sysleuta of units - mechanical parametersSI Imperial【工哩ih m m m cm ft inDen HQ kg/Hi'111? kg/in1IO6 kg/m1tlugt/ft3di bJ in"Foree N kN MN Mdyfies 呵IbfStress內kPa ME bar P点(Iraviij 7m/scc m/sec cnJs*ft/sec^in/sec^m/secwhere 1 bar = 10^ dynes! cm- = 1()^ N/ni- = 1(P Pa;1 atm = 1,013 bars = 14.7 psi = 2U6 lb(/ft2= 1.01325 x 10s Pa;1 slug = 1 Lbj —s- / ft = 14.59 kg;I snail = 1 Ib|- —in: andI gravity = 9.81 ni/s: = 9S1 cm/s2 =32.17 ft/s3.Histhist solve type 1solve rat 1e-5save tun l.savhist ydis 42 -67 sxx 42 -67hist ydis 42 -73 sxx 42 -73 reset disp jdispsolve rat 1e-5save tun 2.savPlotpl syy fill projected pl syy fill bfillJOS 1TTLE :UDEG \VwsiQfiFi却球 X27 Stress ertnfou rs using PLOT syy fill projectedpl ydisp fillpl bl id blackplot pen stress block显示应力pen作用是udec4.0窗口显示pl se ;看裂隙plsig1Lontours of major pnncipal shvss in the 工> -plant; (delunnable blocks only)sig2contours of minor principal stress in the xy-pl tine (deformable b kicks only)slip joints that are at limilin耳fricLioncontours of niuxiTTiuin prin^ipuJ slress in tbe XY-plane (defomiuble blocks only).Because 字亠皿□尸y鰹严.this is themiQ^^rincipiil stress I sm^lc = gig2,.J contours of ininimum pnriLiptd stress in the A v-plLint (deformublt 丿bkeks only).Because compressive stresses nrc negative, (his is the m ajar principal slrciis I smln = $ig1 )・ssi cmil门山*。
统信系统常用命令统信系统是一种常用的企业级通信系统,它提供了一系列的命令来管理和配置系统。
本文将介绍统信系统常用命令的功能和用法。
1. show version:该命令用于显示统信系统的版本信息,包括软件版本、硬件版本等。
通过该命令可以了解系统的当前状态和功能。
2. show status:该命令用于显示统信系统的运行状态,包括CPU利用率、内存利用率、磁盘空间等信息。
通过该命令可以监控系统的性能和资源使用情况。
3. show interface:该命令用于显示统信系统的接口信息,包括网卡接口、电话接口等。
通过该命令可以查看接口的配置和状态。
4. show configuration:该命令用于显示统信系统的配置信息,包括网络配置、电话配置等。
通过该命令可以查看系统的当前配置情况。
5. set configuration:该命令用于修改统信系统的配置信息,包括网络配置、电话配置等。
通过该命令可以对系统进行定制化设置。
6. restart service:该命令用于重启统信系统的服务,包括网络服务、电话服务等。
通过该命令可以重新启动某个服务以解决问题或应用新配置。
7. backup configuration:该命令用于备份统信系统的配置信息。
通过该命令可以将当前配置保存到备份文件中,以便在需要时恢复配置。
8. restore configuration:该命令用于恢复统信系统的配置信息。
通过该命令可以将备份的配置文件导入系统,以恢复之前的配置。
9. show log:该命令用于显示统信系统的日志信息,包括系统日志、事件日志等。
通过该命令可以查看系统的运行记录和故障信息。
10. clear log:该命令用于清除统信系统的日志信息。
通过该命令可以清空日志文件,以释放磁盘空间或保护隐私。
综上所述,统信系统常用命令提供了丰富的功能和灵活的操作方式,可以帮助管理员管理和配置系统。
通过这些命令,管理员可以监控系统状态、修改配置、备份恢复数据等,从而保证系统的稳定运行和高效管理。
Linux提供了大量的命令,利用它可以有效地完成大量的工作,如磁盘操作、文件存取、目录操作、进程管理、文件权限设定等。
所以,在Linux系统上工作离不开使用系统提供的命令。
要想真正理解Linux系统,就必须从Linux命令学起,通过基础的命令学习可以进一步理解Linux系统。
不同Linux发行版的命令数量不一样,但Linux发行版本最少的命令也有200多个。
这里笔者把比较重要和使用频率最多的命令,按照它们在系统中的作用分成下面六个部分一一介绍。
◆ 安装和登录命令:login、shutdown、halt、reboot、install、mount、umount、chsh、exit、last;◆ 文件处理命令:file、mkdir、grep、dd、find、mv、ls、diff、cat、ln;◆ 系统管理相关命令:df、top、free、quota、at、lp、adduser、groupadd、kill、crontab;◆ 网络操作命令:ifconfig、ip、ping、netstat、telnet、ftp、route、rlogin、rcp、finger、mail、 nslookup;◆ 系统安全相关命令:passwd、su、umask、chgrp、chmod、chown、chattr、sudo ps、who;◆ 其它命令:tar、unzip、gunzip、unarj、mtools、man、unendcode、uudecode。
本文以Mandrake Linux 9.1(Kenrel 2.4.21)为例,介绍Linux下的安装和登录命令。
login1.作用login的作用是登录系统,它的使用权限是所有用户。
2.格式login [name][-p ][-h 主机名称]-p:通知login保持现在的环境参数。
-h:用来向远程登录的之间传输用户名。
如果选择用命令行模式登录Linux的话,那么看到的第一个Linux命令就是login:。
精心整理UNIX 系统操作命令大全UNIX 系统操作命令第一章登录和退出用户在登录前,首先要向系统管理员申请一个用户注册名,不论用户从哪台计算机登录到ITPNET 上都将访问相同的文件系统。
1.1登录sun%注:[hide]1.2退出[Ctrl-D]Login :1.3update 等待话xargs第二章目录操作命令2.1显示目录命令:pwd (printworkingdirectory )格式:pwd功能:显示出用户当前工作目录的全路径名。
例如:%pwd2.2建立目录命令:mkdir (makedirectory )格式:mkdirdirectory_name功能:建立新目录例如:%mkdirdirectory2.3改变目录命令:cd(changedirectory)格式:cddirectory_name功能:改变当前工作目录注释:目录名选项有:?(1)目录名,转到当前目录下的子目录。
?(2)..表示转到当前目录的父目录。
?(3)/表示转到根目录。
?(4)没有选项时,转到用户的主目录。
例如:%cddirectory2.4命令:格式:注释:?(1)?(2)例如:3.1命令:格式:??选项:-a??-l?????-r??-g??-t例如:3.2命令:格式:catfilename功能:显示出文件的内容。
注释:当文件内容在屏幕上滚动显示时,可按Ctrl+S键屏幕滚暂停,按Ctrl+Q ??键继续显示。
例如:%catfilename3.3按屏幕显示文件内容命令:more格式:morefilename功能:按屏幕一屏一屏显示出文件内容选项:注释:在屏幕显示过程中有功能键。
空格键继续下一屏的输出;回车键看下??一行的输出;q键跳出more状态;/字符串查旬字符串所在处;a键返??回到上一屏。
例如:%morefilename3.4分页显示文件内容命令:pg(page)格式:pg[option]filename功能:分页显示指定文件的内容。
通用离散元用户指导(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)是一个处理不连续介质的二维离散元程序。
rest fall1.sav
delete range -2,2 -2,2
reset disp ;位移归0 reset hist
hist unbal
hist ydis 0,2
step 2000
plot hold block stress disp
edge
Joints jregion
Systems of Units
Hist
hist solve type 1
solve rat 1e-5
save tun1.sav
hist ydis 42 -67 sxx 42 -67 hist ydis 42 -73 sxx 42 -73 reset disp jdisp
solve rat 1e-5
save tun2.sav
Plot
pl syy fill projected
pl syy fill bfill
pl ydisp fill
pl bl id black
plot pen stress block;显示应力 pen作用是udec4.0窗口显示pl se ;看裂隙
pl
;抓取图像,在autocad内处理set pl dxf 256
copy xxx.dxf
Zone joint 力学属性
zone model elastic density 2E3 bulk 8E9 shear 5E9 range group ’block’group joint ’joint’
joint model area jks 5E10 jkn 2.5E10 jfriction 30 r ange group ’joint’
; new contact default
set jcondf joint model area jks=5E10 jkn=2.5E10 jfriction=30
insitu stress -2500000.0,0.0,-1.0E7
joint model residual jks 1E10 jkn 1E10 jfriction 20 jcohesion 1E6 range &group ’joints’; new contact default
set jcondf joint model residual jks=1E10 jkn=1E10 jfriction=20 jcohesion=1E6
Pline
set plin 0,21 2000,21 76
pl
set log on
print plin 1 ydisp
print plin 1 syy
set log off
Pause key
reset
reset disp hist
Movie
;打开图象抓拍功能,设置图片间隔时间以及窗口大小movie on
movie step 1000
window x1 x2 y1 y2
;movie off
rest seismov.sav
; make a movie of the groundfall ;
wind -12 12 -12 12
set ovtol 0.05
plot block vel max 2.0 blue stress max 50
movie on
movie file = seismic.dcx
movie step 1000
step 4000
;动画,先输出图再开始记录输出图的图,再设定文件名,再设定记录步长plot block; 输出图形
movie on ;动画开始
movie file=kala.dcx;设置动画文件名
movie step 500;动画记录步长
support
; line service tunnel
struct gen xc=42 yc=-70 npoint=16 mat=5 thick=0.37 fang=-11.25 theta 360 prop mat=5 st_d=2400 st_ymod=24.0e9 st_prat=0.20 st_yield=1e10
prop mat=5 if_kn=1e8 if_ks=1e7 if_coh=1.0e10 ; excavate main tunnel
del 28 32 -72 -68
reset disp
solve rat 1e-5
save tun3.sav。