yade-users team mailing list archive

[yangsiyuan <question688938@xxxxxxxxxxxxxxxxxxxxx>]

New question #688938 on Yade:
https://answers.launchpad.net/yade/+question/688938

Hello Everybody, I‘d like to ask How to apply cyclic shear load in periodic cell  , I'm not really understand it .

The following is part of my code, which only applies unidirectional shear stress, and how to modify the code to change the direction of the shear stress periodically.

time=0
count=0
count1=0
zero=1
stop=0
while 1:
  #O.run(time*10,True)
  O.run(100,True)
  count=count+1
  count1=count1+1
  gamma=O.cell.hSize[0,2]*100/O.cell.hSize[2,2]-g0
  #gamma1=gamma0+O.cell.velGrad[0,1]*25*O.dt
  if count==32:
    #saveDat()
    count=0
 
  stress=utils.getStress()
  sigmazx=stress[0,2]/-1000
  sigmazz=stress[2,2]/1000.0
  if sigmazx>26:
    #saveDat()
    defData(cycle)
    count=0
    count1=0
    O.cell.velGrad[0,2]=0
#    O.run(10000,True)
    O.save('Circular-x25-y0-m.xml.bz2')
    break

  print 'Speed:',O.cell.velGrad[0,2],'time:',O.realtime,'cycle:',cycle,'szx:',sigmazx,'stop:',stop
  if count1==4:
    defData(cycle)
    count1=0


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