HELP! This list is far from complete. Some items have no description and others are missing all together. Please send us an email if you can provide info.

Here is the NU table.

Be sure to check the Notes and Obsolete sections at the end.

A

B

C

D

E

F

G

H

I

J

K

L

M

N

O

P

Q

R

S

T

U

V

W

X

Y

Z

UNKNOWN!

A(ns,ns)

?

B(ns,ns) ?

COVA(NEMX,NRMX)

used in data/graphics operations?

EIGV(NOMX,NTMX)

?

ER(nv)

?

MXI(2,NEMX)

used for graphics?

NLCHOR(0:10,NRMX)

used for BEM problems?

NPB(1..6,0)

?

NPB(1..6,np)

?

NTCOVA(NEMX)

used in data/graphics operations?

XA(na,nj,nq)

nodes.global_param_XA

XF(ns,nj)

?

Y

?

ZA(na,nh,ne)

nodes.deform_global_param_ZA

ZC(nj,ne)

used in chnode.f, drelem.f, and moelem.f

ZD2(nj,nd)

used in write.f calcnods.f

ZDD(nd,nj)

data.zdd?

ZF(ns,nh)

?



A


ACOEFF(nactv)

nactv=1,NTACTV are coefficients for linear dynamic terms


ALFA(nactv)

nactv=1,NTACTV are time constants for linear dynamic terms


AXL,AXU
AZL,AZU

are deformed metric tensor components wrt undeformed Nu-coords (KTYP53>1) or wrt undeformed theta coords (KTYP53=1) in zere50.f. SeeRI1,RI2,RI3.

B

colloc.bezier_control_pnts

BCP(4,4,NJMX,NEMX)

Bezier control points for rendering [collocation]

C


C(NYMX,1)

for BEM problems (solve4.f)


CC(NYMX,1)

for BEM problems (solve4.f)


CCOMP(NCYMX,NCYMX)

for BEM problems (solve4.f)

prob.const_mat_param

CE(nm,ne)

is value of piecewise constant material parameters nm in element ne

prob.gauss_mat_param

CG(nm,ng)

is value of material parameter nm at Gauss point ng


CGA(3,ng,nb)

accumulated coordinate transformation matrices [collocation]


CGB(3,3,ng,nb)

accumulated coordinate transformation matrices [collocation]

colloc.integr_scaling_factor

CGDA(ng,nb,ne)

integration scaling factor at the gauss points of basis nb [collocation]

colloc.mat_param

CGMAT(5,ng,nb,ne)

material parameters [collocation]

data.opt_coefs

COYY(noy,ny)

coefficient for optimization DOF NOYY(noy,ny)

prob.linear_mat_const

CP(nm,np)

is value of piecewise linear material constant at node np


CSEG(nosg)

is string defining segment number NOSG as follows:
CSEG(1:4)='NONO' and (5:8) is node number
CSEG(1:4)='ELNO' and (5:8) is element number
CSEG(1:4)='LINE' and (5:8) is line number
CSEG(1:4)='LINO' and (5:8) is line number
CSEG(1:4)='FIBR' and (5:8) is element number
CSEG(1:4)='MATL' and (5:8) is element number
CSEG(1:4)='CPT1' and (5:8) is line number
CSEG(1:4)='CPT2' and (5:8) is line number
CSEG(1:4)='DATA' and (5:8) is data point number
CSEG(1:4)='DANO' and (5:8) is data point number
CSEG(1:4)='CONT' and (5:8) is element number
CSEG(1:4)='CONO' and (5:8) is element number

D


DATYPE(220)

character string of data type drawn in graphics window (set in fe24/drdata.f). Is this ever used anywhere?


DEL_T

is the time step used (taken from load stepping loop in FE07)

elem.face_area

DF

is the logarithmic frequency increment in LPT


DF(nf)

is area of face segment nf

nodes.arc_length_deriv

DL(1..3,nl)

are arc-length derivatives of line segment nl

nodes.old_scaling_factors

DLL(3,nl)

holds old scaling factors DL when rescaling nodes (fe21/upnode.f)


DRDN(NGMX)

for BEM problems (solve4.f)


DRDNS(NGMX)

for BEM problems (solve4.f)


DXIX

are partial derivatives of Xi wrt Xj coords (if KTYP53=1) or wrt Nu coords (if KTYP53>1)


DZDX

are components of the deformation gradient tensor

E

colloc.elem_damping
fit.elem_damping
march1.elem_damping
solve1.elem_damping

ED(nv,nv)

is element damping matrix (1st order time derivatives)

data.diff_XID_ZD

EDD(nd)

difference between element point at XID and data point ZD, weighted by WD (see fe03/fitfld.f)


EG

are physical components of Green's strain

march1.elem_mass
solve1.elem_mass

EM(nv,nv)

is element mass matrix (2nd order time derivatives)

fit.er

ER

is element right hand side; computed in fe03/zder.f


ERRMF

is solution tolerance (epsilon) for the iterative solver

colloc.elem_stiffness
fit.elem_stiffness
march1.elem_stiffness
nonlin.elem_stiffness
solve1.elem_stiffness
solve5.elem_stiffness

ES(nv,nv)

is element stiffness matrix (0th order time derivatives)


ESE(nv,nv,ne)

is element stiffness matrix for element ne (replaces ES when using iterative solver)


ETYP(ie)

ie=1,12 is .TRUE. if element type ie (see TITLE2) used.


EV1RE,EV1IM
EV2RE,EV2IM

are the minimum (EV1*) and maximum (EV2*) Real and Imaginary parts of eigenvalues for the constraint-reduced global system of equations (used in some iterative solvers)


EXR

are extension ratios wrt COORDS

F

prob.fiber_exten

FEXT(8,NGMX,NEMX)

from zgtg5a.f:

FEXT(1)

is current fibre extension ratio

FEXT(2)

is previous fibre extension ratio

FEXT(3)

is fibre extension ratio at the onset of activation

FEXT(4)

is previous value of the hereditary integral

FEXT(5)

is zero when gauss point is inactive, else 1


FLOW_COEFFS(no_coeffs)

no_coeffs=1,nt_coeffs are Fourier coeffs for flow

data.dof_data

FIX(ny,5)

FIX(ny,1)FIX(ny,2)FIX(ny,3)



is used in the assembly of the global system of equations to determine which DOFs are contrained.



defines prescribed force boundary conditions.



is used by subroutine NONLIN to know which DOFs are incremented at each load step.

FIX(ny,1..3) are initialized in ipini5.f

prob.applied_pi

FIXP(2,ne)

is TRUE if a pressure incrementPEis applied to either the Xi_3=0 (1) or Xi_3=1 (2) face of element ne

G

colloc.global_damping
march1.global_damping
solve1.global_damping

GD(ny,ny)

is global damping matrix (1st order time derivatives)

colloc.global_stiffness
fit.global_stiffness
march1.global_stiffness
solve1.global_stiffness

GK(ny,ny)

is global stiffness matrix (0th order time derivatives)


GKC(NCYMX,NCYMX)

for BEM problems (solve4.f)

colloc.reduced_stiffness
fit.reduced_stiffness
march1.reduced_stiffness
nonlin.reduced_stiffness
solve5.reduced_stiffnes_2d

