New question #680397 on Yade:
https://answers.launchpad.net/yade/+question/680397
This was an issue I ran into before with the Cundall-Strack law, and I am now
also seeing with Hertz-Mindlin law. Hence, I am trying to input specific value
for kn and ks for this contact model and the formulation proposed by yade, but
I do not know what I can consider them in my simulation.
Is there a way to use these values? Would a simple
Ip2_FrictMat_FrictMat_MindlinPhys(kn=5e5, ks=5e5)?
You can find herein the necessary part of my script.
Thanks in advance.
Sam
############################## material new parameters
###############################
young =3e5 # elastic modulus [Pa]
poisson = 0.25
finalFricDegree=30
density = 1950
mn,mx=Vector3(-0.52e-3, -1e-3,-1e-3),Vector3(0.52e-3, 1e-3, 1e-3) #size of
particle
radius= 1 # Internal Radius
identifier=str(particleDiameter)
output='./out/'+identifier
rate = strainRate*10000
stabilityThreshold=0.01
iterper=1000 # Periodicity criterion using step number (deactivated if <= 0)
kn=5e1 # Normal stiffness
ks=5e1 # Shear Stiffness
kr=0 # Rotational stiffness
Bn=0
Bs=0
intRadius=1
##########################################################################################
Friction=O.materials.append(FrictMat(young=young, density=density,
frictionAngle=radians(finalFricDegree), poisson=poisson, label='sphere'))
#frictionless walls
wallmat = O.materials.append(FrictMat(young=80e9, poisson=.45,
frictionAngle=radians(finalFricDegree), density=7000, label='Walls'))
walls=aabbWalls([mn,mx], material='Walls')
wallIds=O.bodies.append(walls)
sp = O.bodies.append(ymport.textExt('sandGrain_0-6mmdiameter.spheres',
'x_y_z_r',color=(0,0.2,0.7), material='sphere'))
triax = TriaxialStressController(
thickness=0,
stressMask=7,
internalCompaction=False,
)
O.engines=[
ForceResetter(),
InsertionSortCollider([Bo1_Sphere_Aabb(aabbEnlargeFactor=intRadius),Bo1_Box_Aabb()]),
InteractionLoop(
[Ig2_Sphere_Sphere_ScGeom6D(),Ig2_Box_Sphere_ScGeom()],
[Ip2_FrictMat_FrictMat_MindlinPhys(betan=Bn,betas=Bs,ktwist=0.0,krot=ks,label='ContactModel')],
[Law2_ScGeom_MindlinPhys_Mindlin(neverErase=True,includeMoment=False,includeAdhesion=False,label='Mindlin')]
),
GlobalStiffnessTimeStepper(active=1,timeStepUpdateInterval=100,timestepSafetyCoefficient=0.5),
triax,
VTKRecorder(dead=0,Key=identifier,iterPeriod=iterper,initRun=True,fileName=(output+'-'),recorders=['spheres','intr','stress',
'bstresses','jcfpm','cracks']),
NewtonIntegrator(damping=0.9)
]
O.dt = 0.3 * utils.PWaveTimeStep()
--
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 : yade-users@lists.launchpad.net
Unsubscribe : https://launchpad.net/~yade-users
More help : https://help.launchpad.net/ListHelp
