复杂边坡(pfc).docVIP

复杂颗粒模型的建立方法，我的理解不外乎两种方式;7 U0 b* r# h5 L% n膨胀法：那层状边坡模型来说，首先将每个岩性分层通过wall对象分割开，然后对每个分层区域看做为独立对象，并这些独立区域进行颗粒充填，待每个分层充填完后，对模型执行迭代，最终使得颗粒膨胀、挤密压实，当然如果再使用锁固力控制颗粒体状态自然再好不过，最后的操作只需要把岩性分层间的wall删除即可； 0 z+ g0 f# C. |投掷法：即逐个对每个分层区域进行颗粒充填，这种方法有个要求，即事先需要假定每个区域的孔隙率，就此作为投掷颗粒的终止准则，一旦孔隙率标准达到后，之后的操作就与第一种方法一致了。, U5我建了一个模型$ t- D/ t, k H7 k R r需要运行一段时间2 w* W* @- ?; E R用cycle XXX来设定运行的步数/ Z7 y% K w [5 e- @# y/ l- S5 P8 m- v/ o但是用set dt??XXX 设定步长时老是出错！！. I4 J/ g6 M4 \7 v3 t: i3 t. S Q* K F9 s% n请问高手怎么操作？？: Q [% R F u8 e2 Q$ n5 @* @7 V. z/////////////////////////////////////////////??F- A9 I# j e$ @0 [; 7 ?$ q6 L; [- g3 F6 L# Hset mechanical on6 a D1 @ _2 q6 h7 B: f% O5 xset hist_rep 5 ;数据采集频率 [: w7 b6 f4 B S w: O4 P??|4 Fset dt 2e-8, ! D3 g, K8 T3 [% R, U B 设置视角的问题。有两种办法可以解决: `??j3 o/ L- V??l; ?7 H# s（1）手动。利用键盘上的上、下、左、右、end、home六个键可以对图形的位置与大小调节。; N??W/ @9 p; k: E4 z（2）程序语句。设置plot set dist 5.0e+1中的数字部分，当参数合适时就可以了。new1 J$ T: A6 i U O/ n. [9 mSET random 6 j4 M* M??U! e J F2 lset gen_error off+ F% X* b5 O# _def expand9 |$ G; z! Y8 Y v6 G# qn_stiff 5e6 ; @, [8 @5 @0 u# L8 x$ Us_stiff 5e6 5 b+ z, v7 K ?, M r5 |??twn_stiff 1e73 I$ E8 u, K2 G- t L1 X! P: e3 n7 iws_stiff 1e7 W4 q6 Z, 4 X3 O4 E1 aporos 0.4764 o; `. t b9 X0 `% g; J+ B. Fwidth 6e-4/ M: x+ Y, C9 x; K, Z s9 Gheight 6e-4 e. v: M$ - h ] T6 cdepth 6e-46 P: G5 o9 h3 S# _ Q n$ _: drlo 5e-65 ??E1 l, f; I??Z1 e zrhi 17e-6+ E2 , W s0 ?; I! w/ \ p% Wtot_vol width*height*depth9 l5 S0 y$ c! G4 Frbar 0.5* rlo + rhi ! C- x$ | w0 Anum int 1.0 - poros *tot_vol*3.0/ 4.0*pi*rbar^3 c2 d- ], j9 Z6 \command, R# C/ V ?/ c; w. K8 [wall id 1 ks ws_stiff kn wn_stiff norm 0,1,0 orig 0,0,0 - h1 |* c??a q5 d3 l4 @ N; H+ `wall id 2 ks ws_stiff kn wn_stiff norm -1,0,0 orig width,0,0 3 E y6 t% R. _$ ?wall id 3 ks ws_stiff kn wn_stiff norm 0,-1,0 orig 0,height,0 * f1 E* [; i+ w- Rwall id 4 ks ws_stiff kn wn_stiff norm 1,0,0 orig 0,0,0 1 K: h a; p9 [ s- Kwall id 5