SUMARI 1 I. CLASSES DEL SUBPROGRAMA DE CÀLCUL 3. A. Classe CEstructura...3 B. Classe CElement...74 II. MÒDULS DEL SUBPROGRAMA DE CÀLCUL 127

Transcription

1 1 Sumari SUMARI 1 I. CLASSES DEL SUBPROGRAMA DE CÀLCUL 3 A. Classe CEstructura...3 B. Classe CElement...74 II. MÒDULS DEL SUBPROGRAMA DE CÀLCUL 127 A. Mòdul principal...127

2

3 3 I. Classes del subprograma de càlcul A. Classe CEstructura Option Explicit Const Tiny = 1E-30 Private m_colbarra As CColBarra Private m_colnode As CColNode Private m_nom As String Private m_colaccio As CColAccio Private m_matriuderigidesaglobal As Variant Private m_matriuderigidesareduida As Variant Private m_matriuderigidesareduidacalculada As Boolean Private m_matriuderigidesareduidainvertida As Variant Private m_matriuderigidesareduidainvertidacalculada As Boolean Private m_vectordecarregaglobal As Variant 'EN PROCÉS DE SUBSTITUCIÓ Private m_matriudecarregaglobal As Variant Private m_vectordecarregareduit As Variant 'EN PROCÉS DE SUBSTITUCIÓ Private m_matriudecarregareduida As Variant Private m_vectordecorrimentsnodalsreduit As Variant 'EN PROCÉS DE SUBSTITUCIÓ Private m_matriudecorrimentsnodalsreduit As Variant Private m_vectordecorrimentsnodalsglobal As Variant 'EN PROCÉS DE SUBSTITUCIÓ Private m_matriudecorrimentsnodalsglobal As Variant Private m_vectordereaccionsnodalsglobal As Variant 'EN PROCÉS DE SUBSTITUCIÓ Private m_matriudereaccionsnodalsglobal As Variant Private m_matriudeconnectivitat As Variant Private m_vectorsdecorrimentselementals As Collection 'NO ES FA SERVIR Private m_vectorsdecorrimentselementalscalculats As Boolean 'NO ES FA SERVIR Private m_vectorsdecarregueselementals As Collection 'NO ES FA SERVIR 'VARIABLES DEL CÀLCUL NO LINEAL 'PRIMER IMPLEMENTO EL MÈTODE DE NEWTON-RAPHSON NO MODIFICAT Private m_matriuderigidesaglobalnewton As Variant 'OK Private m_matriuderigidesareduidanewton As Variant 'OK Private m_vectordecorrimentsnodalsreduitparcialnewton As Variant 'EN PROCÉS DE SUBSTITUCIÓ Private m_vectordecorrimentsnodalsglobalparcialnewton As Variant 'EN PROCÉS DE SUBSTITUCIÓ Private m_vectordecorrimentsnodalsreduitnewton As Variant 'EN PROCÉS DE SUBSTITUCIÓ

4 4 Private m_vectordecorrimentsnodalsglobalnewton As Variant 'EN PROCÉS DE SUBSTITUCIÓ Private m_matriudevectorsdecarregaincrementals As Variant 'NOU PER AL CÀLCUL INCREMENTAL Private m_matriudevectorsdecarreganoequilibrats As Variant 'NOU PER AL CÀLCUL INCREMENTAL Private m_matriudevectorsdecorrimentsnodalsreduitnewton As Variant 'NOU PER AL CÀLCUL INCREMENTAL Private m_matriudevectorsdecorrimentsnodalsglobalnewton As Variant 'NOU PER AL CÀLCUL INCREMENTAL Private m_vectordecarregareduitnewton As Variant 'EN PROCÉS DE SUBSTITUCIÓ Private m_matriuderigidesareduidanewtoninvertida As Variant 'OK Private m_matriuderigidesareduidanewtonmodificatinvertida As Variant 'OK Private m_vectorsdecorrimentselementalsnewton As Collection ' Private m_vectorcarreganoequilibradanewton As Variant 'EN PROCÉS DE SUBSTITUCIÓ Private m_vectordecorrimentsnodalsglobalnewtonacumulat As Variant Private m_vectordecarregadereferencianewton As Variant Private m_tin As Double Private m_modalitatdecalcullineal As Integer '1 DESCOMPOSICIÓ LU '2 GAUSS-JORDAN '3 ELIMINACIÓ GAUSSIANA '4 DESCOMPOSICIÓ CHOLESKY Private m_accioactivada As Long Private m_tramsdelcalculincremental As Long Private m_nombremaximditeracions As Long Private m_modalitatdecalculnolineal As Integer '1 MÈTODE INCREMENTAL '2 MÈTODE NEWTON-RAPHSON '3 MÈTODE NEWTON-RAPHSON MODIFICAT '4 MÈTODE NEWTON-RAPHSON INCREMENTAL '5 MÈTODE NEWTON-RAPHSON INCREMENTAL MODIFICAT Sub CalculaMatriuDeVectorsDeCarregaNoEquilibradaNR(Iteracio As Long, AccioActivada As Long) Dim MReal() As Double Dim MNoEq() As Double Dim CIni() As Double Dim i As Long, imax As Long Dim MRealAcum() As Double Dim SumaParcial As Double Dim j As Long

5 5 imax = UBound(Estructura.MatriuDeCarregaReduida, 1) ReDim CIni(1 To imax, 1 To 1) ReDim MRealAcum(i To imax, 1 To 1) For i = 1 To imax CIni(i, 1) = Estructura.MatriuDeCarregaReduida(i, AccioActivada) Next i AssignaVariantAMatriu m_matriudevectorsdecarregaincrementals, MReal For i = 1 To imax SumaParcial = 0 For j = 1 To Iteracio SumaParcial = SumaParcial + MReal(i, j) Next j MRealAcum(i, 1) = SumaParcial Next i AssignaVariantAMatriu m_matriudevectorsdecarreganoequilibrats, MNoEq For i = 1 To imax MNoEq(i, Iteracio) = CIni(i, 1) - MRealAcum(i, 1) Next i m_matriudevectorsdecarreganoequilibrats = MNoEq Sub CalculaMatriuDeVectorsDeCarregaRealNR(Iteracio As Long) Dim MCor() As Double Dim MCorAcum() As Double Dim KR() As Double Dim F() As Double Dim CReal() As Double Dim i As Long, imax As Long Dim j As Long Dim SumaParcial As Double AssignaVariantAMatriu m_matriudevectorsdecorrimentsnodalsreduitnewton, MCor AssignaVariantAMatriu Estructura.MatriuDeRigidesaReduidaNewton, KR imax = UBound(MCor, 1) ReDim MCorAcum(1 To imax, 1 To 1) 'For i = 1 To imax

6 6 ' SumaParcial = 0 ' For j = 1 To Iteracio ' SumaParcial = SumaParcial + MCor(i, j) ' Next j ' MCorAcum(i, 1) = SumaParcial 'Next i For i = 1 To imax SumaParcial = 0 If Iteracio > 1 Then For j = Iteracio - 1 To Iteracio - 1 SumaParcial = SumaParcial + MCor(i, j) Next j Else SumaParcial = 0 MCorAcum(i, 1) = SumaParcial Next i MultiplicaMatrius KR, MCorAcum, F AssignaVariantAMatriu m_matriudevectorsdecarregaincrementals, CReal For i = 1 To imax CReal(i, Iteracio) = F(i, 1) Next i m_matriudevectorsdecarregaincrementals = CReal Sub InicialitzaMatriuDeVectorsDeCarregaNoEquilibradaReduitsNR(AccioActivada As Long, MaxIteracio As Long) Dim A() As Long Dim Max As Long Dim i As Long Max = UBound(Estructura.MatriuDeCarregaReduida, 1) ReDim A(1 To Max, 1 To MaxIteracio) 'For i = 1 To Max 'A(i, 1) = Estructura.MatriuDeCarregaReduida(i, AccioActivada) 'Next i m_matriudevectorsdecarreganoequilibrats = A

