###################Cuerpo del codigo#########################

##########Establecer parametros generales###########

####################Variables del sistema################

########Variables especificas para la simulacion#################

############Definicion de la caja de simulacion###############

###############Definicion del modelo de contacto###############

##############Definicion de las propiedades del material##################

############Carga de la geometria desde .stl Y aplicar modelo de contacto a geometria#################

#################Generacion e insercion de particulas##########################






###############################################################Establecer parametros generales######################################################################

atom_style granular                                           #estilo de particula puede ser sphere o combinaciones 
boundary f f f                                                    #limites no periodicos y arreglados, particulas que salen no pueden volver a entrar 
units si                                                              #unidades sistema internacional
communicate	single vel yes	                        #default
newton off                                                        #tercera ley de newton apagada para relacion entre procesadores
dimension 3                                                      #problema en 3 dimensiones
atom_modify map array                                    #realizar busqueda en modo de matriz





####################################################################Variables del sistema###########################################################################

#Definicion de limites
variable xmin equal -8.0                 #cubo de 10m de lado
variable xmax equal 8.0

variable ymin equal -8.0
variable ymax equal 8.0

variable zmin equal -8.0
variable zmax equal 8.0

#Definicion del paso del tiempo

variable dt equal 1e-4			#paso de tiempo = 0.0001segundos; Cada paso de iteracion 0.0001 segundos








print 1


##########################################################Variables especificas para la simulacion#####################################################################

variable natoms equal 2 	#1:particulas #2:tolva            #cantidad de materiales en la simulacion




####Propiedades mecanicas de materiales####


####Modulos de Young####
variable	youngmodulo1	equal 	2.2e8	                           #modulo de young1/2: particulas / tolva [N/mm²]
variable	youngmodulo2	equal 	2.2e7		                           

####Coeficiente de Poisson####
variable 	poission1	equal	0.35                                      #coeficiente de possion1/2: particulas / tolva
variable 	poission2	equal	0.3


####Propiedades de contacto####


####Coeficiente de restitucion####
variable	cr11	equal	0.186                                           #coefifiente de restitucion entre particulas, particula/tolva, tolva/particula, tolva(entre geometrias)
variable	cr12	equal	0.199
variable	cr21	equal	0.199
variable	cr22	equal	0.186


print 2


####Coeficiente de friccion####
variable	cf11	equal	0.3                                                #coefifiente de friccion entre particulas, particula/tolva, tolva/particula, tolva(entre geometrias)  
variable	cf12	equal	0.2                                                                 
variable	cf21	equal	0.2
variable	cf22	equal	0.1




####Coeficiente k modelo finnie#####                                           #coeficiente k del modelo erosivo de finnie
variable	ck11	equal	1.0                                                                        #BUSCAR VALORES PARA TODOS LOS PARAMETROS                             
variable	ck12	equal	1.0                                              
variable	ck21	equal	1.0
variable        ck22	equal	1.0

####rolling friction coefficient####
variable	rf11	equal	0.21
variable	rf12	equal	0.21
variable	rf21	equal	0.21
variable	rf22	equal	0.21

print 3

####Variables para la particula####

variable	nradio	equal	4                                         #numero de particulas diferentes a agregar

variable	radio1   equal	0.25                                  #medidas de los radios en  [m]                 
variable	radio2	equal	0.27
variable	radio3	equal	0.3



variable	radio4	equal	0.33                                       

variable	frac1     equal	0.25		                                      #Distribucion granulometrica por porcentajes    #%                                    
variable	frac2    	equal	0.25			
variable	frac3    	equal	0.25		
variable	frac4    	equal	0.25		                                                                                                       #%

variable	densidad	equal	2500                                  #densidad de las particulas en [kg/m³]                     




####Parametros de llenado####

variable cara equal 9

variable	tiempollenado	equal	27	                                                              #tiempo de llenado de tolva con material (particulas) en [segundos]                 

variable	tiempollenadoporcara	equal	${tiempollenado}/${cara}

variable	masallenado	    equal	45000                                                             #cantidad de masa que soporta la tolva en [kg]

variable	flujollenado	    equal	${masallenado}/${tiempollenado}	                                 #flujo masico del llenado de la tolva con el material en [kg/s]

