3dec中怎样统计地下工程开挖时不稳定块体体积、数量等信息？

t1996

想统计每次开挖后的不稳定块体体积、数量，以及失稳块体的体积、数量等信息，向高手请教

[ 本帖最后由 t1996 于 2007-6-14 10:13 编辑 ]

oldwell7959

统计破裂后 块体 数量 体积等等
是否需要编制语言

abc8892

seek
;set echo off
def _Mark ;Modify
  ;
  ;token=1---check displacement;
  ;token=2---check stress;
  ;token=3---check velocity;
  
  _disp_e1=0.100 ;units m
  _disp_e2=0.075
  _disp_e3=0.050
  _disp_e4=0.025
  _disp_e5=0.010
  
  bi=block_head
   loop while bi # 0
         zi=b_zone(bi)
       ;;;For each zone of every block, find out average displacement
          loop while zi # 0
            z_state(zi)=0  ; Green
                           _x1=0.0
               _y1=0.0
               _z1=0.0
                  _disp1=0.0
            
             loop _j (1,4)
                gi=z_gp(zi,_j)
                     _x1=gp_xdis(gi)
               _y1=gp_ydis(gi)
               _z1=gp_zdis(gi)  
                           _aa=sqrt(_x1*_x1+_y1*_y1+_z1*_z1)
               
               if _aa > _disp1 then
                            _disp1 = _aa
                     end_if
            
             end_loop
         
              if _disp1 >= _disp_e1 then
                z_state(zi)=1   ;
                 ;command
                         ;mark block bi reg 10001
                         ;change block bi mat 1
                       ;end_command
              end_if
         
              if _disp1 >= _disp_e2 then
                if _disp1 < _disp_e1 then
                 z_state(zi)=2  ;
                  ;command
                          ;mark block bi reg 10002
                          ;change block bi mat 2
                        ;end_command
               end_if
              end_if            
      
              if _disp1 >= _disp_e3 then
                if _disp1 < _disp_e2 then
                  z_state(zi)=4  ;
                   ;command
                           ;mark block bi reg 10003
                           ;change block bi mat 3
                         ;end_command
                end_if
              end_if

              if _disp1 >= _disp_e4 then
                if _disp1 < _disp_e3 then
                  z_state(zi)=8  ;
                   ;command
                           ;mark block bi reg 10004
                           ;change block bi mat 4
                         ;end_command
                end_if
              end_if   
            
              if _disp1 >= _disp_e5 then
                if _disp1 < _disp_e4 then
                  z_state(zi)=16  ;
                   ;command
                           ;mark block bi reg 10005
                           ;change block bi mat 5
                         ;end_command
                end_if
              end_if                 
            
              if _disp1 < _disp_e5 then
                  z_state(zi)=32  ;
                   ;command
                           ;mark block bi reg 10006
                           ;change block bi mat 6
                         ;end_command
              end_if
                              


             zi=z_next(zi)      
        end_loop

       bi=b_next(bi)         
  endloop
end

_Mark

;set echo on

;seek
;pl surf pl

abc8892

请参照上面格式 FISH 修改 即可

t1996

拜读中，感谢

冬之恋4916

请教楼主这个问题怎么解决的，谢谢告知