7 7 Property Get MatriuDeVectorsDeCorrimentsNodalsGlobalNewton() As Long MatriuDeVectorsDeCorrimentsNodalsGlobalNewton = m_matriudevectorsdecorrimentsnodalsglobalnewton Property Let MatriuDeVectorsDeCorrimentsNodalsGlobalNewton(ByVal Val As Long) m_matriudevectorsdecorrimentsnodalsglobalnewton = Val Property Get ModalitatDeCalculNOLineal() As Long ModalitatDeCalculNOLineal = m_modalitatdecalculnolineal Property Let ModalitatDeCalculNOLineal(ByVal Val As Long) m_modalitatdecalculnolineal = Val Property Get AccioActivada() As Long AccioActivada = m_accioactivada Property Let AccioActivada(ByVal Val As Long) m_accioactivada = Val Property Get TramsDelCalculIncremental() As Long TramsDelCalculIncremental = m_tramsdelcalculincremental Property Let TramsDelCalculIncremental(ByVal Val As Long) m_tramsdelcalculincremental = Val Property Get NombreMaximDIteracions() As Long NombreMaximDIteracions = m_nombremaximditeracions Property Let NombreMaximDIteracions(ByVal Val As Long) m_nombremaximditeracions = Val Property Get MatriuDeVectorsDeCarregaIncrementals() As Variant MatriuDeVectorsDeCarregaIncrementals = m_matriudevectorsdecarregaincrementals Property Let MatriuDeVectorsDeCarregaIncrementals(ByVal Val As Variant) m_matriudevectorsdecarregaincrementals = Val Property Get Tin() As Double 'AQUEST VALOR ÉS GAMMA PIVOT Tin = m_tin Property Let Tin(ByVal Val As Double) m_tin = Val

8 8 Property Get ModalitatDeCalculLineal() As Integer ModalitatDeCalculLineal = m_modalitatdecalcullineal Property Let ModalitatDeCalculLineal(ByVal Val As Integer) m_modalitatdecalcullineal = Val Sub CalculaMatriuDeVectorsDeCarregaIncrementalsReduits(NumeroDeIteracions As Long, _ AccioActivada As Long) Dim MVC() As Double 'MATRIUS DE VECTORS DE CÀRREGA PER A CADA ITERACIÓ Dim IteracioEnCurs As Long Dim i As Long, maxi As Long, j As Long, maxj As Long Dim Multiplicador As Double Estructura.CalculaMatriuDeCarregaGlobal Estructura.CalculaMatriuDeCarregaReduida maxi = UBound(Estructura.MatriuDeCarregaReduida, 1) maxj = NumeroDeIteracions ReDim MVC(1 To maxi, 1 To maxj) For j = 1 To maxj Multiplicador = 1 / maxj For i = 1 To maxi MVC(i, j) = Multiplicador * Estructura.MatriuDeCarregaReduida(i, AccioActivada) Next i Next j Estructura.MatriuDeVectorsDeCarregaIncrementals = MVC Sub InicialitzaMatriuDeVectorsDeCarregaNoRepartida(NumeroDeIteracions As Long, _ AccioActivada As Long) Dim MVC() As Double 'MATRIUS DE VECTORS DE CÀRREGA PER A CADA ITERACIÓ Dim IteracioEnCurs As Long Dim i As Long, maxi As Long, j As Long, maxj As Long Dim Multiplicador As Double

9 9 Estructura.CalculaMatriuDeCarregaGlobal Estructura.CalculaMatriuDeCarregaReduida maxi = UBound(Estructura.MatriuDeCarregaReduida, 1) maxj = NumeroDeIteracions ReDim MVC(1 To maxi, 1 To maxj) For j = 1 To maxj Multiplicador = 1 / maxj For i = 1 To maxi MVC(i, j) = Multiplicador * Estructura.MatriuDeCarregaReduida(i, AccioActivada) Next i Next j Estructura.MatriuDeVectorsDeCarregaIncrementals = MVC Public Sub EsborraFilaDUnaMatriu(MatriuOrigen() As Double, MatriuDesti() As Double, Fila As Long) Dim i As Long, j As Long ReDim MatriuDesti(1 To UBound(MatriuOrigen, 1) - 1, 1 To UBound(MatriuOrigen, 2)) For i = 1 To UBound(MatriuOrigen, 1) For j = 1 To UBound(MatriuOrigen, 2) Select Case Fila Case Is > i MatriuDesti(i, j) = MatriuOrigen(i, j) Case Is = i 'NO FAI RES Case Is < i MatriuDesti(i - 1, j) = MatriuOrigen(i, j) End Select Next j Next i Public Sub EsborraColumnaDUnaMatriu(MatriuOrigen() As Double, MatriuDesti() As Double, Columna As Long)

10 10 Dim i As Long, j As Long ReDim MatriuDesti(1 To UBound(MatriuOrigen, 1), 1 To UBound(MatriuOrigen, 2) - 1) For i = 1 To UBound(MatriuOrigen, 1) For j = 1 To UBound(MatriuOrigen, 2) Select Case Columna Case Is > j MatriuDesti(i, j) = MatriuOrigen(i, j) Case Is = j 'NO FAI RES Case Is < j MatriuDesti(i, j - 1) = MatriuOrigen(i, j) End Select Next j Next i Public Function ExisteixElNumeroEnElVector(Vector() As Double, Numero As Double) As Boolean Dim i As Long, imax As Long imax = UBound(Vector, 1) For i = 1 To imax If Vector(i) = Numero Then ExisteixElNumeroEnElVector = True Exit Function ExisteixElNumeroEnElVector = False Next i End Function Sub InversaMatriuLU(MatriuAInvertir() As Double, MatriuInvertida() As Double) Dim d As Double Dim col() As Double Dim i As Long, j As Long, indx() As Long Estructura.LUDescomposicio MatriuAInvertir, indx, d ReDim col(1 To UBound(MatriuAInvertir, 1)) ReDim MatriuInvertida(1 To UBound(MatriuAInvertir, 1), 1 To UBound(MatriuAInvertir, 1)) For j = 1 To UBound(MatriuAInvertir, 1) For i = 1 To UBound(MatriuAInvertir, 1)

11 11 col(i) = 0 Next i col(j) = 1 Estructura.LUbksb MatriuAInvertir, indx, col For i = 1 To UBound(MatriuAInvertir, 1) MatriuInvertida(i, j) = col(i) Next i Principal.Escrollant UBound(MatriuAInvertir, 1), j, "Inversió LU" 'Principal.ProgresText UBound(MatriuAInvertir, 1), j, "LU Etapa ", True Next j Public Property Let MatriuDeReaccionsNodalsGlobal(ByVal Val As Variant) m_matriudereaccionsnodalsglobal = Val Public Property Get MatriuDeReaccionsNodalsGlobal() As Variant MatriuDeReaccionsNodalsGlobal = m_matriudereaccionsnodalsglobal Public Property Let MatriuDeCorrimentsNodalsGlobal(ByVal Val As Variant) m_matriudecorrimentsnodalsglobal = Val Public Property Get MatriuDeCorrimentsNodalsGlobal() As Variant MatriuDeCorrimentsNodalsGlobal = m_matriudecorrimentsnodalsglobal Public Property Let ColAccio(ByVal Val As CColAccio) Set m_colaccio = Val Public Property Get ColAccio() As CColAccio Set ColAccio = m_colaccio Public Property Let VectorDeCarregaDeReferenciaNewton(ByVal Val As Variant) m_vectordecarregadereferencianewton = Val Public Property Get VectorDeCarregaDeReferenciaNewton() As Variant VectorDeCarregaDeReferenciaNewton = m_vectordecarregadereferencianewton Public Property Let VectorDeCorrimentsNodalsGlobalNewtonAcumulat(ByVal Val As Variant) m_vectordecorrimentsnodalsglobalnewtonacumulat = Val Public Property Get VectorDeCorrimentsNodalsGlobalNewtonAcumulat() As Variant VectorDeCorrimentsNodalsGlobalNewtonAcumulat = m_vectordecorrimentsnodalsglobalnewtonacumulat

12 12 Sub InicialitzaMatriuDeVectorsCorrimentsNodalsReduitNewtonIncrementalAcumulats() Dim i As Long, j As Long Dim B() As Double i = UBound(m_MatriuDeVectorsDeCarregaIncrementals, 1) j = UBound(m_MatriuDeVectorsDeCarregaIncrementals, 2) ReDim B(1 To i, 1 To j) m_matriudevectorsdecorrimentsnodalsreduitnewton = B Sub InicialitzaMatriuDeVectorsCorrimentsNodalsReduitNRAcumulats() Dim i As Long, j As Long Dim B() As Double i = UBound(m_MatriuDeVectorsDeCarregaIncrementals, 1) j = UBound(m_MatriuDeVectorsDeCarregaIncrementals, 2) ReDim B(1 To i, 1 To j) m_matriudevectorsdecorrimentsnodalsreduitnewton = B Sub InicialitzaMatriuDeVectorsCorrimentsNodalsGlobalNewtonIncrementalAcumulats() Dim i As Long, j As Long Dim B() As Double i = Estructura.ColNode.Count * 6 j = UBound(m_MatriuDeVectorsDeCarregaIncrementals, 2) ReDim B(1 To i, 1 To j) m_matriudevectorsdecorrimentsnodalsglobalnewton = B Sub InicialitzaMatriuDeVectorsCorrimentsNodalsGlobalNRAcumulats() Dim i As Long, j As Long Dim B() As Double i = Estructura.ColNode.Count * 6 j = UBound(m_MatriuDeVectorsDeCarregaIncrementals, 2) ReDim B(1 To i, 1 To j) m_matriudevectorsdecorrimentsnodalsglobalnewton = B Sub MultiplicaVectorPerNumero(ByVal VectorADividir As Variant, ByVal Numero As Double, ByRef Resultat As Variant) Dim i As Long Dim imax As Long

