how is shear force related to poisson ratio?
Hi all!
I'm stilling studying parameters in "CpmMat" and "FrictMat" material types. In the posted scipt below, a sphere labeled "s1" is fixed at the origin. Another sphere, "s2", approaches and overlaps "s1", and then moves along the z axis. By doing this I want to know how the shear force varies over time between two spheres.
However, I find that changing the poisson ratio from 1 to 20 has no effect on the forces along the z axis in these simulations for both the "CpmMat" and "FrictMat". It seems that they are merely the verticle components of normal stresses.
Another question is that when poisson ratio equals zero, "s2" with "FrictMat" would immediately disappear when the simulation starts. I wonder how this happens.
from yade import plot
concreteId = O.materials.
young = 1,
frictionAngle = 0,
poisson = 0, #defines poisson ratio
relDuctility = 0,
sigmaT = 0,
epsCrackOnset = 1e-5 #low strain limit,
))
frictId = O.materials.
young=1,
frictionAngle = 0,
poisson = 1,
))
s1 = utils.sphere(center = (0,0,0), radius = 0.5, material = concreteId) #assigns the material type
s2 = utils.sphere(center = (1,0,0), radius = 0.5, material = concreteId)
O.bodies.
s1.state.
O.engines=[
ForceResetter(),
InsertionSortC
InteractionLoop(
[Ig2_
[Ip2_
[Law2_
),
NewtonIntegrat
PyRunner(command = 'addPlotData()', iterPeriod = 50),
PyRunner(
]
O.dt=5e-
def addPlotData():
plot.
fz = O.forces.
pz = s2.state.pos[2], # z coordinate
fx = O.forces.
px = s2.state.pos[0], # x coordinate
)
def Pushing():
s2.state.vel = (-0.1,0,0) #lets s2 overlap s1
if s2.state.pos[0] < 0.9:
ctr.command = "Pushing2()"
def Pushing2():
s2.state.vel = (0,0,0.1) #lets s2 move upward
ctr.command = "Stop()"
def Stop():
if s2.state.pos[2] > 0.5:
O.pause()
plot.plots = {'px' : ('fx'), 'pz':('fz')}
plot.plot()
Thanks!
Question information
- Language:
- English Edit question
- Status:
- Solved
- For:
- Yade Edit question
- Assignee:
- No assignee Edit question
- Solved by:
- Chareyre
- Solved:
- Last query:
- Last reply: