不同插值方式下的单元变量值分析

;modeling and assigning groupsmodel newzone create brick size  6 6 6 group 'brick'zone group 'hist' range position-x 0 1  position-y 0 1  position-z 5 6zone face skin 
;assigning cmodel and propertieszone cmodel assign elasticzone prop bulk 1e9 shear 0.5e9 dens 1800
;setting boundary conditionszone face apply velocity-normal 0 range group 'West' or 'East'zone face apply velocity-normal 0 range group 'South' or 'North'zone face apply velocity-normal 0 range group 'Bottom'
;setting history listshistory interval 1zone history name = 'def' stress-zz zoneid 181 zone history name = 'con' stress-zz zoneid 181 method constantzone history name = 'ave' stress-zz zoneid 181 method averagezone history name = 'idw' stress-zz zoneid 181 method i-d-wzone history name = 'poly' stress-zz zoneid 181 method polynomial
;solving initial equilibriummodel gravity 9.81model solve
;setting plot rangezone hide range group 'brick'
;getting hist valueshist list series 'def' begin 260hist list series 'con' begin 260hist list series 'ave' begin 260hist list series 'idw' begin 260hist list series 'poly' begin 260

;getting values of stress-zzfish def z_szz    z = zone.find(181)    pos = zone.pos(z)    zone.field.name = 'stress-zz'    zone.field.method.name = 'constant'    szz_con = zone.field.get(pos)    zone.field.method.name = 'average'    szz_ave = zone.field.get(pos)    zone.field.method.name = 'i-d-w'    szz_idw = zone.field.get(pos)    zone.field.method.name = 'polynomial'    szz_poly = zone.field.get(pos)    szz_default = zone.stress.zz(z)end@z_szzlist @szz_conlist @szz_avelist @szz_idwlist @szz_polylist @szz_default