# Appendix III # Instructions for preparing input data for dynamic transient problems The program DYNPAK has been described in Section 10.6 and MIXDYN in Section 11.5. These programs perform large displacement or viscoplastic or elasto-plastic, transient dynamic analysis of plane stress/strain or axisymmetric problems respectively. The format of the input data is identical for both programs. In this appendix user instructions for preparing input data are provided. CARD SET 1 DYNAMIC DIMENSIONING (4I5)—One card. Cols. 1-5 NPOIN Total number of nodal points. 6-10 NELEM Total number of elements. 11-15 NDOFN Number of degrees of freedom per node $(= 2)$ . 16-20 NMATS Number of different material sets. CARD SET 2 TITLE CARD (10A4)—One card. Cols. 1–40 Title of the problem—limited to 40 alphanumeric characters. CARD SET 3 CONTROL CARD (13I5)—One card. Cols. 1–5 NVFIX Total number of nodal points with fixed degrees of freedom. $$ \begin{array}{l} = 1, \text { Plane stress }, \\ = 2, \text { Plane strain }, \\ = 3, \text { Axisymmetric problem. } \\ \end{array} $$ 6-10 NTYPE Type of problem: 11-15 NNODE Number of nodes per element. 16-20 NPROP Number of material properties $(= 11)$ . 21-25 NGAUS Integration rule for stiffness matrix. 26-30 NDIME Number of coordinate dimensions (=2). 31-35 NSTRE Number of stress components (= 3 for plane stress/strain, = 4 for axisymmetric).
36-40NCRITYield criterion: = 1 — Tresca, = 2 — Von Mises, = 3 — Mohr-Coulomb, = 4 — Drucker-Prager.
41-45NPREVIndicator for the previous state to be read (= 1 for previous state, otherwise, = 0).
46-50NCONMNumber of concentrated masses (≥1 if concentrated mass present, otherwise, = 0).
51-55NLAPSIndicator for large displacement analysis: = 0—Elastic analysis, = 1—Elasto-plastic small displacement analysis, = 2—Elastic large displacement analysis,
56-60NGAUMIntegration rule for mass matrix.
61-65NRADS= 0, Read (r, z) coordinates for nodes, = 1, Read (R, Θ) coordinates for nodes for axisymmetric analysis.
CARD SET 4 ELEMENT CARDS (11I5)—One card for each element, total of NELEM cards. The node numbers are read in anticlockwise sequence. The number of nodes depends upon the type of element. For four and eight noded elements read only four and eight nodes respectively.
Cols.1-5IELEMElement number.
6-10MATNOMaterial identification number.
11-15LNODS(IELEM,1)
16-20LNODS(IELEM,2)
21-25LNODS(IELEM,3)
26-30LNODS(IELEM,4)
31-35LNODS(IELEM,5)Nodal connection numbers.
36-40LNODS(IELEM,6)
41-45LNODS(IELEM,7)
46-50LNODS(IELEM,8)
51-55LNODS(IELEM,9)
CARD SET 5 NODAL COORDINATE CARDS (I5,2F10.5)—One card for each node. Last nodal point (IPOIN=NPOIN) must be read at the end. Only corner and central nodes need to be specified. Midside nodes are interpolated if not specified. For axisymmetric cases, $(R, \Theta)$ values are read for NRADS = 1, and $(r, z)$ coordinates are calculated in the program. Cols. 1-5 IPOIN Current nodal point. 6-15 COORD(IPOIN,1) x-coordinate.\* 16-25 COORD(IPOIN,2) y-coordinate. CARD SET 6 RESTRAINED NODE CARDS (1X,14,3X,211)—One card for each restrained node. Total of NVFIX cards. Cols. 2-5 IPOIN Restrained node number. 9 IFPRE(IVFIX,1) Fixity in $x$ -direction $(= 0$ , Free; $= 1$ , Fixed). 10 IFPRE(IVFIX,2) Fixity in y-direction (= 0, Free; = 1, Fixed). CARD SET 7 MATERIAL CARDS—Three cards for each different material, a total of NMATS\*3 cards. 1st Card MATERIAL IDENTIFICATION CARD (I5) Cols. 1-5 NUMAT Material identification number. 2nd Card MATERIAL PROPERTIES CARD—(a) (8E10.4) Cols. 1-10 PROPS(NUMAT,1) Young's Modulus, $E$ . 11-20 PROPS(NUMAT,2) Poisson's ratio, $\nu$ . 21-30 PROPS(NUMAT,3) Thickness for plane stress problem, t. 31-40 PROPS(NUMAT,4) Mass density per unit volume, $\rho$ . 41-50 PROPS(NUMAT,5) Temperature coefficient, $\alpha_{t}$ . 51-60 PROPS(NUMAT,6) Reference yield value $F_{0}$ : Von Mises, $F_0 = \sigma Y,$ Tresca, $F_{0} = \sigma_{Y},$ Mohr-Coulomb, $F_0 = c\cos \phi$ Drucker-Prager, $F_{0} = 6c\cos \phi /$ $(\sqrt{3}(3 - \sin \phi)).$ 61-70 PROPS(NUMAT,7) Hardening parameter, $H'$ : $$ H ^ {\prime} = \frac {E _ {T}}{1 - E _ {T} / E}, $$ where $E_{T}$ is the hardening tangent modulus, $E$ is the tangent modulus, $\sigma_{Y}$ is the yield stress, $c$ is the cohesion, $\phi$ is the friction angle. 71-80 PROPS(NUMAT,8) Friction angle ‘ $\phi$ ’. 3rd Card MATERIAL PROPERTIES CARD—(b) (3E10.4) Cols. 1–10 PROPS(NUMAT,9) Fluidity parameter, γ. 11-20 PROPS(NUMAT,10)Exponent, $\delta$ . 21-30 PROPS(NUMAT,11)NFLOW code (NFLOW = 1—Power law, NFLOW ≠ 1—Exponential law). CARD SET 8 TIME INTEGRATION CONTROL CARD (1115)—One card. Cols. 1–5 NSTEP Total number of time steps. 6-10 NOUTD Writes displacement and stress history of required points on tapes 10 and 11 respectively at NOUTD timesteps. 11-15 NOUTP Output for displacements and stresses at every NOUTP step (NOUTP $\leqslant 500$ ). 16-20 NREQD Number of nodes for selective output of displacements at NOUTD steps. 21-25 NREQS Number of integration points for selective output of stresses at every NOUTP step. 26-30 NACCE Number of acceleration ordinates (If IFUNC $\neq 0$ , NACCE is not used, then leave blank). 31-35 IFUNC Time function code: IFUNC = 0 Acceleration time history, IFUNC = 1 Heaviside function, $f(t) = 1.0$ , IFUNC = 2 Harmonic excitation, $f(t)$ $= a_{0} + b_{0} \sin \omega t.$ 36-40 IFIXD Indicator for excitation: IFIXD = 0, Horizontal acceleration read from tape 7, Vertical acceleration read from tape 12. IFIXD = 1, Vertical acceleration read from tape 12, IFIXD = 2, Horizontal acceleration read from tape 7. (If IFUNC ≠ 0 IFIXD is not used, then leave blank.) 41-45 MITER Maximum number of iterations. This variable is not used in DYNPAK, so leave blank. 46-50 KSTEP Number of steps after which the stiffness matrix is reformed. Not used in DYN-PAK, leave blank. 51-55 IPRED = 1 Standard algorithm, = 2 Modified algorithm. CARD SET 9 TIME INTEGRATION PARAMETERS CARD (8F10.3)—Two cards. 1st Card Cols. 1-10 DTIME 11-20 DTEND 21-30 DTREC 31-40 AALFA 41-50 BEETA 51-60 DELTA 61-70 GAAMA 71-80 AZERO 2nd Card 1-10 BZERO 11-20 OMEGA 21-30 TOLER Time step length. Time at the end of the excitation force. Time step of acceleration records. $\alpha = \text{Damping} \quad \text{parameter}, \quad C = \alpha M,$ $\alpha = 2\xi_i\omega_i.$ $\beta = \text{Damping parameter}, C = \beta K.$ $(\alpha + \beta \omega_{i}^{2} = 2\omega_{i}\xi_{i}, \text{not used in DYNPAK})$ Newmark's integration parameter $(\delta = 0.25 (\gamma + 0.5)^2$ , not used in DYN-PAK). Newmark's integration parameter ( $\gamma \geqslant 0.5$ for stable solution, not used in DYN-PAK). Constants for harmonic excitation $f(t) = a_{0} + b_{0} \sin \omega t$ . Specified tolerance (Not used in DYN-PAK). CARD SET 10 CARD FOR NODAL POINTS FOR WHICH DISPLACEMENT HISTORY IS REQUIRED (16I5)—Total of NREQD nodes. Cols. 1-5 NPRQD(1) 6-10 NPRQD(2) First nodal point at which displacement history is required. Second nodal point at which displacement history is required. 11-15 . CARD SET 11 CARD FOR INTEGRATION POINTS FOR WHICH STRESS HISTORY IS REQUIRED (16I5)—Total of NREQS integration points. Cols. 1-5 NGRQS(1) First integration point at which stress history is required. 6-10 NGRQS(2) Second integration point at which stress history is required. 11-15 . . . . . . . . . CARD SET 12 IMPLICIT-EXPLICIT ELEMENT INDICATOR CARDS (16I5). Number of cards depends on number of elements. For each 16 elements one card is needed. In DYNPAK, INTGR(IELEM) is 2 for every element. INTGR(IELEM) = 1, Implicit element. INTGR(IELEM) = 2, Explicit element. CARD SET 13 INITIAL DISPLACEMENT CARDS (I5,2F10.5)—One card for each node. If all displacements are zero, read data for last node. Cols. 1–5 NGASH Nodal point. 6-15 XGASH Initial $x$ -displacement. 16-25 YGASH Initial y-displacement. CARD SET 14 INITIAL VELOCITY CARDS (15,2F10.5)—One card for each node. If all velocities are zero, read data for last node. Cols. 1–5 NGASH Nodal point. 6-15 XGASH Initial $x$ -velocity. 16-25 YGASH Initial y-velocity. CARD SET 15 PREVIOUS LOAD STATE CARDS (15,2F10.3)—One card for one node, a total of NNODE cards. Data for the last nodal point should always be read even when it is not loaded. If NPREV = 0 then omit this set of data. Cols. 1-5 NGASH Nodal point. 6-15 XGASH Equivalent nodal load in $x$ direction. 16-25 YGASH Equivalent nodal load in $y$ direction. CARD SET 16 PREVIOUS STRESS STATE CARD (I5,4F10.3)—One card for one integration point. Total of (NELEM\*NGAUS\*NGAUS) cards. If NPREV = 0 omit this set of data.
Cols.1-5KGAUSIntegration point.
6-15STRESS(1)Initial stress, $\sigma_x$ or $\sigma_r$ .
16-25STRESS(2)Initial stress, $\sigma_y$ or $\sigma_z$ .
26-35STRESS(3)Initial stress, $\gamma_{xy}$ or $\gamma_{rz}$ .
36-45STRESS(4)Initial stress, $\sigma_z$ or $\sigma_\theta$ .
CARD SET 17 LOAD TITLE CARD (10A4)—One card. Cols. 1–40 Title of load applied—limited to 40 alphanumeric characters. CARD SET 18 LOAD INDICATOR CARD (415)—One card.
Cols.1-5IPLODPoint load indicator.
6-10IGRAVGravity load indicator.
11-15IEDGEEdge load indicator.
16-20ITEMPTemperature load indicator.
CARD SET 19 POINT LOAD CARD (I5,2F10.3)—One card for each node. Data for the last node must be specified at the end. If IPLOD = 0 then omit this set of data.
Cols.1-5LODPTNode number.
6-15POINT(1)Load in x-direction.
16-25POINT(2)Load in y-direction.
CARD SET 20 GRAVITY LOAD CARD (2F10.3)—One card only. If IGRAV = 0 then omit this set of data.
Cols. 1-10 THETAAngle of gravity axis to the positive y axis.
11-20 GRAVYGravity constant.
CARD SET 21 NUMBER OF PRESSURE EDGE CARD (I5)—One card. If IEDGE = 0, then omit card sets 21 and 22. Cols. 1-5 NEDGE Number of loaded edges. CARD SET 22 PRESSURE CARDS—Two cards for each pressure loaded edge. 1st Card PRESSURE NODES CARD (415)—One card for each edge. Total of NEDGE cards.
Cols. 1-5 NEASSElement number with edge load.
Cols. 6-10 NOPRS(1)Edge nodes read in anticlockwise sequence.
11-15 NOPRS(2)
16-20 NOPRS(3)
2nd Card PRESSURE CARD (6F10.3)—One card for each edge. Total of NEDGE cards. A pressure normal to a face is assumed to be positive if it acts in a direction into the element. A tangential load is assumed to be positive if it acts in an anticlockwise direction with respect to the loanedWW positive if it acts in an anticlockwise direction with respect to the loaded element. Cols. 1-10 PRESS(1,1) } Normal component of edge load for each node. 11-20 PRESS(2,1) 21-30 PRESS(3,1) 31-40 PRESS(1,2) 41-50 PRESS(2,2) 51-60 PRESS(3,2) } Tangential component of edge load for each node. CARD SET 24 TEMPERATURE CARDS (I5, F10.3)—One card for each node. The last card must be for the highest numbered node. If $\text{ITEMP} = 0$ , omit this set of data. Cols. 1–5 NODPT Node number. 6–15 TEMPE Nodal temperature. CARD SET 25 CONCENTRATED MASSES (I5,2F10.3)—One card for each node. Total of NCONM cards. If NCONM = 0, omit this set of data. Cols. 1–5 IPOIN Current nodal point with concentrated mass. 6–15 XCMAS Concentrated mass associated with the x-direction. 16–25 YCMAS Concentrated mass associated with the y-direction. # Appendix IV # Sample input data and line printer output for one – and two-dimensional applications In this appendix input data and line printer output are provided for a selection of the numerical examples presented in the text. This information will be of assistance to readers who wish to implement the programs contained in the book on their own computer. For economy of space, presentation is limited to one example from each area of application. Also in some cases the line printer output is edited for the same reason. # A.4.1 Solution of one-dimensional quasiharmonic problem by direct iteration. Example of Section 3.9.3, Fig. 3.3 Input data
1-D QUASIHARMONIC EXAMPLE, SECTION 3.9.3, FIG. 3.3
11102112111
110.0
1121
2231
3341
4451
5561
6671
7781
8891
99101
1010111
10.0
21.0
32.0
43.0
54.0
65.0
76.0
87.0
98.0
109.0
1110.0
110.0
1111.0
100.00.0
2011.00.5
Line printer output
1-D QUASIHARMONIC EXAMPLE, SECTION 3.9.3, FIG. 3.3
NPOIN = 11 NELEM = 10 NBOUN = 2 NMATS = 1
NPROP = 1 NNODE = 2 NINCS = 1 NALGO = 1
NDOFN = 1
MATERIAL PROPERTIES
1 10.00000
EL NODES MAT.
1 1 2 1
2 2 3 1
3 3 4 1
4 4 5 1
5 5 6 1
6 6 7 1
7 7 8 1
8 8 9 1
9 9 10 1
10 10 11 1
NODE COORD.
1 0.00000
2 1.00000
3 2.00000
4 3.00000
5 4.00000
6 5.00000
7 6.00000
8 7.00000
9 8.00000
10 9.00000
11 10.00000
RES. NODE CODE PRES. VALUES
1 1 0.00000
11 1 1.00000
ELEMENT NODAL LOADS
1 0.00000 0.00000
2 0.00000 0.00000
3 0.00000 0.00000
4 0.00000 0.00000
5 0.00000 0.00000
6 0.00000 0.00000
7 0.00000 0.00000
8 0.00000 0.00000
9 0.00000 0.00000
10 0.00000 0.00000
IINCS = 1 NITER = 20 NOUTP = 1 FACTO = 0.100000E 01 TOLER = 0.500000E 00
CONVERGENCE CODE = 1 NORM OF RESIDUAL SUM RATIO = 0.000000E 00
NODE DISPL. REACTIONS
1 0.000000E 00 -0.100000E 01
2 0.100000E 00 0.000000E 00
3 0.200000E 00 0.000000E 00
4 0.300000E 00 0.000000E 00
5 0.400000E 00 0.000000E 00
6 0.500000E 00 0.000000E 00
7 0.600000E 00 0.000000E 00
8 0.700000E 00 0.000000E 00
9 0.800000E 00 0.000000E 00
10 0.900000E 00 0.000000E 00
11 0.100000E 01 0.100000E 01
ELEMENT STRESSES PL.STRAIN
1 0.000000E 00 0.000000E 00
2 0.000000E 00 0.000000E 00
3 0.000000E 00 0.000000E 00
4 0.000000E 00 0.000000E 00
5 0.000000E 00 0.000000E 00