<div dir="ltr">Hi Omid,<div>Try the attached Block2.zip.</div><p class="separator" style="text-align: center; clear: both;"><span style="font-size: 13px;">Erroll</span></p><p class="separator" style="text-align: left; clear: both;"><a imageanchor="1" href="https://lh5.googleusercontent.com/-EkE2FS91LsI/VAdJvqh3QUI/AAAAAAAAAAw/pfpGUvqDVdQ/s1600/Block2.jpg" style="margin-left: 1em; margin-right: 1em;"><img src="https://lh5.googleusercontent.com/-EkE2FS91LsI/VAdJvqh3QUI/AAAAAAAAAAw/pfpGUvqDVdQ/s320/Block2.jpg" border="0" style="" width="320" height="274"></a></p><div><br><br>On Tuesday, September 2, 2014 9:39:45 PM UTC-7, Omid Mahabadi wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;">
  
    
  
  <div bgcolor="#FFFFFF" text="#000000">
    Hi Errol,<br>
    In the version you sent me, the hole is through-going, in other
    words, it goes through the entire depth of the cube (along z).
    However, my desired geometry requires a hole that doesn't go through
    the entire depth (say it has half of the depth of the exterior
    cube).<br>
    <br>
    Thank you,<br>
    Omid<br>
    <br>
    <br>
    <div>On 14-09-02 06:04 PM,
      <a href="javascript:" target="_blank" gdf-obfuscated-mailto="yRkuJvL7JFcJ" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">erroll...@gmail.com</a> wrote:<br>
    </div>
    <blockquote type="cite">
      <div dir="ltr">
        <div>Hello Omid,</div>
        <div>Solidworks found errors in your IGES file. Try the attached
          STEP, GEO and MSH files.</div>
        <p style="text-align:center;clear:both"><a href="https://lh4.googleusercontent.com/-nxF0eT1U3qQ/VAY-5mA30KI/AAAAAAAAAAg/aQ6Sqk53Fp8/s1600/Block%2BMesh.jpg" style="margin-left:1em;margin-right:1em" target="_blank" onmousedown="this.href='https://lh4.googleusercontent.com/-nxF0eT1U3qQ/VAY-5mA30KI/AAAAAAAAAAg/aQ6Sqk53Fp8/s1600/Block%2BMesh.jpg';return true;" onclick="this.href='https://lh4.googleusercontent.com/-nxF0eT1U3qQ/VAY-5mA30KI/AAAAAAAAAAg/aQ6Sqk53Fp8/s1600/Block%2BMesh.jpg';return true;"><img border="0" height="307" width="320" origsrc="https://lh4.googleusercontent.com/-nxF0eT1U3qQ/VAY-5mA30KI/AAAAAAAAAAg/aQ6Sqk53Fp8/s320/Block%2BMesh.jpg"></a></p>
        <div><br>
        </div>
        <div>Erroll</div>
        <div><br>
          On Tuesday, September 2, 2014 1:28:00 PM UTC-7, Omid Mahabadi
          wrote:
          <blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex">
            <div bgcolor="#FFFFFF" text="#000000"> Hi Oliver,<br>
              Thanks for your response. I had also tried the variation
              that you sent me. The issue is that I don't want to have
              the flat surface that is interfacing the two volumes. I'd
              like the mesh to be continuous across that interface.<br>
              <br>
              Also, when following your suggestion, I still cannot get
              the same exact mesh as you did. Please see the attached
              screenshot. <br>
              <br>
              Is there any other thoughts on how to get the results I'm
              aiming for?<br>
              <br>
              Thank you,<br>
              Omid<br>
              <br>
              <div>On 14-08-28 06:34 PM, Oliver Willekens wrote:<br>
              </div>
              <blockquote type="cite">
                <div dir="ltr">
                  <div>
                    <p style="margin:1.2em 0px!important">You could
                      delete the volume of the cylinder you’ve created
                      by means of extrusion and then delete one of the
                      cylinder’s flat surfaces using the <code style="font-size:0.85em;font-family:Consolas,Inconsolata,Courier,monospace;margin:0px 0.15em;padding:0px 0.3em;white-space:pre-wrap;border:1px solid rgb(234,234,234);background-color:rgb(248,248,248);border-radius:3px;display:inline">Delete</code>
                      command. However, it won’t give you what you want,
                      because the box “holding the cylinder” will be
                      unaware of the presence of the cylindrical hole.</p>
                    <p style="margin:1.2em 0px!important">I tried this:</p>
                    <pre style="font-size:0.85em;font-family:Consolas,Inconsolata,Courier,monospace;font-size:1em;line-height:1.2em;margin:1.2em 0px"><code style="font-size:0.85em;font-family:Consolas,Inconsolata,Courier,monospace;margin:0px 0.15em;padding:0px 0.3em;white-space:pre-wrap;border:1px solid rgb(234,234,234);background-color:rgb(248,248,248);border-radius:3px;display:inline;white-space:pre;overflow:auto;border-radius:3px;border:1px solid rgb(204,204,204);padding:0.5em 0.7em;display:block!important;display:block;padding:0.5em;color:rgb(51,51,51);background:none repeat scroll 0% 0% rgb(248,248,255)">Plane Surface(1) = {1, 2};  // square minus the circle
out1[] = Extrude {0.0, 76.2, 0.0}{Surface{1};};
Plane Surface(2) = {1};
out2[] = Extrude {0.0, -76.2, 0.0} {Surface{2};};
</code></pre>
                    <p style="margin:1.2em 0px!important">Which results
                      in the 2D mesh you find in the attachment. Perhaps
                      this is what you’re looking for?</p>
                    <p style="margin:1.2em 0px!important">I also tried
                      combining those two volumes using <code style="font-size:0.85em;font-family:Consolas,Inconsolata,Courier,monospace;margin:0px 0.15em;padding:0px 0.3em;white-space:pre-wrap;border:1px solid rgb(234,234,234);background-color:rgb(248,248,248);border-radius:3px;display:inline">Compound

                        Volume(4) = {out1[1], out2[1]};</code>, but like
                      you, I got GRegion Compound errors. Probably some
                      basic restriction from triangulations I’m
                      overlooking. This only means your 3D meshes will
                      be clearly stopped at a flat interface about
                      half-way through those two boxes, but you’ll still
                      have a 3D mesh. I was hoping that the Compound
                      Volume would make some tetras across this boundary
                      though… But with the above lines, you’ll have a
                      working mesh.</p>
                    <div title="MDH:PGRpdj48ZGl2PllvdSBjb3VsZCBkZWxldGUgdGhlIHZvbHVtZSBvZiB0aGUgY3lsaW5kZXIgeW91J3ZlIGNyZWF0ZWQgYnkgbWVhbnMgb2YgZXh0cnVzaW9uIGFuZCB0aGVuIGRlbGV0ZSBvbmUgb2Yg
dGhlIGN5bGluZGVyJ3MgZmxhdCBzdXJmYWNlcyB1c2luZyB0aGUgYERlbGV0ZWAgY29tbWFuZC4g
SG93ZXZlciwgaXQgd29uJ3QgZ2l2ZSB5b3Ugd2hhdCB5b3Ugd2FudCwgYmVjYXVzZSB0aGUgYm94
ICJob2xkaW5nIHRoZSBjeWxpbmRlciIgd2lsbCBiZSB1bmF3YXJlIG9mIHRoZSBwcmVzZW5jZSBv
ZiB0aGUgY3lsaW5kcmljYWwgaG9sZS48YnI+Cjxicj48L2Rpdj5JIHRyaWVkIHRoaXM6PGJyPmBg
YDxicj5QbGFuZSBTdXJmYWNlKDEpID0gezEsIDJ9OyZuYnNwOyAvLyBzcXVhcmUgbWludXMgdGhl
IGNpcmNsZTxicj5vdXQxW10gPSBFeHRydWRlIHswLjAsIDc2LjIsIDAuMH17U3VyZmFjZXsxfTt9
Ozxicj5QbGFuZSBTdXJmYWNlKDIpID0gezF9Ozxicj5vdXQyW10gPSBFeHRydWRlIHswLjAsIC03
Ni4yLCAwLjB9IHtTdXJmYWNlezJ9O307PGJyPgpgYGA8YnI+PC9kaXY+V2hpY2ggcmVzdWx0cyBp
biB0aGUgMkQgbWVzaCB5b3UgZmluZCBpbiB0aGUgYXR0YWNobWVudC4gUGVyaGFwcyB0aGlzIGlz
IHdoYXQgeW91J3JlIGxvb2tpbmcgZm9yPzxicj48YnI+SSBhbHNvIHRyaWVkIGNvbWJpbmluZyB0
aG9zZSB0d28gdm9sdW1lcyB1c2luZyBgQ29tcG91bmQgVm9sdW1lKDQpID0ge291dDFbMV0sIG91
dDJbMV19O2AsIGJ1dCBsaWtlIHlvdSwgSSBnb3QgR1JlZ2lvbiBDb21wb3VuZCBlcnJvcnMuIFBy
b2JhYmx5IHNvbWUgYmFzaWMgcmVzdHJpY3Rpb24gZnJvbSB0cmlhbmd1bGF0aW9ucyBJJ20gb3Zl
cmxvb2tpbmcuIFRoaXMgb25seSBtZWFucyB5b3VyIDNEIG1lc2hlcyB3aWxsIGJlIGNsZWFybHkg
c3RvcHBlZCBhdCBhIGZsYXQgaW50ZXJmYWNlIGFib3V0IGhhbGYtd2F5IHRocm91Z2ggdGhvc2Ug
dHdvIGJveGVzLCBidXQgeW91J2xsIHN0aWxsIGhhdmUgYSAzRCBtZXNoLiBJIHdhcyBob3Bpbmcg
dGhhdCB0aGUgQ29tcG91bmQgVm9sdW1lIHdvdWxkIG1ha2Ugc29tZSB0ZXRyYXMgYWNyb3NzIHRo
aXMgYm91bmRhcnkgdGhvdWdoLi4uIEJ1dCB3aXRoIHRoZSBhYm92ZSBsaW5lcywgeW91J2xsIGhh

