Salome HOME
RNC: EDF 1618 Sketcher improvement for arc creation
authorgdd <gdd>
Wed, 24 Nov 2010 15:16:13 +0000 (15:16 +0000)
committergdd <gdd>
Wed, 24 Nov 2010 15:16:13 +0000 (15:16 +0000)
doc/salome/gui/GEOM/images/neo-scetcher1.png
doc/salome/gui/GEOM/images/neo-scetcher2.png
doc/salome/gui/GEOM/input/creating_sketcher.doc
src/GEOMGUI/GEOM_msg_en.ts
src/GEOMGUI/GEOM_msg_fr.ts
src/GEOM_SWIG/GEOM_Sketcher.py
src/GEOM_SWIG/geompyDC.py
src/SKETCHER/Sketcher_Profile.cxx
src/SKETCHER/Sketcher_Profile.hxx

index 2205805bd2b40532ed1167ad781e527648cdc084..1f387e952c5e6624a010e31d72b530e10fdc70b4 100755 (executable)
Binary files a/doc/salome/gui/GEOM/images/neo-scetcher1.png and b/doc/salome/gui/GEOM/images/neo-scetcher1.png differ
index cec5c57b409baaa632d36fcac07a9d9499645c74..bf03122778fee22d7e542d72f2aa8d925c91f120 100755 (executable)
Binary files a/doc/salome/gui/GEOM/images/neo-scetcher2.png and b/doc/salome/gui/GEOM/images/neo-scetcher2.png differ
index c93f5f3ff98eee98becd776196f133ac96bef203..bad3c2c356cb284911ee31c2664156fa3f07d345 100644 (file)
@@ -49,7 +49,7 @@ is predefined and is equal to 90 degrees);</li>
 <li>tangent to the previous segment;</li>
 <li>vector components DX, DY.</li>
 </ul>
-<li>Parameters of an element (for segment : length or target X or Y coordinate value, for arc : radius and angle).</li>
+<li>Parameters of an element (for segment : length or target X or Y coordinate value, for arc : radius, angle or center coordinates).</li>
 </ol>
 
 \b Buttons:
index 21d98dd464bf64cb5aee3e2a35b2e4a2040b90a3..895af0c3f6272858154d9747039757d03c796559 100644 (file)
@@ -1514,6 +1514,10 @@ Please, select face, shell or solid and try again</translation>
         <message>
             <source>GEOM_SKETCHER_ABS</source>
             <translation>Absolute</translation>
+       </message>
+        <message>
+            <source>GEOM_SKETCHER_ADD_PARAMS</source>
+            <translation>Additionnal Parameters</translation>
         </message>
         <message>
             <source>GEOM_SKETCHER_ANGLE</source>
@@ -1530,8 +1534,32 @@ Please, select face, shell or solid and try again</translation>
         <message>
             <source>GEOM_SKETCHER_ARC</source>
             <translation>Arc</translation>
-        </message>
+       </message>
         <message>
+            <source>GEOM_SKETCHER_CENTER</source>
+            <translation>Center</translation>
+       </message>
+       <message>
+            <source>GEOM_SKETCHER_CENTER2</source>
+            <translation>Center :</translation>
+       </message>  
+       <message>
+            <source>GEOM_SKETCHER_CENTER_X</source>
+            <translation>Center X coord. :</translation>
+       </message>
+       <message>
+            <source>GEOM_SKETCHER_CENTER_Y</source>
+            <translation>Center Y coord. :</translation>
+       </message>
+       <message>
+            <source>GEOM_SKETCHER_CENTER_DX</source>
+            <translation>Center DX coord. :</translation>
+       </message>
+       <message>
+            <source>GEOM_SKETCHER_CENTER_DY</source>
+            <translation>Center DY coord. :</translation>
+       </message>
+       <message>
             <source>GEOM_SKETCHER_DEST</source>
             <translation>Destination</translation>
         </message>
@@ -1566,6 +1594,10 @@ Please, select face, shell or solid and try again</translation>
         <message>
             <source>GEOM_SKETCHER_LENGTH2</source>
             <translation>Length :</translation>
