Re: [Yade-users] [Question #706850]: Clumping can't give color to each layer
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?
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
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
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
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
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
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
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
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
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
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
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
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?
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
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
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