qinjianshe 发表于 2004-9-22 11:10:30

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

请教:
关于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 发表于 2004-9-22 17:36:24

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

直接引用 unbal 即可

qinjianshe 发表于 2004-9-22 20:14:30

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

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

qinjianshe 发表于 2004-9-22 20:15:38

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

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

dynamax 发表于 2004-9-22 20:40:02

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

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

qinjianshe 发表于 2004-9-23 09:25:51

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

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

dynamax 发表于 2004-9-23 11:00:30

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

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

qinjianshe 发表于 2004-9-23 19:45:02

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

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

海底冬天 发表于 2010-11-7 15:49:29

页: [1]
查看完整版本: 【求助】请教一个fish语句应用问题