====== 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| |:::|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 |:::| |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}}\] {{ :appendices:a5_1.png |}} 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 -> \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° {{ :appendices:a_2.png?400 |}}| |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 {{ :appendices:a5_3.png?400 |}}| |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] {{ :appendices:a5_4.png?200 |}}| |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.| |:::|≠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|