yade-users team mailing list archive

[ytang <question700475@xxxxxxxxxxxxxxxxxxxxx>]

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

Hi All,

I hope you are doing well.
####################################################
Recently, I am doing a cone penetration simulation. the following is part of the code.
##############################################################
from yade import ymport
import itertools
from numpy import *
import numpy as np
from yade import pack, plot, export, utils,geom
import math
target=0.12
young=4e8
finalcompFricDegree=19.5
frictDegree = 39.5
############################################# sphere particles material properties ########################
sphereMat=O.materials.append(CohFrictMat(young=young,poisson=0.3,frictionAngle=radians(finalcompFricDegree),isCohesive=False,alphaKr=0.2,alphaKtw=0,etaRoll=0.5,momentRotationLaw=True,density=2648))
O.bodies.append(ymport.text("rotation-loose-sample-Dr-03.txt",material = sphereMat))
################################################ create the penetrator and shaft and the cylinder ##############################
facetMat=O.materials.append(FrictMat(young=young,poisson=0.3,frictionAngle=radians(frictDegree),density=0))
O.materials.append(FrictMat(young=young,poisson=0.3,frictionAngle=0,density=0,label='walls'))
facets =  O.bodies.append(geom.facetCylinder((200e-3,200e-3,300e-3),200e-3,600e-3,wallMask=6,material='walls',segmentsNumber=100))
###############################################  coloring #####################################################
############################ cone and shaft geometry parameter #######################
x0=0.2;y0=0.2;z0=0.39665
####################################################################################
cylinderIDS= O.bodies.append(geom.facetCylinder((x0,y0,0.416651),
	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.3308255),
	radiusTop=0.0125,radiusBottom=0.0,height=0.021651,orientation=Quaternion((1, 0, 0), 0),wallMask=6,segmentsNumber=10, angleRange=None,material=facetMat))
conecheck=O.bodies[-1].id
############################ cone and shaft geometry parameter ############################
###########################################################################################
############################## set the global engine for the simulation####################
O.engines=[
	ForceResetter(),
	InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb()]),
	InteractionLoop(
			[Ig2_Sphere_Sphere_ScGeom6D(),Ig2_Facet_Sphere_ScGeom()],
			[Ip2_FrictMat_FrictMat_FrictPhys(),Ip2_CohFrictMat_CohFrictMat_CohFrictPhys()],
					[Law2_ScGeom_FrictPhys_CundallStrack(),Law2_ScGeom6D_CohFrictPhys_CohesionMoment(useIncrementalForm=True,always_use_moment_law=True,label='cohesiveLaw')]
	),
	GlobalStiffnessTimeStepper(active=1,timeStepUpdateInterval=10,timestepSafetyCoefficient=0.8),
	NewtonIntegrator(gravity=(0,0,-9.81),damping=0.0)]
#####################################################################
#####################################################################

Generally speaking, the friction angle for particles is 19.5 degrees, while the friction angle for the penetrator (cylindyer+cone) is 39.5 degrees.
I also set the friction angle for the penetrator as 29.5 degrees. (I just changed the friction angle for the penetrator, the friction angle for the particle remains the same. )

when I compare the penetration resistance on the cone for the 29.5-degree case and the 39.5-degree case. the penetration resistance on the cone is exactly the same. 
###########################################
my question is:
how does Yade choose the interface friction angle between the interactions when the friction angles were not the same for the particles and the penetrator?


best,

Yong



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