13 13 imax = UBound(VectorADividir, 1) ReDim Resultat(1 To imax, 1 To 1) For i = 1 To imax Resultat(i, 1) = VectorADividir(i, 1) * Numero Next i Public Property Let VectorCarregaNoEquilibradaNewton(ByVal Val As Variant) m_vectorcarreganoequilibradanewton = Val Public Property Get VectorCarregaNoEquilibradaNewton() As Variant VectorCarregaNoEquilibradaNewton = m_vectorcarreganoequilibradanewton Public Property Let VectorDeCorrimentsNodalsReduitParcialNewton(ByVal Val As Variant) m_vectordecorrimentsnodalsreduitparcialnewton = Val Public Property Get VectorDeCorrimentsNodalsReduitParcialNewton() As Variant VectorDeCorrimentsNodalsReduitParcialNewton = m_vectordecorrimentsnodalsreduitparcialnewton Sub AssignaResultatsNewtonAVectorOrdinari() Dim NodeEnCurs As Long Dim MaxNode As Long Dim i As Integer 'm_vectordecorrimentsnodalsglobal = _ m_vectordecorrimentsnodalsglobalnewton MaxNode = m_colnode.count For NodeEnCurs = 1 To MaxNode m_colnode.item(nodeencurs).vectorcorriments(1) = m_vectordecorrimentsnodalsglobalnewton((nodeencurs - 1) * 6 + 1, 1) m_colnode.item(nodeencurs).vectorcorriments(2) = m_vectordecorrimentsnodalsglobalnewton((nodeencurs - 1) * 6 + 2, 1) m_colnode.item(nodeencurs).vectorcorriments(3) = m_vectordecorrimentsnodalsglobalnewton((nodeencurs - 1) * 6 + 3, 1) m_colnode.item(nodeencurs).vectorcorriments(4) = m_vectordecorrimentsnodalsglobalnewton((nodeencurs - 1) * 6 + 4, 1) m_colnode.item(nodeencurs).vectorcorriments(5) = m_vectordecorrimentsnodalsglobalnewton((nodeencurs - 1) * 6 + 5, 1) m_colnode.item(nodeencurs).vectorcorriments(6) = m_vectordecorrimentsnodalsglobalnewton((nodeencurs - 1) * 6 + 6, 1) 'For i = 1 To 6 ' If Abs(m_ColNode.Item(NodeEnCurs).VectorCorriments(i)) > # Then ' MsgBox "L'algorisme de Newton-Raphson no ha convergit", vbcritical, "Error" ' Exit Sub ' ' 'Next i Next NodeEnCurs

14 14 'MsgBox " Algorisme de Newton-Raphson realitzat satisfactòriament", vbinformation Sub AssignaResultatsNewtonIncrementalAVectorOrdinari() Dim NodeEnCurs As Long Dim MaxNode As Long Dim i As Integer 'm_vectordecorrimentsnodalsglobal = _ m_vectordecorrimentsnodalsglobalnewton MaxNode = m_colnode.count For NodeEnCurs = 1 To MaxNode m_colnode.item(nodeencurs).vectorcorriments(1) = m_vectordecorrimentsnodalsglobalnewtonacumulat((nodeencurs - 1) * 6 + 1, 1) m_colnode.item(nodeencurs).vectorcorriments(2) = m_vectordecorrimentsnodalsglobalnewtonacumulat((nodeencurs - 1) * 6 + 2, 1) m_colnode.item(nodeencurs).vectorcorriments(3) = m_vectordecorrimentsnodalsglobalnewtonacumulat((nodeencurs - 1) * 6 + 3, 1) m_colnode.item(nodeencurs).vectorcorriments(4) = m_vectordecorrimentsnodalsglobalnewtonacumulat((nodeencurs - 1) * 6 + 4, 1) m_colnode.item(nodeencurs).vectorcorriments(5) = m_vectordecorrimentsnodalsglobalnewtonacumulat((nodeencurs - 1) * 6 + 5, 1) m_colnode.item(nodeencurs).vectorcorriments(6) = m_vectordecorrimentsnodalsglobalnewtonacumulat((nodeencurs - 1) * 6 + 6, 1) For i = 1 To 6 If Abs(m_ColNode.Item(NodeEnCurs).VectorCorriments(i)) > # Then MsgBox "L'algorisme de Newton-Raphson Incremental no ha convergit", vbcritical, "Error" Exit Sub Next i Next NodeEnCurs MsgBox " L' Algorisme Newton-Raphson Incremental ha Convergit Correctament", vbinformation Sub AssignaResultatsIncrementalsAVectorOrdinari() Dim NodeEnCurs As Long Dim MaxNode As Long Dim i As Integer 'm_vectordecorrimentsnodalsglobal = _ m_vectordecorrimentsnodalsglobalnewton MaxNode = m_colnode.count For NodeEnCurs = 1 To MaxNode m_colnode.item(nodeencurs).vectorcorriments(1) = m_vectordecorrimentsnodalsglobalnewton((nodeencurs - 1) * 6 + 1, 1) m_colnode.item(nodeencurs).vectorcorriments(2) = m_vectordecorrimentsnodalsglobalnewton((nodeencurs - 1) * 6 + 2, 1) m_colnode.item(nodeencurs).vectorcorriments(3) =

15 15 m_vectordecorrimentsnodalsglobalnewton((nodeencurs - 1) * 6 + 3, 1) m_colnode.item(nodeencurs).vectorcorriments(4) = m_vectordecorrimentsnodalsglobalnewton((nodeencurs - 1) * 6 + 4, 1) m_colnode.item(nodeencurs).vectorcorriments(5) = m_vectordecorrimentsnodalsglobalnewton((nodeencurs - 1) * 6 + 5, 1) m_colnode.item(nodeencurs).vectorcorriments(6) = m_vectordecorrimentsnodalsglobalnewton((nodeencurs - 1) * 6 + 6, 1) For i = 1 To 6 If Abs(m_ColNode.Item(NodeEnCurs).VectorCorriments(i)) > # Then MsgBox "L'algorisme de Newton-Raphson no ha convergit", vbcritical, "Error" Exit Sub Next i Next NodeEnCurs MsgBox " L'Algorisme Incremental ha Convergit Correctament", vbinformation Sub AssignaResultatsNewtonModificatAVectorOrdinari() Dim NodeEnCurs As Long Dim MaxNode As Long Dim i As Integer 'm_vectordecorrimentsnodalsglobal = _ m_vectordecorrimentsnodalsglobalnewton MaxNode = m_colnode.count For NodeEnCurs = 1 To MaxNode m_colnode.item(nodeencurs).vectorcorriments(1) = m_vectordecorrimentsnodalsglobalnewton((nodeencurs - 1) * 6 + 1, 1) m_colnode.item(nodeencurs).vectorcorriments(2) = m_vectordecorrimentsnodalsglobalnewton((nodeencurs - 1) * 6 + 2, 1) m_colnode.item(nodeencurs).vectorcorriments(3) = m_vectordecorrimentsnodalsglobalnewton((nodeencurs - 1) * 6 + 3, 1) m_colnode.item(nodeencurs).vectorcorriments(4) = m_vectordecorrimentsnodalsglobalnewton((nodeencurs - 1) * 6 + 4, 1) m_colnode.item(nodeencurs).vectorcorriments(5) = m_vectordecorrimentsnodalsglobalnewton((nodeencurs - 1) * 6 + 5, 1) m_colnode.item(nodeencurs).vectorcorriments(6) = m_vectordecorrimentsnodalsglobalnewton((nodeencurs - 1) * 6 + 6, 1) For i = 1 To 6 If Abs(m_ColNode.Item(NodeEnCurs).VectorCorriments(i)) > # Then MsgBox "L'algorisme de Newton-Raphson no ha convergit", vbcritical, "Error" Exit Sub Next i Next NodeEnCurs MsgBox " L'algorisme de Newton-Raphson ha convergit", vbinformation Property Get MatriuDeRigidesaReduidaNewtonInvertida() As Variant

16 16 MatriuDeRigidesaReduidaNewtonInvertida = m_matriuderigidesareduidanewtoninvertida Sub RessetejaMatriuDeRigidesaGlobalNewton(EsElPrimerCop As Long) Dim i As Long, imax As Long Dim j As Long, jmax As Long If EsElPrimerCop = 1 Then Exit Sub i = UBound(m_MatriuDeRigidesaGlobalNewton, 1) j = UBound(m_MatriuDeRigidesaGlobalNewton, 2) For i = 1 To imax For j = 1 To jmax m_matriuderigidesaglobalnewton(i, j) = 0 Next j Next i Property Get VectorDeCarregaReduitNewton() As Variant VectorDeCarregaReduitNewton = m_vectordecarregareduitnewton 'Property Get VectorDeCorrimentsNodalsReduitParcialNewton() As Variant 'VectorDeCorrimentsNodalsReduitParcialNewton = m_vectordecorrimentsnodalsreduitparcialnewton ' Property Get VectorDeCorrimentsNodalsReduitNewton() As Variant VectorDeCorrimentsNodalsReduitNewton = m_vectordecorrimentsnodalsreduitnewton Property Get MatriuDeRigidesaGlobalNewton() As Variant MatriuDeRigidesaGlobalNewton = m_matriuderigidesaglobalnewton Property Get MatriuDeRigidesaReduidaNewton() As Variant MatriuDeRigidesaReduidaNewton = m_matriuderigidesareduidanewton Public Property Get MatriuDeRigidesaReduidaInvertida() As Variant MatriuDeRigidesaReduidaInvertida = m_matriuderigidesareduidainvertida Sub CalculaMatriuDeRigidesaReduidaInvertida() Dim A() As Double Dim B() As Double Dim Omega() As Double Dim Enter As Integer Dim CasEnCurs As Integer

17 17 CasEnCurs = Estructura.ModalitatDeCalculLineal 'CalculaMatriuDeRigidesaReduida AssignaVariantAMatriu m_matriuderigidesareduida, B If m_matriuderigidesareduidainvertidacalculada = False Then Select Case CasEnCurs Case 1 Estructura.InversaMatriuLU B, A Case 2 InversaMatPereSancho B, A 'BY PERE SANCHO Case 3 InversaMatriu B, A 'BY TON (M'HI VAIG PASSAR UN PUTO MATÍ SENCER) Case 4 Estructura.InversaMatriuLU B, A End Select m_matriuderigidesareduidainvertidacalculada = True m_matriuderigidesareduidainvertida = A Sub CalculaMatriuDeRigidesaReduidaMillorada() Dim MaxNodes As Long Dim ElemEnCurs As Long Dim j As Long Dim i As Long Dim MatriuEliminadora() As Double Dim MatriuEliminadoraTrasposta() As Double Dim MatriuAuxiliar() As Double Dim MG() As Double Dim m() As Double Dim FilesiColumnesAEsborrar() As Double Dim FilaEnCurs As Long, ColumnaEnCurs As Long If m_matriuderigidesareduidacalculada = True Then Exit Sub MaxNodes = m_colnode.count 'Primer hem de calcular les dimensions de la matriu de rigidesa _ reduïda j = 0 For ElemEnCurs = 1 To MaxNodes

18 18 'Principal.ProgresText 0, 0, "Etapa 1/3", True If m_colnode.item((elemencurs)).xblock = False Then j = j + 1 If m_colnode.item((elemencurs)).yblock = False Then j = j + 1 If m_colnode.item((elemencurs)).zblock = False Then j = j + 1 If m_colnode.item((elemencurs)).mxblock = False Then j = j + 1 If m_colnode.item((elemencurs)).myblock = False Then j = j + 1 If m_colnode.item((elemencurs)).mzblock = False Then j = j + 1 'ProgresMatriuDeRigidesaReduida1 MaxNodes, ElemEncurs Next ElemEnCurs 'FiProgres ReDim m(1 To j, 1 To j) ReDim FilesiColumnesAEsborrar(1 To MaxNodes * 6 - j) j = 1 For ElemEnCurs = MaxNodes To 1 Step -1 'Principal.ProgresText 0, 0, "Etapa 2/3", True If m_colnode.item((elemencurs)).mzblock = True Then FilesiColumnesAEsborrar(j) = ((ElemEnCurs - 1) * 6 + 6) j = j + 1 If m_colnode.item((elemencurs)).myblock = True Then FilesiColumnesAEsborrar(j) = ((ElemEnCurs - 1) * 6 + 5) j = j + 1 If m_colnode.item((elemencurs)).mxblock = True Then FilesiColumnesAEsborrar(j) = ((ElemEnCurs - 1) * 6 + 4) j = j + 1

19 19 If m_colnode.item((elemencurs)).zblock = True Then FilesiColumnesAEsborrar(j) = ((ElemEnCurs - 1) * 6 + 3) j = j + 1 If m_colnode.item((elemencurs)).yblock = True Then FilesiColumnesAEsborrar(j) = ((ElemEnCurs - 1) * 6 + 2) j = j + 1 If m_colnode.item((elemencurs)).xblock = True Then FilesiColumnesAEsborrar(j) = ((ElemEnCurs - 1) * 6 + 1) j = j + 1 Next ElemEnCurs 'ARA COMENÇA LA CHICHA... i = 0 j = 0 AssignaVariantAMatriu m_matriuderigidesaglobal, MG For FilaEnCurs = 1 To MaxNodes * 6 If Estructura.ExisteixElNumeroEnElVector(FilesiColumnesAEsborrar, (FilaEnCurs)) = False Then i = i + 1 Else GoTo Etiqueta2 For ColumnaEnCurs = 1 To MaxNodes * 6 If Estructura.ExisteixElNumeroEnElVector(FilesiColumnesAEsborrar, (ColumnaEnCurs)) = False Then j = j + 1 m(i, j) = MG(FilaEnCurs, ColumnaEnCurs) Else Etiqueta1: Next ColumnaEnCurs j = 0 'Principal.ProgresText MaxNodes * 6, FilaEnCurs, "Element ", True Etiqueta2: Next FilaEnCurs

20 20 m_matriuderigidesareduida = m m_matriuderigidesareduidacalculada = True Sub CalculaMatriuDeRigidesaReduidaNewtonInvertida() Dim A() As Double Dim B() As Double AssignaVariantAMatriu m_matriuderigidesareduidanewton, B 'InversaMatriu b, a InversaMatPereSancho B, A m_matriuderigidesareduidanewtoninvertida = A 'm_matriuderigidesareduidanewtoninvertida = b Sub CalculaMatriuDeRigidesaReduidaNewtonModificatInvertida() Dim A() As Double Dim B() As Double AssignaVariantAMatriu m_matriuderigidesareduidanewton, B 'InversaMatriu b, a InversaMatPereSancho B, A m_matriuderigidesareduidanewtonmodificatinvertida = A 'm_matriuderigidesareduidanewtonmodificatinvertida = b Public Property Let MatriuDeRigidesaReduidaCalculada(ByVal Val As Boolean) m_matriuderigidesareduidacalculada = Val Public Property Get MatriuDeRigidesaReduidaCalculada() As Boolean MatriuDeRigidesaReduidaCalculada = m_matriuderigidesareduidacalculada Property Get VectorsDeCorrimentsElementals(vntIndClave As Variant) As Variant If m_vectorsdecorrimentselementalscalculats = False Then CalculaVectorsDeCorrimentsElementals m_vectorsdecorrimentselementalscalculats = True VectorsDeCorrimentsElementals = m_vectorsdecorrimentselementals(vntindclave) Property Get MatriuDeConnectivitat() As Variant

21 21 Dim MaxElems As Long Dim ElemEnCurs As Long Dim MC() As Long MaxElems = m_colbarra.count ReDim MC(1 To MaxElems, 1 To 2) 'Construeix la matriu de connectivitat For ElemEnCurs = 1 To MaxElems MC(ElemEnCurs, 1) = _ m_colbarra.item(elemencurs).numnode1 MC(ElemEnCurs, 2) = _ m_colbarra.item(elemencurs).numnode2 Next ElemEnCurs MatriuDeConnectivitat = MC Property Get VectorDeCorrimentsNodalsGlobal() As Variant CalculaVectorDeCorrimentsNodalsGlobal VectorDeCorrimentsNodalsGlobal = m_vectordecorrimentsnodalsglobal Property Get VectorDeCorrimentsNodalsReduit() As Variant CalculaVectorDeCorrimentsNodalsReduit VectorDeCorrimentsNodalsReduit = m_vectordecorrimentsnodalsreduit Property Get MatriuDeRigidesaReduida() As Variant 'CalculaMatriuDeRigidesaReduida MatriuDeRigidesaReduida = m_matriuderigidesareduida Property Get MatriuDeCarregaGlobal() As Variant MatriuDeCarregaGlobal = m_matriudecarregaglobal Property Let MatriuDeCarregaGlobal(Val As Variant) m_matriudecarregaglobal = Val Property Get MatriuDeCarregaReduida() As Variant MatriuDeCarregaReduida = m_matriudecarregareduida Property Let MatriuDeCarregaReduida(Val As Variant) m_matriudecarregareduida = Val Property Get MatriuDeCorrimentsNodalsReduit() As Variant

22 22 MatriuDeCorrimentsNodalsReduit = m_matriudecorrimentsnodalsreduit Property Let MatriuDeCorrimentsNodalsReduit(Val As Variant) m_matriudecorrimentsnodalsreduit = Val Property Get MatriudeRigidesaGlobal() As Variant 'CalculaMatriuDeRigidesaGlobal MatriudeRigidesaGlobal = m_matriuderigidesaglobal Public Property Get ColBarra() As CColBarra 'Propietat per accedir a la col.lecció de barres Set ColBarra = m_colbarra Public Property Let ColBarra(ByVal Val As CColBarra) Set m_colbarra = Val Public Property Get ColNode() As CColNode 'Propietat per accedir a la col.lecció de barres Set ColNode = m_colnode Public Property Let ColNode(ByVal Val As CColNode) Set m_colnode = Val Public Property Let Nom(ByVal Val As String) 'Propietat per assignar la dada nom m_nom = Val Public Property Get Nom() As String Nom = m_nom Private Sub Class_Initialize() Set m_colbarra = New CColBarra Set m_colnode = New CColNode Set m_vectorsdecorrimentselementals = New Collection Set m_vectorsdecarregueselementals = New Collection Set m_colaccio = New CColAccio Private Sub Class_Terminate() Set m_colbarra = Nothing Set m_colnode = Nothing Set m_colaccio = Nothing Public Sub CalculaMatriuDeRigidesaGlobal() Dim K() As Double Dim MaxElemsNodes As Long Dim MaxElemsBarres As Long Dim ElemEnCurs As Long Dim NumNode1 As Long

23 23 Dim NumNode2 As Long Dim i As Integer Dim j As Integer Static Calculat As Boolean If Calculat = True Then Exit Sub MaxElemsNodes = m_colnode.count MaxElemsBarres = m_colbarra.count ReDim K(1 To MaxElemsNodes * 6, 1 To MaxElemsNodes * 6) For ElemEnCurs = 1 To MaxElemsBarres NumNode1 = m_colbarra.item((elemencurs)).numnode1 NumNode2 = m_colbarra.item((elemencurs)).numnode2 For i = 1 To 6 For j = 1 To 6 K((NumNode1-1) * 6 _ + i, (NumNode1-1) * 6 + j) = _ K((NumNode1-1) * 6 _ + i, (NumNode1-1) * 6 + j) + _ m_colbarra.item((elemencurs)).matriuderigidesaglobalelemental(i, j) Next j Next i For i = 1 To 6 For j = 1 To 6 K((NumNode2-1) * 6 _ + i, (NumNode2-1) * 6 + j) = _ K((NumNode2-1) * 6 _ + i, (NumNode2-1) * 6 + j) + _ m_colbarra.item((elemencurs)).matriuderigidesaglobalelemental(i + 6, j + 6) Next j Next i For i = 1 To 6 For j = 1 To 6 K((NumNode1-1) * 6 _ + i, (NumNode2-1) * 6 + j) = _ K((NumNode1-1) * 6 _

24 24 + i, (NumNode2-1) * 6 + j) + _ m_colbarra.item((elemencurs)).matriuderigidesaglobalelemental(i, j + 6) Next j Next i For i = 1 To 6 For j = 1 To 6 K((NumNode2-1) * 6 _ + i, (NumNode1-1) * 6 + j) = _ K((NumNode2-1) * 6 _ + i, (NumNode1-1) * 6 + j) + _ m_colbarra.item((elemencurs)).matriuderigidesaglobalelemental(i + 6, j) Next j Next i 'Principal.ProgresText MaxElemsBarres, ElemEnCurs, "Element :", True Principal.Escrollant MaxElemsBarres, ElemEnCurs, "Construcció Matriu de rigidesa global" Next ElemEnCurs 'FiProgres m_matriuderigidesaglobal = K Calculat = True Sub CalculaMatriuDeRigidesaGlobalNewton() Dim K() As Double Dim MaxElemsNodes As Long Dim MaxElemsBarres As Long Dim ElemEnCurs As Long Dim NumNode1 As Long Dim NumNode2 As Long Dim i As Integer Dim j As Integer MaxElemsNodes = m_colnode.count MaxElemsBarres = m_colbarra.count ReDim K(1 To MaxElemsNodes * 6, 1 To MaxElemsNodes * 6) For ElemEnCurs = 1 To MaxElemsBarres NumNode1 = m_colbarra.item((elemencurs)).numnode1 NumNode2 = m_colbarra.item((elemencurs)).numnode2

25 25 For i = 1 To 6 For j = 1 To 6 j) K((NumNode1-1) * 6 _ + i, (NumNode1-1) * 6 + j) = _ K((NumNode1-1) * 6 _ + i, (NumNode1-1) * 6 + j) + _ m_colbarra.item((elemencurs)).matriuderigidesaglobalelementalnewton(i, Next j Next i For i = 1 To 6 For j = 1 To 6 6, j + 6) K((NumNode2-1) * 6 _ + i, (NumNode2-1) * 6 + j) = _ K((NumNode2-1) * 6 _ + i, (NumNode2-1) * 6 + j) + _ m_colbarra.item((elemencurs)).matriuderigidesaglobalelementalnewton(i + Next j Next i For i = 1 To 6 For j = 1 To 6 + 6) K((NumNode1-1) * 6 _ + i, (NumNode2-1) * 6 + j) = _ K((NumNode1-1) * 6 _ + i, (NumNode2-1) * 6 + j) + _ m_colbarra.item((elemencurs)).matriuderigidesaglobalelementalnewton(i, j Next j Next i For i = 1 To 6 For j = 1 To 6 K((NumNode2-1) * 6 _ + i, (NumNode1-1) * 6 + j) = _ K((NumNode2-1) * 6 _ + i, (NumNode1-1) * 6 + j) + _

26 26 6, j) m_colbarra.item((elemencurs)).matriuderigidesaglobalelementalnewton(i + Next j Next i Next ElemEnCurs 'FiProgres m_matriuderigidesaglobalnewton = K Private Sub CalculaVectorDeCarregaGlobalXXXXXXXXXXXXX() Dim MaxElemsNodes As Long Dim MaxElemsBarres As Long Dim ElemEnCurs As Long Dim NumNode1 As Long Dim NumNode2 As Long Dim V() As Double Static Calculat As Double If Calculat = True Then Exit Sub MaxElemsNodes = m_colnode.count MaxElemsBarres = m_colbarra.count ReDim V(1 To MaxElemsNodes * 6, 1 To 1) 'OWA For ElemEnCurs = 1 To MaxElemsBarres NumNode1 = m_colbarra.item((elemencurs)).numnode1 NumNode2 = m_colbarra.item((elemencurs)).numnode2 V((NumNode1-1) * 6 + 1, 1) = _ V((NumNode1-1) * 6 + 1, 1) + _ m_colbarra.item((elemencurs)).carreguesnodalsequivalentsglobals(1) V((NumNode1-1) * 6 + 2, 1) = _

27 27 V((NumNode1-1) * 6 + 2, 1) + _ m_colbarra.item((elemencurs)).carreguesnodalsequivalentsglobals(2) V((NumNode1-1) * 6 + 3, 1) = _ V((NumNode1-1) * 6 + 3, 1) + _ m_colbarra.item((elemencurs)).carreguesnodalsequivalentsglobals(3) V((NumNode1-1) * 6 + 4, 1) = _ V((NumNode1-1) * 6 + 4, 1) + _ m_colbarra.item((elemencurs)).carreguesnodalsequivalentsglobals(4) V((NumNode1-1) * 6 + 5, 1) = _ V((NumNode1-1) * 6 + 5, 1) + _ m_colbarra.item((elemencurs)).carreguesnodalsequivalentsglobals(5) V((NumNode1-1) * 6 + 6, 1) = _ V((NumNode1-1) * 6 + 6, 1) + _ m_colbarra.item((elemencurs)).carreguesnodalsequivalentsglobals(6) V((NumNode2-1) * 6 + 1, 1) = _ V((NumNode2-1) * 6 + 1, 1) + _ m_colbarra.item((elemencurs)).carreguesnodalsequivalentsglobals(7) V((NumNode2-1) * 6 + 2, 1) = _ V((NumNode2-1) * 6 + 2, 1) + _ m_colbarra.item((elemencurs)).carreguesnodalsequivalentsglobals(8) V((NumNode2-1) * 6 + 3, 1) = _ V((NumNode2-1) * 6 + 3, 1) + _ m_colbarra.item((elemencurs)).carreguesnodalsequivalentsglobals(9) V((NumNode2-1) * 6 + 4, 1) = _ V((NumNode2-1) * 6 + 4, 1) + _ m_colbarra.item((elemencurs)).carreguesnodalsequivalentsglobals(10) V((NumNode2-1) * 6 + 5, 1) = _ V((NumNode2-1) * 6 + 5, 1) + _ m_colbarra.item((elemencurs)).carreguesnodalsequivalentsglobals(11) V((NumNode2-1) * 6 + 6, 1) = _ V((NumNode2-1) * 6 + 6, 1) + _ m_colbarra.item((elemencurs)).carreguesnodalsequivalentsglobals(12) Next ElemEnCurs 'Afegeixo els valors concentrats en els nodes For ElemEnCurs = 1 To MaxElemsNodes V((ElemEnCurs - 1) * 6 + 1, 1) = _ V((ElemEnCurs - 1) * 6 + 1, 1) + _ m_colnode.item((elemencurs)).fx V((ElemEnCurs - 1) * 6 + 2, 1) = _

28 28 V((ElemEnCurs - 1) * 6 + 2, 1) + _ m_colnode.item((elemencurs)).fy V((ElemEnCurs - 1) * 6 + 3, 1) = _ V((ElemEnCurs - 1) * 6 + 3, 1) + _ m_colnode.item((elemencurs)).fz V((ElemEnCurs - 1) * 6 + 4, 1) = _ V((ElemEnCurs - 1) * 6 + 4, 1) + _ m_colnode.item((elemencurs)).momentx V((ElemEnCurs - 1) * 6 + 5, 1) = _ V((ElemEnCurs - 1) * 6 + 5, 1) + _ m_colnode.item((elemencurs)).momenty V((ElemEnCurs - 1) * 6 + 6, 1) = _ V((ElemEnCurs - 1) * 6 + 6, 1) + _ m_colnode.item((elemencurs)).momentz Next ElemEnCurs m_vectordecarregaglobal = V Calculat = False Sub CalculaMatriuDeCarregaGlobal() Dim MaxElemsNodes As Long Dim MaxElemsBarres As Long Dim ElemEnCurs As Long Dim NumNode1 As Long Dim NumNode2 As Long Dim m() As Double Static Calculat As Double Dim AccioEnCurs As Long, MaxAccio As Long If Calculat = True Then Exit Sub MaxElemsNodes = m_colnode.count MaxElemsBarres = m_colbarra.count MaxAccio = m_colaccio.count ReDim m(1 To MaxElemsNodes * 6, 1 To MaxAccio) For AccioEnCurs = 1 To MaxAccio 'OWA For ElemEnCurs = 1 To MaxElemsBarres NumNode1 = m_colbarra.item((elemencurs)).numnode1 NumNode2 = m_colbarra.item((elemencurs)).numnode2 m((numnode1-1) * 6 + 1, AccioEnCurs) = _

29 29 m((numnode1-1) * 6 + 1, AccioEnCurs) + _ m_colbarra.item((elemencurs)).carreguesnodalsequivalentsglobals(1, AccioEnCurs) _ + m_colbarra.item((elemencurs)).carreguesnodalsequivalentsapuntualsglobals(1, AccioEnCurs) m((numnode1-1) * 6 + 2, AccioEnCurs) = _ m((numnode1-1) * 6 + 2, AccioEnCurs) + _ m_colbarra.item((elemencurs)).carreguesnodalsequivalentsglobals(2, AccioEnCurs) _ + m_colbarra.item((elemencurs)).carreguesnodalsequivalentsapuntualsglobals(2, AccioEnCurs) m((numnode1-1) * 6 + 3, AccioEnCurs) = _ m((numnode1-1) * 6 + 3, AccioEnCurs) + _ m_colbarra.item((elemencurs)).carreguesnodalsequivalentsglobals(3, AccioEnCurs) _ + m_colbarra.item((elemencurs)).carreguesnodalsequivalentsapuntualsglobals(3, AccioEnCurs) m((numnode1-1) * 6 + 4, AccioEnCurs) = _ m((numnode1-1) * 6 + 4, AccioEnCurs) + _ m_colbarra.item((elemencurs)).carreguesnodalsequivalentsglobals(4, AccioEnCurs) _ + m_colbarra.item((elemencurs)).carreguesnodalsequivalentsapuntualsglobals(4, AccioEnCurs) m((numnode1-1) * 6 + 5, AccioEnCurs) = _ m((numnode1-1) * 6 + 5, AccioEnCurs) + _ m_colbarra.item((elemencurs)).carreguesnodalsequivalentsglobals(5, AccioEnCurs) _ + m_colbarra.item((elemencurs)).carreguesnodalsequivalentsapuntualsglobals(5, AccioEnCurs) m((numnode1-1) * 6 + 6, AccioEnCurs) = _ m((numnode1-1) * 6 + 6, AccioEnCurs) + _ m_colbarra.item((elemencurs)).carreguesnodalsequivalentsglobals(6, AccioEnCurs) _ + m_colbarra.item((elemencurs)).carreguesnodalsequivalentsapuntualsglobals(6, AccioEnCurs) m((numnode2-1) * 6 + 1, AccioEnCurs) = _ m((numnode2-1) * 6 + 1, AccioEnCurs) + _ m_colbarra.item((elemencurs)).carreguesnodalsequivalentsglobals(7, AccioEnCurs) _ + m_colbarra.item((elemencurs)).carreguesnodalsequivalentsapuntualsglobals(7, AccioEnCurs)

30 30 m((numnode2-1) * 6 + 2, AccioEnCurs) = _ m((numnode2-1) * 6 + 2, AccioEnCurs) + _ m_colbarra.item((elemencurs)).carreguesnodalsequivalentsglobals(8, AccioEnCurs) _ + m_colbarra.item((elemencurs)).carreguesnodalsequivalentsapuntualsglobals(8, AccioEnCurs) m((numnode2-1) * 6 + 3, AccioEnCurs) = _ m((numnode2-1) * 6 + 3, AccioEnCurs) + _ m_colbarra.item((elemencurs)).carreguesnodalsequivalentsglobals(9, AccioEnCurs) _ + m_colbarra.item((elemencurs)).carreguesnodalsequivalentsapuntualsglobals(9, AccioEnCurs) m((numnode2-1) * 6 + 4, AccioEnCurs) = _ m((numnode2-1) * 6 + 4, AccioEnCurs) + _ m_colbarra.item((elemencurs)).carreguesnodalsequivalentsglobals(10, AccioEnCurs) _ + m_colbarra.item((elemencurs)).carreguesnodalsequivalentsapuntualsglobals(10, AccioEnCurs) m((numnode2-1) * 6 + 5, AccioEnCurs) = _ m((numnode2-1) * 6 + 5, AccioEnCurs) + _ m_colbarra.item((elemencurs)).carreguesnodalsequivalentsglobals(11, AccioEnCurs) _ + m_colbarra.item((elemencurs)).carreguesnodalsequivalentsapuntualsglobals(11, AccioEnCurs) m((numnode2-1) * 6 + 6, AccioEnCurs) = _ m((numnode2-1) * 6 + 6, AccioEnCurs) + _ m_colbarra.item((elemencurs)).carreguesnodalsequivalentsglobals(12, AccioEnCurs) _ + m_colbarra.item((elemencurs)).carreguesnodalsequivalentsapuntualsglobals(12, AccioEnCurs) Next ElemEnCurs Next AccioEnCurs 'HI SUMO ELS VALORS DE LES CÀRREGUES NODALS For AccioEnCurs = 1 To MaxAccio For ElemEnCurs = 1 To MaxElemsNodes m((elemencurs - 1) * 6 + 1, AccioEnCurs) = _ m((elemencurs - 1) * 6 + 1, AccioEnCurs) + _ m_colnode.item((elemencurs)).colsolicitacionodal.item(accioencurs).fx m((elemencurs - 1) * 6 + 2, AccioEnCurs) = _

31 31 m((elemencurs - 1) * 6 + 2, AccioEnCurs) + _ m_colnode.item((elemencurs)).colsolicitacionodal.item(accioencurs).fy m((elemencurs - 1) * 6 + 3, AccioEnCurs) = _ m((elemencurs - 1) * 6 + 3, AccioEnCurs) + _ m_colnode.item((elemencurs)).colsolicitacionodal.item(accioencurs).fz m((elemencurs - 1) * 6 + 4, AccioEnCurs) = _ m((elemencurs - 1) * 6 + 4, AccioEnCurs) + _ m_colnode.item((elemencurs)).colsolicitacionodal.item(accioencurs).mx m((elemencurs - 1) * 6 + 5, AccioEnCurs) = _ m((elemencurs - 1) * 6 + 5, AccioEnCurs) + _ m_colnode.item((elemencurs)).colsolicitacionodal.item(accioencurs).my m((elemencurs - 1) * 6 + 6, AccioEnCurs) = _ m((elemencurs - 1) * 6 + 6, AccioEnCurs) + _ m_colnode.item((elemencurs)).colsolicitacionodal.item(accioencurs).mz Next ElemEnCurs Next AccioEnCurs m_matriudecarregaglobal = m Calculat = False Sub CalculaMatriuDeCarregaReduida() Dim MaxNodes As Long Dim MaxElems As Long Dim ElemEnCurs As Long Dim V() As Double Dim B() As Double Dim i As Long Dim j As Long Dim MaxAccio As Long, AccioEnCurs As Long Dim MatriuEliminadora() As Double MaxNodes = m_colnode.count 'MaxElems = m_colbarra.count MaxAccio = m_colaccio.count 'CREACIÓ DE LA MATRIU ELIMINADORA '0 si s'ha d'eliminar i 1 en cas contrari 'Per eliminar, per exemple la segona i tercera fileres 'es tindria (0 1 0) ' (0 0 1) 'Primer hem de calcular quants elements tindrà la matriu eliminadora j = 0 For ElemEnCurs = 1 To MaxNodes

32 32 If m_colnode.item((elemencurs)).xblock = False Then j = j + 1 If m_colnode.item((elemencurs)).yblock = False Then j = j + 1 If m_colnode.item((elemencurs)).zblock = False Then j = j + 1 If m_colnode.item((elemencurs)).mxblock = False Then j = j + 1 If m_colnode.item((elemencurs)).myblock = False Then j = j + 1 If m_colnode.item((elemencurs)).mzblock = False Then j = j + 1 'ProgresMatriuDeRigidesaReduida1 MaxNodes, ElemEncurs Next ElemEnCurs 'FiProgres 'ReDim M(1 To j, 1 To j) ReDim MatriuEliminadora(1 To j, 1 To MaxNodes * 6) j = 1 For ElemEnCurs = 1 To MaxNodes If m_colnode.item((elemencurs)).xblock = False Then MatriuEliminadora(j, (ElemEnCurs - 1) * 6 + 1) = 1 j = j + 1 If m_colnode.item((elemencurs)).yblock = False Then MatriuEliminadora(j, (ElemEnCurs - 1) * 6 + 2) = 1 j = j + 1 If m_colnode.item((elemencurs)).zblock = False Then MatriuEliminadora(j, (ElemEnCurs - 1) * 6 + 3) = 1 j = j + 1 If m_colnode.item((elemencurs)).mxblock = False Then MatriuEliminadora(j, (ElemEnCurs - 1) * 6 + 4) = 1 j = j + 1 If m_colnode.item((elemencurs)).myblock = False Then MatriuEliminadora(j, (ElemEnCurs - 1) * 6 + 5) = 1

33 33 j = j + 1 If m_colnode.item((elemencurs)).mzblock = False Then MatriuEliminadora(j, (ElemEnCurs - 1) * 6 + 6) = 1 j = j + 1 Next ElemEnCurs Principal.AssignaVariantAMatriu m_matriudecarregaglobal, B Principal.MultiplicaMatrius MatriuEliminadora, B, V m_matriudecarregareduida = V Sub CalculaVectorDeCarregaReduitXXXXXXXXXXXXXXXXXX() Dim MaxNodes As Long Dim MaxElems As Long Dim ElemEnCurs As Long Dim V() As Double Dim i As Long Dim j As Long MaxNodes = m_colnode.count MaxElems = m_colbarra.count 'Primer hem de calcular quants elements tindrà el vector de _ càrrega reduït j = 0 For ElemEnCurs = 1 To MaxNodes If m_colnode.item((elemencurs)).xblock = False Then j = j + 1 If m_colnode.item((elemencurs)).yblock = False Then j = j + 1 If m_colnode.item((elemencurs)).zblock = False Then j = j + 1 If m_colnode.item((elemencurs)).mxblock = False Then j = j + 1 If m_colnode.item((elemencurs)).myblock = False Then j = j + 1

34 34 If m_colnode.item((elemencurs)).mzblock = False Then j = j + 1 Next ElemEnCurs ReDim V(1 To j, 1 To 1) For ElemEnCurs = 1 To MaxElems m_colnode.item(str(m_colbarra.item((elemencurs)).numnode1)).fx = _ m_colnode.item(str(m_colbarra.item((elemencurs)).numnode1)).fx + _ m_colbarra.item((elemencurs)).carreguesnodalsequivalentsglobals(1) m_colnode.item(str(m_colbarra.item((elemencurs)).numnode1)).fy = _ m_colnode.item(str(m_colbarra.item((elemencurs)).numnode1)).fy + _ m_colbarra.item((elemencurs)).carreguesnodalsequivalentsglobals(2) m_colnode.item(str(m_colbarra.item((elemencurs)).numnode1)).fz = _ m_colnode.item(str(m_colbarra.item((elemencurs)).numnode1)).fz + _ m_colbarra.item((elemencurs)).carreguesnodalsequivalentsglobals(3) m_colnode.item(str(m_colbarra.item((elemencurs)).numnode1)).momentx = _ m_colnode.item(str(m_colbarra.item((elemencurs)).numnode1)).momentx + _ m_colbarra.item((elemencurs)).carreguesnodalsequivalentsglobals(4) m_colnode.item(str(m_colbarra.item((elemencurs)).numnode1)).momenty = _ m_colnode.item(str(m_colbarra.item((elemencurs)).numnode1)).momenty + _ m_colbarra.item((elemencurs)).carreguesnodalsequivalentsglobals(5) m_colnode.item(str(m_colbarra.item((elemencurs)).numnode1)).momentz = _ m_colnode.item(str(m_colbarra.item((elemencurs)).numnode1)).momentz + _ m_colbarra.item((elemencurs)).carreguesnodalsequivalentsglobals(6) m_colnode.item(str(m_colbarra.item((elemencurs)).numnode2)).fx = _ m_colnode.item(str(m_colbarra.item((elemencurs)).numnode2)).fx + _ m_colbarra.item((elemencurs)).carreguesnodalsequivalentsglobals(7) m_colnode.item(str(m_colbarra.item((elemencurs)).numnode2)).fy = _ m_colnode.item(str(m_colbarra.item((elemencurs)).numnode2)).fy + _ m_colbarra.item((elemencurs)).carreguesnodalsequivalentsglobals(8) m_colnode.item(str(m_colbarra.item((elemencurs)).numnode2)).fz = _ m_colnode.item(str(m_colbarra.item((elemencurs)).numnode2)).fz + _ m_colbarra.item((elemencurs)).carreguesnodalsequivalentsglobals(9) m_colnode.item(str(m_colbarra.item((elemencurs)).numnode2)).momentx = _ m_colnode.item(str(m_colbarra.item((elemencurs)).numnode2)).momentx + _

35 35 m_colbarra.item((elemencurs)).carreguesnodalsequivalentsglobals(10) m_colnode.item(str(m_colbarra.item((elemencurs)).numnode2)).momenty = _ m_colnode.item(str(m_colbarra.item((elemencurs)).numnode2)).momenty + _ m_colbarra.item((elemencurs)).carreguesnodalsequivalentsglobals(11) m_colnode.item(str(m_colbarra.item((elemencurs)).numnode2)).momentz = _ m_colnode.item(str(m_colbarra.item((elemencurs)).numnode2)).momentz + _ m_colbarra.item((elemencurs)).carreguesnodalsequivalentsglobals(12) Next ElemEnCurs j = 1 For ElemEnCurs = 1 To MaxNodes If m_colnode.item((elemencurs)).xblock = False Then V(j, 1) = _ m_colnode.item((elemencurs)).fx j = j + 1 If m_colnode.item((elemencurs)).yblock = False Then V(j, 1) = _ m_colnode.item((elemencurs)).fy j = j + 1 If m_colnode.item((elemencurs)).zblock = False Then V(j, 1) = _ m_colnode.item((elemencurs)).fz j = j + 1 If m_colnode.item((elemencurs)).mxblock = False Then V(j, 1) = _ m_colnode.item((elemencurs)).momentx j = j + 1 If m_colnode.item((elemencurs)).myblock = False Then V(j, 1) = _ m_colnode.item((elemencurs)).momenty j = j + 1 If m_colnode.item((elemencurs)).mzblock = False Then V(j, 1) = _ m_colnode.item((elemencurs)).momentz j = j + 1

36 36 Next ElemEnCurs m_vectordecarregareduit = V Public Sub CalculaMatriuDeRigidesaReduida() Dim MaxNodes As Long Dim ElemEnCurs As Long Dim j As Long Dim i As Long Dim MatriuEliminadora() As Double Dim MatriuEliminadoraTrasposta() As Double Dim MatriuAuxiliar() As Double Dim MG() As Double Dim m() As Double 'CREACIÓ DE LA MATRIU ELIMINADORA '0 si s'ha d'eliminar i 1 en cas contrari 'Per eliminar, per exemple la segona i tercera fileres 'es tindria (0 1 0) ' (0 0 1) If m_matriuderigidesareduidacalculada = True Then Exit Sub MaxNodes = m_colnode.count 'CalculaMatriuDeRigidesaGlobal ' per si no estava calculada 'Primer hem de calcular les dimensions de la matriu de rigidesa _ reduïda j = 0 For ElemEnCurs = 1 To MaxNodes If m_colnode.item((elemencurs)).xblock = False Then j = j + 1 If m_colnode.item((elemencurs)).yblock = False Then j = j + 1 If m_colnode.item((elemencurs)).zblock = False Then j = j + 1

37 37 If m_colnode.item((elemencurs)).mxblock = False Then j = j + 1 If m_colnode.item((elemencurs)).myblock = False Then j = j + 1 If m_colnode.item((elemencurs)).mzblock = False Then j = j + 1 'ProgresMatriuDeRigidesaReduida1 MaxNodes, ElemEncurs Next ElemEnCurs 'FiProgres ReDim m(1 To j, 1 To j) ReDim MatriuEliminadora(1 To j, 1 To MaxNodes * 6) 'Ara posem els elements a la matriu eliminadora j = 1 For ElemEnCurs = 1 To MaxNodes If m_colnode.item((elemencurs)).xblock = False Then MatriuEliminadora(j, (ElemEnCurs - 1) * 6 + 1) = 1 j = j + 1 If m_colnode.item((elemencurs)).yblock = False Then MatriuEliminadora(j, (ElemEnCurs - 1) * 6 + 2) = 1 j = j + 1 If m_colnode.item((elemencurs)).zblock = False Then MatriuEliminadora(j, (ElemEnCurs - 1) * 6 + 3) = 1 j = j + 1 If m_colnode.item((elemencurs)).mxblock = False Then MatriuEliminadora(j, (ElemEnCurs - 1) * 6 + 4) = 1 j = j + 1 If m_colnode.item((elemencurs)).myblock = False Then MatriuEliminadora(j, (ElemEnCurs - 1) * 6 + 5) = 1 j = j + 1 If m_colnode.item((elemencurs)).mzblock = False Then MatriuEliminadora(j, (ElemEnCurs - 1) * 6 + 6) = 1 j = j + 1 'ProgresMatriuDeRigidesaReduida2 MaxNodes, ElemEncurs Next ElemEnCurs

38 38 'FiProgres 'Faig la multiplicació que em donarà el resultat final TrasposaMatriu MatriuEliminadora, MatriuEliminadoraTrasposta AssignaVariantAMatriu m_matriuderigidesaglobal, MG MultiplicaMatrius MatriuEliminadora, MG, _ MatriuAuxiliar MultiplicaMatrius MatriuAuxiliar, MatriuEliminadoraTrasposta, _ m m_matriuderigidesareduida = m m_matriuderigidesareduidacalculada = True Sub CalculaMatriuDeRigidesaReduidaNewton() Dim MaxNodes As Long Dim ElemEnCurs As Long Dim j As Long Dim i As Long Dim MatriuEliminadora() As Double Dim MatriuEliminadoraTrasposta() As Double Dim MatriuAuxiliar() As Double Dim MG() As Double Dim m() As Double '0 si s'ha d'eliminar i 1 en cas contrari 'Per eliminar, per exemple la segona i tercera fileres 'es tindria (0 1 0) ' (0 0 1) MaxNodes = m_colnode.count j = 0 For ElemEnCurs = 1 To MaxNodes If m_colnode.item((elemencurs)).xblock = False Then j = j + 1 If m_colnode.item((elemencurs)).yblock = False Then j = j + 1 If m_colnode.item((elemencurs)).zblock = False Then j = j + 1

39 39 If m_colnode.item((elemencurs)).mxblock = False Then j = j + 1 If m_colnode.item((elemencurs)).myblock = False Then j = j + 1 If m_colnode.item((elemencurs)).mzblock = False Then j = j + 1 ProgresMatriuDeRigidesaReduida1 MaxNodes, ElemEnCurs Next ElemEnCurs 'FiProgres ReDim m(1 To j, 1 To j) ReDim MatriuEliminadora(1 To j, 1 To MaxNodes * 6) 'Ara posem els elements a la matriu eliminadora j = 1 For ElemEnCurs = 1 To MaxNodes If m_colnode.item((elemencurs)).xblock = False Then MatriuEliminadora(j, (ElemEnCurs - 1) * 6 + 1) = 1 j = j + 1 If m_colnode.item((elemencurs)).yblock = False Then MatriuEliminadora(j, (ElemEnCurs - 1) * 6 + 2) = 1 j = j + 1 If m_colnode.item((elemencurs)).zblock = False Then MatriuEliminadora(j, (ElemEnCurs - 1) * 6 + 3) = 1 j = j + 1 If m_colnode.item((elemencurs)).mxblock = False Then MatriuEliminadora(j, (ElemEnCurs - 1) * 6 + 4) = 1 j = j + 1 If m_colnode.item((elemencurs)).myblock = False Then MatriuEliminadora(j, (ElemEnCurs - 1) * 6 + 5) = 1 j = j + 1 If m_colnode.item((elemencurs)).mzblock = False Then MatriuEliminadora(j, (ElemEnCurs - 1) * 6 + 6) = 1 j = j + 1 'ProgresMatriuDeRigidesaReduida2 MaxNodes, ElemEncurs

40 40 Next ElemEnCurs 'FiProgres 'Faig la multiplicació que em donarà el resultat final TrasposaMatriu MatriuEliminadora, MatriuEliminadoraTrasposta AssignaVariantAMatriu m_matriuderigidesaglobalnewton, MG MultiplicaMatrius MatriuEliminadora, MG, _ MatriuAuxiliar MultiplicaMatrius MatriuAuxiliar, MatriuEliminadoraTrasposta, _ m m_matriuderigidesareduidanewton = m Private Sub AssignaVariantAMatriu(Origen As Variant, Desti() As Double) Dim imin As Long Dim imax As Long Dim jmin As Long Dim jmax As Long Dim i As Long Dim j As Long imax = UBound(Origen, 1) imin = LBound(Origen, 1) jmax = UBound(Origen, 2) jmin = LBound(Origen, 2) ReDim Desti(imin To imax, jmin To jmax) For i = imin To imax For j = jmin To jmax Desti(i, j) = Origen(i, j) Next j ' Principal.ProgresText imax, i, " Assigna a Matriu : ", True Next i Public Sub AssignaMatriuAMatriu(Origen() As Double, Desti() As Double) Dim imin As Long Dim imax As Long

41 41 Dim jmin As Long Dim jmax As Long Dim i As Long Dim j As Long imax = UBound(Origen, 1) imin = LBound(Origen, 1) jmax = UBound(Origen, 2) jmin = LBound(Origen, 2) ReDim Desti(imin To imax, jmin To jmax) For i = imin To imax For j = jmin To jmax Desti(i, j) = Origen(i, j) Next j Next i Private Sub TrasposaMatriu(MatriuOrigen() As Double, MatriuDesti() As Double) Dim i As Integer Dim j As Integer Dim Origen() As Double ReDim Origen(1 To UBound(MatriuOrigen, 1), 1 To UBound(MatriuOrigen, 2)) ReDim MatriuDesti(1 To UBound(Origen, 2), 1 To UBound(Origen, 1)) For i = 1 To UBound(Origen, 1) For j = 1 To UBound(Origen, 2) Origen(i, j) = MatriuOrigen(i, j) Next j Next i For i = LBound(Origen, 1) To UBound(Origen, 1) For j = LBound(Origen, 2) To UBound(Origen, 2) MatriuDesti(j, i) = Origen(i, j) Next j Next i Private Sub InversaMatriu(MatriuOrigen() As Double, MatriuDesti() As Double)

42 42 Dim Ix As Long Dim jx As Long Dim ix2 As Long Dim jx2 As Long Dim A As Long Dim Apuja As Long Dim Trobat As Boolean Dim Matraux() As Double Dim Ks As Double Dim Ki As Double Dim aux1 As Double Dim aux2 As Double Dim p As Long Dim Rang As Long Dim i As Long Dim j As Long Dim Origen() As Double 'On Error GoTo Error Rang = UBound(MatriuOrigen, 1) ReDim Matraux(1 To Rang, 1 To Rang) ReDim Origen(1 To Rang, 1 To Rang) ReDim MatriuDesti(1 To Rang, 1 To Rang) 'Assignem a la matriu Origen la matriu MatriuOrigen per _ evitar modificar la matriu d'entrada For i = 1 To Rang For j = 1 To Rang Origen(i, j) = MatriuOrigen(i, j) Next j Next i 'crea la matriu identitat nxn For Ix = 1 To Rang For jx = 1 To Rang If Ix = jx Then Matraux(Ix, jx) = 1 Next jx Next Ix