求助：udf提取力

mirror830

我模拟流固耦合时，边界按正弦方式运动：1、使用F_P提取力，结果力不断增加，到1000以上，并且继续增加。2、用Compute_Force_And_Moment(d,t,CG,F,M,FALSE)提取力发现，x与y方向的力是一样的值。不只是什么原因，希望指教。
另：d=THREAD_DOMAIN(DT_THREAD((Dynamic_Thread*)dt));
是不是可以写成
Domain*domain=Get_Domain(1);
Thread *t=Lookup_Thread(domain,zone_ID)
zone_ID可以从哪里得到？boundary condition还是从gird图上显示呢？
感谢！！

tianwen2009

可以写成，zone_ID从boundary condition可获得，使用F_P提取力只是pressure force，为什么增加和x与y方向的力是一样的值，这个没见具体UDF，不好说。

mirror830

找到错误，是自己程序的问题。谢谢。
现在的问题是，边界使用new-mark法运动，但力不断增加，增的很大。相应udf程序我发到我其他帖子上了。

mirror830

http://forum.simwe.com/thread-958624-1-1.html