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

qiandao

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

xbwei33

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

qiandao

谢谢，我试试

lessse

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

wgy0207

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

wangchengcumt

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

李木木

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

adina81

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

title '屈服接近度计算'
range  name weiyan1  group 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

可惜看不懂

songli269

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

songli269

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

李木木

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

小志

还是很给力呀

adina81

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

李木木

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

wc2010

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

title '屈服接近度计算'

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

nick_h1201

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

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

wc2010

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

title '屈服接近度计算'

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

侯健

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

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

chenpf09

收藏了。。。。。。。。。。。。备用