apply cyclic shear load

Asked by jsonscript on 2020-02-23

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

Question information

Language:
English Edit question
Status:
Answered
For:
Yade Edit question
Assignee:
No assignee Edit question
Last query:
2020-02-23
Last reply:
2020-02-23
Chareyre (bruno-chareyre-9) said : #1

Hello,
It seems to me that you are posting a random collection of lines, I
wouldn't call this a script, in order to justify a vague question.
Please be consistent.
Bruno

Le dim. 23 févr. 2020 11:24, yangsiyuan <
<email address hidden>> a écrit :

> 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 address hidden>
> Unsubscribe : https://launchpad.net/~yade-users
> More help : https://help.launchpad.net/ListHelp
>

Can you help with this problem?

Provide an answer of your own, or ask jsonscript for more information if necessary.

To post a message you must log in.