[etedalgara]

Hi, I need to calculate Force on the blade,
look at the picture:



this is my code but I think this is wrong!
Thanks for your helps

begin_f_loop(f, t)
{
F_CENTROID(x,f,t);
yy = x[0];
F_AREA(area,f,t);
Fpy+=area[1]*F_P(f,t)*1-cos((PI/2)*(yy/L));
}
end_f_loop(f, t)
Fy = Fpy;

[AlexanderZ]

try this code

begin_f_loop(f, t)
{
F_CENTROID(x,f,t);
yy = x[0];
F_AREA(area,f,t);
Fpy+=NV_MAG(area)*F_P(f,t)*1-cos((PI/2)*(yy/L));
}
end_f_loop(f, t)
# if RP_NODE /* Perform node synchronized actions here
Does nothing in Serial */
Fy = PRF_GRSUM1(Fpy);
# endif /* RP_NODE */

you may add
Message0("Fy = %f \n", Fy);
to check the value

[vinerm]

Though the actual equation may be invariant under reflection but the code is not. Therefore, you may have to transform your code if the origin of the coordinate system is not at the fixed end of your beam and also if the free end is not in the positive x. If those conditions are met, then the result would be good. is already defined in Fluent as M_PI.