zoukankan      html  css  js  c++  java
  • UDEC 1

       

    • 在下面这个最简单的模型中,切割块体涉及选择位置和指定裂缝,通过键入这些命令就可以产生10x10个单位的块体,然后劈裂成两个块体
    • crack命令产生一个连续、水平的贯通模型的裂缝
    • 圆角长度指定为0.1
      round 0.1          ;round 指块体与块体之间的圆角半径,默认值是0.5,其值要求小于模型中最小块体的最短那条边长的二分之一。如:round 0.05 set ovtol=0.5   
      block(0,0) (0,10) (10,10) (10,0)       
      Crack(0,5) (10,5)      ;crack x1,y1 x2,y2 两点划一线
      
      ;通过键入如下命令可以产生一个槽口:
      
      crack 2.5,10 5.0,7.5        
      crack 5.0,7.5 7.5,10
      ;通过键入以下命令,就显示出包括块体地址号的块体图形
      plot block num
      ;通过采用delete命令,能从模型中删除一个块体,比如删除槽口
      delete range block 647

        crack命令用于产生块体中的单一直线特诊的裂缝,裂缝由端点坐标(x1,y1)和(x2,y2)所确定 。                                                                                    jset命令则是自动节理组生成器。根据所给定的特征参数(即倾角、迹长、岩桥长度、间距和空间位置)产生一组裂缝。

             

    • 两组连续节理组的产生
    • 第一个jset命令产生一组与x轴方向夹角为45°具有间距为2个单元的连续节理。
    • 第二个jest命令产生与x方向夹角为-10°间距为1.5个单位的连续节理。
    • 圆角长度的选择可能影响节理组的产状。
    • 节理的位置
      new
      round 0.01
      block(0,0) (0,20) (20,20) (20,0)
      jset(45,0) (5,0.5) (0.5,0) (2,0)       ;jset 90,0 4,0 4,0 6,0 0,-50 range jreg 3        jset 倾角,0 线段长,0 线段与线段轴向间隔长,0 垂向间距,0 xm,ym range jregion n
      jset(-10,0)(5,0.5)(0.2,0) (1.5,0)

       遇到了错误却找不对原因/(ㄒoㄒ)/~~ 一个基础案列就在这里卡壳了 

         一开始准备按照教程的基础案例一个一个敲来进行学习,但后来发现这样可能效率不高,掌握的慢,所以就想直接粘贴一个具体案例来搞明白每一块命令行是啥意思,但两三百行的代码从十几行开始就出错,别人在4.0上可以成功运行并且呈现理想的结果,我的7.0就不行T_T。

    下面这图可能是先模拟出不同岩层的模型高度,执行很多次crack命令,接下来就开始出错。目前,找错误原因还是一大难题。

    所以先总结一下用到的各种命令行是啥意思。

      jregion id n x1,y1 x2,y2 x3,y3 x4,y4 delete

      jset 90,0 4,0 4,0 6,0 0,-50 range jreg 3

      ;jset 倾角,0 线段长,0 线段与线段轴向间隔长,0 垂向间距,0 xm,ym range jregion n

      ;其中xm,ym为起始点坐标,n为设置的区域标号

    gen quad 10 range xl xu yl yu
    ; 在指定的区域生成一定宽度的单元
    
    prop mat=1 dens=2000 
    prop mat=2 dens=2650
    prop mat=3 dens=2700 
    ;指定各材料的密度,比如1号材料dens=2000,即1立方米重2吨
    
    
    prop jmat=1 jkn=0.2e8,jks=0.1e7,jcoh=0,jfric=4,jten=0
    prop jmat=2 jkn=8e8,jks=5e7,jcoh=0.1e6,jfric=8,jten=0
    prop jmat=3 jkn=20e8,jks=16e7,jcoh=0.4e6,jfric=15,jten=0  
    ;jkn为节理的法向刚度,jks为节理的切向刚度,jfriction为节理的内摩擦角,
    ;jcohesion为节理的内聚力,jtension为节理的抗拉强度
    
    set gravity 0,-9.81
    ;设置重力加速度,x方向为0,y方向为-9.8
    
    
    bound xvel=0 range -0.1 3.00 -60.1 20.1 
    bound xvel=0 range 97 100.1 -60.1 20.1 
    bound yvel=0 range 0.1 100.1 -60.1 -58 
    ;采用位移法固定边界
    
    
    set plin 1 x1,y1 x2,y2 n
    ;输出观测线的位移与应力数据,默认输出到udec.log文件中,可在excel表中处理数据

     :在摸索了许久之后我似乎明白了UDEC是基于命令驱动格式这句话。

  • 相关阅读:
    HDU1255 覆盖的面积 —— 求矩形交面积 线段树 + 扫描线 + 离散化
    HDU1542 Atlantis —— 求矩形面积并 线段树 + 扫描线 + 离散化
    HDU1540 Tunnel Warfare —— 线段树 区间合并
    HDU3974 Assign the task —— dfs时间戳 + 线段树
    HDU4027 Can you answer these queries? —— 线段树 区间修改
    POJ3264 Balanced Lineup —— 线段树单点更新 区间最大最小值
    ZOJ1610 Count the Colors —— 线段树 区间染色
    HDU1698 Just a Hook —— 线段树 区间染色
    POJ2528 Mayor's posters —— 线段树染色 + 离散化
    POJ3468 A Simple Problem with Integers —— 线段树 区间修改
  • 原文地址:https://www.cnblogs.com/raindawn/p/11923333.html
Copyright © 2011-2022 走看看