Control lamellar direction in multi-scale RVE

[Hx83866]

Hello Everyone,

through the following script I'd like to generate multi-scale RVEs with three hierarchical structure, where the third layer consists of lamellar with specific width.

#!/bin/bash

NEPER="neper --rcfile none"

# RVE Side Length (μm): 60
# Element Size: 60x60x2
# Mesh Size (μm)：1.0

# Definition of Objective Function
PCK="diameq:lognormal(2.527,0.602),aspratio(1.9289,1,1)"
BNT="lamellar(w=2.0)"
T_num=20
iter=1

while [ $iter -le $T_num ]
do
    $NEPER -T -n 2::10::from_morpho\
          -morpho "voronoi::$PCK::$BNT"\
          -o "sheet_08_n$iter"\
          -reg 1\
          -sel 2.0\
          -id $iter\
          -oriformat 'geof'\
          -oridescriptor 'euler-bunge'\
          -statcell "scaleid(1),scaleid(2),scaleid(3)"\
          -domain "cube(60,60,2)" 

    $NEPER -V sheet_08_n$iter.tess -datacellcol 'scaleid(2)' -print sheet_08_n$iter
    $NEPER -M sheet_08_n$iter.tess -statelset 'diameq' -cl 1.0 -dim 3 -order 1 -elttype hex -format inp -o sheet_08_n$iter
    mkdir "sheet_08_n$iter"
    mv sheet_08_n$iter.* ./sheet_08_n$iter/
    iter=$(( iter+1 ))
done
#$NEPER -V sheet_08.tess,sheet_08.msh -dataelsetcol id -print sheet_08_mesh

exit 0

However, the lamellar structure in the generated RVEs still changed in Z-direction (varied in thickness). Reading the manual I knew that the direction of lamellar could be modified by given v=(dir_x, dir_y, dir_z). Nevertheless, the final RVE I got was not what I wanted, even no lamellar existed, despite given lamellar(w=2.0, v=( , , 0)).

So, does anyone could give me some suggestion?

thanks a lot!

Best Regards, Xiang HU

RWTH Aachen Universität