Appendix 1: Printing options

The printings are saved on the file ‘NAME’.OUT. The information relating to the elements are also saved on the file ‘NAME’.STR
The following can be printed:

• IOPT(1) coordinates
• IOPT(2) velocities
• IOPT(3) accelerations
• IOPT(4) stays reaction
• IOPT(5) out-of-balance forces
• IOPT(6) stresses
• IOPT(7) state variables
• IOPT(8) displacements (only works if IMDIS=1, cf. line 2 col. 40 of data file)
• IOPT(9) criterion
• IOPT(10) out-of-balance forces, nodal velocities or actualized coordinates (and displacement) saving for drawing (cf. end of appendix)
• IOPT(11) printing interval
• IOPT(15) computation time at various stages of the time step (≥ 2 to print)
• IOPT(16) print internal energy in file NTWRK (=*.f81 or *.wrk) (≥ 2 to print). In case of a dynamic simulation, this printing is always done, along with the dynamic energy (=0 in case of a static simulation).

One can also print (very large!):

• The various variables in the elements routines IOPT(12)-1 = 1 to 3 (IMPER = IOPT(12)-1)
• The element variables before calling the element routine IOPT(13) > 1 and after the calling of this routine IOPT(13) > 2
• The coordinates and nodal velocities before the resolution IOPT(14) > 1 Remark: For i >10, no list required, regardless of the value of IOPT(i)
  

• IOPT < 2 or > 6 standard option (by default) : no printing
• IOPT = 2a printing at the end of each step for a reduced number (of DOF or elements) (definition of list)
• IOPT = 3a printing at each iteration for a reduced number (of DOF or elements) (definition of list)
• IOPT = 4a printing at the end of each step for all (DOF or elements)
• IOPT = 5a printing at each iteration for a reduced number (DOF or elements) (definition of list) and at the end of each step for all (DOF or elements)
• IOPT = 6a printing at each iteration and before the beginning of the execution for all (DOF or elements)

When IOPT = 2, 3 or 5, one must give the list of the nodes or elements concerned by the printing. If several lists are required, they must be given in the order indicated by IOPT.
Example (IOPT(I),I=1,15) = 2 2 0 4 0 2 0 0 0 0 1 0 0 0 0
List of nodes where the coordinates are printed 10 -21 25 30 0
List of nodes where the velocities are printed 24 10 -20 0
List of nodes where the stresses are printed 3 6 -10 0
Remark The list 3 6 -10 0 means : element 3 + elements 6 to 10
The list 10 -21 25 30 0 means : nodes 10 to 21 + node 25 + node 30
Each new list begins on a new line. 14 numbers can be read on a line (the 0 indicating the end of a list included). If the list contains more than 14 numbers, it continues on the following line.

All the steps whose number is a multiple of IOPT(11) are printed. By default, IOPT(11) = 1.

One can obtain printings of a lot of variables using the parameters IOPT(I) = 12 to 14.
These printings are large and happen at each iteration of each step. They must be used temperately.

If IOPT(10) > 1, Lagamine will save at each iteration of each step (OFTEN !) the out-of-balance forces, the nodal velocities or the actualized coordinates and displacement in the file 14 (*.RES), readable by DESFIN.
How to select which quantity will be printed?
If IOPT(10) = 4 or 6, no list is required. In that case, the out-of-balance forces are saved;
If IOPT(10) = 2, 3 or 5, a list is required:
If the first element of the list = 1, the out-of-balance forces are saved,
If the first element of the list = 2, the nodal velocities are saved,
If the first element of the list = 3, the actualized coordinates (and the displacements in order to be able to scale them by a factor) are saved.
Remark: The choice between the out-of-balance forces or the velocities is NOT based on the IDENT anymore!!!
In order to view the result in Desfin, you have to choose in the general menu :
Drawing type n°7 (“Flux au Noeud”) to observe the out-of-balances forces or the velocities
Drawing type n°2 (“Déformée”) to observe the deformed configuration.
In order to view the content of the RES file (which time step and iterations are saved?), you must first plot any time step. Once the drawing is on the screen, go to the menu “Options” → “Multiplicateurs de charge…”.

IOPT(9) for the printings relating to the criterion

• if IOPT(9) <2 or >6 no printing
• if IOPT(9) >=2 or ⇐6 printing of : - The Darcy nodal speeds
• The Lagrangian concentrations
• The zoomed elements

at the end of each step.