Re: [Yade-users] [Question #706850]: Clumping can't give color to each layer

2023-05-31 Thread Jan Stránský
Question #706850 on Yade changed:
https://answers.launchpad.net/yade/+question/706850

Status: Open => Answered

Jan Stránský proposed the following answer:
Hello,

> sphere_obj.color = color

sphere_obj.shape.color = color # [1]

Cheers
Jan

[1] https://yade-dem.org/doc/yade.wrapper.html#yade.wrapper.Shape.color

-- 
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 #706848]: Constant pressure or Constant volume?

2023-05-31 Thread Jan Stránský
Question #706848 on Yade changed:
https://answers.launchpad.net/yade/+question/706848

Status: Open => Needs information

Jan Stránský requested more information:
Hello,

> I found out that this example program is "constant volume" periodic simple 
> shear.
> So we would like to perform "constant volume" periodic simple shear.
> Is there any program ("constant volume" periodic simple shear) that can help 
> us?

So what is the problem?
How did that example did not help?
What other help do you want (please be more specific)?

Cheers
Jan

-- 
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 #706851]: FreeGlut versus FREEGLUT package

2023-05-31 Thread Luc OGER
New question #706851 on Yade:
https://answers.launchpad.net/yade/+question/706851

hi,
during my cmake process to generate the Makefile for the latest yade version 
2023-05-15 , I got a Warning about freeglut:
CMake Warning (dev) at 
/usr/share/cmake/Modules/FindPackageHandleStandardArgs.cmake:438 (message):
  The package name passed to `find_package_handle_standard_args` (FREEGLUT)
  does not match the name of the calling package (FreeGlut).  This can lead
  to problems in calling code that expects `find_package` result variables
  (e.g., `_FOUND`) to follow a certain pattern.
Call Stack (most recent call first):
  cMake/FindFreeGlut.cmake:14 (FIND_PACKAGE_HANDLE_STANDARD_ARGS)
  CMakeLists.txt:898 (FIND_PACKAGE)
This warning is for project developers.  Use -Wno-dev to suppress it.

-- Found FREEGLUT: /usr/include/GL

I don't know if it is critical but, right now the GUI display is no more 
working (i.e the yade --check option is working!)

for example in debug mode the first thread is :
Thread 1 (Thread 0x7f04fc3a4dc0 (LWP 109253) "yade"):
#0  0x7f04fb7c574f in wait4 () from /lib64/libc.so.6
#1  0x7f04fb73cba7 in do_system () from /lib64/libc.so.6
#2  0x7f04f8a3d73e in crashHandler (sig=11) at 
/home/oger/yade_full_install_local/trunk/core/main/pyboot.cpp:16
#3  
#4  0x7f04fbdb5efb in PyModule_GetState () from 
/usr/lib64/libpython3.6m.so.1.0
#5  0x7f0484192783 in ?? () from 
/usr/lib64/python3.6/lib-dynload/readline.cpython-36m-x86_64-linux-gnu.so
#6  0x7f04895b2e73 in rl_initialize () from /usr/lib64/../lib64/libedit.so.0
#7  0x7f0484192909 in PyInit_readline () from 
/usr/lib64/python3.6/lib-dynload/readline.cpython-36m-x86_64-linux-gnu.so
#8  0x7f04fbe44ba3 in _PyImport_LoadDynamicModuleWithSpec () from 
/usr/lib64/libpython3.6m.so.1.0
#9  0x7f04fbe43031 in ?? () from /usr/lib64/libpython3.6m.so.1.0
#10 0x7f04fbdb4da5 in PyCFunction_Call () from 
/usr/lib64/libpython3.6m.so.1.0
#11 0x7f04fbe2caca in _PyEval_EvalFrameDefault () from 
/usr/lib64/libpython3.6m.so.1.0
#12 0x7f04fbe25fd6 in ?? () from /usr/lib64/libpython3.6m.so.1.0
#13 0x7f04fbe2624f in ?? () from /usr/lib64/libpython3.6m.so.1.0
#14 0x7f04fbe26455 in ?? () from /usr/lib64/libpython3.6m.so.1.0
#15 0x7f04fbe2a9d5 in _PyEval_EvalFrameDefault () from 
/usr/lib64/libpython3.6m.so.1.0
#16 0x7f04fbe2617f in ?? () from /usr/lib64/libpython3.6m.so.1.0
#17 0x7f04fbe26455 in ?? () from /usr/lib64/libpython3.6m.so.1.0
#18 0x7f04fbe2a9d5 in _PyEval_EvalFrameDefault () from 
/usr/lib64/libpython3.6m.so.1.0
#19 0x7f04fbe2617f in ?? () from /usr/lib64/libpython3.6m.so.1.0
#20 0x7f04fbe26455 in ?? () from /usr/lib64/libpython3.6m.so.1.0
#21 0x7f04fbe2a9d5 in _PyEval_EvalFrameDefault () from 
/usr/lib64/libpython3.6m.so.1.0
#22 0x7f04fbe2617f in ?? () from /usr/lib64/libpython3.6m.so.1.0
#23 0x7f04fbe26455 in ?? () from /usr/lib64/libpython3.6m.so.1.0
#24 0x7f04fbe2a9d5 in _PyEval_EvalFrameDefault () from 
/usr/lib64/libpython3.6m.so.1.0
#25 0x7f04fbe2617f in ?? () from /usr/lib64/libpython3.6m.so.1.0
#26 0x7f04fbe26455 in ?? () from /usr/lib64/libpython3.6m.so.1.0
#27 0x7f04fbe2a9d5 in _PyEval_EvalFrameDefault () from 
/usr/lib64/libpython3.6m.so.1.0
#28 0x7f04fbe2617f in ?? () from /usr/lib64/libpython3.6m.so.1.0
#29 0x7f04fbe2e370 in _PyFunction_FastCallDict () from 
/usr/lib64/libpython3.6m.so.1.0
#30 0x7f04fbd6c9ee in _PyObject_FastCallDict () from 
/usr/lib64/libpython3.6m.so.1.0
#31 0x7f04fbd6d8c4 in _PyObject_CallMethodIdObjArgs () from 
/usr/lib64/libpython3.6m.so.1.0
#32 0x7f04fbe43d9c in PyImport_ImportModuleLevelObject () from 
/usr/lib64/libpython3.6m.so.1.0
#33 0x7f04fbe2c79f in _PyEval_EvalFrameDefault () from 
/usr/lib64/libpython3.6m.so.1.0
#34 0x7f04fbe25fd6 in ?? () from /usr/lib64/libpython3.6m.so.1.0
#35 0x7f04fbe2624f in ?? () from /usr/lib64/libpython3.6m.so.1.0
#36 0x7f04fbe26455 in ?? () from /usr/lib64/libpython3.6m.so.1.0
#37 0x7f04fbe2928e in _PyEval_EvalFrameDefault () from 
/usr/lib64/libpython3.6m.so.1.0
#38 0x7f04fbe25fd6 in ?? () from /usr/lib64/libpython3.6m.so.1.0
#39 0x7f04fbe26113 in PyEval_EvalCode () from 
/usr/lib64/libpython3.6m.so.1.0
#40 0x7f04fbe4fa60 in ?? () from /usr/lib64/libpython3.6m.so.1.0
#41 0x7f04fbe51891 in PyRun_FileExFlags () from 
/usr/lib64/libpython3.6m.so.1.0
#42 0x7f04fbe519ed in PyRun_SimpleFileExFlags () from 
/usr/lib64/libpython3.6m.so.1.0
#43 0x7f04fbe68a91 in Py_Main () from /usr/lib64/libpython3.6m.so.1.0
#44 0x562b12a00c90 in main ()
[Inferior 1 (process 109253) detached]
Erreur de segmentation (core dumped)





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   : ht

Re: [Yade-users] [Question #706845]: Young's modulus and normal stiffness

2023-05-31 Thread Jan Stránský
Question #706845 on Yade changed:
https://answers.launchpad.net/yade/+question/706845

Status: Open => Needs information

Jan Stránský requested more information:
Hello,

> What I am gonna ask is a simple question.

Actually it can be very complex :-)

> If I know the normal and shear stiffness kn and ks of spheres, how can
I determine the corresponding Young's modulus E and Poisson's ratio v? E
= kn * D (radius of sphere)?

What is "the normal and shear stiffness kn and ks of spheres"?
What materials and contact law(s) are you using?
What is "the corresponding Young's modulus E and Poisson's ratio v"?

Cheers
Jan

-- 
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 #706843]: VTKRecorder throws Segmentation fault error

2023-05-31 Thread Jan Stránský
Question #706843 on Yade changed:
https://answers.launchpad.net/yade/+question/706843

Status: Open => Needs information

Jan Stránský requested more information:
Thanks for informations.

I have tried it without GUI
yade -n script.py
without error.
Could you please test it?

Cheers
Jan

-- 
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 #706854]: Code availability of Yuan2016, Sweijen2016

2023-05-31 Thread Sebastian Falcioni
New question #706854 on Yade:
https://answers.launchpad.net/yade/+question/706854

Hello,

I noticed the 2 following papers:
- Yuan2016: Pore-scale simulations of drainage in granular materials: Finite 
size effects and the representative elementary volume by C. Yuan et al., in 
Advances in Water Resources (2016)
- Sweijen2016: The Effects of Swelling and Porosity Change on Capillarity: DEM 
Coupled with a Pore-Unit Assembly Method by Sweijen T. et al., in Transport in 
Porous Media (2016)

and would like to replicate the kind of SWCC simulations presented therein (I 
would already be happy with one of the two methods).

Is the following code yet available ? If not, is there some timeline for 
release ? Do you have the link for the code?

Thanks,

Sebastian

-- 
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 #403166]: Code availability of Yuan2016, Sweijen2016

2023-05-31 Thread Sebastian Falcioni
Question #403166 on Yade changed:
https://answers.launchpad.net/yade/+question/403166

Sebastian Falcioni posted a new comment:
Hi.

Do you know if the code is available? Do you have it?

Thanks.
Sebastian.

-- 
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 #705052]: Particle Water Absorption

2023-05-31 Thread Sebastian Falcioni
Question #705052 on Yade changed:
https://answers.launchpad.net/yade/+question/705052

Sebastian Falcioni posted a new comment:
Hi fellas.

Can you access to the code of this article?

Thanks.
Sebastian.

-- 
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 #706854]: Code availability of Yuan2016, Sweijen2016

2023-05-31 Thread Karol Brzezinski
Question #706854 on Yade changed:
https://answers.launchpad.net/yade/+question/706854

Karol Brzezinski posted a new comment:
Hi Sebastian,

Hehe you looked in the examples? I checked the examples in the main repository, 
and
this one seem to be the closest to your problem [1].

Cheers,
Karol

[1] https://gitlab.com/yade-
dev/trunk/-/blob/master/examples/FluidCouplingPFV/drainage-2PFV-
Yuan_and_Chareyre_2017.py

-- 
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 #706854]: Code availability of Yuan2016, Sweijen2016

2023-05-31 Thread Sebastian Falcioni
Question #706854 on Yade changed:
https://answers.launchpad.net/yade/+question/706854

Sebastian Falcioni posted a new comment:
Hi Karol.

Thanks for the quick answer. I'm trying to model hydrogel swelling with
DEM. It's seems to be comparable but if you know another example closest
to mine I'll thank you. Anyway... do you have the video of the example
code?

Thanks a lot.
Sebastian.

-- 
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 #706855]: Segmentation fault (core dumped) when increasing the GridConnection number

2023-05-31 Thread Wang Yaqiong
New question #706855 on Yade:
https://answers.launchpad.net/yade/+question/706855

I have a question when generating the geogrid sample. 

When I increased the number of nodes (e.g., "nbL" from 10 to 60), the 
"Segmentation fault " appears in the terminal. I have confirmed that the 
program woks normally if only unidirectional ribs are generated.

My script is attched below:
___
# -*- encoding=utf-8 -*-
from builtins import zip
from builtins import range
from yade.gridpfacet import *
from yade import pack,ymport,export,geom,bodiesHandling, plot,utils,qt
import pylab
import matplotlib; matplotlib.rc('axes',grid=True)

 tensile control
leftTensilePos = 0.002
rightTensilePos = 0.098
O.engines = [
ForceResetter(),
InsertionSortCollider([
Bo1_GridConnection_Aabb(),
]),
InteractionLoop(
[Ig2_GridNode_GridNode_GridNodeGeom6D()], 
[Ip2_CohFrictMat_CohFrictMat_CohFrictPhys(setCohesionNow=True, 
setCohesionOnNewContacts=False)],
[Law2_ScGeom6D_CohFrictPhys_CohesionMoment()]
),
NewtonIntegrator(gravity=(0, 0, -10), damping=0.1, label='newton'),
#VTKRecorder(iterPeriod=1, 
recorders=['spheres','facets','boxes','stress','colors','intr'], 
fileName='vtkrecorder/p1-'),
PyRunner(command='history()',iterPeriod=1),
PyRunner(command='plotAddData1()',iterPeriod=100),

]


spherematid = O.materials.append(
CohFrictMat(
young=1.74e9,#1.74e9,
alphaKr=2,
alphaKtw=2,
poisson=0.5,
density=800,
frictionAngle=radians(10),
normalCohesion=9e9,
shearCohesion=9e9,
momentRotationLaw=True,
label='spheremat'
)
)
print('spherematid',spherematid)
### Parameters of a rectangular grid ###

L = 2.1  #length [m]   aperture size = 30 mm, grid_gap=35 mm
l = 2.1  #width [m]
nbL =60  #number of grids for the length[#]
nbl =60  #number of grids for the width [#]
nbLL = nbL+1
nbll = nbl+1 #number of nodes for the width
rrr = 0.0025 #L / 100.  #radius
color = [255. / 255., 102. / 255., 0. / 255.]
nodesIds1 = []
#Create all nodes first :
for i in range(0, nbLL):
for j in range(0, nbll):
nodesIds1.append(O.bodies.append(gridNode([i * L / nbL, j * l / 
nbl, 0], rrr, wire=True, fixed=False, material='spheremat', color=color)))

#print(nodesIds[:-1])
ConnectionIds1 = []
#Create connection between the nodes
for i in range(0, len(nodesIds1)):
for j in range(i + 1, len(nodesIds1)):
dist = (O.bodies[i].state.pos - O.bodies[j].state.pos).norm()
if (dist <= l / nbl * 1.01):
ConnectionIds1.append(O.bodies.append(gridConnection(i, 
j, rrr, color=color))) #
print(ConnectionIds1)   
leftNode1 = [O.bodies[s] for s in nodesIds1 if 
O.bodies[s].state.pos[0]rightTensilePos]#right node

