Structure of problem.geof

The structure of the ASCII file describing problem geometries is summarized below:

***geometry
 **node
    n      dim                 number of nodes, space dimension
    x1     y1     [z1]
    ...                        coordinates
    xn     yn     [zn]
 **element
    j                          number of elements
    1   type1   n11  .... n1k  element number, type, list of node numbers
    ...
    j  typej  nj1  ....  njk

 ***group
  **nset name              name of a node set
    n1  n2  ...  nk            list of the k nodes in the node set

 **elset name              name of an element set
    n1  n2  ...  nk            list of the k elements in the element set

 **faset name              surface (3D face) group name
    type1  n11  ...  n1k       type of face, ordered list of nodes in the face
    ...
    typej  nj1  ...  njk

 **liset name              surface (linear) group name
    type1  n11  ...  n1k       type of line, ordered list of nodes in the liset
    ...
    typej  nj1  ...  njk

***return

The position of the nodes must be given in Cartesian coordinates for 2D plane and 3D problems, or in cylindrical coordinates (r,z) for 2D axisymmetric problems.

The groups of nodes, faces, elements, and element lines segments must come after the node and element definitions. These groups are defined between the keywords ***group and ***return.

Example

***geometry
 **node
   349 3
   1 2.000000000000000e+01 0.000000000000000e+00 0.000000000000000e+00
   2 3.000000000000000e+01 0.000000000000000e+00 0.000000000000000e+00
   3 2.121319960000000e+01 0.000000000000000e+00 2.121319960000000e+01
   ...
   349 2.926353200000000e+00 1.471177770000000e+01 2.598075870000000e+01
 **element
   48
   1 c3d20r 1 22 21 41 23 40 8 7 266 278 282 270 54 75 74 94 76 93 61 60
   2 c3d20r 54 75 74 94 76 93 61 60 267 279 283 271 107 128 127 147 129 146 114 113
   3 c3d20r 107 128 127 147 129 146 114 113 268 280 284 272 160 181 180 200 182 199 167 166
   ...
   48 c3d20r 198 212 190 189 162 174 175 210 337 349 345 341 251 265 243 242 215 227 228 263
***group
 **nset planexy
   1 54 107 160 213 266 267 268 269 7 60 113 166 219 8 61 114 167 220 270 271
   272 273 9 62 115 168 221 2 55 108 161 214 274 275 276 277
 **nset planexz
   1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21
   22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42
   43 44 45 46 47 48 49 50 51 52 53
 **nset planeyz
   213 214 215 216 5 6 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233
   234 235 236 237 238 239 240 241 242 243 244 245 246 35 36 37 250 251 252 253 254
   255 256 257 258 259 260 261 262 263 264 265
 **elset inner
   25 30 26 27 28 32 31 29 33 34 35 36 20 19 18 17 9 10 11 12 4
   3 2 1
 **elset outer
   5 13 21 45 41 37 40 44 48 24 16 8 7 6 43 42 46 47 23 22 14
   15 39 38
 **liset internal_line
   quad  4 18 19
   quad  19 20 21
   quad  21 22 1
   quad  5 25 26
   quad  26 27 28
   quad  28 29 4

 **faset internal_face
   q8  4 314 57 71 72 290 19 18
   q8  57 315 110 124 125 291 72 71
   q8  110 316 163 177 178 292 125 124
   ...

***return