Cube_HOG.jou

Go to the documentation of this file.

#!python
#To create RVE for all the three type of boundary conditions, i.e. (Linear dispacemet, traction and periodic)
 
cubit.cmd('new')   
cubit.cmd('brick x 1 y 1 z 1')   
   
str1='block 1 ' +' volume 1'; cubit.cmd(str1)   
str1='block 1 ' +' name "MAT_ELASTIC"'; cubit.cmd(str1)   
cubit.cmd('block 1 attribute count 2')   
Elastic=['200', '0.3']   
for i in range(0, 2):   
    str1='block 1 attribute index ' + str(i+1) +' '+Elastic[i]; cubit.cmd(str1)   
   
interval=1;    
str1='surface 4 interval '+str(interval); cubit.cmd(str1)   
cubit.cmd('surface 4 scheme trimesh')   
cubit.cmd('mesh surface 4')   
cubit.cmd('surface 6 scheme copy source surface 4 source vertex 4 target vertex 1 source curve 9 target curve 10 nosmoothing')   
cubit.cmd('mesh surface 6')   
   
str1='surface 3 interval '+str(interval); cubit.cmd(str1)   
cubit.cmd('surface 3 scheme trimesh')   
cubit.cmd('mesh surface 3')   
cubit.cmd('surface 5 scheme copy source surface 3 source vertex 7 target vertex 8 source curve 9 target curve 11 nosmoothing')   
cubit.cmd('mesh surface 5')   
   
str1='surface 1 interval '+str(interval); cubit.cmd(str1)    
cubit.cmd('surface 1 scheme trimesh')    
cubit.cmd('mesh surface 1')    
cubit.cmd('surface 2 scheme copy source surface 1 source vertex 4 target vertex 7 source curve 3 target curve 7 nosmoothing')    
cubit.cmd('mesh surface 2')    
   
cubit.cmd('volume 1 scheme tetmesh')   
cubit.cmd('mesh volume 1')   
   
cubit.cmd('set duplicate block elements on')    
cubit.cmd('block 2 tet all')         
cubit.cmd('block 2 element type tetra10')      
   
   
cubit.cmd('sideset 101 surface 2 3 4')  # all -ve boundary surfaces   
cubit.cmd('sideset 101 name "NegBoundSurf"')
 
cubit.cmd('sideset 102 surface 1 5 6')  # all +ve boundary surfaces 
cubit.cmd('sideset 102 name "PosBoundSurf"')
  
cubit.cmd('sideset 103 surface 1 2 3 4 5 6')  # all boundary surfaces  
cubit.cmd('sideset 103 name "AllBoundSurf"')
 
  
cubit.cmd('save as "/Users/zahur/Documents/moFEM/mofem-cephas/mofem_v0.2/users_modules/homogenisation/meshes/Cube_HOG.cub" overwrite')  
 
 
 
 
#cubit.cmd('new')  
#cubit.cmd('brick x 1 y 1 z 1')  
#  
#cubit.cmd('create cylinder height 2 radius 0.35')  
#cubit.cmd('subtract volume 2 from volume 1')  
#  
#str1='block 1 ' +' volume 1'; cubit.cmd(str1)  
#str1='block 1 ' +' name "MAT_ELASTIC"'; cubit.cmd(str1)  
#cubit.cmd('block 1 attribute count 2')  
#Elastic=['200', '0.3']  
#for i in range(0, 2):  
#    str1='block 1 attribute index ' + str(i+1) +' '+Elastic[i]; cubit.cmd(str1)  
#  
#interval=4;   
#str1='surface 4 interval '+str(interval); cubit.cmd(str1)  
#cubit.cmd('surface 4 scheme trimesh')  
#cubit.cmd('mesh surface 4')  
#cubit.cmd('surface 6 scheme copy source surface 4 source vertex 4 target vertex 1 source curve 9 target curve 10 nosmoothing')  
#cubit.cmd('mesh surface 6')  
#  
#str1='surface 3 interval '+str(interval); cubit.cmd(str1)  
#cubit.cmd('surface 3 scheme trimesh')  
#cubit.cmd('mesh surface 3')  
#cubit.cmd('surface 5 scheme copy source surface 3 source vertex 7 target vertex 8 source curve 9 target curve 11 nosmoothing')  
#cubit.cmd('mesh surface 5')  
#  
#str1='surface 11 interval '+str(interval); cubit.cmd(str1)   
#cubit.cmd('surface 11 scheme trimesh')   
#cubit.cmd('mesh surface 11')   
#cubit.cmd('surface 12 scheme copy source surface 11 source vertex 4 target vertex 7 source curve 3 target curve 7 nosmoothing')   
#cubit.cmd('mesh surface 12')   
#  
#cubit.cmd('volume 1 scheme tetmesh')  
#cubit.cmd('mesh volume 1')  
#  
#cubit.cmd('set duplicate block elements on')   
#cubit.cmd('block 2 tet all')        
#cubit.cmd('block 2 element type tetra10')     
#  
#  
#cubit.cmd('sideset 101 surface 12 3 4')  # all -ve boundary surfaces  
#cubit.cmd('sideset 102 surface 11 5 6')  # all +ve boundary surfaces  
#cubit.cmd('sideset 103 surface 11 12 3 4 5 6')  # all boundary surfaces  
#  
#cubit.cmd('save as "/Users/zahur/Documents/moFEM/mofem-cephas/mofem_v0.2/users_modules/homogenisation/meshes/Cube_HOG.cub" overwrite')