v0.14.0
torsion.jou
Go to the documentation of this file.
1 reset
2 
3 brick x 100 y 100 z 400
4 
5 brick x 125 y 125 z 25
6 brick x 210 y 25 z 25
7 unite volume 2 3
8 
9 move Volume 2 x 0 y 0 z -187.50 include_merged
10 split curve 41 distance 5 from vertex 18
11 split curve 42 distance 5 from vertex 21
12 split curve 44 distance 5 from vertex 24
13 split curve 45 distance 5 from vertex 19
14 create curve vertex 35 vertex 36 on surface 21
15 create curve vertex 34 vertex 33 on surface 20
16 imprint tolerant surface 21 with curve 65 merge
17 imprint tolerant surface 20 with curve 66 merge
18 delete curve 65 66
19 split curve 39 distance 5 from vertex 17
20 split curve 38 distance 5 from vertex 22
21 split curve 47 distance 5 from vertex 20
22 split curve 48 distance 5 from vertex 23
23 create curve vertex 41 vertex 42 on surface 19
24 create curve vertex 43 vertex 44 on surface 22
25 imprint tolerant surface 19 with curve 77 merge
26 imprint tolerant surface 22 with curve 78 merge
27 delete curve 77 78
28 Volume 2 copy move x 0 y 0 z 375
29 unite volume all
30 
31 brick x 100 y 100 z 400
32 brick x 400 y 400 z 400
33 subtract volume 5 from volume 6
34 Compress Ids All
35 chop volume 1 with volume 2
36 Compress Ids All
37 
38 create Cylinder height 400 radius 50
39 rotate Volume 4 angle 90 about Y include_merged
40 Compress Ids All
41 move Volume 4 location surface 34 include_merged
42 rotate Volume 4 angle -45 about Y include_merged
43 chop volume 3 with volume 4
44 webcut volume 5 with general plane plane x spin axis y angle 45 noimprint nomerge
45 Compress Ids All
46 
47 #Boundary conditions
48 merge volume all group_results
49 split curve 37 38 61 47 48 49 53 57 fraction 0.5
50 create Displacement on vertex 110 dof 1 dof 2 dof 3 fix 0
51 create Displacement on vertex 105 dof 1 dof 2 dof 3 fix 0
52 create Displacement on vertex 107 dof 1 dof 2 dof 3 fix 0
53 create force on vertex 112 force value -1 direction y
54 
55 #Create Block Set with Vertices to ploat load-disp. path
56 block 3 vertex 112 97 98
57 block 3 name 'LoadPath'
58 
59 #Merge Volumes
60 imprint volume all
61 merge volume all
62 
63 #Material Properties
64 block 1 volume 3 4 5
65 block 1 name "MAT_ELASTIC_Concrete"
66 block 1 attribute count 10
67 block 1 attribute index 1 3.5
68 block 1 attribute index 2 0.2
69 
70 block 2 volume 1 2
71 block 2 name "MAT_ELASTIC_Steel"
72 block 2 attribute count 10
73 block 2 attribute index 1 20
74 block 2 attribute index 2 0.3
75 
76 #Create nodesets and sidesets for fracture constraints
77 #Boundary edges
78 sideset 100 curve all
79 sideset 100 curve 156 145 143 146 154 142 144 132 131 remove
80 sideset 100 curve 185 198 remove
81 
82 #Boundary vertices
83 nodeset 101 vertex all
84 nodeset 101 vertex 89 90 99 100 91 92 remove
85 
86 #Boundary surfaces
87 sideset 102 surface all
88 sideset 102 surface 61 53 54 38 39 45 44 36 37 42 43 remove
89 
90 sideset 102 surface 75 remove
91 sideset 10200 surface 75
92 
93 #Crack surface
94 sideset 200 surface 53
95 
96 #Crack front
97 sideset 201 curve 143
98 
99 #!Meshing
100 volume all scheme tetmesh
101 volume all size auto factor 7
102 #!Mesh Crack Surface (those elemenst are streched)
103 surface 53 size auto factor 7
104 
105 #!Mesh All
106 mesh surface 53
107 mesh volume all
108 
109 save as "/Users/likask/MyBuild/mofem-bitbucket/mofem/meshes/torsion.cub" overwrite
110 
Generated by Doxygen 1.8.17 and hosted at University of Glasgow