yade-users team mailing list archive

[Jiannan Wang <question692123@xxxxxxxxxxxxxxxxxxxxx>]

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

    Status: Needs information => Open

Jiannan Wang gave more information on the question:
Hello Robert,

My mistake, I thought you meant just the packing MWE. Really appreciate
your patience. I attached the MWE in the follow. It is modified from
[1]. In order to make it not inconveniently long, I removed the saving
and plotting parts. The packing file is same as #3.

##########
from yade import pack, ymport
from yade import timing

young=5.7e+10
poisson=0.235
density=2630
tensileStrength=1.2e+8
cohesion=3.5e+8
frictionAngle=radians(50)

O.materials.append(JCFpmMat(young=young,poisson=poisson,frictionAngle=frictionAngle,density=density,tensileStrength=tensileStrength,cohesion=cohesion,label='spheres'))
O.materials.append(JCFpmMat(young=young,poisson=poisson,frictionAngle=0,density=density,tensileStrength=tensileStrength,cohesion=cohesion,label='walls'))

walls=aabbWalls([(0,0,0),(0.05,0.05,0.05)],thickness=0,material='walls')
wallIds=O.bodies.append(walls)
sp = O.bodies.append(ymport.textExt('compressedPack.spheres', 'x_y_z_r',material='spheres'))

O.engines=[
	ForceResetter(),
	InsertionSortCollider([Bo1_Sphere_Aabb(aabbEnlargeFactor=1.5,label='is2aabb'),Bo1_Box_Aabb()]),
	InteractionLoop(
		[Ig2_Sphere_Sphere_ScGeom(interactionDetectionFactor=1.5,label='ss2sc'),Ig2_Box_Sphere_ScGeom()],
		[Ip2_JCFpmMat_JCFpmMat_JCFpmPhys(cohesiveTresholdIteration=1,label='interactionPhys')],
		[Law2_ScGeom_JCFpmPhys_JointedCohesiveFrictionalPM(smoothJoint=True,label='interactionLaw'),
		Law2_ScGeom_FrictPhys_CundallStrack()],label="iloop"
	),
	GlobalStiffnessTimeStepper(active=1,timeStepUpdateInterval=100,timestepSafetyCoefficient=0.5),
	FlowEngine(dead=1,label="flow",multithread=False),
	ThermalEngine(dead=1,label='thermal'),
	NewtonIntegrator(damping=0.5)
]
O.step()
ss2sc.interactionDetectionFactor=-1
is2aabb.aabbEnlargeFactor=-1

for b in O.bodies:
	if isinstance(b.shape, Sphere):
		b.dynamic=False

flow.dead=0
flow.defTolerance=-0.003
flow.meshUpdateInterval=-100
flow.boundaryUseMaxMin=[1,1,1,1,1,1]
flow.breakControlledRemesh=True
flow.useSolver=4
flow.permeabilityFactor= 1
flow.viscosity= 0.001
flow.bndCondIsPressure=[0,0,0,0,0,0]
flow.bndCondValue=[0,0,0,0,0,0]
flow.thermalEngine=True
flow.debug=False
flow.fluidRho = 997
flow.fluidCp = 4181.7
flow.bndCondIsTemperature=[1,0,0,0,0,0]
flow.thermalEngine=True
flow.thermalBndCondValue=[343.15,0,0,0,0,0]
flow.tZero=333.15
flow.pZero=0

thermal.dead=0
thermal.debug=False
thermal.fluidConduction=True
thermal.ignoreFictiousConduction=True
thermal.conduction=True
thermal.thermoMech=True
thermal.solidThermoMech = True
thermal.fluidThermoMech = True
thermal.advection=True
thermal.bndCondIsTemperature=[1,0,0,0,0,0]
thermal.thermalBndCondValue=[343.15,0,0,0,0,0]
thermal.fluidK = 0.6069
thermal.fluidConductionAreaFactor=1.
thermal.uniformReynolds=10
thermal.particleT0 = 333.15
thermal.particleDensity=2600.
thermal.particleK = 2.
thermal.particleCp = 710.
thermal.tsSafetyFactor=0
thermal.useKernMethod=True
thermal.useHertzMethod=False

timing.reset()
O.dt=1e-12
O.dynDt=False
O.run(1,1)
#################################
Thank yo again.

Best regards
Jiannan

[1]: https://gitlab.com/yade-
dev/trunk/-/blob/master/examples/ThermalEngine/noFlowScenario.py

-- 
You received this question notification because your team yade-users is
an answer contact for Yade.