YADE Crash - SIGSEGV/SIGABRT handler called
Hi YADE community,
I have a problem as the question's title said.
The procedure that I have done:
1. Packing generation
2. Save the result by using O.save()
3. Reload the saved yade.gz file by using O.load()
4. Working around with the material and interaction and got the crash
Mini working scripts:
--- [Step 1] ---
### INPUT for Material ###
young0= 30e9
poisson0= 0.2
density0 = 2530
compFricDegree=15
### PARAMETERs of particle packing ###
mn=Vector3(
mx=Vector3(
rMean0=4e-4
rRelFuzz0=5./100
porosity0=0.4
distributeMass0
seed0=1
### INPUT for mechanical loading ###
maxVel0=5e-7
confiningS= -2.725e5
### defining MATERIALs ###
O.materials.
O.materials.
### creating the boundary walls ###
walls=aabbWalls
wallIds=
### generating spheres ###
sp=pack.
sp.makeCloud(
O.bodies.
### mechanical loading #### To confer a solid-like nature to the packing
triax=TriaxialS
internalCompac
goal1=confiningS,
goal2=confiningS,
goal3=confiningS,
max_vel=maxVel0,
label="triax"
)
triax.maxMultip
triax.finalMaxM
newton=
O.engines=[
ForceResetter(),
InsertionSortC
InteractionLoop(
[Ig2_
[Ip2_
[Law2_
useIncrement
),
GlobalStiffnes
triax,
newton
]
### --- 1st stage, particle growing while boundary walls are fixed -- ###
NumDEMRun1=100000
while 1:
O.run(
unb=unbalanced
if unb<0.01 and abs(triax.
break
### --- 2nd stage, particle sizes are fixed while boundary walls move very slowly -- ###
finalFricDegree = 30 # contact friction during the deviatoric loading
triax.internalC
setContactFrict
NumDEMRun2=50000
while 1:
O.run(NumDEMRun2 ,True)
unb=unbalanced
if unb<0.001 and abs(triax.
break
*******
--- [Step 4] ---
compFricDegree = 1 # Scholtes
young0= 30e9
poisson0= 0.2
density0 = 2530
tensileStrength = 40e5
shearStrength = 40e6
### setup the material properties, first, Cohesion material
O.materials.
### --- change material of the spheres to create cohesive bond between particles --- ###
pairs = [(i.id1,i.id2) for i in O.interactions]
for b in O.bodies:
if not isinstance(
b.mat = O.materials[2]
### [0] is Friction spheres, [1] is frictionlessWalls, [2] is the newly added CohFrict sphere
O.interactions.
for id1,id2 in pairs: utils.createInt
*******
Any help or suggestion is highly appreciate!
Son Pham
Question information
- Language:
- English Edit question
- Status:
- Expired
- For:
- Yade Edit question
- Assignee:
- No assignee Edit question
- Last query:
- Last reply: