REPOS_ERICCA/magicsld/SuperFace.vb



Public Class SuperFace
    Inherits SuperEntite

    'Public swFace As SldWorks.Face2
    Private Shared compteur As Long
    Private Shared no As Long

    Friend Flag As Integer = 0 ' = 20 si on a une coupeLong


    Friend lst_Faces As New Collections.Generic.List(Of SldWorks.Face2) ' une liste de faces supplémentaires si me.swFace est coupée en 3.
    Friend lst_InterPoutre As New Collection() ' liste des poutres qui intersectionnent
    Friend lst_InterCoqueVolume As New Collections.Generic.List(Of InterCoqueVolume)  ' liste des intersections avec d'autres faces
    Friend lst_InterCoqueCoque As New Collections.Generic.List(Of InterCoqueCoque)


    Friend AttributCL As SldWorks.Attribute ' une l'attribut de condition aux imites qui doit être updaté
    Private swSurface As SldWorks.Surface


    ''' <summary>
    ''' Renvoie le nombre de swFaces qui composent cette superface
    ''' </summary>
    ''' <value></value>
    ''' <returns></returns>
    ''' <remarks></remarks>
    Public ReadOnly Property GetNbFaces() As Long
        Get
            Return lst_Faces.Count
        End Get
    End Property


    ''' <summary>
    ''' New pour un encapsulateur temporaire de face
    ''' </summary>
    ''' <param name="Face"></param>
    ''' <param name="encapsulateur"></param>
    ''' <remarks></remarks>
    Public Sub New(ByRef Face As SldWorks.Face2, ByRef encapsulateur As Boolean)
        lst_Faces.Add(Face)
    End Sub

    Friend Sub New(ByRef face As SldWorks.Face2, Optional ByVal tip As Integer = 0)
        Me.AjouterFace(face)
        Select Case tip
            Case Commun.tipe_e.Volume
                nom = "Face" & compteur
            Case Commun.tipe_e.coque
                nom = "FaceCoque" & compteur
        End Select
        nomOrig = nom
        compteur = compteur + 1
    End Sub


    Friend Overridable Function UpdateApresSplit(ByRef inter As InterPoutreCoque, ByRef poutre As SlyAretePoutre, ByRef x As Double, ByRef y As Double, ByRef z As Double, ByRef Plan As SldWorks.RefPlane, Optional ByRef FI As Boolean = False) As SldWorks.Face2
        MsgBox("La fonction non overridée a été appelée!")
        Return Nothing
    End Function

    Friend Overridable Function UpdateApresSplit(ByRef inter As InterPoutreVolume, ByRef poutre As SlyAretePoutre, ByRef x As Double, ByRef y As Double, ByRef z As Double, ByRef Plan As SldWorks.RefPlane, Optional ByRef FI As Boolean = False) As SldWorks.Face2
        MsgBox("La fonction non overridée a été appelée!")
        Return Nothing
    End Function

    ''' <summary>
    ''' Sub qui ajoute des mini-poutres à la section entre le sommet de la poutre et un point de la face interne
    ''' </summary>
    ''' <param name="poutre">La poutre principale</param>
    ''' <param name="FaceInterne">La face où il faut ajouter UNE mini-poutre</param>
    ''' <param name="x1">Coordonnée en x du point d'intersection</param>
    ''' <param name="y1"></param>
    ''' <param name="z1"></param>
    ''' <remarks>On pourrait ajouter plus d'une mini-poutre.</remarks>
    Friend Sub AjouterMiniPoutresSurFaceInterne(ByRef poutre As SlyAretePoutre, ByRef FaceInterne As SldWorks.Face2, ByVal x1 As Double, ByVal y1 As Double, ByVal z1 As Double)
        Dim x2 As Double, y2 As Double, z2 As Double
        Dim swSketch As SldWorks.Sketch
        Static k As Integer
        Dim swEnt As SldWorks.Entity
        Dim p As SldWorks.Parameter
        Dim refCourbe As SldWorks.ReferenceCurve
        Dim feat As SldWorks.Feature
        Dim attr As SldWorks.Attribute

        '  2 trouver un sommet ou un point sur la face interne.
        Dim swArete As SldWorks.Edge
        Dim swSommet2 As SldWorks.Vertex

        swArete = FaceInterne.GetFirstLoop.GetFirstCoEdge.getedge
        swSommet2 = swArete.GetStartVertex()

        If swSommet2 Is Nothing Then ' un celcle (ou ellipse)
            Dim retval As Object
            retval = swArete.Evaluate(0) ' il va y avor un LC point anyway...
            x2 = retval(0)
            y2 = retval(1)
            z2 = retval(2)

        Else
            Dim retval As Object
            retval = swSommet2.GetPoint()
            x2 = retval(0)
            y2 = retval(1)
            z2 = retval(2)
        End If


        '  3 faire une mini-poutre entre les 2

        swModel.Insert3DSketch2(True)
        swModel.CreateLine2(x1, y1, z1, x2, y2, z2)
        swSketch = swModel.GetActiveSketch2()
        swModel.Insert3DSketch2(True)
        swEnt = swSketch : swEnt.Select2(False, 1)
        swModel.InsertCompositeCurve()

        ' reste à lui mettre les propriétés de mini-poutres
        feat = swModel.FeatureByPositionReverse(0)
        refCourbe = feat.GetSpecificFeature2
        swArete = refCourbe.GetFirstSegment ' y'a juste un segment

        swEnt = swArete
        Dim NomAttr As String
        NomAttr = "miniPoutre" & Me.nom & k

        attr = Intersections.DefAttrRCP1.CreateInstance5(swModel, swArete, NomAttr, 0, 0)
        ' mettre les propriétés de la mini-poutre
        ' attention, elle ne doit pas avoir de masse, sa section est alors Mini...

        p = attr.GetParameter("As")
        p.SetDoubleValue2(-7, 2, "")
        Commun.GererDossiers("Poutres", NomAttr)
        k += 1


    End Sub


    Friend Overridable Function AjouterInterFace(ByRef sPoutre As SlyAretePoutre, ByRef xyz1() As Double, ByVal tipe As Byte) As InterCoqueVolume
        ' sub qui CRÉÉ une instance de la classe InterFaceFace si et seulement si il n'en existe pas avant


        'Dim int As InterFaceFace

        ' créez la routine qui ignore la création de 
        Return Nothing ' pout l'instant
    End Function


    Shared Sub reinitialiser()
        compteur = 0
    End Sub

    Public Overrides Sub SaveNom()
        ' procédure qui enregistre le nom et qui , pour l'instant ne tient pas compte des conditions aux limites
        Dim ent As SldWorks.Entity
        Dim i As Long
        Dim nomtemp As String
        'For i = 1 To Me.lst_Faces.Count
        For Each swFace As SldWorks.Face2 In Me.lst_Faces
            ent = swFace
            Dim retval As Boolean
            If Me.lst_Faces.Count > 1 Then nomtemp = nom & "-" & Chr(i + 96) Else nomtemp = nom
            retval = swPart.SetEntityName(ent, nomtemp)
            If retval = False Then
                Dim nom2 As String
                'Dim swface As SldWorks.Face2
                'swface = ent
                'If Me.SwFace Is swface Then MsgBox("La même face!")
                nom2 = swPart.GetEntityName(ent)
                'MsgBox("Impossibilité d'écrire le nom de l'entité... il y a déjà un nom ( " & nom2 & " ) " & Chr(13) & "Et on veut écrire à la place --> " & nomtemp, MsgBoxStyle.Critical, "Problème dans le setID de SlyFaceVol")
                Debug.Print("Impossibilité d'écrire le nom de l'entité... il y a déjà un nom ( " & nom2 & " ) " & Chr(13) & "Et on veut écrire à la place --> " & nomtemp, MsgBoxStyle.Critical, "Problème dans le setID de SlyFaceVol")
            End If
        Next


    End Sub


    Public Function GetNormale(ByRef x As Double, ByRef y As Double, ByRef z As Double) As Double()
        Dim surf As SldWorks.Surface
        Dim retval As Object
        Dim temp(2) As Double

        surf = Me.lst_Faces.Item(0).GetSurface
        retval = surf.EvaluateAtPoint(x, y, z)

        If retval Is Nothing Then MsgBox("Erreur, dans le Getnormal de la face, le point ne semble pas être sur la face")

        temp(0) = retval(0)
        temp(1) = retval(1)
        temp(2) = retval(2)

        Return temp
    End Function
    Public Function GetNormaleSurface(ByRef x As Double, ByRef y As Double, ByRef z As Double) As Double()
        Dim surf As SldWorks.Surface
        Dim retval As Object = Nothing
        Dim temp(2) As Double

        surf = Me.lst_Faces.Item(0).GetSurface
        retval = surf.EvaluateAtPoint(x, y, z)

        temp(0) = retval(0)
        temp(1) = retval(1)
        temp(2) = retval(2)

        Return temp
    End Function


    Public Function GetRayonCourbureMax(ByRef x As Double, ByRef y As Double, ByRef z As Double) As Double
        Dim surf As SldWorks.Surface
        Dim retval As Object
        Dim temp As Double

        surf = Me.lst_Faces.Item(0).GetSurface
        retval = surf.EvaluateAtPoint(x, y, z)

        temp = retval(9)

        Return temp
    End Function

    Public Function GetRayonCourbureMin(ByRef x As Double, ByRef y As Double, ByRef z As Double) As Double
        Dim surf As SldWorks.Surface
        Dim retval As Object
        Dim temp As Double

        surf = Me.lst_Faces.Item(0).GetSurface
        retval = surf.EvaluateAtPoint(x, y, z)

        temp = retval(10)

        Return temp
    End Function

    ' si la surface est un cylindre, retourne son rayon.  Plante ou valeur aléatoire si pas un cylindre
    Public Function GetRayonCylindre() As Double
        Dim surf As SldWorks.Surface
        Dim retval As Object
        Dim temp As Double

        surf = Me.lst_Faces.Item(0).GetSurface
        retval = surf.CylinderParams

        temp = retval(6)

        Return temp
    End Function


    ''' <summary>
    ''' Retourne vrai si une face est plane
    ''' </summary>
    ''' <returns></returns>
    ''' <remarks></remarks>
    Public Function estPlan() As Boolean
        Dim surf As SldWorks.Surface

        surf = Me.lst_Faces.Item(0).GetSurface()

        If surf.IsPlane Then Return True Else Return False

    End Function


    ''' <summary>
    ''' Retourne vrai si la face est un plan mais pas vu comme tel par solidworks.
    ''' </summary>
    ''' <param name="x"></param>
    ''' <param name="y"></param>
    ''' <param name="z"></param>
    ''' <returns></returns>
    ''' <remarks></remarks>
    Public Function estFauxPlan(ByRef x As Double, ByRef y As Double, ByRef z As Double) As Boolean
        Dim surf As SldWorks.Surface
        surf = Me.lst_Faces.Item(0).GetSurface()
        If Not surf.IsParametric Then Return False

        ' si la normale est la même au point d'intersection et à quelques sommets, alors on a une face plane.
        Dim Normale(2) As Double
        Dim vNormale As Object

        vNormale = surf.EvaluateAtPoint(x, y, z) ' la normale du point d'intersection
        Normale(0) = vNormale(0)
        Normale(1) = vNormale(1)
        Normale(2) = vNormale(2)

        Dim i As Integer
        Dim U As Double, V As Double, Umin As Double, Umax As Double, Vmin As Double, Vmax As Double
        Dim retval As Object

        retval = surf.Parameterization()
        Umin = retval(0)
        Umax = retval(1)
        Vmin = retval(2)
        Vmax = retval(3)

        Randomize()

        For i = 0 To 5
            U = Rnd() * (Umax - Umin) + Umin
            V = Rnd() * (Vmax - Vmin) + Vmin
            vNormale = surf.Evaluate(U, V, 0, 0)
            If Not (Math.Abs(vNormale(3) - Normale(0)) < Epsilon And Math.Abs(vNormale(4) - Normale(1)) < Epsilon And Math.Abs(vNormale(5) - Normale(2)) < Epsilon) Then Return False
        Next i

        Return True


    End Function

    Public Function estCylindre() As Boolean
        Dim surf As SldWorks.Surface

        surf = Me.lst_Faces.Item(0).GetSurface

        If surf.IsCylinder Then Return True Else Return False

    End Function


    Protected Overrides Sub Finalize()
        ' mettre l'effacement des listes ici.

        MyBase.Finalize()
    End Sub


    Public Sub MettreAttributPourConditionLimite()
        Dim swent As SldWorks.Entity
        Dim nom As String
        Dim cond As String

        cond = Me.condition
        If cond = "" Then Exit Sub

        swent = Me.SwFace

        nom = Me.nom & "CL" & CStr(no) & "_" & cond
        Dim Attr As SldWorks.Attribute = Nothing

        Try
            Attr = swent.FindAttribute(Intersections.DefAttrConditionLimite, 0) ' si l'attribut existe déjà on pointe dessus.
        Catch ex As Exception
            MsgBox("N'arrive pas à se lier à l'attribut, erreur: " & ex.Message, MsgBoxStyle.Critical)
        End Try

        If Attr Is Nothing Then Attr = Intersections.DefAttrConditionLimite.CreateInstance5(swModel, Me.SwFace, nom, 0, 2) ' 0 = swThisconfig


        While Attr Is Nothing
            no += 1
            nom = "CL" & Me.nom & "_" & cond
            Attr = Intersections.DefAttrConditionLimite.CreateInstance5(swModel, Me.SwFace, nom, 0, 0)
        End While


        Dim ParamCL As SldWorks.Parameter
        ParamCL = Attr.GetParameter("CL")

        ParamCL.SetStringValue2(cond, 2, "")  ' swAllConfiguration = 2
        Me.AttributCL = Attr
        GererDossiers("Conditions Aux Limites", nom)
        no = no + 1

    End Sub

    ' une fonction qui transforme un attribut en condition aux limites
    Public Sub AttributVersConditionLimite()
        Dim p As SldWorks.Parameter
        Dim ent As SldWorks.Entity
        Dim attr As SldWorks.Attribute

        ent = Me.SwFace
        attr = ent.FindAttribute(Intersections.DefAttrConditionLimite, 0)
        If Not attr Is Nothing Then
            p = attr.GetParameter("CL")
            nom = nomOrig & "@" & p.GetStringValue
        End If

    End Sub


    ' sub qui découpe les bords de la face.
    Friend Sub CoupeCote(ByRef inter As InterPoutreVolume, ByRef poutre As SlyAretePoutre)
        Dim pt3() As Double, pt3Original() As Double
        Dim base(2) As Double, baseOriginal(2) As Double
        Dim swEnt As SldWorks.Entity
        Dim Directionnel As Boolean, Flip As Boolean
        Dim planReference As SldWorks.RefPlane = Nothing
        Dim sketchline As SldWorks.SketchSegment
        Dim swSketch As SldWorks.Sketch
        Dim DemiLargeur As Double
        Dim g As Integer
        Dim Face(1) As SldWorks.Face2
        Dim PlanEntity As SldWorks.Entity = Nothing
        Dim r(2) As Double
        Dim sk(1) As Double
        pt3Original = poutre.GetPoint3


        swModel.SetInferenceMode(False) ' ne pas mettre de contraintes par défaut
        'swModel.SetAddToDB(True)
        'swModel.SetDisplayWhenAdded(False) ' accélérer les performances

        Dim vArete As Object
        Dim cut As Double

        If Me.estPlan Then
            ' la coque est plane, on met une esquisse dessus.#
            PlanEntity = Me.SwFace

        ElseIf Me.estCylindre Then
            ' on doit créer un plan de référence...

        ElseIf Me.estFauxPlan(inter.x, inter.y, inter.z) Then
            Dim vEdge As Object
            Dim i As Integer
            Dim swArete2() As SldWorks.Edge
            Dim swSommet As SldWorks.Vertex

            vEdge = Me.SwFace.GetEdges
            swArete2 = vEdge
            swModel.ClearSelection2(True)

            While planReference Is Nothing
                If UBound(swArete2) - 2 < i Then MsgBox("Dans CoupeLong, problème pour créer un plan avec 3 points.  (La face est un FauxPlan)", MsgBoxStyle.Critical, "Le plan ne sera pas créé") : Exit While
                swSommet = swArete2(i).GetStartVertex()
                swEnt = swSommet
                swEnt.Select4(False, Nothing)
                swArete2(i + 1).GetStartVertex()
                swEnt = swSommet
                swEnt.Select4(True, Nothing)
                swArete2(i + 2).GetStartVertex()
                swEnt = swSommet
                swEnt.Select4(True, Nothing)
                i += 1
                planReference = swModel.CreatePlaneThru3Points3(False)
                PlanEntity = planReference
            End While


        Else ' la face est une spline
            MsgBox("Dans coupeCoté, la face est un type de surface qui n'est pas encore traité")
        End If


        baseOriginal(0) = inter.x : baseOriginal(1) = inter.y : baseOriginal(2) = inter.z


        Dim Psi As Double
        Dim u(2) As Double, v(2) As Double, usketch(2) As Double, vsketch(2) As Double
        Dim Arete As SldWorks.Edge = Nothing
        Dim retval As Object
        u = poutre.GetOrientation(inter.x, inter.y, inter.z)


        vArete = Me.SwFace.GetEdges

        For Each Arete In vArete
            If Commun.Distance(Arete, inter.x, inter.y, inter.z) < Epsilon Then Exit For
        Next

        retval = Arete.GetClosestPointOn(inter.x, inter.y, inter.z)
        retval = Arete.Evaluate(retval(3))
        v(0) = retval(3) : v(1) = retval(4) : v(2) = retval(5)


        For g = 0 To 1

            PlanEntity.Select(False)
            swModel.InsertSketch2(True)
            swSketch = swModel.GetActiveSketch2

            pt3 = Commun.TransfertModelSketch(swSketch, pt3Original)
            usketch = Commun.TransfertModelSketch(swSketch, u)  ' on les met dans le plan du sketch
            vsketch = Commun.TransfertModelSketch(swSketch, v)
            base = Commun.TransfertModelSketch(swSketch, baseOriginal)
            Psi = Outils_Math.cosdir(usketch, vsketch)

            Dim a As Double, b As Double
            'longueur = Math.Sqrt(pt3(0) * pt3(0) + pt3(1) * pt3(1))
            'If pt3(1) = 0 Then a = 999999999999 Else a = Math.Abs(poutre.GetD2() * longueur / pt3(1))
            'If pt3(0) = 0 Then b = 999999999999 Else b = Math.Abs(poutre.GetD1() * longueur / pt3(0))
            ' À revoir. Si le plan est un cylindre ça marche plus. sans compter l'épaisseur de la poutre.
            ' pour l'instant je prend la plus prtite valeur...
            a = poutre.GetD1
            b = poutre.GetD2
            DemiLargeur = Math.Min(a, b)
            cut = DemiLargeur / Math.Sin(Pi / 2 - Psi)


            Dim P1(1) As Double
            Dim P2(1) As Double
            Dim P3(1) As Double
            Dim P4(1) As Double
            Dim Ptest(2) As Double


            If g = 0 Then
                P1(0) = -cut
                P1(1) = -cut '* mult ' 0
                P2(0) = 0
                P2(1) = -cut '* mult ' 0
                P3(0) = 0
                P3(1) = cut   'Intersections.Taille mult
                P4(0) = -cut
                P4(1) = cut   'Intersections.Taille mult
                sk(0) = -Epsilon * 100 + base(0) : sk(1) = 0 + base(1)

            Else
                P1(0) = 0
                P1(1) = -cut '* mult '0
                P2(0) = +cut
                P2(1) = -cut '* mult '0
                P3(0) = +cut
                P3(1) = cut  'Intersections.Taille mult
                P4(0) = 0
                P4(1) = cut  'Intersections.Taille mult
                sk(0) = Epsilon * 100 + base(0) : sk(1) = 0 + base(1)

            End If

            P1 = Outils_Math.Rotation2D(vsketch, P1)
            P2 = Outils_Math.Rotation2D(vsketch, P2)
            P3 = Outils_Math.Rotation2D(vsketch, P3)
            P4 = Outils_Math.Rotation2D(vsketch, P4)
            sk = Outils_Math.Rotation2D(vsketch, sk)

            sketchline = swModel.CreateLine2(P1(0) + base(0), P1(1) + base(1), 0, P2(0) + base(0), P2(1) + base(1), 0)
            sketchline = swModel.CreateLine2(P2(0) + base(0), P2(1) + base(1), 0, P3(0) + base(0), P3(1) + base(1), 0)
            sketchline = swModel.CreateLine2(P3(0) + base(0), P3(1) + base(1), 0, P4(0) + base(0), P4(1) + base(1), 0)
            sketchline = swModel.CreateLine2(P1(0) + base(0), P1(1) + base(1), 0, P4(0) + base(0), P4(1) + base(1), 0)

            swModel.InsertSketch2(True) ' l'option true rebuild le modèle avec les modifications du sketch.
            swModel.ClearSelection2(True)
            swEnt = Me.SwFace : swEnt.Select2(False, 1)
            swEnt = swSketch : swEnt.Select2(True, 4)

            swModel.InsertSplitLineProject(Directionnel, Flip)
            r = Commun.TransfertSketchToModel(swSketch, sk)
            Face(g) = UpdateApresSplit(inter, poutre, r(0), r(1), r(2), planReference)
            'If Face(g) Is Nothing Then
            'swSketch.Select(False)
            'swModel.EditDelete()
            'End If

        Next g


        ' mettre les mini-poutres
        Dim vEdge2 As Object
        Dim swArete As SldWorks.Edge
        Dim vPoint As Object
        Dim Mini1 As SldWorks.Edge = Nothing, Mini2 As SldWorks.Edge = Nothing


        ' 1 - trouver les 2 arrères dont l'orientation est la même (ou l'inverse) que le v
        For g = 0 To 1
            If Not Face(g) Is Nothing Then
                vEdge2 = Face(g).GetEdges()


                ' construire u
                For Each swArete In vEdge2
                    If Commun.Distance(swArete, inter.x, inter.y, inter.z) < Epsilon Then
                        ' l'arête touche à l'intersection, 
                        vPoint = swArete.GetClosestPointOn(inter.x, inter.y, inter.z)
                        vPoint = swArete.Evaluate(vPoint(3))
                        u(0) = vPoint(3) : u(1) = vPoint(4) : u(2) = vPoint(5)

                        If Outils_Math.CompareSens(v, u) Then
                            ' l'arète doit être une mini-poutre
                            If Mini1 Is Nothing Then Mini1 = swArete : Exit For Else Mini2 = swArete : Exit For
                        End If
                    End If

                Next

            End If
        Next

        swEnt = Mini1
        If Not RealConstant.RCCode.CreationAttributPourPoutre(swEnt, "Mini", True) Then MsgBox("A pas marché!")

        If Not Mini2 Is Nothing Then
            swEnt = Mini2
            If Not RealConstant.RCCode.CreationAttributPourPoutre(swEnt, "Mini", True) Then MsgBox("A pas marché!")
        End If

        swModel.SetInferenceMode(True) '
        'swModel.SetAddToDB(False)
        'swModel.SetDisplayWhenAdded(True) ' 
    End Sub


    ' sub qui coupe la face avec une arète qui repose dessus.
    Friend Sub CoupeLong(ByRef inter As InterPoutreVolume, ByVal poutre As SlyAretePoutre)
        Dim swEnt As SldWorks.Entity
        Dim swSketchSegment As SldWorks.SketchSegment
        Dim vSketchSegments As Object
        Dim swSketch As SldWorks.Sketch
        Dim faceinterne(1) As SldWorks.Face2
        Dim swPlan As SldWorks.RefPlane = Nothing
        Dim b As Integer
        Dim swSommet As SldWorks.Vertex
        Dim i As Integer

        swModel.SetInferenceMode(False) ' ne pas mettre de contraintes par défaut

        ' faut découper toutes les faces de la liste si elles ne sont pas des faces internes
        Dim MeFace As SldWorks.Face2
        Dim ListeFace() As SldWorks.Face2
        ReDim ListeFace(Me.lst_Faces.Count - 1)

        For i = 1 To Me.lst_Faces.Count
            ListeFace(i - 1) = Me.lst_Faces.Item(i)
        Next

        For Each MeFace In ListeFace

            If Me.estPlan Then
                swEnt = MeFace
                swEnt.Select(False)
                swPlan = swModel.CreatePlaneAtOffset3(0, False, False)
                swEnt.Select(False)
                swModel.InsertSketch2(True)

                swPlan = swModel.CreatePlaneAtOffset3(0, False, False)

            ElseIf Me.estFauxPlan(inter.x, inter.y, inter.z) Then
                If b = 0 Then
                    Dim vEdge As Object


                    vEdge = MeFace.GetEdges
                    swModel.ClearSelection2(True)
                    While swPlan Is Nothing
                        If UBound(vEdge) - 2 < i Then MsgBox("Dans CoupeLong, problème pour créer un plan avec 3 points.  (La face est un FauxPlan)", MsgBoxStyle.Critical, "Le plan ne sera pas créé") : Exit While
                        swSommet = vEdge(i).GetStartVertex()
                        swEnt = swSommet
                        swEnt.Select4(False, Nothing)
                        swSommet = vEdge(i + 1).GetStartVertex()
                        swEnt = swSommet
                        swEnt.Select4(True, Nothing)
                        swSommet = vEdge(i + 2).GetStartVertex()
                        swEnt = swSommet
                        swEnt.Select4(True, Nothing)
                        i += 1
                        swPlan = swModel.CreatePlaneThru3Points3(False)

                    End While
                End If
                swEnt = swPlan
                swEnt.Select(False)
                swModel.InsertSketch2(True)

            Else
                MsgBox("Dans coupeLong, on a un type de face qui n'est pas encore traité")

            End If


            swSketch = swModel.GetActiveSketch2

            swEnt = poutre.swArete
            swEnt.Select(False)

            ' créer la ligne de  «conversion de entités»
            swModel.SketchUseEdge2(False)

            vSketchSegments = swSketch.GetSketchSegments()
            swSketchSegment = vSketchSegments(0)
            swSketchSegment.Select2(False, 1)  'on sélectionne l'arète de poutre...

            Dim x As Double, y As Double, z As Double
            Commun.GetMidPointSegment(swSketchSegment, x, y, z)


            ' sketchoffset doit avoir un mark de 1 pour l'objet à offsetter.  Une valeur négative inverse la direction
            swModel.SketchManager.SketchOffset(poutre.GetD2, False, 0, 0, 0, 0)
            ' pour rendre le modèle plus beau, on peut enlever la contrainte de offset et laisser solidworks mettre des contraintes automatiques...

            Dim retval As Object
            Dim skPointA1 As SldWorks.SketchPoint = Nothing, skPointA2 As SldWorks.SketchPoint = Nothing, skPointB1 As SldWorks.SketchPoint = Nothing, skPointB2 As SldWorks.SketchPoint = Nothing

            vSketchSegments = swSketch.GetSketchSegments()
            swSketchSegment = vSketchSegments(1)


            Select Case swSketchSegment.GetType()
                Case 0 ' on a une ligne
                    Dim sketchline As SldWorks.SketchLine
                    sketchline = swSketchSegment
                    skPointA1 = sketchline.GetStartPoint2
                    skPointA2 = sketchline.GetEndPoint2()
                Case 1 ' arc
                    Dim arc As SldWorks.SketchArc
                    arc = swSketchSegment
                    skPointA1 = arc.GetStartPoint
                    skPointA2 = arc.GetEndPoint2
                Case 2 ' ellipse
                    Dim sketchEllipse As SldWorks.SketchEllipse
                    sketchEllipse = swSketchSegment
                    skPointA1 = sketchEllipse.GetStartPoint2
                    skPointA2 = sketchEllipse.GetEndPoint2
                Case 3 ' spline
                    Dim spline As SldWorks.SketchSpline
                    Dim pts() As SldWorks.SketchPoint
                    spline = swSketchSegment
                    retval = spline.GetPoints2()
                    pts = retval
                    skPointA1 = pts(0)
                    skPointA2 = pts(UBound(pts))
                Case 5 ' parabole (le 4 est du texte)
                    Dim para As SldWorks.SketchParabola
                    para = swSketchSegment
                    skPointA1 = para.GetStartPoint2
                    skPointA2 = para.GetEndPoint2
            End Select

            swSketchSegment = vSketchSegments(0)
            Select Case swSketchSegment.GetType()
                Case 0 ' on a une ligne
                    Dim sketchline As SldWorks.SketchLine
                    sketchline = swSketchSegment
                    skPointB1 = sketchline.GetStartPoint2
                    skPointB2 = sketchline.GetEndPoint2()
                Case 1 ' arc
                    Dim arc As SldWorks.SketchArc
                    arc = swSketchSegment
                    skPointB1 = arc.GetStartPoint
                    skPointB2 = arc.GetEndPoint2
                Case 2 ' ellipse
                    Dim sketchEllipse As SldWorks.SketchEllipse
                    sketchEllipse = swSketchSegment
                    skPointB1 = sketchEllipse.GetStartPoint2
                    skPointB2 = sketchEllipse.GetEndPoint2
                Case 3 ' spline
                    Dim spline As SldWorks.SketchSpline
                    Dim pts() As SldWorks.SketchPoint
                    spline = swSketchSegment
                    retval = spline.GetPoints2()
                    pts = retval
                    skPointB1 = pts(0)
                    skPointB2 = pts(UBound(pts))
                Case 5 ' parabole (le 4 est du texte)
                    Dim para As SldWorks.SketchParabola
                    para = swSketchSegment
                    skPointB1 = para.GetStartPoint2
                    skPointB2 = para.GetEndPoint2
            End Select

            ' création des 2 lignes pour fermer le sketch.
            swModel.CreateLine2(skPointA1.X, skPointA1.Y, 0, skPointB1.X, skPointB1.Y, 0)
            swModel.CreateLine2(skPointA2.X, skPointA2.Y, 0, skPointB2.X, skPointB2.Y, 0)


            Dim x2 As Double, y2 As Double, z2 As Double ' le midpoint de la poutre
            Dim x3 As Double, y3 As Double, z3 As Double ' le midpoint de la poutre

            swSketchSegment = vSketchSegments(0) ' le midpoint d'une poutre
            Commun.GetMidPointSegment(swSketchSegment, x2, y2, z2)

            swSketchSegment = vSketchSegments(1) ' le midpoint de l'autre poutre
            Commun.GetMidPointSegment(swSketchSegment, x3, y3, z3)

            Dim sk(1) As Double, r(2) As Double
            sk(0) = (x3 + x2) / 2
            sk(1) = (y3 + y2) / 2
            r = Commun.TransfertSketchToModel(swSketch, sk)

            swModel.InsertSketch2(True) ' l'option true rebuild le modèle avec les modifications du sketch.
            swModel.ClearSelection2(True)
            swEnt = MeFace : swEnt.Select2(False, 1)
            swEnt = swSketch : swEnt.Select2(True, 4)

            swModel.InsertSplitLineProject(False, False)

            Me.Flag = 20 ' pour dire que l'on a un coupeLong
            faceinterne(b) = UpdateApresSplit(inter, poutre, r(0), r(1), r(2), swPlan) ' et ça s'occupe de créer la coque... mais je suis pas certain que c'est nécessaire
            Me.Flag = 0

            'If faceinterne(b) Is Nothing Then
            'swEnt = swSketch
            'swEnt.Select(False)
            'swModel.EditDelete()
            'End If


            ' reste à updater, on doit ajouter de 2 à 4 mini-poutres
            'Dim vEdges As Object
            'Dim Arete As SldWorks.Edge
            'Dim vFaces As Object
            'Dim aretePoutre As SldWorks.Edge


            ''For b = 0 To 1
            ''If Not faceinterne(b) Is Nothing Then
            'b = 0
            'While faceinterne(b) Is Nothing
            '    b += 1
            'End While

            'vEdges = faceinterne(b).GetEdges
            'For Each Arete In vEdges
            '    If Distance(Arete.GetStartVertex, inter.x, inter.y, inter.z) < Epsilon Then swSommet = Arete.GetStartVertex : Exit For
            '    If Distance(Arete.GetEndVertex, inter.x, inter.y, inter.z) < Epsilon Then swSommet = Arete.GetEndVertex : Exit For
            'Next

            'Dim T As Double
            'Dim xyz(2) As Double
            'Dim g As Integer
            'Dim courbe As SldWorks.Curve

            'T = poutre.GetT(inter.x, inter.y, inter.z)

            'vEdges = swSommet.GetEdges

            'T -= 10000 * Epsilon

            'For g = 0 To UBound(vEdges) ' boucle pas optimisée en vitesse
            '    Arete = vEdges(g)
            '    If poutre.Evaluer(T, xyz) Then
            '        courbe = Arete.GetCurve
            '        If Distance(courbe, xyz(0), xyz(1), xyz(2)) < Epsilon Then aretePoutre = Arete
            '    ElseIf poutre.Evaluer(T + 20000 * Epsilon, xyz) Then
            '        courbe = Arete.GetCurve
            '        If Distance(courbe, xyz(0), xyz(1), xyz(2)) < Epsilon Then aretePoutre = Arete ' la distance devrait être entre la droite (et non l'arête) et le point.
            '    End If
            'Next g


            'If aretePoutre Is Nothing Then
            '    ' putain d'enfoiré de merde!!!! on trouve pas la courbe de la poutre. alors on sort avec dignité!
            '    Exit Sub
            '    ' anyway on avait dit en réunion de ne pas mettre de minipoutres...
            'End If


            '' on a l'arète de la poutre, avec la boucle de la face on prend une arète avant et une arête après.
            'Dim swBoucle As SldWorks.Loop2
            'Dim AreteAvant As SldWorks.Edge, AreteSuivant As SldWorks.Edge, areteTest As SldWorks.Edge
            'Dim varete As Object
            'Dim Mini1 As SldWorks.Edge, Mini2 As SldWorks.Edge, Mini3 As SldWorks.Edge, Mini4 As SldWorks.Edge

            'Dim j As Integer

            'If Not faceinterne(0) Is Nothing Then
            '    swBoucle = faceinterne(0).GetFirstLoop ' devrait y en avoir juste une...
            '    varete = swBoucle.GetEdges()

            '    For j = 0 To UBound(varete)
            '        areteTest = varete(j)
            '        If areteTest Is aretePoutre Then
            '            If j <> 0 Then Mini1 = varete(j - 1) Else Mini1 = varete(UBound(varete))
            '            If j <> UBound(varete) Then Mini2 = varete(j + 1) Else Mini2 = varete(0)
            '        End If
            '    Next j
            'End If

            'If Not faceinterne(1) Is Nothing Then
            '    swBoucle = faceinterne(1).GetFirstLoop ' devrait y en avoir juste une...
            '    varete = swBoucle.GetEdges()

            '    For j = 0 To UBound(varete)
            '        areteTest = varete(j)
            '        If areteTest Is aretePoutre Then
            '            If j <> 0 Then Mini3 = varete(j - 1) Else Mini3 = varete(UBound(varete))
            '            If j <> UBound(varete) Then Mini4 = varete(j + 1) Else Mini4 = varete(0)
            '        End If
            '    Next j
            'End If
            'If Not Mini1 Is Nothing Then
            '    swEnt = Mini1
            '    If Not RealConstant.RCCode.CreationAttributPourPoutre(swEnt, "Mini", True) Then MsgBox("A pas marché!")
            'End If

            'If Not Mini2 Is Nothing Then
            '    swEnt = Mini2
            '    If Not RealConstant.RCCode.CreationAttributPourPoutre(swEnt, "Mini", True) Then MsgBox("A pas marché!")
            'End If

            'If Not Mini3 Is Nothing Then
            '    swEnt = Mini3
            '    If Not RealConstant.RCCode.CreationAttributPourPoutre(swEnt, "Mini", True) Then MsgBox("A pas marché!")
            'End If

            'If Not Mini4 Is Nothing Then
            '    swEnt = Mini4
            '    If Not RealConstant.RCCode.CreationAttributPourPoutre(swEnt, "Mini", True) Then MsgBox("A pas marché!")
            'End If
        Next MeFace


        swModel.SetInferenceMode(True) ' ne pas mettre de contraintes par défaut
        'swModel.SetAddToDB(False)
        'swModel.SetDisplayWhenAdded(True) ' accélérer les performances

    End Sub


    ' sub qui coupe la face normalement, avec un X.... [cas #1]
    Friend Sub CoupeX(ByRef inter As InterPoutreVolume, ByRef poutre As SlyAretePoutre)

        Dim swEnt As SldWorks.Entity = Nothing
        Dim Directionnel As Boolean, Flip As Boolean
        Dim Faces(3) As SldWorks.Face2
        Dim r(2) As Double
        Dim LaSurface As SldWorks.Surface
        Dim sens As Boolean
        Dim p(2) As Double
        Dim retour() As Double

        swModel.SetInferenceMode(False) ' ne pas mettre de contraintes par défaut
        'swModel.SetAddToDB(True)
        'swModel.SetDisplayWhenAdded(False) ' accélérer les performances


        ' l'idée est de sélectionner le point et l'arète puis d'utiliser CreatePlanePerCurveAndPassPoint3
        Dim planReference As SldWorks.RefPlane
        Dim swsketch As SldWorks.Sketch
        Dim swSommet As SldWorks.Vertex, swSommet2 As SldWorks.Vertex
        Dim pointdeb(2) As Double, pointfin(2) As Double

        'swModel.Extension.SelectByID2("", "POINTREF", inter.x, inter.y, inter.z, False, 0, Nothing, 0)
        ' faut vraiment sélectionner le bon point...
        swSommet = poutre.swArete.GetStartVertex()
        swSommet2 = poutre.swArete.GetEndVertex()
        If swSommet Is Nothing Then
            MsgBox("On a un cercle ou courbe sans sommets, dans coupeX, pas encore traité.  Ne peut pas mettre un plan si pas de sommet")
        Else
            If Distance(swSommet, inter.x, inter.y, inter.z) < Epsilon Then
                swEnt = swSommet
            ElseIf Distance(swSommet2, inter.x, inter.y, inter.z) < Epsilon Then
                swEnt = swSommet2
            Else
                MsgBox("Dans coupeX, l'intersection n'est pas sur un sommet.  Pas encore traité.  Nécessite de créer un point au coordonnées d'intersection")
            End If
        End If

        swEnt.Select4(False, Nothing)
        swEnt = poutre.swArete
        swEnt.Select(True)

        If Me.estPlan Or Me.estFauxPlan(inter.x, inter.y, inter.z) Then
            ' si la coque est plane alors on projette le plan de référence des deux cotés, sinon on doit le décaler vers le bas
            planReference = swModel.CreatePlanePerCurveAndPassPoint3(True, True)
            Directionnel = False
            Flip = False
        ElseIf Me.estCylindre Then
            ' on a un cylindre, on ne projette pas des 2 cotés.  On créé un plan, puis un autre plus bas pour ensuite projeter d'un seul coté.
            Dim PlanDessus As SldWorks.RefPlane
            Dim Rayon As Double, L As Double, B As Double, phi As Double, dist As Double, temp1 As Double, temp2 As Double
            Dim u(2) As Double, v(2) As Double
            PlanDessus = swModel.CreatePlanePerCurveAndPassPoint3(True, True)
            temp1 = poutre.GetD1
            temp2 = poutre.GetD2
            L = Math.Sqrt(temp1 * temp1 + temp2 * temp2)
            Rayon = Me.GetRayonCylindre()
            u = poutre.GetOrientation(inter.x, inter.y, inter.z)
            v = Me.GetNormale(inter.x, inter.y, inter.z)
            phi = -(Math.Acos(Outils_Math.cosdir(u, v)))
            B = Math.Abs(L / 2 * Math.Sin(phi))
            dist = Rayon - Math.Sqrt(Rayon * Rayon - ((L / 2) * (L / 2))) + B
            If dist < 0 Then MsgBox("Gros problème pour couper le cylindre, la poutre est plus grosse!!!!!!", MsgBoxStyle.Critical) : Exit Sub

            swEnt = PlanDessus
            swEnt.Select(False)
            Directionnel = True

            Flip = Flipper(PlanDessus, inter)

            planReference = swModel.CreatePlaneAtOffset3(dist * 2, Flip, True)
        Else
            MsgBox("La coque n'est ni un cylindre, ni un plan" & vbCr & "Le résultat n'est pas certain...", MsgBoxStyle.Information, "Avertissement")
            planReference = swModel.CreatePlanePerCurveAndPassPoint3(True, True)
            Directionnel = False
            Flip = False
        End If


        LaSurface = Me.SwFace.GetSurface()
        sens = Me.SwFace.FaceInSurfaceSense()

        ' skx est la coordonnée du point de ref en coord de sketch, Rx est le point de référence dans le repère global.
        Dim i As Integer, MettreFI As Boolean
        Dim swFeat As SldWorks.Feature

        For i = 0 To 1

            swEnt = planReference
            swEnt.Select(False)
            swModel.InsertSketch2(False)
            swModel.ClearSelection2(True)
            swFeat = swModel.FeatureByPositionReverse(0)
            swModel.SelectByID(swFeat.Name, "SKETCH", 0, 0, 0)
            swModel.EditSketch()
            swsketch = swModel.GetActiveSketch2

            p(0) = inter.x : p(1) = inter.y : p(2) = inter.z
            retour = Commun.TransfertModelSketch(swsketch, p)


            r = DessineSectionPoutre(poutre, retour(0), retour(1), i + 1, swsketch, inter, MettreFI)
            swModel.InsertSketch2(True) ' l'option true rebuild le modèle avec les modifications du sketch.
            swModel.ClearSelection2(True)

            Dim face As SldWorks.Face2
            For Each face In Me.lst_Faces
                swModel.ClearSelection2(True)
                swEnt = face : swEnt.Select2(False, 1)
                swEnt = swsketch : swEnt.Select2(True, 4)
                swModel.InsertSplitLineProject(Directionnel, Flip)
            Next


            Me.SwFace.DetachSurface()
            Me.SwFace.AttachSurface(LaSurface, sens)

            Faces(i) = UpdateApresSplit(inter, poutre, r(0), r(1), r(2), planReference, MettreFI)
            Commun.MettreUnPoint(r(0), r(1), r(2))

            If Faces(i) Is Nothing Then
                swEnt.Select(False)
                swModel.EditDelete()
            End If
            If Flag = 2 Then Flag = 0 : Exit For

        Next i

        swModel.SetInferenceMode(True)
        'swModel.SetAddToDB(False)
        'swModel.SetDisplayWhenAdded(True)

    End Sub


    ' Sub qui appelle le découpage de la face
    Public Sub decouper()

        If lst_InterPoutre.Count = 0 Then Exit Sub ' sortir si on a pas d'intersection


        ' les attributs ne sont pas updatés sur les faces (mais sur les arètes et les sommets c'est OK)
        ' on mémorise l'attribut de la face et on la réapplique à la fin.


        Dim i As Integer
        Dim inter As InterPoutreVolume
        Dim nb1 As Integer, nb2 As Integer, nb3 As Integer, nb5 As Integer
        Dim poutre1 As SlyAretePoutre = Nothing, poutre3 As SlyAretePoutre = Nothing
        Dim lst_poutre2 As New Collection
        Dim aire As Double
        Dim poutreTest As SlyAretePoutre

        Dim lst_coupeXinter As New Collection
        Dim lst_coupeXPoutre As New Collection
        Dim lst_coupeLinter As New Collection
        Dim lst_coupeLPoutre As New Collection
        Dim lst_coupeCinter As New Collection
        Dim lst_coupeCPoutre As New Collection


        For Each inter In lst_InterPoutre

            'pour chaque intersection on peut avoir plusieurs poutres...
            For i = 1 To inter.lst_sPoutre.Count
                poutreTest = inter.lst_sPoutre.Item(i)
                Select Case CInt(inter.lst_type.Item(i))
                    Case 1
                        If poutreTest.GetAireCarree > aire Then poutre1 = poutreTest
                        nb1 += 1
                    Case 2
                        lst_poutre2.Add(poutreTest)
                        nb2 += 1
                    Case 3
                        If poutreTest.GetAireCarree > aire Then poutre3 = poutreTest
                        nb3 += 1
                    Case 5 ' un poutre à faceDeSection
                        nb5 += 1
                    Case 22
                        ' on fait rien, mais c'est pour éviter le msgbox du case else...
                    Case Else
                        MsgBox("Problème dans découper de SlyFaceCoque, le type d'intersection n'est pas reconnu", MsgBoxStyle.Critical)
                End Select
            Next i


            If nb1 > 0 Then 'CoupeX(inter, poutre1) ' on coupe le x en premier
                lst_coupeXinter.Add(inter)
                lst_coupeXPoutre.Add(poutre1)
            End If


            For Each poutreTest In lst_poutre2 ' puis on coupe sur la longueur 'CoupeLong(inter, poutreTest)
                lst_coupeLinter.Add(inter)
                lst_coupeLPoutre.Add(poutreTest)
            Next

            If nb3 > 0 Then 'CoupeCote(inter, poutre3) ' finalement on coupe sur les cotés
                lst_coupeCinter.Add(inter)
                lst_coupeCPoutre.Add(poutre3)
            End If

            If nb5 = 1 And (nb1 > 0 Or nb2 > 0 Or nb3 > 0) Then
                MsgBox("Problème, on a un type d'intersection impossible dans la vraie vie!", MsgBoxStyle.Exclamation, "Design impossible à obtenir en réalité...")
            End If


            lst_poutre2.Clear()
            nb1 = 0 : nb2 = 0 : nb3 = 0


        Next inter


        ' maintenant on a toutes les lists d'intersections.  On les coupe.
        For i = 1 To lst_coupeXinter.Count
            CoupeX(lst_coupeXinter.Item(i), lst_coupeXPoutre.Item(i))
        Next

        For i = 1 To lst_coupeLinter.Count
            CoupeLong(lst_coupeLinter.Item(i), lst_coupeLPoutre.Item(i))
        Next

        For i = 1 To lst_coupeCinter.Count
            CoupeCote(lst_coupeCinter.Item(i), lst_coupeCPoutre.Item(i))
        Next
        If nb5 = 1 Then
            If lst_InterPoutre.Count <> 1 Then MsgBox("Plus d'une intersection du type FacedeSection....")
            CoupeFaceDeSection(lst_InterPoutre(1))
        End If

    End Sub

    Private Sub CoupeFaceDeSection(ByRef inter As InterPoutreVolume)
        Dim swEnt As SldWorks.Entity = Nothing
        Dim Directionnel As Boolean
        Dim Faces(3) As SldWorks.Face2
        Dim r(2) As Double
        Dim p(2) As Double
        Dim planReference As SldWorks.RefPlane = Nothing
        Dim swsketch As SldWorks.Sketch
        Dim pointdeb(2) As Double, pointfin(2) As Double
        Dim sketchline As SldWorks.SketchLine
        Dim swFeat As SldWorks.Feature


        swModel.SetInferenceMode(False) ' ne pas mettre de contraintes par défaut

        swEnt = Me.SwFace
        swEnt.Select(False)
        swModel.InsertSketch2(False)
        swsketch = swModel.GetActiveSketch2

        ' dessin de la forme à faire...

        Dim xyzc() As Double, xyz(2) As Double
        xyz(0) = inter.x : xyz(1) = inter.y : xyz(2) = inter.z
        xyzc = Commun.TransfertModelSketch(swsketch, xyz)

        sketchline = swModel.CreateLine2(xyzc(0), xyzc(1), 0, xyzc(0) + Math.Cos(Pi / 4), xyzc(1) + Math.Sin(Pi / 4), 0)
        sketchline = swModel.CreateLine2(xyzc(0), xyzc(1), 0, xyzc(0) - Math.Cos(Pi / 4), xyzc(1) + Math.Sin(Pi / 4), 0)
        swModel.CreateArc2(xyzc(0), xyzc(1), 0, xyzc(0) + Math.Cos(Pi / 4), xyzc(1) + Math.Sin(Pi / 4), 0, xyzc(0) - Math.Cos(Pi / 4), xyzc(1) + Math.Sin(Pi / 4), 0, 1) ' le dernier param est la direction.  1 ou -1

        swModel.InsertSketch2(True) ' l'option true rebuild le modèle avec les modifications du sketch.
        swModel.ClearSelection2(True)

        swEnt = Me.SwFace : swEnt.Select2(False, 1)
        swEnt = swsketch : swEnt.Select2(True, 4)
        swModel.InsertSplitLineProject(Directionnel, False)


        ' flagger les 2 faces comme faces Internes.
        Dim vface As Object
        Dim face As SldWorks.Face2
        Dim attr As SldWorks.Attribute
        swFeat = swModel.FeatureByPositionReverse(0)
        Try
            vface = swFeat.GetFaces
            For Each face In vface
                '**************
                Dim nom2 As String = "FaceInterne" & no
                swEnt = face
                attr = swEnt.FindAttribute(Intersections.DefAttrFaceInterne, 0) ' si l'attribut existe déjà on pointe dessus.
                If attr Is Nothing Then attr = Intersections.DefAttrFaceInterne.CreateInstance5(swModel, face, nom2, 0, 2) ' 0 = swThisconfig
                While attr Is Nothing
                    no += 1
                    nom2 = "FaceInterne" & CStr(no)
                    attr = Intersections.DefAttrFaceInterne.CreateInstance5(swModel, face, nom2, 0, 2)
                End While
                GererDossiers("FaceInternes", nom2)
                no += 1


                '**************
                Me.AjouterFace(face)
            Next face
        Catch
            ' si on a une poutre dessinée en partie, on aura pas de feature... mais on a la face...
            ' on doit donc le déterminer anyway
        End Try

        ' si ça ne touche pas à la face
        If Not Distance(Me.SwFace, inter.x, inter.y, inter.z) < Epsilon Then
            AjouterMiniPoutresSurFaceInterne(inter.lst_sPoutre.Item(1), face, inter.x, inter.y, inter.z)
            'MsgBox("On ajoute une mini-poutre entre la poutre " & inter.lst_sPoutre.Item(1).nom & vbCr & " et le point ( " & Format(inter.x * 1000, "0000") & " , " & Format(inter.y * 1000, "0000") & " , " & Format(inter.y * 1000, "0000") & " )")
        End If

        swModel.SetInferenceMode(True)

    End Sub


    Friend Overridable Sub chercherAttributs()
        Dim swEnt As SldWorks.Entity
        Dim attr As SldWorks.Attribute

        swEnt = Me.SwFace

        attr = swEnt.FindAttribute(Intersections.DefAttrConditionLimite, 0)
        If Not attr Is Nothing Then Me.AttributCL = attr : attr = Nothing

        'attr = swent.findattribute(intersections.DefAttrFaceInterne,0) ' ne devrait pas s'entrecouper...

    End Sub

    Private Function Flipper(ByRef PlanDessus As SldWorks.RefPlane, ByRef inter As InterPoutreVolume) As Boolean
        ' function qui dit si l'on doit flipper le sens du plan de référence.
        ' calcul de la direction à prendre
        Dim retval As Object
        Dim ret(8) As Double
        Dim ret2(6) As Double
        Dim normalePlan(2) As Double
        Dim OV(2) As Double
        Dim swSurf As SldWorks.Surface

        retval = PlanDessus.GetRefPlaneParams()
        ret = retval
        normalePlan(0) = ret(6) : normalePlan(1) = ret(7) : normalePlan(2) = ret(8)
        swSurf = Me.lst_Faces.Item(1).GetSurface
        retval = swSurf.CylinderParams() ' 7 doubles, les 3 premiers sont l'origine
        ret2 = retval
        OV(0) = ret2(0) - inter.x : OV(1) = ret2(1) - inter.y : OV(2) = ret2(2) - inter.z

        ' l'angle est le produit scalaire divisé par la norme des 2 vecteurs
        Dim temp As Double = Outils_Math.Angle2Vecteurs(OV, normalePlan)
        If (temp < Pi / 2) And (temp > -Pi / 2) Then Return False Else Return True

    End Function


    Private Function DessineSectionPoutre(ByRef Poutre As SlyAretePoutre, ByVal TranslationX As Double, ByVal TranslationY As Double, ByVal numero As Integer, ByRef swSketch As SldWorks.Sketch, ByRef inter As InterPoutreVolume, ByRef MettreFI As Boolean) As Double()
        ' le sketch est déjà inséré, il faut juste mettre des swmodel.line ou autre
        ' doit retourner r() qui est un point situé à l'intérieur de la coupe
        Dim sketchline As SldWorks.SketchSegment
        Dim longueur As Double
        Dim r(2) As Double
        Dim sk(1) As Double
        Dim i As Integer
        Dim Ligne() As Double = Nothing ' liste des lignes (4 valeurs par ligne)
        Dim pt3() As Double
        Dim Nomsection As String

        MettreFI = True

        ' on doit activer le sketch avant d'utiliser la fonction getactivesketch
        pt3 = Poutre.GetPoint3
        longueur = Math.Sqrt(pt3(0) * pt3(0) + pt3(1) * pt3(1))
        Dim IP(2) As Double ' IP est le vecteur directionnel
        IP(0) = pt3(0) - inter.x : IP(1) = pt3(1) - inter.y : IP(2) = pt3(2) - inter.z


        pt3 = Commun.TransfertModelSketch(swSketch, pt3)
        longueur = Math.Sqrt(pt3(0) * pt3(0) + pt3(1) * pt3(1))

        Nomsection = Poutre.GetNomSection
        If Nomsection = "Rectangle" Or Nomsection = " Rectangle générique" Then ' un rectangle
            Select Case numero
                Case 1
                    Dim P(2, 1) As Double
                    P(0, 0) = 0
                    P(0, 1) = 0
                    P(1, 0) = Poutre.GetD1 / 2
                    P(1, 1) = -Poutre.GetD2 / 2
                    P(2, 0) = Poutre.GetD1 / 2
                    P(2, 1) = Poutre.GetD2 / 2

                    ReDim Ligne(11)
                    pt3(0) -= TranslationX
                    pt3(1) -= TranslationY
                    For i = 0 To 2
                        Outils_Math.Rotation2D(pt3, P(i, 0), P(i, 1))
                        P(i, 0) += TranslationX
                        P(i, 1) += TranslationY
                    Next i

                    For i = 0 To 1
                        Ligne(i * 4) = P(i, 0) : Ligne(i * 4 + 1) = P(i, 1) : Ligne(i * 4 + 2) = P(i + 1, 0) : Ligne(i * 4 + 3) = P(i + 1, 1)
                    Next i
                    Ligne(8) = P(2, 0) : Ligne(9) = P(2, 1) : Ligne(10) = P(0, 0) : Ligne(11) = P(0, 1)

                    r(0) = inter.x + 5000 * Epsilon * IP(0)
                    r(1) = inter.y + 5000 * Epsilon * IP(1)
                    r(2) = inter.z + 5000 * Epsilon * IP(2)

                Case 2
                    ReDim Ligne(19)

                    Dim p(4, 1) As Double
                    p(0, 0) = 0
                    p(0, 1) = 0
                    p(1, 0) = Poutre.GetD1 / 2
                    p(1, 1) = Poutre.GetD2 / 2
                    p(2, 0) = -Poutre.GetD1 / 2
                    p(2, 1) = Poutre.GetD2 / 2
                    p(3, 0) = -Poutre.GetD1 / 2
                    p(3, 1) = -Poutre.GetD2 / 2
                    p(4, 0) = Poutre.GetD1 / 2
                    p(4, 1) = -Poutre.GetD2 / 2


                    pt3(0) -= TranslationX
                    pt3(1) -= TranslationY
                    pt3(0) /= longueur : pt3(1) /= longueur
                    For i = 0 To 4
                        Outils_Math.Rotation2D(pt3, p(i, 0), p(i, 1))
                        p(i, 0) += TranslationX
                        p(i, 1) += TranslationY
                    Next i

                    For i = 0 To 3
                        Ligne(i * 4) = p(i, 0) : Ligne(i * 4 + 1) = p(i, 1) : Ligne(i * 4 + 2) = p(i + 1, 0) : Ligne(i * 4 + 3) = p(i + 1, 1)
                    Next i
                    Ligne(16) = p(4, 0) : Ligne(17) = p(4, 1) : Ligne(18) = p(0, 0) : Ligne(19) = p(0, 1)
                    r(0) = inter.x - 5000 * Epsilon * IP(0)
                    r(1) = inter.y - 5000 * Epsilon * IP(1)
                    r(2) = inter.z - 5000 * Epsilon * IP(2)

            End Select
            MettreFI = True
        ElseIf Left(Nomsection, 2) = "ST" Or Nomsection = " Tube carré générique" Then ' tube carré troué
            Dim P(3, 1) As Double
            Select Case numero
                Case 1
                    P(0, 0) = Poutre.GetD1 / 2
                    P(0, 1) = -Poutre.GetD2 / 2
                    P(1, 0) = P(0, 0)
                    P(1, 1) = -P(0, 1)
                    P(2, 0) = -P(0, 0)
                    P(2, 1) = P(1, 1)
                    P(3, 0) = P(2, 0)
                    P(3, 1) = P(0, 1)

                    r(0) = P(0, 0) - 1000 * Epsilon
                    r(1) = 0 : r(2) = 0
                    Outils_Math.Rotation2D(pt3, r(0), r(1))
                    r(0) += TranslationX
                    r(1) += TranslationY
                    r = Commun.TransfertSketchToModel(swSketch, r)

                    pt3(0) -= TranslationX
                    pt3(1) -= TranslationY
                    pt3(0) /= longueur : pt3(1) /= longueur
                    For i = 0 To 3
                        Outils_Math.Rotation2D(pt3, P(i, 0), P(i, 1))
                        P(i, 0) += TranslationX
                        P(i, 1) += TranslationY
                    Next i

                    ReDim Ligne(15)
                    Ligne(0) = P(0, 0) : Ligne(1) = P(0, 1) : Ligne(2) = P(1, 0) : Ligne(3) = P(1, 1)
                    Ligne(4) = P(1, 0) : Ligne(5) = P(1, 1) : Ligne(6) = P(2, 0) : Ligne(7) = P(2, 1)
                    Ligne(8) = P(2, 0) : Ligne(9) = P(2, 1) : Ligne(10) = P(3, 0) : Ligne(11) = P(3, 1)
                    Ligne(12) = P(3, 0) : Ligne(13) = P(3, 1) : Ligne(14) = P(0, 0) : Ligne(15) = P(0, 1)
                    MettreFI = False

                Case 2
                    P(0, 0) = Poutre.GetD1 / 2 - Poutre.GetD3
                    P(0, 1) = -Poutre.GetD2 / 2 + Poutre.GetD3
                    P(1, 0) = P(0, 0)
                    P(1, 1) = -P(0, 1)
                    P(2, 0) = -P(1, 0)
                    P(2, 1) = P(1, 1)
                    P(3, 0) = P(2, 0)
                    P(3, 1) = P(0, 1)

                    r(0) = P(0, 0) + 1000 * Epsilon
                    r(1) = 0 : r(2) = 0
                    Outils_Math.Rotation2D(pt3, r(0), r(1))
                    r(0) += TranslationX
                    r(1) += TranslationY
                    r = Commun.TransfertSketchToModel(swSketch, r)

                    pt3(0) -= TranslationX
                    pt3(1) -= TranslationY
                    pt3(0) /= longueur : pt3(1) /= longueur
                    For i = 0 To 3
                        Outils_Math.Rotation2D(pt3, P(i, 0), P(i, 1))
                        P(i, 0) += TranslationX
                        P(i, 1) += TranslationY
                    Next i

                    ReDim Ligne(15)
                    Ligne(0) = P(0, 0) : Ligne(1) = P(0, 1) : Ligne(2) = P(1, 0) : Ligne(3) = P(1, 1)
                    Ligne(4) = P(1, 0) : Ligne(5) = P(1, 1) : Ligne(6) = P(2, 0) : Ligne(7) = P(2, 1)
                    Ligne(8) = P(2, 0) : Ligne(9) = P(2, 1) : Ligne(10) = P(3, 0) : Ligne(11) = P(3, 1)
                    Ligne(12) = P(3, 0) : Ligne(13) = P(3, 1) : Ligne(14) = P(0, 0) : Ligne(15) = P(0, 1)

                    MettreFI = True  ' lorsque l'on sort on met une face interne

            End Select


        ElseIf Left(Nomsection, 1) = "S" Or Nomsection = " Poutre en I générique" Then ' poutre en I de type S
            Dim P(8, 1) As Double

            Select Case numero
                Case 1
                    Dim d As Double
                    d = Poutre.GetD4 * 0.8660254038 '  section.D4 / (2 * tan(30))

                    P(0, 0) = 0
                    P(0, 1) = 0
                    P(1, 0) = -d
                    P(1, 1) = -Poutre.GetD4 / 2.0R
                    P(2, 0) = (Poutre.GetD1 / 2) - Poutre.GetD3
                    P(2, 1) = -Poutre.GetD4 / 2.0R
                    P(3, 0) = P(2, 0)
                    P(3, 1) = -Poutre.GetD2 / 2
                    P(4, 0) = Poutre.GetD1 / 2
                    P(4, 1) = P(3, 1)
                    P(5, 0) = P(4, 0)
                    P(5, 1) = -P(4, 1)
                    P(6, 0) = P(3, 0)
                    P(6, 1) = -P(3, 1)
                    P(7, 0) = P(2, 0)
                    P(7, 1) = -P(2, 1)
                    P(8, 0) = P(1, 0)
                    P(8, 1) = -P(1, 1)

                    pt3(0) -= TranslationX
                    pt3(1) -= TranslationY
                    pt3(0) /= longueur : pt3(1) /= longueur
                    For i = 0 To 8
                        Outils_Math.Rotation2D(pt3, P(i, 0), P(i, 1))
                        P(i, 0) += TranslationX
                        P(i, 1) += TranslationY
                    Next i

                    ReDim Ligne(35)
                    For i = 0 To 7
                        Ligne(i * 4) = P(i, 0) : Ligne(i * 4 + 1) = P(i, 1) : Ligne(i * 4 + 2) = P(i + 1, 0) : Ligne(i * 4 + 3) = P(i + 1, 1)
                    Next i
                    Ligne(32) = P(8, 0) : Ligne(33) = P(8, 1) : Ligne(34) = P(0, 0) : Ligne(35) = P(0, 1)
                    r(0) = inter.x + 5000 * Epsilon * IP(0)
                    r(1) = inter.y + 5000 * Epsilon * IP(1)
                    r(2) = inter.z + 5000 * Epsilon * IP(2)
                Case 2

                    Dim d As Double
                    d = Poutre.GetD4 * 0.8660254038 '  section.D4 / (2 * tan(30))

                    P(0, 0) = 0
                    P(0, 1) = 0

                    P(1, 0) = -d
                    P(1, 1) = -Poutre.GetD4 / 2.0R
                    P(2, 0) = -((Poutre.GetD1 / 2) - Poutre.GetD3)
                    P(2, 1) = -Poutre.GetD4 / 2.0R
                    P(3, 0) = P(2, 0)
                    P(3, 1) = -Poutre.GetD2 / 2
                    P(4, 0) = -Poutre.GetD1 / 2
                    P(4, 1) = P(3, 1)
                    P(5, 0) = P(4, 0)
                    P(5, 1) = -P(4, 1)
                    P(6, 0) = P(3, 0)
                    P(6, 1) = -P(3, 1)
                    P(7, 0) = P(2, 0)
                    P(7, 1) = -P(2, 1)
                    P(8, 0) = P(1, 0)
                    P(8, 1) = -P(1, 1)

                    pt3(0) -= TranslationX
                    pt3(1) -= TranslationY
                    pt3(0) /= longueur : pt3(1) /= longueur
                    For i = 0 To 8
                        Outils_Math.Rotation2D(pt3, P(i, 0), P(i, 1))
                        P(i, 0) += TranslationX
                        P(i, 1) += TranslationY
                    Next i

                    ReDim Ligne(35)
                    For i = 0 To 7
                        Ligne(i * 4) = P(i, 0) : Ligne(i * 4 + 1) = P(i, 1) : Ligne(i * 4 + 2) = P(i + 1, 0) : Ligne(i * 4 + 3) = P(i + 1, 1)
                    Next i
                    Ligne(32) = P(8, 0) : Ligne(33) = P(8, 1) : Ligne(34) = P(0, 0) : Ligne(35) = P(0, 1)

                    r(0) = inter.x - 5000 * Epsilon * IP(0)
                    r(1) = inter.y - 5000 * Epsilon * IP(1)
                    r(2) = inter.z - 5000 * Epsilon * IP(2)

            End Select
            MettreFI = True

        ElseIf Left(Nomsection, 4) = "Tube" Or Nomsection = " Tuyau (Pipe) générique" Then ' le tube rond
            Dim p(4, 1) As Double
            p(0, 0) = Poutre.GetD1 / 2 - Poutre.GetD3
            p(0, 1) = 0
            p(1, 0) = Poutre.GetD1 / 2
            p(1, 1) = 0
            p(2, 0) = -p(0, 0)
            p(2, 1) = 0
            p(3, 0) = -p(1, 0)
            p(3, 1) = 0
            p(4, 0) = 0
            p(4, 1) = 0

            Select Case numero
                Case 1

                    r(0) = 0
                    r(1) = Poutre.GetD1 / 2 - Poutre.GetD3 / 2 : r(2) = 0
                    Outils_Math.Rotation2D(pt3, r(0), r(1))
                    r(0) += TranslationX
                    r(1) += TranslationY
                    r = Commun.TransfertSketchToModel(swSketch, r)

                    pt3(0) -= TranslationX
                    pt3(1) -= TranslationY
                    pt3(0) /= longueur : pt3(1) /= longueur
                    For i = 0 To 4
                        Outils_Math.Rotation2D(pt3, p(i, 0), p(i, 1))
                        p(i, 0) += TranslationX
                        p(i, 1) += TranslationY
                    Next i

                    ReDim Ligne(7)
                    Ligne(0) = p(0, 0) : Ligne(1) = p(0, 1) : Ligne(2) = p(1, 0) : Ligne(3) = p(1, 1)
                    Ligne(4) = p(2, 0) : Ligne(5) = p(2, 1) : Ligne(6) = p(3, 0) : Ligne(7) = p(3, 1)
                    swModel.CreateArc2(p(4, 0), p(4, 1), 0, p(1, 0), p(1, 1), 0, p(3, 0), p(3, 1), 0, 1)
                    swModel.CreateArc2(p(4, 0), p(4, 1), 0, p(0, 0), p(0, 1), 0, p(2, 0), p(2, 1), 0, 1)
                    MettreFI = True
                    'Flag = 2
                Case 2

                    r(0) = 0
                    r(1) = -Poutre.GetD1 / 2 + Poutre.GetD3 / 2 : r(2) = 0
                    Outils_Math.Rotation2D(pt3, r(0), r(1))
                    r(0) += TranslationX
                    r(1) += TranslationY
                    r = Commun.TransfertSketchToModel(swSketch, r)

                    pt3(0) -= TranslationX
                    pt3(1) -= TranslationY
                    pt3(0) /= longueur : pt3(1) /= longueur
                    For i = 0 To 4
                        Outils_Math.Rotation2D(pt3, p(i, 0), p(i, 1))
                        p(i, 0) += TranslationX
                        p(i, 1) += TranslationY
                    Next i

                    ReDim Ligne(7)
                    Ligne(0) = p(0, 0) : Ligne(1) = p(0, 1) : Ligne(2) = p(1, 0) : Ligne(3) = p(1, 1)
                    Ligne(4) = p(2, 0) : Ligne(5) = p(2, 1) : Ligne(6) = p(3, 0) : Ligne(7) = p(3, 1)
                    swModel.CreateArc2(p(4, 0), p(4, 1), 0, p(1, 0), p(1, 1), 0, p(3, 0), p(3, 1), 0, -1)
                    swModel.CreateArc2(p(4, 0), p(4, 1), 0, p(0, 0), p(0, 1), 0, p(2, 0), p(2, 1), 0, -1)
                    'MettreFI = True  ' lorsque l'on sort on met une face interne
                    MettreFI = False
                    Me.Flag = 2
                    'Case 1 ' le cercle extérieur
                    '    swModel.CreateCircleByRadius2(TranslationX, TranslationY, 0, Poutre.GetD1 / 2)
                    '    MettreFI = False
                    '    r(0) = 0 : r(1) = 0 : r(2) = 0
                    '    r = Commun.TransfertSketchToModel(swSketch, r)
                    'Case 2
                    '    swModel.CreateCircleByRadius2(TranslationX, TranslationY, 0, (Poutre.GetD1 / 2) - Poutre.GetD3)
                    '    r(0) = (Poutre.GetD1 / 2 - (Poutre.GetD3 / 2))
                    '    r(1) = 0 : r(2) = 0
                    '    r = Commun.TransfertSketchToModel(swSketch, r)
                    '    MettreFI = True
            End Select

        ElseIf Left(Nomsection, 1) = "C" Or Nomsection = " Poutre en C générique" Then ' le channel
            Dim P(7, 1) As Double

            Select Case numero
                Case 1 ' le C au complet

                    P(0, 0) = Poutre.GetD1 / 2 - Poutre.GetD3
                    P(0, 1) = Poutre.GetD5
                    P(1, 0) = P(0, 0)
                    P(1, 1) = Poutre.GetD5 + Poutre.GetD4 - Poutre.GetD2
                    P(2, 0) = Poutre.GetD1 / 2
                    P(2, 1) = P(1, 1)
                    P(3, 0) = P(2, 0)
                    P(3, 1) = P(1, 1) + Poutre.GetD2
                    P(4, 0) = -P(3, 0)
                    P(4, 1) = P(3, 1)
                    P(5, 0) = P(4, 0)
                    P(5, 1) = P(1, 1)
                    P(6, 0) = -P(1, 0)
                    P(6, 1) = P(5, 1)
                    P(7, 0) = -P(0, 0)
                    P(7, 1) = P(0, 1)

                    r(0) = P(0, 0) + 1000 * Epsilon
                    r(1) = 0 : r(2) = 0
                    Outils_Math.Rotation2D(pt3, r(0), r(1))
                    r(0) += TranslationX
                    r(1) += TranslationY
                    r = Commun.TransfertSketchToModel(swSketch, r)

                    pt3(0) -= TranslationX
                    pt3(1) -= TranslationY
                    pt3(0) /= longueur : pt3(1) /= longueur
                    For i = 0 To 7
                        Outils_Math.Rotation2D(pt3, P(i, 0), P(i, 1))
                        P(i, 0) += TranslationX
                        P(i, 1) += TranslationY
                    Next i

                    ReDim Ligne(35)
                    For i = 0 To 6
                        Ligne(i * 4) = P(i, 0) : Ligne(i * 4 + 1) = P(i, 1) : Ligne(i * 4 + 2) = P(i + 1, 0) : Ligne(i * 4 + 3) = P(i + 1, 1)
                    Next i
                    Ligne(28) = P(7, 0) : Ligne(29) = P(7, 1) : Ligne(30) = P(0, 0) : Ligne(31) = P(0, 1)

                    MettreFI = False
                    Me.Flag = 2
                Case 2
                    MettreFI = False ' Attention, peut planter à cause de ça.
            End Select

        ElseIf Left(Nomsection, 1) = "L" Or Nomsection = " Poutre en L générique" Then ' l'Angle en L
            Dim P(5, 1) As Double

            Select Case numero
                Case 1 ' le C au complet

                    P(0, 0) = -Poutre.GetD5 + Poutre.GetD1
                    P(0, 1) = -Poutre.GetD6 + Poutre.GetD4
                    P(1, 0) = -Poutre.GetD5 + Poutre.GetD3
                    P(1, 1) = P(0, 1)
                    P(2, 0) = P(1, 0)
                    P(2, 1) = -Poutre.GetD6 + Poutre.GetD2
                    P(3, 0) = -Poutre.GetD5
                    P(3, 1) = P(2, 1)
                    P(4, 0) = P(3, 0)
                    P(4, 1) = -Poutre.GetD6
                    P(5, 0) = P(0, 0)
                    P(5, 1) = P(4, 1)

                    r(0) = P(1, 0) - 1000 * Epsilon
                    r(1) = 0 : r(2) = 0
                    Outils_Math.Rotation2D(pt3, r(0), r(1))
                    r(0) += TranslationX
                    r(1) += TranslationY
                    r = Commun.TransfertSketchToModel(swSketch, r)

                    pt3(0) -= TranslationX
                    pt3(1) -= TranslationY
                    pt3(0) /= longueur : pt3(1) /= longueur
                    For i = 0 To 5
                        Outils_Math.Rotation2D(pt3, P(i, 0), P(i, 1))
                        P(i, 0) += TranslationX
                        P(i, 1) += TranslationY
                    Next i

                    ReDim Ligne(35)
                    For i = 0 To 4
                        Ligne(i * 4) = P(i, 0) : Ligne(i * 4 + 1) = P(i, 1) : Ligne(i * 4 + 2) = P(i + 1, 0) : Ligne(i * 4 + 3) = P(i + 1, 1)
                    Next i
                    Ligne(20) = P(5, 0) : Ligne(21) = P(5, 1) : Ligne(22) = P(0, 0) : Ligne(23) = P(0, 1)

                    MettreFI = False  ' lorsque l'on sort on met une face interne
                    Me.Flag = 2
                Case 2
                    MettreFI = False ' Attention, peut planter à cause de ça.
            End Select


        ElseIf Left(Nomsection, 1) = "T" Or Nomsection = " Poutre en T générique" Then ' le T
            Select Case numero
                Case 1
                    Dim P(8, 1) As Double
                    Dim d As Double
                    d = Poutre.GetD4 * 0.8660254038 '  section.D4 / (2 * tan(30))


                    P(0, 0) = 0
                    P(0, 1) = 0
                    P(1, 0) = -d
                    P(1, 1) = -Poutre.GetD4 / 2.0R
                    P(2, 0) = -(Poutre.GetD1 - Poutre.GetD5 - Poutre.GetD3)
                    P(2, 1) = -Poutre.GetD4 / 2.0R
                    P(3, 0) = P(2, 0)
                    P(3, 1) = -Poutre.GetD2 / 2
                    P(4, 0) = -Poutre.GetD1 + Poutre.GetD5
                    P(4, 1) = P(3, 1)
                    P(5, 0) = P(4, 0)
                    P(5, 1) = -P(4, 1)
                    P(6, 0) = P(3, 0)
                    P(6, 1) = -P(3, 1)
                    P(7, 0) = P(2, 0)
                    P(7, 1) = -P(2, 1)
                    P(8, 0) = P(1, 0)
                    P(8, 1) = -P(1, 1)

                    r(0) = P(0, 0) + 1000 * Epsilon
                    r(1) = 0 : r(2) = 0
                    Outils_Math.Rotation2D(pt3, r(0), r(1))
                    r(0) += TranslationX
                    r(1) += TranslationY
                    r = Commun.TransfertSketchToModel(swSketch, r)

                    pt3(0) -= TranslationX
                    pt3(1) -= TranslationY
                    pt3(0) /= longueur : pt3(1) /= longueur
                    For i = 0 To 8
                        Outils_Math.Rotation2D(pt3, P(i, 0), P(i, 1))
                        P(i, 0) += TranslationX
                        P(i, 1) += TranslationY
                    Next i

                    ReDim Ligne(35)
                    For i = 0 To 7
                        Ligne(i * 4) = P(i, 0) : Ligne(i * 4 + 1) = P(i, 1) : Ligne(i * 4 + 2) = P(i + 1, 0) : Ligne(i * 4 + 3) = P(i + 1, 1)
                    Next i
                    Ligne(32) = P(8, 0) : Ligne(33) = P(8, 1) : Ligne(34) = P(0, 0) : Ligne(35) = P(0, 1)

                Case 2
                    Dim P(4, 1) As Double
                    Dim d As Double
                    d = Poutre.GetD4 * 0.8660254038 '  section.D4 / (2 * tan(30))

                    P(0, 0) = 0
                    P(0, 1) = 0

                    P(1, 0) = -d
                    P(1, 1) = -Poutre.GetD4 / 2.0R
                    P(2, 0) = Poutre.GetD5
                    P(2, 1) = -Poutre.GetD4 / 2.0R
                    P(3, 0) = P(2, 0)
                    P(3, 1) = Poutre.GetD4 / 2
                    P(4, 0) = P(1, 0)
                    P(4, 1) = P(3, 1)

                    r(0) = P(0, 0) - 1000 * Epsilon
                    r(1) = 0 : r(2) = 0
                    Outils_Math.Rotation2D(pt3, r(0), r(1))
                    r(0) += TranslationX
                    r(1) += TranslationY
                    r = Commun.TransfertSketchToModel(swSketch, r)

                    pt3(0) -= TranslationX
                    pt3(1) -= TranslationY
                    pt3(0) /= longueur : pt3(1) /= longueur
                    For i = 0 To 4
                        Outils_Math.Rotation2D(pt3, P(i, 0), P(i, 1))
                        P(i, 0) += TranslationX
                        P(i, 1) += TranslationY
                    Next i

                    ReDim Ligne(19)
                    For i = 0 To 3
                        Ligne(i * 4) = P(i, 0) : Ligne(i * 4 + 1) = P(i, 1) : Ligne(i * 4 + 2) = P(i + 1, 0) : Ligne(i * 4 + 3) = P(i + 1, 1)
                    Next i
                    Ligne(16) = P(4, 0) : Ligne(17) = P(4, 1) : Ligne(18) = P(0, 0) : Ligne(19) = P(0, 1)

            End Select
            MettreFI = True

        ElseIf Left(Poutre.GetNomSection, 1) = "P" Or Nomsection = " Circulaire pleine générique" Then ' Pipe, 
            Dim P(2, 1) As Double
            Dim d As Double, e As Double
            d = Poutre.GetD1 / 4 ' Math.Sin(30)  ( et on doit diviser le diamètre par 2)
            e = Poutre.GetD1 * Math.Sqrt(3) / 4 ' cos (30°)

            P(0, 0) = 0
            P(0, 1) = 0
            P(1, 0) = d
            P(1, 1) = -e
            P(2, 0) = d
            P(2, 1) = e

            Select Case numero
                Case 1
                    r(0) = P(0, 0) + 1000 * Epsilon
                    r(1) = 0 : r(2) = 0
                    Outils_Math.Rotation2D(pt3, r(0), r(1))
                    r(0) += TranslationX
                    r(1) += TranslationY
                    r = Commun.TransfertSketchToModel(swSketch, r)
                    pt3(0) -= TranslationX
                    pt3(1) -= TranslationY
                    pt3(0) /= longueur : pt3(1) /= longueur
                    For i = 0 To 2
                        Outils_Math.Rotation2D(pt3, P(i, 0), P(i, 1))
                        P(i, 0) += TranslationX
                        P(i, 1) += TranslationY
                    Next i
                    ReDim Ligne(7)
                    Ligne(0) = P(0, 0) : Ligne(1) = P(0, 1) : Ligne(2) = P(1, 0) : Ligne(3) = P(1, 1)
                    Ligne(4) = P(0, 0) : Ligne(5) = P(0, 1) : Ligne(6) = P(2, 0) : Ligne(7) = P(2, 1)
                    swModel.CreateArc2(P(0, 0), P(0, 1), 0, P(1, 0), P(1, 1), 0, P(2, 0), P(2, 1), 0, 1) ' le dernier param est la direction.  1 ou -1

                Case 2
                    r(0) = P(0, 0) - 1000 * Epsilon
                    r(1) = 0 : r(2) = 0
                    Outils_Math.Rotation2D(pt3, r(0), r(1))
                    r(0) += TranslationX
                    r(1) += TranslationY
                    r = Commun.TransfertSketchToModel(swSketch, r)

                    pt3(0) -= TranslationX
                    pt3(1) -= TranslationY
                    pt3(0) /= longueur : pt3(1) /= longueur
                    For i = 0 To 2
                        Outils_Math.Rotation2D(pt3, P(i, 0), P(i, 1))
                        P(i, 0) += TranslationX
                        P(i, 1) += TranslationY
                    Next i
                    ReDim Ligne(7)
                    Ligne(0) = P(0, 0) : Ligne(1) = P(0, 1) : Ligne(2) = P(1, 0) : Ligne(3) = P(1, 1)
                    Ligne(4) = P(0, 0) : Ligne(5) = P(0, 1) : Ligne(6) = P(2, 0) : Ligne(7) = P(2, 1)

                    swModel.CreateArc2(P(0, 0), P(0, 1), 0, P(1, 0), P(1, 1), 0, P(2, 0), P(2, 1), 0, -1) ' le dernier param est la direction.  1 ou -1

            End Select
            MettreFI = True

        Else
            MsgBox("Section de poutre non reconnu!", MsgBoxStyle.Critical, "Commun.DessineSectionPoutre")
        End If


        If Not Ligne Is Nothing Then
            For i = 0 To UBound(Ligne) Step 4
                sketchline = swModel.CreateLine2(Ligne(i), Ligne(i + 1), 0, Ligne(i + 2), Ligne(i + 3), 0)
            Next i
        End If

        Return r

    End Function

    Public Function SwFace() As SldWorks.Face2 ' retourne la première face de la liste (dans la partie traitement, ce sera la seule...)
        Return Me.lst_Faces.Item(0)
    End Function

    Public Function IsFaceInterne(ByRef swface As SldWorks.Face2) As Boolean
        Dim attr As SldWorks.Attribute
        Dim SwEnt As SldWorks.Entity
        SwEnt = swface
        attr = SwEnt.FindAttribute(DefAttrFaceInterne, 0)
        If attr Is Nothing Then Return False Else Return True
    End Function

    ''' <summary>
    ''' Fonction qui retourne un tableau de Sldworks.edge (et non slyedges)
    ''' </summary>
    ''' <returns>Un tableau de Edges</returns>
    ''' <remarks></remarks>
    Public Function GetAretes() As SldWorks.Edge()
        Dim face As SldWorks.Face2
        Dim arete As SldWorks.Edge = Nothing
        Dim temp2 As Collections.Generic.List(Of SldWorks.Edge)
        Dim lst As New Collections.Generic.List(Of SldWorks.Edge)

        For Each face In Me.lst_Faces
            temp2 = GetArete1Face(face)
            For Each arete In temp2
                lst.Add(arete)
            Next arete
        Next face

        Return lst.ToArray
    End Function

    Private Function GetArete1Face(ByRef Face As SldWorks.Face2) As Collections.Generic.List(Of SldWorks.Edge)
        Dim vArete As Object
        Dim a As SldWorks.Edge
        Dim arete() As SldWorks.Edge
        Dim lst As New Collections.Generic.List(Of SldWorks.Edge)

        ReDim arete(Face.GetEdgeCount - 1)
        vArete = Face.GetEdges()

        For Each a In vArete
            lst.Add(a)
        Next

        Return lst
    End Function

    Public Overrides Sub Selectionner(Optional ByVal Mark As Integer = 0, Optional ByRef append As Boolean = True)
        Dim swent As SldWorks.Entity
        Dim swface As SldWorks.Face2

        For Each swface In lst_Faces
            swent = swface
            swent.Select2(append, Mark)
        Next swface
    End Sub


    ''' <summary>
    ''' Sélectionne toutes les faces dans la liste de faces
    ''' </summary>
    ''' <param name="Mark"></param>
    ''' <param name="Append"></param>
    ''' <remarks></remarks>
    Public Sub SelectionnerToutes(Optional ByRef Mark As Integer = 0, Optional ByRef Append As Boolean = True)
        Dim swFace As SldWorks.Face2

        Dim swent As SldWorks.Entity
        If Append = False Then swModel.ClearSelection2(True)
        For Each swFace In Me.lst_Faces
            swent = swFace : swent.Select2(True, Mark)
        Next

    End Sub


    Public Function Couleur(ByRef rouge As Double, ByRef Vert As Double, ByRef Bleu As Double, Optional ByVal Ambient As Double = 1, Optional ByVal Diffuse As Double = 1, Optional ByVal Specular As Double = 1, Optional ByVal Shininess As Double = 0.5, Optional ByVal Transparency As Double = 0, Optional ByVal Emission As Double = 0.2) As Integer

        swModel.SelectedFaceProperties(RGB(rouge, Vert, Bleu), Ambient, Diffuse, Specular, Shininess, Transparency, Emission, False, "")

        Return 1
    End Function

    Public Sub AjouterFace(ByRef face As SldWorks.Face2)
        Dim testface As SldWorks.Face2
        Dim faultentity As SldWorks.FaultEntity
        Dim swent As SldWorks.Entity

        If Not Me.lst_Faces.Contains(face) Then Me.lst_Faces.Add(face)

        ' vérifier que les anciennes faces sont toujours ok...
        For Each testface In Me.lst_Faces
            faultentity = testface.Check
            If Not faultentity.Count = 0 Then ' on a un problème avec la face....
                lst_Faces.Remove(testface)
                Dim i As Integer
                For i = 0 To faultentity.Count
                    swent = faultentity.Entity(i)
                    If Not swent Is Nothing Then
                        swent.Select4(True, Nothing)
                    End If
                    Debug.Print("   Fault[" & i & "] = " & swent.errorCode(i))
                Next i
            End If
        Next testface


    End Sub


    'Public Function DonnerFaces() As SldWorks.Face2()
    '    'Dim temp() As SldWorks.Face2
    '    'Dim i As Integer

    '    'ReDim temp(lst_Faces.Count - 1)

    '    'For i = 1 To lst_Faces.Count
    '    '    temp(i - 1) = lst_Faces.Item(i)
    '    'Next

    '    Return lst_Faces.ToArray
    'End Function


    ''' <summary>
    ''' Function qui retourne un pointeur vers la face
    ''' </summary>
    ''' <returns></returns>
    ''' <remarks></remarks>
    Public Function GetFace() As SldWorks.Face2
        Return Me.SwFace
    End Function

    ''' <summary>
    ''' Fonction qui redonne toutes les Faces contenues dans cette face
    ''' </summary>
    ''' <returns></returns>
    ''' <remarks></remarks>
    Public Function GetFaces() As SldWorks.Face2()
        Return Me.lst_Faces.ToArray
    End Function

    ''' <summary>
    ''' Sub qui renvoie les coordonnées min et max des valeurs U et V
    ''' </summary>
    ''' <param name="Umin"></param>
    ''' <param name="UMax"></param>
    ''' <param name="VMin"></param>
    ''' <param name="VMax"></param>
    ''' <remarks></remarks>
    Public Sub UVMinMax(ByRef Umin As Double, ByRef UMax As Double, ByRef VMin As Double, ByRef VMax As Double)
        Dim vBounds As Object
        vBounds = SwFace.GetUVBounds()

        Umin = vBounds(0)
        UMax = vBounds(1)

        VMin = vBounds(2)
        VMax = vBounds(3)
    End Sub

    ''' <summary>
    ''' Calcule la position du point selon les U et V
    ''' </summary>
    ''' <param name="U"></param>
    ''' <param name="V"></param>
    ''' <param name="X"></param>
    ''' <param name="Y"></param>
    ''' <param name="Z"></param>
    ''' <returns>Vrai si l epoint est sur la face, faux sinon</returns>
    ''' <remarks>Retourne X, Y et Z même si le point n'est pas sur la face (mais sur la surface) </remarks>
    Public Function Evaluer(ByRef U As Double, ByVal V As Double, ByRef X As Double, ByRef Y As Double, ByRef Z As Double) As Boolean
        Dim surf As SldWorks.Surface
        Dim vEv As Object, vpoint As Object
        Dim P(2) As Double

        surf = SwFace.GetSurface()

        vEv = surf.Evaluate(U, V, 0, 0)

        X = vEv(0) : Y = vEv(1) : Z = vEv(2)

        vpoint = SwFace.GetClosestPointOn(X, Y, Z)

        If (Math.Abs(vpoint(0) - X) < Epsilon) And (Math.Abs(vpoint(1) - Y) < Epsilon) And (Math.Abs(vpoint(2) - Z) < Epsilon) Then Return True Else Return False

    End Function

    ''' <summary>
    ''' Function qui calcule la normale d'une face au point X,Y,Z.
    ''' </summary>
    ''' <param name="X"></param>
    ''' <param name="Y"></param>
    ''' <param name="Z"></param>
    ''' <returns>Un tableau de 3 doubles correspondant à la normale</returns>
    ''' <remarks></remarks>
    Public Function Normale(ByRef X As Double, ByRef Y As Double, ByRef Z As Double) As Double()
        Dim surf As SldWorks.Surface
        Dim vtemp As Object
        Dim temp() As Double
        Dim sens As Boolean

        surf = SwFace.GetSurface
        If surf.IsPlane Then vtemp = SwFace.Normal : temp = vtemp : Return temp ' si la face est plane alors c'est ok, sinon il faut travailler...


        vtemp = surf.EvaluateAtPoint(X, Y, Z)
        ReDim temp(2)

        ' la normale de la face pointe AWAY from the body

        sens = SwFace.FaceInSurfaceSense()  'TRUE if face normal and surface normal are in the opposite direction and FALSE if they are in the same direction


        If sens Then ' on doit inverser
            temp(0) = -vtemp(0) : temp(1) = -vtemp(1) : temp(2) = -vtemp(2)
        Else
            temp(0) = vtemp(0) : temp(1) = vtemp(1) : temp(2) = vtemp(2)
        End If
        Return temp

    End Function


    Public Sub MettreAttributFaceInterne(Optional ByRef Valeur As Double = 0)
        Dim no As Integer = 0
        Dim nom As String = "FaceInterne" & no
        Dim swent As SldWorks.Entity
        Dim attr As SldWorks.Attribute

        swent = Me.SwFace
        attr = swent.FindAttribute(Intersections.DefAttrFaceInterne, 0) ' si l'attribut existe déjà on pointe dessus.
        If attr Is Nothing Then attr = Intersections.DefAttrFaceInterne.CreateInstance5(swModel, SwFace, nom, 0, 2) ' 0 = swThisconfig
        While attr Is Nothing
            no += 1
            nom = "FaceInterne" & CStr(no)
            attr = Intersections.DefAttrFaceInterne.CreateInstance5(swModel, SwFace, nom, 0, 2)
        End While
        GererDossiers("FaceInternes", nom)

    End Sub


End Class
Revision:	40
Committed:	Mon Aug 20 21:30:28 2007 UTC (17 years, 9 months ago) by bournival
File size:	85633 byte(s)
Log Message:	Projet de these de Sylvain Bournival. Attention projet VB.