[Yade-users] [Question #699044]: 如何设定墙体刚度？
New question #699044 on Yade: https://answers.launchpad.net/yade/+question/699044 What code do I use to set the normal and tangential stiffness of the wall and the parallel bonding parameters? -- You received this question notification because your team yade-users is an answer contact for Yade. ___ Mailing list: https://launchpad.net/~yade-users Post to : yade-users@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp
Re: [Yade-users] [Question #698948]: Stability of advection modeling by ThermalEngine
Question #698948 on Yade changed: https://answers.launchpad.net/yade/+question/698948 Robert Caulk proposed the following answer: Hey Zoheir, I have no way of testing that yadedaily version, and it is not recommended to be compiling old branches of yade. That was a branch that was merged into the master, so it is best to stick to the master as you get updates. But, I can confirm that running the latest master version is stable in all the condtions you asked me to test. Please update your sources and recompile: cd ../trunk git remote update git pull upstream master cd ../build make -jX install # X being the numberof cores on your computer Let me know if this fixes it, sorry for the delay, Robert -- You received this question notification because your team yade-users is an answer contact for Yade. ___ Mailing list: https://launchpad.net/~yade-users Post to : yade-users@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp
Re: [Yade-users] [Question #698948]: Stability of advection modeling by ThermalEngine
Question #698948 on Yade changed: https://answers.launchpad.net/yade/+question/698948 Zoheir Khademian posted a new comment: Hey Robert, I ran it in both Yadedaily 20201215-4510~714a723~bionic1 Using python version: 3.6.9 and "eraseBodyThermalEngine" [1] trunk version. Of 10 times of running the model, I get 8 stable results and 2 unstable results (after 1 iterations) If I change the sphere size to rMean=0.001 from 0.0015, all tries lead to instability. Please reduce the mean sphere size to 0.001 and let me know if you can still get stable results. [1]https://gitlab.com/yade-dev/trunk/-/tree/eraseBodyThermalEngine -- You received this question notification because your team yade-users is an answer contact for Yade. ___ Mailing list: https://launchpad.net/~yade-users Post to : yade-users@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp
Re: [Yade-users] [Question #698970]: Questions about updating yade
Question #698970 on Yade changed: https://answers.launchpad.net/yade/+question/698970 Status: Answered => Solved Ziyu Wang confirmed that the question is solved: Thanks a lot for your help! I will try to solve this problem, or I will learn source-code... Thanks again.That solved my problem. -- You received this question notification because your team yade-users is an answer contact for Yade. ___ Mailing list: https://launchpad.net/~yade-users Post to : yade-users@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp
Re: [Yade-users] [Question #698948]: Stability of advection modeling by ThermalEngine
Question #698948 on Yade changed: https://answers.launchpad.net/yade/+question/698948 Robert Caulk requested more information: Hey zoheir, I just ran your MWE out 160 seconds, 170k iterations. It is stable, temperatures do not exceed the boundary conditions. Am I missing something? Which version of yade are you running? Cheers, Robert -- You received this question notification because your team yade-users is an answer contact for Yade. ___ Mailing list: https://launchpad.net/~yade-users Post to : yade-users@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp
Re: [Yade-users] [Question #698948]: Stability of advection modeling by ThermalEngine
Question #698948 on Yade changed:
https://answers.launchpad.net/yade/+question/698948
Zoheir Khademian posted a new comment:
Thanks Robert!
I went through a trial and error process using a more traditional BC as you
suggested:
low.bndCondIsTemperature=[0,0,0,0,1,1]
flow.thermalBndCondValue=[0,0,0,0,25,45]
thermal.bndCondIsTemperature=[0,0,0,0,1,1]
thermal.thermalBndCondValue=[0,0,0,0,25,45]
and avoid mesh updating by meshUpdateInterval=-1
- First, I turned off the compaction and ran the model. Temperature
ranges seemed to be within expected range 25-45 K. Turning on the
compaction made the model unstable in a way that solid and fluid temp
goes above the assigned boundary condition (45 K).
- I tried to minimize the particle-pore and pore-pore resistivities by
assigning very small numbers to thermal.fluidConductionAreaFactor[1] and
thermal.fluidK[2]. I also made sure that the distance between cells wont be
zero by assigning thermal.minimumFluidCondDist[3] of 0.001. I tried different
numbers but the particle temp still goes up to 1e3 very quickly while cell temp
stays below 45 K.
- I also tried avoiding fictious cells by
thermal.ignoreFictiousConduction[4]=True but still no chance. I also tried
using HertzMethod [5] for area calculation. The cell temp was in the range but
particle temp went to negative.
- I also limited the Reynolds (thermal.uniformReynolds[6]) to 10 to make sure
Nusselt number calc is not the reason.
[1]
https://yade-dem.org/doc/yade.wrapper.html#yade.wrapper.ThermalEngine.fluidConductionAreaFactor
[2] https://yade-dem.org/doc/yade.wrapper.html#yade.wrapper.ThermalEngine.fluidK
[3]
https://yade-dem.org/doc/yade.wrapper.html#yade.wrapper.ThermalEngine.minimumFluidCondDist
[4]
https://yade-dem.org/doc/yade.wrapper.html#yade.wrapper.ThermalEngine.ignoreFictiousConduction
[5]
https://yade-dem.org/doc/yade.wrapper.html#yade.wrapper.ThermalEngine.useHertzMethod
[6]
https://yade-dem.org/doc/yade.wrapper.html#yade.wrapper.ThermalEngine.uniformReynolds
Here is the MWE for reference:
###
from yade import pack, ymport, plot, utils, export, timing
import numpy as np
young=5e6
mn,mx=Vector3(0,0,0),Vector3(0.05,0.05,0.05)
identifier = '-thm_coupling'
if not os.path.exists('VTK'+identifier):
os.mkdir('VTK'+identifier)
if not os.path.exists('txt'+identifier):
os.mkdir('txt'+identifier)
O.materials.append(FrictMat(young=young*100,poisson=0.5,frictionAngle=0,density=2600,label='walls'))
O.materials.append(FrictMat(young=young,poisson=0.5,frictionAngle=radians(30),density=2600,label='spheres'))
walls=aabbWalls([mn,mx],thickness=0,material='walls')
wallIds=O.bodies.append(walls)
sp=pack.SpherePack()
sp.makeCloud(mn,mx,rMean=0.0015,rRelFuzz=.333,num=200,seed=11)
sp.toSimulation(color=(0.752, 0.752, 0.752),material='spheres')
triax=TriaxialStressController(
maxMultiplier=0,
finalMaxMultiplier=0,
thickness = 0,
goal1=-5,
goal2=-5,
goal3=-5,
max_vel=0.1,
stressMask = 7,
internalCompaction=False,
dead=True,
)
O.engines=[
ForceResetter(),
InsertionSortCollider([Bo1_Sphere_Aabb(aabbEnlargeFactor=1,label='is2aabb'),Bo1_Box_Aabb()]),
InteractionLoop(
[Ig2_Sphere_Sphere_ScGeom(interactionDetectionFactor=1,label='ss2sc'),Ig2_Box_Sphere_ScGeom()],
[Ip2_FrictMat_FrictMat_FrictPhys()],
[Law2_ScGeom_FrictPhys_CundallStrack()],label="iloop"
),
#GlobalStiffnessTimeStepper(active=1,timeStepUpdateInterval=100,timestepSafetyCoefficient=0.5),
triax,
FlowEngine(dead=1,label="flow",multithread=False),
ThermalEngine(dead=1,label='thermal'),
VTKRecorder(iterPeriod=500,fileName='VTK'+identifier+'/spheres-',recorders=['spheres','thermal','intr'],dead=0,label='VTKrec'),
NewtonIntegrator(damping=0.5)
]
O.dt=PWaveTimeStep()*.3
O.step()
ss2sc.interactionDetectionFactor=-1
is2aabb.aabbEnlargeFactor=-1
triax.dead=False
while 1:
O.run(1000, True)
unb=unbalancedForce()
print('unbalanced force:',unb,' mean stress: ',triax.meanStress)
if unb<0.1 and abs(triax.goal3-triax.meanStress)/-triax.goal3<0.001:
break
for b in O.bodies:
if isinstance(b.shape,Sphere):
b.dynamic=False
triax.internalCompaction=False
maxY=max([b.state.pos[1]+b.shape.radius for b in O.bodies if
isinstance(b.shape,Sphere)])
maxX=max([b.state.pos[0]+b.shape.radius for b in O.bodies if
isinstance(b.shape,Sphere)])
minY=min([b.state.pos[1]-b.shape.radius for b in O.bodies if
isinstance(b.shape,Sphere)])
minX=min([b.state.pos[0]-b.shape.radius for b in O.bodies if
isinstance(b.shape,Sphere)])
minZ=min([b.state.pos[2]-b.shape.radius for b in O.bodies if
isinstance(b.shape,Sphere)])
maxZ=max([b.state.pos[2]+b.shape.radius for b in O.bodies if
isinstance(b.shape,Sphere)])
print(minX,minY,minZ)
print(maxX,maxY,maxZ)
dz=maxZ-minZ
dy=maxY-minY
dx=maxX-minX
flow.debug=False
# add flow
flow.permeabilityMap = False
flow.pZero = 10
flow.meshUpdateInterval=-1
flow.fluidBulkModulus=2.2e9
flow.useSolver=4
flow.permeabilityFactor=-1e-5
Re: [Yade-users] [Question #698970]: Questions about updating yade
Question #698970 on Yade changed: https://answers.launchpad.net/yade/+question/698970 Status: Open => Answered Robert Caulk proposed the following answer: Sorry to hear that, I am not sure I can do much more to help you. One last thing could be to ensure you are downloading from "Main server" inside the "Software & Updates" menu. Although, I don't think that has anything to do with custom repos like yadedaily. Only stable yade. >>(I have tried the source-code,but because I am not familiar with compilation, there have been more problems..) Only one way to *become* familiar, right? Feel free to post questions with your compilation errors. -- You received this question notification because your team yade-users is an answer contact for Yade. ___ Mailing list: https://launchpad.net/~yade-users Post to : yade-users@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp
