nitsche_test.jou

Go to the documentation of this file.

#!python
autofactor=5; 
cubit.cmd('set duplicate block elements on')    
 
#============================================================= 
#Geometry
#============================================================= 
 
cubit.cmd('reset')
 
cubit.cmd('brick x 1 y 1 z 1')
cubit.cmd('create sphere radius 0.3')
#cubit.cmd('move Volume 2 x 0.14 y 0.14 z 0.14')
cubit.cmd('subtract volume 2 from volume 1 ')
cubit.cmd('delete volume 2')
 
#============================================================= 
#Defining blocks for elastic, transversely-isotropic and potential flow problems  
#============================================================= 
 
vol=['all']  
mat=['MAT_ELASTIC_1']  
for i in range(0, 1):  
    str1='block  ' + str(i+1) +' volume '+vol[i]; cubit.cmd(str1)  
    str1='block  ' + str(i+1) +' name "'+mat[i] + '"'; cubit.cmd(str1)  
 
#============================================================= 
#Surface elements
#============================================================= 
cubit.cmd('sideset 103 surface 1 4 3 2 6 5')  # all boundary surfaces  
cubit.cmd('sideset 102 surface 1 5 6 ')  # all -ve boundary surfaces
cubit.cmd('sideset 101 surface 2 4 3  ')  # all +ve boundary surfaces
 
#============================================================= 
#Material properties for matrix part  
#============================================================= 
  
cubit.cmd('block 1 attribute count 2')  
Em=3.5e3; Enu=0.3;  #giga to mega as we used dimension in mm  
 
Elastic=[str(Em), str(Enu)]  
for i in range(0, 2):  
    str1='block 1 attribute index ' + str(i+1) +' '+Elastic[i]; cubit.cmd(str1)  
  
#============================================================= 
#Defining surfaces for dispacement, traction and periodic boundary conditions  
#============================================================= 
 
cubit.cmd('surface 1 scheme trimesh') 
str1='surface 1  size auto factor '+str(autofactor); cubit.cmd(str1)
cubit.cmd('mesh surface 1') 
cubit.cmd('surface 4 scheme trimesh') 
str1='surface 4  size auto factor '+str(autofactor); cubit.cmd(str1)
cubit.cmd('mesh surface 4') 
cubit.cmd('surface 3 scheme trimesh') 
str1='surface 3  size auto factor '+str(autofactor); cubit.cmd(str1)
cubit.cmd('mesh surface 3') 
 
cubit.cmd('surface 2  scheme copy source surface 1  source vertex 2  target vertex 5  source curve 1  target curve 5  nosmoothing')
cubit.cmd('mesh surface 2')
cubit.cmd('surface 6  scheme copy source surface 4  source vertex 3  target vertex 2  source curve 3  target curve 1  nosmoothing')
cubit.cmd('mesh surface 6')
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')
 
#mesh
str1='volume all size auto factor '+str(autofactor); cubit.cmd(str1)
cubit.cmd('volume all scheme tetmesh')  
cubit.cmd('mesh volume all')  
 
#============================================================= 
#High order geomety approximation
#============================================================= 
cubit.cmd('block 5 tet all')         
cubit.cmd('block 5 element type tetra10')