yade-users team mailing list archive

[Ferenc Safranyik <question678102@xxxxxxxxxxxxxxxxxxxxx>]

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

Dear All,

I try to make a simply simulation of a very cohesive particulate solid. To this CohFrictMat material model was used. After unbalancedForce was decreased I try to turn on cohesion on every contact but I get an erre: "Segmentation fault (core dumped)" and I can not recognise until this time the problem. Can anybody help to me? Thanks in advance!

Some piece of the code:

densT = 2600        
ET =  5e5 * 25       
nuT = .9          
frictAngT = .2      
nCoh = 1e8          
sCoh = 1e8          
eRoll = .1    

talaj = O.materials.append(CohFrictMat(young=ET,poisson=nuT,density=densT,frictionAngle=frictAngT,normalCohesion=nCoh,
                                     shearCohesion=sCoh,isCohesive=True,momentRotationLaw=True,etaRoll=eRoll,label='talaj'))
acel = O.materials.append(FrictMat(young=EA,poisson=nuA,density=densA,frictionAngle=frictAngA,label='acel')) 

doboz = O.bodies.append(geom.utils.facetBox((1.25,0,.05),(.25,.2,.15),wallMask=63, color=(.2,.2,.0),material = 'acel'))

Sugar1 = .02  
felho = pack.SpherePack()
felho.makeCloud((1,-.2,-.1),(1.5,.2,.2),rMean=Sugar1)
particles = O.bodies.append([sphere(c,r,material='talaj',color=(.3,.8,.1)) for c,r in felho])

O.engines=[                   
        ForceResetter(),
        InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Box_Aabb(),Bo1_Facet_Aabb()]),
        InteractionLoop(           
                [Ig2_Sphere_Sphere_ScGeom(),Ig2_Facet_Sphere_ScGeom(),Ig2_Box_Sphere_ScGeom()],
                [Ip2_FrictMat_FrictMat_FrictPhys(),Ip2_CohFrictMat_CohFrictMat_CohFrictPhys(label='coh')],
                [Law2_ScGeom_FrictPhys_CundallStrack(),Law2_ScGeom6D_CohFrictPhys_CohesionMoment(
			useIncrementalForm=True,
			always_use_moment_law = True)]
        ),       
        NewtonIntegrator(damping=.4,gravity=(0,0,-9.81),label='newton'),      
]

O.dt=.01*utils.PWaveTimeStep()

O.engines=O.engines + [PyRunner(command='ment()',iterPeriod = 1,label='ment')]

def ment():
	plot.addData(i=O.iter,uF=utils.unbalancedForce())
	if O.iter > 80000 and utils.unbalancedForce() < .05:
                coh.setCohesionNow=True
      


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