<html>
<head>
<meta content="text/html; charset=UTF-8" http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
Thank you, Errol. I also tried creating the geometry in BricsCAD and
their exported iges file was fine. I guess the issue was with
FreeCAD that I had initially used.<br>
<br>
Anyhow, if anybody is interested, attached is the "manual"
definition of this geometry with a non-through-going hole. As you
can see, it's a lot of work to get the geometry right!<br>
<br>
Down the road, it may be a useful feature to be able to do some
basic 3D solid editing, such as union, subtract, etc. <br>
<br>
Cheers,<br>
Omid<br>
<br>
<br>
<div class="moz-cite-prefix">On 14-09-04 02:06 PM,
<a class="moz-txt-link-abbreviated" href="mailto:erroll.brady@gmail.com">erroll.brady@gmail.com</a> wrote:<br>
</div>
<blockquote
cite="mid:b21d6cd1-05da-4c1e-8448-b3135ab59eee@googlegroups.com"
type="cite">
<div dir="ltr">
<div>Hi Omid,</div>
Solidworks found errors when trying to import your original iges
file; it did not fix the file. The recent object was newly
created in Solidworks and exported to a STEP file. The STEP file
was loaded into Gmsh where it was 3D-meshed and saved as a mesh
and exported as a GEO file.
<div>Erroll<br>
<br>
On Wednesday, September 3, 2014 11:31:25 AM 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>
This is exactly what I wanted. Did Solidworks
automatically find and correct the errors? Do you have any
other suggestions for modifying the iges files?<br>
<br>
Thanks,<br>
Omid<br>
<br>
<div>On 14-09-03 01:03 PM, <a moz-do-not-send="true"
href="javascript:" target="_blank"
gdf-obfuscated-mailto="NqDHbogeD4wJ"
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">Hi Omid,
<div>Try the attached Block2.zip.</div>
<p style="text-align:center;clear:both"><span
style="font-size:13px">Erroll</span></p>
<p style="text-align:left;clear:both"><a
moz-do-not-send="true"
href="https://lh5.googleusercontent.com/-EkE2FS91LsI/VAdJvqh3QUI/AAAAAAAAAAw/pfpGUvqDVdQ/s1600/Block2.jpg"
style="margin-left:1em;margin-right:1em"
target="_blank"
onmousedown="this.href='https://lh5.googleusercontent.com/-EkE2FS91LsI/VAdJvqh3QUI/AAAAAAAAAAw/pfpGUvqDVdQ/s1600/Block2.jpg';return
true;"
onclick="this.href='https://lh5.googleusercontent.com/-EkE2FS91LsI/VAdJvqh3QUI/AAAAAAAAAAw/pfpGUvqDVdQ/s1600/Block2.jpg';return
true;"><img moz-do-not-send="true"
origsrc="https://lh5.googleusercontent.com/-EkE2FS91LsI/VAdJvqh3QUI/AAAAAAAAAAw/pfpGUvqDVdQ/s320/Block2.jpg"
border="0" height="274" width="320"></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
moz-do-not-send="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
moz-do-not-send="true"
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 moz-do-not-send="true"
border="0" height="307" width="320"></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:PGRpdj48ZGl2PllvdSBjb3VsZCBkZWxldGUgdGhlIHZvbHVtZSBvZiB0aGUgY3lsaW5kZXIgeW91J3ZlIGNyZWF0ZWQgYnkgbWVhbnMgb2YgZXh0cnVzaW9uIGFuZCB0aGVuIGRlbGV0ZSBvbmUgb2YgdGhlIGN5bGluZGVyJ3MgZmxhdCBzdXJmYWNlcyB1c2luZyB0aGUgYERlbGV0ZWAgY29tbWFuZC4gSG93ZXZlciwgaXQgd29uJ3QgZ2l2ZSB5b3Ugd2hhdCB5b3Ugd2FudCwgYmVjYXVzZSB0aGUgYm94
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
moz-do-not-send="true">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 moz-do-not-send="true" 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 moz-do-not-send="true" value="+16474789767">+1-647-478-9767</a>
Cell: <a moz-do-not-send="true" value="+14168242408">+1-416-824-2408</a>
</pre>
</font></span></div>
<br>
______________________________<wbr>_________________<br>
gmsh mailing list<br>
<a moz-do-not-send="true">gm...@geuz.org</a><br>
<a moz-do-not-send="true"
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 moz-do-not-send="true"
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
moz-do-not-send="true"
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 moz-do-not-send="true" 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 moz-do-not-send="true" 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 moz-do-not-send="true" 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 class="moz-signature" cols="72">--
Omid Mahabadi, Ph.D.
Geomechanica, Inc.
<a class="moz-txt-link-freetext" href="http://www.geomechanica.com/">http://www.geomechanica.com/</a>
Tel : +1-647-478-9767
Cell: +1-416-824-2408
</pre>
</body>
</html>