自编强度折减法命令

freeflywyx

想问下楼主  
计算得到安全系数后怎么找滑面
请教

diwurou

多谢版主的指点。

diwurou

以计算不收敛来作为破坏标准的？

chunhuo420

原帖由 diwurou 于 2008-3-21 13:23 发表
以计算不收敛来作为破坏标准的？

是的。。

chunhuo420

原帖由 freeflywyx 于 2008-3-17 14:27 发表
想问下楼主  
计算得到安全系数后怎么找滑面
请教

我是用位移等值线来表示滑面的。

boris1949

稳定性算法很多，应该也可以类似对其他方法编程吧

王下长合

请教楼主：在命令流中set mech ratio 9.8e-6是针对这solve step 30000一个循环么？而if mech_ratio<1.0e-5 针对什么的？是整个么？mech_ratio在前面我没有看到对其定义啊？这个怎么回事？问了这么多，麻烦楼主了。

chunhuo420

原帖由 王下长合 于 2008-5-5 15:53 发表
请教楼主：在命令流中set mech ratio 9.8e-6是针对这solve step 30000一个循环么？而if mech_ratio

(1) 是的
(2) 判断计算是否收敛，具体流程见下图
(3) mech_ratio 是系统变量，表示收敛精度

王下长合

感谢楼主，学习了。

rainbowyun

xiexie感谢

jonavy

原帖由 chunhuo420 于 2008-5-6 08:09 发表





(1) 是的
(2) 判断计算是否收敛，具体流程见下图
(3) mech_ratio 是系统变量，表示收敛精度

楼主，set mech ratio 9.8e-6
          solve  step 3000
  前面一句set mech ratio 9.8e-6是针对这solve step 30000一个循环的话，这个step 3000 完了mech还不到9.8e-6怎么办，这个step怎么确定多少步，另外那个上下限到底怎么定的那，要是试算出来的话有时模型太大是不行的，上下限人为性太大了，当考虑局部洞室时，折减一般指针对洞室衬砌可以么，还是还要将围岩全部考虑进去那？

023wangliang

非常感谢版主  谢谢你们的讨论

023wangliang

非常有用
我也正在找安全系数这的算法
顶下

michael02

谢谢楼主分享啊

topccl

程序不收敛，楼主解答！！！

emperor810

谢谢了- -

yc123292

太感谢楼主了，无法表达了！

yc123292

版主，感觉计算过程和求解fos没有什么区别啊，计算简单的模型速度倒更慢！

topccl

我也想问，flac3d中有强度折减法，为何还要自己编写程序来折减呢？

topccl

楼主，您好！
很感谢你编写的这个自定义强度折减，但我觉得还是有些问题请教：
1、强度折减中采用什么作为收敛标准现在有很多人还在讨论这个问题，你这里好像采用的是不平衡力，但
if mech_ratio<1.0e-5
  k11=ks
  k12=k12
else
  k12=ks
  k11=k11
endif
ks=(k11+k12)/2
这段程序不明白是什么意思？举一个简单的例子：假如k11=1，k12=2.5，边坡的安全系数有恰好是ks=(k11+k12)/2=1.75,那上述程序段在一次循环后，mech_ratio<1.0e-5 即满足条件，那么此时k11=ks=1.75,k12=2.5,然后ks=ks=(k11+k12)/2=(1.75+2)/2=1.875,就和先前的1.75不合？？
2、你给的附件程序一直就是计算，好像根本没有结束计算的意思，我算了几十万步了还在计算，所以我也觉得程序有问题。
希望多多交流！