Importing geometry problem (SolidWorks)

Hi all, I am totally new in liggghts...

Until now I worked on tutorials and made necessary changes on them to get worked it...

Right now I want to do something on my own...so I modeled simply box in SolidWorks export as ASCII .STL file and when I import geometry in to liggghts I am getting this message:

"Created orthogonal box = (-0.01 -0.06 -1.22) to (0.51 0.01 0.01)
1 by 1 by 1 MPI processor grid
ERROR on proc 0: Cannot open mesh file /meshes/ohisje.stl (../input_mesh_tri.cpp:105)
(../input_mesh_tri.cpp:105)"

I read something about problems of exporting .STL files from SolidWorks, is this true that mesh prom solid is not OK?

In attachment is my program and .stl file. I just want to fill box with some particles (using inputscript from example of in.meshGran with little changes).
And just to say I tried with "/" in front of meshes, example: /meshes/ohisje.stl....

Any answer will help me very much !!!

Best regards, Neven

And here is program:
# Wall import from CAD
atom_style granular
boundary m m m
newton off
atom_modify sort 0 0

communicate single vel yes
units si

region reg block -0.01 0.26 -0.01 0.1 -0.01 0.1 units box
create_box 1 reg

neighbor 0.02 bin
neigh_modify delay 0

#Material properties required for new pair styles

fix m1 all property/global youngsModulus peratomtype 5.e6
fix m2 all property/global poissonsRatio peratomtype 0.45
fix m3 all property/global coefficientRestitution peratomtypepair 1 0.7
fix m4 all property/global coefficientFriction peratomtypepair 1 0.05

#New pair style
pair_style gran model hooke tangential history
pair_coeff * *

timestep 0.00005

fix 1 all nve/sphere
fix 2 all gravity 9.81 vector 0.0 0.0 -1.0

#import triangular mesh
fix cad all mesh/surface file meshes/ohisje.stl type 1 scale 0.001

#use the imported mesh as granular wall
fix granwalls all wall/gran model hooke tangential history mesh n_meshes 1 meshes cad

region factory cylinder z 0.125 0.1 0.12 0.2 0.1 units box

#distributions for insertion
fix pts1 all particletemplate/sphere 15485863 atom_type 1 density constant 2500 radius constant 0.005
fix pdd1 all particledistribution/discrete 15485867 1 pts1 1.0

group nve_group region reg

#particle insertion
fix ins nve_group insert/stream seed 32452843 distributiontemplate pdd1 &
maxattempt 100 nparticles 10000 particlerate 1000 overlapcheck yes all_in no vel constant 0. 0. -1.0 &
insertion_face inface extrude_length 0.02

fix ts all check/timestep/gran 1000 0.1 0.1
compute 1 all erotate/sphere
thermo_style custom step atoms ke c_1 f_ts[1] f_ts[2] vol
thermo 1000
thermo_modify lost ignore norm no

#dump commands
dump dmp all custom 300 post/dump*.meshGran id type type x y z ix iy iz vx vy vz fx fy fz omegax omegay omegaz radius
dump dmpstl all mesh/stl 300 post/dump*.stl

#insert particles
run 10000 upto
unfix ins

#run
run 40000 upto

Attachment	Size
exercise1.tar_.gz	1.33 KB

Forums:

LIGGGHTS® - User Forum

mschramm | Wed, 04/05/2017 - 21:55

File Name .stl ~= .STL

I took a look and your input script references meshes/ohisje.stl
but the file in your meshes file is ohisje.STL...

There are still problems in the input script though... your next problem will occur in the fix ins nve_group insert_stream... This fix will be looking for a mesh surface with id name inface.

neven.marticnevic | Wed, 04/05/2017 - 22:28

THANKS !!! having new problem

Thanks on answer mschramm :)

I change input script after this problem, so I made geometry one more time and now simulation go throw reading stl, and I change what you mention...
With change input script I drawed new geometry in Rhino, end export as .stl and that is how I solve problem.

And now I get this error: (i google it and nothing comes up)
Reading STL file 'mesh/mesh.stl' (mesh processing step 1/3)

Reading STL file 'mesh/f.stl' (mesh processing step 1/3)
Fix particledistribution/discrete (id pdd1): distribution based on mass%:
pts1: d=1.000000e-02 (max. bounding sphere) mass%=100.000000%
Fix particledistribution/discrete (id pdd1): distribution based on number%:
pts1: d=1.000000e-02 (max. bounding sphere) number%=100.000000%
Setting up run at Wed Apr 5 22:16:44 2017

Import and parallelization of mesh cad containing 3251 triangle(s) successful
ERROR on proc 0: assertion failed (../region_neighbor_list_I.h:156)

What this erro mean?

neven.marticnevic | Wed, 04/05/2017 - 22:34

Here si script!!!

# Wall import from CAD
atom_style granular
boundary f f f
newton off
atom_modify sort 0 0

communicate single vel yes
units si

region reg block -0.4 0.4 -0.4 0.4 -0.4 0.4 units box
create_box 1 reg

neighbor 0.02 bin
neigh_modify delay 0

#Material properties required for new pair styles

#New pair style
pair_style gran model hooke tangential history
pair_coeff * *

timestep 0.00005

fix 1 all nve/sphere
fix 2 all gravity 9.81 vector 0.0 0.0 -1.0

#import triangular mesh
fix cad all mesh/surface file mesh/mesh.stl type 1 scale 0.001

#definition of insertion face
fix inface all mesh/surface file mesh/f.stl type 1

#distributions for insertion
fix pts1 all particletemplate/sphere 15485863 atom_type 1 density constant 2500 radius constant 0.005
fix pdd1 all particledistribution/discrete 15485867 1 pts1 1.0

#particle insertion
fix ins all insert/stream seed 32452843 distributiontemplate pdd1 &
maxattempt 100 nparticles 10000 particlerate 1000 overlapcheck yes all_in no vel constant 0. 0. -1.0 &
insertion_face inface extrude_length 0.02

fix ts all check/timestep/gran 1000 0.1 0.1
compute 1 all erotate/sphere
thermo_style custom step atoms ke c_1 f_ts[1] f_ts[2] vol
thermo 1000
thermo_modify lost ignore norm no

#dump commands
dump dmp all custom 300 post/dump*.meshGran id type type x y z ix iy iz vx vy vz fx fy fz omegax omegay omegaz radius
dump dmpstl all mesh/stl 300 post/dump*.stl

#insert particles
run 10000 upto
unfix ins

#run
run 40000 upto

mschramm | Fri, 04/07/2017 - 00:03

Nothing Jumps out

Could you try using the -echo both command and post what is printed to the screen

$LIGGGHTS -echo both < in.file > out.file

This may provide some insight on where the error is.