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

2022-04-27 Thread xuanshenyu
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

2022-04-27 Thread xuanshenyu
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

2022-04-27 Thread Jan Stránský
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

2022-04-27 Thread Karol Brzezinski
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

2022-04-27 Thread Mithushan Soundaranathan
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

2022-04-27 Thread Janek Kozicki
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

2022-04-27 Thread xuanshenyu
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

2022-04-27 Thread Aoxi Zhang
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

2022-04-27 Thread Karol Brzezinski
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

2022-04-27 Thread Ziyu Wang
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