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