GKK(nz)

is contraint-reduced global stiffness matrix in 1D (vector) form

march1.global_mass
solve1.global_mass

GM(ny,ny)

is global mass matrix (2nd order time derivatives)


GMC(NCYMX,NCYMX)

for BEM problems (solve4.f)

march1.gr1_rhs_vector
nonlin.gr1_rhs_vector
solve1.gr1_rhs_vector

GR1(ny)

is RHS vector for system of global equations (with constraints)

nonlin.gr2_rhs_vector

GR2(ny)

is RHS vector for system of global equations (with constraints) (BFGS method)


GRC(NCYMX)

is the vector resulting from BEM domain integrals (fe90/solve4.f)

fit.grr_rhs_vector
march1.grr_rhs_vector
nonlin.grr_rhs_vector
solve1.grr_rhs_vector
solve5.grr_rhs_vector

GRR(no)

is RHS vector for system of constraint-reduced global equations


GRRC(NCYMX)

for BEM problems (solve4.f)


GXL,GXU

are undeformed metric tensor components wrt Xi-coords (zere50.f).


GZ

is the determinant of GZL


GZL,GZU

are the covariant & contravariant components

H

fit.smoothing_matrix

HK(ny,ny)

smoothing matrix (KTYP12 > 0) for fitting with smoothing constraints

I

elem.index_basis_type

IBT(1,ni,nb)

is index for basis type nb in Xi direction ni:
=1 for Lagrange
=2 for Hermite
=3 for Simplex
=4 for Serendipity
=5 for B-spline
=6 for Infinite
=7 for Transition
=8 for Singular
=9 for Fourier


IBT(2,ni,nb)

=1 for Lagrange linear
=2 for Lagrange quadratic
=3 for Lagrange cubic
=1 for Hermite cubic
=1 for B-Spline linear
=2 for B-Spline quadratic
=3 for B-Spline cubic
=N for Fourier series (number of terms in series incl. constant)

elem.index_deriv_order

IDO(nk,0,nb)

is index for derivative order: NU partial derivative number to derivative nk (seeNU table) of basis nb


IDO(nk,ni,nb)

is an index for derivative order: 1=zeroth order, 2=first order
Thus: IDO(nk,1..,nb)=1,2,1 implies that derivative number nk is a first derivative wrt Xi(2).

colloc.extern_face

IECG(NK,NN,NH,NE)

"external face" tables to figure out which element parameters are handled by the Galerkin method. [collocation]
=1 if the element DOF corresponding to nk,nn,nh,ne is handled with collocation,
=2 if DOF is handled by the boundary galerkin method,
=3 if DOF handled by interior galerkin,
=0 otherwise

colloc.global_equa

IGCG(NYMX)

vector indicates which global equations are from collocation (=1) and which are from Galerkin conditions (>1). [collocation]

elem.index_nodal_pos

INP(nn,ni,nb)

gives the index for element node nn in each Xi direction.
Thus: INP(nn,1..,nb) = 1,2,2 indicates that node nn is the first node in the Xi(1) dir and second in the Xi(2) & Xi(3) directions.


IOTYPE

1 = prompt, 2 = read, 3 = write, 4 = read & list


ISAXES(iw)

is segment number of axes


ISBASE(nb)

is segment number of basis function type nb


ISCONO(nh,ne)

is segment number of contour numbers in element ne


ISCONT(nh,ne,nocont)

is segment number of contour nocont of variable nh in element ne


ISDANO(iw,ne)

is segment number of data point numbers


ISDAPR(iw,ne)

is segment number of data point projections


ISDATA(iw,nodata)

is segment number of data points at set nodata


ISDATR(iw,ne)

is segment number of data point trace in element ne


ISEG(nosg)

is 0,1,2 if segment not yet created / created but not visible / created and visible


ISELNO(iw,ne)

is segment number of element numbers


ISFACE(iw,nf)

is segment number of face nf


ISFANO(iw,nf)

is segment number of face numbers


ISFIBR(iw,ne,nofibr)

is segment number of fibres in element ne at set nofibr


ISGAUS(iw)

is segment number of Gauss points


ISGRID(iw)

is segment number of grid


ISHIST(0)

is segment number of time history axes & labels


ISHIST(np)

is segment number of time history at node np


ISINCR(iw)

is segment number of increments


ISISOC(iw,noisoc)

is segment number of isochrones


ISL2BE(nl)

is segment number of Bezier tangent line 1 on nl


ISL3BE(nl)

is segment number of Bezier tangent line 2 on nl


ISLINE(iw,noline)

is segment number of lines at set noline


ISLINO(iw)

is segment number of line numbers


ISMAP(nomap)

is segment number of map


ISMATL(iw,ne)

is segment number of material in element ne


ISN2BE(nl)

is segment number of Bezier control pt 1 on nl


ISN3BE(nl)

is segment number of Bezier control pt 2 on nl


ISNONO(iw,np)

is segment number of node numbers


ISREAC(iw)

is segment number of reactions


ISSECT(nosect)

is segment number of section


ISSTRE(ne,nostre)

is segment number of principal stresses at set nostre


ISSTRM(ne,nostrm)

is segment number of streamline


ISSURF(ne)

is segment number of surface grid in element ne


ISVELO(ne,novelo)

is segment number of velocity field

data.i_threshold

ITHRES(NGMX,NEMX)

used in threshold modelling activation pattern computations (fe30). Value at Gauss point ng of element ne is...

0

if not active

1

if active

2

if surrounded by active points


ITYP1

is 3,4,5 or 9 for use of FE30,FE40,FE50 or FE90


ITYP2

is equation type [formerly KTYP1 ? FJV]


ITYP3

is equation type qualifier [formerly KTYP11 ? FJV]


IWKDEF(0)

is number of open windows


IWKDEF(noiw)

noiw=1,IWKDEF(0) is list of open windows


IWKG(iw)

is 0 for nongraphics window (eg menu)


IWKS(iw)

= 0 workstation (window) iw is not defined (ie, not open)
= 1 workstation iw is defined but not active
= 2 workstation iw is defined and active


IWKT(iw)

= 1 for GKS workstation
= 2 for PHIGS workstation
= 3 for Frame grabber display screen


IWRIT1

controls output printing frequency.


IWRIT2

is 1,2 for equilibrium solution only / intermediate solutions also


IWRIT3

is 1,2 for solution vectors only / residual vectors also


IWRIT4

is 0..4 for output from linear / nonlinear solver:
= 0 no output
= 1 timing only
= 2 ... and pivots
= 3 ... and element matrices
= 4 ... and global matrices


IWRIT5

is 0..5 for output from iterative solver:
= 0 no output
= 1 total iterations & final residual norm
= 2 iteration number & current residual norm
= 3 ... and pseudo-residual vector Zk
= 4 ... and iteration vector Xk
= 5 ... and true residual vector (Rk = B - A*Xk)

J

data.elem_def
elem.elem_def
march1.elem_def
nodes.elem_def
problem.elem_def
render.elem_def

JTYP1

is 1,2 for elements defined by user / chosen from menu

data.ens_elem_map
elem.ens_elem_map
march1.ens_elem_map
nodes.ens_elem_map
problem.ens_elem_map
render.ens_elem_map

