Пример #1
0
        // ....Compute element Matrix.....................................
        public override double[,] ComputeMatrix()
        {
            double[] gaussCoord = { -1 / Math.Sqrt(3), 1 / Math.Sqrt(3) };
            if (Modell.Material.TryGetValue(ElementMaterialId, out var abstractMaterial))
            {
            }
            material        = (Material)abstractMaterial;
            ElementMaterial = material ?? throw new ArgumentNullException(nameof(material));
            var conductivity = material.Leitfähigkeit[0];

            MatrizenAlgebra.Clear(elementMatrix);
            foreach (var coor1 in gaussCoord)
            {
                foreach (var coor2 in gaussCoord)
                {
                    var z0 = coor1;
                    var z1 = coor2;
                    ComputeGeometry(z0, z1);
                    Sx = ComputeSx(z0, z1);
                    // Ke = C*Sx*SxT*determinant
                    MatrizenAlgebra.MultAddMatrixTransposed(elementMatrix, Determinant * conductivity, Sx, Sx);
                }
            }
            return(elementMatrix);
        }
Пример #2
0
        // ....Compute element matrix.....................................
        public override double[,] ComputeMatrix()
        {
            if (Modell.Material.TryGetValue(ElementMaterialId, out var abstractMaterial))
            {
            }
            Material = (Material)abstractMaterial;
            double[] gCoord  = { -1 / Math.Sqrt(5.0 / 3), 0, 1 / (Math.Sqrt(5.0 / 3)) };
            double[] gWeight = { (5.0 / 9), (8.0 / 9), (5.0 / 9) };                               // gaussian coordinates, weights
            _ = new double[8, 3];
            MatrizenAlgebra.Clear(elementMatrix);

            // material matrix für ebene Verzerrung (plane strain)
            var conduct = ((Material)Material)?.Leitfähigkeit;

            if (conduct != null)
            {
                e[0, 0] = conduct[0];
                e[1, 1] = conduct[1];
                e[2, 2] = conduct[2];
            }

            for (var i = 0; i < gCoord.Length; i++)
            {
                var z0 = gCoord[i]; var g0 = gWeight[i];
                for (var j = 0; j < gCoord.Length; j++)
                {
                    var z1 = gCoord[j]; var g1 = gWeight[j];
                    for (var k = 0; k < gCoord.Length; k++)
                    {
                        var z2 = gCoord[k]; var g2 = gWeight[k];
                        ComputeGeometry(z0, z1, z2);
                        Sx = ComputeSx(z0, z1, z2);
                        // Ke = determinant*g0*g1*g2*Sx*E*SxT
                        var temp = MatrizenAlgebra.Mult(Sx, e);
                        MatrizenAlgebra.MultAddMatrixTransposed(elementMatrix, Determinant * g0 * g1 * g2, temp, Sx);
                    }
                }
            }
            return(elementMatrix);
        }
Пример #3
0
 // ....Compute element matrix.....................................
 public override double[,] ComputeMatrix()
 {
     MatrizenAlgebra.Clear(elementMatrix);
     ComputeMaterial();
     for (var i = 0; i < GCoord.Length; i++)
     {
         z0 = GCoord[i]; g0 = GWeight[i];
         for (var j = 0; j < GCoord.Length; j++)
         {
             z1 = GCoord[j]; g1 = GWeight[j];
             for (var k = 0; k < GCoord.Length; k++)
             {
                 z2 = GCoord[k]; g2 = GWeight[k];
                 ComputeGeometry(z0, z1, z2);
                 Sx = ComputeSx(z0, z1, z2);
                 ComputeStrainDisplacementTransformation();
                 // Ke = determinant*g0*g1*g2*BT*E*B
                 var temp = MatrizenAlgebra.MultTransposedMatrix(Determinant * g0 * g1 * g2, b, e);
                 MatrizenAlgebra.MultAddMatrix(elementMatrix, temp, b);
             }
         }
     }
     return(elementMatrix);
 }