date:20211006

Question #698984 on Yade changed:
https://answers.launchpad.net/yade/+question/698984

Status: Open => Needs information

Jan Stránský requested more information:
Hello,

please read [1] carefully and follow the instructions.

> I know the code to keep the particles stationary, but I don't know how
to use these codes.

then why not to mention this? "How do I keep ... particles ...
motionless?" suggest rather the opposite.

what is "the code" and "these codes"?

If [2] is solved, I see no problem in my answer #1. It is essentially
the same, just using "b.state.blockedDOFs=..." instead of
"O.bodies.erase(b.id)"

Again, to get a good (and quick) answer, first provide a good question.
If you want others to make effort solving your problem, it is much 
better/motivating/... to show your effort (formulating a good question) first.

Cheers
Jan

[1] https://www.yade-dem.org/wiki/Howtoask
[2] https://answers.launchpad.net/yade/+question/698980

-- 
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 #698984]: How do I keep a circle of particles around a circle motionless?

Question #698984 on Yade changed:
https://answers.launchpad.net/yade/+question/698984

Status: Answered => Open

孙灿 is still having a problem:
I know the code to keep the particles stationary, but I don't know how
to use these codes. As in this example, keep the particles stationary on
a circle with a center of (0.05, 1, 1) and a radius of 1

-- 
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 #698985]: Plotting force, unit question.

Question #698985 on Yade changed:
https://answers.launchpad.net/yade/+question/698985

Status: Open => Solved

Szabó Bence confirmed that the question is solved:
Thank you for your advice!
I think i solved the problem. 
I exported the *.stl in meter after that i changed the dimensions of the box 
and the spheres. 
I tought that the dimensions i set is in mm but its in m, my fault. 

Thanks again!

Cheers,
Bence

-- 
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 #698985]: Plotting force, unit question.

Question #698985 on Yade changed:
https://answers.launchpad.net/yade/+question/698985

Status: Answered => Open

Szabó Bence is still having a problem:
Thanks, for your advice.
Just one more question:
O.bodies.append(geom.facetBox((109.14,200,940),(280,320,750),wallMask=55,fixed=True,material='02'))
In my code is this mean the box dimensions are in mm or m?

Thank you in advance!

Cheers,
Bence

-- 
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 #698985]: Plotting force, unit question.

2021-10-06 Thread Karol Brzezinski

Question #698985 on Yade changed:
https://answers.launchpad.net/yade/+question/698985

Status: Open => Answered

Karol Brzezinski proposed the following answer:
You haven't provided the *.stl file, and I cannot run your example.
However, I am quite sure that Yade doesn't exchange unit information
with Solidworks by your *.stl file. Just change units to meters in
Solidworks and export again if you are not convinced.

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 #698985]: Plotting force, unit question.

Question #698985 on Yade changed:
https://answers.launchpad.net/yade/+question/698985

Status: Answered => Open

Szabó Bence is still having a problem:
Hello, 
Thanks for your reply. 
I exactly know the parameters of the tool i designed it in solidworks then 
exported to stl then imported into the simulation (width 220mm, lenght 175, 
height 500mm ). I compared the sizes of the tool and the box so i set these 
parameters and im only guessing the parameters are in mm. 
I tought the sizes of the box are 560mm x 640mm x 1500 mm, and the spheres 
similarly mm. 
Is it true?

Thanks for your helping!

Cheers, 
Bence

-- 
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 #698985]: Plotting force, unit question.

2021-10-06 Thread Karol Brzezinski

Question #698985 on Yade changed:
https://answers.launchpad.net/yade/+question/698985

Status: Open => Answered

Karol Brzezinski proposed the following answer:
Hi,

I think that you should be consistent with units and use meters for the
dimensions of the box, spheres, tool. The forces acting on the tool
originates mostly from gravity. Your dimension is three times too big,
thus mass is overestimated nine times. Mass influences the vertical
load, load influences frictional force, etc. In the end, you get the
results too big 1e9 times.

Best wishes,
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 #698985]: Plotting force, unit question.

New question #698985 on Yade:
https://answers.launchpad.net/yade/+question/698985

Hello Yade friends!

I'm trying to simulate a tillage tool - soil interaction. I plotted the draft 
force, but on the diagram i see 1e11 scale at the force number, but it should 
be around 1e3 in Newton. 
The density [kg/m3], the young modulus [Pa], the dimensions of the box, 
spheres, tool are in [mm].
Am I wrong something with units in my code?

Thanks for your helping!

My code:

from yade import pack, plot, qt

talaj=O.materials.append(FrictMat(density=2600,young=6e7,poisson=0.25,frictionAngle=radians(30),label="01"))
acel=O.materials.append(FrictMat(density=7850,young=2e8,poisson=0.3,frictionAngle=radians(30),label="02"))

O.bodies.append(geom.facetBox((109.14,200,940),(280,320,750),wallMask=55,fixed=True,material='02'))

sp=pack.SpherePack()
c4=pack.SpherePack([((0,0,0),5),((9.8,0,0),4.75),((5.1,7.4,0),4),((5.1,-7.7,0),4.25),((5.1,0,8),4.5)])