JTYP2

is 0,1 for ensemble-to-element map for nodal derivatives is standard / non-standard

data.coor_sys
elem.coor_sys
march1.coor_sys
nodes.coor_sys
problem.coor_sys
render.coor_sys

JTYP3

is 1..5 for coordinate system: rect. cartesian / cylindrical polar / spherical polar / prolate spheroidal / oblate spheroidal

data.sym_type
elem.symm_type
march1.symm_type
nodes.symm_typ
problem.symm_typ
render.symm_typ

JTYP4

is 1..3 for geometry unsymmetric / cyl. symm. / sph. symm.

Note:for cylindrical symmetry:

if JTYP3=1

radius is nj=2 coord.

if JTYP3=2,3

radius is nj=1 coord.

if JTYP3=4

radius=focus*sinhX1*sinX2

data.basis_rep
elem.basis_rep
march1.basis_rep
nodes.basis_rep
problem.basis_rep
render.basis_rep

JTYP5

is 1,2 for basis functions in Lagrange or Hermite / monomial format

data.coor_const
elem.coor_const
march1.coor_const
nodes.coor_const
problem.coor_const
render.coor_const

JTYP6

is 1,2 for global coordinate system constant / specified by elements

data.deformed_coor_sys
elem.deformed_coor_sys
march1.deformed_coor_sys
nodes.deformed_coor_sys
problem.deformed_coor_sys
render.deformed_coor_sys

JTYP7

is 1..5 for dependent variable coordinate system rect. cartesian / cylindrical polar / spherical polar / prolate spheroidal / oblate spheroidal

data.basis_output
elem.basis_output
march1.basis_output
nodes.basis_output
problem.basis_output
rener.basis_output

JTYP8

is 0,1 for no output / output of basis functions

data.fiber_field_type
elem.fiber_field_type
march1.fiber_field_type
nodes.fiber_field_type
problem.fiber_field_type
render.fiber_field_type

JTYP9

is 0,1,2 for fiber direction field not defined / defined / sheet direction field defined
JTYP9 = 1 for fiber directions lying in Xi1-Xi2 plane.
Note:The fiber coordinate Nu(1) lies in the Xi1-Xi2 plane. Nu(2) is orthogonal to the fiber coordinate & lies in the Xi1-Xi2 plane; the Nu(3) coordinate is orthogonal to this plane. The fiber angle eta is eta(1), the angle between the fiber coordinate & the Xi1 coordinate. The Nu coordinates are stress coordinates, and are orthonormal with metric a(i,j)=Kronecker delta.
Note:For JTYP9 = 2, the Nu(1) coordinate is aligned with the local fiber axis and lies in the Xi1-Xi2 plane. Nu(2) lies in the plane of the sheet normal to Nu(1), and Nu(3) is the mutually orthogonal axis normal to the sheet plane. The Nu(2) and Nu(3) axes lie in the (cross-fiber,radial)-plane and therefore the sheet coordinates are obtained by a simple rotation about the fiber axis by an angle of (90-beta) degrees, where beta is the "sheet angle" between the radial axis and Nu(2), and is measured with a positive sense from the radial axis up towards the cross-fiber axis.

data.isochoric_interp
elem.isochoric_interp
march1.isochoric_interp
nodes.isochoric_interp
problem.isochoric_interp
render.isochoric_interp

JTYP10

is 1,2,3 for type of 'radial' interpolation (JTYP3>1 only)

data.num_of_field_vars
elem.num_of_field_vars
march1.num_of_field_vars
nodes.num_of_field_vars
problem.num_of_field_vars
render.num_of_field_vars

JTYP11

is the number of additional geometric / field variables

data.fiber_ref_axis
elem.fiber_ref_axis
march1.fiber_ref_axis
nodes.fiber_ref_axis
problem.fiber_ref_axis
render.fiber_ref_axis

JTYP12

is 1,2 for fibres defined wrt Xi1 or Xi2 coordinates

data.fiber_ang_units
elem.fiber_ang_units
march1.fiber_ang_units
nodes.fiber_ang_units
problem.fiber_ang_units
render.fiber_ang_units

JTYP13

is 1,2 for fibers entered in degrees / radians

data.special_mesh_type
elem.special_mesh_type
march1.special_mesh_type
nodes.special_mesh_type
problem.special_mesh_type
render.special_mesh_type

JTYP14

is mesh type for specialized meshes as follows: 1 & 2) Fractal tree with branch parameters:
Ratio_Angle,Ratio_Length,Ratio_Diameter
B_angle_y(no_gen) : mean angle branch makes with y-axis
B_angle_xy(no_gen): mean angle branch makes with x,y-plane
B_angle_SD(no_gen): std dev of angles
Note:angles are generated with Normal distribution using mean & std dev obtained by dividing previous gen values by Ratio_Angle
B_length(no_gen)
B_diameter(no_gen)
B_volume(no_gen)
B_flow(no_gen) ?? <-- not needed
NW(ne)=generation number
1) Regular fractal tree
2) Stochastic fractal tree

data.extern_stim
elem.extern_stim
march1.extern_stim
nodes.extern_stim
problem.extern_stim
render.extern_stim

JTYP15

0,1 : external stimulus not defined/defined (collocation only)

elem.num_of_mFHN_params
march1.num_of_mFHN_params
nodes.num_of_mFHN_params
problem.num_of_mFHN_params
render.num_of_mFHN_params

JTYP16

number of material parameters for (m)FHN equations (collocation only)

data.scalar_field
elem.scalar_field
march1.scalar_field
nodes.scalar_field
problem.scalar_field
render.scalar_field

JTYP17

0,1 : scalar field (e.g. consistent strains) not defined/defined (collocation only)

K

data.prob_type
elem.prob_type
march1.prob_type
nodes.prob_type
problem.prob_type
render.prob_type

KTYP1

is 1..15 for problem type [or ITYP2 ? -- FJV]

for KTYP3 = 2:

KTYP1 = 1 is Linear elasticity
KTYP1 = 2 is Finite elasticity
KTYP1 = 3 is Advection-diffusion
KTYP1 = 4 is Wave equation
KTYP1 = 5 is Navier-Stokes equations
KTYP1 = 6 is Bio-heat transfer
KTYP1 = 7 is Maxwell equations
KTYP1 = 8 is Cardiac activation
KTYP1 = 9 is Oxygen transport

data.FE_method
elem.FE_method
march1.FE_method
nodes.FE_method
problem.FE_method
render.FE_method

KTYP2

is 1..5 for Galerkin finite elements / direct boundary elements / indirect boundary elements / orthogonal collocation / finite element collocation (formerly arrayITYP4)

data.time_domain_type
elem.time_domain_type
march1.time_domain_type
nodes.time_domain_type
problem.time_domain_type
render.time_domain_type

KTYP3

is 1..5 for static / time integration / modal analysis / Fourier analysis / buckling analysis (formerly arrayITYP5)

data.prob_linear
elem.prob_linear
march1.prob_linear
nodes.prob_linear
problem.prob_linear
render.prob_linear

KTYP4

is 1..2 for linear / nonlinear problem (formerly arrayITYP6)

