|
ALGORITHM = ClosestPoint ElementSizeSearchRadius ConsiderOndulation RVE_Detection | See descriptions below. |
ClosestPoint is the most commonly-used option and accounts for most of the mapping scenarios targeted by the Envyo application. Beginning with the target mesh, the next element or integration point is identified and its data transferred to the target mesh.
The ElementSizeSearchRadius option supports mapping of results from mesoscale process simulations using shell elements onto a target mesh. The target mesh element size is therefore used as a search radius, whereas the source mesh elements are projected onto the target mesh along the source-mesh normal and a check is performed to see if the intersection is within the target element. The element size can be scaled using Scale_SearchRadius (below).
The ConsiderOndulation option considers the through-thickness position of the identified fiber orientation. This information is therefore not considered as provided by user input, but by the position in the source mesh related to the target mesh. This option does not account for the ordering of the input parts - it first maps those that are closest to the target part, then the second closest and so on. No resinuous areas can be assigned with this option. If the MapThickness option (below) is enabled, the offset between the source layers is also considered, which usually leads to varying thicknesses throughout the part.
With RVE Detection, you can define specific criteria which allows the Envyo application to decide if a specific combination of fibers has been identified and therefore assigns a specific, user-defined material ID to the dedicated integration point.
|
MapStress = YES/NO CP AVG SHEPARD NO1 SHEPARD NO2 SHEPARD NO3 SHEPARD NO4 FE FE NODAL AVG | YES/NO activates/deactivates the mapping of stress data (every value stored in the *INITIAL STRESS OPTION fields on the source side). The default algorithm used is the closest point (CP) algorithm (or nearest neighbour search). For explanations of the other options, see below. |
CP - directly activates the closest point algorithm to transfer everything stored in the *INITIAL STRESS OPTION fields on the source side. This is the default algorithm.
AVG - the euclidean average is taken from the values found within a specific search radius. By default, the search radius is based on the average target mesh element length.
SHEPARD NO1 - activates the general Shepard's inverse distance weighting function. You may specify the exponent using the SHEPARD EXPONENT option.
SHEPARD NO2 - activates the Shepard's inverse distance weighting function, with a minimum of four and a maximum of ten points considered for the interpolation function. You may also specify the exponent using the SHEPARD EXPONENT option.
SHEPARD NO3 - activates the Shepard's inverse distance weighting function, with a minimum of four and a maximum of ten points considered for the interpolation function, including direction. You may also specify the exponent using the SHEPARD EXPONENT option.
SHEPARD NO4 - activates the Shepard's inverse distance weighting function, with a minimum of four and a maximum of ten points considered for the interpolation function, including direction and the slope. You may also specify the exponent using the SHEPARD EXPONENT option.
The search radius for the four options listed above is adjusted based on Shepard's proposal, see [22].
FE - activates the interpolation of values found in the *INITIAL STRESS OPTION fields based on finite element shape functions and point projection of the considered integration points into the plane of the source shells. If the source and target through-thickness integration rules are different, through-thickness interpolation is performed based on the method proposed in [24], minimizing the difference in the areas underneath the through-thickness integration rules on source- and target side for all values being transferred. This interpolation result can be visualized using the WriteTTInterpolationFile and TTInterpolationFile options.
FE NODAL AVG - activates the interpolation of values found in the *INITIAL STRESS OPTION fields based on finite element shape functions and point projection of the considered integration points into the plane of the source shells. In addition to the described method above, values of the connected elements on the source side are also considered through extrapolation of the integration point values onto the element nodes, taking the euclidean average and performing a value recovery on the target integration point.
|
MapStrain = YES NO | Specify whether or not strains should be transferred. |
|
MapThickness = YES/NO CP AVG SHEPARD NO1 SHEPARD NO2 SHEPARD NO3 SHEPARD NO4 FE FE NODAL AVG |
Activate or deactivate thickness transfer. This option requires *ELEMENT SHELL THICKNESS fields in the *.dynain file. If TargetThickness is also defined, it will be ignored. The default option is CP (closest point). Thickness mapping is based on the thicknesses stored on the nodes. Thickness is averaged over all elements attached to a node, regardless of the chosen interpolation option. The listed interpolation options work in a similar way as those described under MapStress. |
|
Thck_Avg_Opt = Ele Avg Nodal Avg |
Thickness averaging option. When using *ELEMENT SHELL THICKNESS, nodes may hold different thickness information, depending on the element to which they belong. This is especially criticial for components with ribs. With the default setting (Ele Avg), the Envyo application calculates through-thickness integration points using the nodal thicknesses stored on the element nodes by interpolating the resulting thicknesses at the in-plane integration points via shape functions, allowing for the correct calculation of the local through-thickness integration points. By choosing the Nodal Avg option, the application first collects all thicknesses stored on one node from its attached elements and averages these thicknesses. Interpolation to the in-plane integration points using shape functions is then performed. |
|
TENSORIALINT = AVG INV-N INV-R INV-K INV-U INV-C | This option controls how tensorial values are interpolated. Only stress tensors are considered with the proposed methods. As these methods are still under development, you should use the default AVG method. |
Besides using the standard euclidean average (AVG), which is the default and recommended value, further interpolation methods based on tensor invariants have been investigated due to the fact that swelling effects are often seen when performing standard euclidean average of tensorial data (the tensor shape changes).
Based on the work of [9], several invariant sets have been implemented for the interpolation functions available in MapStress, and linear interpolation is performed for these invariant sets instead of the single tensorial values. These interpolated invariants allow you to recover the eigenvalues based on a closed form solution of the cubic function and together with the eigenvectors determined through linear interpolation of the originating tensors, the resulting tensor can be recovered.
Due to the boundaries of the recovery function, which is based on the arccosine defined between -1.0 ≤ x ≤ 1.0, this might lead to unexpected behavior with the current implementation.
INV-K refers to an orthogonal invariant set based the tensor trace, the Frobenius norm of the tensor's deviatoric (fractional anisotropy of the tensor), and the mode of the deviatoric.
INV-R refers to an orthogonal invariant set based the tensor norm, the relative anisotropy of the tensor, and the mode of the deviatoric as in INV-K.
INV-N is a non-orthogonal invariant set based on the tensor trace, the fractional anisotropy, and the tensor mode.
INV-U is a uniform invariant set based on the tensor trace, the scaled relative anisotropy of the tensor, and the scaled angular mode AM.
These invariant sets have been proposed in [9] for positive definite tensors gained in DT-MRI measurements. Generally, stress tensors are not positive definite. The Cauchy stress tensor invariant set INV-C is therefore added to the interpolation options, even though this might not solve the issue regarding the arccos in the recovery function.
| TargetThickness = DOUBLE | Define the thickness of the target shell mesh. |
| SHEPARD EXPONENT = DOUBLE |
Define the exponent of the Shepard's inverse distance weighting function. Default = 2.0. |
For cases where the target mesh does not contain any information regarding the integration rules being used, such as element formulation and through-thickness integration points, you must define these manually. The input for this is as follows:
| NPLANE = INT |
1 - reduced integrated thick shell elements 4 - fully integrated thick shell elements This option was formerly known as NumberOfTARInPlaneIPs. |
| NTHICK = INT |
Define the number of through-thickness integration points (IPs). If the Algorithm=RVE Detection option is used, ensure that the number of through-thickness IPs refers to the number of fibers being mapped, the number of source shell element stacks, and the number of through-thickness IPs in the source mesh. The number of through thickness IPs can be reduced using the ThroughThicknessAveraging option (see Further Options ). See also Figure 5.3: Number of stacks, integration points and fiber IDs for *MAT 249 and Note 1 under Notes. This option was formerly known as NumberOfTARThroughThicknessIPs. |
|
IntegrationRule = Gauss Lobatto Autoform Moldflow | Define the through-thickness integration rule for the mapping result. This option directly affects the positions of the through-thickness integration points on the target mesh. |
|
BIAX = YES NO |
If fiber orientation is derived from node numbering, usually only the direction from node #1 to node #2 is used. The BIAX flag activates the additional usage of node #1 to node #4 in order to consider two-directional materials (for example, woven materials). The directions are written to two integration points in the same order as described above. |
| INN = INT | This flag is similar to the invariant node numbering flag in LS-DYNA provided in the /CONTROL ACCURACY field, see [19]. To properly calculate orientations with respect to the element coordinate system, the Envyo application requires information about how LS-DYNA calculates the element coordinate system. The default, as with the LS-DYNA application, is off (0). |
|
Search_Radius = SrcEleLen TarEleLen DOUBLE | The search radius declaration for the mapping algorithm. SrcEleLen is the default, which uses average source mesh element size as the search radius. Average target mesh element size can be used by defining TarEleLen, or a positive DOUBLE value can be defined. |
| Scale_SearchRadius = DOUBLE | Coefficient to scale search radius. For the ElementSizeSearchRadius option, the search radius for element projection is scaled with this factor. The default value is 1.0, meaning that the original element size is used as a search radius. |
| Shell_Option = Composite Composite Long |
This option tells the Envyo application to generate *ELEMENT SHELL COMPOSITE( LONG) elements from the mapped orientations. |
| SORT = BUCKET | Using bucket sort is strongly recommended, as it provides a substantial performance improvement for the search algorithm. |
| REPEAT = YES | Enable this option to ensure that all elements and integration points receive mapped data. When there is a significant difference in element sizes between the source and target meshes, the default bucket refinement may be insufficient to cover all points, sometimes by design. In such cases, this flag must be set to guarantee complete data coverage. |