yade-users team mailing list archive

[Rahul R <question701028@xxxxxxxxxxxxxxxxxxxxx>]

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

    Status: Answered => Open

Rahul R is still having a problem:
Hi Jan,
Thank you for your reply.
I tried scaling the packing, but that doesn't seem to eliminate the issue.
As for the first question, I created the specimen by creating a cylindrical
pack and depositing it in a cylindrical facet box of the desired diameter
(5 units), and doing compression and extension cycles manually till the
desired porosity range and height of the specimen was reached. I've shared
the coordinates file in the mail. Here is the code I used to create the
sample:

from yade import pack

#From PSD 1 of PSD curves. Density used is 2530kg/m3
psdSizes,psdCumm=[0.0694, 0.07508, 0.08895, 0.11967, 0.16797, 0.20143,
0.26611, 0.30037, 0.34524, 0.39682, 0.414, 0.44251, 0.46447, 0.49646,
0.52426, 0.60258, 0.70105, 0.73584, 0.80093, 0.89857, 0.99597, 1.23104,
1.64618, 1.99808, 2.57658, 3.44547, 4.69182],[0.00039, 0.01361, 0.02507,
0.04094, 0.08854, 0.12644, 0.20049, 0.23928, 0.29833, 0.38472, 0.41469,
0.45788, 0.49931, 0.55573, 0.59892, 0.69853, 0.77434, 0.80078, 0.83604,
0.87924, 0.90833, 0.94358, 0.96915, 0.98325, 0.99207, 0.99647, 1.]

#Create a cylinder pack
sp=pack.SpherePack();
sp.makeCloud((25,25,55),(29,29,110),psdSizes=psdSizes,psdCumm=psdCumm,distributeMass=True,num=35000);
#This works with f=0.15.
cyl=pack.inCylinder((27,27,55),(27,27,110),radius=2) #Creating a predicate
to fill the pack.
sp1=pack.filterSpherePack(cyl,sp,True) #Create a cylindrical pack using the
predicate "cyl" and the psd sizes and %fines provided in "makeCloud".
sp1.toSimulation();

#Creating a cylindrical box.
O.bodies.append(geom.facetCylinder((27,27,81),radius= 2.5,height=
60,fixed=True));

#Creating and identifying interactions btw particles, walls and facets.
O.engines=[
ForceResetter(),
InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb(),Bo1_Wall_Aabb()]),
InteractionLoop(
[Ig2_Sphere_Sphere_ScGeom(),Ig2_Facet_Sphere_ScGeom(),Ig2_Wall_Sphere_ScGeom()],
[Ip2_FrictMat_FrictMat_FrictPhys()],
[Law2_ScGeom_FrictPhys_CundallStrack()]
              ),
NewtonIntegrator(gravity=(0,0,-9.81),damping=0.5),
#PyRunner(command='compaction()',iterPeriod=50)
]
O.dt=.5*PWaveTimeStep()

On Tue, 22 Mar 2022 at 16:41, Jan Stránský <
question701028@xxxxxxxxxxxxxxxxxxxxx> wrote:

> Your question #701028 on Yade changed:
> https://answers.launchpad.net/yade/+question/701028
>
>     Status: Open => Answered
>
> Jan Stránský proposed the following answer:
> Hello,
>
> > I'm new in Yade.
>
> welcome :-)
>
> > I found a code on github ...
> > http://bazaar.launchpad.net ...
>
> we are using GitLab currently as a main repository platform [1]
>
> > the cylindrical membrane ... is created away from the pack.
>
> >  width = 5,
> > height = 9.85
> > ...
> > sp= ymport.textExt(...)
> > spheres = O.bodies.append(sp)
>
> How did you create the packing? What is its dimensions and position?
> I guess it simply does not fit together with defined cylinder dimensions
> and position..
>
> Shift and/or resize the packing and/or the cylinder.
>
> Cheers
> Jan
>
> [1] https://gitlab.com/yade-dev/trunk
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
> https://answers.launchpad.net/yade/+question/701028/+confirm?answer_id=0
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https://answers.launchpad.net/yade/+question/701028
>
> You received this question notification because you asked the question.
>

-- 
You received this question notification because your team yade-users is
an answer contact for Yade.