Textile_RVE.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')
 
W_wrap=0.3; H_wrap=0.1514; hgap_wrap=0.09 
W_weft=0.3; H_weft=0.0757; hgap_weft=1.2 
L_RVE=3.0;    W_RVE=0.8;        H_RVE=0.3;  
Vgap=0.012 
 
s_weft=hgap_weft+W_weft;    s_wraf=hgap_wrap+W_wrap
T=(H_weft/2+H_wrap/2)/2+Vgap
 
L_RVE=2*s_weft;    W_RVE=2*s_wraf;     H_RVE=0.3;    
 
coordx=-22.5*s_weft;  coordy=T; coordz=0;
for i in range(0, 48):
    if i % 2 == 0:
        str1='create vertex ' + str(coordx) +' '+str(coordy)+' '+str(coordz)+' color' ; coordx=coordx+s_weft; 
        cubit.cmd(str1)
    if i % 2 == 1:
        str1='create vertex ' + str(coordx) +' '+str(-coordy)+' '+str(coordz)+' color' ; coordx=coordx+s_weft; 
        cubit.cmd(str1)
 
coordx=-22.5*s_weft;  coordy=-T; coordz=s_wraf; 
for i in range(0, 48): 
    if i % 2 == 0: 
        str1='create vertex ' + str(coordx) +' '+str(coordy)+' '+str(coordz)+' color' ; coordx=coordx+s_weft;  
        cubit.cmd(str1) 
    if i % 2 == 1: 
        str1='create vertex ' + str(coordx) +' '+str(-coordy)+' '+str(coordz)+' color' ; coordx=coordx+s_weft;  
        cubit.cmd(str1) 
 
coordx=-0.5*s_weft;  coordy=-T; coordz=-22*s_wraf; 
for i in range(0, 48): 
    if i % 2 == 0: 
        str1='create vertex ' + str(coordx) +' '+str(coordy)+' '+str(coordz)+' color' ; coordz=coordz+s_wraf;  
        cubit.cmd(str1) 
    if i % 2 == 1: 
        str1='create vertex ' + str(coordx) +' '+str(-coordy)+' '+str(coordz)+' color' ; coordz=coordz+s_wraf;  
        cubit.cmd(str1) 
 
coordx=0.5*s_weft;  coordy=T; coordz=-22*s_wraf; 
for i in range(0, 48): 
    if i % 2 == 0: 
        str1='create vertex ' + str(coordx) +' '+str(coordy)+' '+str(coordz)+' color' ; coordz=coordz+s_wraf;  
        cubit.cmd(str1) 
    if i % 2 == 1: 
        str1='create vertex ' + str(coordx) +' '+str(-coordy)+' '+str(coordz)+' color' ; coordz=coordz+s_wraf;  
        cubit.cmd(str1) 
 
Tvertices=48;  
vertices=range(1,Tvertices+1); count1=0;   
for i in range(0, 4): 
    str1='create curve spline vertex '
    for j in range(0, 48):
        str1=str1+' '+ str(count1+vertices[j])
    cubit.cmd(str1);  count1=count1+48;
    
cubit.cmd('create planar surface with plane normal to curve 1 distance 0 from vertex 1') 
str1='create curve location vertex 1 direction curve 8 length '+ str(W_wrap/2); cubit.cmd(str1) 
str1='create curve location vertex 1 direction curve 7 length  '+str(H_wrap/2); cubit.cmd(str1) 
cubit.cmd('create surface ellipse vertex 198 200 1')
cubit.cmd('sweep surface 2 along curve 1')
 
cubit.cmd('create planar surface with plane normal to curve 2 distance 0 from vertex 49') 
str1='create curve location vertex 49 direction curve 17 length '+ str(W_wrap/2); cubit.cmd(str1) 
str1='create curve location vertex 49 direction curve 16 length  '+str(H_wrap/2); cubit.cmd(str1) 
cubit.cmd('create surface ellipse vertex 208 210 49')
cubit.cmd('sweep surface 6 along curve 2')
 
cubit.cmd('create planar surface with plane normal to curve 3 distance 0 from vertex 97')
str1='create curve location vertex 97 direction curve 24 length '+ str(W_weft/2); cubit.cmd(str1)
str1='create curve location vertex 97 direction curve 23 length  '+str(H_weft/2); cubit.cmd(str1)
cubit.cmd('create surface ellipse vertex 218 220 97')
cubit.cmd('sweep surface 10 along curve 3')
 
cubit.cmd('create planar surface with plane normal to curve 4 distance 0 from vertex 145')
str1='create curve location vertex 145 direction curve 33 length '+ str(W_weft/2); cubit.cmd(str1)
str1='create curve location vertex 145 direction curve 32 length  '+str(H_weft/2); cubit.cmd(str1)
cubit.cmd('create surface ellipse vertex 228 230 145')
cubit.cmd('sweep surface 14 along curve 4')
 