sp.makeClumpCloud((-170.86,-120,190),(389.14,740,1690),
  [c4],

  num=-1)
sp.toSimulation(color=(0.45,0.28,0.0),material='01')

from yade import ymport

kes=O.bodies.append(ymport.stl('nkapa.STL',wire=False,material='02',color=(0.5,0.5,0.5)))

O.engines=[
ForceResetter(),
InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb()]),
InteractionLoop(

[Ig2_Sphere_Sphere_ScGeom(),Ig2_Facet_Sphere_ScGeom()],
[Ip2_FrictMat_FrictMat_FrictPhys()],
[Law2_ScGeom_FrictPhys_CundallStrack()]
),
TranslationEngine(translationAxis=[0,0,1],velocity=0,ids=kes),
NewtonIntegrator(gravity=(0,-9.81,0),damping=0.2),

PyRunner(command='checkUnbalanced()',realPeriod=0.2),

PyRunner(command='addPlotData()',iterPeriod=100)

]

O.dt=.5*PWaveTimeStep()

def checkUnbalanced():
if unbalancedForce()<0.1:
O.engines[3].velocity=4

def addPlotData():

Fz=sum(O.forces.f(id)[2]for id in kes)
plot.addData(Fz=Fz*-1,t=O.time,unbalanced=unbalancedForce(),i=O.iter)

plot.plots={'t':('Fz')}

plot.plot()

O.saveTmp()


-- 
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 #698984]: How do I keep a circle of particles around a circle motionless?

Question #698984 on Yade changed:
https://answers.launchpad.net/yade/+question/698984

Status: Open => Answered

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

please read [1].

Next time, try to make the question standalone. If it is related to another 
question, please at least reference it, e.g. here it is [2].
The questions and answers are there also for others, who may search them in 
future and find their answers.
As such, the question "How do I keep a circle of particles around a circle 
motionless?" makes no sense.

To keep one particle "b" motionless, block all its degrees of freedom:
b.state.blockedDOFs = "xyzXYZ" # [3,4]

To select relevant particles, use similar code as proposed in [2]

Cheers
Jan

[1] https://www.yade-dem.org/wiki/Howtoask
[2] https://answers.launchpad.net/yade/+question/698980
[3] https://yade-dem.org/doc/yade.wrapper.html#yade.wrapper.State.blockedDOFs
[4] https://yade-dem.org/doc/user.html#motion-constraints

-- 
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 #698984]: How do I keep a circle of particles around a circle motionless?

Question #698984 on Yade changed:
https://answers.launchpad.net/yade/+question/698984

Summary changed to:
How do I keep a circle of particles around a circle motionless?

-- 
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 #698980]: How to remove particles from a packing?

Question #698980 on Yade changed:
https://answers.launchpad.net/yade/+question/698980

Status: Answered => Solved

孙灿 confirmed that the question is solved:
Thank you, Jane

-- 
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 #698980]: How to remove particles from a packing?

Question #698980 on Yade changed:
https://answers.launchpad.net/yade/+question/698980

孙灿 confirmed that the question is solved:
Thanks Jan Stránský, 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

