Unhandled exception in thread started by <function liveUpdate at 0x7f0388a37320>
Hi,
I'm trying to simulation triaxial test by using https:/
In [1]: Unhandled exception in thread started by <function liveUpdate at 0x7f0388a37320>
-------
ValueError Traceback (most recent call last)
/usr/lib/
506 for ax in axes:
507 try:
--> 508 ax.relim() # recompute axes limits
509 ax.autoscale_view()
510 except RuntimeError: pass # happens if data are being updated and have not the same dimension at the very moment
/usr/lib/
1936 for line in self.lines:
1937 if not visible_only or line.get_visible():
-> 1938 self._update_
1939
1940 for p in self.patches:
/usr/lib/
1799 Figures out the data limit of the given line, updating self.dataLim.
1800 """
-> 1801 path = line.get_path()
1802 if path.vertices.size == 0:
1803 return
/usr/lib/
955 """
956 if self._invalidy or self._invalidx:
--> 957 self.recache()
958 return self._path
959
/usr/lib/
665 y = self._y
666
--> 667 self._xy = np.column_
668 self._x, self._y = self._xy.T # views
669
/usr/lib/
247 args = [np.array(_m, copy=False, subok=subok) for _m in args]
248
--> 249 shape = _broadcast_
250
251 if all(array.shape == shape for array in args):
/usr/lib/
182 # use the old-iterator because np.nditer does not handle size 0 arrays
183 # consistently
--> 184 b = np.broadcast(
185 # unfortunately, it cannot handle 32 or more arguments directly
186 for pos in range(32, len(args), 31):
ValueError: shape mismatch: objects cannot be broadcast to a single shape
axial strain reach Limitation, stop loading!
Has anybody met the same problem? How can I solve it ?
regards
Cloud
The whole code :
from yade import pack
num_spheres=5000
key='_My_triax_'
targetPorosity=0.40
compFricDegree=20
finalFricDegree=30
rate=-0.02
damp=0.2
stabilityThresh
young=500e6
mn,mx=Vector3(
confiningStress
O.materials.
O.materials.
walls=aabbWalls
wallIds=
sp=pack.
sp.makeCloud(
O.bodies.
triax=TriaxialS
maxMultipli
finalMaxMul
thickness = 0,
stressMask = 7,
internalCom
)
O.engines=[
ForceResett
InsertionSo
Interaction
),
GlobalStiff
triax,
TriaxialSta
NewtonInteg
]
Gl1_Sphere.
yade.qt.
triax.goal1=
while 1:
O.run(1000, True)
unb=
print 'unbalanced force:',unb,' mean stress: ',triax.meanStress
if unb < stabilityThreshold and abs(-confiningS
break
O.save(
print "### Isotropic state saved ###"
import sys
while triax.porosity > targetPorosity:
compFricDegree = 0.95*compFricDegree
setContactF
print "\r Friction: ", compFricDegree,
sys.
O.run(500,1)
O.save(
print "### Compacted state saved ###"
triax.internalC
setContactFrict
triax.stressMask = 3
triax.goal3 = rate
triax.goal1=
triax.goal2=
newton.damping=0.1
O.saveTmp()
from yade import plot
def history():
plot.addData(
devi = -triax.
def stoploading():
axial_strain = abs(triax.
if axial_strain > 0.2:
print 'axial strain reach Limitation, stop loading!'
O.pause()
if 1:
O.engines=
O.engines=
else:
O.engines[
O.run(100,True)
plot.plots=
}
plot.plot()
plot.saveDataTx
plot.saveGnuplo
Question information
- Language:
- English Edit question
- Status:
- Solved
- For:
- Yade Edit question
- Assignee:
- No assignee Edit question
- Solved by:
- Cloud
- Solved:
- Last query:
- Last reply: