variable | meaning |
REARRANGEMENT OF THE ORDER IN THE INPUT DECK | |
ifreeinp | next blank line in field inp |
ipoinp(1,i) | index of the first column in field inp containing information on a block of lines in the input deck corresponding to fundamental key i; a fundamental key is a key for which the order in the input file matters (the fundamental keys are listed in file keystart.f) |
ipoinp(2,i) | index of the last column in field inp containing information on a block of lines in the input deck corresopnding to fundamental key i; |
inp | a column i in field inp (i.e. inp(1..3,i)) corresponds to a uninterupted block of lines assigned to one and the same fundamental key in the input deck. inp(1,i) is its first line in the input deck, inp(2,i) its last line and inp(3,i) the next column in inp corresponding to the same fundamental key; it takes the value 0 if none other exists. |
MATERIAL DESCRIPTION | |
nmat | # materials |
matname(i) | name of material i |
nelcon(1,i) | # (hyper)elastic constants for material i (negative kode for nonlinear elastic constants) |
nelcon(2,i) | # temperature data points for the elastic constants of material i |
elcon(0,j,i) | temperature at (hyper)elastic temperature point j of material i |
elcon(k,j,i) | (hyper)elastic constant k at elastic temperature point j of material i |
nrhcon(i) | # temperature data points for the density of material i |
rhcon(0,j,i) | temperature at density temperature point j of material i |
rhcon(1,j,i) | density at the density temperature point j of material i |
nshcon(i) | # temperature data points for the specific heat of material i |
shcon(0,j,i) | temperature at specific heat temperature point j of material i |
shcon(1,j,i) | specific heat at the specific heat temperature point j of material i |
nalcon(1,i) | # of expansion constants for material i |
nalcon(2,i) | # of temperature data points for the expansion coefficients of material i |
alcon(0,j,i) | temperature at expansion temperature point j of material i |
alcon(k,j,i) | expansion coefficient k at expansion temperature point j of material i |
ncocon(1,i) | # of conductivity constants for material i |
ncocon(2,i) | # of temperature data points for the conductivity coefficients of material i |
cocon(0,j,i) | temperature at conductivity temperature point j of material i |
cocon(k,j,i) | conductivity coefficient k at conductivity temperature point j of material i |
orname(i) | name of orientation i |
orab(1..6,i) | coordinates of points a and b defining the new orientation |
norien | # orientations |
isotropic hardening | |
---|---|
nplicon(0,i) | # temperature data points for the isotropic hardening curve of material i |
nplicon(j,i) | # of stress - plastic strain data points at temperature j for material i |
plicon(0,j,i) | temperature data point j of material i |
plicon(2*k-1,j,i) | stress corresponding to stress-plastic strain data point k at temperature data point j of material i |
plicon(2*k,j,i) | plastic strain corresponding to stress-plastic strain data point k at temperature data point j of material i |
kinematic hardening | |
nplkcon(0,i) | # temperature data points for the kinematic hardening curve of material i |
nplkcon(j,i) | # of stress - plastic strain data points at temperature j for material i |
plkcon(0,j,i) | temperature data point j of material i |
plkcon(2*k-1,j,i) | stress corresponding to stress-plastic strain data point k at temperature data point j of material i |
plkcon(2*k,j,i) | plastic strain corresponding to stress-plastic strain data point k at temperature data point j of material i |
kode=-1 | Arrudy-Boyce |
-2 | Mooney-Rivlin |
-3 | Neo-Hooke |
-4 | Ogden (N=1) |
-5 | Ogden (N=2) |
-6 | Ogden (N=3) |
-7 | Polynomial (N=1) |
-8 | Polynomial (N=2) |
-9 | Polynomial (N=3) |
-10 | Reduced Polynomial (N=1) |
-11 | Reduced Polynomial (N=2) |
-12 | Reduced Polynomial (N=3) |
-13 | Van der Waals (not implemented yet) |
-14 | Yeoh |
-15 | Hyperfoam (N=1) |
-16 | Hyperfoam (N=2) |
-17 | Hyperfoam (N=3) |
-50 | deformation plasticity |
-51 | incremental plasticity (no viscosity) |
-52 | viscoplasticity |
![]() |
user material routine with -kode-100 user defined constants with keyword *USER MATERIAL |
PROCEDURE DESCRIPTION | |
iperturb | 0 : linear |
1 : second order theory | |
2 : nonlinear geometric | |
3 : nonlinear elastic material (and nonlinear geometric) | |
nmethod | 1 : static (linear or nonlinear) |
2 : frequency(linear) | |
3 : buckling (linear) | |
4 : dynamic (linear or nonlinear) | |
GEOMETRY DESCRIPTION | |
nk | highest node number |
co(i,j) | coordinate i of node j |
inotr(1,j) | transformation number applicable in node j |
inotr(2,j) | a SPC in a node j in which a transformation applies corresponds to a MPC. inotr(2,j) contains the number of a new node generated for the inhomogeneous part of the MPC |
TOPOLOGY DESCRIPTION | |
ne | highest element number |
mint_ | max # of integration points per element (max over all elements) |
kon(i) | field containing the connectivity lists of the elements in successive order |
For element i | |
ipkon(i) | (location in kon of the first node in the element connectivity list of element i)-1 |
lakon(i) | element label |
ielorien(i) | orientation number |
ielmat(i) | material number |
SHELL (2-D) AND BEAM (1-D) VARIABLES (INCLUDING PLANE STRAIN, PLANE STRESS AND AXISYMMETRIC ELEMENTS) | |
iponor(2,i) | two pointers for entry i of kon. The first pointer points to the location in xnor preceding the normals of entry i, the second points to the location in knor of the newly generated dependent nodes of entry i. |
xnor(i) | field containing the normals in nodes on the elements they belong to |
knor(i) | field containing the extra nodes needed to expand the shell and beam elements to volume elements |
thickn(2,i) | thicknesses (one for shells, two for beams) in node i |
thicke(2,i) | thicknesses (one for shells, two for beams) in element nodes. The entries correspond to the nodal entries in field kon |
offset(2,i) | offsets (one for shells, two for beams) in element i |
iponoel(i) | pointer for node i into field inoel, which stores the 1-D and 2-D elements belonging to the node. |
inoel(3,i) | field containing an element number, a local node number within this element and a pointer to another entry (or zero if there is no other). |
inoelfree | next free field in inoel |
rig(i) | integer field indicating whether node i is a rigid node (nonzero value) or not (zero value). In a rigid node or knot all expansion nodes except the ones not in the midface of plane stress, plane strain and axisymmetric elements are connected with a rigid body MPC. If node i is a rigid node rig(i) is the number of the rotational node of the knot; if the node belongs to axisymmetric, plane stress and plane strain elements only, no rotational node is linked to the knot and rig(i)=-1 |
AMPLITUDES | |
nam | # amplitude definitions |
amta(1,j) | time of (time,amplitude) pair j |
amta(2,j) | amplitude of (time,amplitude) pair j |
namtot | total # of (time,amplitude) pairs |
For amplitude i | |
amname(i) | name of the amplitude |
namta(1,i) | location of first (time,amplitude) pair in field amta |
namta(2,i) | location of last (time,amplitude) pair in field amta |
namta(3,i) | in absolute value the amplitude it refers to; if
abs(namta(3,i))=i it refers to itself. If abs(namta(3,i))=j, amplitude i is a
time delay of amplitude j the value of which is stored in
amta(1,namta(1,i)); in the latter case amta(2,namta(1,i)) is without
meaning; If namta(3,i)![]() |
TRANSFORMS | |
ntrans | # transform definitions |
trab(1..6,i) | coordinates of two points defining the transform |
trab(7,i) | =1 for rectangular transformations |
=2 for cylindrical transformations | |
SINGLE POINT CONSTRAINTS | |
nboun | # SPC's |
For SPC i | |
nodeboun(i) | SPC node |
ndirboun(i) | SPC direction |
typeboun(i) | SPC type (SPCs can contain the nonhomogeneous part of MPCs) |
B=prescribed boundary condition | |
M=midplane | |
P=planempc | |
R=rigidbody | |
S=straigthmpc | |
U=usermpc | |
xboun(i) | magnitude of constraint at end of a step |
xbounold(i) | magnitude of constraint at beginning of a step |
xbounact(i) | magnitude of constraint at the end of the present increment |
xbounini(i) | magnitude of constraint at the start of the present increment |
iamboun(i) | amplitude number |
ikboun(i) | ordered array of the DOFs corresponding to the SPC's (DOF=7*(nodeboun(i)-1)+ndirboun(i)) |
ilboun(i) | original SPC number for ikboun(i) |
MULTIPLE POINT CONSTRAINTS | |
j=ipompc(i) | starting location in nodempc and coefmpc of MPC i |
nodempc(1,j) | node of first term of MPC i |
nodempc(2,j) | direction of first term of MPC i |
k=nodempc(3,j) | next entry in field nodempc for MPC i (if zero: no more terms in MPC) |
coefmpc(j) | first coefficient belonging to MPC i |
nodempc(1,k) | node of second term of MPC i |
nodempc(2,k) | direction of second term of MPC i |
coefmpc(k) | coefficient of second term of MPC i |
ikmpc (i) | ordered array of the dependent DOFs corresponding to the MPC's DOF=7*(nodempc(1,ipompc(i))-1)+nodempc(2,ipompc(i)) |
ilmpc (i) | original MPC number for ikmpc(i) |
icascade | 0 : MPC's did not change since the last iteration |
1 : MPC's changed since last iteration : dependency check in cascade.c necessary | |
2 : at least one nonlinear MPC had DOFs in common with a linear MPC or another nonlinear MPC. dependency check is necessary in each iteration | |
POINT LOADS | |
nforc | # of point loads |
For point load i | |
nodeforc(1,i) | node in which force is applied |
nodeforc(2,i) | sector number (only for modal dynamics and steady state dynamics analyses with cyclic symmetry) |
ndirforc(i) | direction of force |
xforc(i) | magnitude of force at end of a step |
xforcold(i) | magnitude of force at start of a step |
xforcact(i) | actual magnitude |
iamforc(i) | amplitude number |
ikforc(i) | ordered array of the DOFs corresponding to the point loads (DOF=7*(nodeboun(i)-1)+ndirboun(i)) |
ilforc(i) | original SPC number for ikforc(i) |
FACIAL DISTRIBUTED LOADS | |
nload | # of facial distributed loads |
For distributed load i | |
nelemload(1,i) | element to which distributed load is applied |
nelemload(2,i) | node for the environment temperature (only for heat transfer analyses); sector number (only for modal dynamics and steady state dynamics analyses with cyclic symmetry) |
sideload(i) | load label; indicated element side to which load is applied |
xload(1,i) | magnitude of load at end of a step or, for heat transfer analyses, the convection (*FILM) or the radiation coefficient (*RADIATE) |
xload(2,i) | the environment temperature (only for heat transfer analyses |
xloadold(1..2,i) | magnitude of load at start of a step |
xloadact(1..2,i) | actual magnitude of load |
iamload(1,i) | amplitude number for xload(1,i) |
iamload(2,i) | amplitude number for xload(2,i) |
MASS FLOW RATE | |
nflow | # of fluid elements |
TEMPERATURE LOADS | |
t0(i) | initial temperature in node i at the start of the calculation |
t1(i) | temperature at the end of a step in node i |
t1old(i) | temperature at the start of a step in node i |
t1act(i) | actual temperature in node i |
iamt1(i) | amplitude number |
MECHANICAL BODY LOADS | |
nbody | # of mechanical body loads |
For body load i | |
ibody(1,i) | code identifying the kind of body load |
1: centrifugal loading | |
2: gravity loading with known gravity vector | |
3: generalized gravity loading | |
ibody(2,i) | amplitude number for load i |
cbody(i) | element number or element set to which load i applies |
xbody(1,i) | size of the body load |
xbody(2..4,i) | for centrifugal loading: point on the axis |
for gravity loading with known gravity vector: normalized gravity vector | |
xbody(5..7,i) | for centrifugal loading: normalized vector on the rotation axis |
xbodyact(1,i) | actual magnitude of load |
xbodyact(2..7,i) | identical to the corresponding entries in xbody |
STRESS AND STRAIN FIELDS | |
eei(i,j,k) | in general : Lagrange strain component i in integration point j of element k (linear strain in linear elastic calculations) |
for elements with &sstarf#star;DEFORMATION PLASTICITY property: Eulerian strain component i in integration point j of element k (linear strain in linear elastic calculations) | |
eeiini(i,j,k) | Lagrange strain component i in integration point of element k at the start of an increment |
een(i,j) | Lagrange strain component i in node j (mean over all adjacent elements linear strain in linear elastic calculations) |
stx(i,j,k) | Cauchy or PK2 stress component i in integration point j of element k at the end of an iteration (linear stress in linear elastic calculations) |
sti(i,j,k) | PK2 stress component i in integration point j of element k at the start of an iteration (linear stress in linear elastic calculations) |
stiini(i,j,k) | PK2 stress component i in integration point j of element k at the start of an increment |
stn(i,j) | Cauchy stress component i in node j (mean over all adjacent elements; "linear" stress in linear elastic calculations) |
THERMAL ANALYSIS | |
ithermal | 0 : no temperatures involved in the calculation |
1 : stress analysis with given temperature field | |
2 : thermal analysis (no displacements) | |
3 : coupled thermal-mechanical analysis : temperatures and displacements are solved for simultaneously | |
v(0,j) | temperature of node j at the end of
an iteration (for ithermal ![]() |
vold(0,j) | temperature of node j at the start
of an iteration (for ithermal ![]() |
vini(0,j) | temperature of node j at the start
of an increment (for ithermal ![]() |
fn(0,j) | actual temperature at node j (for ithermal ![]() |
qfx(i,j,k) | heat flux component i in integration point j of element k at the end of an iteration |
qfn(i,j) | heat flux component i in node j (mean over all adjacent elements) |
DISPLACEMENTS AND SPATIAL/TIME DERIVATIVES | |
v(i,j) | displacement of node j in direction i at the end of an iteration |
vold(i,j) | displacement of node j in direction i at the start of an iteration |
vini(i,j) | displacement of node j in direction i at the start of an increment |
ve(i,j) | velocity of node j in direction i at the end of an iteration |
veold(i,j) | velocity of node j in direction i at the start of an iteration |
veini(i,j) | velocity of node j in direction i at the start of an increment |
accold(i,j) | acceleration of node j in direction i at the start of an iteration |
accini(i,j) | acceleration of node j in direction i at the start of an increment |
vkl(i,j) | (i,j) component of the displacement gradient tensor at the end of an iteration |
xkl(i,j) | (i,j) component of the deformation gradient tensor at the end of an iteration |
xikl(i,j) | (i,j) component of the deformation gradient tensor at the start of an increment |
ckl(i,j) | (i,j) component of the inverse of the deformation gradient tensor |
LINEAR EQUATION SYSTEM | |
ad(i) | element i on diagonal of stiffness matrix |
au(i) | element i in upper triangle of stiffness matrix |
adb(i) | element i on diagonal of mass matrix, or, for buckling, of the incremental stiffness matrix (only nonzero elements are stored) |
aub(i) | element i in upper triangle of mass matrix, or, for buckling, of the incremental stiffness matrix (only nonzero elements are stored) |
neq[0] | # of mechanical equations |
neq[1] | sum of mechanical and thermal equations |
neq[2] | neq[1] + # of single point constraints (only for modal calculations) |
nzl | number of the column such that all columns with
a higher column number do not contain any
(projected) nonzero off-diagonal terms (![]() |
nzs[0] | sum of projected nonzero mechanical off-diagonal terms |
nzs[1] | nzs[0]+sum of projected nonzero thermal off-diagonal terms |
nzs[2] | nzs[1] + sum of nonzero coefficients of SPC degrees of freedom (only for modal calculations) |
nactdof(i,j) | actual degree of freedom (in the system of equations) of DOF i of node j (0 if not active) |
INTERNAL AND EXTERNAL FORCES | |
fext(i) | external mechanical forces in DOF i (due to point loads and distributed loads, including centrifugal and gravity loads, but excluding temperature loading and displacement loading) |
fextini(i) | external mechanical forces in DOF i (due to point loads and distributed loads, including centrifugal and gravity loads, but excluding temperature loading and displacement loading) at the end of the last increment |
finc(i) | external mechanical forces in DOF i augmented by contributions due to temperature loading and prescribed displacements; used in linear calculations only |
f(i) | actual internal forces in DOF i due to : |
actual displacements in the independent nodes; | |
prescribed displacements at the end of the increment in the dependent nodes; | |
temperatures at the end of the increment in all nodes | |
fini(i) | internal forces in DOF i at the end of the last increment |
b(i) | right hand side of the equation system : difference between fext and f in nonlinear calcultions; for linear calculations, b=finc. |
fn(i,j) | actual force at node j in direction i |
INCREMENT PARAMETERS | |
tinc | user given increment size (can be modified by the program if the parameter DIRECT is not activated) |
tper | user given step size |
dtheta | normalized (by tper) increment size |
theta | normalized (by tper) size of all previous increments (not including the present increment) |
reltime | theta+dtheta |
dtime | real time increment size |
time | real time size of all previous increments INCLUDING the present increment |
DIRECT INTEGRATION DYNAMICS | |
alpha,bet,gam | parameter in the alpha-method of Hilber, Hughes and Taylor |
iexpl | =0 : implicit dynamics |
=1 : explicit dynamics | |
FREQUENCY CALCULATIONS | |
mei[0] | number of requested eigenvalues |
mei[1] | number of Lanczos vectors |
mei[2] | maximum number of iterations |
mei[3] | if 1: store eigenfrequencies, eigenmodes, mass matrix and possibly stiffness matrix in .eig file, else 0 |
fei[0] | tolerance (accuracy) |
fei[1] | lower value of requested frequency range |
fei[2] | upper value of requested frequency range |
CYCLIC SYMMETRY CALCULATIONS | |
mcs | number of cyclic symmetry parts |
ics | one-dimensional field; contains all independent nodes, one part after the other, and sorted within each part |
rcs | one-dimensional field; contains the corresponding radial coordinates |
zcs | one-dimensional field; contains the corresponding axial coordinates |
For cyclic symmetry part i | |
cs(1,i) | number of segments in ![]() |
cs(2,i) | minimum nodal diameter |
cs(3,i) | maximum nodal diameter |
cs(4,i) | number of nodes on the independent side |
cs(5,i) | number of sections to be plotted |
cs(6..12,i) | coordinates of two points on the cyclic symmetry axis |
cs(13,i) | number of the element set (for plotting purposes) |
cs(14,i) | total number of independent nodes in all previously defined cyclic symmetry parts |
cs(15,i) | cos(angle) where angle = 2*![]() |
cs(16,i) | sin(angle) where angle = 2*![]() |
cs(17,i) | number of tie constraint |
OUTPUT IN .DAT FILE | |
prset(i) | node or element set corresponding to output request i |
prlab(i) | label corresponding to output request i. It contains 6 characters. The first 4 are reserved for the field name, e.g. 'U ' for displacements, the fifth for the value of the TOTALS parameter ('T' for TOTALS=YES, 'O' for TOTALS=ONLY and ' ' else) and the sixth for the value of the GLOBAL parameter ('G' for GLOBAL=YES and 'L' for GLOBAL=NO). |
nprint | number of print requests |
OUTPUT IN .FIL FILE | |
filab(i) | label corresponding to output field i. It contains 6 characters. The first 4 are reserved for the field name. The order is fixed: filab(1)='U ', filab(2)='NT ',filab(3)='S ',filab(4)='E ', filab(5)='RF ', filab(6)='PE ', filab(7)='ENER', filab(8)='SDV ', filab(9)='HFL ', filab(10)='RFL ', filab(11)='PU ' and filab(12)='PNT '. Results are stored for the complete mesh. A field is not selected if the first 4 characters are blank, e.g. the stress is not stored if filab(3)(1:4)=' '. If the mesh contains 1D or 2D elements, the fifth character takes the value 'I' if the results are to be interpolated, 'M' if the section forces are requested instead of the stresses and 'E' if the 1D/2D element results are to be given on the expanded elements. In all other cases the fifth character is blank: ' '. The sixth character contains the value of the GLOBAL parameter ('G' for GLOBAL=YES and 'L' for GLOBAL=NO). |