[Yade-users] Parallelization tests with more than 8 cores
Hello, I would like to ask, if somebody did recently (or know about) some time testing of Yade parallelization on more than 8 cores? I have found link [1], would it be possible to create some kind of plot time (or speedup) against number of cores? Thanks and have a nice day Jan [1] https://www.yade-dem.org/wiki/Multicore_Performance ___ 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 #235136]: Parallelization tests with more than 8 cores
New question #235136 on Yade: https://answers.launchpad.net/yade/+question/235136 Hello, I would like to ask, if somebody did recently (or know about) some time testing of Yade parallelization on more than 8 cores? I have found link [1], would it be possible to create some kind of plot time (or speedup) against number of cores? Thanks and have a nice day Jan [1] https://www.yade-dem.org/wiki/Multicore_Performance PS: sorry for duplicate, I forgot to use question interface in previous yade-users mail.. -- You received this question notification because you are a member of yade-users, which 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 #235138]: how to define friction at the boundaries in triaxial test
New question #235138 on Yade: https://answers.launchpad.net/yade/+question/235138 hi , please how can one define the friction at the top and bottom of the specimen in a Triaxial test simulation of a granular material. the center ,extension and the wall mask of my model are as below. O.bodies.append(utils.geom.facetBox((0.035,0.035,0.035),(0.035,0.035,0.035),wallMask=15,dynamic=1)) O.bodies.append(utils.geom.facetBox((0.035,0.035,0.035),(0.035,0.035,0.035),wallMask=16,dynamic=0))#bottom of the sample, which is not allowed to move. because of that dynamic=0 [O means False] thanks in advance. best regards stephen -- You received this question notification because you are a member of yade-users, which 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 #235138]: how to define friction at the boundaries in triaxial test
Question #235138 on Yade changed: https://answers.launchpad.net/yade/+question/235138 Status: Open = Answered Christian Jakob proposed the following answer: Hi stephen, Welcome to yade community. In yade you have to define materials (e.g. one material for particles and one for boundaries): ### properties: density = 2650 shear_modulus = 2e5 poisson_ratio = 0.15 young_modulus = 2*shear_modulus*(1+poisson_ratio) ### define materials: SphereMat = FrictMat(young=young_modulus,poisson=poisson_ratio,density=density,frictionAngle=0.5) BoundaryMat = FrictMat(young=young_modulus,poisson=poisson_ratio,frictionAngle=1) O.materials.append(SphereMat) O.materials.append(BoundaryMat) Then you can use the materials in the generation process: O.bodies.append(sphere([0,0,0],radius=R, material = SphereMat, fixed=True)) O.bodies.append(geom.facetBox((0.035,0.035,0.035),(0.035,0.035,0.035),wallMask=15, material = BoundaryMat, dynamic=1)) Hope it helps, Christian -- You received this question notification because you are a member of yade-users, which 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 #235031]: Projected Area of clump
Question #235031 on Yade changed: https://answers.launchpad.net/yade/+question/235031 Status: Open = Answered Jan Stránský proposed the following answer: Hello Sergio, What particle shapes do you use for your clumps? Could you please explain more in detail what do you mean by projected area of several clump? I checked the utils.perpendicularArea function very quickly, but could not find out how it works.. b1,b2=utils.sphere((0,0,0),1), sphere((1,2,3),2) O.bodies.append((b1,b2)) perpendicularArea(0) # result = 30 perpendicularArea(1) # result = 24 perpendicularArea(2) # result = 20 Thanks for the answer and hopefully we can iterate until your question is answered :-) cheers Jan 2013/9/2 Sergio question235...@answers.launchpad.net New question #235031 on Yade: https://answers.launchpad.net/yade/+question/235031 Hello, I am a very new Yade user, and I'd like to post my problem as I know that there are some Yade experts here that could advice me ;). I need to calculate the projected area (respect to the axis) of several clumps of my calculation. But I really don't know if there is an instruction to do that, or how I would have to proceed. I'm looking for doing the same as 'yade.utils.perpendicularArea(axis)', not with spheres, but clumps. Thank you for your support. -- You received this question notification because you are a member of yade-users, which 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 -- You received this question notification because you are a member of yade-users, which 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 #235141]: Simulate the behaviour of a Bingham material (fresh concrete)
New question #235141 on Yade: https://answers.launchpad.net/yade/+question/235141 Hello all, I'm try to simulate the behavior of fresh concrete, this kind of material presents a Bingham behavior (with yield stress to determine the motion). The idea is to use something like the papper Simulation of fresh concrete flow using Discrete Element Method (DEM): theory and applications (Viktor Mechtcherine•Annika Gram•Knut Krenzer•Jörg-Henry Schwabe•Sergiy hyshko•Nicolas Roussel) - DOI 10.1617/s11527-013-0084-7 http://link.springer.com/content/pdf/10.1617%2Fs11527-013-0084-7.pdf#CR14 So the problem is to use a constitutive law like this one. Is ther something similar in Yade where I can begin? Thanks! Ricardo Pieralisi -- You received this question notification because you are a member of yade-users, which 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 #235031]: Projected Area of clump
Question #235031 on Yade changed: https://answers.launchpad.net/yade/+question/235031 Sergio posted a new comment: Hello Jan, Thank you very much for your feedback. My clumps are made from overlapping spheres of different radius (15-30 mm). From what I've understand reading the yade's doc, the 'perpendicularArea' has meaning only with spheres bodies. So, If I made a clump of several overlapping spheres, I understand that the utility is not valid. Appart from that, how can I check how the 'utils.perpendicularArea' works? I have to say that I work with yade-daily, but I downloaded the 'trunk' folder. Am I rigth looking for the the utility in this folder? As you can say, I'm totally lost. But I'll keep trying. Thank you for your time and help ; ) Sergio -- You received this question notification because you are a member of yade-users, which 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 #235138]: how to define friction at the boundaries in triaxial test
Question #235138 on Yade changed:
https://answers.launchpad.net/yade/+question/235138
Status: Answered = Open
stephen is still having a problem:
hi Jakob,
thanks for your prompt response and great contribution.
however, i made some changes and try to run my script again, but the result i
was geting is not what i was expected. here is the detailed copy of the script
for your suggestions and corrections.
thanks in advanced.
#file=open(data1.txt,w)
#file.close()
#stra=open(strain.txt,w)
#stra.write(str(0)+\n) #because we need to add the zero strain. #look at
later, may be we should also add zero force???
#stra.close()
devi=open(deviatoric.txt,w) #deviatoric stresses file
#devi.write(str(0)+\n) #adding zero deviatoric stress
devi.close()
strain=open(strain.txt,w)
#strain.write(str(0)+\n) #adding zero strain
strain.close()
sigma11=open(sigma11.txt,w) # Axial stress file
#strain.write(str(0)+\n) #adding zero strain
sigma11.close()
shear11=open(shear11.txt,w)# shear stress file. shear stress =
((sigma1-sigma2)/2)
#strain.write(str(0)+\n) #adding zero strain
shear11.close()
epsellion=(epsellion.txt,w)
#epsellion.close()
file=open(sigma1.txt,w)
file1=open(sigma2.txt,w)
file.close()
file1.close()
hop=open(leng.txt,w)
#strain.write(str(0)+\n) #adding zero strain
hop.close()
from yade import pack, plot
young=1e6 #young modulus of the spheres [Pa] (1e6 kPa=1e9 Pa)
pr=0.25 #poisson's ratios of the spheres
fa=28.6 #friction angle of the spheres [rad] 0.2 radian=11.46 degree, 0.5
radian= 28.6 degree, 0.4 rad=22.9 degreee, 0.8 rad=45.8 degree
fb=28.6 #friction at the boundaries in degree
number=1500
##create materials for spheres and plates
idSand=O.materials.append(FrictMat(young=young,poisson=pr,frictionAngle=radians(fa),label='idSand'))
#defining sphere properties
BoundaryMat = FrictMat(young=young,poisson=pr,frictionAngle=fb)
O.materials.append(BoundaryMat)
O.bodies.append(geom.facetBox((0.035,0.035,0.035),(0.035,0.035,0.035),wallMask=15,
material = BoundaryMat, dynamic=1))
O.bodies.append(geom.facetBox((0.035,0.035,0.035),(0.035,0.035,0.035),wallMask=15,
dynamic=1))
O.bodies.append(utils.geom.facetBox((0.035,0.035,0.035),(0.035,0.035,0.035),wallMask=16,dynamic=0))#bottom
of the sample, which is not allowed to move. because of that dynamic=0
[O means False]
O.bodies.append(utils.geom.facetBox((0.035,0.035,0.105),(0.035,0.035,0.035),wallMask=15))#
in order to control the particles!!! after the end of the stabilization, we
will remove it
psdSizes,psdCumm=[0.00525,0.00603,0.0105,0.01206,0.02099],[0,0.1,0.5,0.6,1.0]
#d50=10.5mm, cu=2
mincorner,maxcorner=Vector3(0.001,0.001,0.001),Vector3(0.07,0.07,0.14) #min
and max coordinates of the wall
sp=pack.SpherePack() #burda pack modulunun Spherepack modulunu cagiriyo, bi
altinda da o cagirdinin,piramit gibi lan iste!
sp.makeCloud(mincorner,maxcorner,psdSizes=psdSizes,num=number,psdCumm=psdCumm,distributeMass=1)
#sp.makeCloud(mincorner,maxcorner,rMean=0.005,rRelFuzz=0.3)
Gl1_Sphere.stripes=True #showing the spheres in two colors
sp.toSimulation()#color=(0,0,1))
unb=utils.unbalancedForce()
O.engines=[
ForceResetter(),
InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb()]),
InteractionLoop(
[Ig2_Sphere_Sphere_L3Geom(),Ig2_Facet_Sphere_L3Geom()],
[Ip2_FrictMat_FrictMat_FrictPhys()],
[Law2_L3Geom_FrictPhys_ElPerfPl()]
),
GravityEngine(gravity=(0,0,-9.81)),
NewtonIntegrator(damping=0.7), #after the first phase we are going to change
it to 0.5
PyRunner(command='checkUnbalanced()',realPeriod=2,label='checker'),
]
O.dt=0.3*utils.PWaveTimeStep()
D E P O S I T I O N
print '### Deposion is Started ###'
def checkUnbalanced():
print 'iteration number:',O.iter,'unbalanced force:', unbalancedForce()
if O.iter25000: return
# the rest will be run only if unbalanced is .1 (stabilized packing)
if utils.unbalancedForce()0.1: return
if utils.unbalancedForce()0.1:
O.bodies.append(geom.facetBox((0.035,0.035,0.035),(0.035,0.035,0.035),wallMask=15,dynamic=1))
checker.command='erasingwalls()'
print Depostion is Completed and the Top Part of the Membran is Added!
def erasingwalls():
print ### Additional Walls are Erased! ###
O.bodies.erase(10)
O.bodies.erase(11)
O.bodies.erase(12)
O.bodies.erase(13)
O.bodies.erase(14)
O.bodies.erase(15)
O.bodies.erase(16)
O.bodies.erase(17)
checker.command='load()'
print loading is starting!
aim=10 #PASCAL!!
vel=0.03 #applyied velocity
#0.003 calisii
C O N F I N I N G P R E S S U R E
def load():
O.materials[idSand].young=1e9
O.engines=O.engines+[PyRunner(command='writingDatas()',iterPeriod=10)]
print ('CONFINING PRESSURE')
print 'stress at the top
Re: [Yade-users] [Question #235138]: how to define friction at the boundaries in triaxial test
Question #235138 on Yade changed: https://answers.launchpad.net/yade/+question/235138 Status: Open = Answered Christian Jakob proposed the following answer: in your script i see some inconsistencies, that can lead to unexpected results: # 1. problem: idSand=O.materials.append(FrictMat(young=young,poisson=pr,frictionAngle=radians(fa),label='idSand')) #frictionAngle=fa !!! # 2. problem: O.bodies.append(geom.facetBox((0.035,0.035,0.035),(0.035,0.035,0.035),wallMask=15, material = BoundaryMat, dynamic=1)) #this is ok O.bodies.append(geom.facetBox((0.035,0.035,0.035),(0.035,0.035,0.035),wallMask=15, dynamic=1)) #missing material = BoundaryMat !!! O.bodies.append(utils.geom.facetBox((0.035,0.035,0.035),(0.035,0.035,0.035),wallMask=16,dynamic=0)) #missing material = BoundaryMat !!! O.bodies.append(utils.geom.facetBox((0.035,0.035,0.105),(0.035,0.035,0.035),wallMask=15)) #missing material = BoundaryMat # 3.problem: sp=pack.SpherePack() sp.makeCloud(mincorner,maxcorner,psdSizes=psdSizes,num=number,psdCumm=psdCumm,distributeMass=1) sp.toSimulation() you did not set material for spheres. but i do not know how to do this with sp.toSimulation() i prefer to use: O.bodies.append([sphere(c,r,material=idSand) for c,r in sp]) #instead of sp.toSimulation() regards, christian -- You received this question notification because you are a member of yade-users, which 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 #235031]: Projected Area of clump
Question #235031 on Yade changed: https://answers.launchpad.net/yade/+question/235031 Sergio posted a new comment: Hello Jan, Thanks again for your quick response. To clarify what I mean with projected area, I point to a representative image I found in the internet: http://www.real-world-physics- problems.com/images/physics_skydiving_8.png I hope that the figure helps to understand what I mean : ) I will continue to understand a little bit more the function that Yade uses now, beginning by founding it in the trunk version. Thank you Sergio -- You received this question notification because you are a member of yade-users, which 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 #235031]: Projected Area of clump
Question #235031 on Yade changed: https://answers.launchpad.net/yade/+question/235031 Sergio posted a new comment: Thank you again, Jan ; ) With all the clues that you have given to me, I will try to face the problem for myself, but I'm sure that I will continue asking questions referring to it. you can estimate the projected area by defining a grid and checking if the point of the grid is contained by your clump or not. Along what axis you want to make the projection? x/y/z or some generic axis? I think you are considering to do something like is done in mass inertia approax in clumps, right? It could be an option... I will think about it. ; ) I will share my progress referring to the problem. Thanks Sergio -- You received this question notification because you are a member of yade-users, which 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 #225776]: Definition of Yade's frictionAngle for a single material
Question #225776 on Yade changed: https://answers.launchpad.net/yade/+question/225776 Status: Answered = Open Alexander Eulitz [Eugen] is still having a problem: I'd like to reopen this thread. I'm trying to implement a matchmaker for the friction angle. How can I do so? Is the match maker kind of a further engine in the simulation? Thanks, Alex -- You received this question notification because you are a member of yade-users, which 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 #235031]: Projected Area of clump
Question #235031 on Yade changed: https://answers.launchpad.net/yade/+question/235031 Jan Stránský proposed the following answer: Hello Sergio, http://www.real-world-physics- problems.com/images/physics_skydiving_8.png perfect, exactly what I meant :-) I hope that the figure helps to understand what I mean : ) I will continue to understand a little bit more the function that Yade uses now, beginning by founding it in the trunk version. you can find the function in yade/py/utils.py file. It uses aabbExtrema function inside, which you can find in yade/py/_utils.cpp file. Now back to your problem. One thing came to my mind, you can estimate the projected area by defining a grid and checking if the point of the grid is contained by your clump or not. Along what axis you want to make the projection? x/y/z or some generic axis? Would you like to try to implement it, or should we help you by implementing such function in Yade? cheers Jan -- You received this question notification because you are a member of yade-users, which 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
