I am trying to simulate uniaxial compression with DEM on yade. I would like to know how to introduce the hertz mindlin contact model to my model. Also, I need the model to be parallel bonded model, how do I do that in the python script?
I am trying to simulate the uniaxial compression test on the cylinder. I would like to introduce the hertz mindlin contact model, parallel bond bonded model for concrete and also I would like to know how to check the number of elements generated in the geometry in yade.
Also, When I try to save the simulation it appears to save to the directory i provide, but when i open that directory its not saved there. can you also help me with that?
I would really appreciate the response from your side.
Thanks,
Anay
Question information
- Language:
- English Edit question
- Status:
- Solved
- For:
- Yade Edit question
- Assignee:
- No assignee Edit question
- Solved by:
- Luc Scholtès
- Solved:
- Last query:
- Last reply:
|
Revision history for this message
|
#1 |
Hello Anay,
Firstly, please spend some time reading through [1 https:/
To get you started: try starting with an example concrete script [2], and modifying it to what you desire.
>I would like to know how to check the number of elements generated in the geometry in yade.
len(O.bodies)
should do the trick.
>Also, When I try to save the simulation it appears to save to the directory i provide, but when i open that directory >its not saved there. can you also help me with that?
Please be specific with how you save your simulation. In other words, copy and paste the code you use.
Best,
Robert
[1]https:/
[2]https:/
|
Revision history for this message
|
#2 |
Hi Robert,
I am using one of the reference codes I found on github. I am attaching the
code file in this mail. I am really new to yade platform hence, apologies
for asking very basic questions.
I would really appreciate your help with the code. Please have a look at
the code and let me know the necessary changes.
Also, I am trying to implement the parallel bond model for the concrete
cylinder I am using. Would you help with that also?
Thanks,
Anay
On 1 Dec 2017 11:52 am, "Robert Caulk" <email address hidden>
wrote:
> Your question #661250 on Yade changed:
> https:/
>
> Status: Open => Needs information
>
> Robert Caulk requested more information:
> Hello Anay,
>
> Firstly, please spend some time reading through [1
> https:/
> information (e.g. your code, your errors, etc.). Try reformatting the
> question to "I can't seem to [insert problem]. Here is my code [insert
> code] here is the error I am encountering [insert error] here are the ways
> I've tried to solve the problem [insert tries] here are the other launchpad
> questions that seem to pertain to my problem [insert questions]". This way,
> you will likely see far better results.
>
> To get you started: try starting with an example concrete script [2], and
> modifying it to what you desire.
>
> >I would like to know how to check the number of elements generated in
> the geometry in yade.
>
> len(O.bodies)
>
> should do the trick.
>
> >Also, When I try to save the simulation it appears to save to the
> directory i provide, but when i open that directory >its not saved
> there. can you also help me with that?
>
> Please be specific with how you save your simulation. In other words,
> copy and paste the code you use.
>
> Best,
>
> Robert
>
> [1]https:/
> [2]https:/
>
> --
> To answer this request for more information, you can either reply to
> this email or enter your reply at the following page:
> https:/
>
> You received this question notification because you asked the question.
>
|
Revision history for this message
|
#3 |
Hi Robert,
I see that I am using same reference code as you mentioned in the link. I
just have changed geometry and few parameters value. I think this makes it
easier for you to understand my issue.
Regarding the save option, I open the display of the simulation in yade and
hit save button given there. That's doesn't seem to be working in my case.
Again sincere apologies for my silly questions.
Thanks,
Anay
On 1 Dec 2017 11:52 am, "Robert Caulk" <email address hidden>
wrote:
> Your question #661250 on Yade changed:
> https:/
>
> Status: Open => Needs information
>
> Robert Caulk requested more information:
> Hello Anay,
>
> Firstly, please spend some time reading through [1
> https:/
> information (e.g. your code, your errors, etc.). Try reformatting the
> question to "I can't seem to [insert problem]. Here is my code [insert
> code] here is the error I am encountering [insert error] here are the ways
> I've tried to solve the problem [insert tries] here are the other launchpad
> questions that seem to pertain to my problem [insert questions]". This way,
> you will likely see far better results.
>
> To get you started: try starting with an example concrete script [2], and
> modifying it to what you desire.
>
> >I would like to know how to check the number of elements generated in
> the geometry in yade.
>
> len(O.bodies)
>
> should do the trick.
>
> >Also, When I try to save the simulation it appears to save to the
> directory i provide, but when i open that directory >its not saved
> there. can you also help me with that?
>
> Please be specific with how you save your simulation. In other words,
> copy and paste the code you use.
>
> Best,
>
> Robert
>
> [1]https:/
> [2]https:/
>
> --
> To answer this request for more information, you can either reply to
> this email or enter your reply at the following page:
> https:/
>
> You received this question notification because you asked the question.
>
|
Revision history for this message
|
#4 |
Hello Anay,
>I am using one of the reference codes I found on github.
Which code? There are so many on github...feel free to copy and paste links to github, yade documentation, etc. It is very helpful for people trying to help you.
>I am attaching the
>code file in this mail.
You will need to copy and paste your code into your replies. Email attachments do not appear anywhere in launchpad communications.
>I am really new to yade platform hence, apologies
>for asking very basic questions.
Welcome to Yade! Actually, you are very welcome to ask the most basic questions. That is not the problem. I am simply suggesting that it is difficult to help you since you are not providing enough information for anyone to help.
>I would really appreciate your help with the code. Please have a look at
>the code and let me know the necessary changes.
Please also add the problem you are facing. For example "when I try to do [blank], the code does [blank]. Why can't I add/change [blank] without the code throwing the error [blank]?"
|
Revision history for this message
|
#5 |
Hi Robert,
I am using this code to get started with the simulation.
#!/usr/bin/python
# -*- coding: utf-8 -*-
from __future__ import division
from yade import plot,pack,timing
import time, sys, os, copy
#import matplotlib
#matplotlib.
#matplotlib.
"""
A fairly complex script performing uniaxial tension-compression test
on hyperboloid-shaped specimen.
Most parameters of the model (and of the setup) can be read from table
using yade-multi.
After the simulation setup, tension loading is run and stresses are
periodically saved for plotting
as well as checked for getting below the maximum value so far. This
indicates failure (see stopIfDamaged
function). After failure in tension, the original setup is loaded anew
and the sense of loading reversed.
After failure in compression, strain-stress curves are saved via
plot.saveGnuplot and we exit,
giving some useful information like peak stresses in tension/
Running this script for the first time can take long time, as the
specimen is prepared using triaxial
compression. Next time, however, an attempt is made to load
previously-
(from /tmp/triaxPackC
The specimen length can be specified, its diameter is half of the
length and skirt of the hyperboloid is
4/5 of the width.
The particle size is constant and can be specified using the
sphereRadius parameter.
The 3d display has displacement scaling applied, so that the fracture
looks more spectacular. The scale
is 1000 for tension and 100 for compression.
"""
# default parameters or from table
readParamsFromT
young=24e9,
poisson=.2,
sigmaT=3.5e6,
frictionAngle=
epsCrackOnset=
relDuctility=30,
intRadius=1.5,
dtSafety=.8,
damping=0.4,
strainRateTens
strainRateComp
setSpeeds=True,
# 1=tension, 2=compression (ANDed; 3=both)
doModes=3,
specimenLength
sphereRadius=
# isotropic confinement (should be negative)
isoPrestress=0,
)
from yade.params.table import *
if 'description' in O.tags.keys():
O.tags[
# make geom; the dimensions are hard-coded here; could be in param
table if desired
# z-oriented hyperboloid, length 20cm, diameter 10cm, skirt 8cm
# using spheres 7mm of diameter
concreteId=
sps=SpherePack()
sp=pack.
#sp=pack.
sp.toSimulation
bb=uniaxialTest
negIds,
O.dt=dtSafety*
print 'Timestep',O.dt
mm,mx=[pt[axis] for pt in aabbExtrema()]
coord_25,
area_25,
O.engines=[
ForceResetter(),
InsertionSortC
InteractionLoop(
[Ig2_
[Ip2_
[Law2_
),
NewtonIntegrat
CpmStateUpdate
UniaxialStrain
PyRunner(
PyRunner(
]
#O.miscParams=
# plot stresses in ¼, ½ and ¾ if desired as well; too crowded in the
graph that includes confinement, though
plot.plots=
#'sigma.
O.saveTmp(
O.timingEnabled
global mode
mode='tension' if doModes & 1 else 'compression'
def initTest():
global mode
print "init"
if O.iter>0:
O.wait();
O.loadTmp(
print "Reversing plot data"; plot.reverseData()
else: plot.plot(
strainer.
-abs(strainRate
try:
from yade import qt
renderer=
renderer.
except ImportError: pass
print "init done, will now run."
O.step(); # to create initial contacts
# now reset the interaction radius and go ahead
ss2sc.
is2aabb.
O.run()
def stopIfDamaged():
global mode
if O.iter<2 or not plot.data.
the very beginning
sigma,
extremum=
minMaxRatio=0.5 if mode=='tension' else 0.5
if extremum==0: return
import sys; sys.stdout.flush()
if abs(sigma[
if isoPrestress==0 else 5e-2):
if mode=='tension' and doModes & 2: # only if compression is enabled
mode=
O.save(
print "Saved /tmp/uniax-
interaction-
print "Damaged, switching to compression... "; O.pause()
# important! initTest must be launched in a separate thread;
# otherwise O.load would wait for the iteration to finish,
# but it would wait for initTest to return and deadlock would result
import thread; thread.
return
else:
print "Damaged, stopping."
ft,fc=
print 'Strengths fc=%g, ft=%g, |fc/ft|
title=
O.tags['params']
print 'gnuplot'
print 'Bye.'
O.pause()
#sys.exit(0) # results in some threading exception
def addPlotData():
yade.plot.
'sigma.
'sigma.
'sigma.
})
plot.plot(
#O.run()
initTest()
waitIfBatch()
Problem I am facing:
I want to include the parallel bond particle model in this code. How
do I do that?
Also, is there any way to check which contact model and parameters
yade is using in this code?
Thanks,
Anay
On Fri, Dec 1, 2017 at 3:08 PM, Robert Caulk <
<email address hidden>> wrote:
> Your question #661250 on Yade changed:
> https:/
>
> Status: Open => Needs information
>
> Robert Caulk requested more information:
> Hello Anay,
>
> >I am using one of the reference codes I found on github.
>
> Which code? There are so many on github...feel free to copy and paste
> links to github, yade documentation, etc. It is very helpful for people
> trying to help you.
>
> >I am attaching the
> >code file in this mail.
>
> You will need to copy and paste your code into your replies. Email
> attachments do not appear anywhere in launchpad communications.
>
> >I am really new to yade platform hence, apologies
> >for asking very basic questions.
>
> Welcome to Yade! Actually, you are very welcome to ask the most basic
> questions. That is not the problem. I am simply suggesting that it is
> difficult to help you since you are not providing enough information for
> anyone to help.
>
> >I would really appreciate your help with the code. Please have a look at
> >the code and let me know the necessary changes.
>
> Please also add the problem you are facing. For example "when I try to
> do [blank], the code does [blank]. Why can't I add/change [blank]
> without the code throwing the error [blank]?"
>
> --
> To answer this request for more information, you can either reply to
> this email or enter your reply at the following page:
> https:/
>
> You received this question notification because you asked the question.
>
|
Revision history for this message
|
#6 |
>I want to include the parallel bond particle model in this code. How
>do I do that?
Are you referring to [1]? I don't have much experience with CPM, but IINM, CPM is a bonded-particle model, featuring rate dependence, plasticity, and damage. You can investigate CPM further in [2] to see if it is the bonded model you are looking for.
>Also, is there any way to check which contact model and parameters
>yade is using in this code?
These are the model parameters:
young=24e9,
poisson=.2,
sigmaT=3.5e6,
frictionAngle=
epsCrackOnset=
relDuctility=30
concreteId=
more info on these can be found in [3].
>Regarding the save option, I open the display of the simulation in yade and
>hit save button given there. That's doesn't seem to be working in my case.
I recommend using the terminal. The command O.save(
[1]Potyondy, D. O., & Cundall, P. A. (2004). A bonded-particle model for rock. International Journal of Rock Mechanics and Mining Sciences, 41(8 SPEC.ISS.), 1329–1364. http://
[2]http://
[3]https:/
|
Revision history for this message
|
#7 |
Hi Robert,
Apologies for responding so late, I wasn't around the system for the
weekend. Thank you for the help, this will help me with my work. I shall
trouble you if I have further questions. Really appreciate your response
for my queries.
Thanks,
Anay Joshi
On 1 Dec 2017 3:52 pm, "Robert Caulk" <email address hidden>
wrote:
> Your question #661250 on Yade changed:
> https:/
>
> Status: Open => Answered
>
> Robert Caulk proposed the following answer:
> >I want to include the parallel bond particle model in this code. How
> >do I do that?
>
> Are you referring to [1]? I don't have much experience with CPM, but
> IINM, CPM is a bonded-particle model, featuring rate dependence,
> plasticity, and damage. You can investigate CPM further in [2] to see if
> it is the bonded model you are looking for.
>
> >Also, is there any way to check which contact model and parameters
> >yade is using in this code?
>
> These are the model parameters:
>
> young=24e9,
> poisson=.2,
> sigmaT=3.5e6,
> frictionAngle=
> epsCrackOnset=1e-4,
> relDuctility=30
>
> concreteId=
> frictionAngle=
> sigmaT=
> epsCrackOnset,
>
> more info on these can be found in [3].
>
> >Regarding the save option, I open the display of the simulation in yade
> and
> >hit save button given there. That's doesn't seem to be working in my case.
>
> I recommend using the terminal. The command
> O.save(
>
> [1]Potyondy, D. O., & Cundall, P. A. (2004). A bonded-particle model for
> rock. International Journal of Rock Mechanics and Mining Sciences, 41(8
> SPEC.ISS.), 1329–1364. http://
> [2]http://
> [3]https:/
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
> https:/
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https:/
>
> You received this question notification because you asked the question.
>
|
Revision history for this message
|
#8 |
This question was expired because it remained in the 'Open' state without activity for the last 15 days.
|
Revision history for this message
|
#9 |
Hi Robert,
I am running the 'uniax.py' on yade and I am getting the following output
in the terminal. I am seeing a lot of "Damaged, stopping" output. My
question is,
is it normal to see this or am I doing something wrong in the python
script. I am coping my python file as well as the output I am seeing in the
terminal below. Please help me with the issue.
Also, to see how many particles are generated, where in the python code
should I use the syntax len(O.bodies) and how will I print it on the
terminal window?
*Python script:*
#!/usr/bin/python
# -*- coding: utf-8 -*-
from __future__ import division
from yade import plot,pack,timing
import time, sys, os, copy
#import matplotlib
#matplotlib.
#matplotlib.
"""
A fairly complex script performing uniaxial tension-compression test on
hyperboloid-shaped specimen.
Most parameters of the model (and of the setup) can be read from table
using yade-multi.
After the simulation setup, tension loading is run and stresses are
periodically saved for plotting
as well as checked for getting below the maximum value so far. This
indicates failure (see stopIfDamaged
function). After failure in tension, the original setup is loaded anew and
the sense of loading reversed.
After failure in compression, strain-stress curves are saved via
plot.saveGnuplot and we exit,
giving some useful information like peak stresses in tension/
Running this script for the first time can take long time, as the specimen
is prepared using triaxial
compression. Next time, however, an attempt is made to load
previously-
(from /tmp/triaxPackC
The specimen length can be specified, its diameter is half of the length
and skirt of the hyperboloid is
4/5 of the width.
The particle size is constant and can be specified using the sphereRadius
parameter.
The 3d display has displacement scaling applied, so that the fracture looks
more spectacular. The scale
is 1000 for tension and 100 for compression.
"""
# default parameters or from table
readParamsFromT
young=24e9,
poisson=.2,
sigmaT=3.5e6,
frictionAng
epsCrackOns
relDuctilit
intRadius=1.5,
dtSafety=.8,
damping=0.4,
strainRateT
strainRateC
setSpeeds=True,
# 1=tension, 2=compression (ANDed; 3=both)
doModes=2,
specimenLen
sphereRadiu
# isotropic confinement (should be negative)
isoPrestress=0,
)
from yade.params.table import *
if 'description' in O.tags.keys():
O.tags[
rad = 1
# make geom; the dimensions are hard-coded here; could be in param table if
desired
# z-oriented hyperboloid, length 20cm, diameter 10cm, skirt 8cm
# using spheres 7mm of diameter
concreteId=
sps=SpherePack()
sp=pack.
0.1*rad)
#sp=pack.
sp.toSimulation
bb=uniaxialTest
negIds,
O.dt=dtSafety*
print 'Timestep',O.dt
mm,mx=[pt[axis] for pt in aabbExtrema()]
coord_25,
area_25,
O.engines=[
ForceResett
InsertionSortCo
Interaction
[Ig2_Sphere_
),
NewtonInteg
CpmStateUpd
UniaxialStraine
PyRunner(
PyRunner(
]
#O.miscParams=
# plot stresses in ¼, ½ and ¾ if desired as well; too crowded in the graph
that includes confinement, though
plot.plots=
#'sigma.
O.saveTmp(
O.timingEnabled
global mode
mode='tension' if doModes & 1 else 'compression'
def initTest():
global mode
print "init"
if O.iter>0:
O.wait();
print "Reversing plot data"; plot.reverseData()
else: plot.plot(
strainer.
-abs(strainRate
try:
from yade import qt
(100,100,100)
except ImportError: pass
print "init done, will now run."
O.step(); # to create initial contacts
# now reset the interaction radius and go ahead
ss2sc.
is2aabb.
O.run()
def stopIfDamaged():
global mode
if O.iter<2 or not plot.data.
the very beginning
sigma,
extremum=
minMaxRatio=0.5 if mode=='tension' else 0.5
if extremum==0: return
import sys; sys.stdout.flush()
if abs(sigma[
isoPrestress==0 else 5e-2):
if mode=='tension' and doModes & 2: # only if compression is enabled
print "Saved /tmp/uniax-
interaction-
print "Damaged, switching to compression... "; O.pause()
# important! initTest must be launched in a separate thread;
# otherwise O.load would wait for the iteration to finish,
# but it would wait for initTest to return and deadlock would
result
import thread; thread.
return
else:
print "Damaged, stopping."
print 'Strengths fc=%g, ft=%g, |fc/ft|
else O.tags['params']
print 'gnuplot'
print 'Bye.'
def addPlotData():
yade.plot.
'sigma.
'sigma.
'sigma.
})
plot.plot(
#O.run()
initTest()
waitIfBatch()
*Output on the Terminal*
Yade [1]: Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
import yadimport yade.plot; yade.plot.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
-------
ZeroDivisionError Traceback (most recent call last)
/usr/bin/yade in <module>()
/usr/bin/yade in stopIfDamaged()
158 print "Damaged, stopping."
159 ft,fc=max(
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
161 title=O.
'description' in O.tags.keys() else O.tags['params']
162 print
'gnuplot'
ZeroDivisionError: float division by zero
Damaged, stopping.
Strengths fc=-6.43719e+07, ft=133730, |fc/ft|=481.355
gnuplot 20180122T160351
Bye.
Thanks,
-Anay
On Wed, Dec 20, 2017 at 3:38 AM, Launchpad Janitor <
<email address hidden>> wrote:
> Your question #661250 on Yade changed:
> https:/
>
> Status: Open => Expired
>
> Launchpad Janitor expired the question:
> This question was expired because it remained in the 'Open' state
> without activity for the last 15 days.
>
> --
> If you're still having this problem, you can reopen your question either
> by replying to this email or by going to the following page and
> entering more information about your problem:
> https:/
>
> You received this question notification because you asked the question.
>
|
Revision history for this message
|
#10 |
Hello,
Surely, there is no need to post the same error 100s of times to this forum. Please avoid spamming the thread like this in the future.
The error says it all:
--> 160 print 'Strengths fc=%g, ft=%g,
|fc/ft|
ZeroDivisionError: float division by zero
At line 160, it is encountering a float division by zero. You are dividing fc by ft. ft is defined as:
ft = max(sigma)
Since you are doing only compression, your maximum sigma is 0. You can do both modes by changing:
doModes=2 to 3.
Or you can just comment out line 160 if you want to keep doModes=2.
Best,
Robert
|
Revision history for this message
|
#11 |
Hi Robert,
Thanks a lot for the help. I actually couldn't understand why is the error
occuring over and over again so thought it might have some meaning to it. I
will keep it in mind next time.
I really appreciate your help and patience with me.
Thanks,
Anay
On 24 Jan 2018 5:02 pm, "Robert Caulk" <email address hidden>
wrote:
> Your question #661250 on Yade changed:
> https:/
>
> Status: Open => Answered
>
> Robert Caulk proposed the following answer:
> Hello,
>
> Surely, there is no need to post the same error 100s of times to this
> forum. Please avoid spamming the thread like this in the future.
>
> The error says it all:
>
> --> 160 print 'Strengths fc=%g, ft=%g,
> |fc/ft|
>
> ZeroDivisionError: float division by zero
>
> At line 160, it is encountering a float division by zero. You are
> dividing fc by ft. ft is defined as:
>
> ft = max(sigma)
>
> Since you are doing only compression, your maximum sigma is 0. You can
> do both modes by changing:
>
> doModes=2 to 3.
>
> Or you can just comment out line 160 if you want to keep doModes=2.
>
> Best,
>
> Robert
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
> https:/
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https:/
>
> You received this question notification because you asked the question.
>
|
Revision history for this message
|
#12 |
Wow... so many lines.
I would suggest that you use the following first batch of lines to get used to yade first. If you copy paste these lines in a file and name it "uniaxialCompre
Then, once you know how to run and exploit this simulation, you'll be able to modify the script with the contact law that suits you the most. If not, please ask a simple question with minimum lines and clear (short) title.
Luc
--- uniaxial compression:
# -*- coding: utf-8 -*-
from yade import pack, plot
################# SIMULATIONS DEFINED HERE
#### packing (previously constructed)
OUT='compressio
#### Simulation Control
rate=-0.01 #deformation rate
iterMax=10000 # maximum number of iterations
saveVTK=2000 # saving output files for paraview
#### Material microproperties
intR=1.1 # allows near neighbour interaction (can be adjusted for every packing)
DENS=2500 # could be adapted to match material density: dens_DEM=
YOUNG=20e9
FRICT=7
ALPHA=0.1
TENS=1e6
COH=1e6
#### material definition
def sphereMat(): return JCFpmMat(
#### create the specimen
pred=pack.
O.bodies.
R=0
Rmax=0
nbSpheres=0.
for o in O.bodies:
if isinstance(
nbSpheres+=1
R+=o.
if o.shape.
Rmax=
Rmean=R/nbSpheres
print 'nbSpheres=
#### boundary condition (see utils.uniaxialT
bb=utils.
negIds,
################# ENGINES DEFINED HERE
O.engines=[
ForceResetter(),
InteractionLoop(
[Ig2_
[Ip2_
[Law2_
),
UniaxialStrain
GlobalStiffnes
NewtonIntegrat
PyRunner(
]
################# RECORDER DEFINED HERE
def recorder():
yade.
plot.
# if you want to plot during simulation
plot.plots=
#plot.plot()
################# PREPROCESSING
#### manage interaction detection factor during the first timestep and then set default interaction range ((cf. A DEM model for soft and hard rock, Scholtes & Donze, JMPS 2013))
O.step();
### initializes the interaction detection factor
SSgeom.
Saabb.aabbEnlar
#### coordination number verification
numSSlinks=0
numCohesivelinks=0
for i in O.interactions:
if isinstance(
numSSlinks+=1
if i.phys.isCohesive :
numCohesi
print "nblinks=", numSSlinks, " | nbCohesivelinks=", numCohesivelinks, " || Kcohesive=", 2.0*numCohesive
################# SIMULATION REALLY STARTS HERE
strainer.dead=0
O.run(iterMax)
--- plot uniaxial test data:
# -*- coding: utf-8 -*-
from pylab import *
### processing function
def store(var,
data=
it=[]
eps=[]
sig=[]
tc=[]
sc=[]
te=[]
se=[]
ubf=[]
for i in range(0,len(data)):
it.
eps.
sig.
tc.
sc.
te.
se.
ubf.
var.append(it)
var.append(eps)
var.append(sig)
var.append(tc)
var.append(sc)
var.append(te)
var.append(se)
var.append(ubf)
### data input
dataFile1=
a1=[]
store(a1,dataFile1)
rcParams.
lw=2
ms=10
### plots
figure(
ax1=subplot(1,1,1)
xlabel('iteration [-]')
ax1.plot(
ylabel(r'$\sigma_1$ [MPa]')
ax2 = ax1.twinx()
ax2.plot(
ylabel('unbForce [-]')
figure(
ax1=subplot(1,1,1)
xlabel(
plot([x*1e3 for x in a1[1]],[x/1e6 for x in a1[2]],
ylabel(r'$\sigma_1$ [MPa]')
figure(
ax1=subplot(1,1,1)
xlabel(
ax1.plot([x*1e3 for x in a1[1]],[x/1e6 for x in a1[2]],
ylabel(r'$\sigma_1$ [MPa]')
ax2 = ax1.twinx()
ax2.plot([x*1e3 for x in a1[1]],
ax2.plot([x*1e3 for x in a1[1]],
ylabel('cumulative number of microcracks [-]')
legend(
figure(
ax1=subplot(1,1,1)
xlabel(
ax1.plot([x*1e3 for x in a1[1]],[x/1e6 for x in a1[2]],
ylabel(r'$\sigma_1$ [MPa]')
ax2 = ax1.twinx()
ax2.plot([x*1e3 for x in a1[1]],[x+y for x,y in zip(a1[
ylabel('released energy [J]')
### show or save
show()
|
Revision history for this message
|
#13 |
Hi Luc,
Thanks a lot for this descriptive information. I have a few more questions
to ask,
1. I read your paper, ' A DEM model for soft and hard rocks: Role of grain
interlocking on strength' in which you have modeled the granite. However, I
am modeling a concrete cylinder in my work. I have read in several
literatures that they used the bonded particle model (Cundall,2004). My
question is, can I use JCFpm for modeling concrete?
2. I am running the simulation on the remote servers and I would like to
know that is there any way to see the time-dependent animation of the whole
process of my simulation? Also, can you suggest me how to see the
streamlines of the DE with the yade?
Thanks,
Anay
On Wed, Jan 24, 2018 at 9:37 PM, Luc Scholtès <
<email address hidden>> wrote:
> Your question #661250 on Yade changed:
> https:/
>
> Status: Open => Answered
>
> Luc Scholtès proposed the following answer:
> Wow... so many lines.
>
> I would suggest that you use the following first batch of lines to get
> used to yade first. If you copy paste these lines in a file and name it
> "uniaxialCompre
> uniaxialCompres
> with the JCFPM model. Then, if you want to plot some data from the
> simulation, you can copy the second batch of lines, copy paste them in a
> file name "plotUniaxialCo
> plotUniaxialCom
>
> Then, once you know how to run and exploit this simulation, you'll be
> able to modify the script with the contact law that suits you the most.
> If not, please ask a simple question with minimum lines and clear
> (short) title.
>
> Luc
>
> --- uniaxial compression:
>
> # -*- coding: utf-8 -*-
> from yade import pack, plot
>
> ################# SIMULATIONS DEFINED HERE
>
> #### packing (previously constructed)
> OUT='compressio
>
> #### Simulation Control
> rate=-0.01 #deformation rate
> iterMax=10000 # maximum number of iterations
> saveVTK=2000 # saving output files for paraview
>
> #### Material microproperties
> intR=1.1 # allows near neighbour interaction (can be adjusted for every
> packing)
> DENS=2500 # could be adapted to match material density:
> dens_DEM=
> packing porosity as to be computed?
> YOUNG=20e9
> FRICT=7
> ALPHA=0.1
> TENS=1e6
> COH=1e6
>
> #### material definition
> def sphereMat(): return JCFpmMat(
> young=YOUNG,
> tensileStrength
>
> #### create the specimen
> pred=pack.
> O.bodies.
> gap=0.,
>
> R=0
> Rmax=0
> nbSpheres=0.
> for o in O.bodies:
> if isinstance(
> nbSpheres+=1
> R+=o.shape.radius
> if o.shape.
> Rmax=o.shape.radius
> Rmean=R/nbSpheres
>
> print 'nbSpheres=
>
> #### boundary condition (see utils.uniaxialT
> bb=utils.
> negIds,
> bb['posIds'
>
> ################# ENGINES DEFINED HERE
>
> O.engines=[
> ForceResetter(),
> InsertionSortCo
> intR,label=
> InteractionLoop(
> [Ig2_Sphere_
> intR,label=
> [Ip2_JCFpmMat_
> cohesiveTreshol
> [Law2_ScGeom_
> recordCracks=
> ),
> UniaxialStraine
> asymmetry=
> blockDisplaceme
> stopStrain=
> GlobalStiffness
> defaultDt=
> NewtonIntegrato
> PyRunner(
> recorder(
> VTKRecorder(
> fileName=
> Key=OUT,
> ]
>
> ################# RECORDER DEFINED HERE
>
> def recorder():
> yade.plot.
> 'eps':strainer.
> 'sigma'
> 'tc':interactio
> 'sc':interactio
> 'te':interactio
> 'se':interactio
> 'unbF':
> plot.saveDataTx
>
> # if you want to plot during simulation
> plot.plots=
> #plot.plot()
>
> ################# PREPROCESSING
>
> #### manage interaction detection factor during the first timestep and
> then set default interaction range ((cf. A DEM model for soft and hard
> rock, Scholtes & Donze, JMPS 2013))
> O.step();
> ### initializes the interaction detection factor
> SSgeom.
> Saabb.aabbEnlar
>
> #### coordination number verification
> numSSlinks=0
> numCohesivelinks=0
> for i in O.interactions:
> if isinstance(
> isinstance(
> numSSlinks+=1
> if i.phys.isCohesive :
> numCohesivelinks+=1
>
> print "nblinks=", numSSlinks, " | nbCohesivelinks=", numCohesivelinks, "
> || Kcohesive=", 2.0*numCohesive
>
> ################# SIMULATION REALLY STARTS HERE
> strainer.dead=0
> O.run(iterMax)
>
> --- plot uniaxial test data:
>
> # -*- coding: utf-8 -*-
> from pylab import *
>
> ### processing function
> def store(var,
> data=loadtxt(
> it=[]
> eps=[]
> sig=[]
> tc=[]
> sc=[]
> te=[]
> se=[]
> ubf=[]
> for i in range(0,len(data)):
> it.append(
> eps.append(
> sig.append(
> tc.append(
> sc.append(
> te.append(
> se.append(
> ubf.append(
> var.append(it)
> var.append(eps)
> var.append(sig)
> var.append(tc)
> var.append(sc)
> var.append(te)
> var.append(se)
> var.append(ubf)
>
> ### data input
> dataFile1=
>
> a1=[]
> store(a1,dataFile1)
>
> rcParams.
> 20,'axes.
> pad':10,
> lw=2
> ms=10
>
> ### plots
>
> figure(
> ax1=subplot(1,1,1)
> xlabel('iteration [-]')
> ax1.plot(
> ylabel(r'$\sigma_1$ [MPa]')
> ax2 = ax1.twinx()
> ax2.plot(
> ylabel('unbForce [-]')
>
> figure(
> ax1=subplot(1,1,1)
> xlabel(
> plot([x*1e3 for x in a1[1]],[x/1e6 for x in a1[2]],
> ylabel(r'$\sigma_1$ [MPa]')
>
> figure(
> ax1=subplot(1,1,1)
> xlabel(
> ax1.plot([x*1e3 for x in a1[1]],[x/1e6 for x in a1[2]],
> ylabel(r'$\sigma_1$ [MPa]')
> ax2 = ax1.twinx()
> ax2.plot([x*1e3 for x in a1[1]],
> ax2.plot([x*1e3 for x in a1[1]],
> ylabel('cumulative number of microcracks [-]')
> legend(
>
> figure(
> ax1=subplot(1,1,1)
> xlabel(
> ax1.plot([x*1e3 for x in a1[1]],[x/1e6 for x in a1[2]],
> ylabel(r'$\sigma_1$ [MPa]')
> ax2 = ax1.twinx()
> ax2.plot([x*1e3 for x in a1[1]],[x+y for x,y in zip(a1[
> linewidth=lw)
> ylabel('released energy [J]')
>
> ### show or save
> show()
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
> https:/
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https:/
>
> You received this question notification because you asked the question.
>
|
Revision history for this message
|
#14 |
Thanks Luc Scholtès, that solved my question.
