Sidebar
This is an old revision of the document!
Table of Contents
Preprocessor
Preprocessing data for Lagamine program.
This page contains a description of the data to be indicated in the *.lag file for the Prepro.
1. Title
| Format (9A8,I3) | |
|---|---|
| Title(9) | Title, in 72 characters |
| I99999 | Big mesh index: = 1 if NUMNP (total number of nodal points)>99999 = 0 else |
2. Control data of the program
The control data is defined in 4 lines (in exponent, column on which the variable must be fit in with).
| 1st line (14I5) or (2I5, I10, I5, I10, 10I5) if I99999=1 | |
|---|---|
| NTANA5 | Type of analysis (see (1)) |
| IANA10 | Type of problem: |
| =1 - Plane stress state ($EP\sigma$) | |
| =2 - Plane strain state ($EP\varepsilon$) | |
| =3 - Axisymmetric state ($axi$) | |
| =4 - 3D state ($3D$) | |
| =5 - Generalized plane state ($EPG$) see (2)) | |
| NUMNP15 or 20 | Total number of nodal points |
| NELTP20 or 25 | Number of ELement Types used (or groups) |
| MELEM25 or 35 | Total number of elements (all types together) |
| NDISP30 or 40 | Number of degrees of freedom where non zero DISPlacements are imposed |
| NFOUN35 or 45 | Number of foundations |
| NSEGT40 or 50 | Total number of segments defining the NFOUN foundations |
| IPILO45 or 55 | Number of pilot nodes for the foundations |
| NLAW50 or 60 | Number of constitutive LAWs used |
| IRICE55 or 65 | Indicator for the bifurcation calculation (RICE's criterion) – No more used. |
| ICGEO60 or 70 | Indicator for geometrical distortions (%age of distorted elements) |
| IERRO65 or 75 | Indicator for a posteriori error calculations |
| ICRIT70 or 80 | = 0 nil |
| = 1 forging simulation, criteria calculation to detect the too-distorted elements and automatic remeshing if the number of distorted elements exceeds the maximum fixed in the execution data and existence of “EI” (M. Dyduch) | |
| = 2 simulation with shearing bands and automatic remeshing if the number of bifurcated elements exceeds the maximum fixed in the execution data | |
| = 5 simulation where the VILOTTE's criterion is used as indicator | |
| = 6 forging simulation, different from Dyduch | |
| = 7 proportional fatigue computation (see fatigue calculation in Hill_3D_KI law), you must generate the file 61 (.ntfdam) (0) | |
| 2nd line (15I5) | |
| IRENU5 | Index for the DOF renumbering = 0 nil = 1 reading of the RENUM card and renumbering of the DOF |
| IGRAV10 | Index for the analysis with gravity (soil mechanics) = 0 nil = 1 reading of the GRAVI card |
| IDMAT15 | = 0 or 2 reading of the FIXED card = 1 reading of the IDMAT card |
| MPARA20 | Maximum number of real PARAmeters for any constitutive law. (3) Default value : 40 |
| MVARI25 | Maximum number of internal VARiables per integration point. (3) Default value : 40 |
| MSIGP30 | Maximum number of stresses per point. (3) Default value : 18 |
| MPARI35 | Maximum number of integer parameters for any constitutive law. (3) Default value : 20 (Be careful, often not sufficient for thermo-mechanical-metallurgical coupling) |
| IMDIS40 | = 1 displacements can be printed (the CONEC dimension increases) (compulsory for BEM2D) = 0 displacements are not printed NCOL = 1 for static analysis, NCOL = 3 for dynamic analysis. = 2 NCOL = 2 for static analysis useful for OCASFO routine (see data appendix 11) ≥10 presence of a dead load, the NCOL value is incremented by 1. (see Wang) |
| IPIFO45 | = 0 nil = 1 force pilotage by the routine OCASFO (see appendix) |
| ICEVO50 | = 0 nil = 1 checking of the total volume of the piece computed by the Jacobian matrix determinant (ready for BLZ3D) |
| ITREE55 | = 0 nil = 1 search of the 3D contact by arborescence; ITREE = number of foundations using the sub-tree method (4) |
| ICOUP60 | = 0 nil = 1 the bandwidth of the thermomechanical calculation takes into account the splitting between thermics and mechanics |
| ILIEN65 | = 0 nil = 1 the number of the solid element adjacent to each contact element or foundation segment is memorized. |
| NDUPL70 | = 0 nil, no section 7 = n reading of the DUPLI card, n is the number of pair of Duplicata nodes. |
| INCFO75 | = 0 nil, no section INCFO ≥ 1 Read Incremental Forming parameters (SEMPER project – C. Henrard, C. Bouffioux). If INCFO > 1, debug messages will be printed (see module MINCFO → DEBUG_MST=INCFO-1) |
| 3rd line - dynamic analysis (4I5) | |
| This line is left blank if the analysis is not dynamic | |
| NMASC5 | Number of DOF where non zero masses are imposed |
| NDAMC10 | Number of DOF where non zero dampings are imposed |
| NVEL015 | Number of DOF where non zero initial speeds are imposed |
| NACC020 | Number of DOF where non zero initial accelerations are imposed |
| 4th line - metallurgical analysis (4I5) | |
| This line is left blank if the analysis is not metallurgic | |
| NLAWM5 | Number of metallurgic laws (can be equal to 0) Rem. If NLAWM≠0, a x.MET file is required for this type of metallurgical data |
| MTPAM10 | Maximum number of temperatures per table (50 per default) |
| MCPAM15 | Total number of columns for all the metallurgic tables |
| if META = 0 or for META = 4 and IMETH = - 5: 4 (3 Euler's angles + 1 weight) Default value: 25 if NTANA = 3 or 6 85 if NTANA = 4 or 9 |
|
| If META ≠ 0, MCPAM is nil | |
| MDPAM20 | Number of parameters not in the table |
| If META = 0: MDPAM = 108 if NTANA = 3 or 6 MDPAM = 119 if NTANA = 4 or 9 |
|
| If META = 1: MDPAM = 88452 (Leuven, P. Van Houtte's law) imposed in the code | |
| If META = 3 or if (META = 4 and IMETH (cf. law MIPAY3) = ± 2): MDPAM computed according to Leuven files, no default value (214 for a six order series) | |
| If META = 4 and IMETH (cf. law MIPAY3) = 3 or 4: MDPAM = maximum number for crystals necessary to define the texture x 3, no default value |
|
| If META = 4 and IMETH = -5 (Taylor) MDPAM = sliding numbers * 8+1 for the used crystals |
|
| Else: MDPAM imposed to 1 in the code | |
| META25 | =0 Nancy approach (AM Habraken 1989) Phase transformations modelisation |
| =1 Texture approach, Paul Van Houtte's surface (law ANI3VH) | |
| =2 Hill's surface (law ANI3VH) | |
| =3 Surface defined by a 6 order series in the stresses space (Jan Winters' approach, 1996) (law ANI3VH) | |
| =4 Surface defined by pieces and based on texture (Munhoven-Radu-Habraken approach, 1997) (law MIPAY3) Depending on the IMETH variable defined in the law MIPAY3, one can have: - Hill's surface discretized by facets; (IMETH= +1) - Hill's surface non discretized by facets; (IMETH= -1) - Surface defined by a 6 order series in the strains space (Bert van Bael's approach), discretized by facets; (IMETH = +2) - Surface defined by a 6 order series in the stresses space (Jan Winters' approach), non discretized by facets; (IMETH = -2) Equivalent version to META = 3 - Surface discretized by facets and based on the texture (Taylor Ulg A); (IMETH= 3) - Surface discretized by facets and based on the texture with actualisation of this surface (Taylor Ulg B); (IMETH= 4) - No plasticity surface but Taylor ULg A is the mean for the $\sigma_{micro}$ and $\sigma_{macro}$ computation (IMETH = - 5) |
|
(0) NTFDAM file structure
→ for mono-block loading:
| 1st line (I5) | |
|---|---|
| NL | =1 (NL is the number of blocks) |
| 2nd line (I10, G7.5) | |
| ncycle | =0 |
| PERIODE | =input value |
→ for multi-block loading:
| 1st line (I5) | |
|---|---|
| NL | = number of blocks |
| 2nd line (I10, G7.5) - repeated NL times | |
| ncycle | input value |
| PERIODE | input value |
(1) Remark on NTANA
In general, a node has 3 DOF displacements and 1 DOF temperature. In the case of 3D shell elements, it has 3 DOF displacements and 3 DOF rotations. The variable NTANA allows the program to decide which is the maximum number of nodal DOF potentially useful for a given type of analysis (dimension NDOFN, classed as NSPAC, space dimension). It also allows to define the active DOF (i.e. those which can produce an equation) and the passive ones (which do not produce equations).
A negative value of NTANA is used for dynamic analysis.
| NTANA | (2 translations + water pressure + air pressure + temperature) |
| = 1 - Plane mechanical analysis NDOFN = NSPAC = 3 : X,Y,T (no equation for DOF 3) |
|
| = 2 - 3D mechanical analysis NDOFN = NSPAC = 4 : X, Y, Z, T (no equation for DOF 4) |
|
| = 3 - Plane thermal analysis NDOFN = NSPAC = 3 : X, Y, T (no equation for DOF 1 and 2) |
|
| = 4 - Plane thermomech. analysis NDOFN = NSPAC = 3 : X, Y, T |
|
| = 5 - 2D coupled mechanical water – thermal – air flow analysis NDOFN = NSPAC = 5 : X, Y, PW, Pa, T |
|
| = 6 - 3D thermal. analysis NDOFN = NSPAC = 4 : X,Y,Z,T |
|
| = 7 - Plane mechanical analysis with shells NDOFN = NSPAC = 3 (2 translations, 1 rotation) X,Y,ZZ |
|
| = 8 - 3D mechanical analysis with shells NDOFN = NSPAC = 6 : X, Y, Z, ψ1, ψ2, ψ3 (3 translations + 3 coordinates of rotation vector) |
|
| = 9 - 3D thermomechanical. analysis NDOFN = NSPAC = 4 : X, Y, Z, T |
|
| = 10 - Plane mechanical analysis, second gradient method (Ouafa ElHammoumi) NDOFN = NSPAC = 5 : X, Y, λ, $\frac{d\lambda}{dX}$,$\frac{d\lambda}{dY}$ |
|
| = 13 - 3D coupled mechanical water – thermal – air flow analysis NDOFN = NSPAC = 6 : X, Y, Z, PW, Pa, T |
|
| = 14 - Plane mechanical analysis, second gradient method (Grenoble) NDOFN = NSPAC = 6 : X, Y, V11, V12, V21, V22 (ou λ11, λ12, λ21, λ22 for central node) |
|
| = 15 - 2D coupled mechanical water – thermal – air flow – chemical analysis NDOFN = NSPAC = 6 : X, Y, PW, Pa, T, c |
|
| = 16 - Plane coupled mechanical water flow analysis, second gradient method (Grenoble) NDOFN = NSPAC = 7 : X, Y, PW, V11, V12, V21, V22 (ou λ11, λ12, λ21, λ22 for central node) |
|
| = 17 - Plane coupled mechanical water – thermal – air flow analysis, second gradient method (Grenoble) NDOFN = NSPAC = 9 : X, Y, PW, Pa, T, V11, V12, V21, V22 (ou λ11, λ12, λ21, λ22 for central node) |
|
| = 18 : 2D coupled mechanical water – thermal – air flow – chemical analysis NDOFN = NSPAC = 8 : X, Y, PW, Pa, T, c1, c2, c3 |
(2) Remark on IANA=5 (generalized plane state)
- The last node is a special one; it always has 3 DOF, $\alpha_0$, $\alpha_1$ et $\alpha_2$, such that the thickness (in z direction) of a slice of slid varies according to: $b = \alpha_0 + \alpha_1 X + \alpha_2 Y$
- To be consistent with these 3 DOF of the last node, it is necessary to specify NTANA such that there exists at least 3 active DOF per node but may perhaps fix the non-significant ones.
For example, for the mechanical generalized plane state, NTANA = 7, IANA = 5 and fix the 3rd DOF for all the nodes (except for the last one eventually). - Initial values of $\alpha_0$, $\alpha_1$ et $\alpha_2$ can be introduced: they are the “coordinates” of the last node.
2 remarkable particular cases are worth noting :- Plane strain state : $\alpha_0$ = 1, $\alpha_1$ = $\alpha_2$ = 0 (and $\alpha_i$ all fixed)
- Axisymmetric state : $\alpha_0$ = 0, $\alpha_1$ = 1, $\alpha_2$ = 0 (and $\alpha_i$ all fixed).
- The last node appears in all the elements; therefore, the bandwidth is maximum → use the skyline system solver.
(3)
These variables are used in dynamic allocation of memory. In case of thermal analysis, the properties can be described at different temperatures, the default allocations could be insufficient and one must calculate MPARA.
(4)
Sub-tree algorithm for contact searching, developed by OCAS, has been implemented both for solid (CFI3D) and shell element (COQJ4). It is activated by option ITREE≥1 in PREPRO. Details of algorithm can be found in “OCAS Technical Report Hgr/9501, 1995” by H. Bruneel and I. De Rycke.
Currently sub-tree contact searching algorithm only works with tools which are driven by a displacement pilot node. No rotation of tools is allowed. Even if the foundation is fixed, a (fixed) pilot node must be defined. Also the algorithm works with tools described by triangular segments (ICODE=5). (Note by Y. Kim, Jan. 97). Note that the sub-tree method only works with foundations containing more than one triangular segment.
3. Nodes coordinates
| 3.1 Sub-title (A5) | |
|---|---|
| TITLE | “NODES” written in the first 5 columns |
| 3.2 Nodes definition (I5, 8G10.0) or (I10, 8G10.0) if I99999=1 | |
| INODE | Node number |
| XYZ(INODE,J)J=1:NSPAC | Node coordinates |
Remark : convention
The plane or axisymmetric state axis system is direct (right-handed). Conventionally the Y axis is directed toward the top. In the axisymmetric state, the Y axis is the same as the symmetry axis. All the nodes, elements, and local axis lists definitions must be given in the same axis system, implicit in the whole of this handbook.
4. Gravity definition
This section only exists if IGRAV=1.
| 4.1 Sub-title (A5) | |
|---|---|
| TITLE | “GRAVI” written in the first 5 columns |
| 4.2 Gravity definition (3G10.0) | |
| GRAVIX | Acceleration value in the Ox direction |
| GRAVIY | Acceleration value in the Oy direction |
| GRAVIZ | Acceleration value in the Oz direction |
5. Degrees of freedom (DOF) renumbering
This section only exists if IRENU=1.
The renumbering does not change the nodes numbers. It only changes the equations number corresponding to the DOF of these nodes.
File DFIXED.F in Prepro
| 5.1 Sub-title (A5) | |
|---|---|
| TITLE | “RENUM” written in the first 5 columns |
| 5.2 (5I5, G10.0) | |
| ISTAR | For ITYREN=0 or 1: Initial node for renumbering |
| For ITYREN=2: Hundreds: Direction in which the structure has the greatest number of nodes 1 = X, 2 = Y, 3 = Z Tens : Intermediate direction Unity : Direction in which the structure has the smallest number of nodes Possible values: 2-D : 120 or 210 3-D : 123, 132, 231, 213, 312 or 321 |
|
| For ITYREN = 3: not used | |
| KZONE | Default value = 4 = NOT used if ITYREN=2 or 3 |
| IPRIN | Printing level index |
| = 0: bandwidth before and after renum | |
| = 1: idem 0 + new numbering + new ID matrix | |
| = 2: idem 1 + connections matrix + summary of the obtained bandwidths | |
| = 3: idem 2 + list of obtained numberings | |
| = NOT used if ITYREN=2 or 3 | |
| IFOND | = 0: no foundation or if fixed foundation |
| = 1: foundation whose DOF are free | |
| = NOT used if ITYREN=2 or 3 | |
| ITYREN | = 0: Maximum bandwidth optimisation (oil stain) |
| =1 : Mean bandwidth optimisation (oil stain) | |
| =2 : Directional optimisation | |
| =3 : Sloan method based on graph theory (Sloan, IJNME, vol 28, pp 2651-2679, 1989) | |
| BDFAC | Default value = 1.05 NOT used if ITYREN=2 or 3 |
KZONE and BDFAC are factors which have an effect on the new numbering search as the computer time.
6. Duplicata nodes
This section only exists if NDUPL>0.
| 6.1 Sub-title (A5) | |
|---|---|
| TITLE | “DUPLI” written in the first 5 columns |
| 6.2 (2I5) or (2I10) if I99999>1 | |
| IDUPL(I,1), IDUPL(I,2) | I=1, NDUPL Nodes list with IDUPL(I,1)<IDUPL(I,2), and IDUPL(I,2) ≠ IDUPL(II,1) and IDUPL(I,2) ≠ IDUPL(II,2), II=1, I-1 |
7. DOF fixations
This section only exists if IDMAT = 0 or 2. One must define in this section all the stays without structure displacements.
| 7.1 Sub-title (A5) | |
|---|---|
| TITLE | “FIXED” written in the first 5 columns |
| 7.2 (14I5) or (I5,9I10) | |
| In the remeshing case, the nodes have to be given in the trigonometric order and without the sign “-” NB: if IDMAT=0 (default case), the format (14I5) is used if IDMAT=2, the format (I5,9I10) is used. This is useful if NNODE≥10.000 to be able to write a negative number… |
|
| ICOMP | Fixed DOF number |
| NODEF(13) or (9) | List of nodes whose DOF is fixed (max 13 x I5 per line if IDMAT=0: max 9 x I10 per line if IDMAT=2 ) If NODEF(I) < 0, then all the nodes between NODEF(I-1) and NODEF(I) have the DOF ICOMP fixed. |
| 7.3 End of the section – A blank card | |
For the thermal and flow problems, only the DOF nr 3 (2D, axisym) or 4 (3D) can be fixed.
8. ID Matrix
This section only exists if IDMAT = 1.
| 8.1 Sub-title (A5) | |
|---|---|
| TITLE | “IDMAT” written in the first 5 columns |
| 8.2 (14I5) or (I10,13I5) if I99999=1 | |
| INODE(ID,(INODE, IDOFN)) IDOFN=1,NDOFN | ID matrix line by line with the nodes numbers in the increasing order |
9. Setting of the DOF to a non-zero value
This section only exists if NDISP ≠ 0. After a sub-title card, the program reads the cards 9.2 until it founds NDISP imposed displacements to non zero values.
| 9.1 Sub-title (A5) | |
|---|---|
| TITLE | “DISPL” written in the first 5 columns |
| 9.2 Imposed displacements (I5, 6G10.0) or (I10, 6G10.0) if I99999=1 | |
| INODE | Node number |
| DISPL(I) I=1,NDOFN | Imposed displacements at these node's DOF. An imposed to zero displacement is ignored. |
9-bis. Incremental forming
This section only exists if INCFO = 1.
… To be completed …
10. Setting of the initial speeds in dynamics
This section only exists if NVELO ≠ 0 and NTANA < 0.
After a sub-title card, the program reads the cards 10.2. until it founds NVELO initial speeds different from zero.
File VELACC.F in Prepro
| 10.1 Sub-title (A5) | |
|---|---|
| TITLE | “VELO0” written in the first 5 columns |
| 10.2 Speed (I5, 6G10.0) or (I10, 6G10.0) if I99999=1 | |
| INODE | Node number |
| CONEC(INODE,I,3) I=1, NSPAC | Initial speeds at these node's DOF (a speed equal to zero is ignored) |
11. Setting of the initial accelerations in dynamics
This section only exists if NACC0 ≠ 0 and NTANA < 0.
After a sub-title card, the program reads the cards 11.2 until it founds NACC0 initial accelerations different from zero.
| 11.1 Sub-title (A5) | |
|---|---|
| TITLE | “ACCE0” written in the first 5 columns |
| 11.2 Acceleration (I5, 6G10.0) or (I10, 6G10.0) if I99999=1 | |
| INODE | Node number |
| CONEC(INODE,I,3) I=1, NSPAC | Initial speeds at these node's DOF (a speed equal to zero is ignored) |
12. Imposed nodal forces
This section always appears
| 12.1 Sub-title (A5) | |
|---|---|
| TITLE | “FORCE” written in the first 5 columns |
| 12.2 Imposed forces (I5, 7G10.0) or (I10, 7G10.0) if I99999=1 | |
| INODE | Node number |
| FOMAN(I) I=1, NDOFN | Imposed forces at the DOF of this node |
| IMORT | 0 = Ordinary forces |
| 1 = “Dead” forces Rem: The dead forces are forces which do not increase when the spherical step method is used. |
|
| 12.3 End of section – A blank card | |
Remark : In thermal and flow problems, the forces are heat flux or flow, positive when they go in the system. On the other hand, the reactions are the heat flux or flow, positive when they go out the system.
