1.0 INPUT
DATA STRUCTURE
Input data are organized into data blocks by means of corresponding macro commands. The macro commands are used to identify each data block. The macro command is always the first line in the data block and may be entered in uppercase or lowercase. Data associated with the data block must follow the macro command. Only data blocks pertinent to the particular analysis to be performed need be entered.
In many cases, the program will assign default values for input parameters not specified by the user. These default values, when applicable, are tabulated with each variable entry description.
Input data on any given data line are entered in a free format mode. In the free format mode, the data are input as a string of numbers separated by a comma, and/or any number of blanks. Each data line is restricted not to exceed a length of 80 characters.
If needed, data strings may be broken onto several data lines. For that purpose, a backslash (" \ ") or a slash (" / ") need be used at the end of each data line to be continued.
Any line with the letter C in column 1 followed by at least a blank, is treated as a comment line and is ignored by the program. Any statement following the symbol '#' or '!' on any given line is also ignored.
DYNAFLOW main input data are contained in a main input data file. The program creates an output file called: name.results, where name is the name defined by the command DEFINE_PROBLEM name = " " in the main input data file (the default is name = dynaflow). In addition the program creates output files as follows:
- TAPE90.name: contains coordinates and nodal connectivities
- TAPE87.name: contains nodal dumps
- TAPE89.name: contains field dumps
- TAPE88.name: contains nodal and field time histories dumps
- TAPE91.name: contains material data dumps
- TAPE93.name.nstep: restart/backup data for step number nstep
To enter input data and run the program using an input file, type "dynaflow.v02" followed by the name of the input data file, and by the requested program size (size in megawords), and/or the requested length of the pointer(s) array (length in kilowords). Also, an option is available to monitor printing to the logfile by selecting the time step printout frequency (step).
shell prompt > dynaflow.v02 filename [-size = n -length = m -step=k > & logfile &]
By default n = 30, m = 50, and k=1. Alternatively, the program size and pointer length may be defined in the file dynaflow.env. The input data file: filename should be prepared according to the user's manual, and contains input commands for DYNAFLOW input and action.
No system of units is assumed by DYNAFLOW. Rather, it is the user's responsibility to specify the input data in a consistent set of units of his/her choice. Overlooking this requirement is a common source of erroneous results.
In order to guide the user, input parameter dimensions are indicated when appropriate. The notation used follows the SI convention, viz., the symbols [L, M, T] are used to indicate: L = length, M = mass, T = time. Thus, for instance:
Symbol |
Description |
Unit |
Example (SI) |
|
|
|
|
x |
coordinate |
L |
m |
m |
mass |
M |
kg |
t |
time |
T |
sec. |
d |
displacement |
L |
m |
v |
velocity |
L / T |
m / sec. |
a |
acceleration |
L / T2 |
m / sec.2 |
|
mass density (per unit volume) |
M / L3 |
kg / m3 |
b |
body force (per unit mass) |
L / T2 |
m / sec.2 |
g |
acceleration of gravity |
L / T2 |
9.81 m / sec.2 |
f |
force |
ML / T2 |
N |
E |
modulus of elasticity |
M / LT2 |
N / m2 (= Pa) |
|
viscosity |
M / LT |
Pa sec. (= Poiseuille) |
k |
hydraulic conductivity |
L / T |
m / sec. |
|
|
|
|
Note/
The SI unit Poiseuille must be distinguished from the Poise (1 Poiseuille = 10 Poise). For instance, water at 20˚C has a viscosity of almost exactly 10-3 Pa sec.
Note Command Description
PROBLEM DEFINITION
(1) DEFINE_PROBLEM Start New Input Record
(2) OPEN_FILE Open a File
CLOSE_FILE Close a File
PARAMETERS Define Parameters
NODAL_COORDINATES Nodal Coordinates
NODAL_BOUNDARY_CONDITIONS Nodal Boundary Conditions
SLAVE_NODES Slave Nodes
EQUIVALENCE_NODES Equivalence Nodes
ALE_NODES Euler-Lagrange Nodes
INITIAL_D0 Initial Nodal Displacement
INITIAL_V0 Initial Nodal Velocity
BACKGROUND_NODAL_FIELD Background Nodal Field
NODAL_LOADS Prescribed Nodal Forces / Kinematics
SURFACE_LOADS Prescribed Surface Forces
CONVECTIVE_SURFACES Convective/Radiative Surfaces
ALE_FREE_SURFACE Mixed Euler-Lagrange Free Surface
READ_FREE_FIELD_MOTION Read Motion Request
FILTER_FREE_FIELD_MOTION Filter Motion Request
FREE_FIELD_NODES Free-Field Nodal Data
LOAD_TIME_FUNCTION Load-Time Function
FILTER Load-Time Function Filtering
DEFINE_REGION Define a Region
(3) DEFINE_ELEMENT_GROUP Define a Group of Elements
ELEMENT_GROUP Element Group Data
DEFINE_MATERIAL_MODEL Define a Material Model
(4) NODAL_CONNECTIVITY Element Nodal Connectivity
LIST OF MACRO
COMMANDS (Continued)
Note Command Description
OUTPUT / PRINT REQUESTS
ECHO Echo Input
PRINT Print Subsequent Input (Default)
NOPRINT Do Not Print Subsequent Input
PRINT_ITERATIONS Iterations Printout
PRINT_PIVOTS Pivots Printout
PRINTOUT Nodal and/or Spatial Printout
NODAL_PRINT Selective Nodal Printout
NODAL_HISTORY Nodal Time History
NODAL_REACTION Nodal Reaction Time History
TRANSLATOR Target Post-Processor
MESH_DUMP Mesh Dump
NODAL_DUMP Selective Nodal Dump
FIELD_DUMP Selective Field Dump
WRITE_MOTION Write Motion
RECOVER_ERROR Error Recovery
RECOVER_STRAIN_ENERGY Strain-Energy Recovery
SYSTEM_COMPLIANCE System Compliance
LIST OF MACRO
COMMANDS (Continued)
Note Command Description
EXECUTION REQUESTS
TIME_SEQUENCE Specify Global Time Stepping
DEFINE_STAGGER Define a Solution Stagger
STAGGER_CONTROL Staggered Solution Control
(5) INITIALIZE_D0 Initialization Request for D0
INITIALIZE_V0 Initialization Request for V0
INITIALIZE_A0 Initialization Request for A0
(6) TIME_INTEGRATION Time-Integration Parameters
LINEAR_SOLVER Linear Equation Solver Selection
NONLINEAR_ITERATIONS Iteration Requests
EIGENVALUE_SOLUTION Eigenvalue Solution
(7) RECOVER_ERROR Error Recovery
LAYOUT_OPTIMIZATION Layout Optimization
(8) CONSTITUTIVE_EXPERIMENT Constitutive Test Requests
BACKUP Backup Request
(9) RUN_SOLVER Run Solver
(10) STOP Stop Program
UPDATE/CLEAR REQUESTS
UPDATE_COORDINATES Update Nodal Coordinate Array
REMESH Remeshing Request
CLEAR_D Clear Nodal Displacement Array
CLEAR_V Clear Nodal Velocity Array
CLEAR_A Clear Nodal Acceleration Array
CLEAR_STRESS Clear Stress Array
CLEAR_T Reset Time to 0.0
LIST OF MACRO COMMANDS (Continued)
Notes/
(1) Define problem command serves to indicate the start of a new problem record.
(2) The input data segments can be put in any order.
(3) The elements may be read in groups (consult Chapter 9 for details).
(4) Element connectivity data must be entered as part of an element group data block (see Chapter 9 for further details).
(5) For certain problems, it is required that an initialization takes place for D0 (e.g., for pressure dependent materials which require that gravity induced initial stresses be first computed). (See Section 12.3.)
(6) Time-stepping and/or eigenvalue or constitutive experiment requests must be provided in order to direct the code toward an appropriate execution mode.
(7) Determine the magnitude of the error at any given time step.
(8) This option allows testing of the material
constitutive modules by prescribing stress or strain paths.
(9) Run_Solver (mode = data_check or execution)
must be the last command of the problem record to direct the solver toward an
appropriate execution mode. In the data check mode, input data are printed out
and storage requirements are indicated. This mode should be employed prior to making
expensive executions.
(10) STOP must be used to indicate the end of the data file (alternatives are QUIT, EXIT or STOP).
DYNAFLOW includes a general element library. The elements are identified by means of corresponding keywords which are used to identify each element data block.
Keyword Description
QDC_solid Solid Continuum (Lagrangian)
QDC_porous Coupled Solid/Fluid Porous Continuum
QDC_pressure Scalar Diffusion Equation
QDC_Darcy Darcy Flow in Porous Continuum
QDC_Stokes Stokes flow
QDC_fluid Fluid Continuum (Eulerian)
QDC_Helmoltz Helmoltz/Laplace equation
QDC_transport Scalar Convection/Advection-Diffusion Transport Equation
QDC_thermal Coupled Solid/Thermal Continuum
QDC_heat Scalar Heat Equation
QDC_ale Mesh Displacement Field
QDC_reservoir Coupled Porous Solid / Two-phase Flow
QDC_flow Multi-Phase Flow Equations
QDC_charge Electric Charge Equation
Multi_phase_transport Multi-Phase Transport Equation
Multi_phase_heat Multi-Phase Heat Equation
Interface_surface Interface w/ Coulomb Friction
Contact_surface Nodal Contact
Slide_line Slide Line
Slide_coulomb Slide Line w/ Coulomb Friction
Crack_xfem Crack/Joint
Multi_point_constraint Multi Point Constraint
Linear_truss Structural Truss (linear)
Nonlinear_truss Structural Truss (nonlinear)
Linear_beam Structural Beam/Frame (linear)
Nonlinear_beam Structural Beam/Frame (nonlinear)
Plate Structural Plate (linear)
Shell_plate Structural Shell (linear)
Shell_bilinear Structural Bilinear Shell (linear)
Membrane Structural Membrane
Nodal_mass Nodal Mass
Nodal_damping Nodal Damping
Nodal_spring Nodal Stiffness
Nodal_transmitting Transmitting Boundary
Nodal_reaction Boundary Element
Nodal_link Link Element
Nodal_penalty Nodal Penalty (selective)
DYNAFLOW includes a general material library. The materials models are identified by means of corresponding keywords which are used to identify each material data block.
Keyword Description
LINEAR_ELASTIC Isotropic linear elastic
ORTHOTROPIC_ELASTIC Orthotropic linear elastic
HYPERELASTIC Hyperelastic
NEWTONIAN_FLUID Newtonian fluid
MISES Mises elasto-(visco)-plastic
DRUCKER_PRAGER Drucker-Prager elasto-(visco)-plastic
MATSUOKA Matsuoka elasto-(visco)-plastic
CAP Cap model
MULTI_YIELD Multi-yield elasto-plastic
ISHIHARA Multi-mechanism elasto-plastic
MCREEP Mises viscoelastic
PHILLIPS Phillips model
JCR_CHALK JCR chalk model
STRESS_1D 1D multi-yield model
HEAT_CONDUCTION Heat conduction
SCALAR_DIFFUSION Scalar diffusion
ELECTRIC_MODEL Electric models
Notes . .