yade-users team mailing list archive

[xieyihong <question269147@xxxxxxxxxxxxxxxxxxxxx>]

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.