+       </message>
+       <message>
+            <source>GEOM_SKETCHER_NONE</source>
+            <translation>None (Tangential)</translation>
         </message>
         <message>
             <source>GEOM_SKETCHER_PER</source>
@@ -1578,7 +1610,15 @@ Please, select face, shell or solid and try again</translation>
         <message>
             <source>GEOM_SKETCHER_POINT2</source>
             <translation>Point :</translation>
-        </message>
+       </message>
+       <message>
+            <source>GEOM_SKETCHER_END_POINT2</source>
+            <translation>End Point :</translation>
+       </message>
+       <message>
+            <source>GEOM_SKETCHER_RADIUS</source>
+            <translation>Radius</translation>
+        </message> 
         <message>
             <source>GEOM_SKETCHER_RADIUS2</source>
             <translation>Radius :</translation>
index 62dc3b1621415e0f47893554ec61c60b039299cb..71587db9530a1c8c948537b2e20f862e19d480fb 100644 (file)
         <message>
             <source>GEOM_SKETCHER_POINT2</source>
             <translation>Point :</translation>
+       </message>
+       <message>
+            <source>GEOM_SKETCHER_END_POINT2</source>
+            <translation>End Point :</translation>
         </message>
         <message>
             <source>GEOM_SKETCHER_RADIUS2</source>
index 5a8edfacee61903e345249247c3c8b3ff0c7337d..824f758bd2ba180822ea1e92df4bf0614d9cb595 100644 (file)
@@ -49,6 +49,19 @@ import geompy
 #"C radius length" : Create by direction, radius and length(in degree)
 #"AA x y": Create by point at X & Y
 #"A dx dy" : Create by point with DX & DY
+#"UU x y radius flag1": Create by point at X & Y with given radiUs 
+#"U dx dy radius flag1" : Create by point with DX & DY with given radiUs
+#"EE x y xc yc flag1 flag2": Create by point at X & Y with given cEnter 
+#"E dx dy dxc dyc radius flag1 flag2" : Create by point with DX & DY with given cEnter
+
+#Flag1 is 0 or 2 ... 
+# if 0 the drawn arc is the one of lower angle (< Pi)
+# if 2 the drawn arc ius the one of greater angle (> Pi)
+
+#Flag2 is 0 or 1 ...
+# if 0 the specified end point can be at a distance of the arc greater than the tolerance (10^-7)
+# if 1 the wire is built only if the end point is on the arc 
+#   with a tolerance of 10^-7 on the distance else the creation fails
 
 #To finish
 #"WW" : Close Wire
index ce18f9e69b35e556b2f60b4dd5c0a4c8d49b36a3..6ccb27bad9fd6bb25ea6b5f067fbdb5815a9c84b 100644 (file)
@@ -913,11 +913,27 @@ class geompyDC(GEOM._objref_GEOM_Gen):
         #     - "C radius length" : Create arc by direction, radius and length(in degree)
        #     - "AA x y": Create arc by point at X & Y
         #     - "A dx dy" : Create arc by point with DX & DY
-        #     .
+        #     - "A dx dy" : Create arc by point with DX & DY
+       #     - "UU x y radius flag1": Create arc by point at X & Y with given radiUs 
+       #     - "U dx dy radius flag1" : Create arc by point with DX & DY with given radiUs
+       #     - "EE x y xc yc flag1 flag2": Create arc by point at X & Y with given cEnter coordinates 
+       #     - "E dx dy dxc dyc radius flag1 flag2" : Create arc by point with DX & DY with given cEnter coordinates
+       #     .
         #       \n
         #     - "WW" : Close Wire (to finish)
         #     - "WF" : Close Wire and build face (to finish)
-        #
+        #     .
+       #       \n
+       #  - Flag1 (= reverse) is 0 or 2 ... 
+       #     - if 0 the drawn arc is the one of lower angle (< Pi)
+       #     - if 2 the drawn arc ius the one of greater angle (> Pi)
+       #     .
+       #       \n
+       #  - Flag2 (= control tolerance) is 0 or 1 ...
+       #     - if 0 the specified end point can be at a distance of the arc greater than the tolerance (10^-7)
+       #     - if 1 the wire is built only if the end point is on the arc 
+       #       with a tolerance of 10^-7 on the distance else the creation fails
+       #
         #  @param theCommand String, defining the sketcher in local
         #                    coordinates of the working plane.
         #  @param theWorkingPlane Nine double values, defining origin,
