| Cols. | 1-5 NPOIN | Total number of nodal points. |
| 6-10 NELEM | Total number of elements. |
| 11-15 NVFIX | Total number of restrained boundary points—where one or more degrees of freedom are restrained. |
| 16-20 NTYPE | Problem type parameter: |
| | 1—Plane stress, |
| | 2—Plane strain, |
| | 3—Axial symmetry. |
| 21-25 NNODE | Number of nodes per element: |
| | 4—Linear quadrilateral element, |
| | 8—Quadratic Serendipity element, |
| | 9—Quadratic Lagrangian element. |
| 26-30 NMATS | Total number of different materials. |
| 31-35 NGAUS | Order of integration formula for numerical integration: |
| | 2—Two point Gauss quadrature rule, |
| | 3—Three point Gauss quadrature rule. |
36-40 NALGO
Nonlinear solution parameter:
1 Initial stiffness method. The element stiffnesses are calculated at the beginning of the solution process and remain unchanged thereafter.
2 Tangential stiffness method. The element stiffnesses are recalculated for every iteration of each load increment.
3 Combined algorithm (Version I). The element stiffnesses are recalculated for the first iteration of each load increment only.
4 Combined algorithm (Version II). The element stiffnesses are recalculated for the second iteration of each load increment only.
41-45 NCRIT
Yield criterion parameter:
1—Tresca,
2—Von Mises,
3—Mohr-Coulomb,
4—Drucker–Prager.
46-50 NINCS
Number of increments in which the total loading is to be applied.
51-55 NSTRE
Number of stress components at a point:
3—Plane stress or plane strain,
4—Axial symmetry.
CARD SET 3 ELEMENT CARDS (1115)—One card for each element. Total of NELEM cards (See Card Set 2).
Cols. 1-5 NUMEL
Element number.
6-10 MATNO(NUMEL) Material property number.
11-15 LNODS(NUMEL,1) 1st Nodal connection number.
16-20 LNODS(NUMEL,2) 2nd Nodal connection number.
51-55 LNODS(NUMEL,9) 9th Nodal connection number.
Notes: 1) Columns 31–55 remain blank for linear 4-noded elements.
2) Columns 51–55 remain blank for 8-noded elements.
3) The nodal connection numbers must be listed in an anti-clockwise sequence, starting from any corner node.
CARD SET 4 NODE CARDS (I5,2F10.5)—One card for each node whose coordinates are to be input.
Cols. 1-5 IPOIN
Nodal point number.
6-15 COORD(IPOIN,1)
x (or r) coordinate of the node.
16-25 COORD(IPOIN,2)
$y$ (or $z$ ) coordinate of the node.
Notes: 1) The total number of cards in this set will generally differ from NPOIN (see Card Set 2) since for quadratic elements whose sides are linear, it is only necessary to specify data for corner nodes, intermediate nodal coordinates being automatically interpolated if on a straight line.
2) For Lagrangian elements the coordinates of the 9th (central) node are never input.
3) The coordinates of the highest numbered node must be input regardless of whether it is a midside node or not.
CARD SET 5 RESTRAINED NODE CARDS (1X,14,5X,15,5X,2F10.5)—One card for each restrained node. Total of NVFIX cards (See Card Set 2).
Cols. 2-5 NOFIX(IVFIX)
11-15 IFPRE
Restrained node number.
Restraint code:
01 Nodal displacement restrained in the $x$ (or $r$ ) direction,
10 Nodal displacement restrained in the $y$ (or $z$ ) direction,
11 Nodal displacement restrained in both coordinate directions.
21-30 PRESC(IVFIX,1)
The prescribed value of the x (or r) component of nodal displacement.
31-40 PRESC(IVFIX,2)
The prescribed value of the y (or z) component of nodal displacement.
CARD SET 6 MATERIAL CARDS
6(a) CONTROL CARD (I5)—One card.
Cols. 1-5 NUMAT
Material identification number.
6(b) PROPERTIES CARDS (7F10.5)—One card for each different material.
Cols. 1-10 PROPS(NUMAT,I)
Elastic modulus, E.
11-20 PROPS(NUMAT,2)
Poisson's ratio, v.
21-30 PROPS(NUMAT,3)
Material thickness, t (leave blank for plane strain and axisymmetric problems).
31-40 PROPS(NUMAT,4)
Mass density, ρ.
41-50 PROPS(NUMAT,5)
Uniaxial yield stress, $\sigma_{Y}$ (or cohesion $c$ for Mohr-Coulomb or Drucker-Prager materials).
51-60 PROPS(NUMAT,6)
Strain hardening parameter, $H'$ .
61-70 PROPS(NUMAT,7) Friction angle $\phi$ (measured in degrees) for Mohr-Coulomb and Drucker-Prager materials only).
Note: This card set to be repeated for each different material. Total of NMATS card sets (See Card Set 2).
CARD SET 7 LOAD CASE TITLE CARD (12A6)—One card.
Cols. 1–72 TITLE Title of the load case—limited to 72 alphanumeric characters.
CARD SET 8 LOAD CONTROL CARD (315)—One card.
Cols. 1-5 IPLOD Applied point load control parameter: 0 No applied nodal loads to be input, 1 Applied nodal loads to be input.
6-10 IGRAV Gravity loading control parameter:
0 No gravity loads to be considered,
1 Gravity loading to be considered.
11-15 IEDGE Distributed edge load control parameter: 0 No distributed edge loads to be input, 1 Distributed edge loads to be input.
CARD SET 9 APPLIED LOAD CARDS (15,2F10.3)—One card for each loaded nodal point.
Cols. 1-5 LODPT Node number.
6-15 POINT(1) Load component in $x$ (or $r$ ) direction.
16-25 POINT(2) Load component in $y$ (or $z$ ) direction.
Notes: 1) The last card should be that for the highest numbered node whether it is loaded or not.
2) For axisymmetric problems, the loads input should be the total loading on the circumferential ring passing through the nodal point concerned.
3) If IPLOD = 0 in Card Set 8, omit this set.
CARD SET 10 GRAVITY LOADING CARD (2F10.3)—One card.
Cols. 1–10 THETA Angle of gravity axis measured from the positive y axis (see Fig. 6.7).
11-20 GRAVY Gravity constant—specified as a multiple of the gravitational acceleration, g.
Note: If IGRAV = 0 in Card Set 8, omit this set.
CARD SET 11 DISTRIBUTED EDGE LOAD CARDS
11(a) CONTROL CARD (I5)—One card.
Cols. 1-5 NEDGE Number of element edges on which distributed loads are to be applied.
# 11(b) ELEMENT FACE TOPOLOGY CARD (415)
Cols. 1–5 NEASS
$$
\left. \begin{array}{l l} 6 - 1 0 & \text { NOPRS } (1) \\ 1 1 - 1 5 & \text { NOPRS } (2) \\ 1 6 - 2 0 & \text { NOPRS } (3) \end{array} \right\}
$$
The element number with which the element edge is associated.
List of nodal points, in an anticlockwise sequence, of the nodes forming the element face on which the distributed load acts.
Note: For linear 4-noded elements, Cols. 16–20 remain blank.
# 11(c) DISTRIBUTED LOAD CARDS (6F10.3)
Cols. 1-10 PRESS(1,1)
11-20 PRESS(1,2)
21-30 PRESS(2,1)
31-40 PRESS(2,2)
41-50 PRESS(3,1)
51-60 PRESS(3,2)
Value of normal component of distributed load at node NOPRS(1).
Value of tangential component of distributed load at node NOPRS(1).
Value of normal component of distributed load at node NOPRS(2).
Value of tangential component of distributed load at node NOPRS(2).
Value of normal component of distributed load at node NOPRS(3).
Value of tangential component of distributed load at node NOPRS(3).
Notes: 1) For linear 4-noded elements, Cols. 41–60 remain blank.
2) Subsets 11(b) and 11(c) must be repeated in turn for every element edge on which a distributed load acts. The element edges can be considered in any order.
3) If IEDGE = 0 in Card Set 8, omit this card set.
# CARD SET 12 LOAD INCREMENT CONTROL CARDS (2F10.5,3I5)—One card for each load increment. Total of NINCS cards (see Card Set 2).
Cols. 1-10 FACTO
11-20 TOLER
21-25 MITER
26-30 NOUTP(1)
Applied load factor for this increment—specified as a factor of the loading input in Card Sets 8 to 11.
Convergence tolerance factor.—The term TOLER in (3.27).
Maximum number of iterations allowed for the load increment.
Parameter controlling output of results after 1st iteration:
0—No output,
1—Output displacements,
2—Output displacements and reactions,
31-35 NOUTP(2)
3—Output displacements, reactions and stresses.
Parameter controlling output of the converged results:
0—No output,
1—Output displacements,
2—Output displacements and reactions,
3—Output displacements, reactions and stresses.
Note: The applied loading factors are accumulative. If FACTO is specified as 0.6, 0.3, 0.2 for the first three load increments, then the total loading acting during the third increment is 1.1 times that specified in Card Sets 8 to 11.
# A.2.2 Program VISCOUNT for the elasto-viscoplastic analysis of plane and axisymmetric solids
The input data for this application is identical to that described in Section A.2.1, for elasto-plastic problems, with the following exceptions.
# CARD SET 2 CONTROL CARD (1115)
Cols. 36–40 NALGO
Equation solution parameter:
1 Explicit time stepping scheme (i.e. TIMEX = 0—See Card Set 12),
2 Implicit or Semi-implicit schemes (TIMEX ≠ 0).
CARD SET 6(b) PROPERTIES CARDS (8F10.5)—Two cards for each different material.
1st Card
Cols. 1–70
Identical to Card Set 6(b), Section A.2.1.
71-80 PROPS(NUMAT,8)
Fluidity parameter, γ.
2nd Card
Cols. 1-10 PROPS(NUMAT,9)
The constant $M$ in (8.8) or constant $N$ in (8.9).
11-20 PROPS(NUMAT,10)
Parameter controlling choice of the flow function:
0 Expression (8.8) to be used,
1 Expression (8.9) to be used.
CARD SET 12 TIMESTEPPING PARAMETER CARD (4F10.3)—One card.
Cols. 1-10 TIMEX
Timestepping algorithm parameter, $\Theta$ in (8.10).
| Cols.16-20 | NTYPE | Problem type parameter:5—for Heterosis element,0—for 4- or 8-node elements. |
| 21-25 | NNODE | Number of nodes per element:4—Linear 4-node quadrilateral element.8—Quadratic 8-node Serendipity element.9—Quadratic 9-node Lagrangian element or Heterosis element. |
| 31-35 | NGAUS | 2 for 4-node element,3 for 8-, 9-node and Heterosis element.(N.B. This is the integration rule to evaluate the flexural contribution to the element stiffness matrix. Since selective integration is adopted a (NGAUS-1) integration is automatically used to evaluate the transverse shear contribution to the element stiffness matrix.) |
| 41-45 | NCRIT | Yield criterion parameter:1—Tresca,2—Von-Mises.(Mohr-Coulomb and Drucker-Prager yield criteria are not included.) |
| 51-55 | NLAPS | Total number of layers.(for program MINDLAY only—in program MINDLIN leave blank.) |
CARD SET 5 RESTRAINED NODE CARDS (1X, I4, 5X, I5, 5X, 3F10.5) One card for each restrained node. Total of NVFIX cards.
Cols.11-15 IFPRE
Restraint code:
100 Lateral displacement w restrained.
010 Rotation $\theta_{x}$ restrained.
001 Rotation $\theta_{y}$ restrained.
110 Lateral displacement w and rotation $\theta_{x}$ restrained, etc.
21-30 PRESC(IVFIX,1)
The prescribed value of the lateral nodal displacement w.
31-40 PRESC(IVFIX,2)
The prescribed value of the nodal rotation $\theta_{x}$ .
41-50 PRESC(IVFIX.3)
The prescribed value of the nodal rotation $\theta_{y}$ .
# CARD SET 6 MATERIAL CARDS
6(b) PROPERTIES CARDS (7F10.5)—One card for each different material.
Cols.31-40 PROPS(NUMAT,4)
Uniform distributed loading value.
41-50 PROPS(NUMAT,5)
Blank.
51-60 PROPS(NUMAT,6)
Uniaxial yield stress, $\sigma_0$ .
61-70 PROPS(NUMAT,7)
Strain hardening parameter $H'$ .
# CARD SET 6X CONVERGENCE CHECK CARDS
6X(a) DISPLACEMENT CHECK CARD (511)—One card.
Cols. 1 IFDIS
1 The displacement check is to be employed.
2 NCDIS(1)
1 Check based on norm involving w.
3 NCDIS(2)
1 Check based on norm involving $\theta_{x}$ .
4 NCDIS(3)
1 Check based on norm involving $\theta_y$ .
5 NCDIS(4)
1 Check based on $w$ , $\theta_x$ and $\theta_y$ .
6X(b) RESIDUAL FORCE CHECK
CARD (511)—One card.
Cols. 1 IFRES
1 The residual force check is to be employed.
2 NCRES(1)
1 Check based on norm involving residual forces associated with w.
3 NCRES(2)
1 Check based on norm involving residual forces associated with $\theta_{x}$ .
4 NCRES(3)
1 Check based on norm involving residual forces associated with $\theta_{y}$ .
5 NCRES(4)
1 Check based on norm involving residual forces associated with $\dot{w}$ , $\theta_{x}$ and $\theta_{y}$ .
Note: A zero value for any item implies that the check is not being used.
CARD SET 8 LOAD CONTROL CARD (I5)—One card.
Cols. 1–5 IPLOD
Applied point load control parameter:
0 No applied nodal loads to be input.
1 Applied nodal loads to be input.
6-15
Blank.
CARD SET 9 APPLIED LOAD CARDS (I5, 3F10.3)—One card for each loaded nodal point.
Cols. 1-5 LODPT
Node number.
6-15 POINT(1)
Lateral nodal load.
16-25 POINT(2)
Nodal couple in xz plane.
26-35 POINT(3)
Nodal couple in yz plane.
Omit CARD SETS 10, 11(a), 11(b) and 11(c).