yade-users team mailing list archive

[ytang <question694208@xxxxxxxxxxxxxxxxxxxxx>]

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

Hi all,

I'm simulating the CPT tests. I use the geom.facetCylinder and geom.facetCone commands to generate the penetrator.
the code shows below:
#############################################################
x0=0.2;y0=0.2;z0=0.39665
cylinderIDS= O.bodies.append(geom.facetCylinder((x0,y0,0.35665),
	radius=0.0125,height=0.15,orientation=Quaternion((1, 0, 0), 0),wallMask=5,segmentsNumber=10, angleRange=None,material=facetMat))
coneIDS= O.bodies.append(geom.facetCone((x0,y0,0.270825),
	radiusTop=0.0125,radiusBottom=0.0,height=0.02165,orientation=Quaternion((1, 0, 0), 0),wallMask=6,segmentsNumber=10, angleRange=None,material=facetMat))
#######################################
I did two tests. one is 10 segments number, another one is 50 segments number for the cylinder and cone.
###############################
for the 10 segments case: the simulation runs without particles blowing up. while for the 50 segments case, the particles blowing up after the penetrator touches the sample surface. 
#########################################
These two simulations are exactly the same except for the segment number. for these two cases, I set the same timestep. 
#######################################################
GlobalStiffnessTimeStepper(active=1,timeStepUpdateInterval=10,timestepSafetyCoefficient=0.8),
###################
###################

I think the particles blowing up may result from the time step. 

Do you have any advice? 
Thanks!



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