variable	pasosllenado  	equal	${tiempollenadoporcara}/${dt}	                                 #Transformar el tiempo de llenado en cantidad de pasos de llenado

variable    masaporcara    equal ${masallenado}/${cara}

variable	flujollenadoporcara	    equal	${masaporcara}/${tiempollenadoporcara}


print 4

####Parametros de asentamiento####

variable	tiempoespera	    equal	3                                                                  #tiempo de espera entre carga y descarga

variable	pasosespera	        equal	${tiempoespera}/${dt}		                                 #con el paso de tiempo definido se define l cantidad de pasos de asentamiento



####levantamiento de tolva####

variable vellevantar            equal	0.4	                                                                            #velocidad de levantamiento de tolva 

variable	tiempolevantar	    equal	15                                                                            #tiempo de levantamiento de tolva

variable	 pasoslevantar   	equal	${tiempolevantar}/${dt}		                                              #pasos de levantamiento de tolva



print 5



####################################################Definicion de la caja de simulacion################################################################################

region		domain block ${xmin} ${xmax} ${ymin} ${ymax} ${zmin} ${zmax} units box                          #definir limites en el plano        

create_box	2 domain                                                                                                                      #crear caja de simulacion con 2 tipos de materiales y con las dimensiones definidas anteriormente                                 
 
neighbor      0.008 bin
neigh_modify  delay 20




print 6

##################################################Definicion del modelo de contacto#################################################################################

pair_style gran model hertz tangential history rolling_friction cdt                                                                                           #modelo de contacto POR DEFINIR, SI SE DIVIDE LA SIMULACION O NO

pair_coeff	* *				                                                                                                                          #default

timestep ${dt}
fix    integr all nve/sphere

fix		gravi all gravity 9.81 vector 0.0 -1.0 0.0		                                                                                      #gravedad  de 9.81 m/s² en sentido negativo de y




print 7



########################################################Definicion de las propiedades del material######################################################################


fix 		m1 all property/global youngsModulus peratomtype ${youngmodulo1} ${youngmodulo2}                #agregar young y poisson uno por uno     

print 71    

fix 		m2 all property/global poissonsRatio peratomtype ${poission1} ${poission2}                                  #para coeficientes agregar numero de materiales primero y despues valores de coefcientes

print 72
fix 		m3 all property/global coefficientRestitution peratomtypepair ${natoms} ${cr11} ${cr12} ${cr21} ${cr22}   #ID + all property/global + variable + peratomtype/peratomtypepair+ valores
print 73
fix 		m4 all property/global coefficientFriction peratomtypepair ${natoms} ${cf11} ${cf12} ${cf21} ${cf22}
print 74
fix 		m5 all property/global k_finnie peratomtypepair ${natoms} ${ck11} ${ck12} ${ck21} ${ck22}
print 75
fix 		m6 all property/global coefficientRollingFriction peratomtypepair ${natoms} ${rf11} ${rf12} ${rf21} ${rf22}



print 8

###################################Carga de la geometria desde .stl Y aplicar modelo de contacto a geometria#######################################

fix			cad1 	all mesh/surface/stress file tolvajimbo.stl 	type 1  scale 0.001 	wear finnie 	#cargar malla desde archivo stl Tipo 1 para geometria. escala para pasar de mm a m

fix         insercion1 all mesh/surface/planar file cara1.stl type 1 scale 0.001            
fix         insercion2 all mesh/surface/planar file cara2.stl type 1 scale 0.001       
fix         insercion3 all mesh/surface/planar file cara3.stl type 1 scale 0.001       
fix         insercion4 all mesh/surface/planar file cara4.stl type 1 scale 0.001         ####agregar el modelo de contacto a todas las geometrias, es necesario escribir todas las ID en las que se aplica el modelo de contacto
fix         insercion5 all mesh/surface/planar file cara5.stl type 1 scale 0.001       
fix         insercion6 all mesh/surface/planar file cara6.stl type 1 scale 0.001       
fix         insercion7 all mesh/surface/planar file cara7.stl type 1 scale 0.001       
fix         insercion8 all mesh/surface/planar file cara8.stl type 1 scale 0.001       
fix         insercion9 all mesh/surface/planar file cara9.stl type 1 scale 0.001       


fix 		paredes 	all wall/gran model hertz tangential history rolling_friction cdt mesh n_meshes 1 meshes cad1 



print 9




##################################Generacion e insercion de particulas esfericas#########################################

fix		pts1 all particletemplate/sphere 10487 atom_type 1 density constant ${densidad} radius constant ${radio1}	#agregar particulas por medida de radio y despues su porcentaje en pdd1

fix		pts2 all particletemplate/sphere 11887 atom_type 1 density constant ${densidad} radius constant ${radio2}	#pts1,2,3.... + all particletemplate/sphere + numero primo mayor a 10000+ atom_type + densidad + radio (para esferas)

fix		pts3 all particletemplate/sphere 10729 atom_type 1 density constant ${densidad} radius constant ${radio3}

fix		pts4 all particletemplate/sphere 11351 atom_type 1 density constant ${densidad} radius constant ${radio4}

fix 	    pdd1 all particledistribution/discrete 32452867 ${nradio} pts1 ${frac1} pts2 ${frac2} pts3 ${frac3} pts4 ${frac4}                                     #pdd1 + all particledistribution/discrete  + numero primo mayor a 10000 + cantidad de particulas distintas + id particulas1 porcentaje+ id particulas2 porcentaje....





print 10


##################################################region e insercion#################################################################################

run 1
dump   dmp1 all custom/vtk  1000 postesferasgrandes1/particulasesferas_*.vtk id type type x y z ix iy iz vx vy vz fx fy fz omegax omegay omegaz radius 
dump   dumpstress all mesh/gran/VTK 1000 postesferasgrandes1/tolvaesferas_*.vtk stress wear cad1 
compute       1 all erotate/sphere
thermo_style  custom step atoms ke c_1 vol time
thermo        1000
thermo_modify lost ignore norm no


fix    ins1  all insert/stream seed 49979687 distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion3
run	 ${pasosllenado}

fix    ins2  all insert/stream seed 38629 distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion1
run	 ${pasosllenado}

fix    ins3  all insert/stream seed 37619  distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion2
run	 ${pasosllenado}

fix    ins4  all insert/stream seed 36493 distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion4
run	 ${pasosllenado}

fix    ins5  all insert/stream seed 49331 distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion8
run	 ${pasosllenado}

fix    ins6  all insert/stream seed 49747 distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion9
run	 ${pasosllenado}

fix    ins7  all insert/stream seed 49939 distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion5
run	 ${pasosllenado}

fix    ins8  all insert/stream seed 49109 distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion7
run	 ${pasosllenado}

fix    ins9  all insert/stream seed 48311 distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion6
run	 ${pasosllenado}



print 11




######################################################comandos de salida: desgaste y energia termica#######################################################################


print 12





##################################correr tiempo de espera#############################################################
run	 ${pasosespera}
print 15



#################################mover tolva y correr#############################################################
fix		levantar	all move/mesh mesh cad1 	rotate origin -2.4. -0.6. 0. axis 0. 0. 1. period 104
print 16
run	 ${pasoslevantar}
unfix levantar 

run 15000
print 17

fix bajar all move/mesh mesh cad1 	rotate origin -2.4. -0.6. 0. axis 0. 0. -1. period 6.92307692308

run 10000
unfix bajar
fix    ins10  all insert/stream seed 103889 distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion3
run	 ${pasosllenado}

fix    ins11  all insert/stream seed 103903 distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion1
run	 ${pasosllenado}

fix    ins12  all insert/stream seed 103913  distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion2
run	 ${pasosllenado}

fix    ins13  all insert/stream seed 103919   distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion4
run	 ${pasosllenado}

fix    ins14  all insert/stream seed 103951  distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion8
run	 ${pasosllenado}

fix    ins15  all insert/stream seed 103963 distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion9
run	 ${pasosllenado}

fix    ins16  all insert/stream seed 103967 distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion5
run	 ${pasosllenado}

fix    ins17  all insert/stream seed 103969  distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion7
run	 ${pasosllenado}

fix    ins18  all insert/stream seed 103979  distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion6
run	 ${pasosllenado}


print 11




######################################################comandos de salida: desgaste y energia termica#######################################################################


print 12





##################################correr tiempo de espera#############################################################
run	 ${pasosespera}
print 15



#################################mover tolva y correr#############################################################
fix		levantar2	all move/mesh mesh cad1 	rotate origin -2.4. -0.6. 0. axis 0. 0. 1. period 104
print 16
run	 ${pasoslevantar}
unfix levantar2 

run 15000


fix bajar2 all move/mesh mesh cad1 	rotate origin -2.4. -0.6. 0. axis 0. 0. -1. period 6.92307692308

run 10000
unfix bajar2

fix    ins19  all insert/stream seed 103981 distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion3
run	 ${pasosllenado}

fix    ins20  all insert/stream seed 103991  distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion1
run	 ${pasosllenado}

fix    ins21  all insert/stream seed 103993  distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion2
run	 ${pasosllenado}

fix    ins22  all insert/stream seed 103997 distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion4
run	 ${pasosllenado}

fix    ins23  all insert/stream seed 104003   distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion8
run	 ${pasosllenado}

fix    ins24  all insert/stream seed 104009  distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion9
run	 ${pasosllenado}

fix    ins25  all insert/stream seed 104021  distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion5
run	 ${pasosllenado}

fix    ins26  all insert/stream seed 104033  distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion7
run	 ${pasosllenado}

fix    ins27  all insert/stream seed 104047 distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion6
run	 ${pasosllenado}


print 11




######################################################comandos de salida: desgaste y energia termica#######################################################################


print 12





##################################correr tiempo de espera#############################################################
run	 ${pasosespera}
print 15



#################################mover tolva y correr#############################################################
fix		levantar3	all move/mesh mesh cad1 	rotate origin -2.4. -0.6. 0. axis 0. 0. 1. period 104
print 16
run	 ${pasoslevantar}
unfix levantar3 

run 15000


fix bajar3 all move/mesh mesh cad1 	rotate origin -2.4. -0.6. 0. axis 0. 0. -1. period 6.92307692308

run 10000

unfix bajar3



fix    ins28  all insert/stream seed 104053  distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion3
run	 ${pasosllenado}

fix    ins29  all insert/stream seed 104059 distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion1
run	 ${pasosllenado}

fix    ins30 all insert/stream seed 104087 distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion2
run	 ${pasosllenado}

fix    ins31  all insert/stream seed 104089 distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion4
run	 ${pasosllenado}

fix    ins32  all insert/stream seed 104107 distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion8
run	 ${pasosllenado}

fix    ins33  all insert/stream seed 104113 distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion9
run	 ${pasosllenado}

fix    ins34  all insert/stream seed 104119 distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion5
run	 ${pasosllenado}

fix    ins35  all insert/stream seed 104123 distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion7
run	 ${pasosllenado}

fix    ins36  all insert/stream seed 104147 distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion6
run	 ${pasosllenado}


print 11




######################################################comandos de salida: desgaste y energia termica#######################################################################


print 12





##################################correr tiempo de espera#############################################################
run	 ${pasosespera}
print 15



#################################mover tolva y correr#############################################################
fix		levantar4	all move/mesh mesh cad1 	rotate origin -2.4. -0.6. 0. axis 0. 0. 1. period 104
print 16
run	 ${pasoslevantar}
unfix levantar4 

run 15000


fix bajar4 all move/mesh mesh cad1 	rotate origin -2.4. -0.6. 0. axis 0. 0. -1. period 6.92307692308

run 10000
unfix bajar4

fix    ins37  all insert/stream seed 104149 distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion3
run	 ${pasosllenado}

fix    ins38  all insert/stream seed 104161 distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion1
run	 ${pasosllenado}

fix    ins39 all insert/stream seed 104173 distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion2
run	 ${pasosllenado}

fix    ins40  all insert/stream seed 104179 distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion4
run	 ${pasosllenado}

fix    ins41  all insert/stream seed 104183  distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion8
run	 ${pasosllenado}

fix    ins42  all insert/stream seed 104207  distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion9
run	 ${pasosllenado}

fix    ins43  all insert/stream seed 104231 distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion5
run	 ${pasosllenado}

