Question #682161 on Yade changed:
https://answers.launchpad.net/yade/+question/682161
JIPEIQI gave more information on the question:
A more confusing situation would be :
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
--
You received this question notification because your team yade-users is
an answer contact for Yade.
_______________________________________________
Mailing list: https://launchpad.net/~yade-users
Post to : [email protected]
Unsubscribe : https://launchpad.net/~yade-users
More help : https://help.launchpad.net/ListHelp
