private void CreateThicknessShell(int elementDegreeKsi, int elementDegreeHeta, Matrix extractionOperator,
int[] connectivity, ShellElasticMaterial2D shellMaterial, double thickness)
{
Element element = new TSplineKirchhoffLoveShellElementMaterial(elementIDCounter, null,
elementDegreeKsi, elementDegreeHeta, thickness, extractionOperator, shellMaterial)
{
ElementType = new TSplineKirchhoffLoveShellElementMaterial(elementIDCounter, null,
elementDegreeKsi, elementDegreeHeta, thickness, extractionOperator, shellMaterial)
};
foreach (var t in connectivity)
{
element.AddControlPoint(_model.ControlPointsDictionary[t]);
}
_model.ElementsDictionary.Add(elementIDCounter++, element);
_model.PatchesDictionary[0].Elements.Add(element);
}
/// <summary>
/// Create Model from reading an .iga file.
/// </summary>
/// <param name="shellType"><see cref="Enum"/> that specifies the type of the shells that will be generated.</param>
/// <param name="shellMaterial">The material of the shell.</param>
/// <param name="thickness">The tchickness of the shell.</param>
public void CreateTSplineShellsModelFromFile(TSplineShellType shellType = TSplineShellType.Linear, ShellElasticMaterial2D shellMaterial = null, double thickness = 1)
{
char[] delimeters = { ' ', '=', '\t' };
Attributes?name = null;
String[] text = System.IO.File.ReadAllLines(_filename);
_model.PatchesDictionary.Add(0, new Patch());
for (int i = 0; i < text.Length; i++)
{
var line = text[i].Split(delimeters, StringSplitOptions.RemoveEmptyEntries);
if (line.Length == 0)
{
continue;
}
try
{
name = (Attributes)Enum.Parse(typeof(Attributes), line[0].ToLower());
}
catch (Exception exception)
{
throw new KeyNotFoundException($"Variable name {line[0]} is not found. {exception.Message}");
}
switch (name)
{
case Attributes.type:
Types type;
try
{
type = (Types)Enum.Parse(typeof(Types), line[1].ToLower());
}
catch (Exception exception)
{
throw new KeyNotFoundException($"Variable name {line[0]} is not found. {exception.Message}");
}
numberOfDimensions = type == Types.plane ? 2 : 3;
break;
case Attributes.noden:
break;
case Attributes.elemn:
var numberOfElements = int.Parse(line[1]);
break;
case Attributes.node:
var controlPoint = new ControlPoint
{
ID = controlPointIDcounter,
X = double.Parse(line[1], CultureInfo.InvariantCulture),
Y = double.Parse(line[2], CultureInfo.InvariantCulture),
Z = double.Parse(line[3], CultureInfo.InvariantCulture),
WeightFactor = double.Parse(line[4], CultureInfo.InvariantCulture)
};
_model.ControlPointsDictionary.Add(controlPointIDcounter, controlPoint);
((List <ControlPoint>)_model.PatchesDictionary[0].ControlPoints).Add(controlPoint);
controlPointIDcounter++;
break;
case Attributes.belem:
var numberOfElementNodes = int.Parse(line[1]);
var elementDegreeKsi = int.Parse(line[2]);
var elementDegreeHeta = int.Parse(line[3]);
i++;
line = text[i].Split(delimeters);
int[] connectivity = new int[numberOfElementNodes];
for (int j = 0; j < numberOfElementNodes; j++)
{
connectivity[j] = Int32.Parse(line[j]);
}
var extractionOperator = Matrix.CreateZero(numberOfElementNodes,
(elementDegreeKsi + 1) * (elementDegreeHeta + 1));
for (int j = 0; j < numberOfElementNodes; j++)
{
line = text[++i].Split(delimeters);
for (int k = 0; k < (elementDegreeKsi + 1) * (elementDegreeHeta + 1); k++)
{
extractionOperator[j, k] = double.Parse(line[k]);
}
}
if (numberOfDimensions == 2)
{
throw new NotImplementedException("TSpline2D not yet implemented");
}
else
{
switch (shellType)
{
case TSplineShellType.Linear:
CreateLinearShell(elementDegreeKsi, elementDegreeHeta, extractionOperator, connectivity);
break;
case TSplineShellType.Thickness:
CreateThicknessShell(elementDegreeKsi, elementDegreeHeta, extractionOperator, connectivity, shellMaterial, thickness);
break;
}
}
break;
case Attributes.set:
break;
}
}
}
