Question #689234 on Yade changed:
https://answers.launchpad.net/yade/+question/689234
Chien-Cheng Hung gave more information on the question:
Hi Jan,
I would like to export the function in vtk files as below. I adjust the body of
function unlike the above you showed me.
But it shows me that "name i is not defined". How do I modify my code?
Thank you for your help again!
###
def weakeningTemp(interaction):
heatCapacity = 900 # [J/kg*K]
density = 2500 # [kg/m3]
totalHeat = 0
thermal_diff = 1.e-6
shear_modulus = 22e9
tau_c = 0.1*shear_modulus
initialTemp = 20
for i in O.interactions():
if O.bodies[i.id1].isReal and O.bodies[i.id2].isReal:
penetrationDepth = i.geom.penetrationDepth
radius1 = i.geom.refR1
radius2 = i.geom.refR2
effectiveRadius = (i.geom.refR1 * i.geom.refR2 )/(
i.geom.refR1+i.geom.refR2 )
contactRadius = math.sqrt( i.geom.penetrationDepth *effectiveRadius)
Da = 2 * contactRadius
contactArea = np.pi*(contactRadius**2)
relativeVelocity = i.geom.shearInc.norm()/O.dt
weakeningTemp =
math.sqrt((relativeVelocity*Da)/(np.pi*thermal_diff))*tau_c/(density*heatCapacity)
+ initialTemp
return weakeningTemp
builtins.weakeningTemp = weakeningTemp # such that VTKExporter has
access to it
vtk = export.VTKExporter("weakingTemp")
def vtkExport():
vtk.exportSpheres(ids=list(range(2,len(O.bodies))),what=dict(weakeningTemp="weakeningTemp(i)"))
###
Cheers,
Chien-Cheng
--
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 : [email protected]
Unsubscribe : https://launchpad.net/~yade-users
More help : https://help.launchpad.net/ListHelp
