yade-dev team mailing list archive
-
yade-dev team
-
Mailing list archive
-
Message #12954
[Bug 1652529] Re: segmentation fault (core dumped)
Hello,
I can not reproduce the bug: I've tried this example [*] (since it is a
file from "examples" folder, you may just refer directly to it, so that
we're sure it has not been altered) with two recent versions (2016-10-31
.git-1b22b1a and the last one 2016-12-20.git-54c46f3) and did not get
any crash.
I'm marking thus the bug for now as "invalid", before any new
information. You may give us copy of the whole YADE output when you
launch this script, for instance.
Jerome
[*]
https://github.com/yade/trunk/blob/master/examples/capillaryLaplaceYoung
/CapillaryPhys-example.py
** Changed in: yade
Status: New => Invalid
--
You received this bug notification because you are a member of Yade
developers, which is subscribed to Yade.
https://bugs.launchpad.net/bugs/1652529
Title:
segmentation fault (core dumped)
Status in Yade:
Invalid
Bug description:
Hi I am trying to run this code, which is an example given in the
package. But every time i try to execute this file, i get segmentation
fault error. I am using Ubuntu 16.04 LTS and Yade 2016.06a
#!/usr/bin/python
# -*- coding: utf-8 -*-
"""
To run this script you need to have all 10 text files from https://yade-dem.org/wiki/CapillaryTriaxialTest
in the same folder as you run this script in console!
This script shows how to use Law2_ScGeom_CapillaryPhys_Capillarity. The user can switch between hertz and
linear model by setting "model_type"
"""
model_type = 1 #1=Hertz model with capillary forces,
else linear model with capillary model
#some parameters:
shear_modulus = 1e5
poisson_ratio = 0.3
young_modulus = 2*shear_modulus*(1+poisson_ratio)
friction = 0.5
angle = atan(friction)
local_damping = 0.01
viscous_normal = 0.021
viscous_shear = 0.8*viscous_normal
lowercorner = Vector3(0,0,0)
uppercorner = Vector3(0.002,0.002,0.004)
#creating a material (FrictMat):
id_SphereMat=O.materials.append(FrictMat(young=young_modulus,poisson=poisson_ratio,density=2500,frictionAngle=angle))
SphereMat=O.materials[id_SphereMat]
#generate particles:
sp=pack.SpherePack()
sp.makeCloud(lowercorner,uppercorner,.0002,rRelFuzz=.3)
O.bodies.append([sphere(c,r,material=SphereMat) for c,r in sp])
#generate boundary:
O.bodies.append(geom.facetBox(uppercorner/2,uppercorner/2,wire=True,fixed=True,material=SphereMat))
#define engines:
if model_type == 1:#hertz model with capillary forces
O.engines=[
ForceResetter(),
InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb()]),
InteractionLoop(
[Ig2_Sphere_Sphere_ScGeom(),Ig2_Facet_Sphere_ScGeom()],
[Ip2_FrictMat_FrictMat_MindlinCapillaryPhys(label='ContactModel')],#for hertz model only
[Law2_ScGeom_MindlinPhys_Mindlin()]#for hertz model only
),
Law2_ScGeom_CapillaryPhys_Capillarity(capillaryPressure=10000),#for hertz model only
NewtonIntegrator(damping=local_damping,gravity=(0,0,-9.81)),
]
ContactModel.betan=viscous_normal
ContactModel.betas=viscous_shear
ContactModel.useDamping=True
else:
O.engines=[
ForceResetter(),
InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb()]),
InteractionLoop(
[Ig2_Sphere_Sphere_ScGeom(),Ig2_Facet_Sphere_ScGeom()],
[Ip2_FrictMat_FrictMat_CapillaryPhys()], #for linear model only
[Law2_ScGeom_FrictPhys_CundallStrack()], #for linear model only
),
Law2_ScGeom_CapillaryPhys_Capillarity(capillaryPressure=10000),#for linear model only
NewtonIntegrator(damping=local_damping,gravity=(0,0,-9.81)),
]
#set time step and run simulation:
O.dt=0.5*PWaveTimeStep()
from yade import qt
qt.View()
print('Press PLAY button')
#O.run(10000,True)
looking forward to hear soon..!!
Cheers
Amiya
To manage notifications about this bug go to:
https://bugs.launchpad.net/yade/+bug/1652529/+subscriptions
References