3DEC统计剪切和拉伸接触的数量

model restore 'sample.sav'fish define state_statistics    local elastic_count = 0    local slip_n_count = 0 ;1    local tension_n_count = 0 ;2    local slip_p_count = 0 ;4    local tension_p_count = 0 ;8    local slip_n_p_count = 0 ;5     loop foreach local cx block.subcontact.list        b_state = block.subcontact.state(cx)        if b_state = 0 ; none            elastic_count + = 1        else if b_state = 1 ;Failure in shear now            slip_n_count = slip_n_count + 1        else if b_state = 2 ;Failure in tension now             tension_n_count = tension_n_count + 1        else if b_state = 4 ;Failure in shear in the past             slip_p_count = slip_p_count + 1        else if b_state = 8 ;Failure in tension now             tension_p_count = tension_p_count + 1        else if b_state = 5 ;slip_n, slip_p             slip_n_p_count = slip_n_p_count + 1        endif        endloop    io.out('none = '+ string(elastic_count))    io.out('slip-n = '+ string(slip_n_count))    io.out('tension-n = '+ string(tension_n_count))    io.out('slip-p = '+ string(slip_p_count))    io.out('tension-p = '+ string(tension_p_count))    io.out('slip-n, slip-p = '+ string(slip_n_p_count))end[state_statistics]