yade-users team mailing list archive

[JIPEIQI <question682161@xxxxxxxxxxxxxxxxxxxxx>]

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.