public Model CreateModelFromFile()
{
var model = new Model();
model.SubdomainsDictionary[0] = new Subdomain(0);
// Material
char[] delimeters = { ' ', '=', '\t' };
Attributes?name = null;
String[] text = System.IO.File.ReadAllLines(Filename);
elementBoundaries = new List <IList <Element> >();
elementDomains = new List <IList <Element> >();
nodeDomains = new List <IList <Node> >();
nodeBoundaries = new List <IList <Node> >();
quadBoundaries = new List <IList <IList <Node> > >();
triBoundaries = new List <IList <IList <Node> > >();
for (int i = 0; i < 10; i++)
{
elementBoundaries.Add(new List <Element>());
nodeBoundaries.Add(new List <Node>());
quadBoundaries.Add(new List <IList <Node> >());
triBoundaries.Add(new List <IList <Node> >());
}
for (int i = 0; i < 2; i++)
{
elementDomains.Add(new List <Element>());
nodeDomains.Add(new List <Node>());
}
for (int i = 0; i < text.Length; i++)
{
String[] line = text[i].Split(delimeters, StringSplitOptions.RemoveEmptyEntries);
StringBuilder comparisonString = new StringBuilder(null, 50);
if (line.Length == 0)
{
continue;
}
if ((line[0] == "3" || line[0] == "4") && line.Length > 1)
{
try
{
name = (Attributes)Enum.Parse(typeof(Attributes), line[1]);
}
catch (Exception exception)
{
name = null;
}
}
else
{
for (int linePosition = 0; linePosition < line.Length; linePosition++)
{
if (line[linePosition] == "#")
{
for (int ij = linePosition + 1; ij < line.Length; ij++)
{
comparisonString.Append(line[ij]);
}
try
{
name = (Attributes)Enum.Parse(typeof(Attributes), comparisonString.ToString());
}
catch (Exception exception)
{
name = null;
}
}
}
}
switch (name)
{
case Attributes.sdim:
int NumberOfDimensions = Int32.Parse(line[0]);
break;
case Attributes.numberofmeshpoints:
NumberOfNodes = Int32.Parse(line[0]);
break;
case Attributes.Meshpointcoordinates:
IList <Node> nodelist = new List <Node> {
null
};
for (int j = 0; j < NumberOfNodes; j++)
{
i++;
line = text[i].Split(delimeters);
int nodeGlobalID = j;
double x = Double.Parse(line[0], CultureInfo.InvariantCulture);
double y = Double.Parse(line[1], CultureInfo.InvariantCulture);
double z = Double.Parse(line[2], CultureInfo.InvariantCulture);
Node node = new Node(nodeGlobalID, x, y, z);
model.NodesDictionary.Add(nodeGlobalID, node);
nodelist.Add(node);
}
break;
case Attributes.tri:
do
{
i++;
line = text[i].Split(delimeters);
} while (line[0] == "");
i++;
line = text[i].Split(delimeters);
NumberOfTriElements = Int32.Parse(line[0]);
i++;
IList <IList <Node> > triNodesCollection = new List <IList <Node> >();
for (int TriID = 0; TriID < NumberOfTriElements; TriID++)
{
i++;
line = text[i].Split(delimeters);
triNodesCollection.Add(new List <Node>
{
model.NodesDictionary[Int32.Parse(line[2])],
model.NodesDictionary[Int32.Parse(line[1])],
model.NodesDictionary[Int32.Parse(line[0])]
});
}
i = i + 3;
for (int TriID = 0; TriID < NumberOfTriElements; TriID++)
{
i++;
line = text[i].Split(delimeters);
int boundaryID = Int32.Parse(line[0]);
triBoundaries[boundaryID].Add(triNodesCollection[TriID]);
foreach (Node node in triNodesCollection[TriID])
{
nodeBoundaries[boundaryID].Add(node);
}
//int elementBoundaryID = Int32.Parse(line[0]);
//elementBoundaries[elementBoundaryID].Add(model.ElementsDictionary[QuadID]);
}
for (int boundaryID = 0; boundaryID < nodeBoundaries.Count; boundaryID++)
{
nodeBoundaries[boundaryID] = nodeBoundaries[boundaryID].Distinct().ToList();
}
break;
case Attributes.quad:
do
{
i++;
line = text[i].Split(delimeters);
} while (line[0] == "");
i++;
line = text[i].Split(delimeters);
NumberOfQuadElements = Int32.Parse(line[0]);
i++;
IList <IList <Node> > quadNodesCollection = new List <IList <Node> >();
for (int QuadID = 0; QuadID < NumberOfQuadElements; QuadID++)
{
i++;
line = text[i].Split(delimeters);
//IReadOnlyList<Node> nodes = new List<Node>
//{
// model.NodesDictionary[Int32.Parse(line[1])], //0
// model.NodesDictionary[Int32.Parse(line[3])], //1
// model.NodesDictionary[Int32.Parse(line[2])], //2
// model.NodesDictionary[Int32.Parse(line[0])] //3
//};
quadNodesCollection.Add(new List <Node>
{
model.NodesDictionary[Int32.Parse(line[1])], //0
model.NodesDictionary[Int32.Parse(line[3])], //1
model.NodesDictionary[Int32.Parse(line[2])], //2
model.NodesDictionary[Int32.Parse(line[0])] //3
});
//var Quad4 = boundaryFactory3D.CreateElement(CellType.Quad4, nodes);
//var element = new Element();
//element.ID = QuadID;
//element.ElementType = Quad4;
//model.SubdomainsDictionary[0].Elements.Add(element);
//model.ElementsDictionary.Add(QuadID, element);
//foreach (Node node in nodes)
//{
// element.AddNode(node);
//}
}
i = i + 3;
for (int QuadID = 0; QuadID < NumberOfQuadElements; QuadID++)
{
i++;
line = text[i].Split(delimeters);
int boundaryID = Int32.Parse(line[0]);
quadBoundaries[boundaryID].Add(quadNodesCollection[QuadID]);
foreach (Node node in quadNodesCollection[QuadID])
{
nodeBoundaries[boundaryID].Add(node);
}
//int elementBoundaryID = Int32.Parse(line[0]);
//elementBoundaries[elementBoundaryID].Add(model.ElementsDictionary[QuadID]);
}
for (int boundaryID = 0; boundaryID < nodeBoundaries.Count; boundaryID++)
{
nodeBoundaries[boundaryID] = nodeBoundaries[boundaryID].Distinct().ToList();
}
break;
case Attributes.tet:
do
{
i++;
line = text[i].Split(delimeters);
} while (line[0] == "");
i++;
line = text[i].Split(delimeters);
NumberOfTetElements = Int32.Parse(line[0]);
i++;
List <IReadOnlyList <Node> > TetraNodes = new List <IReadOnlyList <Node> >();
for (int TetID = 0; TetID < NumberOfTetElements; TetID++)
{
i++;
line = text[i].Split(delimeters);
TetraNodes.Add(new List <Node>
{
model.NodesDictionary[Int32.Parse(line[3])],
model.NodesDictionary[Int32.Parse(line[2])],
model.NodesDictionary[Int32.Parse(line[1])],
model.NodesDictionary[Int32.Parse(line[0])]
});
}
i = i + 3;
for (int TetID = 0; TetID < NumberOfTetElements; TetID++)
{
i++;
line = text[i].Split(delimeters);
int elementDomainID = Int32.Parse(line[0]);
HyperElasticMaterial3DDefGrad hyperElasticMaterial = new HyperElasticMaterial3DDefGrad()
{
C1 = C1[elementDomainID - 1], C2 = C2[elementDomainID - 1], k_cons = k_cons[elementDomainID - 1]
};
ContinuumElement3DNonLinearDefGradFactory elementFactory;
if (lambdag == null)
{
elementFactory = new ContinuumElement3DNonLinearDefGradFactory(hyperElasticMaterial,
CommonDynamicProperties[elementDomainID - 1]);
}
else
{
elementFactory = new ContinuumElement3DNonLinearDefGradFactory(hyperElasticMaterial,
CommonDynamicProperties[elementDomainID - 1], lambdag[TetID]);
}
IReadOnlyList <Node> nodesTet = TetraNodes[TetID];
var Tet4 = elementFactory.CreateElement(CellType.Tet4, nodesTet);
var element = new Element();
element.ID = TetID;
element.ElementType = Tet4;
foreach (Node node in nodesTet)
{
element.AddNode(node);
nodeDomains[elementDomainID - 1].Add(node);
}
model.SubdomainsDictionary[0].Elements.Add(element);
model.ElementsDictionary.Add(TetID, element);
elementDomains[elementDomainID - 1].Add(model.ElementsDictionary[TetID]);
//int elementBoundaryID = Int32.Parse(line[0]);
//elementBoundaries[elementBoundaryID].Add(model.ElementsDictionary[QuadID]);
}
for (int elementDomainID = 0; elementDomainID < elementDomains.Count; elementDomainID++)
{
nodeDomains[elementDomainID] = nodeDomains[elementDomainID].Distinct().ToList();
}
break;
case Attributes.hex:
do
{
i++;
line = text[i].Split(delimeters);
} while (line[0] == "");
i++;
line = text[i].Split(delimeters);
NumberOfHexElements = Int32.Parse(line[0]);
i++;
List <IReadOnlyList <Node> > HexaNodes = new List <IReadOnlyList <Node> >();
for (int HexID = 0; HexID < NumberOfHexElements; HexID++)
{
i++;
line = text[i].Split(delimeters);
HexaNodes.Add(new List <Node>
{
model.NodesDictionary[Int32.Parse(line[4])],
model.NodesDictionary[Int32.Parse(line[6])],
model.NodesDictionary[Int32.Parse(line[7])],
model.NodesDictionary[Int32.Parse(line[5])],
model.NodesDictionary[Int32.Parse(line[0])],
model.NodesDictionary[Int32.Parse(line[2])],
model.NodesDictionary[Int32.Parse(line[3])],
model.NodesDictionary[Int32.Parse(line[1])],
});
}
i = i + 3;
for (int HexID = 0; HexID < NumberOfHexElements; HexID++)
{
i++;
line = text[i].Split(delimeters);
int elementDomainID = Int32.Parse(line[0]);
HyperElasticMaterial3DDefGrad hyperElasticMaterial = new HyperElasticMaterial3DDefGrad()
{
C1 = C1[elementDomainID - 1],
C2 = C2[elementDomainID - 1],
k_cons = k_cons[elementDomainID - 1]
};
ContinuumElement3DNonLinearDefGradFactory elementFactory;
if (lambdag == null)
{
elementFactory = new ContinuumElement3DNonLinearDefGradFactory(hyperElasticMaterial,
CommonDynamicProperties[elementDomainID - 1]);
}
else
{
elementFactory = new ContinuumElement3DNonLinearDefGradFactory(hyperElasticMaterial,
CommonDynamicProperties[elementDomainID - 1], lambdag[HexID]);
}
IReadOnlyList <Node> nodesHex = HexaNodes[HexID];
var Hexa8 = elementFactory.CreateElement(CellType.Hexa8, nodesHex);
var element = new Element();
element.ID = HexID;
element.ElementType = Hexa8;
foreach (Node node in nodesHex)
{
element.AddNode(node);
nodeDomains[elementDomainID - 1].Add(node);
}
model.SubdomainsDictionary[0].Elements.Add(element);
model.ElementsDictionary.Add(HexID, element);
elementDomains[elementDomainID - 1].Add(model.ElementsDictionary[HexID]);
}
for (int elementDomainID = 0; elementDomainID < elementDomains.Count; elementDomainID++)
{
nodeDomains[elementDomainID] = nodeDomains[elementDomainID].Distinct().ToList();
}
break;
}
}
return(model);
}
internal static int DefineContinuumElement3DNonLinearDefGrad(Model model, IContinuumMaterial3DDefGrad material3D, IDynamicMaterial dynamicMaterial)
{
// Node creation
var node1 = new Node(id: 1, x: 0.0, y: 0.0, z: 0.0);
var node2 = new Node(id: 2, x: 100.0, y: 0.0, z: 0.0);
var node3 = new Node(id: 3, x: 0.0, y: 100.0, z: 0.0);
var node4 = new Node(id: 4, x: 100.0, y: 100.0, z: 0.0);
var node5 = new Node(id: 5, x: 0.0, y: 0.0, z: 100.0);
var node6 = new Node(id: 6, x: 100.0, y: 0.0, z: 100.0);
var node7 = new Node(id: 7, x: 0.0, y: 100.0, z: 100.0);
var node8 = new Node(id: 8, x: 100.0, y: 100.0, z: 100.0);
// Create List of nodes
IList <Node> nodes = new List <Node>();
nodes.Add(node1);
nodes.Add(node2);
nodes.Add(node3);
nodes.Add(node4);
nodes.Add(node5);
nodes.Add(node6);
nodes.Add(node7);
nodes.Add(node8);
numberOfNodes = nodes.Count;
// Add nodes to the nodes dictonary of the model
for (int i = 0; i < numberOfNodes; ++i)
{
model.NodesDictionary.Add(i, nodes[i]);
}
int[][] connectivityMatrix = new int[8][];
connectivityMatrix[0] = new int[] { 1, 1, 2, 4, 3, 5, 6, 8, 7 };
var factory = new ContinuumElement3DNonLinearDefGradFactory(material3D, dynamicMaterial, lambdag);
for (int i = 0; i < 1; i++)
{
List <Node> elementNodeSet = new List <Node>(8);
for (int j = 1; j < 9; j++)
{
elementNodeSet.Add((Node)model.NodesDictionary[connectivityMatrix[i][j] - 1]);
}
var hexa8element = new Element()
{
ID = connectivityMatrix[i][0],
ElementType = factory.CreateElement(CellType.Hexa8, elementNodeSet),
};
for (int j = 1; j < 9; j++)
{
hexa8element.AddNode(model.NodesDictionary[connectivityMatrix[i][j] - 1]);
}
model.ElementsDictionary.Add(hexa8element.ID, hexa8element);
model.SubdomainsDictionary[0].Elements.Add(hexa8element);
}
return(numberOfNodes);
}
