Increase particle radius
Hi,
I am increasing my particle radius at each 100 iteration. But particle radius is not updating when using body.state.radius to update it. The change is not visible in the GUI.
Best,
Mithu
Here is my code:
#!/usr/bin/env python
#encoding: ascii
# Testing of the Deformation Enginge with Luding Contact Law
# Modified Oedometric Test
# The reference paper [Haustein2017]
from __future__ import print_function
from yade import utils, plot, timing
from yade import pack
import pandas as pd
import numpy as np
from scipy.integrate import odeint
o = Omega()
# Physical parameters
fr = 0.54
rho = 1050
Diameter = 0.0012
D=Diameter
r1 = Diameter/2
#r2 = Diameter/2
k1 = 100000
kp = 12.0*k1
kc = k1 * 0.1
ks = k1 * 0.1
DeltaPMax = Diameter/3.0
Chi1 = 0.34
o.dt = 1.0e-7
particleMass = 4.0/3.0*
Vi1 = math.sqrt(
PhiF1=0.999
#PhiF1 = DeltaPMax*
Tab_rad=0.008
Cyl_height=0.045
cross_area=
Comp_press= 0.5e8
Comp_force=
i=0
# swelling_mode function
def model(r,
Q=(
drdt =((D*rho_
return drdt
P=[1.45,
time_save=[]
radius_save=[]
#******
# Add material
mat1 = O.materials.
# Spheres for compression and walls
sp=pack.
sp.makeCloud(
sp.toSimulation
walls=O.
# Add engines
o.engines = [
ForceResetter(),
InsertionSort
InteractionLoop(
[Ig2_
Ig2_
Ig2_
[Ip2_
[Law2_
),
NewtonIntegra
PyRunner(
#DeformContro
]
def ParticleSwelling():
time_
time_
radius_
for b in O.bodies:
if isinstance(b.shape, Sphere):
r= odeint(
radius_
if time_current>8:
O.pause()
Question information
- Language:
- English Edit question
- Status:
- Answered
- For:
- Yade Edit question
- Assignee:
- No assignee Edit question
- Last query:
- Last reply:
Can you help with this problem?
Provide an answer of your own, or ask Mithushan Soundaranathan for more information if necessary.