qiandao 发表于 2008-2-20 11:29:36

FLAC能否显示破坏接近度（或屈服接近度）？

隧道开挖后，想看看围岩的破坏接近度的等值线分布，FLAC该怎么实现？谢谢

xbwei33 发表于 2008-2-20 17:33:33

在FLAC 中，不能显示等值线，只有采用其他的后处理～！！！
tecplot 该软件来实现！！

qiandao 发表于 2008-2-24 10:01:05

谢谢，我试试

lessse 发表于 2008-2-24 13:05:39

屈服接近度的公式有了就可以了，只需编个fish可以显示！

wgy0207 发表于 2008-4-3 15:49:01

请问什么是屈服接近度？？

wangchengcumt 发表于 2010-3-12 10:36:11

4# lessse 请问高手屈服接近度的公司有了，如何编fish语言啊

李木木 发表于 2010-10-19 16:47:55

为什么我编的fish屈服接近度值的范围不在0 1之间

adina81 发表于 2010-11-9 00:05:24

我编写了一个，也不行，大家可以帮我看看吗？

title '屈服接近度计算'
rangename weiyan1group 2 ;对需要输出屈服率的组进行命名

config zextra 1
def yeild1
pnt=zone_head
loop while pnt # null
if inrange('weiyan1',pnt)=1 then
c0=0.1e5;----------输入粘聚力
fai0=20*3.14/180    ;----------输入内摩擦角
st=0.1e7   ;----------输入抗拉强度
s3=z_sig1(pnt)
s2=z_sig2(pnt)
s1=z_sig3(pnt)
s13=(s1+s3)/2.0

;应力罗德角
ylldj=(atan(1.0/(sqrt(3)))*(2.0*s2-s1-s3)/(s1-s3))

i1=s1+s2+s3
j2=(1/6.0)*((s1-s2)*(s1-s2)+(s1-s3)*(s1-s3)+(s3-s2)*(s3-s2))

;GR计算
sigmr=(s1*(2-sin(fai0))-2*c0*cos(fai0))/(2*(1-sin(fai0)))

if s13>sigmr then
z_extra(pnt,1)=(st-s1)/(st-sigmr)
if z_extra(pnt,1)>1 then
z_extra(pnt,1)=1
else
endif
if z_extra(pnt,1)<0 then
z_extra(pnt,1)=0
else
endif
else
z_extra(pnt,1)=(i1*sin(fai0)+(3*cos(ylldj)-sqrt(3)*sin(ylldj)*sin(fai0))*sqrt(j2)-c0*cos(fai0))/(i1*sin(fai0)-3*c0*cos(fai0))
if z_extra(pnt,1)>1 then
z_extra(pnt,1)=1
else
endif
if z_extra(pnt,1)<0 then
z_extra(pnt,1)=0
else
endif
endif

endif
pnt=z_next(pnt)
endloop
end

yeild1
plot contour zextra 1 alias 'Contour of yield rate zone' out on

chen575982 发表于 2010-12-1 10:30:57

可惜看不懂

songli269 发表于 2010-12-2 17:43:21

得到的屈服接近度相反，危险区域屈服接近度反而小。

songli269 发表于 2010-12-2 22:40:25

有谁知道原因啊，告诉一下啊。

李木木 发表于 2011-1-8 23:26:46

adina81 我感觉你得sigmr=(s1*(2-sin(fai0))-2*c0*cos(fai0))/(2*(1-sin(fai0)))
有问题

小志 发表于 2011-1-9 13:51:31

还是很给力呀

adina81 发表于 2011-1-10 01:52:53

sigmr=(s1*(2-sin(fai0))-2*c0*cos(fai0))/(2*(1-sin(fai0)))
这个是在塑性力学里面查到的，但是我用起来就是对的啊

李木木 发表于 2011-1-10 13:31:19

adina81：文章里写的是：sigmr=(st-c0*cos(fai0))/(1-sin(fai0)) ！！你在哪里查的呢？见张传庆博士论文；我有但我还不会在里传东西:$ ！你要是要就给我留e-mail

wc2010 发表于 2011-9-25 15:59:39

adina81 发表于 2010-11-9 00:05 static/image/common/back.gif
我编写了一个，也不行，大家可以帮我看看吗？

title '屈服接近度计算'


您好！请问您这个是用于3D还是2D的啊？最后的图形怎么出来呢？用什么命令？我用2D的话怎么出最后的等值线啊？

nick_h1201 发表于 2011-10-12 13:53:06

adina81 发表于 2011-1-10 01:52 static/image/common/back.gif
sigmr=(s1*(2-sin(fai0))-2*c0*cos(fai0))/(2*(1-sin(fai0)))
这个是在塑性力学里面查到的，但是我用起来就 ...

兄台 你这段程序出来的屈服接近度为何会在0,1区间意外我看已经写了限制性的if语句啊 没搞明白。。还请指点！

wc2010 发表于 2011-12-4 23:43:51

adina81 发表于 2010-11-9 00:05 static/image/common/back.gif
我编写了一个，也不行，大家可以帮我看看吗？

title '屈服接近度计算'


if inrange('weiyan1',pnt)=1 then    不知这句是什么意思，兄弟能不能支教下...谢谢！

侯健 发表于 2014-6-19 11:35:31

李木木 发表于 2011-1-10 13:31
adina81：文章里写的是：sigmr=(st-c0*cos(fai0))/(1-sin(fai0)) ！！你在哪里查的呢？见张传庆博士论文； ...

楼主，能否给我发一份，最近在研究弱胶结软岩的内容，谢谢！这是我的邮箱：kdhoujian@126.com；十分感谢

chenpf09 发表于 2014-6-19 12:47:57

收藏了。。。。。。。。。。。。备用
页: [1] 2
查看完整版本: FLAC能否显示破坏接近度（或屈服接近度）？