for s in leftNode1:
s.shape.color = (1,0,0)
s.dynamic = False

vel = 0

def history():
for s in leftNode1:
s.shape.color = (1,0,0)
s.state.blockedDOFs = 'xyzXYZ'
s.state.vel = (-vel,0,0)


def plotAddData1():
f_leftNode1 = sum(O.forces.f(b.id)[0] for b in leftNode1)
f_rightNode1 = sum(O.forces.f(b.id)[0] for b in rightNode1)

f_both1 = .5*(-f_rightNode1+f_leftNode1)
f1 =f_leftNode1
strain1 = (rightNode1[0].state.displ()[0] - 
leftNode1[0].state.displ()[0]) / (mx-mm)*100 #[%]
height1=rightNode1[0].state.refPos[2]-rightNode1[0].state.pos[2]
plot.addData(
i = O.iter,
f_both1 = f_both1,
height1 = height1,
strain1 = strain1,
)

plot.plots = { 'i':('height1',), }

plot.plot()
O.dt =  0.5*PWaveTimeStep()
O.saveTmp()
qt.View()

-- 
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 #706855]: Segmentation fault (core dumped) when increasing the GridConnection number

2023-05-31 Thread Wang Yaqiong
Question #706855 on Yade changed:
https://answers.launchpad.net/yade/+question/706855

