yade-users team mailing list archive

[chanaka Udaya <question685502@xxxxxxxxxxxxxxxxxxxxx>]

Question #685502 on Yade changed:
https://answers.launchpad.net/yade/+question/685502

chanaka Udaya posted a new comment:
Hi Jan,

I have changed the line below in source code and after compiling ran the
following script.

ConcretePM.hpp,(line 281)

((Real,omegaThreshold,((void)">=1. to deactivate, i.e. never delete any
contacts",0.9999),,"damage after which the contact disappears (<1),
since omega reaches 1 only for strain →+∞"))


here is the my script.

from yade import plot,qt
nh=O.materials.append(CpmMat(young=20e9,frictionAngle=atan(0.5),poisson=.5,density=2600,sigmaT=1.5e6,relDuctility=1.5,epsCrackOnset=5e-5,isoPrestress=0))
s1=utils.sphere((0,0,0),radius=1,color=[1,0,0],material=nh)
s2=utils.sphere((0,0,2),radius=1,color=[0,1,0],material=nh)

O.bodies.append(s1)
O.bodies.append(s2)

Gl1_Sphere.quality=3
# engines definition
O.engines=[
    ForceResetter(),
    InsertionSortCollider([Bo1_Sphere_Aabb(aabbEnlargeFactor=1.3,label='bo1s')]),
    InteractionLoop(
        [Ig2_Sphere_Sphere_ScGeom(interactionDetectionFactor=1.3,label='ig2ss')],
        [Ip2_CpmMat_CpmMat_CpmPhys()],
        [Law2_ScGeom_CpmPhys_Cpm()]
    ),
    TranslationEngine(ids=[s2.id],translationAxis=[0,0,1],velocity=0.01*100), #shear+tension
    TranslationEngine(ids=[s1.id],translationAxis=[0,0,-1],velocity=0.01*100),
    NewtonIntegrator(damping=0.4),
    PyRunner(iterPeriod=1,command='addPlotData()',label='plotData'),
]
O.dt=0.004*PWaveTimeStep()
O.step()
bo1s.aabbEnlargeFactor=ig2ss.interactionDetectionFactor=1.
def addPlotData():
    plot.addData(i=O.iter,
                 displacement_n=O.bodies[1].state.refPos[2]-O.bodies[1].state.pos[2],
                 NF=O.interactions[0,1].phys.normalForce.norm(),
                 SF=O.interactions[0,1].phys.shearForce.norm(),
                 )

plot.plots={'i':('NF'),
           'displacement_n ':('NF'),
	   'i  ':('SF'),
	  }
plot.plot()
O.run()

Thanks

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