public static void TestTriangles(int x = 7, int y = 5)
{
Model model2d = new Cantilever2D(ElementType.TriangleElement, x, y).Model;
FESystem sys = new FESystem(model2d);
sys.Solve();
FEPrint.PrintSystemInfo(sys);
FEPrint.PrintDisplacement(sys);
FEIO.WriteCOOMatrix(sys.GetKG().ToCOO(), "C:/Users/alber/Desktop/matA.mtx");
}
public static void TestHexahedron()
{
List <Node> nodes = new List <Node>()
{
new Node(0.0, 0.0, 0.0),
new Node(0.25, 0.0, 0.0),
new Node(0.25, 0.25, 0.0),
new Node(0.0, 0.25, 0.0),
new Node(0.0, 0.0, 0.25),
new Node(0.25, 0.0, 0.25),
new Node(0.25, 0.25, 0.25),
new Node(0.0, 0.25, 0.25),
new Node(0.5, 0.0, 0.0),
new Node(0.5, 0.0, 0.25),
new Node(0.5, 0.25, 0.25),
new Node(0.5, 0.25, 0.0)
};
var element0 = new Hexahedron(new List <Node>(8)
{
nodes[0], nodes[1], nodes[2], nodes[3], nodes[4], nodes[5], nodes[6], nodes[7]
}, new Material(1, 0.3));
//var element1 = new Hexahedron(new List<Node>(8) { nodes[4], nodes[5], nodes[6], nodes[7], nodes[8], nodes[9], nodes[10], nodes[11] }, new Material(1, 0.3));
var elements = new List <Element>(1)
{
element0
};
//var loads = new List<Load>(4)
//{
// new Load(8, 4.6875, 0.0,0.0),
// new Load(9, 4.6875, 0.0,0.0),
// new Load(10,4.6875, 0.0,0.0),
// new Load(11, 4.6875, 0.0,0.0)
//};
//var supports = new List<Support>(4)
//{
// new Support(0, SupportType.Fixed),
// new Support(1, SupportType.Fixed),
// new Support(2, SupportType.Fixed),
// new Support(3, SupportType.Fixed)
//};
var model = new Model(3, nodes, elements);
FESystem sys = new FESystem(model, Solver.SimplicialLLT);
sys.Initialize();
Console.WriteLine(model.Elements[0].Ke);
Console.WriteLine(model.Elements[0].B);
Console.WriteLine(model.Elements[0].D);
//sys.Solve();
Console.Write(sys.Model.ModelInfo());
Console.Write(sys.MatrixInfo());
Console.WriteLine("Solver: " + sys._Solver.ToString());
Console.WriteLine("Solving: " + sys.TimeCost[3].ToString() + " ms");
Console.WriteLine(sys.DisplacementInfo());
FEIO.WriteKG(sys.GetKG(), "E:\\ALCoding\\ALFE\\topoptTest", false);
Console.WriteLine();
}