FEAP User Forum

FEAP => nurbFEAP => Topic started by: 1user on June 16, 2021, 08:37:37 AM

Title: Closest point projection to NURBS surface (cont_projection06.f)
Post by: 1user on June 16, 2021, 08:37:37 AM

Dear all,

I am using the 3D IGA contact (celmt06.f) and cannot understand what exactly is happening in the subroutine cont_projection06.f.

My question is - why do we not add geometrical part of the matrix while updating the trial point in the Newton loop unless the contribution is very small.

183         ! Add geometric part if small enough
184         if(jgmax.lt.jjmax*1.d-03) then
185            jj = jj + jg
186         endif

Thank you in advance for any hints.

Best regards