Refactoring/MoveSomeFunctionsToBeamLib:

spostate le funzioni GetTunnelDimension e CalcCollisionSafety da LapJoint a BeamLib
2022-09-01 15:24:08 +02:00
parent a9e43984a2
commit b0e10cee87
3 changed files with 90 additions and 163 deletions
@@ -814,5 +814,78 @@ function BeamLib.FindFaceBestOrientedAsAxis( Proc, vtAx, fctExclude)
   return nFaceIndMax, ptC, vtN
 end

+---------------------------------------------------------------------
+function BeamLib.GetTunnelDimension( Proc, nPartId)
+   -- sono necessarie almeno due facce
+   if Proc.Fct < 2 then return 0, 0, 0 end
+   -- recupero l'ingombro della trave
+   local b3Solid = EgtGetBBoxGlob( EgtGetFirstNameInGroup( nPartId, 'Box') or GDB_ID.NULL, GDB_BB.STANDARD)
+   -- recupero centro e normale delle facce
+   local ptC = {}
+   local vtN = {}
+   for i = 1, Proc.Fct do
+      ptC[i], vtN[i] = EgtSurfTmFacetCenter( Proc.Id, i - 1, GDB_ID.ROOT)
+   end
+   -- calcolo l'orientamento del tunnel
+   local vtOrtho
+   local bAdj = EgtSurfTmFacetsContact( Proc.Id, 0, 1)
+   if bAdj then
+      vtOrtho = vtN[1] ^ vtN[2]
+   else
+      if Proc.Fct >= 3 then
+         vtOrtho = vtN[1] ^ vtN[3]
+      else
+         return 0, 0, 0
+      end
+   end
+   -- ottengo il boundingBox e prendo le dimensioni lungo la direzione (Z locale) che rappresenta la profondità della fessura
+   local frFc = Frame3d( ptC[1], vtOrtho) ;
+   local b3BoxLoc = EgtGetBBoxRef( Proc.Id, GDB_BB.STANDARD, frFc)
+   local dDepth = b3BoxLoc:getDimZ()
+   -- recupero gruppo per geometria addizionale
+   local nAddGrpId = BL.GetAddGroup( nPartId)
+   if not nAddGrpId then
+      EgtOutLog( 'Error : missing AddGroup')
+      return 0, 0, 0
+   end
+   -- centro del bounding box locale
+   local ptCen = b3BoxLoc:getCenter()
+   ptCen:toGlob( frFc)
+   -- creo superficie intermedia
+   local nSurfInt = EgtSurfTmPlaneInBBox( nAddGrpId, ptCen, vtOrtho, b3Solid, GDB_ID.ROOT)
+   if not nSurfInt then return 0, 0, 0 end
+   -- ritaglio la superficie con le facce della fessura
+   for i = 1, Proc.Fct do
+      EgtCutSurfTmPlane( nSurfInt, ptC[i], -vtN[i], false, GDB_ID.ROOT)
+   end
+   local frSurfInt, dDimMax, dDimMin = EgtSurfTmFacetMinAreaRectangle( nSurfInt, 0, GDB_ID.ROOT)
+   -- cerco la faccia con larghezza pari a dimensione massima della fessura
+   local nLongIdFace = 0
+   for i = 1, Proc.Fct do
+      if abs( vtN[i] * frSurfInt:getVersX()) < 0.5 then
+         nLongIdFace = i - 1
+         break
+      end
+   end
+   return dDimMin, dDimMax, dDepth, vtOrtho, nLongIdFace, nSurfInt
+end
+
+---------------------------------------------------------------------
+function BeamLib.CalcCollisionSafety( vtDir)
+   local dCollSic = 10 * BD.COLL_SIC
+   if abs( vtDir:getX()) > 0.999 or abs( vtDir:getY()) > 0.999 or abs( vtDir:getZ()) > 0.999 then
+      dCollSic = 0
+   elseif abs( vtDir:getX()) > 0.996 or abs( vtDir:getY()) > 0.996 or abs( vtDir:getZ()) > 0.996 then
+      dCollSic = 1 * BD.COLL_SIC
+   elseif abs( vtDir:getX()) > 0.89 or abs( vtDir:getY()) > 0.89 or abs( vtDir:getZ()) > 0.89 then
+      dCollSic = 2.5 * BD.COLL_SIC
+   elseif abs( vtDir:getX()) > 0.86 or abs( vtDir:getY()) > 0.86 or abs( vtDir:getZ()) > 0.86 then
+      dCollSic = 4 * BD.COLL_SIC
+   elseif abs( vtDir:getX()) > 0.707 or abs( vtDir:getY()) > 0.707 or abs( vtDir:getZ()) > 0.707 then
+      dCollSic = 5.5 * BD.COLL_SIC
+   end
+   return dCollSic
+end
+
 -------------------------------------------------------------------------------------------------------------
 return BeamLib