Re: [Yade-users] [Question #701564]: Install yade in the centos7 system. After executing the make install command, an error is reported when compiling FlowEngine.cpp
Question #701564 on Yade changed: https://answers.launchpad.net/yade/+question/701564 xuanshenyu gave more information on the question: Thank you for your reply. The two posts are not the same, and I closed the other post (question 701378). I have configured it by referring to the link[1] , but an error is reported when compiling FlowEngine.cpp (after executing the make install command). [1]https://gitlab.com/yade-dev/docker-yade/-/issues/5 -- 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 #701378]: install Yade in centos7 system
Question #701378 on Yade changed: https://answers.launchpad.net/yade/+question/701378 Status: Open => Solved xuanshenyu confirmed that the question is solved: Thanks, I will refer to the link to configure and compile. -- 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 #701566]: Two material
Question #701566 on Yade changed: https://answers.launchpad.net/yade/+question/701566 Jan Stránský proposed the following answer: Hello, > My simulation works with single material please provide the code > the difference is the radius and the density. only? what about the makeCloud? > sp_PH101=pack.SpherePack() > sp_PH101.makeCloud(...) > sp_PH101.toSimulation(material=matPH101,color=(0,1,1)) > sp_CCS=pack.SpherePack() > sp_CCS.makeCloud(...) > sp_CCS.toSimulation(material=matCCS,color=(1,0,1)) Why do you use two SpherePack instances? The you get independent clouds possibly with overlapping spheres pushing them with non-realistic speed, possibly causing "flying away". Try with one SpherePack instance, then spheres are not overlapping. Something like: ### sp=pack.SpherePack() sp.makeCloud(...) n1 = len(sp) sp.makeCloud(...) for i,(c,r) in snumerate(sp): mat = matPH101 if i < n1 else matCSS color = (0,1,1) if i < n1 else (1,0,1) O.bodies.append(sphere(c,r,material=mat,color=color) ### 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 #701566]: Two material
Question #701566 on Yade changed: https://answers.launchpad.net/yade/+question/701566 Status: Open => Answered Karol Brzezinski proposed the following answer: Hi, Without running the code I am pretty sure that it is not a MWE (Minimum Working Example). If yoYou need to know that proper dt for stable simulation depends on material density, stiffness, and particle size. You have fixed dt, try to change it as follows: O.dt = 0.8*PWaveTimestep() Cheers Karol -- 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 #701566]: Two material
New question #701566 on Yade: https://answers.launchpad.net/yade/+question/701566 Hi, I working on compaction, with two different material. The material type is the same, the difference is the radius and the density. My simulation works with single materials, however when use two material the particle began to fly away during gravity filling. Best regards, MWE: from __future__ import print_function from yade import utils, plot, timing from yade import pack import pandas as pd import numpy as np from yade import pack, export import os readParamsFromTable(comp_press=0.81e8,h_tab=2.02,m_tab=0.2117, r_tab=5.015,wCCS=0.02,tab_porosity=15,tab_height=1,save=1) from yade.params.table import * o = Omega() save=save # Physical parameters fr = 0.41 rho_PH101 = 1561 rho_CCS =1403 D_PH101 = 7.9e-5 r1_PH101 = D_PH101/2 D_CCS = 5.4e-5 r1_CCS = D_CCS/2 #r2 = Diameter/2 k1 = 1 kp = 14 kc = k1 * 0.1 ks = k1 * 0.1 Chi1 = 0.34 o.dt = 1.0e-8 PhiF1=0.999 Cyl_height=0.008 cross_area=math.pi*(Tab_rad**2) Comp_press_up= comp_press Comp_force_up=Comp_press_up*cross_area Comp_press_lp= comp_press Comp_force_lp=Comp_press_lp*cross_area # Add material matPH101 = O.materials.append(LudingMat(frictionAngle=fr, density=rho_PH101, k1=k1, kp=kp, ks=ks, kc=kc, PhiF=PhiF1, G0 = 0.0)) matCCS = O.materials.append(LudingMat(frictionAngle=fr, density=rho_CCS, k1=k1, kp=kp, ks=ks, kc=kc, PhiF=PhiF1, G0 = 0.0)) # Spheres for compression and walls sp_PH101=pack.SpherePack() sp_PH101.makeCloud((-8*D_PH101,-8*D_PH101,-50*D_PH101),(8*D_PH101,8*D_PH101,40.0*D_PH101), rMean=r1_PH101,rRelFuzz=0.18,num=round(10134)) sp_PH101.toSimulation(material=matPH101,color=(0,1,1)) sp_CCS=pack.SpherePack() sp_CCS.makeCloud((-8*D_PH101,-8*D_PH101,-50*D_PH101),(8*D_PH101,8*D_PH101,40.0*D_PH101), rMean=r1_CCS,rRelFuzz=0.15,num=round(740)) sp_CCS.toSimulation(material=matCCS,color=(1,0,1)) walls=O.bodies.append(yade.geom.facetCylinder((0,0,0),radius=Tab_rad,height=Cyl_height,segmentsNumber=20,wallMask=6,material=matPH101)) # Add engines o.engines = [ ForceResetter(), InsertionSortCollider([Bo1_Sphere_Aabb(aabbEnlargeFactor=1.05), Bo1_Wall_Aabb(), Bo1_Facet_Aabb() ]), InteractionLoop( [Ig2_Sphere_Sphere_ScGeom(interactionDetectionFactor=1.05), Ig2_Facet_Sphere_ScGeom(), Ig2_Wall_Sphere_ScGeom()], [Ip2_LudingMat_LudingMat_LudingPhys()], [Law2_ScGeom_LudingPhys_Basic()] ), NewtonIntegrator(damping=0.1, gravity=[0, 0,-9.81]), PyRunner(command='checkForce()', realPeriod=1, label="fCheck"), #DeformControl(label="DefControl") ] def checkForce(): if O.iter < 300: return if unbalancedForce() > 1: return global upper_punch upper_punch=O.bodies.append(geom.facetCylinder((0,0,((-Cyl_height/2)+0.0001)+utils.aabbDim()[2]),Tab_rad-.1,0,segmentsNumber=50,wallMask=1)) for i in upper_punch: body= O.bodies[i] body.state.vel = (0,0,-0.02) global lower_punch lower_punch= O.bodies.append(geom.facetCylinder((0,0,(-Cyl_height/2)-0.0001),Tab_rad-.1,0,segmentsNumber=50,wallMask=1)) for n in lower_punch: body= O.bodies[n] body.state.vel = (0,0,0.02) O.engines = O.engines + [PyRunner(command='storeData()', iterPeriod=25000)]+ [PyRunner(command='saveData()', iterPeriod=10)] fCheck.command = 'unloadPlate()' -- 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 #701564]: Install yade in the centos7 system. After executing the make install command, an error is reported when compiling FlowEngine.cpp
Question #701564 on Yade changed: https://answers.launchpad.net/yade/+question/701564 Janek Kozicki posted a new comment: That looks like a duplicate of https://answers.launchpad.net/yade/+question/701378 -- 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 #701564]: Install yade in the centos7 system. After executing the make install command, an error is reported when compiling FlowEngine.cpp
New question #701564 on Yade: https://answers.launchpad.net/yade/+question/701564 Hi,everyone I am compiling and installing yade in the centos7 system. Refer to [1]. After executing the make install command,some errors are reported when compiling FlowEngine.cpp. This seems to be related to the design of FlowEngine.cpp. itself . In trunk-2022.01a\pkg\pfv, there are three FlowEngine related header files (FlowEngine.hpp FlowEngine.hpp.in FlowEngine.ipp.in). Please let me know if full file is required or if you find other errors in the script. Thanks in advance for your help! Error Report### [ 3%] Building CXX object CMakeFiles/pkg_pfv.dir/pkg/pfv/DummyFlowEngine.cpp.o [ 3%] Building CXX object CMakeFiles/pkg_pfv.dir/pkg/pfv/FlowEngine.cpp.o In file included from /opt/boost/include/boost/python/detail/copy_ctor_mutates_rhs.hpp:8, from /opt/boost/include/boost/python/detail/value_arg.hpp:7, from /opt/boost/include/boost/python/object/forward.hpp:10, from /opt/boost/include/boost/python/object/pointer_holder.hpp:16, from /opt/boost/include/boost/python/to_python_indirect.hpp:10, from /opt/boost/include/boost/python/converter/arg_to_python.hpp:10, from /opt/boost/include/boost/python/call.hpp:15, from /opt/boost/include/boost/python/object_core.hpp:14, from /opt/boost/include/boost/python/args.hpp:22, from /opt/boost/include/boost/python.hpp:11, from /home/xsy/whut_yade/trunk/lib/high-precision/RealHPConfig.hpp:35, from /home/xsy/whut_yade/trunk/lib/high-precision/Real.hpp:238, from /home/xsy/whut_yade/trunk/lib/base/Math.hpp:12, from /home/xsy/whut_yade/trunk/core/Bound.hpp:10, from /home/xsy/whut_yade/trunk/core/Body.hpp:13, from /home/xsy/whut_yade/trunk/core/PartialEngine.hpp:10, from /home/xsy/whut_yade/build/pkg/pfv/FlowEngine_FlowEngineT.hpp:36, from /home/xsy/whut_yade/trunk/pkg/pfv/FlowEngine.cpp:11: /opt/boost/include/boost/python/detail/is_auto_ptr.hpp:17:105: warning: ‘template class std::auto_ptr’ is deprecated [-Wdeprecated-declarations] 17 | BOOST_PYTHON_IS_XXX_DEF(auto_ptr, std::auto_ptr, 1) | ^ In file included from /opt/rh/devtoolset-10/root/usr/include/c++/10/bits/locale_conv.h:41, from /opt/rh/devtoolset-10/root/usr/include/c++/10/locale:43, from /opt/rh/devtoolset-10/root/usr/include/c++/10/iomanip:43, from /opt/boost/include/boost/multiprecision/cpp_int.hpp:10, from /opt/boost/include/boost/multiprecision/cpp_bin_float.hpp:9, from /home/xsy/whut_yade/trunk/lib/high-precision/RealHP.hpp:51, from /home/xsy/whut_yade/trunk/lib/high-precision/Real.hpp:210, from /home/xsy/whut_yade/trunk/lib/base/Math.hpp:12, from /home/xsy/whut_yade/trunk/core/Bound.hpp:10, from /home/xsy/whut_yade/trunk/core/Body.hpp:13, from /home/xsy/whut_yade/trunk/core/PartialEngine.hpp:10, from /home/xsy/whut_yade/build/pkg/pfv/FlowEngine_FlowEngineT.hpp:36, from /home/xsy/whut_yade/trunk/pkg/pfv/FlowEngine.cpp:11: /opt/rh/devtoolset-10/root/usr/include/c++/10/bits/unique_ptr.h:57:28: note: declared here 57 | template class auto_ptr; |^~~~ In file included from /home/xsy/whut_yade/build/pkg/pfv/FlowEngine_FlowEngineT.hpp:42, from /home/xsy/whut_yade/trunk/pkg/pfv/FlowEngine.cpp:11: /home/xsy/whut_yade/trunk/lib/triangulation/PartialSatLinSolv.hpp:22:90: error: expected template-name before ‘<’ token 22 | elation> class PartialSatLinSolv : public FlowBoundingSphereLinSolv<_Tesselation, FlowBoundingSphere<_Tesselation>> { |^ /home/xsy/whut_yade/trunk/lib/triangulation/PartialSatLinSolv.hpp:22:90: error: expected ‘{’ before ‘<’ token In file included from /home/xsy/whut_yade/trunk/lib/triangulation/PartialSatLinSolv.hpp:193, from /home/xsy/whut_yade/build/pkg/pfv/FlowEngine_FlowEngineT.hpp:42, from /home/xsy/whut_yade/trunk/pkg/pfv/FlowEngine.cpp:11: /home/xsy/whut_yade/trunk/lib/triangulation/PartialSatLinSolv.ipp:30:84: error: invalid use of incomplete type ‘class yade::CGT::PartialSatLinSolv<_Tesselation>’ 30 | ate PartialSatLinSolv<_Tesselation>::~PartialSatLinSolv() | ^ In file included from /home/xsy/whut_yade/build/pkg/pfv/FlowEngine_FlowEngineT.hpp:42, from
Re: [Yade-users] [Question #701438]: Why particle velocity varies with time step
Question #701438 on Yade changed: https://answers.launchpad.net/yade/+question/701438 Aoxi Zhang posted a new comment: Hi Karol, Sorry for not making it clear. I mean the damping is changing properly, but the problem is not solved and I still need help on it. Thanks, Aoxi -- 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 #701438]: Why particle velocity varies with time step
Question #701438 on Yade changed: https://answers.launchpad.net/yade/+question/701438 Karol Brzezinski posted a new comment: Hi Zhang, Do you mean 'work well', so the damping is changing properly, but you still need help or work well and the problem is solved? Cheers, Karol -- 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 #701341]: About ValueError: cannot convert float NaN to integer
Question #701341 on Yade changed: https://answers.launchpad.net/yade/+question/701341 Status: Needs information => Open Ziyu Wang gave more information on the question: Hello, Due to my personal reasons, I can't feed back the problem in time. I'm sorry. >What is your yadedaily version? What it writes when you start it? As described below: Welcome to Yade 20220415-6437~b9e47eb~bionic1 Using python version: 3.6.9 (default, Mar 15 2022, 13:55:28) [GCC 8.4.0] TCP python prompt on localhost:9001, auth cookie `csedak' XMLRPC info provider on http://localhost:21001 [[ ^L clears screen, ^U kills line. F12 controller, F11 3D view (press "h" in 3D view for help), F10 both, F9 generator, F8 plot. ]] >no, sorry, I mean just to change the position of the walls. Yes,according to your suggestion, I have made the following modifications to the wall part: t=0.01 dPos=Vector3.Ones*0.5*t mn,mx=Vector3(0,0,0),Vector3(0.05,0.05,0.05) walls=aabbWalls([(mn-dPos),(mx+dPos)],thickness=t,material='wall') wallIds=O.bodies.append(walls) The other parts are exactly the same as the code I provided before.Although there are no previous problems, new problems have emerged. In previous simulations, the axial strain when the failure condition is reached is very small (0.03 or less),but this time, when the axial strain reaches 0.15 or more, it still does not reach my failure condition (the axial stress does not decrease).. What should I do? Is the wall thickness set too large? Thanks for help. -- 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