yade-users team mailing list archive

[Jiannan Wang <question690820@xxxxxxxxxxxxxxxxxxxxx>]

New question #690820 on Yade:
https://answers.launchpad.net/yade/+question/690820

Hello,

I'm trying to inject water into aJCFpm. It is working fine until I enable the VTKRecorder, I got the following error:

python3.6: /build/vtk6-VHOYAG/vtk6-6.3.0+dfsg1/IO/XML/vtkXMLOffsetsManager.h:142: void OffsetsManagerGroup::Allocate(int, int): Assertion `numElements > 0' failed.
Aborted (core dumped)

My Yade version: try both compiled from the source and yadedaily, same error
Linux version: Ubuntu 18.04.4 LTS

Did I miss anything? Thanks.

Best
Jiannan

Here is the script:
#############
from builtins import range
from yade import pack,plot,utils
import math

intR=0.1 
DENS=2500
YOUNG=1800
FRICT=7
ALPHA=0.1
TENS=100
COH=1000
iterper=100

O.materials.append(JCFpmMat(type=1,density=DENS,young=YOUNG,poisson=ALPHA,frictionAngle=radians(FRICT),tensileStrength=TENS,cohesion=COH,label='mat1'))

mn,mx=Vector3(0,0,0),Vector3(1,2,1) # corners of the initial packing
walls=aabbWalls([mn,mx],thickness=0,material='mat1')
wallIds=O.bodies.append(walls)
sps=SpherePack() sp=pack.randomDensePack(pack.inAlignedBox(mn,mx),spheresInCell=2000,radius=1/20.,returnSpherePack=True)
O.bodies.append([sphere(s[0],s[1],color=(0.6+0.15*rand(),0.5+0.15*rand(),0.15+0.15*rand()),material='mat1') for s in sp])

O.engines=[
    ForceResetter(),
    InsertionSortCollider([Bo1_Sphere_Aabb()]),
    InteractionLoop(
      [Ig2_Sphere_Sphere_ScGeom(interactionDetectionFactor=intR)],
     [Ip2_JCFpmMat_JCFpmMat_JCFpmPhys(cohesiveTresholdIteration=1)],
     [Law2_ScGeom_JCFpmPhys_JointedCohesiveFrictionalPM(smoothJoint=True,label='interactionLaw', recordCracks=True, recordMoments=True,Key='identifier',momentRadiusFactor=5,neverErase=True),Law2_ScGeom_FrictPhys_CundallStrack()]
    ),
    DFNFlowEngine(dead=1,label="flow"),
    GlobalStiffnessTimeStepper(active=1,timeStepUpdateInterval=100,timestepSafetyCoefficient=0.5, defaultDt=utils.PWaveTimeStep()),
    NewtonIntegrator(damping=0.1),
]

flow.dead=0
flow.defTolerance=0.3
flow.meshUpdateInterval=1
flow.updateTriangulation=True
flow.breakControlledRemesh=True
flow.imposePressure(Vector3(0.5,1,0),30000000)
O.dt=0.1e-5
O.dynDt=False

O.engines=O.engines+[VTKRecorder(Key='identifier',label='vtk',iterPeriod=iterper,initRun=True,fileName='output-',recorders=['moments','spheres','intr','stress','velocity','bstresses','jcfpm', 'cracks'])]
######################


-- 
You received this question notification because your team yade-users is
an answer contact for Yade.