<html>
  <head>
    <meta http-equiv="Content-Type" content="text/html; charset=UTF-8">
  </head>
  <body>
    <p>Hello,</p>
    <p>Your Surface is not "in" the volume but on its boundary, your
      surface is counted twice.</p>
    <p>Try BooleanFragments, this code should work:</p>
    <p>SetFactory("OpenCASCADE");<br>
      Box(1) = {0, 0, 0, 100, 10, 10};<br>
      Rectangle(7) = {40,0,0,20,10,0};<br>
      BooleanFragments{Surface{7}; Delete;}{Volume{1}; Delete;}<br>
      <br>
      Physical Volume("beam") = {1};<br>
      Physical Surface("left") = {8};<br>
      Physical Surface("right") = {9};<br>
      Physical Surface("load") = {7};<br>
    </p>
    <p>Best,</p>
    <p><br>
    </p>
    <div class="moz-cite-prefix">Le 13/02/2020 à 12:35, Vincent
      Libertiaux a écrit :<br>
    </div>
    <blockquote type="cite"
cite="mid:CAKQcerjL-OFagHvh0isGo66E4HmmPbDhRePfvrkOw3qtZWybpQ@mail.gmail.com">
      <meta http-equiv="content-type" content="text/html; charset=UTF-8">
      <div dir="ltr">
        <div>Hi, gmsh community,</div>
        <div><br>
        </div>
        <div>I created a simple case (beam fixed at both extremities)
          with a distributed load on a fraction of its upper surface.  I
          reproduce the short *.geo file:</div>
        <div><br>
        </div>
        <div>SetFactory("OpenCASCADE");<br>
          Box(1) = {0, 0, 0, 100, 10, 10};<br>
          Rectangle(7) = {40,0,0,20,10,0};<br>
          Surface{7} In Volume{1};<br>
          surfaces[] = Boundary{Volume{1};};<br>
          Physical Volume("beam") = {1};<br>
          Physical Surface("left") = {surfaces[0]};<br>
          Physical Surface("right") = {surfaces[1]};<br>
          Physical Surface("load") = {7};</div>
        <div><br>
        </div>
        <div>When I try to mesh it, I get the following error:</div>
        <div>Error   : PLC Error:  A segment and a facet intersect at
          point<br>
          Error   : Invalid boundary mesh (segment-facet intersection)
          on surface 7, intersection (40,4.7255,6.14662e-11)</div>
        <div><br>
        </div>
        <div>I think it comes from the fact that I am trying to impose a
          surface on another, but I can't find an easy way (rather than
          defining the geometry with Points, Lines, ... or determining
          the nodes on which I'd like to apply the loading from the msh
          file).</div>
        <div><br>
        </div>
        <div>Could someone please point me the right direction ?</div>
        <div><br>
        </div>
        <div>Thanks,</div>
        <div><br>
        </div>
        <div>kind regards,</div>
        <div><br>
        </div>
        <div>Vincent<br>
        </div>
      </div>
      <br>
      <fieldset class="mimeAttachmentHeader"></fieldset>
      <pre class="moz-quote-pre" wrap="">_______________________________________________
gmsh mailing list
<a class="moz-txt-link-abbreviated" href="mailto:gmsh@onelab.info">gmsh@onelab.info</a>
<a class="moz-txt-link-freetext" href="http://onelab.info/mailman/listinfo/gmsh">http://onelab.info/mailman/listinfo/gmsh</a>
</pre>
    </blockquote>
    <pre class="moz-signature" cols="72">-- 
Bertrand THIERRY
CNRS researcher
Laboratoire Jacques Louis Lions
Office 310, Corridor 15-25
+33 1 44 27 72 01
<a class="moz-txt-link-freetext" href="http://bthierry.pages.math.cnrs.fr/">http://bthierry.pages.math.cnrs.fr/</a></pre>
  </body>
</html>