cannot get contact forces from my code
Hello all and thank you for your help in advanced.
I am using Yade-2019-
I have coded a Trixial test as bellow and have some questions, but the most important one is the 5th question.
1- When I run my code, no plot will show up however I coded plotting commands
2- The strange problem is when I try to modify my code a little, I get nothing (or a filled 1-color page) as the 3D show of the specimen.
3- The timp steps are so low and run will take a long time (more than 1 hour)
4- I cannot save the dada as a text file that I coded its command at the end of my code
5- I do not know how can I get micro information such as contact forces through my code
6- I do not see on the screen the printed test text I used in my code.
7- please let me know if I have any other problem with my code
Also, I get these warnings when I run my code:
This is my code:
In [1]: <WARNING> InsertionSortCo
<WARNING> Shop:445 static Real Shop::PWaveTime
from yade import pack, plot
sp=pack.
mn, mx=Vector3(0,0,0), Vector3(10,10,10)
sp.makeCloud(
sp.toSimulation()
O.materials.
O.materials.
walls=aabbWalls
wallIds=
O.bodies.
triax=TriaxialC
wall_bottom_
wall_top_
wall_left_
wall_right_
wall_back_
wall_front_
internalCompac
sigmaIsoCompac
sigmaLateralCo
max_vel=10,
strainRate=0.01,
label="triax",
)
O.engines=[
ForceResetter(),
InsertionSortC
InteractionLoop(
[Ig2_Sphere_
[Ip2_FrictMat_
[Law2_
),
GlobalStiffnes
triax,
NewtonIntegrat
PyRunner(
PyRunner(
PyRunner(
]
O.dt=.5*
def history():
plot.addData(
e11=-triax.
s11=-triax.
def addPlotData():
plot.addData(
s11=-
e11=-
Etot=
O.trackEnergy=True
def checkUnbalanced():
if unbalancedForce
O.pause()
plot.
print "****** TEST TEXT ******"
plot.plots=
Gl1_Sphere.
O.run()
plot.plot()
O.saveTmp()
plot.saveDataTx
Thank you for your help.
Question information
- Language:
- English Edit question
- Status:
- Invalid
- For:
- Yade Edit question
- Assignee:
- No assignee Edit question
- Solved by:
- Jérôme Duriez
- Solved:
- Last query:
- Last reply: