fishtank--巴西实验与双轴压缩实验 说明.doc

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为例。这个文件看起来很怪，其实很简单。可以用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- Sagc-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. k3 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_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, _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函数外，就是参数设置，正确理解这些参数就可以任意调整实验以达到生成自己特殊模型实验的目的。在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 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到目前位置，重要的建模参数都在这里了，这些