[Yade-users] [Question #698984]: 怎么样让颗粒不动

New question #698984 on Yade:
https://answers.launchpad.net/yade/+question/698984

How do I keep a circle of particles around a circle motionless?

-- 
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 #698980]: How to remove particles from a packing?

Question #698980 on Yade changed:
https://answers.launchpad.net/yade/+question/698980

Status: Open => Answered

Jan Stránský proposed the following answer:
Thanks for clarification.

To remove particles / use only some particles, you have several options.

One is to extract particles and use only those you need in reset
simulation or a new simulation.

Another option is to use O.bodies.erase [1].
###
circleCenter = Vector3(1,2,3)
circleRadius = 123
for b in O.bodies:
d = (b.state.pos - circleCenter).norm() # distance of circleCenter and 
center of "b"
if d < circleRadius:
O.bodies.erase(b.id)   
###

Cheers
Jan

[1] https://yade-
dem.org/doc/yade.wrapper.html#yade.wrapper.BodyContainer.erase

-- 
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 #698980]: How to remove particles from a packing?

Question #698980 on Yade changed:
https://answers.launchpad.net/yade/+question/698980

Status: Needs information => Open

孙灿 gave more information on the question:
First of all, my simulation is that all particles are in a plane, in the
Y-Z plane, that is, a layer structure. I want to remove a circle in the
middle of the plane and make the particles in the middle of the circle
disappear.

-- 
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 #698980]: How to remove particles from a packing?

Question #698980 on Yade changed:
https://answers.launchpad.net/yade/+question/698980

Status: Open => Needs information

Jan Stránský requested more information:
What is "a circle"? How should it be removed (what should be the
result)?

To get a good answer, please provide a good question..

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 #698980]: How to remove particles from a packing?

Question #698980 on Yade changed:
https://answers.launchpad.net/yade/+question/698980

Status: Needs information => Open

孙灿 gave more information on the question:
Remove a circle from the stable structure

-- 
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 #698963]: number of cores to have the max performance with Yade

2021-10-06 Thread ozge

Question #698963 on Yade changed:
https://answers.launchpad.net/yade/+question/698963

Status: Answered => Solved

ozge confirmed that the question is solved:
Thanks Jérôme Duriez, 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 #698948]: Stability of advection modeling by ThermalEngine

2021-10-06 Thread Robert Caulk

Question #698948 on Yade changed:
https://answers.launchpad.net/yade/+question/698948

Status: Open => Answered

Robert Caulk proposed the following answer:
Hey Zoheir,

Thanks for the good feedback. Can you precise to me the example script
that you are basing this on? We have 4 example scripts in the
repository. I guess you are using flowScenario.py since you had
b.dynamic=False before as well as a pressure boundary condition?

>> If so, could you give me some hint on how to avoid this issue.

We've discussed in the past how to estimate the stable timestep
dynamically [1]. Did you explore that already for your case? I started
adding that feature to ThermalEngine [2][3], but it is still
experimental.

Let's try to debug this through some kind of trial and error, I think it
will identify the problem most rapidly. We can consider the factors
discussed in [1] and hold various components constant. Since you know it
is in the fluid conduction, then you just need to consider the
parameters associated with fluid conduction such as
thermal.fluidConductionAreaFactor[4] and thermal.fluidK. As we discuss
in [1], increasing the diffusion coefficient [5] is going to decrease
the maximum allowable timestep. You have pythonic control over all 3
parameters in that diffusion equation which means you should be able to
determine the exact cause of the instability through trial and error. I
would first start by holding them all constant at extraordinary values
i.e. force that diffusion constant to be extremely low to test if the
system is perfectly stable. Then open up each parameter back to your
desired values indivdiually. BTW to be strict about the test you need to
set flow.meshUpdateInterval=-1 and b.dynamic=False - since those two
components adjust geometric quantities used for the diffusion
coefficient and which evolve with time during the simulation.

Another easy test would be to not run compaction (basically just holding
bodies static) and see if the simulation remains stable.

Thanks for testing this stuff, let me know what you come up with.

[1]https://answers.launchpad.net/yade/+question/695542
[2]https://gitlab.com/yade-dev/trunk/-/blob/master/pkg/pfv/Thermal.cpp#L177
[3]https://yade-dem.org/doc/yade.wrapper.html#yade.wrapper.ThermalEngine.tsSafetyFactor
[4]https://yade-dem.org/doc/yade.wrapper.html#yade.wrapper.ThermalEngine.fluidConductionAreaFactor
[5]https://gitlab.com/yade-dev/trunk/-/blob/master/pkg/pfv/Thermal.cpp#L544

Thanks,

Robert

--
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 #698948]: Stability of advection modeling by ThermalEngine

2021-10-06 Thread Zoheir Khademian

Question #698948 on Yade changed:
https://answers.launchpad.net/yade/+question/698948

Status: Answered => Open

Zoheir Khademian is still having a problem:
Hi Robert and thanks for the response.

>>What changes have you made in comparison to the example script? 
1- Turning internal compaction off in the triaxial engine and compacted the 
sample by inwards-moving walls.
2- Reducing the particle density from 2600e10 to 2600 to run the triaixal 
engine. 
3- Setting bodies dynamic false after compaction so there is no movement.
4- Turning off the thermo-mechanical analyses in thermal engine.
5- flow.bndCondIsPressure=[0,0,0,0,1,1] ; flow.bndCondValue=[0,0,0,0,0,10]
6- changing rRelFuzz from 0.333 to 0.33 #Small changes in the packing can make 
the model unstable like changing (mn,mx) a bit.

After making these changes, when I run, the model maximum body and fluid
temp exceed the boundary temp after 2 s.

>> I would apply more traditional boundary conditions a bit to confirm that you 
>> can reach stable states for your setup:
I changed the boundary conditions to the traditional ones as you suggested but 
it did not make any difference in stability.


>>You can take some manual control with thermal.minimumFluidCondDist[1].
I tried assigning different values (0.0015 to 0.015) to minimumFluidCondDist[1] 
but it did not make any difference.


>>However, you should also be able to figure out if that is the issue by simply 
>>turning on and off thermal.fluidConduction.
I turned off thermal.fluidConduction  and it solved the instability. I checked 
the centers of neighboring cells and some of them had pretty close X or Y or Z 
coordinates (less than 1e-5). Is this the issue? If so, could you give me some 
hint on how to avoid this issue. 


[1]https://yade-dem.org/doc/yade.wrapper.html#yade.wrapper.ThermalEngine.minimumFluidCondDist

-- 
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 #698940]: Why particles move along X axis?

2021-10-06 Thread Robert Caulk

Question #698940 on Yade changed:
https://answers.launchpad.net/yade/+question/698940

Summary changed to:
Why particles move along X axis?

-- 
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 #698980]: How to remove particles from a packing?

2021-10-06 Thread Robert Caulk

Question #698980 on Yade changed:
https://answers.launchpad.net/yade/+question/698980

Summary changed to:
How to remove particles from a packing?

-- 
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 #698980]: 结构