REPOS_ERICCA/magicsld/Intersections.vb

Imports SolidWorks.Interop
Imports SolidWorks.Interop.swconst
Imports SolidWorks.Interop.swpublished

Module Intersections
    Public DefAttrInterALAL As SldWorks.AttributeDef
    Public DefAttrConditionLimite As SldWorks.AttributeDef
    Public DefAttrRCP1 As SldWorks.AttributeDef
    Public DefAttrRCCoque As SldWorks.AttributeDef
    Public DefAttrFaceInterne As SldWorks.AttributeDef
    Public DefAttrDoublon As SldWorks.AttributeDef
    Public DefAttrIgnorer As sldworks.AttributeDef

    Public nbMinipoutre As Long

    Public lst_InterCoqueCoque As New Collections.Generic.List(Of InterCoqueCoque)

    Public MettreMiniPoutresSurFaceInternes As Boolean
    Public MultiDecoupageCoques As Boolean = False
    Public SectionSimpleSurPoutre As Boolean

#Region "Enums"
    Public Enum typeInterPoutreVolume
        Centre = 0
        Sur_arete = 1
        Sur_sommet = 2
        Courbe = 3
    End Enum

#End Region


    Public Sub Debuter(Optional ByRef nom2 As String = "C:\Documents and Settings\Sylvain Bournival\Bureau\testMAGiC.sldprt", Optional ByVal nomFichier As String = "", Optional ByRef original As Boolean = True, Optional ByRef modifie As Boolean = True)

        ' *******
        ' quelques options de performance
        ' *******
        swApp.SetUserPreferenceIntegerValue(swconst.swUserPreferenceIntegerValue_e.swAutoSaveInterval, 0)
        swModel.SetAddToDB(True)
        swModel.SetDisplayWhenAdded(False)
        swModel.SetInferenceMode(False)
        ' ******
        ' fin des options de performance
        ' ******


        Memoriser3iemePoint() ' mémorise le coord system car si on découpe, sa coordonnée est perdue.

        CouperPoutres()
        Commun.GenererListes() ' va ignorer les poutres à ignorer... et ajouter les poutres coupées dans la liste.

        lst_InterCoqueCoque.Clear()
        DetectionCoqueCoque()
        DetectionPoutresVolumes() ' doit être avant interpoutreCoque au cas où on aurait une poutre de section
        DetectionPoutresCoques()
        DetectionCoqueVolume()


        ' Traitement des intersection poutres-Volumes
        DecouperPoutreVolume()
        swModel.EditRebuild3()

        ' Traitement des intersection entre poutre et coques
        DecouperPoutreCoque()
        swModel.EditRebuild3()

        ' traitement des coques-volumes
        DécouperCoqueVolume()
        swModel.EditRebuild3()

        ' traitement des coques-coques
        DecouperCoqueCoque()
        swModel.EditRebuild3()


        ' *******
        ' quelques options de performance, remettre à la position initiale
        ' *******
        swApp.SetUserPreferenceIntegerValue(swconst.swUserPreferenceIntegerValue_e.swAutoSaveInterval, 15)
        swModel.SetAddToDB(False)
        swModel.SetDisplayWhenAdded(True)
        swModel.GraphicsRedraw2()
        swModel.SetInferenceMode(True)
        ' ******
        ' fin des options de performance
        ' ******
        ReplacerFolder()


    End Sub


    ''' <summary>
    ''' Sub qui découpe les coques en fonction des informations placées dans le InterCoqueCoque
    ''' </summary>
    ''' <remarks></remarks>
    Private Sub DecouperCoqueCoque()
        Dim rayon As Double

        For Each interCC As InterCoqueCoque In lst_InterCoqueCoque
            'rayon = IIf(Coque Is interCC.sFaceCoque1, interCC.sFaceCoque2.GetEpaisseur, interCC.sFaceCoque1.GetEpaisseur)
            'MsgBox("Traitement de l'intersection de coque - coque # " & interCC.Numero)

            Try
                If interCC.DoitCouperCoque1 Then
                    rayon = interCC.sFaceCoque2.GetEpaisseur / 2
                    Dim sweep As sldworks.Body2 = interCC.GénérerSweep(interCC.sketch, rayon)
                    interCC.DecouperCoque(interCC.sFaceCoque1, sweep)
                    If Intersections.MettreMiniPoutresSurFaceInternes Then interCC.MarquerFacesInternes(interCC.sFaceCoque2, interCC.sFaceCoque1)
                End If

                If interCC.DoitCouperCoque2 Then
                    rayon = interCC.sFaceCoque1.GetEpaisseur / 2
                    Dim sweep As sldworks.Body2 = interCC.GénérerSweep(interCC.sketch, rayon)
                    interCC.DecouperCoque(interCC.sFaceCoque2, sweep)
                    If Intersections.MettreMiniPoutresSurFaceInternes Then interCC.MarquerFacesInternes(interCC.sFaceCoque1, interCC.sFaceCoque2)
                End If
            Catch
                If interCC.DoitCouperCoque2 Then
                    rayon = interCC.sFaceCoque1.GetEpaisseur / 2
                    Dim sweep As sldworks.Body2 = interCC.GénérerSweep(interCC.sketch, rayon)
                    interCC.DecouperCoque(interCC.sFaceCoque2, sweep)
                    If Intersections.MettreMiniPoutresSurFaceInternes Then interCC.MarquerFacesInternes(interCC.sFaceCoque1, interCC.sFaceCoque2)
                End If

                If interCC.DoitCouperCoque1 Then
                    rayon = interCC.sFaceCoque2.GetEpaisseur / 2
                    Dim sweep As sldworks.Body2 = interCC.GénérerSweep(interCC.sketch, rayon)
                    interCC.DecouperCoque(interCC.sFaceCoque1, sweep)
                    If Intersections.MettreMiniPoutresSurFaceInternes Then interCC.MarquerFacesInternes(interCC.sFaceCoque2, interCC.sFaceCoque1)
                End If


            End Try


            If interCC.FaceAPlat Then
                interCC.CoupeAPlat()
            End If

            ' reste à retrouver les faces internes.
            'interCC.MarquerFacesInternes(IIf(Coque Is interCC.sFaceCoque1, interCC.sFaceCoque2, interCC.sFaceCoque1), IIf(Coque Is interCC.sFaceCoque1, interCC.sFaceCoque1, interCC.sFaceCoque2))
        Next


    End Sub

    'Private Sub DecouperCoqueCoque()
    '    Dim rayon As Double


    '    For Each Coque As SlyFaceCoque In Commun.lst_FaceCoque
    '        For Each interCC As InterCoqueCoque In Coque.lst_InterCoqueCoque
    '            rayon = IIf(Coque Is interCC.sFaceCoque1, interCC.sFaceCoque2.GetEpaisseur, interCC.sFaceCoque1.GetEpaisseur)
    '            Dim sweep As SldWorks.Body2 = interCC.GénérerSweep(interCC.sketch, rayon)
    '            interCC.DecouperCoque(Coque, sweep)

    '            ' reste à retrouver les faces internes.
    '            interCC.MarquerFacesInternes(IIf(Coque Is interCC.sFaceCoque1, interCC.sFaceCoque2, interCC.sFaceCoque1), IIf(Coque Is interCC.sFaceCoque1, interCC.sFaceCoque1, interCC.sFaceCoque2))
    '        Next
    '    Next


    'End Sub


    ''' <summary>
    '''  sub qui créé une instance de la classe interCoqueCoque s'il y a une intersection de ce type
    ''' </summary>
    ''' <remarks></remarks>
    Private Sub DetectionCoqueCoque()

        Dim sketch As SldWorks.Sketch = Nothing
        Dim interCC As InterCoqueCoque = Nothing
        Dim Coque1 As SlyFaceCoque, Coque2 As SlyFaceCoque

        For i As Integer = 0 To Commun.lst_FaceCoque.Count - 2 'For Each Coque1 As SlyFaceCoque In Commun.lst_FaceCoque
            Coque1 = Commun.lst_FaceCoque.Item(i)
            For j As Integer = i + 1 To Commun.lst_FaceCoque.Count - 1 '  For Each coque2 As SlyFaceCoque In lst_FaceCoque
                Coque2 = Commun.lst_FaceCoque.Item(j)
                If DetectFaceFace(Coque2.SwFace, Coque1.SwFace, True, sketch) Then
                    ' création de l'instance de interFace-face entre coque et coque

                    interCC = New InterCoqueCoque(Coque1, Coque2)
                    interCC.FaceDeSection = False
                    interCC.sketch = sketch
                    interCC.determineType()
                    lst_InterCoqueCoque.Add(interCC)


                    'Coque1.lst_InterCoqueCoque.Add(interCC)
                    'Coque2.lst_InterCoqueCoque.Add(interCC)
                End If

            Next j
        Next i

    End Sub


    Private Sub ReplacerFolder()
        ' on doit mettre les folders à la fin pour que ça marche dans MAGiC
        'Si on ne met pas les attributs à la fin on est baisé...

        Dim swFeat As SldWorks.Feature = Nothing
        Dim nomdernier As String
        Dim dansFolder As Boolean = False
        Dim i As Integer = 0
        Dim ok As Boolean = False
        Dim SelMgr As SldWorks.SelectionMgr

        'trouver le premier feature qui n'est pas un folder...
        Do Until ok
            swFeat = swModel.FeatureByPositionReverse(i)
            nomdernier = swFeat.GetTypeName

            If (nomdernier = "FtrFolder") Then
                dansFolder = Not dansFolder
            Else
                If Not dansFolder Then ok = True
            End If
            i += 1
        Loop

        nomdernier = swFeat.Name

        SelMgr = swModel.SelectionManager
        If swModel.Extension.SelectByID2("Poutres", "FTRFOLDER", 0, 0, 0, False, 0, Nothing, 0) Then swModel.ReorderFeature("Poutres", nomdernier)
        If swModel.Extension.SelectByID2("Coques", "FTRFOLDER", 0, 0, 0, False, 0, Nothing, 0) Then swModel.ReorderFeature("Coques", nomdernier)
        If swModel.Extension.SelectByID2("Conditions Aux Limites", "FTRFOLDER", 0, 0, 0, False, 0, Nothing, 0) Then swModel.ReorderFeature("Conditions Aux Limites", nomdernier)
        If swModel.Extension.SelectByID2("FaceInternes", "FTRFOLDER", 0, 0, 0, False, 0, Nothing, 0) Then swModel.ReorderFeature("FaceInternes", nomdernier)

    End Sub


    Public Sub RegisterAttribut()
        Static nouveau As Boolean
        'propriétés des intersections entre 2 poutres (Arete-libre Arête-libre)

        If nouveau Then Exit Sub
        Dim nom As String
        Dim retval As Boolean

        nom = "InterALAL"
        DefAttrInterALAL = swApp.DefineAttribute(nom)
        DefAttrInterALAL.AddParameter("X", swconst.swParamType_e.swParamTypeDouble, 0, 0)
        DefAttrInterALAL.AddParameter("Y", swconst.swParamType_e.swParamTypeDouble, 0, 0)
        DefAttrInterALAL.AddParameter("Z", swconst.swParamType_e.swParamTypeDouble, 0, 0)
        DefAttrInterALAL.AddParameter("T", swconst.swParamType_e.swParamTypeDouble, -1, 0)
        retval = DefAttrInterALAL.Register()
        If retval = False Then MsgBox("Enregistrement raté pour le InterALAL")


        nom = "ConditionLimite"
        DefAttrConditionLimite = swApp.DefineAttribute(nom)
        DefAttrConditionLimite.AddParameter("CL", swconst.swParamType_e.swParamTypeString, 0, 0)
        retval = DefAttrConditionLimite.Register()
        If retval = False Then MsgBox("Enregistrement raté pour le COndition Limite")


        DefAttrRCP1 = swApp.DefineAttribute("Poutre1")
        DefAttrRCP1.AddParameter("M", 1, 0, 0) ' 0 = double, 1 = string, 2 = integer
        DefAttrRCP1.AddParameter("S", 1, 0, 0)
        DefAttrRCP1.AddParameter("As", 0, 0, 0)
        DefAttrRCP1.AddParameter("I1", 0, 0, 0)
        DefAttrRCP1.AddParameter("I2", 0, 0, 0)
        DefAttrRCP1.AddParameter("N3", 1, 0, 0) ' le nom du troisième point
        DefAttrRCP1.AddParameter("X3", 0, 0, 0) ' le x du troisième point
        DefAttrRCP1.AddParameter("Y3", 0, 0, 0) ' le y du troisième point
        DefAttrRCP1.AddParameter("Z3", 0, 0, 0) ' le z du troisième point
        DefAttrRCP1.AddParameter("D1", 0, 0, 0)
        DefAttrRCP1.AddParameter("D2", 0, 0, 0)
        DefAttrRCP1.AddParameter("D3", 0, 0, 0)
        DefAttrRCP1.AddParameter("D4", 0, 0, 0)
        DefAttrRCP1.AddParameter("D5", 0, 0, 0)
        DefAttrRCP1.AddParameter("D6", 0, 0, 0)
        DefAttrRCP1.AddParameter("Flag", 0, 0, 0)
        retval = DefAttrRCP1.Register()
        If retval = False Then MsgBox("Enregistrement raté pour le RCPoutre")

        DefAttrRCCoque = swApp.DefineAttribute("Coque")
        DefAttrRCCoque.AddParameter("M", 1, 0, 0) 'le matériau' 0 = double, 1 = string, 2 = integer
        DefAttrRCCoque.AddParameter("Ep", 0, 0, 0) ' épaisseur
        DefAttrRCCoque.AddParameter("Flag", 0, 0, 0) ' un flag pour les coques...
        retval = DefAttrRCCoque.Register()
        If retval = False Then MsgBox("Enregistrement raté pour le RCCoque")

        DefAttrFaceInterne = swApp.DefineAttribute("FaceInterne")
        DefAttrFaceInterne.AddParameter("FI", 0, 0, 0)  ' la taille des éléments
        DefAttrFaceInterne.AddParameter("Po", 0, 0, 0)  ' =0 si poutre, =1 si coque 
        retval = DefAttrFaceInterne.Register()
        If retval = False Then MsgBox("Enregistrement raté pour le FaceInterne")

        DefAttrDoublon = swApp.DefineAttribute("Doublon")
        DefAttrDoublon.AddParameter("Maitre", 1, 0, 0)  ' le nom du maitre...
        DefAttrDoublon.AddParameter("Sens", 0, 0, 0)  ' -1 si dans le sens oposé
        retval = DefAttrDoublon.Register()
        If retval = False Then MsgBox("Enregistrement raté pour le Doublon")

        nom = "Ignorer"
        DefAttrIgnorer = swApp.DefineAttribute(nom)
        DefAttrIgnorer.AddParameter("Rien", 0, 0, 0)
        retval = DefAttrIgnorer.Register()
        If retval = False Then MsgBox("Enregistrement raté pour le Ignorer")


        nouveau = True

    End Sub

    Private Sub TraiteAPAP()
        ' procédure plus utilisée.
        ' procédure qui trouve les intersections entre les poutres et qui créé un attribut contenant les informations

        ' en premier on efface les features attributs, pour ne pas avoir de doubles
        Dim feat As SldWorks.Feature
        feat = swPart.FirstFeature
        Dim nextfeat As SldWorks.Feature

        Do While Not feat Is Nothing
            If Left(feat.Name, 9) = "InterAPAP" Then
                nextfeat = feat.GetNextFeature
                swModel.Extension.SelectByID2(feat.Name, "ATTRIBUTE", 0, 0, 0, False, 0, Nothing, 0)
                swPart.EditDelete()
                feat = nextfeat
            Else
                feat = feat.GetNextFeature
            End If
        Loop
        ' faut aussi détruire les liste de points dans les poutres
        ' ce qui suit n'est pas optimisé en temps d'exécution, mais en temps de programmation  ;-)
        ' j'aurais besoin d'un trouver qui retourne le sly... en fonction d'un feature(edge)
        Dim p As SlyAretePoutre
        For Each p In lst_AretePoutre
            p.EffacerIntersection()
        Next


        '    ' fin de l'effacement des attributs
        Dim swArete1 As SldWorks.Edge
        Dim swArete2 As SldWorks.Edge
        Dim SlyArete1 As SlyAretePoutre
        Dim SlyArete2 As SlyAretePoutre

        Dim xyz() As Double = Nothing
        Dim pt As New InterPoutrePoutre
        'Dim lst_pts As New Collection
        Dim i As Long

        Dim a1 As Integer
        Dim a2 As Integer


        For a1 = 0 To lst_AretePoutre.Count - 2
            SlyArete1 = lst_AretePoutre.Item(a1)
            swArete1 = SlyArete1.swArete

            For a2 = a1 + 1 To lst_AretePoutre.Count - 1
                SlyArete2 = lst_AretePoutre.Item(a2)
                swArete2 = SlyArete2.swArete

                If DetectAreteArete(swArete1, swArete2, xyz) Then ' intersection courbe-courbe
                    For i = 0 To CInt((UBound(xyz) + 1) / 4) Step 4 ' si y'a plus d'une intersection entre les 2 courbes...

                        Select Case isIntersect_milieu(swArete1, swArete2, xyz(i), xyz(i + 1), xyz(i + 2))
                            Case 1 ' la première courbe est coupée
                                pt = New InterPoutrePoutre
                                pt.Arete = swArete1
                                pt.x = xyz(i) : pt.y = xyz(i + 1) : pt.z = xyz(i + 2)
                                SlyArete1.AjouterPointAPAP(pt)
                                pt = Nothing
                            Case 2 ' la seconde courbe est coupée
                                pt = New InterPoutrePoutre
                                pt.Arete = swArete2
                                pt.x = xyz(i) : pt.y = xyz(i + 1) : pt.z = xyz(i + 2)
                                SlyArete2.AjouterPointAPAP(pt)
                                pt = Nothing
                            Case 3 ' les doux courbes sont coupées
                                pt = New InterPoutrePoutre
                                pt.Arete = swArete1
                                pt.x = xyz(i) : pt.y = xyz(i + 1) : pt.z = xyz(i + 2)
                                SlyArete1.AjouterPointAPAP(pt)
                                pt = Nothing
                                pt = New InterPoutrePoutre
                                pt.Arete = swArete2
                                pt.x = xyz(i) : pt.y = xyz(i + 1) : pt.z = xyz(i + 2)
                                SlyArete2.AjouterPointAPAP(pt)
                                pt = Nothing
                        End Select
                    Next i
                End If

            Next a2
        Next a1
        pt = Nothing
    End Sub


    Public Function DetectAreteArete(ByRef swArete1 As SldWorks.Edge, ByRef swArete2 As SldWorks.Edge, ByRef xyz() As Double) As Boolean
        ' function qui détermine si 2 arêtes se touchent, si oui alors la fonction retourne vrai et le tableau XYZ contient le point d'intersection
        xyz = Nothing
        If swArete2 Is Nothing Then Exit Function
        If swArete1 Is Nothing Then Exit Function

        Dim swSommet As SldWorks.Vertex
        Dim point1 As Object  ' point début courbe 1
        Dim point2 As Object  ' point fin,  courbe 1
        Dim point3 As Object  ' point début courbe 2
        Dim point4 As Object  ' point fin,  courbe 2
        Dim swCourbe1 As SldWorks.Curve
        Dim swCourbe2 As SldWorks.Curve
        Dim vIntersectPts As Object

        swSommet = swArete1.GetStartVertex
        If swSommet Is Nothing Then ' cercle fermé
            point1 = swArete1.Evaluate(0)
            point2 = point1
        Else
            point1 = swSommet.GetPoint
            swSommet = swArete1.GetEndVertex
            point2 = swSommet.GetPoint
        End If

        swSommet = swArete2.GetStartVertex
        If swSommet Is Nothing Then
            point3 = swArete2.Evaluate(0)
            point4 = point3
        Else
            point3 = swSommet.GetPoint

            swSommet = swArete2.GetEndVertex
            point4 = swSommet.GetPoint
        End If
        swCourbe1 = swArete1.GetCurve
        swCourbe2 = swArete2.GetCurve


        vIntersectPts = swCourbe1.IntersectCurve(swCourbe2, point1, point2, point3, point4)

        Try  ' je déteste cette façon de procéder, mais il ne détecte pas que le vintersectpts est nothing, il le croit system.dbnull.  Et je ne peut comparer avec ça.
            xyz = vIntersectPts

            If Not xyz Is Nothing Then
                DetectAreteArete = True
                Dim i As Integer
                For i = 0 To ((UBound(xyz) + 1) / 4) - 1
                    Debug.Write("x" & i & " = " & xyz(i * 4) * 1000 & " mm")
                    Debug.Write("y" & i & " = " & xyz(i * 4) * 1000 & " mm")
                    Debug.Write("z" & i & " = " & xyz(i * 4) * 1000 & " mm")
                Next i
            Else
                Return False
            End If
        Catch ex As Exception
            DetectAreteArete = False
        End Try

    End Function

    Public Function DetectAreteArete(ByRef swArete1 As SldWorks.Edge, ByRef swCourbe2 As SldWorks.Curve, ByRef xyz() As Double) As Boolean
        ' function qui détermine si 2 arêtes se touchent, si oui alors la fonction retourne vrai et le tableau XYZ contient le point d'intersection
        xyz = Nothing
        If swArete1 Is Nothing Then Exit Function
        Dim swSommet As SldWorks.Vertex
        Dim point1 As Object  ' point début courbe 1
        Dim point2 As Object  ' point fin,  courbe 1
        Dim point3 As Object  ' point début courbe 2
        Dim point4 As Object  ' point fin,  courbe 2
        Dim swCourbe1 As SldWorks.Curve

        Dim vIntersectPts As Object

        swSommet = swArete1.GetStartVertex
        If swSommet Is Nothing Then ' cercle fermé
            point1 = swArete1.Evaluate(0)
            point2 = point1
        Else
            point1 = swSommet.GetPoint
            swSommet = swArete1.GetEndVertex
            point2 = swSommet.GetPoint
        End If


        If Not swCourbe2.IsTrimmedCurve Then
            Dim vP1 As Object, vP2 As Object
            vP1 = swCourbe2.Evaluate(-9000)
            vP2 = swCourbe2.Evaluate(9000)
            swCourbe2 = swCourbe2.CreateTrimmedCurve2(vP1(0), vP1(1), vP1(2), vP2(0), vP2(1), vP2(2))
        End If

        If Not swCourbe2.IsTrimmedCurve Then
            MsgBox("Ça plante, le trimmage n'a pas fonctionné...")
        End If

        Dim dpoint3 As Double, dpoint4 As Double, isClosed As Boolean, isPeriodic As Boolean
        swCourbe2.GetEndParams(dpoint3, dpoint4, isClosed, isPeriodic)

        point3 = swCourbe2.Evaluate(dpoint3)
        point4 = swCourbe2.Evaluate(dpoint4)

        swCourbe1 = swArete1.GetCurve

        vIntersectPts = swCourbe1.IntersectCurve(swCourbe2, point1, point2, point3, point4)

        Try  ' je déteste cette façon de procéder, mais il ne détecte pas que le vintersectpts est nothing, il le croit system.dbnull.  Et je ne peut comparer avec ça.
            xyz = vIntersectPts

            If Not xyz Is Nothing Then
                DetectAreteArete = True
                Dim i As Integer
                For i = 0 To ((UBound(xyz) + 1) / 4) - 1
                    Debug.Write("x" & i & " = " & xyz(i * 4) * 1000 & " mm")
                    Debug.Write("y" & i & " = " & xyz(i * 4) * 1000 & " mm")
                    Debug.Write("z" & i & " = " & xyz(i * 4) * 1000 & " mm")
                Next i
            Else
                Return False
            End If
        Catch ex As Exception
            DetectAreteArete = False
        End Try

    End Function

    Private Function isIntersect_milieu(ByRef Arete1 As SldWorks.Edge, ByRef Arete2 As SldWorks.Edge, ByRef x As Double, ByRef y As Double, ByRef z As Double) As Byte
        'procédure qui dit si les courbes s'intersectent au milieu ou à leurs point d'extrémité.
        ' 0 si les deux intersections sont sur des points
        ' 1 si la première courbe est coupée
        ' 2 si la deuxième est coupée
        ' 3 si les deux courbes sont coupées


        Dim retv As Object  ' mettre les valeurs de T dans la classe

        retv = Arete1.GetCurveParams2()  ' retv 0,1,2 point de départ, 3,4,5 point final
        If (Math.Abs(retv(0) - x) < Epsilon) And (Math.Abs(retv(1) - y) < Epsilon) And (Math.Abs(retv(2) - z) < Epsilon) Or (Math.Abs(retv(3) - x) < Epsilon) And (Math.Abs(retv(4) - y) < Epsilon) And (Math.Abs(retv(5) - z) < Epsilon) Then isIntersect_milieu = 0 Else isIntersect_milieu = 1

        retv = Arete2.GetCurveParams2()  ' retv 0,1,2 point de départ, 3,4,5 point final
        If (Math.Abs(retv(0) - x) < Epsilon) And (Math.Abs(retv(1) - y) < Epsilon) And (Math.Abs(retv(2) - z) < Epsilon) Or (Math.Abs(retv(3) - x) < Epsilon) And (Math.Abs(retv(4) - y) < Epsilon) And (Math.Abs(retv(5) - z) < Epsilon) Then isIntersect_milieu += 0 Else isIntersect_milieu += 2

    End Function


    ''' <summary>
    ''' Sub qui traite les intersections entre les Poutres et les coques
    ''' </summary>
    ''' <remarks></remarks>
    Private Sub DetectionPoutresCoques()
        '#1 on détecte les intersections, on en a 4 types, 
        ' #1a) intersection au milieu de la coque
        ' #1b) intersection sur une arête
        ' #1c) intersection sur un sommet
        ' #1d) intersection sur une certaine longueur
        ' faut également faire attention, il peut y avoir plus d'une poutre qui se rejoint à la même intersection.
        ' hypothèses: 1- les zones d'influence des poutres ne se croisent pas.
        '             2- une arête ne fait qu'une seule intersection, à un pount ou sur une ligne, mais pas à plusieurs endroits.
        '             3- une arête dont l'intersection est une courbe n'a pas d'autres intersections (pout l'instant)
        ' lorsqu'une intersection est détectée on créé une instance de la classe InterPoutreCoque 

        Dim sPoutre As SlyAretePoutre
        Dim sCoque As SlyFaceCoque
        Dim inter As InterPoutreCoque
        Dim xyz() As Double = Nothing
        Dim tipe As Byte
        Dim i As Integer
        Dim premier2 As Boolean


        For Each sCoque In lst_FaceCoque
            For Each sPoutre In lst_AretePoutre
                ' on cherche entre la coque et la poutre

                If DetectFaceArete(sPoutre.swArete, sCoque.SwFace, xyz) Then
                    For i = 0 To UBound(xyz) - 1 Step 3
                        ' trouver le tipe d'intersection...

                        Dim u() As Double
                        Dim v() As Double

                        u = sPoutre.GetOrientation(xyz(i + 0), xyz(i + 1), xyz(i + 2))
                        v = sCoque.GetNormale(xyz(i + 0), xyz(i + 1), xyz(i + 2))

                        If Math.Abs((Math.Abs(Math.Acos(Outils_Math.cosdir(u, v))) - Pi / 2)) < Epsilon Then
                            ' on est dans le même plan que la coque, on doit déterminer si on sort ou entre dans la coque.
                            Dim T As Double, T1 As Double, T2 As Double
                            Dim PointTest(2) As Double
                            T = sPoutre.GetT(xyz(i + 0), xyz(i + 1), xyz(i + 2))

                            ' on fait T + ou - 15 * epsilon.  si une valeur est dans la coque alors on considère que l'on a le type 2
                            ' on prend 15 fois epsilon car à 10, la fonction de solidworks considère que l'on est tellement près que l'on est sur la face
                            T1 = T + 15 * Epsilon
                            T2 = T - 15 * Epsilon

                            Dim effacer As Double

                            If sPoutre.Evaluer(T1, PointTest) Then
                                ' la valeur de T appartient à la poutre, maintenant on vérifie s'il appartient aussi à la coque
                                If Distance(sCoque.lst_Faces.Item(0), PointTest(0), PointTest(1), PointTest(2)) < Epsilon Then
                                    ' on est dans la coque.
                                    If Not premier2 Then tipe = 2 : premier2 = True Else tipe = 22
                                End If
                            End If

                            If sPoutre.Evaluer(T2, PointTest) Then
                                effacer = Distance(sCoque.lst_Faces.Item(0), PointTest(0), PointTest(1), PointTest(2))
                                If Distance(sCoque.lst_Faces.Item(0), PointTest(0), PointTest(1), PointTest(2)) < Epsilon Then
                                    If Not premier2 Then tipe = 2 : premier2 = True Else tipe = 22
                                End If
                            End If


                            ' si tipe n'est pas à 2 alors on est par défaut à 3
                            If tipe = 0 Then tipe = 3

                        Else
                            tipe = 1 ' on va devoir couper en X
                        End If

                        inter = sCoque.AjouterInterPoutre(sPoutre, xyz, tipe)
                        sPoutre.lst_InterCoque.Add(inter)
                        tipe = 0  ' reset

                    Next i ' autre point d'intersection

                End If
                premier2 = False ' reset
            Next sPoutre
        Next sCoque
    End Sub


    ''' <summary>
    ''' Sub qui procède au découpage des coques en fonction des poutres
    ''' </summary>
    ''' <remarks></remarks>
    Private Sub DecouperPoutreCoque()
        ' #2 on procède au découpage de la face
        Dim sCoque As SlyFaceCoque
        For Each sCoque In lst_FaceCoque
            If sCoque.lst_InterPoutre.Count > 0 Then
                'sVol.chercherAttributs()
                sCoque.decouper()

                ' on met-a-jour l'attribut des conditions aux limites
                Dim attr As SldWorks.Attribute
                Dim swent As SldWorks.Entity
                Dim nom3 As String = Nothing
                Dim p As SldWorks.Parameter
                If Not sCoque.AttributCL Is Nothing Then
                    nom3 = "CL_" & sCoque.nom
                    swent = sCoque.lst_Faces.Item(1)
                    attr = swent.FindAttribute(Intersections.DefAttrConditionLimite, 0) ' si l'attribut existe déjà on pointe dessus.

                    If attr Is Nothing Then attr = Intersections.DefAttrConditionLimite.CreateInstance5(swModel, sCoque.lst_Faces.Item(1), nom3, 0, 0)
                    p = attr.GetParameter("CL")
                    p.SetStringValue(sCoque.condition)

                End If
                GererDossiers("Conditions Aux Limites", nom3)
            End If

        Next
    End Sub


    ''' <summary>
    ''' Sub qui détecte les intersections entre les faces du volume et les poutres
    ''' </summary>
    ''' <remarks></remarks>
    Private Sub DetectionPoutresVolumes()
        Dim sPoutre As SlyAretePoutre

        For Each sPoutre In lst_AretePoutre
            If Not sPoutre.IsFaceDeSection Then
                Call Intersections.GestionPoutreNormaleAvecVolume(sPoutre)
            Else
                ' la poutre est flaggée pour prendre la face de section
                Call Intersections.GestionFace_De_section(sPoutre)
            End If
        Next sPoutre
    End Sub

    ''' <summary>
    ''' Sub qui prend les intersections entre les faces et les volumes et qui coupe les volumes
    ''' </summary>
    ''' <remarks></remarks>
    Private Sub DecouperPoutreVolume()
        ' #2 on procède au découpage de la face
        Dim sFaceVol As SlyFaceVolume

        For Each sFaceVol In lst_FaceVolume
            If sFaceVol.lst_InterPoutre.Count > 0 Then
                sFaceVol.decouper()

                ' on met-a-jour l'attribut des conditions aux limites
                Dim attr As sldworks.Attribute
                Dim swent As sldworks.Entity
                Dim nom3 As String = Nothing
                Dim p As sldworks.Parameter
                If Not sFaceVol.AttributCL Is Nothing Then
                    nom3 = "CL_" & sFaceVol.nom
                    swent = sFaceVol.SwFace
                    attr = swent.FindAttribute(Intersections.DefAttrConditionLimite, 0) ' si l'attribut existe déjà on pointe dessus.

                    If attr Is Nothing Then attr = Intersections.DefAttrConditionLimite.CreateInstance5(swModel, sFaceVol.SwFace, nom3, 0, 0)
                    p = attr.GetParameter("CL")
                    p.SetStringValue(sFaceVol.condition)

                End If
                GererDossiers("Conditions Aux Limites", nom3)
            End If

        Next
    End Sub


    ''' <summary>
    ''' Sub qui gère la poutre qui a été flaggée comme ayant une face de section.  Met à jour les attributs avec les bonnes inerties et autre
    ''' </summary>
    ''' <param name="sPoutre">La poutre avec laquelle il faut travailler</param>
    ''' <remarks></remarks>
    Private Sub GestionFace_De_section(ByRef sPoutre As SlyAretePoutre)

        Dim sVol As SlyFaceVolume
        Dim swExt As SldWorks.ModelDocExtension
        Dim xyz() As Double = Nothing, xyz2() As Double = Nothing
        Dim section As Object
        Dim Face1 As SlyFaceVolume = Nothing, prop1() As Double = Nothing
        Dim Face2 As SlyFaceVolume = Nothing, prop2() As Double = Nothing
        Dim proprietes() As Double
        Dim i As Integer
        swExt = swModel.Extension

        ' 1 - Spotter la ou les faces en question
        xyz = sPoutre.GetStartPoint
        xyz2 = sPoutre.GetEndPoint

        For Each sVol In Commun.lst_FaceVolume
            If Intersections.DetectSurfaceArete(sPoutre.swArete, sVol.SwFace, Nothing) Then
                swModel.ClearSelection2(True)
                sVol.Selectionner()
                section = swExt.GetSectionProperties2(sVol.SwFace) : proprietes = section
                swModel.ClearSelection2(True)
                ' la fonction getsectionproperties renvoie des valeurs dont la précision est très douteuse...
                ' on met alors une beaucoup plus grosse tolérance...
                If Math.Abs(xyz(0) - proprietes(2)) < 0.0001 And Math.Abs(xyz(1) - proprietes(3)) < 0.0001 And Math.Abs(xyz(2) - proprietes(4)) < 0.0001 Then
                    Face1 = sVol : prop1 = proprietes
                ElseIf Math.Abs(xyz2(0) - proprietes(2)) < 0.0001 And Math.Abs(xyz2(1) - proprietes(3)) < 0.0001 And Math.Abs(xyz2(2) - proprietes(4)) < 0.0001 Then
                    Face2 = sVol : prop2 = proprietes
                Else
                    ' on a une intersection où la poutre touche à une face «guide» et une autre partie touche à une face normale
                    ' la face est automatiquement une face de volume.  Si c'est une coque, elle est traitée ailleurs.

                    'Dim interNormale As InterPoutreVolume
                    'Dim point1 As Object = Nothing
                    'swModel.ClosestDistance(sPoutre.swArete, sVol.SwFace, point1, Nothing)
                    'Dim xyzNormale() As Double = point1
                    'interNormale = sVol.AjouterInterPoutre(sPoutre, xyzNormale, 6) ' y'a des cas où ça pourrait ne pas être 1...
                    'sPoutre.lst_InterCoque.Add(interNormale)

                End If

            End If
        Next

        ' si les 2 faces sont nothing, c'est que l'on a un sérieux problème...
        If Face1 Is Nothing AndAlso Face2 Is Nothing Then
            MsgBox("La poutre " & sPoutre.nom & " doit avoir au moins une face pour évaluer sa section.  Cette face n'a pas été trouvée." & vbCr & vbCr & " Vérifiez que: " & vbCr & "    - Un sommet de la poutre repose sur le centre de gravité d'une face" & vbCr & "    - Que l'arète de la poutre est perpendiculaire à la face" & vbCr & "    - Que la face est plane " & vbCr & "    - Je crois avoir tout couvert..." & vbCr & vbCr & "  La poutre a problème a été colorée en rouge.", MsgBoxStyle.Critical, "Impossible de trouver la face représentant la section")
            Dim e As New SuperArete(sPoutre.swArete, True)
            e.Colorer(3, 1, 0, 0)
            sPoutre.Selectionner()
            Err.Raise(520, "Gestion_Face_De_Section", "Une poutre ne peut être traitée car elle n'a pas été correctement définie.")
            Exit Sub
        End If


        ' 2 - s'assurer que:
        ' 2.1 les 2 faces sont identiques (si 2 faces)
        If Face1 IsNot Nothing And Face2 IsNot Nothing Then
            If Not Math.Abs(prop1(1) - prop2(1)) < 0.0000001 Then
                Dim chaine As String
                chaine = "Les 2 faces ne sont pas identiques... "
                For i = 0 To 14
                    chaine = chaine & vbCr & i & " - " & Format(prop1(i), "0.000E+00") & "  " & Format(prop2(i), "0.000E+00")
                Next i
                MsgBox(chaine)
            End If
        End If
        ' 2.2 le centroide est au point d'intersection
        ' déjà fait
        ' 2.3 la face est plane
        Dim prop() As Double = Nothing

        If Face1 IsNot Nothing Then
            prop = prop1
            If Not (Face1.estPlan Xor Face1.estFauxPlan(prop1(2), prop1(3), prop1(4))) Then
                MsgBox("La face 1 n'est pas plane!")
                Err.Raise(513, , "Ne peut pas prendre une face non plane comme source de la section de la poutre")
            End If
        End If
        If Face2 IsNot Nothing Then
            prop = prop2
            If Not (Face2.estPlan Xor Face2.estFauxPlan(prop2(2), prop2(3), prop2(4))) Then
                MsgBox("La face 2 n'est pas plane!")
                Err.Raise(513, , "Ne peut pas prendre une face non plane comme source de la section de la poutre")
            End If
        End If


        ' 2.4 les face sont perpendiculaire aux point d'intersection
        Dim u() As Double
        Dim v() As Double
        Dim angle As Double
        If Face1 IsNot Nothing Then
            u = Face1.GetNormaleSurface(prop1(2), prop1(3), prop1(4))
            v = sPoutre.GetOrientation(prop1(2), prop1(3), prop1(4))
            angle = Outils_Math.Angle2Vecteurs(u, v)
            If Not (Math.Abs(angle - Pi) < (100 * Epsilon)) And Not (angle < 0.005) Then
                MsgBox("La poutre n'est pas perpendiculaire à la face1")
                Err.Raise(514, "GestionFace_de_section - Dans intersection.vb", "Ne peut pas traiter une poutre dont la face de base n'est pas perpendiculaire")
            End If
        End If

        If Face2 IsNot Nothing Then
            u = Face2.GetNormaleSurface(prop2(2), prop2(3), prop2(4))
            v = sPoutre.GetOrientation(prop2(2), prop2(3), prop2(4))
            angle = Outils_Math.Angle2Vecteurs(u, v)
            If Not (Math.Abs(angle - Pi) < (100 * Epsilon)) And Not (angle < Epsilon) Then
                MsgBox("La poutre n'est pas perpendiculaire à la face2")
                Err.Raise(514, "GestionFace_de_section - Dans intersection.vb", "Ne peut pas traiter une poutre dont la face de base n'est pas perpendiculaire")
            End If
        End If

        ' 3 - trouver l'inertie et l'aire, placer le point 3 de façon cohérente, updater l'attribut.
        Dim Nom3 As String

        Commun.MettreUnPoint(prop(2) + prop(18) / 1000, prop(3) + prop(19) / 1000, prop(4) + prop(20) / 1000, True)
        Dim selmgr As sldworks.SelectionMgr = swModel.SelectionManager
        Nom3 = RealConstant.RCCode.Creation3iemePoint(selmgr.GetSelectedObject(1))
        sPoutre.SetAttributsDePoutre(False, Nom3, , , prop(13), prop(14), prop(1), , , , , , , 1)

        ' 4 créer une nouvelle instance de la classe interFacePoutre de tipe 5
        Dim inter As New InterPoutreVolume
        Dim inter2 As New InterPoutreVolume
        If Face1 IsNot Nothing Then
            xyz(0) = prop1(2) : xyz(1) = prop1(3) : xyz(2) = prop1(4)
            inter = Face1.AjouterInterPoutre(sPoutre, xyz, 5)
            sPoutre.lst_InterCoque.Add(inter)
        End If

        If Face2 IsNot Nothing Then
            xyz(0) = prop2(2) : xyz(1) = prop2(3) : xyz(2) = prop2(4)
            inter2 = Face2.AjouterInterPoutre(sPoutre, xyz, 5)
            sPoutre.lst_InterCoque.Add(inter2)
        End If

    End Sub


    ''' <summary>
    ''' Sub qui trouve les intersections entre les poutres et les faces de volumes et qui créé une instance de la classe d'intersection
    ''' </summary>
    ''' <param name="sPoutre">La poutre avec laquelle il faut trouver les intersections</param>
    ''' <remarks></remarks>
    Private Sub GestionPoutreNormaleAvecVolume(ByRef sPoutre As SlyAretePoutre)
        Dim inter As InterPoutreVolume
        Dim xyz() As Double = Nothing
        Dim tipe As Byte
        Dim i As Integer
        Dim premier2 As Boolean
        Dim sVol As SlyFaceVolume
        Dim SurSurface As Boolean = False

        For Each sVol In lst_FaceVolume

            ' on cherche entre la coque et la poutre
            SurSurface = False
            If DetectFaceArete(sPoutre.swArete, sVol, xyz, SurSurface) Then
                For i = 0 To UBound(xyz) - 1 Step 3
                    ' trouver le tipe d'intersection...

                    Dim u() As Double
                    Dim v() As Double

                    u = sPoutre.GetOrientation(xyz(i + 0), xyz(i + 1), xyz(i + 2))
                    v = sVol.GetNormale(xyz(i + 0), xyz(i + 1), xyz(i + 2))

                    ' ***** ici, différencier entre coque et volume ****
                    If Math.Abs((Math.Abs(Math.Acos(Outils_Math.cosdir(u, v))) - Pi / 2)) < Epsilon Then
                        ' on est dans le même plan que la coque, on doit déterminer si on sort ou entre dans la coque.
                        Dim T As Double, T1 As Double, T2 As Double
                        Dim PointTest(2) As Double
                        T = sPoutre.GetT(xyz(i + 0), xyz(i + 1), xyz(i + 2))

                        ' on fait T + ou - 15 * epsilon.  si une valeur est dans la coque alors on considère que l'on a le type 2
                        ' on prend 15 fois epxilon car à 10, la fonction de solidworks considère que l'on est tellement près que l'on est sur la face
                        T1 = T + 15 * Epsilon
                        T2 = T - 15 * Epsilon


                        If sPoutre.Evaluer(T1, PointTest) Then
                            ' la valeur de T appartient à la poutre, maintenant on vérifie s'il appartient aussi à la coque
                            If Distance(sVol.SwFace, PointTest(0), PointTest(1), PointTest(2)) < Epsilon Then
                                ' on est dans la coque.
                                If Not premier2 Then tipe = 2 : premier2 = True Else tipe = 22
                            End If
                        End If

                        If sPoutre.Evaluer(T2, PointTest) Then

                            If Distance(sVol.SwFace, PointTest(0), PointTest(1), PointTest(2)) < Epsilon Then
                                If Not premier2 Then tipe = 2 : premier2 = True Else tipe = 22
                            End If
                        End If

                        ' si tipe n'est pas à 2 alors on est par défaut à 3
                        If tipe = 0 Then tipe = 3

                    Else
                        tipe = 1 ' on va devoir couper en X
                    End If

                    inter = sVol.AjouterInterPoutre(sPoutre, xyz, tipe)
                    sPoutre.lst_InterCoque.Add(inter)
                    tipe = 0  ' reset

                Next i ' autre point d'intersection

            ElseIf SurSurface = True Then
                ' Si les 2 points de la courbe touchent à la surface alors  coupe-long
                If Commun.Distance(sVol.SwFace, sPoutre.GetStartPoint()) < Epsilon AndAlso Commun.Distance(sVol.SwFace, sPoutre.GetEndPoint()) < Epsilon Then
                    inter = sVol.AjouterInterPoutre(sPoutre, xyz, 2)
                    sPoutre.lst_InterCoque.Add(inter)
                End If
            End If
            premier2 = False ' reset

        Next sVol
    End Sub


    'Private Function DetectFaceArete(ByRef swCurve As SldWorks.Curve, ByRef swFace As SldWorks.Face2, ByRef xyz() As Double) As Boolean
    '    ' function qui détecte si une arête coupe une face, si c'est le cas la function retourne true et remplie le tableau xyz avec le point d'intersection
    '    Dim vCurveParam As Object
    '    Dim swSurf As SldWorks.Surface
    '    Dim vCurveBounds As Object, vPointArray As Object, vTArray As Object, vUVArray As Object
    '    Dim nCurveBounds(5) As Double
    '    Dim bRet As Boolean
    '    Dim i As Integer

    '    swSurf = swFace.GetSurface

    '    ' 2- on va chercher les paramètres de la spline avec Curve::ConvertLineToBcurve
    '    If swCurve.IsLine Then
    '        Dim startp(2) As Double
    '        Dim endp(2) As Double
    '        Dim vStart As Object
    '        Dim vEnd As Object

    '        vCurveParam = swCurve.LineParams()
    '        startp(0) = vCurveParam(0)
    '        startp(1) = vCurveParam(1)
    '        startp(2) = vCurveParam(2)
    '        vStart = startp

    '        endp(0) = vCurveParam(3)
    '        endp(1) = vCurveParam(4)
    '        endp(2) = vCurveParam(5)
    '        vEnd = endp

    '        Dim retval As Object
    '        Dim valeur As Double
    '        Dim dimension As Integer
    '        Dim ordre As Integer
    '        Dim nbpoints As Integer
    '        Dim periodique As Integer

    '        retval = swCurve.ConvertLineToBcurve(vStart, vEnd)

    '        ' on suppose que la droite est toujours transformée en spline non-rationelle de dimension 3 et d'ordre 2
    '        Dim knots(3) As Double
    '        Dim ctrlPoints(5) As Double

    '        knots(0) = retval(2)
    '        knots(1) = retval(3)
    '        knots(2) = retval(4)
    '        knots(3) = retval(5)

    '        ctrlPoints(0) = retval(6)
    '        ctrlPoints(1) = retval(7)
    '        ctrlPoints(2) = retval(8)
    '        ctrlPoints(3) = retval(9)
    '        ctrlPoints(4) = retval(10)
    '        ctrlPoints(5) = retval(11)

    '        '3 - on créé une spline dans le modeleur


    '        Dim modeler As SldWorks.Modeler
    '        modeler = swApp.GetModeler
    '        Dim props As Object
    '        Dim dProps(1) As Double
    '        Dim vKnots As Object, vCtrlPoints As Object

    '        dProps(0) = retval(0)
    '        dProps(1) = retval(1)
    '        props = dProps
    '        vKnots = knots
    '        vCtrlPoints = ctrlPoints

    '        swCurve = modeler.CreateBsplineCurve(props, vKnots, vCtrlPoints)
    '        ' 4 - on a une spline, on peut utiliser la fonction IntersectCurve
    '    End If


    '    For i = 0 To 5
    '        nCurveBounds(i) = vCurveParam(i)
    '    Next i

    '    vCurveBounds = nCurveBounds
    '    bRet = swSurf.IntersectCurve(swCurve, vCurveBounds, vPointArray, vTArray, vUVArray)

    '    Dim Pts() As Double
    '    Dim PointsTemp() As Double
    '    Dim j As Integer

    '    Pts = vPointArray

    '    If UBound(Pts) < 1 Then Return False
    '    ' ATTENTION, le point n'est pas nécessairement dans la face
    '    For i = 0 To UBound(Pts) - 1 Step 3
    '        If Commun.Distance(swFace, Pts(i + 0), Pts(i + 1), Pts(i + 2)) < Epsilon Then
    '            ReDim Preserve PointsTemp(j + 2)
    '            PointsTemp(j) = Pts(i + 0)
    '            PointsTemp(j + 1) = Pts(i + 1)
    '            PointsTemp(j + 2) = Pts(i + 2)
    '            j += 3
    '        End If
    '    Next

    '    xyz = PointsTemp
    '    If xyz Is Nothing Then Return False
    '    If UBound(xyz) > 0 Then Return True


    'End Function

    Private Function DetectFaceArete(ByRef swArete As SldWorks.Edge, ByRef swFace As SldWorks.Face2, ByRef xyz() As Double) As Boolean
        ' function qui détecte si une arête coupe une face, si c'est le cas la function retourne true et remplie le tableau xyz avec le point d'intersection

        Dim P1 As Object = Nothing, p2 As Object = Nothing
        If swModel.ClosestDistance(swArete, swFace, P1, p2) > Epsilon Then Return False

        Dim swCurve As SldWorks.Curve
        Dim swSurf As SldWorks.Surface
        Dim vCurveParam As Object

        Dim vCurveBounds As Object, vPointArray As Object = Nothing, vTArray As Object = Nothing, vUVArray As Object = Nothing
        Dim nCurveBounds(5) As Double
        Dim bRet As Boolean
        Dim i As Integer

        '1 - on va chercher la courbe et la surface
        swCurve = swArete.GetCurve
        swSurf = swFace.GetSurface


        ' 2- on va chercher les paramètres de la spline avec Curve::ConvertLineToBcurve
        vCurveParam = swArete.GetCurveParams2
        If swCurve.IsLine Then
            Dim startp(2) As Double
            Dim endp(2) As Double
            Dim vStart As Object
            Dim vEnd As Object


            startp(0) = vCurveParam(0)
            startp(1) = vCurveParam(1)
            startp(2) = vCurveParam(2)
            vStart = startp

            endp(0) = vCurveParam(3)
            endp(1) = vCurveParam(4)
            endp(2) = vCurveParam(5)
            vEnd = endp

            Dim retval As Object

            retval = swCurve.ConvertLineToBcurve(vStart, vEnd)

            ' on suppose que la droite est toujours transformée en spline non-rationelle de dimension 3 et d'ordre 2
            Dim knots(3) As Double
            Dim ctrlPoints(5) As Double

            knots(0) = retval(2)
            knots(1) = retval(3)
            knots(2) = retval(4)
            knots(3) = retval(5)

            ctrlPoints(0) = retval(6)
            ctrlPoints(1) = retval(7)
            ctrlPoints(2) = retval(8)
            ctrlPoints(3) = retval(9)
            ctrlPoints(4) = retval(10)
            ctrlPoints(5) = retval(11)

            '3 - on créé une spline dans le modeleur


            Dim modeler As SldWorks.Modeler
            modeler = swApp.GetModeler
            Dim props As Object
            Dim dProps(1) As Double
            Dim vKnots As Object, vCtrlPoints As Object

            dProps(0) = retval(0)
            dProps(1) = retval(1)
            props = dProps
            vKnots = knots
            vCtrlPoints = ctrlPoints

            swCurve = modeler.CreateBsplineCurve(props, vKnots, vCtrlPoints)
            ' 4 - on a une spline, on peut utiliser la fonction IntersectCurve
        End If


        For i = 0 To 5
            nCurveBounds(i) = vCurveParam(i)
        Next i

        vCurveBounds = nCurveBounds
        bRet = swSurf.IntersectCurve(swCurve, vCurveBounds, vPointArray, vTArray, vUVArray)

        Dim Pts() As Double
        Dim PointsTemp() As Double = Nothing
        Dim j As Integer

        Pts = vPointArray

        If UBound(Pts) < 1 Then Return False
        ' ATTENTION, le point n'est pas nécessairement dans la face
        For i = 0 To UBound(Pts) - 1 Step 3
            If Commun.Distance(swFace, Pts(i + 0), Pts(i + 1), Pts(i + 2)) < Epsilon Then
                ReDim Preserve PointsTemp(j + 2)
                PointsTemp(j) = Pts(i + 0)
                PointsTemp(j + 1) = Pts(i + 1)
                PointsTemp(j + 2) = Pts(i + 2)
                j += 3
            End If
        Next

        xyz = PointsTemp
        If xyz Is Nothing Then Return False
        If UBound(xyz) > 0 Then Return True

    End Function


    Private Function DetectFaceArete(ByRef swArete As sldworks.Edge, ByRef Slyface As SlyFaceVolume, ByRef xyz() As Double, Optional ByRef SurSurface As Boolean = False) As Boolean
        Dim swface As sldworks.Face2

        For Each swface In Slyface.lst_Faces
            If DetectFaceArete(swArete, swface, xyz) Then
                Dim vEdges As Object
                Dim Arete As sldworks.Edge
                Dim inu() As Double = Nothing

                vEdges = swface.GetEdges

                For Each Arete In vEdges
                    If DetectAreteArete(Arete, swArete, inu) Then
                        ' on a un type d'intersection entre une face et une arète...
                        Dim swSurf As sldworks.Surface
                        Dim retval As Object
                        Dim v(2) As Double
                        Dim u(2) As Double
                        Dim angle As Double
                        swSurf = swface.GetSurface
                        retval = swSurf.EvaluateAtPoint(xyz(0), xyz(1), xyz(2))
                        ' les 3 premiers de retval sont la normale...
                        v(0) = retval(0) : v(1) = retval(1) : v(2) = retval(2)

                        retval = swArete.GetParameter(inu(0), inu(1), inu(2))
                        retval = swArete.Evaluate(retval(0))
                        u(0) = retval(3) : u(1) = retval(4) : u(2) = retval(5)
                        angle = Outils_Math.Angle2Vecteurs(u, v)
                        'MsgBox(angle * 180 / Pi & "<-- angle sens --> " & swface.FaceInSurfaceSense & vbCr & "Aire = " & Slyface.Aire)


                        'Si l'angle entre la normale de la face et la normale de la courbe est de 0
                        ' ou de pi/2 alors on a une SurSurface
                        If Math.Abs(angle - Pi / 2) < 100 * Epsilon Or Math.Abs(angle) < 100 * Epsilon Then
                            SurSurface = True
                        Else
                            Return True
                        End If

                        'If ((angle < (Pi / 2)) Xor swface.FaceInSurfaceSense()) Then
                        '    Return True
                        'Else
                        '    SurSurface = True
                        'End If
                    End If
                Next


                If Not SurSurface Then Return True Else Return False
            End If
        Next
    End Function


    Private Function DetectSurfaceArete(ByRef swArete As SldWorks.Edge, ByRef swFace As SldWorks.Face2, ByRef xyz() As Double) As Boolean
        ' function qui détecte si une arête coupe une face, si c'est le cas la function retourne true et remplie le tableau xyz avec le point d'intersection
        Dim swCurve As SldWorks.Curve
        Dim swSurf As SldWorks.Surface
        Dim vCurveParam As Object

        Dim vCurveBounds As Object, vPointArray As Object = Nothing, vTArray As Object = Nothing, vUVArray As Object = Nothing
        Dim nCurveBounds(5) As Double
        Dim bRet As Boolean
        Dim i As Integer

        '1 - on va chercher la courbe et la surface
        swCurve = swArete.GetCurve
        swSurf = swFace.GetSurface


        ' 2- on va chercher les paramètres de la spline avec Curve::ConvertLineToBcurve
        vCurveParam = swArete.GetCurveParams2
        If swCurve.IsLine Then
            Dim startp(2) As Double
            Dim endp(2) As Double
            Dim vStart As Object
            Dim vEnd As Object


            startp(0) = vCurveParam(0)
            startp(1) = vCurveParam(1)
            startp(2) = vCurveParam(2)
            vStart = startp

            endp(0) = vCurveParam(3)
            endp(1) = vCurveParam(4)
            endp(2) = vCurveParam(5)
            vEnd = endp

            Dim retval As Object

            retval = swCurve.ConvertLineToBcurve(vStart, vEnd)

            ' on suppose que la droite est toujours transformée en spline non-rationelle de dimension 3 et d'ordre 2
            Dim knots(3) As Double
            Dim ctrlPoints(5) As Double

            knots(0) = retval(2)
            knots(1) = retval(3)
            knots(2) = retval(4)
            knots(3) = retval(5)

            ctrlPoints(0) = retval(6)
            ctrlPoints(1) = retval(7)
            ctrlPoints(2) = retval(8)
            ctrlPoints(3) = retval(9)
            ctrlPoints(4) = retval(10)
            ctrlPoints(5) = retval(11)

            '3 - on créé une spline dans le modeleur


            Dim modeler As SldWorks.Modeler
            modeler = swApp.GetModeler
            Dim props As Object
            Dim dProps(1) As Double
            Dim vKnots As Object, vCtrlPoints As Object

            dProps(0) = retval(0)
            dProps(1) = retval(1)
            props = dProps
            vKnots = knots
            vCtrlPoints = ctrlPoints

            swCurve = modeler.CreateBsplineCurve(props, vKnots, vCtrlPoints)
            ' 4 - on a une spline, on peut utiliser la fonction IntersectCurve
        End If


        For i = 0 To 5
            nCurveBounds(i) = vCurveParam(i)
        Next i

        vCurveBounds = nCurveBounds
        bRet = swSurf.IntersectCurve(swCurve, vCurveBounds, vPointArray, vTArray, vUVArray)

        Dim Pts() As Double
        Dim PointsTemp() As Double = Nothing
        Dim j As Integer

        Pts = vPointArray

        If UBound(Pts) < 1 Then Return False

        For i = 0 To UBound(Pts) - 1 Step 3
            ReDim Preserve PointsTemp(j + 2)
            PointsTemp(j) = Pts(i + 0)
            PointsTemp(j + 1) = Pts(i + 1)
            PointsTemp(j + 2) = Pts(i + 2)
            j += 3
        Next i

        xyz = PointsTemp
        If xyz Is Nothing Then Return False
        If UBound(xyz) > 0 Then Return True


    End Function


    Private Function DetectSommetArete(ByRef swsommet As SldWorks.Vertex, ByRef swArete As SldWorks.Edge, ByRef xyz() As Double) As Boolean
        Dim x As Double, y As Double, z As Double
        Dim vPoint As Object
        Dim vPoint2 As Object
        vPoint = swsommet.GetPoint

        x = vPoint(0)
        y = vPoint(1)
        z = vPoint(2)

        vPoint2 = swArete.GetClosestPointOn(x, y, z)
        If Math.Abs(vPoint2(0) - x) < Epsilon And Math.Abs(vPoint2(1) - y) < Epsilon And Math.Abs(vPoint2(2) - z) < Epsilon Then
            ReDim xyz(2)
            xyz(0) = x
            xyz(1) = y
            xyz(2) = z
            DetectSommetArete = True
        End If

    End Function


    Private Function DetectSommetArete(ByRef x As Double, ByRef y As Double, ByRef z As Double, ByRef swArete As SldWorks.Edge) As Byte
        Dim vPoint As Object
        Dim vPoint2 As Object

        vPoint2 = swArete.GetClosestPointOn(x, y, z)
        If Math.Abs(vPoint2(0) - x) < Epsilon And Math.Abs(vPoint2(1) - y) < Epsilon And Math.Abs(vPoint2(2) - z) < Epsilon Then
            DetectSommetArete = 1

            ' maintenant on cherche à savoir si ca touche au premier sommet de l'arrête. 
            vPoint = swArete.GetCurveParams2()

            If Math.Abs(vPoint(0) - x) < Epsilon And Math.Abs(vPoint(1) - y) < Epsilon And Math.Abs(vPoint(2) - z) < Epsilon Then DetectSommetArete = 2
            If Math.Abs(vPoint(3) - x) < Epsilon And Math.Abs(vPoint(4) - y) < Epsilon And Math.Abs(vPoint(5) - z) < Epsilon Then DetectSommetArete = 2

        End If

    End Function

    Private Function DetectSommetArete(ByRef swArete2som As SldWorks.Edge, ByRef swArete As SldWorks.Edge, ByRef xyz() As Double, Optional ByRef ret As Byte = 0) As Boolean
        ' ret:  1 = sur premier sommet, 2 = sur dernier sommet
        Dim x As Double, y As Double, z As Double
        Dim vPoint As Object
        Dim vPoint2 As Object
        Dim swSommet As SldWorks.Vertex

        swSommet = swArete2som.GetStartVertex
        If IsNothing(swSommet) Then Exit Function ' on a un cercle...


        vPoint = swSommet.GetPoint

        x = vPoint(0)
        y = vPoint(1)
        z = vPoint(2)

        vPoint2 = swArete.GetClosestPointOn(x, y, z)
        If Math.Abs(vPoint2(0) - x) < Epsilon And Math.Abs(vPoint2(1) - y) < Epsilon And Math.Abs(vPoint2(2) - z) < Epsilon Then
            ReDim xyz(2)
            xyz(0) = x
            xyz(1) = y
            xyz(2) = z
            DetectSommetArete = True
            ret += 1
        End If

        swSommet = swArete2som.GetEndVertex
        vPoint = swSommet.GetPoint

        x = vPoint(0)
        y = vPoint(1)
        z = vPoint(2)

        vPoint2 = swArete.GetClosestPointOn(x, y, z)
        If Math.Abs(vPoint2(0) - x) < Epsilon And Math.Abs(vPoint2(1) - y) < Epsilon And Math.Abs(vPoint2(2) - z) < Epsilon Then
            ReDim xyz(2)
            xyz(0) = x
            xyz(1) = y
            xyz(2) = z
            DetectSommetArete = True
            ret += 2
        End If


    End Function

    Private Sub Memoriser3iemePoint()
        Dim poutre As SlyAretePoutre
        Dim attr As SldWorks.Attribute
        Dim p As SldWorks.Parameter
        Dim swEnt As SldWorks.Entity


        For Each poutre In Commun.lst_AretePoutre
            Dim xyz() As Double
            Dim nom As String = Nothing

            If Not poutre.IsFaceDeSection Then
                xyz = poutre.GetPoint3(nom)

                poutre.X3 = xyz(0)
                poutre.Y3 = xyz(1)
                poutre.Z3 = xyz(2)

                swEnt = poutre.swArete

                attr = swEnt.FindAttribute(Intersections.DefAttrRCP1, 0)
                p = attr.GetParameter("N3")

                p = attr.GetParameter("X3")
                p.SetDoubleValue2(poutre.X3, 2, "")

                p = attr.GetParameter("Y3")
                p.SetDoubleValue2(poutre.Y3, 2, "")

                p = attr.GetParameter("Z3")
                p.SetDoubleValue2(poutre.Z3, 2, "")
            End If
        Next poutre

    End Sub


    Private Sub CouperPoutres()

        ' procédure qui trouve les intersections entre les poutres et qui les coupe.
        ' à la fin on a plus d'intersections entre les poutres (ailleur qu'un sommet)

        Dim swArete1 As SldWorks.Edge
        Dim swArete2 As SldWorks.Edge
        Dim SlyArete1 As SlyAretePoutre
        Dim SlyArete2 As SlyAretePoutre

        Dim xyz() As Double = Nothing
        Dim pt As New InterPoutrePoutre
        'Dim lst_pts As New Collection
        Dim i As Long

        Dim a1 As Integer
        Dim a2 As Integer

        Try
            For a1 = 0 To lst_AretePoutre.Count - 1
                SlyArete1 = lst_AretePoutre.Item(a1)
                swArete1 = SlyArete1.swArete

                If Not a1 = lst_AretePoutre.Count - 1 Then
                    For a2 = a1 + 1 To lst_AretePoutre.Count - 1
                        SlyArete2 = lst_AretePoutre.Item(a2)
                        swArete2 = SlyArete2.swArete

                        If DetectAreteArete(swArete1, swArete2, xyz) Then ' intersection courbe-courbe
                            For i = 0 To CInt((UBound(xyz) + 1) / 4) Step 4 ' si y'a plus d'une intersection entre les 2 courbes...

                                Select Case isIntersect_milieu(swArete1, swArete2, xyz(i), xyz(i + 1), xyz(i + 2))
                                    Case 1 ' la première courbe est coupée
                                        pt = New InterPoutrePoutre
                                        pt.Arete = swArete1
                                        pt.x = xyz(i) : pt.y = xyz(i + 1) : pt.z = xyz(i + 2)
                                        SlyArete1.AjouterPointAPAP(pt)
                                        pt = Nothing
                                    Case 2 ' la seconde courbe est coupée
                                        pt = New InterPoutrePoutre
                                        pt.Arete = swArete2
                                        pt.x = xyz(i) : pt.y = xyz(i + 1) : pt.z = xyz(i + 2)
                                        SlyArete2.AjouterPointAPAP(pt)
                                        pt = Nothing
                                    Case 3 ' les doux poutres sont coupées
                                        pt = New InterPoutrePoutre
                                        pt.Arete = swArete1
                                        pt.x = xyz(i) : pt.y = xyz(i + 1) : pt.z = xyz(i + 2)
                                        SlyArete1.AjouterPointAPAP(pt)
                                        pt = Nothing
                                        pt = New InterPoutrePoutre
                                        pt.Arete = swArete2
                                        pt.x = xyz(i) : pt.y = xyz(i + 1) : pt.z = xyz(i + 2)
                                        SlyArete2.AjouterPointAPAP(pt)
                                        pt = Nothing
                                End Select
                            Next i
                        End If

                    Next a2
                End If

                Dim SlyArete3 As SlyAreteCoque
                For a2 = 0 To lst_AreteCoque.Count - 1
                    SlyArete3 = lst_AreteCoque.Item(a2)
                    swArete2 = SlyArete3.swArete

                    If DetectAreteArete(swArete1, swArete2, xyz) Then ' intersection courbe-courbe
                        For i = 0 To CInt((UBound(xyz) + 1) / 4) Step 4 ' si y'a plus d'une intersection entre les 2 courbes...

                            Select Case isIntersect_milieu(swArete1, swArete2, xyz(i), xyz(i + 1), xyz(i + 2))
                                Case 1 ' la première courbe est coupée
                                    pt = New InterPoutrePoutre
                                    pt.Arete = swArete1
                                    pt.x = xyz(i) : pt.y = xyz(i + 1) : pt.z = xyz(i + 2)
                                    SlyArete1.AjouterPointAPAP(pt)
                                    pt = Nothing
                                    ' seul la coque est coupée, elle sera découpée anyway.
                                Case 3 ' les doux poutres sont coupées, mais on note juste la deuxième
                                    pt = New InterPoutrePoutre
                                    pt.Arete = swArete1
                                    pt.x = xyz(i) : pt.y = xyz(i + 1) : pt.z = xyz(i + 2)
                                    SlyArete1.AjouterPointAPAP(pt)
                                    pt = Nothing
                            End Select
                        Next i
                    End If
                Next a2

                Dim slyarete4 As SlyAreteVol
                For a2 = 0 To lst_AreteVolume.Count - 1
                    slyarete4 = lst_AreteVolume.Item(a2)
                    swArete2 = slyarete4.swArete

                    If DetectAreteArete(swArete1, swArete2, xyz) Then ' intersection courbe-courbe
                        For i = 0 To CInt((UBound(xyz) + 1) / 4) Step 4 ' si y'a plus d'une intersection entre les 2 courbes...

                            Select Case isIntersect_milieu(swArete1, swArete2, xyz(i), xyz(i + 1), xyz(i + 2))
                                Case 1 ' la première courbe est coupée
                                    pt = New InterPoutrePoutre
                                    pt.Arete = swArete1
                                    pt.x = xyz(i) : pt.y = xyz(i + 1) : pt.z = xyz(i + 2)
                                    SlyArete1.AjouterPointAPAP(pt)
                                    pt = Nothing
                                Case 3 ' les doux poutres sont coupées, mais on note juste la deuxième
                                    pt = New InterPoutrePoutre
                                    pt.Arete = swArete1
                                    pt.x = xyz(i) : pt.y = xyz(i + 1) : pt.z = xyz(i + 2)
                                    SlyArete1.AjouterPointAPAP(pt)
                                    pt = Nothing
                            End Select
                        Next i
                    End If
                Next a2


            Next a1
            pt = Nothing

        Catch
            MsgBox(" La première boucle n'a pas marchée")
        End Try


        Dim count As Long = 0, count2 As Long = 0
        Try

            'toutes les poutres ont des points où elles doivent être coupées
            ' il suffit de couper.
            ' en réalité on suppress et on créé 2 ou plus courbes par dessus.
            Dim attr As sldworks.Attribute
            Dim swEnt As sldworks.Entity

            For Each SlyArete1 In lst_AretePoutre
                If SlyArete1.lst_PtsInterAPAP.Count > 0 Then ' on coupe
                    '1 - ordonner les points, de Tmin à Tmax et inclure les 2 extrémités de la poutre

                    count += 1
                    ' si c'est une droite
                    If SlyArete1.IsLine Then
                        Dim Pts(,) As Double
                        Dim swSketch As sldworks.Sketch
                        ReDim Pts(3, SlyArete1.lst_PtsInterAPAP.Count + 1)

                        SlyArete1.GetStartPoint(Pts(0, 0), Pts(1, 0), Pts(2, 0)) ' le premier point
                        Pts(3, 0) = SlyArete1.GetT(Pts(0, 0), Pts(1, 0), Pts(2, 0))

                        For i = 1 To SlyArete1.lst_PtsInterAPAP.Count
                            Pts(0, i) = SlyArete1.lst_PtsInterAPAP.Item(i).x
                            Pts(1, i) = SlyArete1.lst_PtsInterAPAP.Item(i).y
                            Pts(2, i) = SlyArete1.lst_PtsInterAPAP.Item(i).z
                            Pts(3, i) = SlyArete1.GetT(Pts(0, i), Pts(1, i), Pts(2, i))
                        Next i

                        Dim max As Integer = SlyArete1.lst_PtsInterAPAP.Count + 1
                        SlyArete1.GetEndPoint(Pts(0, max), Pts(1, max), Pts(2, max)) ' le dernier point
                        Pts(3, max) = SlyArete1.GetT(Pts(0, max), Pts(1, max), Pts(2, max))

                        ' faut ordonner les points selon T...
                        Dim j As Integer
                        Dim T1 As Double, T2 As Double, T3 As Double, T0 As Double
                        For i = 0 To max - 2
                            For j = 0 To max - i - 1
                                If Pts(3, j) > Pts(3, j + 1) Then
                                    T0 = Pts(0, j) : T1 = Pts(1, j) : T2 = Pts(2, j) : T3 = Pts(3, j)
                                    Pts(0, j) = Pts(0, j + 1) : Pts(1, j) = Pts(1, j + 1) : Pts(2, j) = Pts(2, j + 1) : Pts(3, j) = Pts(3, j + 1)
                                    Pts(0, j + 1) = T0 : Pts(1, j + 1) = T1 : Pts(2, j + 1) = T2 : Pts(3, j + 1) = T3
                                End If
                            Next j
                        Next i


                        For i = 0 To UBound(Pts, 2) - 1
                            swModel.Insert3DSketch2(False)
                            swModel.CreateLine2(Pts(0, i), Pts(1, i), Pts(2, i), Pts(0, i + 1), Pts(1, i + 1), Pts(2, i + 1)) ' et pour chaque segment
                            swSketch = swModel.GetActiveSketch2
                            swModel.Insert3DSketch2(False)

                            swEnt = swSketch : swEnt.Select2(False, 1) ' doit avoir un mark de 1
                            swModel.InsertCompositeCurve()
                            Try
                                count2 += 1
                                UpdateAttributs(SlyArete1, i) 'ajouter les attributs de la vieille poutre sur la nouvelle
                            Catch
                                MsgBox("UpdateAttribut n'a pas marché au compte # " & count2)
                            End Try
                        Next i


                    Else ' If SlyArete1.IsCircle Then  ' si c'est un cercle

                        Dim Pts(,) As Double
                        Dim swSketch As sldworks.Sketch
                        ReDim Pts(3, SlyArete1.lst_PtsInterAPAP.Count + 1)

                        SlyArete1.GetStartPoint(Pts(0, 0), Pts(1, 0), Pts(2, 0)) ' le premier point
                        Pts(3, 0) = SlyArete1.GetT(Pts(0, 0), Pts(1, 0), Pts(2, 0))

                        For i = 1 To SlyArete1.lst_PtsInterAPAP.Count
                            Pts(0, i) = SlyArete1.lst_PtsInterAPAP.Item(i).x
                            Pts(1, i) = SlyArete1.lst_PtsInterAPAP.Item(i).y
                            Pts(2, i) = SlyArete1.lst_PtsInterAPAP.Item(i).z
                            Pts(3, i) = SlyArete1.GetT(Pts(0, i), Pts(1, i), Pts(2, i))
                        Next i

                        Dim max As Integer = SlyArete1.lst_PtsInterAPAP.Count + 1
                        SlyArete1.GetEndPoint(Pts(0, max), Pts(1, max), Pts(2, max)) ' le dernier point
                        Pts(3, max) = SlyArete1.GetT(Pts(0, max), Pts(1, max), Pts(2, max))

                        ' faut ordonner les points selon T...
                        Dim j As Integer
                        Dim T1 As Double, T2 As Double, T3 As Double, T0 As Double
                        For i = 0 To max - 2
                            For j = 0 To max - i - 1
                                If Pts(3, j) > Pts(3, j + 1) Then
                                    T0 = Pts(0, j) : T1 = Pts(1, j) : T2 = Pts(2, j) : T3 = Pts(3, j)
                                    Pts(0, j) = Pts(0, j + 1) : Pts(1, j) = Pts(1, j + 1) : Pts(2, j) = Pts(2, j + 1) : Pts(3, j) = Pts(3, j + 1)
                                    Pts(0, j + 1) = T0 : Pts(1, j + 1) = T1 : Pts(2, j + 1) = T2 : Pts(3, j + 1) = T3
                                End If
                            Next j
                        Next i

                        Dim skSeg As sldworks.SketchSegment
                        Dim x As Double, y As Double, z As Double
                        Dim vretval As Object
                        Dim useEdge As sldworks.SketchSegment
                        Dim m As Integer


                        For i = 0 To UBound(Pts, 2) - 1
                            swModel.Insert3DSketch2(False)
                            ' sélectionner la edge originale
                            swEnt = SlyArete1.swArete
                            swEnt.Select4(False, Nothing)
                            swModel.SketchUseEdge2(False)
                            swSketch = swModel.GetActiveSketch2()

                            ' on créé 2 lignes de construction et on pick de chaque coté... mais on ne le fait pas si on est au premier ou au dernier segment. là on fait juste un pick.
                            If i <> 0 Then ' premier pick, élimine ce qui est avant.
                                skSeg = swModel.CreateLine2(Pts(0, i), Pts(1, i), Pts(2, i), 0.01, 0.01, 0.01)    'pts(0, i - 1) + 10000000 * Epsilon, pts(1, i - 1) + 100000 * Epsilon, pts(2, i - 1) + 100000 * Epsilon)
                                skSeg.ConstructionGeometry = True ' ligne de construction
                                swModel.ClearSelection2(True)
                                ' x et y doit être sélectionné sur un point avec un T inférieur au point d'intersection
                                SlyArete1.Evaluer((Pts(3, i - 1) + Pts(3, i)) / 2, x, y, z)
                                vretval = swSketch.GetSketchSegments
                                useEdge = vretval(0) : m = 0
                                Do Until useEdge.ConstructionGeometry = False
                                    m += 1
                                    useEdge = vretval(m)
                                Loop
                                useEdge.Select4(False, Nothing)
                                swModel.SketchTrim(1, 0, x, y) ' option = 1 pour trim, selEnd est pas utilisé ?, puis un point x et Y pour sélectionner.  et y'a pas de Z????? c'est un sketch3D!!!!!
                                skSeg = swModel.CreateLine2(0, 0, 0, x, y, 0)
                                skSeg.ConstructionGeometry = True
                            End If

                            If i <> UBound(Pts, 2) - 1 Then 'Second(pick)
                                skSeg = swModel.CreateLine2(Pts(0, i + 1), Pts(1, i + 1), Pts(2, i + 1), 0.05, 0, 0.01) 'pts(0, i + 1) + 10000000 * Epsilon, pts(1, i + 1) + 1000000 * Epsilon, pts(2, i + 1) + 100000 * Epsilon)
                                skSeg.ConstructionGeometry = True ' ligne de construction
                                swModel.ClearSelection2(True)

                                ' x et y doit être sélectionné sur un point avec un T inférieur au point d'intersection
                                SlyArete1.Evaluer((Pts(3, i + 1) + Pts(3, i + 2)) / 2, x, y, z)
                                vretval = swSketch.GetSketchSegments
                                useEdge = vretval(0) : m = 0
                                Do Until useEdge.ConstructionGeometry = False
                                    m += 1
                                    useEdge = vretval(m)
                                Loop

                                useEdge.Select4(False, Nothing)
                                swModel.SketchTrim(1, 0, x, y) ' option = 1 pour trim, selEnd est pas utilisé ?, puis un point x et Y pour sélectionner.  et y'a pas de Z????? c'est un sketch3D!!!!!
                                skSeg = swModel.CreateLine2(0, 0.02, 0, x, y, 0)
                                skSeg.ConstructionGeometry = True
                            End If

                            swModel.Insert3DSketch2(False)
                            swEnt = swSketch : swEnt.Select2(False, 1) ' doit avoir un mark de 1
                            swModel.InsertCompositeCurve()
                            UpdateAttributs(SlyArete1, i)
                        Next i

                    End If


                    ' on met un attribut pour ignorer l'arète.  Les sommets devraient donc aussi être ignorés.

                    SlyArete1.MettreAttributIgnorer()

                    'Dim nom As String
                    'Dim no As Integer
                    'Dim arete As sldworks.Edge
                    'arete = SlyArete1.swArete
                    'swEnt = arete
                    'nom = "Ignorer" & SlyArete1.nom & "_" & CStr(no)
                    'attr = swEnt.FindAttribute(Intersections.DefAttrRCP1, 0)
                    ''attr = DefAttrRCP1.CreateInstance5(swModel, arete, nom, 0, 2) ' une deuxième instance du RCPoutre...
                    'If attr Is Nothing Then
                    '    Commun.ColorerAretes()
                    '    swEnt = SlyArete1.swArete
                    '    attr = swEnt.FindAttribute(Intersections.DefAttrRCP1, 0)
                    'End If

                    'Dim p As sldworks.Parameter
                    'p = attr.GetParameter("D1")
                    'p.SetDoubleValue(-9)
                    'p = attr.GetParameter("D2")
                    'p.SetDoubleValue(-9)
                    'p = attr.GetParameter("D3")
                    'p.SetDoubleValue(-9)
                    'p = attr.GetParameter("D4")
                    'p.SetDoubleValue(-9)

                    'If attr Is Nothing Then MsgBox("Pas marché")


                End If
            Next SlyArete1
        Catch
            MsgBox("La seconde boucle n'a pas marchée  Compte: " & count)
        End Try

    End Sub


    Private Sub UpdateAttributs(ByRef slyarete1 As SlyAretePoutre, ByVal i As Integer)
        Dim newArete As sldworks.Edge
        Dim refcurve As sldworks.ReferenceCurve
        Dim attr As sldworks.Attribute = Nothing
        Dim swfeat As sldworks.Feature
        swModel.EditRebuild3()
        swfeat = swModel.FeatureByPositionReverse(0)

        Debug.Print(swfeat.Name & "<-- Nom  Typename--> " & swfeat.GetTypeName)

        refcurve = swfeat.GetSpecificFeature2()
        newArete = refcurve.GetFirstSegment()

        Dim ParamM As sldworks.Parameter
        Dim ParamS As sldworks.Parameter
        Dim ParamI1 As sldworks.Parameter
        Dim ParamI2 As sldworks.Parameter
        Dim ParamD1 As sldworks.Parameter
        Dim ParamD2 As sldworks.Parameter
        Dim ParamD3 As sldworks.Parameter
        Dim ParamD4 As sldworks.Parameter
        Dim ParamD5 As sldworks.Parameter
        Dim ParamD6 As sldworks.Parameter
        Dim ParamAs As sldworks.Parameter
        Dim ParamN3 As sldworks.Parameter

        Do While attr Is Nothing
            attr = Intersections.DefAttrRCP1.CreateInstance5(swModel, newArete, "Nouveau" & i & slyarete1.nom, 0, 2)
            i += 1
        Loop

        ParamM = attr.GetParameter("M")
        ParamS = attr.GetParameter("S")
        ParamI1 = attr.GetParameter("I1")
        ParamI2 = attr.GetParameter("I2")
        ParamD1 = attr.GetParameter("D1")
        ParamD2 = attr.GetParameter("D2")
        ParamD3 = attr.GetParameter("D3")
        ParamD4 = attr.GetParameter("D4")
        ParamD5 = attr.GetParameter("D5")
        ParamD6 = attr.GetParameter("D6")
        ParamAs = attr.GetParameter("As")

        ParamM.SetStringValue2(slyarete1.GetM, 2, "")  ' swAllConfiguration = 2
        ParamS.SetStringValue2(slyarete1.GetSection, 2, "")
        ParamI1.SetDoubleValue2(slyarete1.GetInertieXX, 2, "")
        ParamI2.SetDoubleValue2(slyarete1.GetInertieYY, 2, "")
        ParamD1.SetDoubleValue2(slyarete1.GetD1, 2, "")
        ParamD2.SetDoubleValue2(slyarete1.GetD2, 2, "")
        ParamD3.SetDoubleValue2(slyarete1.GetD3, 2, "")
        ParamD4.SetDoubleValue2(slyarete1.GetD4, 2, "")
        ParamD5.SetDoubleValue2(slyarete1.GetD5, 2, "")
        ParamD6.SetDoubleValue2(slyarete1.GetD6, 2, "")
        ParamAs.SetDoubleValue2(slyarete1.GetAireSection, 2, "")

        Dim p As sldworks.Parameter
        p = attr.GetParameter("N3")
        p.SetStringValue(slyarete1.GetN3)
        p = attr.GetParameter("X3")
        p.SetDoubleValue2(slyarete1.X3, 2, "")
        p = attr.GetParameter("Y3")
        p.SetDoubleValue2(slyarete1.Y3, 2, "")
        p = attr.GetParameter("Z3")
        p.SetDoubleValue2(slyarete1.Z3, 2, "")

        ' Commun.GererDossiers("Poutres", "Nouveau" & i & slyarete1.nom)

    End Sub


    ''' <summary>
    ''' sub qui gère la détection des intersections entre les coques et les volumes.
    ''' </summary>
    ''' <remarks></remarks>
    Private Sub DetectionCoqueVolume()

        ' détection de coque-coque
        Dim Face1 As SlyFaceVolume
        Dim coque2 As SlyFaceCoque
        Dim sketch As SldWorks.Sketch = Nothing
        Dim interFF As InterCoqueVolume

        For Each Face1 In Commun.lst_FaceVolume
            For Each coque2 In lst_FaceCoque

                If DetectFaceFace(coque2, Face1, True, sketch) Then
                    ' création de l'instance de interFace-face entre coque et volume

                    Dim vSeg As Object = sketch.GetSketchSegments()
                    If vSeg Is Nothing Then
                        DetectFaceFace(coque2, Face1, True, sketch)
                        vSeg = sketch.GetSketchSegments()
                        If vSeg Is Nothing Then MsgBox("Problème") : Err.Raise(555)
                    End If


                    interFF = New InterCoqueVolume
                    interFF.sFaceVolume = Face1
                    interFF.sFaceCoque = coque2
                    If coque2.PossedeFaceDeSection = True Then
                        interFF.FaceDeSection = True
                        Verifier_Coque_Section(interFF)
                    Else
                        interFF.FaceDeSection = False
                    End If
                    interFF.sketch = sketch
                    coque2.lst_InterCoqueVolume.Add(interFF)
                End If

            Next coque2
        Next Face1


    End Sub

    ''' <summary>
    ''' Sub qui fait les 6 vérifications pour le type d'intersection où la coque est partiellement dessinée en 3D.
    ''' En profite également pour remettre à jour la proriété de l'épaisseur de la coque 
    ''' </summary>
    ''' <param name="interFF">L'intersection</param>
    ''' <remarks>On fait au total 6 vérifications</remarks>
    Private Sub Verifier_Coque_Section(ByRef interFF As InterCoqueVolume)
        Dim sFaceCoque As SlyFaceCoque
        Dim sFaceVol As SlyFaceVolume
        Dim swFace1 As SldWorks.Face2
        Dim swent As SldWorks.Entity

        sFaceCoque = interFF.sFaceCoque
        sFaceVol = interFF.sFaceVolume

        ' 1 - vérifier que la face du volume est plane
        If Not sFaceVol.estPlan Then
            Err.Raise(513, "Verifier_Coque_Section", "La face représentant la section de la coque n'est pas plane...")
        End If

        ' 1.5 - La face du volume a 4 courbes, 2 petites de même longueur et 2 longues
        ' 1.6 ou qu'elle est une courbe fermée avec 2 loop d'une seule arète chaque.
        Dim vEdges As Object
        Dim Aretes() As SldWorks.Edge
        Dim longueur(3) As Double
        Dim i As Integer
        Dim a As SldWorks.Edge
        ReDim Aretes(0)

        swFace1 = sFaceVol.SwFace

        vEdges = swFace1.GetEdges()
        For Each a In vEdges
            ReDim Preserve Aretes(i)
            Aretes(i) = a
            i += 1
        Next

        If UBound(Aretes) = 3 Then
            For i = 0 To 3
                longueur(i) = Commun.GetLongueurArete(Aretes(i))
            Next i
            longueur = Ordonner(longueur, Aretes)
            If Not Math.Abs(longueur(0) - longueur(1)) < Epsilon Then
                swent = swFace1 : swent.Select2(False, 0)
                Err.Raise(514, "Verifier_Coque_Section", "La face du volume représentant la coque n'a pas 2 petites arètes de la même longueur, ce n'est pas normal.   La face a problème a été sélectionée")
            End If

        ElseIf UBound(Aretes) = 1 Then ' la coque représente un cylindre (ou un autre truc fermé)
            If Not swFace1.GetLoopCount = 2 Then
                swent = swFace1 : swent.Select2(False, 0)
                Err.Raise(514, "Verifier_Coque_Section", "La face du volume représentant la coque n'a pas 2 ou 4 arêtes ou n'est pas correcte...  La face a problème a été sélectionée")
            End If
        Else
            swent = swFace1 : swent.Select2(False, 0)
            Err.Raise(514, "Verifier_Coque_Section", "La face du volume représentant la coque n'a pas 2 ou 4 arêtes.  La face a problème a été sélectionée")
        End If

        ' 2 - Vérifier que la ligne est plane, mais ça je ne sais pas comment

        ' 3 - Perpendiculaire
        ' ok, on va chercher le midpoint de la ligne et on demande la normale à cette position des 2 faces.  Puis on demande(l) 'angle entre les 2 normales...

        Dim swArete As SldWorks.Edge
        Dim vmid As Object
        Dim mid() As Double = Nothing
        Dim T As Double
        Dim u() As Double, v() As Double
        Dim angle As Double

        If interFF.AreteCoque Is Nothing Then interFF.QuelleAreteCoqueToucheVolume() : Dim effacer As New SuperArete(interFF.AreteCoque, True) ': effacer.Colorer(2, 0.7, 0, 0.7)
        swArete = interFF.AreteCoque
        T = Commun.GetTMilieu(swArete)
        vmid = swArete.Evaluate(T)
        mid = vmid

        u = sFaceVol.GetNormale(mid(0), mid(1), mid(2))
        v = sFaceCoque.GetNormale(mid(0), mid(1), mid(2))

        angle = Outils_Math.Angle2Vecteurs(u, v)

        If Not Math.Abs(angle - Pi / 2) < (Epsilon * 10) Then
            sFaceVol.Selectionner(, False)
            sFaceCoque.Selectionner(, True)
            Err.Raise(515, "Verifier_Coque_Section", "La face de la coque n'est pas perpendiculaire à la face représentant sa section.  Les faces fautives ont été sélectionées.")
        End If

        ' 4 - vérifier que les 2 petites arètes sont coupées au milieu


        ' maintenant on update l'attribut de la coque, par chance, j'ai les coordonnées du point milieu...
        Dim epaisseur As Double
        Dim distance1 As Double, distance2 As Double
        Dim point1(2) As Double, point2(2) As Double
        Dim p As Object = Nothing, p2 As Object = Nothing

        If UBound(Aretes) = 3 Then
            p = Aretes(2).GetClosestPointOn(mid(0), mid(1), mid(2))
            point1(0) = p(0) : point1(1) = p(1) : point1(2) = p(2)
            distance1 = Math.Sqrt((point1(0) - mid(0)) ^ 2 + (point1(1) - mid(1)) ^ 2 + (point1(2) - mid(2)) ^ 2)

            p = Aretes(3).GetClosestPointOn(mid(0), mid(1), mid(2))
            point2(0) = p(0) : point2(1) = p(1) : point2(2) = p(2)
            distance2 = Math.Sqrt((point2(0) - mid(0)) ^ 2 + (point2(1) - mid(1)) ^ 2 + (point2(2) - mid(2)) ^ 2)

            If Not Math.Abs(distance1 - distance2) < Epsilon Then
                Err.Raise(516, "Verifier_Coque_Section", "La distance n'est pas la même de chaque coté de la coque")
            Else
                epaisseur = distance1 * 2
            End If

        ElseIf UBound(Aretes) = 1 Then ' on a des cercles
            epaisseur = swModel.ClosestDistance(Aretes(0), Aretes(1), p, p2)
            If epaisseur <= 0 Then MsgBox("Il y a un problème ici...")
        Else
            MsgBox("On ne doit jamias passer ici en principe")
        End If

        sFaceCoque.SetAttributDeCoque(epaisseur)

    End Sub

    ''' <summary>
    '''  sub qui découpe une face en fonction des paramètres qu'on lui passe.
    ''' </summary>
    ''' <remarks></remarks>
    Private Sub DécouperCoqueVolume()
        Static vc As Integer
        ' Algo:
        ' Pour chaque coque
        '   Si on doit couper la coque, alors on la coupe
        ' On découpe la face du volume
        ' ON identifie les faces internes


        Dim coque1 As SlyFaceCoque
        Dim interFF As InterCoqueVolume

        ' 1 - Pour toutes les coques
        Dim lst_face As New Collections.Generic.List(Of SlyFaceVolume)

        For Each coque1 In Commun.lst_FaceCoque

            For Each interFF In coque1.lst_InterCoqueVolume

                'If interFF.Face_A_Plat Then
                '    If Not lst_face.Contains(interFF.sFaceVolume) Then
                '        lst_face.Add(interFF.sFaceVolume)
                '        interFF.DécouperFace_A_Plat()
                '    End If


                If interFF.FaceDeSection Then
                    CoupeCoque1(interFF)
                    If Not Intersections.MultiDecoupageCoques Then interFF.MarquerFacesInternes()
                Else
                    interFF.DécouperVolume()
                    ' s'il y a plus d'une arête sur la coque qui touche au volume alors on doit découper la coque aussi
                    ' on sélectionne les 2 entités, puis on fait le sketch d'intersection et on évalue le nombre de courbes.
                    Dim swBod1 As sldworks.Body2 = interFF.sFaceVolume.SwFace.GetBody() : swBod1.Select2(False, Nothing)
                    Dim swBod2 As sldworks.Body2 = interFF.sFaceCoque.SwFace.GetBody() : swBod2.Select2(True, Nothing)

                    swModel.Sketch3DIntersections()
                    Dim sketch As sldworks.Sketch = swModel.GetActiveSketch2
                    swModel.InsertSketch2(True)
                    Dim swfeatS As sldworks.Feature = swModel.FeatureByPositionReverse(0)
                    swfeatS.Name = "Vérif_decoupage_coque#" & vc : vc += 1

                    If UBound(sketch.GetSketchSegments) > 1 Then
                        interFF.DecouperCoque()
                        interFF.QuelleAreteCoqueToucheVolume()
                    End If

                    If Not Intersections.MultiDecoupageCoques Then interFF.MarquerFacesInternes() ' étape 6

                End If
            Next interFF
        Next coque1


    End Sub


    ''' <summary>
    ''' Lorsque la coque est modélisée en partie par un morceau de volume.  Dessine et coupe.
    ''' </summary>
    ''' <param name="int">L'interFaceFace</param>
    ''' <remarks></remarks>
    Private Sub CoupeCoque1(ByRef int As InterCoqueVolume)
        Dim sVol As SlyFaceVolume
        Dim swFace As SldWorks.Face2
        Dim swent As SldWorks.Entity
        Dim sketch As SldWorks.Sketch
        Dim vEdges As Object
        Dim swArete As SldWorks.Edge
        Dim aretes() As SldWorks.Edge = Nothing
        Dim i As Integer

        ' a - prendre la face du volume  
        sVol = int.sFaceVolume
        swFace = sVol.SwFace
        sVol.Selectionner(0, False)

        '' b - lui mettre une esquisse
        'swModel.InsertSketch2(True)

        '' c - convertir l'esquisse déjà créée
        'sketch = int.sketch
        'swent = sketch
        'swent.Select2(False, 0)

        'swModel.SketchUseEdge2(False)

        '' si la face est carrée

        'If swFace.GetEdgeCount = 4 Then
        '    ' d - ajouter des lignes pour  compléter l'esquisse
        '    Dim longueur(3) As Double


        '    Dim vLine As Object
        '    Dim line1 As SldWorks.SketchLine, line2 As SldWorks.SketchLine
        '    Dim P1 As SldWorks.SketchPoint, P2 As SldWorks.SketchPoint, P3 As SldWorks.SketchPoint, P4 As SldWorks.SketchPoint
        '    Dim skSeg As SldWorks.SketchSegment
        '    vEdges = swFace.GetEdges()
        '    For Each swArete In vEdges
        '        ReDim Preserve aretes(i)
        '        aretes(i) = swArete
        '        i += 1
        '    Next
        '    For i = 0 To 3
        '        longueur(i) = Commun.GetLongueurArete(vEdges(i))
        '    Next i
        '    longueur = Ordonner(longueur, aretes)

        '    swent = aretes(3) : swent.Select2(False, 0) ' on prend la plus grande arète
        '    swModel.SketchUseEdge2(False)
        '    sketch = swModel.GetActiveSketch()
        '    vLine = sketch.GetSketchSegments()
        '    skSeg = vLine(0) : line1 = skSeg
        '    skSeg = vLine(1) : line2 = skSeg

        '    P1 = line1.GetStartPoint2()
        '    P2 = line1.GetEndPoint2()
        '    P3 = line2.GetStartPoint2()
        '    P4 = line2.GetEndPoint2()

        '    If Distance(P1.X, P1.Y, P1.Z, P3.X, P3.Y, P3.Z) < Distance(P1.X, P1.Y, P1.Z, P4.X, P4.Y, P4.Z) Then ' ligne entre 1 et 3
        '        swModel.CreateLine2(P1.X, P1.Y, P1.Z, P3.X, P3.Y, P3.Z)
        '        swModel.CreateLine2(P2.X, P2.Y, P2.Z, P4.X, P4.Y, P4.Z)
        '    Else ' ligne entre 1 et 4 & 2 et 3
        '        swModel.CreateLine2(P1.X, P1.Y, P1.Z, P4.X, P4.Y, P4.Z)
        '        swModel.CreateLine2(P2.X, P2.Y, P2.Z, P3.X, P3.Y, P3.Z)
        '    End If
        'Else
        '    sketch = swModel.GetActiveSketch()
        'End If
        'swModel.InsertSketch2(True)


        '' e - splitter
        'swent = swFace : swent.Select2(False, 1)
        'swent = sketch : swent.Select2(True, 4)
        'swModel.InsertSplitLineProject(False, False)


        sVol.Selectionner(32, False) 'swPart.Extension.SelectByID2("", "FACE", 0.05746341258515, 0.007456177698316, 0.04437034503314, False, 16, Nothing, 0)
        int.sFaceCoque.Selectionner(16, True) 'swPart.Extension.SelectByID2("", "FACE", 0.03921396269902, -0.007016448377556, 0, True, 32, Nothing, 0)
        swPart.FeatureManager.InsertSplitLineIntersect(7)

        ' f - mettre un FaceInterne sur les 2 faces résultantes.
        Dim swFeat As SldWorks.Feature
        Dim vFace As Object
        Dim swFace1 As SldWorks.Face2, swFace2 As SldWorks.Face2
        Dim attr As SldWorks.Attribute
        Dim nom As String
        Dim no As Long
        swFeat = swModel.FeatureByPositionReverse(0)
        vFace = swFeat.GetFaces

        swFace1 = vFace(0)
        swFace2 = vFace(1)

        sVol.AjouterFace(swFace1)
        sVol.AjouterFace(swFace2)

        'Dim eFace As New SlyFaceVolume(swFace1, True)
        'eFace.MettreAttributFaceInterne(swFace1, , False)

        'swent = swFace1
        'attr = swent.FindAttribute(Intersections.DefAttrFaceInterne, 0) ' si l'attribut existe déjà on pointe dessus.
        'nom = "FaceInterneCoque1"
        'If attr Is Nothing Then attr = Intersections.DefAttrFaceInterne.CreateInstance5(swModel, swFace1, nom, 0, 2) ' 0 = swThisconfig

        'While attr Is Nothing
        '    no += 1
        '    nom = "FaceInterneCoque" & CStr(no)
        '    attr = Intersections.DefAttrFaceInterne.CreateInstance5(swModel, swFace1, nom, 0, 2)
        'End While
        'GererDossiers("FaceInternes", nom)

        'swent = swFace2
        'attr = swent.FindAttribute(Intersections.DefAttrFaceInterne, 0) ' si l'attribut existe déjà on pointe dessus.
        'nom = "FaceInterneCoque1"
        'If attr Is Nothing Then attr = Intersections.DefAttrFaceInterne.CreateInstance5(swModel, swFace2, nom, 0, 2) ' 0 = swThisconfig

        'While attr Is Nothing
        '    no += 1
        '    nom = "FaceInterneCoque" & CStr(no)
        '    attr = Intersections.DefAttrFaceInterne.CreateInstance5(swModel, swFace2, nom, 0, 2)
        'End While
        'GererDossiers("FaceInternes", nom)


    End Sub


    ''' <summary>
    ''' Détecte s'il y a une intersection entre une coque et une face d'un volume
    ''' </summary>
    ''' <param name="sFaceCoque"></param>
    ''' <param name="sFaceVolume"></param>
    ''' <param name="dessiner">Si oui, alors on dessine une ligne à l'intersection</param>
    ''' <param name="sketch">Si dessiner est vrai, alors ce sketch contient la ligne d'intersection</param>
    ''' <returns>Vrai s'il y a une intersection</returns>
    ''' <remarks></remarks>
    Private Function DetectFaceFace(ByRef sFaceCoque As SlyFaceCoque, ByRef sFaceVolume As SlyFaceVolume, Optional ByRef dessiner As Boolean = False, Optional ByRef sketch As SldWorks.Sketch = Nothing) As Boolean
        For Each swfc As SldWorks.Face2 In sFaceCoque.lst_Faces
            For Each swFV As SldWorks.Face2 In sFaceVolume.lst_Faces
                If DetectFaceFace(swfc, swFV, dessiner, sketch) Then Return True
            Next
        Next
    End Function


    ''' <summary>
    ''' Function qui retourne vrai ou faux si 2 FACES (pas surfaces) se touchent.  La routine finit par appeler la sub dessinecourbe pour dessiner la courbe.
    ''' </summary>
    ''' <param name="face1">Première face</param>
    ''' <param name="face2">Seconde face</param>
    ''' <param name="dessiner">Si l'on veut dessiner une esquisse contenant la courbe</param>
    ''' <param name="Sketch">Si la précédente option est vrai, ce paramètre redonne l'esquisse qui contient la courbe (peut contenir plusieurs segments) d'intersection</param>
    ''' <returns>Vrai si les faces se touchent</returns>
    ''' <remarks>Si les 2 faces se touchent en un seul point alors ça retourne faux.</remarks>
    Private Function DetectFaceFace(ByRef face1 As SldWorks.Face2, ByRef face2 As SldWorks.Face2, Optional ByRef dessiner As Boolean = False, Optional ByRef Sketch As SldWorks.Sketch = Nothing) As Boolean
        'sub qui détecte si 2 faces se touchent et retourne vrai si c'est le cas.
        Dim Surface1 As sldworks.Surface
        Dim Surface2 As sldworks.Surface
        Dim curveArray As Object = Nothing
        Dim curve As sldworks.Curve
        Dim ret As Boolean
        Dim boundsArray As Object = Nothing
        Dim bounds() As Double
        Dim Point1(2) As Double
        Dim Point2(2) As Double
        Dim P1 As Object = Nothing, P2 As Object = Nothing
        Dim ClosestDist As Double

        Surface1 = face1.GetSurface
        Surface2 = face2.GetSurface

        ClosestDist = swModel.ClosestDistance(face1, face2, P1, P2)
        If ClosestDist > (Epsilon * 10) Then Return False

        ' curvearray est un tableau de curves, boundsarray est un tableau de T des limites de la courbe
        ret = Surface1.IntersectSurface(Surface2, curveArray, boundsArray) ' si c'est une ligne, retourne une ligne infinie...
        ' il faut alors renvoyer vers le detectFaceArete()
        If Not ret Then Return False
        bounds = boundsArray

        Try
            curve = curveArray(0)
        Catch
            Return False  'GoTo Fairesketch 'ne va pas toujours marcher, mais je n'ai rien de mieu pour l'instant...
            'MsgBox("On a une intersection où 2 faces sont sur la même surface...") ' en théorie...
            ' ouch... pas certain.
        End Try


        If Not curve.IsLine Then ' si c'est une ligne, alors le principe de longueur ne fonctionnera pas...
            Dim longueur As Double = curve.GetLength2(bounds(0), bounds(1))
            If longueur < Epsilon Then Return False ' on a juste un point d'intersection
        Else
            Dim vParam As Object
            vParam = curve.GetClosestPointOn(P1(0), P1(1), P1(2))  ' vparam(3) est le U

            P1 = curve.Evaluate(vParam(3) + 100000 * Epsilon) : Point1(0) = P1(0) : Point1(1) = P1(1) : Point1(2) = P1(2)
            P2 = curve.Evaluate(vParam(3) - 100000 * Epsilon) : Point2(0) = P2(0) : Point2(1) = P2(1) : Point2(2) = P2(2)

            ' si p1 et p2 sont identiques alors on a un point D'intersection, ce que l'on ne veut pas
            If Distance(Point1, Point2) < 1000 * Epsilon Then Return False

            If ((Distance(face1, Point1(0), Point1(1), Point1(2)) < (1000 * Epsilon)) AndAlso (Distance(face2, Point1(0), Point1(1), Point1(2)) < (1000 * Epsilon))) OrElse ((Distance(face1, Point2(0), Point2(1), Point2(2)) < (Epsilon * 1000)) AndAlso (Distance(face2, Point2(0), Point2(1), Point2(2)) < (1000 * Epsilon))) Then
                ' return true
            Else

                Dim swent2 As sldworks.Entity
                Dim feat2 As sldworks.Feature
                swModel.Insert3DSketch2(False)
                swModel.ClearSelection2(True)
                'swent2 = face1 : swent2.Select2(False, 0)
                'swent2 = face2 : swent2.Select2(True, 0)
                Dim swBod1 As sldworks.Body2 = face1.GetBody : swBod1.Select2(False, Nothing)
                Dim swBod2 As sldworks.Body2 = face2.GetBody : swBod2.Select2(True, Nothing)

                swModel.Sketch3DIntersections()

                swModel.Insert3DSketch2(False)
                swModel.EditRebuild3()
                feat2 = swModel.FeatureByPositionReverse(0)

                Debug.Print(feat2.Name)

                Sketch = feat2.GetSpecificFeature2
                feat2.Name = "TouchePas" & CStr(Rnd())
                'MsgBox(face1.GetArea)
                Dim vSeg2 As Object = Sketch.GetSketchSegments()
                If vSeg2 Is Nothing Then
                    swent2 = feat2
                    swent2.Select(False)
                    swModel.EditDelete()
                    Return False
                Else
                    If Math.Abs(bounds(0)) > 490 Then Return False ' se touchent à l'infini...
                    Return True
                End If

            End If
        End If

Fairesketch:
        Dim swent As sldworks.Entity
        Dim feat As sldworks.Feature
        swModel.Insert3DSketch2(False)
        swModel.ClearSelection2(True)
        'swent = face1 : swent.Select2(False, 0)
        'swent = face2 : swent.Select2(True, 0)
        Dim swBod1a As sldworks.Body2 = face1.GetBody : swBod1a.Select2(False, Nothing)
        Dim swBod2a As sldworks.Body2 = face2.GetBody : swBod2a.Select2(True, Nothing)

        swModel.Sketch3DIntersections()

        swModel.Insert3DSketch2(False)
        swModel.EditRebuild3()
        feat = swModel.FeatureByPositionReverse(0)
        Debug.Print(feat.Name)
        Sketch = feat.GetSpecificFeature2
        Dim vSeg As Object = Sketch.GetSketchSegments()
        If vSeg Is Nothing Then
            swent = feat
            swent.Select(False)
            swModel.EditDelete()
            Return False
        Else
            Return True
        End If

        Return True

    End Function


    Public Function trier(ByRef lst() As Double) As Double()
        ' sub qui retourne les T min et T max en fonction de la collection.
        ' normalement on enlève les 2 points max et min et il devrait nous rester un nombre pair de valeurs.
        Dim i As Integer, j As Integer
        Dim lst2() As Double

        Dim temp As Double

        If UBound(lst) < 1 Then Return Nothing

        For j = 0 To UBound(lst) - 1
            For i = 0 To UBound(lst) - j - 1
                If lst(i) > lst(i + 1) Then
                    temp = lst(i)
                    lst(i) = lst(i + 1)
                    lst(i + 1) = temp
                End If
            Next i
        Next j

        ReDim lst2(UBound(lst) - 2)
        For i = 1 To UBound(lst) - 1
            lst2(i - 1) = lst(i)
        Next
        Return lst2


    End Function

    ''' <summary>
    ''' Function qui remet en ordre croissant les valeurs d'un tableau
    ''' </summary>
    ''' <param name="lst">Le tableau de valeurs à ordonner</param>
    ''' <param name="liste2">un second tableau que l'on peut ordonner</param>
    ''' <returns></returns>
    ''' <remarks></remarks>
    Public Function Ordonner(ByRef lst() As Double, Optional ByRef liste2() As SldWorks.Edge = Nothing) As Double()

        Dim i As Integer, j As Integer
        Dim temp2 As Object

        Dim temp As Double

        If UBound(lst) < 1 Then Return Nothing

        For j = 0 To UBound(lst) - 1
            For i = 0 To UBound(lst) - j - 1
                If lst(i) > lst(i + 1) Then
                    temp = lst(i)
                    lst(i) = lst(i + 1)
                    lst(i + 1) = temp
                    If Not liste2 Is Nothing Then
                        temp2 = liste2(i)
                        liste2(i) = liste2(i + 1)
                        liste2(i + 1) = temp2
                    End If

                End If
            Next i
        Next j

        Return lst


    End Function


    ''' <summary>
    ''' Sub qui fait comme le tavail de Antoine et Rémi et qui créé un fichier pour la carte de taille.  Si l'option de pré-optimisation est enclenchée, la sub ne fait qu'ajouter des éléments au fichier .txt
    ''' </summary>
    ''' <remarks></remarks>
    Public Sub FairePreCarte()

        Dim ligne_txt As String
        Dim fichier As System.IO.StreamWriter
        Dim EcartNodal As Double = Commun.ÉcartNodal

        ' si le fichier pog n'est pas encore créé, on doit le faire...
        If Commun.NomFichierPog = Nothing Then
            ' la première ligne donne la boite englobante, la seconde la taille ENG

            Dim path As String = Nothing
            Dim CMDialogl As New Windows.Forms.SaveFileDialog
            CMDialogl.DefaultExt = ".txt"
            CMDialogl.Filter = "Fichiers PoG (*.txt)|*.txt|Tout fichiers(*.*)|*.*"

            CMDialogl.OverwritePrompt = True
            CMDialogl.Title = "Sélectionnez le fichier pour enregistrer les points"
            CMDialogl.ShowDialog()
            path = CMDialogl.FileName
            path = Txtpath(path) : Commun.NomFichierPog = path
            If path Is Nothing Or path = "" Then MsgBox("Aucun fichier sélectionné, sortie du programme!", MsgBoxStyle.Critical, "Erreur!")

            fichier = System.IO.File.CreateText(path)

            Dim vBox As Object = swPart.GetPartBox(True)
            Dim box() As Double = vBox

            Dim centre(2) As Double ' le centre de la boite englobante
            Dim longueurs(2) As Double

            centre(0) = (box(3) + box(0)) / 2
            centre(1) = (box(4) + box(1)) / 2
            centre(2) = (box(5) + box(2)) / 2

            longueurs(0) = (box(3) - box(0)) * 1.25
            longueurs(1) = (box(4) - box(1)) * 1.25
            longueurs(2) = (box(5) - box(2)) * 1.25


            ligne_txt = CStr(centre(0) - longueurs(0) / 2) & " " & CStr(centre(1) - longueurs(1) / 2) & " " & CStr(centre(2) - longueurs(2) / 2) & " " & CStr(centre(0) + longueurs(0) / 2) & " " & CStr(centre(1) + longueurs(1) / 2) & " " & CStr(centre(2) + longueurs(2) / 2)
            fichier.WriteLine(Replace(ligne_txt, ",", "."))

            fichier.WriteLine(Replace(CStr(EcartNodal), ",", "."))
        Else
            ' Le fichier pog existe déjà...
            fichier = System.IO.File.AppendText(Commun.NomFichierPog)
            'fichier.WriteLine(Replace(CStr(EcartNodal), ",", "."))
        End If


        ' *** C'est ici que le fun se passe.


        ' 1 - On parcourt toutes les faces, les listes sont déjà RE-crées, donc pas de multifaces
        For Each sFace As SlyFaceVolume In Commun.lst_FaceVolume
            sFace.MettrePointSurPOG(fichier)
        Next

        For Each sFace As SlyFaceCoque In Commun.lst_FaceCoque
            sFace.MettrePointSurPOG(fichier)
        Next

        ' *** Fin
        fichier.Close()
    End Sub

    ''' <summary>
    ''' Fonction qui compare 2 surfaces
    ''' </summary>
    ''' <param name="swSurf1">La première surface</param>
    ''' <param name="swSurf2">La seconde</param>
    ''' <returns>Vrai si les 2 surfaces sont identiques, faux sinon</returns>
    ''' <remarks></remarks>
    Public Function ComparerSurfaces(ByRef swSurf1 As sldworks.Surface, ByRef swSurf2 As sldworks.Surface) As Boolean
        'swSurf1.GetBSurfParams2(True,False,


        If swSurf1.IsPlane And swSurf2.IsPlane Then
            Dim obj1 As Object = swSurf1.PlaneParams
            Dim obj2 As Object = swSurf2.PlaneParams
            For i As Integer = 0 To 2
                If Not Math.Abs(obj1(i) - obj2(i)) < Epsilon Then Return False
            Next
            Return True
        ElseIf swSurf1.IsBlending And swSurf2.IsBlending Then
            ' pas de blendingParams

        ElseIf swSurf1.IsCone And swSurf2.IsCone Then
            Dim obj1 As Object = swSurf1.ConeParams
            Dim obj2 As Object = swSurf2.ConeParams
            For i As Integer = 0 To 7
                If Not Math.Abs(obj1(i) - obj2(i)) < Epsilon Then Return False
            Next
            Return True

        ElseIf swSurf1.IsCylinder And swSurf2.IsCylinder Then
            Dim obj1 As Object = swSurf1.CylinderParams
            Dim obj2 As Object = swSurf2.CylinderParams
            For i As Integer = 0 To 6
                If Not Math.Abs(obj1(i) - obj2(i)) < Epsilon Then Return False
            Next
            Return True
        ElseIf swSurf1.IsForeign And swSurf2.IsForeign Then
            ' ???
        ElseIf swSurf1.IsOffset And swSurf2.IsOffset Then
            ' Pas de offsetParam

        ElseIf swSurf1.IsParametric And swSurf2.IsParametric Then
            ' pas de ParametricParams

        ElseIf swSurf1.IsRevolved And swSurf2.IsRevolved Then
            ' pas de revolvedparams

        ElseIf swSurf1.IsSphere And swSurf2.IsSphere Then
            Dim obj1 As Object = swSurf1.SphereParams
            Dim obj2 As Object = swSurf2.SphereParams

            For i As Integer = 0 To 3
                If Not Math.Abs(obj1(i) - obj2(i)) < Epsilon Then Return False
            Next
            Return True
        ElseIf swSurf1.IsSwept And swSurf2.IsSwept Then
            ' merde, il n'y a pas de Sweptparams

        ElseIf swSurf1.IsTorus And swSurf2.IsTorus Then
            Dim obj1 As Object = swSurf1.TorusParams
            Dim obj2 As Object = swSurf2.TorusParams

            For i As Integer = 0 To 7
                If Not Math.Abs(obj1(i) - obj2(i)) < Epsilon Then Return False
            Next
            Return True
        Else
            ' faut quand même évaluer si on a la bonne chose...
            Try
                Dim pt As Object = swSurf1.Evaluate(0.5, 0.5, 0, 0)
                Dim obj As Object = swSurf1.GetClosestPointOn(pt(0), pt(1), pt(2)) ' x, y , z, U, V
                Dim obj1 As Object = swSurf1.EvaluateAtPoint(obj(0), obj(1), obj(2))

                Dim objb As Object = swSurf2.GetClosestPointOn(pt(0), pt(1), pt(2)) ' x, y , z, U, V
                Dim obj2 As Object = swSurf2.EvaluateAtPoint(objb(0), objb(1), objb(2))

                If Math.Abs(obj1(9) - obj2(9)) < Epsilon And Math.Abs(obj1(10) - obj2(10)) < Epsilon Then
                    ' c'est cheap.... mais
                    Return True
                End If
                Return False
            Catch
                Return False
            End Try

        End If

        Return False


    End Function


End Module
Revision:	205
Committed:	Thu Jul 23 20:53:57 2009 UTC (15 years, 9 months ago) by bournival
File size:	111956 byte(s)
Log Message:	Commit de MAGiC_SLD pendant que j'y pense. Les modifications ne devraient pas concerner personne d'autre que moi. -- Sylvain