data.init_state
elem.init_state
march1.init_state
nodes.init_state
problem.init_state
render.init_state

KTYP5

is 1..3 for initial solution zero / read in / restarted

data.gauss_pnt_fit
elem.gauss_pnt_fit
march1.gauss_pnt_fit
nodes.gauss_pnt_fit
problem.gauss_pnt_fit
render.gauss_pnt_fit

KTYP6

is 1 for Gauss point fitting; is number of boundary integral equation domains

data.equa_params
elem.equa_params
march1.equa_params
nodes.equa_params
problem.equa_params
render.equa_params

KTYP7

is 1..3 for equation parameters constant wrt time / user defined in subroutine USER / read from file IPC at each time step

data.fitting_type
elem.fitting_type
march1.fitting_type
nodes.fitting_type
problem.fitting_type
render.fitting_type

KTYP8

is 1..5 for geometry / fibre or field / motion / optical flow / Fourier motion fitting
=1 with ITYP6=2 is nonlinear geometric fitting problem
=2 with JTYP11 > 0 is field fitting problem (ITYP6=1) (& Gauss)
=4 is image fitting problem
=5 is motion fitting problem with Fourier basis
=6 is stripe fitting problem !new AAY 16 Nov 92

data.nonli_solver
elem.nonli_solver
march1.nonli_solver
nodes.nonli_solver
problem.nonli_solver
render.nonli_solver

KTYP9

is 1..4 for full Newton / modified Newton / BFGS inverse / element-by-element method

data.line_search
elem.line_search
march1.line_search
nodes.line_search
problem.line_search
render.line_search

KTYP10

is 1..2 for solution with no search / linear search

data.opt_prob
elem.opt_prob
march1.opt_prob
nodes.opt_prob
problem.opt_prob
render.opt_prob

KTYP11

is 1..5 for additional options on problem type [or ITYP3 ? -- FJV]

for KTYP1 = 8:

KTYP11 = 1 is Threshold activation
KTYP11 = 2 is Cable equation (1D)
KTYP11 = 3 is FitzHugh-Nagumo equations
KTYP11 = 4 is Noble-DiFrancesco equations
KTYP11 = 5 is unused
KTYP11 = 6 is unused
KTYP11 = 7 is unused
KTYP11 = 8 is Modified FitzHugh-Nagumo equations

data.constraint_type
elem.constraint_type
march1.constraint_type
nodes.constraint_type
problem.constraint_type
render.constraint_type

KTYP12

is 0..3 for fitting without / with constraints:
=0 for no smoothing constraints
=1 for smoothing constraints on displacement ("Sobolev smoothing") -- weights defined in WU(i,ne)
=2 for smoothing constraints on deformation (strain energy smoothing)
=3 for smoothing constraints on deformation wrt fibre coords

data.ode_sovler_type
elem.ode_sovler_type
march1.ode_sovler_type
nodes.ode_sovler_type
problem.ode_sovler_type
render.ode_sovler_type

KTYP13

is 1 if pressure read from file (PRESS.VSAERO)
PROPAGATION PROBLEMS (finite element/collocation): indicates type of ODE integrator:
  =1 for non-stiff equations (RKSUITE)
  =2 for stiff equations (LSODES)

data.inc_params
elem.inc_params
march1.inc_params
nodes.nc_params
problem.nc_params
render.nc_params

KTYP14

is > 0 if material parameter is incremented


KTYP15

is equation type parameter


KTYP16

is 1..2 for lowest / highest eigenvalues required


KTYP17

is number of eigenvalue pairs required


KTYP18

is number of subspace iteration vectors


KTYP19

is number of starting vectors


KTYP22

is 1..3 for time integration algorithm linear / quadratic / cubic


KTYP23

is fixed time step / automatic stepping / read from file


KTYP25

is type of driving function in Fourier analysis 1..3 for impulse / step / sine wave


KTYP26

is 1..2 for optimization of material params / geometric params


KTYP27

is 1..9 for type of minimization objective function
Note:KTYP26 = 3 KTYP27 = 1 is used for stripe intersection calculation


KTYP28

is number of sets of measurements in fit


KTYP31

is 1..2 for Cardaic activation model implemented forwards / backwards


KTYP43

is 0..3 for thermal strains not included / included as fixed initial strain / constrained by displacement b.c.s / fully coupled


KTYP45

is 1..4 for type of beam cross-section

data.elast_prob

KTYP51

is 1..6 for plain stress / plain strain / 3D / membrane/ thin shell / thick shell

data.compressibility

KTYP52

is 1..3 for compressible / incompressible / incomp with fluid perfusate

data.anisotropy

KTYP53

is 1..3 for isotropic / aeleotropic / aeleotropic + active fibres

data.hyperelasticity

KTYP54

is 1..3 for hyperelastic / Cauchy-elasticity / creep

data.kinematic_params

KTYP55

is 1..3 for princ. strain invariants / extension ratios / fibre strains

data.energy_function_form

KTYP56

is 1..3 for polynomial / special function / exponential strain energy function

data.press_bc_type

KTYP57

is no pressure bc / press.incr / press.read / volume incr / vol.computed
Note:If KTYP57>1 a pressure b.c. is applied to the Xi(3)=0,1 faces of elements with NW(ne)=2,3, respectively, or both if NW(ne)=4.

data.isochoric

KTYP58

is 1,2 for conventional / isochoric element

data.active_stress_terms

KTYP59

is number of terms in active fibre stress relation

data.SL_parameter

KTYP5A

is material parameter number of stress-free SL distribution

data.time_delay_param

KTYP5B

is material parameter number of time-delay variable

data.Vcf_param

KTYP5C

is material parameter number of capillary volume fraction


KTYP5D

is unused


KTYP71

is 1 if pressure loads read from file PRESS.VSAERO (ID=14) after flow solution by VSAERO


KTYP90

is saturated-unsaturated / heart-body/ ? coupling

L

data.associated_data_point

LD(nd)

line or face l number associated with data point nd


LDR(nd)

is unused and should be removed


LFR(no)

is unused and should be removed

prob.elem_location

LGE(nv)

is location of element variable nv in global system. LGE(nv) is negative if NE is the last element to reference that variable.


LGEE(nv,ne)

is location of element variable nv in global system for element ne (replaces LGE when using iterative solver)

data.fitting_elem

LN(0)

number of elements in fitting


LN(l)

element number for l=1..LN(0)

prob.boundary_region_var

LRE(nv)

is for the shell/fluid interface case. It is the boundary region variable LRE of the local finite element variable NVE. (fe02/melge.f). This may be unused.


LUMP

is .TRUE. if mass lumping is used

M

prob.last_influenced_elem

ME(np)

is the last element to be influenced by global node NP. WHERE IS THIS SET???


MOTION_TYPE

is 1,2 for Fourier coeffs / Spreadsheet column

N


NA

auxiliary element variable (na=1,NAT(nb))

dims.num_aux_var
elem.num_aux_var
prob.num_aux_var

NAMX

is maximum number of auxiliary parameters

elem.poly_degree

NAN(ni,na,nb)

is polynomial degree in Xi(ni) direction for basis nb for auxiliary variable na


NAT(nb)

is number of auxiliary or spline basis functions


NB

basis function type (nb=1,NBT)


NBC(nb)

is basis function type choice

elem.basis_num

NBH(nh,ne)

is basis number for dependent variable nh in element ne


NBI(nb)

is 1..5 for unit scale factors / elem scale factors read in / global scale factors read in / arc-length / angle-change

elem.basis_type_num

NBJ(nj,ne)

is basis function type number for geometric variable nj in element ne


NBL(0,0,nb)

is number of boundary line / face segments for blending functions
NOT USED ANY MORE!


NBL(0,ls,nb)

is nb number of boundary line / face segment ls


NBL(1..4,ls,nb)

are element node numbers of line / face ls (used for blending function interpolants)

colloc.num_basis
dims.num_basis
elem.num_basis
nodes.num_basis

NBMX

maximum number of basis functions


NBT

number of basis function types


NB_MOTION

Fourier basis number

data.data_in_elem
dims.data_in_elem
fit.data_in_elem

NCMX

maximum number of equations for a given dependent var (the old meaning of NCMX was the max number of data points in an element)

elem.coor_type
nodes.coor_type

NCO(ne)

is the coordinate type used in element ne
NOTE: removed inContinuity 6.0

dims.BE_array_size

NCYMX

maximum size of complex arrays (for BEM problems)


ND

data point (nd=1,NDT)


NDAL(0:nde)

is unused and should be removed

data.global_data_pts

NDDL(ne,nde)

global data pt no. of local data point nde


NDEMX

number of data points in one element

data.num_of_data_pts

NDLT(ne)

number of data points within element ne

data.data_points
dims.data_points

NDMX

number of data points total

data.data_pts_index

NDP(NDMX)

index of data points that project onto an element


NDT

is total number of data points


NE

an element (ne=1,NET)

nodes.shared_corner

NECOR(NPMX,0:4)

stores the number of elements which share the current corner or edge node NP. NECOR(NP,0) is the total number of elements sharing the corner node (BEM problems, fe90/solve4.f)


NEELEM(0)

total number of elements.


NEELEM(noelem)

noelem=1..NEELEM(0) are the element numbers.

colloc.elements
data.elements
dims.elements_nemx
elem.elements
march1.elements
nodes.elements
nonlin.elements
prob.elements

NEMX

maximum number of elements


NEIBS(NI,IN,NE)

Element neighbors of NE on INth face in NIth direction [collocation]

colloc.local_node

NEP(1,IN,np)

the INth local node number coincident with global node np [collocation]


NEP(2,IN,np)

the element the INth local node belongs to. [collocation]


NET

is the highest element number


NF

global face (nf=1,NFT)


NFE(nb)

number of faces for element basis type nb.

elem.global_face_num

NFF(nf,ne)

are the global face numbers of side nf of element ne

dims.global_face_seg
elem.global_face_seg
nodes.global_face_seg

NFMX

maximum number of global face segments


NFT

total number global element faces


NG

Gaussian quadrature point (ng=1,NGT(nb))

elem.num_gauss_pts

NGAP(ni,nb)

is number of Gauss points in Xi direction ni for basis nb

colloc.gaus_per_elem
data.gaus_per_elem
dims.gaus_per_elem
elem.gaus_per_elem
prob.gaus_per_elem

NGMX

maximum number Gauss points per element


NGT(nb)

number of Gauss points per element


NH

dependent variable (nh=1,NHP(np))

elem.num_dep_var_elem

NHE(ne)

number of dependent variables defined in element ne

colloc.depen_var
data.depen_var
dims.depen_var
nodes.depen_var
prob.depen_var

NHMX

maximum number of dependent variablesNote:NHMX must be equal to NJMX for problems where the dependent variable array carries deformed coordinates.


NHO

is Gauss variable to be fitted in Gauss point fitting

nodes.num_dep_var_node

NHP(np)

is number of dependent variables defined at node np


NHT(nje,nve,ie)

is number of global variables required for element type ie when number of local variables is nve & number of dimensions nje.

where:

mh,nh are global element variables, = 1..nhe


mv,nv are local element variables, = 1..nve


NHV(nv,ie)

is global variable number, 1 .. NHE(ne), corresponding to dependent variable nv in element type ie.


NI

Xi-coordinate (ni=1,NIT(nb))

colloc.local_xi
data.local_xi
dims.local_xi
elem.local_xi
fit.local_xi

NIMX

maximum number of local Xi-coordinates


NIT(nb)

is number of local Xi-coordinates for basis nb


NJ

Xj-coordinate (nj=1,NJT)

elem.num_xj_coor_in_elem

NJE(ne)

number of Xj-coords defined in element ne (excluding JTYP9)


NJG

is geometric variable number in linear field fitting.

colloc.global_ref_coor
data.global_ref_coor
dims.global_ref_coor
elem.global_ref_coor
fit.global_ref_coor
nodes.global_ref_coor
prob.global_ref_coor

NJMX

maximum number of global reference coordinates
Note:NHMX must be equal to NJMX for problems where the dependent variable array carries deformed coordinates.


NJO

is field variable number in linear field fitting (eg =NJT+1).

nodes.num_xj_coor

NJP(np)

is number of Xj-coordinates (ie, geometric variables) defined at node np
Node:removed inContinuity 6.0(can use a temp variable when required)


NK

derivative number (nk=1,NKT(nb))

nodes.nodal_derv_dep_var

NKH(nh,np)

is number of nodal derivatives for dependent variable nh at node np (also called nodal_derv_dep_var in some places)

nodes.nodal_derv_geom_var

NKJ(nj,np)

is number of nodal derivatives for geometric variable nj at node np

colloc.deriv_per_var
dims.deriv_per_var
elem.deriv_per_var
nodes.deriv_per_var
prob.deriv_per_var

NKMX

maximum number of derivatives per variable

elem.nodal_spline_term

NKT(nb)

is number of nodal derivatives or spline polynomial terms for basis nb


NL

global line (nl=1,NLT)


NLE(nb)

is number of element line segments for basis nb

elem.line_num_of_arcface

NLF(naf,nf)

are global line numbers of local arc naf of face nf

elem.line_num_of_arcelem

NLL(nae,ne)

are global line numbers of local arc nae of element ne

dims.global_line_seg
nodes.global_line_seg

NLMX

maximum number of global line segments


NLT

total global line segments


NLV(ie,njt)

is number of local variables possible in element type ie.


NM

material parameter (nm=1,NMT)

dims.materials
prob.materials

NMMX

maximum number of material parameters


NMOPTI(noopti)

noopti=NMOPTI(0)+1,NTOPTI is list of displacements in fit


NN

element node (nn=1,NNT(nb))

elem.num_of_elem_nodes_in_face

NNF(0,nf,nb)

is number of element nodes in face NF of element with basis NB


NNF(1 ,nf,nb)

is Xi-direction normal to face


NNF(11.14,nf,nb)

are element derivative numbers in face


NNF(2..10,nf,nb)

are the element node numbers in face

elem.elemnode_face_num

NNL(1..4,nae,nb)

are the element / face node numbers along local arc nae

colloc.nodes_pre_elem
dims.nodes_per_elem
elem.nodes_per_elem
nodes.nodes_per_elem

NNMX

