《fishtank--巴西实验与双轴压缩实验 说明.doc》由会员分享，可在线阅读，更多相关《fishtank--巴西实验与双轴压缩实验 说明.doc（3页珍藏版）》请在金锄头文库上搜索。
1、fishtank巴西实验与双轴压缩实验 说明许多同学可能都对这部分很疑惑，这部分在manual中显得比较神秘，原因很简单，是因为这不风的程序代码比较复杂。仅仅et2.fis一个文件就长达30多页。但实际上我们在调用这些函数的时候，不需要逐一去看。我们要做的是理解每个函数的功能及函数间的参数传递。0 # n3 y. u q3 n2 J5 JY我将部分稍做整理如下，对部分文件和函数进行了初步解读： R9 e4 t$ n5 r- . I4 v- v0 $ Uq( Er# z% L$ R一：文件+ Z6 X, x, X1 d+ k7 N1 f实例中提供了不同类型实验的文件，以AGC.dvr为例。这个文
2、件看起来很怪，其实很简单。可以用txt或其他文本编辑器打开。也可以令存为txt文件。其中A表示材料类型，G表示实验精度，gross；c表示粒子的接触模型，为contactbond。. G1 + P) + |, x/ p; Z+ G# 4 & W 打开这个文件，其文件编辑和普通的建模文件没两样。这个文件开始就调用了三个fish函数文件，分别是md.fis, et2.fis, flt.fis。（稍后详细说明这三个文件）。就是参数设置，包括文件自动运行后保存的文件名的赋值，与具体fish函数中参数的赋值。在agc.dvr文件运行后将生成以下文件：; - W7 F* P( x! Q: D6 m- Sa
3、gc-bal.sav 生成的粒子集合，没有运算7 I9 * H. U# q& r3 b: Q5 Q8 Xagc-pck.sav 接上一步调整模型粒子大小，达到近似无摩擦状态的平衡9 r( TT- J+ U. h3 zagc-iso.sav进一步调整并赋值粘结模型，是模型达到各向同性应力状态! d6 p6 m; y6 qagc-flt.sav 消除模型的中的浮动粒子，默认值为接触数小于2的粒子；2 M2 X# A/ K. Kl8 * Pagc-spc.sav 生成最终模型3 i% b5 o* Q$ M* v! X接下来就是调用agc_ex.dvr 文件进行压缩实验了；2 _% _3 Y d. k
4、3 V0 q这个文件将巴西圆实验与双轴压缩实验一同解决了，都用的同一个样本，即前面生成的文件agc-spc.sav 中的模型。在这个dvr文件中调用了两个dvr文件， 即：_brazil1.dvr 和 _btw.dvr。 这两个文件分别是调用巴西圆测试功能函数与双轴压缩功能函数。生成的文件分别为：3 u/ q$ u% w& n ?7 d4 B: _巴西测试：agc_bt0.sav 巴西测试条件初始化3 j% J1 v9 # L agc_bt1.sav 巴西测试结果 lQ5 $ m) c/ h: 6 n6 x( H5 q$ t9 |$ |0 p8 r3 8 h/ l: h: y双轴压缩：agc_
5、nn_bw0.sav 测试条件初始化3 3 O: M# a% agc_nn_bw1.sav 测试结果 B: I6 f6 G: O 这里的nn表示围呀的大小，在manual中，粒子的围压为0， 所以保存的文件就是agc_00_bw0.sav和agc_00_bw1.sav4 f3 W6 r7 m, f1 s) t- w; F! N; G! o4 M% l, U# V; - z0 z; h到此为止，基本上将运行中相关文件和运行的大致过程说完。补充说明在manual中出现的相关的其他文件：. ( p p& g& _ _2 + B这部分的文件有头三个字母加横短线后的补充说明部分组成。T i7 E1 H,
6、 _7 j6 , q |) w前三个字母分别表示材料类型，计算精度和接触的本构模型。# u& J. + 2 T6 G2 W! r- d材料类型由ABC来表示；计算精度分别与哦那个g,c,f,v表示，分别表示总体，粗糙，精确，非常精确，四类；接触的本构模型分别用C,P,U表示，分别为接触粘结，平行粘结与散体）。因此在mannual中出的相关文件名就不难理解了，如afc表示a行材料，精确计算精度，接触粘结实验。二 参数 （还是以agc.dvr 为例）* l+ D; a- _4 Q! r! e. s z在agc.dvr文件中，除了调用几个fish函数外，就是参数设置，正确理解这些参数就可以任意调整实
7、验以达到生成自己特殊模型实验的目的。在mannual中，提到了两类重要参数的设置，并对这些参数进行了解释。$ k7 / P9 u% 0 S; I( L, l( T第一类：样本模型生成及相关过程中的参数 P8 F& J# h8 j* u2 F+ et2_ylen样本高et2_xlen样本宽et2_rlo最小粒子半径et2_radius_ratio最大与最小粒子半径比例md_wEcfac墙体法相刚度的扩大器tm_req_isostr内锁各向同性应力flt_def非浮动粒子最少接触数flt_remain剩余的浮动粒子所占比例$ O0 W5 R; f# m) p, # L# A( p5 V( - S4
8、 V& b8 f5 d第二类：接触本构模型的微观参数# w0 3 m$ a d$ d o# z* 在这类参数中，不同的接触模型要求的参数不一样，所以理解接触模型的意义很重要。/ R w+ dG5 L8 MD) y! _以接触粘结模型为例，因为接触粘结模型是点接触型，与平行接触比较没有刚度，也没有粘结圆半径等参数。其需要的参数如下表： md_dens密度md_Ec9 o; P# H/ Y7 . Q# D接触模量md_knoverks球法向与切向刚度比例md_fric球摩擦系数cb_sn_mean接触粘结平均法向强度cb_sn_sdev接触粘结偏载法向强度cb_ss_mean接触粘结平均切向强度cb_ss_sdev接触粘结偏载切向强度1 l0 w, R+ M. ( g3 E2 D* t: 7 ?$ H& N& D1 + o到目前位置，重要的建模参数都在这里了，这些弄清楚了，模型基本能建成功并运行。8 G x- ! T/ * E0 i大家是否注意到这些参数的命名很特别前缀表示该参数所在fish函数的函数名。如et2_xlen在et2.fis文件中调用。+ a0 y2 q# s2 r8 P下面就介绍核心部分，fish函数。
