Question #707278 on Yade changed: https://answers.launchpad.net/yade/+question/707278
Ruidong LI gave more information on the question: Thanks for your reply, Jan! I'd also like to ask whether the function for calculating the stress on the loading plate and confining stress is correct or not: def addPlotData(): # axial stress global wAz global wAc global x_cyl global y_cyl fwt = Vector3(0, 0, 0) for i in idTplate: fwt += O.forces.f(i) fwb = Vector3(0, 0, 0) for j in idBplate: fwb += O.forces.f(j) # cylindrical stress fcyl = 0 for n in nodes: fnode = O.forces.f(n.id) x,y,z = n.state.pos dirV = Vector3(x-x_cyl,y-y_cyl,0).normalized() fnodeMagn = fnode.dot(dirV) fcyl += fnodeMagn plot.addData(z=O.iter, swt=fwt[2]/wAz, swb=fwb[2]/wAz, swc = fcyl/wAc) Many thanks Kyle -- 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 : yade-users@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp
