Re: [Yade-users] [Question #696481]: Geometry of hyperboloid
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
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
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
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
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
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
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
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
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
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