maximum number of nodes per element

colloc.coincident_nodes

NNP(np)

the number of element nodes coincident with global node np [collocation]

elem.elem_nodes

NNT(nb)

is number of element nodes for basis nb


NO

a d.o.f in the set of global equations (no=1,NOT)

data.dof
dims.dof
fit.dof
march1.dof
nonlin.dof
solve1.dof
solve5.dof

NOMX

maximum number of optimization DOFs.
Note:NZMX must be > NOMX2 and NYMX2 for some problems and large enough to take the fill-in during sparse matrix solution in others


NORMX

maximum number of size of fractal tree order arrays


NOT

total number of optimization DOFs

data.mesh_opt_dof
dims.mesh_opt_dof

NOYMX

maximum number of optimization DOFs attached to one mesh DOF

data.num_of_optimization_dof

NOYT(ny)

number of optimization DOFs coupled to mesh variable ny

data.optimization_dof_num

NOYY(noy,ny)

optimization DOF number of variable noy for mesh variable ny


NP

global node (np=1,NPT)


NPCOR(0:NPMX,2)

NPCOR(i,1)=0 if node i is not on a corner =column number of GM matrix for second dH/dn contribution if node i is at some corner or edge (3D) (the column number corresponds to the contribution from the element at the corner with the higher element number if it is a 2d corner or 3d edge, otherwise it corresponds to the middle corner element number).
NPCOR(i,2)=0 unless node i is at a 3D corner in which case it gives the column number of the GK matrix for the third dH/dn contribution (corresponding to the highest element number of the elements at the corner).
NPCOR(0,1)=number of corner nodes (2d) and edge nodes (3d).
NPCOR(0,2)=number of corner nodes (3d).
(BEM problems, fe90/solve4.f)

dims.domain_nodes

NPDMX

maximum number of domain nodes for BEM problems

nodes.global_node_num

NPE(nn,nb,ne)

is global node number of local node nn of element ne for basis nb


NPE(nn,nb,nf)

are global node numbers of face nf, where nb is basis type number for first geometric variable.

nodes.nodal_basis_type

NPF(1,nf)

is the 1st Xi-direction of face segment NF


NPF(2,nf)

is the basis fn type for 1st Xi-dir (1,2,3 or 4)


NPF(3,nf)

is the 2nd Xi-direction


NPF(4,nf)

is the basis function type for 2nd Xi-dir


NPF(5..8,nf)

are the basis function type (nb) numbers for nj=1..4


NPF(9,nf)

is number of elements adjoining face (1,2 for external / internal)

nodes.connect

NPL(1,nl)

is the Xi-direction of line segment NL


NPL(2..5,nl)

are the basis function types for nj=1..4 :
=1 linear Lagrange
=2 quadratic Lagrange
=3 cubic Lagrange
=4 cubic Hermite


NPL(6..9,nl)

are the global nodes along line NL in the direction of Xi (8 & 9 are NK numbers of 1st derivatives wrt Xi for any geometric variable which is cubic Hermite)


NPL(10,nl)

is the number of elements adjoining line NL


NPL(11..18,nl)

are the element numbers of elements adjoining line NL

colloc.global_nodes
data.global_nodes
dims.global_nodes
nodes.global_nodes
prob.global_nodes

NPMX

maximum number of global nodes


NPNODE(0)

is the total number of nodes.


NPNODE(nonode)

nonode=1..NPNODE(0) are the node numbers.

data.global_node_fitting

NPO(0)

number of global nodes in data fitting


NPO(np)

global node number corresponding to n=1..NPO(0)


NPT

is the highest node number.


NP_MOTION

is node number for applying motion


NQ

a general DOF (nq=1,NQT)

nodes.global_DOF_num

NQE(ns,nb,ne)

is global DOF number of local DOF ns in element ne for basis nb.
Note:If JTYP2=1 NQE stores the global NK number of local DOFs ns for element ne.

dims.global_dof
nodes.global_dof

NQMX

maximum number of global DOFs


NQT

total number of global DOFs / geometric variables in XA

dims.contours_per_elem

NRMX

maximum number of contours per element & various segs


NS

ns=nk+(nn-1)*NKT(nb)

colloc.elem_dof
data.elem_dof
dims.elem_dof
elem.elem_dof
march1.elem_dof
nodes.elem_dof
prob.elem_dof
solve1.elem_dof
solve5.elem_dof

NSMX

maximum number of element DOFs per variable ( <= NNMX*NKMX )


NST(nb)

total number of element parameters for basis nb [ =NKT(nb)*NNT(nb)+NAT(nb) ]


NTACTV

is number of dynamic terms in the material response function

dims.eigenvalues

NTMX

maximum number of eigenvalues


NTSG

is current number of segments


NTYP1

is 1..16 for type of iterative solver:
= 1 : Conjugate Gradient (CG)
= 2 : Conjugate Gradient Squared (CGS)
= 3 : BiConjugate Gradient Stabilized (Bi-CGSTAB)
= 4 : Bi-CGSTAB with Restart (RBi-CGSTAB)
= 5 : Transpose-free Quasi-Minimal Residual (TFQMR)
= 6 : Generalized Minimal Residual with Restart (RGMRES)
= 7 : RGMRES with eigenvalue estimation (RGMRESEV)NOTE:requires diagonal or left preconditioning.
= 8 : Generalized Conjugate Residual with Restart (RGCR)
= 9 : Chebyshev acceleration (CHEBYSHEV).NOTE:CHEBYSHEV requires min & max eigenvalue estimates.
=10 : unused =11 : Bi-Conjugate Gradient (Bi-CG)
=12 : CG Normal Equations (CGNE)
=13 : CG Normal Residuals (CGNR)
=14 : Quasi-Minimal Residual (QMR)
=15 : unused
=16 : unused


NTYP2

is the number of Krylov subspace basis vectors (used for NTYP1 = 4, 6, 7, or 8)NOTE:max is 250 for RGMRES & RGMRESEV.


NTYP3

is the stopping criterion for the iterative solver


NTYP4

is the maximum number of iterative solver restarts or iterations


NTYP5

is the type of preconditioning for the iterative solver:
= 0 : none
= 1 : left diagonal scaling (Jacobi)
= 2 : right diagonal scaling (Jacobi)
= 3 : symmetric diagonal scaling (Jacobi)
= 4 : left preconditioning (with matrix)
= 5 : right preconditioning (with matrix)
= 6 : symmetric preconditioning (with matrix)
= 7 : unused
= 8 : unused

colloc.deriv_terms
data.deriv_terms
dims.deriv_terms
elem.deriv_terms

NUMX

maximum number of derivative terms up to 2nd order


NU

is index of Xi-coordinate derivative

NUT

derivative

1

u

2

u1

3

u11

4

u2

5

u22

6

u12

7

u3

8

u33

9

u23

10

u13

11

u123

e.g.: u1 = partial(u)/partial(Xi_1)


NUT(nb)

is number of Xi-coordinate derivatives up to 2nd order for basis nb


NVE(ie)

is number of local variables used in element type ie

colloc.total_elem_dof
fit.total_elem_dof
march1.total_elem_dof
nonlin.total_elem_dof
prob.total_elem_dof
solve1.total_elem_dof
solve5.total_elem_dof

