Re: [Yade-users] [Question #707675]: Shear force is lower than MaxShearForce but sliding occurs

2023-09-17 Thread Leonard
Question #707675 on Yade changed:
https://answers.launchpad.net/yade/+question/707675

Status: Answered => Solved

Leonard confirmed that the question is solved:
Thanks Karol Brzezinski, that solved my question.

-- 
You received this question notification because your team yade-users is
an answer contact for Yade.

___
Mailing list: https://launchpad.net/~yade-users
Post to : yade-users@lists.launchpad.net
Unsubscribe : https://launchpad.net/~yade-users
More help   : https://help.launchpad.net/ListHelp


Re: [Yade-users] [Question #707882]: Cylindrical triaxial test using walls

2023-09-17 Thread Jérôme Duriez
Question #707882 on Yade changed:
https://answers.launchpad.net/yade/+question/707882

Status: Open => Answered

Jérôme Duriez proposed the following answer:
Hi,

A Wall [*] shaped-body is an infinite plane (or infinite half-domain
depending on Wall.sense), so no, it can not correspond to a cylindrical
wall.

 If this is really what you want, you have to go through Facet shapes
indeed, but you should note there is no Ig2_Facet_LevelSet_* functor
available at the moment (and interaction with LevelSet-shaped bodies
will not be possible)


[*] https://yade-dem.org/doc/yade.wrapper.html#yade.wrapper.Wall

-- 
You received this question notification because your team yade-users is
an answer contact for Yade.

___
Mailing list: https://launchpad.net/~yade-users
Post to : yade-users@lists.launchpad.net
Unsubscribe : https://launchpad.net/~yade-users
More help   : https://help.launchpad.net/ListHelp


Re: [Yade-users] [Question #698464]: asphalt mixture model, contact, code

2023-09-17 Thread isadjkn123
Question #698464 on Yade changed:
https://answers.launchpad.net/yade/+question/698464

isadjkn123 proposed the following answer:

Hi,

Asphalt mixture modeling is a complex topic, but it is very important
for designing and building durable and safe roads.

To arrange asphalt mortar in the form of code, you can use a variety of
methods, such as discrete element methods (DEM) simulators or continuum
models.

To establish the contact between particles in asphalt mortar, you can
use a variety of contact models, such as the Hertz model or the Johnson-
Kendall-Roberts (JKR) model.

The best method and contact model for your simulations will depend on
the specific problem you are trying to solve.

If you need help with asphalt calculations, you can use an
https://asphaltcalculators.com/

-- 
You received this question notification because your team yade-users is
an answer contact for Yade.

___
Mailing list: https://launchpad.net/~yade-users
Post to : yade-users@lists.launchpad.net
Unsubscribe : https://launchpad.net/~yade-users
More help   : https://help.launchpad.net/ListHelp


Re: [Yade-users] [Question #707955]: Unexpected Termination of Yade Script After Particle Deletion

2023-09-17 Thread Masataka Shiga
Question #707955 on Yade changed:
https://answers.launchpad.net/yade/+question/707955

Masataka Shiga posted a new comment:
I did not notice the #707886 announcement. I will move this question to gitlab.
Many thanks.

-- 
You received this question notification because your team yade-users is
an answer contact for Yade.

___
Mailing list: https://launchpad.net/~yade-users
Post to : yade-users@lists.launchpad.net
Unsubscribe : https://launchpad.net/~yade-users
More help   : https://help.launchpad.net/ListHelp


[Yade-users] [Question #707955]: Unexpected Termination of Yade Script After Particle Deletion

2023-09-17 Thread Masataka Shiga
New question #707955 on Yade:
https://answers.launchpad.net/yade/+question/707955

Hello Yade Community,

I'm working on a simulation to create a sample through natural sedimentation. 
Here's a brief description of my issue:

I'm trying to create a specimen through natural sedimentation.
By setting O.periodic=True, I've established periodic boundary conditions 
laterally and fixed boundary conditions vertically using a Facet whose position 
is fixed at the bottom.
I generate particles at a height well above the target, allow them to settle 
under gravity, and once the average coordination number surpasses 3, I aim to 
delete particles above a specific height to achieve a group of particles at the 
desired height.
However, after deleting the particles using the code below, the program stops 
after a certain number of iterations (ranging from a few tens to thousands). 
There's no error message displayed in the terminal.
Moreover, the iteration at which the program stops varies with each script 
execution.

I'm puzzled by this unexpected behavior, especially given the absence of any 
error messages. Any insights or suggestions would be greatly appreciated.

Thank you in advance for your assistance.

# Environment:

Execution Command: yade -j16 code.py
Yade Version: 2022.01a
Linux Version: Ubuntu 22.04.3 LTS
Architecture: amd64

# Script (MWE):

from yade import pack, utils
import numpy as np

state_index = 0

pack_sp = pack.SpherePack()
pack_sp.makeCloud((0.02, 0.02, 0.02), 
  (0.08, 0.98, 0.08), 
  rMean=0.01)
sphere_id = pack_sp.toSimulation()

O.periodic = True

base_facet_1 = utils.facet([Vector3(-0.02, 0, -0.02), 
Vector3(0.12,  0, -0.02),
Vector3(-0.02, 0, 0.12)],
fixed=True)
base_facet_2 = utils.facet([Vector3(0.12, 0,  0.12), 
Vector3(0.12,  0, -0.02),
Vector3(-0.02, 0, 0.12)],
fixed=True)

base_facet_id = O.bodies.append([base_facet_1, base_facet_2])

O.cell.hSize = Matrix3(0.1, 0, 0,
   0, 1, 0,
   0, 0, 0.1)
O.periodic = True

O.engines = [
ForceResetter(),
InsertionSortCollider(
[Bo1_Sphere_Aabb(), Bo1_Facet_Aabb()], allowBiggerThanPeriod=True),
InteractionLoop(
[Ig2_Sphere_Sphere_ScGeom(), Ig2_Facet_Sphere_ScGeom()],
[Ip2_FrictMat_FrictMat_FrictPhys()],
[Law2_ScGeom_FrictPhys_CundallStrack()]),
NewtonIntegrator(damping=0.2, gravity=(0, -9.81, 0)),
PyRunner(iterPeriod=100, command="checkState()")
]

O.dt = PWaveTimeStep() * 0.5

def checkState():
global sphere_id, state_index

if state_index == 0:

temp_sphere_coord_num = [len(O.interactions.withBody(i)) for i in 
sphere_id]
temp_sphere_Y = [O.bodies[i].state.pos[1] for i in sphere_id]
temp_sphere_data = np.array([sphere_id, temp_sphere_coord_num, 
temp_sphere_Y])

if temp_sphere_data[1, :].mean() >=3:

temp_unused_sphere_id = temp_sphere_data[0, temp_sphere_data[2, :] 
> 0.1]

for temp_unused_sphere_id_each in temp_unused_sphere_id:
O.bodies.erase(int(temp_unused_sphere_id_each))
sphere_id.remove(int(temp_unused_sphere_id_each))

print("some particles are deleted at", O.iter)

state_index = 1

elif state_index == 1:
pass


-- 
You received this question notification because your team yade-users is
an answer contact for Yade.

___
Mailing list: https://launchpad.net/~yade-users
Post to : yade-users@lists.launchpad.net
Unsubscribe : https://launchpad.net/~yade-users
More help   : https://help.launchpad.net/ListHelp