<div dir="ltr"><div>Hello Bertrand,</div><div>thank you for your answer; your solution worked perfectly.</div><div><br></div><div>kind regards,</div><div><br></div><div>Vincent<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Le jeu. 13 févr. 2020 à 16:48, Bertrand Thierry <<a href="mailto:bertrand.thierry@sorbonne-universite.fr">bertrand.thierry@sorbonne-universite.fr</a>> a écrit :<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
  
    
  
  <div>
    <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>Le 13/02/2020 à 12:35, Vincent
      Libertiaux a écrit :<br>
    </div>
    <blockquote type="cite">
      
      <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></fieldset>
      <pre>_______________________________________________
gmsh mailing list
<a href="mailto:gmsh@onelab.info" target="_blank">gmsh@onelab.info</a>
<a href="http://onelab.info/mailman/listinfo/gmsh" target="_blank">http://onelab.info/mailman/listinfo/gmsh</a>
</pre>
    </blockquote>
    <pre cols="72">-- 
Bertrand THIERRY
CNRS researcher
Laboratoire Jacques Louis Lions
Office 310, Corridor 15-25
+33 1 44 27 72 01
<a href="http://bthierry.pages.math.cnrs.fr/" target="_blank">http://bthierry.pages.math.cnrs.fr/</a></pre>
  </div>

_______________________________________________<br>
gmsh mailing list<br>
<a href="mailto:gmsh@onelab.info" target="_blank">gmsh@onelab.info</a><br>
<a href="http://onelab.info/mailman/listinfo/gmsh" rel="noreferrer" target="_blank">http://onelab.info/mailman/listinfo/gmsh</a><br>
</blockquote></div>