Hello. I just reminded that I made a mistake.
Stephane Del Pino a écrit :
Yes. The force on the object can be approximated by the integral over the
object of 1/epsilon*U :
double fx = int[M](1/epsilon*ux);
double fy = int[M](1/epsilon*uy);
double fz = int[M](1/epsilon*uz);
where U = (ux,uy,uz) and the epsilon is used for the penalty term P in the
example for Dominique, 1/epsilon = 10^3.
One should have read
double fx = int[M](ki/epsilon*ux);
where ki is a function that is 0 outside the obstacle and 1 inside. For
instance:
function ki = one(<1,0,0>);
where <1,0,0> is the reference of the object.
Regards,
Stephane.
_______________________________________________
ff3d-users mailing list
[email protected]
http://lists.nongnu.org/mailman/listinfo/ff3d-users
