yade-users team mailing list archive

[JIPEIQI <question682161@xxxxxxxxxxxxxxxxxxxxx>]

Question #682161 on Yade changed:
https://answers.launchpad.net/yade/+question/682161

    Status: Answered => Open

JIPEIQI is still having a problem:
Hi Jan. Sorry about the inaccurate script above.  The script indeed goes
fine. But if you add anything to the last line, then it  just goes
wrong! Like this :

dd='2'
radmult=1
savepath='s'
#Law2_ScGeom_JCFpmPhys_JointedCohesiveFrictionalPM(momentRadiusFactor=5.,smoothJoint=True,recordCracks=False,momentsFileExist=False,cracksFileExist=False,label='jcf',neverErase=True)
#print jcf
dd='2'
radmult=1
savepath='s'
O.engines=[
 ForceResetter(),
 InsertionSortCollider([Bo1_Sphere_Aabb(aabbEnlargeFactor=radmult),Bo1_Facet_Aabb(),Bo1_Box_Aabb()]),
 InteractionLoop(
  [Ig2_Sphere_Sphere_ScGeom(interactionDetectionFactor=radmult),Ig2_Facet_Sphere_ScGeom(),Ig2_Box_Sphere_ScGeom()],
  [Ip2_JCFpmMat_JCFpmMat_JCFpmPhys()],
  [Law2_ScGeom_JCFpmPhys_JointedCohesiveFrictionalPM(momentRadiusFactor=2.,smoothJoint=True,recordCracks=True,recordMoments=True,momentsFileExist=False,cracksFileExist=False,label='jcf',Key=dd,neverErase=True)]
 ),
 #GlobalStiffnessTimeStepper(timestepSafetyCoefficient=0.8),
 newton,
 #VTKRecorder(fileName="%suniax"%(savepath),recorders=['all','cracks','jcfpm','moments'],label='vvtk',iterPeriod=5000,Key=dd),
]
print jcf
triax = ThreeDTriaxialEngine( thickness=0, stressControl_1=False,stressControl_2=True, stressControl_3=False, goal1=0, goal2=-123, goal3=0., max_vel=0.1, internalCompaction=False,)
O.run(1)# This this line is commented, then everything works fine

O.engines=[
 ForceResetter(),
 InsertionSortCollider([Bo1_Sphere_Aabb(aabbEnlargeFactor=radmult),Bo1_Facet_Aabb(),Bo1_Box_Aabb()]),
 InteractionLoop(
  [Ig2_Sphere_Sphere_ScGeom(interactionDetectionFactor=radmult),Ig2_Facet_Sphere_ScGeom(),Ig2_Box_Sphere_ScGeom()],
  [Ip2_JCFpmMat_JCFpmMat_JCFpmPhys()],
  [Law2_ScGeom_JCFpmPhys_JointedCohesiveFrictionalPM(momentRadiusFactor=2.,smoothJoint=True,recordCracks=True,recordMoments=True,momentsFileExist=False,cracksFileExist=False,label='jcf2',Key=dd,neverErase=True)]
 ),
 #GlobalStiffnessTimeStepper(timestepSafetyCoefficient=0.8),
 #triax,
 newton,
 VTKRecorder(fileName="%suniax"%(savepath),recorders=['all','cracks','jcfpm','moments'],label='vvtk',iterPeriod=5000,Key=dd),
 PyRunner(command='addPlotData()',iterPeriod=100),
 PyRunner(command='stopIt()',iterPeriod=1000),
 PyRunner(command='saveplot()',iterPeriod=1000),
]
print jcf2
33#  do anthing

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