- 积分
- 0
- 注册时间
- 2012-5-20
- 仿真币
-
- 最后登录
- 1970-1-1
|
发表于 2014-2-15 10:56:59
|
显示全部楼层
来自 武汉大学
这里有一个例子,仅供参考
new
gen zone brick size 1 1 1
...
def _locations
_gp1 = gp_near(0,0,0)
_gp2 = gp_near(0,1,0)
_gp3 = gp_near(1,0,0)
_gp4 = gp_near(1,1,0)
_iaz1 = z_near(0.5,0.5,0.5)
end
@_locations
def _record_variables
;
_disp_0 = 0.5*(gp_xdis(_gp1) + gp_xdis(_gp2))
_disp_1 = 0.5*(gp_xdis(_gp3) + gp_xdis(_gp4))
_eps_xx = -(_disp_0 - _disp_1)/1.0
;
_disp_0 = 0.5*(gp_ydis(_gp1) + gp_ydis(_gp3))
_disp_1 = 0.5*(gp_ydis(_gp2) + gp_ydis(_gp4))
_eps_yy = -(_disp_0 - _disp_1)/1.0
;
_sig_zz = z_szz(_iaz1)
_sig_xx = z_sxx(_iaz1)
_sig_yy = z_syy(_iaz1)
;
_record_variables = 1.0
;
END...
hist @_record_variables
hist @_eps_xx
hist @_eps_yy
hist @_sig_zz
hist @_sig_xx
hist @_sig_yy
...
|
|