www.pudn.com > flac-in-coal.zip > flac in coal.txt, change:2012-11-24,size:4970b


n  
set logfile coal_41			;log 
title 'sunchun workface model' 		;title 
;create bound------30m 			;创建底层bound组 
macro model_size 'size 50 50'        ;用宏命令  建立 一种网格:  一层网格,size 50 50 1 
gen zone brick model_size 10 p0 0 0 0  p1 142.56 0 0  p2 0 156 0 p3 0 0 30 p4 142.56 156 0 p5 0 156 30 p6 142.56 0 102.546 p7 142.56 156 102.546;用宏的网格建了10层,相当于size 50 50 10 
group bound 				;定义为组:bound是组的名字 
 
inter 1 face range plane dip 27 dd 270 ori 142.56 0 102.546  	;接触面单元在所有在指定的范围里面的格子面上被建立。位置在bound组的上层面建立1号面 
 
;create bottom_4 -----2m 
gen zone brick model_size 1 p0 0 0 30  p1 142.56 0 102.546  p2  0 156 30  p3 0 0 32 p4 142.56 156 102.546 p5 0 156 32  p6 142.56 0 104.546 p7 142.56 156 104.546 
group bottom_4 range group bound not	;除去之前的 bound组,将刚生成的定义为bottom_4组 
inter 2 face range plane dip 27 dd 270 ori 142.56 0 104.546  	;建立2号面,位置在bottom_4组的上层面 
 
 
;create coal_4 ------2m 
gen zone brick model_size 1 p0 0 0 32 p1 142.56 0 104.546 p2 0 156 32 p3 0 0 34 p4 142.56 156 104.546 p5 0 156 34 p6 142.56 0 106.546 p7 142.56 156 106.546  
group coal_4 range group bound not group bottom_4 not		;同上 
inter 3 face range plane dip 27 dd 270 ori 142.56 0 106.546  	;同上 
 
 
;create roof_4 -----15.5m 
gen zone brick model_size 5 p0 0 0 34 p1 142.56 0 106.546 p2 0 156 34 p3 0 0 49.5 p4 142.56 156 106.546 p5 0 156 49.5 p6 142.56 0 122.046 p7 142.56 156 122.046  
group roof_4 range group bound not  group bottom_4 not  group coal_4 not  
inter 4 face range plane dip 27 dd 270 ori 142.56 0 121.546   
 
 
;create roof_3----3m 
 
gen zone brick model_size 1 p0 0 0 49.5  p1 142.56 0 122.046 p2 0 156 49.5 p3 0 0 52.5 p4 142.56 156 122.046 p5 0 156 52.5 p6 142.56 0 125.046 p7 142.56 156 125.046  
group roof_3 range group bound not  group bottom_4 not  group coal_4 not group roof_4 not  
 
inter 5 face range plane dip 27 dd 270 ori 142.56 0 125.046   
 
 
;create coal_2-----3m  
gen zone brick model_size 1  p0 0 0 52.5  p1 142.56 0 125.046 p2 0 156 52.5 p3 0 0 55.5 p4 142.56 156 125.046 p5 0 156 55.5 p6 142.56 0 128.046 p7 142.56 156 128.046  
group coal_2 range group bound not  group bottom_4 not  group coal_4 not group roof_4 not group roof_3 not   
inter 6 face range plane dip 27 dd 270 ori 142.56 0 128.046 
 
 
;create directroof_2 ------8m 
gen zone brick model_size 2  p0 0 0 55.5  p1 142.56 0 128.046 p2 0 156 55.5 p3 0 0 63.5 p4 142.56 156 128.046 p5 0 156 63.5 p6 142.56 0 136.046 p7 142.56 156 136.046  
group directroof_2 range group bound not  group bottom_4 not  group coal_4 not group roof_4 not group roof_3 not    group coal_2 not 
inter 7 face range plane dip 27 dd 270 ori 142.56 0 136.046  
 
 
;create basicroof_2-------30m 
gen zone brick model_size 10 p0 0 0 63.5  p1 142.56 0 136.046 p2 0 156 63.5 p3 0 0 93.5 p4 142.56 156 136.046 p5 0 156 93.5 p6 142.56 0 166.046 p7 142.56 156 166.046  
group basicroof_2 range group bound not  group bottom_4 not  group coal_4 not group roof_4 not  
group roof_3 not    group coal_2 not  group directroof_2 not   
inter 8 face range plane dip 27 dd 270 ori 142.56 0 166.046   
 
 
;create ubound---------30m 
gen zone brick model_size 10 p0 0 0 93.5  p1 142.56 0 166.046 p2 0 156 93.5 p3 0 0 196.046 p4 142.56 156 166.046 p5 0 156 196.046 p6 142.56 0 196.046 p7 142.56 156 196.046  
group ubound range group bound not  group bottom_4 not  group coal_4 not group roof_4 not & 
group roof_3 not    group coal_2 not  group directroof_2 not & 
group basicroof_2 not 
 
model e   						;elastic  isotropic elastic model 弹性模型 
;model m						;mohr   Mohr-Coulomb plasticity model  摩尔库伦模型   前面加了分号,表示不运行。 
  
 
 
prop bu 4.86e9 shear 1.02e9  				;定义属性,体积模量 
 
 
prop den 2000						;定义密度(将整体的密度都简化视为 2000),若参数足够,应按以下定义 
;prop den 1334 range group coal_4   group coal_2   
;prop den 2500 range group bound  group ubound 
;prop den 2591 range group bottom_4 
;prop den 2534 range group roof_4  group roof_3 
;prop den 2662 range group directroof_2  
;prop den 2558 range group basicroof_2  
 
inter 1 pro ks 0.6e9 kn 0.6e9 fri 15			;设置接触面属性:ks(切向刚度),kn(法向刚度)	,fric(摩擦角)	 
inter 2 pro ks 0.6e9 kn 0.6e9 fri 15 
inter 3 pro ks 0.6e9 kn 0.6e9 fri 15 
inter 4 pro ks 0.6e9 kn 0.6e9 fri 15 
inter 5 pro ks 0.6e9 kn 0.6e9 fri 15 
inter 6 pro ks 0.6e9 kn 0.6e9 fri 15 
inter 7 pro ks 0.6e9 kn 0.6e9 fri 15 
inter 8 pro ks 0.6e9 kn 0.6e9 fri 15 
 
 
fix x range x -0.1 0.1					;设置边界约束条件 
fix x range x 142.5 142.6 
fix y range y -0.1 0.1  
fix y range y 155.9 156.1 
fix z range z -0.1 0.1 
 
set grav 0 0 -10					;设置重力加速度 
 
hist unbal 						;设置记录点:unbalance :maximum unbalanced force 
 
 
his gp vel 0 0 196.046 					;记录节点(gridpoint、网格点)vel:velocity magnitude(速度值) 
his gp zdis 0 0 196.046 				;记录该点的  竖向(z轴方向)的位移变化 
 
solve 							;计算 
save coal_42.sav					;保存