yade-users team mailing list archive

[Amrisha Khandelwal <question706072@xxxxxxxxxxxxxxxxxxxxx>]

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

    Status: Needs information => Open

Amrisha Khandelwal gave more information on the question:
My code is running upto print stress at different iterations. It stopped 
when stress<limitstress (i.e. 1000). After that a wall is generated at 
top but this wall doesn't move as per the given velocity gradient ( 0, 0 
,-0.1). Following the is part of script showing the same.

def checkStress():
  # stress tensor as the sum of normal and shear contributions
  # Matrix3.Zero is the intial value for sum(...)
  stress = getStress().trace()/3.
  print('Stress at', O.iter, '=', stress)
  # if mean stress is below (bigger in absolute value) limitMeanStress, 
start shearing
  if stress < limitMeanStress:
   O.bodies.append(wall(max([b.state.pos[2] + b.shape.radius for b in 
O.bodies if isinstance(b.shape, Sphere)]), axis=2, sense=-1))
   global plate
          plate = O.bodies[-1] # the last particles is the plate
   # plate moves downwards by giving velocity gradient
   plate.state.vel = (0, 0, -.1)
   checker.command='stopplate()'

#Apply oedometric compression by applying velocity over a plate/wall

def stopplate():
  force= O.forces.f(plate.id)[2]
  print('Force on plate=', force)
  # if the force on plate exceeds maximum load, start unloading
  if O.forces.f(plate.id)[2]> 1e5:
   plate.state.vel = (0, 0, 0)
   print('Plate stops and Distortion is applied at 0.1 strain rate')
   #Apply Distorsion
   O.cell.velGrad = Matrix3(0, 0, .1, 0, 0, 0, 0, 0, 0)
   checker.command = 'checkDistorsion()'

def checkDistorsion():
  # if the distorsion value is >.5, exit; otherwise do nothing
  if abs(O.cell.trsf[0, 2]) > .5:
   O.pause()
   print('Distortion of cell exceeds 0.5')



-------- Original Message --------
Subject: Re: [Question #706072]: Simple shear with oedometric loading
Date: 14.04.2023 16:25
 From: Jan Stránský <question706072@xxxxxxxxxxxxxxxxxxxxx>
To: cez228139@xxxxxxxxxx

Your question #706072 on Yade changed:
https://answers.launchpad.net/yade/+question/706072

     Status: Open => Needs information

Jan Stránský requested more information:
Please provide a MWE [1], a complete script.
Cheers
Jan

[1] https://www.yade-dem.org/wiki/Howtoask

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