15.建立interface的几种方法首先给出一个简单的例子new
genzonebrickp0000p1111.400p20100p30045size1115group1
genzonebrickp00045p1111.4045p201045p30090p4111.41045&p501090p647090p7471090size1119group3
genzonewedgep066.4090.1p147090.1p266.41090.1p3111.4045.1&p4471090.1p5111.41045.1size519group4
;interface1facerangplanenorm45,0,64.4origin47590interface1facerangeplanenorm45,0,64.4origin47590dist0.1
inizadd-0.1rangegroup4
genzonebrickp0111.400p112800p2111.4100p3111.4045&
size215group2modelmohr
probulk=1e8she=3e7dens2e3coh12e3&friction20dilation0.tension1.0e10rangegroup1probulk=1e8she=3e7dens2e3coh12e3&friction20dilation0.tension1.0e10rangegroup2int1propks2e12kn2e12fric36
probulk=1e8she=3e7dens2e3coh12e3&friction20dilation0.tension1.0e10rangegroup3probulk=1e8she=3e7dens2e3coh12e3&friction20dilation0.tension1.0e10rangegroup4fixzrangez-0.10.1fixxrangex-0.10.1fixxrangex127.9128.1fixy
setgravity00-10step1000inistate0
inixdis0ydis0zdis0;inixvel0yvel0zvel0;setlargestep1000
;solvefosfilecable.savassociated修改上面蓝色字体部分
第一种建立interface的方法（书上的移来移去式）：
;interface1facerangplanenorm45,0,64.4origin47590interface1facerangeplanenorm45,0,64.4origin47590dist=0.1
size215group2
通过计算得到的xdisp云图
plintyelsket
问：为什么第一种方法要加dist0.1？加入节理以后，这个例子就不能计算fos了？
回答：你的第2种，把group2和group4删掉，这样造成的新的模型group2和group4有重复节点，这样是不正确的。方法是在一个实体的外表面全部都设置了接触面，就是可以将两个体完全分开。而你的这个接触问题，在坡脚的位置很难处理。只用把斜坡上的部分删掉，然后导入就可以了。
人为的把group1和group2之间的交界面上多出了节点，下图是两个group的节点图，是分开显示的，可以看到两个共同面上界面标号不同，这样是不对的。计算结果也同样是错误的。
16.callfishcall.fis总是会出错呢？老是给出erroropeningfile，该怎么操作？
错误是：文件名为fishcall.fis.txt(因为在记事本里这个txt的后缀是默认的）。修改方法是：将fishcall.fis.txt里面的内容全部复制到新的记事本里，再按保存，保存时候将保存类型改为所有文件(*,*),这时你便可在文件名栏内输入fishcall.fis，类型也就正确了
17.问题：
fixzrangez-0.10.1
fixxrangex-0.10.1anyx1.92.1any
fixyrangey-0.10.1anyy1.92.1any
和
fixxrangex-0.10.1
fixxrange1.92.1
fixyrangey-0.10.1
fixyrangey1.92.1
问：两者有区别吗？我感觉是一样的
回答：建一个模型，然后实验一下，然后用plotgpfixt看一下。
18.得到初始应力的方法：
方法1、可以先给一些材料参数很大的值，进行初始求解，在计算之前再将材料参数设为正常值，即可。如在手册中给的第一个示例中就是这样做的。下面是例子，Theseareonlyinitialvaluesthatareusedduringthedevelopmentofgravitationalstresseswithinthebody.Ineffect,weareforcingthebodyto
behaveelasticallyduringthedevelopmentoftheinitialin-situstressstate.*Thispreventsanyplasticyieldduringtheinitialloadingphaseoftheanalysis.
Genzonebricksize688
Modemohr
Propbulk1e8shear0.3e8fric35
Propcohesion1e10tens1e10;注意在此这个值给的很大。
Initdens1000
Setgravity00-10
Fixxrangex-0.10.1
Fixxrangex5.96.1
Fixyrangey-0.10.1
Fixyrangey7.98.1
Fixzrangez-0.10.1
setmechforce=50
solve
;----------------------modenull---------------------
Propcoh1e3tens1e3;改为正常值(在此例中我们故意给小值)Modenullrangex2,4y2,6z5,10
Setlarge
Inixdis0ydis0zdis0;清零，不影响结果，为画图方便。
Setp2000
return
方法2：得到初始应力的一个不错办法是设为弹性模型，然后赋予地勘信息，so以后就ok了，得到初始应力了，再把相应的位移清零，下面就可以正式进入所分析的内容了。如上例还可以先设为弹性(MODELelastic)，之后在开挖之前再设为Mohr-Coulombmodel，就行。
19．关于初始应力的问题
Example3.21Initialstressstatewithgravitationalgradient
forexample,a20m20m20mboxof
homogeneousmaterialatadepthof200munderground,withfixedbaseandstressboundariesontheothersides.
genzonebricksize10,10,10p120,0,0p20,20,0p30,0,20modelmohr