Obtaining Node to Element Connectivity

Utility routine GETNODETOELEMCONN can be called from user subroutine UMESHMOTION to retrieve a list of elements connected to a specified node.

This page discusses:

 

See Also
UMESHMOTION
Obtaining Material Point Information Averaged at a Node
Erosion of material (sand production) in an oil wellbore

Utility Routine Interface

PARAMETER ( MAXNELEMS = 100 )
DIMENSION JELEMLIST(MAXNELEMS),JELEMTYPE(MAXNELEMS),JGVBLOCK(*)
...
NELEMS = MAXNELEMS
CALL GETNODETOELEMCONN(NODE, NELEMS, JELEMLIST, JELEMTYPE,
     JRCD, JGVBLOCK)

Variables to Be Provided to the Utility Routine

NODE

User node number.

NELEMS

You must set NELEMS to the maximum allowable length of the JELEMLIST and JELEMTYPE arrays. This value corresponds to the maximum expected number of elements attached to an adaptive mesh constraint node in your model. GETNODETOELEMCONN will assume that your JELEMLIST and JELEMTYPE arrays are NELEMS long. In the event that the actual element connectivity exceeds NELEMS, no result will be returned and the return code JRCD will indicate an error. An NELEMS value of 100 is typically more than adequate for common meshes. NELEMS is modified by GETNODETOELEMCONN and should not be a Fortran parameter-statement constant.

JGVBLOCK

Variable that must be passed into the GETNODETOELEMCONN utility routine. This variable is available in user subroutine UMESHMOTION for this purpose.

Variables Returned from the Utility Routine

JELEMLIST

Array of element numbers for elements connected to NODE. The list will contain elements only in adaptive mesh domains active in the step as well as any contact elements associated with the domain. The number of entries in this array corresponds to the returned value of NELEMS.

JELEMTYPE

Array of element type designators describing the element types corresponding to each element entry in JELEMLIST. The number of entries in this array corresponds to the returned value of NELEMS.

JELEMTYPE entries:

1 indicates a solid element.

2 indicates a contact element.

NELEMS

Actual length of the JELEMLIST and JELEMTYPE arrays.

JRCD

Return code (0 indicates no error, 1 indicates an output request error). An output request error indicates either that the requested variable is not available or that your NELEMS parameter setting is smaller than the element connectivity list at this node.