neperfepx / neper

Polycrystal generation and meshing
http://neper.info
GNU General Public License v3.0
199 stars 53 forks source link

Create an elset for face Z1 #334

Closed mzavard closed 2 years ago

mzavard commented 2 years ago

By defalt in Neper, when you get an .inp you have node sets for all the faces of your cube. I would like to create a set of element for the face Z1 (coo.z=1) of my cube.

After reading neper's documentation i choose to call my elset Y1ELMT and i tried to give her this definition : z>=0.7

Here are my commands lines:

neper -T \ -n from_morpho \ -morpho "diameq:lognormal(0.7,0.05),1-sphericity:lognormal(0.145,0.03)" \ -reg 1 \ -group "id<=25?1:2" \ -domain "cube(1,1,1)" \ -o simulation \ -statpoly id,vol\ -statseed id,z

Generate a finite element mesh

neper -M simulation.tess \ -order 2 \ -elttype tet \ -statelset id,coo\ -statelt3d id,coo\ -part 4 \ -rcl 0.35\ -elset "Y1ELMT:z>=0.7"\ -format vtk,msh,phase,inp exit 0

I got this as an error: Info : [o] Writing file `simulation.msh'... ./generate_mesh.sh: line 27: 2024 Segmentation fault (core dumped) neper -M simulation.tess -order 2 -elttype tet -statelset id,coo -statelt3d id,coo -part 4 -rcl 0.35 -elset "Y1ELMT:z>=0.7" -format vtk,msh,phase,inp

I'm sure that it's a syntax problem when i defined my element set but i can't find documentation about how to do that...

mzavard commented 2 years ago

I could just reuse the nset of z1 and define them as an elset...