Yade

true triaxial test

Asked by xieyihong

Hi
   Can we use the yade to do some research about the true Triaxial test? Thank you!

Question information

Language:
English Edit question
Status:
Solved
For:
Yade Edit question
Assignee:
No assignee Edit question
Solved by:
Jan Stránský
Solved:
Last query:
Last reply:
Revision history for this message
Best Jan Stránský (honzik) said : 		#1

Hi,
yes, we can
cheers
Jan

PS: to get more precise answer, please be more specific on "do some
research" and "true Triaxial test" terms :-)

2015-07-12 4:16 GMT+02:00 xieyihong <email address hidden>:

> New question #269147 on Yade:
> https://answers.launchpad.net/yade/+question/269147
>
> Hi
> Can we use the yade to do some research about the true Triaxial test?
> Thank you!
>
> --
> You received this question notification because you are a member of
> yade-users, which 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
>

Revision history for this message
xieyihong (1501740992-p) said : 		#2

Hi,
Thank you very much. I also think it can be. but I found no thing when I use the "true triaxial test"terms. I wrote the code by the following terms.But the result is not the true triaxial test. Could you help me? thank you very much.

def servo():
 gain=0.000002
 tol=0.00001
 sx_ini=utils.getStress()[0,0]
 sy_ini=utils.getStress()[1,1]
 sz_ini=utils.getStress()[2,2]
 sr_x=-.1
 sr_y=0.
 sr_z=0.
 sx_curr=utils.getStress()[0,0]
 sy_curr=utils.getStress()[1,1]
 sz_curr=utils.getStress()[2,2]
 sy_req=sy_ini-0.28571429*(sx_curr-sx_ini)
 sz_req=sz_ini-0.71428571*(sx_curr-sx_ini)
 sry_curr=O.cell.velGrad[1,1]
 srz_curr=O.cell.velGrad[2,2]
 if fabs(sy_curr-sy_req)>=tol:
  sr_y=sry_curr+gain*(sy_req-sy_curr)
 else:
  sr_y=sry_curr
 if fabs(sz_curr-sz_req)>=tol:
  sr_z=srz_curr+gain*(sz_req-sz_curr)
 else:
  sr_z=srz_curr
 O.cell.velGrad=Matrix3(sr_x,0,0,0,sr_y,0,0,0,sr_z)

Revision history for this message
xieyihong (1501740992-p) said : 		#3

Thanks Jan Stránský, that solved my question.

Revision history for this message
Jan Stránský (honzik) said : 		#4

Hello,
you can open a new question, including a working script and with more
detailed description about what results you get and why they are wrong.
Cheers
Jan

2015-07-12 9:16 GMT+02:00 xieyihong <email address hidden>:

> Question #269147 on Yade changed:
> https://answers.launchpad.net/yade/+question/269147
>
> Status: Open => Solved
>
> xieyihong confirmed that the question is solved:
> Thanks Jan Stránský, that solved my question.
>
> --
> You received this question notification because you are a member of
> yade-users, which 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
>