yade-users team mailing list archive

[Ferenc Safranyik <question254731@xxxxxxxxxxxxxxxxxxxxx>]

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

    Status: Answered => Open

Ferenc Safranyik is still having a problem:
Hi Jerome,

thanks for the quick help. If I set up the cohesion, I get an error message "Segmentation fault".
Here is my script:


from yade import pack                                                       
from yade import qt	                                                    
qt.View()		                                                    


Alsoresz = O.bodies.append(geom.utils.facetBox((.5,.5,.1),(.5,.5,.1),dynamic=False,wallMask=31))        
Felsoresz = O.bodies.append(geom.utils.facetBox((.5,.5,.5),(.5,.5,.3), wallMask = 15))               


densT = 2600        
ET =  15e4         
nuT = .4           
frictAngT = .52   
nCoh = 1e1        
sCoh = 5e30        
eRoll = .1        
aKr = 2
aKtw = 2


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'))      


felho = pack.SpherePack()                                                            
felho.makeCloud((0,0,0),(1,1,.8),rMean=.03, rRelFuzz=.001)                           
particles = O.bodies.append([sphere(c,r,material='talaj') for c,r in felho])                 


O.engines=[                   
        ForceResetter(),
        InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Wall_Aabb(),Bo1_Facet_Aabb()]),
        InteractionLoop(           
                [Ig2_Sphere_Sphere_ScGeom(),Ig2_Wall_Sphere_ScGeom(), Ig2_Facet_Sphere_ScGeom()],
                [Ip2_FrictMat_FrictMat_FrictPhys(),Ip2_CohFrictMat_CohFrictMat_CohFrictPhys(label='coh')],
                [Law2_ScGeom_FrictPhys_CundallStrack(),Law2_ScGeom6D_CohFrictPhys_CohesionMoment()]
        ),             
        NewtonIntegrator(damping=.1,gravity=(0,0,-9.81),label='newton'),      
]

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


def Atlag():
	print utils.unbalancedForce()
	if O.iter > 1000:
		if utils.unbalancedForce() < .01: 
			coh.setCohesionNow = True
                      

O.engines=O.engines + [PyRunner(command='Atlag()',iterPeriod = 1)]                      
O.run()           


Thanks.

-- 
You received this question notification because you are a member of
yade-users, which is an answer contact for Yade.