Appendix 5: Data file for remeshing

This file is used by the remeshing programs CRECOU, CREDAT, INTERP, and REMAIL (see \\Obelix\c\Progs\sources\Remeshing\). These programs were created in the 90's and have not been modified since. It is likely that some of them do not work properly anymore. Some info on these programs can be found in the “REMAIL1.DOC”, “REMAIL2.DOC”, “REMAIL3.DOC” files on the server (\\Obelix\c\Progs\sources\Manuels).
Note: By default, one can generate a remeshing with an empty .REM file, i.e. containing 8 blank lines. In this case, the default options are taken into account.

1st line (12I5)
LCASE⁵	Choice of the transfer method for σ and κ for the internal elements
	= 0 → interpolation between the new IP and the old IP by the1/RN method (default option)
	= 1 → interpolation between the new IP and the old IP by the generalized finite differences method
	= 2 → 2-step method a) in the old mesh, extrapolation of the IP at the nodes by modified interpolation functions and then average at the nodes which belong to the same material; b) interpolation at the IP of the new mesh by the classical interpolation functions of the element.
NSTAR¹⁰	If LCASE = 0 → exponent N of the 1/R^N interpolation (N=2 by default)
	If LCASE = 1 → number of points included in the radius of the central point to be interpolated = 16 or 20 = 100 if one wants all the nodes of the mesh
	If LCASE = 2 → no effect
METH¹⁵	If LCASE = 0 or 2 → no effect
METH¹⁵	If LCASE = 1 → choice of the method to select the points of the patch = 1 inscribed in a square = 0 inscribe in a circle
IANGLI²⁰	Parameter for the automatic computation of the lines (solid angles of intersections) in case of parabolic elements
	= 0 (default value)
	= 1 allows not to consider as line the elements sides topologically parabolic (see parameter ANGPRE for the definition of a topologically parabolic side)
NPT²⁵	Number of points for the smoothing of P and τ of the contact = 1 by default
IEMUL³⁰	Signification of the real parameter EMULT (cf. line 4)
	= 0 → EMULT is the multiplyer of the number of elements of the father mesh (default value)
	= 1 → EMULT is the absolute number of elements of the son mesh
IADP³⁵	Plastic admissibility for the elastoplastic law 72 (PLASOL)
	= 0 → nil (default value)
	= 1 → the stress states transgressing the plastic criterion are corrected by a return on the plasticity surface (recommended if LCASE=2 and/or if hardening). This return is made along the normal to the modified plastic surface (dilatancy equal to zero). In case of hardening, the hardening variable (equivalent plastic strain) is computed by the transfer of remeshing algorithm, and the derived variables are actualized from the hardening rule.
MODETQ⁴⁰	Regrouping method for the triangles used by GMAILL (see manual)
	= 0 → default value
	= 1 → other method
ITYPC⁴⁵	Type of meshing created by GMAILL (see manual)
	= 1 → frontal meshing (default value)
	= 2 → structured meshing
IMARK⁵⁰	Actualization of the passive markers stocked on the file *.MKP
	= 0 → nil (default value)
	= 1 → actualization of the markers
ICONIN⁵⁵	Particular treatment of the new interfaces
	= 0 → nil (default value)
	= 1 → the ε_eq are put to zero for the new CNTCP groups (useful if interfaces for which the elastoplastic law PLASOL is used with hardening are added).
INDMOD⁶⁰	Number of nodes whose coordinates have to be modified according to the user's choice
	= 0 → nil (default value)
	> 0 → Number of nodes whose coordinates have to be modified (to be defined later in the file)
2nd line (7I5, 2G10.0) - Parameters for the computation of crack (interposition of interface elements)
IFIX	For the computation of a crack which develops progressively, number of the imposed node where the crack segment will be located. If < 0 : the progression of the crack proportionally to the appearance of a new crack is privileged.
NODE1	Numbers of nodes which will be joined to define a crack
NODE2	Numbers of nodes which will be joined to define a crack
LMATE	Material number for potential cracks (to be defined a priori in the section COLAW of the file .DAT)
ILEAST	= 1 : the representative segment is computed by LEASTSQ
	= 2 : the representative segment is computed by MEDFIT
	If value < φ, the smoothing is applied
JCHOI	= 1 : only young points are admitted to define the crack
	= 2 : only points which belong to the shear band this time (“belonging” ≠ 0) are admitted
	= 3 : the points of any age (≠ 0) are admitted
	= 4 : the points of age ≤ IAGEM are admitted
	If value < φ, the normal to the previous segment (fixed origin) is computed and the i.p. of the “good” half-plan are taken into account