fix    ins44  all insert/stream seed 104233 distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion7
run	 ${pasosllenado}

fix    ins45  all insert/stream seed 104239 distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion6
run	 ${pasosllenado}


print 11




######################################################comandos de salida: desgaste y energia termica#######################################################################


print 12





##################################correr tiempo de espera#############################################################
run	 ${pasosespera}
print 15



#################################mover tolva y correr#############################################################
fix		levantar5	all move/mesh mesh cad1 	rotate origin -2.4. -0.6. 0. axis 0. 0. 1. period 104
print 16
run	 ${pasoslevantar}
unfix levantar5 

run 15000


fix bajar5 all move/mesh mesh cad1 	rotate origin -2.4. -0.6. 0. axis 0. 0. -1. period 6.92307692308

run 10000
unfix bajar5
fix    ins46  all insert/stream seed 104243 distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion3
run	 ${pasosllenado}

fix    ins47  all insert/stream seed 104281 distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion1
run	 ${pasosllenado}

fix    ins48 all insert/stream seed 104287 distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion2
run	 ${pasosllenado}

fix    ins49  all insert/stream seed 104297 distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion4
run	 ${pasosllenado}

fix    ins50  all insert/stream seed 104309   distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion8
run	 ${pasosllenado}

fix    ins51  all insert/stream seed 104311 distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion9
run	 ${pasosllenado}

fix    ins52  all insert/stream seed 104323 distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion5
run	 ${pasosllenado}

fix    ins53  all insert/stream seed 104327 distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion7
run	 ${pasosllenado}

fix    ins54  all insert/stream seed 104369  distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion6
run	 ${pasosllenado}


print 11




######################################################comandos de salida: desgaste y energia termica#######################################################################


print 12





##################################correr tiempo de espera#############################################################
run	 ${pasosespera}
print 15



#################################mover tolva y correr#############################################################
fix		levantar6	all move/mesh mesh cad1 	rotate origin -2.4. -0.6. 0. axis 0. 0. 1. period 104
print 16
run	 ${pasoslevantar}
unfix levantar6 

run 15000


fix bajar6 all move/mesh mesh cad1 	rotate origin -2.4. -0.6. 0. axis 0. 0. -1. period 6.92307692308

run 10000
unfix bajar6
fix    ins55  all insert/stream seed 104381 distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion3
run	 ${pasosllenado}

fix    ins56  all insert/stream seed 104383   distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion1
run	 ${pasosllenado}

fix    ins57 all insert/stream seed 104393  distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion2
run	 ${pasosllenado}

fix    ins58  all insert/stream seed 104399 distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion4
run	 ${pasosllenado}

fix    ins59  all insert/stream seed 104417 distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion8
run	 ${pasosllenado}

fix    ins60  all insert/stream seed 104459 distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion9
run	 ${pasosllenado}

fix    ins61  all insert/stream seed 104471 distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion5
run	 ${pasosllenado}

fix    ins62  all insert/stream seed 104473 distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion7
run	 ${pasosllenado}

fix    ins63  all insert/stream seed 104479 distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion6
run	 ${pasosllenado}


print 11




######################################################comandos de salida: desgaste y energia termica#######################################################################


print 12





##################################correr tiempo de espera#############################################################
run	 ${pasosespera}
print 15



#################################mover tolva y correr#############################################################
fix		levantar7	all move/mesh mesh cad1 	rotate origin -2.4. -0.6. 0. axis 0. 0. 1. period 104
print 16
run	 ${pasoslevantar}
unfix levantar7 

run 15000


fix bajar7 all move/mesh mesh cad1 	rotate origin -2.4. -0.6. 0. axis 0. 0. -1. period 6.92307692308

run 10000
unfix bajar7
fix    ins64  all insert/stream seed 104491 distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion3
run	 ${pasosllenado}

fix    ins65  all insert/stream seed 104513 distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion1
run	 ${pasosllenado}

fix    ins66 all insert/stream seed 104527 distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion2
run	 ${pasosllenado}

fix    ins67  all insert/stream seed 104537 distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion4
run	 ${pasosllenado}

fix    ins68  all insert/stream seed 104543 distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion8
run	 ${pasosllenado}

