UDEC实例翻译与命令解析.pdf

《UDEC实例翻译与命令解析.pdf》由会员分享，可在线阅读，更多相关《UDEC实例翻译与命令解析.pdf（85页珍藏版）》请在得力文库 - 分享文档赚钱的网站上搜索。

1、UDEC实例翻译与命令解析中铁隧道集团科研所珠穆浪玛UDEC实例翻译与命令解析翻 译：珠穆朗玛UDEC Version 4.0UDEC实例翻译与命令解析 中铁隧道集团科研所珠穆浪玛1 地震诱发地层坍塌 Seismic-Induced Groundfall1.1问题描述本例展示使用UDEC模拟分析地震诱发地层坍塌的一类的问题，模 型 见 图1.1,该模型基于加拿大安大略省萨德伯里市鹰桥公司弗雷则矿34-1-554切割断面的一个剖面图的结构和尺寸.用二维平面应变模型代表垂直于超采轴向方向的平面效应，超采面高5m,宽10m.Tiine-vaiyiug Stress AppliedFigure 1.1

2、 I DEC model for seismic-induced groimdfall假定两个连续节理交叉平面分析:一个角度为4 5度，另一个为-9度，两者节理间距均为5m,为了演示的目的，一个近似垂直的虚拟节理”也被添加到块体内开挖面顶部以增强不稳定性。围岩参数来自试验室平均测试数值，假定岩石块体参数如下：density 3000 kg/m3Ybiuigs modulus 75,000 MPaPoissons ratio 0.18假定块体仅具有弹性行为，节理假定符合库伦滑动准则，选择典型的教课书数值作为节理 参 数，如 下：joint noiiual stiffiiessjoint shea

3、r stiffiiessfriction anglecohesion20.000 MPa/m20.000 MPa/m300UDEC Version 4.0-2-UDEC实例翻译与命令解析 中铁隧道集团科研所珠穆浪玛初始应力状态按各向同性估计为24Mpa（假定垂直荷载由覆盖深度大约800m的岩层产生卜1.2 U D EC 分析UDEC模拟顺序分三个阶段，首先，模型在初始应力状态下进行无超采固结.其次，进行开挖并且模型循环至平衡状态.本阶段超采面周围的应力分布见图1 2超采正上方和下方的块体滑动后稳定.在第三阶段.估计了两个不同的峰值速度的地震事件.对所有地震模拟,在问题域的外周边界引入粘滞边界用

4、以消除波的反射.从而模拟有限的岩体，地震事件用施加到模型顶部y方向的正弦应力波表现.应力波被叠加到已存在的初始地应力上.在第一个模拟中，施加1.25Mpa的峰值应力,应当注意的是，由于粘滞边界条件实际是在模型顶部，施加的有效影响应力应该是1.25 MPa/2,or 0.625 MPa.0.02秒后的开挖面拱顶的应力分布见图1.3,两点的位移被监测,1点位于开挖面的左角，点2位于拱顶块体的右角，图1.4的位移时间曲线显示两点本质上是弹性反应Figure 1.2 Stress distribution around excavation at end of excavation stageUDEC

