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.
_______________________________________________
Mailing list: https://launchpad.net/~yade-users
Post to : [email protected]
Unsubscribe : https://launchpad.net/~yade-users
More help : https://help.launchpad.net/ListHelp
