Knot insertion in order to split patch

[alina.i]

Dear FEAP Team,

I was testing knot insertion command in order to split my patch into separated patches (obtain C-1 continuity). So, for second-order functions, I insert knots 3 times. However, after the last insertion, I obtain wrong values of coordinates and weights.  Is there a way to solve this problem? The body is small, see input and two meshes attached.

I am using Version 8.4. License belongs to TU Dresden.
I would really appreciate any hint.

Best regards,
Alina Israfilova

[Prof. R.L. Taylor]

The knot insertion algorithm only makes sense to C_0.  After that there are no control points at the ends.  Once you are at C_0, to split the patch would be done the same way you would split standard C_0 isoparametric elements.  To understand this think of what you would have for control points for the knot vectors (0,1) and (0,0,1,1).

[alina.i]

I see. Thank you for your help, Professor!

Best regards,
Alina Israfilova

[Yaakov]

 In pre-process, knot insertion ('not local') can be done very easily in NURBS Toolbox ('h,p,k', fully 'free', Octave), which is still being developed by Prof. Alessandro Reali's group at the University of  Pavia.

https://octave.sourceforge.io/nurbs/

If the code works in Octave, the code could also work in Matlab. ('openoffice <-> microsoft office')

But NURBS Toolbox is only used for pre-processing, the main job (99,999999%) must be done in FEAP85/84.

[Yaakov]

@Alina Israfilova, compare with IgA in FEAP84, IgA in FEAP85 includes many new features, more robust, much more fast, etc.

[alina.i]

In pre-process, knot insertion ('not local') can be done very easily in NURBS Toolbox ('h,p,k', fully 'free', Octave), which is still being developed by Prof. Alessandro Reali's group at the University of  Pavia.

https://octave.sourceforge.io/nurbs/

If the code works in Octave, the code could also work in Matlab. ('openoffice <-> microsoft office')

But NURBS Toolbox is only used for pre-processing, the main job (99,999999%) must be done in FEAP85/84.

Thank you!! I will take a look!

[Yaakov]

You are welcome!

NURBS Toolbox contains very limit number of NURBS GEOs. But for h-p-k global refinement, NURBS Toolbox can be a very good assist for pre-processing in  IGA85/FEAP85. 

P.S. if you need k refinement in  NURBS Toolbox , first step:  p refinement,  then next step h refinement.

In pre-process, knot insertion ('not local') can be done very easily in NURBS Toolbox ('h,p,k', fully 'free', Octave), which is still being developed by Prof. Alessandro Reali's group at the University of  Pavia.

https://octave.sourceforge.io/nurbs/

If the code works in Octave, the code could also work in Matlab. ('openoffice <-> microsoft office')

But NURBS Toolbox is only used for pre-processing, the main job (99,999999%) must be done in FEAP85/84.

Thank you!! I will take a look!