ISAAR.MSolve.FEM.Tests/Elements/Tri3.cs
/// <summary>
/// To reproduce the reference solution, run the Abaqus job tri3_test0.inp and look at 
/// TRI3_TEST0_STIFFNESS_MATRICES.mtx. 
/// </summary>
To reproduce the reference solution, run the Abaqus job tri3_test0.inp and look at
TRI3_TEST0_STIFFNESS_MATRICES.mtx.
To reproduce the reference solution run the Abaqus job tri3_test0 . inp and look at TRI3_TEST0_STIFFNESS_MATRICES . mtx .
[Fact]
        private static void TestStiffness0()
        {
            var factory = new ContinuumElement2DFactory(thickness, material0, dynamicMaterial);
            ContinuumElement2D tri3 = factory.CreateElement(CellType.Tri3, nodeSet0);
            IMatrix K = tri3.BuildStiffnessMatrix();
            double[,] expectedK = new double[,]
            {
                { 43303.16742081400, 5294.11764705880, -43778.28054298600, -44796.38009049800, 475.11312217194, 39502.26244343900 },
                { 5294.11764705880, 122361.99095023000, -39027.14932126700, -104660.63348416000, 33733.03167420800, -17701.35746606300 },
                { -43778.28054298600, -39027.14932126700, 151866.51583710000, 66176.47058823500, -108088.23529412000, -27149.32126696800 },
                { -44796.38009049800, -104660.63348416000, 66176.47058823500, 127601.80995475000, -21380.09049773800, -22941.17647058800 },
                { 475.11312217194, 33733.03167420800, -108088.23529412000, -21380.09049773800, 107613.12217195000, -12352.94117647100 },
                { 39502.26244343900, -17701.35746606300, -27149.32126696800, -22941.17647058800, -12352.94117647100, 40642.53393665200 }
            }; 
            Assert.True(K.Equals(Matrix.CreateFromArray(expectedK), 1e-10));
        }
[ Fact ] private static void TestStiffness0 ( ) { var factory = new ContinuumElement2DFactory ( thickness , material0 , dynamicMaterial ) ; ContinuumElement2D tri3 = factory . CreateElement ( CellType . Tri3 , nodeSet0 ) ; IMatrix K = tri3 . BuildStiffnessMatrix ( ) ; double [ , ] expectedK = new double [ , ] { { 43303.16742081400 , 5294.11764705880 , - 43778.28054298600 , - 44796.38009049800 , 475.11312217194 , 39502.26244343900 } , { 5294.11764705880 , 122361.99095023000 , - 39027.14932126700 , - 104660.63348416000 , 33733.03167420800 , - 17701.35746606300 } , { - 43778.28054298600 , - 39027.14932126700 , 151866.51583710000 , 66176.47058823500 , - 108088.23529412000 , - 27149.32126696800 } , { - 44796.38009049800 , - 104660.63348416000 , 66176.47058823500 , 127601.80995475000 , - 21380.09049773800 , - 22941.17647058800 } , { 475.11312217194 , 33733.03167420800 , - 108088.23529412000 , - 21380.09049773800 , 107613.12217195000 , - 12352.94117647100 } , { 39502.26244343900 , - 17701.35746606300 , - 27149.32126696800 , - 22941.17647058800 , - 12352.94117647100 , 40642.53393665200 } } ; Assert . True ( K . Equals ( Matrix . CreateFromArray ( expectedK ) , 1e-10 ) ) ; }
To reproduce the reference solution, run the Abaqus job tri3_test0.inp and look at
TRI3_TEST0_STIFFNESS_MATRICES.mtx.
To reproduce the reference solution run the Abaqus job tri3_test0 . inp and look at TRI3_TEST0_STIFFNESS_MATRICES . mtx .
// from Abaqus