cubit.cmd('delete vertex all')
cubit.cmd('delete curve all')
cubit.cmd('delete body 7 3 5 1')
 
str1='brick x '+str(L_RVE)+' y '+str(H_RVE)+' z '+str(W_RVE); cubit.cmd(str1) 
str1='move Volume '+str(9)+' x '+str(0)+' y '+str(0)+' z '+str(s_wraf/2); cubit.cmd(str1) 
 
cubit.cmd('intersect volume all keep') 
cubit.cmd('delete volume 6 2 8 4') 
cubit.cmd('subtract volume 10 11 12 13 from volume 9 keep') 
cubit.cmd('delete volume 9') 
cubit.cmd('imprint volume 10 11 12 13 14')   
cubit.cmd('merge volume 10 11 12 13 14') 
 
 
### =============== Meshing ================
cubit.cmd('group 1 add surface 42 27 24')  
cubit.cmd('group "g2" add surface 44 25 28')  
cubit.cmd('group "g3" add surface 39 30 33')  
cubit.cmd('group "g4" add surface 31 40 34')  
cubit.cmd('group "g5" add surface 43')  
cubit.cmd('group "g6" add surface 41')  
 
 
cubit.cmd('surface 42 27 24 size auto factor 6')       
cubit.cmd('surface 42 27 24 scheme trimesh')      
cubit.cmd('mesh surface 42 27 24')      
cubit.cmd('surface 28 25 44 scheme mirror source surface 27 24 42 source vertex 263 target vertex 262 nosmoothing')  
cubit.cmd('mesh surface 28 25 44')  
 
cubit.cmd('surface 39 30 33 size auto factor 5')     
cubit.cmd('surface 39 30 33 scheme trimesh')     
cubit.cmd('mesh surface 39 30 33')     
cubit.cmd('surface 31 34 40 scheme mirror source surface 30 33 39 source vertex 257 target vertex 262 nosmoothing') 
cubit.cmd('mesh surface 31 34 40') 
  
cubit.cmd('surface 43 scheme trimesh')    
cubit.cmd('mesh surface 43')    
cubit.cmd('surface 41 scheme mirror source surface 43 source vertex 259 target vertex 260 nosmoothing')
cubit.cmd('mesh surface 41')
 
##============================================================= 
##Mesh volume 
##============================================================= 
cubit.cmd('volume all scheme Tetmesh')  
cubit.cmd('mesh volume all')  
  
vol=['14', '10,11,12,13']         
mat=['MAT_ELASTIC_MATRIX','MAT_ELASTIC_FIBRES']     
Elastic=['5.0e3', '0.3', '5.0e3', '0.3']   #E in N/mm^2 for geometry dimensions in mm
count=0; 
for i in range(0, 2):         
    cubit.cmd('set duplicate block elements on')         
    str1='block ' + str(i+1) +' volume '+vol[i]; cubit.cmd(str1)         
    str1='block ' + str(i+1) +' name "'+mat[i] + '"'; cubit.cmd(str1)      
    str1='block ' + str(i+1) +' attribute count 2'; cubit.cmd(str1)      
    str1='block ' + str(i+1) +' attribute index 1 '+str(Elastic[count]); cubit.cmd(str1)      
    str1='block ' + str(i+1) +' attribute index 2 '+str(Elastic[count+1]); cubit.cmd(str1)      
    count=count+2
 
vol=['14', '10,11,12,13']         
mat=['MAT_MOISTURE','MAT_MOISTURE_LOW_DIFF']     
Diff=[0.01, 0.01]   
block_id=3
for i in range(0, 2):         
    cubit.cmd('set duplicate block elements on')         
    str1='block ' + str(block_id) +' volume '+vol[i]; cubit.cmd(str1)         
    str1='block ' + str(block_id) +' name "'+mat[i] + '"'; cubit.cmd(str1)      
    str1='block ' + str(block_id) +' attribute count 1'; cubit.cmd(str1)      
    str1='block ' + str(block_id) +' attribute index 1 '+str(Diff[i]); cubit.cmd(str1)
    block_id=block_id+1      
 
#Surfaces         
cubit.cmd('sideset 101 surface 41 31 34 40 42 27 24')  # all -ve boundary surfaces for periodic boundary conditions         
cubit.cmd('sideset 102 surface 43 30 33 39 25 44 28')  # all +ve boundary surfaces  for periodic boundary conditions            
cubit.cmd('sideset 103 surface 41 31 34 40 42 27 24 43 30 33 39 25 44 28')  # all boundary surfaces         
 
cubit.cmd('save as "/Users/zahur/Documents/moFEM/mofem-cephas/mofem_v0.2/users_modules/homogenisation/meshes/Textile_RVE.cub" overwrite')