dmUgYSB3b3JraW5nIG1lc2guPGJyPgo8ZGl2Pjxicj48L2Rpdj4=" style="min-height:0;font-size:0em;padding:0;margin:0">​</div>
                  </div>
                </div>
                <div><br>
                  <br>
                  <div class="gmail_quote">2014-08-28 17:08 GMT+02:00
                    Omid Mahabadi <span dir="ltr"><<a>omid.m...@geomechanica.com</a>></span>:<br>
                    <blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
                      <div bgcolor="#FFFFFF" text="#000000"> Hi
                        Christophe and Gmsh team,<br>
                        I've been trying to model a simple cube with a
                        cylindrical hole that is NOT through-going
                        (i.e., its depth is smaller than the depth of
                        the cube), as shown in the attached picture. The
                        Extrude command doesn't seem to work since it
                        will extrude both surfaces (exterior and
                        interior -- hole) at the same time. I know I can
                        define all the surfaces, surface loops, and
                        volumes manually, but is there a better way of
                        defining something like this in Gmsh? <br>
                        <br>
                        I also tried to use the Compound Volume command
                        by first defining two volumes from Extrude and
                        then trying to combine them but I'm getting
                        errors for the Compound Volume visualization
                        (Error: Cannot evaluate bounds on GRegion
                        Compound) and when I mesh the geometry, the
                        actual shared interfaces are still existing,
                        although by the notion of compound from the
                        documentation, the internal interfaces should be
                        neglected. Here is the geometry file:<br>
                        <blockquote>// Characteristic length (==>
                          element size)<br>
                          cl_external   = 25;<br>
                          cl_excavation = 5;<br>
                          <br>
                          // External boundaries<br>
                          Point(1)      = {-127, 0.0, -127,
                          cl_external};<br>
                          Point(2)      = {+127, 0.0, -127,
                          cl_external};<br>
                          Point(3)      = {+127, 0.0, +127,
                          cl_external};<br>
                          Point(4)      = {-127, 0.0, +127,
                          cl_external};<br>
                          Line(1)       = {1, 2};<br>
                          Line(2)       = {2, 3};<br>
                          Line(3)       = {3, 4};<br>
                          Line(4)       = {4, 1};<br>
                          Line Loop(1)  = {1, 2, 3, 4};<br>
                          <br>
                          // Excavation boundaries<br>
                          Point(5)      = {0.0, 0.0, 0.0,
                          cl_excavation};<br>
                          Point(6)      = {19.05, 0.0, 0.0,
                          cl_excavation};<br>
                          Point(7)      = {0.0, 0.0, 19.05,
                          cl_excavation};<br>
                          Point(8)      = {-19.05, 0.0, 0.0,
                          cl_excavation};<br>
                          Point(9)      = {0.0, 0.0, -19.05,
                          cl_excavation};<br>
                          Circle(5)     = {6, 5, 7};<br>
                          Circle(6)     = {7, 5, 8};<br>
                          Circle(7)     = {8, 5, 9};<br>
                          Circle(8)     = {9, 5, 6};<br>
                          Line Loop(2)  = {5, 6, 7, 8};<br>
                          <br>
                          //Using compound volumes<br>
                          Plane Surface(1) = {1};<br>
                          out1[] = Extrude {0.0, 76.2,
                          0.0}{Surface{1};};<br>
                          <br>
                          Plane Surface(2) = {1};<br>
                          out2[] = Extrude {0.0, -76.2,
                          0.0}{Surface{2};};<br>
                          <br>
                          Compound Volume(3) = {1, 2};<br>
                        </blockquote>
                        <br>
                        I also tried to create the geometry in CAD
                        software and imported it as iges, brep or step
                        formats but for reason the hole is not there
                        completely. See for instances the iges file
                        attached.<br>
                        <br>
                        Can you kindly shed some light here? Am I doing
                        something wrong? Or are there Gmsh
                        tricks/commands that I can use to achieve my
                        goal?<br>
                        <br>
                        Thank you,<br>
                        Omid<span><font color="#888888"><br>
                            <br>
                            <pre cols="72">-- 
