【求助】请教一个fish语句应用问题

qinjianshe

请教：
  关于fish语句获得计算过程中不平衡力的问题，应用plot hist id=1 unbal可以显示系统不平衡力。但是在fish语句中如何获得该值？
本人用如下方法：
  
V_unbal=0
p_gp = gp_head
loop while p_gp # null
V_unbal=V_unbal+abs(gp_xfunbal(p_gp))+abs(gp_yfunbal(p_gp))+abs(gp_zfunbal(p_gp))
p_gp = gp_next(p_gp)
end_loop
  
我用节点三个方向不平衡力的绝对值总和计算，不知到对否？
通过上述计算我获得的不平衡力与hist unbal显示的结果不一致，是否与外荷载或其它因素有关，请熟悉fish的大侠指教。

dynamax

直接引用 unbal 即可

qinjianshe

还请dynamax进一步指点：
  如何通过fish语句实现，引用unbal函数怎么写。
多谢

qinjianshe

我想获得该计算时间步中的unbal值。

dynamax

def aaa
whilestepping
ii = out(string(unbal))
end

qinjianshe

谢谢dynamax的帮助，深表感谢  
     计算的不平衡力大小与所有结点不平衡力总和相差较大，问题在那儿还不态清楚
  另外：clock变量指的是flac开始计算的时间吧（但是好像比我的计算机计时要快些），如何得到计算机时间，该用那个函数，请指教。

dynamax

unbal好像不是一个总和
而是某个节点的某个方向上的正不平衡力
我也不知道flac3d中具体是怎么取
clock is number of hundredths-of-a-second from midnight
好像没有相关函数
根据clock自己算或许可以
但得不到绝对的时间
计计时还勉强

qinjianshe

unbal与单元体力、面力及外力等有关，很复杂，只好等等有时间再研究了。
clock is number of hundredths-of-a-second from midnight，好像不是吧，我试过好像是程序开始运行其计时。