yade-users team mailing list archive

[Luis Barbosa <question403846@xxxxxxxxxxxxxxxxxxxxx>]

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

Hi,

I'm using Law2_ScGeom_JCFpmPhys_JointedCohesiveFrictionalPM in a simple simulation of two spheres in contact [1].

Increasing the values of poisson (0.2,0.3,0.5,0.7,1) I realized a reduction in max Normal force. My doubt is why it happens, once poisson is supposed to affect Ks only [2].

[1]

#!/usr/bin/python
# -*- coding: utf-8 -*-

from yade import plot

O.materials.append(JCFpmMat(type=1,young=1e8,poisson=0.7,frictionAngle=radians    (30),density=3000,tensileStrength=1e6,cohesion=1e6,jointNormalStiffness=1e6,jointShearStiffness=1e6,jointCohesion=1e6,jointFrictionAngle=radians(20),jointDilationAngle=0.0,label='spheres'))
O.bodies.append([
utils.sphere(center=(0,0,0),radius=1,material='spheres',fixed=True),
utils.sphere(center=(0,0,2),radius=1,material='spheres',fixed=False),
])

O.engines=[
  ForceResetter(),
  InsertionSortCollider([Bo1_Sphere_Aabb()]),
  InteractionLoop(
		[Ig2_Sphere_Sphere_ScGeom()],
[Ip2_JCFpmMat_JCFpmMat_JCFpmPhys(cohesiveTresholdIteration=1,alpha=0.3,tensileStrength=1e6,cohesion=1e6,jointNormalStiffness=1e7,jointShearStiffness=1e7,jointCohesion=1e6,jointFrictionAngle=radians(20),jointDilationAngle=0.0)],
		[Ip2_JCFpmMat_JCFpmMat_JCFpmPhys(cohesiveTresholdIteration=1)],
		[Law2_ScGeom_JCFpmPhys_JointedCohesiveFrictionalPM(smoothJoint=True)]
  ),
  TranslationEngine(ids=[1],translationAxis=[0,1,0],velocity=0.1),
  NewtonIntegrator(damping=0.0),
  PyRunner(iterPeriod=100,command='myAddPlotData()')
]
O.dt=1e-4*utils.PWaveTimeStep()
from yade import qt
qt.View()
qt.Controller()

from yade import plot
plot.labels=dict(Fn='$Normal\,force$',un='$Overlapping$',time='$Time$',time_='$Time$',Fs='$Shear\,force$',t='$Time$')
plot.plots={'time_':('Fn',),'t':('Fs')}
def myAddPlotData():
	if O.interactions[0,1].isReal:
		i=O.interactions[0,1]
		plot.addData(Fn=i.phys.normalForce[2],Fs=i.phys.shearForce[1],time=O.time,time_=O.time,t=O.time)
	else:
		plot.addData(Fn=0,Fs=0,un=0,time=O.time,time_=O.time,t=O.time)
plot.plot(subPlots=False)


[2]http://www.sciencedirect.com/science/article/pii/S1365160912000391


Thanks,
Luis

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