User Tools

Site Tools


appendices:a5

This is an old revision of the document!


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 realize 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)
LCASE5Choice 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.
NSTAR10 If LCASE = 0 → exponent N of the 1/RN 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
METH15If 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
IANGLI20Parameter 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)
NPT25Number of points for the smoothing of P and τ of the contact
= 1 by default
IEMUL30Signification 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
IADP35Plastic 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.
MODETQ40 Regrouping method for the triangles used by GMAILL (see manual)
= 0 → default value
= 1 → other method
ITYPC45Type of meshing created by GMAILL (see manual)
= 1 → frontal meshing (default value)
= 2 → structured meshing
IMARK50Actualization of the passive markers stocked on the file *.MKP
= 0 → nil (default value)
= 1 → actualization of the markers
ICONIN55Particular 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).
INDMOD60 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)
IFIXFor 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
LMATEMaterial 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
IAGEMdefined if SCHOI = 4 only (otherwise, they are not used)
VALFIXvalue 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
FMULTHigh density coefficient of the mesh (low density = 1)
PINSDinfluence parameters of the high density points
PINWDinfluence parameters of the low density points
POIDSponderation factor of the points according to their age (= 0 by default)
If < φ: weighting of the points according to their indic (= 0 by default)
IDIWinteger divisor of the diagonal width of the mesh L for the suppression of the weak points (W = Weak)
IDISidem 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)
appendices/a5.1561454819.txt.gz · Last modified: 2020/08/25 15:33 (external edit)