定义单元体能量为微单元体的应变比能，即当应力和应变满足线性关系时，微单元体在三向应力状态下的应变比能为：
（3.1）
下面代码为用户自定义云图显示变量。
Flac3d Code
1 new 2 gen zone brick size 10 10 10 group lay-1 3 mo el 4 prop bulk 3e8 shear 1e8 5 ini den 20e3 6 set gravity 0 0 -9.8 7 fix z range z 0 8 solve 9 config zextra 2 10 def _udf_zextra 11 p_z = zone_head 12 _u=0.3 13 _E=2.7e8 ;_E=9*bulk*shear/3*bulk+shear) 14 loop while p_z # null 15 if z_groupp_z)='lay-1';只计算在组内的单元体 16 _sig1=z_sig3p_z) 17 _sig2=z_sig2p_z) 18 _sig3=z_sig1p_z) 19 _a=_sig1^2+_sig2^2+_sig3^2 20 _b=2*_u*_sig1*_sig2+_sig2*_sig3+_sig3*_sig1) 21 z_extrap_z,1) =_a-_b)/2*_E);微单元体应变比能 22 z_extrap_z,2) = z_sxxp_z);自定义为单元体sxx应力 23 endif 24 p_z = z_nextp_z) 25 endloop 26 end 27 _udf_zextra 28 plot add axes 29 plot set plane nor 0 1 0 ori 0 5 0 30 plot con zextra 1 plane ;range 5 7,不能加range 31 ;plot add con zextra 2 plane
第二个变量取为sxx只是为了显示在一段代码中如何同时自定义两个变量。多个变量与此类似，不再祥述。
 |
 |
|
a Szz应力云图 a The szz stress cloud |
b 应变比能应力云图 b Strain specific stress cloud |
图3.3 自定义单元应力有关云图
Fig.3.3 The custom cloud map with zone stress
下面是节点变量：
Flac3d Code
1 new 2 gen zone brick size 10 10 10 group lay-1 3 mo el 4 prop bulk 3e8 shear 1e8 5 ini den 20e3 6 set gravity 0 0 -9.8 7 fix z range z 0 8 solve 9 config gpextra 2 10 def _udf_gpextra 11 p_gp = gp_head 12 loop while p_gp # null 13 gp_extrap_gp,1) = sqrtgp_xposp_gp)-xc)^ 2+gp_zposp_gp)-zc)^ 2) 14 gp_extrap_gp,2) = atan2gp_xposp_gp)-xc),zc-gp_zposp_gp))) 15 p_gp = gp_nextp_gp) 16 endloop 17 end 18 set xc=5 zc=5 19 _udf_gpextra 20 plot set background white 21 ;plot block gr 22 plot con gpextra 1 23 plot con gpextra 2
 |
 |
|
a 自定义节点云图1 a custom node cloud 1 |
b 自定义节点云图2 b custom node cloud 2
|
图3.4 自定义节点运动量有关云图
Fig.3.4 The custom cloud map with gridpoint movement
想看详细的.sav和fish代码请关注下面的二维码（flac3d自定义变量输出云图显示）。
