mixed_modeI_III.jou

Go to the documentation of this file.

reset
#{angle=45}
create vertex -130  -5 0
create vertex  -130   5 0
 
create vertex  130   -5 0
create vertex  130   5 0
 
#{gama=angle*PI/180}
#{x=5/tan(gama)}
 
create vertex  {x}      -5 0
create vertex  {-1*x}  5 0 
 
 
create curve vertex 6 vertex 2
create curve vertex 2 vertex 1
create curve vertex 1 vertex 5
create curve vertex 5 vertex 6
 
 
create curve vertex 5 vertex 3
create curve vertex 3 vertex 4
create curve vertex 4 vertex 6
 
create surface curve 1 2 3 4
create surface curve 5 6 7 4
 
sweep surface 1 vector 0 0 1 distance 20
sweep surface 2 vector 0 0 1 distance 20
 
 
brick x 260 y 10 z 40 
move volume 3 x 0 y 0 z 40
split surface 13 across location position 0 -5 60 location position 0 5 60
 
split surface 7   across location position -120 -5 0 location position -120 5 0
split surface 12 across location position   120 -5 0 location position   120 5 0
 
 
brick x 50 y 10 z 60 
move volume 4 x 0 y 0 z 30
volume 4 copy
volume 4 copy
 
chop volume 1 with volume 4
chop volume 2 with volume 5
chop volume 3 with volume 6
 
imprint volume all  
merge volume all  
  
  
block 1 volume all  
block 1 name 'MAT_ELASTIC'  
block 1 attribute count 2  
block 1 attribute index 1 2800 
block 1 attribute index 2 0.38  
  
#surface 6 size auto factor 7 
#surface 6 scheme trimesh  
#mesh surface 6  
  
volume 7 9 11 scheme Tetmesh  
volume 7 9 11  size auto factor 7
mesh volume 7 9 11   
 
 
volume 8 10 12 13 scheme Tetmesh  
volume 8 10 12 13  size auto factor 8
mesh volume 8 10 12 13   
 
 
create displacement on curve 42 dof 1 dof 2 dof 3 fix 0 
create displacement on curve 47 dof 3 dof 2 fix 0 
create force on curve 37 vector 0 0 10 0 0 0 
 
#actual geometry boundries curves on which you don't allow nodes to move   
sideset 100  curve 92 100 101 44 43 2 94 108 110 118 117 51 50 6 21 18 36 34 25 10 13 35 33 27 37     
  
#actual geometry boundries vertices on which you don't allow nodes to move (we also include vertices of the force curve because we don't allow it move)  
nodeset 101 vertex 25 26 4 3 5 6 2 1 28 27 
  
#actual geometry boundries surfaces from which you don't allow nodes to move out (nodes can only move on the surface)  
sideset 102  surface 16 4 50 78 67 44 56 61 73 18 9 71 66 75 51 46 54 59 24 60 55 45 49 21       
  
#The edges between 102 and 10200,10201,.... will allow nodes to move out from surface 102  
Sideset 10200  surface 77 65 69 74  
  
#crack surface  
Sideset 200  surface 6  
  
#crack front  
Sideset 201 curve 16  
 
#color Volume 7 mediumturquoise 
#color Volume 11 mediumturquoise 
#color Volume 12 mediumturquoise 
#color Volume 13 mediumturquoise 
#color Volume 8 mediumturquoise 
#color Volume 10 mediumturquoise 
#color Volume 9 mediumturquoise 
 
color Volume all white 
color Surface 6 firebrick
 
#ADJUST FOR THE NEW VERSION OF MESH CUTTING
 
block 2 curve 100 118 110 6 94 35 117 44 51 43 34 13 2 101 21 10 36 50 25 18 27 33 148 144 128 133 141 138 131 129 
block 3 vertex 4 3 26 25 27 28 2 1
 
delete sideset 102 200 201 10200 100
 
 
create sheet extended from surface 6
move curve 157 midpoint location curve 16 include_merged
surface 82 scale x 1 y 1 z 0.5 about vertex 17
 
surface 82 size auto factor 4
surface 82 scheme trimesh
mesh surface 82
#move surface 82 z 0.1 include_merged
 
sideset 400 surface 82
 
unite volume 9 7  
 
imprint volume 9  with volume 11
merge volume all 
 
volume 9 scheme Tetmesh  
volume 9  size auto factor 7
mesh volume 9   
 
 
save as "/Users/zahur/Desktop/mixed_modeI_III_75.cub" overwrite 
 
hardcopy "/Users/zahur/Desktop/mesh_45deg1.bmp" bmp