Re: [Yade-users] [Question #695332]: Triaxial test with rolling resistance effect
Question #695332 on Yade changed: https://answers.launchpad.net/yade/+question/695332 Status: Answered => Open Alireza Sadeghi is still having a problem: Hello, Thank you very much for your response. I added some new comments in the Ip2 and Law2 sections. It seems to work now. Is it correct? O.engines=[ ForceResetter(), InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Box_Aabb(),Bo1_Facet_Aabb()]), InteractionLoop([Ig2_Sphere_Sphere_ScGeom6D(),Ig2_Box_Sphere_ScGeom6D(),Ig2_Facet_Sphere_ScGeom6D()], [Ip2_CohFrictMat_CohFrictMat_CohFrictPhys(setCohesionNow=True)], [Law2_ScGeom_FrictPhys_CundallStrack(),Law2_ScGeom6D_CohFrictPhys_CohesionMoment(useIncrementalForm=True, always_use_moment_law=True, label='cohesiveLaw')]), TriaxialStateRecorder(iterPeriod=10,file='WallStresses'), VTKRecorder(fileName='3d-vtk-',recorders=['all'],iterPeriod=100), NewtonIntegrator(damping=0.4,gravity=[0,0,-9.81]), PyRunner(iterPeriod=1,command='addPlotData()'), PyRunner(command='calm()',iterPeriod=10,label='calmEngine') ] -- 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 #695332]: Triaxial test with rolling resistance effect
New question #695332 on Yade:
https://answers.launchpad.net/yade/+question/695332
Hello,
I am working on the behavior of granular material under the triaxial test. I
want to find the effect of rolling resistance on the shear behavior of the
material. I simulated as the instruction, but I have no idea why the rolling
resistance does not affect the shear response. A very short version of my
script is here. I was wondering if you could give me a suggestion if I am doing
wrong in my code.
I really appreciate your help in advance.
Best Regards,
Alireza
from yade import utils, plot
from yade import pack, qt
from datetime import datetime
import math
import pylab
from numpy import *
#==
#= define the materials ===
#==
O.materials.append(CohFrictMat(normalCohesion= 1e20, shearCohesion= 1e20,
isCohesive= True, young=6.81e8, density=1377.5, poisson=0.3, frictionAngle=
0.31, fragile=False, alphaKr=1.5, momentRotationLaw=True, etaRoll=0.5,
label='aggregate-48'))
O.materials.append(CohFrictMat(normalCohesion= 1e20, shearCohesion= 1e20,
isCohesive= True, young=6.81e8, density=1532.2, poisson=0.3, frictionAngle=
0.31, fragile=False, alphaKr=1.5, momentRotationLaw=True, etaRoll=0.5,
label='aggregate-814'))
O.materials.append(CohFrictMat(normalCohesion= 1e20, shearCohesion= 1e20,
isCohesive= False, young=2e11, density=0, poisson=0.3, frictionAngle= 0,
fragile=False, alphaKr=0, momentRotationLaw=True, etaRoll=0,
label='wall-frictionless'))
#===
#=== define walls
#===
walls=aabbWalls([(0,0,0),(9e-2,9e-2,9e-2)],thickness=0.001,oversizeFactor=1.5,material='wall-frictionless')
wallIds=O.bodies.append(walls)
#===
#=== define packing
#===
nums=['t','t','t']
mats=['aggregate-48','aggregate-814','aggregate-814']
coke=((1.875e-3,500),(0.9475e-3,1074),(0.50125e-3,2090))
color=((0,0,1),(0,1,0),(1,0,0))
tolerance=[(0.4752e-3),(0.2055e-3),(0.1539e-3)]
for i in range(len(nums)):
nums[i]=pack.SpherePack()
nums[i].makeCloud((5e-3,5e-3,5e-3),(85e-3,85e-3,85e-3),rMean=coke[i][0],rRelFuzz=tolerance[i],num=coke[i][1])
O.bodies.append([utils.sphere(c,r,material=mats[i],color=color[i]) for c,r
in nums[i]])
#===
#=== define Engine =
#===
O.engines=[
ForceResetter(),
InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Box_Aabb(),Bo1_Facet_Aabb()]),
InteractionLoop([Ig2_Sphere_Sphere_ScGeom6D(),Ig2_Box_Sphere_ScGeom6D(),Ig2_Facet_Sphere_ScGeom6D()],
[Ip2_CohFrictMat_CohFrictMat_CohFrictPhys()],
[Law2_ScGeom6D_CohFrictPhys_CohesionMoment()]
),
TriaxialStateRecorder(iterPeriod=10,file='WallStresses'),
VTKRecorder(fileName='3d-vtk-',recorders=['all'],iterPeriod=10),
#qt.SnapshotEngine(fileBase='3d-',iterPeriod=200,label='snapshot'),
NewtonIntegrator(damping=0.4,gravity=[0,0,-9.81]),
PyRunner(iterPeriod=1,command='addPlotData()'),
PyRunner(command='calm()',iterPeriod=10,label='calmEngine')
]
O.dt=6e-7
def addPlotData():
plot.addData(exx=-triax.strain[0], eyy=-triax.strain[1],
ezz=-triax.strain[2], ev=(triax.strain[0]+triax.strain[1]+triax.strain[2]),
sxx=-triax.stress(triax.wall_right_id)[0],
syy=-triax.stress(triax.wall_top_id)[1],
szz=-triax.stress(triax.wall_front_id)[2],
q=-(triax.stress(triax.wall_front_id)[2]-triax.stress(triax.wall_right_id)[0]),
p=triax.meanStress, R=triax.stress(triax.wall_front_id)[2]/sigmaIso,
por=triax.porosity, i=O.iter,**O.energy,total=O.energy.total()),
#=
#= APPLYING CONFINING PRESSURE
#=
sigmaIso=-1e5
triax=TriaxialStressController(
maxMultiplier=1.000,
finalMaxMultiplier=1.000,
thickness = 0,
stressMask = 7,
internalCompaction=False,
max_vel=0.05,
)
O.engines=O.engines+[triax]
triax.goal1=sigmaIso
triax.goal2=sigmaIso
triax.goal3=sigmaIso
triax.wall_back_activated=True
while 1:
O.run(100, True)
unb=unbalancedForce()
meanS=(triax.stress(triax.wall_right_id)[0]+triax.stress(triax.wall_top_id)[1]+triax.stress(triax.wall_front_id)[2])/3
print('mean Stress:',triax.meanStress,'porosity:',
triax.porosity,'Unbalancing:',unb)
if abs(triax.meanStress)<100:
[Yade-users] [Question #694273]: Using clump and micro-strain field
New question #694273 on Yade: https://answers.launchpad.net/yade/+question/694273 Hello, Is it possible to get a microstrain field inside a sample when we are using a clump instead of a sphere? I want to have a microstrain field of a sample during the Triaxial test but I got a Segmentation fault. Thank you very much for your help in advance. Best Regards, Alireza -- 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 #694158]: Segmentation fault (core dumped) when I am using clump in Tiaxial test
New question #694158 on Yade: https://answers.launchpad.net/yade/+question/694158 Hello, I use Yadedaily on ubuntu 18 and 20. I want to simulate a triaxial test. When I use clump instead of spherical particles, I received this error: Segmentation fault (core dumped) I don't know how to fix it. I was wondering if you could help me. An example of my script is: from yade import utils, plot from yade import pack, qt from datetime import datetime import math import pylab import sys from numpy import * qtr=qt.Renderer() qtr.bgColor=(1,1,1) #== #= define the materials === #== O.materials.append(CohFrictMat(normalCohesion= 1e20, shearCohesion= 1e20, isCohesive= True, young=6.81e8, density=1377.5, poisson=0.3, frictionAngle= 0.31, fragile=False, label='aggregate-48')) O.materials.append(CohFrictMat(normalCohesion= 1e20, shearCohesion= 1e20, isCohesive= False, young=2e11, density=0, poisson=0.3, frictionAngle= 0, fragile=False, label='wall-frictionless')) #=== #=== define walls #=== walls=aabbWalls([(0,0,0),(4e-2,4e-2,4e-2)],thickness=0.001,oversizeFactor=1.5,material='wall-frictionless') wallIds=O.bodies.append(walls) #=== #===clump shape= #=== radz1=[0.355155e-3,0.505113e-3,0.397713e-3,0.465286e-3,0.484395e-3,0.394534e-3,0.493151e-3,0.487328e-3,0.613619e-3,0.413455e-3] poz1= [[1.13418e-3,-0.703895e-3,-1.20338e-3],[-0.390408e-3,0.476061e-3,-0.150612e-3],[-0.556545e-3,0.451341e-3,1.1495e-3],[-0.633942e-3,0.498253e-3,0.348231e-3],[0.0256934e-3,0.388855e-3,-0.733445e-3],[-0.218563e-3,0.504478e-3,1.54117e-3],[-0.319601e-3,0.104778e-3,0.742895e-3],[0.650678e-3,-0.76675e-3,-0.289908e-3],[0.0113115e-3,-0.207684e-3,0.00255944e-3],[0.594902e-3,-0.301473e-3,-0.878654e-3]] template1= [] template1.append(clumpTemplate(relRadii=radz1,relPositions=poz1)) #=== #=== define packing #=== coke=((1.875e-3,100),(1.875e-3,100),(1.875e-3,100)) nums=['t','t','t'] temps=[template1,template1,template1] mats=['aggregate-48','aggregate-48','aggregate-48'] tolerance=[(0.4752e-3),(0.4752e-3),(0.4752e-3)] for i in range(len(nums)): nums[i]=pack.SpherePack() nums[i].makeCloud((0.005,0.005,0.005),(0.035,0.035,0.035),rMean=coke[i][0],rRelFuzz=tolerance[i],num=coke[i][1]) O.bodies.append([utils.sphere(c,r,material=mats[i]) for c,r in nums[i]]) O.bodies.replaceByClumps(temps[i],[1.0],discretization=5) #=== #=== define Engine = #=== O.engines=[ ForceResetter(), InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Box_Aabb(),Bo1_Facet_Aabb()]), InteractionLoop([Ig2_Sphere_Sphere_ScGeom6D(),Ig2_Box_Sphere_ScGeom6D(),Ig2_Facet_Sphere_ScGeom6D()], [Ip2_CohFrictMat_CohFrictMat_CohFrictPhys()], [Law2_ScGeom_FrictPhys_CundallStrack(),Law2_ScGeom6D_CohFrictPhys_CohesionMoment()] ), TriaxialStateRecorder(iterPeriod=1,file='WallStresses'), NewtonIntegrator(damping=0.4,gravity=[0,0,-9.81]), PyRunner(iterPeriod=1,command='addPlotData()'), ] O.dt=5e-7 def addPlotData(): plot.addData(exx=-triax.strain[0], eyy=-triax.strain[1], ezz=-triax.strain[2], ev=(triax.strain[0]+triax.strain[1]+triax.strain[2]), sxx=-triax.stress(triax.wall_right_id)[0], syy=-triax.stress(triax.wall_top_id)[1], szz=-triax.stress(triax.wall_front_id)[2], q=-(triax.stress(triax.wall_front_id)[2]-triax.stress(triax.wall_right_id)[0]), p=triax.meanStress, R=triax.stress(triax.wall_front_id)[2]/sigmaIso, por=triax.porosity, i=O.iter,**O.energy,total=O.energy.total()), #= #= APPLYING CONFINING PRESSURE #= sigmaIso=-1e5 triax=TriaxialStressController( maxMultiplier=1.000, finalMaxMultiplier=1.000, thickness = 0, stressMask = 7, internalCompaction=False, max_vel=1, ) O.engines=O.engines+[triax] triax.goal1=sigmaIso triax.goal2=sigmaIso triax.goal3=sigmaIso triax.wall_back_activated=True while 1: O.run(100, True) unb=unbalancedForce() meanS=(triax.stress(triax.wall_right_id)[0]+triax.stress(t
[Yade-users] [Question #694002]: initial unbalancing and WallStresses
New question #694002 on Yade: https://answers.launchpad.net/yade/+question/694002 Hello, I wish you and your family perfect health. I am modeling a triaxial test. The problem that I have is that at the beginning of the simulation, some particles come out from the walls. Therefore, it destroys the shape of micro-strain contours. In addition, the initial stresses which are sensed at the walls are so high (more than 1 MPa). I am using YADE 2020.01a, and Ubuntu 20.04. My python version is 3.8.5. Attached, you can find my script. I know it should not be too long, but this is the problem that I have. I was wondering if you could give me a solution. Thank you very much. Best Regards, Alireza from yade import utils, plot from yade import pack, qt from datetime import datetime import math import pylab from numpy import * qtr=qt.Renderer() qtr.bgColor=(1,1,1) #== #= define the materials === #== O.materials.append(CohFrictMat(normalCohesion= 1e20, shearCohesion= 1e20, isCohesive= True, young=6.81e8, density=1377.5, poisson=0.3, frictionAngle= 0.31, fragile=False, label='aggregate-48')) O.materials.append(CohFrictMat(normalCohesion= 1e20, shearCohesion= 1e20, isCohesive= True, young=6.81e8, density=1532.2, poisson=0.3, frictionAngle= 0.31, fragile=False, label='aggregate-814')) O.materials.append(CohFrictMat(normalCohesion= 1e20, shearCohesion= 1e20, isCohesive= False, young=2e11, density=0, poisson=0.3, frictionAngle= 0, fragile=False, label='wall-frictionless')) #=== #=== define walls #=== walls=aabbWalls([(0,0,0),(10e-2,10e-2,10e-2)],thickness=0.001,oversizeFactor=1.5,material='wall-frictionless') wallIds=O.bodies.append(walls) #=== #=== define packing #=== nums=['t','t','t'] mats=['aggregate-48','aggregate-814','aggregate-814'] coke=((1.875e-3,3*1000),(0.9475e-3,3*2734),(0.50125e-3,3*24970)) color=((0,0,1),(0,1,0),(1,0,0)) tolerance=[(2e-4),(1e-4),(5e-5)] for i in range(len(nums)): nums[i]=pack.SpherePack() nums[i].makeCloud((5e-3,5e-3,5e-3),(9.95e-2,9.95e-2,9.95e-2),rMean=coke[i][0],rRelFuzz=tolerance[i],num=coke[i][1]) O.bodies.append([utils.sphere(c,r,material=mats[i],color=color[i]) for c,r in nums[i]]) #=== #=== define Engine = #=== O.engines=[ ForceResetter(), InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Box_Aabb(),Bo1_Facet_Aabb()]), InteractionLoop([Ig2_Sphere_Sphere_ScGeom6D(),Ig2_Box_Sphere_ScGeom6D(),Ig2_Facet_Sphere_ScGeom6D()], [Ip2_CohFrictMat_CohFrictMat_CohFrictPhys()], [Law2_ScGeom_FrictPhys_CundallStrack(),Law2_ScGeom6D_CohFrictPhys_CohesionMoment()] ), TriaxialStateRecorder(iterPeriod=1,file='WallStresses'), #VTKRecorder(fileName='3d-vtk-',recorders=['all'],iterPeriod=1000), #qt.SnapshotEngine(fileBase='3d-',iterPeriod=200,label='snapshot'), NewtonIntegrator(damping=0.4,gravity=[0,0,-9.81]), PyRunner(iterPeriod=1,command='addPlotData()'), ] #yade.qt.Controller(), yade.qt.View() O.dt=5e-8 #O.dt=.2*PWaveTimeStep() def addPlotData(): plot.addData(exx=-triax.strain[0], eyy=-triax.strain[1], ezz=-triax.strain[2], ev=(triax.strain[0]+triax.strain[1]+triax.strain[2]), sxx=-triax.stress(triax.wall_right_id)[0], syy=-triax.stress(triax.wall_top_id)[1], szz=-triax.stress(triax.wall_front_id)[2], mass=sum(b.state.mass for b in O.bodies if isinstance(b.shape,Sphere)), q=-(triax.stress(triax.wall_front_id)[2]-triax.stress(triax.wall_right_id)[0]), p=triax.meanStress, R=triax.stress(triax.wall_front_id)[2]/sigmaIso, por=triax.porosity, i=O.iter,**O.energy,total=O.energy.total()), #= #= APPLYING CONFINING PRESSURE #= sigmaIso=-5e4 triax=TriaxialStressController( maxMultiplier=1.000, finalMaxMultiplier=1.000, thickness = 0, stressMask = 7, internalCompaction=False, max_vel=0.05, ) O.engines=O.engines+[triax] triax.goal1=sigmaIso triax.goal2=sigmaIso triax.goal3=sigmaIso triax.wall_back_activated=True while 1: O.run(100, True) unb=unbalancedForce() meanS=(triax.stress(triax.wall_right_id)[0]+triax.stress(triax.wall_top_id)[1]+triax.stress(triax.wall_front_id)
Re: [Yade-users] [Question #688771]: run with Yadedaily
Question #688771 on Yade changed: https://answers.launchpad.net/yade/+question/688771 Status: Answered => Solved Alireza Sadeghi confirmed that the question is solved: Thank you, dear Jerome. -- 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 #688771]: run with Yadedaily
Question #688771 on Yade changed: https://answers.launchpad.net/yade/+question/688771 Alireza Sadeghi posted a new comment: Hi Ehsan, Thank you for your reply. My problem is "when I run a simulation with 12,000 particles with 8 CPUs, the software just uses one CPU and the other seven CPUs do not participate in the simulation. Therefore, it takes a long time. But, when I am using, for example, 1000 particles, all of my CPUs work 100%. Do you have a similar problem? Best Regards, Alireza -- 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 #688771]: run with Yadedaily
Question #688771 on Yade changed: https://answers.launchpad.net/yade/+question/688771 Alireza Sadeghi posted a new comment: Dear Robert, You are absolutely right. Sure, I will do those simulations in the three days later. Now, I am running a similar code (1000 particles+using clumps+periodic boundary) but all the particles are in the same size (with different templates but the overall size is the same). This code is running parallel very well. The problem arises when I add small particles in the aforementioned mixture. Best Regards, Alireza -- 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 #688771]: run with Yadedaily
Question #688771 on Yade changed: https://answers.launchpad.net/yade/+question/688771 Status: Needs information => Open Alireza Sadeghi gave more information on the question: Dear Jan, Thank you very much for your response. I applied your comments and I have this in my terminal: ForceResetter49 54037.627us 0.02% InsertionSortCollider49 269965135.677us 96.36% InteractionLoop 495101713.404us 1.82% "triax" 49 647611.067us 0.23% NewtonIntegrator 491206346.451us 0.43% forces sync 49 126156.097us 10.46% motion integration 491080081.651us 89.53% sync max vel 49 25.183us 0.00% terminate49 17.333us 0.00% TOTAL 1961206280.264us 99.99% "recorder" 47 3182985.26us 1.14% TOTAL 280157829.486us 100.00% In addition, when I used 8 CPU for the larger simulation, after two hours, the program used one CPU 100% and other seven CPUs worked at almost 50%. It means that, although all the cpus did not participate in the simulation equally, they have effect on the speed. Do you know how can I receive to faster calculation with my code? Thank you again. Best Regards, Alireza -- 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 #688771]: run with Yadedaily
Question #688771 on Yade changed: https://answers.launchpad.net/yade/+question/688771 Status: Answered => Open Alireza Sadeghi is still having a problem: No, I have no warning. In addition, I can use parallel when I use 500 particles. Now I have more than 12000 particles. Therefore, the number of bodies is not too small for parallel calculation. Thank you. Alireza -- 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 #688771]: run with Yadedaily
Question #688771 on Yade changed: https://answers.launchpad.net/yade/+question/688771 Status: Needs information => Open Alireza Sadeghi gave more information on the question: Hello Robert, Thank you very much. How can I fix my problem? I mean how can I use parallel calculation for this simulation? Is it not possible to run this code in parallel? Thanks a lot. Best Regards, Alireza -- 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 #688771]: run with Yadedaily
Question #688771 on Yade changed: https://answers.launchpad.net/yade/+question/688771 Status: Answered => Open Alireza Sadeghi is still having a problem: Hello, Thank you very much for your answers. I want to simulate a PeriTriaxController test with 12485 particles. I am using clump instead of spheres. Yes, you are right, when I used 500 particles in the same code, I can sun it in parallel. What should I do now? Thank you very much. Best Regards, Alireza P.S. my code is: from yade import utils, plot from yade import pack, qt from datetime import datetime qtr=qt.Renderer() qtr.bgColor=(1,1,1) #=== #==set a periodic boundary== #=== O.periodic=True O.cell.refSize=(5e-2,5e-2,5e-2) #== #= define the materials === #== O.materials.append(CohFrictMat(normalCohesion= 1e20, shearCohesion= 1e20, isCohesive= True, young=6.81e8, density=1377.5, poisson=0.3, frictionAngle= 0.31, fragile=False, label='aggregate-48')) O.materials.append(CohFrictMat(normalCohesion= 1e20, shearCohesion= 1e20, isCohesive= True, young=6.81e8, density=1532.2, poisson=0.3, frictionAngle= 0.31, fragile=False, label='aggregate-814')) O.materials.append(CohFrictMat(normalCohesion= 1e20, shearCohesion= 1e20, isCohesive= False, young=2e11, density=1523.6, poisson=0.3, frictionAngle= 0.28, fragile=False, label='wall')) #=== #===clump shape= #=== radz1=[0.355155e-3,0.505113e-3,0.397713e-3,0.465286e-3,0.484395e-3,0.394534e-3,0.493151e-3,0.487328e-3,0.613619e-3,0.413455e-3] poz1= [[1.13418e-3,-0.703895e-3,-1.20338e-3],[-0.390408e-3,0.476061e-3,-0.150612e-3],[-0.556545e-3,0.451341e-3,1.1495e-3],[-0.633942e-3,0.498253e-3,0.348231e-3],[0.0256934e-3,0.388855e-3,-0.733445e-3],[-0.218563e-3,0.504478e-3,1.54117e-3],[-0.319601e-3,0.104778e-3,0.742895e-3],[0.650678e-3,-0.76675e-3,-0.289908e-3],[0.0113115e-3,-0.207684e-3,0.00255944e-3],[0.594902e-3,-0.301473e-3,-0.878654e-3]] template1= [] template1.append(clumpTemplate(relRadii=radz1,relPositions=poz1)) radz2=[0.330164e-3,0.504115e-3,0.399587e-3,0.614205e-3,0.466444e-3,0.495302e-3,0.394324e-3,0.486898e-3,0.444037e-3,0.489396e-3] poz2= [[1.16386e-3,-0.70706e-3,-1.25222e-3],[-0.39596e-3,0.482385e-3,-0.144097e-3],[-0.554928e-3,0.448475e-3,1.14685e-3],[-0.00361766e-3,-0.198211e-3, 0.00106559e-3],[-0.633362e-3,0.490424e-3,0.357391e-3],[0.0434148e-3,0.367924e-3,-0.736319e-3],[-0.218749e-3,0.504703e-3,1.54165e-3],[-0.311777e-3, 0.101954e-3,0.750235e-3],[0.621565e-3,-0.779387e-3,-0.179029e-3],[0.74e-3,-0.503202e-3,-0.795602e-3]] template2= [] template2.append(clumpTemplate(relRadii=radz2,relPositions=poz2)) radz3=[0.959101e-3,0.774782e-3,0.924682e-3,0.882986e-3,0.572207e-3,0.875338e-3,0.499054e-3,0.582586e-3,1.03184e-3,0.561428e-3] poz3=[[-0.742455e-3,-0.539002e-3,0.030705e-3],[0.800775e-3,1.00778e-3,0.366095e-3],[-0.217508e-3,-0.423924e-3,-0.671114e-3],[0.65016e-3,0.741405e-3,0.305558e-3],[-0.770717e-3,0.423006e-3,0.33259e-3],[0.20475e-3,-0.707669e-3,-0.332745e-3],[-0.0356616e-3,1.34122e-3,0.703809e-3],[0.64889e-3,0.340976e-3,-0.456228e-3],[0.258817e-3,0.13357e-3,0.178886e-3],[0.35466e-3,1.36934e-3,0.583507e-3]] template3= [] template3.append(clumpTemplate(relRadii=radz3,relPositions=poz3)) radz4=[1.18842e-3,1.0381e-3,0.664711e-3,0.531753e-3,0.853808e-3,0.789023e-3,0.717061e-3,1.0081e-3,0.967001e-3,0.535189e-3] poz4=[[-0.0354242e-3,-0.245642e-3,-0.0514299e-3],[0.444228e-3,-0.342353e-3,-0.327639e-3],[-0.815227e-3,1.41647e-3,-0.246128e-3],[1.38493e-3,-0.405625e-3,-0.69207e-3],[0.608272e-3,0.0171296e-3,0.863712e-3],[-0.191982e-3,-1.05752e-3,-0.341411e-3],[-1.08258e-3,-0.336001e-3,-0.482846e-3],[-0.396508e-3,0.47029e-3,0.0820364e-3],[0.486211e-3,0.345623e-3,0.492724e-3],[-0.0502768e-3,0.977767e-3,0.467401e-3]] template4= [] template4.append(clumpTemplate(relRadii=radz4,relPositions=poz4)) radz5= [0.734571e-3,0.933485e-3,0.483025e-3,0.563064e-3,0.457405e-3,0.955996e-3,0.722849e-3,0.71635e-3,0.486014e-3,0.735338e-3] poz5= [[-0.452026e-3,-0.455606e-3,-0.0619758e-3],[-0.126025e-3,0.168901e-3,0.0762994e-3],[-1.6211e-3,-0.316016e-3,1.23288e-3],[-0.715097e-3,0.497059e-3,0.0664344e-3],[0.220643e-3,0.879083e-3,0.344862e-3],[0.394206e-3,-0.220853e-3,-0.189581e-3],[-0.998749e-3,0.0346725e-3,0.758223e-3],[0.518915e-3,0.743702e-3,-0.383334e-3],[0.877614e-3,0.598898e-3,-0.809004e-3],[1.00845e-3,-0.510125e-3,-0.532033e-3]] template5= [] template5.append(clumpTemplate(relRadii=radz5,relPositions=poz5))
[Yade-users] [Question #688771]: run with Yadedaily
New question #688771 on Yade: https://answers.launchpad.net/yade/+question/688771 Hello All, I am using yadedaily and I want to use parallel computation. I typed "yadedaily -j8 test.py", but the yade just used one CPU. I checked it in the System Monitor in Ubunto 18. My yadedaily is the version of 20200213-3366~7bb55ad~bionic1. Could you please help me if you have some ideas? Thank you very much. Best Regards, Alireza -- 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 #688073]: using clump in PeriTriaxController
Question #688073 on Yade changed: https://answers.launchpad.net/yade/+question/688073 Status: Needs information => Open Alireza Sadeghi gave more information on the question: Hello Jan, Thank you very much for your reply. When I work with YADE (2018.02b), I just type "yade -j4 test.py" in the test.py file directory. Now, I am using yade-2020-01-git-146d6de. In the directory of "yade-2020-01-git-146d6de", I typed "./yade-2020-01-02.git-146d6de". I have: alsac8@GMN-HA-02LX:~/yade-dev/install/bin$ ./yade-2020-01-02.git-146d6de Welcome to Yade 2020-01-02.git-146d6de Using python version: 3.6.9 (default, Nov 7 2019, 10:44:02) [GCC 8.3.0] TCP python prompt on localhost:9000, auth cookie `dasksc' XMLRPC info provider on http://localhost:21000 [[ ^L clears screen, ^U kills line. F12 controller, F11 3D view (press "h" in 3D view for help), F10 both, F9 generator, F8 plot. ]] In [1]: Now, how should I run my code (for example test.py which is located in a directory different with directory of yade-2020-01-02.git-146d6de) ? Thank you very much. Best Regards, Alireza -- 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 #688073]: using clump in PeriTriaxController
Question #688073 on Yade changed: https://answers.launchpad.net/yade/+question/688073 Status: Answered => Open Alireza Sadeghi is still having a problem: Dear Jan, Thank you very much for your response. I installed yade-2020-01-git- 146d6de. Now I don't know how should I run my code. Because in YADE, I just type for example (yade -j2 test.py), but now I have to go through the folder that yade-2020-01-git-146d6de is located. Thank you very much. Best Regards, Alireza -- 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 #688073]: using clump in PeriTriaxController
Question #688073 on Yade changed: https://answers.launchpad.net/yade/+question/688073 Status: Answered => Open Alireza Sadeghi is still having a problem: Dear Jan, I see that yade 2018.02b is newest one for Ubuntu 18.04. 1. Is it possible to give me the exact version of YADE you are using? 2. How about yadedaily?I installed it on my computer but I cannot run my code with it. How does it work?Does it work as the same as YADE with writing "yade test.py" in the terminal?? 3. For the last version, Should I install Ubuntu 19??If yes, do you think that it could solve my problem? Thank you very much. best regards, Alireza -- 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 #688073]: using clump in PeriTriaxController
Question #688073 on Yade changed: https://answers.launchpad.net/yade/+question/688073 Alireza Sadeghi posted a new comment: Dear Jan, Thank you very much for your response. I will install the newest version of YADE and see what will be happened. I will definitely let you know about the problem. Thank you very much. Best Regards, Alireza -- 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 #688073]: using clump in PeriTriaxController
Question #688073 on Yade changed: https://answers.launchpad.net/yade/+question/688073 Status: Needs information => Open Alireza Sadeghi gave more information on the question: Hello Jan, Thank you very much for your reply. Do you mean that you get zero stress at the boundaries for 1M initial steps? I am using YADE 2018.02b. Does the problem come from my YADE version? Thank you very much for your help. Best Regards, Alireza -- 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 #688073]: using clump in PeriTriaxController
New question #688073 on Yade: https://answers.launchpad.net/yade/+question/688073 Hello All, I want to use clump shape particles instead of spherical particles in the PeriTriaxController. As you can see in the following code, there is only one particle in the middle of a cell. There is no interactions between the boundaries and the particle at the begininng. Why the stress at the boundaries are not zero? I asked this question several times, and nobody answered me. Therfore, the questions were expired and removed. I have stopped in this position for two months. I will be so grateful if you could help me. Thank you very much for your help in advance. Best Regards, Alireza ##code## from yade import utils, plot from yade import pack, qt from datetime import datetime #==set a periodic boundary== O.periodic=True O.cell.refSize=(10e-3,10e-3,10e-3) #= define the materials === #O.interactions.clear() O.materials.append(CohFrictMat(normalCohesion= 1e20, shearCohesion= 1e20, isCohesive= True, young=6.81e8, density=1377.5, poisson=0.3, frictionAngle= 0.31, fragile=False, label='aggregate-48')) #===clump shape= radz1=[3e-3,3e-3] poz1= [[0,0,-1e-3],[1e-3,0,0]] template1= [] template1.append(clumpTemplate(relRadii=radz1,relPositions=poz1)) nums=['t'] mats=['aggregate-48'] coke=(1e-3,1) color=((0,0,1)) tolerance=[(2e-4)] temps=[template1] nums=pack.SpherePack() nums.makeCloud((3e-3,3e-3,3e-3),(7e-3,7e-3,7e-3),rMean=coke[0],rRelFuzz=tolerance[0],num=coke[1]) O.bodies.append([utils.sphere(c,r,material=mats[0],color=color[0]) for c,r in nums]) O.bodies.replaceByClumps(temps[0],[1.0],discretization=5) #O.interactions.clear() #Engine sigmaIso=-1e5 O.engines=[ ForceResetter(), InsertionSortCollider([Bo1_Sphere_Aabb()]), InteractionLoop( [Ig2_Sphere_Sphere_ScGeom()], [Ip2_FrictMat_FrictMat_FrictPhys()], [Law2_ScGeom_FrictPhys_CundallStrack()] ), PeriTriaxController(label='triax', goal=(sigmaIso,sigmaIso,sigmaIso),stressMask=7, dynCell=True,maxStrainRate=(10,10,10), maxUnbalanced=.01,relStressTol=1e-3, doneHook='compactionFinished()' ), NewtonIntegrator(damping=0.4,gravity=[0,0,0]), PyRunner(command='history()',iterPeriod=1,label='recorder'), ] O.dt=5e-8 def history(): plot.addData(unbalanced=unbalancedForce(), sxx=-triax.stress[0],syy=-triax.stress[1],szz=-triax.stress[2], exx=-triax.strain[0],eyy=-triax.strain[1],ezz=-triax.strain[2], q=-(triax.stress[2]-triax.stress[1]), p=-(triax.stress[2]+triax.stress[1]+triax.stress[0])/3, por=porosity(),i=O.iter,) print 'porosity:', porosity(),'Unbalancing:',unbalancedForce(),'Strain',-triax.strain[2], 'Stress',-triax.stress[2] def compactionFinished(): O.pause() -- 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 #687396]: Initial stress when I use clump instead of spherical particles
New question #687396 on Yade: https://answers.launchpad.net/yade/+question/687396 Hello All, I want to use clumps instead of spherical particles in PeriTriaxController. My problem is at the initial time, the stresses on the wall is more that my goal stress. I don't know where this stresses come from. I will be gratefule if you could help me. Thank you very much. Best Regards Alireza #=== from yade import utils, plot from yade import pack, qt from datetime import datetime #==set a periodic boundary== O.periodic=True O.cell.refSize=(2e-2,2e-2,2e-2) #= define the materials === O.materials.append(CohFrictMat(normalCohesion= 1e20, shearCohesion= 1e20, isCohesive= True, young=6.81e8, density=1377.5, poisson=0.3, frictionAngle= 0.31, fragile=False, label='aggregate-48')) #===clump shape= radz1=[3e-3,3e-3] poz1= [[-1e-3,0,0],[1e-3,0,0]] template1= [] template1.append(clumpTemplate(relRadii=radz1,relPositions=poz1)) nums=['t'] mats=['aggregate-48'] coke=(3e-3,4) color=((0,0,1)) tolerance=[(2e-4)] temps=[template1] nums=pack.SpherePack() nums.makeCloud((0,0,0),(2e-2,2e-2,2e-2),rMean=coke[0],rRelFuzz=tolerance[0],num=coke[1]) O.bodies.append([utils.sphere(c,r,material=mats[0],color=color[0]) for c,r in nums]) O.bodies.replaceByClumps(temps[0],[1.0],discretization=5) #Engine sigmaIso=-1e5 O.engines=[ ForceResetter(), InsertionSortCollider([Bo1_Sphere_Aabb()]), InteractionLoop( [Ig2_Sphere_Sphere_ScGeom()], [Ip2_FrictMat_FrictMat_FrictPhys()], [Law2_ScGeom_FrictPhys_CundallStrack()] ), PeriTriaxController(label='triax', goal=(sigmaIso,sigmaIso,sigmaIso),stressMask=7, dynCell=True,maxStrainRate=(0.1,0.1,0.1), maxUnbalanced=.01,relStressTol=1e-3, doneHook='compactionFinished()' ), NewtonIntegrator(damping=0.4,gravity=[0,0,-10]), PyRunner(command='history()',iterPeriod=10,label='recorder'), ] O.dt=5e-8 def history(): plot.addData(unbalanced=unbalancedForce(), sxx=-triax.stress[0],syy=-triax.stress[1],szz=-triax.stress[2], exx=-triax.strain[0],eyy=-triax.strain[1],ezz=-triax.strain[2], q=-(triax.stress[2]-triax.stress[1]), p=-(triax.stress[2]+triax.stress[1]+triax.stress[0])/3, por=porosity(),i=O.iter,) print 'porosity:', porosity(),'Unbalance', unbalancedForce(), 'Strain',-triax.strain[2], 'Stress',-triax.stress[2] -- 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 #686343]: initial mean stress in Triaxial compaction when I used clumps
Question #686343 on Yade changed: https://answers.launchpad.net/yade/+question/686343 Status: Open => Solved Alireza Sadeghi confirmed that the question is solved: Hello Jerome, Thank you very much for your reply. I checked it and my problem is solved. Thank you again. Best Regards, Alireza -- 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 #686343]: initial mean stress in Triaxial compaction when I used clumps
Question #686343 on Yade changed: https://answers.launchpad.net/yade/+question/686343 Alireza Sadeghi gave more information on the question: Dear Jerome, Thank you very much for your answer. For particles (instead of clumps) the initial stress is zero, but when I use clumps, the initial stress is zo high. The clumps did not touch the boundaries at the first time. why does it show 1 kPa mean stress? =>if I were you I would check your initial creation of sample to check things are as you want them to be. What does it mean?do you mean that I should check each clumps individually? Thank you very much. Best Regards, Alireza -- 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 #686343]: initial mean stress in Triaxial compaction when I used clumps
Question #686343 on Yade changed: https://answers.launchpad.net/yade/+question/686343 Alireza Sadeghi posted a new comment: Also, I have the same problem when I am using PeriTriaxController. -- 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 #686343]: initial mean stress in Triaxial compaction when I used clumps
Question #686343 on Yade changed: https://answers.launchpad.net/yade/+question/686343 Description changed to: Hello All, I have a problem when I am using triaxial compaction on clumps particles. The initial mean stress is not zero, it means when I run the code, at the first steps in which porosity is near 1, the mean stress is not zero. For example in my code [1], the initial mean stress is near 1 kPa. This can cause problems when I increase the number of particles. For example, when I used 1 particles, the initial mean stress is more that 2 MPa, whereas my target confining pressure is 100 kPa. Do you have any experiences like this? Thank you very much. Best Regards, Alireza [1] #code# O.reset() from yade import utils, plot from yade import pack, qt from datetime import datetime #==COKE AGGREGATE CLUMP TEMPLATES radz1=[0.355155e-3,0.505113e-3,0.397713e-3,0.465286e-3,0.484395e-3,0.394534e-3,0.493151e-3,0.487328e-3,0.613619e-3,0.413455e-3] poz1= [[1.13418e-3,-0.703895e-3,-1.20338e-3],[-0.390408e-3,0.476061e-3,-0.150612e-3],[-0.556545e-3,0.451341e-3,1.1495e-3],[-0.633942e-3,0.498253e-3,0.348231e-3],[0.0256934e-3,0.388855e-3,-0.733445e-3],[-0.218563e-3,0.504478e-3,1.54117e-3],[-0.319601e-3,0.104778e-3,0.742895e-3],[0.650678e-3,-0.76675e-3,-0.289908e-3],[0.0113115e-3,-0.207684e-3,0.00255944e-3],[0.594902e-3,-0.301473e-3,-0.878654e-3]] template1= [] template1.append(clumpTemplate(relRadii=radz1,relPositions=poz1)) radz2=[0.330164e-3,0.504115e-3,0.399587e-3,0.614205e-3,0.466444e-3,0.495302e-3,0.394324e-3,0.486898e-3,0.444037e-3,0.489396e-3] poz2= [[1.16386e-3,-0.70706e-3,-1.25222e-3],[-0.39596e-3,0.482385e-3,-0.144097e-3],[-0.554928e-3,0.448475e-3,1.14685e-3],[-0.00361766e-3,-0.198211e-3, 0.00106559e-3],[-0.633362e-3,0.490424e-3,0.357391e-3],[0.0434148e-3,0.367924e-3,-0.736319e-3],[-0.218749e-3,0.504703e-3,1.54165e-3],[-0.311777e-3, 0.101954e-3,0.750235e-3],[0.621565e-3,-0.779387e-3,-0.179029e-3],[0.74e-3,-0.503202e-3,-0.795602e-3]] template2= [] template2.append(clumpTemplate(relRadii=radz2,relPositions=poz2)) #= define the materials === O.materials.append(CohFrictMat(normalCohesion= 1e20, shearCohesion= 1e20, isCohesive= True, young=1.95e7, density=1532.2, poisson=0.3, frictionAngle= 0.0, fragile=False, label='aggregate-814')) O.materials.append(CohFrictMat(normalCohesion= 1e20, shearCohesion= 1e20, isCohesive= False, young=4e9, density=1523.6, poisson=0.3, frictionAngle= 0.0, fragile=False, label='wall')) #= creating walls == walls=aabbWalls([(-15e-3,-15e-3,-15e-3),(15e-3,15e-3,15e-3)],thickness=0.0003,oversizeFactor=1.0,material='wall') wallIds=O.bodies.append(walls) ### DEFINING ENGINES ### O.engines=[ ForceResetter(), InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Box_Aabb(),Bo1_Facet_Aabb()]), InteractionLoop([Ig2_Sphere_Sphere_ScGeom6D(),Ig2_Box_Sphere_ScGeom6D(),Ig2_Facet_Sphere_ScGeom6D()], [Ip2_CohFrictMat_CohFrictMat_CohFrictPhys()], [Law2_ScGeom_FrictPhys_CundallStrack(),Law2_ScGeom6D_CohFrictPhys_CohesionMoment()] ), NewtonIntegrator(damping=0.4,gravity=[0,0,-10]) ,PyRunner(command='calm()',iterPeriod=10,label='calmEngine') ] #=Generating the aggregates== coke=((4.303694e-04,500),(3.971601e-04,500)) nums=['t','t'] temps=[template1,template2] mats=['aggregate-814','aggregate-814'] for i in range(len(nums)): nums[i]=pack.SpherePack() nums[i].makeCloud((-0.0145,-0.0145,-0.0145),(0.0145,0.0145,0.0145),rMean=coke[i][0],rRelFuzz=0.0,num=coke[i][1]) O.bodies.append([utils.sphere(c,r,material=mats[i]) for c,r in nums[i]]) O.bodies.replaceByClumps(temps[i],[1.0],discretization=5) #=== #=== Compaction #=== triax=TriaxialStressController( maxMultiplier=1.000, finalMaxMultiplier=1.000, thickness = 0, stressMask = 7, internalCompaction=False, ) O.engines=O.engines+[triax] triax.goal1=-1.0e5 triax.goal2=-1.0e5 triax.goal3=-1.0e5 triax.wall_back_activated=True O.dt=2e-6 while 1: O.run(100, True) unb=unbalancedForce() meanS=(triax.stress(triax.wall_right_id)[0]+triax.stress(triax.wall_top_id)[1]+triax.stress(triax.wall_front_id)[2])/3 print 'mean Stress:',triax.meanStress,'porosity:', triax.porosity,'meanS:',unb if triax.porosity<0.4 and abs(-1e5-triax.meanStress)/abs(-1e5)<0.01: print 'Isotropic strain1:',-triax.strain[0], 'Isotropic strain 2:',-triax.strain[1], 'Isotropic strain 3:',-triax.strain[2] break print "### Isotropic state saved ###" triax.depth0=triax.depth triax.height0=tri
[Yade-users] [Question #686343]: initial mean stress in Triaxial compaction when I used clumps
New question #686343 on Yade: https://answers.launchpad.net/yade/+question/686343 Hello All, I incounter to a problem when I was using triaxial compaction on clumps particles. The initial mean stress is not zero, it means when I run the code, at the first steps in which porosity is near 1, the mean stress is not zero. For example in my code, the initial mean stress is near 1 kPa. This can cause problem when I increase the number of particles. For example, when I used 1 particles, the initial mean stress is more that 2 MPa, whereas my target confining pressure is 100 kPa. Do you have some experiences like that? Thank you very much. Best Regards, Alireza #code# O.reset() from yade import utils, plot from yade import pack, qt from datetime import datetime #==COKE AGGREGATE CLUMP TEMPLATES radz1=[0.355155e-3,0.505113e-3,0.397713e-3,0.465286e-3,0.484395e-3,0.394534e-3,0.493151e-3,0.487328e-3,0.613619e-3,0.413455e-3] poz1= [[1.13418e-3,-0.703895e-3,-1.20338e-3],[-0.390408e-3,0.476061e-3,-0.150612e-3],[-0.556545e-3,0.451341e-3,1.1495e-3],[-0.633942e-3,0.498253e-3,0.348231e-3],[0.0256934e-3,0.388855e-3,-0.733445e-3],[-0.218563e-3,0.504478e-3,1.54117e-3],[-0.319601e-3,0.104778e-3,0.742895e-3],[0.650678e-3,-0.76675e-3,-0.289908e-3],[0.0113115e-3,-0.207684e-3,0.00255944e-3],[0.594902e-3,-0.301473e-3,-0.878654e-3]] template1= [] template1.append(clumpTemplate(relRadii=radz1,relPositions=poz1)) radz2=[0.330164e-3,0.504115e-3,0.399587e-3,0.614205e-3,0.466444e-3,0.495302e-3,0.394324e-3,0.486898e-3,0.444037e-3,0.489396e-3] poz2= [[1.16386e-3,-0.70706e-3,-1.25222e-3],[-0.39596e-3,0.482385e-3,-0.144097e-3],[-0.554928e-3,0.448475e-3,1.14685e-3],[-0.00361766e-3,-0.198211e-3, 0.00106559e-3],[-0.633362e-3,0.490424e-3,0.357391e-3],[0.0434148e-3,0.367924e-3,-0.736319e-3],[-0.218749e-3,0.504703e-3,1.54165e-3],[-0.311777e-3, 0.101954e-3,0.750235e-3],[0.621565e-3,-0.779387e-3,-0.179029e-3],[0.74e-3,-0.503202e-3,-0.795602e-3]] template2= [] template2.append(clumpTemplate(relRadii=radz2,relPositions=poz2)) #= define the materials === O.materials.append(CohFrictMat(normalCohesion= 1e20, shearCohesion= 1e20, isCohesive= True, young=1.95e7, density=1532.2, poisson=0.3, frictionAngle= 0.0, fragile=False, label='aggregate-814')) O.materials.append(CohFrictMat(normalCohesion= 1e20, shearCohesion= 1e20, isCohesive= False, young=4e9, density=1523.6, poisson=0.3, frictionAngle= 0.0, fragile=False, label='wall')) #= creating walls == walls=aabbWalls([(-15e-3,-15e-3,-15e-3),(15e-3,15e-3,15e-3)],thickness=0.0003,oversizeFactor=1.0,material='wall') wallIds=O.bodies.append(walls) ### DEFINING ENGINES ### O.engines=[ ForceResetter(), InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Box_Aabb(),Bo1_Facet_Aabb()]), InteractionLoop([Ig2_Sphere_Sphere_ScGeom6D(),Ig2_Box_Sphere_ScGeom6D(),Ig2_Facet_Sphere_ScGeom6D()], [Ip2_CohFrictMat_CohFrictMat_CohFrictPhys()], [Law2_ScGeom_FrictPhys_CundallStrack(),Law2_ScGeom6D_CohFrictPhys_CohesionMoment()] ), NewtonIntegrator(damping=0.4,gravity=[0,0,-10]) ,PyRunner(command='calm()',iterPeriod=10,label='calmEngine') ] #=Generating the aggregates== coke=((4.303694e-04,500),(3.971601e-04,500)) nums=['t','t'] temps=[template1,template2] mats=['aggregate-814','aggregate-814'] for i in range(len(nums)): nums[i]=pack.SpherePack() nums[i].makeCloud((-0.0145,-0.0145,-0.0145),(0.0145,0.0145,0.0145),rMean=coke[i][0],rRelFuzz=0.0,num=coke[i][1]) O.bodies.append([utils.sphere(c,r,material=mats[i]) for c,r in nums[i]]) O.bodies.replaceByClumps(temps[i],[1.0],discretization=5) #=== #=== Compaction #=== triax=TriaxialStressController( maxMultiplier=1.000, finalMaxMultiplier=1.000, thickness = 0, stressMask = 7, internalCompaction=False, ) O.engines=O.engines+[triax] triax.goal1=-1.0e5 triax.goal2=-1.0e5 triax.goal3=-1.0e5 triax.wall_back_activated=True O.dt=2e-6 while 1: O.run(100, True) unb=unbalancedForce() meanS=(triax.stress(triax.wall_right_id)[0]+triax.stress(triax.wall_top_id)[1]+triax.stress(triax.wall_front_id)[2])/3 print 'mean Stress:',triax.meanStress,'porosity:', triax.porosity,'meanS:',unb if triax.porosity<0.4 and abs(-1e5-triax.meanStress)/abs(-1e5)<0.01: print 'Isotropic strain1:',-triax.strain[0], 'Isotropic strain 2:',-triax.strain[1], 'Isotropic strain 3:',-triax.strain[2] break print "### Isotropic state saved ###" triax.depth0=triax.depth triax.height0=triax.height triax.width0=tri
Re: [Yade-users] [Question #686273]: micro strain in clumps
Question #686273 on Yade changed: https://answers.launchpad.net/yade/+question/686273 Status: Needs information => Open Alireza Sadeghi gave more information on the question: Dear Robert, Thank you very much for your answer. It does not show any errors in the terminal. Just I have something similar this: . . . 10398 : Vh==NULL!! id=10398 Point=-0.0135973 0.0130332 0.00710087 rad=0 10409 : Vh==NULL!! id=10409 Point=0.00605096 -0.00551408 -0.00619671 rad=0 Triangulated Grains : 8687 sym_grad_u_total_g (wrong averaged strain): 0.451431 -0.0578464 -0.166245 -0.0578464 0.562565 -0.0421092 -0.166245 -0.0421092 0.16872 Total volume = 0.168629, grad_u = 0.451431 -0.0917767 -0.254586 -0.0239162 0.562565 -0.0416519 -0.0779029 -0.0425665 0.16872 sym_grad_u (true average strain): 0.451431 -0.0578464 -0.166245 -0.0578464 0.562565 -0.0421092 -0.166245 -0.0421092 0.16872 Macro strain : 0 0 0 0 0 0 0 0 0 I can not open the vtk files in paraview. Thanks a lot. Best Regards, Alireza -- 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 #686273]: micro strain in clumps
New question #686273 on Yade: https://answers.launchpad.net/yade/+question/686273 Hello All, I want to simulate a triaxial test of sample of granular material with real shape. I am using clumps instead of spheres. I want to plot micro strain of sampel in paraview. Is it possible to get micro strain when the particles are clumps? Because my code does not work. Thank you very much for your help in advance. Best Regards, Alireza my code# O.reset() from yade import utils, plot from yade import pack, qt from datetime import datetime #==COKE AGGREGATE CLUMP TEMPLATES radz1=[0.355155e-3,0.505113e-3,0.397713e-3,0.465286e-3,0.484395e-3,0.394534e-3,0.493151e-3,0.487328e-3,0.613619e-3,0.413455e-3] poz1= [[1.13418e-3,-0.703895e-3,-1.20338e-3],[-0.390408e-3,0.476061e-3,-0.150612e-3],[-0.556545e-3,0.451341e-3,1.1495e-3],[-0.633942e-3,0.498253e-3,0.348231e-3],[0.0256934e-3,0.388855e-3,-0.733445e-3],[-0.218563e-3,0.504478e-3,1.54117e-3],[-0.319601e-3,0.104778e-3,0.742895e-3],[0.650678e-3,-0.76675e-3,-0.289908e-3],[0.0113115e-3,-0.207684e-3,0.00255944e-3],[0.594902e-3,-0.301473e-3,-0.878654e-3]] template1= [] template1.append(clumpTemplate(relRadii=radz1,relPositions=poz1)) radz2=[0.330164e-3,0.504115e-3,0.399587e-3,0.614205e-3,0.466444e-3,0.495302e-3,0.394324e-3,0.486898e-3,0.444037e-3,0.489396e-3] poz2= [[1.16386e-3,-0.70706e-3,-1.25222e-3],[-0.39596e-3,0.482385e-3,-0.144097e-3],[-0.554928e-3,0.448475e-3,1.14685e-3],[-0.00361766e-3,-0.198211e-3, 0.00106559e-3],[-0.633362e-3,0.490424e-3,0.357391e-3],[0.0434148e-3,0.367924e-3,-0.736319e-3],[-0.218749e-3,0.504703e-3,1.54165e-3],[-0.311777e-3, 0.101954e-3,0.750235e-3],[0.621565e-3,-0.779387e-3,-0.179029e-3],[0.74e-3,-0.503202e-3,-0.795602e-3]] template2= [] template2.append(clumpTemplate(relRadii=radz2,relPositions=poz2)) radz3=[0.959101e-3,0.774782e-3,0.924682e-3,0.882986e-3,0.572207e-3,0.875338e-3,0.499054e-3,0.582586e-3,1.03184e-3,0.561428e-3] poz3=[[-0.742455e-3,-0.539002e-3,0.030705e-3],[0.800775e-3,1.00778e-3,0.366095e-3],[-0.217508e-3,-0.423924e-3,-0.671114e-3],[0.65016e-3,0.741405e-3,0.305558e-3],[-0.770717e-3,0.423006e-3,0.33259e-3],[0.20475e-3,-0.707669e-3,-0.332745e-3],[-0.0356616e-3,1.34122e-3,0.703809e-3],[0.64889e-3,0.340976e-3,-0.456228e-3],[0.258817e-3,0.13357e-3,0.178886e-3],[0.35466e-3,1.36934e-3,0.583507e-3]] template3= [] template3.append(clumpTemplate(relRadii=radz3,relPositions=poz3)) radz4=[1.18842e-3,1.0381e-3,0.664711e-3,0.531753e-3,0.853808e-3,0.789023e-3,0.717061e-3,1.0081e-3,0.967001e-3,0.535189e-3] poz4=[[-0.0354242e-3,-0.245642e-3,-0.0514299e-3],[0.444228e-3,-0.342353e-3,-0.327639e-3],[-0.815227e-3,1.41647e-3,-0.246128e-3],[1.38493e-3,-0.405625e-3,-0.69207e-3],[0.608272e-3,0.0171296e-3,0.863712e-3],[-0.191982e-3,-1.05752e-3,-0.341411e-3],[-1.08258e-3,-0.336001e-3,-0.482846e-3],[-0.396508e-3,0.47029e-3,0.0820364e-3],[0.486211e-3,0.345623e-3,0.492724e-3],[-0.0502768e-3,0.977767e-3,0.467401e-3]] template4= [] template4.append(clumpTemplate(relRadii=radz4,relPositions=poz4)) #= define the materials === O.materials.append(CohFrictMat(normalCohesion= 1e20, shearCohesion= 1e20, isCohesive= True, young=1.95e7, density=1532.2, poisson=0.3, frictionAngle= 0.0, fragile=False, label='aggregate-814')) O.materials.append(CohFrictMat(normalCohesion= 1e20, shearCohesion= 1e20, isCohesive= False, young=4e9, density=1523.6, poisson=0.3, frictionAngle= 0.0, fragile=False, label='wall')) #= creating walls == #fc=yade.geom.facetCylinder((0,0,0.05),12.7e-3, 0.1, orientation=Quaternion((0, 0, 1), 0), segmentsNumber=100, wallMask=7, angleRange=None, closeGap=False, radiusTopInner=-1, radiusBottomInner=-1,material='wall') #O.bodies.append(fc) walls=aabbWalls([(-0.05,-0.05,-0.05),(0.05,0.05,0.05)],thickness=0.0003,oversizeFactor=1.0,material='wall') wallIds=O.bodies.append(walls) ### DEFINING ENGINES ### O.engines=[ ForceResetter(), InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Box_Aabb(),Bo1_Facet_Aabb()]), InteractionLoop([Ig2_Sphere_Sphere_ScGeom6D(),Ig2_Box_Sphere_ScGeom6D(),Ig2_Facet_Sphere_ScGeom6D()], [Ip2_CohFrictMat_CohFrictMat_CohFrictPhys()], [Law2_ScGeom_FrictPhys_CundallStrack(),Law2_ScGeom6D_CohFrictPhys_CohesionMoment()] ), # TriaxialStateRecorder(iterPeriod=1000,file='WallStresses'), NewtonIntegrator(damping=0.4,gravity=[0,0,-10]) ,PyRunner(command='calm()',iterPeriod=10,label='calmEngine') ] O.dt=1e-7 #Generating the aggregates=== coke=((1.875e-3,100),(0.9475e-3,367),(0.50125e-3,500),(0.50125e-3,500)) nums=['t','t','t','t'] t
Re: [Yade-users] [Question #686024]: micro strain-paraview
Question #686024 on Yade changed: https://answers.launchpad.net/yade/+question/686024 Status: Answered => Open Alireza Sadeghi is still having a problem: Dear Bruno, Thank you very much for your reply. I am not sure that the problem is in post processing. Because in low particles sample, the strain on the boundary is occurding to prescribed deformation, but when I increased the number of particles, the strain on the boundary does not follow the prescribed deformation. That is my problem. What do you mean about prescribed motion to all bodies?Do you mean that I set a motion to each particle?What can it show? thank you very much. Best Regards, Alireza -- 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 #686024]: micro strain-paraview
Question #686024 on Yade changed: https://answers.launchpad.net/yade/+question/686024 Status: Answered => Open Alireza Sadeghi is still having a problem: Hello Bruno, When I increase the number of particles, as I said in the previeus massage, the strain tensor is not correct. In the code, I applied compaction strian rate (negetive) in z direction in the deviatoric part. Therefore, It should have negetive strain in the z direction. But now I have: Triangulated Grains : 12491 sym_grad_u_total_g (wrong averaged strain): 0.0271423 -0.00121828 -0.00169358 -0.00121828 -0.00108962 -5.17863e-05 -0.00169358 -5.17863e-05 0.00136745 Total volume = 0.018195, grad_u = 0.0271423 6.67822e-06 -4.62223e-06 -0.00244324 -0.00108962 0.000223831 -0.00338253 -0.000327403 0.00136745 sym_grad_u (true average strain): 0.0271423 -0.00121828 -0.00169358 -0.00121828 -0.00108962 -5.17863e-05 -0.00169358 -5.17863e-05 0.00136745 Macro strain : 0 0 0 0 0 0 0 0 0 In addition, When I opened the micro strain in Paraview, there is no specific shape. However, when I open view in the controller, the shape of sample is correct. I don't know what should I do. All the conditions are the same, just I increased the number of particles. Thank you very much. best Regards, Alireza -- 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 #686024]: micro strain-paraview
Question #686024 on Yade changed: https://answers.launchpad.net/yade/+question/686024 Status: Answered => Open Alireza Sadeghi is still having a problem: Hello Dear Robert, Thank you very much for your help. Actually the recent code (less particles) works great, and I can open the micro strain in paraview without any problem. But when I increased the number of particles as it was in the question (coke=((1.875e-3,500),(0.9475e-3,1367),(0.50125e-3,10618))), it does not work. My problem is "why the code does not work when I increase the number of particles?" I tried it several times. Without any changes, when I increase the number of particles, I can't reach to micro strain. Thank you very much again. Best Regards, Alireza -- 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 #686024]: micro strain-paraview
Question #686024 on Yade changed:
https://answers.launchpad.net/yade/+question/686024
Status: Needs information => Open
Alireza Sadeghi gave more information on the question:
Dear Robert,
Thank you very much. I reduced the number of particles and now it takes time
around 1 minute.
Sorry for the long code. In other cases, I don't have similar problem.
Thanks a lot.
Best Regards,
Alireza
code#
from yade import utils, plot
from yade import pack, qt
from datetime import datetime
qtr=qt.Renderer()
qtr.bgColor=(1,1,1)
#==
#= define the materials ===
#==
O.materials.append(CohFrictMat(normalCohesion= 1e20, shearCohesion=
1e20, isCohesive= True, young=6.81e8, density=1377.5, poisson=0.3,
frictionAngle= 0.31, fragile=False, label='aggregate-48'))
O.materials.append(CohFrictMat(normalCohesion= 1e20, shearCohesion=
1e20, isCohesive= True, young=6.81e8, density=1532.2, poisson=0.3,
frictionAngle= 0.31, fragile=False, label='aggregate-814'))
O.materials.append(CohFrictMat(normalCohesion= 1e20, shearCohesion= 1e20,
isCohesive= False, young=2e11,
density=1523.6e2, poisson=0.3, frictionAngle= 0.28, fragile=False,
label='wall'))
#===
#=== define walls
#===
walls=aabbWalls([(0,0,0),(2e-1,2e-1,2e-1)],thickness=0.0003,oversizeFactor=1.0,material='wall')
wallIds=O.bodies.append(walls)
#===
#=== define packing
#===
nums=['t','t','t']
mats=['aggregate-48','aggregate-814','aggregate-814']
coke=((1.875e-3,50),(0.9475e-3,136),(0.50125e-3,106))
color=((0,0,1),(0,1,0),(1,1,0))
tolerance=[(2e-4),(1e-4),(1e-5)]
for i in range(len(nums)):
nums[i]=pack.SpherePack()
nums[i].makeCloud((0,0,0),(1e-1,1e-1,1e-1),rMean=coke[i][0],rRelFuzz=tolerance[i],num=coke[i][1])
O.bodies.append([utils.sphere(c,r,material=mats[i],color=color[i]) for c,r
in nums[i]])
#===
#=== define Engine =
#===
O.engines=[
ForceResetter(),
InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Box_Aabb(),Bo1_Facet_Aabb()]),
InteractionLoop([Ig2_Sphere_Sphere_ScGeom6D(),Ig2_Box_Sphere_ScGeom6D(),Ig2_Facet_Sphere_ScGeom6D()],
[Ip2_CohFrictMat_CohFrictMat_CohFrictPhys()],
[Law2_ScGeom_FrictPhys_CundallStrack(),Law2_ScGeom6D_CohFrictPhys_CohesionMoment()]
),
TriaxialStateRecorder(iterPeriod=1000,file='WallStresses'),
#VTKRecorder(fileName='3d-vtk-',recorders=['all'],iterPeriod=1000),
#qt.SnapshotEngine(fileBase='3d-',iterPeriod=200,label='snapshot'),
NewtonIntegrator(damping=0.4,gravity=[0,0,-9.81])
]
yade.qt.Controller(), yade.qt.View()
O.dt=5e-7
#=
#= APPLYING CONFINING PRESSURE
#=
sigmaIso=-1e5
triax=TriaxialStressController(
maxMultiplier=1.000,
finalMaxMultiplier=1.000,
thickness = 0,
stressMask = 7,
internalCompaction=False,
)
O.engines=O.engines+[triax]
triax.goal1=sigmaIso
triax.goal2=sigmaIso
triax.goal3=sigmaIso
triax.wall_back_activated=True
while 1:
O.run(100, True)
unb=unbalancedForce()
meanS=(triax.stress(triax.wall_right_id)[0]+triax.stress(triax.wall_top_id)[1]+triax.stress(triax.wall_front_id)[2])/3
print 'mean Stress:',triax.meanStress,'porosity:', triax.porosity,'meanS:',unb
if triax.porosity<0.48 and abs(sigmaIso-triax.meanStress)/abs(sigmaIso)<0.01:
print 'Isotropic strain1:',-triax.strain[0], 'Isotropic strain
2:',-triax.strain[1], 'Isotropic strain 3:',-triax.strain[2]
break
print "### Isotropic state saved ###"
triax.depth0=triax.depth
triax.height0=triax.height
triax.width0=triax.width
O.save('RVE-sizeDis-solid-Isoe5-Isopart.yade')
#
#= DEVIATORIC LOADING ===
#
triax.stressMask = 3
triax.goal1=sigmaIso
triax.goal2=sigmaIso
triax.goal3=-1
O.saveTmp()
#===
Re: [Yade-users] [Question #686024]: micro strain-paraview
Question #686024 on Yade changed: https://answers.launchpad.net/yade/+question/686024 Status: Needs information => Open Alireza Sadeghi gave more information on the question: Dear Robert, Thank you very much for your reply. When I open the micro strain.vtk files, the paraview give me these errors: vtkVolumeTextureMapper3D was deprecated for VTK 7.0 and will be removed in a future version. vtkOpenGLVolumeTextureMapper3D was deprecated for V7X 7.0 will be removed in a future version. and after it there is nothing in it. Does it mean that the problem is from my Paraview?I am using ParaView 5.4.1. Thank you very much. Best Regards, Alireza -- 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 #686024]: micro strain-paraview
New question #686024 on Yade:
https://answers.launchpad.net/yade/+question/686024
Hello all,
I want to simulate a sample in triaxial test. The sample has particles with
three different sizes. The problem is that when I have only mono-size paricles,
I can see the micro strain in paraview, but when I have three different size
particles, the paraview does not show me the correct result. I was wondering if
you could check my code.
Thank you very much.
Best Regards,
Alireza
### code##
from yade import utils, plot
from yade import pack, qt
from datetime import datetime
#==
#= define the materials ===
#==
O.materials.append(CohFrictMat(normalCohesion= 1e20, shearCohesion= 1e20,
isCohesive= True, young=6.81e8, density=1377.5, poisson=0.3, frictionAngle=
0.31, fragile=False, label='aggregate-48'))
O.materials.append(CohFrictMat(normalCohesion= 1e20, shearCohesion= 1e20,
isCohesive= True, young=6.81e8, density=1532.2, poisson=0.3, frictionAngle=
0.31, fragile=False, label='aggregate-814'))
O.materials.append(CohFrictMat(normalCohesion= 1e20, shearCohesion= 1e20,
isCohesive= False, young=2e11,
density=1523.6e2, poisson=0.3, frictionAngle= 0.28, fragile=False,
label='wall'))
#===
#=== define walls
#===
walls=aabbWalls([(0,0,0),(2e-1,2e-1,2e-1)],thickness=0.0003,oversizeFactor=1.0,material='wall')
wallIds=O.bodies.append(walls)
#===
#=== define packing
#===
nums=['t','t','t']
mats=['aggregate-48','aggregate-814','aggregate-814']
coke=((1.875e-3,500),(0.9475e-3,1367),(0.50125e-3,10618))
color=((0,0,1),(0,1,0),(1,1,0))
tolerance=[(2e-4),(1e-4),(1e-5)]
for i in range(len(nums)):
nums[i]=pack.SpherePack()
nums[i].makeCloud((0,0,0),(2e-1,2e-1,2e-1),rMean=coke[i][0],rRelFuzz=tolerance[i],num=coke[i][1])
O.bodies.append([utils.sphere(c,r,material=mats[i],color=color[i]) for c,r
in nums[i]])
#===
#=== define Engine =
#===
O.engines=[
ForceResetter(),
InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Box_Aabb(),Bo1_Facet_Aabb()]),
InteractionLoop([Ig2_Sphere_Sphere_ScGeom6D(),Ig2_Box_Sphere_ScGeom6D(),Ig2_Facet_Sphere_ScGeom6D()],
[Ip2_CohFrictMat_CohFrictMat_CohFrictPhys()],
[Law2_ScGeom_FrictPhys_CundallStrack(),Law2_ScGeom6D_CohFrictPhys_CohesionMoment()]
),
TriaxialStateRecorder(iterPeriod=1000,file='WallStresses'),
#VTKRecorder(fileName='3d-vtk-',recorders=['all'],iterPeriod=1000),
#qt.SnapshotEngine(fileBase='3d-',iterPeriod=200,label='snapshot'),
NewtonIntegrator(damping=0.4,gravity=[0,0,-9.81])
]
yade.qt.Controller(), yade.qt.View()
O.dt=5e-8
#=
#= APPLYING CONFINING PRESSURE
#=
sigmaIso=-1e5
triax=TriaxialStressController(
maxMultiplier=1.000,
finalMaxMultiplier=1.000,
thickness = 0,
stressMask = 7,
internalCompaction=False,
)
O.engines=O.engines+[triax]
triax.goal1=sigmaIso
triax.goal2=sigmaIso
triax.goal3=sigmaIso
triax.wall_back_activated=True
while 1:
O.run(100, True)
unb=unbalancedForce()
meanS=(triax.stress(triax.wall_right_id)[0]+triax.stress(triax.wall_top_id)[1]+triax.stress(triax.wall_front_id)[2])/3
print 'mean Stress:',triax.meanStress,'porosity:', triax.porosity,'meanS:',unb
if triax.porosity<0.38 and abs(sigmaIso-triax.meanStress)/abs(sigmaIso)<0.01:
print 'Isotropic strain1:',-triax.strain[0], 'Isotropic strain
2:',-triax.strain[1], 'Isotropic strain 3:',-triax.strain[2]
break
print "### Isotropic state saved ###"
triax.depth0=triax.depth
triax.height0=triax.height
triax.width0=triax.width
O.save('RVE-sizeDis-solid-Isoe5-Isopart.yade')
#
#= DEVIATORIC LOADING ===
#
triax.stressMask = 3
triax.goal1=sigmaIso
triax.goal2=sigmaIso
triax.goal3=-1
O.saveTmp()
#
# Micro Strain===
#
TW=TesselationWrapper()
TW.comp
Re: [Yade-users] [Question #684547]: A proper PC for DEM or FEMxDEM simulation
Question #684547 on Yade changed: https://answers.launchpad.net/yade/+question/684547 Status: Answered => Open Alireza Sadeghi is still having a problem: Thank you very much for your response. Is there any paper about the hardware which is proper for DEM simulation?how about parallel calculation and using PC for YADE simulation? Thank you very much. -- 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 #684547]: A proper PC for DEM or FEMxDEM simulation
New question #684547 on Yade: https://answers.launchpad.net/yade/+question/684547 Hello, Do have any suggestions for a proper system? My project is about DEM simulation and I have problem with high run time. My professor asked me if I need faster PC. I have a PC now with this configuration: Intel Core i7-4790@3.6GHz x 8. It is really fast but it is not enough. The price is not important. Thank you for your help. Best Regards Alireza -- 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 #684030]: example of FEMxDEM simulation
Question #684030 on Yade changed: https://answers.launchpad.net/yade/+question/684030 Status: Answered => Solved Alireza Sadeghi confirmed that the question is solved: Thanks Yaniv Fogel, 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 #684030]: example of FEMxDEM simulation
Question #684030 on Yade changed: https://answers.launchpad.net/yade/+question/684030 Status: Answered => Open Alireza Sadeghi is still having a problem: Dear Huanran, Thank you very much for your reply. Is it possible to share the link of the example?I couldn't find it. Thank you again. Best Regards Alireza -- 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 #684030]: example of FEMxDEM simulation
New question #684030 on Yade: https://answers.launchpad.net/yade/+question/684030 Hello All, Is it possible to share the python code of the example of [1]? I am new in multi scale simulation and I want to have a complete example of it. In addition, I install both the YADE and Escript on my PC. I know how to work with both softwares. Am I have to install another software except YADE and Escript? Thank you very much for your help in advance. Best Regards Alireza [1] https://yade-dem.org/doc/FEMxDEM.html#example-tests -- 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 #683375]: uncertainty in the shear behavior of material
Question #683375 on Yade changed: https://answers.launchpad.net/yade/+question/683375 Status: Answered => Solved Alireza Sadeghi confirmed that the question is solved: Thanks Robert Caulk, 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 #683375]: uncertainty in the shear behavior of material
New question #683375 on Yade:
https://answers.launchpad.net/yade/+question/683375
Hello All,
I was encountered to a problem in my dem modeling. I want to simulate the bulk
behavior of material. Due to reach this goal, I have to choose a number of
particles which could represent the bulk behavior as similar as experiment. I
did two type of simulation.
in the first one, I put the particles in the box and the the box faces compress
the paricles, after in the x and y direction, stress remain constant and in the
z direction, we increase stress and plot the strain-stress in z direction. Due
to initial conditions (randomness of initial position of particles) we will
reach to different results in each realization. The code is as below:
from yade import utils, plot
from yade import pack, qt
from datetime import datetime
from yade import export
O.periodic=True
O.cell.refSize=(2e-1,2e-1,2e-1)
O.materials.append(CohFrictMat(normalCohesion= 1e20, shearCohesion= 1e20,
isCohesive= True, young=6.81e8, density=1377.5, poisson=0.3, frictionAngle=
0.31, fragile=False, label='Coke'))
nums=['t']
mats=['Coke']
coke=(5e-3,3000)
nums=pack.SpherePack()
nums.makeCloud((0,0,0),(2e-1,2e-1,2e-1),rMean=coke[0],rRelFuzz=0,num=coke[1])
O.bodies.append([utils.sphere(c,r,material=mats[0],color=(0,0,1)) for c,r in
nums])
sigmaIso=-1e5
O.engines=[
ForceResetter(),
InsertionSortCollider([Bo1_Sphere_Aabb()]),
InteractionLoop(
[Ig2_Sphere_Sphere_ScGeom()],
[Ip2_FrictMat_FrictMat_FrictPhys()],
[Law2_ScGeom_FrictPhys_CundallStrack()]
),
PeriTriaxController(label='triax',
goal=(sigmaIso,sigmaIso,sigmaIso),stressMask=7,
dynCell=True,maxStrainRate=(1,1,1),
maxUnbalanced=.01,relStressTol=1e-3,
doneHook='compactionFinished()'
),
NewtonIntegrator(damping=0.4),
PyRunner(command='history()',iterPeriod=1000,label='recorder'),
]
O.dt=1e-6
O.saveTmp()
def compactionFinished():
O.cell.trsf=Matrix3.Identity
triax.goal=(sigmaIso,sigmaIso,-0.5)
triax.stressMask=3
triax.maxStrainRate=(0.1,0.1,0.1)
triax.doneHook='triaxFinished()'
triax.maxUnbalanced=1
dataOld = plot.data
plot.saveDataTxt('RVE3000-iso_part-e5-5.txt')
plot.data = {}
plot.plot()
def history():
plot.addData(unbalanced=unbalancedForce(),
sxx=-triax.stress[0],syy=-triax.stress[1],szz=-triax.stress[2],
exx=-triax.strain[0],eyy=-triax.strain[1],ezz=-triax.strain[2],
sxz=-0.5*(triax.stress[2]+triax.stress[0]),
ev=(triax.strain[2]+triax.strain[1]+triax.strain[0]),
q=-(triax.stress[2]-triax.stress[1]),
p=-(triax.stress[2]+triax.stress[1]+triax.stress[0])/3 ,
D=(triax.stress[2]-triax.stress[1])/sigmaIso,
R=3*(triax.stress[2]-triax.stress[1])/(triax.stress[2]+triax.stress[1]+triax.stress[0]+1e-2),
por=porosity(),i=O.iter,)
print 'porosity:', porosity()
print 'iteration', O.iter
O.run(800,True)
plot.saveDataTxt('RVE3000-dev_part-e5-5.txt')
plot.plots={'ezz':('szz')}
plot.plot()
def triaxFinished():
print 'Compaction finished'
O.pause()
In the second one, I made the sample first with the code below:
from yade import utils, plot
from yade import pack, qt
from datetime import datetime
from yade import export
O.periodic=True
O.cell.refSize=(2e-1,2e-1,2e-1)
O.materials.append(CohFrictMat(normalCohesion= 1e20, shearCohesion= 1e20,
isCohesive= True, young=6.81e8, density=1377.5, poisson=0.3, frictionAngle=
0.31, fragile=False, label='Coke'))
nums=['t']
mats=['Coke']
coke=(5e-3,3000)
nums=pack.SpherePack()
nums.makeCloud((0,0,0),(2e-1,2e-1,2e-1),rMean=coke[0],rRelFuzz=0,num=coke[1])
O.bodies.append([utils.sphere(c,r,material=mats[0],color=(0,0,1)) for c,r in
nums])
sigmaIso=-1e5
O.engines=[
ForceResetter(),
InsertionSortCollider([Bo1_Sphere_Aabb()]),
InteractionLoop(
[Ig2_Sphere_Sphere_ScGeom()],
[Ip2_FrictMat_FrictMat_FrictPhys()],
[Law2_ScGeom_FrictPhys_CundallStrack()]
),
PeriTriaxController(label='triax',
goal=(sigmaIso,sigmaIso,sigmaIso),stressMask=7,
dynCell=True,maxStrainRate=(1,1,1),
maxUnbalanced=.01,relStressTol=1e-3,
doneHook='compactionFinished()'
),
NewtonIntegrator(damping=0.4),
PyRunner(command='history()',iterPeriod=1000,label='recorder'),
]
O.dt=1e-6
O.saveTmp()
def compactionFinished():
print 'Finished'
O.pause()
plot.saveDataTxt('Uncertainty-Pre-compact3000.txt')
plot.data = {}
O.save('Pre-compact3000.yade')
O.run(150,True)
and after that, I compress the sample of "Pre-compact3000.yade" in the z
direction and plot the stress-strain in the z direction with code below:
from yade import utils, plot
from yade
Re: [Yade-users] [Question #681474]: force chain in paraview
Question #681474 on Yade changed: https://answers.launchpad.net/yade/+question/681474 Status: Needs information => Solved Alireza Sadeghi confirmed that the question is solved: YES, Thank you Jan, Thank you very much. It solved. -- 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 #681474]: force chain in paraview
Question #681474 on Yade changed:
https://answers.launchpad.net/yade/+question/681474
Status: Needs information => Open
Alireza Sadeghi gave more information on the question:
Hello Jan,
Thank you very much for your response. I think the vtk export express just the
initial interaction, because at the beginning there is no interaction between
the particles.
This is my complete code:
from yade import utils, plot
from yade import pack, qt
from datetime import datetime
from yade import export
qtr=qt.Renderer()
qtr.bgColor=(1,1,1)
#===
#==set a periodic boundary==
#===
O.periodic=True
O.cell.refSize=(2e-1,2e-1,2e-1)
#==
#= define the materials ===
#==
O.materials.append(CohFrictMat(normalCohesion= 1e20, shearCohesion=
1e20, isCohesive= True, young=6.81e8, density=1377.5, poisson=0.3,
frictionAngle= 0.31, fragile=False, label='Coke'))
#===
#=== define packing
#===
nums=['t']
mats=['Coke']
coke=(5e-3,150)
nums=pack.SpherePack()
nums.makeCloud((0,0,0),(2e-1,2e-1,2e-1),rMean=coke[0],rRelFuzz=0,num=coke[1])
O.bodies.append([utils.sphere(c,r,material=mats[0],color=(0,0,1)) for
c,r in nums])
sigmaIso=-1e8
#===
#=VTK===
#===
vtkExporter = export.VTKExporter('VTK_simulatione')
#===
#=== define Engine =
#===
O.engines=[
ForceResetter(),
InsertionSortCollider([Bo1_Sphere_Aabb()]),
InteractionLoop(
[Ig2_Sphere_Sphere_ScGeom()],
[Ip2_FrictMat_FrictMat_FrictPhys()],
[Law2_ScGeom_FrictPhys_CundallStrack()]
),
PeriTriaxController(label='triax',
goal=(sigmaIso,sigmaIso,sigmaIso),stressMask=7,
dynCell=True,maxStrainRate=(1,1,1),
maxUnbalanced=.01,relStressTol=1e-3,
doneHook='compactionFinished()'
),
NewtonIntegrator(damping=0.4),
PyRunner(iterPeriod=100,
command="vtkExporter.exportInteractions(ids='all',what=[('forceN','i.phys.normalForce.norm()')])"),
PyRunner(command='history()',iterPeriod=10,label='recorder'),
]
O.dt=1e-6
O.saveTmp()
def history():
plot.addData(unbalanced=unbalancedForce(),
szz=-triax.stress[2],ezz=-triax.strain[2],i=O.iter,)
print 'porosity:', porosity()
print 'iteration', O.iter
O.run()
plot.plots={'ezz':('szz')}
plot.plot()
def compactionFinished():
print 'Finished'
O.pause()
--
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 #681474]: force chain in paraview
Question #681474 on Yade changed:
https://answers.launchpad.net/yade/+question/681474
Alireza Sadeghi gave more information on the question:
Dear Jerome and Dear Jan,
Thank you very much for your help. I increased the confining pressure to
-1e8 Pa. the particles go through each other. but in the ParaView there
is no interaction. in addition, I added
PyRunner(command='history()',iterPeriod=10,label='recorder'),
to my engine and added
def history():
plot.addData(szz2=-triax.stress[2], ezz=-triax.strain[2], i=O.iter,)
O.run()
plot.plots={'ezz':('szz2')}
plot.plot()
to my code, and the diagram shows the stress versus strain and it is not zero.
All in all, my problem does bot solve.
Thank you very much for your help.
Best Regards
Alireza
--
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 #681474]: force chain in paraview
Question #681474 on Yade changed: https://answers.launchpad.net/yade/+question/681474 Status: Answered => Open Alireza Sadeghi is still having a problem: Dear Jan, Thank you very much. You comment works great. But my question doesn't solve completely. Now, the program run successfully, but I have a warning "Warning (yade.export.VTKExporter): no interactions" The particles hit to each other and compaction process is completed. but the VTKExporter says that no interactions. I don't know what should I do now. Thank you very much for your help. Best Regards Alireza -- 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 #681474]: force chain in paraview
Question #681474 on Yade changed: https://answers.launchpad.net/yade/+question/681474 Status: Answered => Open Alireza Sadeghi is still having a problem: Hello Jerome, Thank you very much for your response. I added the sentences which you told. Now, I have this error "TypeError: 'Vector3' object is not callable". The program gave me a file with name "VTK_simulatione- intrs-.vtk". I don't know where did the zeros come from. However, when I deleted the zeroes I could open the file with Paraview. In the Paraview, I can't select Cell Data to Point Data. In addition, all the tube has the same radius and I can't set the color by forceN. I was wondering if you could help me. If there is an example which is similar to my problem, I will be grateful if you could share it with me. Thank you very much Alireza -- 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 #681474]: force chain in paraview
New question #681474 on Yade:
https://answers.launchpad.net/yade/+question/681474
Hello All,
I want to model force chain in the paraview, but I have problem to do it. I
have the error that "NameError: name 'VTKExporter' is not defined". I was
wondering if you could give me an example to solve my problem. Thank you very
much.
Best Regards
Alireza
P.S. my code is :
from yade import utils, plot
from yade import pack, qt
from datetime import datetime
qtr=qt.Renderer()
qtr.bgColor=(1,1,1)
#===
#==set a periodic boundary==
#===
O.periodic=True
O.cell.refSize=(2e-1,2e-1,2e-1)
#==
#= define the materials ===
#==
O.materials.append(CohFrictMat(normalCohesion= 1e20, shearCohesion= 1e20,
isCohesive= True, young=6.81e8, density=1377.5, poisson=0.3, frictionAngle=
0.31, fragile=False, label='Coke'))
#===
#=== define packing
#===
nums=['t']
mats=['Coke']
coke=(5e-3,150)
nums=pack.SpherePack()
nums.makeCloud((0,0,0),(2e-1,2e-1,2e-1),rMean=coke[0],rRelFuzz=0,num=coke[1])
O.bodies.append([utils.sphere(c,r,material=mats[0],color=(0,0,1)) for c,r in
nums])
sigmaIso=-1e5
#===
#=VTK===
#===
vtkExporter = VTKExporter('VTK_simulatione')
#===
#=== define Engine =
#===
O.engines=[
ForceResetter(),
InsertionSortCollider([Bo1_Sphere_Aabb()]),
InteractionLoop(
[Ig2_Sphere_Sphere_ScGeom()],
[Ip2_FrictMat_FrictMat_FrictPhys()],
[Law2_ScGeom_FrictPhys_CundallStrack()]
),
PeriTriaxController(label='triax',
goal=(sigmaIso,sigmaIso,sigmaIso),stressMask=7,
dynCell=True,maxStrainRate=(1,1,1),
maxUnbalanced=.01,relStressTol=1e-3,
doneHook='compactionFinished()'
),
NewtonIntegrator(damping=0.4),
PyRunner(iterPeriod=50,
command="vtkExporter.exportInteractions(ids='all',what=[('forceN','i.phys.normalForce().norm()')])")
]
O.dt=1e-6
O.saveTmp()
O.run(20,True)
def compactionFinished():
print 'Finished'
O.pause()
--
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 #680144]: plot the deviatoric part of TriaxialStressController simulation
Question #680144 on Yade changed: https://answers.launchpad.net/yade/+question/680144 Alireza Sadeghi 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 #680144]: plot the deviatoric part of TriaxialStressController simulation
Question #680144 on Yade changed: https://answers.launchpad.net/yade/+question/680144 Status: Answered => Solved Alireza Sadeghi confirmed that the question is solved: Thank you dear Jerome, This solved my problem. Best Regards Alireza -- 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 #680144]: plot the deviatoric part of TriaxialStressController simulation
Question #680144 on Yade changed:
https://answers.launchpad.net/yade/+question/680144
Status: Answered => Open
Alireza Sadeghi is still having a problem:
Hello Jerome,
Thank you for your response and help. I will use
https://gitlab.com/yade-dev/trunk. Thank you.
about my question, I used plot.resetData() but it did not solve my problem.
The simulation of triaxial test has two parts. The particles compact with
isometric pressure and in the second step we have deviatoric compaction. The
problem is, strain in the first step is added to the strain in the second step.
but I want to draw strain-Stress diagram just for the second step of triaxial
test simulation. my code is in the below of the mail. Thank you very much for
your help.
Best Regards
Alireza
#Triaxial compresion for finding a proper size for RVE (stress control
boundary condition)
from yade import utils, plot
from yade import pack, qt
from datetime import datetime
qtr=qt.Renderer()
qtr.bgColor=(1,1,1)
#==
#= define the materials ===
#==
O.materials.append(CohFrictMat(normalCohesion= 1e20, shearCohesion=
1e20, isCohesive= True, young=6.81e8, density=1377.5, poisson=0.3,
frictionAngle= 0.31, fragile=False, label='Coke'))
O.materials.append(CohFrictMat(normalCohesion= 1e20, shearCohesion= 1e20,
isCohesive= False, young=4e9,
density=1523.6, poisson=0.3, frictionAngle= 0.0, fragile=False, label='wall'))
#===
#=== define walls
#===
walls=aabbWalls([(0,0,0),(2e-1,2e-1,2e-1)],thickness=0.0003,oversizeFactor=1.0,material='wall')
wallIds=O.bodies.append(walls)
#===
#=== define packing
#===
nums=['t']
mats=['Coke']
coke=(5e-3,2000)
nums=pack.SpherePack()
nums.makeCloud((0,0,0),(2e-1,2e-1,2e-1),rMean=coke[0],rRelFuzz=1e-4,num=coke[1])
O.bodies.append([utils.sphere(c,r,material=mats[0],color=(0,0,1)) for
c,r in nums])
#===
#=== define Engine =
#===
O.engines=[
ForceResetter(),
InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Box_Aabb(),Bo1_Facet_Aabb()]),
InteractionLoop([Ig2_Sphere_Sphere_ScGeom6D(),Ig2_Box_Sphere_ScGeom6D(),Ig2_Facet_Sphere_ScGeom6D()],
[Ip2_CohFrictMat_CohFrictMat_CohFrictPhys()],
[Law2_ScGeom_FrictPhys_CundallStrack(),Law2_ScGeom6D_CohFrictPhys_CohesionMoment()]
),
# TriaxialStateRecorder(iterPeriod=100,file='WallStresses'+table.key),
#VTKRecorder(fileName='3d-vtk-',recorders=['all'],iterPeriod=100),
#qt.SnapshotEngine(fileBase='3d-',iterPeriod=200,label='snapshot'),
NewtonIntegrator(damping=0.4,gravity=[0,0,0])
]
O.dt=2e-6
#=
#= APPLYING CONFINING PRESSURE
#=
stabilityThreshold=0.2
sigmaIso=-1e5
triax=TriaxialStressController(
maxMultiplier=1.000,
finalMaxMultiplier=1.000,
thickness = 0,
stressMask = 7,
internalCompaction=False,
)
O.engines=O.engines+[triax]
triax.goal1=sigmaIso
triax.goal2=sigmaIso
triax.goal3=sigmaIso
triax.wall_back_activated=True
while 1:
O.run(100, True)
unb=unbalancedForce()
meanS=(triax.stress(triax.wall_right_id)[0]+triax.stress(triax.wall_top_id)[1]+triax.stress(triax.wall_front_id)[2])/3
print 'mean Stress:',triax.meanStress,'porosity:', triax.porosity,'meanS:',unb
if triax.porosity < 0.497 and ((sigmaIso-triax.meanStress)/sigmaIso) < 0.001:
print 'Isotropic strain1:',-triax.strain[0], 'Isotropic strain
2:',-triax.strain[1], 'Isotropic strain 3:',-triax.strain[2]
break
print "### Isotropic state saved ###"
#
#= DEVIATORIC LOADING ===
#
plot.reset()
triax.stressMask = 3
triax.goal1=sigmaIso
triax.goal2=sigmaIso
triax.goal3=-10
#dataOld = plot.data
#plot.saveDataTxt('isotropic_part_deformableRVE.txt')
plot.resetData()
#plot.plot()
O.saveTmp()
#=
#=== data c
[Yade-users] [Question #680144]: plot the deviatoric part of TriaxialStressController simulation
New question #680144 on Yade: https://answers.launchpad.net/yade/+question/680144 Hello All, I am using TriaxialStressController to simulate a triaxial test. I want to plot stress in the z-direction vs strain in the z-direction. my code is similar to [1], except without growing of particles. The strain which is plotted in the code has isotropic compaction part too. I was wondering if I could remove an isotropic compaction part and reset the strain and plot just deviatoric part. Thank you very much for your help in advance. Best Regards Alireza [1] http://bazaar.launchpad.net/~yade-pkg/yade/git-trunk/view/head:/examples/triax-tutorial/script-session1.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 #679661]: plot diagram in PeriTriaxController
Question #679661 on Yade changed: https://answers.launchpad.net/yade/+question/679661 Status: Answered => Solved Alireza Sadeghi confirmed that the question is solved: Dear Jan, Thank you for your help. Your help works great. Best Regards Alireza -- 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 #679661]: plot diagram in PeriTriaxController
New question #679661 on Yade:
https://answers.launchpad.net/yade/+question/679661
Hello all,
I encountered to a problem in plotting diagram in PeriTriaxController. I want
to simulate a triaxial test with below code. The simulation consist of tow
part. 1. isotropic compaction and 2. deviatoric compaction. but when I plotted
strain vs stress in compaction direction (szz-ezz), it plotted both isotropic
and deviatoric part, but I want just the deviatoric part. Could you please help
me to find how can I plot the diagram just for deviatoric part?
Thank you very much for your help in advance.
Best regards
Alireza
P.S. the code is:
from yade import utils, plot
from yade import pack, qt
from datetime import datetime
qtr=qt.Renderer()
qtr.bgColor=(1,1,1)
#===
#==set a periodic boundary==
#===
O.periodic=True
O.cell.refSize=(2e-1,2e-1,2e-1)
#==
#= define the materials ===
#==
O.materials.append(CohFrictMat(normalCohesion= 1e20, shearCohesion= 1e20,
isCohesive= True, young=6.81e8, density=1377.5e2, poisson=0.3, frictionAngle=
0.31, fragile=False, label='Coke'))
#===
#=== define packing
#===
nums=['t']
mats=['Coke']
coke=(5e-3,200)
nums=pack.SpherePack()
nums.makeCloud((0,0,0),(2e-1,2e-1,2e-1),rMean=coke[0],rRelFuzz=1e-4,num=coke[1])
O.bodies.append([utils.sphere(c,r,material=mats[0],color=(0,0,1)) for c,r in
nums])
#===
#=== define Engine =
#===
sigmaIso=-1e6
O.engines=[
ForceResetter(),
InsertionSortCollider([Bo1_Sphere_Aabb()]),
InteractionLoop(
[Ig2_Sphere_Sphere_ScGeom()],
[Ip2_FrictMat_FrictMat_FrictPhys()],
[Law2_ScGeom_FrictPhys_CundallStrack()]
),
PeriTriaxController(label='triax',
goal=(sigmaIso,sigmaIso,sigmaIso),stressMask=7,
dynCell=True,maxStrainRate=(1,1,1),
maxUnbalanced=.1,relStressTol=1e-3,
doneHook='compactionFinished()'
),
NewtonIntegrator(damping=0.4),
PyRunner(command='history()',iterPeriod=10,label='recorder'),
]
O.dt=0.5*PWaveTimeStep()
O.saveTmp()
def history():
plot.addData(unbalanced=unbalancedForce(),
sxx=triax.stress[0],syy=triax.stress[1],szz=-triax.stress[2],
exx=triax.strain[0],eyy=triax.strain[1],ezz=-triax.strain[2],
sxz=-0.5*(triax.stress[2]+triax.stress[0]),
ev=-(triax.strain[2]-triax.strain[0]),
q=-(triax.stress[2]-triax.stress[1]),
p=-(triax.stress[2]+triax.stress[1]+triax.stress[0])/3 ,
R=3*(triax.stress[2]-triax.stress[1])/(triax.stress[2]+triax.stress[1]+triax.stress[0]+1e-2),
por=porosity(),i=O.iter,)
print 'stress difference:', -(triax.stress[2]-triax.stress[1])
print 'mean
stress:',-(triax.stress[2]+triax.stress[1]+triax.stress[0])/3
print 'porosity:', porosity()
def compactionFinished():
O.cell.trsf=Matrix3.Identity
triax.goal=(sigmaIso,sigmaIso,-0.2)
triax.stressMask=3
triax.maxStrainRate=(.1,.1,.1)
triax.doneHook='triaxFinished()'
triax.maxUnbalanced=10
O.run(40,True)
plot.plots={'ezz':('szz')}
plot.plot()
def triaxFinished():
O.run(3,True)
print 'Finished'
O.pause()
--
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 #679585]: Using CPU cores for a simulation
Question #679585 on Yade changed: https://answers.launchpad.net/yade/+question/679585 Status: Answered => Solved Alireza Sadeghi confirmed that the question is solved: Hello dear Jerome, Thank you very much for your help. Best Regards Alireza -- 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 #679585]: Using CPU cores for a simulation
New question #679585 on Yade: https://answers.launchpad.net/yade/+question/679585 Hello I want to use all of the CPU cores which I have for my simulation, but I don't know how can I do this?Should I write something in the python code?or in the run section, should I add something? Thank you very much for your help in advance. Bests, Alireza -- 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 #678982]: about TriaxialStressController
Question #678982 on Yade changed: https://answers.launchpad.net/yade/+question/678982 Status: Answered => Solved Alireza Sadeghi confirmed that the question is solved: Dear Robert, Thank you very much for your help. It solved my problem. Best Regards Alireza -- 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 #678982]: about TriaxialStressController
Question #678982 on Yade changed: https://answers.launchpad.net/yade/+question/678982 Status: Answered => Open Alireza Sadeghi is still having a problem: Dear Robert, Thank you very much for your help. I used the sample, but the code didn't stop. I think the threshold for instability force was too small. I run the below code, and in the second section of it (deviatoric part) it doesn't stop when it reach to -0.2 strain in z direction. What can I do for it? Thank you very much for your help. Best Regards Alireza #=== from yade import utils, plot from yade import pack, qt from datetime import datetime qtr=qt.Renderer() qtr.bgColor=(1,1,1) stabilityThreshold=0.01 #= define the materials === O.materials.append(CohFrictMat(normalCohesion= 1e20, shearCohesion= 1e20, isCohesive= True, young=6.81e8, density=1377.5e2, poisson=0.3, frictionAngle= 0.31, fragile=False, label='Coke')) O.materials.append(CohFrictMat(normalCohesion= 1e20, shearCohesion= 1e20, isCohesive= False, young=4e9, density=1523.6, poisson=0.3, frictionAngle= 0.0, fragile=False, label='wall')) #=== define walls == walls=aabbWalls([(0,0,0),(2e-1,2e-1,2e-1)],thickness=0.0003,oversizeFactor=1.0,material='wall') wallIds=O.bodies.append(walls) #=== define packing nums=['t'] mats=['Coke'] coke=(5e-3,1000) nums=pack.SpherePack() nums.makeCloud((0,0,0),(2e-1,2e-1,2e-1),rMean=coke[0],rRelFuzz=1e-4,num=coke[1]) O.bodies.append([utils.sphere(c,r,material=mats[0],color=(0,0,1)) for c,r in nums]) #=== define Engine = O.engines=[ ForceResetter(), InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Box_Aabb(),Bo1_Facet_Aabb()]), InteractionLoop([Ig2_Sphere_Sphere_ScGeom6D(),Ig2_Box_Sphere_ScGeom6D(),Ig2_Facet_Sphere_ScGeom6D()], [Ip2_CohFrictMat_CohFrictMat_CohFrictPhys()], [Law2_ScGeom_FrictPhys_CundallStrack(),Law2_ScGeom6D_CohFrictPhys_CohesionMoment()] ), NewtonIntegrator(damping=0.4,gravity=[0,0,0]) ] O.dt=2e-6 #= APPLYING CONFINING PRESSURE triax=TriaxialStressController( maxMultiplier=1.000, finalMaxMultiplier=1.000, thickness = 0, stressMask = 7, internalCompaction=False, ) O.engines=O.engines+[triax] triax.goal1=-1.0e5 triax.goal2=-1.0e5 triax.goal3=-1.0e5 triax.wall_back_activated=True while 1: O.run(1000, True) unb=unbalancedForce() print 'unbalanced force:',unb,' mean stress: ',triax.meanStress if unbhttps://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 #678982]: about TriaxialStressController
Question #678982 on Yade changed: https://answers.launchpad.net/yade/+question/678982 Status: Answered => Open Alireza Sadeghi is still having a problem: Hello Jérôme, Thank you for your help. I saw the code but it couldn't solve my problem. I want to do a simulation of triaxial test similar to periodic triaxial test [1] but its boundary is not priodic. Thank you for your help. Best Regards Alireza [1] https://yade-dev.gitlab.io/trunk/tutorial-examples.html#periodic-triaxial-test -- 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 #678982]: about TriaxialStressController
New question #678982 on Yade: https://answers.launchpad.net/yade/+question/678982 Hello all, I have a problem. I want to compress particles with an isotropic compaction (for example with 10 kpa compression in the x, y, and z direction) and after that put the stress control in x and y directions (for example 10 kpa) and put strain (for example -0.2) in the z direction. I don't know how should I do this simulation with TriaxialStressController. Could you please give me a hint? A sample of my code is: #=== from yade import utils, plot from yade import pack, qt from datetime import datetime O.materials.append(CohFrictMat(normalCohesion= 1e20, shearCohesion= 1e20, isCohesive= True, young=1.95e7, density=1532.2, poisson=0.3, frictionAngle= 0.0, fragile=False, label='Coke')) O.materials.append(CohFrictMat(normalCohesion= 1e20, shearCohesion= 1e20, isCohesive= False, young=4e9, density=1523.6, poisson=0.3, frictionAngle= 0.0, fragile=False, label='wall')) walls=aabbWalls([(0,0,0),(2e-1,2e-1,2e-1)],thickness=0.0003,oversizeFactor=1.0,material='wall') wallIds=O.bodies.append(walls) nums=['t'] mats=['Coke'] coke=(5e-3,1000) nums=pack.SpherePack() nums.makeCloud((0,0,0),(2e-1,2e-1,2e-1),rMean=coke[0],rRelFuzz=1e-4,num=coke[1]) O.bodies.append([utils.sphere(c,r,material=mats[0],color=(0,0,1)) for c,r in nums]) O.engines=[ ForceResetter(), InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Box_Aabb(),Bo1_Facet_Aabb()]), InteractionLoop([Ig2_Sphere_Sphere_ScGeom6D(),Ig2_Box_Sphere_ScGeom6D(),Ig2_Facet_Sphere_ScGeom6D()], [Ip2_CohFrictMat_CohFrictMat_CohFrictPhys()], [Law2_ScGeom_FrictPhys_CundallStrack(),Law2_ScGeom6D_CohFrictPhys_CohesionMoment()] ), NewtonIntegrator(damping=0.4), PyRunner(command='calm()',iterPeriod=10,label='calmEngine') ] triax=TriaxialStressController( maxMultiplier=1.000, finalMaxMultiplier=1.000, thickness = 0, stressMask = 7, internalCompaction=False, ) O.engines=O.engines+[triax] triax.goal1=-1.0e5 triax.goal2=-1.0e5 triax.goal3=-1.0e5 triax.wall_back_activated=True O.dt=2e-6 calmEngine.dead=True triax2=TriaxialStressController( maxMultiplier=1.000, finalMaxMultiplier=1.000, thickness = 0, stressMask = 3, internalCompaction=False, ) O.engines=O.engines+[triax2] triax2.goal1=-1.0e5 triax2.goal2=-1.0e5 triax2.goal3=-0.2 triax2.wall_back_activated=True # Thank you very much for your help in advance. Best Regards Alireza -- 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 #678796]: Periodic triaxial test
Question #678796 on Yade changed: https://answers.launchpad.net/yade/+question/678796 Status: Answered => Solved Alireza Sadeghi confirmed that the question is solved: Thanks Robert Caulk, 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 #678796]: Periodic triaxial test
Question #678796 on Yade changed: https://answers.launchpad.net/yade/+question/678796 Status: Solved => Open Alireza Sadeghi is still having a problem: Dear Robert, Thank you very much for your help. My script does not work and it can not plot any thing. I don't know how should I fix it. Thank you for your help. Best Regards Alireza -- 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 #678796]: Periodic triaxial test
Question #678796 on Yade changed: https://answers.launchpad.net/yade/+question/678796 Status: Answered => Solved Alireza Sadeghi confirmed that the question is solved: Thanks Robert Caulk, 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 #678796]: Periodic triaxial test
Question #678796 on Yade changed: https://answers.launchpad.net/yade/+question/678796 Status: Needs information => Open Alireza Sadeghi gave more information on the question: Dear Robert, Thank you very much for your help. I will obey the law (one question with one object). Thanks for your help. for [1], is there any example for it? for question number 3, it doesn't have any lines in the diagram. the diagram is empty. without any error. I run the sample which I put in this question, but I could not draw the diagram. Thank you for your help. Best Regards Alireza -- 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 #678796]: Periodic triaxial test
New question #678796 on Yade:
https://answers.launchpad.net/yade/+question/678796
hello all,
1. How can I control strain in the "PeriTriaxController"?
2.Is it possible to add shear on the boundary condition with this engine?
3. another question is that I can not draw diagram of Stress in compaction
direction vs strain in compaction direction. what is the problem in your
opinion?
my code is:
from yade import utils, plot
from yade import pack, qt
from datetime import datetime
O.periodic=True
O.cell.refSize=(2e-1,2e-1,2e-1)
O.materials.append(CohFrictMat(normalCohesion= 1e20, shearCohesion= 1e20,
isCohesive= True, young=6.81e8, density=1377.5, poisson=0.3, frictionAngle=
0.31, fragile=False, label='Coke'))
nums=['t']
mats=['Coke']
coke=(5e-3,1000)
nums=pack.SpherePack()
nums.makeCloud((0,0,0),(2e-1,2e-1,2e-1),rMean=coke[0],rRelFuzz=1e-4,num=coke[1])
O.bodies.append([utils.sphere(c,r,material=mats[0]) for c,r in nums])
sigmaIso=-1e5
O.engines=[
ForceResetter(),
InsertionSortCollider([Bo1_Sphere_Aabb()]),
InteractionLoop(
[Ig2_Sphere_Sphere_ScGeom()],
[Ip2_FrictMat_FrictMat_FrictPhys()],
[Law2_ScGeom_FrictPhys_CundallStrack()]
),
PeriTriaxController(label='triax',
goal=(sigmaIso,sigmaIso,sigmaIso),stressMask=7,
dynCell=True,maxStrainRate=(10,10,10),
maxUnbalanced=.1,relStressTol=1e-3,
doneHook='compactionFinished()'
),
NewtonIntegrator(damping=0.4),
PyRunner(command='addPlotData()',iterPeriod=100),
]
O.dt=0.5*PWaveTimeStep()
def addPlotData():
plot.addData(unbalanced=unbalancedForce(),i=O.iter,
sxx=triax.stress[0],syy=triax.stress[1],szz=triax.stress[2],
exx=triax.strain[0],eyy=triax.strain[1],ezz=triax.strain[2])
#print 'normal stress', triax.stress[2]
#print 'normal strain', triax.strain[2]
plot.plots={('ezz'):('szz')}
plot.plot()
def compactionFinished():
O.cell.trsf=Matrix3.Identity
triax.goal=(sigmaIso,sigmaIso,-.2)
triax.stressMask=3
triax.maxStrainRate=(10,10,10)
triax.doneHook='triaxFinished()'
triax.maxUnbalanced=10
def triaxFinished():
print 'Finished'
O.pause()
The code does not draw "plot.plots={('ezz'):('szz')}". I was wondering if you
could help me to fix this problem. Thank you very much for your help in advance.
Best Regards
Alireza
--
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 #678645]: Shrinkage of particles
Question #678645 on Yade changed: https://answers.launchpad.net/yade/+question/678645 Status: Answered => Solved Alireza Sadeghi confirmed that the question is solved: Dear Jan and Dear Bruno, Thank you very much for your help. My problem is solved with your comments. Best Regards Alireza -- 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 #678645]: Shrinkage of particles
Question #678645 on Yade changed: https://answers.launchpad.net/yade/+question/678645 Status: Needs information => Open Alireza Sadeghi gave more information on the question: Hello Jan, This is a simple case similar to my simulation. O.reset() from yade import utils, plot from yade import pack, qt from datetime import datetime qtr=qt.Renderer() qtr.bgColor=(1,1,1) #== #= define the materials === #== O.materials.append(CohFrictMat(normalCohesion= 1e20, shearCohesion= 1e20, isCohesive= True, young=6.81e8, density=1377.5, poisson=0.3, frictionAngle= 0.31, fragile=False, label='Coke')) O.materials.append(CohFrictMat(normalCohesion= 1e20, shearCohesion= 1e20, isCohesive= False, young=6.81e8, density=1523.6, poisson=0.3, frictionAngle= 0.28, fragile=False, label='mold')) #=== #=== define wall #=== mld=yade.geom.facetCylinder((0,0,0),0.6e-2, 10e-2, orientation=Quaternion((0, 0, 1), 0), segmentsNumber=30, wallMask=6, angleRange=None, closeGap=False, radiusTopInner=-1, radiusBottomInner=-1, wire=True, material='mold') O.bodies.append(mld) #=== #=== define packing #=== nums=['t'] mats=['Coke'] coke=(1e-3,2000) nums=pack.SpherePack() nums.makeCloud((-4.5e-3,-4.5e-3,-5e-2),(4.5e-3,4.5e-3,0.08),rMean=coke[0],rRelFuzz=1e-4,num=coke[1]) O.bodies.append([utils.sphere(c,r,material=mats[0]) for c,r in nums]) #=== #=== define Engine = #=== O.engines=[ ForceResetter(), InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb(),Bo1_Box_Aabb()]), InteractionLoop( [Ig2_Sphere_Sphere_ScGeom(),Ig2_Facet_Sphere_ScGeom(),Ig2_Box_Sphere_ScGeom()], [Ip2_FrictMat_FrictMat_FrictPhys()], [Law2_ScGeom_FrictPhys_CundallStrack()] ), NewtonIntegrator(gravity=(0,0,-9.81),damping=0.4), ] O.dt=1.00e-6 #= #=define Shrinkage O.engines=O.engines+[PyRunner(command='Shrinkage()', label='Shrink',iterPeriod=70)] def Shrinkage(): for i in O.bodies[Sphere]: if i.state.pos[2]<-40e-3: utils.growParticles(0.5,updateMass=True) #for i in O.bodies: # if isinstance(i.shape,Sphere) and i.state.pos[2]<-20e-3 and i.state.pos[2]>-30e-3: #r1=i.shape.radius # i.shape.radius=(1-2.3e-2)*r1 # i.material=O.materials['Coke'] Thank you very much for your help. Best Regards Alireza -- 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 #678646]: tracking a particle
Question #678646 on Yade changed: https://answers.launchpad.net/yade/+question/678646 Alireza Sadeghi posted a new comment: Dear Vasileios and Dear Mina, Thank you very much for your help. Best Regards Alireza -- 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 #678645]: Shrinkage of particles
Question #678645 on Yade changed: https://answers.launchpad.net/yade/+question/678645 Status: Needs information => Open Alireza Sadeghi gave more information on the question: Dear Jan, Thank you very much for your consideration. Actually , I want to reduce the particles' size in specific area for each iteration (time step), for example: r1= (1-5e-10) r0 where r0 is initial radius of particle in each iteratition and r1 is the final radius at the end of specific iteratition. > do you mean growParticle(), or really grow.Particle()? growParticle() doesn't work. It was a mistake in dictation :) > please again, always be as specific as possible. Did you get an error? If yes, which? ? The particles did not grow/shrink? ...? with growParticle() the particles disappeared. I don't know why. In addition, via the code for i in O.bodies: if isinstance(i.shape,Sphere) and i.state.pos[2]<-46e-3 and i.state.pos[2]>-48e-3: r2=i.shape.radius i.shape.radius=(1-2.25e-2)*r2 i.material=O.materials['Coke'] the radius will be reduced instantly. but I want their reduction in each step. Thank you very much for your help. Best Regards Alireza -- 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 #678646]: tracking a particle
Question #678646 on Yade changed: https://answers.launchpad.net/yade/+question/678646 Status: Answered => Solved Alireza Sadeghi confirmed that the question is solved: Hello Jan, Thank you very much for your response and help. Best Regards Alireza -- 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 #678646]: tracking a particle
New question #678646 on Yade: https://answers.launchpad.net/yade/+question/678646 Hello All, I have another question. How can I track an specific particle during the deformation? Thank you. Best Regards Alireza -- 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 #678645]: Shrinkage of particles
New question #678645 on Yade: https://answers.launchpad.net/yade/+question/678645 Hello All, I want to simulate shrinkage of particles during the deformation. the summary of my problem is: there are some particles in the cylinder and gravity are applied to them. The particles which are in the specific area will be reduced their radii. I did the below simulation but it is not the thing that I want to do. for i in O.bodies: if isinstance(i.shape,Sphere) and i.state.pos[2]<-46e-3 and i.state.pos[2]>-48e-3: r2=i.shape.radius i.shape.radius=(1-2.25e-2)*r2 i.material=O.materials['Coke'] I want to reduce the size of particles in a specific position. I used the grow.Particle() but it doesn't work. I was wondering if you could give me a solution. the reduction of particles' radius must be a function of time and the position of particles. Thank you very much for your help in advance. Best Regards Alireza -- 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 #677511]: access to yade-dem.org
Question #677511 on Yade changed: https://answers.launchpad.net/yade/+question/677511 Status: Needs information => Solved Alireza Sadeghi confirmed that the question is solved: After asking I saw them :) -- 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 #677511]: access to yade-dem.org
New question #677511 on Yade: https://answers.launchpad.net/yade/+question/677511 Hello all, I can't access to site yade-dem.org. Is there any problem with the website or the problem is from my PC? Thank you for your help. Best Regards Alireza -- 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 #676026]: definition of stress-strain in a sample
Question #676026 on Yade changed: https://answers.launchpad.net/yade/+question/676026 Status: Answered => Solved Alireza Sadeghi confirmed that the question is solved: Thanks Bruno Chareyre, 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 #676026]: definition of stress-strain in a sample
New question #676026 on Yade: https://answers.launchpad.net/yade/+question/676026 Hello All, I am new on YADE. I want to find stress and strain in a very simple sample. I don't know how can I use TesselationWrapper as described in [1]. could you please give me a simple example to find how can find stress or strain in the YADE. Thank you very much for your help in advance. Best Regards Alireza [1]https://www.yade-dem.org/doc/user.html#micro-stress-and-micro-strain -- 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 #673380]: VTKRecorder
New question #673380 on Yade: https://answers.launchpad.net/yade/+question/673380 Hello All, I have tow systems and I run my code on those. One of them, give me this error: elcome to Yade 2016-05-04.git-040f642 TCP python prompt on localhost:9000, auth cookie `akseys' XMLRPC info provider on http://localhost:21000 Running script T22.py WARN /Applications/Yade/trunk/pkg/dem/SpherePack.cpp:248 makeCloud: Exceeded 1000 tries to insert non-overlapping sphere to packing. Only 75 spheres were added, although you requested 80. Traceback (most recent call last): File "/usr/bin/yade", line 182, in runScript execfile(script,globals()) File "T22.py", line 96, in VTKRecorder(fileName='/Home/T22/3d-vtk-',recorders=['all'],iterPeriod=1000), NameError: name 'VTKRecorder' is not defined Do you know what should I do? Thank you very much for your help. Best Regards Alireza -- 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 #671132]: tow engines in a program
Question #671132 on Yade changed: https://answers.launchpad.net/yade/+question/671132 Status: Answered => Open Alireza Sadeghi is still having a problem: Hello Jan, Thank you for your response. I can't run it. it gives me the error: TypeError: No registered converter was able to produce a C++ rvalue of type int from this Python object of type Body when I run the below engines: #== engine for movement of particles = O.engines=[ ForceResetter(), InsertionSortCollider([Bo1_Box_Aabb(),Bo1_Sphere_Aabb()]), InteractionLoop( [Ig2_Sphere_Sphere_ScGeom(),Ig2_Facet_Sphere_ScGeom()], [Ip2_FrictMat_FrictMat_FrictPhys()], [Law2_ScGeom_FrictPhys_CundallStrack()] ), NewtonIntegrator(damping=0.4,gravity=[0,0,-9.8],label="newton"), ] #== engine for movement of piston === O.engines=[ ForceResetter(), InsertionSortCollider([Bo1_Box_Aabb(),Bo1_Sphere_Aabb()]), InteractionLoop( [Ig2_Sphere_Sphere_ScGeom(),Ig2_Box_Sphere_ScGeom()], [Ip2_JCFpmMat_JCFpmMat_JCFpmPhys()], [Law2_ScGeom_JCFpmPhys_JointedCohesiveFrictionalPM()] ), NewtonIntegrator(damping=0,gravity=[0,0,-9.8],label="newton",), TranslationEngine(translationAxis=[0,0,-5],velocity=0.5,ids=pst) GlobalStiffnessTimeStepper(active=1,timeStepUpdateInterval=10,timestepSafetyCoefficient=0.4, defaultDt=0.1*utils.PWaveTimeStep()), ] O.dt=.5*PWaveTimeStep() O.saveTmp() Bests Alireza -- 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 #671132]: tow engines in a program
Question #671132 on Yade changed: https://answers.launchpad.net/yade/+question/671132 Status: Needs information => Open Alireza Sadeghi gave more information on the question: Hello Jan, Thank you very much for your response. I want to simulate a simple model which have a mold and piston and particles in there. I applied a engine for movement of particles caused by gravity, after that I want to compress particles by the piston to find force chain in the system during compaction. Bests Alireza -- 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 #671132]: tow engines in a program
New question #671132 on Yade: https://answers.launchpad.net/yade/+question/671132 Hello All, I am a beginner in YADE. I want to simulate a simple compaction test which a piston will apply compaction force to particles. I want to use an engine for moving the particles and use another engine for movement of the piston. Could you please help me which engine should be the first?my engine for piston movement is: O.engines=[ ForceResetter(), InsertionSortCollider([Bo1_Box_Aabb(),Bo1_Sphere_Aabb()]), InteractionLoop( [Ig2_Sphere_Sphere_ScGeom(),Ig2_Box_Sphere_ScGeom()], [Ip2_JCFpmMat_JCFpmMat_JCFpmPhys()], [Law2_ScGeom_JCFpmPhys_JointedCohesiveFrictionalPM()] ), TranslationEngine(dead=False,translationAxis=[0,0,-5],velocity=0.5,ids=pst), # piston motion GlobalStiffnessTimeStepper(active=1,timeStepUpdateInterval=10,timestepSafetyCoefficient=0.4, defaultDt=0.1*utils.PWaveTimeStep()), NewtonIntegrator(damping=0,gravity=[0,0,-9.8]), ] Is it correct? Thank you very much for your consideration in advance. Best Regards Alireza -- 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 #665529]: about indentation error in example of yade website
Question #665529 on Yade changed: https://answers.launchpad.net/yade/+question/665529 Status: Answered => Solved Alireza Sadeghi confirmed that the question is solved: Hello Gentlemen, Thank you very much for your help. I check it again and I found my mistakes. Thanks a lot. Best Regards Alireza -- 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 #665529]: about indentation error in example of yade website
New question #665529 on Yade: https://answers.launchpad.net/yade/+question/665529 Hello, I have a question about indentation error in loading example number two (Gravity deposition) of https://yade-dem.org/doc/tutorial-examples.html When I run it, it say that unbalancedForce() is not defined! why does this error appear?how can I fix it? Does yade know how to define unbalancedForce() or I have to define for it? Thank you very much for your help in advance. Best Regards Alireza -- 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 #665440]: how can make a path for yade to load python files
Question #665440 on Yade changed: https://answers.launchpad.net/yade/+question/665440 Status: Answered => Solved Alireza Sadeghi confirmed that the question is solved: Dear Jan, Thank you very much for your help. It was about the folder that my program was in it. it was solved. Thanks a lot. Best Regards Alireza -- 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 #665440]: how can make a path for yade to load python files
Question #665440 on Yade changed: https://answers.launchpad.net/yade/+question/665440 Alireza Sadeghi 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 #665440]: how can make a path for yade to load python files
New question #665440 on Yade: https://answers.launchpad.net/yade/+question/665440 Hello, I am new user of YADE. I can't load any python files in the yade. when I write the file (for example: yade Hello.py) it shows me an error like this: Welcome to Yade 2016-05-04.git-040f642 TCP python prompt on localhost:9000, auth cookie `uaydks' XMLRPC info provider on http://localhost:21000 Running script Hello.py Traceback (most recent call last): File "/usr/bin/yade", line 182, in runScript execfile(script,globals()) IOError: [Errno 2] No such file or directory: 'Hello.py' [[ ^L clears screen, ^U kills line. F12 controller, F11 3d view (use h-key for showing help), F10 both, F9 generator, F8 plot. ]] Could you please help me?I don't know what should I do :( Thanks a lot -- 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
