Cutting a triangular mesh and keeping the pieces

-
When we use a plane to cut a triangular mesh, several cases can be considered. If there is a non-null intersection between the plane and a given triangle, there are mostly five different cases:


  1. The plane contains one vertex of the triangle: not very useful for cutting purposes if the rest of the triangle lays above the cutting plane and we want to keep all that is below.
  2. The plane contains two vertices of the triangle, then the intersection is a line segment but we will only keep that triangle if the other vertex is below the cutting plane.
  3. The plane contains three vertices of the triangle: that is a keeper.
  4. The plane does not contain any vertices of the triangle, but there is an intersection with two of the edges: this is what we will talk about today here.
  5. A fifth case can be considered when the plane contains one vertex but the two other are on opposite sides of the plane, and it can easily be incorporated to the previous case.
The following picture shows us a couple of triangles of a mesh that is going to be cut by a plane:
 The dotted horizontal line represents the cutting plane. If we want to create a new mesh of all the model that lays below the cutting plane, the triangles shown will need to be divided in a top and a bottom part.

But depending on how many vertices of the triangle are above the cutting plane, the resultant polygon of cutting the triangle by the cutting plane may be a smaller triangle, as shown in the left in pale orange or a quad as shown on the right case.
But if we want to keep the resulting mesh as a triangular mesh, the quads that might be obtained will need to be converted into triangles too.

We can split the quad into two triangles easily, but what may not be as easy to do is to maintain the proper winding order of the triangle nodes. As that is the base of many programs for determining the normal direction of the triangle, failing to restore the correct order for these new triangles will result in a messy mesh.

What we have to remember is that these new triangles obtained by dividing a quad share the same orientation (normal) as the original triangle they are coming from. So a simple thing to do can be to calculate the new normal for the new triangles and, if found too far away from the original normal, two vertices of the triangle can be swapped to restore the right winding order. That approach worked for me nicely. 

Of course, another task that needs to be achieved to obtain a 2-manifold mesh is to cover the holes created by the cutting plane in the remaining mesh below, but that will be covered in an upcoming entry.

Comments

Post a Comment

Popular posts from this blog

VFD control with Arduino using RS485 link

-
Image
Variable Frequency Drives ( VFD ) allow the control of spindles so speed can accurately be controlled and a detailed acceleration profile for the spindle and reverse rotation can all be handled. In essence a VFD is an three phase inverter for three-phase AC motor. I am using a popular (I mean cheap) Chinese VFD and though the reference manual is not great, I could see there is a built-in RS-485 port. I usually control the start, stop and speed selection using the keyboard on the unit but I thought it will be more useful if I could control everything from the same Arduino is doing our CNC table control. Some cheap RS-485 off eBay and some lines of code later, I can start, stop and change the speed from an Arduino. What a cool thing to have! Some configuration of your VFD are needed before you can use it like that.  You need to set PD163=1 (I am using address 1 in the code). PD=164=1 (for setting serial to 9600bps) and PD165=0 (for using ASCII and 8N1 character format)
Read more

How to get sinusoidal s-curve for a stepper motor

-
Image
After some experimentation, and delving a bit into the math of motion, I realized there can be a significant difference on the motion of a machine by just making small changes. A common and popular approach to smooth one axis of motion is to use a trapezoidal speed pattern, where acceleration and deceleration are uniformly accelerated movements (aka constant acceleration motion) separated by a portion happening at cruising speed. The problem comes with the fact that, while that approach is simple to understand and to code, it does contain sudden changes in the acceleration, that show as spikes on the acceleration derivative, called jerk. These sudden changes in the acceleration translate into undesired vibration in our machines. Most of that can be eliminated by selecting a motion pattern that does not contain sudden changes of acceleration. If we chose a mathematical function for our speed whose second derivative is continuous, then we reduce system vibration quite a lot. One ma
Read more

Importing OpenSCAD designs into Onshape

-
Image
After using Onshape CAD software for a while I am really fond of it. So one natural thing to do, specially once they have released the Instructor's kit is to start using it in the classroom too. But one thing I was not so sure how to do was to use old designs I made using OpenSCAD. It does not mean I am completely quitting using OpenSCAD but suddenly I can see I can do assemblies in Onshape that will look quite nice for showing and documentation purposes if only I could easily import from OpenSCAD. I am glad to report that the route OpenSCAD --> FreeCAD --> Onshape worked brilliantly. Even more when a few tricks were applied. While some people suggest to import CSG files from OpenSCAD to FreeCAD, I have had more success by directly importing SCAD files into FreeCAD. One trick that worked well for me was to get rid of $fn references, so now circles become real circles (or cylinders) in FreeCAD. If a part keeps some $fn references then your cylinders or circles wi
Read more