www.pudn.com > dichouhangkaiwa.rar > dichouhangkaiwa.txt, change:2013-11-04,size:5140b


restore ʼƽ.sav 
 
plot con szz 
 
; 
 
;ɰҳ 
  
;prop bulk=1.35e9 shear=0.71e9 c=0.87e6 f=19 ten 0.27e6  range group ɰҳ3                     
;ini density 2492 range  group ɰҳ3 
 
;prop bulk=1.35e9 shear=0.71e9 c=0.87e6 f=19 ten 0.27e6  range group 1                     
;ini density 2492 range  group 1 
 
;prop bulk=1.35e9 shear=0.71e9 c=0.87e6 f=19 ten 0.27e6  range group 2                     
;ini density 2492 range  group 2 
 
 
;prop bulk=2.0e9 shear=1.05e9 c=1.67e6 f=25 ten 0.47e6  range group ɰҳ                   
;ini density 2492 range  group ɰҳ 
 
 
 
;ɰ 
;prop bulk=0.47e9 shear=0.21e9 c=2.9e6 f=26 ten 0.82e6 range  group ɰ 
;ini density 2644 range  group ɰ 
 
 
 
 ;Ұú 
 
; prop bulk=0.23e9 shear=0.14e9 c=0.58e6 f=17 ten 0.08e6 range  group Ұú 
; ini density 1350 range  group Ұú 
 
 
;Ұ 
  
; prop bulk=4.08e9 shear=1.63e9 c=4.1e6 f=33 ten 1.87e6 range  group Ұ 
; ini density 2716  range group Ұ 
 
 
;ɰ 
;  prop bulk=2.60e9 shear=1.04e9 c=2.6e6 f=31 ten 0.67e6 range group ɰ 
;  ini density 2597 range  group ɰ 
 
 
 ;2ú 
 
; prop bulk=0.23e9 shear=0.14e9 c=0.58e6 f=17 ten 0.08e6 range  group 2ú 
; ini density 1350 range  group 2ú 
 
 
 ;ϸɰ 
 
; prop bulk=2.48e9 shear=1.05e9 c=2.75e6 f=29 ten 0.71e6 range  group ϸɰ 
; ini density 2613 range  group ϸɰ 
 
 
;ϸɰ 
 
; prop bulk=2.48e9 shear=1.05e9 c=2.8e6 f=28 ten 0.78e6 range  group ϸɰ 
; ini density 2628 range  group ϸɰ 
 
 
; prop bulk=1.72e9 shear=0.73e9 c=2.0e6 f=22 ten 0.48e6 range  group ϸɰ1 
; ini density 2628 range  group ϸɰ1 
 
 
 
;̿ҳ 
  
; prop bulk=0.73e9 shear=0.37e9 c=1.18e6 f=21 ten 0.21e6 range  group ̿ҳ 
; ini density 2451 range  group ̿ҳ  
 
 
;ϵͳλٶ 
 
ini xdis=0 ydis=0 zdis=0 xvel=0 yvel=0 zvel=0 
   
model null range group 1 y 50 60 
 
set large 
 
;step 50 
;save 50.sav 
 
 
def  ins_cab   
 
    y_t=0 
  
    loop t(0,9) 
            y_t=50.5+t*1 
;======================================================================= 
;ê2.4 ê6/5 15 54 
;=======================================================================    
	command 
 
	    ;ê 
	    sel cable id 1 begin 139.99 y_t 48.4  end 137.67     y_t  47.78    nseg=6 ;ê 
            sel cable id 1 begin 139.99 y_t 49.3  end 137.59     y_t  49.3     nseg=6 ;ê 
             
	    ;Ұê 
	    sel cable id 1 begin 144.01 y_t 48.4  end 146.33     y_t  47.78    nseg=6 ;Ұê 
	    sel cable id 1 begin 144.01 y_t 49.3  end 146.41     y_t  49.3     nseg=6 ;Ұê 
	     
	    ;ê 
	    sel cable id 1 begin 140.268   y_t 50.31   end 138.1896   y_t 51.5      nseg=6  ;һê 
            sel cable id 1 begin 141       y_t 51.042  end 139.8      y_t 53.1104   nseg=6  ;ڶê 
            sel cable id 1 begin 142       y_t 51.31   end 142,01     y_t 53.7      nseg=6  ;ê 
            sel cable id 1 begin 143       y_t 51.042  end 144.2      y_t 53.1104   nseg=6  ;ĸê 
            sel cable id 1 begin 143.732   y_t 50.31   end 145.8104   y_t 51.5      nseg=6  ;ê 
 
 
	end_command 
     end_loop 
 
    loop t(0,4) 
          
	   y_t=t*2+51 
 
	command 
 
	    ;ê 
	    sel cable id 2 begin 141.01  y_t 51.032  end 138.01   y_t 56.228  nseg=10 ;ê 
	    sel cable id 2 begin 143.01  y_t 51.032  end 146.01   y_t 56.228  nseg=10 ;ê 
 
	    sel cable id 2 begin 139.99   y_t 49.31    end 134.99  y_t 49.31  nseg=10 ;ê 
	    sel cable id 2 begin 144.01   y_t 49.31    end 149.01  y_t 49.31  nseg=10 ;Ұê 
 
	end_command 
 
     end_loop 
end 
ins_cab 
 
;êֱ22êֱ21.8ê̺8mm 
 
sel cable id 1 prop emod=200e9   ytens=190e7 xcarea=1.3e-3 gr_k=5.35e9 gr_coh=4.2e5 gr_per=0.22  gr_fric=25    
sel cable id 2 prop emod=200e9   ytens=190e7 xcarea=1.3e-3 gr_k=5.35e9 gr_coh=4.2e5 gr_per=0.218 gr_fric=25 
;============================================================= 
 
;pause;********************************************************************************************************** 
 
;êlink 
def camd 
    new_id=100000000 
    loop k(0,89) 
	id_node=1+k*7 
	new_id=new_id+1 
	command 
	sel delete link range id id_node 
	sel link id=new_id id_node target zone 
	sel link attach xdir=rigid ydir=rigid zdir=rigid xrdir=rigid yrdir=rigid zrdir=rigid range id new_id 
	sel cable id=1 pretension 60e3  
	end_command 
    end_loop 
end 
camd 
;=================================================================================================================  
;êlink 
def cdma 
    new_id=200000000 
    loop k(0,39) 
	id_node=631+k*11 
	new_id=new_id+1 
	command 
	sel delete link range id id_node    
	sel link id=new_id id_node target zone 
	sel link attach xdir=rigid ydir=rigid zdir=rigid xrdir=rigid yrdir=rigid zrdir=rigid range id new_id 
	sel cable id=2 pretension 60e3 
	end_command 
    end_loop 
end 
cdma 
 
 
;plot con xdisp 
 
plot con zdisp 
pause 
 
 
set large 
solve  
save ê˵׳↑110m.sav