update

Author: epsi1on

BriefFiniteElementNet.Validation/Case_03/Validator.cs

@@ -4,19 +4,20 @@
 using System.Linq;
 using System.Text;
 using System.Threading.Tasks;
+using BriefFiniteElementNet.Common;
 using BriefFiniteElementNet.Elements;
 
 namespace BriefFiniteElementNet.Validation.Case_03
 {
-    //[ValidationCase("Console beam with tetrahedron", typeof(TetrahedronElement))]
+    [ValidationCase("Console beam with tetrahedron", typeof(TetrahedronElement))]
     public class Validator : IValidationCase
     {
         public ValidationResult Validate()
         {
 
             /**/
             {
-                var model = StructureGenerator.Generate3DTetrahedralElementGrid(3, 3, 100);
+                var model = StructureGenerator.Generate3DTetrahedralElementGrid(4, 4, 100);
 
                 var e = 210e9;
 
@@ -40,7 +41,7 @@ public ValidationResult Validate()
                 var cnt = model.Nodes.Where(i => i.Location.Z == l);
 
                 var f = 1e7;
-                var I = dx * dy * dy * dy / 12;
+                var I = dy * dx * dx * dx / 12;
                 var rigid = new MpcElements.RigidElement_MPC() { UseForAllLoads = true };
 
                 foreach (var node in cnt)
@@ -49,8 +50,10 @@ public ValidationResult Validate()
                     rigid.Nodes.Add(node);
                 }
 
+                //model.MpcElements.Add(rigid);
+                model.Trace.Listeners.Add(new ConsoleTraceListener()); 
                 model.Solve_MPC();
-
+                
 
                 var delta = f * l * l * l / (3 * e * I);
 

BriefFiniteElementNet/ElementHelpers/IElementHelper.cs

@@ -19,7 +19,7 @@ public interface IElementHelper
         /// Gets the B matrix at defined isometric coordinates (B is derivation of N regarding to local x or y or z, not regarding to ξ, η or γ -- it is ∂N/∂x or ..., it is not ∂N/∂ξ or ...).
         /// </summary>
         /// <param name="targetElement">The target element.</param>
-        /// <param name="isoCoords">The isometric coordinations.</param>
+        /// <param name="isoCoords">The isometric coordination.</param>
         /// 
         /// <returns></returns>
         Matrix GetBMatrixAt(Element targetElement, params double[] isoCoords);

BriefFiniteElementNet/ElementHelpers/TetrahedronHelper.cs

@@ -104,7 +104,7 @@ public Matrix GetBMatrixAt(Element targetElement, params double[] isoCoords)
                 buf.SetRow(4, new[] { 0, d1, c1, 0, d2, c2, 0, d3, c3, 0, d4, c4 });
                 buf.SetRow(5, new[] { d1, 0, b1, d2, 0, b2, d3, 0, b3, d4, 0, b4 });
 
-                //buf.Scale(1 / (6.0 * v));
+                buf.Scale(1 / (6.0 * v));
             }
 
             return buf.AsMatrix();