5、 Version 4.0-3-UDEC实例翻译与命令解析中铁隧道集团科研所珠穆浪玛Figure 1.3 Stress distribution in roof of excavation after 0.02 seconds(applied stress=1.25 x cos(In 100/)JOB TITLE:SEISMIC INDUCED ROOF COLLAPSEUDEC(Version 4.00)LEGEND27-Aug-04 23:10cycle 2871tme 1.9&9E-02sechistory plot-1.02E-04 1.14E-04-2.73E-O4 1.64E-04-

6、Vs.2.41E-O4 2.00E-02Itasca Consultng Group.Inc.Minneapolis.Minnesota USA(002)Figure 1.4 y-displacetnent histories for nvo points on excavation boundary(applied stress-1.25 x cos(2兀 100/)UDEC Version 4.0-4-UDEC实例翻译与命令解析 中铁隧道集团科研所珠穆浪玛本例关心的问题是在模型顶部施加的速度和计算速度的对比，下面的公式可以用以估计施加的波速.2(p Cp)(1.1)where=(K+(4/

7、3)G)/p 二使用这个方程，施加的最大波速大概是0.04m/sec,图1.5显示的峰值波速小于0.06m/sec.估计的波速和监测波速的不同在于使用的围岩模量.而是没有考虑节理变形的相等变形模量.Figure 1.5 Plot of y-velocity at top of model(applied stress=1.25 x co$(27r 100r)在第二个案例中，施加应力波峰值12.5 Mpa(有效应力6.25Mpa).0.02秒后的开挖拱顶应力 分 布 见 图1 6该图显示出拱顶岩体不受力，表面该块体已经松散并正在下落.对于关心的问题，后来三个时间的几何体和应力分布见图1.8至 图

8、1.10.在问题的顶部预测的波速(从上面的方程)是0.4m/sec.从模型中计算的波速见图1.11,再次，由于使用的是原岩弹性模量而不是岩体的变形模量导致预测和监测的波速之间的差异.1.3节包含了该模型的数据列表，该列表包含了一个FISH函数(show)被用来创建坍塌的动画文件，每隔0.02秒俘获一个显示的图片.通过改变FISH参数timejnt可以改变动画帧的UDEC Version 4.0-5-UDEC实例翻译与命令解析 中铁隧道集团科研所珠穆浪玛间隔.视图的总数也可以通过改变snap_shot的数值进行改变.为了显示80帧的显示图片而创建的该电影文件需要大概13MB的硬盘空间.Figur

9、e 1.6 Stress distribution in roof of excavation after 0.02 seconds(applied stress=12.5 x COS(2TT 100z)JOB TITLE :SE ISMIC INDUCE D ROOF COLLAPSEUDEC(Version 4.00)LE GE ND27-Aug-04 23:10cycle 2871time 1.d99E-02 sechtstory plot-124E-03 8.6&E-O4-9.03E-03-2.47E-08-Vs.2.41E-O4me 2.00E-02Itasca Ccnsultng

10、Group.Inc.Minneapolts.Minnesota USA(e-003)M2)Figure 1.7 y-displacement histories for two points on excavation boundary(applied stress=12.5 x COS(2TT 100，)，UDEC Version 4.0-6-UDEC实例翻译与命令解析中铁隧道集团科研所珠穆浪玛Figure 1.8 Stress distribution around excavation after 0.25 seconds(appliedstress=12.5 x COS(2T T 10

11、0r)JFigure 1.9 Stress distribution around excavation after 0.50 seconds(appliedstress=12.5 x cos(2/r 100/)UDEC Version 4.0-7-UDEC实例翻译与命令解析中铁隧道集团科研所珠穆浪玛Figure 1.10 Stress distribution around excavation after 0.75 seconds(appliedstress=12.5 x COS(2TT 100/)Figure 1.11 Plot of y-velocity at top of model

12、(applied stress=12.5 x COS(2TT i00r)；UDEC Version 4.0-8-UDEC实例翻译与命令解析 中铁隧道集团科研所珠穆浪玛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

13、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

14、,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 mode

15、l under field stressesbound stress=-24.0,0.0,-24.0 ygrad=-.O3 0-.03insitu stress=-24.0,0.0,-24.0 ygrad=-.O3 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 1 e-5;save consolidated statesave

16、 seismicl.sav;make excavationdelete range-5,5-2.5,2.5solve rat 1 e-5;save excavated statesave seismic2.savrest seismic2.sav;apply seismic load from top(peak velocity=0.04 m/sec)UDEC Version 4.0-9-UDEC实例翻译与命令解析;set up nonreflecting boundarybound mat=1bound xvisc range-26-23-21 21bound xvisc range 23

17、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 21reset 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)yve

18、l(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.eye time 0.02save seismic3.savrest seismic2.sav;apply seismic load from top(peak velocity=0.4 m/sec);set up nonreflecting boundarybound mat=1bound xvisc range-2

19、6-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)hi

20、st 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.中铁隧道集团科研所珠穆浪玛UDEC Version 4.0-10UDEC实例翻译与命令解析 中铁隧道集团科研所珠穆浪玛eye time 0.02save seismic4.sav;0.25 sec.eye time 0.23save

21、seismic5.sav;0.50 sec.eye time 0.25save seismic6.sav;0.75 seceye time 0.25save seismic7.savrest 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

22、问题陈述本例模拟展示了 UDEC在检查衬砌隧道方面的应用，着重强调了荷载在混凝土衬砌中的发展，本例也解释了模拟连续建造操作中独立阶段的模拟程序.隧道系统的理想几何体见图3.1.系统包含在海床下大约70m（中线）深度，中 线 间 距12m的两个隧道，初始水位在隧道中线上方110m处.服务隧道直径5.24m,衬 砌 厚 度37cm.主隧道 直 径8.22m,衬砌厚度46cm.服务隧道先于主隧道开挖和衬砌.随后设置主隧道衬砌,水位上升增加到100m.UDEC Version 4.0-11UDEC实例翻译与命令解析中铁隧道集团科研所珠穆浪玛Figure 3.1 Idealized geometry o

23、f service tunnel and main tunnel施工顺序是：(1)开挖服务隧道 excavation of the service tunnel;(2)衬砌服务隧道 lining of the service tunnel;开挖主隧道 excavation of the main tunnel;(4)衬砌主隧道 lining of the main tunnel;and(5)升高水位raising of the water level.分析的目的是评价每个施工阶段服务隧道和主隧道支护状况.本例的材料参数见下：岩体一开 挖 隧 道 的 围 岩 参 数 为：弹性模量 elastic

24、 modulus 0.89 GPa泊松比 Poisson s ratio 0.35单轴抗压强度 uniaxial compressive strength 3.5 MPa粘聚力 cohesion 1 MPa密度 density 1340 kg/rn3混凝土衬砌弹性模量为24 GPa,泊松比为0.19.假定衬砌为线弹性材料。3.2 UDE C 分析本问题创建的UDEC模型见图3.2,隧道中心线在Y=-70的 位 置，注 意，模型边界距隧UDEC Version 4.0-12UDEC实例翻译与命令解析 中铁隧道集团科研所珠穆浪玛道开挖面很近，本模型为了演示的目的只是提供一个快速计算，实际计算时应采

25、用大的模型。底边和侧边界采用滚轴固定，海底之上海水的重量按30m的水头施加在模型表面上，假定隧道为线性防水衬砌，因此不需要进行瞬态地下水流分析。忽略孔隙水压力，而将岩体单位重量设置为浮容重。垂直对水平应力比率假定为静水压力0.3o对本例而言，上 面 列 的5个施工顺序按三个阶段进行模拟，主隧道开挖和服务隧道衬砌按一个瞬时的活动进行模拟。主隧道衬砌和水位上升也假定瞬时产生，这些活动可以通过线性（梯 度）降低隧道周围反力或者设定产生一些松弛后安装支护而被分开模拟。在模拟的第一阶段，在重力应力被初始化到块体中后，服务隧道被开挖，UDEC循环直到达到平衡状态，弹性位移结果见图3.3.JOB TITLE

26、：UDEC(Version 4.00)LE GE ND27-Aug-04 2325cycle 1980time 7.365E-01 seczones in fdef blocksblock plotItasca Consulting Group.Inc.Minneapolis.Minnesota USAIO*!)UDEC model zoning with service tunnel excavatedUDEC Version 4.0-13UDEC实例翻译与命令解析中铁隧道集团科研所珠穆浪玛Figure 3.3 Elastic displacements due to excavation

27、of service tunnel在第二阶段，衬砌服务隧道和开挖主隧道.采用1 6个梁单元来模拟服务隧道混凝土衬砌,图3.4和 图3.5显示了主隧道开挖后的位移和主应力分布状况.注意，整个服务隧道向主隧道转移.UDEC Version 4.0-14UDEC实例翻译与命令解析中铁隧道集团科研所珠穆浪玛JOB TITLE:UDEC(Version 4.00)LEGEND27-Aug-04 23:27cyde 3440time 1.281E+00 secblock plotpnnapal stressesminimum=-1.882E+06maximum=-1.142E*040 1E 7struct

28、ural etements plottedItasca Consisting Group,Inc.Minneapolis,Minnesota USAro*iXXXx xx xGZ/Jr xLxr xlxxlxxrX或+.axra我a心wAxaXXEkwxxxx九rxxaxuxx/*x/*x/xx.会Figure 3.5 Principal stress distribution after mining of main tunnel在第三个阶段，一个额外的荷载被施加在模型的顶部用来模拟增加的1 0 0 m 的水的重量，使 用 S T R U C T apply pressure命令来施加静水压

29、力荷载（施加到隧道衬砌上）,用以代表隧道中心线上2 1 0 m 的水位.图3.6和 图 3.7显示了施加额外荷载后相应的轴力和弯矩分布，服务隧道最大轴力为5.59X106N,主隧道最大轴力为9.06X106N.Figure 3.6 Thrust 加 liners after water level is raisedUDEC Version 4.0-15UDEC实例翻译与命令解析中铁隧道集团科研所珠穆浪玛Figure 3.7 Moment in liners after water level is raised3.3 数据文件列表Example 3.1 TUNNEL.DAT;tunnel s

30、upport loadingblock 0-30 60-30 60-90 0-90round 0.1crack 0-70 60-70crack 30 0 30-90crack 42 0 42-90tun 30-70 4.11 12tun 42-70 2.62 8tun 30-70 5.5 12tun 42-70 5.5 12gen edge 2.0;initial stress statebound stress 1.02e5 0 1.02e5 ygrad 1.34e4 0 1.34e4;bound stress sxxo sxyo syyo ygrad sxxy sxyy syyy;sxx=

31、sxxo+(sxxx-x)+(sxxy-y)insit stres 1.02e5 0 1.02e5 szz 1.02e5 ygra 1.34e4 0 1.34e4 zgra 0 1.34e4gravity 0-10;rock propertiesprop mat=1 d=1340zone model mohrzone shear=.33e9 bulk=.99e9 coh=1e6 fric=30.0UDEC Version 4.0-16UDEC实例翻译与命令解析_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _;prop mat=1 d=1340

32、g=.33e9 k=.99e9 coh=1e6 fric=30.0;change con=3;elastic joint propertiesprop jmat=1 jkn=1 e9 jks=1e9 jcoh=1e10 jtens=1e10;cycle to initial equilibriumhist solve type 1solve rat 1 e-5save tun1.sav;excavate service tunneldel 40 44-72-68bound-1 1 -91 0 xvel=0.0bound-1 90-91-89 yvel=0.0bound 59 61-91 0 x

33、vel=0.0;histories around tunnel 1hist yd is 42-67 sxx 42-67hist ydis 42-73 sxx 42-73hist xdis 39-70 syy 39-70hist xdis 45-70 syy 39-70;histories around tunnel 2hist ydis 30-65.0 sxx 30-65.0hist ydis 30-75.0 sxx 30-75.0hist xdis 25.0-70 syy 25.0-70hist xdis 35.0-70 syy 35.0-70reset disp jdispsolve ra

34、t 1e-5save tun2.sav;line service tunnelstruct gen xc=42 yc=-70 npoint=16 mat=5 thick=0.37 fang=-11.25 theta 360prop mat=5 st_d=2400 st_ymod=24.0e9 st_prat=0.20 st_yield=1e10prop mat=5 if_kn=1e8 if_ks=1e7 if_coh=1.0e10;excavate main tunneldel 28 32-72-68reset dispsolve rat 1 e-5save tun3.sav;line mai

35、n tunnelstruc gen xc=30 yc=-70 npoint=8 mat=5 thick=0.46 fang=22.5 theta-360;add additional load representing raised water levelbound stress 0.0 0.0-1.0e6 range-1 91-31-29;add hydrostatic loads to tunnel linersstruct apply press 0.0 2.06e6reset dispsolve rat 1 e-5save tun4.sav中铁隧道集团科研所珠穆浪玛UDEC Versi

36、on 4.0-17UDEC实例翻译与命令解析 中铁隧道集团科研所珠穆浪玛ret4流体和动态荷载下的重力坝分析4.1 问题描述本例展示了一个坐落在节理岩体基础上的100m高混凝土重力坝建造问题，节理平均间距 是50m,节理方向为20度和-70度，主要研究了两个荷载条件,首先是研究包含流体在岩体裂隙内渗流的水库蓄水效应分析.其次，为了研究潜在地震类型的荷载,一个动态波被施加在模型的基础上进行分析.4.2 UDE C 分析本问题的UDEC模型见图4.10模型理想化的展示了一个坐落在节理岩体基础上的重力坝，目的是展示这种类型问题的推荐求解过程。数据文件见4.3节，按下面的顺序进行分析。Figure 4

37、.1 U D E C model of gravity dam with principal stresses plotted atStage 1第 一 步：重力荷载空的水库假定初始应力状态下的有效应力比率为oH/oV=0.69,水位假定在y=0,由大坝重量产生的初始应力状态和未蓄水的水库见图4.1o注 意，使 用INSITU命令设定的应力是总应力，被加载在块体上。对节理裂隙，UDEC计算有效应力并将域压力设置到静水压力。第 二 步：水库蓄水本阶段，假定水位上升到大坝的顶部，在大坝上游侧和岩石基础处施加静水压力作用，U D E C Version 4.0-18UDEC实例翻译与命令解析 中铁隧

38、道集团科研所珠穆浪玛在侧边界由于施加在大坝上的荷载产生的水平反作用力被认为是滚轴支撑。对于流体渗流假定下面的条件：1.假定沿着模型底部和侧边的节理裂隙连接渗透性为零。2.大坝上游的岩体面，通过使用B O U N D pp命令将水头固定在100m(0.98Mpa卜下游面，水头设置为零。3.大坝和岩体基础交界面假定有较低的渗透性。4.使用稳定流体算法(SET flow steady)选定的第二步结果见图4.2至 图4.6,蓄水后的位移结果见图4.20大坝顶部的X方向和Y方向位移历史见图4.3O后来的图形显示蓄水后模型处于平衡状态。图4.4的流体速率显示了大多数流体被集中在大坝基础下面的裂隙内。图4

39、.5显示了沿着节理(at x=-33.42,y=-30.37)的法向和剪切位移。正法向位移表明在本阶段节理张开。在位置x=-22.1,y=-26.3处沿着裂隙的流体压力历史见图4.6oFigure 4.2 Principal stress state and displacemetns ar Stage 2UDEC Version 4.0-19UDEC实例翻译与命令解析中铁隧道集团科研所珠穆浪玛Figure 4.3 x-andyconiponents of displacement at crest of da tn during Stage2 reservoir fillingFigure

40、4.4 Flow rates at Stage 2UDEC Version 4.0-20UDEC实例翻译与命令解析中铁隧道集团科研所珠穆浪玛Figure 4.5 Shear and normal displacements along joint at x=-33.42,j =-3037JOB TITLE DAM2 WATER LOADS ANO FLOWUDEC(Version 4.00)LE GE ND27-AugXM 22:52cy比 1790ume 3.890E 4-00 secflow nme-3.890E+00 sechtftofy plotFluid PressureV8.1.2

41、8E+O0 3 89E+00nasca ConsuRing Group.Inc.Minneapotis.Minnesota USAFigure 4.6 Fluid pressure in domain at location x=-22.1,j =-26.3第 三 步：动态荷载UDEC Version 4.0-21UDEC实例翻译与命令解析 中铁隧道集团科研所珠穆浪玛本阶段中，一个垂直传播的正弦波（frequency=5 Hz）被施加在模型基础上10秒种。下面的边界条件被采用。1.假定底部边界在水平（剪 切）方向为无反射边界且固定垂直方向，动态输入以剪切应力历史的形式被施加。2.类似的条件施加

42、在自由区域的的底部。3.基础一边的节点与自由区域的计算相关联。自由区域被离散化成20个 单 元，并设定为弹塑性行为。在自由区域使用FFIELD in it命令初始化当前静力水平以匹配模型中的网格。4.在动态阶段，由渗流阶段产生的静力作用仍旧被施加块体上。因 此，在无动态荷载作 用 时，块体仍处于平衡状态。在本阶段，渗流计算被关闭（S E T flow off）o这是近似假定当10秒动态荷载时无流体产生。图4.7至4.9显示了本阶段在1.5秒后动态荷载作用下的分析结果。每个荷载下大坝顶部X和Y方向的位移历史见图4.70在大坝下方节理的剪切和方向位移（图4.9）显示了剪切位移的累积情况。注 意，在

43、 图4.10中，沿着大坝基础下的第一个节理产生了大量滑移；上 游70度节理张开（无有效应力），这个边产生了较大位移，连同累计的位移显示在图4.8和4.9.表明大坝很可能破坏。Figure 4.7 x-andy-velocity histories at dam crest at Stage 3UDEC Version 4.0-22UDEC实例翻译与命令解析中铁隧道集团科研所珠穆浪玛Figure 4.8 x-andy-displaceinents at dam crest at Stage 3Figure 4.9 Joint shear and normal displacement histo

44、ries at joint beneathdam foundation at Stage 3UDEC Version 4.0-23UDEC实例翻译与命令解析中铁隧道集团科研所珠穆浪玛Figure 4.10 Displacements,principal stresses and joint shear ciisplacenienr atStage 34.3 数据文件列表Exanwle 4.1 DAM.DAT;dam discontinuous joints:20 and-70 deg.大坝-不连续节理20和-70度;insitu stresses(k=0.5)初始应力(K=0.5)?与文中叙述

45、不符，应为0.69否;free-field(20 nodes)applied only in phase 3 自由区域(20节点)，仅在第三阶段应用;第一阶段一重力荷载;固定X方向边界;phase 1 -gravity loading;x-fixed boundariestitle(DAM1)IN-SITU STRESSES;GRAVITY*;rock blocks;joints(no cohesion)岩块 节理无黏聚力prop mat=1 d=0.00265 k=33333 g=20000prop mat=1 jkn=1000 jks=1000 jf=30.0;节理参数为啥不用prop j

46、mat命令呢？下同;dam and foundation joint 大坝和基础节理prop mat=2 d=0.00240 k=16667 g=12500prop mat=2 jkn=1000 jks=1000 jf=30.0 jcoh=2 jtens=2;above y=-150:free field(low cohesion);joints(no cohesion);y=150上部自由区域(低黏聚力)节理无黏聚力prop mat=3 d=0.00265 k=33333 g=20000 f=30.0 coh=0.2UDEC Version 4.0-24UDEC实例翻译与命令解析 中铁隧道集

47、团科研所珠穆浪玛prop mat=3 jkn=1000 jks=1000 jf=30.0;below y=-150:free field(cohesion);joints(cohesion);y=-150下的自由区域(黏聚力)，节理(黏聚力)prop mat=4 d=0.00265 k=33333 g=20000 f=30.0 coh=2prop mat=4 jkn=1000 jks=1000 jf=30.0 jcoh=2;viscous boundaries(equivalent elastic properties);joints(cohesion);粘结边界(等价弹性参数)，节 理(黏聚

48、力)prop mat=5 d=0.00265 k=11680 g=11111prop mat=5 jkn=1000 jks=1000 jf=30.0 jcoh=2round 0.5;set minimum edge length 设置最小边长set edge 8.0;set minimum contact length 设置最小连接长set clemin=5.0block-200,-200-200,100 200,100 200,-200 模型块体;structure:gravity dam 结构 重力坝crack-201,0 201,0;根据几何形体对块体进行分割crack-40,-1-40

49、,101crack 40,0-40,100delete range-200,-40 0,100 删除多余部分，形成坝体delete range 40,200 0,100jregion id=1-200,-200-200,0 200,0 200,-200jset 20,0 800,0 0,0 50,0(50,0)range jreg 1jset-70,0 50,0 50,0 50,0(-50.99,16.45)range jreg 1jset-70,0 50,0 50,0 50,0(-10.40,-21.98)range jreg 1;节理生成区域-基础区；20度；-70 度；-70 度chan

50、ge mat=1 range-200,200-200,0change mat=2 range-40,40 0,100generate edge 60 range-200,200-200,0generate edge 30 range-40,40 0,100;all joints所有节理change jmat=1 range-210,210-210,10;cohesion below y=-150change jmat=4 range-210,210-210,-150;foundation joint 基础下节理change jmat=2 range-41,41-1,1 ang=-5,5;坝下基