Wang Yaqiong posted a new comment:
I found that if I change "O.saveTmp()" into “O.run()”, no "segmentation
fault" is printed and the program works well.

The problem seems to be fixed.

But I don't know that's why.

-- 
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 #706843]: VTKRecorder throws Segmentation fault error

2023-05-31 Thread Yuri Bezmenov
Question #706843 on Yade changed:
https://answers.launchpad.net/yade/+question/706843

Status: Needs information => Open

Yuri Bezmenov gave more information on the question:
hi Jan, 
I am still facing the same issue. 

baadalvm@baadalvm:~/file/triax$ yade -n new.py
Welcome to Yade 2022.01a 
Using python version: 3.10.6 (main, Mar 10 2023, 10:55:28) [GCC 11.3.0]
TCP python prompt on localhost:9000, auth cookie `ucseks'
XMLRPC info provider on http://localhost:21000
Running script new.py
unbalanced force: 0.001  mean stress:  -13.73470581507172
unbalanced force: 0.0954468397211678  mean stress:  -24701.39828766909
unbalanced force: 0.0006870880148087554  mean stress:  -9.34510748321
###  Isotropic state reached  ###
Segmentation fault (core dumped)

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


Re: [Yade-users] [Question #706848]: Constant pressure or Constant volume?

2023-05-31 Thread 内山康太郎
Question #706848 on Yade changed:
https://answers.launchpad.net/yade/+question/706848

内山康太郎 posted a new comment:
Thanks for the reply.

I'm sorry, Mr Jan.
I made a mistake in my description.

What I am looking for is a "constant pressure“ periodic simple shear.

So we would like to perform "constant pressure" periodic simple shear

-- 
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 #706106]: Load the data to plot

2023-05-31 Thread alaia woods
Question #706106 on Yade changed:
https://answers.launchpad.net/yade/+question/706106

alaia woods posted a new comment:
To plot the results from a file and compare them with DEM (Discrete
Element Method) simulations in real-time, you can follow these steps:

Load the data from your file: You'll need to determine the format of
your data file (e.g., CSV, TXT) and use appropriate libraries or
functions to read the data into your program. For examplehttps://laptopspapa.com/";>, if you're using Python, you can
use libraries like NumPy or Pandas to load and manipulate the data.

Prepare the data for comparison: Depending on the structure of your
data, you may need to preprocess or reformat it to align with the DEM
simulation data. This may involve organizing the data into appropriate
variables or arrays.

Plot the data: Utilize a plotting library such as Matplotlib or Plotly
to visualize the loaded data. These libraries provide a range of
functions and customization options for creating plots, such as line
plots, scatter plots, or 3D plots. You can plot the loaded data
alongside the DEM simulation results to compare them in real-time.

Update the plot in real-time: If you want to continuously compare the
loaded data with the evolving DEM simulation, you'll need to update the
plot dynamically. This can be achieved by periodically loading new data
from the file and updating the plot accordingly. You may need to
incorporate a loop or an event-driven approach to enable real-time
updates

-- 
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 #706843]: VTKRecorder throws Segmentation fault error

2023-05-31 Thread Yuri Bezmenov
Question #706843 on Yade changed:
https://answers.launchpad.net/yade/+question/706843

Yuri Bezmenov posted a new comment:
Hi Jan, 
I just update Ubuntu and again run 
  baadalvm@baadalvm:~/file/triax$ yade -n new.py
this time I didnt got ant error. 


just one more thing, could you please explain this to me. I didn't got why it 
doesn't work earlier. 

Thank You

-- 
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