tutorials/solidMechanics/elasticAcpSolidFoam/crackingBiMatDcbLinear/constant/polyMesh/blockMeshDict.simpleDCB.m4

   1 /*--------------------------------*- C++ -*----------------------------------*\
   2 | =========                 |                                                 |
   3 | \\      /  F ield         | foam-extend: Open Source CFD                    |
   4 |  \\    /   O peration     | Version:     3.2                                |
   5 |   \\  /    A nd           | Web:         http://www.foam-extend.org         |
   6 |    \\/     M anipulation  |                                                 |
   7 \*---------------------------------------------------------------------------*/
   8 FoamFile
   9 {
  10     version     2.0;
  11     format      ascii;
  12     class       dictionary;
  13     location    "constant/polyMesh";
  14     object      blockMeshDict;
  15 }
  16 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  17
  18 // Model Description
  19 // Plane strain 2-D model of simple DCB
  20 // two beams adhered together
  21 // thickness of each beam can be set
  22 // or the material of each beam could be set using setFields
  23
  24 // Setup m4 stuff
  25 changecom(//)changequote([,])
  26 define(calc, [esyscmd(perl -e 'printf ($1)')])
  27 //define(VCOUNT, 0)
  28 //define(vlabel, [[// ]Vertex $1 = VCOUNT define($1, VCOUNT)define([VCOUNT], incr(VCOUNT))])
  29
  30 // define geometry in mm
  31 define(bl, 60) // beams length
  32 define(bh, 10) // individual beam height
  33 define(pn, 10) // pre-notch length
  34 define(z, 1) // out of plane thickness
  35
  36 // define mesh density
  37 define(meshPn, 10) // number of cells in beam length direction before notch
  38 define(meshBl, 60) // number of cells in beam length direction after notch
  39 define(meshBh, 10)  // number of cells in beam height direction
  40
  41 // start of blockMeshDict
  42
  43 convertToMeters 0.001;
  44
  45 vertices
  46 (
  47  //- dimension in mm
  48  (0 -bh 0)
  49  (pn -bh 0)
  50  (bl -bh 0)
  51  (bl 0 0)
  52  (bl bh 0)
  53  (pn bh 0)
  54  (0 bh 0)
  55  (0 0 0)
  56  (pn 0 0)
  57  (0 0 0) // 9
  58
  59  (0 -bh z)
  60  (pn -bh z)
  61  (bl -bh z)
  62  (bl 0 z)
  63  (bl bh z)
  64  (pn bh z)
  65  (0 bh z)
  66  (0 0 z)
  67  (pn 0 z)
  68  (0 0 z) // 19
  69
  70 );
  71
  72 blocks
  73 (
  74 //  hex (0 1 8 9 10 11 18 19) sheet (meshPn meshBh 1) simpleGrading (1 1 1)
  75 //  hex (1 2 3 8 11 12 13 18) sheet (meshBl meshBh 1) simpleGrading (1 1 1)
  76 //  hex (7 8 5 6 17 18 15 16) sheet (meshPn meshBh 1) simpleGrading (1 1 1)
  77 //  hex (8 3 4 5 18 13 14 15) sheet (meshBl meshBh 1) simpleGrading (1 1 1)
  78  hex (0 1 8 9 10 11 18 19) (meshPn meshBh 1) simpleGrading (1 1 1)
  79  hex (1 2 3 8 11 12 13 18) (meshBl meshBh 1) simpleGrading (1 1 1)
  80  hex (7 8 5 6 17 18 15 16) (meshPn meshBh 1) simpleGrading (1 1 1)
  81  hex (8 3 4 5 18 13 14 15) (meshBl meshBh 1) simpleGrading (1 1 1)
  82  );
  83
  84 edges
  85 (
  86 );
  87
  88 patches
  89 (
  90  patch tractionFree
  91  (
  92   (15 16 6 5)
  93   (14 15 5 4)
  94
  95   (0 1 11 10)
  96   (1 2 12 11)
  97
  98   (2 3 13 12)
  99   (3 4 14 13)
 100
 101   (7 8 18 17)
 102   (19 18 8 9)
 103   )
 104
 105  patch topLoading
 106  (
 107   (17 16 6 7)
 108   )
 109
 110  patch bottomLoading
 111  (
 112   (10 19 9 0)
 113   )
 114
 115  empty back
 116  (
 117   (7 6 5 8)
 118   (8 5 4 3)
 119   (0 9 8 1)
 120   (1 8 3 2)
 121   )
 122
 123  empty front
 124  (
 125   (10 11 18 19)
 126   (11 12 13 18)
 127   (17 18 15 16)
 128   (18 13 14 15)
 129   )
 130
 131  cohesive crack ()
 132  );
 133
 134 mergePatchPairs
 135 (
 136 );
 137
 138 // ************************************************************************* //