yade-users team mailing list archive

[Yor1 <question276604@xxxxxxxxxxxxxxxxxxxxx>]

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

Yor1 posted a new comment:
Hello Jan,

I use the formula of plastic dissipation used in sources.
In fact, i integrate the formula in the recorder like this :

tensCks=shearCks=cks=cks0=0
e10=e20=e30=0

def recorder():
    E=0
    Ub=0
    Uc=0
    Wp=0
    Wpla=0

    global tensCks, shearCks, e10,e20,e30
    tensCks=0
    shearCks=0
    for o in O.bodies:
	tensCks+=o.state.tensBreak
	shearCks+=o.state.shearBreak

    for i in O.interactions:
       #Ft=i.phys.shearForce
       #du=i.geom.shearInc
       #kt=i.phys.ks
       #Fn=i.phys.normalForce
       #Ftt=Ft-kt*du
       #maxFs=(Fn.norm()*tan(math.pi*FRICT/180)

      if((i.phys.shearForce-i.phys.ks*i.geom.shearInc).norm() > i.phys.normalForce.norm()*tan(math.pi*FRICT/180) ):
                Ftt=i.phys.shearForce-i.phys.ks*i.geom.shearInc
	 	ratio = i.phys.normalForce.norm()*tan(math.pi*FRICT/180) / Ftt.norm()
   	 	trialForce=Ftt
   	 	Ftt *= ratio
   	 	e=((1/i.phys.ks)*(trialForce-Ftt)).dot(Ftt)
         	Wpla=Wpla+e
      if not i.isReal : continue
      if isinstance(O.bodies[i.id1].shape,Sphere) and isinstance(O.bodies[i.id2].shape,Sphere):
       	E+=0.5*(i.phys.normalForce.squaredNorm()/i.phys.kn + i.phys.shearForce.squaredNorm()/i.phys.ks)
		
    yade.plot.addData( t=O.time
			,i=O.iter
			,e1=triax.strain[0]-e10
			,e2=triax.strain[1]-e20
			,e3=triax.strain[2]-e30
			,s1=triax.stress(triax.wall_right_id)[0]
			,s2=triax.stress(triax.wall_top_id)[1]
			,s3=triax.stress(triax.wall_front_id)[2]
			,tc=0.5*tensCks,sc=0.5*shearCks,unbF=utils.unbalancedForce()
			,E=E   #Energie elastique
			,Wp=triax.externalWork #Travail des contraintes appliquees aux limites de echantillon
			,Wk=utils.kineticEnergy() #Energie cinetique
			,Wd=2*utils.kineticEnergy()*DAMP*0.1*utils.PWaveTimeStep() #Energie dissipee par amortissement
    			,Wpla=Wpla #energie plastique
    )
    plot.saveDataTxt(OUT)

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