Dear FEAP community,
The OUTMesh command in nurbfeap does not produce the mesh output properly (it lacks correct side information) most of the times. For instance, the side data for the OUTMesh outcome of the attached input file is:
NSIDes 6
side 1 4 1 0 0 0 0
side 2 4 1 0 0 0 0
side 3 4 1 0 0 0 0
side 4 3 2 0 0 0
side 5 3 2 0 0 0
side 6 3 2 0 0 0
which lacks control points data.
The problem appears to be in poutm.f file. I think
include 'cdata.h'
logical setvar,palloc
...
if(np(79).eq.0) then
setvar = palloc( 79,'IPOS ',numnp, 1)
call pseqn(mr(np(79)),numnp)
endif
...
must be added before calling poutsid subroutine in outm_nurbs subroutine.
NSIDes 6
side 1 4 1 1 2 3 4
side 2 4 1 5 6 7 8
side 3 4 1 9 10 11 12
side 4 3 2 1 5 9
side 5 3 2 2 6 10
side 6 3 2 3 7 11
or, the user must use a command which activates the np(79) array (e.g. TIE) before calling OUTMesh.
Kind regards,
RM
