Simualtion for the Workprocess of a Shovel
Dear all,
I have tried a workprocess of a shovel. In first 5800 steps work the Simulation stable, but after 6000 steps comes it a worning 'Segmentation fault (core dumped)
wuxin@37l4247e2
and brocken the simultion.
It is the Code:
Part 1,
from yade import utils
from numpy import linspace
from numpy import arange
import gts
import itertools
from yade import geom,pack,qt,plot
## PhysicalParameters
Density=2400
frictionAngle=
frictionAngleSt
frictionAngleBo
tc = 0.001
en = 0.3
es = 0.3
rMean=.0096
rRelFuzz=.0016
rRelFuzz=0.3
VelocityX=0
spheremat=
steel=O.
kwMeshes=
from yade import ymport
sp=pack.
sp.makeCloud(
sp.toSimulation
wallmat=
O.bodies.
O.engines=[
ForceResetter(),
InsertionSortC
InteractionLoop(
[Ig2_
[Ip2_
[Law2_
),
NewtonIntegrato
PyRunner(
qt.SnapshotEngi
]
O.dt=.3*
def Save():
if utils.unbalance
if max([b.
checker.
def Save_Teil1():
O.save(
O.saveTmp()
O.pause()
qt.Controller()
qt.View()
O.run()
Part 2,
O.load(
import math
from yade import utils
from numpy import linspace
from numpy import arange
import gts
import itertools
from yade import geom,pack,qt,plot
## PhysicalParameters
Density=2650
frictionAngle=
frictionAngleSt
frictionAngleBo
tc = 0.001
en = 0.3
es = 0.3
rMean=.0096
rRelFuzz=.0016
rRelFuzz=0.3
VelocityX=0
spheremat=
steel=O.
kwMeshes=
from yade import ymport
ori = Quaternion(
fctIds=
O.engines=[
ForceResetter(),
InsertionSortC
InteractionLoop(
[Ig2_
[Ip2_
[Law2_
),
NewtonIntegrato
PyRunner(
CombinedKinemat
PyRunner(
PyRunner(
qt.SnapshotEngi
]
def Kinematik_
for n in fctIds:
if PX4>0.28:
O.pause()
O.engines=
def RTranslation():
for j in fctIds:
if PX1>0.35:
else:
return
def T1():
T=O.iter
Vx=O.dt*(float (T**(0.5)))*2000.0
Engin.
for l in fctIds:
if PX2<0.3:
Vy=0
return
else:
return
def RT():
Vr=O.dt*(float (T**(0.5)))*900.0
Engin.
def T3():
Engin.
Engin.
O.dt=.3*
yade.plot.
def addplotdata():
Fx1 = 0
Fy1 = 0
Fz1 = 0
for i in fctIds:
Fx1 += O.forces.f(i)[0]
Fy1 += O.forces.f(i)[1]
Fz1 += O.forces.f(i)[2]
plot.
def stopSimulation():
if Engin.comb[
plot.
O.pause()
else:
return
plot.liveInterv
plot.plots=
plot.plot()
O.saveTmp()
qt.Controller()
qt.View()
O.run()
Thanks in advance
Question information
- Language:
- English Edit question
- Status:
- Answered
- For:
- Yade Edit question
- Assignee:
- No assignee Edit question
- Last query:
- Last reply:
Can you help with this problem?
Provide an answer of your own, or ask Xin Wu for more information if necessary.