NVMX

maximum number of element DOFs ( = NHMX*NSMX )

elem.type_num

NW(ne)

is the type number (= 1,2..12) for element ne.

dims.workstations

NWMX

maximum number of workstations for segment arrays (e.g. ISAXES, ISFIBR...)

elem.elem_adjacency

NXI(-ni:ni,0:ne)

Filled in fe02/nenxi.f. Element number adjacent to element ne:
NXI(-3,ne) is element in -Xi(3) direction
NXI(-2,ne) is element in -Xi(2) direction
NXI(-1,ne) is element in -Xi(1) direction
NXI( 1,ne) is element in +Xi(1) direction
NXI( 2,ne) is element in +Xi(2) direction
NXI( 3,ne) is element in +Xi(3) direction
NXI( ., 0) is 0

dims.image_cell

NXMX

maximum number of image cell array dimension (ASSERTed in deimag.f)


NY

a finite element mesh DOF (ny=1,NYT)

colloc.mesh_dof
data.mesh_dof
dims.mesh_dof
march1.mesh_dof
nodes.mesh_dof
nonlin.mesh_dof
prob.mesh_dof
solve1.mesh_dof

NYMX

maximum number of mesh DOFs.
Note:NZMX must be > NOMX2 and NYMX2 for some problems and large enough to take the fill-in during sparse matrix solution in others

prob.ny_ne_mapping

NYNE(na,nh,ne)

is mapping from na,nh,ne to ny (set in ipini5.f)

prob.ny_np_mapping

NYNP(nk,nh,np)

is mapping from nk,nh,np to ny (set in ipini5.f)


NYT

total finite element mesh DOFs (initialized in ipini5.f)


NZ

a stiffness matrix component in a 1D vector (nz=1,NZT)

colloc.stiffness_vector
dims.stiffness_vector
fit.stiffness_vector
march1.stiffness_vector
nonlin.stiffness_vector

NZMX

maximum number of coefficients in 1D global stiffness vector GKK.
Note:NZMX must be > NOMX2 and NYMX2 for some problems and large enough to take the fill-in during sparse matrix solution in others.


NZT

total number of nonzero coefficients in 1D global stiffness vector GKK.

O

P

prob.pressure_incr

PE(1..2,ne)

the pressure increment applied to the Xi_3=0 face (1) or Xi_3=1 face (2) of element ne

prob.pressure_load

PF(1..2,ne)

the pressure load applied to the Xi_3=0 face (1) or Xi_3=1 face(2) of element ne


PG(ns,nu,ng,nb)

are basis function values for element DOF ns at Gauss point ng, basis number nb. Set in fe02/gauss1.f when bases are defined.
nu=     1 is zeroth order derivative,
nu= 2..10 are derivatives up to 2nd order
nu=    11 is third order derivative


PGF(ns,nu,ng,nf,nb)

basis function values (as in PG) on face nf. [collocation]


PHI

are the Euler angles wrt COORDS of the principal extensions


PMAX(noopti)

maximum parameter values allowed


PMIN(noopti)

minimum parameters value allowed


PROMPT

is .TRUE. if computations await prompt after IWRIT1 steps.


PST

are principal strains / stresses

Q

R


R

is the orthogonal rotation tensor


RAD(NGMX)

for BEM problems (solve4.f)


RADS(NGMX)

for BEM problems (solve4.f)

prob.elem_resid_1

RE1(ns,nh)

element residual array

prob.elem_resid_2

RE2(ns,nh)

element residual array

elem.jacobian

RG(ng)

is the Jacobian for a length, area, or volume integral (computed as RG and RGX in fe02/xgmg.f)


RGS(NGMX)

for BEM problems (solve4.f)


RHO(NNMX,NJMX,2*NNMX)

are the (rho,phi) coordinates of a 2D boundary element which has been subdivided with a local polar coordinate system. NNMIN is set up in solve4.f and is the local node at which the origin of the local polar coordinate scheme is located. BEM
NOT USED ANY MORE!


RI1,RI2,RI3

are principal invariants of AZL


RM

is the modal matrix whose columns are the eigenvectors assoc with PST

prob.global_resid_vector

RP1(ny)

global residual vector


RWG(ng)

is the square root of the determinant of the undeformed metric tensor multiplied by the quadratrure weight at gauss point ng (this variable isRGXin fe50/zeex50.f)

S

data.dof_opti_scale

SCALE(no)

scale associated with optimization DOF no for geometric fitting (see fe03/ipfit.f and fitgeo.f). Currently unused and should be removed.

elem.scale_factor

SE(ns,nb,ne)

are scaling factors for basis nb of element ne. Arc-length scale factors are set in fe02/lincal.f and fe02/dlse.f; face scaling factors are set in fe02/facseg.f


SNLPA

is static nonlinearity parameter "a"

data.squared_dist

SQ(nd)

square of distance from mesh to data point nd


SS

total sum of squared distances SQ(nd)

T


TC

are physical components of Cauchy stress


TG

are tensor components of 2nd Piola-Kirchhoff stresses

data.threshold

THRES(3,ng,ne)

used in threshold modelling activation pattern computations (fe30). Value at Gauss point ng of element ne is...

1

time since Gauss point became active

2

membrane potential at Gauss point

3

recovery variable at Gauss point


TN

are physical components of Nominal stresses


TV_SLO

is slope of the force / velocity relation in stretching (before yield)

U


U

is the right stretch tensor

V

elem.b_spline

VE(ns,nk,ne)

are B-spline polynomial coefficients for each DOF ns of element ne

W

data.weighting_factor

WD(nj,nd)

weighting factor for data point nd (this is gamma in Hashima et.al. 1993)

fit.line_waiting_factor

WDL(nj,nde)

weighting factor for line data point nde

elem.gaus_weight

WG(ng,nb)

are Gauss point weights at Gauss point ng for basis number nb


WGF(ng,nf,nb)

gaussian integration weights (as in WG) on face nf. [collocation]


WS(5*NYMX)

is unused and should be removed

data.Sobolev_smoothing_weights

WU(nu,ne)

Sobolev smoothing weights for DOF nu (these are alpha & beta in Hashima et.al. 1993)

X


XB(1..2,NJMX,NLMX)

is first/second Bezier slope control point for line NL (fe21/crherm.f)

elem.xj_geom_pos_XE

XE(ns,nj)

is Xj geometric position or derivative for element DOF ns, coordinate nj


XD(NKMX,NPDMX)

for BEM problems (solve4.f)


XDC(NKMX,NPDMX)

for BEM problems (solve4.f)

elem.xj_geom_pos_XG

XG

are undeformed theta coords and derivatives wrt Xi


XG(nj,nu)

is Xj geometric position (nu=1) or derivative (nu>1) for coord nj within an element


XG1(NJMX,NUMX,NGMX)

for BEM problems (solve4.f)


XGS(NJMX,NUMX,NGMX)

for BEM problems (solve4.f)

data.xi_coor_data

XID(ni,nd)

Xi-coordinate of data point nd

fit.xi_of_line

XIDL(ni,nde)

Xi-coordinate of line data point nde


XIDR(ni,nd)

is unused and should be removed

elem.xi_gaus

XIG(ni,ng,nb)

is Xi-coordinate (ni) at Gauss point ng for basis nb


XIGF(ni,ng,nf,nb)

gauss point locations (as in XIG) on face nf. [collocation]


XMGF(nfg,ng,njmat)

Material parameter njmat evaluated at gauss point ng of face nfg [collocation]


XN(NJMX,NGMX)

for BEM problems (solve4.f)


XNS(NJMX,NGMX)

for BEM problems (solve4.f)

data.opti_var

XO(no)

optimization variable


XOC(NCYMX)

for BEM problems (solve4.f)

nodes.global_param_XP

XP(nk,nj,np)

is Xj geometric position (nk=1) or derivative (nk>1) for coordinate nj at global node np


XPD(NKMX,NJMX,NPDMX)

for BEM problems (solve4.f)

Y


YD(NKMX)

calculated BEM solution at the point XPFP (fe90/domsol.f)


YDC(NKMX)

calculated BEM solution at the point XPFP (fe90/domsol.f)


YIELDR

is ratio of yield tension to isometric tension

data.data_at_guass

YG(ng,nj,ne)

is data defined at Gauss point ng (used in fe03/yger.f and fe03/fitgau.f)
- or -
in fe40/opst40.f:

YG(ng,1,ne)

1st principal stresses

YG(ng,2,ne)

2nd principal stresses

YG(ng,3,ne)

principal angles

YG(ng,4,ne)

strain energy

YG(ng,5,ne)

material density

- or -
is updated with ZG, principal stresses, activation timing, material parameters, or data proximity in fe21/upgaus.f.

nodes.mesh_DOF

YP(ny,1..16)

is value of mesh DOF ny

for FE30 & linear (KTYP4=1) FE40 problems:

YP(ny,1)

essential b.c.s as defined by FIX(ny,1)

YP(ny,2)

natural b.c.s as defined by FIX(ny,2)

YP(ny,4)

initial solution

for FE50 & nonlinear (KTYP4=2) FE40 problems:

YP(ny,1)

incremental solution and is used in NONLIN

YP(ny,2)

current force b.c.s which are updated in NONLIN

YP(ny,3)

prescribed b.c. increments defined by FIX(NY,2..3) (initialized in ipini5.f)

YP(ny,4)

initial equilibrium solution (initialized in ipini5.f)

YP(ny,5)

stores residual vector (RP1) of previous iteration within NONLIN; contains RP1 on exit

YP(ny,6..16)

used for update information in BFGS & Broyden schemes, or for previous time solutions in Newton & modified Newton schemes (in NONLIN)

YP(ny,17)

applied force b.c.s for elastica problem (NW=13)

Z

data.rect_cart_coord_data

ZD(nj,nd)

rectangular cartesian coords of data point nd

data.rect_cart_coord

ZDL(nj,nde)

rectangular cartesian coords of line data point nde

data.deformed_param_array

ZE(ns,nh)

local deformed element parameter array [FJV 9 April 1996]


ZEC(ns,nh)

used for BEM problems. This is currently unused and should be removed.


ZFC(ns,nh)

used for BEM problems. This is currently unused and should be removed.

data.deformed_theta_coor

ZG(nh,nu)

are deformed theta coordinates and derivatives wrt undeformed coordinates (computed in fe02/zezg.f):

1

is value of dependent variable nh at current Gauss point

2,4,7

are 1st derivatives wrt Xi material coords (JP=1)

or:

are 1st derivatives wrt Nu or Theta coords (JP=2)


ZK

is unused and should be removed

nodes.deform_global_param_ZP

ZP(nk,nh,np)

global deformed element parameter array [FJV 9 April 1996]

Notes

  1. Material Nu-coordinates (reference for stresses): are orthogonal and (Nu1,Nu2) lie in the (Xi1-Xi2) plane such that Nu(1) is aligned with the `fibres' to which material aeolotropy is referred; The undeformed base vectors are defined such that the undeformed metric tensors wrt the Nu are delta(i,j).
  2. The fibre coord Nu(1) lies in the Xi1-Xi2 plane Nu(2) is orthog to the fibre coord & lies in the Xi1-Xi2 plane & the remaining Nu coord is orthog to this plane. The fibre angle eta is eta(1), the angle between the fibre coord & the Xi1 coord. The Nu coords are stress coords & are orthonormal with metric a(i,j)=Kronecker delta.

  3. KTYP26 = 3 KTYP27 = 1 is used for stripe intersection calculation.
  4. If KTYP57>1 a pressure b.c. is applied to the Xi(3)=0,1 faces of elements with NW(ne)=2,3, respectively, or both if NW(ne)=4.

  5. NHMX must be equal to NJMX for problems where the dependent variable array carries deformed coordinates.
  6. NZMX must be > NOMX2 and NYMX2 for some problems and large enough to take the fill-in during sparse matrix solution in others.

  7. If JTYP2=1 NQE stores the global NK number of local DOFs ns for element ne.

    For cylindrical symmetry (JTYP4=2):

    if JTYP3=1

    radius is nj=2 coord.

    if JTYP3=2,3

    radius is nj=1 coord.

    if JTYP3=4

    radius=focus*sinhX1*sinX2

  8. Dependent and Geometric Variables counterparts

    Dep Var Fortran Name

    Dep Var Python Name

    Geometric Var Fortran Name

    Geom Var Python Name

    nh

    nj

    zp

    xp

    ze

    xe

    zg

    xg

    nbh

    elem.basis_num

    nbj

    elem.basis_type_num

    nhe

    nje

    nhp

    njp

    nht

    njt

    nkh

    nkj



Obsolete

GB(no)

is RHS vector for system of reduced global equations

GS(nz)

is reduced global solution matrix in 1D form

INDEX_PLIN

is the index for a particular polyline.

INDEX_PLIN_TYPE(index_plin)

indicates the type (piecewise linear / Bezier).

ITYP4

is Galerkin FEM / Direct bem / Indirect bem / Orthogonal collocation (replaced by KTYP2)

ITYP5

is static / time integration / modal / Fourier / Laplace / buckling (replaced by KTYP3)

ITYP6

is linear / nonlinear problem (replaced by KTYP4)

LXI( 0,nl)
LXI( ni,nl)
LXI(-ni,nl)

replaced by NXI

NCNP(nh,np)

is the mapping between nc and nh, and np which is set up in IPEQUA or DEEQUA and possibly modified in DECORN.

NKE(nk,nb,ne)

is global derivative number of local derivative nk of element ne for basis nb

NO_BRANCH_PATTERN(no_gen,inumber_of_branches)

is number of elements which have inumber_of_branches at generation no_gen.

NT_PLIN

is the total number of polylines.

NT_PLIN_SECTIONS(index_plin)

is the no of sections to the polyline.

PLIN_DATA(nj,no_point,index_plin)

are the coordinates of each point of the polyline.

RP2(ny)

formerly passed to GENSOL, now unused.

YNP(nk,nh,nc,np,0)

is the global mapping between nk,nh,nc,np and ny or use in coupled problems (yet to be set up).



A

B

C

D

E

F

G

H

I

J

K

L

M

N

O

P

Q

R

S

T

U

V

W

X

Y

Z