Re: [Yade-users] [Question #688400]: Energy of the body is very high, defying gravity and normal stress - cpm material
Question #688400 on Yade changed: https://answers.launchpad.net/yade/+question/688400 Status: Needs information => Open Akm gave more information on the question: Hi, I downloaded the source code from gitlab - it mentions the yade as Yade unknown while running. I have checked everything that you said - but something is not right and I am unable to figure it out. -- 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 #688400]: Energy of the body is very high, defying gravity and normal stress - cpm material
Question #688400 on Yade changed: https://answers.launchpad.net/yade/+question/688400 Status: Needs information => Open Akm gave more information on the question: Hi Jan, I increased the normal load by 10 times, but still, the residual strength drops to zero. If I increase by 100x or something, the entire sample crushes under the load and there is no point in doing the simulation for that high normal stress. Please help me understand the cause and solution to this problem. Thanks, Arun -- 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 #688400]: Energy of the body is very high, defying gravity and normal stress - cpm material
Question #688400 on Yade changed: https://answers.launchpad.net/yade/+question/688400 Status: Solved => Open Akm is still having a problem: Sorry to bring this problem up again. I thought the 'flying up' phenomenon was due to very high tensile strength. But even when I use strengths of 0.8Mpa, I still can see that the sample rebounds and flys up, thus losing contact with the bottom box. I could verify it in two ways : 1. Seeing the images in ParaView - the upper box was literally up in the air without any contact with the lower box. 2. Mathematically, the residual shear strength should be "Normal stress*tan(phi)". But the residual went to zero. Even when I do the same code for just Direct shear where I have a plain interface, there is the same behaviour. The top sample is losing contact with the bottom sample at any cost. Can't we model the direct shear of two rigid bodies using Cpmmat? I tried using JCFPm as well but could see the same results. Can someone please guide me into this? Thanks, Arun -- 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 #688696]: Calculate the number of broken bonds in Cpm mat
Question #688696 on Yade changed: https://answers.launchpad.net/yade/+question/688696 Status: Answered => Solved Akm confirmed that the question is solved: Thanks Jan Stránský, that solved my question. -- You received this question notification because your team yade-users is an answer contact for Yade. ___ Mailing list: https://launchpad.net/~yade-users Post to : yade-users@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp
Re: [Yade-users] [Question #688696]: Calculate the number of broken bonds in Cpm mat
Question #688696 on Yade changed: https://answers.launchpad.net/yade/+question/688696 Status: Needs information => Open Akm gave more information on the question: Hi Jan and Robert, My script is pretty lengthy and has a lot of imports. So let me just stick to those commands where I calculate the bond breakage statistics. I do not get any change in the cohesive bond numbers even when I can see particles flying after bond breakage. ###At the start of the program global cohesive_count cohesive_count=0 for i in O.interactions: if hasattr(i.phys, 'isCohesive'): if i.phys.isCohesive == True: cohesive_count+=1 ###Inside a Pyrunner: presentcohesive_count = 0 for i in O.interactions: if hasattr(i.phys, 'isCohesive'): if i.phys.isCohesive == 1: presentcohesive_count+=1 noncohesive_count= cohesive_count-presentcohesive_count print (noncohesive_count) -- 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 #688696]: Calculate the number of broken bonds in Cpm mat
New question #688696 on Yade: https://answers.launchpad.net/yade/+question/688696 Hi, I am trying to calculate the number of bonds broken at each iteration in my simulation. I am using the cpm material and tried to plot it using the concept where I check if 'iscohesive==True ' but it doesn't work. Is there any other way to get these numbers? Many thanks, Arun -- 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 #688400]: Energy of the body is very high, defying gravity and normal stress - cpm material
Question #688400 on Yade changed: https://answers.launchpad.net/yade/+question/688400 Status: Needs information => Solved Akm confirmed that the question is solved: Sure many thanks. -- 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 #688400]: Energy of the body is very high, defying gravity and normal stress - cpm material
Question #688400 on Yade changed: https://answers.launchpad.net/yade/+question/688400 Status: Needs information => Open Akm gave more information on the question: I think that might have been the problem. I had given the tensile strength as too high and that has created a lot of problems for me. How do I know which version of Yade I am using? I had downloaded the package from https://github.com/yade/trunk. So it shows the version as Yade unknown. Should I be downloading it from somewhere else? Many thanks, Arun -- 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 #688400]: Energy of the body is very high, defying gravity and normal stress - cpm material
Question #688400 on Yade changed: https://answers.launchpad.net/yade/+question/688400 Status: Needs information => Open Akm gave more information on the question: Hi Jan, It's 1.20.0. Were you able to figure out any discrepancy in my code that led to the overall behaviour? The conditions are valid and experiments have been performed under those conditions. I was looking for an approximate behaviour like that. Please let me know if you found any issues out here. Many thanks, Arun -- 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 #688400]: Energy of the body is very high, defying gravity and normal stress - cpm material
Question #688400 on Yade changed: https://answers.launchpad.net/yade/+question/688400 Status: Needs information => Open Akm gave more information on the question: Hi Jan, Thanks for looking into the code. >what version of Yade you use? It is 1.2.0. For the sake of modifyng the source code, I had downloaded the Yade/trunk files - https://github.com/yade/trunk and then installed Yade using three folders as outlined in the installation manual: https://www.yade- dem.org/doc/installation.html?highlight=sudo%20add%20apt%20repository%20ppa%20yade%20users%20external So when I run it, I get it as Yade Unknown. >next time please make the question and code consistent Sorry, next time, will do. >I could not reproduce this (probably because there is cohesive connection between top and bottom) I have only changed three lines in the source code. It will be not good to paste the entire thing here. Let me put that source code in the drive as well so that you can download it and compile it. ''https://drive.google.com/drive/folders/1M2rKs- vnR3uHC_ybv2vliYLhvO2g0RTm?usp=sharing'' >realistic? Maybe a real specimen would also jump in reality moving horizontally 10 cm per second We know that the physical nature of the experiment cannot be simulated here in the simulations. This is a 1 particle sample. 0.1mps is a 125% strain rate which covers like 5% strain in 20 minutes. I tried to use a 12.5% strain rate as well, it covered 5% of the strain in 13 hours. But the sample behaved the same way like it did with higher strain rates. There is absolutely no difference. -- 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 #688400]: Energy of the body is very high, defying gravity and normal stress - cpm material
Question #688400 on Yade changed:
https://answers.launchpad.net/yade/+question/688400
Akm gave more information on the question:
I will post the full script here. I will attach the drive links for the
particles because I exported those particles using a gts surface. I did
the sample preparation in a separate file using radial expansion. So the
particle positions are here-kindly download those text files.
''https://drive.google.com/drive/folders/1M2rKs-
vnR3uHC_ybv2vliYLhvO2g0RTm?usp=sharing''
The complete script is as follows:
from yade import pack, qt, export, ymport, plot
#Default parameters or from table
readParamsFromTable(noTableOk=True, # unknownOk=True,
intfactor = 1.2, #Max value of 1.2 is preferred as it can be reverted back to
1 after a step.
Shear_Vel=0.1,
S_length=0.060,
S_area=0.060*0.030,
Sig0=300e3, #300kN/m2
kn=600e3, #300kPa/mm
plotperiod=500,
topname='/home/akm/Documents/Yade/install/bin/Full-10k-Top.txt',
botname='/home/akm/Documents/Yade/install/bin/Full-10k-Bot.txt',
spname='/home/akm/Documents/Yade/install/bin/snapshots-trial/plot',
vtkname='/home/akm/Documents/Yade/install/bin/vtk-trial/s',
)
from yade.params.table import *
F_init=Sig0*S_area #Initial Normal force
#MATERIAL SPECIFICATION
mat1=CpmMat(
young = 32e9,
density=2400,
poisson = 0.30,
frictionAngle = radians(35),
epsCrackOnset = 5e-1,
relDuctility = 1.1,
sigmaT = 1.2e6,
)
mat2=CpmMat(
young = 5e9,
density=1940,
poisson = 0.30,
frictionAngle = radians(35),
epsCrackOnset = 5e-2,
relDuctility = 1.0001,
sigmaT = 1.2e1,
)
O.materials.append(mat1)
O.materials.append(mat2)
#UPPER PARTICLE IMPORT:
upbox_particles=ymport.text(topname,material=mat2, color=(1,0,0))
O.bodies.append(upbox_particles)
#LOWER PARTICLE IMPORT:
lowbox_particles=ymport.text(botname,material=mat1, color=(0,1,0))
O.bodies.append(lowbox_particles)
###
##GEOMETRY CREATION :
O.materials.append(FrictMat(frictionAngle=0,density=0,label='walls'))
#Dimension of Packing
d1,d2=Vector3(.01,.01,.01),Vector3(.07,.0783,.04) # corners of the initial
packing
xmin=d1[0]
xmax=d2[0]
ymin=d1[1]
ymax=d2[1]
zmin=d1[2]
zmax=d2[2]
ht_lower=0.03
ht_upper=0.03
asp_ht=0.0083
#Shear plate for shearing in the +X direction
P_Shear=utils.box((xmin,ymin+(ht_lower/2)-(ht_lower/20),zmin+(zmax-zmin)/2 ),
(0,(ht_lower/2),(zmax-zmin)/2 ),fixed=True, wire=True, color=(0,0,1))
O.bodies.append(P_Shear)
P_Shear.state.blockedDOFs = "xyzXYZ"
#Surrounding setup geometry
U_box_1=geom.facetBox((xmin+(xmax-xmin)/2,
ymin+ht_lower+asp_ht+(asp_ht+ht_upper)/2,zmin+(zmax-zmin)/2 ),
((xmax-xmin)/2,(asp_ht+ht_upper)/2,(zmax-zmin)/2 ), wallMask=2) # +X Boundary
constraints
U_box_2=geom.facetBox((xmin+(xmax-xmin)/2,
ymin+ht_lower+(asp_ht*2+ht_upper)/2,zmin+(zmax-zmin)/2 ),
((xmax-xmin)/2,(asp_ht*2+ht_upper)/2,(zmax-zmin)/2 ), wallMask=1) # -X Boundary
constraints
L_box=geom.facetBox((xmin+(xmax-xmin)/2,
ymin+(asp_ht+ht_lower)/2,zmin+(zmax-zmin)/2 ),
((xmax-xmin)/2,(asp_ht+ht_lower)/2,(zmax-zmin)/2 ), wallMask=6)
Cover_box=geom.facetBox(((xmin+xmax)/2,(ymin+ymax)/2,(zmin+zmax)/2
),((xmin+xmax+0.06)/2,(ymin+ymax+0.02)/2,(zmin+zmax+0.02)/2 ) ,wallMask=63)
O.bodies.append(U_box_1)
O.bodies.append(U_box_2)
O.bodies.append(L_box)
O.bodies.append(Cover_box)
idlist=[]
idlist.append(P_Shear.id)
for i in L_box:
idlist.append(i.id) #ids for the translation motion of the box
for i in O.bodies:
if isinstance(i.shape,Sphere):
if i.mat==mat1:
i.state.blockedDOFs = "yzXYZ" #blocking the mat1 movement to ensure that
there is no moment at the end point
#Adding a plate at the top:
global plate
plate_pos=max([b.state.pos[1]+(b.shape.radius) for b in O.bodies if
isinstance(b.shape,Sphere)])
plate=utils.box((xmin+(xmax-xmin)/2,plate_pos,zmin+(zmax-zmin)/2),
((xmax-xmin)/2 - xmin/60, 0, (zmax-zmin)/2),fixed=True,
wire=False,color=(0,0,1))
O.bodies.append(plate)
print('Plate added')
plate.state.blockedDOFs = "xyzXYZ"
O.engines=[
ForceResetter(),
InsertionSortCollider([Bo1_Sphere_Aabb(aabbEnlargeFactor=intfactor,label='bo1s'),Bo1_Facet_Aabb(),Bo1_Wall_Aabb(),Bo1_Box_Aabb()]),
InteractionLoop(
[Ig2_Sphere_Sphere_ScGeom(interactionDetectionFactor=intfactor,label='ig2s'),Ig2_Facet_Sphere_ScGeom(),Ig2_Wall_Sphere_ScGeom(),Ig2_Box_Sphere_ScGeom()],
[Ip2_FrictMat_FrictMat_FrictPhys(),
Ip2_CpmMat_CpmMat_CpmPhys(cohesiveThresholdIter=1)],
[Law2_ScGeom_FrictPhys_CundallStrack(), Law2_ScGeom_CpmPhys_Cpm()]
),
NewtonIntegrator(gravity=(0,-9.81,0),damping=0.5),
GlobalStiffnessTimeStepper(active=True,timeStepUpdateInterval=1,timestepSafetyCoefficient=0.8),
PyRunner(iterPeriod=1,command='addPlotData()'),
]
O.step()
bo1s.aabbEnlargeFactor=1
ig2s.interactionDetectionFactor=1
def addPlotData():
Sf = -O.forces.f(P_Shear.id)[0] #Negative sign since the force applied is in
-X direction
Hor_dspl = 1000*P_Shear.state.displ()[0] #Hor disp in mm
Shear_stress
Re: [Yade-users] [Question #688400]: Energy of the body is very high, defying gravity and normal stress - cpm material
Question #688400 on Yade changed:
https://answers.launchpad.net/yade/+question/688400
Akm gave more information on the question:
I will post the full script here. I will attach the drive links for the
particles because I exported those particles using a gts surface. I did
the sample preparation in a separate file using radial expansion. So the
particle positions are here-kindly download those text files.
''https://drive.google.com/drive/folders/1M2rKs-
vnR3uHC_ybv2vliYLhvO2g0RTm?usp=sharing''
The complete script is as follows:
from yade import pack, qt, export, ymport, plot
#Default parameters or from table
readParamsFromTable(noTableOk=True, # unknownOk=True,
intfactor = 1.2, #Max value of 1.2 is preferred as it can be reverted
back to 1 after a step.
Shear_Vel=0.1,
S_length=0.060,
S_area=0.060*0.030,
Sig0=300e3, #300kN/m2
kn=600e3, #300kPa/mm
F_init=Sig0*S_area,
#Initial Normal force
plotperiod=500,
topname='/home/akm/Documents/Yade/install/bin/Full-10k-Top.txt',
botname='/home/akm/Documents/Yade/install/bin/Full-10k-Bot.txt',
spname='/home/akm/Documents/Yade/install/bin/snapshots-trial/plot',
vtkname='/home/akm/Documents/Yade/install/bin/vtk-trial/s',
)
from yade.params.table import *
#MATERIAL SPECIFICATION
mat1=CpmMat(
young = 32e9,
density=2400,
poisson = 0.30,
frictionAngle = radians(35),
epsCrackOnset = 5e-1,
relDuctility = 1.1,
sigmaT = 1.2e6,
)
mat2=CpmMat(
young = 5e9,
density=1940,
poisson = 0.30,
frictionAngle = radians(35),
epsCrackOnset = 5e-2,
relDuctility = 1.0001,
sigmaT = 1.2e1,
)
O.materials.append(mat1)
O.materials.append(mat2)
#UPPER PARTICLE IMPORT:
upbox_particles=ymport.text(topname,material=mat2, color=(1,0,0))
O.bodies.append(upbox_particles)
#LOWER PARTICLE IMPORT:
lowbox_particles=ymport.text(botname,material=mat1, color=(0,1,0))
O.bodies.append(lowbox_particles)
###
##GEOMETRY CREATION :
O.materials.append(FrictMat(frictionAngle=0,density=0,label='walls'))
#Dimension of Packing
d1,d2=Vector3(.01,.01,.01),Vector3(.07,.0783,.04) # corners of the initial
packing
xmin=d1[0]
xmax=d2[0]
ymin=d1[1]
ymax=d2[1]
zmin=d1[2]
zmax=d2[2]
ht_lower=0.03
ht_upper=0.03
asp_ht=0.0083
#Shear plate for shearing in the +X direction
P_Shear=utils.box((xmin,ymin+(ht_lower/2)-(ht_lower/20),zmin+(zmax-zmin)/2 ),
(0,(ht_lower/2),(zmax-zmin)/2 ),fixed=True, wire=True, color=(0,0,1))
O.bodies.append(P_Shear)
P_Shear.state.blockedDOFs = "xyzXYZ"
#Surrounding setup geometry
U_box_1=geom.facetBox((xmin+(xmax-xmin)/2,
ymin+ht_lower+asp_ht+(asp_ht+ht_upper)/2,zmin+(zmax-zmin)/2 ),
((xmax-xmin)/2,(asp_ht+ht_upper)/2,(zmax-zmin)/2 ), wallMask=2) # +X
Boundary constraints
U_box_2=geom.facetBox((xmin+(xmax-xmin)/2,
ymin+ht_lower+(asp_ht*2+ht_upper)/2,zmin+(zmax-zmin)/2 ),
((xmax-xmin)/2,(asp_ht*2+ht_upper)/2,(zmax-zmin)/2 ), wallMask=1) # -X
Boundary constraints
L_box=geom.facetBox((xmin+(xmax-xmin)/2,
ymin+(asp_ht+ht_lower)/2,zmin+(zmax-zmin)/2 ),
((xmax-xmin)/2,(asp_ht+ht_lower)/2,(zmax-zmin)/2 ), wallMask=6)
Cover_box=geom.facetBox(((xmin+xmax)/2,(ymin+ymax)/2,(zmin+zmax)/2
),((xmin+xmax+0.06)/2,(ymin+ymax+0.02)/2,(zmin+zmax+0.02)/2 ) ,wallMask=63)
O.bodies.append(U_box_1)
O.bodies.append(U_box_2)
O.bodies.append(L_box)
O.bodies.append(Cover_box)
idlist=[]
idlist.append(P_Shear.id)
for i in L_box:
idlist.append(i.id) #ids for the translation motion
of the box
for i in O.bodies:
if isinstance(i.shape,Sphere):
if i.mat==mat1:
i.state.blockedDOFs = "yzXYZ" #blocking the mat1
movement to ensure that there is no moment at the end point
#Adding a plate at the top:
global plate
plate_pos=max([b.state.pos[1]+(b.shape.radius) for b in O.bodies if
isinstance(b.shape,Sphere)])
plate=utils.box((xmin+(xmax-xmin)/2,plate_pos,zmin+(zmax-zmin)/2),
((xmax-xmin)/2 - xmin/60, 0, (zmax-zmin)/2),fixed=True,
wire=False,color=(0,0,1))
O.bodies.append(plate)
print('Plate added')
plate.state.blockedDOFs = "xyzXYZ"
O.engines=[
ForceResetter(),
InsertionSortCollider([Bo1_Sphere_Aabb(aabbEnlargeFactor=intfactor,label='bo1s'),Bo1_Facet_Aabb(),Bo1_Wall_Aabb(),Bo1_Box_Aabb()]),
InteractionLoop(
[Ig2_Sphere_Sphere_ScGeom(interactionDetecti
Re: [Yade-users] [Question #688399]: How to simulate a single particle crushing test?
Question #688399 on Yade changed: https://answers.launchpad.net/yade/+question/688399 Status: Open => Answered Akm proposed the following answer: Hi Vincent, 1. There are two options to provide a constant velocity to the plate. You can use the translation engine or use the body.state.vel=(0,0,0,) command for whatever direction you want it in. 2. You can check the force on the wall by the command "O.forces.f(body.id)" Hope it helps. -- 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 #688400]: Energy of the body is very high, defying gravity and normal stress - cpm material
New question #688400 on Yade: https://answers.launchpad.net/yade/+question/688400 Hi all, I am trying to simulate a shear test where I have two CPM materials with different properties - top_mat(lower strength) and bot_mat(higher strength). top_mat rests on bot_mat at an angle <45. Try to imagine a direct shear test where there is an asperity at the interface. There is only friction between the two materials-no cohesion. There is a plate at the top which applies a load on top_mat in Constant Normal stiffness condition(CNS). top_mat has its boundaries free only in the Y direction where it is allowed to move vertically up or down. bot_mat can only move in the X direction(horizontal). When bot_mat(bottom box) is moved using a plate in the horizontal direction(X), the top_mat should ''SLIDE'' along the asperity in the Y direction and then be subjected to damage when the shear strength threshold is reached for top_mat. But even though I have a plate with a CNS load applied on top, during the shearing process, top_mat starts to ''FLY'' up instead of sliding and loses contact with bot_mat by pushing the CNS plate up, thereby defying both the CNS and the gravitational forces in the system. I am not sure why this is happening. P.S: I had modified the source code to have zero cohesion ONLY AT THE INTERFACE: iscohesive=0, undamagedcohesion=0 Sorry that I am unable to think of a way to provide an MWE because I cut out that interface using a gts surface and this platform doesn't support the inclusion of external links. If it is alright, I can provide the code and the external links here. Kindly help me out as I have been stuck in this issue for many days. Help as soon as possible will be much appreciated. Many thanks in advance. -Arun -- 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 #688321]: Not able to apply DOF constraints
Question #688321 on Yade changed: https://answers.launchpad.net/yade/+question/688321 Status: Needs information => Open Akm gave more information on the question: The box started to disappear after a while when I entered the specific degrees of freedom to be blocked. After your advice to include all the DOFs as blocked, this works perfectly fine. Many thanks. -- 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 #688321]: Not able to apply DOF constraints
New question #688321 on Yade: https://answers.launchpad.net/yade/+question/688321 Hi all, I tried to move a sphere(Cpmmat) using a utilsbox acting as a wall(Frictmat). I want to restrain the degrees of freedom in the other two directions but when I do that, the wall tends to fly away after bouncing on the sphere. Here is a MWE of that phenomenon. Kindly help me to understand where I have gone wrong. from yade import pack, qt, export, ymport, plot O.reset() concrete_material=CpmMat( young = 1e10, density=1400, poisson = 0.25, frictionAngle = radians(40), epsCrackOnset = 1e-4, relDuctility = 10, sigmaT = 1e18, label='concrete_mat' ) mat1=O.materials.append(concrete_material) #O.materials.append(FrictMat(frictionAngle=0,density=0,label='walls')) O.bodies.append([ sphere(center=(0,0,0),radius=.5,fixed=False,material=concrete_material), ]) O.materials.append(FrictMat(frictionAngle=0,density=0,label='walls')) O.engines=[ ForceResetter(), InsertionSortCollider([Bo1_Sphere_Aabb(aabbEnlargeFactor=1,label='bo1s'),Bo1_Facet_Aabb(),Bo1_Wall_Aabb(),Bo1_Box_Aabb()]), InteractionLoop( [Ig2_Sphere_Sphere_ScGeom(interactionDetectionFactor=1,label='ig2s'),Ig2_Facet_Sphere_ScGeom(),Ig2_Wall_Sphere_ScGeom(),Ig2_Box_Sphere_ScGeom()], [Ip2_FrictMat_FrictMat_FrictPhys(), Ip2_CpmMat_CpmMat_CpmPhys(cohesiveThresholdIter=1)], [Law2_ScGeom_FrictPhys_CundallStrack(), Law2_ScGeom_CpmPhys_Cpm()] ), NewtonIntegrator(gravity=(0,0,0),damping=0.5), #CpmStateUpdater(iterPeriod=1), GlobalStiffnessTimeStepper(active=True,timeStepUpdateInterval=1,timestepSafetyCoefficient=0.8), ] O.step() U_box= utils.box((.5,1,-.5),(1,0,.5),fixed=True,wire=False,color=(1,0,0)) O.bodies.append(U_box) U_box.state.mass=1 U_box.state.inertia = (1,1,1) U_box.state.blockedDOFs = "xzXYZ" U_box.state.vel=(0,-0.01,0) yade.qt.View() O.saveTmp() ## Many thanks in advance - Arun -- 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 #687294]: Interaction issue with permanent force function
Question #687294 on Yade changed: https://answers.launchpad.net/yade/+question/687294 Status: Answered => Solved Akm confirmed that the question is solved: Thanks Jan Stránský, that solved my question. -- You received this question notification because your team yade-users is an answer contact for Yade. ___ Mailing list: https://launchpad.net/~yade-users Post to : yade-users@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp
Re: [Yade-users] [Question #687294]: Interaction issue with permanent force function
Question #687294 on Yade changed: https://answers.launchpad.net/yade/+question/687294 Akm posted a new comment: Hi Jerome, Thanks for the call. I checked and there are no interactions between the two. How do I solve this issue? -- 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 #687294]: Interaction issue with permanent force function
New question #687294 on Yade: https://answers.launchpad.net/yade/+question/687294 I am trying to apply a permanent force on a sphere with the help of a utils box(only one side as a plate). But when I apply a permanent force on the plate, the plate passes through the sphere as if there is no interaction instead of moving the sphere. Kindly help me with this issue. Yade 20191217-2900~6d5e1cb~xenial1 Using python version: 3.5.2 Here is the Minimal working script: from yade import pack, qt, export, ymport, plot O.reset() concrete_material=CpmMat( young = 1e10, density=1400, poisson = 0.25, frictionAngle = radians(40), epsCrackOnset = 1e-4, relDuctility = 10, sigmaT = 1e18, label='concrete_mat' ) mat1=O.materials.append(concrete_material) #O.materials.append(FrictMat(frictionAngle=0,density=0,label='walls')) O.bodies.append([ sphere(center=(0,0,0),radius=.5,fixed=False,material='concrete_mat'), ]) O.engines=[ ForceResetter(), InsertionSortCollider([Bo1_Sphere_Aabb(aabbEnlargeFactor=1,label='bo1s'),Bo1_Facet_Aabb(),Bo1_Wall_Aabb(),Bo1_Box_Aabb()]), InteractionLoop( [Ig2_Sphere_Sphere_ScGeom(interactionDetectionFactor=1,label='ig2s'),Ig2_Facet_Sphere_ScGeom(),Ig2_Wall_Sphere_ScGeom(),Ig2_Box_Sphere_ScGeom()], [Ip2_CpmMat_CpmMat_CpmPhys(cohesiveThresholdIter=1)], [Law2_ScGeom_CpmPhys_Cpm()] ), NewtonIntegrator(gravity=(0,0,0),damping=0.5), GlobalStiffnessTimeStepper(active=True,timeStepUpdateInterval=1,timestepSafetyCoefficient=0.8), ] O.step() U_box= utils.box((.5,1,-.5),(1,0,.5),fixed=False,wire=False,color=(1,0,0)) O.bodies.append(U_box) U_box.state.mass=10 O.forces.setPermF(U_box.id,(0,-20,0)) yade.qt.View() O.saveTmp() Thanks in advance Arun -- 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 #683934]: Bond broken and non movement of cohesive particles
Question #683934 on Yade changed: https://answers.launchpad.net/yade/+question/683934 Status: Needs information => Open Akm gave more information on the question: Hi Jan, Thanks for the call. But if I do not use the gts surface, then the essence of the question will be lost. I will try to upload an MWE though. -- 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 #683934]: Bond broken and non movement of cohesive particles
Question #683934 on Yade changed:
https://answers.launchpad.net/yade/+question/683934
Description changed to:
Hi all,
Hope you all are doing well. I have a question with respect to the
CohFric model. I am trying to implement a shearing mechanism between two
body masses which are of different material properties. It is kind of
like a direct shear test. The setup consists of two boxes - top and
bottom. The material is contained in them. I created the particles
through radial expansion, used a gts surface predicate to cut them
appropriately and then imported them in another file with different
material properties.
My yade version is Yade 2018.02b-142067caa7~xenial
Issue 1:
I tried to do the shearing through a translation engine after applying a
cohesive bond between the particles of the same material. But when I
apply a force on the upper box , the particles inside should move as
well with the box. But only the box moves and the cohesive network of
particles just stay there without any movement. Can anyone please tell
me where I went wrong with the code?
Issue 2:
The particles in the lower box tend to break their cohesive bond after a
particular iteration. I am unable to comprehend the reason why this is
happening.
I have attached the code below for your reference. Since I have created
the particles separately using a gts predicate, I have attached the
drive link for the particle import here as well.
https://drive.google.com/drive/folders/1bCS_w917aP0Tr_Jfy6YfN20c9W0cdJYM?usp=sharing
Note: Kindly download the particle file using the download link above to
impor the particles.
from yade import pack, qt, export, ymport, plot
#MATERIAL SPECIFICATION:
rock_material=CohFrictMat(
young=6e10 ,poisson=0.21 ,density=3400 ,frictionAngle=radians(45)
,normalCohesion=1e20,shearCohesion=1e8 ,isCohesive=True, label='rock_mat')
mat1=O.materials.append(rock_material)
concrete_material=CohFrictMat(
young=1e10 ,poisson=0.25 ,density=1400 ,frictionAngle=radians(30)
,normalCohesion=1e20,shearCohesion=1e6 ,isCohesive=True, label='concrete_mat')
mat2=O.materials.append(concrete_material)
O.materials.append(FrictMat(frictionAngle=0,density=0,label='walls'))
#UPPER PARTICLE IMPORT:
upbox_particles=ymport.text('Upperbox_particles-22.5.txt',material='rock_mat',
color=(1,0,0)) #Arguement passed to utils.sphere
O.bodies.append(upbox_particles)
shiftlist=[]
for i in O.bodies:
if (i.shape, Sphere):
shiftlist.append(i.id)
shiftBodies(shiftlist, (0,0.005,0))
#LOWER PARTICLE IMPORT:
lowbox_particles=ymport.text('Lowerbox_particles-22.5.txt',material='concrete_mat',
color=(0,1,0)) #Arguement passed to utils.sphere
O.bodies.append(lowbox_particles)
#GEOMETRY DESIGN:
#Creating the upper box
U_box=geom.facetBox((.13,.11,0.09),(.08,.02,.04),wallMask=51)
O.bodies.append(U_box)
#Creating the lower box
L_box=geom.facetBox((.13,.07,0.09),(.08,.02,.04),wallMask=55)
O.bodies.append(L_box)
#Creating the extension for sliding plane so particles don't fall out
#Give extra 1e-3 m in Y dir. so that the two facet walls...
#...at the shear interface do not rub against each other
U_box_extension1=geom.facetBox((.01,.111,0.09),(.04,.02,.04),wallMask=4)
O.bodies.append(U_box_extension1)
U_box_extension2=geom.facetBox((.25,.111,0.09),(.04,.02,.04),wallMask=4)
O.bodies.append(U_box_extension2)
L_box_extension1=geom.facetBox((.01,.11,0.09),(.04,.02,.04),wallMask=4)
O.bodies.append(L_box_extension1)
L_box_extension2=geom.facetBox((.25,.11,0.09),(.04,.02,.04),wallMask=4)
O.bodies.append(L_box_extension2)
intfactor = 1.2
T_Vel=10
#SIMULATION PROCESSES:
O.engines=[
ForceResetter(),
InsertionSortCollider([Bo1_Sphere_Aabb(aabbEnlargeFactor=intfactor,label='bo1s'),Bo1_Facet_Aabb(),Bo1_Wall_Aabb(),Bo1_Box_Aabb()]),
InteractionLoop(
[Ig2_Sphere_Sphere_ScGeom6D(interactionDetectionFactor=intfactor,label='ig2s'),Ig2_Facet_Sphere_ScGeom(),Ig2_Wall_Sphere_ScGeom(),Ig2_Box_Sphere_ScGeom6D()],
[Ip2_FrictMat_FrictMat_FrictPhys(),
Ip2_CohFrictMat_CohFrictMat_CohFrictPhys(setCohesionNow=True,label="cohesiveIp",
normalCohesion=MatchMaker(matches=((mat1,mat1,1e20),(mat2,mat2,1e20),(mat1,mat2,0))),
shearCohesion=MatchMaker(matches=((mat1,mat1,1e9),(mat2,mat2,1e9),(mat1,mat2,0],
[Law2_ScGeom_FrictPhys_CundallStrack(),
Law2_ScGeom6D_CohFrictPhys_CohesionMoment()]
),
NewtonIntegrator(gravity=(0,0,0),damping=0.3),
GlobalStiffnessTimeStepper(active=True,timeStepUpdateInterval=1,timestepSafetyCoefficient=0.8),
#using global timestep for faster simulation
#PyRunner(command='check()',iterPeriod=10,label='runner')
]
O.step()
#for i in O.interactions:
#i.phys.normalAdhesion = 1e9
#i.phys.shearAdhesion = 1e3
bo1s.aabbEnlargeFactor=1
ig2s.interactionDetectionFactor=1
#Setting default material values:
rock_material.normalCohesion=1e9
concrete_material.normalCohesion=1e9
[Yade-users] [Question #683934]: Bond broken and non movement of cohesive particles
New question #683934 on Yade:
https://answers.launchpad.net/yade/+question/683934
Hi all,
Hope you all are doing well. I have a question with respect to the CohFric
model. I am trying to implement a shearing mechanism between two body masses
which are of different material properties. It is kind of like a direct shear
test. The setup consists of two boxes - top and bottom. The material is
contained in them. I created the particles through radial expansion, used a
gts surface predicate to cut them appropriately and then imported them in
another file with different material properties.
My yade version is Yade 2018.02b-142067caa7~xenial
Issue 1:
I tried to do the shearing through a translation engine after applying a
cohesive bond between the particles of the same material. But when I apply a
force on the upper box , the particles inside should move as well with the box.
But only the box moves and the cohesive network of particles just stay there
without any movement. Can anyone please tell me where I went wrong with the
code?
Issue 2:
The particles in the lower box tend to break their cohesive bond after a
particular iteration. I am unable to comprehend the reason why this is
happening.
I have attached the code below for your reference. Since I have created the
particles separately using a gts predicate, I have attached the drive link for
the particle import here as well.
https://drive.google.com/drive/folders/1bCS_w917aP0Tr_Jfy6YfN20c9W0cdJYM?usp=sharing
Note: Kindly download the particle file using the download link above to impor
the particles.
from yade import pack, qt, export, ymport, plot
#MATERIAL SPECIFICATION:
rock_material=CohFrictMat(
young=6e10 ,poisson=0.21 ,density=3400 ,frictionAngle=radians(45)
,normalCohesion=1e20,shearCohesion=1e8 ,isCohesive=True, label='rock_mat')
mat1=O.materials.append(rock_material)
concrete_material=CohFrictMat(
young=1e10 ,poisson=0.25 ,density=1400 ,frictionAngle=radians(30)
,normalCohesion=1e20,shearCohesion=1e6 ,isCohesive=True, label='concrete_mat')
mat2=O.materials.append(concrete_material)
O.materials.append(FrictMat(frictionAngle=0,density=0,label='walls'))
#UPPER PARTICLE IMPORT:
upbox_particles=ymport.text('Upperbox_particles-22.5.txt',material='rock_mat',
color=(1,0,0)) #Arguement passed to utils.sphere
O.bodies.append(upbox_particles)
shiftlist=[]
for i in O.bodies:
if (i.shape, Sphere):
shiftlist.append(i.id)
shiftBodies(shiftlist, (0,0.005,0))
#LOWER PARTICLE IMPORT:
lowbox_particles=ymport.text('Lowerbox_particles-22.5.txt',material='concrete_mat',
color=(0,1,0)) #Arguement passed to utils.sphere
O.bodies.append(lowbox_particles)
#GEOMETRY DESIGN:
#Creating the upper box
U_box=geom.facetBox((.13,.11,0.09),(.08,.02,.04),wallMask=51)
O.bodies.append(U_box)
#Creating the lower box
L_box=geom.facetBox((.13,.07,0.09),(.08,.02,.04),wallMask=55)
O.bodies.append(L_box)
#Creating the extension for sliding plane so particles don't fall out
#Give extra 1e-3 m in Y dir. so that the two facet walls...
#...at the shear interface do not rub against each other
U_box_extension1=geom.facetBox((.01,.111,0.09),(.04,.02,.04),wallMask=4)
O.bodies.append(U_box_extension1)
U_box_extension2=geom.facetBox((.25,.111,0.09),(.04,.02,.04),wallMask=4)
O.bodies.append(U_box_extension2)
L_box_extension1=geom.facetBox((.01,.11,0.09),(.04,.02,.04),wallMask=4)
O.bodies.append(L_box_extension1)
L_box_extension2=geom.facetBox((.25,.11,0.09),(.04,.02,.04),wallMask=4)
O.bodies.append(L_box_extension2)
intfactor = 1.2
T_Vel=10
#SIMULATION PROCESSES:
O.engines=[
ForceResetter(),
InsertionSortCollider([Bo1_Sphere_Aabb(aabbEnlargeFactor=intfactor,label='bo1s'),Bo1_Facet_Aabb(),Bo1_Wall_Aabb(),Bo1_Box_Aabb()]),
InteractionLoop(
[Ig2_Sphere_Sphere_ScGeom6D(interactionDetectionFactor=intfactor,label='ig2s'),Ig2_Facet_Sphere_ScGeom(),Ig2_Wall_Sphere_ScGeom(),Ig2_Box_Sphere_ScGeom6D()],
[Ip2_FrictMat_FrictMat_FrictPhys(),
Ip2_CohFrictMat_CohFrictMat_CohFrictPhys(setCohesionNow=True,label="cohesiveIp",
normalCohesion=MatchMaker(matches=((mat1,mat1,1e20),(mat2,mat2,1e20),(mat1,mat2,0))),
shearCohesion=MatchMaker(matches=((mat1,mat1,1e9),(mat2,mat2,1e9),(mat1,mat2,0],
[Law2_ScGeom_FrictPhys_CundallStrack(),
Law2_ScGeom6D_CohFrictPhys_CohesionMoment()]
),
NewtonIntegrator(gravity=(0,0,0),damping=0.3),
GlobalStiffnessTimeStepper(active=True,timeStepUpdateInterval=1,timestepSafetyCoefficient=0.8),
#using global timestep for faster simulation
#PyRunner(command='check()',iterPeriod=10,label='runner')
]
O.step()
#for i in O.interactions:
#i.phys.normalAdhesion = 1e9
#i.phys.shearAdhesion = 1e3
bo1s.aabbEnlargeFactor=1
ig2s.interactionDetectionFactor=1
#Setting default material values:
rock_material.normalCohesion=1e9
concrete_material.normalCohesion=1e9
Re: [Yade-users] [Question #683526]: Interaction problem between two bodies
Question #683526 on Yade changed: https://answers.launchpad.net/yade/+question/683526 Akm posted a new comment: Thanks, Jan. I guess that the problem is with the version of Yade. I will upgrade it appropriately. -- 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 #683526]: Interaction problem between two bodies
Question #683526 on Yade changed: https://answers.launchpad.net/yade/+question/683526 Status: Open => Invalid Akm rejected the question: The question contains sample data -- 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 #683526]: Interaction problem between two bodies - Direct Shear
Question #683526 on Yade changed: https://answers.launchpad.net/yade/+question/683526 Akm posted a new comment: Can anyone please tell me where I have gone wrong with the code? The 'dynamic' property is deprecated as listed in the Yade documentation. So I am using the 'Fixed' parameter to specify the body's movement. I have tried several ways but when I apply a permanent force on the main plate, the plate just passes through the spheres without any interaction. I am unsure about why this is happening. Help would be very much appreciated. -- 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 #683526]: Interaction problem between two bodies - Direct Shear
Question #683526 on Yade changed:
https://answers.launchpad.net/yade/+question/683526
Status: Needs information => Open
Akm gave more information on the question:
from yade import pack, plot, qt
O.reset()
#Creating the upper box
U_box=geom.facetBox((.2,.1,.8),(.2,.1,.2),wallMask=51)
O.bodies.append(U_box)
#Creating the lower box
L_box=geom.facetBox((.2,-.1,.8),(.2,.1,.2),wallMask=55)
O.bodies.append(L_box)
'''
#Creating the upper box
U_left=utils.box((0,.1,.8),(0,.1,.2),dynamic=False,fixed=True,wire=True,color=None)
O.bodies.append(U_left)
U_right=utils.box((.4,.1,.8),(0,.1,.2),dynamic=False,fixed=True,wire=True,color=None)
O.bodies.append(U_right)
U_front=utils.box((.2,.1,1),(.2,.1,0),dynamic=False,fixed=True,wire=True,color=None)
O.bodies.append(U_front)
U_back=utils.box((.2,.1,.6),(.2,.1,0),dynamic=False,fixed=True,wire=True,color=None)
O.bodies.append(U_back)
#Creating the lower box
L_left=utils.box((0,-.1,.8),(0,.1,.2),dynamic=False,fixed=True,wire=True,color=None)
O.bodies.append(L_left)
L_right=utils.box((.4,-.1,.8),(0,.1,.2),dynamic=False,fixed=True,wire=True,color=None)
O.bodies.append(L_right)
L_front=utils.box((.2,-.1,1),(.2,.1,0),dynamic=False,fixed=True,wire=True,color=None)
O.bodies.append(L_front)
L_back=utils.box((.2,-.1,.6),(.2,.1,0),dynamic=False,fixed=True,wire=True,color=None)
O.bodies.append(L_back)
L_bottom
=utils.box((.2,-0.05,.8),(.2,0,.2),dynamic=False,fixed=True,wire=True,color=None)
O.bodies.append(L_bottom)
#Creating the extension for sliding plane so particles don't fall out
#Give extra .01m in Y dir. so that the two facet walls...
#...at the shear interface do not rub against each other
U_box_extension1=geom.facetBox((-.2,.5,.8),(.2,.5,.2),wallMask=4)
O.bodies.append(U_box_extension1)
U_box_extension2=geom.facetBox((.6,.5,.8),(.2,.5,.2),wallMask=4)
O.bodies.append(U_box_extension2)
#Creating the extension for sliding plane so particles don't fall out
L_box_extension1=geom.facetBox((.6,-.5001,.8),(.2,.5,.2),wallMask=8)
O.bodies.append(L_box_extension1)
L_box_extension2=geom.facetBox((-.2,-.5001,.8),(.2,.5,.2),wallMask=8)
O.bodies.append(L_box_extension2)
#Material assignment- normalcohesion = Tensile strength, shearCohesion = Shear
strength
sample_material=FrictMat(
young=50e6 ,poisson=0.25 ,density=1800
,frictionAngle=radians(45),label='sample_mat')
O.materials.append(sample_material)
wall_material= FrictMat(
young=4e9 ,poisson=0.25 ,density=5000 ,frictionAngle=radians(50)
,label='wall_mat')
O.materials.append(wall_material)
#Speed of the translation engine
Vel=1
#Sphere Pack Creation
sp=pack.SpherePack()
sp.makeCloud((0,-0.05,0.6),(0.4,0.2,1.0),rMean=0.008,rRelFuzz=0.3,porosity=0.5,
num = 6000)
sp.toSimulation(material=sample_material) #already one cloud of particles here
#Collecting the ids for translation engine run
id_list=[]
#for i in L_box: id_list.append(i.id)#for facet box
id_list.append(L_left.id)
id_list.append(L_right.id)
id_list.append(L_front.id)
id_list.append(L_back.id)
id_list.append(L_bottom.id)
for i in L_box_extension1: id_list.append(i.id)
for i in L_box_extension2: id_list.append(i.id)
#Simulation
O.engines=[
ForceResetter(),
InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb(),Bo1_Wall_Aabb(),Bo1_Box_Aabb()]),
InteractionLoop(
[Ig2_Sphere_Sphere_ScGeom6D(),Ig2_Facet_Sphere_ScGeom(),Ig2_Wall_Sphere_ScGeom(),Ig2_Box_Sphere_ScGeom6D()],
[Ip2_FrictMat_FrictMat_FrictPhys(),
Ip2_CohFrictMat_CohFrictMat_CohFrictPhys(setCohesionNow=False,label="cohesiveIp")
],
[Law2_ScGeom_FrictPhys_CundallStrack(),
Law2_ScGeom6D_CohFrictPhys_CohesionMoment()]
),
NewtonIntegrator(gravity=(0,-1,0),damping=0.3),
GlobalStiffnessTimeStepper(active=True,timeStepUpdateInterval=1,timestepSafetyCoefficient=0.8),
#using global timestep for faster simulation
# (command='checkUnbalanced()',realPeriod=0.5,label='checker'),
PyRunner(command='addplate()',iterPeriod=10,label='runner') #real period is
actually not good
]
maxLoad=1e5
minLoad=1e1
inst_position_lower=min([b.state.pos[1]-(b.shape.radius) for b in
O.bodies if isinstance(b.shape,Sphere)])
global plate_b
plate_b=wall(min([b.state.pos[1]-(b.shape.radius) for b in O.bodies if
isinstance(b.shape,Sphere)]),axis=1,sense=+1)
#O.bodies.append(plate_b)
print('Plate in -Y added')
#Plate is being added to the simulation to consolidate the sample
def addplate():
#if unbalancedForce()<0.5 and O.iter>15000:
global plate
inst_position1=max([b.state.pos[1]+(b.shape.radius) for b in O.bodies
if isinstance(b.shape,Sphere)])
plate=utils.box((.2,inst_position1,.8),(.199,0,.199),dynamic=False,fixed=False,wire=False,color=(0,0,1))#,material=wall_material)
O.bodies.append(plate)
print('Plate added')
plate.state.vel=(0,-2,0)
runner.command='unloadPlate()'
#Plate rebounding after consolidation upto a maximum load
[Yade-users] [Question #683526]: Interaction problem between two bodies - Direct Shear
New question #683526 on Yade: https://answers.launchpad.net/yade/+question/683526 Hi, I have been trying to simulate the direct shear test in YADE version 1.20.0. I am able to construct the code upto the point where the box moves to create shearing in the sample. But there is a problem in the simulation at the last step. When I try to apply a permanent force on the plate to induce a constant normal load condition on the sample, the plate passes through the sample as if there is no interaction between the bodies at all. The expected behaviour is that the plate applying load on the sample permanently. But it doesn't seem to happen. Kindly help me with this issue. I have attached the link to the code below. https://drive.google.com/file/d/1jiS_ndotegcLJJsVmaaTnFvdhcZfkMeV/view?usp=sharing Thanks in advance. -- 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 #683390]: Constant Normal Load for Direct Shear Test
Question #683390 on Yade changed: https://answers.launchpad.net/yade/+question/683390 Akm posted a new comment: Hi Jerome, Thanks for the tip. I did not mean to use triaxial to simulate the direct shear experiment, I only wanted to impose a constant normal load on the sample using the engine and then set the upper box with a prescribed constant horizontal velocity. -- 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 #683390]: Constant Normal Load for Direct Shear Test
New question #683390 on Yade: https://answers.launchpad.net/yade/+question/683390 Hello, I am a beginner in Yade trying to create the simulation for a direct shear test. I am using a faceted box to create the sample. But I am unclear on how to provide a constant load at the top plate to simulate the condition. I tried to use the triaxial stress controller but couldn't do it as the top box of the sample has to move horizontally and Triaxial stress controller doesn't work that way. Kindly let me know your suggestions on how to achieve this. I can upload the code in launchpad if you need to take a look. If someone could upload an example code for direct shear, I would be more than happy. -- 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 #683113]: Clarification required for Ques 402952: Using a gts surface as a shell and filling it using radius expansion method
Question #683113 on Yade changed: https://answers.launchpad.net/yade/+question/683113 Description changed to: Hi, I looked into the answer for that question hoping that I would get clarified but it seems that the author found out another way of doing the work and the answer got deviated to polyhedral particles. Can you please let me know how to prepare a sample with a cluster of spheres having an irregular surface as a predicate? Is it possible to do with the Triax engine? If so, kindly let me know the steps. Thanks. -- 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