fix    ins69  all insert/stream seed 104549 distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion9
run	 ${pasosllenado}

fix    ins70  all insert/stream seed 104551 distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion5
run	 ${pasosllenado}

fix    ins71  all insert/stream seed 104561 distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion7
run	 ${pasosllenado}

fix    ins72  all insert/stream seed 104579 distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion6
run	 ${pasosllenado}


print 11




######################################################comandos de salida: desgaste y energia termica#######################################################################


print 12





##################################correr tiempo de espera#############################################################
run	 ${pasosespera}
print 15



#################################mover tolva y correr#############################################################
fix		levantar8	all move/mesh mesh cad1 	rotate origin -2.4. -0.6. 0. axis 0. 0. 1. period 104
print 16
run	 ${pasoslevantar}
unfix levantar8 

run 15000


fix bajar8 all move/mesh mesh cad1 	rotate origin -2.4. -0.6. 0. axis 0. 0. -1. period 6.92307692308

run 10000
unfix bajar8
fix    ins73  all insert/stream seed 104593  distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion3
run	 ${pasosllenado}

fix    ins74  all insert/stream seed 104597 distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion1
run	 ${pasosllenado}

fix    ins75 all insert/stream seed 104623  distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion2
run	 ${pasosllenado}

fix    ins76  all insert/stream seed 104639  distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion4
run	 ${pasosllenado}

fix    ins77  all insert/stream seed 104651  distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion8
run	 ${pasosllenado}

fix    ins78  all insert/stream seed 104659 distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion9
run	 ${pasosllenado}

fix    ins79  all insert/stream seed 104677 distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion5
run	 ${pasosllenado}

fix    ins80  all insert/stream seed  104681  distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion7
run	 ${pasosllenado}

fix    ins81  all insert/stream seed 104683   distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion6
run	 ${pasosllenado}


print 11




######################################################comandos de salida: desgaste y energia termica#######################################################################


print 12





##################################correr tiempo de espera#############################################################
run	 ${pasosespera}
print 15



#################################mover tolva y correr#############################################################
fix		levantar9	all move/mesh mesh cad1 	rotate origin -2.4. -0.6. 0. axis 0. 0. 1. period 104
print 16
run	 ${pasoslevantar}
unfix levantar9

run 15000


fix bajar9 all move/mesh mesh cad1 	rotate origin -2.4. -0.6. 0. axis 0. 0. -1. period 6.92307692308

run 10000
unfix bajar9
fix    ins82  all insert/stream seed 104693 distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion3
run	 ${pasosllenado}

fix    ins83 all insert/stream seed 104701 distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion1
run	 ${pasosllenado}

fix    ins84 all insert/stream seed 104707 distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion2
run	 ${pasosllenado}

fix    ins85  all insert/stream seed 104711  distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion4
run	 ${pasosllenado}

fix    ins86  all insert/stream seed 104717 distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion8
run	 ${pasosllenado}

fix    ins87  all insert/stream seed 104723  distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion9
run	 ${pasosllenado}

fix    ins88  all insert/stream seed 104729   distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion5
run	 ${pasosllenado}

fix    ins89  all insert/stream seed 100333  distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion7
run	 ${pasosllenado}

fix    ins90  all insert/stream seed 100343   distributiontemplate pdd1 mass ${masaporcara} massrate ${flujollenadoporcara} insert_every 10000 overlapcheck yes all_in no vel constant 0. -1. 0.  insertion_face insercion6
run	 ${pasosllenado}


print 11




######################################################comandos de salida: desgaste y energia termica#######################################################################


print 12





##################################correr tiempo de espera#############################################################
run	 ${pasosespera}
print 15



#################################mover tolva y correr#############################################################
fix		levantar10	all move/mesh mesh cad1 	rotate origin -2.4. -0.6. 0. axis 0. 0. 1. period 104
print 16
run	 ${pasoslevantar}
unfix levantar10 

run 15000


fix bajar10 all move/mesh mesh cad1 	rotate origin -2.4. -0.6. 0. axis 0. 0. -1. period 6.92307692308

run 10000
unfix bajar10