Omid Mahabadi, Ph.D.
Geomechanica, Inc.
<a href="http://www.geomechanica.com/" target="_blank" onmousedown="this.href='http://www.google.com/url?q\75http%3A%2F%2Fwww.geomechanica.com%2F\46sa\75D\46sntz\0751\46usg\75AFQjCNF6W3NuiF8qST2qDZz4P5VNTQsiVw';return true;" onclick="this.href='http://www.google.com/url?q\75http%3A%2F%2Fwww.geomechanica.com%2F\46sa\75D\46sntz\0751\46usg\75AFQjCNF6W3NuiF8qST2qDZz4P5VNTQsiVw';return true;">http://www.geomechanica.com/</a>
Tel : <a value="+16474789767">+1-647-478-9767</a>
Cell: <a value="+14168242408">+1-416-824-2408</a>
</pre>
                          </font></span></div>
                      <br>
                      ______________________________<wbr>_________________<br>
                      gmsh mailing list<br>
                      <a>gm...@geuz.org</a><br>
                      <a href="http://www.geuz.org/mailman/listinfo/gmsh" target="_blank" onmousedown="this.href='http://www.google.com/url?q\75http%3A%2F%2Fwww.geuz.org%2Fmailman%2Flistinfo%2Fgmsh\46sa\75D\46sntz\0751\46usg\75AFQjCNEJ3QLvAdHWrw4QzBDqKw4PAU6GzA';return true;" onclick="this.href='http://www.google.com/url?q\75http%3A%2F%2Fwww.geuz.org%2Fmailman%2Flistinfo%2Fgmsh\46sa\75D\46sntz\0751\46usg\75AFQjCNEJ3QLvAdHWrw4QzBDqKw4PAU6GzA';return true;">http://www.geuz.org/mailman/<wbr>listinfo/gmsh</a><br>
                      <br>
                    </blockquote>
                  </div>
                  <br>
                  <br clear="all">
                  <br>
                  -- <br>
                  Oliver Willekens<br>
                  PhD Student <br>
                  <a href="https://lcp.elis.ugent.be/" target="_blank" onmousedown="this.href='https://www.google.com/url?q\75https%3A%2F%2Flcp.elis.ugent.be%2F\46sa\75D\46sntz\0751\46usg\75AFQjCNEL0HUpZZ3TcsvAIBMVl6cseoC1Gw';return true;" onclick="this.href='https://www.google.com/url?q\75https%3A%2F%2Flcp.elis.ugent.be%2F\46sa\75D\46sntz\0751\46usg\75AFQjCNEL0HUpZZ3TcsvAIBMVl6cseoC1Gw';return true;"><img alt="LCP group
                      logo" style="float:left;margin-right:20px;margin-top:20px"></a>
                  <p style="margin-top:20px"> Liquid Crystals &
                    Photonics Group <br>
                    Sint- Pietersnieuwstraat 41 <br>
                    9000 Gent <br>
                    Phone: +32 9 264.89.51 </p>
                </div>
              </blockquote>
              <br>
              <pre cols="72">-- 
Omid Mahabadi, Ph.D.
Geomechanica, Inc.
<a href="http://www.geomechanica.com/" target="_blank" onmousedown="this.href='http://www.google.com/url?q\75http%3A%2F%2Fwww.geomechanica.com%2F\46sa\75D\46sntz\0751\46usg\75AFQjCNF6W3NuiF8qST2qDZz4P5VNTQsiVw';return true;" onclick="this.href='http://www.google.com/url?q\75http%3A%2F%2Fwww.geomechanica.com%2F\46sa\75D\46sntz\0751\46usg\75AFQjCNF6W3NuiF8qST2qDZz4P5VNTQsiVw';return true;">http://www.geomechanica.com/</a>
Tel : +1-647-478-9767
Cell: +1-416-824-2408
</pre>
            </div>
          </blockquote>
        </div>
      </div>
    </blockquote>
    <br>
    <pre cols="72">-- 
Omid Mahabadi, Ph.D.
Geomechanica, Inc.
<a href="http://www.geomechanica.com/" target="_blank" onmousedown="this.href='http://www.google.com/url?q\75http%3A%2F%2Fwww.geomechanica.com%2F\46sa\75D\46sntz\0751\46usg\75AFQjCNF6W3NuiF8qST2qDZz4P5VNTQsiVw';return true;" onclick="this.href='http://www.google.com/url?q\75http%3A%2F%2Fwww.geomechanica.com%2F\46sa\75D\46sntz\0751\46usg\75AFQjCNF6W3NuiF8qST2qDZz4P5VNTQsiVw';return true;">http://www.geomechanica.com/</a>
Tel : +1-647-478-9767
Cell: +1-416-824-2408
</pre>
  </div>

</blockquote></div></div>