Which subroutine/variable contains information about the type of load?

[serkuz]

Dear Prof. Taylor, dear FEAP community,

For my algorithm one must be able to distinguish between the load types (e.g. FORCe or DISPlacement) and I would like to know where in the program this information is stored.

Best Regards
serkuz

[Prof. R.L. Taylor]

Nodal values of forces FORCe and displacements DISPl are stored in the array "F": f(ndf,numnp,2) --> 1 = force; 2= displacements. 

However, these are not used directly as they need to be multiplied by the proportional load.    You can look in pmacr1.f and follow what happens with the call to "pload".

[FEAP_Admin]

To tell what type of condition is set look at ID(ndf,numnp,2) == mr(np(31)+nneq).

For example, for degree of freedom 3 at node 10, ID(3,10,2) .eq. 0 means a force boundary condition, ID(3,10,2) .ne. 0 means a displacement boundary condition.

[serkuz]

Nodal values of forces FORCe and displacements DISPl are stored in the array "F": f(ndf,numnp,2) --> 1 = force; 2= displacements.

However, these are not used directly as they need to be multiplied by the proportional load. You can look in pmacr1.f and follow what happens with the call to "pload".

Thank you very much Professor Taylor.

[serkuz]

To tell what type of condition is set look at ID(ndf,numnp,2) == mr(np(31)+nneq).

For example, for degree of freedom 3 at node 10, ID(3,10,2) .eq. 0 means a force boundary condition, ID(3,10,2) .ne. 0 means a displacement boundary condition.

Thank you very much, this is exactly what I'm looking for.

I wrote the following code:

Code: [Select]

...
      integer    id(ndf,numnp,2)
...
      id(ndf,numnp,2) = mr(np(31)+nneq)
...
open(unit = ios, file = "test.txt")
        write(ios,*) id(2,1,2)
close(ios)
...
However, at node 1, is the clamping and the code gives me a zero value, so the value should be non-zero.
Have I misunderstood your post or wrongly coded?

Thanks in advance.

Best regards
serkuz

[Prof. S. Govindjee]

You have misunderstood.
mr(np(31)+nneq) is the same location in memory as id(1,1,2).  Thus a normal way to access id would be

Code: [Select]

include 'comblk.h'
include 'sdata.h'
call suba(mr(np(31)+nneq)
then

Code: [Select]

subroutine suba(id)
include 'sdata.h'
integer id(ndf,*)
at this stage inside subroutine suba you have access to ID.  The other way to access ID would be to directly index into MR( ) at location
np(31)+nneq and beyond, which is not advised.

[serkuz]

You have misunderstood.
mr(np(31)+nneq) is the same location in memory as id(1,1,2).  Thus a normal way to access id would be

Code: [Select]

include 'comblk.h'
include 'sdata.h'
call suba(mr(np(31)+nneq)
then

Code: [Select]

subroutine suba(id)
include 'sdata.h'
integer id(ndf,*)
at this stage inside subroutine suba you have access to ID.  The other way to access ID would be to directly index into MR( ) at location
np(31)+nneq and beyond, which is not advised.

Sorry for last answer. Thank you very much Prof. S. Govindjee.