对于一个简化成二维的圆柱流体，如何通过UDF计算其底面受力

求是创新

如题，在三维中可以通过F_area函数提取底面压力向量，简化为二维的时候该宏经测试计算的就不是一个圆底面的压力了，三维中UDF如下，请问变为二维后应该怎么改呢？
begin_f_loop(f,t)
      {
                  if (PRINCIPAL_FACE_P(f,t))
         {
                         F_AREA(A,f,t);
             force += F_P(f,t) * NV_DOT(A,vector);
                  }
      }
    end_f_loop(f,t)