VECFEM3 Reference Manual: vempat

Type: FORTRAN routine

NAME

vempat - creates a PATRAN neutral file

SYNOPSIS

CALL VEMPAT(: TEXT1, LIVEM, IVEM, LNEK, NEK, LRPARM, RPARM, LIPARM, IPARM, LDNOD, DNOD, LRDPRM, RDPARM, LIDPRM, IDPARM, LNODN, NODNUM, LNOD, NOD, LNOPRM, NOPARM, LBIG, RBIG, IBIG)
INTEGER: LIVEM, LNEK, LRPARM, LIPARM, LDNOD, LRDPRM, LIDPRM, LNODN, LNOPRM, LBIG
CHARACTER*80: TEXT1
INTEGER: IVEM(LIVEM), NEK(LNEK), IPARM(LIPARM), DNOD(LDNOD), IDPARM(LIDPRM), NODNUM(LNODN), IBIG(*)
DOUBLE PRECISION: RPARM(LRPARM), RDPARM(LRDPRM), NOD(LDNOD), NOPARM(LNOPRM), RBIG(LBIG)

PURPOSE

PATRAN is a program package for the pre- and postprocessing of the FEM (in structural analysis). The generation of the FEM mesh is supported by graphical equipment. The solutions can be plotted. The interface to FEM programs is the neutral file.

vempat writes a PATRAN neutral file from the geometrical mesh of the FEM mesh. The elements (except for the nodal element with CLASS=0) must have a unique element id number stored as the first integer vector parameter. The solution is handed over to PATRAN by the nodal result file, which is created by vepa97. Element results (e.g. stresses computed by vemu03) are written to the element result file by vepa99. vempat can be called without a preceding vemdis call.

ARGUMENTS

TEXT1 character*80, scalar, input, local

Title of neutral file.

LIVEM integer, scalar, input, local

Length of the integer information vector, LIVEM>= MESH+ NINFO.

IVEM integer, array: IVEM(LIVEM), input/output, local/global

Integer information vector.

(1)=MESH, input, local

Start address of the mesh informations in IVEM, MESH>203+ NPROC.

(2)=ERR, output, global

Error number.

0	program terminated without error.
90	`LBIG` is too small.
96	unexpected element type.
99	fatal error.

(5)=NIVEM, output, local

Used length of IVEM.

(120)=LOUT, input, local

Unit number of the standard output file, normally 6.

(121)=OUTCNT, input, local

Output control flag, normally 1.

0	only error messages are printed.
>0	a protocol is printed.

(124)=COMP6, input, local

If COMP6=0, the Dirichlet conditions create one constraint set. In the other cases the Dirichlet conditions of component d create one constraint set with constraint set id d.

(125), input, local

Unit of the neutral file.

(200)=NPROC, input, global

Number of processes, see combgn.

(201)=MYPROC, input, local

Logical process id number, see combgn.

(202)=NMSG, input/output, global

Message counter. The difference of the input and the output values gives the number of communications during the vempat call.

(204)=TIDS(1), input, global

Begin of the list TIDS which defines the mapping of the logical process ids to the physical process ids. See combgn.

(MESH), input, local

Start of mesh informations, see mesh.

LNEK integer, scalar, input, local

Length of the element array.

NEK integer, array: NEK(LNEK), input, local

Array of the elements, see mesh.

LRPARM integer, scalar, input, local

Length of the real parameter array.

RPARM double precision, array: RPARM(LRPARM), input, local

Real parameter array, see mesh.

LIPARM integer, scalar, input, local

Length of the integer parameter array.

IPARM integer, array: IPARM(LIPARM), input, local

Integer parameter array, see mesh.

LDNOD integer, scalar, input, local

Length of the array of the Dirichlet nodes.

DNOD integer, array: DNOD(LDNOD), input, local

Array of the Dirichlet nodes, see mesh.

LRDPRM integer, scalar, input, local

Length of the real Dirichlet parameter array.

RDPARM double precision, array: RDPARM(LRDPRM), input, local

Array of the real Dirichlet parameters, see mesh.

LIDPRM integer, scalar, input, local

Length of the integer Dirichlet parameter array.

IDPARM integer, array: IDPARM(LIDPRM), input, local

Array of the integer Dirichlet parameters, see mesh.

LNODN integer, scalar, input, local

Length of the array of the id numbers of the geometrical nodes.

NODNUM integer, array: NODNUM(LNODN), input, local

Array of the id numbers of the geometrical nodes, see mesh.

LNOD integer, scalar, input, local

Length of the array of the coordinates of the geometrical nodes.

NOD double precision, array: NOD(LNOD), input, local

Array of the coordinates of the geometrical nodes, see mesh.

LNOPRM integer, scalar, input, local

Length of the array of the node parameters.

NOPARM double precision, array: NOPARM(LNOPRM), input, local

Array of the node parameters, see mesh.

LBIG integer, scalar, input, local

Length of the real work array. The needed length of LBIG depends on the given mesh. A minimal length of LBIG cannot be given. It should be as large as possible.

RBIG double precision, array: RBIG(LBIG), work array, local

Real work array.

IBIG integer, array: IBIG(*), work array, local

Integer work array, RBIG and IBIG have to be defined by the EQUIVALENCE statement.

EXAMPLE

See vemexamples.

METHOD

vempat writes to the neutral file on the first process only.

The Nodes

If DIM=3, VEMIDE writes the x-,y- and z-coordinates from the array NOD to the PATRAN neutral file (if DIM=2 or DIM=1, missing coordinates are filled by zeroes).

The Elements

The first integer parameter is the element id number (except for the nodal elements (CLASS=0)). It must be unique, i.e. there must exist no elements with the same element number. The elements with CLASS=0 create the nodal forces in a load set. The real vector parameters are written to the components of the forces. If NRVP is greater than 6, only the first six vector parameters are considered. Missing values are set to zero.

Dirichlet Conditions

The Dirichlet conditions generate the displacements in special constraint sets. If COMP6=0, the Dirichlet conditions create one constraint set. If NK is greater than 6, only the first six components are considered. If COMP6<>0, every component d creates one constraint set with constraint set id d, which contains all Dirichlet conditions of the component d. The first real vector parameter defines the value of the displacements. Missing values are set to zero.

REFERENCES

[FAQ], [DATAMAN], [DATAMAN2], [P_MPI], [PATRAN].

COPYRIGHTS

by L. Grosz, Auckland , 11. June, 2000.