IAGEM	defined if SCHOI = 4 only (otherwise, they are not used)
VALFIX	value of the indicator for ICRIT = 8 value of COMP = percentage of the indic if ICRIT = 5
POIDS	ponderation factor of the point age (= 0∞ by default) \[WEIGHT=\frac{1}{Age^{POIDS}} \text{ or } \frac{1}{g^{POIDS}} \text{ if } POIDS<0\]
3rd line (4G10.0,4I5) - Density parameters
FMULT	High density coefficient of the mesh (low density = 1)
PINSD	influence parameters of the high density points
PINWD	influence parameters of the low density points
POIDS	ponderation factor of the points according to their age (= 0 by default) If < φ: weighting of the points according to their indic (= 0 by default)
IDIW	integer divisor of the diagonal width of the mesh L for the suppression of the weak points (W = Weak)
IDIS	idem for the strong points
ITDEP	= 0 mechanic forging case
	= 1 partially uncoupled thermomechanic forging case (to be reexamined 7/94, AMH)
	= 2 soil mechanics case
	= 3 mechanic case with introduction of shear bands (in the form of “cracks”)
	= 4 mechanic case with simplified cracks (J-D. Barnichon)
IGENER	= 0 only the density points which belong to the foundation will be taken into account
	= 1 only the density points situated in the solid and those situated on the foundation will be taken into account
	= 2 the high density points will be replaced by the minor squares line
	= 3 the density data of the father mesh will be taken into account (read on the file name_old.gml)

Remarks
IDIW is smaller than IDIS in all cases. The low density points situated at a square distance lower than $L^2/IDIW$ of the strong points and the strong points situated at a distance lower than $L^2/IDIS$ of the other strong points will be eliminated. In other words, if the user wants density points at least separated by a distance $D_{critical}$ that he chose, he will have to choose the corresponding factor IDI(W or S) such that: \[IDI = \left(\frac{L}{D}\right)^2\] Geometrically, if one considers each point successively as the center of a circle, one suppresses all the points (for example, the weak points) which are in the circle of radius R at this point with: \[R=\frac{L}{\sqrt{IDIW}}\] For more information details relating to the computation mode of the density points, see the description of the routine RECDEN given in appendix.
The POIDS will balance FMULT as a function of age according to the following rule: \[FMULT \rightarrow \frac{FMULT}{POIDS^{age-1}}\] It will be chosen:

> 1 if the density has to decrease when the point becomes older
< 1 if the density has to increase when the point becomes older

4th line (5G10.0)
EMULT	If IEMUL = 0 it is the multiplyer of the number of elements of the father mesh If IEMUL=1, it is the absolute number of the elements of the son mesh
ANGLIM	value of the angle α∈[0,π] under which there is a line (solid angle of intersection) Default value = 160°
ANGPRE	Useful if and only if ANGLI=1 for parabolic element Limit value of the angle β between the two segments going through the common node of two parabolic elements, for which a parabolic side is considered topologically linear Default value =5° A line will not be defined if β > ANGPRE
MAXELE	Maximum number of elements which may not be exceeded
RIGM	Penalty value of the new contact elements when a solid law is used
5th line (I5)
IFDENS	Number of foundation nodes to which one wants to give a density (constant during the whole simulation) Max 20
5th bis line (I5, 2G10,0) - repeated IFDENS times if IFDENS ≠ 0
NF	Number of the foundation node
DPL	Associated density coefficient
CPL	Influence reduction parameter
6th line (I5)
IFFAIBL	Number of density points generated between 2 different foundations nodes
6th bis line (2I5, 3G10.0) - repeated IFFAIBL times if IFFAIBL ≠ 0
N1	N° of the node belonging to the foundation 1
N2	N° of the node belonging to the foundation 2
RATIO	Ratio a/b ∈ [0,1]
DPL	Density coefficient of the generated points
CPL	Influence parameters of the generated points
7th line (2I5, 1G10.0) - Memorization parameters of the successive remeshings to limit the number of remeshings
The number of the step where the latest meshing has been realized is stocked in the file .LOA and the remeshing is effectuated only if the usual criteria are achieved and if the number of steps realized since the last remeshing > NDSTEPM
ILOA	≠ 0 if one wants to use a file .LOA
NDSTEPM	Minimum number of steps between two remeshings
DIV	Divisor of the time step (1.0 per default)
8th line (I5)
ILIC	=0 : no file .LIC. ETAGM represents then the absolute percentage of localized elements which provokes the stop of the solicitation for the remeshing.
ILIC	≠0 : one wants to print in the file .LIC the number of elements to be remeshed at each solicitation step. The remeshing will be effectuated only if the increment of the number of elements to be remeshed from step to step is sufficient, that is > ETAGM. In this case, ETAGM becomes the relative percentage of the localized elements since the latest converged step.
9th line (I5, 2G10.0) - repeated INDMOD times if and only if INDMOD>0 (cf 1st line)
INOD	Number of the node
XNEW	New abscissa for this node
YNEW	New ordinate for this node