Question #269147 on Yade changed:
https://answers.launchpad.net/yade/+question/269147
Status: Answered => Open
xieyihong is still having a problem:
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)
--
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 protected]
Unsubscribe : https://launchpad.net/~yade-users
More help : https://help.launchpad.net/ListHelp
