New question #677437 on Yade: https://answers.launchpad.net/yade/+question/677437
Hi， I want to save the result to txt.I learned the code here,https://answers.launchpad.net/yade/+question/209832.But I found the data in txt.bz2 is just 5000 iter at the beginning.How do I get all the data？ Thanks， Feng Here's my code： ###################### # encoding: utf-8 from yade import pack, qt, plot import matplotlib; matplotlib.rc('axes',grid=True) from yade import pack import pylab import yade.timing; O.timingEnabled=True sigmaIso=-25000 O.periodic=True spheres=O.materials.append(FrictMat(young=64e9,poisson=0.12,density=2650,frictionAngle=0.0005)) s=O.materials.append(FrictMat(young=64e9,poisson=0.12,density=2650,frictionAngle=0.24)) psdSizes,psdCumm=[0.1,0.12,0.15,0.19,0.20,0.23,0.27,0.28,0.30,0.32,0.33,0.50],[0.0001,0.003,0.03,0.11,0.25,0.43,0.70,0.85,0.95,0.97,0.98,1] sp=pack.SpherePack() sp.makeCloud((0,0,0),(3,3,3),psdSizes=psdSizes,psdCumm=psdCumm,distributeMass=True,periodic=True) n=sp.toSimulation(material=spheres) O.engines=[ ForceResetter(), InsertionSortCollider([Bo1_Sphere_Aabb()]), InteractionLoop( [Ig2_Sphere_Sphere_ScGeom()], [Ip2_FrictMat_FrictMat_FrictPhys()], [Law2_ScGeom_FrictPhys_CundallStrack()] ), PeriTriaxController(label='triax', goal=(sigmaIso,sigmaIso,sigmaIso),stressMask=7, dynCell=True,maxStrainRate=(0.5,0.5,0.5), maxUnbalanced=1,relStressTol=0.05, doneHook='compactionFinished()' ), NewtonIntegrator(damping=.1), PyRunner(command='addPlotData1()',iterPeriod=200), ] O.dt=1*PWaveTimeStep() O.trackEnergy=True def addPlotData1(): plot.addData(s0=triax.stress[0],s1=triax.stress[1],s2=triax.stress[2],q=triax.stress[2]-0.5*(triax.stress[0]+triax.stress[1]), p=(triax.stress[0]+triax.stress[1]+triax.stress[2])/3) print('stress1',triax.stress[0],triax.stress[1],triax.stress[2]) print('strain',triax.strain[0],triax.strain[1],triax.strain[2]) print('q',triax.stress[2]-0.5*(triax.stress[0]+triax.stress[1]),'p',(triax.stress[0]+triax.stress[1]+triax.stress[2])/3) pylab.semilogx(*sp.psd(bins=30,mass=True),label='Mass PSD of (free) %d random spheres'%len(sp)) pylab.legend() pylab.show() O.run(5000,True) plot.plots={'q':('p',)} plot.plot() plot.saveDataTxt('pq.txt.bz2') qt.Controller() O.saveTmp() def compactionFinished(): O.cell.trsf=Matrix3.Identity triax.goal=(sigmaIso,sigmaIso,-800000) triax.stressMask=4 triax.maxStrainRate=(0,0,0.5) triax.maxUnbalanced=10 triax.relStressTol=0.05 triax.doneHook='triaxFinished()' def triaxFinished(): print 'Finished' print(yade.timing.stats()) O.pause() -- 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
