自编强度折减法命令

hql1225 · 发表于 2008-10-21 23:25:06

算是能够算完的，只是收敛的很慢，时间很长，最后算出来结果是1.849609375000e＋000，来慢慢研究下这个咯

Lyhn · 发表于 2008-10-22 16:04:00

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

rocky5678 · 发表于 2008-10-23 15:40:31

学习学习

chengjia02 · 发表于 2008-10-24 15:26:40

楼主我采用此多层赋值
运行一个循环后弹出：bad type conversion
不知道是什么原因呢~

mifan · 发表于 2008-10-26 09:53:47

我也出现ls的问题，请大牛解释下！另外建议针对常出现的error notification 能够专门地整理，让新手能够更系统的学习下，o(∩_∩)o...，先谢谢！

nosil007 · 发表于 2008-12-2 21:41:11

楼上的兄弟说的很对，我算过了可行

liuyongseu · 发表于 2008-12-2 21:45:58

支持支持支持支持

yuanpengbo · 发表于 2008-12-13 13:43:07

谢谢楼主。谢谢分享。

csuroylin · 发表于 2009-2-15 05:40:56

楼主，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怎么确定多少步，另外 ...
jonavy 发表于 2008-6-5 23:32

到了30000步还没完，就是说明不收敛，将降低折减系数。
30000和折减系数的上下限都是人为定出来的。
对于洞室不太清楚，边坡是所有的岩土体都折减。

csuroylin · 发表于 2009-2-15 05:43:27

我也想问，flac3d中有强度折减法，为何还要自己编写程序来折减呢？
topccl 发表于 2008-9-19 09:56

比如有的时候，我们考虑锚杆参数对边坡稳定性的影响，
我们也能知道变化过程中，安全系数的变化大致范围，
则用这种方法能够有效的提高计算效率。
另外，目前solve fos 只是在mc准则中应用，在其他准则中，
就要自己编写命令流了。

csuroylin · 发表于 2009-2-15 05:45:38

楼主，您好！
很感谢你编写的这个自定义强度折减，但我觉得还是有些问题请教：
1、强度折减中采用什么作为收敛标准现在有很多人还在讨论这个问题，你这里好像采用的是不平衡力，但
if mech_ratio
topccl 发表于 2008-9-19 10:48

loop while (k12-k11)>ait1 是计算的标志，
(k12-k11)<=ait1 就结束了
另外，如果你计算的够多的话，最终是1.75的

ulzeng · 发表于 2009-2-23 20:57:30

谢谢。学习中。。。

li_2001 · 发表于 2009-3-21 19:33:37

好啊好支持

li_2001 · 发表于 2009-3-24 10:45:25

可以可以支持支持

whlxpp521 · 发表于 2009-3-26 14:54:35

luhuan · 发表于 2009-4-9 00:20:34

学习，强顶此贴

umengzhu · 发表于 2009-5-5 19:00:47

好东西啊！
学习中！

东门下雨 · 发表于 2009-5-13 14:07:26

谢啦~~
下来看看学习一下！

ezheng_08hi · 发表于 2009-5-14 22:20:27

目前，强度折减法的研究方兴未艾，FLAC3D只能对于Mohr-coulomb准则实施强度折减技术，通过solve fos，
而且计算速度比较慢，为了提高计算速度，即自定义强度折减法折减系数的上下限值，编了个程序，请大家
斧正。 ...
chunhuo420 发表于 2007-5-21 16:56

根据陈育民、徐顶平编的FLAC3D基础与实例中14章有折减法。与楼主的基本一致。
只是那里在
;==================! B" p9 e' j( L- s8 c3 L. _; y% R9 {
command
8 R. f+ K/ S" W: p7 \model null    ;挖空# i6 y7 h! V# H; h1 [* V9 E
model mohr  ;采用莫尔库仑模型* o1 \% V1 J9 T( ~
pro bulk K1 she G1 dens dens1 coh coh1 &
friction fri1 dilation 0. tens ten1
fix x range x -0.1 0.1
fix x range x 104.9 105.1( f' g: p1 N( h  C" T1 q
fix y1 u$ q& V4 Q* w# c, Z  @
fix x y z range z -0.1 0.1. l) Q- O3 v0 f
set grav 0 0 grav0; }8 u: p! ?" o0 i
set mech ratio 9.8e-68 t0 p, t! o- g3 f2 I
solve step 30000! m  h( i% c4 I. ^! M7 w3 E7 A
endcommand
===================
部分
那书里循环里每次都得先用弹性求初始应力场，然后再塑性求解。
感觉计算起来也时间挺长的。

想问与这里的二者有区别吗？

@' Z7 T1 b6 R3 D5 l

flyeagle009 · 发表于 2009-6-5 10:05:40

model null ;挖空
有何作用啊？

账号		自动登录	找回密码
密码			注册

[命令/FISH] 自编强度折减法命令

回复 20# chunhuo420 的帖子

浏览过的版块