private static Model CreateModel()
{
IList <IList <Node> > nodeBoundaries;
IList <IList <Element> > elementBoundaries;
double density = 1.0;
double c = 1.0;
double k = 1.0;
double[] U = { 1.0, 1, 1 };
double h = .0;
var elementFactory3D = new ConvectionDiffusionElement3DFactory(new ConvectionDiffusionMaterial(k, U, 0));
var boundaryFactory3D = new SurfaceBoundaryFactory3D(0, new ConvectionDiffusionMaterial(k, U, 0));
var model = new Model();
model.SubdomainsDictionary[0] = new Subdomain(0);
elementBoundaries = new List <IList <Element> >();
nodeBoundaries = new List <IList <Node> >();
for (int i = 0; i < 10; i++)
{
elementBoundaries.Add(new List <Element>());
nodeBoundaries.Add(new List <Node>());
}
double[,] nodeData = new double[, ]
{
{ 0, 0, 0 }, { 0, .1, 0 }, { 0, .1, .1 }, { 0, 0, .1 }, { .1, 0, 0 }, { .1, .1, 0 }, { .1, .1, .1 }, { .1, 0, .1 },
{ 0, 0.05, 0 }, { 0, 0, 0.05 }, { .05, 0, 0 }, { 0, .1, 0.05 }, { .05, .1, 0 }, { 0, 0.05, .1 }, { .05, .1, .1 },
{ .05, 0, .1 }, { .1, 0.05, 0 }, { .1, 0, 0.05 }, { .1, .1, 0.05 }, { .1, 0.05, .1 }, { 0, 0.05, 0.05 },
{ .05, 0.05, 0 }, { .05, 0, 0.05 }, { .05, .1, 0.05 }, { .05, 0.05, .1 }, { .1, 0.05, 0.05 },
{ .05, 0.05, 0.05 },
};
for (int nNode = 0; nNode < nodeData.GetLength(0); nNode++)
{
model.NodesDictionary.Add(nNode, new Node(id: nNode, x: nodeData[nNode, 0], y: nodeData[nNode, 1], z: nodeData[nNode, 2]));
}
IReadOnlyList <Node> nodes = new List <Node>
{
model.NodesDictionary[0], model.NodesDictionary[1], model.NodesDictionary[2],
model.NodesDictionary[3], model.NodesDictionary[4], model.NodesDictionary[5],
model.NodesDictionary[6], model.NodesDictionary[7], model.NodesDictionary[8],
model.NodesDictionary[9], model.NodesDictionary[10], model.NodesDictionary[11],
model.NodesDictionary[12], model.NodesDictionary[13], model.NodesDictionary[14],
model.NodesDictionary[15], model.NodesDictionary[16], model.NodesDictionary[17],
model.NodesDictionary[18], model.NodesDictionary[19], model.NodesDictionary[20],
model.NodesDictionary[21], model.NodesDictionary[22], model.NodesDictionary[23],
model.NodesDictionary[24], model.NodesDictionary[25], model.NodesDictionary[26],
};
IList <IReadOnlyList <Node> > face = new List <IReadOnlyList <Node> >();
//face = new List<IReadOnlyList<Node>>();
for (int i = 0; i < 6; i++)
{
face.Add(new List <Node>());
}
face[0] = new List <Node>
{
model.NodesDictionary[0],
model.NodesDictionary[1],
model.NodesDictionary[2],
model.NodesDictionary[3],
model.NodesDictionary[4],
model.NodesDictionary[5],
model.NodesDictionary[6],
model.NodesDictionary[7],
model.NodesDictionary[8],
};
face[1] = new List <Node>
{
model.NodesDictionary[0],
model.NodesDictionary[1],
model.NodesDictionary[2],
model.NodesDictionary[9],
model.NodesDictionary[10],
model.NodesDictionary[11],
model.NodesDictionary[18],
model.NodesDictionary[19],
model.NodesDictionary[20],
};
face[2] = new List <Node>
{
model.NodesDictionary[18],
model.NodesDictionary[19],
model.NodesDictionary[20],
model.NodesDictionary[21],
model.NodesDictionary[22],
model.NodesDictionary[23],
model.NodesDictionary[24],
model.NodesDictionary[25],
model.NodesDictionary[26],
};
face[3] = new List <Node>
{
model.NodesDictionary[6],
model.NodesDictionary[7],
model.NodesDictionary[8],
model.NodesDictionary[15],
model.NodesDictionary[16],
model.NodesDictionary[17],
model.NodesDictionary[24],
model.NodesDictionary[25],
model.NodesDictionary[26],
};
face[4] = new List <Node>
{
model.NodesDictionary[0],
model.NodesDictionary[3],
model.NodesDictionary[6],
model.NodesDictionary[9],
model.NodesDictionary[12],
model.NodesDictionary[15],
model.NodesDictionary[18],
model.NodesDictionary[21],
model.NodesDictionary[24],
};
face[5] = new List <Node>
{
model.NodesDictionary[2],
model.NodesDictionary[5],
model.NodesDictionary[8],
model.NodesDictionary[11],
model.NodesDictionary[14],
model.NodesDictionary[17],
model.NodesDictionary[20],
model.NodesDictionary[23],
model.NodesDictionary[26],
};
var Hexa27 = elementFactory3D.CreateElement(CellType.Hexa27, nodes);
var element = new Element();
element.ID = 0;
element.ElementType = Hexa27;
foreach (Node node in nodes)
{
element.AddNode(node);
}
model.SubdomainsDictionary[0].Elements.Add(element);
model.ElementsDictionary.Add(0, element);
for (int i = 1; i < 7; i++)
{
var Quad = boundaryFactory3D.CreateElement(CellType.Quad9, face[i - 1]);
var faceElement = new Element();
faceElement.ID = i;
faceElement.ElementType = Quad;
foreach (Node node in face[i - 1])
{
faceElement.AddNode(node);
}
model.SubdomainsDictionary[0].Elements.Add(faceElement);
model.ElementsDictionary.Add(i, faceElement);
}
//constraints
foreach (int i in new int[] { 2, 6, 7 })
{
model.NodesDictionary[i].Constraints.Add(new Constraint()
{
Amount = 0,
DOF = ThermalDof.Temperature
});
}
// fortish korufhs
Load load1;
foreach (int i in new int[] { 0, 1, 3, 4, 5 })
{
load1 = new Load()
{
Node = model.NodesDictionary[i],
DOF = ThermalDof.Temperature,
Amount = 100
};
model.Loads.Add(load1);
}
return(model);
}
private static Model CreateModel()
{
var model = new Model();
// Subdomains
model.SubdomainsDictionary.Add(0, new Subdomain(subdomainID));
// Material
double density = 1.0;
double k = 1.0;
double c = 1.0;
double h = 0;
// Nodes
int numNodes = 27;
var nodes = new Node[numNodes];
nodes[0] = new Node(id: 0, x: 2.0, y: 2.0, z: 2.0);
nodes[1] = new Node(id: 1, x: 2.0, y: 1.0, z: 2.0);
nodes[2] = new Node(id: 2, x: 2.0, y: 0.0, z: 2.0);
nodes[3] = new Node(id: 3, x: 2.0, y: 2.0, z: 1.0);
nodes[4] = new Node(id: 4, x: 2.0, y: 1.0, z: 1.0);
nodes[5] = new Node(id: 5, x: 2.0, y: 0.0, z: 1.0);
nodes[6] = new Node(id: 6, x: 2.0, y: 2.0, z: 0.0);
nodes[7] = new Node(id: 7, x: 2.0, y: 1.0, z: 0.0);
nodes[8] = new Node(id: 8, x: 2.0, y: 0.0, z: 0.0);
nodes[9] = new Node(id: 9, x: 1.0, y: 2.0, z: 2.0);
nodes[10] = new Node(id: 10, x: 1.0, y: 1.0, z: 2.0);
nodes[11] = new Node(id: 11, x: 1.0, y: 0.0, z: 2.0);
nodes[12] = new Node(id: 12, x: 1.0, y: 2.0, z: 1.0);
nodes[13] = new Node(id: 13, x: 1.0, y: 1.0, z: 1.0);
nodes[14] = new Node(id: 14, x: 1.0, y: 0.0, z: 1.0);
nodes[15] = new Node(id: 15, x: 1.0, y: 2.0, z: 0.0);
nodes[16] = new Node(id: 16, x: 1.0, y: 1.0, z: 0.0);
nodes[17] = new Node(id: 17, x: 1.0, y: 0.0, z: 0.0);
nodes[18] = new Node(id: 18, x: 0.0, y: 2.0, z: 2.0);
nodes[19] = new Node(id: 19, x: 0.0, y: 1.0, z: 2.0);
nodes[20] = new Node(id: 20, x: 0.0, y: 0.0, z: 2.0);
nodes[21] = new Node(id: 21, x: 0.0, y: 2.0, z: 1.0);
nodes[22] = new Node(id: 22, x: 0.0, y: 1.0, z: 1.0);
nodes[23] = new Node(id: 23, x: 0.0, y: 0.0, z: 1.0);
nodes[24] = new Node(id: 24, x: 0.0, y: 2.0, z: 0.0);
nodes[25] = new Node(id: 25, x: 0.0, y: 1.0, z: 0.0);
nodes[26] = new Node(id: 26, x: 0.0, y: 0.0, z: 0.0);
for (int i = 0; i < numNodes; ++i)
{
model.NodesDictionary[i] = nodes[i];
}
// Elements
int numElements = 8;
var elementFactory = new ConvectionDiffusionElement3DFactory(new ConvectionDiffusionMaterial(k, new double[] { 0, 0, 0 }, 0));
var elements = new ConvectionDiffusionElement3D[8];
elements[0] = elementFactory.CreateElement(CellType.Hexa8, new Node[] { nodes[13], nodes[4], nodes[3], nodes[12], nodes[10], nodes[1], nodes[0], nodes[9] });
elements[1] = elementFactory.CreateElement(CellType.Hexa8, new Node[] { nodes[14], nodes[5], nodes[4], nodes[13], nodes[11], nodes[2], nodes[1], nodes[10] });
elements[2] = elementFactory.CreateElement(CellType.Hexa8, new Node[] { nodes[16], nodes[7], nodes[6], nodes[15], nodes[13], nodes[4], nodes[3], nodes[12] });
elements[3] = elementFactory.CreateElement(CellType.Hexa8, new Node[] { nodes[17], nodes[8], nodes[7], nodes[16], nodes[14], nodes[5], nodes[4], nodes[13] });
elements[4] = elementFactory.CreateElement(CellType.Hexa8, new Node[] { nodes[22], nodes[13], nodes[12], nodes[21], nodes[19], nodes[10], nodes[9], nodes[18] });
elements[5] = elementFactory.CreateElement(CellType.Hexa8, new Node[] { nodes[23], nodes[14], nodes[13], nodes[22], nodes[20], nodes[11], nodes[10], nodes[19] });
elements[6] = elementFactory.CreateElement(CellType.Hexa8, new Node[] { nodes[25], nodes[16], nodes[15], nodes[24], nodes[22], nodes[13], nodes[12], nodes[21] });
elements[7] = elementFactory.CreateElement(CellType.Hexa8, new Node[] { nodes[26], nodes[17], nodes[16], nodes[25], nodes[23], nodes[14], nodes[13], nodes[22] });
for (int i = 0; i < numElements; ++i)
{
var elementWrapper = new Element()
{
ID = i, ElementType = elements[i]
};
foreach (var node in elements[i].Nodes)
{
elementWrapper.AddNode(node);
}
model.ElementsDictionary[i] = elementWrapper;
model.SubdomainsDictionary[subdomainID].Elements.Add(elementWrapper);
}
// Dirichlet BC
model.NodesDictionary[0].Constraints.Add(new Constraint()
{
DOF = ThermalDof.Temperature, Amount = 100.0
});
model.NodesDictionary[1].Constraints.Add(new Constraint()
{
DOF = ThermalDof.Temperature, Amount = 100.0
});
model.NodesDictionary[2].Constraints.Add(new Constraint()
{
DOF = ThermalDof.Temperature, Amount = 100.0
});
model.NodesDictionary[9].Constraints.Add(new Constraint()
{
DOF = ThermalDof.Temperature, Amount = 100.0
});
model.NodesDictionary[10].Constraints.Add(new Constraint()
{
DOF = ThermalDof.Temperature, Amount = 100.0
});
model.NodesDictionary[11].Constraints.Add(new Constraint()
{
DOF = ThermalDof.Temperature, Amount = 100.0
});
model.NodesDictionary[18].Constraints.Add(new Constraint()
{
DOF = ThermalDof.Temperature, Amount = 100.0
});
model.NodesDictionary[19].Constraints.Add(new Constraint()
{
DOF = ThermalDof.Temperature, Amount = 100.0
});
model.NodesDictionary[20].Constraints.Add(new Constraint()
{
DOF = ThermalDof.Temperature, Amount = 100.0
});
// Neumann BC
double q = 100;
model.Loads.Add(new Load()
{
Amount = q, Node = model.NodesDictionary[6], DOF = ThermalDof.Temperature
});
model.Loads.Add(new Load()
{
Amount = q, Node = model.NodesDictionary[24], DOF = ThermalDof.Temperature
});
return(model);
}
public Model UpdateModel(Model structuralModel = null, Dictionary <int, IVector> displacements = null, bool updateNodes = true)
{//only works for tet4 now
this.displacements = displacements;
if (updateNodes && displacements != null)
{
for (int id = 0; id < displacements.Count; id++)
{
Model.NodesDictionary.Clear();
int count = 0;
foreach (Node n in structuralModel.Nodes)
{
double x = n.X;
double y = n.Y;
double z = n.Z;
if (structuralModel.GlobalDofOrdering.GlobalFreeDofs.Contains(n, StructuralDof.TranslationX))
{
x += displacements[id][count];
count += 1;
}
if (structuralModel.GlobalDofOrdering.GlobalFreeDofs.Contains(n, StructuralDof.TranslationY))
{
y += displacements[id][count];
count += 1;
}
if (structuralModel.GlobalDofOrdering.GlobalFreeDofs.Contains(n, StructuralDof.TranslationZ))
{
z += displacements[id][count];
count += 1;
}
Node node = new Node(n.ID, x, y, z);
Model.NodesDictionary.Add(n.ID, node);
}
}
}
ConvectionDiffusionMaterial[] CDMaterial = new ConvectionDiffusionMaterial[diffusionCoeff.Length];
ConvectionDiffusionElement3DFactory[] elementFactory3D = new ConvectionDiffusionElement3DFactory[diffusionCoeff.Length];
for (int i = 0; i < diffusionCoeff.Length; i++)
{
CDMaterial[i] = new ConvectionDiffusionMaterial(capacityCoeff[i], diffusionCoeff[i], convectionCoeff[i], loadFromUnknownCoeff[i]);
elementFactory3D[i] = new ConvectionDiffusionElement3DFactory(CDMaterial[i]);
}
for (int domain = 0; domain < elementDomains.Count; domain++)
{
foreach (Element elem in elementDomains[domain])
{
var oldElemNodes = elem.Nodes.ToList();
var newElemNodes = new List <Node>();
foreach (Node n in oldElemNodes)
{
newElemNodes.Add(Model.Nodes[n.ID]);
}
elem.NodesDictionary.Clear();
elem.ElementType = elementFactory3D[domain].CreateElement(CellType.Tet4, newElemNodes);
elem.AddNodes(newElemNodes);
}
}
//ConvectionDiffusionMaterial[] CDMaterial = new ConvectionDiffusionMaterial[diffusionCoeff.Length];
//ConvectionDiffusionElement3DFactory[] elementFactory3D = new ConvectionDiffusionElement3DFactory[diffusionCoeff.Length];
//for (int i = 0; i < diffusionCoeff.Length; i++)
//{
// CDMaterial[i] = new ConvectionDiffusionMaterial(diffusionCoeff[i], convectionCoeff[i], loadFromUnknownCoeff[i]);
// elementFactory3D[i] = new ConvectionDiffusionElement3DFactory(CDMaterial[i]);
//}
//for (int domain = 0; domain < elementDomains.Count; domain++)
// foreach (Element elem in elementDomains[domain])
// elem.ElementType = elementFactory3D[domain].CreateElement(CellType.Tet4, elem.Nodes.ToList());
Model.BodyLoads.Clear();
return(Model);
}
public Model CreateModelFromFile(bool AllDomains = true)
{
ConvectionDiffusionMaterial[] CDMaterial = new ConvectionDiffusionMaterial[diffusionCoeff.Length];
ConvectionDiffusionElement3DFactory[] elementFactory3D = new ConvectionDiffusionElement3DFactory[diffusionCoeff.Length];
for (int i = 0; i < diffusionCoeff.Length; i++)
{
CDMaterial[i] = new ConvectionDiffusionMaterial(capacityCoeff[i], diffusionCoeff[i], convectionCoeff[i], loadFromUnknownCoeff[i]);
elementFactory3D[i] = new ConvectionDiffusionElement3DFactory(CDMaterial[i]);
}
//var boundaryFactory3D = new SurfaceBoundaryFactory3D(0, new ConvectionDiffusionMaterial(diffusionCoeff, new double[] {0,0,0}, 0));
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]);
IReadOnlyList <Node> nodesTet = TetraNodes[TetID];
var Tet4 = elementFactory3D[elementDomainID - 1].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]);
IReadOnlyList <Node> nodesHex = HexaNodes[HexID];
var Hexa8 = elementFactory3D[elementDomainID - 1].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;
}
}
Model = model;
return(model);
}
private static Model CreateModel()
{
IList <IList <Node> > nodeBoundaries;
IList <IList <Element> > elementBoundaries;
double density = 1.0;
double[] U = { 2, 2, 2 };
double k = 1.0;
double L = .0;
var elementFactory3D = new ConvectionDiffusionElement3DFactory(new ConvectionDiffusionMaterial(k, U, L));
var boundaryFactory3D = new SurfaceBoundaryFactory3D(0, new ConvectionDiffusionMaterial(k, U, L));
var model = new Model();
model.SubdomainsDictionary[0] = new Subdomain(0);
elementBoundaries = new List <IList <Element> >();
nodeBoundaries = new List <IList <Node> >();
for (int i = 0; i < 10; i++)
{
elementBoundaries.Add(new List <Element>());
nodeBoundaries.Add(new List <Node>());
}
double[,] nodeData = new double[, ] {
{ 0, 0, 0 },
{ 0, 1, 0 },
{ 0, 1, 1 },
{ 0, 0, 1 },
{ 1, 0, 0 },
{ 1, 1, 0 },
{ 1, 1, 1 },
{ 1, 0, 1 }
};
for (int nNode = 0; nNode < nodeData.GetLength(0); nNode++)
{
model.NodesDictionary.Add(nNode, new Node(id: nNode, x: nodeData[nNode, 0], y: nodeData[nNode, 1], z: nodeData[nNode, 2]));
}
IReadOnlyList <Node> nodes = new List <Node>
{
model.NodesDictionary[0],
model.NodesDictionary[1],
model.NodesDictionary[2],
model.NodesDictionary[3],
model.NodesDictionary[4],
model.NodesDictionary[5],
model.NodesDictionary[6],
model.NodesDictionary[7],
};
IList <IReadOnlyList <Node> > face = new List <IReadOnlyList <Node> >();
//face = new List<IReadOnlyList<Node>>();
for (int i = 0; i < 6; i++)
{
face.Add(new List <Node>());
}
face[0] = new List <Node>
{
model.NodesDictionary[0],
model.NodesDictionary[1],
model.NodesDictionary[2],
model.NodesDictionary[3],
};
face[1] = new List <Node>
{
model.NodesDictionary[4],
model.NodesDictionary[5],
model.NodesDictionary[6],
model.NodesDictionary[7],
};
face[2] = new List <Node>
{
model.NodesDictionary[0],
model.NodesDictionary[1],
model.NodesDictionary[5],
model.NodesDictionary[4],
};
face[3] = new List <Node>
{
model.NodesDictionary[1],
model.NodesDictionary[2],
model.NodesDictionary[6],
model.NodesDictionary[5],
};
face[4] = new List <Node>
{
model.NodesDictionary[2],
model.NodesDictionary[6],
model.NodesDictionary[7],
model.NodesDictionary[3],
};
face[5] = new List <Node>
{
model.NodesDictionary[0],
model.NodesDictionary[4],
model.NodesDictionary[7],
model.NodesDictionary[3],
};
int[] elementData = new int[] { 0, 1, 2, 3, 4, 5, 6, 7 };// the last line will not be used. We assign only one element
var Hexa8 = elementFactory3D.CreateElement(CellType.Hexa8, nodes);
var element = new Element();
element.ID = 0;
element.ElementType = Hexa8;
for (int j = 0; j < 8; j++)
{
element.NodesDictionary.Add(elementData[j], model.NodesDictionary[elementData[j]]);
}
model.SubdomainsDictionary[0].Elements.Add(element);
model.ElementsDictionary.Add(0, element);
var flux1 = new FluxLoad(1);
var flux2 = new FluxLoad(10);
var dir1 = new DirichletDistribution(list => {
return(Vector.CreateWithValue(list.Count, 0));
});
var dir2 = new DirichletDistribution(list => {
return(Vector.CreateWithValue(list.Count, 10));
});
var weakDirichlet1 = new WeakDirichlet(dir1, k);
var weakDirichlet2 = new WeakDirichlet(dir2, k);
var dirichletFactory1 = new SurfaceLoadElementFactory(weakDirichlet1);
var dirichletFactory2 = new SurfaceLoadElementFactory(weakDirichlet2);
var fluxFactory1 = new SurfaceLoadElementFactory(flux1);
var fluxFactory2 = new SurfaceLoadElementFactory(flux2);
//foreach (int i in new int[] { 1 })
//{
// var Quad = boundaryFactory3D.CreateElement(CellType.Quad4, face[i - 1]);
// var faceElement = new Element();
// faceElement.ID = i;
// faceElement.ElementType = Quad;
// foreach (Node node in face[i - 1])
// {
// faceElement.AddNode(node);
// }
// model.SubdomainsDictionary[0].Elements.Add(faceElement);
// model.ElementsDictionary.Add(i, faceElement);
// var dirichletElement = dirichletFactory1.CreateElement(CellType.Quad4, face[i - 1]);
// //var fluxElement = fluxFactory.CreateElement(CellType.Quad4, face[i - 1]);
// model.SurfaceLoads.Add(dirichletElement);
// //model.SurfaceLoads.Add(fluxElement);
//}
foreach (int i in new int[] { 0 })
{
//var Quad = boundaryFactory3D.CreateElement(CellType.Quad4, face[i]);
//var faceElement = new Element();
//faceElement.ID = 1;
//faceElement.ElementType = Quad;
//foreach (Node node in face[i])
//{
// faceElement.AddNode(node);
//}
//model.SubdomainsDictionary[0].Elements.Add(faceElement);
//model.ElementsDictionary.Add(1, faceElement);
//var dirichletElement = dirichletFactory2.CreateElement(CellType.Quad4, face[i]);
//model.SurfaceLoads.Add(dirichletElement);
var fluxElement = fluxFactory1.CreateElement(CellType.Quad4, face[i]);
model.SurfaceLoads.Add(fluxElement);
}
foreach (int i in new int[] { 1 })
{
//var Quad = boundaryFactory3D.CreateElement(CellType.Quad4, face[i]);
//var faceElement = new Element();
//faceElement.ID = 2;
//faceElement.ElementType = Quad;
//foreach (Node node in face[i])
//{
// faceElement.AddNode(node);
//}
//model.SubdomainsDictionary[0].Elements.Add(faceElement);
//model.ElementsDictionary.Add(2, faceElement);
//var dirichletElement = dirichletFactory2.CreateElement(CellType.Quad4, face[i]);
//model.SurfaceLoads.Add(dirichletElement);
var fluxElement = fluxFactory2.CreateElement(CellType.Quad4, face[i]);
model.SurfaceLoads.Add(fluxElement);
}
//constraints
//foreach (int i in new int[] { 4, 5, 6, 7 })
//{
// model.NodesDictionary[i].Constraints.Add(new Constraint()
// {
// Amount = 0,
// DOF = ThermalDof.Temperature
// });
//}
//// loads
//Load load1;
//foreach (int i in new int[] {0, 1, 2, 3})
//{
// load1 = new Load()
// {
// Node = model.NodesDictionary[i],
// DOF = ThermalDof.Temperature,
// Amount = 12.5
// };
// model.Loads.Add(load1);
//}
return(model);
}
public Model CreateModelFromFile()
{
double c = 1.0;
double[] U = { 1.0, 1, 1 };
double k = 1.0;
double L = .0;
diffusionCeoff = k;
var elementFactory3D = new ConvectionDiffusionElement3DFactory(new ConvectionDiffusionMaterial(c, k, U, L));
var boundaryFactory3D = new SurfaceBoundaryFactory3D(0, new ConvectionDiffusionMaterial(c, k, U, L));
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> >();
nodeBoundaries = new List <IList <Node> >();
for (int i = 0; i < 10; i++)
{
elementBoundaries.Add(new List <Element>());
nodeBoundaries.Add(new List <Node>());
}
for (int i = 0; i < 2; i++)
{
elementDomains.Add(new List <Element>());
}
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.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 <Node> nodesCollection = new List <Node>();
for (int TriID = 0; TriID < NumberOfTriElements; TriID++)
{
i++;
line = text[i].Split(delimeters);
for (int j = 0; j < (line.Length - 1); j++)
{
nodesCollection.Add(model.NodesDictionary[Int32.Parse(line[j])]);
}
IReadOnlyList <Node> nodes = new List <Node>
{
model.NodesDictionary[Int32.Parse(line[2])],
model.NodesDictionary[Int32.Parse(line[1])],
model.NodesDictionary[Int32.Parse(line[0])]
};
var Tri3 = boundaryFactory3D.CreateElement(CellType.Tri3, nodes);
var element = new Element();
element.ID = TriID;
element.ElementType = Tri3;
model.SubdomainsDictionary[0].Elements.Add(element);
model.ElementsDictionary.Add(TriID, element);
double r3 = 0;
foreach (Node node in nodes)
{
element.AddNode(node);
r3 += Math.Sqrt(Math.Pow(node.X, 2) + Math.Pow(node.Y, 2) + Math.Pow(node.Z, 2));
}
if (element.Nodes[0].X == 0 && r3 < 3 * 5e-4)
{
elementBoundaries[0].Add(model.ElementsDictionary[TriID]);
//model.NodesDictionary[element.Nodes[0].ID].Constraints.Add(new Constraint() { DOF = ThermalDof.Temperature, Amount = 10 });
}
else if (element.Nodes[0].Y == 0 && r3 < 3 * 5e-4)
{
elementBoundaries[1].Add(model.ElementsDictionary[TriID]);
}
else if (element.Nodes[0].Z == 0 && r3 < 3 * 5e-4)
{
elementBoundaries[2].Add(model.ElementsDictionary[TriID]);
}
else if (element.Nodes[0].X == 0 && r3 > 3 * 5e-4)
{
elementBoundaries[3].Add(model.ElementsDictionary[TriID]);
}
else if (element.Nodes[0].Y == 0 && r3 > 3 * 5e-4)
{
elementBoundaries[4].Add(model.ElementsDictionary[TriID]);
}
else if (r3 == 3 * 5e-4)
{
elementBoundaries[5].Add(model.ElementsDictionary[TriID]);
}
else if (element.Nodes[0].Z == 0.1 && element.Nodes[1].Z == 0.1 && element.Nodes[2].Z == 0.1)
{
elementBoundaries[6].Add(model.ElementsDictionary[TriID]);
}
else if (element.Nodes[0].Z == 0 && r3 > 3 * 5e-4)
{
elementBoundaries[7].Add(model.ElementsDictionary[TriID]);
}
else if (element.Nodes[0].Y == 0.1 && element.Nodes[1].Y == 0.1 && element.Nodes[2].Y == 0.1)
{
elementBoundaries[8].Add(model.ElementsDictionary[TriID]);
}
else if (element.Nodes[0].X == 0.1 && element.Nodes[1].X == 0.1 && element.Nodes[2].X == 0.1)
{
elementBoundaries[9].Add(model.ElementsDictionary[TriID]);
}
}
nodesCollection = nodesCollection.Distinct().ToList();
foreach (Node node in nodesCollection)
{
double r = Math.Sqrt(Math.Pow(node.X, 2) + Math.Pow(node.Y, 2) + Math.Pow(node.Z, 2));
if (node.X == 0 && node.Y <= 5e-4 && node.Z <= 5e-4)
{
nodeBoundaries[0].Add(model.NodesDictionary[node.ID]);
//model.NodesDictionary[node.ID].Constraints.Add(new Constraint() { DOF = ThermalDof.Temperature, Amount = 10 });
}
else if (node.X <= 5e-4 && node.Y == 0 && node.Z <= 5e-4)
{
nodeBoundaries[1].Add(model.NodesDictionary[node.ID]);
}
else if (node.X <= 5e-4 && node.Y <= 5e-4 && node.Z == 0)
{
nodeBoundaries[2].Add(model.NodesDictionary[node.ID]);
}
else if (node.X == 0 && r > 5e-4)
{
nodeBoundaries[3].Add(model.NodesDictionary[node.ID]);
}
else if (node.Y == 0 && r > 5e-4)
{
nodeBoundaries[4].Add(model.NodesDictionary[node.ID]);
}
else if (r == 5e-4)
{
nodeBoundaries[5].Add(model.NodesDictionary[node.ID]);
}
else if (node.Z == 0.1)
{
nodeBoundaries[6].Add(model.NodesDictionary[node.ID]);
}
else if (node.Z == 0 && r > 5e-4)
{
nodeBoundaries[7].Add(model.NodesDictionary[node.ID]);
}
else if (node.Y == 0.1)
{
nodeBoundaries[8].Add(model.NodesDictionary[node.ID]);
}
else if (node.X == 0.1)
{
nodeBoundaries[9].Add(model.NodesDictionary[node.ID]);
}
//else
//{
// model.Loads.Add(new Load() { Node = model.NodesDictionary[node.ID], DOF = ThermalDof.Temperature, Amount = 10 });
//}
}
break;
case Attributes.quad:
do
{
i++;
line = text[i].Split(delimeters);
} while (line[0] == "");
i++;
line = text[i].Split(delimeters);
NumberOfTriElements = Int32.Parse(line[0]);
i++;
IList <Node> quadNodesCollection = new List <Node>();
for (int QuadID = 0; QuadID < NumberOfTriElements; QuadID++)
{
i++;
line = text[i].Split(delimeters);
for (int j = 0; j < (line.Length - 1); j++)
{
quadNodesCollection.Add(model.NodesDictionary[Int32.Parse(line[j])]);
}
IReadOnlyList <Node> nodes = 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])]
};
var Quad = boundaryFactory3D.CreateElement(CellType.Quad4, nodes);
var element = new Element();
element.ID = QuadID;
element.ElementType = Quad;
model.SubdomainsDictionary[0].Elements.Add(element);
model.ElementsDictionary.Add(QuadID, element);
//double r3 = 0;
//foreach (Node node in nodes)
//{
// element.AddNode(node);
// r3 += Math.Sqrt(Math.Pow(node.X, 2) + Math.Pow(node.Y, 2) + Math.Pow(node.Z, 2));
//}
//if (element.Nodes[0].X == 0 && r3 < 3 * 5e-4)
//{
// elementBoundaries[0].Add(model.ElementsDictionary[QuadID]);
// //model.NodesDictionary[element.Nodes[0].ID].Constraints.Add(new Constraint() { DOF = ThermalDof.Temperature, Amount = 10 });
//}
//else if (element.Nodes[0].Y == 0 && r3 < 3 * 5e-4)
//{
// elementBoundaries[1].Add(model.ElementsDictionary[QuadID]);
//}
//else if (element.Nodes[0].Z == 0 && r3 < 3 * 5e-4)
//{
// elementBoundaries[2].Add(model.ElementsDictionary[QuadID]);
//}
//else if (element.Nodes[0].X == 0 && r3 > 3 * 5e-4)
//{
// elementBoundaries[3].Add(model.ElementsDictionary[QuadID]);
//}
//else if (element.Nodes[0].Y == 0 && r3 > 3 * 5e-4)
//{
// elementBoundaries[4].Add(model.ElementsDictionary[QuadID]);
//}
//else if (r3 == 3 * 5e-4)
//{
// elementBoundaries[5].Add(model.ElementsDictionary[QuadID]);
//}
//else if (element.Nodes[0].Z == 0.1 && element.Nodes[1].Z == 0.1 && element.Nodes[2].Z == 0.1)
//{
// elementBoundaries[6].Add(model.ElementsDictionary[QuadID]);
//}
//else if (element.Nodes[0].Z == 0 && r3 > 3 * 5e-4)
//{
// elementBoundaries[7].Add(model.ElementsDictionary[QuadID]);
//}
//else if (element.Nodes[0].Y == 0.1 && element.Nodes[1].Y == 0.1 && element.Nodes[2].Y == 0.1)
//{
// elementBoundaries[8].Add(model.ElementsDictionary[QuadID]);
//}
//else if (element.Nodes[0].X == 0.1 && element.Nodes[1].X == 0.1 && element.Nodes[2].X == 0.1)
//{
// elementBoundaries[9].Add(model.ElementsDictionary[QuadID]);
//}
}
quadNodesCollection = quadNodesCollection.Distinct().ToList();
foreach (Node node in quadNodesCollection)
{
double r = Math.Sqrt(Math.Pow(node.X, 2) + Math.Pow(node.Y, 2) + Math.Pow(node.Z, 2));
if (node.X == 0 && node.Y <= 5e-4 && node.Z <= 5e-4)
{
nodeBoundaries[0].Add(model.NodesDictionary[node.ID]);
//model.NodesDictionary[node.ID].Constraints.Add(new Constraint() { DOF = ThermalDof.Temperature, Amount = 10 });
}
else if (node.X <= 5e-4 && node.Y == 0 && node.Z <= 5e-4)
{
nodeBoundaries[1].Add(model.NodesDictionary[node.ID]);
}
else if (node.X <= 5e-4 && node.Y <= 5e-4 && node.Z == 0)
{
nodeBoundaries[2].Add(model.NodesDictionary[node.ID]);
}
else if (node.X == 0 && r > 5e-4)
{
nodeBoundaries[3].Add(model.NodesDictionary[node.ID]);
}
else if (node.Y == 0 && r > 5e-4)
{
nodeBoundaries[4].Add(model.NodesDictionary[node.ID]);
}
else if (r == 5e-4)
{
nodeBoundaries[5].Add(model.NodesDictionary[node.ID]);
}
else if (node.Z == 0.1)
{
nodeBoundaries[6].Add(model.NodesDictionary[node.ID]);
}
else if (node.Z == 0 && r > 5e-4)
{
nodeBoundaries[7].Add(model.NodesDictionary[node.ID]);
}
else if (node.Y == 0.1)
{
nodeBoundaries[8].Add(model.NodesDictionary[node.ID]);
}
else if (node.X == 0.1)
{
nodeBoundaries[9].Add(model.NodesDictionary[node.ID]);
}
//else
//{
// model.Loads.Add(new Load() { Node = model.NodesDictionary[node.ID], DOF = ThermalDof.Temperature, Amount = 10 });
//}
}
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++;
for (int TetID = NumberOfTriElements + 1; TetID < NumberOfTriElements + NumberOfTetElements; TetID++)
{
i++;
line = text[i].Split(delimeters);
IReadOnlyList <Node> nodes = 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])]
};
var Tet4 = elementFactory3D.CreateElement(CellType.Tet4, nodes);
var element = new Element();
element.ID = TetID;
element.ElementType = Tet4;
double r4 = 0;
foreach (Node node in nodes)
{
element.AddNode(node);
r4 += Math.Sqrt(Math.Pow(node.X, 2) + Math.Pow(node.Y, 2) + Math.Pow(node.Z, 2));
}
model.SubdomainsDictionary[0].Elements.Add(element);
model.ElementsDictionary.Add(TetID, element);
if (r4 < 4 * 5e-4)
{
elementDomains[0].Add(element);
}
else
{
elementDomains[1].Add(element);
}
}
name = null;
break;
}
}
return(model);
}
