public ODEElement3D(IReadOnlyList <Node> nodes, IIsoparametricInterpolation3D interpolation,
IQuadrature3D quadratureForStiffness, IQuadrature3D quadratureForMass,
IGaussPointExtrapolation3D gaussPointExtrapolation, ODEMaterial material)
{
this.material = material;
this.GaussPointExtrapolation = gaussPointExtrapolation;
this.Nodes = nodes;
this.Interpolation = interpolation;
this.QuadratureForConsistentMass = quadratureForMass;
this.QuadratureForStiffness = quadratureForStiffness;
dofTypes = new IDofType[nodes.Count][];
for (int i = 0; i < interpolation.NumFunctions; ++i)
{
dofTypes[i] = new IDofType[] { ThermalDof.Temperature }
}
;
}
public ODEDomainLoad(ODEMaterial material, double load, IDofType dofType)
{
_material = material;
_load = load;
_dofType = dofType;
}
public Model CreateModelFromFile()
{
var ODEMaterial = new ODEMaterial(C1, C0);
var elementFactory3D = new ODEElement3DFactory(ODEMaterial);
//var elementFactory3D = new ConvectionDiffusionElement3DFactory(new ConvectionDiffusionMaterial(diffusionCoeff, convectionCoeff, loadFromUnknownCoeff));
//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> >();
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>());
}
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]);
//int elementBoundaryID = Int32.Parse(line[0]);
//elementBoundaries[elementBoundaryID].Add(model.ElementsDictionary[QuadID]);
}
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]);
//int elementBoundaryID = Int32.Parse(line[0]);
//elementBoundaries[elementBoundaryID].Add(model.ElementsDictionary[QuadID]);
}
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 = 0; TetID < 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;
foreach (Node node in nodes)
{
element.AddNode(node);
}
model.SubdomainsDictionary[0].Elements.Add(element);
model.ElementsDictionary.Add(TetID, element);
}
i = i + 3;
for (int TetID = 0; TetID < NumberOfTetElements; TetID++)
{
i++;
line = text[i].Split(delimeters);
int elementDomainID = Int32.Parse(line[0]);
elementDomains[elementDomainID - 1].Add(model.ElementsDictionary[TetID]);
//int elementBoundaryID = Int32.Parse(line[0]);
//elementBoundaries[elementBoundaryID].Add(model.ElementsDictionary[QuadID]);
}
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++;
for (int HexID = 0; HexID < NumberOfHexElements; HexID++)
{
i++;
line = text[i].Split(delimeters);
IReadOnlyList <Node> nodes = 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])],
};
var Hexa8 = elementFactory3D.CreateElement(CellType.Hexa8, nodes);
var element = new Element();
element.ID = HexID;
element.ElementType = Hexa8;
foreach (Node node in nodes)
{
element.AddNode(node);
}
model.SubdomainsDictionary[0].Elements.Add(element);
model.ElementsDictionary.Add(HexID, element);
}
i = i + 3;
for (int HexID = 0; HexID < NumberOfHexElements; HexID++)
{
i++;
line = text[i].Split(delimeters);
int elementDomainID = Int32.Parse(line[0]);
elementDomains[elementDomainID - 1].Add(model.ElementsDictionary[HexID]);
}
break;
}
}
return(model);
}
public ODEElement3DFactory(ODEMaterial commonMaterial)
{
this.commonMaterial = commonMaterial;
}
