SET md_Ec=62e9 md_knoverks=2.5 SET md_fric=0.50 ;
; Specify parallel-bond parameters. ;
SET md_add_pbonds=1
SET pb_radmult=1.0 pb_Ec=62e9 pb_knoverks=2.5 SET pb_sn_mean=157e6 pb_sn_sdev=36e6 SET pb_ss_mean=157e6 pb_ss_sdev=36e6 ;
; Specify contact-bond parameters. ;
SET md_add_cbonds=0
SET cb_sn_mean=0.0 cb_sn_sdev=0.0 SET cb_ss_mean=0.0 cb_ss_sdev=0.0 ;
=========================================================================== return
;EOF: mA-param.DAT
之后et2_prep命令开始调用上面的md.FIS；et2.FIS；和flt.FIS文件，生成模型。
et2_prep命令是et2.FIS中的函数，也是整个模型生成的总命令。 1. et2_prep首先调用et2_setup函数
et2_setup函数用来创建试验试样的函数 内部分别调用bg_fillbox函数
bg_fillbox：根据mA-param.DAT文件生成最小粒径为et2_rlo和最大粒径为et2_radius_ratio×et2_rlo均匀分布的颗粒
生成的区域为宽为et2_xlen=31.7e-3，高为et2_ylen=63.4e-3的矩形区域 然后调用et2_install_meas_ciecles函数
et2_install_meas_ciecles：生成三个上中下meascircles。（如下图所示）
试验及测量圆
2. et2_prep接着调用et2_plot_assembly函数（根据mA-param.DAT文件，以后省略）
et2_plot_assembly：建立所要显示的内容（球、墙和坐标如上图所示）。
3. et2_prep接着调用md_save_state函数
md_save_state：保存现阶段的所获结果 4. et2_prep接着调用tm_assemble_isopack函数
tm_assemble_isopack：主要用来使试样获得四周相同的压力，即消除试样之间的不平衡力
tm_assemble_isopack分别一次调用一下内部函数：
md_balldeform：根据md_Ec=62e9 md_knoverks=2.5定义球的刚度（切向刚度和法向刚度）
法向刚度：kn = 2.0 * md_Ec * md2_thick 切向刚度：ks = _kn / md_knoverks md_wallkn： 根据md_wEcfac=1.1定义墙的刚度。
即为求得刚度的1.1倍。（此只设置了墙的法向刚度，并未设置切向刚度和摩擦系数）
之后cyc消除不平衡力。
5. et2_prep接着调用md_save_state函数（md_save_state：保存现阶段的所获结果）
6. et2_prep接着调用tm_assemble_isostr函数
tm_assemble_isostr函数为了获得说要达到的初始应力调用md_expand_radii函数，通过适当的改变球的半径，满足abs(_diso/tm_req_isostr)
<= tm_req_isostr_tol的要求，即为所要求达到的初始的应力状态。
7 et2_prep接着调用flt_eliminate函数
flt_eliminate函数就是消除悬浮颗粒（及颗粒的接触数小于给定的flt_def则视为悬浮颗粒），由于说明起来较为复杂，在此省略。
8. et2_prep接着调用md_pbprops函数
md_pbprops函数：设置parallel-bond参数
9. et2_prep最后调用md_ballfric函数
md_ballfric函数设置球的摩擦系数
则et2_prep函数运行完毕，试样也建好。
下面讲讲 call sCx_mA_tAy-bt.DVR文件，文件的主要内容是 _btw.DVR的调用。 _btw.DVR源文件为： ; Filename: _btw.DVR ;
; PURPOSE: Perform biaxial test using existing four walls and ; extract responses upon test completion. If et2_ucs=1, ; then move sidewalls away from the specimen and lock in place ; to perform fully unconfined test. ;
; INPUT: 'md_run_name' 'bnd' bonded assembly
; OUTPUT: 'md_run_name' 'bw0' stage-0 during testing, start of test ; 'md_run_name' 'bw1' stage-1, test complete ;
; INPUT(Fish): et2_knxfac, et2_knyfac
; et2_ucs, et2_wsxx_req, et2_wsyy_req, et2_ws_tol ; p_vel, p_cyc, p_stages
; md_granular, et2_peakfac, et2_wexx_lim ; pk_ci_fac ;
; Itasca Consulting Group, Inc. ;
===================================================================
SET echo off ; load support functions
call %itascaFishTank%\\FishPfc\\md\\fishcall.FIS call %itascaFishTank%\\FishPfc\\md\\crk.FIS SET echo on ;
=================================================================== ; Specify in calling routine: ; SET et2_knxfac= et2_knyfac=
; SET et2_ucs= et2_wsxx_req= et2_wsyy_req= et2_ws_tol= p_vel=