index 5e04c891e2ce15c5fbaa4389f9ff39a1bfb6df6c..d7c46982fd93aa55f92f38f5fd89df4b5ed0a8a0 100644 (file)
@@ -67,6 +67,7 @@ Sketcher_Profile::Sketcher_Profile()
 //=======================================================================
 Sketcher_Profile::Sketcher_Profile(const char* aCmd)
 {
+  myErrMsg = "";
   enum {line, circle, point, none} move;
 
   Standard_Integer i = 1;
@@ -78,6 +79,9 @@ Sketcher_Profile::Sketcher_Profile(const char* aCmd)
   first = Standard_True;
   stayfirst = face = close = Standard_False;
 
+  Standard_Integer reversed = 0;
+  Standard_Integer control_Tolerance = 0;
+
   TopoDS_Shape S;
   TopoDS_Vertex MP;
   BRepBuilderAPI_MakeWire MW;
@@ -185,7 +189,7 @@ Sketcher_Profile::Sketcher_Profile(const char* aCmd)
         }
       case 'T':
         {
-          if (n1 != 3) goto badargs;
+         if (n1 != 3) goto badargs;
           Standard_Real vx = a(1).RealValue();
           Standard_Real vy = a(2).RealValue();
           if (a(0) == "TT") {
@@ -245,7 +249,7 @@ Sketcher_Profile::Sketcher_Profile(const char* aCmd)
             move = none;
           break;
         }    
-      case 'A':
+      case 'A':                                // TAngential arc by end point   
         { 
           if (n1 != 3) goto badargs;
          Standard_Real vx = a(1).RealValue();
@@ -254,16 +258,14 @@ Sketcher_Profile::Sketcher_Profile(const char* aCmd)
             vx -= x;
             vy -= y;
           }
-         Standard_Real det = (dx * vy - dy * vx);
+         Standard_Real det = dx * vy - dy * vx;
          if ( Abs(det) > Precision::Confusion()) {
-           Standard_Real c = (dx * vx + dy * vy)   // Cosine of alpha = arc angle / 2
-                             / Sqrt((dx * dx + dy * dy)    
-                                  * (vx * vx + vy * vy));  
-            radius = (vx * vx + vy * vy)       // radius = distance between start  
-                    * Sqrt(dx * dx + dy * dy) // and end point / 2 * sin(alpha)
-                    / (2.0 * det);            // radius is > 0 or < 0
+           Standard_Real c = (dx * vx + dy * vy)                                            
+                             / Sqrt((dx * dx + dy * dy) * (vx * vx + vy * vy));                 // Cosine of alpha = arc of angle / 2 , alpha in [0,Pi]
+            radius = (vx * vx + vy * vy)* Sqrt(dx * dx + dy * dy)                                // radius = distance between start and end point / 2 * sin(alpha)  
+                    / (2.0 * det);                                                              // radius is > 0 or < 0
            if (Abs(radius) > Precision::Confusion()) {
-             angle = 2.0 * acos(c);           // angle in [0,Pi]  
+             angle = 2.0 * acos(c);                                                             // angle in [0,2Pi]  
              move = circle;
             }
            else
@@ -274,6 +276,83 @@ Sketcher_Profile::Sketcher_Profile(const char* aCmd)
             move = none;
           break;
         }       
+      case 'U':                                // Arc by end point and radiUs
+        { 
+         if (n1 != 5) goto badargs;
+         Standard_Real vx = a(1).RealValue();
+          Standard_Real vy = a(2).RealValue();
+         radius  = a(3).RealValue();
+         reversed = a(4).IntegerValue();
+         if (a(0) == "UU") {                 // Absolute
+            vx -= x;
+            vy -= y;
+          }
+         Standard_Real length = Sqrt(vx * vx + vy * vy);
+         if ( (4.0 - (vx * vx + vy * vy) / (radius * radius) >= 0.0 ) && (length > Precision::Confusion()) ) {
+           Standard_Real c = 0.5 * Sqrt(4.0 - (vx * vx + vy * vy) / (radius * radius));        // Cosine of alpha = arc angle / 2 , alpha in [0,Pi/2]
+           angle = 2.0 * acos(c);                                                              // angle in [0,Pi]
+           if ( reversed == 2 )
+             angle = angle - 2 * PI; 
+           dx =    0.5 * (  vy * 1.0/radius 
+                          + vx * Sqrt(4.0  / (vx * vx + vy * vy) - 1.0 / (radius * radius)));    
+            dy = -  0.5 * (  vx * 1.0/radius 
+                          - vy * Sqrt(4.0  / (vx * vx + vy * vy) - 1.0 / (radius * radius)));    
+           move = circle;
+          }
+         else{
+           move = none;
+         }
+         break;
+        }       
+      case 'E':                                // Arc by end point and cEnter
+        { 
+         myErrMsg = "";
+         if (n1 != 7) goto badargs;
+         Standard_Real vx = a(1).RealValue();
+          Standard_Real vy = a(2).RealValue();
+         Standard_Real vxc  = a(3).RealValue();
+         Standard_Real vyc  = a(4).RealValue();
+         reversed = a(5).IntegerValue();
+         control_Tolerance = a(6).IntegerValue();
+
+         if (a(0) == "EE") {                 // Absolute
+            vx -= x;
+            vy -= y;
+           vxc -= x;
+           vyc -= y; 
+          }
+         radius = Sqrt( vxc * vxc + vyc * vyc );
+         Standard_Real det = vx * vyc - vy * vxc;
+         Standard_Real length = Sqrt(vx * vx + vy * vy);
+         Standard_Real length2 = Sqrt((vx-vxc) * (vx-vxc) + (vy-vyc) * (vy-vyc));
+         Standard_Real length3 = Sqrt(vxc * vxc + vyc * vyc);
+         Standard_Real error = Abs(length2 - radius);
+         if ( error > Precision::Confusion() ){
+           MESSAGE("Warning : The specified end point is not on the Arc, distance = "<<error);
+           myErrMsg = "Warning : The specified End Point is not on the Arc";
+         }
+          if ( error > Precision::Confusion() &&                                                // Don't create the arc if the end point 
+              control_Tolerance == 1)                                                          // is too far from it
+           move = none;
+         else if ( (length > Precision::Confusion()) && 
+                   (length2 > Precision::Confusion()) && 
+                   (length3 > Precision::Confusion()) ) {
+           Standard_Real c = ( radius * radius - (vx * vxc + vy * vyc) ) 
+                             / ( radius * Sqrt((vx-vxc) * (vx-vxc) + (vy-vyc) * (vy-vyc)) ) ;  // Cosine of arc angle 
+           angle = acos(c);                                                                    // angle in [0,Pi] 
+           if ( reversed == 2 )
+             angle = angle - 2 * PI;
+           if (det < 0)
+             angle = -angle; 
+           dx =  vyc / radius;
+           dy = -vxc / radius; 
+           move = circle;
+         }
+         else {
+           move = none;
+         }
+         break;
+        }      
       case 'I':
         {
           if (n1 != 2) goto badargs;
index 29801915f08550d706de9d41f0f5d1fb58bfa484..e68e2ebfed2f4e04c791fdea3c929912d7c51f4f 100644 (file)
@@ -29,6 +29,7 @@
 #include <gp_Pnt.hxx>
 #include <gp_Dir.hxx>
 #include <TopoDS_Shape.hxx>
+#include <string>
 
 class Sketcher_Profile
 {
@@ -43,6 +44,7 @@ private:
 
   TopoDS_Shape myShape;
   bool myOK;
+  std::string myErrMsg;
 
 public:
   Standard_EXPORT gp_Pnt GetLastPoint(){return myLastPoint;};
@@ -50,5 +52,6 @@ public:
 
   Standard_EXPORT const TopoDS_Shape& GetShape(){return myShape;};
   Standard_EXPORT bool IsDone(){return myOK;};
+  Standard_EXPORT std::string ErrMsg(){return myErrMsg;};
 
 };