Re: [Yade-users] [Question #696481]: Geometry of hyperboloid

2021-04-09 Thread Faqih Maarif
Question #696481 on Yade changed:
https://answers.launchpad.net/yade/+question/696481

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


Re: [Yade-users] [Question #696481]: Geometry of hyperboloid

2021-04-09 Thread Faqih Maarif
Question #696481 on Yade changed:
https://answers.launchpad.net/yade/+question/696481

Status: Answered => Solved

Faqih Maarif confirmed that the question is solved:
Dear Jan,

But please provide information, what is the problem.
I looking for pack.Hyperboloid but I cannot find it. Thank you for the solution 
in [4].
Furthermore, I am following [4] for new geometry of Hyperboloid in my model, it 
can run properly.

it is just your [1] with #L81 at the end.
Thank you, I got it

Regards,
Faqih

-- 
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 #696481]: Geometry of hyperboloid

2021-04-09 Thread Faqih Maarif
Question #696481 on Yade changed:
https://answers.launchpad.net/yade/+question/696481

Faqih Maarif posted a new comment:
Dear Jan,

But please provide information, what is the problem.
I looking for pack.Hyperboloid but I cannot find it.  Thank you for the 
solution in [4].
Furthermore, I am following [4] for new geometry of Hyperboloid in my model, it 
can run properly.

it is just your [1] with #L81 at the end.
Thank you, I got it
 

Regards,
Faqih

-- 
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 #696481]: Geometry of hyperboloid

2021-04-09 Thread Jan Stránský
Question #696481 on Yade changed:
https://answers.launchpad.net/yade/+question/696481

Status: Needs information => Answered

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

thanks for clarification, now it is clear what you want.
But please provide information, what is the problem.

Everything what you need is specify suitable numbers to pack.inHyperboloid [4] 
- centerBottom, centerTop, radius, skirt.
E.g.:
diameter = 0.15 # 15 cm
diameterSkirt = 0.10 # 10 cm, or other value
length = 0.30 # 30 cm
pack.inHyperboloid((0,0,-.5*length),(0,0,.5*length),.5*diameter,.5*diameterSkirt)

> but I cannot open the [3]

it is just your [1] with #L81 at the end.
Note that yade main repository is now on gitlab, so [1] and [3] should be as 
below.

cheers
Jan

[1] https://gitlab.com/yade-dev/trunk/-/blob/master/examples/concrete/uniax.py
[3] 
https://gitlab.com/yade-dev/trunk/-/blob/master/examples/concrete/uniax.py#L81

-- 
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 #696481]: Geometry of hyperboloid

2021-04-09 Thread Faqih Maarif
Question #696481 on Yade changed:
https://answers.launchpad.net/yade/+question/696481

Faqih Maarif posted a new comment:
Dear Jan,
Thank you for your quick response.

Sorry for not clear on the questions.
I would like to create the geometry of hyperboloids in a uniaxial test of 
concrete. Based on the example [1] the dimensions of the model of 10cm (in 
diameter) x 20cm (in length), I want to change the new geometry of Hyperboloid 
to 15cm (in diameter) x 30cm (in length).

**Guessing** that links to pack.inHyperboloid ([3] in the script, [4] 
documentation) could be helpful
Thank you for the script, but I cannot open the [3]

Regards,
Faqih

[1] https://gitlab.com/yade-
dev/trunk/-/blob/master/examples/concrete/uniax.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 #696481]: Geometry of hyperboloid

2021-04-09 Thread Faqih Maarif
Question #696481 on Yade changed:
https://answers.launchpad.net/yade/+question/696481

Faqih Maarif posted a new comment:
Dear Jan,
Thank you for your quick response.

Sorry for not clear on the questions.
I would like to create the geometry of hyperboloids in a uniaxial test of 
concrete. Based on the example [1] the dimensions of the model of 10cmx15cm, I 
want to change the new geometry of Hyperboloid be 15cmx30cm. 

**Guessing** that links to pack.inHyperboloid ([3] in the script, [4] 
documentation) could be helpful
Thank you for the script, but I cannot open the [3]

Regards,
Faqih


[1] https://gitlab.com/yade-dev/trunk/-/blob/master/examples/concrete/uniax.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 #696376]: Contact between sphere-gridconnection not detected in a pfacet edge

2021-04-09 Thread Bruno Chareyre
Question #696376 on Yade changed:
https://answers.launchpad.net/yade/+question/696376

Status: Open => Answered

Bruno Chareyre proposed the following answer:
Hi,
Additional bugs could be discussed directly hre: 
https://gitlab.com/yade-dev/trunk/-/issues/210

Bruno

-- 
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 #696481]: Geometry of hyperboloid

2021-04-09 Thread Jan Stránský
Question #696481 on Yade changed:
https://answers.launchpad.net/yade/+question/696481

Status: Open => Needs information

Jan Stránský requested more information:
Hello,

> Could anyone help me with the problem?

firstly, please describe clearly the problem [2].
This way, the question is too confusing and ambiguous.

> I will create the diameter from 10cm to 15cm?

like here. 
Is this the problem how to create such hyperboloid? Why is then there the 
question mark?
Or is it a question if the dimensions are OK?
Or ... ?
What does "create" mean?
What does "the diameter" mean?
What does "from ... to" mean? One fixed value from that interval? Parameter 
study? ... ?
... ?

**Guessing** that links to pack.inHyperboloid ([3] in the script, [4]
documentation) could be helpful

cheers
Jan

[2] https://www.yade-dem.org/wiki/Howtoask
[3] https://github.com/yade/trunk/blob/master/examples/concrete/uniax.py#L81
[4] 
https://yade-dem.org/doc/yade.pack.html#yade._packPredicates.inHyperboloid.__init__

-- 
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 #696481]: Geometry of hyperboloid

2021-04-09 Thread Faqih Maarif
New question #696481 on Yade:
https://answers.launchpad.net/yade/+question/696481

Dear All,
I would like to model the geometry of hyperboloids for compression stress of 
concrete [1]. I have a question about these sentences in the code:
“make geom; the dimensions are hard-coded here; could be in param table if 
desired”. In the model [1], z-oriented hyperboloid, length 20cm, diameter 10cm, 
skirt 8cm. In my case, I will create the diameter from 10cm to 15cm? 

Could anyone help me with the problem?
Thank you very much for your attention.

Regards,
Faqih

[1] https://github.com/yade/trunk/blob/master/examples/concrete/uniax.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


[Yade-users] [Question #696480]: PFacet model - contact data mining

2021-04-09 Thread Paul Pircher
New question #696480 on Yade:
https://answers.launchpad.net/yade/+question/696480

Hello there!
Recently I've been working a lot with the PFacet approach to obtain an elastic 
model for bigger but simple bodies in DEM simulations. 
I managed to obtain correct mechanical behaviour of my elastic body. Now I am 
researching the contact with external particles.

I made up a simple example of a PFacet cuboid that is penetrated by a sphere. 
The goal is to do that at every different position on the cuboid (edge, corner, 
intersection...) to verify that at all positions the same contact behaviour 
(penetration depth, force etc.) appears. I do this in order to get a better 
understanding of how exactly the model works (coding wise) as I want to change 
the contact model to an adapted hertz-mindlin model afterwards. 

When I perform this penetration test e.g. within the area of a PFacet and on 
the intersection of 2 PFacets (so directly on a cylinder connection) I can see 
that the behaviour is according to the Cundall Strack law, which was expected, 
and that I only have 1 interaction on the sphere. 
However, if I perform this test on the very edge of the cuboid, like the 
following script shows, I get unexpected results. 

Mostly there is no interaction with the sphere, but according to the position 
and velocity of the sphere it behaves correctly and identical to other sphere 
positions (as mentioned before). 

If I go through the simulation step by step, at some timesteps 3 interactions 
are found. Those 3 interactions are between the sphere and 3 different PFacets, 
one of them should be not even close to the sphere. If 3 interactions are 
found, I would have expected that there are 2 interactions with PFacets and 1 
with the cylindrical connection. 
Still, this isn't in accordance with the other tests where the sphere is at a 
different position and I am still not able to extract the wanted data. 

So my question is: How can I correctly mine the contact parameters like 
penetration depth, force on the sphere and number of interactions for my 
penetration tests?
Am I using the wrong functions or am I looking at the wrong spot?

Thanks in advance!

P.S.: Since it will be needed in the future: If there is an "easy" way (only 
changing python code, not c++ source code) to change the contact law for the 
PFacet - external particles interactions, I would be happy to receive hints. 

Working with YADE 2020.01a on Ubuntu 20.04.2 LTS

##

from yade import plot, qt
import os, sys, time
from yade.gridpfacet import *
import numpy as np, math

""" All manual entries for the simulation """
O.dt = 1e-8 # timestep

node_mat_name = 'gridNodeMat'
node_r = 10 * 1e-3 # radius for the nodes, cylinders and pfacets
node_E = 52.*1e9 # (concrete)
node_poisson = 0.167
node_rho = 2750
node_phi = 35. 
node_normal_coh = 3e1000 # high values to not lose node interactions. total 
elasticity, no plasticity
node_shear_coh = 3e1000

pfacet_mat_name = 'pFacetMat' # this material has no influence on the beam 
stiffness, only for contacts
pfacet_E = node_E
pfacet_poisson = node_poisson
pfacet_rho = node_rho
pfacet_phi = node_phi

color = [0, 0, 1]   # color for certain bodies (PFacets)

sphere_mat_name = 'sphereMat'   
sphere_poisson = 0.2
sphere_phi = 24.   
sphere_E = 30.*1e6  
sphere_rho = 2660 

sphere_r = 40 * 1e-3 # radius of the penetrationg sphere
sphere_v = -0.1*100 # velocity of the penetrating sphere

dim_x = 0.5 # length in x-direction of the cuboid
dim_y = 0.6 # width of the cuboid
dim_z = 0.7 # height of the cuboid

# initial position of the penetrating sphere
sphere_x = (1/2) * dim_x 
sphere_y = 0
sphere_z = dim_z + sphere_r + node_r

engine_gravity = (0,0,0)# gravity acceleration
engine_damping = 0.0# general damping factor

plot_label = "plotter"
plotData_period = 50# iterations

 The simulation engine """

O.engines=[
ForceResetter(),
InsertionSortCollider([Bo1_Sphere_Aabb(), Bo1_Node_Aabb(), 
Bo1_GridConnection_Aabb(), Bo1_PFacet_Aabb(), ]),
InteractionLoop(
[Ig2_GridNode_GridNode_GridNodeGeom6D(), 
Ig2_GridConnection_GridConnection_GridCoGridCoGeom(),
Ig2_Sphere_PFacet_ScGridCoGeom(), Ig2_PFacet_PFacet_ScGeom(),
Ig2_Sphere_GridConnection_ScGridCoGeom(), Ig2_Sphere_Sphere_ScGeom(), 
],   

[Ip2_CohFrictMat_CohFrictMat_CohFrictPhys(setCohesionNow=True,setCohesionOnNewContacts=True),
Ip2_FrictMat_FrictMat_FrictPhys() ],
[Law2_ScGeom6D_CohFrictPhys_CohesionMoment(), 
Law2_ScGeom_FrictPhys_CundallStrack(), 
Law2_ScGridCoGeom_FrictPhys_CundallStrack(), 
Law2_GridCoGridCoGeom_FrictPhys_CundallStrack(), ]
),
NewtonIntegrator(  gravity=engine_gravity,  damping=engine_damping),
PyRunner( command = 'addPlotData()', iterPeriod = plotData_period),
]

""" Materials """
O.materials.append( 
CohFrictMat( young = node_E, poisson = node_poisson, density = node_rho,
frictionAngle = r